Issue 
A&A
Volume 519, September 2010



Article Number  A44  
Number of page(s)  14  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201014416  
Published online  09 September 2010 
Testing magnetofrictional extrapolation with the TitovDémoulin model of solar active regions
G. Valori^{1}  B. Kliem^{2,3,4}  T. Török^{1}  V. S. Titov^{5}
1  LESIA, Observatoire de Paris, CNRS, UPMC, Université Paris Diderot, 5 place Jules Janssen, 92190 Meudon, France
2 
Universität Potsdam, Institut für Physik und Astronomie, 14482 Potsdam, Germany
3 
University College London, Mullard Space Science Laboratory, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK
4 
Naval Research Laboratory, Space Science Division, Washington, DC 20375, USA
5 
Predictive Science, Inc., 9990 Mesa Rim Road, Suite 170, San Diego, CA 921212910, USA
Received 12 March 2010 / Accepted 28 April 2010
Abstract
We examine the nonlinear magnetofrictional extrapolation scheme
using the solar active region model by Titov and Démoulin as test
field. This model consists of an arched, linetied current channel held
in forcefree equilibrium by the potential field of a bipolar flux
distribution in the bottom boundary. A modified version with a
parabolic current density profile is employed here.
We find that the equilibrium is reconstructed with very high accuracy
in a representative range of parameter space, using only the vector
field in the bottom boundary as input. Structural features formed in
the interface between the flux rope and the surrounding arcade 
``hyperbolic flux tube'' and ``bald patch separatrix surface''  are
reliably reproduced, as are the flux rope twist and the energy and
helicity of the configuration. This demonstrates that forcefree fields
containing these basic structural elements of solar active regions can
be obtained by extrapolation. The influence of the chosen initial
condition on the accuracy of reconstruction is also addressed,
confirming that the initial field that best matches the external
potential field of the model quite naturally leads to the best
reconstruction. Extrapolating the magnetogram of a TitovDémoulin
equilibrium in the unstable range of parameter space yields a sequence
of two opposing evolutionary phases, which clearly indicate the
unstable nature of the configuration: a partial buildup of the flux
rope with rising free energy is followed by destruction of the rope,
losing most of the free energy.
Key words: magnetic fields  Sun: corona  Sun: surface magnetism
1 Introduction
The computation of the coronal magnetic field from boundary data is the only available technique to obtain fully threedimensional information, i.e., the complete field vector in the whole coronal volume bounded below by the magnetogram. Deductions of the field from spectral lines formed in the coronal plasma, from radio maps across a range of microwave frequencies, or from observations of specific structures, such as prominences, are restricted to twodimensional projections or to cuts through the volume, or else they do not yield the full vector. The threedimensional information is required to advance a wide array of issues in coronal physics, for example, the structure of prominences, the onset of eruptions, and the way the coronal field is connected to the Sun's interior and the solar wind.
In equilibrium, the coronal magnetic field is nearly forcefree
in the major part of active regions, except their upper periphery (see,
e.g., Gary 2001; Metcalf et al. 1995; Moon et al. 2002).
To a very good approximation it can therefore be described by
Here is a scalar function that is constant along each individual field line (which follows directly from Eqs. (1)) but in general has different values on different field lines.
The computation of forcefree coronal fields can be formulated as a boundaryvalue problem for Eqs. (1). Such extrapolation can be performed at different levels of approximation for the field. If is assumed to be constant, then the first of Eqs. (1) is linear, and only the knowledge of the normal field component is needed as boundary condition for the extrapolation. The solution of the linear problem was given in closed form in Nakagawa & Raadu (1972), Chiu & Hilton (1977), and (). Linear methods have several intrinsic limitations, the most severe one being that they cannot match most observed magnetograms closely because is often found to vary strongly across active regions (e.g., Régnier et al. 2002). The particular case of yields a potential field, which is inadequate for many purposes, since it does not contain any free energy to power coronal activity. Therefore, it is necessary to proceed beyond the linear approximation and to solve the nonlinear extrapolation problem, permitting to vary across the field (Amari et al. 1999; Sakurai 1981; McClymont et al. 1997). The nonlinearity means that such extrapolation requires a numerical approach. Various numerical schemes and implementations have been developed to construct a solution; for a recent review, see Schrijver et al. (2006).
A partially alternative strategy, the flux rope insertion method (van Ballegooijen 2004; van Ballegooijen et al. 2007), permits structures to be modeled that are supposed to contain a flux rope, particularly filament channels. Here an extrapolated potential field is replaced by a flux rope in part of the volume, and the resulting configuration is numerically relaxed to find an equilibrium. The parameters of the inserted rope are varied until the field lines of the resulting equilibrium match observed features, e.g., threads in a filament or overlying loops.
A convenient tool for both nonlinear forcefree extrapolation and the flux rope insertion method is the socalled magnetofrictional relaxation of the field to forcefree numerical equilibria (Craig & Sneyd 1986; Roumeliotis 1996; Chodura & Schlüter 1981; Yang et al. 1986). In a previous paper (Valori et al. 2005), we described our implementation of magnetofrictional nonlinear extrapolation and its application to reconstructing a relatively simple forcefree magnetic field, which contains an approximately currentneutralized flux rope rooted in the flux concentrations of a simple bipolar magnetogram (Török & Kliem 2003). That paper included a thorough study of the dependence of the reconstruction quality on the numerical resolution. In Valori et al. (2007), an improved implementation of the magnetofrictional method allowed us to significantly increase the accuracy in the reconstruction of the well known forcefree equilibria by Low & Lou (1990) in comparison to similar extrapolations summarized in Schrijver et al. (2006). Applications of this code to measured magnetograms are included in Metcalf et al. (2008), Schrijver et al. (2008a), and DeRosa et al. (2009).
Here we extend the testing of our magnetofrictional code in the reconstruction of solarrelevant solutions of Eqs. (1) by applying it to the flux rope equilibrium constructed by Titov & Démoulin (1999, hereafter TD). This configuration describes structures that emerge from the solar interior already twisted and that carry a net current (e.g., Wheatland 2000; Lites et al. 1995; Leka et al. 1996), and it is in this sense complementary to the one in Valori et al. (2005), which represents the case of an originally potential field twisted solely by photospheric displacements. The relevance of the TD equilibrium as a model for solar active regions and eruptive configurations has been demonstrated by many investigations (e.g., Zarro et al. 2000; McKenzie & Canfield 2008; Török & Kliem 2005; Schrijver et al. 2008b; Kliem et al. 2004). In particular, it has been shown that the TD flux rope can be subject to ideal MHD instabilities (Török & Kliem 2007; Török et al. 2004; Isenberg & Forbes 2007; Roussev et al. 2003). Depending on its parameters, the equilibrium can include a hyperbolic flux tube (HFT) or bald patches (BPs) (see Sect. 2 for detail). Such structural features are thought to play an important role in the initiation and evolution of coronal mass ejections (CMEs) and their associated flares. Moreover, the occurrence of BPs changes the topology of the field.
Our goal in this paper is to provide a detailed analysis of the reconstruction capabilities of the magnetofrictional code for this test field. We address the flux rope morphology, including the twist and the structural features HFT and BP, the energy and helicity contents, and the influence of the initial condition used by the extrapolation code. We also applied the code to an unstable case, expecting this to be useful in judging extrapolations of magnetograms that were taken during the initial phase of an eruption or immediately prior to its onset. Finally, we compare our results to the previous reconstruction of the TD field by Wiegelmann et al. (2006a), which employs a different extrapolation scheme.
Using a forcefree test field implies that the magnetogram is compatible with the forcefree hypothesis of the extrapolation method. Therefore, in the following we do not need to consider deviations from forcefreeness in the photosphere, the ambiguity of the transverse field direction, noise, and insufficient resolution, all of which are associated with measured magnetograms. Attempts to deal with some of these problems are discussed, e.g., in McClymont et al. (1997), Amari et al. (1997), Wiegelmann et al. (2006b), Fuhrmann et al. (2007), Metcalf et al. (2008), Wiegelmann et al. (2008), Schrijver et al. (2008a), and DeRosa et al. (2009).
The paper is organized as follows. The TD equilibria used below are described in Sect. 2, the extrapolation method is briefly outlined in Sect. 3, and the results are presented and discussed in Sect. 4. Our conclusions are summarized in Sect. 5.
2 Test equilibria
Titov & Démoulin (1999, TD) constructed a forcefree equilibrium of a toroidal current channel to obtain an analytical model of a solar active region that contains free magnetic energy. Because it is force free and carries a net current, the equilibrium belongs to the class of tokamak equilibria, also known as Shafranov equilibria (Shafranov 1966), which require an external poloidal field, i.e., one due to sources other than the current channel, to balance the always outwardly directed Lorentz selfforce (hoop force) of the bent current channel. TD used a pair of magnetic sources, placed symmetrically to the torus plane at the symmetry axis of the torus, to model the external poloidal field, , of a solar active region. They also added a line current running along the symmetry axis to include an external toroidal field, , which permits modeling of the magnetic shear typically present in active regions and prevents the twist profile from becoming infinite at the surface of the torus. The torus is arranged vertically, with its symmetry axis submerged below the plane , which represents the photosphere, resulting in a smooth, divergencefree field in the coronal volume . Inside the torus, the field is obtained by modifying the forcefree field of a locally straight current channel of prescribed radial profile to match the external potential field smoothly at the surface. The original version of the equilibrium uses a uniform radial profile of the current density, resulting, after the matching, in a profile that increases moderately toward the surface at , where it drops to zero [ ]. The requirement of forcefreeness implies that a toroidal field component is always present inside the torus, regardless of the value of . The forcefreeness within the current channel (i.e., the equilibrium condition) is satisfied only approximately, but numerical relaxation runs verify that the approximation is very good for a high toroidal aspect ratio, (Kliem et al. 2004), and tolerable down to moderate and likely more realistic aspect ratios of about (Schrijver et al. 2008b).
If characterized by its field lines, the equilibrium represents an arched, linetied flux rope that is embedded in a (generally sheared) potential field of arcade topology. The flux rope has a somewhat larger crosssection than the current channel. The ``lineofsight'' magnetogram, B_{z}(x,y,0), possesses four flux concentrations, a pair of ``sunspots'' above the point sources, which are located at , and a pair of ``satellite polarities'' in the intersection between the current channel and the ``photosphere'' (Fig. 1). Depending on the parameters chosen for the particular realization of the equilibrium, the satellite polarities can actually encompass more flux than the sunspots. The neutral line of the normal field component in the magnetogram, given by B_{z}(x,y,0)=0, runs approximately parallel to the magnetic flux rope axis beneath its apex and runs relatively straight between the two positive and the two negative flux concentrations; i.e., the ``active region'' is to be characterized as essentially bipolar^{}.
Figure 1: ``Lineofsight'' magnetogram B_{z}(x,y,0) with superimposed arrows representing the transverse field ( top), and forcefree parameter ( bottom), for the reference fields of the cases High_HFT ( left) and BP ( right). 

