Showing 1–36 of 36 results for author: Linton, M G

Assessing the Performance of the ADAPT and AFT Flux Transport Models Using In-Situ Measurements From Multiple Satellites

Authors: Kalman J. Knizhnik, Micah J. Weberg, Elena Provornikova, Harry P. Warren, Mark G. Linton, Shaheda Begum Shaik, Yuan-Kuen Ko, Samuel J. Schonfeld, Ignacio Ugarte-Urra, Lisa A. Upton

Abstract: The launches of Parker Solar Probe (Parker) and Solar Orbiter (SolO) are enabling a new era of solar wind studies that track the solar wind from its origin at the photosphere, through the corona, to multiple vantage points in the inner heliosphere. A key ingredient for these models is the input photospheric magnetic field map that provides the boundary condition for the coronal portion of many hel… ▽ More The launches of Parker Solar Probe (Parker) and Solar Orbiter (SolO) are enabling a new era of solar wind studies that track the solar wind from its origin at the photosphere, through the corona, to multiple vantage points in the inner heliosphere. A key ingredient for these models is the input photospheric magnetic field map that provides the boundary condition for the coronal portion of many heliospheric models. In this paper, we perform steady-state, data-driven magnetohydrodynamic (MHD) simulations of the solar wind during Carrington rotation 2258 with the GAMERA model. We use the ADAPT and AFT flux transport models and quantitatively assess how well each model matches in-situ measurements from Parker, SolO, and Earth. We find that both models reproduce the magnetic field components at Parker quantitatively well. At SolO and Earth, the magnetic field is reproduced relatively well, though not as well as at Parker, and the density is reproduced extremely poorly. The velocity is overpredicted at Parker, but not at SolO or Earth, hinting that the Wang-Sheeley-Arge (WSA) relation, fine-tuned for Earth, misses the deceleration of the solar wind near the Sun. We conclude that AFT performs quantitatively similarly to ADAPT in all cases and that both models are comparable to a purely WSA heliospheric treatment with no MHD component. Finally, we trace field lines from SolO back to an active region outflow that was observed by Hinode/EIS, and which shows evidence of elevated charge state ratios. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Solar Eruptions Triggered by Flux Emergence Below or Near a Coronal Flux Rope

Authors: T. Török, M. G. Linton, J. E. Leake, Z. Mikić, R. Lionello, V. S. Titov, C. Downs

Abstract: 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. He… ▽ More 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. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 18 pages, 8 figures, accepted for publication by The Astrophysical Journal

Simulating the photospheric to coronal plasma using magnetohydrodyanamic characteristics I: data-driven boundary conditions

Authors: Lucas A. Tarr, N. Dylan Kee, Mark G. Linton, Peter W. Schuck, James E. Leake

Abstract: We develop a general description of how information propagates through a magnetohydrodynamic (MHD) system based on the method of characteristics and use that to formulate numerical boundary conditions (BCs) that are intrinsically consistent with the MHD equations. Our formulation includes two major advances for simulations of the Sun. First, we derive data-driven BCs that optimally match the state… ▽ More We develop a general description of how information propagates through a magnetohydrodynamic (MHD) system based on the method of characteristics and use that to formulate numerical boundary conditions (BCs) that are intrinsically consistent with the MHD equations. Our formulation includes two major advances for simulations of the Sun. First, we derive data-driven BCs that optimally match the state of the plasma inferred from a time series of observations of a boundary (e.g., the solar photosphere). Second, our method directly handles random noise and systematic bias in the observations, and finds a solution for the boundary evolution that is strictly consistent with MHD and maximally consistent with the observations. We validate the method against a Ground Truth (GT) simulation of an expanding spheromak. The data-driven simulation can reproduce the GT simulation above the photosphere with high fidelity when driven at high cadence. Errors progressively increase for lower driving cadence until a threshold cadence is reached and the driven simulation can no longer accurately reproduce the GT simulation. However, our characteristic formulation of the BCs still requires adherence of the boundary evolution to the MHD equations even when the driven solution departs from the true solution in the driving layer. That increasing departure clearly indicates when additional information at the boundary is needed to fully specify the correct evolution of the system. The method functions even when no information about the evolution of some variables on the lower boundary is available, albeit with a further decrease in fidelity. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: Submitted for publication to the Astrophysical Journal; 77 pages, 12 figures

Deflection of Coronal Mass Ejections in Unipolar Ambient Magnetic Fields

Authors: Michal Ben-Nun, Tibor Török, Erika Palmerio, Cooper Downs, Viacheslav S. Titov, Mark G. Linton, Ronald M. Caplan, Roberto Lionello

Abstract: The trajectories of coronal mass ejections (CMEs) are often seen to substantially deviate from a purely radial propagation direction. Such deviations occur predominantly in the corona and have been attributed to "channeling" or deflection of the eruptive flux by asymmetric ambient magnetic fields. Here, we investigate an additional mechanism that does not require any asymmetry of the pre-eruptive… ▽ More The trajectories of coronal mass ejections (CMEs) are often seen to substantially deviate from a purely radial propagation direction. Such deviations occur predominantly in the corona and have been attributed to "channeling" or deflection of the eruptive flux by asymmetric ambient magnetic fields. Here, we investigate an additional mechanism that does not require any asymmetry of the pre-eruptive ambient field. Using magnetohydrodynamic numerical simulations, we show that the trajectory of CMEs through the solar corona can significantly deviate from a radial direction when propagation takes place in a unipolar radial field. We demonstrate that the deviation is most prominent below ~15 solar radii and can be attributed to an "effective IxB force" that arises from the intrusion of a magnetic flux rope with a net axial electric current into a unipolar background field. These results are important for predictions of CME trajectories in the context of space weather forecasts, as well as for reaching a deeper understanding of the fundamental physics underlying CME interactions with the ambient fields in the extended solar corona. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 15 pages, 12 figures

Disentangling the Entangled Linkages of Relative Magnetic Helicity

Authors: Peter W. Schuck, Mark G. Linton

Abstract: Magnetic helicity, $H$, measures magnetic linkages in a volume. The early theoretical development of helicity focused on magnetically closed systems in $\mathcal{V}$ bounded by $\mathcal{S}$. For magnetically closed systems, $\mathcal{V}\in\mathbb{R}^3=\mathcal{V}+\mathcal{V}^*$, no magnetic flux threads the boundary, $\hat{\boldsymbol{n}}\cdot\boldsymbol{B}|_\mathcal{S}=0$. Berger and Field (1984… ▽ More Magnetic helicity, $H$, measures magnetic linkages in a volume. The early theoretical development of helicity focused on magnetically closed systems in $\mathcal{V}$ bounded by $\mathcal{S}$. For magnetically closed systems, $\mathcal{V}\in\mathbb{R}^3=\mathcal{V}+\mathcal{V}^*$, no magnetic flux threads the boundary, $\hat{\boldsymbol{n}}\cdot\boldsymbol{B}|_\mathcal{S}=0$. Berger and Field (1984) and Finn and Antonsen (1985) extended the definition of helicity to relative helicity, $\mathcal{H}$, for magnetically open systems where magnetic flux may thread the boundary. Berger (1999,2003) expressed this relative helicity as two gauge invariant terms that describe the self helicity of magnetic field that closes inside $\mathcal{V}$ and the mutual helicity between the magnetic field that threads the boundary $\mathcal{S}$ and the magnetic field that closes inside $\mathcal{V}$. The total magnetic field that permeates $\mathcal{V}$ entangles magnetic fields that are produced by current sources $\boldsymbol{J}$ in $\mathcal{V}$ with magnetic fields that are produced by current sources $\boldsymbol{J}^*$ in $\mathcal{V}^*$. Building on this fact, we extend Berger's expressions for relative magnetic helicity to eight gauge invariant quantities that simultaneously characterize both of these self and mutual helicities and attribute their origins to currents $\boldsymbol{J}$ in $\mathcal{V}$ and/or $\boldsymbol{J}^*$ in $\mathcal{V}^*$, thereby disentangling the domain of origin for these entangled linkages. We arrange these eight terms into novel expressions for internal and external helicity (self) and internal-external helicity (mutual) based on their domain of origin. The implications of these linkages for interpreting magnetic energy is discussed and new boundary observables are proposed for tracking the evolution of the field that threads the boundary. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 27 pages, 2 figures

Characterizing the Umbral Magnetic Knots of $δでるた$-Sunspots

Authors: A. A. Norton, P. J. Levens, K. J. Knizhnik, M. G. Linton, Y. Liu

Abstract: Delta ($δでるた$)-spots are active regions (ARs) in which positive and negative umbrae share a penumbra. They are known to be the source of strong flares. We introduce a new quantity, the degree of $δでるた$ (Do$δでるた$), to measure the fraction of umbral flux participating in the $δでるた$-configuration and to isolate the dynamics of the magnetic knot, i.e. adjacent umbrae in the $δでるた$-configuration. Using Helioseismic a… ▽ More Delta ($δでるた$)-spots are active regions (ARs) in which positive and negative umbrae share a penumbra. They are known to be the source of strong flares. We introduce a new quantity, the degree of $δでるた$ (Do$δでるた$), to measure the fraction of umbral flux participating in the $δでるた$-configuration and to isolate the dynamics of the magnetic knot, i.e. adjacent umbrae in the $δでるた$-configuration. Using Helioseismic and Magnetic Imager data, we analyze 19 $δでるた$-spots and 11 $βべーた$-spots in detail, and 120 $δでるた$-spots in less detail. We find that $δでるた$-regions are not in a $δでるた$-configuration for the entire time but spend 55$\%$ of their observed time as $δでるた$-spots with an average, maximum Do$δでるた$ of 72$\%$. Compared to $βべーた$-spots, $δでるた$-spots have 2.6$\times$ the maximum umbral flux, 1.9$\times$ the flux emergence rate, 2.6$\times$ the rotation, and 72$\times$ the flare energy. On average, the magnetic knots rotate 17$^{\circ}$ day$^{-1}$ while the $βべーた$-spots rotate 2$^{\circ}$ day$^{-1}$. Approximately 72$\%$ of the magnetic knots present anti-Hale or anti-Joy tilts, contrasting starkly with only 9$\%$ of the $βべーた$-spots. A positive correlation exists between $φふぁい_{Doδでるた}$ and the flare energy emitted by that region. The $δでるた$-spots obey the hemispheric current helicity rule 64$\%$ of the time. 84$\%$ of the $δでるた$-spots are formed by single flux emergence events and 58$\%$ have a quadrupolar magnetic configuration. The $δでるた$-spot characteristics are consistent with the formation mechanism signatures as follows: 42$\%$ with the kink instability or Sigma effect, 32$\%$ with multi-segment buoyancy, 16$\%$ with collisions and two active regions that are unclassified but consistent with a rising O-ring. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: Accepted by the Astrophysical Journal. 28 pages, 12 figures

Continued PSP/WISPR Observations of a Phaethon-related Dust Trail

Authors: Karl Battams, Angel J. Gutarra-Leon, Brendan M. Gallagher, Matthew M. Knight, Guillermo Stenborg, Sarah Tanner, Mark G. Linton, Jamey R. Szalay, Michael S. P. Kelley, Russell A. Howard

Abstract: We present an update to the first white-light detections of a dust trail observed closely following the orbit of asteroid (3200) Phaethon, as seen by the Wide-field Imager for Parker Solar Probe (WISPR) instrument on the NASA Parker Solar Probe (PSP) mission. Here we provide a summary and analysis of observations of the dust trail over nine separate mission encounters between October 2018 and Augu… ▽ More We present an update to the first white-light detections of a dust trail observed closely following the orbit of asteroid (3200) Phaethon, as seen by the Wide-field Imager for Parker Solar Probe (WISPR) instrument on the NASA Parker Solar Probe (PSP) mission. Here we provide a summary and analysis of observations of the dust trail over nine separate mission encounters between October 2018 and August 2021 that saw the spacecraft approach to within 0.0277 au of the orbit of Phaethon. We find the photometric and estimated dust mass properties to be inline with those in the initial publication, with a visual (V) magnitude of approximately 16.1$\pm$0.3 per pixel, corresponding to a surface brightness of 26.1 mag arcsec$^{-2}$, and an estimated mass of dust within the range $10^{10}$ kg - $10^{12}$ kg depending on the assumed dust properties. However, the key finding of this survey is the discovery that the dust trail does not perfectly follow the orbit of Phaethon, with a clear separation noted between them that increases as a function of true anomaly, though the trail may differ from Phaethon's orbit by as little as 1-degree in periapsis. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 25 pages, 14 figures; Accepted for publication in Astrophysical Journal (Jul 25, 2022)

Overview of the remote sensing observations from PSP solar encounter 10 with perihelion at 13.3 Rsun

Authors: Russell A. Howard, Guillermo Stenborg, Angelos Vourlidas, Brendan M. Gallagher, Mark G. Linton, Phillip Hess, Nathan B. Rich, Paulett C. Liewer

Abstract: The closest perihelion pass of Parker Solar Probe (PSP), so far, occurred between 16 and 26 of November 2021 and reached ~13.29 Rsun from Sun center. This pass resulted in very unique observations of the solar corona by the Wide-field Instrument for Solar PRobe (WISPR). WISPR observed at least ten CMEs, some of which were so close that the structures appear distorted. All of the CMEs appeared to h… ▽ More The closest perihelion pass of Parker Solar Probe (PSP), so far, occurred between 16 and 26 of November 2021 and reached ~13.29 Rsun from Sun center. This pass resulted in very unique observations of the solar corona by the Wide-field Instrument for Solar PRobe (WISPR). WISPR observed at least ten CMEs, some of which were so close that the structures appear distorted. All of the CMEs appeared to have a magnetic flux rope (MFR) structure and most were oriented such that the view was along the axis orientation, revealing very complex interiors. Two CMEs had a small MFR develop in the interior, with a bright circular boundary surrounding a very dark interior. Trailing the larger CMEs were substantial outflows of small blobs and flux-rope like structures within striated ribbons, lasting for many hours. When the heliophysics plasma sheet (HPS) was inclined, as it was during the days around perihelion on November 21, 2021, the outflow was over a very wide latitude range. One CME was overtaken by a faster one, with a resultant compression of the rear of the leading CME and an unusual expansion in the trailing CME. The small Thomson Surface creates brightness variations of structures as they pass through the field of view. In addition to this dynamic activity, a brightness band from excess dust along the orbit of asteroid/comet 3200 Phaethon is also seen for several days. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 28 pages, 15 figures, accepted to ApJ

The Dynamic Evolution of Solar Wind Streams Following Interchange Reconnection

Authors: Roger B. Scott, Stephen J. Bradshaw, Mark G. Linton

Abstract: Interchange reconnection is thought to play an important role in determining the dynamics and material composition of the slow solar wind that originates from near coronal hole boundaries. To explore the implications of this process we simulate the dynamic evolution of a solar wind stream along a newly-opened magnetic flux tube. The initial condition is composed of a piecewise continuous dynamic e… ▽ More Interchange reconnection is thought to play an important role in determining the dynamics and material composition of the slow solar wind that originates from near coronal hole boundaries. To explore the implications of this process we simulate the dynamic evolution of a solar wind stream along a newly-opened magnetic flux tube. The initial condition is composed of a piecewise continuous dynamic equilibrium in which the regions above and below the reconnection site are extracted from steady-state solutions along open and closed field lines. The initial discontinuity at the reconnection site is highly unstable and evolves as a Riemann problem, decomposing into an outward-propagating shock and inward-propagating rarefaction that eventually develop into a classic N-wave configuration. This configuration ultimately propagates into the heliosphere as a coherent structure and the entire system eventually settles to a quasi-steady wind solution. In addition to simulating the fluid evolution we also calculate the time-dependent non-equilibrium ionization of oxygen in real time in order to construct in situ diagnostics of the conditions near the reconnection site. This idealized description of the plasma dynamics along a newly-opened magnetic field line provides a baseline for predicting and interpreting the implications of interchange reconnection for the slow solar wind. Notably, the density and velocity within the expanding N-wave are generally enhanced over the ambient wind, as is the O7+/O6+ ionization ratio, which exhibits a discontinuity across the reconnection site that is transported by the flow and arrives later than the propagating N-wave. △ Less

Submitted 13 June, 2022; originally announced June 2022.

The Rise and Emergence of Untwisted Toroidal Flux Ropes on the Sun

Authors: Kalman J. Knizhnik, James E. Leake, Mark G. Linton, Sally Dacie

Abstract: Magnetic flux ropes (MFRs) rising buoyantly through the Sun's convection zone are thought to be subject to viscous forces preventing them from rising coherently. Numerous studies have suggested that MFRs require a minimum twist in order to remain coherent during their rise. Furthermore, even MFRs that get to the photosphere may be unable to successfully emerge into the corona unless they are at le… ▽ More Magnetic flux ropes (MFRs) rising buoyantly through the Sun's convection zone are thought to be subject to viscous forces preventing them from rising coherently. Numerous studies have suggested that MFRs require a minimum twist in order to remain coherent during their rise. Furthermore, even MFRs that get to the photosphere may be unable to successfully emerge into the corona unless they are at least moderately twisted, since the magnetic pressure gradient needs to overcome the weight of the photospheric plasma. To date, however, no lower limit has been placed on the critical minimum twist required for an MFR to rise coherently through the convection zone or emerge through the photosphere. In this paper, we simulate an untwisted toroidal MFR which is able to rise from the convection zone and emerge through the photosphere as an active region that resembles those observed on the Sun. We show that untwisted MFRs can remain coherent during their rise and then pile-up near the photosphere, triggering the undular instability, allowing the MFR to emerge through the photosphere. We propose that the toroidal geometry of our MFR is critical for its coherent rise. Upon emerging, a pair of lobes rises into the corona which interact and reconnect, resulting in a localized high speed jet. The resulting photospheric magnetogram displays the characteristic salt-and-pepper structure often seen in observations. Our major result is that MFRs need not be twisted to rise coherently through the convection zone and emerge through the photosphere. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: Accepted to ApJ

Analysis of the Helical Kink Stability of Differently Twisted Magnetic Flux Ropes

Authors: Marta Florido-Llinas, Teresa Nieves-Chinchilla, Mark G. Linton

Abstract: Magnetic flux ropes (MFRs) are usually considered to be the magnetic structure that dominates the transport of helicity from the Sun into the heliosphere. They entrain a confined plasma within a helically organized magnetic structure and are able to cause geomagnetic activity. The formation, evolution and twist distribution of MFRs are issues subject to strong debate. Although different twist prof… ▽ More Magnetic flux ropes (MFRs) are usually considered to be the magnetic structure that dominates the transport of helicity from the Sun into the heliosphere. They entrain a confined plasma within a helically organized magnetic structure and are able to cause geomagnetic activity. The formation, evolution and twist distribution of MFRs are issues subject to strong debate. Although different twist profiles have been suggested so far, none of them has been thoroughly explored yet. The aim of this work is to present a theoretical study of the conditions under which MFRs with different twist profiles are kink stable and thereby shed some light on the aforementioned aspects. The magnetic field is modeled according to the circular-cylindrical analytical flux rope model in Nieves-Chinchilla et al. (Astrophys. J. 823, 27, 2016) as well as the Lundquist and Gold-Hoyle models, and the kink stability is analyzed with a numerical method that has been developed based on Linton, Longcope, and Fisher (Astrophys. J. 469, 954, 1996). The results are discussed in relation to MFR rotations, magnetic forces, the reversed chirality scenario, and the expansion throughout the heliosphere, among others, providing a theoretical background to improve the current understanding of the internal magnetic configuration of coronal mass ejections (CMEs). The data obtained by new missions like Parker Solar Probe or Solar Orbiter will give the opportunity to explore these results and ideas by observing MFRs closer than ever to the Sun. △ Less

Submitted 11 August, 2020; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 26 pages, 5 figures, 3 tables, accepted for publication in Solar Physics; typos and references corrected

Journal ref: Sol. Phys. 295, 118 (2020)

The Role of Twist in Kinked Flux Rope Emergence and Delta-Spot Formation

Authors: Kalman J. Knizhnik, Mark G. Linton, C. Richard DeVore

Abstract: It has been observationally well established that the magnetic configurations most favorable for producing energetic flaring events reside in delta-spots, a class of sunspots defined as having opposite polarity umbrae sharing a common penumbra. They are frequently characterized by extreme compactness, strong rotation and anti-Hale orientation. Numerous studies have shown that nearly all of the lar… ▽ More It has been observationally well established that the magnetic configurations most favorable for producing energetic flaring events reside in delta-spots, a class of sunspots defined as having opposite polarity umbrae sharing a common penumbra. They are frequently characterized by extreme compactness, strong rotation and anti-Hale orientation. Numerous studies have shown that nearly all of the largest solar flares originate in delta-spots, making the understanding of these structures a fundamental step in predicting space weather. Despite their important influence on the space environment, surprisingly little is understood about the origin and behavior of delta-spots. In this paper, we perform a systematic study of the behavior of emerging flux ropes to test a theoretical model for the formation of delta-spots: the kink instability of emerging flux ropes. We simulated the emergence of highly twisted, kink-unstable flux ropes from the convection zone into the corona, and compared their photospheric properties to those of emerged weakly twisted, kink-stable flux ropes. We show that the photospheric manifestations of the emergence of highly twisted flux ropes closely match the observed properties of delta-spots, and we discuss the resulting implications for observations. Our results strongly support and extend previous theoretical work that suggested that the kink instability of emerging flux ropes is a promising candidate to explain delta-spot formation, as it reproduces their key characteristics very well. △ Less

Submitted 16 August, 2018; originally announced August 2018.

Comments: Accepted for publication in ApJ, but comments welcome

The Grad-Shafranov Reconstruction of Toroidal Magnetic Flux Ropes: First Applications

Authors: Qiang Hu, M. G. Linton, B. E. Wood, P. Riley, T. Nieves-Chinchilla

Abstract: This article completes and extends a recent study of the Grad-Shafranov (GS) reconstruction in toroidal geometry, as applied to a two and a half dimensional configurations in space plasmas with rotational symmetry. A further application to the benchmark study of an analytic solution to the toroidal GS equation with added noise shows deviations in the reconstructed geometry of the flux rope configu… ▽ More This article completes and extends a recent study of the Grad-Shafranov (GS) reconstruction in toroidal geometry, as applied to a two and a half dimensional configurations in space plasmas with rotational symmetry. A further application to the benchmark study of an analytic solution to the toroidal GS equation with added noise shows deviations in the reconstructed geometry of the flux rope configuration, characterized by the orientation of the rotation axis, the major radius, and the impact parameter. On the other hand, the physical properties of the flux rope, including the axial field strength, and the toroidal and poloidal magnetic flux, agree between the numerical and exact GS solutions. We also present a real event study of a magnetic cloud flux rope from \textit{in situ} spacecraft measurements. The devised procedures for toroidal GS reconstruction are successfully executed. Various geometrical and physical parameters are obtained with associated uncertainty estimates. The overall configuration of the flux rope from the GS reconstruction is compared with the corresponding morphological reconstruction based on white-light images. The results show overall consistency, but also discrepancy in that the inclination angle of the flux rope central axis with respect to the ecliptic plane differs by about 20-30 degrees in the plane of the sky. We also compare the results with the original straight-cylinder GS reconstruction and discuss our findings. △ Less

Submitted 28 July, 2017; originally announced July 2017.

Comments: submitted to Sol. Phys.; under revision

Magnetic Flux Emergence and Decay Rates for Preceder and Follower Sunspots Observed with HMI

Authors: A. A. Norton, E. H. Jones, M. G. Linton, J. E. Leake

Abstract: We quantify the emergence and decay rates of preceder (p) and follower (f) sunspots within ten active regions from 2010-2014 using Space-weather Helioseismic Magnetic Imager Active Region Patch data. The sunspots are small- to mid-sized regions and contain a signed flux within a single polarity sunspot of $(1.1-6.5)\times 10^{21}$ Mx. The net unsigned flux within the regions, including plage, rang… ▽ More We quantify the emergence and decay rates of preceder (p) and follower (f) sunspots within ten active regions from 2010-2014 using Space-weather Helioseismic Magnetic Imager Active Region Patch data. The sunspots are small- to mid-sized regions and contain a signed flux within a single polarity sunspot of $(1.1-6.5)\times 10^{21}$ Mx. The net unsigned flux within the regions, including plage, ranges from $(5.1-20)\times 10^{21}$ Mx. Rates are calculated with and without intensity contours to differentiate between sunspot formation and flux emergence. Signed flux emergence rates, calculated with intensity contours, for the p (f) spots average $6.8$ (4.9) $\times 10^{19}$ Mx h$^{-1}$, while decay rates are $-1.9 (-3.4)\times 10^{19}$ Mx h$^{-1}$. The mean, signed flux emergence rate of the regions, including plage, is $7.1 \times 10^{19}$ Mx h$^{-1}$ for a mean peak flux of $5.9 \times 10^{21}$ Mx. Using a synthesis of these results and others reported previously, there is a clear trend for larger flux regions to emerge faster than smaller ones. Observed emergence rates ($dφふぁい/dt$, Mx h$^{-1}$) scale with total signed peak flux, $\tildeφふぁい_{max}$, as a power law with an exponent of 0.36, i.e., $dφふぁい/dt = A \tildeφふぁい_{max}^{0.36}$. The observed rates may assist in constraining the boundary and initial conditions in simulations which already demonstrate increased rates for flux tubes with higher buoyancy and twist, or in the presence of a strong upflow. Overall, the observed emergence rates are smaller than those in simulations, which may indicate a slower rise of the flux in the interior than captured in simulations. △ Less

Submitted 4 May, 2017; originally announced May 2017.

Comments: Accepted by ApJ, 17 pages, 6 figures

Relative magnetic helicity as a diagnostic of solar eruptivity

Authors: E. Pariat, J. E. Leake, G. Valori, M. G. Linton, F. P. Zuccarello, K. Dalmasse

Abstract: The discovery of clear criteria that can deterministically describe the eruptive state of a solar active region would lead to major improvements on space weather predictions. Using series of numerical simulations of the emergence of a magnetic flux rope in a magnetized coronal, leading either to eruptions or to stable configurations, we test several global scalar quantities for the ability to disc… ▽ More The discovery of clear criteria that can deterministically describe the eruptive state of a solar active region would lead to major improvements on space weather predictions. Using series of numerical simulations of the emergence of a magnetic flux rope in a magnetized coronal, leading either to eruptions or to stable configurations, we test several global scalar quantities for the ability to discriminate between the eruptive and the non-eruptive simulations. From the magnetic field generated by the three-dimensional magnetohydrodynamical simulations, we compute and analyse the evolution of the magnetic flux, of the magnetic energy and its decomposition into potential and free energies, and of the relative magnetic helicity and its decomposition. Unlike the magnetic flux and magnetic energies, magnetic helicities are able to markedly distinguish the eruptive from the non-eruptive simulations. We find that the ratio of the magnetic helicity of the current-carrying magnetic field to the total relative helicity presents the highest values for the eruptive simulations, in the pre-eruptive phase only. We observe that the eruptive simulations do not possess the highest value of total magnetic helicity. In the framework of our numerical study, the magnetic energies and the total relative helicity do not correspond to good eruptivity proxies. Our study highlights that the ratio of magnetic helicities diagnoses very clearly the eruptive potential of our parametric simulations. Our study shows that magnetic-helicity-based quantities may be very efficient for the prediction of solar eruptions. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Journal ref: A&A 601, A125 (2017)

Testing the Accuracy of Data Driven MHD Simulations of Active Region Evolution

Authors: James E. Leake, Mark G. Linton, Peter W. Schuck

Abstract: Models for the evolution of the solar coronal magnetic field are vital for understanding solar activity, yet the best measurements of magnetic field lie at the photosphere, necessitating the development of coronal models which are \textit{"data-driven"} at the photosphere. We present an investigation to determine the feasibility and accuracy of such methods. Our validation framework uses a simulat… ▽ More Models for the evolution of the solar coronal magnetic field are vital for understanding solar activity, yet the best measurements of magnetic field lie at the photosphere, necessitating the development of coronal models which are \textit{"data-driven"} at the photosphere. We present an investigation to determine the feasibility and accuracy of such methods. Our validation framework uses a simulation of active region (AR) formation, modeling the emergence of magnetic flux from the convection zone to the corona, as a ground-truth dataset, to supply both the photospheric information, and to perform the validation of the data-driven method. We focus our investigation on how the data-driven model accuracy depends on the temporal frequency of the driving data. The Helioseismic and Magnetic Imager on NASA's Solar {\jl Dynamics} Observatory produces full-disk vector magnetic field measurements at a 12 minute cadence. Using our framework we show that ARs that emerge over 25 hours can be modeled by the data-driving method with only $\sim$1\% error in the free magnetic energy, assuming the photospheric information is specified every 12 minutes. However, for rapidly evolving features, under-sampling of the dynamics at this cadence leads to a strobe effect, generating large electric currents and incorrect coronal morphology and energies. We derive a sampling condition for the driving cadence based on the evolution of these small-scale features, and show that higher-cadence driving can lead to acceptable errors. Future work will investigate the source of errors associated with deriving plasma variables from the photospheric magnetograms as well as other sources of errors such as reduced resolution and instrument bias and noise. △ Less

Submitted 22 February, 2017; originally announced February 2017.

A STEREO Survey of Magnetic Cloud Coronal Mass Ejections Observed at Earth in 2008--2012

Authors: Brian E. Wood, Chin-Chun Wu, Ronald P. Lepping, Teresa Nieves-Chinchilla, Russell A. Howard, Mark G. Linton, Dennis G. Socker

Abstract: We identify coronal mass ejections (CMEs) associated with magnetic clouds (MCs) observed near Earth by the Wind spacecraft from 2008 to mid-2012, a time period when the two STEREO spacecraft were well positioned to study Earth-directed CMEs. We find 31 out of 48 Wind MCs during this period can be clearly connected with a CME that is trackable in STEREO imagery all the way from the Sun to near 1 AUえーゆー… ▽ More We identify coronal mass ejections (CMEs) associated with magnetic clouds (MCs) observed near Earth by the Wind spacecraft from 2008 to mid-2012, a time period when the two STEREO spacecraft were well positioned to study Earth-directed CMEs. We find 31 out of 48 Wind MCs during this period can be clearly connected with a CME that is trackable in STEREO imagery all the way from the Sun to near 1 AUえーゆー. For these events, we perform full 3-D reconstructions of the CME structure and kinematics, assuming a flux rope morphology for the CME shape, considering the full complement of STEREO and SOHO imaging constraints. We find that the flux rope orientations and sizes inferred from imaging are not well correlated with MC orientations and sizes inferred from the Wind data. However, velocities within the MC region are reproduced reasonably well by the image-based reconstruction. Our kinematic measurements are used to provide simple prescriptions for predicting CME arrival times at Earth, provided for a range of distances from the Sun where CME velocity measurements might be made. Finally, we discuss the differences in the morphology and kinematics of CME flux ropes associated with different surface phenomena (flares, filament eruptions, or no surface activity). △ Less

Submitted 6 January, 2017; originally announced January 2017.

Comments: 32 pages, 11 figures, to appear in The Astrophysical Journal Supplement

The Contribution of Coronal Jets To The Solar Wind

Authors: R. Lionello, T. TÖrÖk, V. S. Titov, J. E. Leake, Z. MikiĆ, J. A. Linker, M. G. Linton

Abstract: 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… ▽ More 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. △ Less

Submitted 10 October, 2016; originally announced October 2016.

Imaging Prominence Eruptions Out to 1 AUえーゆー

Authors: Brian E. Wood, Russell A. Howard, Mark G. Linton

Abstract: Views of two bright prominence eruptions trackable all the way to 1AUえーゆー are here presented, using the heliospheric imagers on the Solar TErrestrial RElations Observatory (STEREO) spacecraft. The two events first erupted from the Sun on 2011 June 7 and 2012 August 31, respectively. Only these two examples of clear prominence eruptions observable this far from the Sun could be found in the STEREO imag… ▽ More Views of two bright prominence eruptions trackable all the way to 1AUえーゆー are here presented, using the heliospheric imagers on the Solar TErrestrial RElations Observatory (STEREO) spacecraft. The two events first erupted from the Sun on 2011 June 7 and 2012 August 31, respectively. Only these two examples of clear prominence eruptions observable this far from the Sun could be found in the STEREO image database, emphasizing the rarity of prominence eruptions this persistently bright. For the 2011 June event, a time-dependent 3-D reconstruction of the prominence structure is made using point-by-point triangulation. This is not possible for the August event due to a poor viewing geometry. Unlike the coronal mass ejection (CME) that accompanies it, the 2011 June prominence exhibits little deceleration from the Sun to 1 AUえーゆー, as a consequence moving upwards within the CME. This demonstrates that prominences are not necessarily tied to the CME's magnetic structure far from the Sun. A mathematical framework is developed for describing the degree of self-similarity for the prominence's expansion away from the Sun. This analysis suggests only modest deviations from self-similar expansion, but close to the Sun the prominence expands radially somewhat more rapidly than self-similarity would predict. △ Less

Submitted 21 December, 2015; originally announced December 2015.

Comments: to appear in The Astrophysical Journal

Modeling Jets in the Corona and Solar Wind

Authors: T. Torok, R. Lionello, V. S. Titov, J. E. Leake, Z. Mikic, J. A. Linker, M. G. Linton

Abstract: 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 bri… ▽ More 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. △ Less

Submitted 30 November, 2015; originally announced November 2015.

Comments: 12 pages, 4 figures, to appear in ASP proceedings

The Evolution of Writhe in Kink-Unstable Flux Ropes and Erupting Filaments

Authors: Tibor Torok, Bernhard Kliem, Mitchell A. Berger, Mark G. Linton, Pascal Demoulin, Lidia van Driel-Gesztelyi

Abstract: The helical kink instability of a twisted magnetic flux tube has been suggested as a trigger mechanism for solar filament eruptions and coronal mass ejections (CMEs). In order to investigate if estimations of the pre-eruptive twist can be obtained from observations of writhe in such events, we quantitatively analyze the conversion of twist into writhe in the course of the instability, using numeri… ▽ More The helical kink instability of a twisted magnetic flux tube has been suggested as a trigger mechanism for solar filament eruptions and coronal mass ejections (CMEs). In order to investigate if estimations of the pre-eruptive twist can be obtained from observations of writhe in such events, we quantitatively analyze the conversion of twist into writhe in the course of the instability, using numerical simulations. We consider the line tied, cylindrically symmetric Gold-Hoyle flux rope model and measure the writhe using the formulae by Berger and Prior which express the quantity as a single integral in space. We find that the amount of twist converted into writhe does not simply scale with the initial flux rope twist, but depends mainly on the growth rates of the instability eigenmodes of higher longitudinal order than the basic mode. The saturation levels of the writhe, as well as the shapes of the kinked flux ropes, are very similar for considerable ranges of initial flux rope twists, which essentially precludes estimations of pre-eruptive twist from measurements of writhe. However, our simulations suggest an upper twist limit of $\sim 6πぱい$ for the majority of filaments prior to their eruption. △ Less

Submitted 6 March, 2014; originally announced March 2014.

Comments: 16 pages, 4 figures, to appear in Plasma Physics and Controlled Fusion

On Flux Rope Stability and Atmospheric Stratification in Models of Coronal Mass Ejections Triggered by Flux Emergence

Authors: E. Lee, V. S. Lukin, M. G. Linton

Abstract: Flux emergence is widely recognized to play an important role in the initiation of coronal mass ejections. The Chen-Shibata (2000) model, which addresses the connection between emerging flux and flux rope eruptions, can be implemented numerically to study how emerging flux through the photosphere can impact the eruption of a pre-existing coronal flux rope. The model's sensitivity to the initial co… ▽ More Flux emergence is widely recognized to play an important role in the initiation of coronal mass ejections. The Chen-Shibata (2000) model, which addresses the connection between emerging flux and flux rope eruptions, can be implemented numerically to study how emerging flux through the photosphere can impact the eruption of a pre-existing coronal flux rope. The model's sensitivity to the initial conditions and reconnection micro-physics is investigated with a parameter study. In particular, we aim to understand the stability of the coronal flux rope in the context of X-point collapse and the effects of boundary driving in both unstratified and stratified atmospheres. In the absence of driving, we assess the behavior of waves in the vicinity of the X-point. With boundary driving applied, we study the effects of reconnection micro-physics and atmospheric stratification on the eruption. We find that the Chen-Shibata equilibrium can be unstable to an X-point collapse even in the absence of driving due to wave accumulation at the X-point. However, the equilibrium can be stabilized by reducing the compressibility of the plasma, which allows small-amplitude waves to pass through the X-point without accumulation. Simulations with the photospheric boundary driving evaluate the impact of reconnection micro-physics and atmospheric stratification on the resulting dynamics: we show the evolution of the system to be determined primarily by the structure of the global magnetic fields with little sensitivity to the micro-physics of magnetic reconnection; and in a stratified atmosphere, we identify a novel mechanism for producing quasi-periodic behavior at the reconnection site behind a rising flux rope as a possible explanation of similar phenomena observed in solar and stellar flares. △ Less

Submitted 14 August, 2014; v1 submitted 2 March, 2014; originally announced March 2014.

Comments: Submitted Feb 28, 2014 to, accepted Aug 14, 2014 by Astronomy & Astrophysics. 13 pages, 10 figures, 2 tables

Journal ref: A&A 569, A94 (2014)

Simulations of Emerging Magnetic Flux. II: The formation of Unstable Coronal Flux Ropes and the Initiation of CMEs

Authors: J. E. Leake, M. G. Linton, S. K. Antiochos

Abstract: We present results from 3D magnetohydrodynamic (MHD) simulations of the emergence of a twisted convection zone flux tube into a pre-existing coronal dipole field. As in previous simulations, following the partial emergence of the sub-surface flux into the corona, a combination of vortical motions and internal magnetic reconnection forms a coronal flux rope. Then, in the simulations presented here,… ▽ More We present results from 3D magnetohydrodynamic (MHD) simulations of the emergence of a twisted convection zone flux tube into a pre-existing coronal dipole field. As in previous simulations, following the partial emergence of the sub-surface flux into the corona, a combination of vortical motions and internal magnetic reconnection forms a coronal flux rope. Then, in the simulations presented here, external reconnection between the emerging field and the pre-existing dipole coronal field allows further expansion of the coronal flux rope into the corona. After sufficient expansion, internal reconnection occurs beneath the coronal flux rope axis, and the flux rope erupts up to the top boundary of the simulation domain ($\sim$ 36 Mm above the surface). We find that the presence of a pre-existing field, orientated in a direction to facilitate reconnection with the emerging field, is vital to the fast rise of the coronal flux rope. The simulations shown in this paper are able to self-consistently create many of the surface and coronal signatures used by coronal mass ejection (CME) models. These signatures include: surface shearing and rotational motions; quadrupolar geometry above the surface; central sheared arcades reconnecting with oppositely orientated overlying dipole fields; the formation of coronal flux ropes underlying potential coronal field; and internal reconnection which resembles the classical flare reconnection scenario. This suggests that proposed mechanisms for the initiation of a CME, such as "magnetic breakout", are operating during the emergence of new active regions. △ Less

Submitted 26 March, 2014; v1 submitted 7 February, 2014; originally announced February 2014.

Distribution of Electric Currents in Solar Active Regions

Authors: Tibor Török, James E. Leake, Viacheslav S. Titov, Vasilis Archontis, Zoran Mikić, Mark G. Linton, Kévin Dalmasse, Guillaume Aulanier, Bernhard Kliem

Abstract: There has been a long-lasting debate on the question of whether or not electric currents in solar active regions are neutralized. That is, whether or not the main (or direct) coronal currents connecting the active region polarities are surrounded by shielding (or return) currents of equal total value and opposite direction. Both theory and observations are not yet fully conclusive regarding this q… ▽ More There has been a long-lasting debate on the question of whether or not electric currents in solar active regions are neutralized. That is, whether or not the main (or direct) coronal currents connecting the active region polarities are surrounded by shielding (or return) currents of equal total value and opposite direction. Both theory and observations are not yet fully conclusive regarding this question, and numerical simulations have, surprisingly, barely been used to address it. Here we quantify the evolution of electric currents during the formation of a bipolar active region by considering a three-dimensional magnetohydrodynamic simulation of the emergence of a sub-photospheric, current-neutralized magnetic flux rope into the solar atmosphere. We find that a strong deviation from current neutralization develops simultaneously with the onset of significant flux emergence into the corona, accompanied by the development of substantial magnetic shear along the active region's polarity inversion line. After the region has formed and flux emergence has ceased, the strong magnetic fields in the region's center are connected solely by direct currents, and the total direct current is several times larger than the total return current. These results suggest that active regions, the main sources of coronal mass ejections and flares, are born with substantial net currents, in agreement with recent observations. Furthermore, they support eruption models that employ pre-eruption magnetic fields containing such currents. △ Less

Submitted 13 January, 2014; originally announced January 2014.

Comments: 6 pages, 5 figures, to appear in Astrophysical Journal Letters

A Plasma βべーた Transition Within a Propagating Flux Rope

Authors: Neel P. Savani, A. Vourlidas, D. Shiota, M. G. Linton, K. Kusano, N. Lugaz, A. P. Rouillard

Abstract: We present a 2.5D MHD simulation of a magnetic flux rope (FR) propagating in the heliosphere and investigate the cause of the observed sharp plasma beta transition. Specifically, we consider a strong internal magnetic field and an explosive fast start, such that the plasma beta is significantly lower in the FR than the sheath region that is formed ahead. This leads to an unusual FR morphology in t… ▽ More We present a 2.5D MHD simulation of a magnetic flux rope (FR) propagating in the heliosphere and investigate the cause of the observed sharp plasma beta transition. Specifically, we consider a strong internal magnetic field and an explosive fast start, such that the plasma beta is significantly lower in the FR than the sheath region that is formed ahead. This leads to an unusual FR morphology in the first stage of propagation, while the more traditional view (e.g. from space weather simulations like Enlil) of a `pancake' shaped FR is observed as it approaches 1 AUえーゆー. We investigate how an equipartition line, defined by a magnetic Weber number, surrounding a core region of a propagating FR can demarcate a boundary layer where there is a sharp transition in the plasma beta. The substructure affects the distribution of toroidal flux, with the majority of the flux remaining in a small core region which maintains a quasi-cylindrical structure. Quantitatively, we investigate a locus of points where the kinetic energy density of the relative inflow field is equal to the energy density of the transverse magnetic field (i.e. effective tension force). The simulation provides compelling evidence that at all heliocentric distances the distribution of toroidal magnetic flux away from the FR axis is not linear; with 80% of the toroidal flux occurring within 40% of the distance from the FR axis. Thus our simulation displays evidence that the competing ideas of a pancaking structure observed remotely can coexist with a quasi-cylindrical magnetic structure seen in situ. △ Less

Submitted 27 November, 2013; v1 submitted 17 October, 2013; originally announced October 2013.

Comments: 11 pages of text + 6 figures. Accepted to ApJ on 16 Oct 2013

Ionized Plasma and Neutral Gas Coupling in the Sun's Chromosphere and Earth's Ionosphere/Thermosphere

Authors: J. E. Leake, C. R. DeVore, J. P. Thayer, A. G. Burns, G. Crowley, H. R. Gilbert, J. D Huba, J. Krall, M. G. Linton, V. S. Lukin, W. Wang

Abstract: We review our understanding of ionized plasma and neutral gas coupling in the weakly ionized, stratified, electromagnetically-permeated regions of the Sun's chromosphere and Earth's ionosphere/thermosphere. Using representative models for each environment we derive fundamental descriptions of the coupling of the constituent parts to each other and to the electric and magnetic fields, and we examin… ▽ More We review our understanding of ionized plasma and neutral gas coupling in the weakly ionized, stratified, electromagnetically-permeated regions of the Sun's chromosphere and Earth's ionosphere/thermosphere. Using representative models for each environment we derive fundamental descriptions of the coupling of the constituent parts to each other and to the electric and magnetic fields, and we examine the variation in magnetization of the ionized component. Using these descriptions we compare related phenomena in the two environments, and discuss electric currents, energy transfer and dissipation. We present a coupled theoretical and numerical study of plasma instabilities in the two environments that serves as an example of how the chromospheric and ionospheric communities can further collaborate. We also suggest future collaborative studies that will help improve our understanding of these two different atmospheres which share many similarities, but have large disparities in key quantities. △ Less

Submitted 29 August, 2016; v1 submitted 1 October, 2013; originally announced October 2013.

Simulations of Emerging Magnetic Flux. I: The Formation of Stable Coronal Flux Ropes

Authors: James E. Leake, Mark G. Linton, Tibor Török

Abstract: We present results from 3D visco-resistive magnetohydrodynamic (MHD) simulations of the emergence of a convection zone magnetic flux tube into a solar atmosphere containing a pre-existing dipole coronal field, which is orientated to minimize reconnection with the emerging field. We observe that the emergence process is capable of producing a coronal flux rope by the transfer of twist from the conv… ▽ More We present results from 3D visco-resistive magnetohydrodynamic (MHD) simulations of the emergence of a convection zone magnetic flux tube into a solar atmosphere containing a pre-existing dipole coronal field, which is orientated to minimize reconnection with the emerging field. We observe that the emergence process is capable of producing a coronal flux rope by the transfer of twist from the convection zone as found in previous simulations. We find that this flux rope is stable, with no evidence of a fast rise, and that its ultimate height in the corona is determined by the strength of the pre-existing dipole field. We also find that although the electric currents in the initial convection zone flux tube are almost perfectly neutralized, the resultant coronal flux rope carries a significant net current. These results suggest that flux tube emergence is capable of creating non-current-neutralized stable flux ropes in the corona, tethered by overlying potential fields, a magnetic configuration that is believed to be the source of coronal mass ejections. △ Less

Submitted 28 August, 2013; originally announced August 2013.

Comments: 14 figures

Magnetic Reconnection in a Weakly Ionized Plasma

Authors: James E Leake, Vyacheslav S Lukin, Mark G Linton

Abstract: Magnetic reconnection in partially ionized plasmas is a ubiquitous phenomenon spanning the range from laboratory to intergalactic scales, yet it remains poorly understood and relatively little studied. Here, we present results from a self-consistent multi-fluid simulation of magnetic reconnection in a weakly ionized reacting plasma with a particular focus on the parameter regime of the solar chrom… ▽ More Magnetic reconnection in partially ionized plasmas is a ubiquitous phenomenon spanning the range from laboratory to intergalactic scales, yet it remains poorly understood and relatively little studied. Here, we present results from a self-consistent multi-fluid simulation of magnetic reconnection in a weakly ionized reacting plasma with a particular focus on the parameter regime of the solar chromosphere. The numerical model includes collisional transport, interaction and reactions between the species, and optically thin radiative losses. This model improves upon our previous work in Leake et al. 2012 \cite{Leake2012} by considering realistic chromospheric transport coefficients, and by solving a generalized Ohm's law that accounts for finite ion-inertia and electron-neutral drag. We find that during the two dimensional reconnection of a Harris current sheet with an initial width larger than the neutral-ion collisional coupling scale, the current sheet thins until its width becomes less than this coupling scale, and the neutral and ion fluids decouple upstream from the reconnection site. During this process of decoupling, we observe reconnection faster than the single-fluid Sweet-Parker prediction, with recombination and plasma outflow both playing a role in determining the reconnection rate. As the current sheet thins further and elongates it becomes unstable to the secondary tearing instability, and plasmoids are seen. The reconnection rate, outflows and plasmoids observed in this simulation provide evidence that magnetic reconnection in the chromosphere could be responsible for jet-like transient phenomena such as spicules and chromospheric jets. △ Less

Submitted 13 February, 2013; originally announced February 2013.

Comments: 11 Figures

Journal ref: Phys. Plasmas 20, 061202 (2013)

Effect of ion-neutral collisions in simulations of emerging active regions

Authors: James E. Leake, Mark G. Linton

Abstract: We present results of 2.5D numerical simulations of the emergence of sub-surface magnetic flux into the solar atmosphere, with emerging flux regions ranging from $10^{18}$ to $10^{21}$ Mx, representing both ephemeral and active regions. We include the presence of neutral Hydrogen in the governing equations, improve upon previous models by including the ionization in the equation of state, and use… ▽ More We present results of 2.5D numerical simulations of the emergence of sub-surface magnetic flux into the solar atmosphere, with emerging flux regions ranging from $10^{18}$ to $10^{21}$ Mx, representing both ephemeral and active regions. We include the presence of neutral Hydrogen in the governing equations, improve upon previous models by including the ionization in the equation of state, and use a more realistic convection zone model. We find that ionization and recombination of plasma during the rise of a convection zone flux tube reduces the rise speed of the tube's axis. The presence of neutral Hydrogen allows the effective flow of mass across fieldlines, by the addition of a Pedersen resistivity to the generalized Ohm's law, which dissipates current perpendicular to the magnetic field. This causes an increase of up to 10% in the amount of magnetic in-plane flux supplied to the corona and a reduction of up to 89% in the amount of sub-surface plasma brought up into the corona. However, it also reduces the amount of free magnetic energy supplied to the corona, and thus does not positively affect the likelihood of creating unstable coronal structures. △ Less

Submitted 17 December, 2012; originally announced December 2012.

Comments: 18 figures

Multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma

Authors: James E. Leake, Vyacheslav S. Lukin, Mark G. Linton, Eric T. Meier

Abstract: We present results from the first self-consistent multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma. We simulate two dimensional magnetic reconnection in a Harris current sheet with a numerical model which includes ion-neutral scattering collisions, ionization, recombination, optically thin radiative loss, collisional heating, and thermal conduction… ▽ More We present results from the first self-consistent multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma. We simulate two dimensional magnetic reconnection in a Harris current sheet with a numerical model which includes ion-neutral scattering collisions, ionization, recombination, optically thin radiative loss, collisional heating, and thermal conduction. In the resulting tearing mode reconnection the neutral and ion fluids become decoupled upstream from the reconnection site, creating an excess of ions in the reconnection region and therefore an ionization imbalance. Ion recombination in the reconnection region, combined with Alfvénic outflows, quickly removes ions from the reconnection site, leading to a fast reconnection rate independent of Lundquist number. In addition to allowing fast reconnection, we find that these non-equilibria partial ionization effects lead to the onset of the nonlinear secondary tearing instability at lower values of the Lundquist number than has been found in fully ionized plasmas.These simulations provide evidence that magnetic reconnection in the chromosphere could be responsible for jet-like transient phenomena such as spicules and chromospheric jets. △ Less

Submitted 5 October, 2012; originally announced October 2012.

Comments: 8 Figures, 32 pages total

Journal ref: ApJ 760, 109 (2012)

Tests of Dynamical Flux Emergence as a Mechanism for CME Initiation

Authors: James E. Leake, Mark G. Linton, Spiro K. Antiochos

Abstract: Current coronal mass ejection (CME) models set their lower boundary to be in the lower corona. They do not calculate accurately the transfer of free magnetic energy from the convection zone to the magnetically dominated corona because they model the effects of flux emergence using kinematic boundary conditions or simply assume the appearance of flux at these heights. We test the importance of incl… ▽ More Current coronal mass ejection (CME) models set their lower boundary to be in the lower corona. They do not calculate accurately the transfer of free magnetic energy from the convection zone to the magnetically dominated corona because they model the effects of flux emergence using kinematic boundary conditions or simply assume the appearance of flux at these heights. We test the importance of including dynamical flux emergence in CME modeling by simulating, in 2.5D, the emergence of sub-surface flux tubes into different coronal magnetic field configurations. We investigate how much free magnetic energy, in the form of shear magnetic field, is transported from the convection zone to the corona, and whether dynamical flux emergence can drive CMEs. We find that multiple coronal flux ropes can be formed during flux emergence, and although they carry some shear field into the corona, the majority of shear field is confined to the lower atmosphere. Less than 10% of the magnetic energy in the corona is in the shear field, and this, combined with the fact that the coronal flux ropes bring up significant dense material, means that they do not erupt. Our results have significant implications for all CME models which rely on the transfer of free magnetic energy from the lower atmosphere into the corona but which do not explicitly model this transfer. Such studies of flux emergence and CMEs are timely, as we have new capabilities to observe this with Hinode and SDO, and therefore to test the models against observations. △ Less

Submitted 3 August, 2010; v1 submitted 30 July, 2010; originally announced July 2010.

Reconnection of a kinking flux rope triggering the ejection of a microwave and hard X-ray source. II. Numerical Modeling

Authors: B. Kliem, M. G. Linton, T. Török, M. Karlický

Abstract: Numerical simulations of the helical ($m\!=\!1$) kink instability of an arched, line-tied flux rope demonstrate that the helical deformation enforces reconnection between the legs of the rope if modes with two helical turns are dominant as a result of high initial twist in the range $Φふぁい\gtrsim6πぱい$. Such reconnection is complex, involving also the ambient field. In addition to breaking up the origina… ▽ More Numerical simulations of the helical ($m\!=\!1$) kink instability of an arched, line-tied flux rope demonstrate that the helical deformation enforces reconnection between the legs of the rope if modes with two helical turns are dominant as a result of high initial twist in the range $Φふぁい\gtrsim6πぱい$. Such reconnection is complex, involving also the ambient field. In addition to breaking up the original rope, it can form a new, low-lying, less twisted flux rope. The new flux rope is pushed downward by the reconnection outflow, which typically forces it to break as well by reconnecting with the ambient field. The top part of the original rope, largely rooted in the sources of the ambient flux after the break-up, can fully erupt or be halted at low heights, producing a "failed eruption." The helical current sheet associated with the instability is squeezed between the approaching legs, temporarily forming a double current sheet. The leg-leg reconnection proceeds at a high rate, producing sufficiently strong electric fields that it would be able to accelerate particles. It may also form plasmoids, or plasmoid-like structures, which trap energetic particles and propagate out of the reconnection region up to the top of the erupting flux rope along the helical current sheet. The kinking of a highly twisted flux rope involving leg-leg reconnection can explain key features of an eruptive but partially occulted solar flare on 18 April 2001, which ejected a relatively compact hard X-ray and microwave source and was associated with a fast coronal mass ejection. △ Less

Submitted 13 July, 2010; originally announced July 2010.

Comments: Solar Physics, in press

Gas-dynamic shock heating of post-flare loops due to retraction following localized, impulsive reconnection

Authors: D. W. Longcope, S. E. Guidoni, M. G. Linton

Abstract: We present a novel model in which shortening of a magnetic flux tube following localized, three-dimensional reconnection generates strong gas-dynamic shocks around its apex. The shortening releases magnetic energy by progressing away from the reconnection site at the Alfven speed. This launches inward flows along the field lines whose collision creates a pair of gas-dynamic shocks. The shocks ra… ▽ More We present a novel model in which shortening of a magnetic flux tube following localized, three-dimensional reconnection generates strong gas-dynamic shocks around its apex. The shortening releases magnetic energy by progressing away from the reconnection site at the Alfven speed. This launches inward flows along the field lines whose collision creates a pair of gas-dynamic shocks. The shocks raise both the mass density and temperature inside the newly shortened flux tube. Reconnecting field lines whose initial directions differ by more that 100 degrees can produce a concentrated knot of plasma hotter that 20 MK, consistent with observations. In spite of these high temperatures, the shocks convert less than 10% of the liberated magnetic energy into heat - the rest remains as kinetic energy of bulk motion. These gas-dynamic shocks arise only when the reconnection is impulsive and localized in all three dimensions; they are distinct from the slow magnetosonic shocks of the Petschek steady-state reconnection model. △ Less

Submitted 3 November, 2008; v1 submitted 20 October, 2008; originally announced October 2008.

Patchy Reconnection in a Y-Type Current Sheet

Authors: M. G. Linton, C. R. DeVore, D. W. Longcope

Abstract: We study the evolution of the magnetic field in a Y-type current sheet subject to a brief, localized magnetic reconnection event. The reconnection produces up- and down-flowing reconnected flux tubes which rapidly decelerate when they hit the Y-lines and underlying magnetic arcade loops at the ends of the current sheet. This localized reconnection outflow followed by a rapid deceleration reprodu… ▽ More We study the evolution of the magnetic field in a Y-type current sheet subject to a brief, localized magnetic reconnection event. The reconnection produces up- and down-flowing reconnected flux tubes which rapidly decelerate when they hit the Y-lines and underlying magnetic arcade loops at the ends of the current sheet. This localized reconnection outflow followed by a rapid deceleration reproduces the observed behavior of post-CME downflowing coronal voids. These simulations support the hypothesis that these observed coronal downflows are the retraction of magnetic fields reconnected in localized patches in the high corona. △ Less

Submitted 1 March, 2010; v1 submitted 7 December, 2007; originally announced December 2007.

Comments: 4 pages, 3 figures

Journal ref: Earth Planets, and Space (2009), 61 p573

A Model for Patchy Reconnection in Three Dimensions

Authors: M. G. Linton, D. W. Longcope

Abstract: We show, theoretically and via MHD simulations, how a short burst of reconnection localized in three dimensions on a one-dimensional current sheet creates a pair of reconnected flux tubes. We focus on the post-reconnection evolution of these flux tubes, studying their velocities and shapes. We find that slow-mode shocks propagate along these reconnected flux tubes, releasing magnetic energy as i… ▽ More We show, theoretically and via MHD simulations, how a short burst of reconnection localized in three dimensions on a one-dimensional current sheet creates a pair of reconnected flux tubes. We focus on the post-reconnection evolution of these flux tubes, studying their velocities and shapes. We find that slow-mode shocks propagate along these reconnected flux tubes, releasing magnetic energy as in steady-state Petschek reconnection. The geometry of these three-dimensional shocks, however, differs dramatically from the classical two-dimensional geometry. They propagate along the flux tube legs in four isolated fronts, whereas in the two-dimensional Petschek model, they form a continuous, stationary pair of V-shaped fronts. We find that the cross sections of these reconnected flux tubes appear as teardrop shaped bundles of flux propagating away from the reconnection site. Based on this, we argue that the descending coronal voids seen by Yohkoh SXT, LASCO, and TRACE are reconnected flux tubes descending from a flare site in the high corona, for example after a coronal mass ejection. In this model, these flux tubes would then settle into equilibrium in the low corona, forming an arcade of post-flare coronal loops. △ Less

Submitted 13 September, 2005; originally announced September 2005.

Comments: 27 pages plus 16 figures

Journal ref: Astrophys.J.642:1177-1192,2006

Magnetic Flux Tube Reconnection: Tunneling Versus Slingshot

Authors: M. G. Linton, S. K. Antiochos

Abstract: The discrete nature of the solar magnetic field as it emerges into the corona through the photosphere indicates that it exists as isolated flux tubes in the convection zone, and will remain as discrete flux tubes in the corona until it collides and reconnects with other coronal fields. Collisions of these flux tubes will in general be three dimensional, and will often lead to reconnection, both… ▽ More The discrete nature of the solar magnetic field as it emerges into the corona through the photosphere indicates that it exists as isolated flux tubes in the convection zone, and will remain as discrete flux tubes in the corona until it collides and reconnects with other coronal fields. Collisions of these flux tubes will in general be three dimensional, and will often lead to reconnection, both rearranging the magnetic field topology in fundamental ways, and releasing magnetic energy. With the goal of better understanding these dynamics, we carry out a set of numerical experiments exploring fundamental characteristics of three dimensional magnetic flux tube reconnection. We first show that reconnecting flux tubes at opposite extremes of twist behave very differently: in some configurations, low twist tubes slingshot while high twist tubes tunnel. We then discuss a theory explaining these differences: by assuming helicity conservation during the reconnection one can show that at high twist, tunneled tubes reach a lower magnetic energy state than slingshot tubes, whereas at low twist the opposite holds. We test three predictions made by this theory. 1) We find that the level of twist at which the transition from slingshot to tunnel occurs is about two to three times higher than predicted on the basis of energetics and helicity conservation alone, probably because the dynamics of the reconnection play a large role as well. 2) We find that the tunnel occurs at all flux tube collision angles predicted by the theory. 3) We find that the amount of magnetic energy a slingshot or a tunnel reconnection releases agrees reasonably well with the theory, though at the high resistivities we have to use for numerical stability, a significant amount of magnetic energy is lost to diffusion, independent of reconnection. △ Less

Submitted 21 January, 2005; originally announced January 2005.

Comments: 21 pages, 15 figures, submitted to ApJ

Journal ref: Astrophys.J. 625 (2005) 506-521