Search for: All records

The Wang–Sheeley–Arge (WSA) model has been in use for decades and remains a popular, economical approach to modeling the solar coronal magnetic field and forecasting conditions in the inner heliosphere. Given its usefulness, it is unsurprising that a number of WSA implementations have been developed by various groups with different computational approaches. While the WSA magnetic field model has traditionally been calculated using a spherical harmonic expansion of the solar magnetic field, finite-difference potential field solutions can offer speed and/or accuracy advantages. However, the creation of new versions of WSA requires that we ensure the solutions from these new models are consistent with established versions and that we quantify for the user community to what degree and in what ways they differ. In this paper, we present side-by-side comparisons of WSA models produced using the traditional, spherical harmonic–based implementation developed by Wang, Sheeley, and Arge with WSA models produced using a recently open-sourced finite-difference code from the CORHEL modeling suite called POT3D. We present comparisons of the terminal solar wind speed and magnetic field at the outer boundaries of the models, weighing these against the variation of the WSA model in the presence of small perturbations in the computational procedure, parameters, and inputs. We also compare the footpoints of magnetic field lines traced from the outer boundaries and the locations of open field in the models. We find that the traced field-line footpoints show remarkable agreement, with the greatest differences near the magnetic neutral line and in the polar regions.

 

Unconditionally stable time stepping schemes are useful and often practically necessary for advancing parabolic operators in multi-scale systems. However, serious accuracy problems may emerge when taking time steps that far exceed the explicit stability limits. In our previous work, we compared the accuracy and performance of advancing parabolic operators in a thermodynamic MHD model using an implicit method and an explicit super time-stepping (STS) method. We found that while the STS method outperformed the implicit one with overall good results, it was not able to damp oscillatory behavior in the solution efficiently, hindering its practical use. In this follow-up work, we evaluate an easy-to-implement method for selecting a practical time step limit (PTL) for unconditionally stable schemes. This time step is used to ‘cycle’ the operator-split thermal conduction and viscosity parabolic operators. We test the new time step with both an implicit and STS scheme for accuracy, performance, and scaling. We find that, for our test cases here, the PTL dramatically improves the STS solution, matching or improving the solution of the original implicit scheme, while retaining most of its performance and scaling advantages. The PTL shows promise to allow more accurate use of unconditionally stable schemes for parabolic operators and reliable use of STS methods.

 

To address Objective II of the National Space Weather Strategy and Action Plan “Develop and Disseminate Accurate and Timely Space Weather Characterization and Forecasts” and US Congress PROSWIFT Act 116–181, our team is developing a new set of open-source software that would ensure substantial improvements of Space Weather (SWx) predictions. On the one hand, our focus is on the development of data-driven solar wind models. On the other hand, each individual component of our software is designed to have accuracy higher than any existing SWx prediction tools with a dramatically improved performance. This is done by the application of new computational technologies and enhanced data sources. The development of such software paves way for improved SWx predictions accompanied with an appropriate uncertainty quantification. This makes it possible to forecast hazardous SWx effects on the space-borne and ground-based technological systems, and on human health. Our models include (1) a new, open-source solar magnetic flux model (OFT), which evolves information to the back side of the Sun and its poles, and updates the model flux with new observations using data assimilation methods; (2) a new potential field solver (POT3D) associated with the Wang–Sheeley–Arge coronal model, and (3) a new adaptive, 4-th order of accuracy solver (HelioCubed) for the Reynolds-averaged MHD equations implemented on mapped multiblock grids (cubed spheres). We describe the software and results obtained with it, including the application of machine learning to modeling coronal mass ejections, which makes it possible to improve SWx predictions by decreasing the time-of-arrival mismatch. The tests show that our software is formally more accurate and performs much faster than its predecessors used for SWx predictions.

 

The arrival time prediction of coronal mass ejections (CMEs) is an area of active research. Many methods with varying levels of complexity have been developed to predict CME arrival. However, the mean absolute error (MAE) of predictions remains above 12 hr, even with the increasing complexity of methods. In this work we develop a new method for CME arrival time prediction that uses magnetohydrodynamic simulations involving data-constrained flux-rope-based CMEs, which are introduced in a data-driven solar wind background. We found that for six CMEs studied in this work the MAE in arrival time was ∼8 hr. We further improved our arrival time predictions by using ensemble modeling and comparing the ensemble solutions with STEREO-A and STEREO-B heliospheric imager data. This was done by using our simulations to create synthetic J-maps. A machine-learning (ML) method called the lasso regression was used for this comparison. Using this approach, we could reduce the MAE to ∼4 hr. Another ML method based on the neural networks (NNs) made it possible to reduce the MAE to ∼5 hr for the cases when HI data from both STEREO-A and STEREO-B were available. NNs are capable of providing similar MAE when only the STEREO-A data are used. Our methods also resulted in very encouraging values of standard deviation (precision) of arrival time. The methods discussed in this paper demonstrate significant improvements in the CME arrival time predictions. Our work highlights the importance of using ML techniques in combination with data-constrained magnetohydrodynamic modeling to improve space weather predictions.

 

