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.

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.

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.

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.

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.

Using solar vector magnetograms of the highest available spatial resolution and signal-to-noise ratio we perform a detailed study of electric current patterns in two solar active regions: a flaring/eruptive, and a flare-quiet one. We aim to determine whether active regions inject non-neutralized (net) electric currents in the solar atmosphere, responding to a debate initiated nearly two decades ago that remains inconclusive. We find that well-formed, intense magnetic polarity inversion lines (PILs) within active regions are the only photospheric magnetic structures that support significant net current. More intense PILs seem to imply stronger non-neutralized current patterns per polarity. This finding revises previous works that claim frequent injections of intense non-neutralized currents by most active regions appearing in the solar disk but also works that altogether rule out injection of non-neutralized currents. In agreement with previous studies, we also find that magnetically isolated active regions remain globally current-balanced. In addition, we confirm and quantify the preference of a given magnetic polarity to follow a given sense of electric currents, indicating a dominant sense of twist in active regions. This coherence effect is more pronounced in more compact active regions with stronger PILs and must be of sub-photospheric origin. Our results yield a natural explanation of the Lorentz force, invariably generating velocity and magnetic shear along strong PILs, thus setting a physical context for the observed pre-eruption evolution in solar active regions.

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.

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 coming decade will see the routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness. To capitalize on the new (or soon to be available) facilities such as SDO, ATST and FASR, and the challenges they present in the visualization and synthesis of multi-wavelength datasets, we propose that realistic, sophisticated, 3D active region and flare modeling is timely and critical, and will be a forefront of coronal studies over the coming decade. To make such modeling a reality, a broad, concerted effort is needed to capture the wealth of information resulting from the data, develop a synergistic modeling effort, and generate the necessary visualization, interpretation and model-data comparison tools to accurately extract the key physics.

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.

Many stormy events in astrophysics occur due to the sudden magnetic energy release. This is possible if a magnetic configuration abruptly changes its topology, an event usually referred to as magnetic reconnection. It is known that pure Ohmic decay is inefficient, occurring during cosmological times (due to the huge characteristic scales $L$). It is recognized that the presence of current sheets speeds up the process, but still insufficiently$^{1,2,3,4,5}$. We show that, in highly compressible and substantially gravitational media, the reconnection is fast enough to account for stormy events. Thus, highly compressible situations offer exiting opportunities in explanations of violent events, although full-scale compressible and gravitational simulations proved to be quite challenging.

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.

The solar magnetic field is the primary agent that drives solar activity and couples the Sun to the Heliosphere. Although the details of this coupling depend on the quantitative properties of the field, many important aspects of the corona - solar wind connection can be understood by considering only the general topological properties of those regions on the Sun where the field extends from the photosphere out to interplanetary space, the so-called open field regions that are usually observed as coronal holes. From the simple assumptions that underlie the standard quasi-steady corona-wind theoretical models, and that are likely to hold for the Sun, as well, we derive two conjectures on the possible structure and dynamics of coronal holes: (1) Coronal holes are unique in that every unipolar region on the photosphere can contain at most one coronal hole. (2) Coronal holes of nested polarity regions must themselves be nested. Magnetic reconnection plays the central role in enforcing these constraints on the field topology. From these conjectures we derive additional properties for the topology of open field regions, and propose several observational predictions for both the slowly varying and transient corona/solar wind.