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.
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- Award ID(s):
- 1651459
- PAR ID:
- 10518805
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Space Weather
- Volume:
- 21
- Issue:
- 9
- ISSN:
- 1542-7390
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract In this study, we present ionospheric observations of field‐aligned currents from AMPERE and the ESA Swarm A satellite, in conjunction with high‐resolution thermospheric density measurements from accelerometers on board Swarm C and GRACE‐FO, for the third and 4 February 2022 geomagnetic storms that led to the loss of 38 Starlink internet satellites. We study the global storm time response of the thermospheric density enhancements, including their decay and latitudinal distribution. We find that the thermospheric density enhances globally in response to high‐latitude energy input from the magnetosphere‐solar wind system and takes at least a full day to recover to pre‐storm density levels. We also find that the greatest density perturbations occur at polar latitudes consistent with the magnetosphere‐ionosphere dayside cusp, and that there appeared to be a saturation of the thermospheric density during the geomagnetic storm on the fourth. Our results highlight the critical importance of high‐latitude ionospheric observations when diagnosing potentially hazardous conditions for low‐Earth‐orbit satellites.
-
more » « less
Abstract The geomagnetic storm on February 3, 2022 caused the loss of 38 Starlink satellites of Space‐X. The Global‐scale Observations of the Limb and Disk (GOLD) observations and Multi‐Scale Atmosphere Geospace Environment (MAGE) model simulations are utilized to investigate the thermospheric composition responses to the Space‐X storm. The percentage difference of the GOLD observed thermospheric O and N2column density ratio (∑O/N2) between the storm time (February 3, Day‐of‐Year [DOY] 34) and quiet time (DOY 32) shows a depletion region in the local noon sector mid‐high latitudes in the southern hemisphere, which corresponds to the east side of GOLD field‐of‐view (FOV). This is different from the classic theory of thermospheric composition disturbance during geomagnetic storms, under which the ∑O/N2depletion is usually generated at local midnight and high latitudes, and thus, appear on the west side of GOLD FOV. MAGE simulations reproduce the observations qualitatively and indicate that the ∑O/N2depletion is formed due to strong upwelling in the local morning caused by strong Joule heating. Interestingly, enhanced equatorward winds appear near local midnight, but also in the local morning sector, which transports ∑O/N2depletion equatorward. The depletion corotates toward the local afternoon and is observed in the GOLD FOV. The equatorward winds in the local morning are due to the ion‐neutral coupling under the conditions of a dominant positive interplanetary magnetic field east‐west component (
B y) during the storm. -
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 Effects of Hurricane Grace in August 2021 are studied in the thermosphere and ionosphere, using data from the COSMIC‐2, ICON, and GOLD satellites. Significant impacts on electron density, thermospheric winds, and temperature are observed after the onset of the hurricane, compared to the pre‐hurricane phase. Comparison of the observations during the hurricane with the ones during a non‐hurricane year clearly provides further evidence for substantial hurricane‐induced thermospheric and ionospheric changes. We reveal an enhancement in electron density during the hurricane's rapid intensification and pronounced changes in thermospheric winds. Additionally, the low‐latitude thermosphere exhibits considerable warming of up to 70 K around 150 km during this period. These changes highlight the long‐range vertical coupling mechanisms between hurricanes and the upper atmosphere, and provide valuable insights into the profound influence of meteorological events on upper atmospheric dynamics, emphasizing the need for further exploration.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1029/2023SW003521
