Atomic oxygen (O) in the mesosphere and lower thermosphere (MLT) results from a balance between production via photo‐dissociation in the lower thermosphere and chemical loss by recombination in the upper mesosphere. The transport of O downward from the lower thermosphere into the mesosphere is preferentially driven by the eddy diffusion process that results from dissipating gravity waves and instabilities. The motivation here is to probe the intra‐annual variability of the eddy diffusion coefficient (k
NRLMSIS® 2.0 is an empirical atmospheric model that extends from the ground to the exobase and describes the average observed behavior of temperature, eight species densities, and mass density via a parametric analytic formulation. The model inputs are location, day of year, time of day, solar activity, and geomagnetic activity. NRLMSIS 2.0 is a major, reformulated upgrade of the previous version, NRLMSISE‐00. The model now couples thermospheric species densities to the entire column, via an effective mass profile that transitions each species from the fully mixed region below ~70 km altitude to the diffusively separated region above ~200 km. Other changes include the extension of atomic oxygen down to 50 km and the use of geopotential height as the internal vertical coordinate. We assimilated extensive new lower and middle atmosphere temperature, O, and H data, along with global average thermospheric mass density derived from satellite orbits, and we validated the model against independent samples of these data. In the mesosphere and below, residual biases and standard deviations are considerably lower than NRLMSISE‐00. The new model is warmer in the upper troposphere and cooler in the stratosphere and mesosphere. In the thermosphere, N2and O densities are lower in NRLMSIS 2.0; otherwise, the NRLMSISE‐00 thermosphere is largely retained. Future advances in thermospheric specification will likely require new in situ mass spectrometer measurements, new techniques for species density measurement between 100 and 200 km, and the reconciliation of systematic biases among thermospheric temperature and composition data sets, including biases attributable to long‐term changes.
more » « less- Award ID(s):
- 1829138
- NSF-PAR ID:
- 10374826
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Earth and Space Science
- Volume:
- 8
- Issue:
- 3
- ISSN:
- 2333-5084
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract zz ) and eddy velocity in the MLT based on the climatology of the region, initially accomplished by Garcia and Solomon (1985,https://doi.org/10.1029/JD090iD02p03850 ). In the current study, the intra‐annual cycle was divided into 26 two‐week periods for each of three zones: the northern hemisphere (NH), southern hemisphere (SH), and equatorial (EQ). Both 16 years of SABER (2002–2018) and 10 years of SCIAMACHY (2002–2012) O density measurements, along with NRLMSIS®2.0 were used for calculation of atomic oxygen eddy diffusion velocities and fluxes. Our prominent findings include a dominant annual oscillation below 87 km in the NH and SH zones, with a factor of 3–4 variation between winter and summer at 83 km, and a dominant semiannual oscillation at all altitudes in the EQ zone. The measured global average kzz at 96 km lacks the intra‐annual variability of upper atmosphere density data deduced by Qian et al. (2009,https://doi.org/10.1029/2008JA013643 ). The very large seasonal (and hemispherical) variations in kzz and O densities are important to separate and isolate in satellite analysis and to incorporate in MLT models. -
more » « less
Abstract The ultraviolet‐imaging spectrograph that comprises Global‐scale Observations of the Limb and Disk (GOLD) mission in geostationary orbit at 47.5°W longitude has taken full disk images at high cadence throughout the deep solar minimum period of 2019–2020. Synoptic (i.e., concurrent and spatially unified and resolved) observations of thermospheric temperature and composition at ∼150 km altitude are made for the first time, allowing GOLD to disambiguate temporal and spatial variations. Here we analyze the daytime effective temperature and column integrated O and N2density ratio (ΣO/N2) data simultaneously observed by GOLD over 120°W–20°E longitude and 60°S–60°N latitude from 13 October 2019 to 12 October 2020. Daily zonal mean values are calculated for each latitude and compared with NRLMSIS 2.0 and simulations from the Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension (WACCM‐X). On average, the GOLD observations show higher temperatures than Mass Spectrometer Incoherent Scatter radar (MSIS) and WACCM‐X by ∼20–60 K (5%–10%) and 80–120 K (12%–18%), respectively. The ΣO/N2ratios observed by GOLD are larger than the MSIS results by ∼0.4 (40%) but smaller than the WACCM‐X simulations by ∼0.3 (30%). The observed and modeled results are correlated at most latitudes (
r = 0.4–0.8), and GOLD, MSIS, and WACCM‐X all display a similar seasonal variation and change with latitude. WACCM‐X simulates a larger annual variation in ΣO/N2, suggesting that the thermospheric circulation is overestimated and atmospheric waves and turbulence transport are not properly represented in the model. -
more » « less
Abstract The space weather research community relies heavily on thermospheric density data to understand long‐term thermospheric variability, construct assimilative, empirical, and semiempirical global atmospheric models and validate model performance. One of the challenges in resolving accurate thermospheric density data sets from satellite orbital drag measurements is modeling appropriate physical aerodynamic drag force coefficients. The drag coefficient may change throughout the thermosphere due to model dependencies on composition and altitude. As such, existing drag coefficient model errors and corresponding errors in orbit‐derived density data sets and models may be altitude and solar cycle dependent with greater errors at higher altitudes around 500 km near the oxygen‐to‐helium transition region. In this paper, inter‐satellite observed‐to‐modeled density comparisons at ∼500 km are evaluated to constrain drag coefficient modeling assumptions. Observed densities are derived from accelerometer data for the Gravity Recovery and Climate Experiment (GRACE) satellites and Two‐Line Element data for a set of compact satellites, while the NRLMSISE‐00 atmospheric model is used to obtain modeled densities and composition information. Density consistency results indicate that drag coefficient models with incomplete energy and momentum accommodation produce the most consistent densities, while the standard diffuse modeling approach may not be appropriate at these altitudes. Models with momentum accommodation between 0.5 and 0.9 and energy accommodation between 0.83 and 0.96 may be most appropriate at upper thermospheric altitudes. Modeling drag coefficients with diffuse gas‐surface interactions for the GRACE satellites could lead to errors in derived density of ∼25% and in‐track satellite orbit prediction uncertainty during solar maximum conditions on the order of kilometers.
-
more » « less
Abstract The Starlink satellites launched on 3 February 2022 were lost before they fully arrived in their designated orbits. The loss was attributed to two moderate geomagnetic storms that occurred consecutively on 3–4 February. We investigate the thermospheric neutral mass density variation during these storms with the Multiscale Atmosphere‐Geospace Environment (MAGE) model, a first‐principles, fully coupled geospace model. Simulated neutral density enhancements are validated by Swarm satellite measurements at the altitude of 400–500 km. Comparison with standalone TIEGCM and empirical NRLMSIS 2.0 and DTM‐2013 models suggests better performance by MAGE in predicting the maximum density enhancement and resolving the gradual recovery process. Along the Starlink satellite orbit in the middle thermosphere (∼200 km altitude), MAGE predicts up to 150% density enhancement near the second storm peak while standalone TIEGCM, NRLMSIS 2.0, and DTM‐2013 suggest only ∼50% increase. MAGE also suggests altitudinal, longitudinal, and latitudinal variability of storm‐time percentage density enhancement due to height dependent Joule heating deposition per unit mass, thermospheric circulation changes, and traveling atmospheric disturbances. This study demonstrates that a moderate storm can cause substantial density enhancement in the middle thermosphere. Thermospheric mass density strongly depends on the strength, timing, and location of high‐latitude energy input, which cannot be fully reproduced with empirical models. A physics‐based, fully coupled geospace model that can accurately resolve the high‐latitude energy input and its variability is critical to modeling the dynamic response of thermospheric neutral density during storm time.
-
more » « less
Abstract We examine characteristics of the seasonal variation of thermospheric composition using column number density ratio ∑
O/N 2 observed by the NASA Global Observations of Limb and Disk (GOLD) mission from low‐mid to mid‐high latitudes. We also use ∑O/N 2 derived from the Global Ultraviolet Imager (GUVI) limb measurements onboard the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite and estimated by the NRLMSISE‐00 empirical model to aid our investigation. We found that the∑O/N 2 seasonal variation is hemispherically asymmetric: in the southern hemisphere, it exhibits the well‐known annual and semiannual pattern, with highs near the equinoxes, and primary and secondary lows near the solstices. In the northern hemisphere, it is dominated by an annual variation, with a minor semiannual component with the highs shifting toward the wintertime. We also found that the durations of the December and June solstice seasons in terms of∑O/N 2 are highly variable with longitude. Our hypothesis is that ion‐neutral collisional heating in the equatorial ionization anomaly region, ion drag, and auroral Joule heating play substantial roles in this longitudinal dependency. Finally, the rate of change in ∑O/N 2 from one solstice season to the other is dependent on latitude, with more dramatic changes at higher latitudes.
