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]:
where
is the signal power spectrum,
is the transformed travel time,
is the migration velocity,
is the wavenumber, and
is the circular frequency of the echo signal. The exponential term in Equation (1) contributes to the phase change when
; on the contrary, it is an attenuation term and the signals are migrated to the evanescent region during de-migration, when
. 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]
where
and
are the F2 layer peak height and critical frequency, respectively.
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,
,
,
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]
where
is the plasma frequency,
is the electron plasma frequency, which is also called Langmuir frequency,
is the electron cyclotron frequency,
is the Boltzmann constant,
is the electron temperature,
is electron mass,
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
. Therefore, if the plasma frequency can be accurately recovered, then the absolute electron density profile can be calculated.