Next Article in Journal
Assessment of the Capability of Landsat and BiodivMapR to Track the Change of Alpha Diversity in Dryland Disturbed by Mining
Previous Article in Journal
Impacts of Shape Assumptions on Z–R Relationship and Satellite Remote Sensing Clouds Based on Model Simulations and GPM Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Ion Line Calibration by Plasma Line in ISR Measurements

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(6), 1553; https://doi.org/10.3390/rs15061553
Submission received: 12 February 2023 / Accepted: 10 March 2023 / Published: 12 March 2023

Abstract

:
The radar constant calibration in incoherent scatter radar ion line processing is essential for the data quality and was not paid enough attention in previous studies. In this investigation, based on several experiments made by the newly built Sanya incoherent scatter radar (SYISR), we made and evaluated the ion line calibration by plasma line both in case study and statistically. The calibration factor had local time and altitude variations, due to the corresponding variations of the transmitted power, the radar gain, and the noise temperature. We obtained a mean calibration factor of 1.35 by the simultaneous measured plasma line and ion line electron density and applied it to a one-month ion line observation calibration. Through a co-located ionosonde measured foF2 evaluation, the calibration decreased the mean deviation from −1.92 MHz (−18%) to −0.33 MHz (−3%), which resulted in much better agreement between the ion line foF2 after calibration and the ionosonde results. The existed deviations between after calibration and ionosonde results were due to the uncertainties either in the used calibration factor or the ionosonde measurements. An empirical Te/Ti usage in raw electron density estimation and ignoring the seasonal and short-term variations of the effecting factors might influence the calibration performance. Using the to-be-completed SYISR Tristatic System, the performance of plasma line calibration technique is expected to be improved in the future.

1. Introduction

The incoherent scatter radar (ISR) has been one of the most powerful ionospheric monitoring techniques since it was invented [1]. Based on radar equation, the received signal-to-noise ratio (SNR) of monostatic ISR scattered from electrons with respect to the electron number density (Ne) under the condition of ignoring the Debye effect is expressed as follows [2]:
S N R = ( P t G T s ) ( c τ 2   r 0 2 λ 2 K b B w ( 4 π ) R 2 ) ( N e 1 + T e / T i )
where the parameters in the first bracket, including the transmitted power (Pt), the radar gain (G), and the noise temperature (Ts), are the variables that are difficult to accurately determine and, therefore, need calibration; the parameters in the second bracket, including the bandwidth (Bw), the wave length ( λ ), effective pulse width ( τ ), light speed (c), electron radius (r0), Boltzmann constant (Kb), and range (R), have known values for specific radar and experimental setup; the parameters in the third bracket are related to the ionosphere and need derivation from the received SNR. Using Equation (1), the raw electron density can be derived based on the measured power by ISR and the predetermined ratio Te/Ti by an empirical model. In this study, the international reference ionosphere (IRI) model was used to predefine the ratio Te/Ti, regardless which type of waveform was used in the experiment [3]. However, in the actual operation of ISR, the transmitted power (Pt), the radar gain (G), and the noise temperature (Ts) usually have some deviations from the theoretical value, and, therefore, the derived raw electron density is biased to some extent, and needs further calibration for some usage requiring absolute electron density value. For the traditional centralized Klystron ISR, the factors that result in these deviations include the transmitter frequency switch, the antenna scanning, the long-term stability of the transmitter and receiver performance [4], and even external environment such as weather condition. For example, Rietveld et al. reported that during some extreme situations such as wet snow, the calibration factor might change to some extent [5]. For the phased array radar, there exist additional factors such as the stability of usable transmitting/receiving (T/R) units and the beam direction, that can result in such kind of deviation due to its emission mechanism [6,7].
In previous reports, several independent measurements, including ionosonde measured critical frequency and ISR plasma line-derived electron density, were used to calibrate the ISR ion line measurements [4,6,7]. In comparison with ionosonde calibration, the plasma line calibration has the following advantages: (1) the plasma line spectra can be used to derive the absolute electron density profiles with a high accuracy and resolution, which can give more accurate calibration [8,9]; (2) the plasma line could be observed by ISR simultaneously, which enable the calibration independent of an additional radar. However, the amplitude of plasma line spectra is usually several orders lower than that of the ion line; therefore, its detectability depends on the radar sensitivity, local time, and radar geometry very much [10].
In EISCAT UHF radar measurement, the ion line data was calibrated through separately determining the transmitted power, the received power, and the radar constant in the radar equation [4]. Specifically, the transmitted power was obtained from a monitor of the applied high voltage given the known transmitter performance. The received power was measured by comparison with a reference noise signal injected into the receiver system before the first preamplifier. The radar constant was derived from the simultaneous ion line and plasma line observations. Through independent evaluation, their absolute accuracy after calibration was ~10%. The newly built Resolute ISR (RISR) tried both the plasma line calibration and independent co-located ionosonde calibration in their multiple beams’ measurements and found reasonably good agreement between both type of observations [6,7]. In addition, Setov et al. and Alsatkin et al. reported that they used the radar measured sky noise to calibrate the ion line data [11,12].
Although there exist several publications mentioned above that involved ISR ion line calibration, this issue was not given enough attention. Especially, there is still a lack of a special articles dedicated to introduce the method and evaluate its effect systematically. In this paper, we attempted to use the simultaneous measured plasma line to calibrate the ion line electron density for the newly built Sanya incoherent scatter radar (SYISR) [13,14]. We used independent measurement to evaluate the performance of the calibration statistically. Such investigation is important and necessary for a newly built radar, which can help to understand and use the produced data product properly. Furthermore, on the basis of SYISR, we are currently developing the SYISR Tristatic System (SYISR-TS) [15]. Specifically, the SYISR-TS will double the current SYISR transmitter antenna array and build two new receivers. Once the SYISR-TS is completed, its performance, especially the plasma line detection capability, will be greatly enhanced given the doubled radar gain and aperture. We are considering to run the plasma line mode as a normal operational mode in the future by SYISR-TS, although this will require much more computational and storage resources. By doing that, we can potentially calibrate the ion line measurements by plasma line in real time. Therefore, this paper could be a necessary demonstration for that purpose. In Section 2, we will briefly describe the experimental data and calibration methodology. The main results will be given in Section 3. We finally present the discussion in Section 4 and conclude the whole investigation in Section 5, respectively.

