skip to main content


Title: NRLMSIS 2.0: A Whole‐Atmosphere Empirical Model of Temperature and Neutral Species Densities
Abstract

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
PAR ID:
10374826
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
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
  1. 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
  2. Abstract

    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 (kzz) 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 kzzat 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 kzzand O densities are important to separate and isolate in satellite analysis and to incorporate in MLT models.

     
    more » « less
  3. Abstract

    The Mass Spectrometer and Incoherent Scatter radar (MSIS) model family has been developed and improved since the early 1970's. The most recent version of MSIS is the Naval Research Laboratory (NRL) MSIS 2.0 empirical atmospheric model. NRLMSIS 2.0 provides species density, mass density, and temperature estimates as function of location and space weather conditions. MSIS models have long been a popular choice of thermosphere model in the research and operations community alike, but—like many models—does not provide uncertainty estimates. In this work, we develop an exospheric temperature model based in machine learning that can be used with NRLMSIS 2.0 to calibrate it relative to high‐fidelity satellite density estimates directly through the exospheric temperature parameter. Instead of providing point estimates, our model (called MSIS‐UQ) outputs a distribution which is assessed using a metric called the calibration error score. We show that MSIS‐UQ debiases NRLMSIS 2.0 resulting in reduced differences between model and satellite density of 25% and is 11% closer to satellite density than the Space Force's High Accuracy Satellite Drag Model. We also show the model's uncertainty estimation capabilities by generating altitude profiles for species density, mass density, and temperature. This explicitly demonstrates how exospheric temperature probabilities affect density and temperature profiles within NRLMSIS 2.0. Another study displays improved post‐storm overcooling capabilities relative to NRLMSIS 2.0 alone, enhancing the phenomena that it can capture.

     
    more » « less
  4. 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
  5. Abstract

    Previous work suggested that the peak response time of the mass densities of atomic oxygen (O) and molecular nitrogen (N2) in the thermosphere had more than a 1‐day difference relative to the peak of the 27‐day periodic variation of solar extreme ultraviolet (EUV) flux. In this study, we used the Thermosphere Ionosphere Electrodynamic General Circulation Model (TIEGCM) to explore the physical mechanisms responsible for the different peak response times of the daytime thermospheric neutral species. It was found that the peak response time of O or N2mass density corresponds to the time of equilibrium between the contributions from the barometric effect and the change in its abundance. The peak response time of O is shorter than that of thermospheric temperature Tn, due to a dynamic change in the circulation that acts to cancel out the contribution from the barometric process prior to the peak of Tn. On the contrary, the change of N2abundance contributes further to a decrease of N2mass density on a constant pressure surface when the thermosphere is expanding. The change of chemical loss leads to a longer peak response time of N2abundance than that due to barometric motion. Therefore, an equilibrium is reached after the barometric effect turns from expansion (contraction) to contraction (expansion), so that the peak response time of N2is longer than that of Tn. Moreover, the meridional circulation in the thermosphere modulates the latitudinal dependence of the peak response time of thermospheric neutral species.

     
    more » « less