Next Article in Journal
Evaluation of the Accuracy of Spectral Calibration Light Source on Spectral Radiance Acquired by the Greenhouse-Gases Absorption Spectrometer-2 (GAS-2)
Next Article in Special Issue
Preliminary Results of the Three-Dimensional Plasma Drift Velocity at East Asian Low-Latitudes Observed by the Sanya Incoherent Scattering Radar (SYISR)
Previous Article in Journal
Editorial on New Challenges in Solar Radiation, Modeling and Remote Sensing
Previous Article in Special Issue
Wavenumber-4 Structure in COSMIC-2 Observations: Vertical Plane Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Method of Extracting Plasma Lines Based on Sanya Incoherent Scatter Radar

1
Key Laboratory of Earth and Planetary Physics, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2
College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
3
Beijing National Observatory of Space Environment, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(10), 2634; https://doi.org/10.3390/rs15102634
Submission received: 18 April 2023 / Revised: 15 May 2023 / Accepted: 16 May 2023 / Published: 18 May 2023

Abstract

:
The plasma lines observed by Sanya incoherent scatter radar (SYISR) are dependent on the enhancement of Langmuir waves due to superthermal photoelectrons generated by solar EUV radiation. The plasma line power spectrum can be obtained using long-pulse and alternating-code transmission signals during the period from sunrise to noon almost every day. For the power spectrum of the long pulse, the CLEAN algorithm that has been applied in this field is used to verify the feasibility of this method for SYISR in only a few cases. However, it is difficult to deal with alternating code with such a low SNR using the general deconvolution method. The irreversible-migration filtering (IMF) method has been developed to separate signal noise from the measurements of the alternating code. Some experimental results from the SYISR measurements validate the excellent performance of the IMF method for alternating code. Additionally, an example observation of the electron density with a high time and range resolution is derived. The results show that plasma line detection can be a powerful new observational capability for SYISR as an ionospheric experimental mode for ionospheric calibration, when possible, which can be simultaneously measured with the ion line for constant radar calibration in the standard fitting of the ion line.

1. Introduction

Incoherent scatter radar is used for measuring the properties of ionospheric plasma by receiving the radar backscatter echo signal modulated by the thermal fluctuation of electrons and ions in the ionosphere. Following the contribution of electrons and ions to the fluctuation of electron density, the incoherent scatter spectrum can be divided into the ion line near the zero frequency and the plasma line near the plasma frequency. Its spectral peak is related to the solution of the dispersion relationship in cold, magnetized plasma, namely, the ion acoustic wave and the Langmuir wave.
Because of strong echo power, the ion lines with a classical double-humped spectra shape are relatively easy to routinely detect, so they provide a wealth of information for obtaining the ionospheric parameter estimates of electron density, electron temperature, ion temperature and plasma drift velocity. Compared with ion lines, the echo signal amplitude of plasma lines is very weak, lower than the measured power spectrum noise level, which makes it difficult for most ISRs to detect the echo signal of plasma lines.
In the presence of suprathermal photoelectrons, the velocity of the superthermal charged particles matches the phase velocity of the Langmuir wave, the amplitude of the Langmuir wave can be enhanced, and the plasma lines near the plasma frequency can also be correspondingly enhanced and detected by ISR [1]. The first observation of plasma lines enhanced by the photoelectrons produced by solar UV radiation in the daytime was successfully performed with the Arecibo Observatory radar [2]. In the absence of direct solar EUV during the night, plasma lines enhanced by the injection of energetic particles into the ionosphere, such as auroral electron precipitation, can also be detected by EISCAT UHF and Sondrestrom ISR [3,4]. In addition, the technology of the artificially enhanced plasma line has also been developed by injecting high-power, high-frequency electromagnetic waves into the ionosphere through ground heating facilities. The Arecibo ISR [5,6,7] and Tromsø ISR [8,9] radar observed the enhanced plasma line after the artificial heating of the ionosphere.
The ionospheric electron density with a high time and range resolution can be determined based on the offset frequency of the plasma line from the transmitter frequency, which so far is the most accurate ground measurement method to obtain the absolute electron density [4,10,11]. In addition, plasma lines appear in pairs on both sides of the ion lines, which are regarded as up- and down-shifted plasma lines, respectively. If the plasma line with a high spectral resolution can be measured, then the high-precision electron temperature can be obtained by determining the asymmetrical frequency difference between the up- and down-shifted plasma lines [12,13,14]. In the study of Bjørnå and Kirkwood [15], the electron density and electron temperature measured by the plasma line were used as the prior knowledge of the ion line fitting to reduce the number of fitting parameters and eliminate the ambiguity of temperature–ion composition. Moreover, the offset frequency of the plasma lines was added into the ion line least-square fitting, and the ionospheric ion composition was successfully extracted from the measured data obtained by the EISCAT UHF and VHF radars [16]. Similarly, if the measurement of the plasma line is combined in the ion line fitting, then the ion neutral collision frequency in the lower height E region can also be obtained [17].
Although plasma lines have been studied by many scholars for a long time, the work on plasma lines in the low latitudes of East Asia has not yet been carried out. Sanya incoherent scattering radar (SYISR) is China’s first phased-array incoherent scatter radar located in Sanya (18.34°N, 109.62°E). Preliminary ion line results and evaluations are presented in Hao et al. [18]. A prominent capability of SYISR is the ability to detect weak plasma line spectra, except for ion lines, in the regular mode. In this paper, we focus on the natural enhanced plasma line power spectrum during the daytime measured by a long pulse and alternating code from SYISR. The signal-to-noise ratio of the echo from the long pulse is much stronger than that from the alternating code. So, we use the CLEAN method with iterative automatic extraction to process long-pulse data and also use a novel method, irreversible-migration filtering (IMF) for alternating-code data by mainly taking advantage of migration and de-migration. In Section 2, we describe the experimental details for obtaining the measured plasma power spectrum. We apply the CLEAN and IMF algorithms to extract the plasma frequency in Section 3. Further, the electron density estimation from the plasma lines is presented. The applicable requirements for the two different methods are summarized in Section 4.

