Stability of CubeSat Clocks and Their Impacts on GNSS Radio Occultation

Global Navigation Satellite Systems’ radio occultation (GNSS-RO) provides the upper troposphere-lower stratosphere (UTLS) vertical atmospheric profiles that are complementing radiosonde and reanalysis data. Such data are employed in the numerical weather prediction (NWP) models used to forecast global weather as well as in climate change studies. Typically, GNSS-RO operates by remotely sensing the bending angles of an occulting GNSS signal measured by larger low Earth orbit (LEO) satellites. However, these satellites are faced with complexities in their design and costs. CubeSats, on the other hand, are emerging small and cheap satellites; the low prices of building them and the advancements in their components make them favorable for the GNSS-RO. In order to be compatible with GNSS-RO requirements, the clocks of the onboard receivers that are estimated through the precise orbit determination (POD) should have short-term stabilities. This is essential to correctly time tag the excess phase observations used in the derivation of the GNSS-RO UTLS atmospheric profiles. In this study, the stabilities of estimated clocks of a set of CubeSats launched for GNSS-RO in the Spire Global constellation are rigorously analysed and evaluated in comparison to the ultra-stable oscillators (USOs) onboard the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC-2) satellites. Methods for improving their clock stabilities are proposed and tested. The results (i) show improvement of the estimated clocks at the level of several microseconds, which increases their short-term stabilities, (ii) indicate that the quality of the frequency oscillator plays a dominant role in CubeSats’ clock instabilities, and (iii) show that CubeSats’ derived UTLS (i.e., tropopause) atmospheric profiles are comparable to those of COSMIC-2 products and in situ radiosonde observations, which provided external validation products. Different comparisons confirm that CubeSats, even those with unstable onboard clocks, provide high-quality RO profiles, comparable to those of COSMIC-2. The proposed remedies in POD and the advancements of the COTS components, such as chip-scale atomic clocks and better onboard processing units, also present a brighter future for real-time applications that require precise orbits and stable clocks.

On the quality of RS41 radiosonde descent data

Radiosonde descent profiles have been available from tens of stations for several years now – mainly from Vaisala RS41 radiosondes. They have been compared with the ascent profiles, with ECMWF short-range forecasts and with co-located radio occultation retrievals. Over this time, our understanding of the data has grown, and the comparison has also shed some light on radiosonde ascent data. The fall rate is very variable and is an important factor, with high fall rates being associated with temperature biases, especially at higher altitudes. Ascent winds are affected by pendulum motion; on average, descent winds are less affected by pendulum motion and are smoother. It is plausible that the true wind variability in the vertical lies between that shown by ascent and descent profiles. This discrepancy indicates the need for reference wind measurements. With current processing, the best results are for radiosondes with parachutes and pressure sensors. Some of the wind, temperature and humidity data are now assimilated in the ECMWF forecast system.

Sounding of sporadic E layers from CSES radio occultation and comparing with ionosonde measurements

GNSS radio occultation (RO) plays an important role in ionospheric electron density inversion and sounding of sporadic E layers. As the China's first electromagnetic satellite, China Seismo Electromagnetic Satellite (CSES) has collected the RO data from both GPS and BDS-2 satellites since March 2018. In this study, we extracted the carrier to noise density ratio (CNR) data of CSES and calculated the standard deviation of normalized CNR. A new criterion is developed to determine the Es events, that is when the mean value of the absolute value of the difference between the normalized CNR is greater than 3 times of the standard deviation. The statistics show that sporadic E layers have strong seasonal variations with highest occurrence rates in summer season at middle latitudes. It is also found that the occurrence height of Es is mainly located at 90–110 km, and the period of local time 15:00–18:00 is the high incidence period of Es. In addition, the geometric altitudes of a sporadic E layer detected in CSES radio occultation profiles and the virtual heights of a sporadic E layer obtained by the Wuhan Zuo Ling Tai (ZLT) ionosonde show four different space-time matching criterions. Our results reveal that there is a good agreement between both parameters which is reflected in the significant correlation.

A simple model of ozone–temperature coupling in the tropical lower stratosphere

