Radio Occultation Data Research Articles

Assessment of satellite images terrestrial surface temperature and WVP using GNSS radio occultation data

Abstract Land Surface Temperature (LST) and Water Vapor Pressure (WVP) contour maps can be produced using cameras aboard satellites, for instance, under the name “Remote Sensing (RS)”. Satellite image observations should be verified before using based on a reliable data. Global Navigation Satellite System Radio Occultation (GNSS-RO) method is observing accurate Earth atmosphere parameters continuously. In the present research, LST and WVP differences between Landsat 8 (LC08), Sentinel-3 (S3), and MODIS (Terra and Aqua) images and GNSS-RO are assessed in Egypt depending on the satellites operating periods and data availability during the years from 2015 to 2020. Statistically, S3 and Terra have insignificant differences with RO temperature with an average bias of 3.48 °C and 1.47 °C, respectively, but LC08 and Aqua have significant differences with it. For WVP, Aqua and LC08 have insignificant differences with an average bias of 0.02 kg/m2 and 2.31 kg/m2, respectively, but S3 and Terra have significant differences with RO observations. When comparing LC08 LST data to other satellites, it was found that there were insignificant differences between LC08 and S3 as well as Terra. However, significant differences were observed when comparing LC08 LST data to Aqua. Additionally, significant differences were noted when comparing LC08 WVP data to other satellites. In response to these differences, improvement models have been developed to enhance the estimation of terrestrial data through remote sensing, particularly for satellites that exhibited significant disparities when compared to reference observations (RO).

Analysis of Two Selected Solar Events in 2011 and 2015 With Mars Express Radio Occultation Data

AbstractThe temporal behavior of the Martian ionosphere is highly variable due to various dynamic processes including space weather events. Here, we study the effect of solar flares and coronal mass ejections (CMEs) on the Martian ionosphere for two selected solar events in 2011 and 2015, using the publicly available Mars EXpress (MEX) radio occultation (RO) data (MaRS). We developed a data processing software that converts the calibrated Radio Occultation (RO) Doppler data to scientifically valuable atmospheric profiles. Using this software and previously unexplored MaRS observations, the variations in ionospheric parameters (electron density profiles and total electron content (TEC)) are calculated in order to evaluate the ionospheric changes due to solar flares and CMEs. The RO measurements mostly available 1–4 days apart from the peak events, showed no evident change in the TEC nor in the shape of electron density profiles except for a possible gradual decrease in altitude of M2 (main layer) peak density following the arrival of CMEs. To better quantify the effect of solar events on electron density profiles, RO observations near the time of arrival of solar flares and CMEs are crucial. This can be achieved by frequent RO measurements by various Mars orbiters including spacecraft‐to‐spacecraft measurements assisted by multi‐instrument monitoring of the ionosphere.

Estimation of dusk time F-region electron density vertical profiles using LSTM neural networks: A preliminary investigation

The vertical profile of the ionosphere density plays a significant role in the development of low-latitude Equatorial Plasma Bubbles (EPBs), that in turn lead to ionospheric scintillation which can severely degrade precision and availability of critical users of the Global Navigation Satellite System (GNSS). Accurate estimation of ionospheric delays through vertical electron density profiles is vital for mitigating GNSS errors and enhancing location-based services. The objective of this study is to propose a neural network, trained with radio occultation data from the COSMIC-1 mission, that generates average ionospheric electron density profiles during dusk, focusing on the pre-reversal enhancement of the zonal electric field. Results show that the estimated profiles exhibit a clear seasonal pattern, and reproduce adequately the climatological behavior of the ionosphere, thus presenting strong appeal on ionospheric error attenuation.

A Comparison of Atmospheric Boundary Layer Height Determination Methods Using GNSS Radio Occultation Data

