The National Science Foundation (NSF) Daniel K. Inouye Solar Telescope (DKIST) is about to start operations at the summit of Haleakala (Hawaii). DKIST will join the early science phases of the NASA and ESA Parker Solar Probe and Solar Orbiter encounter missions. By combining in-situ measurements of the near-sun plasma environment and detail remote observations of multiple layers of the Sun, the three observatories form an unprecedented multi-messenger constellation to study the magnetic connectivity inside the solar system. This white paper outlines the synergistic science that this multi-messenger suite enables.

The topology of the coronal magnetic field produces a strong impact on the properties of the solar corona and presumably on the origin of the slow solar wind. To advance our understanding of this impact, we revisit the concept of so-called slip-back mapping (Titov et al. 2009) and adapt it for determining open, closed, and disconnected flux systems that are formed in the solar corona by magnetic reconnection during a given time interval. The developed method allows us, in particular, to describe the magnetic flux transfer between open and closed flux regions via so-called interchange reconnection with an unprecedented level of details. We illustrate the application of this method to the analysis of the global MHD evolution of the solar corona that is driven by an idealized differential rotation of the photospheric plasma.

We simulate a coronal mass ejection (CME) using a three-dimensional magnetohydrodynamic (MHD) code that includes coronal heating, thermal conduction, and radiative cooling in the energy equation. The magnetic flux distribution at 1 R$_s$ is produced by a localized subsurface dipole superimposed on a global dipole field, mimicking the presence of an active region within the global corona. Transverse electric fields are applied near the polarity inversion line to introduce a transverse magnetic field, followed by the imposition of a converging flow to form and destabilize a flux rope, producing an eruption. We examine the quantities responsible for plasma heating and cooling during the eruption, including thermal conduction, radiation, adiabatic effects, coronal heating, and ohmic heating. We find that ohmic heating is an important contributor to hot temperatures in the current sheet region early in the eruption, but in the late phase adiabatic compression plays an important role in heating the plasma there. Thermal conduction also plays an important role in the transport of thermal energy away from the current sheet region throughout the reconnection process, producing a ``thermal halo'' and widening the region of high temperatures. We simulate emission from solar telescopes for this eruption and find that there is evidence for emission from heated plasma above the flare loops late in the eruption, when the adiabatic heating is the dominant heating term. These results provide an explanation for hot supra-arcade plasma sheets that are often observed in X-rays and extreme ultraviolet wavelengths during the decay phase of large flares.

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.

Ion fractional charge states, measured in situ in the heliosphere, depend on the properties of the plasma in the inner corona. As the ions travel outward in the solar wind and the electron density drops, the charge states remain essentially unaltered or "frozen in". Thus they can provide a powerful constraint on heating models of the corona and acceleration of the solar wind. We have implemented non-equilibrium ionization calculations into a 1D wave-turbulence-driven (WTD) hydrodynamic solar wind model and compared modeled charge states with the Ulysses 1994-5 in situ measurements. We have found that modeled charge state ratios of $C^{6+}/C^{5+}$ and $O^{7+}/O^{6+}$, among others, were too low compared with Ulysses measurements. However, a heuristic reduction of the plasma flow speed has been able to bring the modeled results in line with observations, though other ideas have been proposed to address this discrepancy. We discuss implications of our results and the prospect of including ion charge state calculations into our 3D MHD model of the inner heliosphere.

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.

Observations of coronal loops have identified several common loop characteristics, including that loops appear to cool and have higher than expected densities. Two potential heating scenarios have been suggested to explain these observations. One scenario is that the loops are heated by a series of small-scale impulsive heating events, or nanoflares. Another hypothesis is that the heating is quasi-steady and highly-stratified, i.e., ``footpoint heating'. The goal of this paper is to identify observables that can be used to differentiate between these two heating scenarios. For footpoint heating, we vary the heating magnitude and stratification, for impulsive heating, we vary the heating magnitude. We use one-dimensional hydrodynamic codes to calculate the resulting temperature and density evolution and expected lightcurves in four channels of AIA and one channel of XRT. We consider two principal diagnostics: the time lag between the appearance of the loop in two different channels, and the ratio of the peak intensities of the loop in the two channels. We find that 1) footpoint heating can predict longer time lags than impulsive heating, 2) footpoint heating can predict zero or negative time lags, 3) the intensity ratio expected from impulsive heating is confined to a narrow range, while footpoint heating predicts a wider range of intensity ratios, and 4) the range of temperatures expected in impulsive heating is broader than the range of temperatures expected in footpoint heating. This preliminary study identifies observables that may be useful in discriminating between heating models in future work.

