Ionospheric Electron Density Distribution Research Articles

Investigation on Chasing and Interaction of Traveling Ionospheric Disturbances Based on Multi‐Instrument

AbstractIn this study, we use multi‐instrument observations (all‐sky imager (ASI), global navigation satellite system (GPS) receivers, digisonde) to study the interaction of nighttime medium‐scale traveling ionospheric disturbances (MSTIDs) on 13 November 2018. The most attractive aspect of this event is that the interaction appeared between two dark bands both propagated southwestward. The airglow observations show that the latter band moved faster and caught up with the former, and these two bands merged into a new one. The propagating characteristics and morphology of the MSTIDs changed during the interaction process. The simulations from the Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model (TIEGCM) suggested that the ionospheric background zonal winds and electron density distributions could play essential roles in the interaction of the MSTIDs. Moreover, the merging process might be associated with the electrostatic reconnection.

Role of Martian Crustal Fields in Ionospheric Electron Density Distribution and Subsequent South‐North Asymmetry: Insights From Multi‐Year MAVEN Observations During (MYs 33–36)

AbstractThis study uses Mars Atmosphere and Volatile EvolutioN observations of electron density and magnetic field for a period of four Martian years (MYs 33–36) (∼8 Earth years) to investigate the effects of Martian crustal magnetic fields on the distribution and variability of Mars' ionosphere. The results show a clear enhancement in electron density in the southern hemisphere in the region where the strong crustal magnetic fields are present with the longitudes between 120° and 240° (i.e., the central longitude), which is in agreement with previous studies. On the contrary, the corresponding northern hemisphere region in the central longitudes shows an exactly opposite behavior that the electron density is lower compared to the surrounding longitude regions. These effects are found to be primarily dayside phenomena. As opposed to dayside, the nightside electron density in the central longitudes are slightly reduced at altitudes below 200 km, compared to longitudes on its western and eastern sides. Above 200 km, the nightside effects are not very clear. Significant hemispheric asymmetry is observed in the longitude regions of enhanced crustal magnetic fields compared to other longitude regions during the daytime. This dayside south‐north asymmetry in the central longitude region is observed to be a constant feature across all seasons. However, on the nightside, the south‐north asymmetry remains more or less similar across all longitude regions, during all seasons implying a weakened control of the crustal fields over the nightside ionosphere. Even then, the northern hemisphere retains a stronger nightside ionosphere during all seasons except summer.

Application of generalized aurora computed tomography to the EISCAT_3D project

Abstract. EISCAT_3D is a project to build a multi-site phased-array incoherent scatter radar system in northern Fenno-Scandinavia. We demonstrate via numerical simulation how useful monochromatic images taken by a multi-point imager network are for auroral research in the EISCAT_3D project. We apply the generalized aurora computed tomography (G-ACT) method to modelled observational data from real instruments, such as the Auroral Large Imaging System (ALIS) and the EISCAT_3D radar. G-ACT is a method for reconstructing the three-dimensional (3D) distribution of auroral emissions and ionospheric electron density (corresponding to the horizontal two-dimensional (2D) distribution of energy spectra of precipitating electrons) from multi-instrument data. It is assumed that the EISCAT_3D radar scans an area of 0.8° in geographic latitude and 3° in longitude at an altitude of 130 km with 10 × 10 beams from the radar core site at Skibotn (69.35° N, 20.37° E). Two neighboring discrete arcs are assumed to appear in the observation region of the EISCAT_3D radar. The reconstruction results from G-ACT are compared with those from the normal ACT as well as the ionospheric electron density from the radar. It is found that G-ACT can interpolate the ionospheric electron density at a much higher spatial resolution than that observed by the EISCAT_3D radar. Furthermore, the multiple arcs reconstructed by G-ACT are more precise than those by ACT. In particular, underestimation of the ionospheric electron density and precipitating electrons' energy fluxes inside the arcs is significantly improved by G-ACT including the EISCAT_3D data. Even when the ACT reconstruction is difficult due to the unsuitable locations of the imager sites relative to the discrete arcs and/or a small number of available images, G-ACT allows us to obtain better reconstruction results.

