Radio Occultation Measurements Research Articles

Development of the Ionospheric E‐Region Prompt Radio Occultation Based Electron Density (E‐PROBED) Model

AbstractThis work reports the development of the first version of the E‐region Prompt Radio Occultation Based Electron Density (E‐PROBED) Model. This is an empirical model of E‐region electron density (Ne) between 90 and 120 km developed using radio occultation measurements from the COSMIC‐1 mission. This first version captures more than 80% of the observed variability in monthly‐mean latitude‐local time‐altitude E‐region Ne profiles but it does not account for longitudinal variability at constant local‐time. This work also reports a validation of E‐PROBED simulations through comparisons with ionosondes and incoherent scatter radar (ISR) E‐region Ne profiles. E‐PROBED generally agrees with these ground‐based observations during day‐time. During night‐time, there is a large disparity between E‐PROBED and ISR values. Finally, this work compares E‐PROBED with E‐region Ne simulated by the International Reference Ionosphere (IRI) and the Specified Dynamics—Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD‐WACCM‐X). One of the main differences amongst these models is on the simulation of variabilities that cannot be attributed to photoionization. IRI barely simulates any variability not driven by photoionization. Both E‐PROBED and SD‐WACCM‐X simulates variability not driven by photoionization. Another main difference is in the absolute magnitude of night‐time E‐region Ne values. Both IRI and SD‐WACCM‐X are substantially lower than E‐PROBED. This work first concludes that E‐PROBED can conveniently provide E‐region Ne latitude—local time variabilities and structures that COSMIC‐1 observes. This work also concludes that E‐region Ne have significant non‐photoionization driven variabilities.

Investigation of the global climatologic performance of ionospheric models utilizing in-situ Swarm satellite electron density measurements

The Swarm constellation is a triplet of satellites, flying, since their final configuration reached in April 2014, at altitudes of about 420 to 490 km (the lower pair) and about 500 to 530 km (the upper satellite). All the three satellites provide in-situ measurements of the plasma density in the topside ionosphere using Langmuir Probe sensors onboard the Electrical Field Instrument. The present study is a comprehensive investigation into the climatologic performance of three ionospheric models when compared to the Swarm satellite in-situ measurements. The models are the International Reference Ionosphere (IRI) model, a quick run ionospheric electron density model (NeQuick), and a 3-dimensional electron density model based on artificial neural network training of COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate) satellites radio occultation measurements (3D-NN). The mean monthly quiet-time latitudinal profile of Swarm measurements was computed by binning the Swarm electron density measurements in 15-degree longitudes starting from longitude −180° in steps of 15° to 180°, and corresponding model predictions were obtained. The data used in the study covers the years 2014, 2016, 2019, and 2022, capturing various phases of the solar activity cycle. Results from the study show that modelled electron density predictions from all three climatologic models are fairly good representations of the Swarm satellite measurements, with some exceptions in which the models underestimate or overestimate the Swarm satellite values. The IRI model performed best at the northern hemisphere mid latitude, and it overestimated the Swarm measurements at altitudes of ∼ 450 km, especially at the southern hemisphere mid and high latitudes. The NeQuick performed best during the night times, and it overestimated the Swarm measurements, especially at the mid latitudes. The NeQuick was also observed to overestimate the Swarm measurements during the winter solstices at both hemispheres, which is June solstice in the southern hemisphere and December solstice in the northern hemisphere. Overall, the 3D-NN model most often performed better than the IRI model and the NeQuick, especially during the day times and during the high solar activity year (2014), but it underestimated the Swarm measurements, especially at the low and mid latitudes. For all categories explored in the study, the 3D-NN consistently performed better than the other two models. The NeQuick performed better than the IRI model at altitude of satellite B, while the IRI model performed slightly better than the NeQuick at altitude of satellites A and C. The NeQuick also performed better than the IRI model in the local time category, whereas the IRI model performed better than the NeQuick in categories of season, solar activity, and longitudinal sector.

Characterizing global equatorial sporadic-E layers through COSMIC GNSS radio occultation measurements

Characterizing global equatorial sporadic-E layers through COSMIC GNSS radio occultation measurements

Variations of Kelvin Waves around the Tropical Tropopause Inversion Layer from GNSS RO Measurements