The accurate determination of Atmospheric Boundary Layer Height (ABLH) is crucial in various atmospheric studies and practical applications. In this study, we present a comprehensive comparative analysis of five distinct methods for estimating ABLH using Global Navigation Satellite System (GNSS) Radio Occultation (RO) data. These methods encompass the use of bending angle and refractivity profiles, namely Minimum Gradient methods of the Bending Angle (MGBA) and Refractivity (MGR) profiles, breaking point, Wavelet Covariance Transform (WCT), and Double-Parameter Model Function (DPMF). GNSS-RO data from COSMIC-2 and Spire are used. To establish robust validation, radiosonde data are employed as a reference, ensuring the reliability of our findings. The results reveal notable variations in the performances of these ABLH estimation methods. Specifically, the MGBA, MGR, breaking point, and DPMF methods exhibit strong correlations with the reference data. Conversely, the WCT method displays weaker correlations, higher biases, and elevated root-mean-square-errors, suggesting limitations in capturing the true ABLH. Furthermore, we remove outlier screening to facilitate a comparison of the differences among the five methods. The WCT and DPMF methods can detect strong variations in the profiles near the Earth’s surface and consider them as ABLH. However, these variations are caused by errors. The MGBA method emerges as a reliable and stable option, while the WCT and DPMF methods should be used with caution due to the lower quality of the GNSS-RO profiles near the Earth’s surface.

Determination of the Effective Collision Frequency of Electrons in the E and D Regions of the High-Latitude Ionosphere from Analysis of Radio Occultation Measurements

Collisions between electrons and neutral molecules are of special interest for the physics of the Earth’s ionosphere, in particular, for determining the ionospheric conductivity and current systems in the lower ionosphere of the planet, as well as elucidating the role they play in attenuating radio waves propagating inside the D and E regions of the ionosphere. The effective collision frequency of electrons can be estimated from laboratory studies of electron mobility in atmospheric gases in combination with rocket measurements of temperature and particle density in the Earth’s upper atmosphere, or it can be determined independently from analysis of radio occultation data. We have developed a method for reconstructing the vertical profiles of the absorption coefficient of decimeter (wavelength ~19 cm) radio waves by solving the inverse problem of signal absorption in the D and E regions of the Earth’s ionosphere. Based on the analysis of radio occultation data from the FORMOSAT-3/COSMIC satellites, the altitude profiles of the absorption coefficient of decimeter (DM) radio waves in the planet’s ionosphere during the geomagnetic storm on June 22–23, 2015, were determined. It is known that the absorption coefficient at a given fixed frequency is directly proportional to both the electron density and the collision frequency of electrons with ions and neutrals. Using the data on the vertical profiles of the absorption coefficient of DM radio waves and the electron density reconstructed from the analysis of FORMOSAT-3/COSMIC radio occultations, we estimated the effective collision frequency of electrons in the D and E regions of the Earth’s high-latitude ionosphere. The practical significance of studying the frequency of electron collisions and the effects of radio wave absorption in the D and E regions of the planet’s ionosphere is associated with maintaining the uninterrupted operation of space radio communication and navigation systems.

Estimating Individual Radio Occultation Uncertainties Using the Observations and Environmental Parameters

Abstract Estimation of uncertainties (random error statistics) of radio occultation (RO) observations is important for their effective assimilation in numerical weather prediction (NWP) models. Average uncertainties can be estimated for large samples of RO observations and these statistics may be used for specifying the observation errors in NWP data assimilation. However, the uncertainties of individual RO observations vary, and so using average uncertainty estimates will overestimate the uncertainties of some observations and underestimate those of others, reducing their overall effectiveness in the assimilation. Several parameters associated with RO observations or their atmospheric environments have been proposed to estimate individual RO errors. These include the standard deviation of bending angle (BA) departures from either climatology in the upper stratosphere and lower mesosphere (STDV) or the sample mean between 40 and 60 km (STD4060), the local spectral width (LSW), and the magnitude of the horizontal gradient of refractivity (|∇HN|). In this paper we show how the uncertainties of two RO datasets, COSMIC-2 and Spire BA, as well as their combination, vary with these parameters. We find that the uncertainties are highly correlated with STDV and STD4060 in the stratosphere, and with LSW and |∇HN| in the lower troposphere. These results suggest a hybrid error model for individual BA observations that uses an average statistical model of RO errors modified by STDV or STD4060 above 30 km, and LSW or |∇HN| below 8 km. Significance Statement These results contribute to the understanding of the sources of uncertainties in radio occultation observations. They could be used to improve the effectiveness of these observations in their assimilation into numerical weather prediction and reanalysis models by improving the estimation of their observational errors.

