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
Thermospheric Neutral Density Variation During the “SpaceX” Storm: Implications From Physics‐Based Whole Geospace Modeling
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
- Award ID(s):
- 2033843
- PAR ID:
- 10505016
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- Space Weather
- Volume:
- 20
- Issue:
- 12
- ISSN:
- 1542-7390
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Mesoscale high‐latitude electric fields are known to deposit energy into the ionospheric and thermospheric system, yet the energy deposition process is not fully understood. We conduct a case study to quantify the energy deposition from mesoscale high‐latitude electric fields to the thermosphere. For the investigation, we obtain the high‐latitude electric field with mesoscale variabilities from Poker Flat Incoherent Scatter Radar measurements during a moderate geomagnetic storm, providing the driver for the Global Ionosphere and Thermosphere Model (GITM) via the High‐latitude Input for Mesoscale Electrodynamics framework. The HIME‐GITM simulation is compared with GITM simulations driven by the large‐scale electric field from the Weimer model. Our modeling results indicate that the mesoscale electric field modifies the thermospheric energy budget primarily through enhancing the Joule heating. Specifically, in the local high‐latitude region of interest, the mesoscale electric field enhances the Joule heating by up to five times. The resulting neutral temperature enhancement can reach up to 50 K above 200 km altitude. Significant increase in the neutral density above 250 km altitude and in the neutral wind speed are found in the local region as well, lagging a few minutes after the Joule heating enhancement. We demonstrate that the energy deposited by the mesoscale electric field transfers primarily to the gravitational potential energy in the thermosphere.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 (By) during the storm.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 During the second recovery phase of the 13–14 March 2022 storm, intense high‐latitude neutral mass density spikes are detected by satellites at ∼500 km in both hemispheres. These density spikes, accurately modeled by the Global Ionospheric Thermosphere Model (GITM), are identified as high‐latitude neutral mass density anomalies (HDAs). The GITM simulation indicates that these HDAs, which extends over the polar region, are important features in high‐latitude neutral density. Furthermore, GITM reveals that these HDAs are manifestations of transpolar traveling atmospheric disturbances triggered on the dawn side. Moreover, GITM also reveals significant interhemispheric asymmetries (IHAs) in the magnitude, propagation speed, and propagation direction of HDAs, which are linked to the IHAs in the distribution and magnitude of Joule heating deposited as well as the thermospheric background conditions. This study presents a dynamic perspective on the IHA of storm‐time high‐latitude neutral density variations that is particularly helpful to the proper interpretation of satellite observations.more » « less
An official website of the United States government

