Global Non-Potential Magnetic Models of the Solar Corona During the March 2015 Eclipse

Seven different models are applied to the same problem of simulating the Sun's coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models.


Introduction
In recent years, a number of different approaches have been developed for modelling non-potential magnetic fields in the Sun's atmosphere, based on measurements of the magnetic field on the solar surface.Non-potential means that electric currents are allowed to be present within the coronal volume, in contrast to traditional potential-field extrapolations.However, different modellers use not only different approximations for the coronal magnetic field but also different input data and boundary conditions.To date, a direct comparison of these various approaches has been lacking in the literature.With this in mind, we convened a scientific team at the International Space Science Institute in Bern, Switzerland, and we report here on the results.
In this paper, we present a number of non-potential models side-by-side in a way that enables a direct comparison.Since the problem of coronal modelling is  (Druckmüller 2009).This is overlaid with a combination of 171, 193 and 211 Å (red, green, blue respectively) channels from the AIA instrument on the Solar Dynamics Observatory satellite.Each of the channels were processed using Multiscale Gaussian Normalization (MGN, Morgan and Druckmüller 2014) prior to combination.
very much an active area of research, there are inevitable disagreements between the models.By highlighting the similarities and differences between different modelling approaches, and why they occur, we hope to assist the community in moving toward a single non-potential modelling approach for the corona.We would like to highlight what factors should be taken into account, and what are the biggest uncertainties in existing models.
We have fixed a single date and time, coinciding with the total solar eclipse of 2015 March 20.Fixing a single time allows us to compare both static and time-dependent models, and we chose an eclipse date during the Solar Dynamics Observatory (SDO) era so as to maximize the available observations of the real coronal structure.Figure 1 shows the structure of the observed corona on that day.There are no large active regions on the visible solar disk, although there is some activity on the far side and visible above both East and West limbs.The streamer structure is relatively complex, consistent with the fact that this eclipse occurred shortly after Solar Maximum.The streamers also show an asymmetry between the north and south poles, consistent with the formation of a polar coronal hole in the south but not in the north (Petrie 2017).In terms of non-potential structure, there is a clear polar crown prominence on the north-east limb, with suggestions of a coronal cavity surrounding it.Although this is not a time of particularly high solar activity, the coronal structure is nevertheless more complex than that typically found around Solar Minimum.It is therefore quite challenging to model all of the aspects of this observed corona within the context of a single global magnetic field model.
The compared models are summarized in Section 2. All cover the full global solar corona, except for the polar regions in some cases.The outer boundaries of the models vary; for the comparisons in this paper we will therefore focus on the region R ⊙ ≤ r ≤ 2.5R ⊙ that is covered by all of the models.Grid resolution and boundary conditions were chosen by each individual modeller.In particular, different photospheric boundary conditions were used for each model, and these differences in input data are highlighted in Section 2. All models used magnetic data from the same SDO/HMI instrument, but different methods were used to reconstruct boundary maps of the full solar surface, leading to quite different photospheric magnetic fields.Owing to these differences in input data, our study should not be viewed as a direct side-by-side comparison of numerical codes for modelling the coronal magnetic field, such as (for example) the active region exercises by Schrijver et al. (2008) and DeRosa et al. (2015).Rather, the approach is more qualitative: to compare and contrast the results from different non-potential modelling approaches, where these approaches include the different methods of processing magnetogram input data that are characteristic of present-day coronal modelling.Through this we can assess the strengths and weaknesses of each approach.
Because we only observe (usually) one side of the Sun, it is difficult to make detailed comparisons with observations using only a single snapshot such as this.Moreover, the lack of direct magnetic measurements above the photosphere is a major driver for the development of these models in the first place.Nevertheless, we are able to compare with indirect observations of the coronal magnetic field, and do so qualitatively in this paper.It is hoped that by combining all available information, we can ultimately build a "best guess" picture of the coronal magnetic structure during the eclipse.Following our summary of the models in Section 2, the resulting coronal magnetic fields are compared in Section 3, both with each other and with the indirect observations.The overall findings are discussed in Section 4.

Non-Potential Models
The models presented in this paper are summarized in Table 1, and their computational domains are summarized in Table 2.The latter table also shows magnetic energies, which will be discussed in Section 3.For purposes of illustration we have selected a single representative model of each type, although additional runs were carried out in many cases.These models are still under development, and not yet freely available, but readers interested in their use are encouraged to contact the individual modellers, as listed in Table 1.
Space precludes a detailed description of each model, so the following subsections concentrate on their similarities and differences.For a review of the various non-potential modelling approaches, see Mackay and Yeates (2012) or Wiegelmann, Petrie, and Riley (2017), or follow the references given in Table 1.Some additional details are given  for the mhd-mas and mf models in Appendices A and B respectively, since these differ somewhat in their setup from previously published simulations.
We stress that, although we consider only a single model of each type, it would be possible to "mix and match" aspects of different models in future.For example, the techniques used to energize filament channels in the mhd-mas model could equally be applied to any of the other static models.Or thermodynamics effects could be included in the mhd-mas model.Thus our comparison should be treated only as a representative selection, aiming to illustrate the consequences of including different model features.

