Next Article in Journal
Structural Characteristics of Endorheic Rivers in the Tarim Basin
Next Article in Special Issue
Analysis of the Ionospheric Irregularities and Phase Scintillation at Low and Middle Latitudes Based on Swarm Observations
Previous Article in Journal
Validation of Remote-Sensing Algorithms for Diffuse Attenuation of Downward Irradiance Using BGC-Argo Floats
Previous Article in Special Issue
A Novel Method for Improving Quality of Oblique Incidence Sounding Ionograms Based on Eigenspace-Based Beamforming Technology and Capon High-Resolution Range Profile
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Comparison of Electron Density between CSES In Situ and Ground-Based Observations in China

1
Institute of Earthquake Forecasting, China Earthquake Administration, Beijing 100036, China
2
National Key Laboratory of Electromagnetic Environment, China Research Institute of Radio Wave Propagation, Qingdao 266107, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(18), 4498; https://doi.org/10.3390/rs14184498
Submission received: 27 July 2022 / Revised: 29 August 2022 / Accepted: 6 September 2022 / Published: 9 September 2022

Abstract

:
As the observation accuracy of parameters in the ionosphere cannot be directly checked, the comparison with other observations is the main way to evaluate the data quality of satellite measurements. Through the comparative analysis between the in situ electron density (Ne) observed by the China Seismo-Electromagnetic Satellite (CSES) and Ne at about 500 km altitude detected by Qujing Incoherent Scatter Radar (ISR), it was found that the pattern of CSES Ne is consistent with that of ISR Ne, and the correlation coefficient between the two sets of data is above 0.88 for different groups according to the magnitude. The value of CSES Ne is lower than that of ISR Ne, and the median value of the ratio for the difference between the conjugate data is 84.04%. Based on the comparison in the daytime between CSES Ne and ionosonde observations in China, it was found that the trend of the two datasets is mostly similar, and the correlation coefficient in some locations can reach up to 0.7. The distribution of CSES Ne and correlation coefficients at different latitudes show that the relationship is relatively better around the peak of the equatorial ionization anomaly (EIA). The differences in the value between CSES Ne and ionosonde data also exist, the relative change of which is about 80–95% in the daytime.

1. Introduction

The China Seismo-Electromagnetic Satellite (CSES), also called ZhangHeng-1 (ZH-1), was successfully launched on 2 February 2018. There are eight scientific payloads on the CSES, including a high-precision magnetometer (HPM), an electric field detector (EFD), a search coil magnetometer (SCM), a plasma analyzer package (PAP), a Langmuir probe (LAP), a high energetic particle package (HEPP) and high energetic particle detector (HEPD), a GNSS occultation receiver (GOR) and a tri-band beacon (TBB), which can provide ionospheric observations of electromagnetic fields, plasma and energetic particles [1].
As the accuracy of physical parameters in the ionosphere cannot be directly checked, the comparison with other observations or models is the main way to evaluate the quality of the satellite data. For the CSES data, Zhou et al. [2] found that the three HPM sensors have good consistency after in-orbit linear correction. The comparison of HPM data between CSES and Swarm shows that both the vector and the scalar differences are quite consistent, and the responses to a geomagnetic storm were similar for the two satellites. Diego et al. [3] reported that the observation of EFD highlights a very good agreement with the theoretical value, and can respond to both natural and artificial sources. The results of Zhao et al. [4] show that the electromagnetic field recorded by CSES EFD and SCM has a good correlation with the simulated result. The correlation coefficients of electron density (Ne) observed by the LAP payload between CSES and the Detection of Electro-Magnetic Emissions Transmitted from Earthquake Regions (DEMETER) satellite, Swarm constellation and the International Reference Ionosphere (IRI) model are all high, and the variation patterns also have similar trends between conjugate data [5,6,7]. By analyzing the data in the commissioning phase, Picozza et al. [8] attested that the particle rates of the HEPD proton and electron are consistent with radiation models, as well as their distribution along the orbit. Ambrosi et al. [9] found that not only the satellite housekeeping data, but also the scientific data of HEPD show very good stability. Wang et al. [10] found that the measurements of CSES GOR are consistent with those of ionosonde and incoherent scatter radar (ISR) in Millstone Hill, including F2 layer peak electron density (NmF2), its height (hmF2) and electron density profile.
The study of data validation for other satellite missions have also been carried out by scientists, which gives confidence for scientific research using these data. The electron density and temperature of the Swarm constellation were compared with ISR, ionosonde and COSMIC data, and the sufficiently high correlation coefficients were reported by Lomidze et al. [11]. The study of Cai et al. [12] shows that the data from Global-Scale Observations of the Limb and Disk (GOLD) match well with the GPS TEC in the morphology of equatorial ionization anomaly (EIA) and its seasonal variability. By comparing the thermospheric neutral wind observed by the Ionospheric Connection Explorer (ICON) spacecraft, and four ground-based specular meteor radars, Harding et al. [13] found the agreement between space-based and ground-based winds. Pedatella et al. [14] reported that the morphology of the total electron content (TEC) from FORMOSAT-7/COSMIC-2 and GPS TEC is similar, and the absolute TEC of MOSAT-7/COSMIC-2 shows an overall consistency with Swarm-B observations.
The data quality is worth being researched for a long time. Eight years after Swarm launched in 2013, scientists are still comparing the Swarm data with other observations for future improvement [15,16]. Although the comparisons between LAP Ne and other observations/models have been carried out in previous studies [5,6,7], the data analysis is still needed from different aspects to meticulously and extensively evaluate the observations. CSES payloads operate in two different modes, called survey and burst modes. In most regions of the world, the satellite works in survey mode. When the satellite flies over China, the Circum-Pacific and Eurasia seismic belts, the burst mode is automatically triggered with a higher resolution. Fortunately, there is an incoherent scatter radar (ISR) station and several ionosonde stations in China, which support data to carry out the comparison between the in situ Ne of CSES with burst mode and the ground-based remote sensing of the ionosphere. In the article, Section 2 gives the description of the data. Section 3 and Section 4 illustrate the comparison between LAP Ne and ISR data according to the ionosonde observations. The discussion and conclusions are in Section 5 and Section 6, respectively.