Open with DEXTER 
Figure 2: Contours of , with superimposed arrows representing the transverse field in the vicinity of the neutral line of B_{z}, for the reference ( top) and reconstructed ( bottom) BP cases. Note the inverse crossing of the neutral line, the characteristic signature of a BP, in the middle of the plots. 

Open with DEXTER 
The bipolar structure of the external field implies the absence of null points. These relevant topological features (e.g., Antiochos 1998) can be included by generalizing the external field; however, we leave this step for future work.
With regard to the reconstruction task, it is important to note that TD equilibria that contain a separatrix are structurally quite different from any potential field that can be constructed from their corresponding lineofsight magnetogram. The potential field will not include a flux rope. In the majority of cases, HFT or BP(s) will not be present either, and in case they happen to exist, they will have very different locations than in the TD field. At the same time, these equilibria appear to be the simplest forcefree configurations currently considered in coronal physics, which are structurally (and, in the presence of BP(s), even topologically) different from a shearedarcade field. They consist of a single current channel and the minimum two flux systems required: the flux in the rope and the external poloidal flux. These are the generic elements of more complex forcefree configurations that carry a net current. The requirement of a structural or even topological change makes the reconstruction qualitatively different from reconstructions of Low & Lou fields or of fields obtained by shearing or twisting the photospheric flux with smooth profiles of substantial symmetry (e.g., Török & Kliem 2003). Therefore, the reconstruction of the TD field is a major and necessary step in verifying the capabilities of nonlinear extrapolation schemes.
Figure 3: Selected field lines for the Low_HFT case; a) reference, b) extrapolated, c) potential, and d) linear field (with ). In a) and b) the flux rope is depicted by two groups of field lines (in redorange and green), each one starting from six uniformly spaced points on a circle of radius 0.1, centered on the intersections of the magnetic axis of the flux rope (CFL) with the z=0 plane. Gray field lines are started from the zaxis and visualize the potential field above and below the current channel. In c) and d) the field lines are started from the same set of footpoints as in a) and b). Contour plots of B_{z}(x,y,0) and its neutral line B_{z}(x,y,0)=0 (dotdashes) are shown in the bottom plane of each panel. 

