Impact of geomagnetic storm intensity on thermospheric density and orbital decay of LEO satellites: a G1–G5 per-event analysis

Technische Samenvatting: Impact van de Intensiteit van Geomagnetische Stormen op de Thermosferische Dichtheid en Orbitale Verval van LEO-satellieten

Probleemstelling
De snelle expansie van constellaties in een lage aardbaan (LEO) heeft de operationele behoefte aan nauwkeurige karakterisering van de thermosferische dichtheid vergroot. Atmosferische weerstand is de dominante niet-gravitatieve perturbatie voor satellieten onder de 600 km, en geomagnetische stormen verhogen de thermosferische dichtheid aanzienlijk, waardoor het orbitale verval wordt versneld. Terwijl historische gebeurtenissen, zoals het verlies van 38 Starlink-satellieten in februari 2022 en de voortijdige terugkeer in de atmosfeer van een Starlink-satelliet in oktober 2024, de ernstige operationele risico's aantonen, vertrouwen bestaande statistische methoden vaak op superposed epoch analysis (SEA). SEA middelt responsen over vele gebeurtenissen, wat potentieel de scherpe piekresponsen afvlakt die cruciaal zijn voor directe baanvoorspelling en botsingsvermijding. Bovelijk worstelen empirische thermosferische modellen er vaak mee om dichtheidsverhogingen tijdens intense stormgebeurtenissen nauwkeurig te reproduceren.

Methodologie
Deze studie hanteert een "per-event peak" analyse-aanpak om de instantane piekrespons van 22 geomagnetische stormen van G1 tot G5 (maximale $K_p$ van 5 tot 9) te isoleren, verspreid over de zonnecycli 23, 24 en 25 (2003–2024). De dataset omvat de mei 2024 Gannon G5-storm.

• Databronnen: In-situ neutrale massadichtheid en precieze baangegevens werden afgeleid van versnellingsmeter-metingen van de GRACE (2002–2017) en GRACE-FO (2018–heden) tweelingsatellieten. Geomagnetische indices ($K_p$, $Dst$, $SYM-H$, $AU$, $AL$) en zonnewindparameters werden verkregen van NASA/GSFC OMNIWeb.
• Selectie van Gebeurtenissen: 22 gebeurtenissen werden geselecteerd op basis van drie criteria: (1) maximale $K_p \ge 5$ (G1 of hoger); (2) continue beschikbaarheid van versnellingsmetergegevens van 48 uur vóór de Storm Sudden Commencement (SSC) tot 48 uur na de stormpiek; en (3) geen overlappende voorafgaande stormfases. De selectie bevat telkens vijf gebeurtenissen voor de categorieën G1–G4 en de twee beschikbare G5-gebeurtenissen (oktober 2003 en mei 2024).
• Verwerking:
  1. Baanselectie: Voor elke gebeurtenis werd de "stormbaan" gedefinieerd als de baan die het dichtst bij de stormpiek (minimum $Dst$) ligt, en de "rustige referentiebaan" als de baan die het dichtst bij de SSC ligt die plaatsvond na ten minste 24 opeenvolgende uren van $K_p \le 3$.
  2. Berekening van de Vervalrate: Ruwe 10-seconden hoogtegegevens werden per baan gemiddeld, en vervolgens gladgestreken met een 24-uurs voortschrijdend gemiddelde om semi-diurnale modulatie te verwijderen. Lineaire fits werden toegepast op zowel de rustige als de stormperiodes om de hoogte-vervalrates ($R_{quiet}$ en $R_{storm}$) te bepalen.
  3. Verhogingsmetriek: De verhoging van de orbitale vervalrate tijdens de storm werd berekend als de procentuele toename: $\frac{R_{storm} - R_{quiet}}{R_{quiet}} \times 100$.
  4. Statistische Modellering: Een empirisch multi-lineair regressiemodel (MLR) werd ontwikkeld om storm-tijd vervalrates te voorspellen met behulp van drie toegankelijke voorspellers: de vervalrate tijdens rustige periodes, de operationele hoogte van de satelliet en de maximale $K_p$.

Belangrijkste Resultaten

• Monotone Schaling met Intensiteit: De gemiddelde orbitale vervalrate-verhoging neemt monotoon toe met de stormintensiteit. De geobserveerde verhogingen zijn: 77,2% (G1), 98,5% (G2), 138,1% (G3), 149,6% (G4) en 166,5% (G5).
• Hoogte als Kritieke Modulator: De operationele hoogte van de satelliet fungeert als een onafhankelijke en significante factor. Gebeurtenissen op lage hoogte (bijv. ~350 km) vertoonden sterke weerstandsresponsen zelfs onder relatief zwakke geomagnetische forcering, waarbij ze vaak de gemiddelde verhoging van hoog-intensiteit stormen op grotere hoogtes overtroffen.
• Gebeurtenis-specifieke Variabiliteit: Er bestaat aanzienlijke spreiding binnen de intensiteitscategorieën. Opvallend genoeg vertoonden de twee G5-gebeurtenissen (oktober 2003 en mei 2024) een ~2,6-voudig verschil in vervalverhoging, ondanks identieke $K_p$-waarden en vergelijkbare hoogtes. De auteurs schrijven de hogere verhoging in de 2024 Gannon-storm toe aan een lager achtergrondniveau van de zonneactiviteit (een lagere noemer voor de rustige-tijd vervalrate) en een langdurige duur gedreven door opeenvolgende CMEs (>30 uur), vergeleken met de kortere, intensere Halloween-storm van 2003.
• Prestaties van het MLR-model: Het empirische model, $D_s = 1.55D_q - 0.29(h_0 - 450) + 7.39Kp_{max} - 35.60$, bereikte een determinatiecoëfficiënt ($R^2$) van 0,91 en een wortelmiddelkwadraatfout (RMSE) van 12,3 m/dag over alle 22 gebeurtenissen. Het model vangt de afhankelijkheid van stormintensiteit en hoogte succesvol op zonder dat satelliet-specifieke fysieke parameters (massa, dwarsdoorsnede, weerstandscoëfficiënt) nodig zijn.

Betekenis en Claims
Het artikel positioneert dit werk als een complementaire benchmark voor bestaande SEA-gebaseerde analyses. Door te focussen op per-event piekresponsen, behoudt het de scherpe signalen die door gemiddelde procedures worden afgevlakt, wat een directere validatie-referentie biedt voor empirische thermosferische modellen. De studie levert een dataset en een eenvoudig, operationeel toegankelijk MLR-instrument waarmee satellietoperators de vervalrates tijdens stormen kunnen schatten met enkel standaard trackinggegevens (hoogte, rustige-tijd verval) en ruimteweerindices ($K_p$). De auteurs stellen dat deze resultaten de ontwikkeling van verbeterde baanvoorspellings- en botsingsvermijdingsstrategieën ondersteunen voor de expanderende LEO-populatie, met name in de context van de huidige Zonnecyclus 25.

De auteurs hanteren een bescheiden reikwijdte met betrekking tot de generaliseerbaarheid van het model, waarbij zij opmerken dat de MLR is afgeleid van twee specifieke satellietplatforms (GRACE/GRACE-FO) met vergelijkbare fysieke kenmerken. Zij erkennen dat toekomstig werk vereist is om massa, dwarsdoorsnede en snelheid te integreren om de toepasbaarheid van het instrument uit te breiden naar de bredere en meer diverse LEO-satellietpopulatie.