skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Comparison of a Neutral Density Model With the SET HASDM Density Database
Abstract The EXospheric TEMperatures on a PoLyhedrAl gRid (EXTEMPLAR) method predicts the neutral densities in the thermosphere. The performance of this model has been evaluated through a comparison with the Air Force High Accuracy Satellite Drag Model (HASDM). The Space Environment Technologies (SET) HASDM database that was used for this test spans the 20 years 2000 through 2019, containing densities at 3 hr time intervals at 25 km altitude steps, and a spatial resolution of 10° latitude by 15° longitude. The upgraded EXTEMPLAR that was tested uses the newer Naval Research Laboratory MSIS 2.0 model to convert global exospheric temperature values to neutral density as a function of altitude. The revision also incorporated time delays that varied as a function of location, between the total Poynting flux in the polar regions and the exospheric temperature response. The density values from both models were integrated on spherical shells at altitudes ranging from 200 to 800 km. These sums were compared as a function of time. The results show an excellent agreement at temporal scales ranging from hours to years. The EXTEMPLAR model performs best at altitudes of 400 km and above, where geomagnetic storms produce the largest relative changes in neutral density. In addition to providing an effective method to compare models that have very different spatial resolutions, the use of density totals at various altitudes presents a useful illustration of how the thermosphere behaves at different altitudes, on time scales ranging from hours to complete solar cycles.  more » « less
Award ID(s):
2019465
PAR ID:
10375156
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Space Weather
Volume:
19
Issue:
12
ISSN:
1542-7390
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The Iterative Driver Estimation and Assimilation (IDEA) data assimilation technique was used with the Whole Atmosphere Model (WAM) to improve neutral density specification in the upper thermosphere. Two different neutral density data sources were examined to enhance the capability of simulating the global thermospheric state. The first were accelerometer estimates of neutral density from the Challenging Mini‐Satellite Payload (CHAMP) satellite. The second were neutral density estimates from the Global Ultraviolet Imager (GUVI) limb‐scan airglow observations aboard the Thermosphere Ionosphere Mesosphere Energetics and Dynamics satellite. Due to the intensity of the November 2003 storm, two changes were necessary in WAM. The first was allowing the Kp geomagnetic index to exceed 9 and the second was changing the relationship between Kp and the solar wind parameters used to drive the model. With these changes, results show that IDEA effectively captures the thermospheric neutral density at the CHAMP satellite altitude and follows the time‐dependence through the November 2003 storm period. Furthermore, a cross‐comparison was conducted with the GUVI dayside limb scan measurements. GUVI neutral densities within 270–320 km show the closest agreement with WAM when CHAMP data was assimilated by IDEA. We speculate on the potential for observations from GUVI at 300 km to be used as a data source in the IDEA‐WAM simulations. These simulations demonstrate the utility of the IDEA data assimilation technique with physical models and that using either accelerometer observations or ultraviolet airglow limb measurement during extreme storm periods could be used. 
    more » « less
  2. Abstract The largest obstacle to managing satellites in low Earth orbit (LEO) is accurately forecasting the neutral mass densities that appreciably impact atmospheric drag. Empirical thermospheric models are often used to estimate neutral densities but they struggle to forecast neutral densities during geomagnetic storms when they are highly variable. Physics‐based models are thus increasingly turned to for their ability to describe the dynamical evolution of neutral densities. However, these models require observations to constrain dynamical state variables to be able to forecast mass densities with adequate fidelity. The LEO environment has scarce neutral state observations. Here, we demonstrate, in simulated experiments, a reduction in orbit errors and neutral densities using a physics‐based, data assimilation approach with ionospheric observations. Using a coupled thermosphere‐ionosphere model, the Thermosphere Ionosphere Electrodynamics General Circulation Model, we assimilate Constellation Observing System for Meterology, Ionosphere, and Climate electron density profiles (EDPs) derived from radio occultation (RO) observations. We use the EDPs to directly update neutral states, improving errors for neutral temperature by 70% and neutral winds by 20%. Updated neutral temperature and neutral winds additionally improve helium composition errors by 60% and 40%, respectively. Improved neutral density estimates correspond to a reduction in orbit errors of 1.2 km over 2 days, a 70% reduction over a no‐assimilation control, and a 29 km improvement over 9 days. This study builds on the results of our earlier work to further develop and demonstrate the potential of using a vast and growing RO data source, with a physics‐based model, to overcome our limited number of neutral observations. 
    more » « less
  3. Observations of solar terminator waves (STWs) in thermospheric mass density (TMD) measurements above 500 km reveal STWs as dominant features of the upper thermosphere. While previous investigations have shown that STWs in the middle‐lower thermosphere have amplitudes on the order of 6%–8% of the background TMD in that region, this study shows that STWs exhibit a striking amplification with altitude, producing density perturbations of up to a factor of two near 500 km. The study analyzes STWs in TMD data across altitude, solar cycle, and both solstices, leveraging a direct comparative methodology with Challenging Minisatellite Payload, Gravity Recovery and Climate Experiment Follow‐On, and High‐Accuracy Satellite Drag Model. The first detection of a dawn STW with a magnitude comparable to dusk is made near 500 km, contrasting the asymmetry seen at lower altitudes. The newfound prominence of STWs highlights the need for further research into their generation mechanisms, role in geophysical variability, and broader implications for thermospheric modeling and spacecraft operations. 
    more » « less
  4. Abstract On 3 February 2022, at 18:13 UTC, SpaceX launched and a short time later deployed 49 Starlink satellites at an orbit altitude between 210 and 320 km. The satellites were meant to be further raised to 550 km. However, the deployment took place during the main phase of a moderate geomagnetic storm, and another moderate storm occurred on the next day. The resulting increase in atmospheric drag led to 38 out of the 49 satellites reentering the atmosphere in the following days. In this work, we use both observations and simulations to perform a detailed investigation of the thermospheric conditions during this storm. Observations at higher altitudes, by Swarm‐A (∼438 km, 09/21 Local Time [LT]) and the Gravity Recovery and Climate Experiment Follow‐On (∼505 km, 06/18 LT) missions show that during the main phase of the storms the neutral mass density increased by 110% and 120%, respectively. The storm‐time enhancement extended to middle and low latitudes and was stronger in the northern hemisphere. To further investigate the thermospheric variations, we used six empirical and first‐principle numerical models. We found the models captured the upper and lower thermosphere changes, however, their simulated density enhancements differ by up to 70%. Further, the models showed that at the low orbital altitudes of the Starlink satellites (i.e., 200–300 km) the global averaged storm‐time density enhancement reached up to ∼35%–60%. Although such storm effects are far from the largest, they seem to be responsible for the reentry of the 38 satellites. 
    more » « less
  5. 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