Coronal Multi-channel Polarimeter (CoMP)

• CoMP is a ground-based imaging spectropolarimeter that systematically probes the solar corona's magnetic fields, plasma structures, and MHD wave dynamics using infrared diagnostic lines.
• It employs a 20 cm Lyot coronagraph, tunable birefringent filters, and full-Stokes polarimetry with high cadence (30 s) to capture detailed spectral and polarimetric images.
• CoMP's advanced diagnostics facilitate coronal seismology, turbulence analysis, and space weather forecasting by integrating observational data with numerical models and inversion techniques.

The Coronal Multi-channel Polarimeter (CoMP) is a ground-based, imaging spectropolarimeter and coronagraph designed to systematically probe the magnetic field, wave dynamics, and plasma structure of the solar corona. Deployed at the Mauna Loa Solar Observatory, CoMP combines high-cadence polarimetric and spectroscopic imaging in infrared forbidden lines of Fe xiii to deliver unprecedented diagnostics of coronal magnetism, MHD wave energy, and eruptive phenomena in the region 1.05–1.4 R⊙. Its full-Stokes vector capability and multi-channel architecture have transformed the empirical mapping and theoretical modeling of the coronal magnetic field, enabling advances in coronal seismology, turbulence, and space weather forecasting.

1. Instrument Architecture and Core Observational Capabilities

CoMP operates as a 20 cm Lyot-type coronagraph coupled to a narrowband, tunable birefringent filter and a high-sensitivity polarimetric imaging system. Incoming solar light is filtered to select forbidden Fe xiii emission lines at 1074.7 nm and 1079.8 nm, and optionally He i 1083.0 nm. This light is passed through a modulation system (rotating waveplates or liquid-crystal variable retarders) and a polarization analyzer, allowing sequential measurement of the full Stokes vector S = [I, Q, U, V]^T across each spectral channel. Standard operation involves three to five wavelength steps per line (Δλ ≈ 0.13 nm spectral resolution), spatial sampling of 4.4–4.5″ per pixel, and a temporal cadence as fast as 30 s for full-field imaging (Gibson, 2015, Tian et al., 2013).

The key instrument metrics are:

Parameter	Value	Reference
Aperture	20 cm Lyot coronagraph	(Gibson, 2015)
Spectral Lines	Fe xiii 1074.7 nm, 1079.8 nm; He i 1083.0 nm	(Gibson, 2015)
Spectral Resolution	0.13 nm (typical), up to 0.01 nm	(Bak-Steslicka et al., 2013)
Spatial Sampling	4.4–4.5″/pixel (~3–3.2 Mm at Sun)	(Tian et al., 2013)
Field of View	1.05–1.40 R⊙ (extensible to 1.5 R⊙)	(Bak-Steslicka et al., 2013)
Temporal Cadence	30 s (full Stokes sequence per channel)	(Gibson, 2015)
Polarimetric Sens.	Q, U: ~3×10⁻⁴ I after integration	(Kramar et al., 2015)

Calibration involves dark and flat-field corrections, Mueller-matrix demodulation for Stokes retrieval, wavelength registration, and correction for instrumental polarization and seeing-induced cross-talk (French et al., 2019, Liu et al., 2014).

2. Spectropolarimetric Diagnostics of the Coronal Magnetic Field

CoMP’s polarimetric measurements directly access the vector components of the coronal magnetic field through atomic alignment and magneto-optical effects (Gibson, 2015, Kramar et al., 2015, Judge et al., 2013). Hanle-effect physics governs the linear polarization (Q, U) of the forbidden Fe xiii lines, encoding the plane-of-sky field orientation. The fundamental diagnostic relations are:

• Linear polarization degree: $P_L = \sqrt{Q^2 + U^2}/I$
• Azimuth angle (POS field): $\chi = (1/2)\arctan(U/Q)$
• Hanle effect sensitivity to field inclination, subject to Van Vleck (θ_VV=54.74°) 90° ambiguities
• Circular polarization (V): weak Zeeman effect, providing line-of-sight field, but at 10⁻⁴ I levels requiring long integration (Gibson, 2015)

Emission coefficients are synthesized by integrating the formal Stokes transport equation along each line of sight:

