Modeling of Coronal Mass Ejection Originated from a Sheared Arcade of Realistic Active-Region Scale and Its Propagation in the Heliosphere: Methodology

TL;DR

This paper introduces a comprehensive methodology for simulating coronal mass ejections from active regions to Earth, combining advanced adaptive mesh refinement and semi-relativistic corrections for realistic, efficient modeling.

Contribution

The authors develop an end-to-end CME modeling approach with high-resolution AR scale resolution and realistic magnetic fields, enabling detailed simulation of CME initiation and propagation.

Findings

Simulated CME shows a three-part structure in synthetic coronagraph images.

Model captures shock formation, density compression, and southward Bz at 1 au.

Simulation predicts CME arrival at 1 au with a two-day forecast lead time.

Abstract

Simulating coronal mass ejections (CMEs) from their origin in active regions (ARs) to their propagation to Earth remains challenging, particularly when aiming to resolve AR scales and employ realistic magnetic field strengths without compromising computational efficiency. Here we present a methodology for end-to-end CME modeling that addresses these challenges. Three nested magnetohydrodynamic simulations are coupled to jointly cover the heliosphere from solar surface to beyond $1.5$ au. A block-structured adaptive mesh refinement scheme is employed to achieve $\sim 700$ km resolution in the low corona, allowing AR scales to be resolved while maintaining the total grid count below $10^8$ across the entire computational domain. A semi-relativistic Boris correction combined with a relativistic mass-density factor is used to handle magnetic field strengths up to $10^3$ G without prohibitively small time steps. Using this model, we simulate the emergence of a bipolar AR into the corona, the initiation of a CME by shearing of the AR core field and the subsequent evolution. Our simulation captures its pre-eruption energy buildup, triggering by magnetic reconnection, rapid acceleration, and propagation to 1 au and beyond. The simulated CME exhibits a three-part structure in synthetic coronagraph images and a torus-shaped flux rope in the heliosphere, with synthetic in-situ observations showing shock formation, density compression, and a prolonged southward $B_z$ component at 1 au. The entire simulation requires about one day on a moderately sized cluster (e.g., $600$ processors), while the simulated CME takes three days to arrive at $1$ au, offering a lead time of two days if used for forecasting.