Modelling approaches
In the low corona, on a global scale, magnetic forces dominate over gravity and pressure gradients.Thus the large-scale magnetic field configuration is expected to be approximately force-free at the heights r < 2.5 R ⊙ of interest here (Wiegelmann and Sakurai 2012).This means that the Lorentz force vanishes and the magnetic field B satisfies where j = c∇ × B/(4π) is the electric current density.In a non-potential magnetic field with free energy, we must have j = 0 at least in some regions of the corona.All of the models in Table 1 aim to approximate force-free equilibria of the form (1), except for mhs and mhd-cese where plasma pressure and gravity also influence the equilibria.Namely, the mhs model uses a special class of magnetohydrostatic equilibria where the electric current density flows on spherical shells.Here we set the parameter a = 1.0 (for more details and the definition of a, see Bogdan and Low 1986).In this model, the radial magnetic field B r on r = R ⊙ is used as a boundary condition for spherical harmonic decomposition.The resulting magnetic field contains a non-vanishing Lorentz force, which is compensated by pressure gradients and gravity.The density and plasma pressure distribution are computed as a superposition of these terms and a stratified solar atmosphere model.The mhd-cese model solves the full magnetohydrodynamic (MHD) equations including gravity, centrifugal and Coriolis terms, as well as a volumetric heating term, relaxing to a static equilibrium that extends out to r = 30 R ⊙ , albeit at a lower resolution than the other models.All of the other models have an outer boundary at or near r = 2.5 R ⊙ , except for ffe where the computational domain extends to 5.8 R ⊙ , although a perfectly-matched (absorbing) layer is implemented above 2.5 R ⊙ .
Whilst the nlf-op, nlf-gr , ffe, mf and mhd-mas models all use some form of iteration toward j × B = 0 -constrained by the imposed photospheric boundary conditions -they differ significantly in how the electric current structures in the final magnetic field are built up.In the first three models, the electric current distribution is fixed by the photospheric boundary conditions.By contrast, in the mf model the coronal magnetic field is driven by large-scale horizontal flows in the photosphere, at the same time as relaxing the coronal magnetic field toward j × B = 0.This leads to an evolving photospheric B r distribution.The advantage of this approach is that it can build up large-scale current structures in a manner mimicking that in the real Sun.But the disadvantage is that a whole time series of input data extending over several months is required on the photospheric surface, both for B r and the horizontal flows.Moreover, the final B r distribution in the mf model matches that of the original magnetogram only on a coarse scale.The mhdmas model is in some sense a hybrid: although it uses a single observed map of B r on the solar surface, additional electric fields are imposed to drive the formation of filament channels at specified locations in the corona.The locations of these channels, and in particular their chirality (direction of the transverse magnetic field) were guided by the results of the mf model and by EUV observations (see Appendix A).