Seven different models are applied to the same problem of simulating the Sun's coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models.

Both direct observations and reconstructions from various datasets, suggest that conditions were radically different during the Maunder Minimum (MM) than during the space era. Using an MHD model, we develop a set of feasible solutions to infer the properties of the solar wind during this interval. Additionally, we use these results to drive a global magnetospheric model. Finally, using the 2008/2009 solar minimum as an upper limit for MM conditions, we use results from the International Reference Ionosphere (ILI) model to speculate on the state of the ionosphere. The results describe interplanetary, magnetospheric, and ionospheric conditions that were substantially different than today. For example: (1) the solar wind density and magnetic field strength were an order of magnitude lower; (2) the Earth's magnetopause and shock standoff distances were a factor of two larger; and (3) the maximum electron density in the ionosphere was substantially lower.

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.

Using a potential field source surface (PFSS) model, we recently analyzed the global topology of the background coronal magnetic field for a sequence of coronal mass ejections (CMEs) that occurred on 2010 August 1-2. Here we repeat this analysis for the background field reproduced by a magnetohydrodynamic (MHD) model that incorporates plasma thermodynamics. As for the PFSS model, we find that all three CME source regions contain a coronal hole that is separated from neighboring coronal holes by topologically very similar pseudo-streamer structures. However, the two models yield very different results for the size, shape, and flux of the coronal holes. We find that the helmet-streamer cusp line, which corresponds to a source-surface null line in the PFSS model, is structurally unstable and does not form in the MHD model. Our analysis indicates that generally, in MHD configurations, this line rather consists of a multiple-null separator passing along the edge of disconnected flux regions. Some of these regions are transient and may be the origin of so-called streamer blobs. We show that the core topological structure of such blobs is a three-dimensional "plasmoid", consisting of two conjoined flux ropes of opposite handedness, which connect at a spiral null point of the magnetic field. Our analysis reveals that such plasmoids appear also in pseudo-streamers on much smaller scales. These new insights into the coronal magnetic topology provide some intriguing implications for solar energetic particle events and for the properties of the slow solar wind.

We study the link between an expanding coronal shock and the energetic particles measured near Earth during the Ground Level Enhancement (GLE) of 17 May 2012. We developed a new technique based on multipoint imaging to triangulate the 3-D expansion of the shock forming in the corona. It uses images from three vantage points by mapping the outermost extent of the coronal region perturbed by the pressure front. We derive for the first time the 3-D velocity vector and the distribution of Mach numbers, $M_{FM}$, of the entire front as a function of time. Our approach uses magnetic field reconstructions of the coronal field, full magneto-hydrodynamic simulations and imaging inversion techniques. We find that the highest $M_{FM}$ values appear near the coronal neutral line within a few minutes of the Coronal Mass Ejection (CME) onset; this neutral line is usually associated with the source of the heliospheric current and plasma sheet. We illustrate the variability of the shock speed, shock geometry and Mach number along different modeled magnetic field lines. Despite the level of uncertainty in deriving the shock Mach numbers, all employed reconstruction techniques show that the release time of GeV particles occurs when the coronal shock becomes super-critical ($M_{FM}>3$). Combining in-situ measurements with heliospheric imagery, we also demonstrate that magnetic connectivity between the accelerator (the coronal shock of 17 May 2012) and the near-Earth environment is established via a magnetic cloud that erupted from the same active region roughly five days earlier.

