Search

Besides the total electron content (TEC) and the peak electron density of the F2 layer (NmF2) also the equivalent slab thickness (τ) is considered as a key parameter of the ionosphere. The equivalent slab thicknessτ is a shape parameter of the vertical electron density profile defined by the ratio TEC/NmF2. Because the equivalent slab thickness is closely related to the neutral gas scale height under steady-state conditions at daytime, τ is correlated with thermospheric temperature changes. Based on earlier studies that have revealed a continuous decrease of the slab thickness over European ionosonde stations since 1996, the presentation analyses trends of τ in relation to TEC and NmF2 for an extended data base up to the end of 2022. A continuous TEC data base is available from DLR and the International GNSS Service (IGS). The peak electron density NmF2 is obtained from vertical sounding stations located over Europe covering solar cycles 23 and 24. In addition to the ground-based observations of NmF2 also radio occultation measurements are used to derive the peak electron density for equivalent slab thickness computations. It is found that the slab thickness behavior continues decreasing, e.g. by 50-60 km over the Juliusruh ionosonde station since 1996. The geophysical and solar activity background will be considered to explain the significant contraction of the equivalent slab thickness associated with a long-term cooling of the thermosphere. It is concluded that the slab thickness analysis has a high potential for estimating long-term changes in the thermosphere/ionosphere systems., The 28th IUGG General Assembly (IUGG2023) (Berlin 2023)

High-resolution 3D ionosphere modelling is crucial to monitor and understand space weather. Prol et al. (2021) is one of the most recent developments in advanced global-scale tomography that can represent the ionospheric electron density in a relatively high spatial and temporal resolution. The study was applied to the March 17, 2015 geomagnetic storm, showing that global-scale tomography could be useful to reproduce the system dynamics during a severe geomagnetic storm. Our results show good agreement with several ground- and space-based measurements of the ionosphere. The dataset used as reference is based on electron density observations from worldwide ionosondes, Millstone Hill incoherent scatter radar and in-situ measurements from the DMSP, GRACE and SWARM missions. The obtained results presented more accurate electron density profiles when compared to the Neustrelitz Electron Density Model, which is the model used as background, and physics-based thermosphere-ionosphere-electrodynamics general circulation model (TIE-GCM). Here we will present the main results in Prol et al. (2021), as well as provide guidelines and recommendations for future improvements in the development of 3D ionospheric models. Reference: Prol, F. S., Kodikara, T., Hoque, M. M., & Borries, C. (2021). Global-scale ionospheric tomography during the March 17, 2015 geomagnetic storm. Space Weather, 19, e2021SW002889. https://doi.org/10.1029/2021SW002889

GNSS single-frequency applications are affected by the interaction of the radio signals with the free electrons of the ionosphere, introducing range errors of up to 100 m in the L-band. Besides the GPS Klobuchar and the Galileo NeQuick ionospheric models, also the Neustrelitz Total Electron Content Model (NTCM) has been proposed as a practicable solution to mitigate such propagation errors.In this investigation we present a global statistical validation of the NTCM version driven by Galileo ionization coefficients (NTCM-GlAzpar) by comparing its performance with the performance achieved by the NeQuickG and Klobuchar models in the position domain. For this aim, we used the ESA analysis tool “gLAB” in the Standard Point Positioning (SPP) mode and GNSS data over two different time periods. The first period covers one month of perturbed solar activity (December 2014) and the second period corresponds to a month of quiet conditions (December 2019). We achieved a worldwide coverage with data from up to 73 receivers of the International GNSS Service (IGS).Our statistical analysis allows us to conclude that the NTCM-GlAzpar model clearly outperforms the results achieved with the Klobuchar model and is slightly better, or at least comparable, to the performance shown by NeQuickG. Indeed, the root mean squared (RMS) values of the hourly mean 3D position errors obtained for the global dataset are 4.36, 4.61 and 6.71 meters for perturbed conditions and 2.32, 2.35 and 2.75 meters for the quiet period, respectively for the NTCM-GlAzpar, NeQuickG and Klobuchar models. Nevertheless, through a geographic- and diurnal-specific analysis, we identify also that the performance of NTCM-GlAzpar slightly decreases for conditions of reduced solar activity – at night time, higher latitudes and low perturbations.We further discuss the applicability of the NTCM-GlAzpar ionospheric model in GNSS single-frequency applications motivated by its simple software adaptations and low computational cost.