2. Data and Methodology

The SYISR was a newly developed ISR located in Sanya (18.3°N, 109.6°E), led by the Institute of Geology and Geophysics, Chinese Academy of Sciences, with the main features of an active digital phased array, with all solid-state transmission and digital receiving. The readers can refer to Yue et al. for detailed information regarding the construction process, the hardware composition, radar index and evaluation, signal processing, and preliminary experimental results of SYISR [13,14]. Several detection modes, including zenith stare, perpendicular to the geomagnetic field, meridian scan, and all sky scan, were developed for ionospheric monitoring. The most common waveforms used by ISR, such as Barker code, long pulse, and alternating code, were implemented in SYISR. The detailed ion line processing and evaluation were given by Hao et al. [16]. The hard target, which occurred frequently over Sanya and challenged the data processing significantly, was eliminated by the method described in Wang et al. [17].
The plasma line spectra are induced by the electrostatic Langmuir waves in the plasma [8]. It is usually too weak to be observed by ISR. However, its intensity can be enhanced significantly due to the increased photoelectron flux around sunrise and sunset time or in the aurora region, which enables the detectability of plasma line by ISR. However, it is still much weaker in comparison to that of ion line [18]. Therefore, the detectability of plasma line challenges the sensitivity of ISR [10]. One outstanding feature of SYISR is that it can detect the weak plasma line spectra in addition to the ion line spectra [14]. The maximum bandwidth of SYISR was 4 MHz, which was not sufficient to cover the entire plasma line spectra. To resolve this problem, we used three parallel channels to splice a maximum bandwidth of 12 MHz in one side of ion line to cover the plasma line spectra dynamical range through the optimized configuration of the center frequency for each channel [14]. By our test, SYISR could observe plasma line over Sanya around sunrise and sunset time, and occasionally during the whole daytime. Ignoring the temperature related term, which is insignificant, in the dispersion relation of plasma line, the electron density (Ne) can be calculated from plasma line frequency ( ω P L ) using the following formula [8,19]:
N e = e 2 ε 0 m e ( ω P L 2 ( e B m e ) 2 sin 2 ( α ) )
where e, ε 0 , me, B, and α are charge constant, electrical constant, electron static mass, geomagnetic field intensity, and the angle between scattered beam direction and geomagnetic field, respectively. For SYISR, the α is ~65° when the beam points to the zenith direction.
In this study, five different experimental data were used based on SYISR, the configurations of which are summarized in Table 1. Using experimental #1, we obtained simultaneous ion line and plasma line, which were used to calculate the calibration factor. A 100 μs long pulse waveform selected here was a balanced choice of range resolution and spectra intensity of plasma line. The experimental #2 data was used to evaluate the effect of calibration on multiple parameters of alternating code measured ion line. In experimental #3, a one-month long experimental data were used to statistically analyze the effect of calibration. The experimental #4–5 were used in the Section 4. In all five experiments, the radar was run with a single beam pointing to the zenith direction over Sanya to ensure maximized gain and SNR.
As an example, Figure 1 shows a simultaneously observed ion line spectra and plasma line spectra by SYISR at 28 December 2021 10:14:17 local time. The ion line spectra shows enhanced intensity around E and F2 layers, which is very common in ISR measurements [14]. However, the usual double crests structure of ion line spectra was not identified. This was due to the relatively lower frequency resolution of 100 μs long pulse waveform used in the experiment. However, this did not influence the calibration effect because the calibration factor was calculated using the raw electron density derived from received power using Formula (1) and the plasma line-derived electron density using Formula (2). In the plasma line spectra, although its intensity was not very significant, it still showed obvious altitude variation of the ionosphere. Furthermore, according to Formula (1), the altitude variation of the signal power was determined by both range and ionosphere status (electron density and Te/Ti). We can see that the plasma line intensity shown in Figure 1b had obvious enhancement in the F2 peak altitude, corresponding to the highest ionospheric electron density. However, the effect of range variation was not very obvious, which was probably due to the relative narrow altitude interval in this experiment. We also investigated other experiments with higher altitude coverage and found that the plasma line spectra intensity evidently decreased with altitude increase. Before using Formula (2) to calculate the electron density profile from the spectra, we first denoised the spectra image and identified the plasma line profile. In the bottom panel of Figure 1, the corresponding ion line raw electron density and the plasma line-derived electron density are also displayed. We can see that the plasma line-derived electron density was generally larger than that of the ion line. This is understandable since Formula (1) assumed a 100% radiation efficiency, which is usually lower than 100%. This was expected to result in underestimation of electron density. In addition, the profile shape especially in the lower altitude showed difference between two results. This might be due to the uncertainty of the used empirical Te/Ti in Formula (1).
Furthermore, observations made by a co-located ionosonde over Sanya, named portable digital ionosonde (PDI), was also used for evaluation purpose in the study. The readers can refer to Lan et al. for detailed information regarding PDI [20]. In experiment #1, the corresponding original ionograms were manually scaled, while in experiment #3, the automatically scaled peak density and height were used.