Transient collimated plasma eruptions in the solar corona, commonly known as coronal (or X-ray) jets, are among the most interesting manifestations of solar activity. It has been suggested that these events contribute to the mass and energy content of the corona and solar wind, but the extent of these contributions remains uncertain. We have recently modeled the formation and evolution of coronal jets using a three-dimensional (3D) magnetohydrodynamic (MHD) code with thermodynamics in a large spherical domain that includes the solar wind. Our model is coupled to 3D MHD flux-emergence simulations, i.e, we use boundary conditions provided by such simulations to drive a time-dependent coronal evolution. The model includes parametric coronal heating, radiative losses, and thermal conduction, which enables us to simulate the dynamics and plasma properties of coronal jets in a more realistic manner than done so far. Here we employ these simulations to calculate the amount of mass and energy transported by coronal jets into the outer corona and inner heliosphere. Based on observed jet-occurrence rates, we then estimate the total contribution of coronal jets to the mass and energy content of the solar wind to (0.4-3.0) % and (0.3-1.0) %, respectively. Our results are largely consistent with the few previous rough estimates obtained from observations, supporting the conjecture that coronal jets provide only a small amount of mass and energy to the solar wind. We emphasize, however, that more advanced observations and simulations are needed to substantiate this conjecture.

To simulate the energy balance of coronal plasmas on macroscopic scales, we often require the specification of the coronal heating mechanism in some functional form. To go beyond empirical formulations and to build a more physically motivated heating function, we investigate the wave-turbulence-driven (WTD) phenomenology for the heating of closed coronal loops. Our implementation is designed to capture the large-scale propagation, reflection, and dissipation of wave turbulence along a loop. The parameter space of this model is explored by solving the coupled WTD and hydrodynamic evolution in 1D for an idealized loop. The relevance to a range of solar conditions is also established by computing solutions for over one hundred loops extracted from a realistic 3D coronal field. Due to the implicit dependence of the WTD heating model on loop geometry and plasma properties along the loop and at the footpoints, we find that this model can significantly reduce the number of free parameters when compared to traditional empirical heating models, and still robustly describe a broad range of quiet-sun and active region conditions. The importance of the self-reflection term in producing relatively short heating scale heights and thermal nonequilibrium cycles is also discussed.

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.

The light curves of solar coronal loops often peak first in channels associated with higher temperatures and then in those associated with lower. The time delays between the different narrowband EUV channels have been measured for many individual loops and recently for every pixel of an active region observation. Time delays between channels for an active region exhibit a wide range of values, with maxima $>$ 5,000\,s. These large time delays make up 3-26\% (depending on the channel pair) of the pixels where a significant, positive time delay is measured. It has been suggested that time delays can be explained by impulsive heating. In this paper, we investigate whether the largest observed time delays can be explained by this hypothesis by simulating a series of coronal loops with different heating rates, loop lengths, abundances, and geometries to determine the range of expected time delays between a set of four EUV channels. We find that impulsive heating cannot address the largest time delays observed in two of the channel pairs and that the majority of the large time delays can only be explained by long, expanding loops with photospheric abundances. Additional observations may rule out these simulations as an explanation for the long time delays. We suggest that either the time delays found in this manner may not be representative of real loop evolution, or that the impulsive heating and cooling scenario may be too simple to explain the observations and other heating scenarios must be explored.

Coronal jets are transient, collimated eruptions that occur in regions of predominantly open magnetic field in the solar corona. Our understanding of these events has greatly evolved in recent years but several open questions, such as the contribution of coronal jets to the solar wind, remain. Here we present an overview of the observations and numerical modeling of coronal jets, followed by a brief description of "next-generation" simulations that include an advanced description of the energy transfer in the corona ("thermodynamic MHD"), large spherical computational domains, and the solar wind. These new models will allow us to address some of the open questions.