Comparative Assessment of Spire and COSMIC-2 Radio Occultation Data Quality

In this study, we investigate the performances of a commercial Global Navigation Satellite System (GNSS) Radio Occultation (RO) mission and a new-generation RO constellation, i.e., Spire and Constellation Observing System for Meteorology, Ionosphere, and Climate 2 (COSMIC-2), respectively. In the statistical comparison between Spire and COSMIC-2, the results indicate that although the average signal-to-noise ratio (SNR) of Spire is far weaker than that of COSMIC-2, the penetration of Spire is comparable to, and occasionally even better than, that of COSMIC-2. In our analysis, we find that the penetration depth is contingent upon various factors including SNR, GNSS, RO modes, topography, and latitude. With the reanalysis of the European Centre for Medium-Range Weather Forecasts and Radiosonde as the reference data, the identical error characteristics of Spire and COSMIC-2 reveal that overall, the accuracy of Spire’s neutral-atmosphere data products was found to be comparable to that of COSMIC-2.

Influences of Sudden Stratospheric Warming Events on Tropopause Based on GNSS Radio Occultation Data

Sudden Stratospheric Warming (SSW) events have a strong impact on the tropospheric weather and climate. Past researchers have carried out extensive studies investigating the theories of interactions between the stratosphere and the troposphere. However, detailed studies on the influences of the global tropopause are rarely shown. This study uses Global Navigation Satellite System (GNSS) Radio Occultation (RO) data from the years 2007 to 2013 to investigate the influences of different types of SSW events on the tropopause over latitude bands from 30° S to 90° N. It was found that SSW events have strong influences on the tropopause over 60° N–90° N and over 20° N–30° N regions. In 60° N–90° N, SSW events cause a tropopause temperature increase and, therefore, a tropopause height decrease. The increment in the tropopause temperature are more than 10 K and the decrement in the tropopause height is about to 2 km during strong events. Such influences last for about 1.5 months for strong split events and about 10 days for weaker and/or displacement type events. The influences of SSW events on 20° N–30° N are weaker. Only the January 2009 SSW event shows a visible influence on the tropopause layer with a tropopause temperature decrease of about 4 K and a tropopause height increase of about 1 km. Other SSW events share no common characteristics on the tropical tropopause. This is mainly because SSW events are not strong enough to dominate the tropopause variations and other factors, especially the planetary waves in the troposphere, have stronger impacts on the tropopause layer.

Using the Commercial GNSS RO Spire Data in the Neutral Atmosphere for Climate and Weather Prediction Studies

Recently, the NOAA has included GNSS (Global Navigation Satellite System) Radio Occultation (RO) data as one of the crucial long-term observables for weather and climate applications. To include more GNSS RO data in its numerical weather prediction systems, the NOAA Commercial Weather Data Pilot program (CWDP) started to explore the commercial RO data available on the market. After two rounds of pilot studies, the CWDP decided to award the first Indefinite Delivery Indefinite Quantity (IDIQ) contract to GeoOptics and Spire Incs. in 2020. This study examines the quality of Spire RO data products for weather and climate applications. Spire RO data collected from commercial CubeSats are carefully compared with data from Formosa Satellite Mission 7–Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2), the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis (ERA5), and high-quality radiosonde data. The results demonstrate that, despite their generally lower Signal-Noise-Ratio (SNR), Spire RO data show a pattern of lowest penetration height similar to that of COSMIC-2. The Spire and COSMIC-2 penetration heights are between 0.6 and 0.8 km altitude over tropical oceans. Although using different GNSS RO receivers, the precision of Spire STRATOS receivers is of the same quality as those of the COSMIC-2 TriG (Global Positioning System—GPS, GALILEO, and GLObal NAvigation Satellite System—GLONASS) RO Receiver System (TGRS) receivers. Furthermore, the Spire and COSMIC-2 retrieval accuracies are quite comparable. We validate the Spire temperature and water vapor profiles by comparing them with collocated radiosonde observation (RAOB) data. Generally, over the height region between 8 km and 16.5 km, the Spire temperature profiles match those from RS41 RAOB very well, with temperature biases of <0.02 K. Over the height range from 17.8 to 26.4 km, the temperature differences are ~−0.034 K, with RS41 RAOB being warmer. We also estimate the error covariance matrix for Spire, COSMIC-2, and KOMPSAT-5. The results show that the COSMIC-2 estimated error covariance values are slightly more significant than those from Spire over the oceans at the mid-latitudes (45°N–30°N and 30°S–45°S), which may be owing to COSMIC-2 SNR being relatively lower at those latitudinal zones.