Photospheric boundary conditions
The different modelling techniques each require particular input data at the solar photosphere.All of the models except mf use only a single magnetic map covering the full photospheric surface, notionally at the time of the eclipse.Vector magnetic data with all three components of B are required for nlf-op, nlf-gr and ffe, whereas only the radial component B r is utilized by all of the other models.
Since the observed vector magnetograms are taken at photospheric heights where the magnetic field is not generally force free, these are generally not consistent with (1).The three models that use vector magnetic input deal with this in different ways.The nlf-gr model uses the Grad-Rubin scheme.By imposing the force-free function α = j r /B r in only one polarity (here, B r > 0), as well as B r everywhere, it matches the transverse components of B in locations of this polarity.Some cut-offs on the B r and on the transverse components of B are applied so as to generate regular values of α, and in particular it is assumed that j r = 0 in regions of weak B r .The nlf-op method uses an optimization that penalizes differences in the simulated and observed vector B over the whole map, but allows differences, particularly in pixels with smaller observed transverse field.The ffe method imposes the vector B over the whole map, although this leads to a boundary layer just above the photosphere where B becomes more force free.
The mhs, mf , mhd-cese and mhd-mas models require only the observed radial magnetic field B r on the photospheric boundary rather than the full vector B, although as mentioned above the mf model requires a whole time sequence of B r measurements.In the Bogdan-Low solution used by the mhs model, the transverse components of B on the boundary are uniquely determined by the form of the solution.However, in the methods used by mf , mhd-cese and mhd-mas, a transverse electric field is required to specify a unique evolution from the initial potential field.
This transverse electric field is set either by boundary flows or through adhoc relations.In mhd-cese, the transverse electric field is simply set to zero as the corona is evolved toward equilibrium (the initial plasma parameters are taken from Parker's hydrodynamic isothermal solar wind solution).In mhd-mas, the transverse electric field is set to zero except where it is imposed to create magnetic flux ropes in filament channels (see Appendix A).The mf model is rather different since transverse electric fields are imposed continually throughout the 200-day evolution, so as to evolve the photospheric B r distribution as well as the coronal magnetic field.This is described in more detail in Appendix B.

Magnetogram input data
All models used photospheric magnetogram data from the HMI instrument on Solar Dynamics Observatory (Schou et al. 2012).However, it is important to note that these data were processed in different ways by each modeller, leading to different distributions of B r on r = R ⊙ .Rather than standardize this processing and do a direct model comparison, we have decided to highlight the fact that the input data themselves represent a significant difference between present models of the coronal magnetic field.The differences are summarized in this section.
Fig. 2 The Br, B θ , and B φ components of the custom SDO/HMI vector map, with color tables saturated at ±50 G (red positive, blue negative).The horizontal axis shows Carrington longitude and the vertical dashed line marks the location of central meridian at the time of eclipse (99 • ).The vertical solid line marks the time discontinuity between data from the synoptic maps for CR2161 (to the left) and CR2162 (to the right).This is at 279 • rather than the typical 360 • .
Several of the models used a custom SDO/HMI vector map that was created specially for this project, shown in Figure 2.Here data from the successive Carrington rotations CR2161 and CR2162 were combined so as to extract a synoptic map centered on the eclipse Central Meridian (Carrington longitude ≈ 99 • during CR2161), with the discontinuity located 180 • out of phase from this (at Carrington longitude 279 • ).The data were processed using the standard HMI pipeline (Hoeksema et al. 2014), but with azimuth-angle disambiguated full-disk images.Strong fields (|B| ≥ 50 G) were left unsmoothed, but boxcar smoothing was applied to low-latitude quiet-Sun fields (|B| < 50 G), in order to improve the signal-tonoise ratio.High latitude fields (poleward of ±70 • latitude) were filled using a combination of smoothed vector B r and radially-corrected line-of-sight field.
For the ffe model, this custom vector map (with dimensions 1440 × 3600) was simply boxcar binned to a lower resolution.To generate the initial potential fields for the nlf-op and mhs models, the map was initially smoothed by retaining only spherical harmonics of degree ℓ ≤ 25.The full map was then used in the ensuing iteration for nlf-op, while the mhs input map remains smoothed.For nlf-op, note that no pre-processing of the vector data (as described by Tadesse, Wiegelmann, and Inhester 2009) was applied.
The nlf-gr model took the basic map as a starting point but made two modifications.Firstly, B r on the visible disk (centred at Carrington longitude 99 • ) was replaced by an HMI longitudinal 720-second full-disk magnetogram from 09:47:58UT on 2015 March 20.Secondly, six active regions on the visible disk were replaced by HMI vector Space-Weather HMI Active Region Patches (HARPs; Bobra et al. 2014).Specifically, patches number 5337, 5342, 5345, 5347, 5348, and 5350, again observed at 09:47:58UT.These additions were made on a higher resolution grid (3804 × 5784), which was subsequently interpolated to a final mesh of 250 × 500 for the computation.This final mesh is non-uniform with more grid points within active regions.
The mhd-cese model used a standard SDO/HMI synoptic magnetogram for CR2161, smoothed with a simple Gaussian kernel (as a function of heliocentric angle).The mhd-mas model used a B r map constructed from HMI data in a similar way to the custom map, combining synoptic maps from CR2161 (from 0 • to 279 • ) and CR2162 (from 279 • to 360 • ).The high-latitude fields were filled in from previous Carrington rotations when the poles were visible, and the maps were suitably smoothed for the MHD calculation.The mf model did not use HMI data directly, being driven instead by a flux transport model coupled with the emergence of bipolar magnetic regions starting from 2014 September 1 -this is described in Appendix B.
It is important to note that the models using synoptic input data from both CR2161 and CR2162 do not represent what could have been achieved in an advance prediction of the eclipse, since the maps incorporate magnetogram data taken after the eclipse time.The exceptions to this are mf and mhd-cese, which use only information from the visible face of the Sun up to but not after the eclipse.
Figure 3 shows the distribution of radial magnetic field at the photosphere r = R ⊙ in each model (at the closest simulated time to the eclipse for the mf model).It is clear that there are differences between these input data for the different models.Most obvious are the variations in resolution and in the amount of smoothing applied to the data, as described above.But important also are the timing variations, namely how recently different regions of the maps have been updated.There are two notable differences that affect the model results presented below.Firstly, a small active region (NOAA 12304) is located at approximately 140 • Carrington longitude, in the Northern hemisphere.This region is visible in the EUV eclipse image (Figure 1), in the Northern hemisphere between disk center and the West limb.However, it is missing from the HMI vector synoptic map in Figure 2 because it emerged after that Carrington longitude had rotated past Central Meridian.Accordingly, it is omitted from models nlf-op, mhs, ffe, mhd-cese and mhd-mas.It is included, however, in nlf-gr thanks to the replacement of the visible disk by an HMI full-disk magnetogram taken just prior to the eclipse.The region is also included in the mf model owing to the method of bipolar magnetic region (BMR) insertion (Appendix B).The second difference is potentially more significant: all models except for mf and mhd-cese include the activity complex centred around 10 • Carrington longitude in the Southern hemisphere, and appearing at the East limb in Figure 1.In the mf model the corresponding BMRs have not yet been assimilated at the time of eclipse, since they were not yet fully visible at that time and their properties could not be determined.Similarly they do not appear in mhd-cese since the corresponding longitudes were taken from earlier in CR2161, rather than CR2162 as in the other models.Such differences are typical of coronal magnetic field modelling based on synoptic magnetic field observations, but must be taken into account when assessing our results, which we consider next.We reiterate that the inclusion of this activity complex in the other models was possible only by using observations taken after the eclipse; if the comparison had been carried out in real time then this correction would not have been possible.

Model Comparison
A number of aspects of the models are compared in the following subsections, with the aim of highlighting the similarities and differences between them.Since many of the models do not include realistic thermodynamics (or indeed any plasma at all), we limit our comparison to the magnetic field, rather than plasma density or temperature.Some of the models, such as the mhd-mas model, have already been developed extensively to describe the flow of energy in the corona, as well as the acceleration of the solar wind, and can deduce the coronal temperature and density.In fact, we also studied this eclipse using this more sophisticated version of the mhd-mas model, but we do not describe those results here since they exceed the scope of the comparisons presented in this paper.Ground truth observations of the real corona are included, although these can provide only indirect, qualitative information about the coronal magnetic field.

Magnetic Flux and Energy
We begin with some overall diagnostics.Firstly, Figure 4(a) shows the total unsigned magnetic flux, for each model at a sequence of heights in the corona.At r = R ⊙ there are significant differences in the photospheric flux arising from the difference in resolution between the models (these different resolutions are stated in Table 2 and evident in Figure 3).The model photospheric fluxes span the range 3 − 5.5 × 10 23 Mx, from mhd-cese with the lowest flux to ffe with the highest.For comparison, note that the SDO/HMI synoptic map in Figure 2 and plotted against the corresponding energy E p for a potential field source surface extrapolation with source surface at r = 2.5R ⊙ and the same B r (R ⊙ , θ, φ) for each particular model.These energies are also listed in Table 2.There are considerable variations in both E and E p , which depend to some extent on the model resolution.
The distance of each symbol above the solid line indicates the excess of nonpotential energy E above E p , or the relative free energy.This varies from very little free energy for nlf-gr to around 40% for nlf-op, mhs, mf and mhd-mas.Whilst the latter models all produce a similar percentage of free energy, the absolute values differ significantly, and we will see in Section 3.2 that the distributions of electric current within the volume also differ significantly.In the case of the mhs model, the free energy depends on the choice of the a parameter; for example, taking a = 2 would lead to E = 1.73 E p rather than 1.42 E p for a = 1.The particularly low free energy for nlf-gr arises because the electric currents are strongly concentrated within active regions, with the magnetic field being close to potential throughout most of the volume.

Electric Currents
Figures 5(a) and (b) show measures of the average perpendicular and parallel electric currents in each model, as a function of radius.Specifically we plot the quantities and The perpendicular measure F (r) is essentially a measure of the Lorentz force j×B.
Figure 5(c) shows the ratio J (r)/F (r), which is a measure of "force-freeness" for each model.Below about r = 1.5 R ⊙ , the ratio J/F divides the models into three broad classes: those models that are relatively force-free with J/F ≫ 1 (mhd-mas, nlfgr , nlf-op and mf ), those that are not force-free, with J/F < 1 (mhd-cese and mhs), and the ffe model which is not force-free at the photosphere but becomes rather more so above about 1.1 R ⊙ .These differences arise from the physics of the models: the mhd-cese and mhs models include non-magnetic terms in their force balances, while the ffe model matches to a boundary condition at r = R ⊙ that does not satisfy j × B = 0.The model with highest ratio J/F is mhd-mas, owing to the relatively smooth boundary data and level of numerical relaxation applied.
The actual amount of current near the photosphere varies between models, as seen by J (r) in Figure 6(b).This is highest for nlf-gr , because it includes the most fine-scale structure in the HARP patches.On the other hand, J (r) falls off rapidly with height in this model, consistent with the field being closest to potential overall  (Table 2 and Figure 4b).By contrast, the mhs model has relatively low J (r) near the photosphere, owing to its lower input resolution, and the lowest J (r) is found for the mhd-cese model, which lacks a mechanism for energizing the magnetic field there; its free magnetic energy is located at larger radii where the non-magnetic terms become important.
Both the mf and mhd-mas models have lower J (r) very close to the photosphere, because they do not resolve such fine-scale structures as nlf-gr and nlf-op.However, these two models have larger current in the low corona (1.1 to 1.2 R ⊙ ), arising from the incorporation of low-lying filament channels.These form selfconsistently in the mf model, while their locations (for imposed currents) were chosen by design in the mhd-mas model.Between about 1.2 and 1.5 R ⊙ , the mf model has the largest J (r), owing to the significant volume currents that have been ejected over time by flux emergence and surface motions.As in the nlf-op model, which has almost as high J (r), these currents are not limited to active regions, and reach greater heights in the corona.
Above about 1.7 R ⊙ , the ffe model has the largest J (r), although this likely results from incomplete relaxation to equilibrium in the model and deserves further investigation.The mhd-cese model also has a more gradual fall-off of J (r) with height than many of the other models, again reflecting the influence of non forcefree effects at larger heights in the model as the magnetic field is opened out by the solar wind.Similarly, the mf model becomes non force-free (J/F < 1) above about 2 R ⊙ , caused by the imposed outflow at the upper boundary.This outflow is used to simulate the effect of the solar wind (Mackay and van Ballegooijen 2006), so that the magnetic field is no longer force-free above 2 R ⊙ but is in a steady-state balance between the outflow and the Lorentz force.
To illustrate the spatial distributions of currents, Figure 6 shows the magnitude of vertical current density j r (1.02 R ⊙ , θ, φ) in each model, with the same logarithmic color scale used for all plots.For comparison, panel (h) shows the magnitude of j r (R ⊙ , θ, φ) computed directly from the SDO/HMI vector synoptic map (Figure 2).
In all models, the strongest currents are located in the active regions, as might be expected since the magnetic field is strongest there.The magnitudes of these currents differ significantly between the models, consistent with the variations seen in Figure 5. Outside the active regions, the models differ in their distributions of current.The nlf-op, nlf-gr , mhs and ffe models all show current distributions that correlate with the locations of strongest |j r | in the photospheric input map (Figure 6h).These locations are essentially those with strongest large-scale B r in the input map (Figure 2).Compared to the other models, nlf-op shows a higher background level of current in the quiet Sun, consistent with its values of J (r) and F (r) at this height in Figure 5.By contrast, nlf-gr shows a much lower level of background current outside of active regions, due to the assumption that j r = 0 in those footpoints.The mhd-cese model has small electric currents at low heights, as previously discussed.
The mf and mhd-mas models have additional current concentrations outside the locations of strong observed photospheric |j r |.These take the form of concentrated filament channels, seen in Figure 6 as parallel lines of |j r |, and lying above polarity inversion lines in the photospheric B r .A good example lies between about 60 • and 125 • Carrington longitude in the Southern hemisphere in both models.This current concentration is absent from the other models, and is not seen in the observations at the photospheric level.It is a concentration of electric current density in the coronal volume.We will return to this below in Section 3.4.

Plane-of-Sky Magnetic Structure
Figures 7(a) to (g) show visualizations of each model with field lines selected to show the plane-of-sky magnetic structure, as viewed from Earth at approximately the eclipse time.We must bear in mind that, with magnetogram observations presently available only from the Earth's viewpoint, all of the models use primarily synoptic observations built up from central meridian data.These do not include cotemporal information near the limbs, so any comparison can only be approximate.Recall also that the nlf-op and mhs models do not include the region poleward of ±70 • latitude.The model images can be compared to observations of the real corona in EUV (Figure 7h) and the Fe XIV line (Figure 7i), as well as in white light out to a larger radius (Figure 1).These particular observations have been chosen and processed to bring out as clearly as possible the structure of coronal streamers above the limb.References describing the processing are given in the figure caption.When comparing with the models it must be remembered that the observations see total emission along the line-of-sight, which includes structures in front of or behind the sky plane.It is beyond the scope of this project to do forward modelling of coronal emission, particularly since most of the models are purely magnetic.Nevertheless, it is clear that there are significant differences between the streamer structure in the different models, as well as similarities.
Firstly, we observe that streamers (closed field regions) are often too high in many models, particularly mhd-cese, mhs, ffe and nlf-gr .The white light observations (Figure 1) suggest their cusps to lie below 2 R ⊙ in most cases, although there could be larger closed loops that cannot be seen due to signal-to-noise issues.In the nlf-gr model, these closed-field regions are close to potential, so the cusp height is set by the potential field source surface at 2.5 R ⊙ .In the mhd-mas model, the outer boundary condition is one of zero velocity, tending also to keep the streamer cusps at the outer boundary.In ffe and mhd-cese there is no source surface at all, allowing the closed regions to extend even further out.Similarly,  2014) (with parameter h = 0.9) to bring out the streamer structure.Several angular regions are labelled in panel (h) for ease of reference in the text.
the mhs model is an infinite-space solution with no outer boundary imposed.In the nlf-op and mf models the electric currents allow the streamers to have lower cusps.This is further helped in the mf model by the radial outflow which pulls out the field lines.These differences highlight the importance of the outer boundary conditions for this kind of modelling.Some models show unusual field line behavior near the outer boundary.For example, mf has disconnected U-loops, which are the result of the ejection of a magnetic flux rope (Mackay and van Ballegooijen 2006;Yeates 2014).Also, some of the field lines for nlf-op and ffe are less smooth near the outer boundary.In the case of ffe this is consistent with the large currents seen at larger heights, and the fact that the relaxation has not reached a stable equilibrium.Although these field lines are strongly non-potential, they contribute little to E because the magnetic field strength is weak.
We can also compare the angular positions of particular streamers.These are sensitive to both the input data and the locations of electric currents in the corona.For ease of reference, eight approximate locations are labelled in Figure 7(h), although these do not necessarily correspond to individual streamers.Looking first at the West limb (A to D), the most prominent streamers in the observations are at B and C. All models show evidence of closed field near these locations, although their morphology differs significantly between models.Except for mhd-cese, the models also tend to agree that there is a narrow pseudo-streamer structure at D. At location A, the observed streamer is less clear in EUV, although there are indications of closed field in Figures 1 and 7(i).Many of the models show at least some closed field at this location, although its structure and orientation are quite varied between models.As mentioned above, the mf model shows U-loops here resulting from a recent flux rope ejection.Interestingly, there are indications in SWAP of a high-altitude EUV cavity at this location, although it is clearest on the day before that shown in Figure 7(h).Such a cavity is often associated with an underlying filament or filament channel.
On the East limb, the EUV observations indicate prominent structures at E, F and G around the equator, and there is a filament channel at H on the polar crown.All models indicate the presence of closed field around E, F and G, associated with the several active regions spread around Carrington longitude 0 • .The mf model has rather less closed field here than the other models, owing to the fact that two of the main bipolar regions have not yet reached the visible face of the Sun and been incorporated into the time-evolving mf simulation.This is an example of when additional magnetogram observations around the East limb, such as might be obtained from a satellite at the L5 Lagrange point, would be beneficial if the simulation were to be carried out in real time (Mackay, Yeates, and Bocquet 2016).
At H, all models agree that there is a closed field arcade overlying the East-West polarity inversion line.This is supported by the white light and EUV observations which show a filament at this location.Sheared magnetic field at low heights corresponding to this filament channel has formed naturally in the mf model, and has correspondingly been added to the mhd-mas model.This sheared field was created in mf by differential rotation and flux cancellation at the polarity inversion line, and is not present in the other static models.In these other models, this arcade is closer to potential.