$\frac{dS}{ds} = -K\cdot S + \varepsilon,$

where $K$ is the propagation matrix and $\varepsilon$ is the emissivity vector with contributions from radiative excitation (Q/U: Hanle), density, temperature, and the local magnetic field (Gibson, 2015, Judge et al., 2013).

Vector magnetic tomography combines CoMP’s Stokes I, Q, U with EUV-derived density/temperature maps to solve for the 3D coronal field via inversion algorithms regularized by smoothness and divergence-free constraints. This procedure leverages the fact that the observed linear polarization is a LOS integral of Hanle-modified emissivities:

$$P_\text{model,i}(B) = \int_{\text{LOS}} \varepsilon(I, Q, U; B(r))\,d\ell,$$

with inversion of the misfit metric $$\chi^2 = \sum_{i} \|P_{\text{obs},i} - P_\text{model,i}(B)\|^2/\sigma_i^2$$ under additional physical constraints (Kramar et al., 2015).

Direct comparison of observed polarization signatures to forward-modeled MHD or force-free extrapolations is routine, including within established frameworks such as the FORWARD IDL package (Gibson, 2015).

3. Coronal Seismology, Wave Kinematics, and Magnetoseismology

CoMP uniquely enables routine measurement of propagating and standing transverse MHD waves (kink/Alfvénic modes) through high-cadence Doppler and polarimetric imaging (Threlfall et al., 2013, Liu et al., 2014, Bak-Steslicka et al., 2013, Yang et al., 2020, Duckenfield et al., 7 Mar 2025). Its methodology encompasses:

• Phase-speed measurement by building time–distance diagrams along field-guided paths and performing cross-correlations or 2D FFTs to extract coherent, outward-propagating wave packets (Yang et al., 2020, Yang et al., 15 Jan 2026).
• Combined density mapping using the Fe xiii 1079.8/1074.7 nm intensity ratio, analyzed via CHIANTI-based calibration, yielding 2D electron density maps.
• Determination of the magnetic field strength using the coronal kink speed in the thin-tube, low-β limit:

$$B = v_\mathrm{ph} \sqrt{\mu_0 \langle \rho \rangle}$$

where $v_\mathrm{ph}$ is the measured phase speed and $\langle \rho \rangle$ is the local emissivity-weighted density; mapping is performed pixelwise to produce 2D coronal “magnetograms” (Yang et al., 2020, Yang et al., 15 Jan 2026).

• Recovery of the field direction both via polarization azimuth (Hanle) and via wave-propagation coherence, resolving the plane-of-sky vector up to Van Vleck-induced 90° ambiguities.

Statistical seismology of hundreds of loops identifies relationships between damping rate, loop length, and density contrast (e.g., equilibrium parameter ξ versus L; longer loops display weaker damping due to declining ζ = ρ_i/ρ_e) (Tiwari et al., 2021).

Recent advances have enabled single-viewpoint inference of kink-mode polarization states—critical for wave energy flux estimation—by combining time-series analysis of the plane-of-sky displacement and LOS Doppler velocity at loop apices (Duckenfield et al., 7 Mar 2025). Periodic line width enhancements in phase with kink oscillations further trace mode coupling and unresolved turbulent flows.

4. Turbulence, Wavefront Dislocations, and Energy Cascade

Comprehensive CoMP wave observations reveal the presence and onsets of Alfvénic turbulence in the corona (Liu et al., 2014). High-frequency Doppler wave power (periods <3 min) preferentially accumulates at the apex of sufficiently long loops (L/λ ≳3), consistent with the emergence of nonlinear interactions and the stochastic cascade characteristic of turbulence. This is quantified via the “ratio difference” (RD) of FFT power between apex and footpoints across spectral bands, and correlated with line width enhancements as signatures of unresolved velocity structures (Liu et al., 2014).

Analysis of Doppler time–distance diagrams has demonstrated robust detection of wavefront dislocations—phase singularities arising from interference between kink (m=1) and sausage (m=0) modes of comparable but distinct frequencies. Such topological structures, identified via phase monodromy integrals, demand at least two co-propagating modes and a loop cross-section below CoMP’s spatial resolution (Ariste et al., 2015). The distribution and statistics of dislocations critically constrain the modal content of coronal waves, supporting the view that Alfvénic and compressible modes are generically excited and interact nonlinearly.