Recent Impact of COSMIC-2 with Improved Radio Occultation Data Assimilation Algorithms

Abstract A Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) follow-on constellation, COSMIC-2, was successfully launched into equatorial orbit on 24 June 2019. With an increased signal-to-noise ratio due to improved receivers and digital beam-steering antennas, COSMIC-2 is producing about 5000 high-quality radio occultation (RO) profiles daily over the tropics and subtropics. The initial evaluation of the impact of assimilating COSMIC-2 into NOAA’s Global Forecast System (GFS) showed mixed results, and adjustments to quality control procedures and observation error characteristics had to be made prior to the assimilation of this dataset in the operational configuration in May 2020. Additional changes in the GFS that followed this initial operational implementation resulted in a larger percentage of rejection (∼90%) of all RO observations, including COSMIC-2, in the mid- to lower troposphere. Since then, two software upgrades directly related to the assimilation of RO bending angle observations were developed. These improvements were aimed at optimizing the utilization of COSMIC-2 and other RO observations to improve global weather analyses and forecasts. The first upgrade was implemented operationally in September 2021 and the second one in November 2022. This study describes both RO software upgrades and evaluates the impact of COSMIC-2 with this most recently improved configuration. Specifically, we show that the assimilation of COSMIC-2 observations has a significant impact in improving temperature and winds in the tropics, though benefits also extend to the extratropical latitudes.

Global Ionosphere Modeling Based on GNSS, Satellite Altimetry, Radio Occultation, and DORIS Data Considering Ionospheric Variation

AbstractThe accuracy of ionospheric models estimated by ground‐based multiple global navigation satellite system ionospheric data over regions with sparse tracking stations is not ideal. To improve the accuracy of the estimated ionospheric model, different types of ionospheric data with different combinations were employed for previous studies. However, the ionospheric observational ranges for different types of ionospheric data are not the same. In this study, the accuracy of ionospheric maps generated by ground‐based ionospheric data (ground‐based strategy) and ground‐based ionospheric data combined with data provided by other geodetic measurements normalized by the single‐layer normalization method (multi‐source strategy) were studied. The results showed that the main differences between the ionospheric models estimated by the two strategies occur for data taken over the ocean, which mainly range from −1 to 0 total electron content unit (TECU). When assessed using Jason‐3 vertical total electron content data, the mean root mean square (RMS) value of the ionospheric model estimated by the multi‐source strategy was 5.03 TECU, which is approximately 15% smaller than that estimated by the ground‐based strategy. The maximum reduction in results using the multisource strategy was approximately 25% over different latitudes compared with that of the ground‐based strategy. Furthermore, the self‐consistency evaluation method was employed for evaluation. The results showed that the RMS of the ionospheric model estimated by the multi‐source strategy was 2.41 TECU, which is 3.60% better than that of the ground‐based strategy. The maximum reduction was 15% on different days.

Investigation of the aberrant record from bus GPS receiver onboard FORMOSAT-7/COSMIC-2 satellite constellation in low Earth orbit

The FORMOSAT-7/COSMIC-2 (F7/C2) satellite constellation, a joint Taiwan-US collaboration mission, was launched on 25 June 2019. Following the FORMOSAT-3/COSMIC, the F7/C2 mission aims to provide global users the next-generation global navigation satellite system (GNSS) radio occultation data. The mission consists of six identical small satellites with an inclination angle of 24°, launched to the parking orbit of 720 km, and then transferred to the mission orbit of 550 km. Each satellite bus is equipped with two commercial off-the-shelf global positioning system receivers manufactured by Surrey Satellite Technology Ltd. for positioning, navigation, and timing functions. This receiver generates “GPS Receiver (bus GPSR) Warm Start Requested (WSR)” messages after commissioning, which could be considered as a kind of single event effect (SEE). More than 7000 WSR events have been recorded in 3 years after launch, providing another opportunity to investigate the SEE with the commercial off-the-shelf bus GPSR. In addition, the F7/C2 mission is the first satellite constellation that allows scientists to investigate the SEE using the same constellation at different altitudes of Low-Earth orbit. This study investigates and discusses the cumulative rate, distribution of WSRs, and their relationships between altitude and geomagnetic field. The statistical method is applied to define a warning area of the high probability occurring the SEE of this kind of component based on the geomagnetic field, whose corresponding strength is around 21,200/20,550 nT and close to the South Atlantic Anomaly area. Moreover, the relationships with space weather, including the disturbance of the Earth's magnetic field (Kp index) and the solar activity (F10.7 index), are also analyzed and reported.