Under certain conditions, the ionization of the E layer can dominate over that of the F2 layer. This phenomenon is called the E layer dominated ionosphere (ELDI) and occurs mainly in the auroral regions. In the present work, we model the variation of the ELDI for the Northern and Southern Hemispheres. Our proposed Neustrelitz ELDI Event Model (NEEM) is an empirical, climatological model that describes ELDI characteristics by means of four submodels for selected model observables, considering the dependencies on appropriate model drivers. The observables include the occurrence probability of ELDI events and typical E layer parameters that are important to describe the propagation medium for High Frequency (HF) radio waves. The model drivers are the geomagnetic latitude, local time, day of year, solar activity and the convection electric field. During our investigation, we found clear trends for the model observables depending on the drivers, which can be well represented by parametric functions. In this regard, the submodel NEEM-N characterizes the peak electron density NmE of the E layer, while the submodels NEEM-H and NEEM-W describe the corresponding peak height hmE and the vertical width wvE of the E layer electron density profile, respectively. Furthermore, the submodel NEEM-P specifies the ELDI occurrence probability %ELDI. The dataset underlying our studies contains more than two million vertical electron density profiles covering a period of almost 13 years. These profiles were derived from ionospheric GPS radio occultation observations on board the six COSMIC/FORMOSAT-3 satellites (Constellation Observing System for Meteorology, Ionosphere and Climate/Formosa Satellite Mission 3). We divided the dataset into a modeling dataset for determining the model coefficients and a test dataset for subsequent model validation. The normalized root mean square deviation (NRMS) between the original and the predicted model observables yields similar values across both datasets and both hemispheres. For NEEM-N, we obtain an NRMS varying between 36.1% and 47.1% and for NEEM-H, between 6.1% and 6.3%. In the case of NEEM-W, the NRMS varies between 38.5% and 41.1%, while it varies between 56.5% and 60.3% for NEEM-P. In summary, the proposed NEEM utilizes primary relationships with geophysical and solar wind observables, which are useful for describing ELDI occurrences and the associated changes of the E layer properties. In this manner, the NEEM paves the way for future prediction of the ELDI and of its characteristics in technical applications, especially from the fields of telecommunications and navigation.

High-precision GNSS (global navigation satellite e system) measurements can be used for remote sensing and nowadays play a significant role in atmospheric sounding (station data, radio occultation observations) and sea surface altimetry based on reflectometry. A limiting factor of high-precision reflectometry is the loss of coherent phase information due to sea-state-induced surface roughness. This work studies airborne reflectometry observations recorded over coastal waters to examine the sea-state influence on Doppler distribution and the coherent residual phase retrieval. From coherent observations, the possibility of zenith total delay inversion is also investigated, considering the hydrostatic mapping factor from the Vienna mapping function and an exponential vertical decay factor depending on height receiver changes. The experiment consists of multiple flights performed along the coast between the cities of Calais and Boulogne-sur-Mer, France, in July 2019. Reflected signals acquired in a right-handed circular polarization are processed through a model-aided software receiver and passed through a retracking module to obtain the Doppler and phase-corrected signal. Results from grazing angle observations (elevation < 15°) show a high sensitivity of Doppler spread with respect to sea state with correlations of 0.75 and 0.88 with significant wave height and wind speed, respectively. An empirical Doppler spread threshold of 0.5 Hz is established for coherent reflections supported by the residual phase observations obtained. Phase coherence occurs in 15% of the observations; however, the estimated zenith total delay for the best event corresponds to 2.44 m, which differs from the typical zenith total delay (2.3 m) of 5%.

The correct representation of global-scale electron density is crucial for monitoring and exploring the space weather. This study investigates whether the ground-based Global Navigation Satellite System (GNSS) tomography can be used to reflect the global spatial and temporal responses of the ionosphere under storm conditions. A global tomography of the ionosphere electron density is constructed based on data from over 2,700 GNSS stations. In comparison to previous techniques, advances are made in spatial and temporal resolution, and in the assessment of results. To demonstrate the capabilities of the approach, the developed method is applied to the March 17, 2015 geomagnetic storm. The tomographic reconstructions show good agreement with electron density observations from worldwide ionosondes, Millstone Hill incoherent scatter radar and in-situ measurements from satellite missions. Also, the results show that the tomographic technique is capable of reproducing plasma variabilities during geomagnetically disturbed periods including features such as equatorial ionization anomaly enhancements and depletion. Validation results of this brief study period show that the accuracy of our tomography is better than the Neustrelitz Electron Density Model, which is the model used as background, and physics-based thermosphere-ionosphere-electrodynamics general circulation model. The results show that our tomography approach allows us to specify the global electron density from ground to ∼900 km accurately. Given the demonstrated quality, this global electron density reconstruction has potential for improving applications such as assessment of the effects of the electron density on radio signals, GNSS positioning, computation of ray tracing for radio-signal transmission, and space weather monitoring.

