2D Coronal Seismology Technique

• 2D Coronal Seismology is a diagnostic technique that maps the coronal magnetic field and plasma properties by analyzing the propagation of MHD waves.
• It employs advanced wave tracking, density mapping, and inversion algorithms to derive spatially resolved parameters from high-cadence imaging and spectropolarimetric data.
• The technique overcomes traditional 1D limitations and has been validated against simulations and complementary methods for robust field strength and direction mapping.

A two-dimensional (2D) coronal seismology technique enables the spatially resolved mapping of coronal plasma parameters—most notably the magnetic field—across extended regions by analyzing the propagation of magnetohydrodynamic (MHD) waves observed in imaging or spectropolarimetric data. This methodology overcomes the limitations of classical one-dimensional or event-based coronal seismology by exploiting the spatial ubiquity of propagating MHD waves, thus facilitating the construction of coronal “magnetograms” based on wave diagnostics. It serves as a key tool for the global measurement of magnetic field strength and orientation, as well as related plasma properties, by leveraging observations from instruments such as the Coronal Multi-channel Polarimeter (CoMP), UCoMP, SDO/AIA, and complementary extreme ultraviolet (EUV) imagers (Yang et al., 2020, Jess et al., 2016, Yang et al., 15 Jan 2026, Anfinogentov et al., 2019).

1. Theoretical Foundations of 2D Coronal Seismology

The core of 2D coronal seismology is the direct inversion of measured MHD wave properties into plasma parameters. For transverse (kink-mode) oscillations, the thin-tube, low-β approximation yields the local kink speed as:

$$c_k = \sqrt{\frac{\rho_i v_{A,i}^2 + \rho_e v_{A,e}^2}{\rho_i + \rho_e}}$$

where $v_{A} = B/\sqrt{\mu_0 \rho}$, $B$ is the magnetic field strength, and $\rho_{i,e}$ represent internal and external densities. In the unresolved ensemble limit relevant for CoMP/UCoMP imaging, the measured phase speed approaches $$v_{\rm ph} \approx B/\sqrt{\mu_0 \langle \rho \rangle}$$, with $\langle \rho \rangle$ as the local, emissivity-weighted mean density.

For slow magneto-acoustic tube waves, as in sunspot fan regions, the phase speed is set by

$c_t^2 = \frac{c_s^2 v_A^2}{c_s^2 + v_A^2}$

where $c_s$ is the sound speed derived from DEM-based temperature mapping (Jess et al., 2016). These analytic relationships are the basis for mapping measured wave phase speeds and local density into $B$.

2. Observational and Analytical Methodologies

2D seismology is enabled by high-cadence, multi-pixel imaging or spectropolarimetric data. The methodologies can be broadly summarized as:

• Wave Tracking: Motion magnification algorithms (e.g., Dual-Tree Complex Wavelet Transform, DT-CWT) enhance low-amplitude transverse oscillations (including sub-pixel motions), allowing Bayesian or sinusoidal fitting of oscillatory displacements in time-distance slices. For CoMP/UCoMP, Dopplergram time series are cross-correlated spatially to extract local phase speed $v_{\rm ph}(x,y)$ and propagation direction (Yang et al., 2020, Anfinogentov et al., 2019).
• Density Mapping: For lines such as Fe XIII, the intensity ratio $R(x,y) = I_{10798}/I_{10747}$ is inverted via CHIANTI-based lookup curves, accounting for photoexcitation, to map electron density $n_e(x,y)$. EUV-based DEM inversions similarly provide $n_e(x,y)$ and temperature (Jess et al., 2016, Yang et al., 15 Jan 2026).
• Seismological Inversion: The field strength is extracted algebraically,

$$B(x,y) = v_{\rm ph}(x,y) \sqrt{ \mu_0 \langle \rho \rangle(x,y) }$$

where $\langle \rho \rangle = 1.2 m_p n_e$ (factor 1.2 for abundance of He) (Yang et al., 2020, Yang et al., 15 Jan 2026). For slow-mode diagnostics, slow and Alfvén speeds, and ultimately $B$, are computed via the tube mode relationships.