Filament Channels
An important observable signature of non-potential magnetic structure in the corona is the presence of filament channels and filaments.These are located above polarity inversion lines in the photospheric B r , and are understood to have a On the solar disk there are a number of filaments visible in the observations, both in Hα (Figure 8h) and in EUV (Figure 8i).One filament extends around the North-East limb as a prominence in the white-light eclipse image (Figure 1).This is the structure labelled H in Figure 7(h).The presence of filaments at these particular locations indicates that sheared magnetic field is present.Only the mf and mhd-mas models have significantly sheared magnetic field along polarity inversion lines.In the case of the mf model, this has built up naturally over time due to surface motions and flux cancellation, whereas in the mhd-mas model it has been imposed based on the mf model results (Appendix A).
In principle, it would be possible for the nlf-op and nlf-gr models to contain sheared fields along these polarity inversion lines.Indeed the nlf-gr model does recover a sigmoidal structure within the active region located around Carrington longitude 190 • (on the far side of the Sun during the eclipse).However, the transverse magnetic field measurements, from which the currents are determined, have a poor signal-to-noise ratio outside of active regions.This arises because the linear polarization scales with |B| 2 , in contrast to the circular polarization which is linear in |B|.Moreover, by creating synoptic maps, these horizontal magnetic fields and resulting currents are further smeared out.As a result, the corresponding electric currents outside of active regions are not accurate enough to allow the nlf-op and nlf-gr extrapolations to recover the sheared magnetic fields in filament channels.In fact, as we have seen, the nlf-gr model assumed j r = 0 in the weak field regions of the photosphere, owing to this uncertainty.
The resolution of the BBSO Hα image in Figure 8 is sufficient to show filaments themselves, as dark structures, but is not sufficient to show empty filament channels.Thus absence of a filament can not be taken to mean absence of sheared magnetic fields.Filament channels may also be identified from alignment of coronal cells in AIA images, particularly in the 193 Å channel (Sheeley and Warren 2012).Careful inspection of SDO/AIA movies of EUV emission, especially using composite images processed by the Morgan and Druckmüller (2014) technique, was used here to identify which filament channels to energize in the mhd-mas model (see Appendix A).This method identifies many filament channels that do not show up clearly in Figure 8(h) and (i).As a case in point, consider the location labelled A in Figure 8(g), where a sheared field was inserted in the mhd-mas model.There is only a small amount of Hα filament material visible at the west end of the channel in Figure 8(h), but analysis of the EUV observations at higher resolution suggests the possible presence of a long East-West filament channel all along this polarity inversion line, bending northward at its eastern end as in the mhd-mas model.A sheared filament channel is also present here in the mf model.

