Ground-based Global Navigation Satellite System Research Articles

The Zenith Total Delay Combination of International GNSS Service Repro3 and the Analysis of Its Precision

Currently, ground-based global navigation satellite system (GNSS) techniques have become widely recognized as a reliable and effective tool for atmospheric monitoring, enabling the retrieval of zenith total delay (ZTD) and precipitable water vapor (PWV) for meteorological and climate research. The International GNSS Service analysis centers (ACs) have initiated their third reprocessing campaign, known as IGS Repro3. In this campaign, six ACs conducted a homogeneous reprocessing of the ZTD time series spanning the period from 1994 to 2022. This paper primarily focuses on ZTD products. First, the data processing strategies and station conditions of six ACs were compared and analyzed. Then, formal errors within the data were examined, followed by the implementation of quality control processes. Second, a combination method is proposed and applied to generate the final ZTD products. The resulting combined series was compared with the time series submitted by the six ACs, revealing a mean bias of 0.03 mm and a mean root mean square value of 3.02 mm. Finally, the time series submitted by the six ACs and the combined series were compared with VLBI data, radiosonde data, and ERA5 data. In comparison, the combined solution performs better than most individual analysis centers, demonstrating higher quality. Therefore, the advanced method proposed in this study and the generated high-quality dataset have considerable implications for further advancing GNSS atmospheric sensing and offer valuable insights for climate modeling and prediction.

Evaluation of GOCE/GRACE and combined global geopotential models using GNSS/levelling data over Nigeria

AbstractGlobal Geopotential Models (GGMs) provide valuable information about Earth’s gravity field functionals, such as geoid heights and gravity anomalies. However, ground-based datasets are required to validate these GGMs at the regional and local scales. In this study, we validated the accuracy of GGMs by comparing them with ground-based Global Navigation Satellite System (GNSS)/levelling data for the first time in Nigeria. We employed two validation scenarios: with and without considering spectral consistency using the spectral enhancement method (SEM) to incorporate high and very high frequencies of the gravity field spectrum from the combined global gravity field model (XGM2019e_2159) and the residual terrain model (RTM) derived from the Shuttle Radar Topography Mission (SRTM) data, respectively. The results of this evaluation confirmed that the application of SEM improved the assessment of the GGM solutions in an unbiased manner. Integrating XGM2019e_2159 and SRTM data to constrain the high-frequency component of geoid heights in Gravity Field and Steady-State Ocean Circulation Explorer (GOCE)-based GGMs led to an improvement of approximately 10% in reducing the standard deviation (STD) relative to when SEM was not applied. GO_CONS_GCF_2_TIM_R6 at spherical harmonics (SH) of up to degree and order 260 demonstrated the lowest STD when compared to GO_CONS_GCF_2_DIR_R6 and GO_CONS_GCF_2_SPW_R5, with a reduction from 0.380 m without SEM application to 0.342 m with SEM implementation. In addition, four transformation models, namely, linear, four-parameter, five-parameter, and seven-parameter models, were evaluated. The objective is to mitigate the reference system offsets between the GNSS/levelling data and the GGMs and to identify the particular parametric model with the smallest STD across all GGMs. This effort reduced the GGMs misfits to GNSS/levelling to 0.30 m, representing a 15.3% decrease in STD. Notably, the XGM2019e_2159 model provides this improvement.

An Observed Relationship Between Satellite‐Estimated Transmittance and Ground‐Estimated Water Vapor: Implications for High‐Temporal‐Resolution Water Vapor Retrieval From Non‐Geostationary Satellite Measurements

AbstractThe magnitude of integrated water vapor (IWV) varies considerably in the spatial‐temporal domain, which demonstrates the significance of high‐spatiotemporal‐resolution IWV observations for atmospheric water vapor distribution monitoring, both locally and globally. Unlike previously published algorithms based on data fusion, an empirical retrieval algorithm is for the first time proposed to directly retrieve high‐temporal‐resolution IWV estimates from near‐infrared radiance observations of the non‐geostationary Ocean and Land Color Instrument (OLCI). The retrieval algorithm is developed based on an observed regression relationship between satellite‐based OLCI‐estimated transmittance and ground‐based Global Navigation Satellite System (GNSS)‐estimated IWV in the temporal domain. The results show that all newly retrieved IWV estimates have an overall good consistency with ground‐based IWV from additional GNSS‐sensed measurements, indicating the feasibility of the retrieval approach. The performance of the retrieval algorithm is acceptable and satisfactory when compared with that of IWV retrievals listed in previous studies.