3. Results

3.1. The Calibration Factor

The calibration factor was defined as the ratio between plasma line derived electron density to the ion line raw electron density in this study. In the actual operation of SYISR, the plasma line was not a regular monitoring mode. In the occasional experiment, the plasma line spectra were only obvious around sunrise and sunset time over Sanya. During the experiment on 28 December 2021, we fortunately measured clear plasma line spectra during 09:00–16:30 local time, the data of which, therefore, were used to determine the calibration factor in this study. Figure 2 shows the local time and altitude variation of electron density derived from the SYISR plasma line (a), ion line (b), and the co-located PDI (c). Three observations showed generally similar time variations. The plasma line electron density was systematically larger than that of the ion line. The PDI electron density profiles showed obvious fluctuations above the peak height. This was probably because the profiles by ionosondes were extrapolated under the assumption of constant scale height in the topside ionosphere rather than directly measured [21].
Figure 3 further shows the electron density ratio between the plasma line and ion line (a) and between the plasma line and PDI (b) during the same day as Figure 2. Regarding the ratio of plasma line to ion line, its value ranged from ~0.8 to ~2.5, with obvious local time and altitude variations, which might be explained by the corresponding variations of the uncertain terms in Formula (1). The altitudinal mean value showed no significant local time variation, with the values oscillating around ~1. The local time mean value had evident altitude variation, with values ranging from ~1 to ~2. For the ratio of plasma line to ionosonde, the main feature was the relatively larger value above the peak height region. Again, this was due to the extrapolated density profile in ionosonde inversion under the assumption of constant scale height in the topside ionosphere, which deviated from the real condition. During the altitude range of ~200 km to the peak height, the mean value of this ratio was around ~1, which indicated good consistency between the plasma line and ionosonde in this altitude interval.
Since few plasma line experimental data were available, it was difficult to derive accurate altitude and local time variations of the calibration factor statistically. Therefore, we calculated the mean value of plasma line to ion line ratio around peak height in Figure 3a, whose value was 1.35, to represent the calibration factor in the experiment hereafter. The corresponding uncertainty of this factor, represented by the standard deviation in the mean value calculation, was ~0.24.

3.2. Calibration Effect on Multiple Ionospheric Parameters

Theoretically, power calibration will not influence the temperature and ion velocity determination. However, in the least square fitting of the auto correlation function (ACF), the electron density, temperature, and ion velocity were derived simultaneously numerically. To determine the calibration effect on other parameters, Figure 4 shows the retrieved multiple ionospheric parameters using alternating code waveform during 6 May 2022 and the difference between before and after calibration. We observed that, in this experiment, the raw ion line value underestimated the electron density up to ~0.5 × 1012 m−3, depending on the amplitude of the background ionosphere. However, the difference of electron temperature, ion temperature, and line of sight ion bulk velocity, between before and after calibration value, was negligible, which meant that the calibration had no effect on those parameters. Some occasional relatively larger value around lower altitude of nighttime might have been due to the retrieval noise because of the low SNR.

3.3. Statistical Effect of the Calibration on Ionospheric F2 Peak Results