There is a growing appreciation that the environmental conditions that we call space weather impact the technological infrastructure that powers the coupled economies around the world. With that comes the need to better shield society against space weather by improving forecasts, environmental specifications, and infrastructure design. [...] advanced understanding of space weather requires a coordinated international approach to effectively provide awareness of the processes within the Sun-Earth system through observation-driven models. This roadmap prioritizes the scientific focus areas and research infrastructure that are needed to significantly advance our understanding of space weather of all intensities and of its implications for society. Advancement of the existing system observatory through the addition of small to moderate state-of-the-art capabilities designed to fill observational gaps will enable significant advances. Such a strategy requires urgent action: key instrumentation needs to be sustained, and action needs to be taken before core capabilities are lost in the aging ensemble. We recommend advances through priority focus (1) on observation-based modeling throughout the Sun-Earth system, (2) on forecasts more than 12 hrs ahead of the magnetic structure of incoming coronal mass ejections, (3) on understanding the geospace response to variable solar-wind stresses that lead to intense geomagnetically-induced currents and ionospheric and radiation storms, and (4) on developing a comprehensive specification of space climate, including the characterization of extreme space storms to guide resilient and robust engineering of technological infrastructures. The roadmap clusters its implementation recommendations by formulating three action pathways, and outlines needed instrumentation and research programs and infrastructure for each of these. [...]

Although it is widely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process, the detailed mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still not fully understood. Cranmer et al. (2007) developed a sophisticated, 1D, time-steady model of the solar wind with turbulence dissipation. By varying the coronal magnetic field, they obtain, for a single choice of wave properties, a realistic range of slow and fast wind conditions with a sharp latitudinal transition between the two streams. Using a 1D, time-dependent model of the solar wind of Lionello et al. (2014), which incorporates turbulent dissipation of Alfv\'en waves to provide heating and acceleration of the plasma, we have explored a similar configuration, obtaining qualitatively equivalent results. However, our calculations suggest that the rapid transition between slow and fast wind suggested by this 1D model may be disrupted in multidimensional MHD simulations by the requirement of transverse force balance.

Force-free equilibria containing two vertically arranged magnetic flux ropes of like chirality and current direction are considered as a model for split filaments/prominences and filament-sigmoid systems. Such equilibria are constructed analytically through an extension of the methods developed in Titov & D\'emoulin (1999) and numerically through an evolutionary sequence including shear flows, flux emergence, and flux cancellation in the photospheric boundary. It is demonstrated that the analytical equilibria are stable if an external toroidal (shear) field component exceeding a threshold value is included. If this component decreases sufficiently, then both flux ropes turn unstable for conditions typical of solar active regions, with the lower rope typically being unstable first. Either both flux ropes erupt upward, or only the upper rope erupts while the lower rope reconnects with the ambient flux low in the corona and is destroyed. However, for shear field strengths staying somewhat above the threshold value, the configuration also admits evolutions which lead to partial eruptions with only the upper flux rope becoming unstable and the lower one remaining in place. This can be triggered by a transfer of flux and current from the lower to the upper rope, as suggested by the observations of a split filament in Paper~I (Liu et al. 2012). It can also result from tether-cutting reconnection with the ambient flux at the X-type structure between the flux ropes, which similarly influences their stability properties in opposite ways. This is demonstrated for the numerically constructed equilibrium.

Although the mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still being actively investigated, it is largely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process. \citet{2010ApJ...708L.116V} presented a model for heating and accelerating the solar wind based on the turbulent dissipation of Alfv\'en waves. We first use a time-dependent model of the solar wind to reproduce one of \citeauthor{2010ApJ...708L.116V}'s solutions; then we extend its application to the case when the energy equation includes thermal conduction and radiation losses, and the upper chromosphere is part of the computational domain. Using this model, we explore parameter space and describe the characteristics of a fast-solar-wind solution. We discuss how this formulation may be applied to a 3D MHD model of the coron a and solar wind \citep{2009ApJ...690..902L}.

