This work estimates global‐mean K_zz using Sounding of the Atmosphere using Broadband Emission Radiometry/Thermosphere‐Ionosphere‐Mesosphere Energetics and Dynamics monthly global‐mean CO₂ profiles and a one‐dimensional transport model. It is then specified as a lower boundary into the Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model (TIE‐GCM). Results first show that global‐mean CO₂ in the mesosphere and lower thermosphere region has annual and semiannual oscillations (AO and SAO) with maxima during solstice seasons along with a primary maximum in boreal summer. Our calculated AO and SAO in global‐mean CO₂ are then modeled by AO and SAO in global‐mean K_zz. It is then shown that our estimated global‐mean K_zz is lower in magnitude than the suggested global‐mean K_zz from Qian et al. (2009) that can model the observed AO and SAO in the ionosphere/thermosphere (IT) region. However, our estimated global‐mean K_zz is similar in magnitude with recent suggestions of global‐mean K_zz in models with explicit gravity wave parameterization. Our work therefore concludes that global‐mean K_zz from global‐mean CO₂ profiles cannot model the observed AO and SAO in the IT region because our estimated global‐mean K_zz may only be representing eddy diffusion due to gravity wave breaking. The difference between our estimated global‐mean K_zz and the global‐mean K_zz from Qian et al. (2009) thus represents diffusion and mixing from other nongravity wave sources not directly accounted for in the TIE‐GCM lower boundary conditions. These other sources may well be the more dominant lower atmospheric forcing behind the AO and SAO in the IT region.

Quantifying the eddy diffusion coefficient profile in the mesosphere and lower thermosphere (MLT) is critical to the constituent density distributions in the upper mesosphere and thermosphere. Previous work by Swenson et al. (2018, https://doi.org/10.1016/j.jastp.2018.05.014) estimated the global mean eddy diffusion (k_zz) values in the upper mesosphere using atomic oxygen (O), derived from Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) hydroxyl (OH). In this study, vertical eddy diffusive transport velocities of O were determined from continuity of mass in the mesopause region (80–97 km), primarily via the HO_x chemistry. Global average constituent climatology from previously deduced SABER ozone (O₃) and atomic hydrogen (H) was applied. Furthermore, we extended the global mean eddy transport velocities to new heights (105 km) in the MLT using the newly available global mean Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) data. The combined method of determining O₃ loss and O density climatology from SCIAMACHY, as well as an improved global mean background atmosphere from SABER, provides new information for eddy diffusion determination in the MLT. Three prominent results to emerge from this study include (i) global mean k_zz profiles between 80 and 105 km derived from MLT constituent climatologies, SABER, and SCIAMACHY global mean O density profiles averaged for approximately one solar cycle, (ii) determination of O eddy diffusion velocities in the MLT consistent between two satellite measurements and the thermosphere‐ionosphere‐mesosphere‐electrodynamics general circulation model, and (iii) resolution of historically large differences between deduced k_zz determined from O versus CO₂ by analysis of SABER and SCIAMACHY measurements.

Global Navigation Satellite System radio occultation signals often show extremely strong levels of scintillation when passing through the ionospheric E region. This is related to sporadic E ( E s )—dense layers of metallic ions that can form in the E region, influencing terrestrial and satellite radio propagation. In our report on the 2007–2014 variation of E region S 4 amplitude fluctuation indices measured by the FORMOSAT‐3/Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) satellite constellation, we find that the spatial and temporal variation of the maximum S 4 index in the E region is proportionate to the occurrence rate of extreme scintillation and by extension, sporadic E . We also find that the monthly median extreme S 4 amplitude fluctuation index in the E region midlatitudes shows a dependence on variation of the El Niño–Southern Oscillation (ENSO) in the troposphere that has not been previously reported. The ENSO‐related variation of the E region median extreme S 4 indices varies closely with the tropopause height, with both parameters lagging the Oceanic Niño Index by roughly 1 to 2 months, while also displaying a similar spectrum of periodicities. This similarity is especially strong in the southern midlatitudes. These results indicate that ENSO signatures can be transmitted to E s formation mechanisms, potentially through modulation of vertically propagating atmospheric tides that alter lower thermospheric wind shears. The end result is the modulation of the interannual variation of extreme E s values by ENSO.

This work uses 17 years of upper mesospheric carbon monoxide (CO) and temperature observations by the microwave limb sounder (MLS) on-board the Aura satellite to present and explain the seasonal and interannual variability of the migrating diurnal tide (DW1) component of upper mesospheric CO. This work then compares these observations to simulations by the specified dynamics – whole atmosphere community climate model with ionosphere/thermosphere extension (SD-WACCM-X). Results show that, for all seasons, MLS CO local-time perturbations peaks above 85 km and has a latitude structure resembling the (1,1) mode in temperature. On the other hand, SD-WACCM-X DW1 also peaks above 85 km and has a latitude structure resembling the (1,1) mode, but it simulates two local maximum of the (1,1) mode between 85 and 92 km. Despite the differences in altitude structure, a tendency analysis and the adiabatic displacement method revealed that, on seasonal and interannual timescales, observed and modeled CO's (1,1) component can be reproduced solely using vertical advection. It was also found that both observed and modeled CO's (1,1) component contains interannual oscillations with periodicities close to that of the quasi-biennial oscillation and the El Niño–Southern Oscillation. From these results, this work concludes that on seasonal and interannual timescales, the observed and modeled (1,1) mode affects the global structure of upper mesospheric CO primarily through vertical advection.

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in characterizing the global 3D distribution and variability of ionospheric electron density (N_e). In this study, we developed an optimal estimation (OE) method to retrieve N_e profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles (ε < 0°). Although both OE and onion-peeling (OP) methods use the Abel weighting function in the N_e inversion, they are significantly different in terms of performance in the lower ionosphere. The new OE results can overcome the large N_e oscillations, sometimes negative values, seen in the OP retrievals in the E-region ionosphere. In the companion paper in this Special Issue, the HmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006–2019), COSMIC-2 (2019–present), and Spire (2019–present) constellations, showing a consistent ionospheric N_e morphology. The unprecedented spatiotemporal sampling of the ionosphere from these constellations now allows a detailed analysis of the frequency–wavenumber spectra for the N_e variability at different heights. In the lower ionosphere (~150 km), we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2, and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4, and SW4. The co-existence of eastward- and westward-propagating wave4 components implies the presence of a stationary wave4 (SPW4), as suggested by other earlier studies. Further improvements to the OE method are proposed, including a tomographic inversion technique that leverages the asymmetric sampling about the tangent point associated with GNSS-LEO links.