To observe the calibration effect statically, we further used the above calibration factor to calibrate the continuously SYISR measurements for one whole month. Then, the co-located ionosonde critical frequency (foF2, in MHz) and peak height (hmF2, in km) were used to make an independent evaluation.
Figure 5 shows the statistical comparison of ionosonde foF2 and SYISR ion line-derived foF2 before and after calibration during the whole month. Before calibration, the ion line foF2 underestimated the ionosonde foF2 systematically with a mean deviation of −1.92 MHz. While after calibration, this value decreased to −0.33 MHz, although underestimation still existed.
To demonstrate the local time variation of the calibration effect, Figure 6 shows the mean local time variation of foF2 before and after calibration from SYISR ion line, and that from the ionosonde. We clearly observed that the after calibration foF2 showed much better agreement with the ionosonde result than before calibration. Before calibration, the relative deviation ranged from ~−25% to ~−4%, with the mean value of −18%. After calibration, the mean deviation decreased to −3%, ranging from ~−11% to ~+13%. The overall agreement after calibration between SYISR and PDI appeared better during the interval of 04:00 to 16:00. This was probably because the calibration factor was calculated based on the observations during the similar time interval. Furthermore, in both Figure 5 and Figure 6, the shown results were derived in the interval of March to April of 2022, during which the mean F10.7 index was ~112 sfu (solar flux unit, 1 sfu = 10−22 W m−2 Hz−1), which represented a middle-low solar activity level. Therefore, the absolute deviation between the after calibration and ionosonde measurement should be larger in higher solar activity since the corresponding background ionospheric electron density was higher. However, in terms of relative deviation, the results should be similar, given a single calibration factor was used in our calibration processing.

4. Discussions

4.1. Potential Reason of the Residuals after Calibration

As indicated in Figure 5 and Figure 6, the calibration improved the ion line raw foF2 significantly. However, deviations between calibrated foF2 and ionosonde foF2 results, including a mean deviation of −0.33 MHz and the local time variation of the calibration performance, still existed. We think the probable reasons for this include:
(1) The current used calibration factor was calculated from one experiment, which might not represent the radar status statistically. In addition, we used one single calibration factor rather than considering its altitude and local time variation. This will also result in deviation to some extent. Furthermore, we used the calculated calibration factor from the measurements in 2021 to calibrate the results in 2022. Since the radar circumstance might change versus times, this will also result in deviations. For example, the radar noise temperature might change with seasons; the usable T/R units might be different between two experiments; the radar stability might also change in a longer time period [5]. In the future, more plasma line experiments are needed to consider the altitude and local time variation of the calibration factor with statistical significance. In addition, an operational mode with observing both ion line and plasma line simultaneously might be very useful for the ion line calibration, although realizing this is very challenging from the perspective of radar resources and computation resources.
(2) In the evaluation, we used a co-located ionosonde measured foF2 and hmF2, both of which were direct measurements and were believed to be accurate enough. Scotto and Sabbagh used to compare the MillStone Hill ISR measured hmF2 with that from a co-located ionosonde during ~20 years [22]. They found that the ionosonde hmF2 tended to underestimate the ISR hmF2 with ~10 km, regardless of which kind of inversion was used in the ionogram processing. This indicated that the deviations in Figure 5 and Figure 6 might also have been due to the bias of ionosonde measurements. To confirm this issue, we also made a statistical comparison between the SYISR ion line hmF2 and ionosonde hmF2, as shown in Figure 7. The ion line hmF2 were derived from the original power profile and were not influenced by the calibration, therefore, should be accurate. As indicated, the ionosonde hmF2 underestimated the SYISR hmF2 by ~25.4 km on average, although both had a high correlation coefficient of 0.97. The larger underestimation in our results than that of Scotto and Sabbagh might have been due to different seasons and locations [22]. Furthermore, the deviation of ionosonde hmF2 from that of SYISR was most significant around the middle noon time, with the mean deviation potentially reaching up to ~80 km.

4.2. The Effect of the Used Te/Ti in the Formula (1)

When deriving the raw electron density using Formula (1), we used the empirical IRI model to provide the Te/Ti ratio [3]. However, this was a climatological model, which might have deviated from the realistic situation; therefore, this resulted in the uncertainty of the calculated calibration factor. Especially, the derived calibration factor shown in Figure 3 had altitude dependency, which might be related to the IRI Te/Ti error. To clarify this point, we made another simultaneous plasma line and ion line experiment with combined 100 μs LP waveform and 32-bit AC waveform during 6 January 2022. Therefore, the AC measurements could be used to derive the Te/Ti self-consistently through the ACF fitting [16]. Unfortunately, the visible PL spectra, that can be used to derive electron density profile successfully, only lasted for ~2 h in this experiment. However, it was enough for our clarification purpose here, especially because this valid time interval covered the typical Te/Ti enhancement around the dawn sector. Figure 8 shows the comparison of Te/Ti and the corresponding ratio of PL electron density to LP electron density between IRI model calculation and SYISR AC measurements.
We observed that the overall local time and altitude variation of Te/Ti between IRI model and SYISR AC measurement appeared similar. For example, both Te/Ti peaked around ~200–300 km, and had higher value around 07:00 due to the photoelectron heating of the low-density thermal electrons [23]. However, the AC measured Te/Ti value and its peak height was more or less lower than that of IRI model after 07:30 local time. This is understandable since the IRI model was a climatological model, while the realistic ionosphere showed significant day-to-day variability, whose effect was most significant around morning hours [24]. In addition, the AC measured Te/Ti was much nosier than the model result because the electron density was much lower around the dawn sector, which resulted in lower measured SNR and, therefore, larger retrieval uncertainty of Te/Ti [14].
After replacing the IRI Te/Ti with the AC-measured Te/Ti in deriving the LP electron density using the radar equation, the electron density ratio of PL to LP measurements showed some differences. The variability among different local times became smaller, especially in the higher altitudes, which can be indicated from the degree of dispersion among different color lines. The mean ratio with AC-measured Te/Ti generally increased with the altitude increase, while for the IRI Te/Ti, it showed a decreasing trend during the 250–300 km altitude interval. However, the altitude dependency of the ratio still existed with AC measured Te/Ti, which implied that the altitude variation of the calibration factor not only resulted from the Te/Ti factor. Furthermore, in terms of mean profile of the ratio, the difference between both was insignificant, with IRI Te/Ti based ranges from 1.13 to 2.26 and AC Te/Ti ranges from 1.35 to 2.72. The corresponding mean value was ~1.5 and ~1.6, respectively. Please note that the mean altitude dependency here was slightly different from that shown in Figure 3, which was probably due to local time coverage difference in both experiments. In summary, replacing the empirical Te/Ti ratio with the AC measured Te/Ti ratio, the calibration factor will be slightly different in terms of both mean value and altitude variation. However, the altitude dependency still exists.

