Formation of a Magnetic Flux Rope Prior to the Eruption: Insight from a Radiative MHD Simulation of Active Region Emergence

Magnetic flux ropes (MFRs) are fundamental magnetic structures in solar eruptions, whose formation is generally attributed to (1) the emergence of subsurface flux tubes or (2) flux cancellation driven by photospheric horizontal flows and magnetic reconnection. Both mechanisms can operate simultaneously during active region evolution, making their relative contributions challenging to quantify. Here, we analyze the formation of a flux rope in a MURaM radiative magnetohydrodynamic (RMHD) simulation, which formed and evolved for approximately three hours before an M-class flare. The formation process is quantified by magnetic helicity flux, which drives the non-potential evolution of magnetic field, with its advection and shear terms on the photosphere corresponding to the emergence and photospheric horizontal flows, respectively. Examining the helicity injected into the flux rope through the photosphere, we find both terms increase significantly as the eruption approaches, with the shear term prevailing overall. Height-dependent analysis of helicity flux, together with magnetic field and velocity distributions, further reveals a gradual transition from the shear to the advection term with an increasing altitude, which is driven by magnetic reconnection above the photosphere. Our results provide quantitative evidence that flux cancellation governs flux rope formation, arising naturally from magnetic field reorganization during active region evolution: as flux emergence transports magnetic flux upward, photospheric shearing motions adjust magnetic field and inject helicity into solar atmosphere, and magnetic reconnection ultimately assembles the main body of flux ropes.

💡 Research Summary

**
This paper investigates how a magnetic flux rope (MFR) forms prior to a solar eruption using a state‑of‑the‑art radiative magnetohydrodynamic (RMHD) simulation performed with the MURaM code. The simulation domain spans 196 × 196 × 123 Mm³ with a uniform grid of 192 km (horizontal) and 64 km (vertical). A realistic bottom boundary supplies velocity and magnetic‑field data from a solar convective dynamo simulation, allowing self‑consistent flux emergence and photospheric flows. Over ~48 h of simulated time the active region produces more than 100 eruptive events; the authors focus on an M‑class flare and its associated pre‑eruptive flux rope, tracking its evolution for about three hours before eruption.

The flux rope is identified by two topological diagnostics: the twist number (Tw) and the squashing factor (Q). In a horizontal plane at β≈1 (z≈0.83 Mm) a coherent region of negative Tw is surrounded by a high‑Q quasi‑separatrix layer (QSL). Using a region‑growing algorithm on the smoothed Tw map and a log Q threshold of 2–3, the authors isolate the rope’s cross‑section and reconstruct its three‑dimensional field lines. The onset of a well‑defined rope is set at t = −186.3 min relative to the GOES peak.

Magnetic helicity flux through the photosphere is decomposed into a shear term (associated with horizontal motions) and an advection term (associated with vertical emergence). Both terms increase markedly as the eruption approaches, but the shear term dominates the total helicity budget throughout the formation phase, contributing roughly twice as much as the advection term. Height‑dependent analysis reveals a transition: near the photosphere the shear term is dominant, while at heights of 5–10 Mm the advection term becomes comparable, indicating that reconnection above the photosphere transfers helicity upward.

The study also documents persistent flux cancellation at the polarity inversion line, driven by converging/shearing horizontal flows. This cancellation creates current sheets that reconnect, converting arcade‑like fields into the twisted rope. Thus, flux emergence supplies magnetic flux upward, photospheric shearing injects helicity, and reconnection assembles the rope’s core—supporting the view that flux cancellation is the primary mechanism governing flux‑rope formation in this realistic AR evolution.

By coupling a realistic radiative MHD simulation with quantitative helicity‑flux diagnostics, the paper demonstrates that (1) the MURaM framework can self‑consistently reproduce the full life cycle of an active region, (2) helicity‑flux decomposition provides a robust tool to separate emergence and cancellation contributions, and (3) the identified transition from shear‑ to advection‑dominated helicity with height offers new insight into how reconnection mediates rope assembly. These results have direct implications for interpreting observations of helicity injection and for improving eruption‑forecasting models.