Open with DEXTER 
We note that the configuration considered in Metcalf et al. (2008) represents an intermediate case between purely sheared configurations (e.g., the family of Low & Lou fields) and the TD field. It was obtained by means of the flux rope insertion method such that the magnetic axis of the rope passes at a low height of only a few Mm above the photosphere along much of the length of the rope (van Ballegooijen 2004). This corresponds to a strongly submerged TD rope with Ra<d<R+a and contains a separatrix that may or may not form an HFT or a BP. The configuration used in Metcalf et al. (2008) includes a BP, but the extrapolations presented in that paper were not checked for the reconstruction of this feature.
We have searched published extrapolations of solar magnetogram data for the occurrence of BPs (inverse field direction at the neutral line) and of the flux penetration feature characteristic of an HFT. There appear to only be two cases of a BP (Guo et al. 2010; Canou et al. 2009) and one case of flux penetration (Régnier & Amari 2004).
For the purposes of the present paper, the original TD equilibrium is modified in two ways. First, in place of the line current, we use a pair of dipoles as sources of the external toroidal field (Török & Kliem 2005). This yields a compact distribution of magnetic flux in the magnetogram plane and a steeper decrease of the external field with height than for the original TD field, both of which are more realistic with regard to solar active regions. The dipoles are placed below the footpoints of the flux rope and point nearly vertically, such that the field lines connecting them through the torus section in match the geometry of the torus section closely. The conservation of toroidal flux in the rope now requires the minor radius of the torus to vary along its length. Since still points in the approximate direction of the current, its strength remains a free parameter to a good approximation for parameters of interest here, similar to the original equilibrium. The break of the toroidal symmetry also causes the magnetic axis of the rope and the HFT to bulge out of the plane , but the effect remains weak, since varies only moderately along the coronal flux rope for the chosen considerable depth of the dipoles. Second, the radial profile of the current density is changed such that it no longer peaks at the surface but in the middle of the channel, and decreases monotonically towards zero at the surface. This removes unrealistically steep gradients of the current density at the surface of the torus, including the magnetogram plane. The expressions for the equilibrium with this parabolic current density profile will be provided in a future publication (Titov, in preparation).
Table 1: Normalized parameters defining the TD equilibria.
We consider the modified TD equilibrium for several parameter sets, four stable configurations and an unstable one. The parameters defining each equilibrium are given in Table 1. All variables are normalized by a characteristic length, taken here to be the apex height of the geometrical torus axis, Rd, by the magnetic field strength at the apex point, and by quantities derived from them.Three of the four stable equilibria have identical geometric flux rope parameters and an almost identical lefthanded, average endtoend twist of in the center of the current channel (see Sect. 4 for the averaging procedure), but different properties of the ambient field. In the first two cases an HFT is present, staying close to the photosphere in one (Low_HFT) and reaching significantly into the volume in the other (High_HFT). A third case (No_HFT) has no HFT above the photospheric plane, so it has a simpler magnetic structure than the cases with an HFT. Owing to the relatively strong toroidal field component chosen, these three configurations do not have bald patches. The fourth stable case (BP) has a lefthanded average twist of , close to the twist of the first three equilibria, but the minor radius of the torus is enlarged (and correspondingly reduced) to introduce a bald patch in the resulting field.
Figure 1 illustrates for the High_HFT and BP cases that the field strength in the magnetogram decreases rapidly with distance from the four flux concentrations. The distribution of in the same figure shows that the field is highly nonlinear (see below for the occurrence of currents outside the torus in these plots). The transverse field near the section of the neutral line between the flux concentrations is strongly aligned with the neutral line for the BP configuration, as in a filament channel on the Sun. A zoom into this area is plotted in Fig. 2, demonstrating by the inverse crossing of the neutral line in the central part of the magnetogram that a BP is indeed present.
To obtain an unstable case, we increase the average twist to by reducing the minor torus radius. This configuration is unstable with respect to the helical kink instability.
The analytical expressions in TD and their modifications described above represent only approximate equilibria. Moreover, the numerical determination of the current density introduces additional (``discretization'') errors. Therefore, the individual magnetic configurations were first relaxed to numerical equilibria, using the MHD code described in Török & Kliem (2003) on uniform Cartesian grids with a resolution in all directions. This introduces very little change in the shape of the flux rope, but new layers of current density, not present in the analytical model, form in the separatrix and its vicinity, primarily below the current channel. These are visible in the plots of Fig. 1 as narrow strips of forward and reverse current between the footpoints of the channel.
From each numerically relaxed stable TD equilibrium, a data cube of size was extracted, containing the nonlinear magnetic field to be used as the test field for the extrapolation code (the somewhat larger cube used for the unstable case is given in Sect. 4.3). We refer to these data cubes as the reference field, , for the given set of parameters. The reference fields were then reconstructed by magnetofrictional extrapolation, employing the same grid as in the MHD relaxation. The High_HFT reference field was also constructed and reconstructed using a uniform grid resolution of 0.03 in each direction. Since the achieved accuracy of reconstruction turned out to be comparable to the one obtained with a resolution of 0.06, we selected the latter value for the present study. Although the nonlinear solution is known in the whole discretized volume, and specifically on its lateral and top boundaries, the vector field in the bottom plane, , is the only information used in the reconstruction of all test fields, because this is the challenge that the extrapolation codes have to face when they are applied to observed magnetograms.
Since the reference TD fields satisfy Eqs. (1) only to the accuracy of the finitedifference scheme used in the relaxation, it is necessary to quantify the compatibility of their magnetograms with the forcefree hypothesis implicit in the extrapolation method. For this purpose, we use normalized integral averages of flux, force, and torque in the lower boundary as defined in Wiegelmann et al. (2006b), which vanish for perfect compatibility. These metrics, reported in Table 2, show that the flux imbalance is negligible and that the residual force and torque values remain small, even for the unstable case.
3 Extrapolation method
The implementation of magnetofrictional extrapolation used in this paper was described in detail in Valori et al. (2005,2007), so that we can restrict ourselves here to a short outline of the method and its recent improvements. Our formulation of the method consists in a pseudotemporal evolution of an initial field subject to the equation
and c_{Y} and
are numerical factors fixed by suitable stability criteria.
The initial field is usually chosen to be a potential field derived from the ``lineof sight'' magnetogram
,
with the full vector magnetogram
overwritten in the photospheric layer.
The rightmost term in Eq. (2) is the
cleaner (see, e.g., Keppens et al. 2003, and references therein).
By applying the divergence operator, it is readily seen that Eq. (2) reduces to a diffusion equation for the scalar field
.
Since also the first term on the righthand side of Eq. (2) is of parabolic type (Craig & Sneyd 1986), Lorentz forces and errors of
diffuse out of the computation domain through the boundaries in the course of the extrapolation.
This is best seen by computing the time derivative of the magnetic energy
.
Multiplying Eq. (2) with
and integrating by parts, it is found that
which implies that the magnetic energy decreases as long as there are Lorentz forces and magnetic sources in the domain. Here it has been assumed that the boundaries are perfectly conducting (as in Craig & Sneyd 1986) and that the field is solenoidal at them, so that surface integrals do not contribute. During the pseudotemporal evolution, the decrease in the perpendicular current density also reduces the relaxation velocity (Eq. (3)). If a static state is reached, then the field is divergencefree and implies that it is also forcefree. By construction, such a static solution extends smoothly from the vector magnetogram in the bottom boundary and represents the extrapolated nonlinear forcefree field. By overwriting the vector magnetogram in , we actually introduce a finite jump in the field close to the photosphere, yielding Lorentz forces and errors into the solenoidal property; consequently, surface integral terms in Eq. (4) are important initially. This results in an initial rise in the energy not described by Eq. (4). A discussion of this effect will be given in a future paper.
Table 2: Normalized flux, force, and torque imbalance of the five analyzed magnetograms.
Equation (2) can be formally derived from the MHD equations (as in Yang et al. 1986) by assuming that viscosity dominates the evolution (hence the characterization of the method as ``magnetofrictional''). However, since such a strong viscosity has no immediate physical justification for the nearly collisionless coronal plasma, one can equivalently consider Eq. (2) to be a numerical tool that finds a forcefree magnetic field for given boundary conditions, with only the initial potential field (before the vector magnetogram is overwritten) and the final nonlinear forcefree field having a welldefined physical meaning.
Our implementation of the magnetofrictional code employs a fourthorder centraldifference scheme for the space discretization. The pseudotime discretization is based on a general MHDstability prescription (see details in Valori et al. 2007). To reduce the execution time, the time discretization was recently complemented by a RungeKuttaChebyshev acceleration technique (Evje et al. 2000; Alexiades et al. 1996). The detailed description and testing of this upgraded pseudotime discretization is beyond the scope of the present paper and will be provided in a future publication.
The field in the photospheric boundary is fixed to keep the input magnetogram values throughout the whole extrapolation. On the side and top boundaries of the extrapolation volume, which is identical to the volume of the reference field given in Sect. 2, the normal component is derived from the inner field values so as to have vanishing in the boundary. The transverse component is given by a fourthorder polynomial extrapolation of the inner field values onto the boundaries. These boundary conditions were tested in Valori et al. (2007), where more detail can be found. Finally, all extrapolations presented here start from the potential field obtained by Schmidt's method (Sakurai 1989; Schmidt 1964), except for a potentialfield and a linearfield case discussed in Sect. 4.2, for which the method by Seehafer (1978) was used.
4 Analysis of the extrapolation results
The reconstructed fields, , are compared to their respective reference field, , in an ``analysis domain'' given by clipping 20 grid layers from the top and lateral boundaries of the extrapolation volume. This is done in recognition of the fact that the magnetograms have considerable margins of currentfree field around the flux rope footpoints. The currents remain relatively small in the corresponding parts of the extrapolation volume, where the relaxation is slowed, so the quality of reconstruction is reduced. Similar practice, with larger margins, was adopted in Schrijver et al. (2006) and Metcalf et al. (2008). The quality of the reconstructed equilibria is quantified in the following in several ways.
The degree of forcefreeness is measured by the currentweighted, average sineangle between current and magnetic field,
and the sums are taken over all grid nodes i in the analysis domain. This same quantity is sometimes referred to as CWsin (e.g., Schrijver et al. 2008a). The solenoidal property is quantified, following Wheatland et al. (2000), by the average over the grid nodes of the fractional flux
through the surface of a small volume including the node i. For uniform grid spacing , this reduces to and , where N is the number of grid points^{}. The influence of the residual spurious magnetic charges implied by can also be quantified by the nonsolenoidal part of the magnetic field, . Using the Gauss theorem, we have
which is evaluated on an auxiliary grid displaced by in each direction. All three quantities should be as small as possible. We have demonstrated in Valori et al. (2007) that our code reduces the value of to machine precision in some cases; however, due to the parabolic nature of Eq. (2), this may require long execution times. For the purposes of the present work, we terminated the relaxation when the forcefreeness metric tended to reach a plateau of sufficiently low value (10^{2}) and was within the same order of magnitude as the corresponding value for the initial potential field or less.
Since the reference field
is known, we can directly compare
and the reconstructed field
obtained from the magnetogram.
Several comparison metrics were defined in Schrijver et al. (2006), out of which we adopt the mean vector error
The smaller , the better the reconstruction. As for previous reconstructions of test fields (Schrijver et al. 2006; Valori et al. 2007; Metcalf et al. 2008), this metric is the one most sensitive to reconstruction errors in the TD field. The further metrics of direct comparison defined in Schrijver et al. (2006), , , , yield even smaller deviations from the ideal value (1 or 0) than for all cases analyzed here.
To judge geometry and topology, we first determine the apex height of the flux rope's magnetic axis, also referred to as the ``central field line'', CFL, in the following. (The CFL differs slightly from the geometrical axis of the torus, since the external potential field drops with height, thereby contributing differently to the bottom and top parts of the flux rope.) Next we check for the presence of an HFT and for its apex height. Both apex heights are given by inversion points of the poloidal component of the field at the linesymmetric z axis. This is approximated well by B_{x}(0,0,z) because the flux rope writhes only rather weakly out of the plane in the equilibria considered in this paper. Finally, we check whether the BP in the fourth stable test field is recovered at the correct section of the neutral line.
The TD equilibrium can be susceptible to the helical kink instability (Török et al. 2004) and to the torus instability (Kliem & Török 2006).
The parameters controlling these instabilities should be reliably and accurately reproduced by the extrapolation.
The main parameter that regulates the stability with respect to the helical kink mode is the field line twist.
This quantity,
,
is defined with respect to the CFL of length l using local cylindrical coordinates
.
We use its average in a circular cross section of normalized radius 0.3 centered on the apex of the CFL.
Within this radius, the average twist is
for three of the four stable equilibria (High_HFT, Low_HFT, and No_HFT) and
for the fourth stable case (BP).
For the unstable case discussed in Sect. 4.3, the average twist in
is .
The threshold of the torus instability is given in terms of the ``decay index'' of the external poloidal field,
,
where R is the major torus radius.
For
,
with the threshold lying in the range
(Kliem & Török 2006), an expansion of the current ring out of equilibrium
cannot be balanced by the external poloidal field.
Computing the index n requires knowledge of the external
poloidal field, which is not immediately available for our numerical
solutions, since it cannot be easily separated from the poloidal field
produced by the ring current.
Therefore, we compute the decay index of the total poloidal field at
the z axis, given to a good approximation by
B_{x}(0,0,z),
The average is taken in a small interval in z, which is located just above the current channel in all equilibria considered here. The quantity is not suitable to determining the onset of the torus instability, but it serves the purpose of quantifying how well the stabilizing potential field above the current channel is reproduced by the extrapolation.
The magnetic energy, , and relative magnetic helicity, , should also be accurately reproduced by the reconstructed field. We use the formulae by DeVore (2000) to obtain an approximation of for both the reference and the reconstructed fields.
In order to reliably compare the effects that topological differences have on the extrapolation quality, all computations presented in this paper for stable TD equilibria were performed using identical numerical parameters, rather than being individually optimized for best reconstruction.
Table 3: Reconstruction accuracy of the stable TD equilibria.
4.1 Stable equilibria
We consider the reconstruction of four stable TD equlibria, two of which contain an HFT (of different height) and one contains a BP (Table 1), using the potential field by Schmidt (1964) as initial condition. The main results are summarized in Table 3, where measures of the reconstruction quality of the four reference fields are listed. (See also the comprehensive compilation of the various metrics for all our extrapolations in Table 5.) The table shows that the errors are very small and of comparable magnitude for all cases (within fractions of a percent up to a few percent). In particular, we point out that the mean vector error of all four reconstructed fields is 0.054 or smaller, a result far better than the accuracy achieved in any previous similar reconstruction study, i.e., in those that used only the vector magnetogram as input (Schrijver et al. 2006; Wiegelmann et al. 2006a; Valori et al. 2007; Metcalf et al. 2008), although a smaller margin between the analysis and extrapolation volumes is used here. This mean vector error corresponds to a mean angular deviation between and of at most 3 degrees (if the field strengths agree on all grid nodes) and to a mean deviation of field strength of at most 5% (if all directions agree)  values lower than the noise level for most pixels in observed magnetograms. It is also worth noting that the reconstruction accuracy is now comparable to the one achieved previously in ``extrapolations'' that have provided the reference field vector on all six boundaries of the reconstruction volume. The accuracy reached here is clearly superior to all but one of those results reported in Tables 1 and 2 in Schrijver et al. (2006) and is approaching the best one in this paper, as well as the results in Amari et al. (2006) and Wiegelmann et al. (2006a). The field line plots in Fig. 3 show that such small errors have practically no visible influence on the reconstruction, which is virtually indistinguishable from the reference field.
The best reconstruction is obtained for the case with a bald patch, while the least accurate one is the No_HFT case, where neither an HFT nor a BP are present. There is apparently no relation between the quality of the reconstruction (Table 3) and the compatibility of the corresponding magnetogram with the forcefree constraints (Table 2). This shows that the level of inconsistencies in the vector magnetograms of our stable reference fields is small enough to have essentially no influence on the extrapolation quality.
Figure 4: Plane averages of the relative solenoidal error in the No_HFT case, as a function of position. 

