NASA's Parker Solar Probe (Parker) spacecraft reached its first perihelion of 35.7 solar radii on November 5th, 2018. To aid in mission planning, and in anticipation of the unprecedented measurements to be returned, in late October, we developed a three-dimensional magnetohydrodynamic (MHD) solution for the solar corona and inner heliosphere, driven by the then available observations of the Sun's photospheric magnetic field. Our model incorporates a wave-turbulence-driven (WTD) model to heat the corona. Here, we present our predictions for the structure of the solar corona and the likely {\it in situ} measurements that Parker will be returning over the next few months. We infer that, in the days prior to 1st Encounter, Parker was immersed in wind emanating from a positive-polarity equatorial coronal hole. During the encounter, however, field lines from the spacecraft mapped to another, negative-polarity equatorial coronal hole. Following the encounter, Parker was magnetically connected to the large, positive-polarity northern polar coronal hole. When the Parker data become available, these model results can be used to assist in their calibration and interpretation, and, additionally, provide a global context for interpreting the localized {\it in situ} measurements. In particular, we can identify what types of solar wind Parker encountered, what the underlying magnetic structure was, and how complexities in the orbital trajectory can be interpreted within a global, inertial frame. Ultimately, the measurements returned by Parker can be used to constrain current theories for heating the solar corona and accelerating the solar wind.

We describe the initial version of the Solar Particle Event (SPE) Threat Assessment Tool or STAT. STAT relies on elements of Corona-Heliosphere (CORHEL) and the Earth-Moon-Mars Radiation Environment Module (EMMREM), and allows users to investigate coronal mass ejection (CME) driven SPEs using coupled magnetohydrodynamic (MHD) and focused transport solutions. At the present time STAT focuses on modeling solar energetic particle (SEP) acceleration in and transport from the low corona, where the highest energy SEP events are generated. We illustrate STAT's capabilities with a model of the July 14, 2000 "Bastille Day" event, including innovative diagnostics for understanding the three-dimensional distribution of particle fluxes and their relation to the structure of the underlying CME driver. A preliminary comparison with NOAA GOES measurements is shown.

GPU accelerators have had a notable impact on high-performance computing across many disciplines. They provide high performance with low cost/power, and therefore have become a primary compute resource on many of the largest supercomputers. Here, we implement multi-GPU acceleration into our Solar MHD code (MAS) using OpenACC in a fully portable, single-source manner. Our preliminary implementation is focused on MAS running in a reduced physics "zero-beta" mode. While valuable on its own, our main goal is to pave the way for a full physics, thermodynamic MHD implementation. We describe the OpenACC implementation methodology and challenges. "Time-to-solution" performance results of a production-level flux rope eruption simulation on multi-CPU and multi-GPU systems are shown. We find that the GPU-accelerated MAS code has the ability to run "zero-beta" simulations on a single multi-GPU server at speeds previously requiring multiple CPU server-nodes of a supercomputer.

Solar eruptions are the main driver of space-weather disturbances at the Earth. Extreme events are of particular interest, not only because of the scientific challenges they pose, but also because of their possible societal consequences. Here we present a magnetohydrodynamic (MHD) simulation of the 14 July 2000 Bastille Day eruption, which produced a very strong geomagnetic storm. After constructing a thermodynamic MHD model of the corona and solar wind, we insert a magnetically stable flux rope along the polarity inversion line of the eruption's source region and initiate the eruption by boundary flows. More than 10^33 ergs of magnetic energy are released in the eruption within a few minutes, driving a flare, an EUV wave, and a coronal mass ejection (CME) that travels in the outer corona at about 1500 km/s, close to the observed speed. We then propagate the CME to Earth, using a heliospheric MHD code. Our simulation thus provides the opportunity to test how well in situ observations of extreme events are matched if the eruption is initiated from a stable magnetic-equilibrium state. We find that the flux-rope center is very similar in character to the observed magnetic cloud, but arrives about 8.5 hours later and about 15 degrees too far to the North, with field strengths that are too weak by a factor of about 1.6. The front of the flux rope is highly distorted, exhibiting localized magnetic-field concentrations as it passes 1 AU. We discuss these properties with regard to the development of space-weather predictions based on MHD simulations of solar eruptions.

