Signatures of Damping Nonlinear Oscillations by KHI-induced Turbulence in Synthetic Observations

Large-amplitude decaying kink oscillations of coronal loops are strongly influenced by nonlinear processes, such as Kelvin-Helmholtz instability (KHI) and turbulence, though comprehensive theory and observational confirmation remain limited. Building on the recently developed theory on nonlinear damping by KHI-induced turbulence in impulsively driven transverse loop oscillations, we investigate its observational signatures using 3D magnetohydrodynamic simulations and forward-modelled EUV images. The simulated oscillations exhibit time-varying frequency shifts and damping rates, which are broadly consistent with nonlinear turbulence-damping theory. Additionally, they exhibit excitation of higher-order modes, slightly increased periods relative to the linear kink period, and reduced displacement amplitudes. These features are generally preserved in synthetic observations, though resolving higher-order modes requires higher spatial resolution than currently available. For loops embedded in a hotter background, hotter channels (e.g., 193 Angstroms) are more sensitive to boundary dynamics, thus their oscillations decay faster with smaller displacements and larger phase shifts than those in cooler channels (e.g., 171 Angstroms). Comparisons of simulated and synthetic oscillations show close agreement at the early stage. At later times, synthetic oscillations exhibit smaller displacements and larger phase shifts, due to turbulence-induced asymmetry in the loop cross-section. Bayesian fitting shows that the initial oscillation amplitude and kink period are robustly constrained, whereas parameters controlling the damping profile are degenerate, indicating that additional observables would aid reliable seismological inference. These results provide a quantitative basis for identifying nonlinear damping and detecting KHI-driven turbulence in transverse loop oscillations.

💡 Research Summary

This paper investigates the observational signatures of nonlinear damping of large‑amplitude transverse (kink) oscillations in coronal loops, focusing on the role of Kelvin‑Helmholtz instability (KHI)‑driven turbulence. Building on the analytic framework introduced by Zhong et al. (2025), which predicts a time‑varying damping rate and a frequency drift as the loop transitions from a KHI‑developing stage to a fully turbulent stage, the authors use three‑dimensional ideal MHD simulations together with forward‑modeled EUV images to test these predictions.

The simulations model a straight, density‑enhanced loop (radius R = 0.5 loop‑diameters) embedded in a uniform background. The density contrast ζ = ρᵢ/ρₑ is varied (0.3–10), and either a temperature ratio Tᵢ/Tₑ = 0.5 or a pressure ratio Pᵢ/Pₑ = 3/2 is imposed to explore the effect of thermal imbalance. An impulsive transverse velocity pulse is applied such that the non‑linearity parameter V₀L/(CₖR) ≥ 1, guaranteeing large‑amplitude motions. The loop is half‑modeled (foot‑to‑apex) and mirrored to obtain the full cross‑section, saving computational cost while preserving symmetry.

Key diagnostics extracted from the simulations include the centre‑of‑mass velocity V_CoM, the displacement ξ(t), the initial velocity Vᵢ (after the impulsive kick), the shear ΔV across the loop boundary, and the mixing‑layer growth coefficient C₁, which quantifies the efficiency of KHI‑induced mixing. The authors confirm the two‑stage turbulence evolution reported in earlier works (Hillier et al. 2023, 2024): an early growth phase where the mixing‑layer thickness h ∝ t, followed by a decay phase where turbulent energy decays as ∝ t⁻². The growth coefficient C₁ ranges from ~0.1 to 1, with larger ζ yielding higher C₁ and faster damping.

The simulated oscillations display several hallmark features of nonlinear damping: (i) a frequency drift proportional to the square of the amplitude (Δω ∝ A²), consistent with Ruderman & Goossens (2014); (ii) a time‑dependent damping rate that accelerates as the KHI develops; (iii) excitation of higher‑order azimuthal modes (m ≥ 2) at twice the kink frequency, although their amplitudes remain modest. The fundamental kink period is slightly longer than the linear prediction, reflecting the additional inertia contributed by the turbulent mixing layer.

To assess observability, the authors forward‑model the simulation data with the FoMo code, generating synthetic AIA images in the 131 Å, 171 Å, 193 Å, and 211 Å channels. The loop radius is scaled to 1 Mm, external density to 10⁹ cm⁻³, and internal temperature to 1 MK. The synthetic data are rebinned to AIA’s 440 km pixel⁻¹ resolution, convolved with the instrument PSF, and Poisson noise is added. The analysis shows clear channel‑dependent behaviour: the cooler 171 Å channel, whose temperature response peaks near 1 MK, captures the loop core and exhibits relatively large displacement amplitudes with slower decay. In contrast, the hotter 193 Å and 211 Å channels are more sensitive to the loop boundary and the surrounding hotter plasma; they display faster amplitude decay, smaller displacements, and larger phase shifts (up to 0.3–0.6 π) after a few periods. This reflects the fact that KHI‑generated roll‑ups and turbulence preferentially affect the boundary region, which contributes more strongly to the emission in hotter lines.

Higher‑order modes are present in the synthetic intensity time‑series but remain below AIA’s spatial resolution (≈440 km), implying that current instruments can only hint at their existence through subtle changes in the apparent period or damping profile.

Finally, a Bayesian inference framework is applied to the synthetic displacement curves, comparing three damping models: (a) linear exponential damping, (b) nonlinear damping with a constant rate, and (c) nonlinear damping with a time‑varying rate as predicted by the turbulence theory. The posterior distributions robustly constrain the initial amplitude A₀ and the linear kink period Pₖ, but the damping‑related parameters (τ, the characteristic damping time, and α, the nonlinearity exponent) are strongly degenerate. This degeneracy arises because a single channel’s intensity integrates over the entire cross‑section, mixing contributions from the turbulent boundary and the relatively undisturbed core. The authors argue that additional observables—such as multi‑channel phase differences, line‑width broadening, or high‑resolution imaging from DKIST or future missions—are required to break this degeneracy and enable reliable seismological inference of KHI‑driven turbulence.

In summary, the study provides a quantitative bridge between theory, numerical simulation, and synthetic observation for KHI‑induced turbulent damping of coronal loop kink oscillations. It identifies three observable signatures—time‑varying damping rate, frequency drift, and channel‑dependent phase shift—that can be used to detect nonlinear damping in existing AIA data, while emphasizing the need for higher spatial and spectral resolution to fully exploit these diagnostics for coronal seismology.