Predicting Atmospheric Re-entries Using Segmented M/A Calibration of Public TLE Data

Technische Samenvatting: Voorspellen van Atmosferische Re-entries middels Gesegmenteerde M/A-Kalibratie van Publieke TLE-data

Probleemstelling
Het voorspellen van de exacte timing van atmosferische re-entries voor objecten op lage hoogte wordt bemoeilijkt door de intrinsieke ruis en de cadans van publiekelijk beschikbare Two-Line Element (TLE) data. Traditionele long-arc drag fitting divergeert vaak omdat het kortstondige catalogusruis accumuleert, waardoor de werkelijke dynamische evolutie van het object wordt vertroebeld. Bovendien kan de fysieke staat van een herintredend object (attitude, tuimelen, thermische condities en ballistische coëfficiënt) variëren op tijdsschalen van minuten, waardoor gedetailleerde ab-initio fysieke modellering verouderd is voordat deze voltooid kan worden. De uitdaging ligt in het extraheren van een stabiel dynamisch signaal uit ruisige, discrete TLE-observaties zonder dat gedetailleerde kennis vereist is van de geometrie of de aerodynamische staat van het object.

Methodologie: Gesegmenteerde M/A-Kalibratie
Het artikel stelt een empirische, reproduceerbare methode voor die het terminale vervalpad reconstrueert door de effectieve massa-tot-oppervlakteverhouding ($M/A$) te kalibreren over opeenvolgende korte propagatiebogen. De kernprincipes en stappen zijn:

• Gesegmenteerde aanpak: In plaats van een enkele long-arc fit, gebruikt de methode paren van opeenvolgende TLE's om onafhankelijke segmenten te definiëren. De eerdere TLE wordt gepropageerd naar de epoch van de latere TLE.
• Parameterkalibratie: Het effectieve drag-oppervlak ($A$) wordt aangepast zodat de gesimuleerde semi-major-axis afname ($\Delta a_{sim}$) overeenkomt met de geobserveerde afname ($\Delta a_{obs}$) tussen de twee TLE's. Het optimale oppervlak wordt bepaald via lineaire interpolatie, wat een effectieve $M/A$-verhouding oplevert voor dat specifieke interval.
• Ruisonderdrukking: Om de seculaire afname te isoleren van catalogusruis, wordt een minimale scheiding van 12 uur tussen TLE's afgedwongen voor standaardkalibratie. In de terminale fase, waar de daling van de semi-major axis groot genoeg wordt om fluctuaties te domineren, worden segmenten geaccepteerd als de semi-major axis met ten minste 10 km afneemt, zelfs bij kortere tussenpozen.
• Force Model: Propagaties worden uitgevoerd met NASA's GMAT met een vast force model (EGM96 zwaartekracht, NRLMSISE-00 dichtheid, zonnestralingsdruk en aardoriëntatieparameters). De enige aangepaste parameter is het effectieve oppervlak $A$, dat fungeert als een proxy voor het netto drag-gedrag en variaties in geometrie, attitude en thermische staat absorbeert.
• Databron: De methode vertrouwt uitsluitend op publieke data: opeenvolgende TLE's en de uiteindelijke Tracking and Impact Prediction (TIP) berichten van Space-Track.

Belangrijkste Resultaten
De methode werd toegepast op recente natuurlijke re-entries van raketlichamen, puin en niet-manoeuvrerende satellieten, evenals specifieke Starlink-gevallen:

• Nauwkeurigheid: In de terminale fase reproduceert de methode consistent gepubliceerde TIP-tijden met een nauwkeurigheid van ongeveer $\pm 1$ uur. De afwijking ten opzichte van de uiteindelijke TIP's neemt in de laatste 6–12 uur van het baanleven vaak met een factor twee tot vier af vergeleken met eerdere schattingen.
• Identificatie van de laatste baan: De aanpak identificeert betrouwbaar de laatste volledige baan voorafgaand aan de atmosferische intreding, waardoor alle regio's buiten de laatste grondtrack met zekerheid kunnen worden uitgesloten.
• Casestudies:
  • Natuurlijk verval: Voor objecten zoals NORAD 43802 en 61440 convergeerde de methode soepel naar de uiteindelijke TIP.
  • Starlink Variaties: De methode volgde succesvol zowel natuurlijk verval (STARLINK-2171) als thrusting-events (STARLINK-1363). Echter, voor snelle deorbit-manoeuvres (STARLINK-34110) was real-time voorspelling onmogelijk tijdens de thrusting-fase; de gekalibreerde $M/A$-waarden piekten (wat duidt op voortstuwing in plaats van drag), en een nauwkeurige reconstructie was pas retrospectief mogelijk nadat de manoeuvre was voltooid.
• Onafhankelijke Validatie: Verschillende voorspellingen werden bevestigd door onafhankelijke observaties, waaronder een ISS-waarneming van de PROGRESS MS-34 re-entry, getuigenrapporten in Sri Lanka voor ICEYE-X26, en 91 rapporten van de American Meteor Society voor CLARITY-1.

Beperkingen
Het artikel definieert expliciet de operationele reikwijdte en beperkingen van de methode:

• Catalogusbeperkingen: De belangrijkste beperkende factor is de intrinsieke ruis, kwantisatie en publicatiecadans van de TLE-catalogus, niet de atmosferische of zonneflux-modellering.
• Manoeuvres: Gecontroleerde manoeuvres of thrusting-episoden maken kortetermijnvoorspellingen tijdelijk ongeldig omdat de evolutie van de semi-major axis voortstuwing reflecteert in plaats van drag. Voorspellingen worden pas weer betrouwbaar zod{%} het object terugkeert naar een puur ballistisch regime.
• Precisieplafond: Zelfs onder gunstige omstandigheden kan de uiteindelijke vervaltijd doorgaans worden gereconstrueerd met een nauwkeurigheid van ongeveer $\pm 1$ uur. Dit is een fundamentele limiet opgelegd door de informatiehouding van de TLE-data.

Betekenis en Claims
Het artikel claimt dat gesegmenteerde $M/A$-kalibratie een transparant, lichtgewicht en operationeel robuust alternatief biedt voor long-arc drag fitting wanneer alleen publieke trackingdata beschikbaar zijn. De betekenis ligt in:

1. Empirische Robuustheid: Door elk segment te verankeren aan een nieuwe catalogusmeting, voorkomt de methode long-arc divergentie en isoleert het catalogusruis als de dominante onzekerheid.
2. Praktische Bruikbaarheid: Het vereist geen gedetailleerde kennis van de aerodynamische staat of geometrie van het object, waardoor het toepasbaar is op een breed scala aan objecten (puin, raketlichamen, satellieten) waarover dergelijke gegevens onbekend zijn.
3. Reproduceerbaarheid: De methode biedt een reproduceerbaar kader voor eventidentificatie en snelle vervalreconstructie, geschikt voor onafhankelijke verificatie van publiekelijk gerapporteerde re-entries.

De auteur concludeert dat hoewel de methode de intrinsieke beperkingen van de TLE-catalogus niet kan overwinnen, het effectief het maximale dynamische signaal uit publieke data extraheert, wat zorgt voor nauwkeurigheid op uurschaal in de kritieke terminale fase van re-entry.