Enlil (heliospheric model)

Enlil is a time-dependent three-dimensional magnetohydrodynamic model of the inner heliosphere that simulates the background solar wind and the interplanetary propagation of coronal mass ejections, CMEs. It is a core component of the operational WSA–Enlil modeling system used for space weather forecasting and research at the National Oceanic and Atmospheric Administration Space Weather Prediction Center, at NASA's Community Coordinated Modeling Center, and at other national forecasting centers.^[1][2][3]

NOAA states the system provides "1-4 day advance warning" of solar-wind structures and Earth-directed CMEs that can drive geomagnetic storms.^[1] The code is named for Enlil, the Sumerian god of wind and storms, reflecting the model's focus on solar-wind dynamics.^[4][5]

Enlil emerged from work by Czech-American space physicist Dusan Odstrcil on three-dimensional simulations of the solar wind and CMEs in the late 1990s. Foundational Journal of Geophysical Research papers modeled CME propagation in a structured solar wind, one case for events launched within the streamer belt and a companion for events launched adjacent to the streamer belt.^[6][7] A 2003 review set out the numerical approach for modeling three-dimensional solar-wind structure and transient disturbances that later informed the operational code.^[8]

The coupled prediction system that is now known as WSA–Enlil was developed when American space weather researcher Charles N. Arge and Odstrcil began coupling the semi-empirical Wang–Sheeley–Arge, WSA, coronal model to Enlil's heliospheric domain in the early 2000s.^[9][10] A parallel line of work demonstrated the coupling of coronal MHD outputs from the Magnetohydrodynamics Around a Sphere, MAS, model into Enlil for Sun-to-Earth event simulations, which American space physicist Janet G. Luhmann described as "Sun-to-Earth space weather."^[11][12]

NOAA began routine operations of WSA–Enlil in 2011 and continues to archive and distribute model outputs.^[13] In April 2023 NOAA deployed WSA–Enlil v3.0, combining WSA-5.4 with Enlil-2.9e and updating magnetogram processing, with the goal of improving ambient-wind fidelity and CME arrival forecasts.^[14]

Enlil solves the ideal MHD equations for plasma mass, momentum, and energy density, together with the interplanetary magnetic field, using a flux-corrected transport scheme on a spherical grid.^[2] The inner radial boundary is placed beyond the solar sonic point, typically at 12.5 R_☉ solar radii when driven by WSA or 30 R_☉ solar radii when driven by MAS, and the computational domain is often configured from about 0.1 to about 2 AU in radius for operational runs.^[2][15][16] At the inner boundary Enlil commonly assumes equal electron and proton temperatures and reconstructs the tangential IMF from the radial component and solar-wind speed under Parker-spiral geometry.^[17]

To model CME transients, forecasters or automated tools estimate CME parameters from coronagraph imagery and insert a hydrodynamic ejecta at the Enlil inner boundary using a cone geometry. Inputs include onset time, direction, speed, and angular width, sometimes with density and thickness scalings.^[18][19] American space physicist V. J. Pizzo emphasized that, in this configuration, there is "no magnetic cloud component," which limits predictions of the CME's internal magnetic structure and geoeffectiveness, although arrival time and bulk properties can be forecast.^[19]

WSA–Enlil is operated by NOAA's Space Weather Prediction Center to provide global predictions of solar wind conditions and CME arrivals at Earth and other locations throughout the heliosphere. SWPC publishes real-time plots and maintains an archive of model outputs, with continuously available data dating back to November 2011.^[1][13] NASA's Community Coordinated Modeling Center operates WSA–Enlil for research purposes, mission support, and event-driven requests. The model is also used operationally by the UK Met Office and other national forecasting centers. Odstrcil describes the model as providing "near-real-time prediction" of heliospheric space weather conditions across multiple spacecraft domains.^[3] Model outputs are distributed as an open dataset to support research and applications.^[20]

Operational runs typically use the WSA model to specify the background solar-wind speed and radial magnetic field structure at the inner boundary of 21.5 R_☉ solar radii. The WSA model calculates wind speed based on two key parameters: the magnetic-flux-tube expansion factor and the distance from the magnetic footpoint to the nearest coronal hole boundary. These parameters are derived from synoptic magnetograms that map the Sun's magnetic field structure.^[10][17] In NOAA's 2023 upgrade to WSA–Enlil v3.0, the system adopted zero-point-corrected GONG synoptic magnetograms and implemented a retuned velocity relationship in WSA-5.4, with the goal of improving ambient solar wind predictions at Earth's L1.^[14]

Researchers have used WSA–Enlil to maintain continuous upstream solar wind conditions at Mars for the MAVEN mission and to support analyses involving Solar Orbiter, Parker Solar Probe, STEREO, and other spacecraft.^[17][3] WSA–Enlil outputs also serve as inputs for solar energetic particle forecasting tools, including the SEPMOD approach developed by space physicist Janet G. Luhmann and colleagues.^[21]

An evaluation of 25 events during the first year of SWPC operations found an average absolute error of approximately 7.5 hours.^[18] A larger verification study conducted at NASA's CCMC, covering March 2010 to December 2016 and analyzing more than 1,800 CMEs, found a mean absolute error of 10.4 ± 0.9 hours, with a tendency toward early arrival predictions.^[22] Ensemble forecasting techniques applied to WSA–Enlil have been used to estimate uncertainty ranges, with reported mean absolute errors near 12 hours in early real-time trials and smaller errors for selected cases.^[23] Comparative studies of community submissions to NASA's CME Scoreboard indicate typical arrival-time uncertainties of approximately 10 hours for both physics-based and analytical models.^[24]

The operational "cone" insertion method treats CME ejecta hydrodynamically, which means Enlil cannot predict the internal magnetic field structure of a CME at 1 AU. Space physicist V. J. Pizzo describes this limitation as having "no magnetic cloud component," which restricts forecasts of southward IMF Bz and geomagnetic storm intensity, focusing predictions instead on arrival time and shock properties.^[19] Forecast accuracy depends on the quality of the modeled ambient solar wind and the precision of CME input parameters. WSA–Enlil is a mature, widely used tool and a "workhorse" of both research and operational forecasting, but could better support specification of inner boundary conditions and multi-viewpoint constraints on CME geometry.^[9][3][24]

• Space weather
• Magnetohydrodynamics
• Wang–Sheeley–Arge model

• NOAA Space Weather Prediction Center, WSA–Enlil product page
• NASA CCMC, Enlil model page
• NOAA announcement of WSA–Enlil v3.0
• NOAA WSA–Enlil Open Data registry