Observations of ionospheric irregularities during a geomagnetic storm based on the C-band sentinel-1

ABSTRACT Disturbances in the ionosphere can severely affect synthetic aperture radar (SAR) acquisition and also the accuracy of differential SAR interferometry (D-InSAR) techniques. Focusing on the mid-latitudes, we report for the first time the effect of small-scale ionospheric irregularities on the interferograms generated by the C-band Sentinel-1 satellite. The interferograms obtained during the geomagnetic storm on 4 November 2021 showed finger-like stripes of phase artefacts, while the interferograms obtained on non-storm days did not show such anomalies. Such phase artefacts introduce an error greater than ± 3π in the interferogram. By further checking the gridded total electron content (TEC) from the ground-based global navigation satellite system (GNSS) network, as well as the in-situ electron density measured by ESA’s Swarm satellite, the finger-like stripes in the interferograms are found to be correlated with ionospheric bulge structures generated during the geomagnetic storm. Since both the interferogram stripes and the ionospheric bulges exhibit a comparable east-west extension, it is believed that the interferograms can also be used to reflect the small-scale structures of ionospheric bulge fringes. For example, in this case, the latitudinal width of each finger-like stripe reaches about 8 km. Our result implies that the C-band interferogram has the potential to probe the two-dimensional ionospheric structures with unprecedented resolution, which cannot be achievable by the gridded GNSS-TEC.

Influence of the October 14, 2023, Ring of Fire Annular Eclipse on the Ionosphere: A Comparison Between GNSS Observations and SAMI3 Model Prediction

AbstractCelestial phenomena such as solar eclipses disrupt the ionosphere's inherent photochemical, dynamic, and electrodynamic processes and can be viewed as a natural experiment that provides a unique opportunity to study ionospheric perturbations. We investigate the spatiotemporal ionospheric response induced by the 14 October 2023 annular eclipse using ground‐based Global Navigation Satellite System (GNSS) receivers located over the continental region of North and South America. The largest total electron content (TEC) change (∼22 TECU decrease) is observed over the path of the eclipse at around 16°–18°N which lies just before the greatest eclipse location. However, the percentage change in TEC here is ∼44% (of the background value), which is less than that observed at midlatitudes around 30°–35°N (∼18.7 TECU, ∼50%). We observed a latitudinal dependency of TEC variation and time delay in ionospheric response to the eclipse with midlatitudes experiencing greater TEC changes and longer time delays compared to low‐latitude and equatorial regions. The SAMI3 model used to simulate the impact of the eclipse, captures the large decrease in VTEC along the eclipse path. Interestingly, increases in the TEC are also observed in several GNSS sites and it varied from 4.5 to 6.5 TECU in the Northern Hemisphere and 7 to 12 TECU in the Southern Hemisphere. The SAMI3 model results also show an enhancement in TEC but significantly less than that observed. Moreover, the SAMI3 model simulations for individual GNSS site coordinates within the conjugate regions failed to predict accurately, the TEC enhancement recorded by the GNSS sites.

Differences among the total electron content derived by radio occultation, global ionospheric maps and satellite altimetry

In recent years, significant progress has been in ionospheric modeling research through data ingestion and data assimilation from a variety of sources, including ground-based global navigation satellite systems, space-based radio occultation and satellite altimetry (SA). Given the diverse observing geometries, vertical data coverages and intermission biases among different measurements, it is imperative to evaluate their absolute accuracies and estimate systematic biases to determine reasonable weights and error covariances when constructing ionospheric models. This study specifically investigates the disparities among the vertical total electron content (VTEC) derived from SA data of the Jason and Sentinel missions, the integrated VTEC from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) and global ionospheric maps (GIMs). To mitigate the systematic bias resulting from differences in satellite altitudes, the vertical ranges of various VTECs are mapped to a standardized height. The results indicate that the intermission bias of SA-derived VTEC remains relatively stable, with Jason-1 serving as a benchmark for mapping other datasets. The mean bias between COSMIC and SA-derived VTEC is minimal, suggesting good agreement between these two space-based techniques. However, COSMIC and GIM VTEC exhibit remarkable seasonal discrepancies, influenced by the solar activity variations. Moreover, GIMs demonstrate noticeable hemispheric asymmetry and a degradation in accuracy ranging from 0.7 to 1.7 TECU in the ocean-dominant Southern Hemisphere. While space-based observations effectively illustrate phenomena such as the Weddell Sea anomaly and longitudinal ionospheric characteristics, GIMs tend to exhibit a more pronounced mid-latitude electron density enhancement structure.