D-region ionospheric disturbances due to the December 2019 solar eclipse observed using multi-station VLF radio network

We present the low-latitude D-region (60–90 km) ionospheric disturbances due to the December 2019 annular solar eclipse observed using the multi-station Very Low Frequency (VLF) radio network in West Bengal, India. VLF signals from communication transmitters were received from ten places using low-cost VLF radio receivers. The receivers were capable of recording the amplitudes of VLF signals mainly from the transmitters VTX, India at 18.2 kHz and NWC, Australia at 19.8 kHz. During the solar eclipse, the VTX signal experienced considerable positive, negative, and mixed-type amplitude variations depending on the propagation lengths between the transmitter and receivers, whereas the NWC signal displayed a relatively small amplitude response. We have also found that the highest VLF signal disturbances in each propagation path occurred after the time when the solar obscuration was maximum over the entire path. We used the observations of the VTX signal to formulate the relationship between solar obscuration and D-region electron density reduction during the solar eclipse. On the basis of this, we presented the electron density distribution of the D-region ionosphere during the solar eclipse across the Indian subcontinent. This study shows a useful method for integrating VLF observations and simulations to estimate the D-region electron density distribution over a vast area during a solar eclipse.

Prediction of Ionospheric Electron Density Distribution Based on CNN-LSTM Model

Prediction of Ionospheric Electron Density Distribution Based on CNN-LSTM Model

Computerized Ionospheric Tomography Based on the ADS-B System

The broadcast automatic dependent surveillance (ADS-B) system is a new-generation air traffic control system designed to avoid the waste of resources in secondary radars. The establishment of the spaceborne ADS-B system provides a broad prospect for ionospheric tomography. In this paper, the external observation information of the ionosphere is obtained by measuring the Faraday rotation angle, that is, the total electron content (TEC). Tomography research can be carried out all over the world to conduct large-scale ionospheric electron density research. The experiment selected two different regions and had a time resolution of two hours, a height resolution of 200 km, a latitude resolution of 2°, and a longitude resolution of 5°. Based on the simulated spaceborne ADS-B signal to invert the regional ionospheric electron density, the latitude, longitude, and height distributions of inversion result are basically consistent with those of the actual ionospheric electron density.

Optimal Estimation Inversion of Ionospheric Electron Density from GNSS-POD Limb Measurements: Part I-Algorithm and Morphology

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in characterizing the global 3D distribution and variability of ionospheric electron density (Ne). In this study, we developed an optimal estimation (OE) method to retrieve Ne profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles (ε < 0°). Although both OE and onion-peeling (OP) methods use the Abel weighting function in the Ne inversion, they are significantly different in terms of performance in the lower ionosphere. The new OE results can overcome the large Ne oscillations, sometimes negative values, seen in the OP retrievals in the E-region ionosphere. In the companion paper in this Special Issue, the HmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006–2019), COSMIC-2 (2019–present), and Spire (2019–present) constellations, showing a consistent ionospheric Ne morphology. The unprecedented spatiotemporal sampling of the ionosphere from these constellations now allows a detailed analysis of the frequency–wavenumber spectra for the Ne variability at different heights. In the lower ionosphere (~150 km), we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2, and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4, and SW4. The co-existence of eastward- and westward-propagating wave4 components implies the presence of a stationary wave4 (SPW4), as suggested by other earlier studies. Further improvements to the OE method are proposed, including a tomographic inversion technique that leverages the asymmetric sampling about the tangent point associated with GNSS-LEO links.

Ionospheric Electron Density Model by Electron Density Grid Deep Neural Network (EDG-DNN)