Open Magnetic Flux
The open magnetic field lines are the source of the solar wind, so represent the output of the models as far as the heliosphere is concerned.In Section 3.1 (Figure 4a), we compared the open magnetic flux Φ(2.5 R ⊙ ), as defined in Equation (2).The numerical values for each model are given explicitly in Table 3, along with an observational estimate based on in situ OMNI data.To make this estimate, daily averages of the basic hourly OMNI data were obtained from the GSFC/SPDF OMNIWeb interface at http://omniweb.gsfc.nasa.gov,giving B r (R E ).Averaging over 27 days centered on the eclipse time, and assuming a uniform distribution of magnetic flux over latitude at R E = 1 AU, we estimate the equivalent open flux shown in Table 3 as Mx.As has previously been found during relatively active periods of solar activity (see Linker et al. 2017), this observed value is much higher than that predicted by models -either potential or non-potential.Our results are consistent with their findings.The additional coronal currents in the mf and mhd-cese models, in particular, enhance the open flux significantly, but it still remains below the level inferred from OMNI observations.The reasons for this discrepancy are not yet understood, and are discussed by Linker et al. (2017).It is possible that this may be partially heliospheric in origin, if there are regions in the solar wind where the interplanetary magnetic field folds back on itself and is thus over-counted in the B r (R E ) measurements (Owens, Crooker, and Lockwood 2013).Finally, we consider the spatial distribution of open and closed magnetic field in the models.In addition to indicating the possible source regions of the solar wind, this is important because the open or closed nature of the magnetic field at a given location depends sensitively on both the input data and the distribution of electric currents in the model.Comparing the open field footpoints to the locations of observed coronal holes therefore provides a further observational constraint.
Figure 9 shows the locations of open magnetic field line footpoints on r = R ⊙ for the different models, along with -in panel (h) -a persistence map of observed coronal holes in EUV.This map was built up synoptically over Carrington rotation 2161 (2015 February 28 to 2015 March 27) by thresholding full-disk SDO/AIA 193 Å EUV images, at a cadence of 12 hours.The procedure is described in more detail by Lowder, Qiu, and Leamon (2017), although during this particular time period no far-side EUV data were available.Of course, this should be taken only as a lower bound on the areas of observed open field rather than a true measure.This is because, while a significant fraction of open flux originates from coronal holes, the footpoints of open field can also be bright, particularly if they lie in active regions.Indeed, quantifying the amount of open flux not located in coronal holes is an important observational challenge for reconciling the models and observations in future.
Given the sensitivity of the open/closed footpoint regions as well as the variety of boundary conditions used, there is reasonable agreement between most of the models and with the observed persistence map.Robust features are labelled in Figure 9(h) and include: a negative-polarity polar coronal hole in the Southern Hemisphere (A) but no corresponding hole in the Northern hemisphere; a narrow equatorward extension of the Southern polar hole (B and C); a long positivepolarity hole in the Northern hemisphere (D); a north-south oriented positive-Fig.9 Footpoints of open magnetic field lines shown on r = R ⊙ , and their polarity (red positive, blue negative).Panel (h) shows the observed coronal hole persistence map (see text) using AIA data.In this map, the colors again indicate polarity, with density indicating persistence, i.e., the number of days for which a coronal hole was observed at that location.
polarity hole (E); and a more compact positive-polarity hole near the equator in the Southern Hemisphere (F).
The ffe model differs significantly from the others, as was evident from its magnetic field structure in Figure 7.The mf model reproduces many of the observed coronal holes, although some (particularly E) are shifted in position.However, it has in number of additional open field regions, with significantly more open magnetic field of both polarities at low and medium latitudes.Indeed this additional open flux was seen in Figure 4(a).Partly this additional open field arises from the opening up of streamers discussed in Section 3.3 and seen in Figure 7. Another reason for enhancement is the ejection of magnetic flux ropes in the mf model (one such eruption is responsible for the U-loops in Figure 7e).But the difference from other models is also partly due to the different photospheric distribution of B r (Figure 3), owing to the model being driven by a flux transport model rather than directly from observed magnetograms (Appendix B).A good example is the large extension of the negative-polarity hole into the Northern hemisphere around Carrington longitude 300 • to 360 • (labelled G in Figure 9e).This hole is likely to reduce in size once the new active regions emerge that are already present in the other models.
Figure 10 shows the distribution of B r on the outer boundary r = 2.5 R ⊙ in each model.Here all of the field lines are "open" according to our definition.At this radius, the magnetic field is insensitive to small-scale differences in the input magnetic maps, and there is reasonable agreement between nlf-op, nlf-gr , mhs and mhd-mas.These models all inherit the same basic open/closed topology from the potential field used as their initial condition before the injection of electric currents at the base.In the ffe, mf and mhd-cese models, the heliospheric current sheet (boundary between positive and negative B r ) has a more complex shape, with a disconnected loop in the mf model.The reasons for these differences are the same as in Figure 9.As mentioned above, the ffe model has a complex magnetic structure that likely reflects the incomplete relaxation of the model.These results do suggest that we should not take the potential field topology for granted, although some of these differences arise from the different photospheric B r distribution in the mf and mhd-cese models; indeed, the coronal hole map comparison suggests that the difference may be over-emphasized in these models.The additional open flux in the mhs, mhs and mhd-cese models, compared to the others, is also clear in this plot.Note that this could be reduced in the mhs model by reducing the parameter a, so may not be particularly significant in that case.