Open with DEXTER 
The geometry of the reference field is perfectly recovered, as expected from the very low value of and illustrated in Fig. 3 for the case Low_HFT. The flux rope, as the entire extrapolated field, reproduces the zaxis line symmetry of the TD equilibrium exactly. Therefore, the flux rope location is correctly represented by the apex height of the CFL, found to match the height in the reference field nearly exactly for all stable equilibria. Since includes a normalization to the local field value on the grid, it gives the weakfield areas in the outer parts of the reconstruction volume the same weight as the strongfield areas near the bottom and in the flux rope. From the very low value of , we can therefore be sure that the high overlying field loops are reconstructed with essentially the same accuracy as the flux rope, and this is apparent from Fig. 3 as well. The comparison of such loops with soft Xray or EUV loops has proven insatisfactory in recent nonlinear extrapolations of observed vector magnetograms (DeRosa et al. 2009), although the observed loops very likely outline forcefree or even currentfree field line bundles. The results given here demonstrate that this is not due to limitations of the extrapolation scheme(s).
Figure 5: Isosurface of the reconstructed current channel at 31% of peak current density for the cases High_HFT ( top) and No_HFT ( bottom). Field lines crossing the z axis are shown in gray if below the HFT and in green if above the HFT. 

Open with DEXTER 
Considering the parameters that characterize the relevant instabilities, we find that the relatively high twist of is recovered with very high accuracy, the errors being below 3% in all four cases. Thus, as in Valori et al. (2005), it is confirmed that the magnetofrictional method can reconstruct field lines that wind around a flux rope axis more than one time. In this regard the nonlinear extrapolation differs in principle from linear extrapolations, which limit the number of field line turns in the computed field to about 1/2 for realistic magnetogram sizes, due to the upper limit on the forcefree parameter inherent in linear fields (Kliem et al., in preparation; see also Leka et al. 2005). This means that estimates of twist in fields extrapolated with the magnetofrictional code are by far more reliable than estimates using a linear code.
Figure 6: Poloidal field at the z axis, B_{x}(0,0,z), in the No_HFT case, for the reference (red), extrapolated (green), and TD_ex potential fields (blue). Left: entire z axis extension; right: upper part (above the current channel) in logarithmic scale. 