Knowledge about the background solar wind plays a crucial role in the framework of space weather forecasting. In-situ measurements of the background solar wind are only available for a few points in the heliosphere where spacecraft are located, therefore we have to rely on heliospheric models to derive the distribution of solar wind parameters in interplanetary space. We test the performance of different solar wind models, namely Magnetohydrodynamic Algorithm outside a Sphere/ENLIL (MAS/ENLIL), Wang-Sheeley-Arge/ENLIL (WSA/ENLIL), and MAS/MAS, by comparing model results with in-situ measurements from spacecraft located at 1 AU distance to the Sun (ACE, Wind). To exclude the influence of interplanetary coronal mass ejections (ICMEs), we chose the year 2007 as a time period with low solar activity for our comparison. We found that the general structure of the background solar wind is well reproduced by all models. The best model results were obtained for the parameter solar wind speed. However, the predicted arrival times of high-speed solar wind streams have typical uncertainties of the order of about one day. Comparison of model runs with synoptic magnetic maps from different observatories revealed that the choice of the synoptic map significantly affects the model performance.

A sequence of apparently coupled eruptions was observed on 2010 August 1-2 by SDO and STEREO. The eruptions were closely synchronized with one another, even though some of them occurred at widely separated locations. In an attempt to identify a plausible reason for such synchronization, we study the large-scale structure of the background magnetic configuration. The coronal field was computed from the photospheric magnetic field observed at the appropriate time period by using the potential field source-surface model. We investigate the resulting field structure by analyzing the so-called squashing factor calculated at the photospheric and source-surface boundaries, as well as at different coronal cross-sections. Using this information as a guide, we determine the underlying structural skeleton of the configuration, including separatrix and quasi-separatrix surfaces. Our analysis reveals, in particular, several pseudo-streamers in the regions where the eruptions occurred. Of special interest to us are the magnetic null points and separators associated with the pseudo-streamers. We propose that magnetic reconnection triggered along these separators by the first eruption likely played a key role in establishing the assumed link between the sequential eruptions. The present work substantiates our recent simplified magnetohydrodynamic model of sympathetic eruptions and provides a guide for further deeper study of these phenomena. Several important implications of our results for the S-web model of the slow solar wind are also addressed.

We study an active-region dextral filament which was composed of two branches separated in height by about 13 Mm. This "double-decker" configuration sustained for days before the upper branch erupted with a GOES-class M1.0 flare on 2010 August 7. Analyzing this evolution, we obtain the following main results. 1) During hours before the eruption, filament threads within the lower branch were observed to intermittently brighten up, lift upward, and then merge with the upper branch. The merging process contributed magnetic flux and current to the upper branch, resulting in its quasi-static ascent. 2) This transfer might serve as the key mechanism for the upper branch to lose equilibrium by reaching the limiting flux that can be stably held down by the overlying field or by reaching the threshold of the torus instability. 3) The erupting branch first straightened from a reverse S shape that followed the polarity inversion line and then writhed into a forward S shape. This shows a transfer of left-handed helicity in a sequence of writhe-twist-writhe. The fact that the initial writhe is converted into the twist of the flux rope excludes the helical kink instability as the trigger process of the eruption, but supports the occurrence of the instability in the main phase, which is indeed indicated by the very strong writhing motion. 4) A hard X-ray sigmoid, likely of coronal origin, formed in the gap between the two original filament branches in the impulsive phase of the associated flare. This supports a model of transient sigmoids forming in the vertical flare current sheet. 5) Left-handed magnetic helicity is inferred for both branches of the dextral filament. 6) Two types of force-free magnetic configurations are compatible with the data, a double flux rope equilibrium and a single flux rope situated above a loop arcade.

Sympathetic eruptions on the Sun have been observed for several decades, but the mechanisms by which one eruption can trigger another one remain poorly understood. We present a 3D MHD simulation that suggests two possible magnetic trigger mechanisms for sympathetic eruptions. We consider a configuration that contains two coronal flux ropes located within a pseudo-streamer and one rope located next to it. A sequence of eruptions is initiated by triggering the eruption of the flux rope next to the streamer. The expansion of the rope leads to two consecutive reconnection events, each of which triggers the eruption of a flux rope by removing a sufficient amount of overlying flux. The simulation qualitatively reproduces important aspects of the global sympathetic event on 2010 August 1 and provides a scenario for so-called twin filament eruptions. The suggested mechanisms are applicable also for sympathetic eruptions occurring in other magnetic configurations.

