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

Technische Zusammenfassung: Einfluss der Intensität geomagnetischer Stürme auf die thermosphärische Dichte und den orbitalen Zerfall von LEO-Satelliten

Problemstellung
Die rasche Expansion von Satellitenkonstellationen im niedrigen Erdorbit (LEO) hat die operative Notwendigkeit einer präzisen Charakterisierung der thermosphärischen Dichte erhöht. Der atmosphärische Widerstand ist die dominierende nicht-gravitative Störung für Satelliten unterhalb von 600 km, und geomagnetische Stürme erhöhen die thermosphärische Dichte signifikant, was den orbitalen Zerfall beschleunigt. Während historische Ereignisse, wie der Verlust von 38 Starlink-Satelliten im Februar 2022 und der vorzeitige Wiedereintritt eines Starlink-Satelliten im Oktober 2024, die schweren operativen Risiken verdeutlichen, verlassen sich bestehende statistische Methoden oft auf die Superposed Epoch Analysis (SEA). Die SEA mittelt die Reaktionen über viele Ereignisse hinweg, was potenziell die scharfen Spitzenreaktionen glättet, die für die unmittelbare Orbitvorhersage und Kollisionsvermeidung entscheidend sind. Darüber hinaus haben empirische thermosphärische Modelle oft Schwierigkeiten, Dichtezunahmen während intensiver Sturmereignisse genau zu reproduzieren.

Methodik
Diese Studie verwendet einen „Per-Event-Peak“-Analyseansatz, um die instantane Spitzenreaktion von 22 geomagnetischen Stürmen der Kategorien G1 bis G5 (maximales $K_p$ von 5 bis 9) zu isolieren, die die Sonnenzyklen 23, 24 und 25 (2003–2024) abdecken. Der Datensatz umfasst den Gannon-G5-Sturm vom Mai 2024.

• Datenquellen: In-situ-Neutralmasse-Dichte und präzise Orbitdaten wurden aus Beschleunigermessungen der GRACE- (2002–2017) und GRACE-FO- (2018–gegenwärtig) Zwillingssatelliten abgeleitet. Geomagnetische Indizes ($K_p$, $Dst$, $SYM-H$, $AU$, $AL$) und Sonnenwindparameter wurden aus NASA/GSFC OMNIWeb bezogen.
• Ereignisauswahl: 22 Ereignisse wurden basierend auf drei Kriterien ausgewählt: (1) maximales $K_p \ge 5$ (G1 oder höher); (2) kontinuierliche Verfügbarkeit von Beschleunigerdaten von 48 Stunden vor dem Storm Sudden Commencement (SSC) bis 48 Stunden nach dem Sturmmaximum; und (3) keine überlappenden vorangehenden Hauptphasen des Sturms. Die Auswahl umfasst jeweils fünf Ereignisse für die Kategorien G1–G4 sowie die zwei verfügbaren G5-Ereignisse (Oktober 2003 und Mai 2024).
• Verarbeitung:
  1. Orbit-Auswahl: Für jedes Ereignis wurde der „Sturm-Orbit“ als der Orbit definiert, der dem Sturmmaximum (Minimum $Dst$) am nächsten liegt, und der „ruhige Referenzorbit“ als der Orbit, der dem SSC am nächsten liegt, welcher nach mindestens 24 aufeinanderfolgenden Stunden mit $K_p \le 3$ auftrat.
  2. Berechnung der Zerfallsrate: Die rohen 10-Sekunden-Höhendaten wurden pro Orbit gemittelt und anschließend mittels eines 24-Stunden-gleitenden Durchschnitts geglättet, um die semidiurnale Modulation zu entfernen. Lineare Fits wurden auf die ruhigen und die Sturmperioden angewendet, um die orbitalen Zerfallsraten ($R_{quiet}$ und $R_{storm}$) zu bestimmen.
  3. Verstärkungsmetrik: Die durch den Sturm verursachte Steigerung des orbitalen Zerfalls wurde als prozentuale Zunahme berechnet: $\frac{R_{storm} - R_{quiet}}{R_{quiet}} \times 100$.
  4. Statistische Modellierung: Ein empirisches Multi-Lineares Regressionsmodell (MLR) wurde entwickelt, um die Zerfallsraten während der Sturmzeit unter Verwendung von drei zugänglichen Prädiktoren vorherzusagen: Zerfallsrate in der Ruhephase, Betriebshöhe des Satelliten und maximales $K_p$.