Many existing models assume that magnetic flux ropes play a key role in solar flares and coronal mass ejections (CMEs). It is therefore important to develop efficient methods for constructing flux-rope configurations constrained by observed magnetic data and the morphology of the pre-eruptive source region. For this purpose, we have derived and implemented a compact analytical form that represents the magnetic field of a thin flux rope with an axis of arbitrary shape and circular cross-sections. This form implies that the flux rope carries axial current $I$ and axial flux $F$, so that the respective magnetic field is the curl of the sum of axial and azimuthal vector potentials proportional to $I$ and $F$, respectively. We expressed the vector potentials in terms of modified Biot-Savart laws whose kernels are regularized at the axis in such a way that, when the axis is straight, these laws define a cylindrical force-free flux rope with a parabolic profile for the axial current density. For the cases we have studied so far, we determined the shape of the rope axis by following the polarity inversion line of the eruptions' source region, using observed magnetograms. The height variation along the axis and other flux-rope parameters are estimated by means of potential field extrapolations. Using this heuristic approach, we were able to construct pre-eruption configurations for the 2009 February 13 and 2011 October 1 CME events. These applications demonstrate the flexibility and efficiency of our new method for energizing pre-eruptive configurations in simulations of CMEs.

The heliospheric magnetic field is of pivotal importance in solar and space physics. The field is rooted in the Sun's photosphere, where it has been observed for many years. Global maps of the solar magnetic field based on full disk magnetograms are commonly used as boundary conditions for coronal and solar wind models. Two primary observational constraints on the models are (1) the open field regions in the model should approximately correspond to coronal holes observed in emission, and (2) the magnitude of the open magnetic flux in the model should match that inferred from in situ spacecraft measurements. In this study, we calculate both MHD and PFSS solutions using fourteen different magnetic maps produced from five different types of observatory magnetograms, for the time period surrounding July, 2010. We have found that for all of the model/map combinations, models that have coronal hole areas close to observations underestimate the interplanetary magnetic flux, or, conversely, for models to match the interplanetary flux, the modeled open field regions are larger than coronal holes observed in EUV emission. In an alternative approach, we estimate the open magnetic flux entirely from solar observations by combining automatically detected coronal holes for Carrington rotation 2098 with observatory synoptic magnetic maps. This approach also underestimates the interplanetary magnetic flux. Our results imply that either typical observatory maps underestimate the Sun's magnetic flux, or a significant portion of the open magnetic flux is not rooted in regions that are obviously dark in EUV and X-ray emission.

We describe a method for the automatic mapping of coronal holes (CH) using simultaneous multi-instrument EUV imaging data. Synchronized EUV images from STEREO/EUVI A&B 195A and SDO/AIA 193A are preprocessed, including PSF deconvolution and the application of data-derived intensity corrections that account for center-to-limb variations (limb brightening) and inter-instrument intensity normalization. We systematically derive a robust limb-brightening correction that takes advantage of unbiased long-term averages of data and respects the physical nature of the problem. The new preprocessing greatly assists in CH detection, allowing for the use of a simplified variable-connectivity two-threshold region growing image segmentation algorithm to obtain consistent detection results. We generate synchronic EUV and CH maps, and show a preliminary analysis of CH evolution. Several data and code products are made available to the community (www.predsci.com/chd): For the period of this study (06/10/2010 to 08/18/14) we provide synchronic EUV and coronal hole map data at 6-hour cadence, data-derived limb-brightening corrections for STEREO/EUVI A&B 195A and SDO/AIA 193A, and inter-instrument correction factors to equate their intensities. We also provide the coronal hole image segmentation code module (ezseg) implemented in both FORTRAN OpenMP and GPU-accelerated C-CUDA. A complete implementation of our coronal hole detection pipeline in the form of a ready-to-use MATLAB driver script euv2chm utilizing ezseg is also made available.

Late on July 23, 2012, the STEREO-A spacecraft encountered a fast forward shock driven by a coronal mass ejection launched from the Sun earlier that same day. The estimated travel time of the disturbance ($\sim 20$ hrs), together with the massive magnetic field strengths measured within the ejecta ($> 100$nT), made it one of the most extreme events observed during the space era. In this study, we examine the properties of the shock wave. Because of an instrument malfunction, plasma measurements during the interval surrounding the CME were limited, and our approach has been modified to capitalize on the available measurements and suitable proxies, where possible. We were able to infer the following properties. First, the shock normal was pointing predominantly in the radial direction (${\bf n} = 0.97 {\bf e}_r -0.09 {\bf e}_t -0.23 {\bf e}_n$). Second, the angle between ${\bf n}$ and the upstream magnetic field, $\theta_{Bn}$, was estimated to be $\approx 34^{\circ}$, making the shock "quasi-parallel," and supporting the idea of an earlier "preconditioning" ICME. Third, the shock speed was estimated to be $\approx 3300$ km s$^{-1}$. Fourth, the sonic Mach number, $M_s$, for this shock was $\sim 28$. We support these results with an idealized numerical simulation of the ICME. Finally, we estimated the change in ram pressure upstream of the shock to be $\sim 5$ times larger than the pressure from the energetic particles, suggesting that this was not a standard "steady-state" cosmic-ray modified shock (CRMS). Instead it might represent an early, transient phase in the evolution of the CRMS.