🔭 astrophysics

Global Coronal Equilibria with Solar Wind Outflow II -- Optimizing the Outflow Model

This paper introduces and optimizes a new "Outflow Fields" model for the Sun's coronal magnetic field that significantly reduces discrepancies with in-situ measurements and eclipse observations compared to traditional PFSS extrapolations, with the resulting methodology and code made available to the community via the Python package "outflowpy."

Original authors: Oliver Rice, Anthony Yeates

Published 2026-03-24

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Oliver Rice, Anthony Yeates

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the Sun as a giant, glowing lighthouse. It doesn't just shine light; it constantly blows a powerful wind made of charged particles (the Solar Wind) that sweeps across our entire solar system, reaching Earth and beyond. To understand how this wind behaves, scientists need to map the invisible "lanes" or magnetic field lines that guide it, much like how a river's current follows the shape of its banks.

For decades, scientists have used a standard map-making tool called PFSS (Potential Field Source Surface) to draw these lanes. Think of PFSS as a very simple, rigid ruler. It assumes the magnetic lanes are straight and smooth, and it forces them to become perfectly straight at a specific, arbitrary height (like saying, "At exactly 2.5 miles up, all rivers must flow straight down").

The Problem with the Old Ruler
The trouble is, the Sun isn't simple. The solar wind pushes and stretches these magnetic lanes, making them curve and twist in ways the old ruler can't predict. Because of this, the old maps often get two things wrong:

The "Open Flux" Problem: They underestimate how much magnetic "stuff" is actually escaping into space. It's like trying to measure the flow of a river with a ruler that only measures the width of the stream at the bank, missing the water rushing in the middle.
The Shape Problem: When we look at the Sun during a total eclipse (when the Moon blocks the bright face, revealing the glowing halo), the real magnetic loops look different from the straight lines the old ruler predicts.

The New Solution: "Outflow Fields"
In this paper, the authors introduce a new, smarter tool called Outflow Fields. Instead of a rigid ruler, imagine a flexible, stretchy rubber band that naturally bends when you pull on it.

This new model adds a "wind factor" to the math. It acknowledges that the solar wind is constantly pulling on the magnetic field, stretching it out.

The Analogy: Imagine holding a piece of string (the magnetic field) in a gentle breeze (the solar wind). The old model (PFSS) assumes the string stays straight until you reach a specific height, then it snaps straight. The new model (Outflow) realizes the breeze is already bending the string as soon as you let go, and it calculates exactly how much the wind is pulling.

How They Tested It
The authors didn't just guess; they played a game of "match the photo."

The Eclipse Photos: They took 12 high-quality photos of solar eclipses from 2006 to 2024. These photos show the actual shape of the Sun's magnetic "hair" (streamers) as seen from Earth.
The Evolutionary Algorithm: They used a computer program that acts like a digital evolution. It tried thousands of different "wind speeds" and "stretching factors" automatically. It kept the versions that looked most like the eclipse photos and threw away the ones that looked weird.
The Result: They found the "Goldilocks" settings. The new model's magnetic lanes matched the real eclipse photos much better than the old model.

Why This Matters

Better Weather Forecasts: Space weather (solar storms) can knock out satellites and power grids. To predict these storms, we need to know exactly how the Sun's magnetic field is shaped. The new model gives a much more accurate starting point for these predictions.
Less Guesswork: The old model required scientists to guess a "magic height" where the field becomes straight. The new model figures this out naturally based on physics, making it more reliable.
The Trade-off: Interestingly, the authors found that if you force the model to perfectly match the amount of magnetic wind escaping (the "Open Flux"), the shape of the lanes looks a bit fake. So, they chose to prioritize the shape (matching the eclipse photos) over the exact number, because a realistic shape is more important for understanding how the Sun actually works.

The Takeaway
The authors have released their new tool as a free software package called 'outflowpy'. It's like upgrading from a paper map to a GPS that accounts for traffic and wind. It helps scientists see the Sun's magnetic atmosphere more clearly, leading to better predictions of space weather that affects our daily lives on Earth.

1. Problem Statement

The paper addresses the limitations of the Potential Field Source Surface (PFSS) model, the standard tool for global coronal magnetic field extrapolation. While PFSS is computationally efficient, it suffers from three critical flaws:

The Open Flux Problem: PFSS significantly underestimates the total Open Solar Flux (OSF) compared to in-situ spacecraft measurements (often by ~45%).
Arbitrary Source Surface Height: PFSS results are highly sensitive to the chosen source surface height ( $r_{ss}$ ), typically set arbitrarily at $2.5 R_\odot$ . Changing this height drastically alters the magnetic topology to match OSF data, often resulting in unrealistic field structures.
Streamer Shape Discrepancy: PFSS models fail to accurately reproduce the shapes of coronal streamers observed during solar eclipses, often extending closed-field regions unrealistically high.

The authors aim to quantify and resolve these issues by optimizing their previously introduced "Outflow Fields" model, which incorporates solar wind effects into the magnetic equilibrium.

2. Methodology

Theoretical Framework

The Outflow model seeks equilibrium solutions to the magneto-frictional (MF) equation, including a term for solar wind outflow, without requiring expensive time-stepping.