A comparison of FORMOSAT-3/COSMIC radio occultation and ionosonde measurements in sporadic E detection over mid- and low-latitude regions

The investigation of sporadic E or Es layers typically relies on ground-based or satellite data. This study compares the Es layers recorded in ionograms with those detected using GNSS L1 signal-to-noise ratio data from FORMOSAT-3/COSMIC radio occultation at mid and low latitudes. GPS radio occultation measurements of Es layers, during an 11-year time span of 2007–2017, within a 2° latitude × 5° longitude grid around each ionosonde site are compared to the Es recordings of the ionosonde. By comparing multi-year radio occultation data with recordings from six ionosonde stations at mid and low latitudes, it was discovered that at least 20% of the Es layer detection results between each ionosonde and its crossing GPS radio occultation measurements did not agree. The results show that the agreement between the two methods in Es detection is highly dependent on the season and local time. This study suggests that Es layer recordings from ground-based ionosonde observations have the best agreement with the Es layers detected by radio occultation data during daytime and local summers. The difference in the Es detection mechanisms between the two methods can explain the inconsistency between Es events measured by these two methods. The detection of Es layers in ionograms relies on the high plasma concentration in the E region, whereas signal scintillations caused by a large vertical gradient of the plasma density in the E region are considered a sign of Es occurrence in satellite techniques.

Spire RO Thermal Profiles for Climate Studies: Initial Comparisons of the Measurements from Spire, NOAA-20 ATMS, Radiosonde, and COSMIC-2

Global Navigation Satellite System (GNSS) Radio Occultation (RO) data play an essential role in improving numerical weather prediction (NWP) and monitoring climate change. The NOAA Commercial RO Purchase Program (CDP) purchased RO data provided by Spire Global Inc. To ensure the data quality from Spire Global Inc. is consistent with other RO missions, we need to quantify their accuracy and retrieval uncertainty carefully. In this work, Spire Wet Profile (wet temperature profile) data from 7 September 2021 to 31 October 2022, processed by the University Corporation for Atmospheric Research (UCAR), and COSMIC-2 (Constellation Observing System for Meteorology, Ionosphere, and Climate-2/Formosa Satellite Mission 7) data are evaluated through comparison with NOAA-20 Advanced Technology Microwave Sounder (ATMS) microwave sounder measurements and collocated RS41 radiosonde measurements. Through the Community Radiative Transfer Model (CRTM) simulation, we convert the Spire and COSMIC-2 RO retrievals to ATMS brightness temperature (BT) at sounding channels CH07 to CH14 (temperature channels), with weighting function peak heights from 8 km to 35 km, and CH19 to CH22 (water vapor channels), with weighting function peak heights ranging from 3.2 km to 6.7 km, and compare the simulations with the collocated NOAA-20 ATMS measurements over ocean. Using ATMS observations as references, Spire and COSMIC-2 BTs agree well with ATMS within 0.07 K for CH07-14 and 0.20 K for CH19-22. The trends between Spire and COSMIC-2 are consistent within 0.07 K/year over the oceans for ATMS CH07-CH13 and CH19-22, indicating that Spire/COSMIC-2 wet profiles are, in general, compatible with each other over oceans. The RO retrievals and RS41 radiosonde observation (RAOB) comparison shows that above 0.2 km altitude, RS41 RAOB matches Spire/COSMIC-2 temperature profiles well with a temperature difference of <0.13 K, and the trends between Spire and COSMIC-2 are consistent within 0.08 K/year over land, indicating that Spire/COSMIC-2 wet profiles are overall compatible with each other through RS41 RAOB measurements over land. In addition, the consistency of Spire and COSMIC-2 based on different latitude intervals, local times, and signal-to-noise ratios (SNRs) through ATMS was evaluated. The results show that the performance of Spire is comparable to COSMIC-2, even though COSMIC-2 has a higher SNR. The high quality of RO profiles from Spire is expected to improve short- and medium-range global numerical weather predictions and help construct consistent climate temperature records.