Drawing connections between heliospheric spacecraft and solar wind sources is a vital step in understanding the evolution of the solar corona into the solar wind and contextualizing in situ timeseries. Furthermore, making advanced predictions of this linkage for ongoing heliospheric missions, such as Parker Solar Probe (Parker), is necessary for achieving useful coordinated remote observations and maximizing scientific return. The general procedure for estimating such connectivity is straightforward (i.e., magnetic field line tracing in a coronal model) but validating the resulting estimates is difficult due to the lack of an independent ground truth and limited model constraints. In its most recent orbits, Parker has reached perihelia of 13.3R_⊙and moreover travels extremely fast prograde relative to the solar surface, covering over 120° longitude in 3 days. Here we present footpoint predictions and subsequent validation efforts for Parker Encounter 10, the first of the 13.3R_⊙orbits, which occurred in November 2021. We show that the longitudinal dependence of in situ plasma data from these novel orbits provides a powerful method of footpoint validation. With reference to other encounters, we also illustrate that the conditions under which source mapping is most accurate for near‐ecliptic spacecraft (such as Parker) occur when solar activity is low, but also require that the heliospheric current sheet is strongly warped by mid‐latitude or equatorial coronal holes. Lastly, we comment on the large‐scale coronal structure implied by the Encounter 10 mapping, highlighting an empirical equatorial cut of the Alfvèn surface consisting of localized protrusions above unipolar magnetic separatrices.

 

Abstract Flux-rope-based magnetohydrodynamic modeling of coronal mass ejections (CMEs) is a promising tool for prediction of the CME arrival time and magnetic field at Earth. In this work, we introduce a constant-turn flux rope model and use it to simulate the 2012 July 12 16:48 CME in the inner heliosphere. We constrain the initial parameters of this CME using the graduated cylindrical shell (GCS) model and the reconnected flux in post-eruption arcades. We correctly reproduce all the magnetic field components of the CME at Earth, with an arrival time error of approximately 1 hr. We further estimate the average subjective uncertainties in the GCS fittings by comparing the GCS parameters of 56 CMEs reported in multiple studies and catalogs. We determined that the GCS estimates of the CME latitude, longitude, tilt, and speed have average uncertainties of 5.°74, 11.°23, 24.°71, and 11.4%, respectively. Using these, we have created 77 ensemble members for the 2012 July 12 CME. We found that 55% of our ensemble members correctly reproduce the sign of the magnetic field components at Earth. We also determined that the uncertainties in GCS fitting can widen the CME arrival time prediction window to about 12 hr for the 2012 July 12 CME. On investigating the forecast accuracy introduced by the uncertainties in individual GCS parameters, we conclude that the half-angle and aspect ratio have little impact on the predicted magnetic field of the 2012 July 12 CME, whereas the uncertainties in longitude and tilt can introduce relatively large spread in the magnetic field predicted at Earth. 

The NASA-NSF sponsored Space Weather with Quantified Uncertainty (SWQU) project's main objective is to develop a data-driven, time-dependent, open source model of the solar corona and heliosphere. One key component of the SWQU effort is using a data-assimilation flux transport model to generate an ensemble of synchronic radial magnetic field maps as boundary conditions for the coronal field model. To accomplish this goal, we are developing a new Open-source Flux Transport (OFT) software suite. While there are a number of established flux transport models in the community, OFT is distinguished from many of these efforts in 3 key attributes: (1) It is based on modern computing techniques that will allow many realizations to be rapidly computed on multi-core systems and/or GPUs, (2) it is designed to be easily extensible, and (3) OFT will be released as an open source project. OFT consists of three software packages: 1) OFTpy: a python package for data acquisition, database organization, and Carrington map processing, 2) ConFlow: a Fortran code that generates super granular convective flows, and 3) High-Performance Flux Transport (HipFT): a modular, GPU-accelerated Fortran code for modeling surface flux transport with data assimilation. Here, we present the current state of the OFT project, key features and methods of OFTpy, ConFlow, and HipFt, and real-world examples of data-assimilation and flux transport with HipFT. Validation and performance tests are shown, including generating an ensemble of OFT maps.