Optimization Algorithm for Determining the Source Surface Radius Based on Parker Solar Probe in situ Measurements from Encounters 1 to 19

Technische Samenvatting: Optimalisatiealgoritme voor het bepalen van de bronoppervlakradius op basis van in situ-metingen van de Parker Solar Probe

Probleemstelling
De Potential Field Source Surface (PFSS) extrapolatie is een standaardmethode voor het schatten van het grootschalige coronaal magnetisch veld vanuit fotosferische magnetogrammen. Een kritische parameter in dit model is de bronoppervlakradius ($R_{ss}$), de sferische grens waar het magnetisch veld wordt aangenomen als radiaal te worden. Hoewel $R_{ss} = 2,5 R_s$ een conventionele referentie is, geven eerdere studies aan dat één enkele waarde niet universeel aan alle observationele beperkingen (bijv. omvang van coronaal gaten, open magnetische flux en magnetische polariteit) kan voldoen. De optimale $R_{ss}$ staat bekend als variabel met de zonneactiviteit, maar het bepalen van een rigoureuze, observationeel geconstante waarde blijft een uitdaging. Dit artikel behandelt het inverse probleem van het bepalen van de optimale $R_{ss}$ door de discrepantie tussen PFSS-extrapolaties en in situ magnetische veldmetingen te minimaliseren, waarbij specifiek gebruik wordt gemaakt van gegevens van de Parker Solar Probe (PSP) tijdens de ontmoetingen (encounters) 1–19.

Methodologie
De auteurs benaderen de bepaling van $R_{ss}$ als een vormoptimalisatieprobleem binnen een free boundary framework.

1. Wiskundige Formulering:

   • Het voorwaartse PFSS-probleem wordt geformuleerd als een tweede-orde elliptisch randwaardeprobleem met gemengde Dirichlet- en Neumann-condities. De auteurs bewijzen de welgesteldheid (existentie en uniciteit) van de voorwaartse oplossing in de Sobolev-ruimte $H^1(\Omega)$.
   • Het inverse probleem wordt gegoten in het minimaliseren van een objectieve functionaal over een toegestane verzameling domeinen. De auteurs bewijzen de compactheid van de toegestane verzameling en de continuïteit van de objectieve functionaal, waarmee zij de existentie van een globaal minimum voor het inverse probleem vaststellen. Deze bewijzen zijn van toepassing op zowel sferische als niet-sferische bronoppervlakgeometrieën.

2. Data en Voorverwerking:

   • Inputs: Synoptische fotosferische magnetogrammen van de Global Oscillation Network Group (GONG) dienen als de binnenste randvoorwaarde.
   • Observaties: Radiaal magnetisch veld ($B_r$) en radiale zonne wind snelheid ($v_r$) van PSP (FIELDS en SWEAP instrumenten) voor de ontmoetingen 1–19.
   • Verwerking: PSP-data ondergaan uitschieterverwijdering (Interquartile Range criterium) en temporele middeling (20-minuten cadans). Gaten in de zonne wind snelheid worden opgevuld via lineaire interpolatie om continue Parker-spiraal backmapping mogelijk te maken.
   • Validatie: Onafhankelijke in situ metingen van de Advanced Composition Explorer (ACE) nabij 1 AU worden gebruikt voor cross-validatie.

3. Optimalisatiealgoritme:

   • Objectieve Functionaal: Het algoritme minimaliseert de gemiddelde kwadratische fout (MSE) tussen de geobserveerde PSP $B_r$ en de door PFSS voorspelde $B_r$. De PFSS-voorspelling wordt afgeleid door de positie van het ruimtevaartuig terug te traceren naar het bronoppervlak langs de Parker-spiralen, het PFSS-veld te samplen, en het radiaal te schalen naar de afstand van het ruimtevaartuig ($r^2 B_r = \text{const}$).
   • Normalisatie: Om encounter-afhankelijke amplitudeverschillen te verwijderen, ondergaan zowel de geobserveerde als de gemodelleerde velden een Z-score normalisatie voordat de MSE wordt berekend.
   • Solver: Een gradiëntgebaseerd iteratief schema (Algoritme 1) wordt toegepast. De gradiënt wordt benaderd met behulp van gecentreerde eindverschillen. Het algoritme iteratert tot convergentie of een maximaal aantal iteraties is bereikt.
   • Validatie: Het algoritme wordt gevalideerd tegen een analytische potentieel veld oplossing en synthetische data met toegevoegde Gaussische ruis, wat het vermogen aantoont om een voorgeschreven $R_{ss}$ te herstellen.