Abstract In the present study, the tropical tropopause inversion layer (TIL) Kelvin waves are extracted from the Global Navigation Satellite System (GNSS) radio occultation (RO) temperature data of multiple missions from January 2007 to December 2020. We focus on the variations of TIL Kelvin waves in two longitude regions, the Maritime Continent (MC; 90°–150°E) and the Pacific Ocean (PO; 170°–230°E). The results show that over both regions, ENSO leads to the opposite variations of TIL Kelvin wave temperature amplitude during different ENSO phases. Specifically, during La Niña, the strong (weak) deep convection over MC (PO) leads to strengthened (weakened) static stability. With enhanced easterly (westerly) winds and strengthened (weakened) static stability, the TIL Kelvin wave temperature amplitudes are stronger (weaker) over MC (PO). The opposite phenomenon occurs during El Niño. The zonal-mean zonal winds affect TIL Kelvin wave temperature amplitudes by two mechanisms. First, the prevalence of easterlies (westerlies) in the upper troposphere affects the upward propagation of Kelvin waves, resulting in stronger (weaker) TIL Kelvin wave temperature amplitudes over MC (PO). Second, the TIL Kelvin wave temperature amplitude peaks about 2 months before the zero-wind line of the descending westerly QBO phase occurs, due to dissipation on the critical line. Additionally, the rapid increase of zonal-mean static stability significantly affects the annual variation of TIL Kelvin wave temperature amplitudes. They both reach maxima during DJF and minima during JJA, which should be related to the annual cycles of temperature and ozone mixing ratio in the TIL. Significance Statement Recent studies indicate that the Kelvin wave temperature amplitudes in the tropical tropopause inversion layer (TIL) exhibit distinct characteristics compared with those in other height levels, while the modulation mechanisms of the TIL Kelvin waves need further investigation. The present study aims to study the differences in the variabilities and the modulation factors of TIL Kelvin waves over two longitude regions. Our findings suggest that the different responses of background conditions during ENSO phases influence the spatiotemporal distribution of the TIL Kelvin waves. Besides, the zonal winds and the static stability significantly affect the temporal variations of TIL Kelvin waves. Our work fills the research gap of TIL Kelvin waves and contributes to understanding the dynamics of tropical tropopause variations.

Data Assimilation Impact of GNSS RO Measurements from Cube Satellites on Arctic Weather Forecasts

Abstract In this study, the effects of assimilating Global Navigation Satellite System (GNSS) radio occultation (RO) observations from existing and recently added commercial cube satellites on analyses and forecasts over the Arctic region were investigated by conducting observing system experiments (OSEs). Profiles of refractivity were assimilated with a local observation operator using the three-dimensional variational method. The analyses and forecasts from the OSEs were verified against ERA5 reanalysis, radiosonde observations, and buoy observations. In addition to the averaged impact on forecast skill, the impact of GNSS RO observations was further examined for an individual Arctic cyclone case, focusing on the added value of the cube satellite data. The effects of GNSS RO observations from existing satellites on analyses and forecasts over the Arctic region are positive, and the assimilation of GNSS RO observations from cube satellites leads to additional improvements, particularly for temperature in the upper troposphere and lower stratosphere (UTLS). Temperature biases in the UTLS are significantly reduced in the analyses, and the improved analyses result in better forecasts of upper-level potential vorticity and cyclone development when GNSS RO observations from cube satellites are assimilated. This result demonstrates the potential of GNSS RO data from cube satellites to enhance forecasts over the Arctic region.

Forward operator for polarimetric radio occultation measurements

Abstract. Global Navigation Satellite System (GNSS) polarimetric radio occultation (PRO) observations sense the presence of hydrometeor particles along the ray path by measuring the difference of excess phases in horizontally and vertically polarised carrier waves. As a first step towards using these observations in data assimilation and model diagnostics, a forward operator for the GNSS-PRO observable ΦDP (polarimetric differential phase shift) has been implemented by extending the existing two-dimensional forward operator for radio occultation bending-angle observations. Evaluation of heavy-precipitation cases showed that the implemented forward operator can simulate the observed ΦDP in synoptic-scale atmospheric river (AR) cases very accurately. For tropical cyclone cases it is more challenging to produce reasonable ΦDP simulations, due to the high sensitivity of ΦDP with respect to displacement of the position of the tropical cyclones. It was also found that snow is the dominant contributor to the simulated ΦDP and that the ability to compute the ray paths in two dimensions is essential to accurately simulate ΦDP.

Global zenith wet delay modeling with surface meteorological data and machine learning

