PFSS Model for Solar Coronal Magnetic Fields

• PFSS is a current-free model that reconstructs the Sun’s global coronal magnetic field from photospheric magnetograms by enforcing a radial source surface boundary.
• Its solution uses Laplace’s equation with spherical-harmonic expansion to efficiently derive magnetic topology and support solar wind mapping.
• The method underpins applications in coronal hole diagnostics, solar wind source mapping, and CME kinematics, despite limitations in capturing non-potential dynamics.

The Potential Field Source Surface (PFSS) model is the foundational analytic framework for reconstructing the global coronal magnetic field of the Sun from line-of-sight photospheric magnetograms. PFSS models are built on the assumption that in the low-β corona below a so-called “source surface” (typically at $R_{ss}=2.5R_\odot$), the coronal field is current-free and can be expressed as the gradient of a scalar potential. This approach constrains the open and closed magnetic topology and is a critical component of operational and research pipelines for solar wind source mapping, coronal structure analysis, and predictive space weather modeling.

1. Mathematical Formulation and Boundary Conditions

In the PFSS model, the coronal magnetic field $\mathbf{B}$ between the photosphere ($r=R_\odot$) and a spherical source surface ($r=R_{ss}$) is assumed to satisfy: $$\nabla \times \mathbf{B} = 0 \implies \mathbf{B} = -\nabla \Phi$$ with the scalar potential $\Phi$ satisfying Laplace’s equation: $\nabla^2 \Phi(r,\theta,\phi) = 0$ Boundary conditions are specified as:

• Photospheric boundary ($r=R_\odot$): The observed radial field is imposed: $$B_r(R_\odot, \theta, \phi) = B_{phot}(\theta, \phi)$$, usually from synoptic magnetograms (e.g., GONG, HMI, ADAPT-GONG) (Ervin et al., 2024, Jaffarove et al., 26 Nov 2025).
• Source-surface boundary ($r=R_{ss}$): The field is forced to be purely radial: $\mathbf{B}$0, interpreted as the onset of the high-β heliospheric regime (Amenomori et al., 2013, Benavitz et al., 2024).

The general solution is expressed via a spherical-harmonic expansion: $\mathbf{B}$1 The coefficients are determined from the two boundary conditions.

2. Source Surface Height $\mathbf{B}$2: Physical Interpretation and Optimization

$\mathbf{B}$3 acts as the sole free parameter in the PFSS formalism, demarcating the interface where closed coronal loops are forced to open into the solar wind. The canonical $\mathbf{B}$4 has historical precedent and offers a first-order match to open magnetic flux measured at 1 AU and to heliospheric field observations (Amenomori et al., 2013, Lin et al., 2023). However, both historical empirical (Shoda et al., 7 Oct 2025, Arden et al., 2016) and modern data-driven approaches (Benavitz et al., 2024, Huang et al., 2024) show that $\mathbf{B}$5 should be adjusted dynamically:

• Near solar minimum: Optimal values range from 1.3–1.7 $\mathbf{B}$6; the coronal field opens closer to the Sun.
• Near solar maximum: Values extend up to 3.0 $\mathbf{B}$7, with high-lying loops persisting further outward (Benavitz et al., 2024, Huang et al., 2024).

Empirical relations link $\mathbf{B}$8 to averaged unsigned photospheric field and the global dipolarity (multipole decomposition) (Shoda et al., 7 Oct 2025): $\mathbf{B}$9 at maximum, and

$r=R_\odot$0

at minimum, with $r=R_\odot$1 the normalized dipole content.

3. Computational Solution and Model Implementation

The PFSS solution is implemented by projecting the photospheric boundary field onto spherical harmonics, enforcing the radial boundary at $r=R_\odot$2, and reconstructing $r=R_\odot$3 and hence $r=R_\odot$4 everywhere in $r=R_\odot$5 (Stansby et al., 2022, Jaffarove et al., 26 Nov 2025, Ervin et al., 2024). Standard libraries include pfsspy and the IDL SolarSoft PFSS package. For high-resolution synoptic maps (e.g., 360×180 grid), truncation degrees $r=R_\odot$6–$r=R_\odot$7 are common, with finer grids required for high-fidelity coronal loop tracing and sharp polarity inversion line recovery (Heinemann et al., 16 Jan 2026, Benavitz et al., 2024). Numerical finite-difference solvers complement spherical-harmonic approaches for irregular domain mapping or high-latitude fields.

4. Applications: Magnetic Connectivity, Solar Wind, and Topology

PFSS extrapolations are central to several research and operational applications:

• Solar wind source mapping: PFSS, often coupled with ballistic or Parker-spiral backmapping, links in-situ spacecraft data to coronal and photospheric source regions (e.g., coronal holes, streamers) (Jaffarove et al., 26 Nov 2025, Ervin et al., 2024, Badman et al., 2019).
• Coronal hole and streamer diagnostics: PFSS-derived open field footprints are compared to EUV observations for coronal hole identification, although the Jaccard index evidences poor spatial correspondence in many cases (Heinemann et al., 16 Jan 2026).
• Event forecasting: The model is used as the magnetic module in solar wind speed predictors (e.g., WSA, CNN-based predictors) (Lin et al., 2023).
• Kinematic modeling of CMEs: Parametric studies show that CME deflections and impacts are sensitive to the choice of $r=R_\odot$8 and the magnetogram input (Ledvina et al., 2023).
• Particle propagation: In studies of cosmic-ray shadows and solar gamma rays, the PFSS coronal field acts as a testbed for high-energy particle tracing and secondary production modeling (Amenomori et al., 2013, Li et al., 2020).

5. Model Limitations and Critical Comparisons

Analytical Strengths

• Simplicity and computational efficiency: The linear, static, current-free formulation allows for rapid, analytic solutions, enabling batch processing for time-dependent studies and real-time forecasting (Amenomori et al., 2013).
• Direct data assimilation: PFSS can ingest arbitrary synoptic magnetograms, preserving observed large-scale flux distribution.

Physical and Observational Limitations

• Neglect of currents and plasma forces: PFSS cannot represent non-potentiality (currents, free magnetic energy), dynamic structures (CMEs, reconnection), or pressure-driven phenomena (Antonio et al., 20 Mar 2026, Rice et al., 23 Mar 2026).
• Open flux deficit: Classic PFSS models systematically underestimate open solar magnetic flux by 30–45% at 1 AU (Amenomori et al., 2013, Huang et al., 2024, Rice et al., 23 Mar 2026), particularly at solar minimum.
• Rigid global source surface: A fixed, spherical $r=R_\odot$9 cannot recover spatially or temporally varying footpoint topologies, leading to frequent mismatches with observed coronal hole boundaries and high-latitude field expansions (Heinemann et al., 16 Jan 2026, Ma et al., 24 Oct 2025).
• Limited accuracy in slow-wind/ecliptic regimes: Backmapping uncertainty for photospheric source identification can reach tens of degrees in the ecliptic, undermining rigorous compositional or origin studies in the slow solar wind (Ma et al., 24 Oct 2025).

Comparative and Extended Frameworks

• CSSS and HCCSSS models: These introduce currents and current sheets to recover more open flux, streamer-topology, and helmet streamer structures (Amenomori et al., 2013, Arden et al., 2016). The CSSS model, for example, achieves a better solar cycle correlation with the cosmic-ray shadow than PFSS (Amenomori et al., 2013).
• Outflow Fields/Magnetofrictional extensions: Incorporate effects of coronal wind and weak coronal currents, narrowing the model–observation flux gap and better representing streamer inclination (Rice et al., 23 Mar 2026).
• Data-constrained or multi-constraint models: Hybrid schemes can incorporate coronal loop imaging and 3D stereoscopy into PFSS frameworks, preserving efficiency while improving topological realism (Antonio et al., 20 Mar 2026).

6. Validation, Benchmarking, and Best Practices

Rigorous benchmark tests use single-harmonic analytic solutions to quantify numerical errors in field solution and footpoint tracing (typically <1% and <0.5°, respectively, for well-resolved cases) (Stansby et al., 2022). For operational applications:

• Tuning $r=R_{ss}$0: Best-fit $r=R_{ss}$1 should be empirically optimized by matching modeled open flux to in-situ spacecraft measurements or by minimizing the misalignment angle between model and eclipse-observed streamer directions (Benavitz et al., 2024, Shoda et al., 7 Oct 2025, Huang et al., 2024).
• Empirical parameterization: $r=R_{ss}$2 adjustment as a function of global field measures (e.g., mean unsigned Br, dipolarity) is recommended for time-varying operational pipelines (Shoda et al., 7 Oct 2025).

7. Physical Interpretation and Research Impact

PFSS remains the backbone of coronal topology research, solar-wind mapping, and event-driven heliophysics. The fundamentally potential, current-free description provides global context, but comparison with eclipse coronagraphy, in-situ flux, and dynamic events continually demonstrates both the strengths and the regime-dependent limitations of the approach (Benavitz et al., 2024, Huang et al., 2024, Shoda et al., 7 Oct 2025). The model's role as a baseline for quantifying electric currents and non-potentiality through mismatch analysis is well established, motivating ongoing research into hybrid models, data-driven boundary optimization, and wind-inclusive extensions (Rice et al., 23 Mar 2026).

In practical terms, best practice is to employ PFSS (preferably with a solar-cycle–dependent $r=R_{ss}$3 and ensemble modeling) for high-latitude, coronal-hole, and fast-wind studies (Ma et al., 24 Oct 2025). For low-latitude, streamer-belt, or slow-wind analyses, mapping uncertainties and missing current-driven structure require ensemble, multi-model comparisons or more sophisticated MHD/nonlinear-force-free approaches (Ma et al., 24 Oct 2025, Amenomori et al., 2013).