5. Magnetic Topology, Cavities, and CME Precursors

CoMP’s polarimetric coverage has enabled routine identification of complex coronal magnetic topologies. Application of Hanle-polarimetry to prominence cavities and polar-crown filaments reveals robust “lagomorphic” linear polarization signatures—a rabbit-head pattern diagnosing twisted flux-rope morphology (Bak-Steslicka et al., 2013, Gibson, 2017). Such signatures, forward-modeled from MHD equilibrium, unambiguously signal the presence of flux-rope cores and provide maps of field orientation and twist.

Linear polarization mapping directly identifies coronal magnetic null points in pseudostreamers, visible as the intersection of three Van Vleck lines (L/I minima) and abrupt azimuth flips (Gibson et al., 2017). This capability allows direct localization of reconnection-prone topologies independent of photospheric extrapolation. Furthermore, CoMP yields a polarization-based expansion factor (LPF) quantifying super-radial magnetic flux tube expansion, which is crucial for empirical solar wind speed prediction (Gibson et al., 2017).

CME initiation and propagation is tracked through high-cadence Doppler shift and line-width imaging. CMEs are routinely observed as large, outward-propagating perturbations in v_LOS and line width up to 20–30 km/s, with onsets preceding white-light CME detection by up to 20 minutes (Tian et al., 2013).

6. Data–Model Comparison, Inversion, and Automated Analysis

Community toolchains such as FORWARD integrate CoMP data with analytic and numerical MHD models, providing direct comparison between observed polarimetric and Doppler maps and synthetic observables (Gibson, 2015). Inverse approaches range from direct single-point fitting (solving for local B, θ_B, γ_B, σ_0^2 using Newton-Raphson methods) to global vector-tomographic and forward-model optimization using Radial-basis-function surrogates (ROAM method) (Dalmasse et al., 2016). The latter enables high-dimensional optimization of magnetic models against observed Stokes profiles with dramatic computational speedup, providing orientation recovery to within ≲1° and field strength uncertainties of tens of Gauss when Stokes V is available.

Limitations associated with optically thin LOS integration, ambiguous polarization signatures (notably Van Vleck flips and 180° azimuth ambiguity), and limited S/N at large heights require future efforts in multi-line and multi-viewpoint tomography, higher spatial resolution, and the integration of time-dependent coronal models. Ongoing improvement of inversion routines now incorporates Bayesian techniques and machine learning approaches for rapid, robust field recovery and error quantification (Gibson, 2015, Kramar et al., 2015).

7. Future Directions, Upgrades, and Complementary Facilities

Improvements in instrument design—including the upgrade to UCoMP with finer spatial resolution (∼2.2 Mm/pixel), increased FOV (up to 1.95 R⊙), and extended spectral coverage (530–1083 nm)—are underway or operational (Duckenfield et al., 7 Mar 2025, Yang et al., 15 Jan 2026). Next-generation coronagraphs such as COSMO and DKIST will offer multi-wavelength, full-Stokes polarimetric input with spatial resolution down to ∼0.03″ and polarimetric sensitivity ≲10⁻⁴, enabling routine, high-fidelity 3D coronal magnetometry and direct resolution of sub-Mm wave and turbulence scales (French et al., 2019, Kramar et al., 2015).

Recommendations include:

• Increased spectral sampling across diagnostic lines for improved line-shape analysis
• Higher temporal cadence (≲10 s) to capture eruptive and wave phenomena
• Inclusion of time-varying boundary conditions and real-time boundary forcings in MHD modeling
• Integration of complementary diagnostics (Faraday-rotation, radio, visible lines) for full vector field inversion (Gibson, 2015)

By operating as a multi-channel diagnostic platform for coronal plasma, magnetic field, and wave dynamics, CoMP, and its successors, are positioned as central pillars in observational solar physics, uniquely bridging empirical measurement, data-driven modeling, and physics-based theory across all spatial and temporal scales relevant to the solar corona and heliosphere.