Observations show strong correlations between large-scale ozone and temperature variations in the tropical lower stratosphere across a wide range of timescales. We quantify this behavior using monthly records of ozone and temperature data from Southern Hemisphere Additional Ozonesonde (SHADOZ) tropical balloon measurements (1998–2016), along with global satellite data from Aura microwave limb sounder and GPS radio occultation over 2004–2018. The observational data demonstrate strong in-phase ozone–temperature coherence spanning sub-seasonal, annual and interannual timescales, and the slope of the temperature–ozone relationship (T / O3) varies as a function of timescale and altitude. We compare the observations to idealized calculations based on the coupled zonal mean thermodynamic and ozone continuity equations, including ozone radiative feedbacks on temperature, where both temperature and ozone respond in a coupled manner to variations in the tropical upwelling Brewer–Dobson circulation. These calculations can approximately explain the observed (T / O3) amplitude and phase relationships, including sensitivity to timescale and altitude, and highlight distinct balances for “fast” variations (periods < 150 d, controlled by transport across background vertical gradients) and “slow” coupling (seasonal and interannual variations, controlled by radiative balances).

Determination of Tropical Belt Widening Using Multiple GNSS Radio Occultation Measurements

In the last decades, Global navigation satellite systems (GNSS) have provided an exceptional opportunity to retrieve atmospheric parameters globally through GNSS radio occultation (GNSS-RO). In this paper, data of 12 GNSS-RO missions from June 2001 to November 2020 with high resolution were used to investigate the possible widening of the tropical belt along with the probable drivers and impacts in both hemispheres. Applying both lapse rate tropopause (LRT) and cold point tropopause (CPT) definitions, the global tropopause height shows increase of approximately 36 m/decade and 60 m/decade, respectively. Moreover, the tropical edge latitude (TEL) estimated based on two tropopause height metrics, in the northern hemisphere (NH) and southern hemisphere (SH), are different from each other. For the first metric, subjective method, the tropical width from GNSS has expansion behavior in NH with ~ 0.41°/decade and a minor expansion in SH with ~ 0.08°/decade. In case of ECMWF Reanalysis v5 (ERA5) there is no significant contraction in both NH and SH. For Atmospheric Infrared Sounder (AIRS), there are expansion behavior in NH with ~ 0.34°/decade and strong contraction in SH with ~ −0.48°/decade. Using the second metric, objective method, the tropical width from GNSS has expansion in NH with ~ 0.13°/decade, and no significant expansion in SH. In case of ERA5, there is no significant signal in NH while SH has a minor contraction. AIRS has an expansion with ~ 0.13°/decade in NH, and strong contraction in SH with ~ −0.37°/decade. The variability of tropopause parameters (temperature and height) is maximum around the TEL locations at both hemispheres. The total column ozone (TCO) shows increasing rates globally, and the rate of increase at the SH is higher than that of the NH. There is a good agreement between the spatial and temporal patterns of TCO variability and the TEL location estimated from GNSS LRT height. Carbon dioxide (CO2), and Methane (CH4), the most important greenhouse gases (GHGs) and the main drivers of global warming, have a global increasing rate and the increasing rate of the NH is similar to that of the SH. The spatial pattern in the NH is located more pole ward than its equivalent at the SH. Both surface temperature and precipitation increase in time and have strong correlation with GNSS LRT height. Both show higher increasing rates at the NH, while the precipitation at the SH has slight decrease and the surface temperature increases. The surface temperature shows a spatial pattern with strong variability, which broadly agrees with the TEL locations. The spatial pattern of precipitation shows northward occurrence. In addition, Standardized Precipitation Evapotranspiration Index (SPEI) has no direct connection with the TEL behavior.

Sporadic E morphology based on COSMIC radio occultation data and its relationship with wind shear theory