Open with DEXTER 
Similarly, the magnetic energy and relative magnetic helicity are obtained with high accuracy; their errors remain below 1 and 5 percent, respectively, for all four stable reference fields.
Figure 3 also shows that the potential and linear extrapolations entirely miss the flux rope, typically leading to an arcadetype structure. Despite the failure of the potential field to reproduce any of the salient features of the reference field, still some of the comparison metrics score decently (e.g., ; see Table 5 for more detail). This is because a large fraction of the analysis volume is occupied by potential field. Because it is the minimum for a given normal field distribution at the boundaries, the energy of the potential field reconstruction is of course lower than the energy of the reference field. However, it still accounts for a large part of it, about 80% in the selected analysis volume. The linear reconstruction suffers from similar limitations but, since it fills the entire volume with currents, it is affected by even larger errors than the potential field. The value adopted for the computation of the linear field is the one that gives the best rms matching of the transverse magnetogram components, often referred to as (see Pevtsov et al. 1995).
Finally, we briefly comment on the reconstruction of the TD equilibrium in Wiegelmann et al. (2006a). The original TD equilibrium with a line current was considered, and a different extrapolation scheme, the optimization method, was used. The paper includes a reconstruction based only on the information from the bottom boundary as ``Case II.'' The corresponding field line plot in Fig. 2 of that article (done by one of the authors of the present paper, GV) indicates that the extrapolated field contains two partially twisted flux tubes rather than one. That plot is actually misleading because the appearance of two separate flux tubes was obtained by calculating field lines starting around the nominal footpoint positions of the geometrical torus axis; however, as described above, the footpoints of the flux rope's magnetic axis do not coincide with those of the geometrical torus axis. A proper field line plot (not included here) is obtained by placing the field line start points around the magnetic axis, and it shows a single, higher, and slightly writhed flux rope, similar to the cases analyzed here. Due to the many differences between the employed test fields and reconstruction volumes, we do not pursue the comparison with the results in Wiegelmann et al. (2006a) any further. Those are far less accurate than the ones given here when only the vector magnetogram is used as input.
4.2 Dependence on the initial field
In this section we describe the link between the potential or linear forcefree field that is used as initial condition for the extrapolation and the quality of the obtained reconstruction. A key consideration here is the balance between the hoop force and the Lorentz force provided by the external poloidal field. We consider the Low_HFT case in combination with four different initial fields: two potential fields computed by the methods of Schmidt (1964) and Seehafer (1978), respectively; a linear forcefree field computed by the Seehafer method; and the external potential field of the TD equilibrium (generated by the magnetic charges and dipoles in the modified version of the equilibrium used in this paper). The last, referred to as TD_ex in the following, is, of course, not available for observed magnetograms.
Figure 7: Poloidal field at the z axis, B_{x}(0,0,z), in the Low_HFT case, for the reference (red solid line), initial (blue), and corresponding extrapolated fields (green), for different type of the initial field: potential Schmidt (short dashes), potential Seehafer (dotdash), linear alphabest Seehafer (triple dotsdash), and potential TD_ex (long dashes). 

Open with DEXTER 
Table 4: Reconstruction accuracy of the Low_HFT equilibrium for different initial conditions.
The Seehafer potential and linear fields are weaker at the nominal CFL apex position, with a ratio of the poloidal components to of the TD_ex field of 0.9 and 0.48, respectively, so that the flux rope finds equilibrium at greater heights. However, the apex heights differ by only 14% and 18%, respectively, which shows that the initial field does have an influence on the extrapolation, but also that the magnetofrictional relaxation can tolerate a substantial mismatch between the initial field and the true external field while still producing relatively small errors. The information flowing from the magnetogram into the extrapolation volume can largely, but not completely, modify the properties of the initial field, although open boundaries are implemented at the top and sides of the box, which allow the field to change in response to changes in the interior. This has to do with the fact that large parts of the volume are currentfree, both in the initial and in the final configuration; hence, currents must be built up and subsequently removed, in order to modify the initial field in the large outer parts of the box. This is a slow process, eventually hindered by numerical diffusion.Apart from the CFL apex height and the related decay index , the figures of merit of the four extrapolations in Table 4 show excellent values. The tendency for the Seehafer potential and linear fields to score worse is visible but rather weak, indicating that the initial field has a weaker influence on quantities like flux rope twist, total energy, and helicity. These depend primarily on the information contained in the magnetogram and only to a smaller extent on the flux rope height. It is also clear that the figures of merit for the Schmidt potential field mostly beat those of TD_ex because the former is more consistent with the magnetogram, which includes the contribution by the current channel. All four extrapolations recover the very low lying HFT apex of the test field, essentially at the correct height (Table 5).
These comparisons clearly suggest using the Schmidt potential field as initial condition for fluxbalanced magnetograms of isolated active regions (two conditions that go together in tendency). A flux imbalance can often be reduced by embedding the vector magnetogram in a larger magnetogram, which is typically of lower resolution and (with current instrumentation) only a lineofsight magnetogram. The transverse components must then be modeled, for example by assuming a potential field, and this can strongly influence the outcome of the extrapolation. Alternatively, the use of the imbalanced vector magnetogram joint with open boundary conditions that allow for flux and currents to leave the numerical box through the top and side boundaries, as in our code, can lead to better extrapolations (see DeRosa et al. 2009; and Fuhrmann et al. 2010, A&A, submitted) for a discussion). Finding the best tradeoff between these options will require further study.
4.3 Unstable equilibrium
As mentioned in the Introduction, the parameters defining the TD equilibrium can be chosen such that the configuration is unstable with respect to the helical kink or to the torus instability. Such equilibria offer the opportunity to test the extrapolation code in the reconstruction of unstable forcefree configurations, addressing the question of which signatures of the unstable situation are produced. In the following we compare the extrapolation of the magnetogram provided by an unstable TD equilibrium with the equilibrium itself and with the posteruption configuration obtained by a numerical MHD evolution of the perturbed equlibrium with the magnetogram kept fixed.
The configuration we consider is the same as the eruptive case in Török & Kliem (2005), except that a smaller initial twist of is used. (The twist average over the whole cross section of the current channel, as calculated by Török & Kliem, is for this configuration.) The MHD simulation is performed in a Cartesian box of size , which consists of the uniform grid used for the extrapolation in its central part and a streched grid outside the extrapolation domain. The smaller twist allows us to relax the configuration, including the magnetogram plane, for a few Alfvén times, reducing the initial spurious Lorentz force densities, before the helical displacement due to the developing kink instability becomes visible.
Figure 8: Selected field lines for the reference fields of the unstable case at the initial ( top) and final ( bottom) stages of the MHD evolution (see Fig. 3 for the color coding of the field lines). 