Electron density (or electron concentration) is a critical metric for characterizing the ionosphere’s mobility. Shortwave technologies, remote sensing systems, and satellite communications—all rely on precise estimations of electron density in the ionosphere. Using electron density profiles from FORMOSAT-3/COSMIC (Constellation Observation System for Meteorology, Ionosphere, and Climate) from 2006 to 2013, a four-dimensional physical grid model of ionospheric electron density was created in this study. The model, known as EDG-DNN, utilizes a DNN (deep neural network), and its output is the electron density displayed as a physical grid. The preprocessed electron density data are used to construct training, validation, and test sets. The International Reference Ionosphere model (IRI) was chosen as the reference model for the validation procedure since it predicts electron density well. This work used the IRI-2016 version. IRI-2016 produced more precise results of electron density when time and location parameters were input. This study compares the electron density provided by IRI-2016 to the EDG-DNN to assess the merits of the latter. The final results reveal that EDG-DNN has low-error and strong stability, can represent the global distribution structure of electron density, has some distinctive features of ionospheric electron density distribution, and predicts electron density well during quiet periods.

A Double-Adaptive Adjustment Algorithm for Ionospheric Tomography

A double-adaptive adjustment algorithm (DAAA) is proposed to reconstruct three-dimensional ionospheric electron density (IED) distribution. In the DAAA method, the relaxation factor of the multiplicative algebraic reconstruction technique (MART) is first adaptively adjusted by introducing adaptive MART (AMART). To avoid the voxels without any rays traversing them becoming dependent on the initial IED values, smoothing constraints are generally imposed on the adaptive multiplicative algebraic reconstruction technique (AMART). In general, the elements of the smoothing matrices are invariant in the iterative process. They affect the accuracy and efficiency of the IED inversion. To overcome the above limitation, the adaptive adjustments of the constrained matrix elements are subsequently carried out. Both numerical simulation and actual global navigation satellite system (GNSS) experimental results validate that the accuracy and efficiency of ionospheric tomography have been improved by the DAAA method. Finally, the new algorithm is applied to reconstruct the three-dimensional structure of the ionosphere during different geomagnetic activities. The comparisons show that the vertical profiles of the DAAA method are in agreement with those recorded from the ionosonde, and the inverted vertical total electron content (VTEC) of the DAAA method also agrees with the ionospheric products of center for orbit determination in Europe (CODE) during geomagnetic quiet and geomagnetic storms. The comparisons confirm the reliability and superiority of the DAAA method.

Three-Step Tomographic Algorithm for Ionospheric Electron Density Reconstruction

In this work, a three-step algorithm for the tomographic reconstruction of ionospheric electron density (IED) is proposed. In the new algorithm, the Taylor series expansion of the multiplicative algebraic reconstruction technique is first introduced to perform an adaptive adjustment for the relaxation parameter vector. Second, the horizontal and vertical constraints are imposed on the tomography system aiming to reduce the dependence on the initial iteration values for those voxels without any rays traversing them. Third, the uniformed voxel size is replaced with the variable voxel size by adopting the unequal grid interval in the altitudinal direction. The feasibility and superiority of the three-step algorithm are validated by devising a numerical simulation scheme. Finally, the new algorithm is successfully used to obtain the 3-D IED distribution under magnetically quiet conditions. The comparisons of IED profiles and the statistics of the reconstruction error of different algorithms further validate the superiority of the three-step algorithm.

A Tomographic Method for the Reconstruction of the Plasmasphere Based on COSMIC/ FORMOSAT-3 Data

A tomographic method has been developed for reconstructing the topside ionospheric and plasmaspheric electron density distribution using total electron content (TEC) measurements from global positioning system (GPS) receivers aboard the constellation observing system for meteorology, ionosphere, and climate&#x002F;Formosa Satellite Mission 3 (COSMIC&#x002F;FORMOSAT-3). Since the COSMIC&#x002F;FORMOSAT-3 constellation has an orbit altitude of about 800 km, the integral TEC measurements obtained from the topside GPS navigation data are rather small and imposed relevant challenges to obtaining stable electron density reconstructions. However, the developed method can represent the natural variability of the plasma ambient in terms of latitude, altitude, solar activity, season, and local time when analyzing electron density reconstructions during 2008&#x2013;2013. The method employs independent spatial grids for satellite rising and setting geometries and imposes a set of constraints to stabilize the solution in the presence of noise and ill-conditioned geometry. We further consider background ionosphere and plasmasphere models for electron density initialization and filling data gaps. The quality assessment using TEC and <i>in-situ</i> electron density measurements has shown that the proposed method performs better than the background model, with improvements of about 26&#x0025; in TEC and 20&#x0025; in terms of electron density. Our investigation also reveals the necessity of more accurate background electron density representations and precise TEC measurements in order to have better plasmaspheric specifications at high altitudes.