The S4max data retrieved from the Constellation Observing System for the Meteorology, Ionosphere, and Climate (COSMIC) radio occultation (RO) measurements during 2007 to 2015 is adopted to investigate the global distribution and seasonal variation of the sporadic E (Es) layers in the present work. The long-term and short-term global Es occurrence maps are presented and the spatial and temporal distributions of Es occurrence rates (ORs) are further confirmed and studied. The International Geomagnetic Reference Field model (IGRF12) is used to calculate the horizontal intensity and inclination of the Earth’s magnetic field. The analysis shows that the Earth’s magnetic field is one of the fundamental reasons for the global distribution of the Es layers. In addition, the Horizontal Wind Field model HWM14 and the IGRF12 model were employed to calculate the vertical ion convergence (VIC) to examine the role of neutral wind shear in the global distribution of the Es ORs. The results reveal that the middle latitude distribution of simulated vertical concentration of Fe+ is similar to that of Es ORs, which indicates that the VIC induced by the neutral wind shear is an important factor in determining the geographical distribution, summer maximum (or winter minimum) and diurnal characteristics of Es ORs in middle latitudes. The new findings mainly include the following two aspects: (1) in summer over mid-latitudes, VIC peaks in the morning and afternoon to evening, which explains the semidiurnal behavior of Es ORs; (2) VIC reaches its minimum value in low-altitude (100 km) areas, which is the reason for the significant decrease in Es ORs in low-altitude areas. The disagreements between the VIC and Es ORs indicate that other processes, such as the meteor influx rate, the ionospheric electric fields and atmospheric tides, should also be considered as they may have an important impact on the variation of Es layers. Graphical Abstract

Variations in Stratospheric Gravity Waves Derived from Temperature Observations of Multi-GNSS Radio Occultation Missions

The spatial–temporal distribution of the global gravity wave (GW) potential energy (Ep) at the lower stratosphere of 20–35 km is studied using the dry temperature profiles from multi- Global Navigation Satellite System (GNSS) radio occultation (RO) missions, including CHAMP, COSMIC, GRACE, and METOP-A/B/C, during the 14 years from 2007 to 2020, based on which the linear trends of the GW Ep and the responses of GW Ep to solar activity, quasi biennial oscillation (QBO), and El Niño-Southern Oscillation (ENSO) are analyzed using the multivariate linear regression (MLR) method. It is found that the signs and the magnitudes of the trends of GW Ep during each month vary at different altitude ranges and over different latitudes. At 25–35 km of the middle and high latitudes, GW Ep values generally show significant negative trends in almost all months, and the values of the negative trends become smaller in the regions closer to the poles. The distribution of the deseasonalized trends in the monthly zonal-mean GW Ep demonstrates that the GW activities are generally declining from 2007 to 2020 over the globe. The responses of GW Ep to solar activity are found to be mostly positive at 20–35 km over the globe, and the comparison between the distribution pattern of the deseasonalized trends in the GW activities and that of the responses of GWs to solar activity indicates that the sharp decline in solar activity from 2015 to 2017 might contribute to the overall attenuation of gravity wave activity during the 14 years. Significant negative responses of GW Ep to QBO are found at 30–35 km over 30° S–25° N, and the negative responses extend to the mid and high latitudes in the southern hemisphere at 20–30 km. The responses of GW Ep to QBO change to be significantly positive at 20–30 km over 15° S–15° N, which demonstrates that the zonal wind field should be the main factor affecting the GW activities at 20–30 km over the tropics. The responses of GW Ep at 20–35 km to ENSO are found to be positive over 15° S–15° N, while at 30–35 km over 15° N–30° N and at 20–35 km near 50° N, significant negative responses of GW Ep to ENSO exist.

Comparison of Seasonal foEs and fbEs Occurrence Rates Derived from Global Digisonde Measurements

A global climatology of sporadic-E occurrence rates (ORs) based on ionosonde measurements is presented for the peak blanketing frequency, fbEs, and the ordinary mode peak frequency of the layer, foEs. ORs are calculated for a variety of sporadic-E frequency thresholds: no lower limit, 3, 5, and 7 MHz. Seasonal rates are calculated from 64 Digisonde sites during the period 2006–2020 using ionograms either manually or automatically scaled with ARTIST-5. Both foEs and fbEs ORs peak in the Northern Hemisphere during the boreal summer, with a decrease by roughly a factor of 2–3 in fbEs rates relative to foEs rates without a lower threshold on the sporadic-E intensity. This ratio of foEs to fbEs OR increases with increasing sporadic-E intensity, up to a factor of 5 for the 7 MHz threshold. An asymmetry is observed with the Southern Hemisphere peaks during the austral summer, with slightly lower rates compared with the Northern Hemisphere during the boreal summer. A drastic decrease in ORs is observed for the higher intensity thresholds, such that the fbEs occurrence rates for 7 MHz are nearly zero during most locations and seasons. These updated occurrence rates can be used for future statistical comparisons with GPS radio occultation-based sporadic-E occurrence rates.