2. Experimental Setup

The plasma line echo signal is the backscattered signal from the high-frequency Langmuir wave in the scatterer volume, whose Doppler shifts from transmitter frequency are approximately equal to the plasma frequency. In most cases, the plasma frequency near Sanya is in the range of several MHz to twenty MHz. Therefore, in order to detect the complete plasma line, SYISR adopts three channels for simultaneous digital down-conversion and filtering. Different digital oscillators for the three channels are presented according to the plasma frequency requirements. The local oscillation frequency of one channel is used for both transmission and reception in ion line measurement, and the other two only for different receiver channels in plasma line measurement. Each channel obtains the power spectrum at a different frequency band, and the spectra at these three frequency bands are further synthesized into the full incoherent scatter spectrum, including the ion line and plasma lines. The digital filter bandwidth and output data rate for each channel are 4 MHz, so the maximum coverage of the spectrum that can be detected at once during the plasma line experiments is 12 MHz. If the plasma line is to be continuously observed for a long time, then the digital local oscillator frequency of the plasma line channel needs to be flexibly adjusted, so that the receiver frequency covered by the filter window exactly matches the plasma frequency.
Currently, coded and uncoded long-pulse transmission schemes have been developed for incoherent scatter radar to obtain the measurements of plasma lines [10,19]. Over SYISR, the 430~450 MHz radar frequency is operated at a peak power near 2 MW, and plasma lines enhanced by the natural superthermal electrons can be routinely detected during the daytime with an uncoded long pulse and alternating code. In order to ensure the maximum antenna gain and a better signal-to-noise ratio (SNR) of the plasma line echo, the radar beams in the following experimental results are all pointed vertically. The processing of the observations is offline. For the same code type, the solution of the plasma line power spectrum is the same as that of the ion line spectrum, except that more data samples need to be calculated on a range gate and it is more time-consuming. The examples of the full incoherent scattering spectrum of the ion line and up- and downshifted plasma lines detected at noon are shown in Figure 1. The left panel shows the spectrum measured by 32-bit alternating code with a 20 μs baud length. The right panel represents the spectrum measured by a long pulse with a 100 μs pulse length.
The results spectra displayed in Figure 1 are the spliced power spectra of the three respective channels after removing the background. Due to a different channel background noise for each channel, the intensity shown in the combined power spectrum is not the real relative intensity of the ion line and plasma lines. Here, it is mainly demonstrated that the ion line and plasma lines can be simultaneously obtained by a long pulse or alternating code. The range resolution of the 100 μs long pulse is theoretically 15 km so that the plasma lines can be visibly detected. Then, the corresponding frequency resolution is 5 kHz. The alternating code has a higher range resolution of 3 km and a frequency resolution of 0.781 kHz. Note that we use a 3 km range step for two radar transmit codes in Figure 1, but this measurement of the long pulse is in fact clearly smeared in this range. Based on the characteristic that the plasma line has a strong frequency dependence, alternating code is used to determine the plasma frequency with good accuracy.

3. Method

The obtained echo power spectrum can be regarded as the two-dimensional convolution of the spectral ambiguity function and the incoherent scatter power spectrum in terms of frequency and range. Due to its detectability depending on the intensity of the SNR, the intensity of the decoded plasma line power spectrum is very noisy under the lower SNR condition. In addition, it may also be affected by the response of the radar receiver, which may cause the calculated plasma line power spectrum to contain some unwanted signal contributions, referred to as background noise signal. As shown in Figure 1, there exists obvious convolutional distortion in the range of the long-pulse power spectrum. In comparison, the range resolution of the alternating code is much better, but the SNR is much worse. Therefore, the following different methods are used to extract plasma line profiles as a function of altitude and frequency for the long pulse and alternating code, respectively.