4. Multiobjectieve Analyse:

   • Om trade-offs te begrijpen, introduceren de auteurs twee aanvullende metrieken: de magnetische veldschaalfactor ($\alpha_k$, gerelateerd aan open flux) en de polariteitsvoorspellingsnauwkeurigheid ($P^k_a$).
   • Een Pareto-analyse (Algoritme 2) wordt uitgevoerd om de Pareto-front te identificeren tussen de MSE-objectief en de polariteitsvoorspellingsfout. Dit kwantificeert of de optimalisatie primair wordt gedreven door amplitude-overeenkomst of polariteitsnauwkeurigheid.

Belangrijkste Resultaten

• Optimale $R_{ss}$ Trend: De geoptimaliseerde $R_{ss}$-waarden voor PSP ontmoetingen 1–19 nemen over het algemeen toe van zonneminimum naar de stijgende fase van Cyclus 25. Deze trend wordt bevestigd door onafhankelijke ACE-data.
• Vergelijking met Referentie: Vergeleken met de vaste referentie $R_{ss} = 2,5 R_s$, tonen de geoptimaliseerde oplossingen over het algemeen een betere overeenstemming met geobserveerde magnetische veldamplitudes (het verminderen van de open flux discrepantie), terwijl de polariteitsvoorspellingsnauwkeurigheid behouden blijft of verbetert.
• Zonneminimum vs. Stijgende Fase:
  • Tijdens het zonneminimum (vroege ontmoetingen) wordt de optimalisatie primair gedreven door amplitude-overeenkomst (open flux). Lagere $R_{ss}$-waarden worden gekozen om de gemodelleerde open flux te verhogen om overeen te komen met observaties.
  • Tijdens de stijgende fase worden de Pareto-fronten meer geclusterd, wat aangeeft dat de polariteitsvoorspellingsnauwkeurigheid de dominante metriek wordt. De optimale $R_{ss}$ stijgt om de complexe sectorstructuren en polariteitsomkeringen die worden waargenomen beter te vangen.
• Lokale Afwijkingen: Specifieke ontmoetingen (bijv. 9 en 13) vertonen lokale afwijkingen van de algemene trend. Voor ontmoeting 13, die een complexe Interplanetaire Coronaal Massa Ejectie (ICME) kruiste, koos de optimalisatie een zeer lage $R_{ss}$ ($1,27 R_s$) om complexe polariteitsinversies te vangen, wat de beperkingen van de Parker-spiraal aanname in verstoorde omgevingen benadrukt.

Betekenis en Claims
Het artikel beweert een wiskundig rigoureus fundament te bieden voor het PFSS inverse probleem, waarbij de existentie van een unieke minimizer voor de schatting van de bronoppervlakradius wordt bewezen. Door een specifiek optimalisatiealgoritme te construeren dat gevalideerd is tegen analytische oplossingen en in situ data, toont de studie aan dat $R_{ss}$ geen statische constante is, maar een dynamische parameter die evolueert met de zonneactiviteit.

De auteurs benadrukken dat er niet één enkele "globaal optimale" $R_{ss}$ is; de optimale waarde hangt in plaats daarvan af van de specifieke observationele beperkingen en geëvalueerde metrieken die prioriteit krijgen (bijv. open flux versus polariteit). De studie belicht een verschuiving van de dominante metriek van amplitude-overeenkomst tijdens het zonneminimum naar polariteit-overeenkomst tijdens de stijgende fase. Bovendien erkent het papier het "open flux probleem", waarbij wordt opgemerkt dat hoewel de optimalisatie de overeenstemming met in situ flux verbetert, het de fundamentele discrepantie tussen coronaal gat-gebaseerde fluxschattingen en in situ metingen niet oplost, wat suggereert dat toekomstige verbeteringen mogelijk betere polaire magnetische veldobservaties of meer geavanceerde coronaal modellen vereisen. Het werk legt een kader voor toekomstige vormoptimalisatie van niet-sferische bronoppervlakken naarmate gegevens van missies zoals Solar Orbiter beschikbaar komen.