Performance of GPS Positioning in the Presence of Irregularities in the Auroral and Polar Ionospheres during EISCAT UHF/ESR Measurements

Irregularities in the spatial distribution of ionospheric electron density introduce temporal fluctuations in the intensity and phase of radio signals received from Global Navigation Satellite Systems (GNSS). The impact of phase fluctuations originating from irregularities in the auroral and polar ionospheres on GPS positioning was investigated on three days in March 2018 in the presence of quiet-to-moderately disturbed magnetic conditions by combining measurements from GPS and EISCAT UHF/ESR incoherent scatter radars. Two different positioning solutions were analysed: broadcast kinematic (BK) and precise static (PS). The results show that the propagation through irregularities induced residual errors on the observables leading to an increase in the positioning error, in its variability, and in the occurrence of gaps. An important aspect emerging from this study is that the variability of the 3-D positioning error was reduced, and the presence of gaps disappeared when the positioning solutions were evaluated at a 1 s rate rather than at a 30 s rate. This is due to the transient nature of residual errors that are more significant over 30 s time intervals in the presence of irregularities with scale size between few kilometres in the E region to few tens of kilometres in the F region.

Analysis of different propagation models for the estimation of the topside ionosphere and plasmasphere with an ensemble Kalman filter

Abstract. The accuracy and availability of satellite-based applications, like Global Navigation Satellite System (GNSS) positioning and remote sensing, crucially depend on the knowledge of the ionospheric electron density distribution. The tomography of the ionosphere is one of the major tools for providing links to specific ionospheric corrections and studying and monitoring physical processes in the ionosphere and plasmasphere. In this work, we apply an ensemble Kalman filter (EnKF) approach for the 4D electron density reconstruction of the topside ionosphere and plasmasphere, with the focus on the investigation of different propagation models, and compare them with the iterative reconstruction technique of simultaneous multiplicative column normalized method plus (SMART+). The slant total electron content (STEC) measurements of 11 low earth orbit (LEO) satellites are assimilated into the reconstructions. We conduct a case study on a global grid with altitudes between 430 and 20 200 km, for two periods of the year 2015, covering quiet to perturbed ionospheric conditions. Particularly the performance of the methods for estimating independent STEC and electron density measurements from the three Swarm satellites is analysed. The results indicate that the methods of EnKF, with exponential decay as the propagation model, and SMART+ perform best, providing, in summary, the lowest residuals.

A new three-dimensional computerized ionospheric tomography model based on a neural network

Computerized ionospheric tomography (CIT) is an ill-posed inverse problem owing to insufficient data acquisition. Therefore, the ionospheric electron density (IED) distributions cannot be reconstructed accurately. Although many attempts have been made to deal with this issue, there is still a long way to go before it can be completely overcome. Specifically, the inverted IEDs of voxels without observational information show a strong dependence on initial values, which affects the overall accuracy of CIT. Taking this into account, a new three-dimensional CIT model is developed, based on a backpropagation neural network. The neural network model is trained using the characteristics and inverted IEDs of voxels with observational information, and then, the IEDs of voxels without observational information are predicted again. Careful validation of the proposed model is performed by conducting numerical experiments with GPS simulation and real data under both quiet and disturbed ionospheric conditions. Compared with the traditional non-neural network method in the simulation experiment, the proposed method offers improvements of 62.0 and 56.89% in root mean square error and the mean absolute error for those voxels without observational information, respectively, while it offers improvements of 30.98 and 26.67% for all voxels of the whole region. In the real data experiment, the IEDs of the control groups obtained by the proposed method are compared with the target IEDs for all periods. The result presents correlation coefficient greater than 0.96 between this predicted IEDs and the target IEDs for all periods, and this further certifies the feasibility of the proposed method. Additionally, the latitude–longitude maps and profiles of the ionospheric electron density also show that the ill-posedness problem has a significantly weaker effect for those voxels without observational information.