The tropospheric delay is a major error source for space geodetic techniques, and the performance of its modeling is significantly limited due to the high spatiotemporal variability of the moisture in the lower atmosphere. In this study, global modeling of the tropospheric zenith wet delay (ZWD) was realized based on surface meteorological data obtained from radiosondes and Global Positioning System (GPS) radio occultation (RO) measurements through the random forest (RF) and backpropagation neural network (BPNN) regression analysis. The modeling performance was further validated based on two kinds of global atmospheric profiles for the year 2020. Our results show that the ZWD modeling accuracy gained by two machine learning regression approaches is significantly improved by taking into account surface meteorological parameters, especially the surface water vapor pressure when compared to the Global Pressure and Temperature 3 (GPT3) model. When surface meteorological data are available, the RF-B model yields ZWD estimations with an overall agreement of 3.1 cm in comparison with the sounding profiles and 2.4 cm in contrast to the GPS RO atmospheric profiles. The RF-B is superior to other models based on surface meteorological parameters for ZWD calculation, e.g., the accuracy improves by 21.8–23.8% against the approach by Saastamoinen and 7–12.2% against the formula by Askne and Nordius.

One‐Dimensional Variational Ionospheric Retrieval Using Radio Occultation Bending Angles: 2. Validation

AbstractCulverwell et al. (2023, https://doi.org/10.1029/2023SW003572) described a new one‐dimensional variational (1D‐Var) retrieval approach for ionospheric GNSS radio occultation (GNSS‐RO) measurements. The approach maps a one‐dimensional ionospheric electron density profile, modeled with multiple “Vary‐Chap” layers, to bending angle space. This paper improves the computational performance of the 1D‐Var retrieval using an improved background model and validates the approach by comparing with the COSMIC‐2 profile retrievals, based on an Abel Transform inversion, and co‐located (within 200 km) ionosonde observations using all suitable data from 2020. A three or four layer Vary‐Chap in the 1D‐Var retrieval shows improved performance compared to COSMIC‐2 retrievals in terms of percentage error for the F2 peak parameters (NmF2 and hmF2). Furthermore, skill in retrieval (compared to COSMIC‐2 profiles) throughout the bottomside (∼90–300 km) has been demonstrated. With a single Vary‐Chap layer the performance is similar, but this improves by approximately 40% when using four‐layers.

One‐Dimensional Variational Ionospheric Retrieval Using Radio Occultation Bending Angles: 1. Theory

AbstractA new one‐dimensional variational (1D‐Var) retrieval method for ionospheric GNSS radio occultation (GNSS‐RO) measurements is described. The forward model implicit in the retrieval calculates the bending angles produced by a one‐dimensional ionospheric electron density profile, modeled with multiple “Vary‐Chap” layers. It is demonstrated that gradient based minimization techniques can be applied to this retrieval problem. The use of ionospheric bending angles is discussed. This approach circumvents the need for Differential Code Bias (DCB) estimates when using the measurements. This new, general retrieval method is applicable to both standard GNSS‐RO retrieval problems, and the truncated geometry of EUMETSAT's Metop Second Generation (Metop‐SG), which will provide GNSS‐RO measurements up to about 600 km above the surface. The climatological a priori information used in the 1D‐Var is effectively a starting point for the 1D‐Var minimization, rather than a strong constraint on the final solution. In this paper the approach has been tested with 143 COSMIC‐1 measurements. We find that the method converges in 135 of the cases, but around 25 of those have high “cost at convergence” values. In the companion paper (Elvidge et al., 2023), a full statistical analysis of the method, using over 10,000 COSMIC‐2 measurements, has been made.

Sub‐Hourly Observations of Dust Storm Growth, Lee Waves, and Lyot Crater, by the EMM Camera EXI

AbstractWe explore a sequence of 13 unique high‐cadence images of a dust storm, from the Emirates Mars Mission (EMM). The Emirates eXploration Imager camera took these images in less than 8 hr on 18 December 2022 (Martian Year 36, solar longitude 356°). Most of these images are separated by a time difference of half an hour. The region of interest extends from Lyot crater to the east. During the morning, the EMM images show lee waves (atmospheric gravity waves). In the late morning, the lee waves rapidly change into clearly distinct dust storm texture/convective features. We track the evolution of both lee waves and a local dust storm between sunrise and mid‐afternoon. Also, we relate our observations to atmospheric dynamics. Our analysis is supported by the Mars Climate Database and radio occultation measurement data.

Convective tropopause over the tropics: Climatology, seasonality, and inter-annual variability inferred from long-term FORMOSAT-3/COSMIC-1 RO data

The present study aims to provide a comprehensive analysis of the convective tropopause height (COTH) and temperature (COT-T) over the tropics using long-term radio occultation (RO) measurements from the Formosa Satellite Mission 3–Constellation Observing System for Meteorology, Ionosphere and Climate (FORMOSAT-3/COSMIC) from July 2006 to December 2018. The study reveals that the COT-H (COT-T) exhibits a significant longitudinal structure with a maximum (minimum) over the Western Pacific and two secondary maxima (minima) over the American and African regions. The seasonal cycles of COT-H and COT-T are also prominent, with higher and colder values over the Western Pacific in boreal winter and over the Indian region in boreal summer. This seasonality in COT-H and COT-T is consistent with the seasonal process of tropical convective activity and is explained by outgoing longwave radiation (OLR) data. The correlation analysis shows that COT-H (COT-T) is negatively (positively) correlated with OLR, with the highest correlation found in the Indian region. Additionally, the study investigates the relationship between COT-H (COT-T) and upper tropospheric water vapor (UTWV) over the various regions in the tropics. Over the Indian region, the COT-H (COT-T) exhibits a positive (negative) correlation with UTWV, while over the Western Pacific, American, and African regions, it displays a negative (positive) correlation. Finally, the spatial distribution of the convective tropopause in El Niño and La Niña winters was examined. Results show that the COT-H (COT-T) is higher (lower) over the central and eastern Pacific Ocean and lower (higher) over the maritime continent during El Niño compared with La Niña.

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.

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

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

Radio Occultation Measurements of Europa's Ionosphere From Juno's Close Flyby

AbstractOn 29 September 2022 the Juno spacecraft flew within 354 km of Europa's surface while several instruments probed the moon's surroundings. During the close flyby, radio occultations were performed by collecting single‐frequency Doppler measurements. These investigations are essential to the study of Europa's ionosphere and represent the first repeat sampling of any set of conditions since the Galileo era. Ingress measurements resulted in a marginal detection with a peak ionospheric density of 4,000 ± 3,700 cm−3 (3σ) at 22 km altitude. A more significant detection emerged on egress, with a peak density of 6,000 ± 3,000 cm−3 (3σ) at 320 km altitude. Comparison with Galileo measurements reveals a consistent picture of Europa's ionosphere, and confirms its dependence on illumination conditions and position within Jupiter's magnetosphere. However, the overall lower densities measured by Juno suggest a dependence on time of observation, with implications for the structure of the neutral atmosphere.

GNSS radio occultation excess-phase processing for climate applications including uncertainty estimation

Abstract. Earth observation from space provides a highly valuable basis for atmospheric and climate science, in particular also through climate benchmark data from suitable remote sensing techniques. Measurements by global navigation satellite system (GNSS) radio occultation (RO) qualify to produce such benchmark data records as they globally provide accurate and long-term stable datasets for essential climate variables (ECVs) such as temperature. This requires a rigorous processing of the raw RO measurements to ECVs, with narrow uncertainties. In order to fully exploit this potential, Wegener Center's Reference Occultation Processing System (rOPS) Level 1a (L1a) processing subsystem includes uncertainty estimation in both precise orbit determination (POD) and excess-phase profile derivation. Here we introduce the new rOPS L1a excess-phase processing, the first step in the RO profiles retrieval down to atmospheric profiles, which extracts the atmospheric excess phase from raw SI-traceable RO measurements. This excess-phase processing, for itself algorithmically concise, includes integrated quality control and uncertainty estimation, requiring a complex framework of various subsystems that we first introduce before describing the implementation of the core algorithms. The quality control and uncertainty estimation, computed per RO event, are supported by reliable forward-modeled excess-phase profiles based on the POD orbit arcs and collocated short-range forecast profiles of the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA5). The quality control removes or alternatively flags excess-phase profiles of insufficient or degraded quality. The uncertainty estimation accounts both for relevant random- and systematic-uncertainty components, and the resulting (total) uncertainty profiles serve as a starting point for the subsequent uncertainty propagation through the retrieval processing chain down to the atmospheric ECV profiles. We also evaluated the quality and reliability of the resulting excess-phase profiles based on Metop-A/B/C (Meteorological Operational) RO datasets for three 3-month periods in 2008, 2013, and 2020 by way of a sensitivity analysis for three representative atmospheric layers (tropo-, strato-, mesosphere), investigating consistency with ERA5-derived profiles, influences of different orbit and clock inputs, and consistency across the different Metop satellites. These consistencies range from centimeter to submillimeter levels, indicating that the new processing can provide highly accurate and robust excess-phase profiles. Furthermore, cross-evaluation and intercomparison with excess-phase data from the established data providers EUMETSAT (European Organisation for the Exploitation of Meteorological Satellites) and UCAR (University Corporation for Atmospheric Research) revealed subtle discrepancies but overall very close agreement, with larger differences compared to UCAR in the boundary layer. The new rOPS L1a processing can hence be considered capable of producing reliable long-term data records including uncertainty estimation for the benefit of climate applications.

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.

Undifferenced processing of COSMIC-2 radio occultation measurements with 1-Hz LEO clock corrections

Undifferenced processing of COSMIC-2 radio occultation measurements with 1-Hz LEO clock corrections

Ionospheric irregularity reconstruction using multisource data fusion via deep learning

Abstract. Ionospheric sporadic E layers (Es) are intense plasma irregularities between 80 and 130 km in altitude and are generally unpredictable. Reconstructing the morphology of sporadic E layers is not only essential for understanding the nature of ionospheric irregularities and many other atmospheric coupling systems, but is also useful for solving a broad range of demands for reliable radio communication of many sectors reliant on ionosphere-dependent decision-making. Despite the efforts of many empirical and theoretical models, a predictive algorithm with both high accuracy and high efficiency is still lacking. Here we introduce a new approach for Sporadic E Layer Forecast using Artificial Neural Networks (SELF-ANN). The prediction engine is trained by fusing observational data from multiple sources, including a high-resolution ERA5 reanalysis dataset, Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) radio occultation (RO) measurements, and integrated data from OMNIWeb. The results show that the model can effectively reconstruct the morphology of the ionospheric E layer with intraseasonal variability by learning complex patterns. The model obtains good performance and generalization capability by applying multiple evaluation criteria. The random forest algorithm used for preliminary processing shows that local time, altitude, longitude, and latitude are significantly essential for forecasting the E-layer region. Extensive evaluations based on ground-based observations demonstrate the superior utility of the model in dealing with unknown information. The presented framework will help us better understand the nature of the ionospheric irregularities, which is a fundamental challenge in upper-atmospheric and ionospheric physics. Moreover, the proposed SELF-ANN can make a significant contribution to the development of the prediction of ionospheric irregularities in the E layer, particularly when the formation mechanisms and evolution processes of the Es layer are not well understood.