Performance of Ground-Based Global Navigation Satellite System Precipitable Water Vapor Retrieval in Beijing with the BeiDou B2b Service

Performance of Ground-Based Global Navigation Satellite System Precipitable Water Vapor Retrieval in Beijing with the BeiDou B2b Service

Assessing GNSS hydrological monitoring capability across different climatic settings in China

Vertical position changes of ground-based Global Navigation Satellite System (GNSS) stations have been used to study regional terrestrial water storage (TWS) changes. However, the feasibility is still unclear in many areas due to diverse local effects. This study aims to evaluate the capability of GNSS vertical displacements in monitoring hydrological variations in four climate settings over Chinese mainland. The spatial and temporal variations of hydrological load-induced (HYDL) vertical displacements at 208 GNSS sites during 2011–2020 were analyzed by comparing with Gravity Recovery and Climate Experiment (GRACE)/GRACE Follow-On (GFO) and Global Land Data Assimilation System (GLDAS) derived TWS changes. The results indicate that GNSS vertical positions show different capabilities in capturing seasonal and non-seasonal hydrological dynamics in different climate regions. Among the four climatic settings, the subtropical monsoon climate (SMC) region, with the largest deformation fluctuation (the regional mean root mean square (RMS) is 7.97 mm), has the highest regional mean HYDL-GRACE and HYDL-GLDAS anti-correlation coefficients (CCs) of −0.47 and −0.45 at the seasonal scale, respectively. For the individual GNSS site, the number of the sites with CC <−0.40 between HYDL and GRACE/GLDAS-derived TWS changes accounts for 55.1% and 55.1% (SMC), 13.0% and 7.4% (temperate monsoon climate, TMC), 6.7% and 13.3% (temperate continental climate, TCC), 32.3% and 38.7% (plateau climate, PC), respectively. For the non-seasonal term, although the proportion with CC <−0.40 in each climate type decreases mainly due to the influence of local geodynamic and human activities, especially in the SMC and PC regions, GNSS site vertical deformations still show good capability in monitoring hydrological extremes. The results provide valuable information for better application of GNSS to hydrology.

CORS station for synergistic monitoring of multivariate surface parameters in expansive soils

Expansive soils cause frequent surface deformation due to their expansion and contraction, which is a serious engineering hazard, and long-term subsidence monitoring is a prerequisite for preventing and controlling expansive soil disasters. Currently, the conventional monitoring methods for the above issue include Interferometric Synthetic Aperture Radar (InSAR) technology, but InSAR is not suitable for uninterrupted monitoring of surface deformation and has low sensitivity. Meanwhile, it can’t obtain multiple surface environmental parameters around the station. The Global Navigation Satellite System (GNSS), a system that can directly acquire surface deformation, has been widely used in landslide disaster monitoring, and in recent years, this technology has also been applied to the field of expansive soil disaster monitoring. At the same time, GNSS can also provide a constant stream of L-band microwave signals to obtain ground environmental information such as precipitable rainfall and soil moisture around the station. In previous studies of expansive soil hazards, GNSS technology has been mainly used to provide surface deformation information without exploring its potential to invert ground environmental information around stations. This paper proposes a ground-based GNSS remote sensing integrated monitoring system that integrates expanding land surface parameters such as “precipitable rainfall, soil moisture, and three-dimensional deformation” and analyses the ability of ground-based GNSS to be used for integrated monitoring of expanding soil hazards by combining ten years of consecutive observational data from GNSS stations along the coastal area of Houston. The experimental results show that the GNSS is capable of providing highly accurate time-series characterization of deformation, and inelastic subsidence in recent years has resulted in a cumulative permanent elevation loss of 2 cm along the Houston coast. The correlation coefficient between soil moisture extracted by the fifth-generation European reanalysis data (ERA5) and soil moisture inverted by ground-based GNSS is 0.514. At the same time, the GNSS was also able to monitor the zenithal precipitable water vapor (PWV) and soil moisture changes around the GNSS station and further analyze the response relationship among the three parameters, which could comprehensively evaluate the stability of expansive soils, avoiding the unreliability of relying on a single piece of monitoring information to assess the stability of expansive soils. We hope to construct a more comprehensive ground-based GNSS remote sensing monitoring system to better monitor expansive soil hazards.