Ionospheric response to the 17 March and 22 June 2015 geomagnetic storms over Wuhan region using GNSS-based tomographic technique

By using the data of GNSS (Global Navigation Satellite System) observation from Crustal Movement Observation Network of China (CMONOC), ionospheric electron density (IED) distributions reconstructed by using computerized ionospheric tomography (CIT) technique are used to investigate the ionospheric storm effects over Wuhan region during 17 March and 22 June 2015 geomagnetic storm periods. F-region critical frequency (foF2) at Wuhan ionosonde station shows an obvious decrease during recovery phase of the St. Patrick’s Day geomagnetic storm. Moreover, tomographic results present that the decrease in electron density begins at 12:00 UT on 17 March during the storm main phase. Also, foF2 shows a long-lasting negative storm effect during the recovery phase of the 22 June 2015 geomagnetic storm. Electron density chromatography presents the evident decrease during the storm day in accordance with the ionosonde observation. These ionospheric negative storm effects are probably associated with changes of chemical composition, PPEF and DDEF from high latitudes.

Virtual reference station-based computerized ionospheric tomography

In computerized ionospheric tomography (CIT) with ground-based GNSS, the voxels without satellite-receiver ray traversing cannot be reconstructed directly. We present a CIT algorithm based on virtual reference stations (VRSs), called VRS–CIT, to decrease the number of unilluminated voxels and improve the precision of the estimated ionospheric electron density (IED). The VRSs are set at the nodes of grids with a 0.5° × 0.5° resolution in longitude and latitude. We generate the virtual observations with the observations from nearby six or three stations selected according to azimuths and distances. The generation utilizes multi-quadric surface fitting with six stations and triangular linear interpolation with three stations. With the virtual observations added, the IED distribution is reconstructed by the multiplicative algebraic reconstruction technique with the initial values obtained from IRI-2016. The performance of VRS–CIT is examined by using the data from 127 GNSS stations located in 24–46° N and 122–146° E to derive the IED every 30 min. The study focuses on April 29, 2014, with the adaptability of VRS–CIT analyzed by 12 days, evenly distributed around equinoxes and solstices of 2014. The accuracy of the virtual observation is about 1 TECU. Comparing to that derived from CIT with only real observations, the unsolvability of VRS–CIT declined by 4–12% for the whole region, and for the main area, the improvement can be up to 70%. Taking two IED profiles from radio occultation as reference measurements, the mean absolute error (MAE) of IED by VRS–CIT decreases by 6.88% and 8.43%, respectively. Comparing with slant total electron content (STEC) extracted from five additional GNSS stations, the MAE and the root mean square error of the estimated STEC can be reduced up to 17.24% and 33.81%, respectively.

Imaging the Three-Dimensional Ionospheric Structure with a Blob Basis Functional Ionospheric Tomography Model.

A new ionospheric tomography model is presented in this work. In the new model, the traditional voxel basis function is replaced by the blob basis function. Due to the overlapping nature of their rotational symmetric basis functions, the new model introduces certain weighting from nearby tomographic spherical blobs. To confirm the feasibility of the new tomography model, a numerical simulation scheme is devised, and the simulation demonstrates that the reconstructed quality of the blob basis tomographic model is higher than that of the voxel basis tomographic model. Meanwhile, the variable blob radius is adopted in order to improve the efficiency of the new model. Finally, the new ionospheric tomography model is applied to reconstruct the temporal-spatial distribution of ionospheric electron density using actual global navigation satellite system observations. The comparisons between the tomographic profiles and those obtained from ionosonde data further demonstrate the reliability and the superiority of the new ionospheric tomography model.

3D ionospheric tomography over southeast China using a new scheme of constrained least squares with parameter weight matrix