Open with DEXTER 
To be compatible with fixing the magnetogram, we set the plasma velocities to zero in the grid layers and . Therefore, although the simulation is run for a long time, the Lorentz force densities do not relax completely, so that the posteruption configuration, while largely relaxed, is not perfectly forcefree in the lower layers of the box. Nevertheless, the field line connectivities do not change significantly anymore, so that a qualitative comparison with the forcefree extrapolation is possible. We refer to this posteruption configuration as the final reference field for the extrapolation of the unstable case (see again Fig. 8).
Table 5: Figures of merit (absolute values) for all reference, initial, and extrapolated fields.
The extrapolation of the unstable case is obtained in a similar way to what is described in Sects. 4.1 and 4.2, except that here the extrapolation grid discretizes a larger volume with the same uniform resolution . The extrapolation employs the z=0 magnetogram of the initial/final reference field and starts from the Schmidt potential field in . Table 2 verifies that the magnetogram includes residual forces, but it has a lower value of the applied torque than the stable cases. According to the metrics in the table, the unstable nearequilibrium TD configuration obtained from the initial partial MHD relaxation (before the instabilities become dominant) provides a magnetogram that is indeed acceptably forcefree.Figure 9: Evolution of (solid line) and (dashes) during the extrapolation runs of the cases BP (left) and Unstable (right). The energy is normalized to its value at the first iteration. 

Open with DEXTER 
Figure 10: Selected field lines for the extrapolation of the unstable case at different iterations. From left to right: the initial potential field at it = 1, the maximumenergy configuration at 6700 iterations, and the final lowenergy configuration at iterations. See Fig. 3 for the color coding of the field lines. 