2. Data

The CSES is a circular sun-synchronous satellite, with the descending and ascending nodes at 14:00 LT (local time) and 02:00 LT, respectively. The satellite is three-axes stable to set the altitude at about 507 km with an orbital inclination of 97.4°. As the orbit of the CSES is strictly revisited every five days, the global electromagnetic environment of the ionosphere can be obtained every five days [17].
In this study, we pay attention to the in situ electron density data observed by LAP. The LAP payload is installed on the front of the satellite with two spherical sensors and one electronics box. Sensor 1 is a larger one with a diameter of 5 cm, and the diameter of sensor 2 is 1 cm, designed as a backup. There is a long stub 50 cm in length for each sensor to reduce the interference from the satellite surface [18]. The data we used in this study were observed by sensor 1. The plasma parameters can be obtained by fitting the I–V characteristic curve of LAP based on the sweep voltage and sensor current, which is introduced in detail by Yan et al. [19]. The Ne measurement range of CSES LAP is 5 × 102–1 × 107 cm−3 with a relative accuracy of 10% [20]. For the LAP payload, the time resolutions of the burst and survey modes are 1.5 and 3 s, which means the latitudinal interval of each sample is about 0.1° and 0.2°, respectively. The electron density data observed by the CSES can be downloaded from the mission website (http://www.leos.ac.cn/, accessed on 10 August 2021) after registration. Currently, the CSES mission does not stamp the quality flag to LAP Ne. In our analysis, Ne samples with a value higher than 1 × 107 cm−3 or lower than 5 × 102 cm−3 were omitted, as it is considered that there may be bias in the fitting process if the data exceed the measurement range.
The accuracy of the ionosonde and ISR observations is widely accepted, and these data are usually applied to validate the satellite observations [11,21,22]. In China, there is an ISR in operation at the Qujing site and some ionosondes at different latitudes, the data of which can be applied to compare with the electron density of the CSES in the China region. The Qujing ISR (QJISR) was constructed by the China Research Institute of Radio Wave Propagation (CRIRP) at the beginning of 2012 [23], and the station is situated in Yunnan province with the latitude 25.6°N and longitude 103.7°E, which is plotted in Figure 1 with a green star. There are two groups of ionosonde data to compare with the CSES Ne. One is from the Institute of Earthquake Forecasting (IEF), China Earthquake Administration, including three ionosonde stations in the southwest China, which is marked by a blue point in Figure 1. Leshan station is located in Sichuan province, with the longitude 103.8°E and latitude 30.0°N, and Yunnan province has two stations, located in Tengchong (98.5°E, 25.0°N) and Puer (101.1°E, 22.7°N) [24]. The other group is from CRIRP. Since the 1940s, CRIRP has established several ionosondes from lower (20.0°N) to middle latitudes (49.6°N) in the mainland of China [25], and the locations of the data we used are marked by red triangles in Figure 1.

3. Comparison with ISR Data

QJISR works several hours in each day, almost all during the daytime, with a time resolution of about 5 min. Corresponding to the available data of the CSES Ne, 25 events were found for QJISR from 2018 to 2021. The ISR can give an Ne profile with an altitude resolution of about 5 km. The principles of the selected comparison dataset are as follows. (1) The mean value (M) and standard deviation (std) of QJISR Ne at the altitude 450–550 km from 13:00 to 15:00 (LT) were obtained to represent the observation at the satellite altitude in the daytime for each day. There are differences for the working time of QJISR each day; thus, the time range covers 2 h to ensure sufficient samples are obtained for the comparison between the conjugate data. (2) Corresponding to the date of QJISR observations in 25 days, the M and std of the CSES Ne in the region between 93.8–113.8°E and 15.6–35.6°N were calculated for the daytime data. The purpose of the selected ±10° region is to provide a sufficient number of observations. With the burst mode of CSES, the resolution of latitude is almost 0.1°; thus, there are about 200 data points to obtain the M and std.
Figure 2 gives the Ne time series for the CSES (red line) and QJISR (blue line) in the comparison days, which is not continuous and the detailed date is marked on the X axis. Every day, the circle and bar represent the mean value and standard deviation calculated according to the above principles. The two time series are displayed with different amplitude scales to better underline their own variations. The Ne variation pattern between the CSES and QJISR is almost similar, with a high value in 2018 and a relative low value in 2020 and 2021. Furthermore, some trends in a short time frame fit perfectly, for example, the data from 25 May to 7 June 2018.
To check the difference in the value of electron density, the relative change of CSES Ne corresponding to QJISR data was calculated according to Equation (1).
R e l = N e C S E S N e Q J I S R N e Q J I S R × 100  
where Ne_CSES represents the data of CSES Ne, and Ne_QJISR is the Ne at satellite altitude calculated based on the above principle. A histogram of the relative change value is shown in Figure 3. It can be seen that all of the ratios are negative with a median value −84.04%, which means the average value of CSES Ne is lower than that of QJISR Ne at the satellite altitude.
To further analyze the correlation between the two datasets, panel a of Figure 4 gives the Ne scatters of CSES and QJISR with the detailed date for each point. It seems that these data can be separated into two groups, which are marked by the red oval. It can be seen that except four points in 2018, the time range of the group 1 data is in 2020 and 2021. The data of group 2 belong to 2018, except one point in 2020. Combined with the curves in Figure 2, the magnitude of the group 2 data is mostly higher than that of the group 1 data. Since the available dataset for the comparison of CSES and QJISR Ne is limited, the characteristics of ionospheric variation cannot be completely considered, and the magnitude of observation may play a leading role in the correlation analysis. Based on the two groups, the data were separately fitted in panel b of Figure 4, with the blue and black points denoting the data of group 1 and 2, respectively. The relationship of the Ne data for CSES and QJISR is better represented by a linear correlation in each group with the empirical formula in Equation (2).
Y = 19.94 × X − 2.45 × 1011, Y = 11.49 × X − 4.3 × 1011
Furthermore, the correlation coefficient between the two datasets was calculated according to Equation (3).
R ( X , Y ) = C o v ( X , Y ) V a r [ X ] V a r [ Y ]  
where X and Y represent the Ne values of the CSES and QJISR. Cov means covariance between the two datasets, and Var is the variance of each dataset. The correlation coefficients are 0.9250 and 0.8884 for the conjugate data of groups 1 and 2, respectively, which means the Ne of the CSES and QJISR exhibits a strong correlation.

4. Comparing with Ionosonde Data

Based on the two groups of ionosonde data, the continuous-time observations of three stations from IEF were used to carry out time series analysis, and the data of 10 ionosonde stations from CRIRP were applied to analyze the correlation at different latitudes.
In the ionosphere, the electron density above the F2 region is determined by chemical production and loss, transport due to neutral wind, E × B drift, ambipolar diffusion and so on. In the daytime, the production, loss and transportation of electrons are mainly affected by the sun. The ionospheric electron density is mainly produced by the photoionization of solar EUV and X-ray radiation, which is represented by the fact that NmF2, TEC and Ne at the topside of the ionosphere linearly change with the solar proxies at low and medium solar activity levels [26,27,28]. Solar activity and photochemical equilibrium play a leading role in the variability of Ne in the daytime, and the electron density at different altitudes may synchronously variate with the solar effect. In the nighttime, the compound and transport actions induced by chemical reaction, dynamic and electrodynamic processes are the main effect factors to make the ionosphere express local characteristics at the latitude, longitude and altitude [29,30]. The linear dependence between electron density and altitude may be broken down, which means some biases will be brought when the NmF2 are inferred to Ne at satellite altitude based on the empirical model. Therefore, in the Ne comparison of the CSES and ionosonde, we just paid attention to the data analysis in the daytime.
The detailed data process of the selected comparison dataset is as follows. On the one hand, the mean value (M) and standard deviation (std) of the CSES Ne falling in the region of ±10° around the ionosonde station were calculated to represent the Ne data above the ionosonde station. As the CSES is a circular sun-synchronous satellite, the longitudinal distance of each orbit in a day is about 24°. The reason for selecting ±10° as the range is to provide a sufficient number of observations and temporally continuous data to make the comparison statistically meaningful, and to avoid the ionospheric variation on a large scale. On the other hand, the IRI model was used to infer the Ne data at satellite altitude based on the ionosonde observations. The topside density profile can be obtained by fitting data to a Chapman-α layer [31], while the selection of ionospheric scale height (H) can also affect the fitting result [32,33]. As the electron density is assumed to be mostly linear-dependent at altitudes above the hmF2 [34,35], the ratio outputted by the IRI model was applied to obtain the Ne at satellite orbit in this study. The IRI model is an empirical standard model of the ionosphere, which includes the main ionospheric functions and available data sources [36], and is updated yearly according to the analysis of the IRI working group (http://irimodel.org/, accessed on 1 September 2021). There are three steps to obtaining the comparison dataset for ionosonde data. Firstly, the ratio of Ne at 507 km and NmF2 at the ionosonde location was outputted with a 15-min time resolution based on the IRI 2016 model. Secondly, the NmF2 observed by the ionosonde was multiplied by the above ratio to obtain the Ne data at the satellite altitude. Finally, the mean value and standard deviation of the calculated Ne from 13:00 to 15:00 (LT) were acquired to represent the ionosonde data in the daytime.

4.1. Time Series Analysis

As the first three months after the CSES launch is within the commissioning phase, the data provided by the satellite mission are only taken from 8 May 2018. Since the date of the available data is different, the time range of the comparative dataset is different for each ionosonde—Leshan station from 8 May 2018 to 13 October 2020, Tengchong station from 1 January 2019 to 31 July 2021, and Puer station from 8 May 2018 to 31 July 2021. Figure 5 shows the time series of the CSES Ne and Ne at the satellite altitude inferred by ionosonde observations and the IRI model, which is displayed with the log scale. The patterns of two curves are almost same, and the annual and semiannual variations can both be found for the two datasets, with high values in the summer and low values in the winter. Furthermore, the correlation coefficient between the two datasets was calculated according to Equation (3), where X and Y represent the Ne of the CSES and ionosonde. The scatters of the two datasets and linear fitting curves are plotted in Figure 6. The calculated values for the Leshan, Tengchong, and Puer stations are 0.4858, 0.5973 and 0.4944, respectively, which means the Ne of the CSES and ionosonde represents a medium correlation and can be well linearly fitted. In Figure 5, it can be seen that there is a difference in the Ne absolute value between the CSES and the ionosonde. The relative change between the CSES and ionosonde Ne was calculated according to Equation (1), where Ne_QJISR was changed to Ne_ionosonde. Figure 7 shows a histogram of the relative change between the two datasets for the Leshan Figure 7a), Tengchong (Figure 7b) and Puer (Figure 7c) stations. All of the ratios are negative, which means the electron density of the CSES is lower than the Ne at the satellite altitude inferred by ionosonde observations and the IRI model. Based on the median values of the ratios in each station, the value of CSES Ne is lower by 80–90% in the daytime, corresponding to the ionosonde Ne.

4.2. Latitudinal Analysis

According to the above data process illustrated at the beginning of Section 4, the comparison dataset between the Ne of the CSES and ionosonde data from CRIRP was constructed in March, June, September and December of 2019. Taking Lhasa station as an example, Figure 8 gives the time series of the two data, and the median value of relative change (calculated by Equation (1)) and correlation coefficient (obtained by Equation (3)) between the CSES and ionosonde Ne are also marked in each panel for different months. The variation patterns of the two datasets are almost consistent, and the peaks or troughs coincide with each other at times, for example, the similar trend seen in mid-September and the single peak in December. The correlation coefficients all exceed 0.4, and the value exceeds 0.7 in June. The value of the CSES Ne is lower than that of the ionosonde Ne with a ratio of 80–90%.
In the same way as for the Lhasa station, the correlation coefficients between the CSES and ionosonde Ne for each station from CRIRP were calculated according to Equation (3) for March, June, September and December of 2019, which was summarized in Table 1. For a more intuitive representation, the differential correlation coefficients by 1° × 1° were obtained based on all of the data in Table 1; the contour lines with colors are shown in Figure 9, where the name of each station is used the abbreviation. It can be seen that the values at the low to mid-low latitudes are mostly higher than those at the middle latitude, which is more obvious in June, September and December. In general, the correlation between the two datasets in June and September is relatively better than that in other months, with the highest value in June for the CC, LZ, LS and CQ stations, and in September for the BJ, XA and KM stations.
The median value of the relative change of the CSES Ne corresponding to the ionosonde observations (calculated by Equation (1)) for each station is summarized in Table 2. All of the data are negative, which means that the value of the CSES Ne is lower than that of the Ne at the satellite altitude inferred by ionosonde observations and the IRI model at different latitudes. The values are from −80% to −95%. In general, the difference between the two datasets is a little larger in December compared to other months, when the correlation coefficient is relatively low.

5. Discussion

Figure 5 shows that the day-to-day variation can be found in the observations of both the CSES and ionosonde Ne, and the scaled ionograms of the ionosonde can also be affected by the phenomena of spread F, the stratification of the F2 layer, sporadic E layer (Es) and so on [37,38,39], which can induce differences between the CSES and ionosonde Ne. Meanwhile, the disturbances of the running mean with 15 days (shown by red and blue lines in Figure 5) are not frequent. The correlation coefficients of the running mean values between the two datasets were calculated according to Equation (3), and are listed in Table 3 for these three stations. The coefficients increase compared to the former analysis, with values up to 0.75, especially for the Leshan and Tengchong stations. Furthermore, the ionospheric disturbances can also be induced by solar activity, magnetic storms, and other space events [30]. In order to analyze the correlation on a geomagnetically quiet day, the coefficients between the CSES and ionosonde Ne were only calculated for the data in the days with 10.7 cm solar radio flux (F10.7) < 160, Kp ≤ 3, |Dst| < 30 nT and AE < 500 nT, as shown in Table 3. Compared to the coefficient calculated by all data, the correlation between the CSES and ionosonde Ne is better on geomagnetically quiet days, which indicates that the consistency of the different observations for the ionosphere can be affected by geomagnetic events. It can be seen that the day-to-day variation and disturbance events could result in differences in the correlation analysis.
In terms of time, the correlation coefficient in June and September is generally better than that in March and December. In Figure 5, the values of the CSES Ne in June and September of 2019 are all higher than those in March and December for three ionosonde stations. It seems that the CSES Ne in the month with a relative high value may exhibit good correlation with the ionosonde observations. In general, the signal-to-noise ratio (SNR) of the payload will be better under the conditions of the observation parameters with a relatively high value. In terms of latitude, the correlation between the CSES and ionosonde Ne is better at the low latitude than that at the middle latitude during the daytime. The spatial distributions of the monthly mean for the CSES Ne in March, June, September and December of 2019 are exhibited in Figure 10. The patterns of the CSES Ne in these four months all show the high values at low latitude and low values at middle latitude. The region of the Ne peak around EIA occupies a wider area with stronger strength in March, June and September, and the values at mid to low latitude are relatively high in June. Compared to the pattern of CSES Ne distribution, it can be seen that the high value of the correlation coefficient may have a relationship with the EIA phenomenon. Affected by the fountain effect during the daytime, the phenomenon of EIA can be found at the low latitude for both the CSES and ionosonde; thus, the correlation coefficient is high around the peak of EIA. Compared to other locations, Lomidze et al. [11] also found that the initial Ne observed by Swarm agrees relatively well with COSMIC data near the magnetic equator in the daytime. Yan et al. [6] reported that the phenomena of EIA and longitudinal wave number around the magnetic equator can be observed by both the CSES and the Swarm B satellite. In the comparison between the CSES Ne and the IRI model, Liu et al. [7] detected that the correlation coefficient in the daytime is higher around the peak of EIA than that in other regions.
Compared with the QJISR data and ionosonde observations, the average values of the CSES Ne are relatively lower, with the relative change about 80–95% in the daytime corresponding to the ground-based data. Although the IRI model can exactly exhibit the ionospheric characteristics, some researchers found that the value of Ne at the topside ionosphere are overestimated by the IRI model [32,40]. The systematic biases are also found in the comparison of electron density between the CSES and Swarm constellation [5,6], DEMETER satellite and the IRI model [7], which may be caused by the different designs of the payload, different data processing and the overestimation of the IRI model at the topside ionosphere. Furthermore, similar to the above illustrations, there may be some biases in the interpretation of ionograms when the phenomena of spread F, stratification of F2 layer and Es occur in the ionosphere, which may induce the discrepancies of Ne between the CSES and ionosonde observations. Furthermore, the reliability of ISR data is widely accepted, although it is somewhat model-dependent. Different from ground-based observations for the topside ionosphere, the in situ datum of the satellite is the direct measurement of ionospheric parameters. The payload of LAP operates as an absolute instrument on the platform of the CSES, and is able to identify the plasma potential, thus retrieving the actual electron density, which is a valuable measurement for ionospheric in situ observations. Through the comparison, scientists found that the accuracy of LAP installed on a satellite could be different from mission to mission [21,41]. Although there are differences in the Ne values between ground-based remote sensing and satellite in situ measurements, similar variations and relatively good correlation between the two datasets can certify the authenticity of in situ Ne observed by the CSES.

6. Conclusions

Based on the comparison between the CSES Ne and ground-based observations in China, some conclusions were summarized from the variation trend, correlation and values.
  • The pattern of the CSES Ne is consistent with that of the ISR Ne at the satellite orbit. The two datasets can be linearly fitted, seemly according to the magnitude, and the correlation coefficient is above 0.88. However, more data are needed to more deeply understand the relationship between the two sets of data.
  • The trend of the CSES Ne is mostly similar to the ionosonde observations no matter whether the data are continuous in time or in a certain month; some correlation coefficients between the two datasets exceeded 0.7. On the China mainland, the correlation coefficient is higher around the peak of EIA compared to other locations.
  • The value of the CSES Ne is low at 80–95% in the daytime, corresponding to the ISR and ionosonde observations in China.

Author Contributions

J.L. analyzed the CSES Ne data and some ionosonde observations, and wrote the manuscript. T.X. collected and analyzed the ionosonde data at different latitudes. Z.D. collected and analyzed the ISR data. X.Z. proposed the topic and conceived the study. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Specialized Research Funds of the National Key Laboratory of Electromagnetic Environment (202001002), APSCO Earthquake Research Project Phase II, National Basic Scientific Research Program of China (JCKY2021210C614240302), Civil Project (ZH-1-DMYZ-02) and ISSI-Beijing (2019-IT33).

Data Availability Statement

The CSES Ne data can be downloaded from the website http://www.leos.ac.cn/, accessed on 10 August 2021, after registration.

Acknowledgments

This work is supported by the Specialized Research Funds of the National Key Laboratory of Electromagnetic Environment (202001002), APSCO Earthquake Research Project Phase II, National Basic Scientific Research Program of China (JCKY2021210C614240302), Civil Project (ZH-1-DMYZ-02) and ISSI-Beijing (2019-IT33). The authors acknowledge the CSES mission, and the Ne data can be downloaded from the website http://www.leos.ac.cn/, accessed on 10 August 2021, after registration. The authors acknowledge the IRI website (https://ccmc.gsfc.nasa.gov/modelweb/models/iri2016_vitmo.php, accessed on 1 September 2021) for providing the electron density data. The ionosonde observations from southwest China were provided by the Institute of Earthquake Forecasting, China Earthquake Administration. The ISR data and ionosonde observations at different latitudes were supported by the China Research Institute of Radio Wave Propagation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shen, X.; Zhang, X.; Yuan, S.; Wang, L.; Cao, J.; Huang, J.; Zhu, X.; Piergiorgio, P.; Dai, J. The state-of-the-art of the China Seismo-Electromagnetic Satellite mission. Sci. China Technol. Sci. 2018, 61, 634–642. [Google Scholar] [CrossRef]
  2. Zhou, B.; Cheng, B.; Gou, X.; Li, L.; Zhang, Y.; Wang, J.; Magnes, W.; Lammegger, R.; Pollinger, A.; Ellmeier, M.; et al. First in-orbit results of the vector magnetic field measurement of the High Precision Magnetometer onboard the China Seismo-Electromagnetic Satellite. Earth Planets Space 2019, 71, 119. [Google Scholar] [CrossRef]
  3. Diego, P.; Huang, J.; Piersanti, M.; Badoni, D.; Zeren, Z.; Yan, R.; Rebustini, G.; Ammendola, R.; Candidi, M.; Guan, Y.-B.; et al. The Electric Field Detector on Board the China Seismo Electromagnetic Satellite—In-Orbit Results and Validation. Instruments 2020, 5, 1. [Google Scholar] [CrossRef]
  4. Zhao, S.; Zhou, C.; Shen, X.; Zhima, Z. Investigation of VLF Transmitter Signals in the Ionosphere by ZH-1 Observations and Full-Wave Simulation. J. Geophys. Res. Space Phys. 2019, 124, 4697–4709. [Google Scholar] [CrossRef]
  5. Wang, X.; Cheng, W.; Yang, D.; Liu, D. Preliminary validation of in situ electron density measurements onboard CSES using observations from Swarm Satellites. Adv. Space Res. 2019, 64, 982–994. [Google Scholar] [CrossRef]
  6. Yan, R.; Zhima, Z.; Xiong, C.; Shen, X.; Huang, J.; Guan, Y.; Zhu, X.; Liu, C. Comparison of Electron Density and Temperature From the CSES Satellite With Other Space-Borne and Ground-Based Observations. J. Geophys. Res. Space Phys. 2020, 125, e2019JA027747. [Google Scholar] [CrossRef]
  7. Liu, J.; Guan, Y.; Zhang, X.; Shen, X. The data comparison of electron density between CSES and DEMETER satellite, Swarm constellation and IRI model. Earth Space Sci. 2021, 8, e2020EA001475. [Google Scholar] [CrossRef]
  8. Picozza, P.; Battiston, R.; Ambrosi, G.; Bartocci, S.; Basara, L.; Burger, W.J.; Campana, D.; Carfora, L.; Casolino, M.; Castellini, G.; et al. Scientific Goals and In-orbit Performance of the High-energy Particle Detector on Board the CSES. Astrophys. J. Suppl. Ser. 2019, 243, 16. [Google Scholar] [CrossRef]
  9. Ambrosi, G.; Bartocci, S.; Basara, L.; Battiston, R.; Burger, W.J.; Campana, D.; Caprai, M.; Carfora, L.; Castellini, G.; Cipollone, P.; et al. The electronics of the High-Energy Particle Detector on board the CSES-01 satellite. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2021, 1013, 165639. [Google Scholar] [CrossRef]
  10. Wang, X.; Yang, D.; Zhou, Z.; Cheng, W.; Xu, S.; Guo, F. Validation of CSES RO measurements using ionosonde and ISR observations. Adv. Space Res. 2020, 66, 2275–2288. [Google Scholar] [CrossRef]
  11. Lomidze, L.; Knudsen, D.J.; Burchill, J.; Kouznetsov, A.; Buchert, S.C. Calibration and Validation of Swarm Plasma Densities and Electron Temperatures Using Ground-Based Radars and Satellite Radio Occultation Measurements. Radio Sci. 2018, 53, 15–36. [Google Scholar] [CrossRef]
  12. Cai, X.; Burns, A.G.; Wang, W.; Coster, A.; Qian, L.; Liu, J.; Solomon, S.C.; Eastes, R.W.; Daniell, R.E.; McClintock, W.E. Comparison of GOLD Nighttime Measurements With Total Electron Content: Preliminary Results. J. Geophys. Res. Space Phys. 2020, 125, e2019JA027767. [Google Scholar] [CrossRef]
  13. Harding, B.J.; Chau, J.L.; He, M.; Englert, C.R.; Harlander, J.M.; Marr, K.D.; Makela, J.J.; Clahsen, M.; Li, G.; Ratnam, M.V.; et al. Validation of ICON-MIGHTI thermospheric wind observations: 2. Green-line comparisons to specular meteor radars. J. Geophys. Res. Space Phys. 2021, 126, e2020JA028947. [Google Scholar] [CrossRef]
  14. Pedatella, N.M.; Zakharenkova, I.; Braun, J.J.; Cherniak, I.; Hunt, D.; Schreiner, W.S.; Straus, P.R.; Valant-Weiss, B.L.; Vanhove, T.; Weiss, J.; et al. Processing and Validation of FORMOSAT-7/COSMIC-2 GPS Total Electron Content Observations. Radio Sci. 2021, 56, e2021RS007267. [Google Scholar] [CrossRef]
  15. Catapano, F.; Buchert, S.; Qamili, E.; Nilsson, T.; Bouffard, J.; Siemes, C.; Coco, I.; D’Amicis, R.; Tøffner-Clausen, L.; Trenchi, L.; et al. Swarm Langmuir Probes’ data quality and future improvements. Geosci. Instrum. Methods Data Syst. 2021. In discussions. [Google Scholar]
  16. Singh, A.K.; Maltseva, O.; Panda, S.K. Comparison between Swarm measured and IRI-2016, IRI-Plas 2017 modeled electron density over low and mid latitude region. Acta Astronaut. 2021, 189, 476–482. [Google Scholar] [CrossRef]
  17. Shen, X.H.; Zong, Q.G.; Zhang, X.M. Introduction to special section on the China Seismo-Electromagnetic Satellite and initial results. Earth Planet. Phys. 2018, 2, 439–443. [Google Scholar] [CrossRef]
  18. Guan, Y.-B.; Wang, S.-J.; Liu, C.; Feng, Y.-B. The design of the Langmuir probe onboard a seismo-electromagnetic satellite. In International Symposium on Photoelectronic Detection and Imaging 2011: Space Exploration Technologies and Applications; SPIE: Bellingham, WA, USA, 2011; pp. 631–639. [Google Scholar]
  19. Yan, R.; Guan, Y.; Shen, X.; Huang, J.; Zhang, X.; Liu, C.; Liu, D. The Langmuir Probe Onboard CSES: Data inversion analysis method and first results. Earth Planet. Phys. 2018, 2, 479–488. [Google Scholar] [CrossRef]
  20. Liu, C.; Guan, Y.; Zheng, X.; Zhang, A.; Piero, D.; Sun, Y. The technology of space plasma in-situ measurement on the China Seismo-Electromagnetic Satellite. Sci. China Technol. Sci. 2018, 62, 829–838. [Google Scholar] [CrossRef]
  21. Mcnamara, L.F.; Cooke, D.L.; Valladares, C.E.; Reinisch, B.W. Comparison of CHAMP and Digisonde plasma frequencies at Jicamarca, Peru. Radio Sci. 2007, 42, RS2005. [Google Scholar] [CrossRef]
  22. Hu, L.; Ning, B.; Liu, L.; Zhao, B.; Li, G.; Wu, B.; Huang, Z.; Hao, X.; Chang, S.; Wu, Z. Validation of COSMIC ionospheric peak parameters by the measurements of an ionosonde chain in China. Ann. Geophys. 2014, 32, 1311–1319. [Google Scholar] [CrossRef] [Green Version]
  23. Ding, Z.; Yu, L.; Dai, L.; Xu, Z.; Wu, J. The preliminary measurement and analysis of the power profiles by the Qujing incoherent scatter radar. Chin. J. Geophys. 2014, 57, 3564–3569, (In Chinese with English abstract). [Google Scholar]
  24. Liu, J.; Jiang, C.; Deng, C.; Yang, G.; Zhang, X.; Lou, W.; Yang, C. Vertical ionosonde net and its data application in southwestern China. Acta Seismol. Sin. 2016, 38, 399–407, (In Chinese with English abstract). [Google Scholar]
  25. Xu, T.; Wu, J.; Zhao, Z.; Liu, Y.; He, S.; Li, J.; Wu, Z.; Hu, Y. Monitoring ionospheric variations before earthquakes using the vertical and oblique sounding network over China. Nat. Hazards Earth Syst. Sci. 2011, 11, 1083–1089. [Google Scholar] [CrossRef]
  26. Bhonsle, R.V.; da Rosa, A.V.; Garriott, O.K. Meas- urements of the total electron content and equivalent slab thickness of the mid-latitude ionospher. Radio Sci. 1965, 69D, 929. [Google Scholar]
  27. McNamara, L.F.; Smith, D.H. Total electron content of the ionosphere at 30°S, 1967-1974. J. Atmos. Terr. Phys. 1982, 44, 227. [Google Scholar] [CrossRef]
  28. Liu, L.; Wan, W.; Yue, X.; Zhao, B.; Ning, B.; Zhang, M.L. The dependence of plasma density in the topside ionosphere on the solar activity level. Ann. Geophys. 2007, 25, 1337–1343. [Google Scholar] [CrossRef]
  29. Balan, N.; Rao, P.B. Latitudinal variations of nighttime enhancements in total electron content. J. Geophys. Res. 1987, 92, 3436–3440. [Google Scholar] [CrossRef]
  30. Robert, S.; Andrew, N. Ionospheres: Physics, Plasma Physics, and Chemistry; Cambridge University: New York, NY, USA, 2009. [Google Scholar]
  31. Rishbeth, H.; Garriott, O.K. Introduction to Ionospheric Physics; Elsevier: New York, NY, USA, 1969. [Google Scholar]
  32. Lei, J.; Liu, L.; Wan, W.; Zhang, S.-R. Variations of electron density based on long-term incoherent scatter radar and ionosonde measurements over Millstone Hill. Radio Sci. 2005, 40, RS2008. [Google Scholar] [CrossRef]
  33. Liu, L.; Wan, W.; Zhang, M.-L.; Ning, B.; Zhang, S.-R.; Holt, J.M. Variations of topside ionospheric scale heights over Millstone Hill during the 30-day incoherent scatter radar experiment. Ann. Geophys. 2007, 25, 2019–2027. [Google Scholar] [CrossRef]
  34. Lei, J.; Liu, L.; Wan, W.; Zhang, S.-R.; Van Eyken, A.P. Comparison of the first long-duration IS experiment measurements over Millstone Hill and EISCAT Svalbard radar with IRI2001. Adv. Space Res. 2006, 37, 1102–1107. [Google Scholar] [CrossRef]
  35. Liu, L.; Le, H.; Wan, W.; Sulzer, M.P.; Lei, J.; Zhang, M.-L. An analysis of the scale heights in the lower topside ionosphere based on the Arecibo incoherent scatter radar measurements. J. Geophys. Res. Space Phys. 2007, 112, A06307. [Google Scholar] [CrossRef]
  36. Bilitza, D. (Ed.) International Reference Ionosphere 1990, NSSDC 90-22; Science Applications Research: Greenbelt, MD, USA, 1990. [Google Scholar]
  37. Wei, L.; Jiang, C.; Hu, Y.; Aa, E.; Huang, W.; Liu, J.; Yang, G.; Zhao, Z. Ionosonde Observations of Spread F and Spread Es at Low and Middle Latitudes during the Recovery Phase of the 7–9 September 2017 Geomagnetic Storm. Remote Sens. 2021, 13, 1010. [Google Scholar] [CrossRef]
  38. Jiang, C.; Yang, G.; Liu, J.; Yokoyama, T.; Komolmis, T.; Song, H.; Lan, T.; Zhou, C.; Zhang, Y.; Zhao, Z. Ionosonde observations of daytime spread F at low latitudes. J. Geophys. Res. Space Phys. 2016, 121, 12093–12103. [Google Scholar] [CrossRef]
  39. Jiang, C.; Hu, H.; Yang, G.; Liu, J.; Zhao, Z. A statistical study of the F2 layer stratification at the northern equatorial ionization anomaly. Adv. Space Res. 2019, 63, 3167–3176. [Google Scholar] [CrossRef]
  40. Lühr, H.; Xiong, C. IRI-2007 model overestimates electron density during the 23/24 solar minimum. Geophys. Res. Lett. 2010, 37. [Google Scholar] [CrossRef]
  41. Pedatella, N.M.; Yue, X.; Schreiner, W.S. Comparison between GPS radio occultation electron densities and in situ satellite observations. Radio Sci. 2015, 50, 518–525. [Google Scholar] [CrossRef]
Figure 1. Locations of ground-based observation stations in China. Blue points represent ionosonde stations of IEF with continuous-time data. Red triangles indicate ionosonde stations of CRIRP at different latitudes. Green star means the location of the ISR station. Green lines are the footprints of the CSES in one day.
Figure 1. Locations of ground-based observation stations in China. Blue points represent ionosonde stations of IEF with continuous-time data. Red triangles indicate ionosonde stations of CRIRP at different latitudes. Green star means the location of the ISR station. Green lines are the footprints of the CSES in one day.
Remotesensing 14 04498 g001
Figure 2. Ne time series for CSES (red line) and QJISR (blue line). Each day, the circle means mean value, and the bar represents standard deviation. The two time series are displayed with different amplitude scales to better underline their own variations. The comparison days are not continuous and the detailed dates are marked on the X axis.
Figure 2. Ne time series for CSES (red line) and QJISR (blue line). Each day, the circle means mean value, and the bar represents standard deviation. The two time series are displayed with different amplitude scales to better underline their own variations. The comparison days are not continuous and the detailed dates are marked on the X axis.
Remotesensing 14 04498 g002
Figure 3. Histogram for relative change of CSES Ne corresponding to QJISR Ne. The Y axis represents empirical cumulative distribution function (cdf).
Figure 3. Histogram for relative change of CSES Ne corresponding to QJISR Ne. The Y axis represents empirical cumulative distribution function (cdf).
Remotesensing 14 04498 g003
Figure 4. Ne scatters of CSES and QJISR with two groups (panel (a)) and the fitting lines (panel (b)). In panel (a), the number means the detailed date for each point with the format of YYMMDD. In panel (a), the data of groups 1 and 2 are represented by blue and black dots, respectively.
Figure 4. Ne scatters of CSES and QJISR with two groups (panel (a)) and the fitting lines (panel (b)). In panel (a), the number means the detailed date for each point with the format of YYMMDD. In panel (a), the data of groups 1 and 2 are represented by blue and black dots, respectively.
Remotesensing 14 04498 g004
Figure 5. Time series of Ne for CSES and ionosonde. Panels (ac) are for Leshan, Tengchong and Puer stations, respectively. In each panel, the red dotted line represents CSES Ne, and the blue dotted line denotes Ne at the satellite altitude inferred by ionosonde observations and IRI model. Red and blue lines are the 15-day running mean.
Figure 5. Time series of Ne for CSES and ionosonde. Panels (ac) are for Leshan, Tengchong and Puer stations, respectively. In each panel, the red dotted line represents CSES Ne, and the blue dotted line denotes Ne at the satellite altitude inferred by ionosonde observations and IRI model. Red and blue lines are the 15-day running mean.
Remotesensing 14 04498 g005
Figure 6. Ne scatters for CSES and ionosonde. Panels (ac) are for Leshan, Tengchong and Puer stations, respectively. In each panel, the red line is the linear fitting line.
Figure 6. Ne scatters for CSES and ionosonde. Panels (ac) are for Leshan, Tengchong and Puer stations, respectively. In each panel, the red line is the linear fitting line.
Remotesensing 14 04498 g006
Figure 7. Relative change of CSES Ne corresponding to ionosonde Ne. Panels (ac) are for Leshan, Tengchong and Puer, respectively. In each panel, the Y axis represents empirical cumulative distribution function (cdf).
Figure 7. Relative change of CSES Ne corresponding to ionosonde Ne. Panels (ac) are for Leshan, Tengchong and Puer, respectively. In each panel, the Y axis represents empirical cumulative distribution function (cdf).
Remotesensing 14 04498 g007
Figure 8. Time series of Ne for CSES and Lhasa ionosonde in different months. In each panel, the red curve represents CSES Ne and the blue curve denotes Ne at the satellite altitude inferred by Lhasa ionosonde observations and IRI model. The two time series are displayed with a log scale.
Figure 8. Time series of Ne for CSES and Lhasa ionosonde in different months. In each panel, the red curve represents CSES Ne and the blue curve denotes Ne at the satellite altitude inferred by Lhasa ionosonde observations and IRI model. The two time series are displayed with a log scale.
Remotesensing 14 04498 g008
Figure 9. Contour map of correlation coefficient between CSES and ionosonde Ne in different months. Acronyms of station names are used in each panel.
Figure 9. Contour map of correlation coefficient between CSES and ionosonde Ne in different months. Acronyms of station names are used in each panel.
Remotesensing 14 04498 g009
Figure 10. Spatial distribution of CSES Ne in different months for 2019.
Figure 10. Spatial distribution of CSES Ne in different months for 2019.
Remotesensing 14 04498 g010
Table 1. Correlation coefficient between CSES and ionosonde Ne at different ionosonde stations.
Table 1. Correlation coefficient between CSES and ionosonde Ne at different ionosonde stations.
StationMarchJuneSeptemberDecember
Manzhouli (MZL)0.53310.21920.1460−0.2932
Changchun (CC)0.24100.52430.4167−0.0542
Urumuchi (UC)0.30170.2121−0.0263−0.0101
Beijing (BJ)0.48950.21930.58640.1550
Lanzhou (LZ)0.40180.64450.55850.4069
Xi’an (XA)0.39860.51190.66650.0989
Lhasa (LS)0.44760.71240.51530.6356
Chongqing (CQ)0.54440.75330.53660.3109
Kunming (KM)0.39390.57860.73110.3447
Haikou (HK)0.35520.55240.52280.6494
Table 2. Median value for relative change of CSES Ne corresponding to ionosonde Ne at different stations.
Table 2. Median value for relative change of CSES Ne corresponding to ionosonde Ne at different stations.
StationMarchJuneSeptemberDecember
Manzhouli (MZL)−93.51%−90.35%−94.48%−94.22%
Changchun (CC)−85.52%−80.47%−84.62%−91.08%
Urumuchi (UC)−93.81%−88.35%−93.46%−94.29%
Beijing (BJ)−93.13%−87.11%−93.50%−94.10%
Lanzhou (LZ)−90.48%−85.00%−89.99%−93.45%
Xi’an (XA)−89.72%−84.45%−90.10%−92.86%
Lhasa (LS)−87.09%−82.18%−86.02%−91.17%
Chongqing (CQ)−85.52%−80.47%−84.62%−91.08%
Kunming (KM)−88.72%−84.31%−89.52%−91.91%
Haikou (HK)−85.29%−81.85%−87.13%−93.21%
Table 3. Correlation coefficient between CSES and ionosonde Ne for Leshan, Tengchong and Puer stations in different conditions.
Table 3. Correlation coefficient between CSES and ionosonde Ne for Leshan, Tengchong and Puer stations in different conditions.
StationRaw DataRunning MeanGeomagnetically Quiet Day
Leshan0.48580.75570.5785
Tengchong0.59730.75060.6334
Puer0.49440.64250.5431
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, J.; Xu, T.; Ding, Z.; Zhang, X. The Comparison of Electron Density between CSES In Situ and Ground-Based Observations in China. Remote Sens. 2022, 14, 4498. https://doi.org/10.3390/rs14184498

AMA Style

Liu J, Xu T, Ding Z, Zhang X. The Comparison of Electron Density between CSES In Situ and Ground-Based Observations in China. Remote Sensing. 2022; 14(18):4498. https://doi.org/10.3390/rs14184498

Chicago/Turabian Style

Liu, Jing, Tong Xu, Zonghua Ding, and Xuemin Zhang. 2022. "The Comparison of Electron Density between CSES In Situ and Ground-Based Observations in China" Remote Sensing 14, no. 18: 4498. https://doi.org/10.3390/rs14184498

APA Style

Liu, J., Xu, T., Ding, Z., & Zhang, X. (2022). The Comparison of Electron Density between CSES In Situ and Ground-Based Observations in China. Remote Sensing, 14(18), 4498. https://doi.org/10.3390/rs14184498

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