Equation: Instead of the potential field solution $\mathbf{B} = \nabla\Phi$ , the Outflow field is defined as $\mathbf{B} = f(r)\nabla\Phi$ .
Equilibrium Condition: The field satisfies $(\nabla \times \mathbf{B}) \times \mathbf{B} / (\nu_0 |\mathbf{B}|^2) + v(r)\mathbf{e}_r \times \mathbf{B} = 0$ , where $\nu_0$ is the magnetofrictional constant and $v(r)$ is the solar wind speed profile.
Key Parameter: The model relies on a dimensionless speed $V(r) = R_\odot \nu_0 v(r)$ , which represents the ratio of outflow strength to frictional relaxation.

Data Processing

Input: Lower boundary magnetic field data from HMI (SDO) and MDI (SOHO) instruments (2000–2022).
Preprocessing: Data is downsampled, smoothed using discrete Legendre Polynomial decomposition to remove noise and artifacts, and flux-balanced to ensure zero net flux.
Solar Wind Profile: The authors move beyond the explicit approximation used in their previous work (Paper I) to solve the implicit isothermal Parker solar wind solution exactly. This allows for non-zero wind speeds near the solar surface, accommodating higher sound speeds.

Optimization Strategy

To determine the optimal parameters ( $\nu_0$ and $v(r)$ ), the authors employ two distinct validation approaches:

OSF Comparison: Comparing model predictions against in-situ OSF measurements from the OMNI database (Wind/ACE spacecraft), corrected for switchbacks.
Topological Optimization (Eclipse Data):
- Data Source: 12 total solar eclipse photographs (2006–2024) processed by Emmanoulidis and Druckmüller.
- Feature Extraction: Using edge detection (Canny algorithm) and filtering, the authors identify magnetic field line angles relative to the radial direction up to $2.5 R_\odot$ .
- Algorithm: A Covariance Matrix Adaptation Evolution Strategy (CMA-ES) is used to optimize the outflow profile $V(r)$ . The objective function minimizes the squared error between the model's field line angles and the observed eclipse angles.
- Synthetic Imaging: A new algorithm generates synthetic white-light eclipse images by tracing field lines and weighting emissions based on Thomson scattering, field strength, and footpoint properties.

3. Key Contributions

Quantification of PFSS Limitations: The paper rigorously demonstrates that PFSS models cannot simultaneously match OSF measurements and observed streamer topologies. Reducing $r_{ss}$ to match OSF creates unrealistic field shapes, while keeping $r_{ss}=2.5 R_\odot$ yields poor OSF agreement.
Robustness of Outflow Fields: The authors show that Outflow Fields are largely independent of the source surface height ( $r_{ss}$ ) above $\approx 2.5 R_\odot$ . The solar wind outflow naturally forces field lines to become radial, stabilizing the OSF regardless of the arbitrary $r_{ss}$ choice.
Optimization via Evolutionary Algorithms: The paper introduces a novel method to empirically derive the optimal solar wind speed profile by matching eclipse field line angles, rather than just fitting scalar OSF values.
Open Source Tool (outflowpy): The authors released a Python package designed to be a drop-in replacement for the existing pfsspy package, facilitating community adoption.

4. Key Results

OSF Discrepancy Reduction:
- PFSS: Underestimates OSF by ~45% (at $r_{ss}=2.5 R_\odot$ ) compared to in-situ data.
- Optimized Outflow: Reduces the discrepancy to ~24%.
- Sensitivity: Increasing $r_{ss}$ to $5.0 R_\odot$ worsens PFSS performance (74% underestimation) but has a much smaller negative impact on Outflow fields (34% underestimation).
Topological Accuracy:
- Outflow fields optimized for topology match observed eclipse field line angles significantly better than PFSS.
- Error Metrics: For the 2017 eclipse, the standard deviation of error for PFSS was 7–21°, whereas the optimized Outflow model reduced this to **1–4°** (depending on $r_{ss}$ ).
- The optimized model correctly captures the non-radial nature of field lines at $2.5 R_\odot$ , a feature PFSS fails to reproduce.
Solar Wind Profile Characteristics:
- The optimized profiles generally show low speeds near the surface, rising to dimensionless speeds of 2.5–9 units by $2.5 R_\odot$ .
- There is no strong correlation between the optimized speed and the solar cycle phase, suggesting that local magnetic topology drives the required outflow speed more than global cycle activity.
Trade-off: The paper highlights a fundamental trade-off: tuning the model to match in-situ OSF exactly requires unrealistically high wind speeds that force field lines to be purely radial (contradicting eclipse photos). Conversely, tuning for eclipse topology yields a more realistic field structure but leaves a residual OSF underestimation (~24%).

5. Significance and Conclusion

This paper establishes Outflow Fields as a superior alternative to PFSS for global coronal modeling. By incorporating a physically motivated solar wind term and optimizing parameters against observational topology (eclipse data), the model achieves:

Realistic Streamer Shapes: Better agreement with white-light observations.
Parameter Stability: Reduced sensitivity to the arbitrary choice of source surface height.
Improved OSF Prediction: A significant reduction in the "Open Flux Problem," though a residual discrepancy remains.

The authors conclude that the remaining OSF discrepancy likely stems from unmodeled dynamic processes (e.g., flux ropes, CMEs, helicity condensation) or calibration issues in the lower boundary magnetogram data (HMI/MDI). They assert that for most applications requiring a global magnetic field approximation, matching the field line topology is more physically accurate than forcing a match to OSF values, and the outflowpy package provides a robust tool for this purpose.