3.1. Deconvolution for Long-Pulse Spectrogram

In order to recover the true spectrogram from the convolutional distortion, the well-known CLEAN algorithm [20] is utilized, which is widely used to reconstruct images in radio astronomy. The algorithm is effective against point sources, that is, based on the fact that the map to be deconvolved is itself a convolution of multiple point sources and the point spread function (PSF) that describes the features of the radar receiver. Then, the measured plasma line spectrogram can be regarded as a dirty map to be deconvoluted, the ambiguity function generated from the radar transmit waveform and the receiver impulse response as the PSF, and the plasma frequency at each range gate as a point source. Therefore, the CLEAN algorithm can be used to deconvolve the plasma line spectrogram. It has been proven that this method can remove range-Doppler smearing for the long-pulse measurement with a sufficient SNR, and thus extract the plasma line profiles [4,21].
The CLEAN algorithm is implemented by repeatedly subtracting the dirty beam from the dirty map until the maximum intensity of the residual map is less than the specified noise level or reaches a preset number of iterations, thus removing the smearing in the range. The PSF required is the receiver response to a point source in plasma, and its accurate formation is very important. To obtain complete information, the PSF should be constructed by the real transmitted radar sampling. This allows a more realistic model of the ambiguity of the plasma echo at each height in the range direction. However, the receiver of SYISR does not support the sampling of transmitted signals. Here, we instead use an idealized rectangular pulse as the transmission waveform and receiver pulse response, respectively, to create a matrix of the PSF where the transmitted pulse sampling is at the center. It is achieved by multiplying the transmitted pulse sampling by the shifted version of its complex conjugate signal, where the shifted operation is performed at one sampling interval each time. Then, the non-overlapping part of each row of the matrix undergoes zero-padding, so that the total length of the product is consistent with the transmitted pulse. In this study, we performed the deconvolution procedures as described above for the long-pulse plasma line spectrogram.
Figure 2(left) shows the distorted upshifted plasma line spectrum measured by a long pulse at noon. At this time, the electron density is very high, and the plasma line spectrum with a sufficiently strong SNR can be obtained as a detectable and recognizable feature. The power spectrum of the echo signal with a 100 μs transmitting pulse width are obtained by averaging 36,037 pulses with an integration time of 2 min, so that a thin plasma line can be indistinctly seen, even though it is smeared by about 15 km in the altitude direction. In this case, we can use the CLEAN algorithm to remove range smearing; otherwise, this method is invalid for recovering the true plasma line profiles from the convolutional distortion under a very low SNR. To better estimate the plasma line frequency, it is necessary to remove the background noise from the deconvolved spectrum. We select the median of the spectrum at the last few heights as the background noise, and then subtract it from the spectrum after deconvolution. The result of the deconvolved plasma line spectrum with the background noise removed after 300,000 iterations is shown in Figure 2 (right). It can be seen that the plasma line is extracted below 200 km, but not above 200 km because of contamination from the strong background noise level. The parts that can be recognized by the human eye are actually discrete plasma frequency points.
This method has also been used to handle plasma line observations with lower SNRs at other local times. Figure 3 shows an example of the detectable plasma line profile. Although the signal intensity near the maximum electron density is strong enough to clearly exhibit a distorted plasma line profile, it also includes strong background noise, causing the true plasma line to be completely submerged in noise. That is to say, the point source strength is very weak compared to the extended structures in the range. Under such circumstances, convergence instability in the CLEAN process may occur when the PSF is inaccurately generated with limited information. With more iterations or a smaller loop gain, it is also impossible to find the location of the true plasma line. Therefore, larger error bars are produced after deconvolution, and we do not show the deconvolved result here.

3.2. Separating Signal–Noise for Alternating-Code Spectrogram