High-resolution 3-D imaging of electron density perturbations using ultra-dense GNSS observation networks in Japan: an example of medium-scale traveling ionospheric disturbances

For the first time using computerized ionospheric tomography (CIT) and leveraging ultra-dense slant total electron content (STEC) measurements derived from two ground-based Global Navigation Satellite System (GNSS) receiver networks in Japan, we have reconstructed the 3-D field-aligned structure of nighttime medium-scale traveling ionospheric disturbances (MSTIDs) with high spatiotemporal resolution. The CIT algorithm focuses on electron density perturbation components, allowing for the imaging of disturbances with small amplitudes and scales. Slant TECs used for CIT are setup to consist of two components: the background derived from IRI-2016 model and TEC perturbations obtained by subtracting a 30-min running average from observations. The resolution is set to 0.25º in latitude and longitude, 10 km in altitude, 30 s in time. Simulations were conducted to assess the performance of the CIT algorithm, revealing that this technique has good fidelity by accurately reconstructing more than 80% of the electron density perturbations. The focus is on the nighttime event of July 4, 2022, when data were accessible. The reconstruction results show that the MSTIDs initially form at lower altitudes and subsequently develop to exhibit large amplitudes and scales that extend to higher altitudes, characterized by a well-defined frontal structure with electrodynamic signatures. These results are consistent with theories and snippets of observational evidence regarding electromagnetic-influenced MSTIDs, hence affirming the effectiveness of the developed CIT technique in probing of the variations in the 3-D structure of ionospheric electron density. This is expected to contribute to a compressive understanding of the underlying mechanisms of ionospheric inhomogeneities.Graphical

Comparative analysis of precipitable water vapour data in the Tarim Basin, China

Atmospheric water vapour is an important contributor to global energy and water cycles and its detection is necessary for precipitation forecasting and climate change research. To explore the applicability of reanalysis data for determining precipitable water vapour (PWV) content in the Tarim Basin, China, this study compared the reliability of two reanalysis PWV products, ERA-5 (PWERA) and MERRA-2 (PWMER), using data from four ground-based global navigation satellite system stations representing basin ecosystems (Oasis, Gobi Desert, Central Desert, and Alpine; dataset: PWGNSS) as reference values. We conducted correlation analysis between PWGNSS and both reanalysis PWV products at different time scales to verify the products’ accuracy. The results showed that the reanalysis PWV products were mostly overestimated at the Desert station, while they were mostly underestimated at the Oasis station. The applicability of PWERA and PWMER varied by season and location, with better applicability at the Oasis station from spring to autumn. The applicability of reanalysis PWV products was lower during precipitation periods than during non-precipitation, and varied by location. During non-precipitation situation, the PWMER had better applicability at the Oasis station and at the Central Desert station when the PWGNSS exceeded 25 mm. Meanwhile, PWERA had better applicability in precipitation situation at the Gobi Desert station from April to June and at Oasis station from August to September when PWGNSS was above 30 mm, at the Central Desert station from May to August, and when the PWGNSS exceeded 30 mm.

A Statistical Analysis of the Morphology of Storm‐Enhanced Density Plumes Over the North American Sector