Conclusions
Having analysed seven different non-potential models of the coronal magnetic field on 2015 March 20, we now draw some overall conclusions.
The initial impression from Figure 7, for example, is one of significant disagreement between different models.This is true particularly in regard to the overall magnetic structure.This disagreement arises from several sources: the input data used by each modeller, the coronal modelling techniques themselves, and the outer boundary conditions.
There are two fundamental limitations with currently available magnetogram input, not including the differences between magnetograms from different instruments and observatories (Riley et al. 2014), which we avoid here by focusing on SDO/HMI.One limitation is the lack of co-temporal full surface coverage, which necessitates the use of synoptic maps that combine observations taken at different times.The different models achieve this in different ways, as described in Section 2, and this leads to different maps of the magnetic field across the whole solar surface (Figure 3), with corresponding differences in coronal magnetic topology.One consequence of the use of synoptic maps is that the majority of models (except for mf and mhd-cese) use data taken after the eclipse time, meaning that this exercise is not equivalent to a prediction of the eclipse magnetic field which could have been made in advance.The B r distribution in the mf model differs particularly from the others because it uses a surface flux transport model to evolve the photospheric magnetic field over a period of months, rather than inserting magnetogram data more directly.This has the advantage of allowing the build-up of free magnetic energy in weak-field regions, but does mean that there are resulting differences in the magnetic structure.The use of synoptic data also means that the magnetogram input near the limbs -particularly the East limb -is out-of-date, and this can particularly affect comparisons with eclipse images.The proposed space mission to the L5 Lagrange point would greatly improve the situation at the East limb, provided that a magnetograph were included on board.
The second limitation is with vector magnetograms.Firstly, the signal-to-noise ratio remains low in the transverse components of B, so that the nlf-op, nlf-gr and ffe models do not reproduce the sheared magnetic fields in filament channels.And the fact that the photospheric vector magnetograms do not satisfy j × B = 0 leads to problems in particular with the ffe model, preventing it from reaching equilibrium.Further refinement will be needed before this model can be of practical use for the solar corona.In future, it is hoped that this issue will improve if and when upper-chromospheric magnetogram observations become available.
In the coronal volume itself, we have shown that the different models vary in their degree of non-potentiality, as measured by either electric currents or free energy.The models with greater free energy achieve this with a variety of different current distributions: currents may be concentrated in active regions (nlf-gr ), in filament channels (mhd-mas and mf ), or may be more distributed throughout the corona (mf and mhs).This aspect of the models is perhaps the most difficult to calibrate against observations, since no direct observations of coronal electric currents are available.The presence of filament channels, for example, illustrates the importance of the gradual build-up of coronal electric currents over time.Although these are not reproduced in the nlf-gr , nlf-op, mhs and mhd-cese models, this is not due to a fundamental limitation of the equations used to model the magnetic field in the corona, but due to limitations in the available boundary data.In particular, we have shown with the mhd-mas model how static models can be further energized by the injection of electric currents in filament channels, and this could be applied in future with the other models.This hybrid approach of informing static models with the results either from simplified time-evolving models (like mf ) or from additional coronal observations could well be a useful one in the absence of improved magnetogram input data.
The outer boundary conditions also have a significant impact, and deserve greater attention.With the exception of mhd-cese, none of the models presented here couples physically to a solar wind solution beyond 2.5 R ⊙ .This is because the increasing plasma-β in that region no longer allows for the use of a purely magnetic model.Several of the models impose an artificial source surface at 2.5 R ⊙ , as in the common PFSS model, and this leads to inaccuracies in the streamer structure, and potentially also the open flux.The mf model instead uses a radial outflow boundary condition to mimic the effect of the solar wind, but this seems to be inflating the field too much in this particular case, at least in terms of open field footpoints (though interestingly not in the height of the streamers).The corresponding currents higher in the corona change the open-closed magnetic topology significantly.Overall, this outer boundary is an important problem that requires more sophisticated MHD modelling that includes plasma thermodynamics.Such simulations have been performed with the full-MHD version of the mhd-mas model (e.g., Lionello, Linker, and Mikić 2009;Downs et al. 2013), which can describe the solar wind.
All of this being said, there are also areas of broad agreement between many of the models.For example, while they have very different input grid resolutions on r = R ⊙ , this (in itself) does not affect the estimated open flux and topology of the heliospheric current sheet on r = 2.5 R ⊙ , which arises rather from the differences between the coronal modelling approaches.On r = R ⊙ , the footpoint regions of open magnetic field show similarities in all models, and there is agreement that the strongest currents are within the active regions.Among those models with significant free magnetic energy, there is general agreement on the ratio E/E p ≈ 1.4 to 1.5.And the locations of closed field streamers are broadly in agreement, though not their height and shape.
From this study it is clear that all of the models have positive aspects that agree with observations, but other aspects that do not match so well.Much of this can be related to the distribution of electric current both within and outside of active regions.At present, nonlinear force-free extrapolations such as nlf-op or nlf-gr are best at representing the structure of active regions where reliable vector magnetic field input data are available.But accounting for the free energy outside of active regions is currently possible only with time-evolving models such as mf .Yet, in these models, it is too computationally expensive to account for the full plasma thermodynamics, something that has not been considered here but is already possible in state-of-the-art full-MHD models, albeit for static configurations.
What is clear is that this is an exciting time for coronal magnetic field modelling, with progress on several fronts but much still to do.Rather than simply waiting for better magnetogram data, our comparisons with currently available observations -though qualitative -do suggest that these indirect observational constraints could be better used to optimize the models.The challenge is to do this systematically.The ideal model would match EUV observations of filament chan-nels and coronal loops, the positions of white-light streamers, and the locations of observed coronal holes.A more sophisticated approach would involve forward modelling of actual observed emission, but we suggest that much can already be learned from morphological comparisons.A Details of the mhd-mas model The MAS time-dependent magnetohydrodynamic (MHD) model (Mikic and Linker 1994;Linker and Mikic 1995;Mikić et al. 1999) of the global solar corona has been in development for over two decades.
The simulation mhd-mas illustrated in this paper uses a simplified zero-beta version of the model, in which pressure and gravity forces are neglected, in order to emphasize the nonpotential aspects of the corona, in particular the sheared fields expected in filament channels.
An energization mechanism was applied to create sheared and twisted magnetic fields and flux ropes in the observed filament channels.Unlike the more sophisticated version of the MHD model, this zero-beta version is not able to include explicitly the effect of the solar wind, but in the future it is hoped to include this energization capability into the full MHD model.The filaments (visible for example in Figure 8h) are thought to be supported by twisted flux-rope-like magnetic fields.These were introduced in an ad hoc manner as follows: 1. Start with a PFSS solution but add additional flux to Br on the boundary r = R ⊙ , to compensate for that lost by the cancellation in Step 3. 2. Emerge transverse magnetic field by imposing a transverse electric field Et = ∇tΦ at the boundary r = R ⊙ .Here Φ is chosen to be localized to the filament channels, and reverses sign across the PILs, and ∇t is the transverse gradient.Note that this does not change Br. 3.After the transverse magnetic field emergence is complete, cancel flux at the PILs by applying a transverse electric field of the form Et = ∇t × ψr, so as to form more flux-rope like field lines out of highly sheared arcades.The function ψ is driven by the change in Br, By the end of this process, Br will once again match the magnetogram, thanks to the additional flux added in Step 1.
The total amount of added flux, which was chosen by trial and error, corresponds to 20% of the flux in the magnetogram (over the whole Sun).Since the added flux is localized to filament channels, it corresponds to a greater fraction of the flux in the filament channels (on the order of 40%).The filament channel regions at which the field was twisted in this manner were chosen in an ad hoc manner, guided by both the results of the mf model and the locations of filament channels inferred from SDO/AIA EUV movies during the evolution surrounding the eclipse.Since the mf model assimilated the emerging active regions from HMI observations over time, it was able to give the best estimate of the chirality of the fields in the filament channels.By chirality here we mean the direction of the transverse field introduced.According to the usual definition, a sinistral filament has an axial field that points to the left when looking towards the PIL from the side with Br > 0, and vice versa for a dextral filament.In the technique described above, choosing Φ with the same sign as Br produces a dextral filament, whereas Φ with opposite sign to Br produces a sinistral filament.In some cases, when the chirality changed within a single filament, a single chirality was used to simplify matters.Figure 11 shows the filament channels where the magnetic fields were energized, with the chosen chirality (as deduced from the mf model), and the corresponding Φ potential used.The magnetic energy of the initial PFSS field (with the observed Br plus the added flux) was 2.13 × 10 33 ergs.The energy in the PFSS field corresponding to the observed Br was 1.72 × 10 33 ergs.The energy in the final field, with twisted fields in the filament channels, matching the observed Br, was 2.52 × 10 33 ergs, so that the final field had an energy 47% above the corresponding PFSS field.