When transmitting pulses with a longer pulse width for sufficient radar efficiency, alternating code can be used to achieve a higher range resolution of plasma line measurements, up to one baud width. However, there have been difficulties with using general deconvolution methods to extract the plasma line profiles of alternating code. We use an irreversible-migration filtering (IMF) method [22] to perform signal–noise separation, which can achieve the purpose of background noise removal and SNR improvement. Then, some maximum-intensity points on some range gates are found, which are the plasma frequency values at the corresponding heights.
To illustrate the effectiveness of this method of plasma line extraction, we take the plasma line echo data of 64-bit alternating code with a 10 μs code width as an example. The sampling interval is 0.25 μs. For the alternating code with a 640 μs pulse width, the maximum lag number of the decoded autocorrelation function (ACF) is 2560. In order to improve the SNR of the plasma line power spectrum, we utilize the first 256 complex-valued lags of the ACF with a Hamming window. The resulting frequency resolution is 7.8125 kHz. Figure 4 shows the raw upshifted plasma line power spectrum contaminated by background noise on 8 March 2021. The range step is 10 μs, which results in a 1.5 km altitude resolution, and the integration time is 2 min. Note that there exist some obvious vertical stripe background signals, but a strong intensity of the plasma line at the peak of the F layer is obtained in Figure 4. Thus, under such strong background noise, it is very difficult to obtain the plasma frequency at other range gates, except the peak of the F layer.
The IMF method is predominantly employed for wavefield separation in seismic data processing and has demonstrated its efficacy in various applications, including diffraction separation, structural noise elimination, and surface wave separation. The underlying principle of this technique involves partitioning the migration domain into evanescent and propagating regions; wavefield components residing within the evanescent region undergo automatic attenuation, whereas those within the propagating region are inherently preserved. Consequently, adjusting the migration velocities enables the separation of wavefield components with distinct slopes throughout the migration and de-migration procedures. Given the dissimilar slopes of the background noise and signals, this approach can also be extended to signal–noise separation within radar signal processing. Hence, we can convert the signal in f–k space into the migrated domain as [22,23]:
A ( k , ω , z ) = e j z 4 ω 2 v 2 k 2 A ( k , ω , 0 )
where A ( k , ω , 0 ) is the signal power spectrum, z is the transformed travel time, v is the migration velocity, k is the wavenumber, and ω = 2 π f is the circular frequency of the echo signal. The exponential term in Equation (1) contributes to the phase change when 4 ω 2 / v 2 k 2 > 0 ; on the contrary, it is an attenuation term and the signals are migrated to the evanescent region during de-migration, when 4 ω 2 / v 2 k 2 < 0 . To do this, we tune the migration velocity, although it might not have any obvious physical meaning, so that the background noise with large absolute slope values is irreversible after de-migration, as far as possible. However, the other signals are returned to the original region. Thus, the background signal can be separated from the raw plasma line power spectrum.
To eliminate the background noise more effectively with the IMF method, we need to distinguish the data near the peak of the F layer from the data at the other heights. As can be seen from Figure 4, compared with other range gates, the slope value at the peak of the F layer is relatively large, close to vertical. For example, if we choose a larger migration velocity, then the raw effective plasma line will fall in the irreversible region in the migrated domain, and vice versa. In other words, this method may mistakenly remove the signal at the peak of the F layer as the background if an inappropriate migration velocity is used. Hence, the crucial feature of the IMF method is to choose an optimal velocity. Figure 5 shows an example of signal–noise separation from the raw power spectrum in Figure 4 using the IMF method. The power spectrum results after background removal are shown in Figure 5a,c,e. Figure 5b,d,f shows the separated background noise power spectrum. For range gates less than 200 km and more than 250 km, we use the same relatively large migration velocity (e.g., v = 8.5 m/s), while the smaller migration velocity (e.g., v = 0.8 m/s) is used between 200 and 250 km. It can be seen that the vertical stripe background has been removed so well that we can find the plasma frequency more accurately at each altitude in the subsequent processing.
The white circles in Figure 5a,c,e indicate the plasma frequency extracted at some heights where the plasma line is enhanced. The plasma frequency estimation mainly includes three steps. First, the frequency value corresponding to the maximum intensity at each height is obtained. Second, the second-order difference values of the plasma frequency between different heights are calculated. Finally, these difference values are combined with the threshold value to eliminate the falsely selected points greater than the threshold value. An appropriate threshold value should be selected against the background signal. A constant threshold value is usually enough, even for complex data. For the data of the example discussed here, we set the threshold value as 0.03.
Figure 6 shows sections of the power spectrum before and after using the IMF method at various altitudes. The data in Figure 6a,c,e are all from the raw power spectrum in Figure 4. The data in Figure 6b,d,f are selected from the results after using the IMF method in Figure 5a,c,e respectively. The red circles indicate the true plasma frequency position at the corresponding height. The green circles represent obvious strong background signals, some of which are even greater than the intensity of the true plasma line, as shown in Figure 6a,e. Note that the desired plasma frequency values are well-retained without damaging the valid information, and the unwanted strong signals are clearly removed. Consequently, the signal after background removal by this method is also the result of filtering and smoothing the input signal.

3.3. Plasma Line Profile Fitting