AbstractThe storm‐enhanced density (SED) is a large‐scale midlatitude ionospheric electron density enhancement in the local afternoon sector, which exhibits substantial spatial gradients and thus can impose detrimental effects on modern navigation and communication systems, causing potential space weather hazards. This study has identified a comprehensive list of 49 SED events over the continental US and adjacent regions, by examining strong geomagnetic storms occurring between 2000 and 2023. The ground‐based Global Navigation Satellite System (GNSS) total electron content and data from a new TEC‐based ionospheric data assimilation system were used to analyze the characteristics of SED. For each derived SED events, we have quantified its morphology by employing a Gaussian function to parameterize key characteristics of the SED, such as the plume intensity, central longitude, and half‐width. A statistical analysis of SEDs was conducted for the first time to characterize their climatological features. We found that the SED distribution exhibits a higher peak intensity and a narrower width as geomagnetic activity strengthens. The peak intensity of SED has maximum values around the equinoxes in their seasonal distribution. Additionally, we observed a solar cycle dependence in the SED distribution, with more events occurring during the solar maximum and declining phases compared to the solar minimum. SED plumes exhibit a sub‐corotation feature with respect to the Earth, characterized by a westward drift speed between 50 and 400 m/s and a duration of 3–10 hr. These information advanced the current understanding of the spatial‐temporal variation of SED characteristics.

Crustal deformation from GNSS measurement and earthquake mechanism along Pieniny Klippen Belt, Southern Poland

The dynamic geological features of the Pieniny Klippen Belt (PKB) in southern Poland are nowadays a focal point of researchers as it is recognised as an active zone of crustal discontinuity. In the present study, we employed long-term analysis of ground-based global navigation satellite system (GNSS) measurements (from 2004 to 2020) to probe the crustal deformation, strain rates, and rotational rates in the PKB unit and the surrounding region, i.e., Magura Nappe (MN) and Podhale Flysch (PF) units. Measured velocities from GNSS observables are modelled by the auto-regressive integrated moving average (ARIMA) method to comprehend the long-term tectonic deformation. Our results showed that the ARIMA-modelled velocity varied from ~ 0.15 to ~ 8.86 mm/yr, indicating about 8.71 mm/yr difference along all units. Such differences suggest that crustal slip along the active thrusts and folds is the major factor causing regional deformation. The strain rates in PKB are also varying from the western to the eastern part. The rotational rates in PKB show a counterclockwise (CCW) pattern similar to the strain rates. These patterns suggesting that the PKB was rotated in the CCW direction with a large angle during the Miocene period. Finally, we analysed the seismicity for a period from 2004 to 2020 by using Bayesian moment tensor inversion and multivariate Bayesian inversion. The Bayesian inversion was applied based on bootstrapping chain analysis to figure out the earthquake mechanism using moment tensor inversion for the mainshock that occurred in Poland on 20 July 2018. The inversion results for the 2018 earthquake resolved a thrusting mechanism with nodal plane-1 having a strike of 346°, dip of 32°, and rake of 92° and a nodal plane-2 with a strike of 163°, dip of 58°, and rake of 89°. Since the seismicity in the Poland region has experienced less significant earthquakes in the last century, it is reasonable to attribute this lower seismic activity to the correspondingly low slip rates discerned through geodetic monitoring efforts.

A Case Study of Ionospheric Storm‐Time Altitudinal Differences at Low Latitudes During the May 2021 Geomagnetic Storm

AbstractPrevious studies paid little attention to the ionospheric storm‐time altitudinal differences due to insufficiency of ionospheric measurements. In this work, multiple instrumental observations were used to investigate the ionospheric storm‐time response at low latitudes in the American and Asian‐Australian sectors during the May 2021 geomagnetic storm. The ground‐based Global Navigation Satellite Systems (GNSS) total electron content (TEC) and Low Earth Orbit (LEO) satellite topside TEC presented opposite (positive/negative) variations in the low‐latitude and equatorial region of both sectors during this storm. The electron density profiles from the Constellation Observing System for Meteorology, Ionosphere, and Climate‐2 (COSMIC‐2) and the Sanya Incoherent Scatter Radar showed a good agreement and well explained the opposite variations between GNSS TEC and LEO satellite topside TEC. The F2‐layer peak height (hmF2) and peak density (NmF2) displayed inverse variations, and the feature was present mostly in the regions between equatorial ionization anomaly crests. The combined modulation effects of the storm‐time zonal electric fields and the field‐aligned transports possibly resulted in the contrary variations of hmF2 and NmF2 in the low‐latitude and equatorial region, leading to the storm‐time altitudinal differences during this storm. Relatively, the storm‐time thermospheric composition disturbances might be a minor factor responsible for these differences.

GNSS Differential Code Bias Determination Using Rao‐Blackwellized Particle Filtering