4.3. Shortcomings of the Current Calibration Method

Most other ISRs, especially those equipped with parabolic antenna, such as the EISCAT radar [25], usually employ a noise injection technique to monitor the noise temperature. In addition, they can measure the transmitted power by a separate channel. Therefore, they can separate the effect of transmitter power and noise temperature in their calibration process. However, in our phased array SYISR radar, we did not design a mechanism to monitor the noise temperature and the transmitted power, given that we had 4096 channels in total and a limited budget [14]. Therefore, our calibration accuracy was determined by a mixed effect of the variations of the transmitted power, the noise temperature, and radar constant.
In the current calibration investigation, we used a constant calibration factor to make the correction. However, most effecting factors, including the transmitted power, the gain and the noise temperature, might have time variations. Generally, some components, especially the T/R unit, might break down during the operation for a phased radar. By our monitoring, this was a very small probability event primarily because SYISR is a brand new developed radar. So, the breakdown of the T/R unit resulted variations of the transmitted power and gain could be ignored in the short term. The noise temperature roughly included the system noise temperature and the external noise temperature [14]. In the factory test stage, we measured the system noise coefficient and gain for each T/R unit, amplifier, digital receiver, and the SNR loss of the feeder network and cable. Then, we calculated the noise temperature of the system using those measured value [14]. In the field test stage, we re-tested the noise coefficient of the sampled 1% T/R units and found that it was relatively stable. The external noise temperature was determined by many factors, such as atmospheric absorption, cosmic noise, and even anthropogenic noise. Most of these factors can vary from time to time and, therefore, result in the short-term variability of external noise temperature [26]. Using one constant calibration factor definitely cannot indicate this short-term variability. To demonstrate the short-term time variation of the noise, we made a sampling experiment without transmitting power using SYISR on 26 November 2022. As shown in Figure 9, the measured noise showed short-time variation, with peaks around 07:00 and 17:30. The variability (maximum minus minimum) was ~6% relatively to the mean value. In real signal processing, we usually prepared the range gate as high as possible (~1200 km in this article). Then, the echo from highest altitude was considered as the background noise. Doing this would produce the defect of non-noise injection mechanism in the radar to some extent.

4.4. Applicability to Multiple Beams Experiment

In this article, we only calculated the calibration factor and evaluated its performance in single vertical beam experiments. Actually, the SYISR can be configured to make measurements in arbitrary beam direction within its effective detection area [14]. In realistic operation, we also ran the radar in multiple beams mode for the purpose of monitoring the regional variation of the ionosphere. For those beams that deviated from the vertical direction, the current calibration factor was not valid anymore, since the corresponding beam width, radar gain, and the environmental noise were different. However, the current method could also be used to calculate the corresponding calibration factor if simultaneous plasma line measurements are available. We actually made some multiple beams experiments with simultaneous ion line and plasma line measurements. However, due to the decreased gain and, therefore, SNR in non-vertical direction, the effect, especially the visibility of the plasma line spectra in comparison with that of the vertical direction under the same waveforms, was not good enough. As stated previously, once the SYISR-TS is completed, we will conduct another experiment to carry out the multiple beam’s calibration investigation using the plasma line. Considering that both the transmitted power and gain will be increased, we expect this plasma line calibration to work well.

5. Conclusions