An evaluation of GNSS radio occultation atmospheric profiles from Sentinel-6

The Global Navigation Satellite System (GNSS) radio occultation (RO) significantly impacts climate change monitoring. GNSS RO is based mainly on the refraction of the GNSS signals transmitted by satellites and received by a receiver on a Low Earth Orbit (LEO) satellite. Sentinel-6 (S6) is the first satellite from the Sentinel missions to retrieve atmospheric parameters using the RO technique. Temperature, pressure, and water vapor profiles can be retrieved using S6 based on the refraction of GNSS signals. This research evaluated the distribution and coverage of GNSS RO data from S6. This research also aims to evaluate the results of GNSS RO products, especially temperature from S6, from January 2022 to June 2022. The results of S6 are compared with the results from the Meteorological Operational Polar Satellite series (MetOp-B and MetOp-C) and the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5) dataset. Our results indicate that S6 allows for many RO observations in space by receiving signals from the Global Positioning System (GPS) and Russia's Global Navigation Satellite System (GLONASS). S6 provides monthly RO profiles from GPS that are greater than monthly RO profiles from GLONASS. Results show that S6 provides more GNSS RO data than MetOp-B and MetOp-C. The findings indicate that the distribution of GNSS RO events from S6 over the globe is uniform in longitude and non-uniform at latitude, with the minimum number at the poles. This paper makes a significant contribution to exploring the atmosphere with the S6 as a promising satellite for weather forecasting and monitoring climate change.

Prospects for meteotsunami detection in earth’s atmosphere using GNSS observations

We study, for the first time, the physical coupling and detectability of meteotsunamis in the earth’s atmosphere. We study the June 13, 2013 event off the US East Coast using Global Navigation Satellite System (GNSS) radio occultation (RO) measurements, Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) temperatures, and ground-based GNSS ionospheric total electron content (TEC) observations. Hypothesizing that meteotsunamis also generate gravity waves (GWs), similar to tsunamigenic earthquakes, we use linear GW theory to trace their dynamic coupling in the atmosphere by comparing theory with observations. We find that RO data exhibit distinct stratospheric GW activity at near-field that is captured by SABER data in the mesosphere with increased vertical wavelength. Ground-based GNSS-TEC data also detect a far-field ionospheric response 9 h later, as expected by GW theory. We conclude that RO measurements could increase understanding of meteotsunamis and how they couple with the earth’s atmosphere, augmenting ground-based GNSS TEC observations.

Efficient collocation of global navigation satellite system radio occultation soundings with passive nadir microwave soundings

Abstract. Radio occultation (RO) using the global navigation satellite system (GNSS) can be used to infer atmospheric profiles of microwave refractivity in the Earth's atmosphere. GNSS RO data are now assimilated into numerical weather prediction models and used for climate monitoring. New remote sensing applications are being considered that fuse GNSS RO soundings and passive nadir-scanned radiance soundings. Collocating RO soundings and nadir-scanned radiance soundings, however, is computationally expensive, especially as new commercial GNSS RO constellations greatly increase the number of global daily RO soundings. This paper develops a new and efficient technique, called the “rotation–collocation method”, for collocating RO and nadir-scanned radiance soundings in which all soundings are rotated into the time-dependent reference frame in which the nadir sounder's scan pattern is stationary. Collocations with RO soundings are then found when the track of an RO sounding crosses the line corresponding to the nadir sounder's scan pattern. When applied to finding collocations between RO soundings from COSMIC-2, Metop-B-GRAS, and Metop-C-GRAS and the passive microwave (MW) soundings of the Advanced Technology Microwave Sounder (ATMS) on NOAA-20 and Suomi-NPP and the Advanced Microwave Sounding Unit (AMSU-A) on Metop-B and Metop-C for the month of January 2021, the rotation–collocation method proves to be 99.0 % accurate and is hundreds to thousands of times faster than traditional approaches to finding collocations.

