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

Résumé technique : Impact de l'intensité des tempêtes géomagnétiques sur la densité thermosphérique et la désintégration orbitale des satellites en orbite basse (LEO)

Énoncé du problème
L'expansion rapide des constellations de satellites en orbite terrestre basse (LEO) a accentué le besoin opérationnel d'une caractérisation précise de la densité thermosphérique. La traînée atmosphérique est la perturbation non gravitationnelle dominante pour les satellites situés en dessous de 600 km, et les tempêtes géomagnétiques augmentent considérablement la densité thermosphérique, accélérant ainsi la désintégration orbitale. Alors que des événements historiques, tels que la perte de 38 satellites Starlink en février 2022 et la rentrée prématurée d'un satellite Starlink en octobre 2024, démontrent les risques opérationnels sévères, les méthodes statistiques existantes reposent souvent sur l'analyse de l'époque superposée (SEA). La méthode SEA moyenne les réponses à travers de nombreux événements, ce qui peut lisser les pics de réponse aigus, pourtant critiques pour la prédiction orbitale immédiate et l'évitement de collisions. De plus, les modèles thermosphériques empiriques peinent souvent à reproduire avec précision les augmentations de densité lors d'événements de tempête intenses.

Méthodologie
Cette étude emploie une approche d'analyse par « pic par événement » afin d'isoler la réponse de pointe instantanée de 22 tempêtes géomagnétiques allant de G1 à G5 (valeur maximale de $K_p$ de 5 à 9) couvrant les cycles solaires 23, 24 et 25 (2003–2024). L'ensemble de données inclut la tempête G5 de Gannon de mai 2024.

• Sources de données : Les mesures de masse neutre in situ et les données d'orbite précises ont été dérivées des mesures d'accéléromètre des satellites jumeaux GRACE (2002–2017) et GRACE-FO (2018–présent). Les indices géomagnétiques ($K_p$, $Dst$, $SYM-H$, $AU$, $AL$) et les paramètres du vent solaire ont été obtenus via NASA/GSFC OMNIWeb.
• Sélection des événements : 22 événements ont été sélectionnés selon trois critères : (1) un $K_p$ maximal $\ge 5$ (G1 ou supérieur) ; (2) la disponibilité de données d'accéléromètre continues de 48 heures avant le commencement soudain de la tempête (SSC) jusqu'à 48 heures après le pic de la tempête ; et (3) l'absence de phases principales de tempêtes précédentes chevauchantes. La sélection comprend cinq événements pour chaque catégorie G1–G4 et les deux événements G5 disponibles (octobre 2003 et mai 2024).
• Traitement :
  1. Sélection de l'orbite : Pour chaque événement, l'« orbite de tempête » a été définie comme l'orbite la plus proche du pic de la tempête (minimum $Dst$), et l'« orbite de référence calme » comme l'orbite la plus proche du SSC survenant après au moins 24 heures consécutives de $K_p \le 3$.
  2. Calcul du taux de désintégration : Les données d'altitude brutes toutes les 10 secondes ont été moyennées par orbite, puis lissées à l'aide d'une moyenne mobile de 24 heures pour éliminer la modulation semi-diurne. Des ajustements linéaires ont été appliqués aux périodes calmes et de tempête pour déterminer les taux de désintégration d'altitude ($R_{calme}$ et $R_{tempête}$).
  3. Métrique d'augmentation : L'augmentation de la désintégration orbitale pendant la tempête a été calculée comme l'augmentation en pourcentage : $\frac{R_{tempête} - R_{calme}}{R_{calme}} \times 100$.
  4. Modélisation statistique : Un modèle de régression linéaire multiple (MLR) empirique a été développé pour prédire les taux de désintégration en période de tempête à l'aide de trois prédicteurs accessibles : le taux de désintégration en période calme, l'altitude de fonctionnement du satellite et le $K_p$ maximal.

Résultats clés

• Mise à l'échelle monotone avec l'intensité : Le taux moyen d'augmentation de la désintégration orbitale augmente de manière monotone avec l'intensité de la tempête. Les augmentations observées sont : 77,2 % (G1), 98,5 % (G2), 138,1 % (G3), 149,6 % (G4) et 166,5 % (G5).
• L'altitude comme modulateur critique : L'altitude de fonctionnement du satellite agit comme un facteur indépendant et significatif. Les événements à basse altitude (par exemple, ~350 km) ont présenté de fortes réponses de traînée même sous un forçage géomagnétique relativement faible, dépassant souvent l'augmentation moyenne des tempêtes de plus haute intensité à des altitudes plus élevées.
• Variabilité spécifique aux événements : Une dispersion significative existe au sein des catégories d'intensité. Notamment, les deux événements G5 (octobre 2003 et mai 2024) ont montré une différence d'environ 2,6 fois dans l'augmentation de la désintégration malgré des valeurs de $K_p$ identiques et des altitudes similaires. Les auteurs attribuent l'augmentation plus élevée lors de la tempête Gannon de 2024 à un niveau d'activité solaire de fond plus faible (dénominateur du taux de désintégration en période calme plus bas) et à une durée prolongée pilotée par des CME successives (>30 heures), comparé à la tempête d'Halloween de 2003, plus courte et plus intense.
• Performance du modèle MLR : Le modèle empirique, $D_s = 1,55D_q - 0,29(h_0 - 450) + 7,39Kp_{max} - 35,60$, a atteint un coefficient de détermination ($R^2$) de 0,91 et une erreur quadratique moyenne (RMSE) de 12,3 m/jour sur l'ensemble des 22 événements. Le modèle capture avec succès la dépendance à l'intensité de la tempête et à l'altitude sans nécessiter de paramètres physiques spécifiques au satellite (masse, section efficace, coefficient de traînée).

Signification et revendications
L'article positionne ce travail comme une référence complémentaire aux analyses existantes basées sur la méthode SEA. En se concentrant sur les réponses de pic par événement, il préserve les signaux aigus que les procédures de moyennage masquent, offrant ainsi une référence de validation directe pour les modèles thermosphériques empiriques. L'étude fournit un ensemble de données et un outil MLR simple et opérationnellement accessible qui permet aux opérateurs de satellites d'estimer les taux de désintégration en période de tempête en utilisant uniquement les données de suivi standard (altitude, désintégration en période calme) et les indices de la météo spatiale ($K_p$). Les auteurs affirment que ces résultats soutiennent le développement de stratégies améliorées de prédiction d'orbite et d'évitement de collision pour la population croissante de satellites LEO, particulièrement dans le contexte de l'actuel cycle solaire 25.

Les auteurs maintiennent une portée modeste concernant la généralisabilité du modèle, notant que le MLR a été dérivé de deux plateformes satellitaires spécifiques (GRACE/GRACE-FO) présentant des caractéristiques physiques similaires. Ils reconnaissent que des travaux futurs seront nécessaires pour intégrer la masse, la section efficace et la vitesse du satellite afin d'étendre l'applicabilité de cet outil à la population plus large et plus diversifiée de satellites LEO.