Accessibility navigation

A computationally efficient, time-dependent model of the solar wind for use as a surrogate to three-dimensional numerical magnetohydrodynamic simulations

Owens, M., Lang, M., Barnard, L., Riley, P., Ben-Nun, M., Scott, C. J., Lockwood, M., Reiss, M. A., Arge, C. N. and Gonzi, S. (2020) A computationally efficient, time-dependent model of the solar wind for use as a surrogate to three-dimensional numerical magnetohydrodynamic simulations. Solar Physics, 295 (3). 43. ISSN 0038-0938

Text (Open Access) - Published Version
· Available under License Creative Commons Attribution.
· Please see our End User Agreement before downloading.


It is advisable to refer to the publisher's version if you intend to cite from this work. See Guidance on citing.

To link to this item DOI: 10.1007/s11207-020-01605-3


Near-Earth solar-wind conditions, including disturbances generated by coronal mass ejections (CMEs), are routinely forecast using three-dimensional, numerical magnetohydrodynamic (MHD) models of the heliosphere. The resulting forecast errors are largely the result of uncertainty in the near-Sun boundary conditions, rather than heliospheric model physics or numerics. Thus ensembles of heliospheric model runs with perturbed initial conditions are used to estimate forecast uncertainty. MHD heliospheric models are relatively cheap in computational terms, requiring tens of minutes to an hour to simulate CME propagation from the Sun to Earth. Thus such ensembles can be run operationally. However, ensemble size is typically limited to 101 to 102 members, which may be inadequate to sample the relevant high-dimensional parameter space. Here, we describe a simplified solar-wind model that can estimate CME arrival time in approximately 0.01 seconds on a modest desktop computer and thus enables significantly larger ensembles. It is a one-dimensional, incompressible, hydrodynamic model, which has previously been used for the steady-state solar wind, but it is here used in time-dependent form. This approach is shown to adequately emulate the MHD solutions to the same boundary conditions for both steady-state solar wind and CME-like disturbances. We suggest it could serve as a “surrogate” model for the full three-dimensional MHD models. For example, ensembles of 105 to 106 members can be used to identify regions of parameter space for more detailed investigation by the MHD models. Similarly, the simplicity of the model means it can be rewritten as an adjoint model, enabling variational data assimilation with MHD models without the need to alter their code. The model code is available as an Open Source download in the Python language.

Item Type:Article
Divisions:Science > School of Mathematical, Physical and Computational Sciences > Department of Meteorology
ID Code:89616


Downloads per month over past year

University Staff: Request a correction | Centaur Editors: Update this record

Page navigation