Stratospheric gravity wave potential energy and tropospheric parameters relationships over South America: a study using COSMIC-2 and METOP radio occultation measurements

Using COSMIC-2 and METOP radio occultation measurements during the years 2020 and 2021, the study presents the first direct and independent relationship between the potential energy (Ep) in the stratosphere, precipitable water vapour (PWV), tropopause heights (TPH), and cold-point heights (CPH) over South America. The South American continent comprises the tropical region, the Andes Mountain range, and mid-latitude climates. The seasonal mean of the potential energy (Ep), the PWV, and the tropopause parameters height (TPH and CPH) were obtained to investigate the relationship between the stratospheric gravity wave (SGW) Ep and the tropospheric parameters (PWV, TPH, and CPH). Around the Andes Mountains to the east, there is significantly less water vapour (PWV < 10 mm) and a relatively high gravity wave Ep (Ep > 8 kJ kg−1). A good correlation of variability was found between the PWV and the lower SGW Ep in summer over the tropical region (± 20◦). Generally, good and strong correlations were observed in the summer and spring, with negative/no correlations in the winter in 2020 and 2021. Also, good and strong correlations between SGW, PWV, and TPH were observed in the summer at {20}^{o}N - {10}^{o}N in 2020 and 2021. Our result demonstrated the possibility that convective activity was a major driver of the tropical gravity waves over South America. In the subtropical (30^circ–40^circ) region, especially in the winter, the tropospheric parameters make little or no contribution to gravity wave activity in the region. The CPH generally showed a no/negative with SGW over the South American tropics. The SGW activities in the tropical region showed an impact on the structure of the tropopause parameters, which could be a result of the convective activity in this region.Graphical

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.