Open with DEXTER 
We recall that the extrapolation starts from a potential field that obviously has no flux rope (as shown in the left panel of Fig. 10), and therefore the flux rope in the maximumenergy reconstruction is formed by the extrapolation process. The extrapolation does not succeed to form the flux rope completely because the underlying TD configuration is in the unstable domain, so that the Lorentz forces tend to propagate the information about the current channel from the magnetogram into the volume and tend to move the system away from that configuration at the same time. The partial buildup of the current channel is also reflected by the maximumenergy field having reached 77% of the energy in the initial reference field. In light of this number, the very close agreement of the magnetic helicities must be considered a coincidence or a consequence of the approximate nature of its calculation (using the expressions in DeVore 2000).
Conversely, the minimumenergy solution has no flux rope: field lines started at the original location of the flux rope connect to the sunspots (rather than to the flux rope footpoints in the TD and initial reference fields), in a manner largely similar to the final reference field (compare the bottom panel in Figs. 8 with the right panel in 10). The magnetic energy has a smaller relative error (of 16%) than the maximumenergy reconstruction, while is a rather poor match. Although both the MHD and magnetofrictional runs eventually proceed to a lowenergy state, the closed box used in the MHD run prevents the system from relaxing as deep as in the magnetofrictional run. Furthermore, as argued above, we cannot expect the reconstruction of the final reference field to be accurate, since this field presents relatively strong forces close to the magnetogram due to the MHD eruption process.
5 Conclusions
This investigation demonstrates that a forcefree model of solar active regions, which includes relevant structural (i.e., topological or geometrical) features and is free of noise, can be reconstructed to a very high accuracy based only on its vector magnetogram. The reconstructed configuration consists of a flux rope with a nonneutralized current channel in its core. A bald patch separatrix surface or two quasiseparatrix layers combined into a hyperbolic flux tube form the interface between the flux rope and the ambient potential field. These elements are regarded as generic for many active regions, so that the model realistically captures their structural properties, in particular for essentially bipolar flux distributions.
Previous reconstructions of nonlinear forcefree test fields considered the simpler configurations of an approximately currentneutralized flux rope (Valori et al. 2005) and of a sheared arcade (Schrijver et al. 2006; Valori et al. 2007). A forcefree model of a filament channel containing a hollow current channel could also be reconstructed, albeit at significantly lower accuracy, which was likely caused by the combination with a noisy lineofsight magnetogram in the exterior of the channel (Metcalf et al. 2008). Thus, the reconstruction capability has now been demonstrated, mostly at excellent accuracy, for all basic types of forcefree equilibria currently considered in coronal physics, with the one treated here providing the highest degree of structural complexity. This confirms and substantiates previous conclusions (DeRosa et al. 2009; Metcalf et al. 2008) that insufficient accuracy of nonlinear forcefree extrapolation originates from the violation of the forcefree assumption by the input vector magnetogram, not from an insufficient intrinsic accuracy of the extrapolation codes.
The accuracy achieved in the reconstruction of stable modified TD equilibria clearly permits configurations containing an HFT, a BP, or none of these features to be differentiated. The energy, flux rope twist and position, and helicity are obtained at relative error levels of 1, 3, 4, and , respectively, even for a high twist near the threshold of the helical kink instability, slightly exceeding one full turn.
These quantities are found to depend only weakly on the initial field chosen, except for the apex height of the flux rope, which shows a displacement from the true value of up to (if a linear forcefree field is used). Starting the extrapolation from the conventional potential field, which assumes vanishing flux outside the magnetogram, yields the most accurate results for the considered equilibria. Other potential fields result in somewhat lower accuracy, and a linear forcefree field scores worst but still produces a rather acceptable reconstruction.
The high accuracy is reached using vector magnetograms that have a realistic lowflux margin of width 1/4 the size of the strongflux area and are exactly flux balanced. The effects that result from violating the latter condition require further study.
Attempting to extrapolate the magnetogram of an unstable modified TD equilibrium, the magnetofrictional relaxation produces an evolution in two phases of opposing tendency that clearly reflect the unstable nature of the equilibrium. A phase of rising free energy yielding a partial buildup of the flux rope is superseded by a decrease in the energy back to nearly the initial (potentialfield) value, destroying the flux rope.
AcknowledgementsWe thank G. Aulanier and the referee for constructive comments. This work was supported by the DFG, an STFC Rolling Grant, and by NASA grants NNH06AD58I and NNX08AG44G. The research leading to these results has received funding from the European Commission's Seventh Framework Programme (FP7/20072013) under the grant agreement No. 218816 (SOTERIA project, www.soteriaspace.eu). Funding from the European Comission through the SOLAIRE network (MTRMCT2006035484) is also gratefully acknowledged. The contribution of V. S. Titov was supported by NASA's Heliophysics Theory, Living With a Star, and SR&T programs, and the Center for Integrated Space Weather Modeling (an NSF Science and Technology Center).
References
 Alexiades, V., Amiez, G., & Gremaud, P. 1996, in IMACS'94, 2 [Google Scholar]
 Amari, T., Aly, J. J., Luciani, J. F., Boulmezaoud, T. Z., & Mikic, Z. 1997, Sol. Phys., 174, 129 [NASA ADS] [CrossRef] [Google Scholar]
 Amari, T., Boulmezaoud, T. Z., & Mikic, Z. 1999, A&A, 350, 1051 [NASA ADS] [Google Scholar]
 Amari, T., Boulmezaoud, T. Z., & Aly, J. J. 2006, A&A, 446, 691 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Antiochos, S. K. 1998, ApJ, 502, L181 [NASA ADS] [CrossRef] [Google Scholar]
 Bungey, T. N., Titov, V. S., & Priest, E. R. 1996, A&A, 308, 233 [NASA ADS] [Google Scholar]
 Canou, A., Amari, T., Bommier, V., et al. 2009, ApJ, 693, L27 [NASA ADS] [CrossRef] [Google Scholar]
 Chiu, Y. T., & Hilton, H. H. 1977, ApJ, 212, 873 [NASA ADS] [CrossRef] [Google Scholar]
 Chodura, R., & Schlüter, A. 1981, J. Comp. Phys., 41, 68 [NASA ADS] [CrossRef] [Google Scholar]
 Craig, I. J. D., & Sneyd, A. D. 1986, ApJ, 311, 451 [NASA ADS] [CrossRef] [Google Scholar]
 Demoulin, P., Henoux, J. C., Priest, E. R., & Mandrini, C. H. 1996, A&A, 308, 643 [NASA ADS] [Google Scholar]
 DeRosa, M. L., Schrijver, C. J., Barnes, G., et al. 2009, ApJ, 696, 1780 [Google Scholar]
 DeVore, C. R. 2000, ApJ, 539, 944 [NASA ADS] [CrossRef] [Google Scholar]
 Evje, S., Karlsen, K., Lie, K., & Risebro, N. 2000, in IMA Volumes in Mathematics and its Applications, Parallel Solution of Partial Differential Equations, ed. P. Bjorstad, & M. Luskin (Springer Science + Business Media, LLC), 120, 209 [Google Scholar]
 Fuhrmann, M., Seehafer, N., & Valori, G. 2007, A&A, 476, 349 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Galsgaard, K., Titov, V. S., & Neukirch, T. 2003, ApJ, 595, 506 [NASA ADS] [CrossRef] [Google Scholar]
 Gary, G. A. 2001, Sol. Phys., 203, 71 [NASA ADS] [CrossRef] [Google Scholar]
 Gibson, S. E., Fan, Y., Mandrini, C., Fisher, G., & Demoulin, P. 2004, ApJ, 617, 600 [NASA ADS] [CrossRef] [Google Scholar]
 Green, L. M., & Kliem, B. 2009, ApJ, 700, L83 [NASA ADS] [CrossRef] [Google Scholar]
 Guo, Y., Schmieder, B., Démoulin, P., et al. 2010, ApJ, 714, 343 [NASA ADS] [CrossRef] [Google Scholar]
 Isenberg, P. A., & Forbes, T. G. 2007, ApJ, 670, 1453 [NASA ADS] [CrossRef] [Google Scholar]
 Keppens, R., Nool, M., Tóth, G., & Goedbloed, J. 2003, Comp. Phys. Comm., 153, 317 [Google Scholar]
 Kliem, B. & Török, T. 2006, Phys. Rev. Lett., 96, 255002 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Kliem, B., Titov, V. S., & Török, T. 2004, A&A, 413, L23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Leka, K. D., Canfield, R. C., McClymont, A. N., & van DrielGesztelyi, L. 1996, ApJ, 462, 547 [NASA ADS] [CrossRef] [Google Scholar]
 Leka, K. D., Fan, Y., & Barnes, G. 2005, ApJ, 626, 1091 [NASA ADS] [CrossRef] [Google Scholar]
 Lites, B. W., Low, B. C., Martinez Pillet, V., et al. 1995, ApJ, 446, 877 [NASA ADS] [CrossRef] [Google Scholar]
 Low, B. C., & Lou, Y. Q. 1990, ApJ, 352, 343 [NASA ADS] [CrossRef] [Google Scholar]
 McClymont, A. N., Jiao, L., & Mikic, Z. 1997, Sol. Phys., 174, 191 [NASA ADS] [CrossRef] [Google Scholar]
 McKenzie, D. E., & Canfield, R. C. 2008, A&A, 481, L65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Metcalf, T. R., Jiao, L., McClymont, A. N., Canfield, R. C., & Uitenbroek, H. 1995, ApJ, 439, 474 [NASA ADS] [CrossRef] [Google Scholar]
 Metcalf, T. R., Derosa, M. L., Schrijver, C. J., et al. 2008, Sol. Phys., 247, 269 [Google Scholar]
 Moon, Y., Choe, G. S., Yun, H. S., Park, Y. D., & Mickey, D. L. 2002, ApJ, 568, 422 [NASA ADS] [CrossRef] [Google Scholar]
 Nakagawa, Y., & Raadu, M. A. 1972, Sol. Phys., 25, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Pevtsov, A. A., Canfield, R. C., & Metcalf, T. R. 1995, ApJ, 440, L109 [NASA ADS] [CrossRef] [Google Scholar]
 Priest, E. R., & Démoulin, P. 1995, J. Geophys. Res., 100, 23443 [NASA ADS] [CrossRef] [Google Scholar]
 Régnier, S., & Amari, T. 2004, A&A, 425, 345 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Régnier, S., Amari, T., & Kersalé, E. 2002, A&A, 392, 1119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roumeliotis, G. 1996, ApJ, 473, 1095 [NASA ADS] [CrossRef] [Google Scholar]
 Roussev, I. I., Forbes, T. G., Gombosi, T. I., et al. 2003, ApJ, 588, L45 [NASA ADS] [CrossRef] [Google Scholar]
 Sakurai, T. 1981, Sol. Phys., 69, 343 [NASA ADS] [CrossRef] [Google Scholar]
 Sakurai, T. 1989, Space Sci. Rev., 51, 11 [NASA ADS] [Google Scholar]
 Schmidt, H. U. 1964, NASA Spec. Publ., 50, 107 [NASA ADS] [Google Scholar]
 Schrijver, C. J., DeRosa, M. L., Metcalf, T. R., et al. 2006, Sol. Phys., 235, 161 [Google Scholar]
 Schrijver, C. J., DeRosa, M. L., Metcalf, T., et al. 2008a, ApJ, 675, 1637 [Google Scholar]
 Schrijver, C. J., Elmore, C., Kliem, B., Török, T., & Title, A. M. 2008b, ApJ, 674, 586 [NASA ADS] [CrossRef] [Google Scholar]
 Seehafer, N. 1978, Sol. Phys., 58, 215 [NASA ADS] [CrossRef] [Google Scholar]
 Seehafer, N. 1986, Sol. Phys., 105, 223 [NASA ADS] [CrossRef] [Google Scholar]
 Shafranov, V. D. 1966, Rev. Plasma Phys., 2, 103 [NASA ADS] [Google Scholar]
 Titov, V. S. 2007, in Reconnection of magnetic fields: magnetohydrodynamics and collisionless theory and observations, ed. J. Birn, & E. R. Priest (Cambridge, UK: Cambridge University Press), 250 [Google Scholar]
 Titov, V. S., & Démoulin, P. 1999, A&A, 351, 707 [NASA ADS] [Google Scholar]
 Titov, V. S., Priest, E. R., & Demoulin, P. 1993, A&A, 276, 564 [NASA ADS] [Google Scholar]
 Titov, V. S., Démoulin, P., & Hornig, G. 1999, in Magnetic Fields and Solar Processes, ed. A. Wilson, et al., ESA SP, 448, 715 [Google Scholar]
 Titov, V. S., Hornig, G., & Démoulin, P. 2002, J. Geophys. Res. (Space Physics), 107, 1164 [Google Scholar]
 Titov, V. S., Galsgaard, K., & Neukirch, T. 2003, ApJ, 582, 1172 [NASA ADS] [CrossRef] [Google Scholar]
 Török, T., & Kliem, B. 2003, A&A, 406, 1043 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Török, T., & Kliem, B. 2005, ApJ, 630, L97 [Google Scholar]
 Török, T., & Kliem, B. 2007, Astron. Nachr., 328, 743 [NASA ADS] [CrossRef] [Google Scholar]
 Török, T., Kliem, B., & Titov, V. S. 2004, A&A, 413, L27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valori, G., Kliem, B., & Keppens, R. 2005, A&A, 433, 335 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valori, G., Kliem, B., & Fuhrmann, M. 2007, Sol. Phys., 245, 263 [NASA ADS] [CrossRef] [Google Scholar]
 van Ballegooijen, A. A. 2004, ApJ, 612, 519 [NASA ADS] [CrossRef] [Google Scholar]
 van Ballegooijen, A. A., DeLuca, E. E., Squires, K., & Mackay, D. H. 2007, JASTP, 69, 24 [NASA ADS] [Google Scholar]
 Wheatland, M. S. 2000, ApJ, 532, 616 [NASA ADS] [CrossRef] [Google Scholar]
 Wheatland, M. S., Sturrock, P. A., & Roumeliotis, G. 2000, ApJ, 540, 1150 [Google Scholar]
 Wiegelmann, T., Inhester, B., Kliem, B., Valori, G., & Neukirch, T. 2006a, A&A, 453, 737 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wiegelmann, T., Inhester, B., & Sakurai, T. 2006b, Sol. Phys., 233, 215 [NASA ADS] [CrossRef] [Google Scholar]
 Wiegelmann, T., Thalmann, J. K., Schrijver, C. J., Derosa, M. L., & Metcalf, T. R. 2008, Sol. Phys., 247, 249 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, W. H., Sturrock, P. A., & Antiochos, S. K. 1986, ApJ, 309, 383 [NASA ADS] [CrossRef] [Google Scholar]
 Zarro, D. M., Canfield, R. C., Nitta, N., et al. 2000, BAAS, 32, 817 [NASA ADS] [Google Scholar]