Wichtigste Ergebnisse

• Monotonische Skalierung mit der Intensität: Die mittlere Steigerung der orbitalen Zerfallsrate nimmt monoton mit der Sturmintensität zu. Die beobachteten Verstärkungen liegen bei: 77,2 % (G1), 98,5 % (G2), 138,1 % (G3), 149,6 % (G4) und 166,5 % (G5).
• Höhe als kritischer Modulator: Die Betriebshöhe des Satelliten fungiert als unabhängiger und signifikanter Faktor. Ereignisse in niedrigen Höhen (z. B. ~350 km) zeigten selbst unter relativ schwacher geomagnetischer Forcierung starke Drag-Reaktionen, die oft die mittlere Verstärkung höherer Intensitätsstürme in größeren Höhen übertrafen.
• Ereignisspezifische Variabilität: Es besteht eine signifikante Streuung innerhalb der Intensitätskategorien. Bemerkenswert ist, dass die beiden G5-Ereignisse (Oktober 2003 und Mai 2024) trotz identischer $K_p$-Werte und ähnlicher Höhen einen Unterschied in der Zerfallsverstärkung von etwa dem 2,6-fachen aufwiesen. Die Autoren führen die höhere Verstärkung im Gannon-Sturm 2024 auf eine geringere Hintergrund-Sonnenaktivität (niedrigerer Nenner der Zerfallsrate in der Ruhephase) und eine längere Dauer zurück, die durch aufeinanderfolgende CMEs (>30 Stunden) getrieben wurde, im Vergleich zum kürzeren, intensiveren Halloween-Sturm von 2003.
• MLR-Modellleistung: Das empirische Modell, $D_s = 1,55D_q - 0,29(h_0 - 450) + 7,39Kp_{max} - 35,60$, erreichte ein Bestimmtheitsmaß ($R^2$) von 0,91 und einen quadratischen Mittelwertfehler (RMSE) von 12,3 m/Tag über alle 22 Ereignisse hinweg. Das Modell erfasst erfolgreich die Abhängigkeit von der Sturmintensität und der Höhe, ohne satellitenspezifische physikalische Parameter (Masse, Querschnitt, Drag-Koeffizient) zu benötigen.

Bedeutung und Ansprüche
Das Paper positioniert diese Arbeit als ergänzende Benchmark zu bestehenden, auf SEA basierenden Analysen. Durch die Konzentration auf die Per-Event-Peak-Reaktionen werden die scharfen Signale bewahrt, die durch Mittelungsverfahren verborgen werden, was eine direktere Validierungsreferenz für empirische thermosphärische Modelle bietet. Die Studie liefert einen Datensatz und ein einfaches, operativ zugängliches MLR-Werkzeug, das es SatBetreibern ermöglicht, die Zerfallsraten während eines Sturms unter Verwendung von Standard-Trackingdaten (Höhe, Zerfallsrate in der Ruhephase) und Weltraumwetter-Indizes ($K_p$) zu schätzen. Die Autoren behaupten, dass diese Ergebnisse die Entwicklung verbesserter Orbitvorhersage- und Kollisionsvermeidungsstrategien für die expandierende LEO-Population unterstützen, insbesondere im Kontext des aktuellen Sonnenzyklus 25.

Die Autoren halten den Umfang der Arbeit hinsichtlich der Verallgemeinerbarkeit moderat und merken an, dass das MLR von zwei spezifischen Satellitenplattformen (GRACE/GRACE-FO) mit ähnlichen physikalischen Eigenschaften abgeleitet wurde. Sie räumen ein, dass zukünftige Arbeiten erforderlich sind, um Masse, Querschnittsfläche und Geschwindigkeit des Satelliten einzubeziehen, um die Anwendbarkeit des Werkzeugs auf die breitere und vielfältigere LEO-Satellitenpopulation zu erweitern.