AbstractThe Assimilative Canadian High Arctic Ionospheric Model (A‐CHAIM) is a near‐real‐time data assimilation model of the high latitude ionosphere, incorporating measurements from many instruments, including slant Total Electron Content measurements from ground‐based Global Navigation Satellite System (GNSS) receivers. These measurements have receiver‐specific Differential Code Biases (DCB) which must be resolved to produce an absolute measurement, which are resolved simultaneously with the ionospheric state using Rao‐Blackwellized particle filtering. These DCBs are compared to published values and to DCBs determined using eight different Global Ionospheric Maps (GIM), which show small but consistent systematic differences. The potential cause of these systematic biases is investigated using multiple experimental A‐CHAIM test runs, including the effect of plasmaspheric electron content. By running tests using the GIM‐derived DCBs, it is shown that using A‐CHAIM DCBs produces the lowest overall error, and that using GIM DCBs causes an overestimation of the topside electron density which can exceed 100% when compared to in situ measurements from DMSP.

2‐D Total Electron Content and 3‐D Ionospheric Electron Density Variations During the 14 October 2023 Annular Solar Eclipse

AbstractThis study investigates the ionospheric total electron content (TEC) responses in the 2‐D spatial domain and electron density variations in the 3‐D spatial domain during the annular solar eclipse on 14 October 2023, using ground‐based Global Navigation Satellite System (GNSS) observations, a novel TEC‐based ionospheric data assimilation system (TIDAS), ionosonde measurements, and satellite in situ data. The main results are summarized as follows: (a) The 2‐D TEC responses exhibited distinct latitudinal differences. The mid‐latitude ionosphere exhibited a more substantial TEC decrease of 25%–40% along with an extended recovery time of 3–4 hr. In contrast, the equatorial and low‐latitude ionosphere experienced a smaller TEC reduction of 10%–25% and a faster recovery time of 20–50 min. The minimal eclipse effect was observed near the northern equatorial ionization anomaly crest region. (b) The ionospheric electron density variations during the eclipse were effectively reconstructed by TIDAS data assimilation in the 3‐D domain, providing important altitude information with validity. (c) The ionospheric electron density variations showed a notable altitude‐dependent feature. The eclipse led to a substantial electron density reduction of 30%–50%, with the maximum depletion occurring around the ionospheric F2‐layer peak height (hmF2) of 250–350 km. The post‐eclipse recovery of electron density exhibited a relatively slower pace near the F2‐layer peak height than that at lower and higher altitudes.

3‐D Ionospheric Imaging Over the South American Region With a New TEC‐Based Ionospheric Data Assimilation System (TIDAS‐SA)

AbstractThis study has developed a new TEC‐based ionospheric data assimilation system for 3‐D regional ionospheric imaging over the South American sector (TIDAS‐SA) (45°S–15°N, 35°–85°W, and 100–800 km). The TIDAS‐SA data assimilation system utilizes a hybrid Ensemble‐Variational approach to incorporate a diverse set of ionospheric data sources, including dense ground‐based Global Navigation Satellite System (GNSS) line‐of‐sight Total Electron Content (TEC) data, radio occultation data from the Constellation Observing System for Meteorology, Ionosphere, and Climate‐2 (COSMIC‐2), and altimeter TEC data from the JASON‐3 satellite. TIDAS‐SA can produce a reanalyzed three‐dimensional (3‐D) electron density spatial variation with a high time cadence, yielding spatial‐temporal resolution of 1° (latitude) × 1° (longitude) × 20 km (altitude) × 5 min. This allows us to reconstruct and study the 3‐D ionospheric morphology with multi‐scale structures. The performance of the data assimilation system is validated against independent ionosonde and in situ measurements through an experiment for a strong geomagnetic storm event on 03–04 November 2021. The results demonstrate that TIDAS‐SA can provide detailed and altitude‐resolved information that accurately characterizes the storm‐time ionospheric disturbances in vertical and horizontal domains over the equatorial and low‐latitude regions of South America.

Ionospheric response to the 2020 Samos earthquake and tsunami