Footnotes
 ... bipolar^{}
 For a quadrupolar active region in the classical sense one would require the neutral line to meander substantially between opposite flux concentrations or even to be split, with one or several parts forming a closed loop.
 ... points^{}
 The derivatives involved in and are computed to fourth order accuracy, except for the reference fields and for the potential and linear fields, where a second order differencing is used, and onesided derivatives of corresponding order are employed at, and where necessary also just above, the plane .
All Tables
Table 1: Normalized parameters defining the TD equilibria.
Table 2: Normalized flux, force, and torque imbalance of the five analyzed magnetograms.
Table 3: Reconstruction accuracy of the stable TD equilibria.
Table 4: Reconstruction accuracy of the Low_HFT equilibrium for different initial conditions.
Table 5: Figures of merit (absolute values) for all reference, initial, and extrapolated fields.
All Figures
Figure 1: ``Lineofsight'' magnetogram B_{z}(x,y,0) with superimposed arrows representing the transverse field ( top), and forcefree parameter ( bottom), for the reference fields of the cases High_HFT ( left) and BP ( right). 

Open with DEXTER  
In the text 
Figure 2: Contours of , with superimposed arrows representing the transverse field in the vicinity of the neutral line of B_{z}, for the reference ( top) and reconstructed ( bottom) BP cases. Note the inverse crossing of the neutral line, the characteristic signature of a BP, in the middle of the plots. 

Open with DEXTER  
In the text 
Figure 3: Selected field lines for the Low_HFT case; a) reference, b) extrapolated, c) potential, and d) linear field (with ). In a) and b) the flux rope is depicted by two groups of field lines (in redorange and green), each one starting from six uniformly spaced points on a circle of radius 0.1, centered on the intersections of the magnetic axis of the flux rope (CFL) with the z=0 plane. Gray field lines are started from the zaxis and visualize the potential field above and below the current channel. In c) and d) the field lines are started from the same set of footpoints as in a) and b). Contour plots of B_{z}(x,y,0) and its neutral line B_{z}(x,y,0)=0 (dotdashes) are shown in the bottom plane of each panel. 

Open with DEXTER  
In the text 
Figure 4: Plane averages of the relative solenoidal error in the No_HFT case, as a function of position. 

Open with DEXTER  
In the text 
Figure 5: Isosurface of the reconstructed current channel at 31% of peak current density for the cases High_HFT ( top) and No_HFT ( bottom). Field lines crossing the z axis are shown in gray if below the HFT and in green if above the HFT. 

Open with DEXTER  
In the text 
Figure 6: Poloidal field at the z axis, B_{x}(0,0,z), in the No_HFT case, for the reference (red), extrapolated (green), and TD_ex potential fields (blue). Left: entire z axis extension; right: upper part (above the current channel) in logarithmic scale. 

Open with DEXTER  
In the text 
Figure 7: Poloidal field at the z axis, B_{x}(0,0,z), in the Low_HFT case, for the reference (red solid line), initial (blue), and corresponding extrapolated fields (green), for different type of the initial field: potential Schmidt (short dashes), potential Seehafer (dotdash), linear alphabest Seehafer (triple dotsdash), and potential TD_ex (long dashes). 

Open with DEXTER  
In the text 
Figure 8: Selected field lines for the reference fields of the unstable case at the initial ( top) and final ( bottom) stages of the MHD evolution (see Fig. 3 for the color coding of the field lines). 

Open with DEXTER  
In the text 
Figure 9: Evolution of (solid line) and (dashes) during the extrapolation runs of the cases BP (left) and Unstable (right). The energy is normalized to its value at the first iteration. 

Open with DEXTER  
In the text 
Figure 10: Selected field lines for the extrapolation of the unstable case at different iterations. From left to right: the initial potential field at it = 1, the maximumenergy configuration at 6700 iterations, and the final lowenergy configuration at iterations. See Fig. 3 for the color coding of the field lines. 

Open with DEXTER  
In the text 
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.