Search

Observations have shown a clear association of filament/prominence eruptions with the emergence of magnetic flux in or near filament channels. Magnetohydrodynamic (MHD) simulations have been employed to systematically study the conditions under which such eruptions occur. These simulations to date have modeled filament channels as two-dimensional (2D) flux ropes or 3D uniformly sheared arcades. Here we present MHD simulations of flux emergence into a more realistic configuration consisting of a bipolar active region containing a line-tied 3D flux rope. We use the coronal flux-rope model of Titov et al. (2014) as the initial condition and drive our simulations by imposing boundary conditions extracted from a flux-emergence simulation by Leake et al. (2013). We identify three mechanisms that determine the evolution of the system: (i) reconnection displacing foot points of field lines overlying the coronal flux rope, (ii) changes of the ambient field due to the intrusion of new flux at the boundary, and (iii) interaction of the (axial) electric currents in the pre-existing and newly emerging flux systems. The relative contributions and effects of these mechanisms depend on the properties of the pre-existing and emerging flux systems. Here we focus on the location and orientation of the emerging flux relative to the coronal flux rope. Varying these parameters, we investigate under which conditions an eruption of the latter is triggered., Comment: 18 pages, 8 figures, accepted for publication by The Astrophysical Journal

The so-called regularized Biot-Savart laws (RBSLs) provide an efficient and flexible method for modeling pre-eruptive magnetic configurations of coronal mass ejections (CMEs) whose characteristics are constrained by observational images and magnetic-field data. This method allows one to calculate the field of magnetic flux ropes (MFRs) with small circular cross-sections and an arbitrary axis shape. The field of the whole configuration is constructed as a superposition of (1) such a flux-rope field and (2) an ambient potential field derived, for example, from an observed magnetogram. The RBSL kernels are determined from the requirement that the MFR field for a straight cylinder must be exactly force-free. For a curved MFR, however, the magnetic forces are generally unbalanced over the whole path of the MFR. To minimize these forces, we apply a modified Gauss-Newton method to find optimal MFR parameters. This is done by iteratively adjusting the MFR axis path and axial current. We then try to relax the resulting optimized configuration in a subsequent line-tied zero-beta magnetohydrodynamic simulation toward a force-free equilibrium. By considering two models of the sigmoidal pre-eruption configuration for the 2009 February 13 CME, we demonstrate how this approach works and what it is capable of. We show, in particular, that the building blocks of the core magnetic structure described by these models match to morphological features typically observed in such type of configurations. Our method will be useful for both the modeling of particular eruptive events and theoretical studies of idealized pre-eruptive MFR configurations., Comment: 21 pages, 11 figures, accepted to ApJS

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., Comment: To appear in Ap.J

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.

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.

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., Comment: 12 pages, 4 figures, to appear in ASP proceedings

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., Comment: 32 pages, 10 figures, accepted to ApJ

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., Comment: 25 pages, 11 figures, to appear in ApJ 2008

The precise nature of the heating mechanism (location, duration) in coronal loops is still a matter of enormous research. We present results from a one-dimensional (1D) hydrodynamic loop simulation of a coronal loop which was run using different parameters such as loops length (50, 200 and 500 light-seconds), maximum temperature reached (3 million degrees Kelvin and 10 million degrees Kelvin), and abundances. For each scenario the model outputs were used to calculate the corresponding lightcurves as seen by X-ray telescope/Be-thin filter and various Extreme Ultra Violet Atmospheric Imaging Assembly channels. The lag time between the peak of these lightcurves was computed using cross-correlation and plotted as a function of loop length. Additional results were computed using the zero-dimensional Enthalpy-Based Thermal Evolution of Loops (EBTEL) code in order to test the compatibility of the two codes and to investigate additional loop lengths. Initial results indicate that the long (greater than 5000 seconds) lags observed in the approximately 100 light-seconds loops of active regions can only be reproduced using photospheric abundances and much longer loop lengths. This result suggests that the observed time lags cannot be completely explained by impulsive heating.

Time correlation analysis has been used to show widespread cooling in the solar corona; this cooling has been interpreted as a result of impulsive (nanoflare) heating. In this work, we investigate wide-spread cooling using a 3D model for a solar active region which has been heated with highly stratified heating. This type of heating drives thermal non-equilibrium solutions, meaning that though the heating is effectively steady, the density and temperature in the solution are not. We simulate the expected observations in narrowband EUV images and apply the time correlation analysis. We find that the results of this analysis are qualitatively similar to the observed data. We discuss additional diagnostics that may be applied to differentiate between these two heating scenarios.

Many different techniques have been used to characterize the plasma in the solar corona: density-sensitive spectral line ratios are used to infer the density, the evolution of coronal structures in different passbands is used to infer the temperature evolution, and the simultaneous intensities measured in multiple passbands are used to determine the emission measure. All these analysis techniques assume that the intensity of the structures can be isolated through background subtraction. In this paper, we use simulated observations from a 3D hydrodynamic simulation of a coronal active region to verify these diagnostics. The density and temperature from the simulation are used to generate images in several passbands and spectral lines. We identify loop structures in the simulated images and calculate the loop background. We then determine the density, temperature and emission measure distribution as a function of time from the observations and compare with the true temperature and density of the loop. We find that the overall characteristics of the temperature, density, and emission measure are recovered by the analysis methods, but the details of the true temperature and density are not. For instance, the emission measure curves calculated from the simulated observations are much broader than the true emission measure distribution, though the average temperature evolution is similar. These differences are due, in part, to inadequate background subtraction, but also indicate a limitation of the analysis methods.

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 ux 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 eld 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. Such topologies are extremely favorable for magnetic reconnection, because it can now occur over the entire length of the separators rather than being con ned to a small region around the nulls. Finally, the coronal holes are not connected by an open- eld 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 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.

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 has large angular width, up to approximately 60 degrees, 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 front the heliospheric current sheet. We then use an MHD code and MIDI/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.

A simple model of the coronal magnetic field prior to the coronal mass ejection (CME) eruption on 1997 May 12 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. We demonstrate 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. Such a partition is approximate since the entire flux distribution is perfectly continuous. The vortex motions applied to the AR inter-lock 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 describe the use of a three-dimensional MHD model to predict the s tructure of the corona prior to the total solar eclipse of March 29, 2006. The calculation uses the observed photospheric radial magnetic f ield as a boundary condition. We use a new version of our model that has an improved description of energy transport in the corona. The mo del allows us to predict the emission of X-ray and EUV radiation in t he corona. We compare the predicted polarization brightness in the co rona with four observations of the eclipse from Greece, Egypt, and Li bya, and we demonstrate that the model accurately predicts the largescale structure of the corona. We also compare X-ray emission from the model with GOES/SXI images.

The Sentinels mission promises to provide a unique view of the acceleration and transport processes of energetic particles as well as the initiation and evolution of coronal mass ejections (CMEs). A crucial component in understanding the physics associated with these processes lies in the large-scale structure of the corona and heliosphere, particularly during the eruption and propagation of fast CMEs. In this talk we review the current status of our MHD modeling efforts, focusing on Sentinels-specific science, and project forward to envisage what capabilities we may have developed by the time that the Sentinels launch (2012). In conjunction with Solar Orbiter, the farside and near-Earth Sentinels spacecraft will provide simultaneous photospheric magnetograph measurements at multiple longitudes, which will lead to major improvements in our ability to prescribe accurate, time-dependent global boundary conditions. Data returned from the inner heliospheric Sentinels will be used to validate these results. The model results can be used in a variety of ways to interpret the observations. For example, products from the model, such as the properties of CME-driven shocks and CME-associated reconnection sites can be used to interpret complex energetic particle profiles. Also, by tracing along magnetic field lines, the inferred sites of the energetic particles can be connected directly with the in situ measurements at each spacecraft. Perhaps more so than any previous mission, sophisticated models will be required to unravel the broad and disparate measurements returned by the suite of Sentinels spacecraft.