In this study, based on the plasma line experiment by the newly built SYISR, we calculated the calibration factor by the simultaneous plasma line and ion line electron density. We then applied this calibration factor to several SYISR ion line experiments and made a statistical evaluation on the calibration effect. The main conclusions are summarized as follows:
(1) The calibration factor had local time and altitude variations, which might be due to the corresponding variations of the transmitted power, the radar gain, and the noise temperature in the soft target radar equation. Then, a mean calibration factor of 1.35 was obtained by averaging the corresponding value around peak height during one specific plasma line experiment.
(2) We then applied this calibration factor to calibrate one-month SYISR measurements and made a statistical comparison between SYISR result and a co-located ionosonde measured foF2 and hmF2. The results showed that the calibration decreased the mean deviation from −1.92 MHz (−18%) to −0.33 MHz (−3%), which resulted in much better agreement between the ion line foF2 after calibration and the ionosonde results. The existed deviations between after calibration and ionosonde results were due to the uncertainties either in the used calibration factor or the ionosonde measurements.
(3) Considering that the SYISR had no noise injection mechanism to monitor the radar’s noise temperature and no additional channel to measure the transmitted power in operation, the current calibration cannot distinguish the effects of those factors. We used an empirical Te/Ti in the raw electron density estimation using the radar equation, which might deviate from the real situation to some extent. In addition, we did not consider the seasonal and short-term variations of those effecting factors such as noise temperature and transmitted power in this calibration. All the above facts should influence the calibration performance. Furthermore, the current calculated calibration factor is inapplicable to other non-vertical beam experiments, given that the corresponding beam width, radar gain, and the environmental noise were different. Considering that the developing SYISR-TS will have higher gain and power, we can use AC code to measure more obvious plasma line and self-consistently determine Te/Ti. More accurate calibration using plasma line is expected for both single vertical beam and multiple beams experiments.

Author Contributions

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

Funding

This research was funded by the Project of Stable Support for Youth Team in Basic Research Field of CAS (YSBR-018), and the National Natural Science Foundation of China (41427901). This research made use of the data from the Meridian Project.

Data Availability Statement

The experimental results in the figures could be obtained upon request (X.Y.).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gordon, W.E. Incoherent scattering of radio waves by free electrons with applications to space exploration by radar. Proc. IEEE 1958, 46, 1824–1829. [Google Scholar] [CrossRef]
  2. Lehtinen, M.S. Statistical Theory of Incoherent Scatter Radar Measurements. Ph.D. Thesis, University of Helsinki, Helsinki, Finland, 1986. [Google Scholar]
  3. Bilitza, D.; Pezzopane, M.; Truhlik, V.; Altadill, D.; Reinisch, B.W.; Pignalberi, A. The International Reference Ionosphere model: A review and description of an ionospheric benchmark. Rev. Geophys. 2022, 60, e2022RG000792. [Google Scholar] [CrossRef]
  4. Kirkwood, S.; Collis, P.N.; Schmidt, W. Calibration of electron densities for the EISCAT UHF radar. J. Atmos. Terr. Phys. 1986, 48, 773–775. [Google Scholar] [CrossRef]
  5. Rietveld, M.T.; Isham, B.; Häggström, I. Calibration of EISCAT Incoherent Scatter Radar Electron Densities and the Anomaly of 23–25 October 2003. In EISCAT Scientific Association Report; EISCAT: Kiruna, Sweden, 2005. [Google Scholar]
  6. Bahcivan, H.; Tsunoda, R.; Nicolls, M.; Heinselman, C. Initial ionospheric observations made by the new Resolute incoherent scatter radar and comparison to solar wind IMF. Geophys. Res. Lett. 2010, 37, L15103. [Google Scholar] [CrossRef]
  7. Gillies, R.G.; van Eyken, A.; Spanswick, E.; Nicolls, M.; Kelly, J.; Greffen, M.; Knudsen, D.; Connors, M.; Schutzer, M.; Valentic, T.; et al. First observations from the RISR-C incoherent scatter radar. Radio Sci. 2016, 51, 1645–1659. [Google Scholar] [CrossRef] [Green Version]
  8. Akbari, H.; Bhatt, A.; Hoz, C.L.; Semeter, J.L. Incoherent Scatter Plasma Lines: Observations and Applications. Space Sci. Rev. 2017, 212, 249–294. [Google Scholar] [CrossRef]
  9. Showen, R.L. The spectral measurement of plasma lines. Radio Sci. 1979, 14, 503–508. [Google Scholar] [CrossRef]
  10. Vierinen, J.; Bhatt, A.; Hirsch, M.A.; Strømme, 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]
  11. Setov, A.G.; Medvedev, A.V.; Lebedev, V.P.; Kushnarev, D.S.; Alsatkin, S.S.; Tashlykov, V.P. Calibration methods for absolute measurements at the Irkutsk incoherent scatter radar. In Proceedings of the 23rd International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics, Irkutsk, Russia, 3–7 July 2017; Volume 10466, p. 104665Y. [Google Scholar] [CrossRef]
  12. Alsatkin, S.S.; Medvedev, A.V.; Ratovsky, K.G. Features of Ne recovery at the Irkutsk Incoherent Scatter Radar. Sol.-Terr. Phys. 2020, 6, 77–88. [Google Scholar] [CrossRef]
  13. Yue, X.; Wan, W.; Ning, B.; Jin, L. An active phased array radar in China. Nat. Astron. 2022, 6, 619. [Google Scholar] [CrossRef]
  14. Yue, X.; Wan, W.; Ning, B.; Jin, L.; Ke, C.; Ding, F.; Zhao, B.; Zeng, L.; Deng, X.; Wang, J.; et al. Development of the Sanya Incoherent Scatter Radar and Preliminary Results. J. Geophys. Res. Space Phys. 2022, 127, e2022JA030451. [Google Scholar] [CrossRef]
  15. Li, M.; Yue, X.; Zhao, B.; Zhang, N.; Wang, J.; Zeng, L.; Hao, H.; Ding, F.; Ning, B.; Wan, W. Simulation of the Signal-to-Noise Ratio of Sanya Incoherent Scatter Radar Tristatic System. IEEE Trans. Geosci. Remote Sens. 2021, 59, 2982–2993. [Google Scholar] [CrossRef]
  16. 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]
  17. Wang, J.; Yue, X.; Ding, F.; Ning, B.; Jin, L.; Ke, C.; Zhang, N.; Luo, J.; Wang, Y.; Yin, H.; et al. The Effect of Space Objects on Ionospheric Observations: Perspective of SYISR. Remote Sens. 2022, 14, 5092. [Google Scholar] [CrossRef]
  18. Feng, T.; Liu, M.; Xu, B.; Xu, T.; Gao, S.; Wang, X.; Zhou, C. Auroral-enhanced plasma lines by suprathermal electrons observed by EISCAT. J. Geophys. Res. Space Phys. 2021, 126, e2020JA028495. [Google Scholar] [CrossRef]
  19. Salpeter, E.E. Electron density fluctuations in a plasma. Phys. Rev. 1961, 120, 1528. [Google Scholar] [CrossRef]
  20. Lan, J.; Ning, B.; Li, G.; Zhu, Z.; Hu, L.; Sun, W. Observation of short-period ionospheric disturbances using a portable digital ionosonde at Sanya. Radio Sci. 2018, 53, 1521–1532. [Google Scholar] [CrossRef]
  21. Reinisch, B.W.; Huang, X. Deducing topside profiles and total electron content from bottomside ionograms. Adv. Space Res. 2001, 27, 23–30. [Google Scholar] [CrossRef]
  22. Scotto, C.; Sabbagh, D. The Accuracy of Real-Time hmF2 Estimation from Ionosondes. Remote Sens. 2020, 12, 2671. [Google Scholar] [CrossRef]
  23. Dalgarno, A.; McElrow, M.B. Ionospheric electron temperatures near dawn. Planet. Space Sci. 1965, 13, 143–145. [Google Scholar] [CrossRef]
  24. Zhou, X.; Liu, H.L.; Lu, X.; Zhang, R.; Maute, A.; Wu, H.; Yue, X.; Wan, W. Quiet-time day-to-day variability of equatorial vertical E × B drift from atmosphere perturbations at dawn. J. Geophys. Res. Space Phys. 2020, 125, e2020JA027824. [Google Scholar] [CrossRef]
  25. Lehtinen, M.; Markkanen, J.; Väänänen, A.; Huuskonen, A.; Damtie, B.; Nygrén, T.; Rahkola, J. A new incoherent scatter technique in the EISCAT Svalbard Radar. Radio Sci. 2002, 37, 3-1–3-14. [Google Scholar] [CrossRef]
  26. Skolnik, M.I. Introduction to Radar System, 3rd ed.; McGraw-Hill Education: New Delhi, India, 2001; pp. 482–535. [Google Scholar]