At some locations, especially in the auroral regions, the ionization of the E layer can dominate over that of the F2 layer, which is called the E layer dominated ionosphere (ELDI). In the present work we investigate the spatiotemporal variation of the ELDI depending on the season, solar activity, geomagnetic activity, interplanetary magnetic field, convection electric field, and solar wind energy. We specify each distribution of ELDI events by the values of four parameters. In this regard, we compute the height, full width at half maximum, and position of a Gaussian function relative to a precomputed reference ellipse as parameters to describe the spatial distribution of ELDI events in geocentric latitude/longitude coordinates. To study the temporal variation of the ELDI events, we estimate the weighted mean local time of the distribution as the fourth parameter. The database used for our investigations contains more than 3.5 million vertical electron density profiles derived from ionospheric GPS radio occultation observations on board the COSMIC/FORMOSAT-3 (Constellation Observing System for Meteorology, Ionosphere, and Climate/Formosa Satellite Mission 3) mission, covering a period of almost 13 years. The analysis of observations representing changing geophysical conditions results in clear trends for all ELDI parameters. In this context, the mean local time varies mostly between 01:00 and 02:00 local time, while the probability of ELDI occurrence is increased in local winter and in the case of low solar activity. Likewise, an increase in the solar wind parameters increases the number of ELDI events and leads to an equatorward shift of their position. The relationships found in our investigations can serve as a basis for future modeling studies addressing ELDI occurrences as a function of selected geophysical quantities.

The next generation of satellite-based augmentation systems (SBAS) will support aviation receivers that take advantage of the ionosphere-free dual-frequency combination. By combining signals of the L1 and L5 bands, about 99% of the ionospheric refraction effects on the GNSS (Global Navigation Satellite Systems) signals can be removed in the user receivers without additional SBAS corrections. Nevertheless, even if most of the negative impacts on GNSS signals are removed by the ionospheric-free combination, some residuals remain and have to be taken into account by overbounding models in the integrity computation conducted by safety-of-live (SoL) receivers in airplanes. Such models have to overbound residuals as well, which result from the most rare extreme ionospheric events, e.g., such as the famous “Halloween Storm”, and should thus include the tails of the error distribution. Their application shall lead to safe error bounds on the user position and allow the computation of protection levels for the horizontal and vertical position errors. Here, we propose and justify such an overbounding model for residual ionospheric delays that remain after the application of the ionospheric-free linear combination. The model takes into account second- and third-order ionospheric refraction effects, excess path due to ray bending, and increased ionospheric total electron content (TEC) along the signal path due to ray bending.

Ionospheric delay is a significant error source in multi-GNSS positioning. We present different processing strategies to fully exploit the ionospheric delay effects on multi-frequency and multi-GNSS positioning performance, including standard point positioning (SPP) and precise point positioning (PPP) scenarios. Datasets collected from 10 stations over thirty consecutive days provided by multi-GNSS experiment (MGEX) stations were used for single-frequency SPP/PPP and dual-frequency PPP tests with quad-constellation signals. The experimental results show that for single-frequency SPP, the Global Ionosphere Maps (GIMs) correction achieves the best accuracy, and the accuracy of the Neustrelitz TEC model (NTCM) solution is better than that of the broadcast ionospheric model (BIM) in the E and U components. Eliminating ionospheric parameters by observation combination is equivalent to estimating the parameters in PPP. Compared with the single-frequency uncombined (UC) approach, the average convergence time of PPP with the external ionospheric models is reduced. The improvement in BIM-, NTCM- and GIM-constrained quad-constellation L2 single-frequency PPP was 15.2%, 24.8% and 28.6%, respectively. The improvement in convergence time of dual-frequency PPP with ionospheric models was different for different constellations and the GLONASS-only solution showed the least improvement. The improvement in the convergence time of BIM-, NTCM- and GIM-constrained quad-constellation L1/L2 dual-frequency PPP was 5.2%, 6.2% and 8.5%, respectively, compared with the UC solution. The positioning accuracy of PPP is slightly better with the ionosphere constraint and the performance of the GIM-constrained PPP is the best. The combination of multi-GNSS can effectively improve the positioning performance

To facilitate F2-layer peak density (NmF2) modeling, a nonlinear polynomial model approach based on global NmF2 observational data from ionospheric radio occultation (IRO) measurements onboard the CHAMP, GRACE, and COSMIC satellites, is presented in this paper. We divided the globe into 63 slices from 80°S to 80°N according to geomagnetic latitude. A Nonlinear Polynomial Peak Density Model (NPPDM) was constructed by a multivariable least squares fitting to NmF2 measurements in each latitude slice and the dependencies of NmF2 on solar activity, geographical longitude, universal time, and day of year were described. The model was designed for quiet and moderate geomagnetic conditions (Ap ≤ 32). Using independent radio occultation data, quantitative analysis was made. The correlation coefficients between NPPDM predictions and IRO data were 0.91 in 2002 and 0.82 in 2005. The results show that NPPDM performs better than IRI2016 and Neustrelitz Peak Density Model (NPDM) under low solar activity, while it undergoes performance degradation under high solar activity. Using data from twelve ionosonde stations, the accuracy of NPPDM was found to be better than that of NPDM and comparable to that of IRI2016. Additionally, NPPDM can well simulate the variations and distributions of NmF2 and describe some ionospheric features, including the equatorial ionization anomaly, the mid-latitude trough, and the wavenumber-four longitudinal structure.