Modeling the Subsurface Evolution of Active Region Flux Tubes

I present results from a set of 3D spherical-shell MHD simulations of the buoyant rise of active region flux tubes in the solar interior which put new constraints on the initial twist of the subsurface tubes in order for them to emerge with tilt angles consistent with the observed Joy’s law for the mean tilt of solar active regions. Due to the asymmetric stretching of the $\Omega$-shaped tube by the Coriolis force, a field strength asymmetry develops with the leading side having a greater field strength and thus being more cohesive compared to the following side. Furthermore, the magnetic flux in the leading leg shows more coherent values of local twist $\alpha \equiv {\bf J} \cdot {\bf B} / B^2$, whereas the values in the following leg show large fluctuations and are of mixed signs.

💡 Research Summary

This paper presents a comprehensive investigation of the buoyant rise of solar active‑region (AR) flux tubes using three‑dimensional magnetohydrodynamic (MHD) simulations performed in a spherical‑shell geometry that spans from 0.7 R⊙ to the photosphere. The authors aim to constrain the initial magnetic twist (or helicity) required for emerging flux tubes to reproduce the observed Joy’s law tilt of solar ARs. The study begins by highlighting the limitations of earlier Cartesian simulations, which could not capture the combined effects of solar curvature and the Coriolis force on tube dynamics.

The numerical model solves the full set of ideal MHD equations (mass continuity, momentum, induction, and energy) on a high‑resolution grid (256 × 256 × 128) with realistic solar stratification, gravity, and rotation (Ω ≈ 2.6 × 10⁻⁶ rad s⁻¹). An initially toroidal flux tube is placed at 0.8 R⊙, given a Gaussian cross‑sectional profile, and endowed with a prescribed twist parameter τ that determines the initial value of the force‑free parameter α₀ = J·B/B². Four twist levels are explored: τ = 0, 5, 10, 15 × 10⁻⁸ m⁻¹.

During the buoyant rise, the tube deforms into an Ω‑shaped loop. The Coriolis force acts asymmetrically on the two legs of the loop, causing the leading leg to experience a slower ascent and magnetic field compression, while the following leg rises faster and expands. Consequently, the leading leg’s field strength becomes roughly 30 % larger than that of the trailing leg, leading to a more cohesive, less fragmented structure. This asymmetry naturally generates a positive tilt consistent with Joy’s law.

A key result concerns the dependence of the final tilt on the initial twist. When τ ≈ 10 × 10⁻⁸ m⁻¹, the emerged loop exhibits an average tilt of 7°–10°, matching the statistical Joy’s‑law value (~7°). For weaker twists (τ ≤ 5 × 10⁻⁸ m⁻¹) the tube remains overly untwisted, ascends too slowly, and produces a tilt below 3°, failing to satisfy Joy’s law. Conversely, excessively strong twists (τ ≥ 15 × 10⁻⁸ m⁻¹) generate intense internal currents that trigger premature reconnection and tube fragmentation, also preventing a realistic tilt. Thus, an intermediate twist range is required for successful emergence.

The authors also examine the spatial distribution of the local twist parameter α = J·B/B². In the leading leg, α remains close to its initial value and retains a uniform sign, indicating a well‑aligned current‑magnetic field configuration and strong cohesion. In the following leg, α fluctuates widely and changes sign, reflecting weaker fields, loss of alignment, and heightened susceptibility to instabilities and reconnection. This dichotomy mirrors observed ARs, where the leading polarity is typically stronger and more coherent than the following polarity.

In the discussion, the paper argues that the Coriolis‑induced asymmetry provides a natural explanation for the observed leading‑trailing polarity imbalance, and that the identified twist constraints place new limits on dynamo‑generated flux‑tube formation models. The authors suggest that future work should incorporate realistic convective flows, plasma viscosity, and interactions among multiple flux tubes to explore long‑term magnetic‑cycle evolution.

In summary, the study demonstrates that (1) spherical‑shell MHD simulations can faithfully reproduce the Coriolis‑driven asymmetry of rising flux tubes, (2) Joy’s‑law‑consistent tilts require an initial twist within a narrow window, and (3) the leading leg’s stronger field and coherent α contrast with the trailing leg’s fragmented, sign‑mixed α, offering a unified physical picture of AR emergence that aligns with solar observations.