Figure 1. An example of simultaneously measured ion line (a) spectra and plasma line (b) spectra (top row) and the corresponding derived electron density profiles (bottom row). The experiment was produced by SYISR at 28 December 2021 10:14:17 local time.
Figure 1. An example of simultaneously measured ion line (a) spectra and plasma line (b) spectra (top row) and the corresponding derived electron density profiles (bottom row). The experiment was produced by SYISR at 28 December 2021 10:14:17 local time.
Remotesensing 15 01553 g001
Figure 2. Local time variation of electron density profiles derived from simultaneous plasma line spectra (a) and ion line power (uncalibrated, (b)) obtained by SYISR and a collocated ionosonde (c) during 28 December 2021.
Figure 2. Local time variation of electron density profiles derived from simultaneous plasma line spectra (a) and ion line power (uncalibrated, (b)) obtained by SYISR and a collocated ionosonde (c) during 28 December 2021.
Remotesensing 15 01553 g002
Figure 3. The ratio of plasma line-derived electron density to ion line derived raw electron density (a) and plasma line-derived electron density to ionosonde electron density (b) during 28 December 2021. The white line in each colormap represents the corresponding ionospheric peak height. The mean ratio variations versus local time and altitude are, respectively, shown in the bottom and right subplots of each colormap.
Figure 3. The ratio of plasma line-derived electron density to ion line derived raw electron density (a) and plasma line-derived electron density to ionosonde electron density (b) during 28 December 2021. The white line in each colormap represents the corresponding ionospheric peak height. The mean ratio variations versus local time and altitude are, respectively, shown in the bottom and right subplots of each colormap.
Remotesensing 15 01553 g003
Figure 4. Local time and altitude variations of the retrieved electron density (a), electron temperature (b), ion temperature (c), and line-of-sight ion bulk velocity (d) during the entire day of 6 May 2022. The left panels are the results after calibration, while the right panels are the differential results between parameters before-calibration and after-calibration.
Figure 4. Local time and altitude variations of the retrieved electron density (a), electron temperature (b), ion temperature (c), and line-of-sight ion bulk velocity (d) during the entire day of 6 May 2022. The left panels are the results after calibration, while the right panels are the differential results between parameters before-calibration and after-calibration.
Remotesensing 15 01553 g004
Figure 5. Versus plot between ionosonde measured ionospheric critical frequency (foF2) and SYISR ion line-derived foF2 (top panels) and the corresponding statistical histogram of their difference (bottom panels) for both before calibration (left) and after calibration (right) during the one-month long experiment of 12 March 2022 to 10 April 2022. The corresponding correlation coefficient, mean deviation, and root mean square error between SYISR and Ionosonde foF2 were also calculated in both conditions. The color in the top panel represents the data number intensity in each pixel.
Figure 5. Versus plot between ionosonde measured ionospheric critical frequency (foF2) and SYISR ion line-derived foF2 (top panels) and the corresponding statistical histogram of their difference (bottom panels) for both before calibration (left) and after calibration (right) during the one-month long experiment of 12 March 2022 to 10 April 2022. The corresponding correlation coefficient, mean deviation, and root mean square error between SYISR and Ionosonde foF2 were also calculated in both conditions. The color in the top panel represents the data number intensity in each pixel.
Remotesensing 15 01553 g005
Figure 6. Local time variation of mean foF2 during 12 March 2022 to 10 April 2022 of ionosonde measurements (red), SYISR ion line before calibration (green) and after calibration (blue).
Figure 6. Local time variation of mean foF2 during 12 March 2022 to 10 April 2022 of ionosonde measurements (red), SYISR ion line before calibration (green) and after calibration (blue).
Remotesensing 15 01553 g006
Figure 7. Comparison of ionospheric peak height (hmF2) between ionosonde measurements and SYISR ion line results during 12 March 2022 to 10 April 2022 in terms of versus plot (left panel) and mean local time variation (right panel).
Figure 7. Comparison of ionospheric peak height (hmF2) between ionosonde measurements and SYISR ion line results during 12 March 2022 to 10 April 2022 in terms of versus plot (left panel) and mean local time variation (right panel).
Remotesensing 15 01553 g007
Figure 8. Local time and altitude variation of the Te/Ti ratio made by the IRI model calculation (a) and SYISR alternating code (AC) experiment (b) during 6 January 2022. The white lines in subplots (a,b) are the ionospheric peak heights derived from PL measurements. Subplots (c,d) are the electron density ratio of plasma line (PL) direct measurements to the 100 μs long pulse (LP) derivations before calibration using radar equation on the basis of Te/Ti corresponding to subplots (a,b), respectively. The thick black lines are the corresponding mean profiles, while the thin color lines represent the profiles of different local times.
Figure 8. Local time and altitude variation of the Te/Ti ratio made by the IRI model calculation (a) and SYISR alternating code (AC) experiment (b) during 6 January 2022. The white lines in subplots (a,b) are the ionospheric peak heights derived from PL measurements. Subplots (c,d) are the electron density ratio of plasma line (PL) direct measurements to the 100 μs long pulse (LP) derivations before calibration using radar equation on the basis of Te/Ti corresponding to subplots (a,b), respectively. The thick black lines are the corresponding mean profiles, while the thin color lines represent the profiles of different local times.
Remotesensing 15 01553 g008
Figure 9. Local time and altitude variation of the measured noise in dBw in a non-transmitting experiment made by SYISR on 26 November 2022. The signal was integrated within 1 min.
Figure 9. Local time and altitude variation of the measured noise in dBw in a non-transmitting experiment made by SYISR on 26 November 2022. The signal was integrated within 1 min.
Remotesensing 15 01553 g009
Table 1. Experimental setup of SYISR observations used in the figures. LP: long pulse; AC: alternating code; BC: Barker code.
Table 1. Experimental setup of SYISR observations used in the figures. LP: long pulse; AC: alternating code; BC: Barker code.
#WaveformSampling RateTime IntervalFigures
1LP, 100 μs4 MHz (3 channels)28 December 20211–3
216-bit AC, 480 μs0.1 MHz (1 channel)6 May 20224
313-bit BC, 390 μs0.1 MHz (1 channel)12 March 2022 to 10 April 20225–7
432-bit AC, 100 μs LP4 MHz (3 channels)6 January 20228
5N/A0.1 MHz26 November 20229
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

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. https://doi.org/10.3390/rs15061553

AMA Style

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 Sensing. 2023; 15(6):1553. https://doi.org/10.3390/rs15061553

Chicago/Turabian Style

Yue, Xinan, Fanyu Liu, Junyi Wang, Yonghui Wang, Yihui Cai, Feng Ding, Baiqi Ning, Mingyuan Li, Ning Zhang, Zhongqiu Wang, and et al. 2023. "On the Ion Line Calibration by Plasma Line in ISR Measurements" Remote Sensing 15, no. 6: 1553. https://doi.org/10.3390/rs15061553

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