B Details of the mf model
For the mf calculation presented here, the coronal magnetic field was evolved for 200 days up to the eclipse date, allowing the self consistent build-up of electric currents and free magnetic energy, driven by the photospheric evolution (both footpoint motions and flux emergence).The simulation was initialized with a potential field extrapolation for 2014 September 1, and a fully non-potential corona was arrived at after 6-8 weeks (as measured by average current density), well before the eclipse date.
Compared to previously published simulations, a different evolution of Br at the solar surface r = R ⊙ was used.Namely, this was derived from the Advective Flux Transport (AFT) model of Upton and Hathaway (2014b,a) and Ugarte-Urra et al. (2015), which assimilates HMI magnetograms on the visible disk and advects the distribution over the full solar surface using differential rotation, meridional flow and small-scale convective flows.For driving the magnetofrictional model, it is actually the horizontal electric field that is required.Since this is unknown in the regions assimilated from HMI, the AFT model could not be used directly to drive the 3D simulations.Instead, it was used to determine the locations of new bipolar magnetic regions (BMRs), which were then inserted into a simpler flux transport model which also includes differential rotation and meridional circulation, but approximates the small-scale convective flows by a simple supergranular diffusion term.The resulting smoother distribution of the photospheric field is evident in Figure 3(e).
The BMR properties were derived from the AFT model by a three-stage automated procedure: 1. Compare successive Br maps (once every 24 hrs) to identify new BMRs (leading to 197 distinct BMRs).2. Compute the BMR properties (location, size, magnetic flux and tilt angle) once per hour, from 48 hrs before to 120 hrs after initial identification, and determine the time of maximum flux.Insert the (3D) BMR in the new simulation at this time with the corresponding properties.3. Check the resulting Br evolution manually, and correct the BMR properties to best reproduce their times of emergence and also the structure of multiple-BMR activity complexes.
In the run shown here, each BMR was given a twist when inserted into the 3D simulation -so as to model the active region helicity -with magnitude β = ±0.4(see Yeates, Mackay, and van Ballegooijen 2008).