The study concerns the streamer belt observed at high spectral resolution during the minimum of solar cycle 22 with the Ultraviolet Coronagraph Spectrometer (UVCS) onboard SOHO. On the basis of a spectroscopic analysis of the O VI doublet, the solar wind plasma parameters are inferred in the extended corona. The analysis accounts for the coronal magnetic topology, extrapolated through a 3D magneto-hydrodynamic model, in order to define the streamer boundary and to analyse the edges of coronal holes. The results of the analysis allow an accurate identification of the source regions of the slow coronal wind that are confirmed to be along the streamer boundary in the open magnetic field region.

Models for the origin of the slow solar wind must account for two seemingly contradictory observations: The slow wind has the composition of the closed field corona, implying that it originates from the continuous opening and closing of flux at the boundary between open and closed field. On the other hand, the slow wind also has large angular width, up to ~ 60{\circ}, suggesting that its source extends far from the open-closed boundary. We propose a model that can explain both observations. The key idea is that the source of the slow wind at the Sun is a network of narrow (possibly singular) open-field corridors that map to a web of separatrices and quasi-separatrix layers in the heliosphere. We compute analytically the topology of an open-field corridor and show that it produces a quasi-separatrix layer in the heliosphere that extends to angles far from the heliospheric current sheet. We then use an MHD code and MDI/SOHO observations of the photospheric magnetic field to calculate numerically, with high spatial resolution, the quasi-steady solar wind and magnetic field for a time period preceding the August 1, 2008 total solar eclipse. Our numerical results imply that, at least for this time period, a web of separatrices (which we term an S-web) forms with sufficient density and extent in the heliosphere to account for the observed properties of the slow wind. We discuss the implications of our S-web model for the structure and dynamics of the corona and heliosphere, and propose further tests of the model.

In recent work, Antiochos and coworkers argued that the boundary between the open and closed field regions on the Sun can be extremely complex with narrow corridors of open flux connecting seemingly disconnected coronal holes from the main polar holes, and that these corridors may be the sources of the slow solar wind. We examine, in detail, the topology of such magnetic configurations using an analytical source surface model that allows for analysis of the field with arbitrary resolution. Our analysis reveals three important new results: First, a coronal hole boundary can join stably to the separatrix boundary of a parasitic polarity region. Second, a single parasitic polarity region can produce multiple null points in the corona and, more important, separator lines connecting these points. It is known that such topologies are extremely favorable for magnetic reconnection, because they allow this process to occur over the entire length of the separators rather than being confined to a small region around the nulls. Finally, the coronal holes are not connected by an open-field corridor of finite width, but instead are linked by a singular line that coincides with the separatrix footprint of the parasitic polarity. We investigate how the topological features described above evolve in response to the motion of the parasitic polarity region. The implications of our results for the sources of the slow solar wind and for coronal and heliospheric observations are discussed.

We consider a three-dimensional bipolar force-free magnetic field with non zero magnetic helicity, occupying a half-space, and study the problem of its evolution driven by an imposed photospheric flux decrease. For this specific setting of the Flux Cancellation Model describing coronal mass ejections occuring in active regions, we address the issues of the physical meaning of flux decrease, of the influence on field evolution of the size of the domain over which this decrease is imposed, and of the existence of an energetic criterion characterizing the possible onset of disruption of the configuration. We show that: (1) The imposed flux disappearance can be interpreted in terms of transport of positive and negative fluxes towards the inversion line, where they get annihilated. (2) For the particular case actually computed, in which the initial state is quite sheared, the formation of a twisted flux rope and the subsequent global disruption of the configuration are obtained when the flux has decreased by only a modest amount over a limited part of the whole active region. (3) The disruption is produced when the magnetic energy becomes of the order of the decreasing energy of a semi-open field, and then before reaching the energy of the associated fully open field. This suggests that the mechanism leading to the disruption is nonequilibrium as in the case where flux is imposed to decrease over the whole region.