• Field Direction Diagnostics: The wave-propagation direction from coherence mapping is combined with the linear polarization azimuth derived from Stokes $Q,U$ (with Van Vleck correction) to resolve the POS field vector orientation (Yang et al., 2020, Yang et al., 15 Jan 2026).

3. Key Implementations and Parameter Inversion

Distinct 2D coronal seismology approaches have been implemented for both spectropolarimetric and EUV imaging data:

• CoMP/UCoMP 2D Seismology: Yang et al. (2020)(Yang et al., 2020) and subsequent validation (Yang et al., 15 Jan 2026) construct field strength and direction maps over global coronal regions by mapping phase speed and density pixel-by-pixel. LOS integration and spatial resolution are addressed via emissivity-weighted inversion and by quantifying the optimal path length for reliable phase lag fitting through the dimensionless parameter $\alpha = npt \cdot dx \cdot f / v_{\rm ph}$, with an operational regime $\alpha \gtrsim 0.5$ for minimal error (Yang et al., 15 Jan 2026).
• Sunspot Wave-Based Mapping: Jess et al. (2016) (Jess et al., 2016) use propagating slow-mode waves in sunspot fan regions with 2D DEM analysis and pixel-by-pixel phase speed measurements, providing high-resolution, rapidly updated $B(x,y)$ maps. Loop inclination from NLFFF extrapolation or stereoscopy is critical for deprojection.
• Decayless Kink Oscillation Mapping: Tracking ubiquitous low-amplitude kink oscillations (including via motion magnification (Anfinogentov et al., 2016, Anfinogentov et al., 2021)) allows determination of kink and Alfvén speeds in active regions, even in quiet periods (Anfinogentov et al., 2019).

Parameter inference commonly relies on Bayesian methods (e.g., MCMC) for robust uncertainty quantification, with typical accuracies for local field strength of $10$–$20\%$ and for direction of $\sim 5^\circ$ under optimal parameter regimes (Yang et al., 2020, Yang et al., 15 Jan 2026, Anfinogentov et al., 2019).

4. Results, Validation, and Accuracy

2D coronal-seismology–derived field strength and plasma maps exhibit strong consistency with complementary diagnostics (e.g., Zeeman, Hanle, shock standoff, and radio methods) and with forward-modeled MHD simulations:

• Typical measured values in the low corona ($h = 1.05$–$1.35\,R_\odot$): $n_e = 10^{7.5}–10^{8.5}\,{\rm cm}^{-3}$, $v_{\rm ph} = 300–700\,{\rm km\,s}^{-1}$, $B_{\rm POS} = 1–5$ G (Yang et al., 2020).
• Radial power-law decline in $B(r)$, with $\alpha_{\rm QS} \approx 2.9$, $\alpha_{\rm AR} \approx 4.1$ for quiet Sun and active regions, respectively (Yang et al., 2020).
• Forward modeling in volumetric MURaM simulations demonstrates RMSE errors of $\sim 22\%$ in $B_{\rm POS}$, negligible bias, and $\sigma_{\rm angle} \lesssim 5^\circ$ in direction (Yang et al., 15 Jan 2026).
• For sunspot fans, field strength falls from $\sim 32$ G above the umbra to $\sim 1$ G at the edge over $7\,{\rm Mm}$, exceeding the spatial resolution and cadence of prior techniques (Jess et al., 2016).
• For decayless kink oscillations, local kink speeds and inferred Alfvén speeds are mapped along entire loop systems, with typical uncertainties of $10$–$20\%$ for $c_{A0}$ (Anfinogentov et al., 2019).

5. Algorithmic and Data Analysis Advances

Modern 2D coronal seismology incorporates several advanced methodologies for signal extraction, noise mitigation, and inversion:

• Motion Magnification (DT-CWT-based): Enables detection and measurement of sub-pixel oscillations, acting linearly for $A_\mathrm{in} \lesssim 2.5$ pixel and $A_\mathrm{out} \lesssim 10$ pixel, period-independent amplification for $6 \leq P \leq w$ frames (Anfinogentov et al., 2016, Anfinogentov et al., 2021). Advised to verify linear regime and filter windowing to avoid artifacts.
• Empirical Mode Decomposition (EMD): 2D EMD and ensemble (EEMD) variants decompose image sequences into intrinsic mode functions and mitigate mode-mixing or spurious detections, supporting significance assessment against power-law noise spectra (Anfinogentov et al., 2021).
• Spectral and Wavelet Analyses: Multidimensional Fourier/wavelet transforms isolate spatial and temporal oscillations; significance thresholds are set via $\chi^2$-confidence against modeled noise backgrounds (Anfinogentov et al., 2021).
• Bayesian and Forward Modeling: MCMC sampling and synthetic forward-modeling (e.g., FoMo) provide robust uncertainty propagation and validation against realistic topologies (Yang et al., 15 Jan 2026, Anfinogentov et al., 2021).
• LOS Integration and Bias Quantification: Utilize emissivity-weighted diagnostics to minimize ambiguity from overlapping structures along the line of sight, with systematic uncertainties from LOS integration of $\lesssim 12\%$ in $B$ (Yang et al., 2020).

6. Strengths, Limitations, and Operational Best Practices

Strengths:

• Delivers spatially continuous 2D maps of magnetic field strength and direction in the corona, overcoming the isolated-event and single-loop constraints of earlier seismology (Yang et al., 2020, Jess et al., 2016).
• Applicable over active regions, quiet Sun, and sunspot fans, with potential for real-time or high-cadence mapping (down to $\sim$1 min) (Jess et al., 2016).
• Validated by MHD simulations and cross-checked with direct measurements and independent remote sensing (Yang et al., 15 Jan 2026).

Limitations:

• Moderate spatial (e.g., $\sim$9″ for CoMP; $\sim$1.2″ for AIA) and temporal resolution; finite pixel and cadence must satisfy Nyquist limits to avoid spatial or temporal aliasing (Anfinogentov et al., 2021).
• Density diagnostics may be poorly constrained in regions with insufficient signal (e.g., coronal holes) (Yang et al., 2020).
• Systematic uncertainties arise from LOS integration, density inhomogeneities, inclination angle errors, and non-idealized loop geometry (Jess et al., 2016, Yang et al., 2020, Anfinogentov et al., 2019).
• Ambiguities in polarization azimuth (90° Van Vleck effect) and 180° wave-direction require careful multi-diagnostic resolution (Yang et al., 2020, Yang et al., 15 Jan 2026).

Best Practices:

• Optimize the dimensionless control parameter $\alpha$ in phase lag measurements; $\alpha \gtrsim 0.5$ is required to resolve phase lags while minimizing path-averaging bias (Yang et al., 15 Jan 2026).
• Calibrate filtering parameters (frequency, window width) and magnification factors to align with expected wave periods and amplitudes (Anfinogentov et al., 2016).
• Employ robust error propagation and, where possible, validate with forward modeling using synthetic observables.

7. Future Prospects and Directions

Forthcoming instrumentation and methodological advances will enable substantial progress:

• UCoMP and DKIST will extend spatial and spectral capabilities, allowing denser spatial coverage, access to additional coronal lines, and direct LOS field (Zeeman) measurements (Yang et al., 2020).
• Joint analysis combining 2D POS seismology, LOS field measurements, and advanced forward/MHD modeling will enable reconstruction of full 3D coronal vector magnetograms (Yang et al., 2020).
• Further automation and integration of motion magnification, EMD, and Bayesian approaches will facilitate routine and unbiased production of high-fidelity coronal B-field and plasma parameter maps.
• Quantitative assessment of inversion biases using realistic, forward-modeled, multi-structure scenarios remains an active area, guiding the operational regimes and error budgets of 2D seismology (Yang et al., 15 Jan 2026).

2D coronal seismology has thus matured into a robust diagnostic framework at the intersection of remote-sensing, MHD theory, and advanced data analytics, with an increasingly central role in solar coronal physics, space weather diagnosis, and the empirical validation of coronal heating and structuring models (Yang et al., 2020, Anfinogentov et al., 2016, Yang et al., 15 Jan 2026, Anfinogentov et al., 2021, Jess et al., 2016, Anfinogentov et al., 2019).