Real-Time Tomographic Inversion of Truncated Ionospheric GNSS Radio Occultations

This paper presents a new way of combining Abel inversion and the Chapman model with a linearly increasing scale height to retrieve ionospheric electron density vertical profiles from truncated-sounding radio-occultation data. A linear Vary–Chap model is used to cover the blind region due to data truncation, with parameters estimated by enumeration of the possible values in a grid centered around a set of parameters compatible with ionospheric physics. The resulting electron density is estimated with its corresponding error from the linear least-squares solution presenting the smaller post-fit residual on the input GNSS carrier-phase measurements. The results, tested on a set of representative GNSS RO measurements obtained by COSMIC/FORMOSAT-3, show that this method can retrieve EDVPs with a predominant absolute and relative error of 1010e−m−3 and 5%, respectively, and in less than 10 s per profile, which makes this method suitable for near real-time applications in upcoming missions such as EUMETSAT Polar System-Second Generation.

A Disturbance Frequency Index in Earthquake Forecast Using Radio Occultation Data

Earthquake forecasting is the process of forecasting the time, location, and magnitude of an earthquake, hoping to gain some time to prepare to reduce the disasters caused by earthquakes. In this paper, the possible relationship between the maximum electron density, the corresponding critical frequency, and the occurrence of earthquakes is explored by means of radio occultation data based on mechanism analysis and actual earthquake-nearby data. A new disturbance frequency index is proposed in this paper as a novel method to help forecast earthquakes. Forecasting of the location and timing of earthquakes is based on the connection between proven new frequency distributions and earthquakes. The effectiveness of this index is verified by backtracking observation around the 2022 Ya’an earthquake. Using this index, occultation data can forecast the occurrence of earthquakes five days ahead of detection, which can help break the bottleneck in earthquake forecasting.

Modeling the Topside Ionosphere Effective Scale Height through In Situ Electron Density Observations by Low-Earth-Orbit Satellites

In this work, we aim to characterize the effective scale height at the ionosphere F2-layer peak (H0) by using in situ electron density (Ne) observations by Langmuir Probes (LPs) onboard the China Seismo-Electromagnetic Satellite (CSES—01). CSES—01 is a sun-synchronous satellite orbiting at an altitude of ~500 km, with descending and ascending nodes at ~14:00 local time (LT) and ~02:00 LT, respectively. Calibrated CSES—01 LPs Ne observations for the years 2019–2021 provide information in the topside ionosphere, whereas the International Reference Ionosphere model (IRI) provides Ne values at the F2-layer peak altitude for the same time and geographical coordinates as CSES—01. CSES—01 and IRI Ne datasets are used as anchor points to infer H0 by assuming a linear scale height in the topside representation given by the NeQuick model. COSMIC/FORMOSAT—3 (COSMIC—1) radio occultation (RO) data are used to constrain the vertical gradient of the effective scale height in the topside ionosphere in the linear approximation. With the CSES—01 dataset, we studied the global behavior of H0 for daytime (~14:00 LT) and nighttime (~02:00 LT) conditions, different seasons, and low solar activity. Results from CSES—01 observations are compared with those obtained through Swarm B satellite Ne-calibrated measurements and validated against those from COSMIC—1 RO for similar diurnal, seasonal, and solar activity conditions. H0 values modeled by using CSES—01 and Swarm B-calibrated observations during daytime both agree with corresponding values obtained directly from COSMIC—1 RO profiles. Differently, H0 modeling for nighttime conditions deserves further investigation because values obtained from both CSES—01 and Swarm B-calibrated observations show remarkable and spatially localized differences compared to those obtained through COSMIC—1. Most of the H0 mismodeling for nighttime conditions can probably to be attributed to a sub-optimal spatial representation of the F2-layer peak density made by the underlying IRI model. For comparison, H0 values obtained with non-calibrated CSES—01 and Swarm B Ne observations are also calculated and discussed. The methodology developed in this study for the topside effective scale height modeling turns out to be applicable not only to CSES—01 satellite data but to any in situ Ne observation by low-Earth-orbit satellites orbiting in the topside ionosphere.