Search

Transient collimated plasma ejections (jets) occur frequently throughout the solar corona, in active regions, quiet Sun, and coronal holes. Although magnetic reconnection is generally agreed to be the mechanism of energy release in jets, the factors that dictate the location and rate of reconnection remain unclear. Our previous studies demonstrated that the magnetic breakout model explains the triggering and evolution of most jets over a wide range of scales, through detailed comparisons between our numerical simulations and high-resolution observations. An alternative explanation, the resistive-kink model, invokes breakout reconnection without forming and explosively expelling a flux rope. Here we report direct observations of breakout reconnection and plasmoid formation during two jets in the fan-spine topology of an embedded bipole. For the first time, we observed the formation and evolution of multiple small plasmoids with bidirectional flows associated with fast reconnection in 3D breakout current sheets (BCSs) in the solar corona. The first narrow jet was launched by reconnection at the BCS originating at the deformed 3D null, without significant flare reconnection or a filament eruption. In contrast, the second jet and release of cool filament plasma were triggered by explosive breakout reconnection when the leading edge of the rising flux rope formed by flare reconnection beneath the filament encountered the preexisting BCS. These observations solidly support both reconnection-driven jet models: the resistive kink for the first jet, and the breakout model for the second explosive jet with a filament eruption.

Jets (transient/collimated plasma ejections) occur frequently throughout the solar corona and contribute mass/energy to the corona and solar wind. By combining numerical simulations and high-resolution observations, we have made substantial progress recently on determining the energy buildup and release processes in these jets. Here we describe a study of 27 equatorial coronal-hole jets using Solar Dynamics Observatory/Atmospheric Imaging Assembly and Helioseismic and Magnetic Imager observations on 2013 June 27–28 and 2014 January 8–10. Out of 27 jets, 18 (67%) are associated with mini-filament ejections; the other nine (33%) do not show mini-filament eruptions but do exhibit mini-flare arcades and other eruptive signatures. This indicates that every jet in our sample involved a filament-channel eruption. From the complete set of events, six jets (22%) are apparently associated with tiny flux-cancellation events at the polarity inversion line, and two jets (7%) are associated with sympathetic eruptions of filaments from neighboring bright points. Potential-field extrapolations of the source-region photospheric magnetic fields reveal that all jets originated in the fan-spine topology of an embedded bipole associated with an extreme ultraviolet coronal bright point. Hence, all our jets are in agreement with the breakout model of solar eruptions. We present selected examples and discuss the implications for the jet energy buildup and initiation mechanisms.

Solar coronal jets are small, transient, collimated ejections most easily observed in coronal holes (CHs). The upcoming Parker Solar Probe (PSP) mission provides the first opportunity to encounter CH jets in situ near the Sun and examine their internal structure and dynamics. Using projected mission orbital parameters, we have simulated PSP encounters with a fully three-dimensional magnetohydrodynamic (MHD) model of a CH jet. We find that three internal jet regions, featuring different wave modes and levels of compressibility, have distinct identifying signatures detectable by PSP. The leading Alfvén wave front and its immediate wake are characterized by transAlfvénic plasma flows with mild density enhancements. This front exhibits characteristics of a fast switch-on MHD shock, whose arrival is signaled by the sudden onset of large-amplitude transverse velocity and magnetic-field oscillations highly correlated in space and time. The trailing portion is characterized by supersonic but subAlfvénicout flows of dense plasma with uncorrelated velocity and magnetic-field oscillations. This compressible region contains most of the jet's mass. The volume between the immediate wake and dense jet, the remote wake,mixes and transitions the characteristics of the two other regions. In addition to probing each region separately, we also simulate a corotational PSP-jet encounter. In this scenario, the simulated spacecraft hovers over the jet producing CH, as may occur during the mission's corotational phases, sampling each jet region in turn. We estimate that PSP will encounter numerous CH jets over the lifetime of the mission.

Small, impulsive jets commonly occur throughout the solar corona, but are especially visible in coronal holes. Evidence is mounting that jets are part of a continuum of eruptions that extends to much larger coronal mass ejections and eruptive flares. Because coronal-hole jets originate in relatively simple magnetic structures, they offer an ideal testbed for theories of energy buildup and release in the full range of solar eruptions. We analyzed an equatorial coronal-hole jet observed by the Solar Dynamics Observatory (SDO)/AIA (Atmospheric Imaging Assembly)) on 2014 January 9 in which the magnetic-field structure was consistent with the embedded-bipole topology that we identified and modeled previously as an origin of coronal jets. In addition, this event contained a mini-filament, which led to important insights into the energy storage and release mechanisms. SDO/HMI (Solar Dynamics Observatory/Helioseismic and Magnetic Imager) magnetograms revealed footpoint motions in the primary minority-polarity region at the eruption site, but show negligible flux emergence or cancellation for at least 16 hours before the eruption. Therefore, the free energy powering this jet probably came from magnetic shear concentrated at the polarity inversion line within the embedded bipole. We find that the observed activity sequence and its interpretation closely match the predictions of the breakout jet model, strongly supporting the hypothesis that the breakout model can explain solar eruptions on a wide range of scales.

One of the leading models for electron acceleration in solar flares is stochastic acceleration by weakly turbulent fast magnetosonic waves ("fast waves"). In this model, large-scale flows triggered by magnetic reconnection excite large-wavelength fast waves, and fast-wave energy then cascades from large wavelengths to small wavelengths. Electron acceleration by large-wavelength fast-waves is weak, and so the model relies on the small-wavelength waves produced by the turbulent cascade. In order for the model to work, the energy cascade time for large-wavelength fast waves must be shorter than the time required for the waves to propagate out of the solar-flare acceleration region. To investigate the effects of wave escape, we solve the wave kinetic equation for fast waves in weak turbulence theory, supplemented with a homogeneous wave-loss term.We find that the amplitude of large-wavelength fast waves must exceed a minimum threshold in order for a significant fraction of the wave energy to cascade to small wavelengths before the waves leave the acceleration region.We evaluate this threshold as a function of the dominant wavelength of the fast waves that are initially excited by reconnection outflows.

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 60 degrees, suggesting that its source extends far from the open-closed boundary. We describe 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 (the S-Web) and quasi-separatrix layers in the heliosphere. We discuss the dynamics of the S-Web model and its implications for present observations and for the upcoming observations from Solar Orbiter and Solar Probe Plus.

We have developed a comprehensive, multi-threaded, three-dimensional model of the plasma dynamics and energetics of a prominence and its overlying arcade (Luna et al. 2012). In this model, the basic magnetic structure is that of two interacting sheared arcades, while the cool condensations composing the prominence are formed by the well-studied thermal nonequilibrium mechanism. In a given filament-channel flux tube, the mass is evaporated from the chromosphere by heating localized near the footpoints, and condenses in the form of transient blobs or a persistent thread. Our previous studies of thermal nonequilibrium used steady or impulsive heating functions with no dependence on local physical. However, parametric active-region models with steady heating proportional to B/L, where B is the flux-tube magnetic field strength at each footpoint and L is the flux-tube length, yield the best agreement with observations (e.g., Schrijver et al. 2008). We have determined the effects of this active-region heating function on our model for the formation and evolution of prominence mass. We have also expanded the range of our computational domain to include more ofthe overlying arcade (the so-called "cavity") than in Luna et al. (2012), and have increased the number of selected flux tubes from 125 to 533. We will illustrate the time-dependent plasma behavior produced by the B/L heating function with synthetic images in several ALA passbands, and compare the resulting prominence properties with those predicted by our model with steady heating.

We have investigated the onset and acceleration of coronal mass ejections (CMEs) and eruptive flares. To isolate the eruption physics, our study uses the breakout model, which is insensitive to the energy buildup process leading to the eruption. We performed 2.5D simulations with adaptive mesh refinement that achieved the highest overall spatial resolution to date in a CME/eruptive flare simulation. The ultra-high resolution allows us to separate clearly the timing of the various phases of the eruption. Using new computational tools, we have determined the number and evolution of all X- and O-type nulls in the system, thereby tracking both the progress and the products of reconnection throughout the computational domain. Our results show definitively that CME onset is due to the start of fast reconnection at the breakout current sheet. Once this reconnection begins, eruption is inevitable; if this is the only reconnection in the system, however, the eruption will be slow. The explosive CME acceleration is triggered by fast reconnection at the flare current sheet. Our results indicate that the explosive eruption is caused by a resistive instability, not an ideal process. Moreover, both breakout and flare reconnections begin first as a form of weak tearing characterized by a slowly evolving plasmoids, but eventually transition to a fast form with well-defined Alfvenic reconnection jets and rapid flux transfer. This transition to fast reconnection is required for both CME onset and explosive acceleration. We discuss the key implications of our results for CME/flare observations and for theories of magnetic reconnection.

Magnetic reconnection is widely believed to be the primary process by which the magnetic field releases energy to plasma in the Sun's corona. For example, in the breakout model for the initiation of coronal mass ejections/eruptive flares, reconnection is responsible for the catastrophic destabilizing of magnetic force balance in the corona, leading to explosive energy release. A critical requirement for the reconnection is that it have a "switch-on' nature in that the reconnection stays off until a large store of magnetic free energy has built up, and then it turn on abruptly and stay on until most of this free energy has been released. We discuss the implications of this requirement for reconnection in the context of the breakout model for CMEs/flares. We argue that it imposes stringent constraints on the properties of the flux breaking mechanism, which is expected to operate in the corona on kinetic scales. We present numerical simulations demonstrating how the reconnection and the eruption depend on the effective resistivity, i.e., the effective Lundquist number, and propose a model for incorporating kinetic flux-breaking mechanisms into MHO calculation of CMEs/flares.

Cancellation of magnetic flux in the solar photosphere and chromosphere has been linked observationally and theoretically to a broad range of solar activity, from filament channel formation to CME initiation. Because this phenomenon is typically measured at only a single layer in the atmosphere, in the radial (line of sight) component of the magnetic field, the actual processes behind this observational signature are ambiguous. It is clear that reconnection is involved in some way, but the location of the reconnection sites and associated connectivity changes remain uncertain in most cases. We are using numerical modeling to demystify flux cancellation, beginning with the simplest possible configuration: a subphotospheric Lundquist flux tube surrounded by a potential field, immersed in a gravitationally stratified atmosphere, spanning many orders of magnitude in plasma beta. In this system, cancellation is driven slowly by a 2-cell circulation pattern imposed in the convection zone, such that the tops of the cells are located around the beta= 1 level (Le., the photosphere) and the flows converge and form a downdraft at the polarity inversion line; note however that no flow is imposed along the neutral line. We will present the results of 2D and 3D MHD-AMR simulations of flux cancellation, in which the flux at the photosphere begins in either an unsheared or sheared state. In all cases, a lOW-lying flux rope is formed by reconnection at the polarity inversion line within a few thousand seconds. The flux rope remains stable and does not rise, however, in contrast to models which do not include the presence of significant mass loading.

Transition to turbulence in magnetohydrodynamic (MHD) tearing jets has been invoked as a mechanism underlying some of the complex behavior observed in solar surges, including deceleration of the upflowing plasma and temporal correlations with types I and III radio bursts. In this paper we investigate a possible mechanism for this transition: three-dimensional secondary instabilities on two-dimensional saturated states. We find through linear analysis that these MHD configurations -- in particular, the tearing jet -- are secondarily unstable, with the dominant energy transfer from the one-dimensional field into the 3-dimensional fields. Using nonlinear simulations, we also investigate the system evolution after the secondary modes attain finite amplitude. When the tearing jet transitions to turbulence, the total kinetic energy drops rapidly corresponding to the deceleration of the jet. The electric field grows rapidly as the primary mode saturates and the three-dimensional secondary mode develops, and then decays quickly as the tearing jet becomes turbulent, providing a possible explanation for the finite duration of the associated meter-wave bursts. The electric field decays as the magnetic and velocity fields both decay. The system is dominated at late times by spanwise modes, which strongly resemble the magnetic field-aligned filamentary flows characteristic of many surges.

One of the long-standing uncertainties in the wave-resonance theory of coronal heating is the stability of the resonance layer. The wave motions in the resonance layer produce highly localized shear flows which vary sinusoidally in time with the resonance period. This configuration is potentially susceptible to the Kelvin-Helmholtz instability (KHI), which can enhance small-scale structure and turbulent broadening of shear layers on relatively rapid ideal timescales. We have investigated numerically the response of a characteristic velocity profile, derived from resonance absorption models, to finite fluid perturbations comparable to photospheric fluctuations. We find that the KHI primarily should affect long (approximately greater than 6 x 10(exp 4) km) loops where higher velocity flows (M approximately greater than 0.2) exist in resonance layers of order 100 km wide. There, the Kelvin-Helmholtz growth time is comparable to or less than the resonance quarter-period, and the potentially stabilizing magnetic effects are not felt until the instability is well past the linear growth stage. Not only is the resonance layer broadened by the KHI, but also the convective energy transport out of the resonance layer is increased, thus adding to the efficiency of the wave-resonance heating process. In shorter loops, e.g., those in bright points and compact flares, the stabilization due to the magnetic field and the high resonance frequency inhibit the growth of the Kelvin-Helmholtz instability beyond a minimal level.

A series of hydrodynamic numerical simulations has been used to investigate the nonlinear evolution of driven, subsonic velocity shears under a range of typical photospheric conditions. These calculations show that typical photospheric flows are susceptible to the Kelvin-Helmholtz instability (KHI), with rapid nonlinear growth times that are approximately half of a typical granule lifetime. The KHI produces vortical structures in intergranule lanes comparable to a typical fluxule radius; this is precisely the correct scale for maximum power transfer to the corona.

A 2.5D numerical code is used to investigate the results of an asymmetric shear imposed on a potential quadrupolar magnetic field under two sets of atmospheric boundary conditions - a low-beta plasma with line tying at the base, similar to the line-tied analytic model, and a hydrostatic-equilibrium atmosphere with solar gravity, typical of the observed photosphere-chromosphere interface. The low-beta simulation confirms the crucial role of the line-tying assumption in producting current sheets. The effects of a symmetric shear on the same hydrostatic-equilibrium atmosphere is examined, using more grid points to improve the resolution of the current structures which form along the flux surfaces. It is found that true current sheets do not form in the corona when a more realistic model is considered. The amount of Ohmic dissipation in the thick currents is estimated to be two to four orders of magnitude below that required to heat the corona. It is concluded that magnetic topologies of the type examined here do not contribute significantly to coronal heating.

Physical parameters obtained for a flaring solar atmosphere in an earlier paper are used here to predict time-dependent emission-line profiles and integrated intensities as a function of position for two spectral lines commonly observed during solar flares: the X-ray resonance lines of Ca XIX and Mg XI. Considerations of ionization nonequilibrium during the rise phase of the flare are addressed, and the effects on the predicted spectral-line characteristics are discussed. It is concluded that some spectroscopic diagnostics favor the nonlocal model, but other long-standing discrepancies between the numerical models and the observations remain unresolved.

Ways in which the radiation-driven thermal instability might affect the structure of the solar atmosphere are considered. It is found that the ultimate state of the medium is highly sensitive to the evolving modal content of the perturbation in that both the initial modal composition and the extent of mode coupling determine the final structure of the atmosphere. It is also found that the condensation process generates highly rotational flows during and after the transition to a new stable state.

A flaring solar atmosphere is modeled assuming classical thermal transport, locally limited thermal transport, and nonlocal thermal transport. The classical, local, and nonlocal expressions for the heat flux yield significantly different temperature, density, and velocity profiles throughout the rise phase of the flare. Evaporation of chromospheric material begins earlier in the nonlocal case than in the classical or local calculations, but reaches much lower upward velocities. Much higher coronal temperatures are achieved in the nonlocal calculations owing to the combined effects of delocalization and flux limiting. The peak velocity and momentum are roughly the same in all three cases. A more impulsive energy release influences the evolution of the nonlocal model more than the classical and locally limited cases.