Fig. 1
Fig. 1 An extreme-ultraviolet and white-light composite of the corona on 2015 March 20.The white light image is a combination of 29 exposures made from Longyearbyen, Svalbard, and aligned with sub-pixel precision using the Phase Correlation technique (Druckmüller 2009).This is overlaid with a combination of 171, 193 and 211 Å (red, green, blue respectively) channels from the AIA instrument on the Solar Dynamics Observatory satellite.Each of the channels were processed using Multiscale Gaussian Normalization (MGN, Morgan and Druckmüller 2014) prior to combination.

Fig. 3
Fig. 3 Radial magnetic field Br(R ⊙ , θ, φ) in each model, over the full solar surface.Here the horizontal axis shows Carrington longitude, with Central Meridian at eclipse time located at 99 • .The color tables are saturated at ±50 G (red positive, blue negative).Panel (h) repeats Br from the SDO/HMI map in Figure 2, for ease of comparison.

Fig. 4
Fig. 4 Magnetic flux (a) and energy (b) of the models.Shown are (a) unsigned radial magnetic flux Φ(r) as a function of height, and (b) total magnetic energy E versus potential magnetic energy Ep.The solid line indicates E = Ep.
has a flux Φ(R ⊙ ) = 5.8 × 10 23 Mx at its original resolution.The difference in resolution between the models becomes insignificant above 1.5 R ⊙ , where the behaviour is dominated by low-order spherical harmonics.Accordingly, most of the models predict a similar open magnetic flux of around 3 × 10 22 Mx.The exceptions are mhs, mf and mhd-cese, where the open flux is inflated due to the presence of significant volumetric currents in the corona -up to 6 × 10 22 Mx in the case of mf .The mf model includes the ejection of magnetic flux ropes which also increases the open flux.These open flux values, and the open magnetic field distribution, will be further discussed in Section 3.5.

Figure 4
Figure 4(b) shows the total magnetic energy of each model between r = R ⊙ and r = 2.5R ⊙ , given by

Fig. 5
Fig. 5 Measures of the (a) perpendicular and (b) parallel currents for each model, as functions of height.Panel (c) shows the ratio of parallel to perpendicular measures.

Fig. 6
Fig.6Magnitude of the vertical electric current density, log 10 |jr|, at r = 1.02R ⊙ in each model, over the whole spherical surface.The color tables are saturated at 10 −10 G cm −1 (black) and 10 −6.5 G cm −1 (white).Panel (h) shows the same quantity computed at r = R ⊙ directly from the SDO/HMI synoptic map in Figure2(boxcar smoothed to reduce noise).The same color scale is used.

Fig. 7
Fig. 7 Magnetic field structure of each model in the plane-of-sky, as viewed from the Earth at the (approximate) eclipse time.(The density of magnetic field lines is arbitrary and does not correspond to field strength.)Red/blue shading shows Br on r = R ⊙ (saturated at ±50 G).Panel (h) shows an EUV 174 Å image of the real corona from PROBA2/SWAP (with median stacking and stray-light correction as described by Seaton et al. 2013).Panel (i) shows an image of the corona in the Fe XIV 5305 Å line from the Habbal eclipse expedition, processed by the MGN algorithm of Morgan and Druckmüller (2014) (with parameter h = 0.9) to bring out the streamer structure.Several angular regions are labelled in panel (h) for ease of reference in the text.

Fig. 8
Fig. 8 Magnetic field lines for each model chosen to show the low coronal magnetic structure above each polarity inversion line on the visible disk.Red/blue shading shows Br on r = R ⊙ (saturated at ±50 G).Panel (h) shows an Hα image from Big Bear Solar Observatory (taken at 1733UT on the eclipse day).Panel (i) shows a composite of the 171, 193 and 211 Å channels (red, green, blue respectively) from AIA, processed as in Figure 1.

Fig. 10
Fig. 10 Radial magnetic field at r = 2.5R ⊙ in each model, over the full spherical surface.The color tables are saturated at ±0.4 G.

Fig. 11
Fig. 11 Energization of filament channels in the mhd-mas model.Panel (a) shows radial magnetic field Br at r = R ⊙ for the observed HMI field with the added flux.The colored contours are the PILs in the PFSS model at r = 1.02R ⊙ , with thick colored line segments showing the filament channels in which the field was energized.Panel (b) shows the energizing Φ potential.The chirality selected for the magnetic field in these filament channels is shown in (b), with "S" denoting sinistral and "D" denoting dextral filaments.

Table 1
Summary of the non-potential models presented.Mikić et al. (1999) -MAS * See Appendix B. † Filament channel locations based on mf model (see Appendix A).

Table 2
Domains, resolutions and energies * of the non-potential models.

Table 3
Unsigned open magnetic flux Φ(2.5 R ⊙ ) of the non-potential models.