Computerized ionospheric tomography (CIT) technique allows reconstructing the 3-dimensional and even 4-dimensional state of the ionosphere in terms of electron content. It is typically denoted as an inverse problem. Due to the unevenly distributed measurements collected from ground receivers, together with limited cut-off elevation angle and density of observations, the normal matrix of the inverse problem becomes sparse, which makes it difficult to solve and potentially unstable. The main methods to resolve the problem are the constraint algorithms during the process of the electron density inversion. A new approach is proposed for CIT that aims to mitigate the ill-posed problem and improve the precision of ionospheric electron density (IED) resolution by using constrained least squares. In this method, a new scheme of constructing parameter weight matrices in both horizontal and vertical directions is designed using the correlation of IED among those neighboring voxels, which helps to add the information needed in CIT, and the ill-posed problem is efficiently resolved. The accuracy and feasibility of the additional constraints in CIT algorithms are verified by numerical simulation tests. Finally, this newly developed method of constrained least squares is used to effectively perform the tomographic reconstruction of IED distribution over southeast China.

Comparative Investigation of the Digisonde‐Derived Electron Density and IRI Profiles at a Station Near the African Magnetic Equator During LowSolar Activity

AbstractInformation about the altitudinal distribution of ionospheric electron density (Ne) in the African equatorial region is very scanty. Such information is vital in modeling of key ionospheric parameters and in solving problems associated with ionospheric effect on the applications that depend on space‐based positioning and navigation systems. For the first time in the region, the entire Ne profile estimated by IRI is evaluated with profile deduced from Digisonde measurements at Ilorin (geog. 8.50°N, 4.50°E; dip. −7.9°) during low solar activity condition. Results obtained indicate marked discrepancies between the two profiles with magnitudes revealing local time and altitudinal and seasonal trends. The largest discrepancies are found between the post‐sunset and sunrise hours regardless of the seasons. IRI underestimates the Ne around the F1 region and is most obvious in September equinox and December solstice with deviation up to −42%. Further, IRI overestimates the Ne around the F2 peak at daytime, mostly in solstices with maximum deviation of 26% at F2 peak. All the three IRI topside models generally underestimate Ne between ~F2 peak height and ~800 km in all seasons except in June solstice. Beyond 800 km, the Digisonde, IRI‐NeQuick, and IRI‐2001 Corr profiles tend to converge, whereas the IRI‐2001 profiles still show marked divergence from the Digisonde profile in all seasons. Overall, the IRI‐NeQuick gives the best topside representation at this station at daytime/nighttime. The comparative analysis of Digisonde‐TEC and GPS‐TEC suggests that the Digisonde alone may not give a reliable estimate of the nighttime topside altitudinal distribution of Ne at this location.

A new parameterized approach for ionospheric tomography

GNSS ionospheric tomography technique is capable to reconstruct the high-quality 3D ionospheric electron density (IED) images with a relatively low cost. We present a new parameterized approach for refining the voxel-based ionospheric tomography modeling. This approach is different to most other voxel-based techniques as they assume a homogeneous IED distribution in each voxel that is unreasonable for the tomography modeling. In this method, IED of any point within a voxel is determined via vertically exponential interpolation and horizontally inverse distance weighted interpolation from the IED values at the eight corners of that voxel. The parameterized tomography is tested with real data collected over the period of June 1–30, 2015 from 45 GPS stations in south China. The superiority of the new parameterized method is verified by comparison with the traditional nonparametric method. The new parameterized method outperforms the traditional method by 12%, 10%, 5% and 2% for vertical resolutions of 25 km, 50 km, 75 km and 100 km, respectively, in the self-consistency validation by GPS data. Such improvements are 20%, 24%, 22% and 16%, respectively, when assessed by the Swarm in situ IEDs. In terms of the vertical layer discretization, configurations using the resolution of 25 km generally performs better than the other three vertical resolutions. Overall, the parameterized method using a vertical resolution of 25 km achieves the best performance from the comprehensive comparisons with ionospheric data derived by GPS, ionosonde and Swarm satellites.