Low and mid-latitude coronal holes (CHs) observed on the Sun during the current solar activity minimum (from September 21, 2006, Carrington rotation (CR) 2048, until June 26, 2009 (CR 2084)) were analyzed using {\it SOHO}/EIT and STEREO-A SECCHI EUVI data. From both the observations and Potential Field Source Surface (PFSS) modeling, we find that the area occupied by CHs inside a belt of $\pm 40^\circ$ around the solar equator is larger in the current 2007 solar minimum relative to the similar phase of the previous 1996 solar minimum. The enhanced CH area is related to a recurrent appearance of five persistent CHs, which survived during 7-27 solar rotations. Three of the CHs are of positive magnetic polarity and two are negative. The most long-lived CH was being formed during 2 days and existed for 27 rotations. This CH was associated with fast solar wind at 1 AU of approximately 620$\pm 40$ km s$^{-1}$. The 3D MHD modeling for this time period shows an open field structure above this CH. We conclude that the global magnetic field of the Sun possessed a multi-pole structure during this time period. Calculation of the harmonic power spectrum of the solar magnetic field demonstrates a greater prevalence of multi-pole components over the dipole component in the 2007 solar minimum compared to the 1996 solar minimum. The unusual large separation between the dipole and multi-pole components is due to the very low magnitude of the dipole component, which is three times lower than that in the previous 1996 solar minimum.

A general method for describing magnetic reconnection in arbitrary three-dimensional magnetic configurations is proposed. The method is based on the field-line mapping technique previously used only for the analysis of magnetic structure at a given time. This technique is extended here so as to analyze the evolution of magnetic structure. Such a generalization is made with the help of new dimensionless quantities called "slip-squashing factors". Their large values define the surfaces that border the reconnected or to-be-reconnected magnetic flux tubes for a given period of time during the magnetic evolution. The proposed method is universal, since it assumes only that the time sequence of evolving magnetic field and the tangential boundary flows are known. The application of the method is illustrated for simple examples, one of which was considered previously by Hesse and coworkers in the framework of the general magnetic reconnection theory. The examples help us to compare these two approaches; they reveal also that, just as for magnetic null points, hyperbolic and cusp minimum points of a magnetic field may serve as favorable sites for magnetic reconnection. The new method admits a straightforward numerical implementation and provides a powerful tool for the diagnostics of magnetic reconnection in numerical models of solar-flare-like phenomena in space and laboratory plasmas.

A simple model of the coronal magnetic field prior to the CME eruption on May 12 1997 is developed. First, the magnetic field is constructed by superimposing a large-scale background field and a localized bipolar field to model the active region (AR) in the current-free approximation. Second, this potential configuration is quasi-statically sheared by photospheric vortex motions applied to two flux concentrations of the AR. Third, the resulting force-free field is then evolved by canceling the photospheric magnetic flux with the help of an appropriate tangential electric field applied to the central part of the AR. To understand the structure of the modeled configuration, we use the field line mapping technique by generalizing it to spherical geometry. It is demonstrated that the initial potential configuration contains a hyperbolic flux tube (HFT) which is a union of two intersecting quasi-separatrix layers. This HFT provides a partition of the closed magnetic flux between the AR and the global solar magnetic field. The vortex motions applied to the AR interlock the field lines in the coronal volume to form additionally two new HFTs pinched into thin current layers. Reconnection in these current layers helps to redistribute the magnetic flux and current within the AR in the flux-cancellation phase. In this phase, a magnetic flux rope is formed together with a bald patch separatrix surface wrapping around the rope. Other important implications of the identified structural features of the modeled configuration are also discussed.

We use direct numerical simulations to study the evolution, or relaxation, of magnetic configurations to an equilibrium state. We use the full single-fluid equations of motion for a magnetized, non-resistive, but viscous fluid; and a Lagrangian approach is used to obtain exact solutions for the magnetic field. As a result, the topology of the magnetic field remains unchanged, which makes it possible to study the case of topological nonequilibrium. We find two cases for which such nonequilibrium appears, indicating that these configurations may develop singular current sheets.