On 30 October 2020 at 11:51 UT, a magnitude 7.0 earthquake occurred in the Dodecanese sea (37.84°N, 26.81°E, 10 km depth) and generated a tsunami with an observed run-up of more than 1 m on the Turkish coasts. Both the earthquake and the tsunami produced acoustic and gravity waves that propagated upward, triggering co-seismic and co-tsunamic ionospheric disturbances. This paper presents a multi-instrumental study of the ionospheric impact of the earthquake and related tsunami based on ionosonde data, ground-based Global Navigation Satellite Systems (GNSS) data and data from DORIS beacons received by Jason3 in the Mediterranean region. Our study focuses on the Total Electron Content to describe the propagation of co-seismic and co-tsunami ionospheric disturbances (CSID, CTID), possibly related to gravity waves triggered by the earthquake and tsunami. We use simultaneous vertical ionosonde soundings to study the interactions between the upper and lower atmosphere, highlighting the detection of acoustic waves generated by the seismic Rayleigh waves reaching the ionosonde locations and propagating vertically up to the ionosphere. The results of this study provide a detailed picture of the Lithosphere-Atmosphere–Ionosphere coupling in the scarcely investigated Mediterranean region and for a relatively weak earthquake.Graphical abstract

Performance of the double-thin-shell approach for studying nighttime medium-scale traveling ionospheric disturbances using two dense GNSS observation networks in Japan

Electrodynamic coupling between the ionospheric E and F regions is widely recognized as the underlying mechanism for generating medium-scale traveling ionospheric disturbances (MSTIDs) during nighttime at midlatitudes. Recently, the double-thin-shell approach has proven to be a useful tool for studying the E–F coupling. By using total electron content (TEC) measurements, this approach enables the simultaneous reconstruction of electron density perturbations in both the E and F regions with broad and continuous coverage. However, the current reconstruction performance is limited when using only GPS-TEC measurements from GEONET, a dense network of ground-based Global Navigation Satellite System (GNSS) receivers over Japan. The expansion of available data sources and the integration of multi-GNSS observation data are considered important to enhance the double-thin-shell model. Fortunately, SoftBank Corp., a Japanese telecommunications provider, has recently developed a dense independent GNSS observation network to improve positioning services. In this paper, we analyze the potential of the improved double-thin-shell approach and emphasize the importance of incorporating multi-GNSS observation data from both GEONET and SoftBank networks. The solvability analysis, simulation, and observation results collectively indicate a substantial improvement in the spatiotemporal resolution. Specifically, the longitudinal and latitudinal resolution is improved from 0.15° to 0.1° in the E region, and from 0.5° to 0.3° in the F region. The temporal resolution is also improved from 2 to 1 min. In addition, significant improvements have been achieved in the reconstruction performance, particularly for the E region under complex background conditions. Based on these assessments, we conclude that the incorporation of GEONET and SoftBank GNSS observation data holds significant potential for improving the double-thin-shell model and advancing our understanding of MSTIDs.Graphical abstract

Statistical models of the variability of plasma in the topside ionosphere: 2. Performance assessment

Statistical models of the variability of plasma in the topside ionosphere based on the Swarm data have been developed in the “Swarm Variability of Ionospheric Plasma” (Swarm-VIP) project within the European Space Agency’s Swarm+4D-Ionosphere framework. The models can predict the electron density, its gradients for three horizontal spatial scales – 20, 50 and 100 km – along the North-South direction and the level of the density fluctuations. Despite being developed by leveraging on Swarm data, the models provide predictions that are independent of these data, having a global coverage, fed by various parameters and proxies of the helio-geophysical conditions. Those features make the Swarm-VIP models useful for various purposes, which include the possible support for already available ionospheric models and proxy of the effect of ionospheric irregularities of the medium scales that affect the signals emitted by Global Navigation Satellite Systems (GNSS). The formulation, optimisation and validation of the Swarm-VIP models are reported in Paper 1 (Wood et al. 2024. J Space Weather Space Clim. in press). This paper describes the performance assessment of the models, by addressing their capability to reproduce the known climatological variability of the modelled quantities, and the ionospheric weather as depicted by ground-based GNSS, as a proxy for the ionospheric effect on GNSS signals. Additionally, we demonstrate that, under certain conditions, the model can better reproduce the ionospheric variability than a physics-based model, namely the Thermosphere-Ionosphere-Electrodynamics General Circulation Model (TIE-GCM).