In most cases, only a limited number of plasma frequency values can be obtained from the observed plasma line power spectra. This is not good for reproducing a plasma line profile containing the complete altitude information. Based on the finite number of effective plasma frequency values extracted in the previous process, the α-Chapman function is used to fit the change of the enhanced plasma line to the height. Here, we use the varied scale height to fit the profile of the plasma line that is similar to the fitted form of the electron density profiles used in [24]
{ f p ( h ) = f o F 2 e x p ( 0.5 ( 1 z e z ) ) z = h h m F 2 H m + α ( h h m F 2 )
where h m F 2 and f o F 2 are the F2 layer peak height and critical frequency, respectively. H m is the scale height at the F2 layer peak, α is the ratio of altitude variations of the scale height. In those four parameters, the peak height and density are known. Thus, f o F 2 , h m F 2 , H m and α can be determined by using the least-squares fitting approach.
The comparison of the fitted results with the raw enhanced plasma lines from Figure 4 is shown in Figure 7 (left). It can be seen that the white dotted line basically coincides with the enhanced plasma spectral line above 220 km, and the fitting result is the best between 220 km and 250 km. The deviation below 220 km may be caused by the lack of effective plasma frequency points for the fitting between 180–220 km due to a poor SNR.
The power spectrum containing the background signal as white noise is shown in Figure 7 (right). Compared with the results in Figure 6 (left), there is no obvious large slope for the background signal. The IMF method can also deal with this situation. If a migration velocity is selected, then the smoothing effect of the noisy power spectrum can be achieved, which also improves the SNR of the plasma line power spectrum. Further, the plasma frequency extraction method described above is used to accurately extract the plasma line profile.

3.4. Computation of Plasma Electronic Density

Based on the incoherent scatter theory of magnetized plasma, the linear dispersion relationship of the Langmuir wave with a magnetic field is expressed as [25]
f r ± 2 = f p e 2 + 3 k ± 2 4 π 2 k b T e m e + f c e 2 s i n 2 α
where f r ± is the plasma frequency, f p e = n e e 2 / 4 π 2 ε 0 m e is the electron plasma frequency, which is also called Langmuir frequency, f c e = e B / 4 π 2 m e is the electron cyclotron frequency, k b is the Boltzmann constant, T e is the electron temperature, m e is electron mass, k ± is the scatter wave number, and α is the angle between the backscattered wave vector and the magnetic field. For SYISR, the angle between the zenith direction and the magnetic field is about 65° in the F region over Sanya. Due to the main dependence of the plasma frequency on the Langmuir frequency, which is only related to the electron density and some fundamental constants, the inaccuracies associated with using the simplified Langmuir dispersion relation are usually negligible when determining the electron density from plasma line measurements [26]. Simplistically, the derived electron density can be approximately expressed as n e = e 2 4 π 2 ε 0 m e f r 2 . Therefore, if the plasma frequency can be accurately recovered, then the absolute electron density profile can be calculated.

4. Discussion and Conclusions

In this paper, the methods of extracting the plasma frequency with a long pulse and alternating code on the basis of SYISR plasma line experiment data are introduced. Each method has its own advantages and limitations. The CLEAN algorithm does not have the ability to deal with weaker signals under the condition of a lower SNR, since the algorithm is largely dependent on the intensity of the true signal. In other words, the feature of the true signal to be deconvoluted must be clearly identified by human experts and is merely smeared through the convolution process, but without losing its inherent characteristics. For SYISR, this method is used for the plasma line spectrum, which is greatly enhanced by photoelectrons produced by solar UV irradiation over a very short period at noon, when the electron density is at its highest. In addition, it is a time-consuming calculation task to use the highly iterative algorithm to remove the smearing and background noise and obtain the true plasma line profiles. Consequently, the CLEAN algorithm is only suitable for long-pulse deconvolution with a sufficiently strong SNR, not for alternating code with a high range resolution but weak intensity.
The proposed IMF method describes a non-iterative algorithm for separating the signal–noise of the alternating-code spectrogram, and we have conducted a large amount of statistical analyses using this method and presented some examples on SYISR data in this paper to verify its capability. Although the SNR of the alternating-code echo is worse than that of the long pulse, the obtained power spectrum has the great advantage that the features of the plasma line profile are prominent due to the small ambiguity in the range direction. As long as the appropriate migration velocity is selected for the plasma line spectrogram with obvious characteristics, the background noise can be removed well and the plasma line profiles can be obtained after fitting. Further, the electron density profiles can be obtained by the detailed method outlined above for processing such alternating-code data.
A comparison of the electron density derived from the ion line, plasma line and ionosonde analyses is shown in Figure 8. Figure 8a,b uses the same alternating-code echo data and the integration time of 2 min. The electron density in Figure 8a is derived from the standard ion line analysis and is the result of further calibration with Figure 8c after fitting. The result of Figure 8b is derived from the plasma line analysis with the IMF method. We can see that the magnitude of the electron density measured by the ion line is not very smooth in terms of time and height compared with that by the plasma line. This is because there is parameter ambiguity in the multiparameter fitting of the ion line, where the incoherent scatter model is a highly nonlinear function of its plasma parameters. The high-time and range resolution of the electron density obtained by the plasma line measurement can also be used as a calibration method for the electron density of the ion line. For the plasma line measurement, the SNR only needs to be sufficient to obtain the corresponding plasma frequency, and then the electron density with a high range resolution from the alternating-code data can be directly obtained.
The plasma line technique is the most accurate approach to obtain the absolute electron density with a high time and range resolution, but it also has some drawbacks. For the present SYISR method, only enhanced plasma line profiles above and below the peak height can be detected, such as from sunrise to noon, and then the complete electron density profile can be derived through the IMF method. Hence, its detectability largely depends on the radar power and receiver sensitivity, as well as the local time. Especially, the developing SYISR Tristatic System will double the transmitter power of the current SYISR method. It is expected that the simultaneous measurement of the ion line and plasma line can be routinely carried out when the radar energy is strong enough. In this case, the electron density measured by the plasma line can be used for the real-time calibration of a constant radar system instead of ionosonde. Considering that a large number of effective plasma line spectra will be obtained in the future, we will also develop a deep-learning approach to achieve the automatic recognition and extraction of plasma lines. Then, the IMF algorithm will be used as a potential method to extract plasma lines with different scales of the ionosphere structure as its training dataset, which also solves the manual scaling problem in the processing of deep learning. This is beneficial for the real-time efficiency of obtaining calibration factors that vary with local time and altitude [27]. In addition, due to the superiority of SYISR’s geographical location, the radar scanning detection angle can almost meet the requirements of any angle of the magnetic field. Therefore, while conducting plasma line detection, plasma line splitting may be found from a lower-SNR echo when possible, which will greatly improve the global contribution of SYISR to space plasma experiments and basic plasma scatter theory.

Author Contributions

Conceptualization, H.H.; methodology, H.H.; formal analysis, H.H.; investigation, H.H. and B.Z.; software, H.H.; validation, H.H.; data curation, H.H.; visualization, H.H.; supervision, B.Z. and X.Y.; project administration, X.Y.; funding acquisition, B.Z. writing—original draft preparation, H.H.; writing—review and editing, B.Z., X.Y., F.D., B.N. and L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (42174206). The ionogram data were provided by Beijing National Observatory of Space Environment, Institute of Geology and Geophysics Chinese Academy of Sciences through the Geophysics Center, National Earth System Science Data Center (http://wdc.geophys.ac.cn, accessed on 17 April 2023). This research was funded by the Young Talent Program of Chinese Academy of Science (CAS) and Youth Innovation Promotion Association CAS. This work was also supported by Stable-Support Scientific Project of China Research Institute of Radiowave Propaga-tion (Grant No. A13XXXXWXX). The authors acknowledge the use of data from the Chinese Meridian Project.

Data Availability Statement

The data presented in this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Perkins, F.; Salpeter, E.E. Enhancement of Plasma Density Fluctuations by Nonthermal Electrons. Phys. Rev. 1965, 139, A55–A62. [Google Scholar] [CrossRef]
  2. Perkins, F.W.; Salpeter, E.E.; Yngvesson, K.O. Incoherent Scatter from Plasma Oscillations in the Ionosphere. Phys. Rev. Lett. 1965, 14, 579–581. [Google Scholar] [CrossRef]
  3. Haggstrom, I.; Hedin, M.; Aso, T.; Pellinen-Wannberg, A.; Westman, A. Auroral field-aligned currents by incoherent scatter plasma line observations in the E region. Adv. Polar Up. Atmos. Res. 2000, 14, 103–121. [Google Scholar]
  4. Vierinen, J.; Bhatt, A.; Hirsch, M.A.; Stromme, A.; Semeter, J.L.; Zhang, S.R.; Erickson, P.J. High temporal resolution observations of auroral electron density using superthermal electron enhancement of Langmuir waves. Geophys. Res. Lett. 2016, 43, 5979–5987. [Google Scholar] [CrossRef]
  5. Carlson, H.C.; Djuth, F.T.; Zhang, L.D. Creating space plasma from the ground. J. Geophys. Res. Space Phys. 2017, 122, 978–999. [Google Scholar] [CrossRef]
  6. Carlson, H.C.; Wickwar, V.B.; Mantas, G.P. Observations of fluxes of suprathermal electrons accelerated by HF excited instabilities. J. Atmos. Terr. Phys. 1982, 44, 1089–1100. [Google Scholar] [CrossRef]
  7. Muldrew, D.B.; Showen, R.L. Height of the HF-enhanced plasma line at Arecibo. J. Geophys. Res. 1977, 82, 4793–4804. [Google Scholar] [CrossRef]
  8. Rietveld, M.T.; Isham, B.; Kohl, H.; La Hoz, C.; Hagfors, T. Measurements of HF-enhanced plasma and ion lines at EISCAT with high-altitude resolution. J. Geophys. Res. Space Phys. 2000, 105, 7429–7439. [Google Scholar] [CrossRef]
  9. Rietveld, M.T.; Kohl, H.; Kopka, H.; Stubbe, P. Introduction to ionospheric heating at Tromsø—I. Experimental overview. J. Atmos. Terr. Phys. 1993, 55, 577–599. [Google Scholar] [CrossRef]
  10. Djuth, F.T.; Sulzer, M.P.; Elder, J.H. Application of the Coded Long-Pulse Technique to Plasma Line Studies of the Ionosphere. Geophys. Res. Lett. 1994, 21, 2725–2728. [Google Scholar] [CrossRef]
  11. Hagfors, T.; Birkmayer, W.; Sulzer, M. A new method for accurate ionospheric electron density measurements by incoherent scatter radar. J. Geophys. Res. Space Phys. 1984, 89, 6841–6845. [Google Scholar] [CrossRef]
  12. Hagfors, T.; Lehtinen, M.S. Electron temperature derived from incoherent scatter radar observations of the plasma line frequency. J. Geophys. Res. 1981, 86, 119–124. [Google Scholar] [CrossRef]
  13. Kofman, W.; Lejeune, G.; Hagfors, T.; Bauer, P. Electron temperature measurements by the plasma line technique at the French Incoherent Scatter Radar Facilities. J. Geophys. Res. Space Phys. 1981, 86, 6795–6801. [Google Scholar] [CrossRef]
  14. Nicolls, M.J.; Sulzer, M.P.; Aponte, N.; Seal, R.; Nikoukar, R.; González, S.A. High-resolution electron temperature measurements using the plasma line asymmetry. Geophys. Res. Lett. 2006, 33, L18107. [Google Scholar] [CrossRef]
  15. Bjørnå, N.; Kirkwood, S. Derivation of ion composition from a combined ion line/plasma line incoherent scatter experiment. J. Geophys. Res. 1988, 93, 5787–5793. [Google Scholar] [CrossRef]
  16. Fredriksen, Å. High latitude quiet summer ion composition profiles derived from a combined ion line/plasma line incoherent scatter experiment. J. Atmos. Terr. Phys. 1990, 52, 77–84. [Google Scholar] [CrossRef]
  17. Bjørnå, N. Derivation of ion-neutral collision frequencies from a Combined Ion Line/Plasma Line Incoherent Scatter Experiment. J. Geophys. Res. 1989, 94, 3799–3804. [Google Scholar] [CrossRef]
  18. Hao, H.; Zhao, B.; Wan, W.; Yue, X.; Ding, F.; Ning, B.; Zeng, L.; Jin, Y.; Wang, J.; Zhang, N. Initial Ionospheric Ion Line Results and Evaluation by Sanya Incoherent Scatter Radar (SYISR). J. Geophys. Res. Space Phys. 2022, 127, e2022JA030563. [Google Scholar] [CrossRef]
  19. Guio, P.; Bjørnå, N.; Kofman, W. Alternating-code experiment for plasma-line studies. Ann. Geophys. 1997, 14, 1473–1479. [Google Scholar] [CrossRef]
  20. Högbom, J.A. Aperture synthesis with a non-regular distribution of interferometer baselines. Astron. Astrophys. Suppl. Ser. 1974, 500, 417–426. [Google Scholar]
  21. Vierinen, J.; Gustavsson, B.; Hysell, D.L.; Sulzer, M.P.; Perillat, P.; Kudeki, E. Radar observations of thermal plasma oscillations in the ionosphere. Geophys. Res. Lett. 2017, 44, 5301–5307. [Google Scholar] [CrossRef]
  22. Li, C.; Zhang, J.H. Wavefield separation using irreversible-migration filtering. Geophysics 2022, 87, A43–A48. [Google Scholar] [CrossRef]
  23. Li, C.; Wang, Y.; Zhang, J.H.; Geng, J.H.; You, Q.Y.; Hu, Y.X.; Liu, Y.Z.; Hao, T.Y.; Yao, Z.X. Separating Scholte Wave and Body Wave in OBN Data Using Wave-Equation Migration. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5914213. [Google Scholar] [CrossRef]
  24. Lei, J.H.; Liu, L.B.; Wan, W.X.; Zhang, S.R. Variations of electron density based on long-term incoherent scatter radar and ionosonde measurements over Millstone Hill. Radio Sci. 2005, 40, 1–10. [Google Scholar] [CrossRef]
  25. Yngvesson, K.O.; Perkins, F.W. Radar Thomson scatter studies of photoelectrons in the ionosphere and Landau damping. J. Geophys. Res. 1968, 73, 97–110. [Google Scholar] [CrossRef]
  26. Akbari, H.; Bhatt, A.; Hoz, C.L.; Semeter, J. Incoherent Scatter Plasma Lines: Observations and Applications. Space Sci. Rev. 2017, 212, 249–294. [Google Scholar] [CrossRef]
  27. Yue, X.; Liu, F.; Wang, J.; Wang, Y.; Cai, Y.; Ding, F.; Ning, B.; Li, M.; Zhang, N.; Wang, Z.; et al. On the Ion Line Calibration by Plasma Line in ISR Measurements. Remote Sens. 2023, 15, 1553. [Google Scholar] [CrossRef]
Figure 1. Full incoherent scatter radar spectrum as a function of altitude and Doppler shift, measured by 32−bit alternating code on 1 March 2021 13:06 LT (left) and a long pulse on 24 February 2021 13:05 LT (right). The center straight line is the ion line while the curve lines are the up− and downshifted plasma lines.
Figure 1. Full incoherent scatter radar spectrum as a function of altitude and Doppler shift, measured by 32−bit alternating code on 1 March 2021 13:06 LT (left) and a long pulse on 24 February 2021 13:05 LT (right). The center straight line is the ion line while the curve lines are the up− and downshifted plasma lines.
Remotesensing 15 02634 g001
Figure 2. Distorted long-pulse plasma line spectrum (left) and deconvolved plasma line spectrogram with background noise removed using the CLEAN algorithm (right).
Figure 2. Distorted long-pulse plasma line spectrum (left) and deconvolved plasma line spectrogram with background noise removed using the CLEAN algorithm (right).
Remotesensing 15 02634 g002
Figure 3. Detectable plasma line profile measured by long pulse.
Figure 3. Detectable plasma line profile measured by long pulse.
Remotesensing 15 02634 g003
Figure 4. The raw upshifted plasma line power spectrum on 8 March 2021.
Figure 4. The raw upshifted plasma line power spectrum on 8 March 2021.
Remotesensing 15 02634 g004
Figure 5. An example of signal–noise separation from the raw power spectrum using the IMF method in Figure 4. (a,c,e) The pure plasma line power spectrum obtained using the IMF method. (b,d,f) The noise separated using the IMF method. The white circles represent the extracted discrete plasma frequency points.
Figure 5. An example of signal–noise separation from the raw power spectrum using the IMF method in Figure 4. (a,c,e) The pure plasma line power spectrum obtained using the IMF method. (b,d,f) The noise separated using the IMF method. The white circles represent the extracted discrete plasma frequency points.
Remotesensing 15 02634 g005
Figure 6. The sections of the power spectrum before and after using the IMF method at various altitudes. (a,c,e) They are all from the raw power spectrum at altitudes of 260.18, 240.70, 171.74 km in Figure 4. (b,d,f) They are the results after using the IMF method at corresponding heights. The red circles indicate the true plasma frequency position. The green circles represent obvious strong background signal.
Figure 6. The sections of the power spectrum before and after using the IMF method at various altitudes. (a,c,e) They are all from the raw power spectrum at altitudes of 260.18, 240.70, 171.74 km in Figure 4. (b,d,f) They are the results after using the IMF method at corresponding heights. The red circles indicate the true plasma frequency position. The green circles represent obvious strong background signal.
Remotesensing 15 02634 g006
Figure 7. Comparison of fitting results with raw enhanced plasma lines. The white dotted line is the plasma line profile after α-Chapman function fitting.
Figure 7. Comparison of fitting results with raw enhanced plasma lines. The white dotted line is the plasma line profile after α-Chapman function fitting.
Remotesensing 15 02634 g007
Figure 8. Comparison of electron density derived from ion line, plasma line and ionosonde analyses. Subplot (a) is fitted from ion line. Subplot (b) is derived from plasma line. Subplot (c) is derived from ionosonde.
Figure 8. Comparison of electron density derived from ion line, plasma line and ionosonde analyses. Subplot (a) is fitted from ion line. Subplot (b) is derived from plasma line. Subplot (c) is derived from ionosonde.
Remotesensing 15 02634 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hao, H.; Zhao, B.; Yue, X.; Ding, F.; Ning, B.; Zeng, L. Study on the Method of Extracting Plasma Lines Based on Sanya Incoherent Scatter Radar. Remote Sens. 2023, 15, 2634. https://doi.org/10.3390/rs15102634

AMA Style

Hao H, Zhao B, Yue X, Ding F, Ning B, Zeng L. Study on the Method of Extracting Plasma Lines Based on Sanya Incoherent Scatter Radar. Remote Sensing. 2023; 15(10):2634. https://doi.org/10.3390/rs15102634

Chicago/Turabian Style

Hao, Honglian, Biqiang Zhao, Xinan Yue, Feng Ding, Baiqi Ning, and Lingqi Zeng. 2023. "Study on the Method of Extracting Plasma Lines Based on Sanya Incoherent Scatter Radar" Remote Sensing 15, no. 10: 2634. https://doi.org/10.3390/rs15102634

APA Style

Hao, H., Zhao, B., Yue, X., Ding, F., Ning, B., & Zeng, L. (2023). Study on the Method of Extracting Plasma Lines Based on Sanya Incoherent Scatter Radar. Remote Sensing, 15(10), 2634. https://doi.org/10.3390/rs15102634

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop