We present imaging and spectroscopic observations from the Interface Region Imaging Spectrograph (IRIS) of the evolution of the flare ribbon in the SOL2014-04-18T13:03 M-class flare event, at high spatial resolution and time cadence. These observatio
ns reveal small-scale substructure within the ribbon, which manifests as coherent quasi-periodic oscillations in both position and Doppler velocities. We consider various alternative explanations for these oscillations, including modulation of chromospheric evaporation flows. Among these we find the best support for some form of wave localized to the coronal current sheet, such as a tearing mode or Kelvin-Helmholtz instability.
Magnetic energy released in the corona by solar flares reaches the chromosphere where it drives characteristic upflows and downflows known as evaporation and condensation. These flows are studied here for the case where energy is transported to the c
hromosphere by thermal conduction. An analytic model is used to develop relations by which the density and velocity of each flow can be predicted from coronal parameters including the flares energy flux $F$. These relations are explored and refined using a series of numerical investigations in which the transition region is represented by a simplified density jump. The maximum evaporation velocity, for example, is well approximated by $v_esimeq0.38(F/rho_{co,0})^{1/3}$, where $rho_{co,0}$ is the mass density of the pre-flare corona. This and the other relations are found to fit simulations using more realistic models of the transition region both performed in this work, and taken from a variety of previously published investigations. These relations offer a novel and efficient means of simulating coronal reconnection without neglecting entirely the effects of evaporation.
In the present work we study evolution of magnetic helicity in the solar corona. We compare the rate of change of a quantity related to the magnetic helicity in the corona to the flux of magnetic helicity through the photosphere and find that the two
rates are similar. This gives observational evidence that helicity flux across the photosphere is indeed what drives helicity changes in solar corona during emergence. For the purposes of estimating coronal helicity we neither assume a strictly linear force-free field, nor attempt to construct a non-linear force-free field. For each coronal loop evident in Extreme Ultraviolet (EUV) we find a best-matching line of a linear force-free field and allow the twist parameter alpha to be different for each line. This method was introduced and its applicability was discussed in Malanushenko et. al. (2009). The object of the study is emerging and rapidly rotating AR 9004 over about 80 hours. As a proxy for coronal helicity we use the quantity <alpha_i*L_i/2> averaged over many reconstructed lines of magnetic field. We argue that it is approximately proportional to flux-normalized helicity H/Phi^2, where H is helicity and Phi is total enclosed magnetic flux of the active region. The time rate of change of such quantity in the corona is found to be about 0.021 rad/hr, which is compatible with the estimates for the same region obtained using other methods Longcope et. al. (2007), who estimated the flux of normalized helicity of about 0.016 rad/hr.
The Minimum Current Corona (MCC) model provides a way to estimate stored coronal energy using the number of field lines connecting regions of positive and negative photospheric flux. This information is quantified by the net flux connecting pairs of
opposing regions in a connectivity matrix. Changes in the coronal magnetic field, due to processes such as magnetic reconnection, manifest themselves as changes in the connectivity matrix. However, the connectivity matrix will also change when flux sources emerge or submerge through the photosphere, as often happens in active regions. We have developed an algorithm to estimate the changes in flux due to emergence and submergence of magnetic flux sources. These estimated changes must be accounted for in order to quantify storage and release of magnetic energy in the corona. To perform this calculation over extended periods of time, we must additionally have a consistently labeled connectivity matrix over the entire observational time span. We have therefore developed an automated tracking algorithm to generate a consistent connectivity matrix as the photospheric source regions evolve over time. We have applied this method to NOAA Active Region 11112, which underwent a GOES M2.9 class flare around 19:00 on Oct.16th, 2010, and calculated a lower bound on the free magnetic energy buildup of ~8.25 x 10^30 ergs over 3 days.
We show, through a simple patchy reconnection model, that retracting reconnected flux tubes may present elongated regions relatively devoid of plasma, as well as long lasting, dense central hot regions. Reconnection is assumed to happen in a small pa
tch across a Syrovatskii (non-uniform) current sheet (CS) with skewed magnetic fields. The background magnetic pressure has its maximum at the center of the CS plane, and decreases toward the edges of the plane. The reconnection patch creates two V-shaped reconnected tubes that shorten as they retract in opposite directions, due to magnetic tension. One of them moves upward toward the top edge of the CS, and the other one moves downward toward the top of the underlying arcade. Rotational discontinuities (RDs) propagate along the legs of the tubes and generate parallel super-sonic flows that collide at the center of the tube. There, gas dynamics shocks that compress and heat the plasma are launched outwardly. The descending tube moves through the bottom part of the CS where it expands laterally in response to the background magnetic pressure. This effect may decrease plasma density by 30 % to 50 % of background levels. This tube will arrive at the top of the arcade that will slow it down to a stop. Here, the perpendicular dynamics is halted, but the parallel dynamics continues along its legs; the RDs are shut down, and the gas is rarified to even lower densities. The hot postshock regions continue evolving, determining a long lasting hot region on top of the arcade. We provide an observational method based on total emission measure and mean temperature, that indicates where in the CS the tube has been reconnected.
Non-linear force-free fields are the most general case of force-free fields, but the hardest to model as well. There are numerous methods of computing such fields by extrapolating vector magnetograms from the photosphere, but very few attempts have s
o far made quantitative use of coronal morphology. We present a method to make such quantitative use of X-Ray and EUV images of coronal loops. Each individual loop is fit to a field line of a linear force-free field, allowing the estimation of the field lines twist, three-dimensional geometry and the field strength along it. We assess the validity of such a reconstruction since the actual corona is probably not a linear force-free field and that the superposition of linear force-free fields is generally not itself a force-free field. To do so, we perform a series of tests on non-linear force-free fields, described in Low & Lou (1990). For model loops we project field lines onto the photosphere. We compare several results of the method with the original field, in particular the three-dimensional loop shapes, local twist (coronal alpha), distribution of twist in the model photosphere and strength of the magnetic field. We find that, (i) for these trial fields, the method reconstructs twist with mean absolute deviation of at most 15% of the range of photospheric twist, (ii) that heights of the loops are reconstructed with mean absolute deviation of at most 5% of the range of trial heights and (iii) that the magnitude of non-potential contribution to photospheric field is reconstructed with mean absolute deviation of at most 10% of the maximal value.
In this paper we propose that additive self helicity, introduced by Longcope and Malanushenko (2008), plays a role in the kink instability for complex equilibria, similar to twist helicity for thin flux tubes (Hood and Priest (1979), Berger and Field
(1984)). We support this hypothesis by a calculation of additive self helicity of a twisted flux tube from the simulation of Fan and Gibson (2003). As more twist gets introduced, the additive self helicity increases, and the kink instability of the tube coincides with the drop of additive self helicity, after the latter reaches the value of $H_A/Phi^2approx 1.5$ (where $Phi$ is the flux of the tube and $H_A$ is additive self helicity). We compare additive self helicity to twist for a thin sub-portion of the tube to illustrate that $H_A/Phi^2$ is equal to the twist number, studied by Berger and Field (1984), when the thin flux tube approximation is applicable. We suggest, that the quantity $H_A/Phi^2$ could be treated as a generalization of a twist number, when thin flux tube approximation is not applicable. A threshold on a generalized twist number might prove extremely useful studying complex equilibria, just as twist number itself has proven useful studying idealized thin flux tubes. We explicitly describe a numerical method for calculating additive self helicity, which includes an algorithm for identifying a domain occupied by a flux bundle and a method of calculating potential magnetic field confined to this domain. We also describe a numerical method to calculate twist of a thin flux tube, using a frame parallelly transported along the axis of the tube.
