Issue 
A&A
Volume 644, December 2020



Article Number  A57  
Number of page(s)  7  
Section  The Sun and the Heliosphere  
DOI  https://doi.org/10.1051/00046361/202039079  
Published online  01 December 2020 
Preprocessing of vector magnetograms for magnetohydrostatic extrapolations
MaxPlanckInstitut für Sonnensystemforschung, JustusvonLiebigWeg 3, 37077 Göttingen, Germany
email: zhu@mps.mpg.de
Received:
31
July
2020
Accepted:
14
September
2020
Context. Understanding the 3D magnetic field as well as the plasma in the chromosphere and transition region is important. One way is to extrapolate the magnetic field and plasma from the routinely measured vector magnetogram on the photosphere based on the assumption of the magnetohydrostatic (MHS) state. However, photospheric data may be inconsistent with the MHS assumption. Therefore, we must study the restriction on the photospheric magnetic field, which is required by the MHS system. Moreover, the data should be transformed accordingly before MHS extrapolations can be applied.
Aims. We aim to obtain a set of surface integrals as criteria for the MHS system and use this set of integrals to preprocess a vector magnetogram.
Methods. By applying Gauss’ theorem and assuming an isolated active region on the Sun, we related the magnetic energy and forces in the volume to the surface integral on the photosphere. The same method was applied to obtain restrictions on the photospheric magnetic field as necessary criteria for a MHS system. We used an optimization method to preprocess the data to minimize the deviation from the criteria as well as the measured value.
Results. By applying the virial theorem to the active region, we find the boundary integral that is used to compute the energy of a forcefree field usually underestimates the magnetic energy of a large active region. We also find that the MHS assumption only requires the x, ycomponent of net Lorentz force and the zcomponent of net torque to be zero. These zero components are part of Aly’s criteria for a forcefree field. However, other components of net force and torque can be nonzero values. According to new criteria, we preprocess the magnetogram to make it more consistent with the MHS system and, at the same time close, to the original data.
Key words: Sun: magnetic fields / Sun: photosphere / magnetohydrodynamics (MHD)
© X. Zhu et al. 2020
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Open Access funding provided by Max Planck Society.
1. Introduction
The forcefree assumption ∇ × B = α(r)B has long been the basis for a magnetic field extrapolation from the solar photosphere to the corona (Wiegelmann & Sakurai 2012). As a result, a number of forcefree field extrapolation techniques have been developed over half a century. The simplest case is the potential field or currentfree field extrapolation in which α = 0 (Schmidt 1964; Schatten et al. 1969). A next step with a spatially constant α is the linear forcefree field extrapolation (Chiu & Hilton 1977; Seehafer 1978; Alissandrakis 1981). At last, a model with nonconstant α is called the nonlinear forcefree field (NLFFF) extrapolation. The NLFFF extrapolations include the (1) upward integration method (Nakagawa 1974; Wu et al. 1985, 1990; Cuperman et al. 1991; Demoulin & Priest 1992; Song et al. 2006), (2) GradRubin method (Grad et al. 1958; Sakurai 1981; Amari et al. 1997, 1999, 2006; Wheatland 2004, 2006), (3) relaxation method (Mikic et al. 1988; Roumeliotis 1996; Valori et al. 2005; Jiang & Feng 2012; Fan et al. 2011; Guo et al. 2016), (4) optimization method (Wheatland et al. 2000; Wiegelmann 2004; Wiegelmann & Inhester 2010), (5) boundaryelement method (Yan & Sakurai 2000; Yan & Li 2006; He & Wang 2006) and (6) forwardfitting method (Malanushenko et al. 2012; Aschwanden 2013, 2016).
In a forcefree model, some necessary conditions of the magnetic field on the boundary have to be fulfilled. Molodenskii (1969), Molodensky (1974), and Aly (1984, 1989) obtained several integral relations on the boundary corresponding to the vanishing net magnetic force and vanishing torque, respectively. If these integral relations are not fulfilled then the boundary is not consistent with a forcefree field. Based on this, Wiegelmann et al. (2006) proposed a preprocessing algorithm to modify the magnetogram within the error margins of the measurements to minimize the net magnetic force and torque. The resulting vector magnetogram is more suitable for a forcefree extrapolation. Further developments include using the method of simulated annealing (Fuhrmann et al. 2007), extending to the spherical geometry (Tadesse et al. 2009), adding a new term concerning chromospheric longitudinal fields (Yamamoto & Kusano 2012), and dealing with the potential and nonpotential components separately (Jiang & Feng 2014).
With increasing spatial resolution of the measured vector magnetogram, we can study the magnetic field in the lower solar atmosphere in detail. However, the forcefree assumption is not valid anymore in this layer as the plasma β is much larger than it is in the corona (Gary 2001). A straightforward approach is to take into account the plasma in the extrapolation. Zhu et al. (2013, 2016), and Miyoshi et al. (2020) proposed the magnetohydrodynamic relaxation method to obtain a magnetohydrostatic (MHS) solution with a nonforcefree layer near the bottom boundary. Gilchrist & Wheatland (2013) and Gilchrist et al. (2016) extended the Grad et al. (1958) method to compute the MHS equilibria. Wiegelmann & Neukirch (2006) used the optimization method to solve the MHS equations without gravitational force in the Cartesian coordinate system. Wiegelmann et al. (2007) extended the code with gravitational force in spherical coordinate system. Recently, we extended the optimization method by introducing the gravitational force (Zhu & Wiegelmann 2018) in the Cartesian coordinate system, tested the code with an realistic radiative MHD simulation (Zhu & Wiegelmann 2019), and also applied the code to a SUNRISE/IMaX dataset (Zhu et al. 2020). It is worth noting that the new algorithm ensures the positive definiteness of gas pressure and mass density.
As a result of the MHS assumption, we can also define several integral relations just as we did for the forcefree field. A preprocessing algorithm (extended from the NLFFF case) is proposed to modify the vector magnetogram within the error margins of the measurement. The resulting magnetogram is expected to be more consistent with the assumption of a MHS extrapolation.
The remainder of the paper is organized as follows: in Sect. 2 we define the integral relations for a MHS system and also apply the virial theorem to an active region. In Sect. 3 we describe the optimization algorithm to derive consistent boundary conditions for a MHS extrapolation. In Sect. 4 we apply the algorithm to an example of the observed vector magnetogram. In Sect. 5 we draw conclusions.
2. Boundary integrals of a MHS system
A MHS equilibrium can be described as
where B, p, and ρ are the magnetic field, plasma pressure, and plasma density, respectively. We note that B, p, and ρ have been appropriately normalized to simplify the equation. For example, for a study that includes photosphere, the following normalization constants are convenient: ρ_{0} = 2.7 × 10^{−1} g cm^{−3} (density), T_{0} = 6 × 10^{3} K (temperature), g = 2.7 × 10^{4} cm s^{−2} (gravitational acceleration), cm (length), dyn cm^{−2} (plasma pressure), and G (magnetic field), where ℛ is the ideal gas constant.
As discussed in many papers (e.g., Chandrasekhar 1961; Molodenskii 1969; Molodensky 1974; Aly 1984), the above force balance equation may be written in the form analogous to equations of elasticity as follows:
where is the Maxwellian tensor and f includes forces such as pressure gradient ( − ∇p) and gravitational force ().
2.1. Equations of momenta
Assuming a Cartesian coordinate, integrating Eq. (3) over the volume V with surface S, we get
We apply these three relations to an isolated active region with the photosphere as the bottom boundary. As to a closed current system, the magnetic field falls off as 1/r^{3}. Therefore as r → ∞, the lateral and top surface integrals with magnetic terms approach zero. We note that in ∫ − ∇pdv = −∮_{S} pds the transverse surfaces make zero net contribution since, without magnetic fields, the forces acting on the opposite side boundaries by the external plasma pressure are equal and opposite. Then we have net Lorentz force restrictions written as
where S_{1} is the bottom boundary. We note that Eqs. (7) and (8) are exactly the same as those for a forcefree field. However, that is not the case vertically. Equation (9) implies that the difference between the pressure at the bottom boundary and the weight of the plasma above is compensated by the Lorentz force.
The net Lorentz force torque restrictions can also be obtained in a similar way with the cross product of Eq. (3) and r as follows:
We see Eq. (12) is exactly the same as that for a forcefree field, while that is not the case in two transverse directions. Equations (10) and (11) imply that the plasma may induce rotational moments relative to the x or y axes.
2.2. Virial theorem
Multiplying Eq. (3) by and integrate over volume V with a surface S, we derive the virial theorem
which relates the magnetic energy and forces inside the volume to the integral on the surface. If the electric currents are closed, the magnetic field falls off as 1/r^{3}, and as r → ∞ the surface integral approaches zero. Therefore the magnetic energy vanishes if f = 0. That means a forcefree field does not exist in an isolated current system (Chandrasekhar 1961).
Plugging into Eq. (13), we have
where is the total energy including magnetic energy , internal energy and gravitational potential energy −ρz. If then the lefthand side (LHS) of Eq. (14) becomes exactly the total energy. Equation (14) relates the total energy and pressure in the volume to the integral on the surface.
2.3. Apply to an active region
Let us focus on an active region. Assume a semiinfinite halfspace with the bottom boundary on the photopshere, the magnetic field of an isolated active region decreases with distance as 1/r^{3}, thus the surface integral in Eq. (13) on the side and top boundaries will approach zero. Then we get:
where S_{1} is the bottom boundary.
2.3.1. Effect of translation on the integration
In the forcefree case, the integral at the righthand side (RHS) of Eq. (15) is invariant under translation. However, it seems from the LHS of Eq. (15) that this nice property is not fulfilled any more in a MHS system. Suppose a translation: (x, y, z)→(x′ + Δx, y′ + Δy, z′ + Δz). Then B′(r′) = B(r), dx′ = dx, and dy′ = dy. Thus we have the RHS of Eq. (15) in the new coordinate system as follows:
where Δr = (Δx, Δy, Δz). We note that Eqs. (7) and (8) are used to infer formula (18) from (17). Thus the integration on the bottom boundary is invariant under the translation in the same horizontal plane with Δz = 0.
2.3.2. Estimate the magnetic energy
Suppose r = 0 on the bottom boundary, Eq. (15) becomes
As with a forcefree field, the integral at the RHS of Eq. (19) is often used to estimate the magnetic energy in the semiinfinite volume. However, the integral is not equal to the magnetic energy in a MHS equilibrium. In the following analysis, we try to judge, in a MHS equilibrium, whether the integral is an underestimate or an overestimate of the magnetic energy.
Consider the case of an axisymmetric monopole sunspot (see Fig. 1). The field line approaches a radial direction since the gas pressure becomes small. Near the spot, however, the field line has to bend to compensate the pressure gradient. According to Sect. 2.3.1, the boundary integration of Eq. (19) is invariant under the translation on the bottom plane. For the sake of convenience, we chose r = 0 at the center of the spot. As gas pressure is depleted in the spot, the negative pressure gradient usually points to the interior of the spot. To compensate this force, a reversed Lorentz force is generated transversely. Then a vertically downward component of the Lorentz force is required to make sure that the Lorentz force is perpendicular to the magnetic field. The negative net Lorentz force has been confirmed in most active regions by statistical studies (Moon et al. 2002; Tiwari 2012; Liu et al. 2013; Liu & Hao 2015). According to Fig. 1, far away from the spot, the contribution of f ⋅ r vanishes as the field line becomes radial. However, near the spot, the contribution of f ⋅ r is always negative. Thus we get, in a monopole spot case,
Fig. 1. Typical interaction of the magnetic force and plasma forces in an active region with single compact polarity, in a MHS state. 
The above inequality also holds in an active region with multiple spots. See Fig. 2 with two spots, the f ⋅ r of the left spot is equal to that of the monopole spot case. As to the right spot, with the coordinate transformation r = r_{1} + R, we have
Fig. 2. Typical interaction of the magnetic force and plasma forces in an active region with double compact polarities, in a MHS state. 
where R is transverse while ∫_{V} fdv is vertical. Term R⋅∫_{V} fdv vanishes is because the horizontal vector R is perpendicular to the vector ∫_{V} fdv, which is vertical due to the axisymmetry assumption. It is worth noting that, in the multiple spots case, the magnetic field line deviates from the radial direction because of the attraction from other spots. This effect, however, is minor near the spot. In regions far away from spots where field lines are not radial any more, the magnitude of f ⋅ r decreases with altitude rapidly because f decreases exponentially (depending on temperature) while r can only grow linearly. Thus the analysis expression above still works. Therefore, we conclude that, as in an active region with sunspots, the surface integral of Eq. (19) is usually an underestimate of the magnetic energy of the active region.
3. Preprocessing method
To see if a vector magnetogram can be served as the boundary condition for a MHS system, three parameters are proposed (similar to those for a forcefree field Wiegelmann et al. 2006) as follows:
A vector magnetogram is suitable for a MHS extrapolation at least if ϵ_{flux} ≪ 1, ϵ_{force} ≪ 1 and ϵ_{torque} ≪ 1.
Since the vector magnetogram sometimes do not fulfill the aforementioned criteria, a preprocessing procedure is required to modify the data within the freedom of the noise. The algorithm is based on the preprocessing method that was developed by Wiegelmann et al. (2006). As the restrictions on the boundary values for a MHS equilibrium are weaker compared with those for a forcefree field, we need to change the preprocessing procedure accordingly. Nonzero values of zcomponent of net Lorentz force and x, ycomponent of net torque are allowed to exist in a MHS equilibrium. Therefore these three numbers should be retained during the preprocessing process.
To do so, we define the functional
where
are the zcomponent of net Lorentz force and the x, ycomponent of net Lorentz torque, respectively. The summation is over all p grid nodes on the photosphere. The weighting factors μ_{n} are as yet undetermined. The terms L_{1} and L_{2} correspond to the net force and net torque constraints. The term L_{3} measures the deviation between the original data and the preprocessed data. The term L_{4} controls the smoothing. For computational reasons, sufficiently smooth data are necessary for an optimization to obtain a good solution. The new algorithm differs from that for the forcefree field (Wiegelmann et al. 2006) by introducing three quantities: a_{0}, a_{1}, and a_{2}, which ensure that the preprocessing does not change the corresponding integration values.
The strategy of preprocessing is to use the gradient descent method to minimize L, and meanwhile make all L_{n} small as well. The magnetic field is optimized as follows:
with a gradient descent method
The three functional derivatives at node (q) are defined as
Smoothing was performed for all three components. Effects from terms that have mixed products of vertical and transverse magnetic field components in functional Eqs. (30–33) are not included when evaluating B_{z}. This is designed for the fact that B_{z} is measured with much higher accuracy than B_{x} and B_{y} (Martínez Pillet et al. 2011).
4. Application to SUNRISE/IMaX data
For the test we used a combined vector magnetogram in which the SUNRISE/IMaX data (Martínez Pillet et al. 2011; Solanki et al. 2017) is embedded in the HMI data (Scherrer et al. 2012). We have used this dataset to extrapolate the magnetic field as well as the plasma using two approaches (Wiegelmann et al. 2017; Zhu et al. 2020). The top panels of Fig. 3 show the original vector magnetogram within the IMaX field of view (FOV). The combined data have a flux imbalance of ϵ_{flux} = 0.013 (almost balanced). The MHS criteria are not fulfilled to some extent with ϵ_{force} = 0.091 and ϵ_{torque} = 0.066. For comparison, Aly’s criteria for forcefree field are largely violated with ϵ_{force} = 0.29 and ϵ_{torque} = 0.32.
Fig. 3. Vector magnetogram of IMaX measured on 2013 June 12 at 23:48 UT. Top: original data. Bottom: preprocessed data. 
Before the magnetogram is preprocessed, we first need to choose the appropriate μ_{n}. There are four parameters of μ_{n}. Only three of them are independent since only the ratio of the parameters really counts. We further assume μ_{1} = μ_{2}D^{2} ≡ μ_{12}, where D is the average edge length of the magnetogram. Thus we could define L_{12} = L_{1} + D^{2}L_{2}. This assumption gives the same weight to the momentum and torque constraints. As only the ratio of the parameters counts, we specify μ_{12} = 1/B_{ave}, where B_{ave} is the average magnetic field strength in the magnetogram. Then only two independent parameters remain (μ_{3} and μ_{4}). A survey of the two parameters for the combined magnetogram shows that, as was also found in Wiegelmann et al. (2006), log(L_{3}) and log(L_{4}) are almost determined by the ratio of μ_{3} to μ_{4}, while log(L_{12}) depends on the magnitude of μ_{3} and μ_{4}.
With a deviation of the magnetic field value (i.e., a finite L_{3}), we obtain a smoothed solution that satisfies the criteria of a MHS system. The deviation is tolerable as long as L_{3} does not exceed noise level of the magnetogram, that is, L_{3} = 2.2 × 10^{−9} in this case. The noise is retrieved from the HMI dataset (Hoeksema et al. 2014). Figure 4 shows optimal μ_{3} and μ_{4} combinations at which log(L_{3}) equals to noise level of the magnetogram. As μ_{3} and μ_{4} decrease, L_{12} converges. Table 1 shows Lvalues as well as three summations of the initial data and the preprocessed data. The three summations are not changed during the preprocessing owing to the special design of the algorithm for fixing them. An improvement on L_{12} is obvious in that L_{12} is reduced by five orders of magnitude, which enforce a good compliance with the MHS criteria. The quantity L_{3} has to be finite because we allow the field values to deviate from the observed values. Even though the smoothness is hardly visible in Fig. 3, L_{4} after preprocessing is 1/20 smaller than it was before preprocessing. Now we get ϵ_{force} = 7.0 × 10^{−5} and ϵ_{torque} = 8.7 × 10^{−5} with an optimal values of μ_{3} = 1.0 × 10^{−3} and μ_{4} = 8.9 × 10^{−4}, which are more suitable for a MHS extrapolation.
Fig. 4. Test runs with various μ_{3} and μ_{4} combinations. The quantity L_{3} of each test reaches the noise level of the magnetogram. The optimal combination is denoted with a rhombus, where μ_{3} = 1.0 × 10^{−3} and μ_{4} = 8.9 × 10^{−4}. 
Lvalues and other indexes before and after preprocessing.
Figure 5 shows field lines of different models. The original lineofsight magnetogram was used for the potential field modeling, the forcefree preprocessing was applied for the NLFFF modeling, and the MHS preprocessing was applied for the MHS modeling. All extrapolations were done in a 2336 × 1824 × 128 box. A central box with 936 × 936 × 128 dimensions above the IMaX FOV was cut to display the result. It is clear from Fig. 5 that field lines that share the same footpoints have different structures. A relatively large difference (e.g., the connectivity) can be found between panels a and b due to the forcefree currents in the NLFFF. The difference between panels b and c as a consequence of the perpendicular component of the current are much smaller. Field lines in both panels b and c have quite a similar pattern, but still have a different connectivity as well as length of field lines.
Fig. 5. Panel a: potential field reconstruction. Panel b: NLFF reconstruction. Panel c: MHS reconstruction. The field lines originate from the same footpoints in the negative polarity of each panels. 
In Fig. 6 we compare the simultaneous chromospheric observation by SUNRISE/SuFI 3968 Å (Gandorfer et al. 2011; Solanki et al. 2017) with field lines within subvolumes spanning the 600–1400 km height range. We note that the extrapolation data were cut according to the SuFI FOV. The observed slender fibrils are generally believed to outline the magnetic fields in this layer. We find that most field lines trace the fibrils nicely. However, deviations can also be observed in some regions.
Fig. 6. MHS field lines within the heights [600, 1400] km overplotted on the image observed in Ca II H core line with 1.1 Å wide filter. Footpoints of the field lines are uniformly selected on the photosphere. 
5. Conclusions
We have investigated the virial theorem for a MHS system. By applying it to a large active region, we find that the surface integral that was often used to compute the energy of the forcefree magnetic field is a lower bound of the magnetic energy in the semiinfinite volume. We also find a set of surface integrals as necessary criteria for a MHS system. These integrals are equivalent to Aly’s criteria for a forcefree field. According to the new set of criteria, we proposed an optimization algorithm to preprocess the vector magnetogram with the aim to use the result as a suitable input for a MHS extrapolation. The optimization strategy is similar to that designed in Wiegelmann et al. (2006) for the forcefree modeling, which is to force the data to fulfill the criteria for the MHS system and be sufficiently smooth within the freedom of the measurement noise. We also show an application of the preprocessing method to a combined vector magnetogram in which the IMaX data are embedded in the HMI data.
Acknowledgments
We thank the referee for helpful comments and suggestions. The German contribution to SUNRISE and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. This work was supported by DFGgrant WI 3211/41.
References
 Alissandrakis, C. E. 1981, A&A, 100, 197 [NASA ADS] [Google Scholar]
 Aly, J. J. 1984, ApJ, 283, 349 [NASA ADS] [CrossRef] [Google Scholar]
 Aly, J. J. 1989, Sol. Phys., 120, 19 [NASA ADS] [CrossRef] [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]
 Aschwanden, M. J. 2013, Sol. Phys., 287, 323 [NASA ADS] [CrossRef] [Google Scholar]
 Aschwanden, M. J. 2016, ApJS, 224, 25 [Google Scholar]
 Chandrasekhar, S. 1961, Hydrodynamic and Hydromagnetic Stability (Clarendon Press) [Google Scholar]
 Chiu, Y. T., & Hilton, H. H. 1977, ApJ, 212, 873 [NASA ADS] [CrossRef] [Google Scholar]
 Cuperman, S., Demoulin, P., & Semel, M. 1991, A&A, 245, 285 [NASA ADS] [Google Scholar]
 Demoulin, P., & Priest, E. R. 1992, A&A, 258, 535 [Google Scholar]
 Fan, Y. L., Wang, H. N., He, H., & Zhu, X. S. 2011, ApJ, 737, 39 [CrossRef] [Google Scholar]
 Fuhrmann, M., Seehafer, N., & Valori, G. 2007, A&A, 476, 349 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gandorfer, A., Grauf, B., Barthol, P., et al. 2011, Sol. Phys., 268, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Gary, G. A. 2001, Sol. Phys., 203, 71 [NASA ADS] [CrossRef] [Google Scholar]
 Gilchrist, S. A., & Wheatland, M. S. 2013, Sol. Phys., 282, 283 [NASA ADS] [CrossRef] [Google Scholar]
 Gilchrist, S. A., Braun, D. C., & Barnes, G. 2016, Sol. Phys., 291, 3583 [NASA ADS] [CrossRef] [Google Scholar]
 Grad, H., Rubbin, H., et al. 1958, in Theoretical and Experimental Aspects of Controlled Nuclear Fusion, eds. J. H. Martens, L. Ourom, W. M. Barss, et al., 31, 190 [Google Scholar]
 Guo, Y., Xia, C., Keppens, R., & Valori, G. 2016, ApJ, 828, 82 [NASA ADS] [CrossRef] [Google Scholar]
 He, H., & Wang, H. 2006, MNRAS, 369, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Hoeksema, J. T., Liu, Y., Hayashi, K., et al. 2014, Sol. Phys., 289, 3483 [Google Scholar]
 Jiang, C., & Feng, X. 2012, ApJ, 749, 135 [NASA ADS] [CrossRef] [Google Scholar]
 Jiang, C., & Feng, X. 2014, Sol. Phys., 289, 63 [CrossRef] [Google Scholar]
 Liu, S., & Hao, J. 2015, Adv. Space Res., 55, 1563 [CrossRef] [Google Scholar]
 Liu, S., Su, J. T., Zhang, H. Q., et al. 2013, PASA, 30, e005 [CrossRef] [Google Scholar]
 Malanushenko, A., Schrijver, C. J., DeRosa, M. L., Wheatland, M. S., & Gilchrist, S. A. 2012, ApJ, 756, 153 [NASA ADS] [CrossRef] [Google Scholar]
 Martínez Pillet, V., Del Toro Iniesta, J. C., ÁlvarezHerrero, A., et al. 2011, Sol. Phys., 268, 57 [NASA ADS] [CrossRef] [Google Scholar]
 Mikic, Z., Barnes, D. C., & Schnack, D. D. 1988, ApJ, 328, 830 [NASA ADS] [CrossRef] [Google Scholar]
 Miyoshi, T., Kusano, K., & Inoue, S. 2020, ApJS, 247, 6 [CrossRef] [Google Scholar]
 Molodenskii, M. M. 1969, Sov. Ast., 12, 585 [Google Scholar]
 Molodensky, M. M. 1974, Sol. Phys., 39, 393 [NASA ADS] [CrossRef] [Google Scholar]
 Moon, Y. J., Choe, G. S., Yun, H. S., Park, Y. D., & Mickey, D. L. 2002, ApJ, 568, 422 [NASA ADS] [CrossRef] [Google Scholar]
 Nakagawa, Y. 1974, ApJ, 190, 437 [CrossRef] [Google Scholar]
 Roumeliotis, G. 1996, ApJ, 473, 1095 [NASA ADS] [CrossRef] [Google Scholar]
 Sakurai, T. 1981, Sol. Phys., 69, 343 [NASA ADS] [CrossRef] [Google Scholar]
 Schatten, K. H., Wilcox, J. M., & Ness, N. F. 1969, Sol. Phys., 6, 442 [Google Scholar]
 Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, Sol. Phys., 275, 207 [Google Scholar]
 Schmidt, H. U. 1964, On the Observable Effects of Magnetic Energy Storage and Release Connected With Solar Flares, 50, 107 [Google Scholar]
 Seehafer, N. 1978, Sol. Phys., 58, 215 [NASA ADS] [CrossRef] [Google Scholar]
 Solanki, S. K., Riethmüller, T. L., Barthol, P., et al. 2017, ApJS, 229, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Song, M. T., Fang, C., Tang, Y. H., Wu, S. T., & Zhang, Y. A. 2006, ApJ, 649, 1084 [NASA ADS] [CrossRef] [Google Scholar]
 Tadesse, T., Wiegelmann, T., & Inhester, B. 2009, A&A, 508, 421 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tiwari, S. K. 2012, ApJ, 744, 65 [NASA ADS] [CrossRef] [Google Scholar]
 Valori, G., Kliem, B., & Keppens, R. 2005, A&A, 433, 335 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wheatland, M. S. 2004, Sol. Phys., 222, 247 [NASA ADS] [CrossRef] [Google Scholar]
 Wheatland, M. S. 2006, Sol. Phys., 238, 29 [NASA ADS] [CrossRef] [Google Scholar]
 Wheatland, M. S., Sturrock, P. A., & Roumeliotis, G. 2000, ApJ, 540, 1150 [Google Scholar]
 Wiegelmann, T. 2004, Sol. Phys., 219, 87 [NASA ADS] [CrossRef] [Google Scholar]
 Wiegelmann, T., & Inhester, B. 2010, A&A, 516, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wiegelmann, T., & Neukirch, T. 2006, A&A, 457, 1053 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wiegelmann, T., & Sakurai, T. 2012, Liv. Rev. Sol. Phys., 9, 5 [Google Scholar]
 Wiegelmann, T., Inhester, B., & Sakurai, T. 2006, Sol. Phys., 233, 215 [NASA ADS] [CrossRef] [Google Scholar]
 Wiegelmann, T., Neukirch, T., Ruan, P., & Inhester, B. 2007, A&A, 475, 701 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wiegelmann, T., Neukirch, T., Nickeler, D. H., et al. 2017, ApJS, 229, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Wu, S. T., Chang, H. M., & Hagyard, M. J. 1985, in Measurements of Solar Vector Magnetic Fields, ed. M. J. Hagyard, NASA Conf. Publ., 2374, 17 [Google Scholar]
 Wu, S. T., Sun, M. T., Chang, H. M., Hagyard, M. J., & Gary, G. A. 1990, ApJ, 362, 698 [NASA ADS] [CrossRef] [Google Scholar]
 Yamamoto, T. T., & Kusano, K. 2012, ApJ, 752, 126 [NASA ADS] [CrossRef] [Google Scholar]
 Yan, Y., & Li, Z. 2006, ApJ, 638, 1162 [NASA ADS] [CrossRef] [Google Scholar]
 Yan, Y., & Sakurai, T. 2000, Sol. Phys., 195, 89 [NASA ADS] [CrossRef] [Google Scholar]
 Zhu, X., & Wiegelmann, T. 2018, ApJ, 866, 130 [NASA ADS] [CrossRef] [Google Scholar]
 Zhu, X., & Wiegelmann, T. 2019, A&A, 631, A162 [Google Scholar]
 Zhu, X. S., Wang, H. N., Du, Z. L., & Fan, Y. L. 2013, ApJ, 768, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Zhu, X., Wang, H., Du, Z., & He, H. 2016, ApJ, 826, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Zhu, X., Wiegelmann, T., & Solanki, S. K. 2020, A&A, 640, A103 [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
All Figures
Fig. 1. Typical interaction of the magnetic force and plasma forces in an active region with single compact polarity, in a MHS state. 

In the text 
Fig. 2. Typical interaction of the magnetic force and plasma forces in an active region with double compact polarities, in a MHS state. 

In the text 
Fig. 3. Vector magnetogram of IMaX measured on 2013 June 12 at 23:48 UT. Top: original data. Bottom: preprocessed data. 

In the text 
Fig. 4. Test runs with various μ_{3} and μ_{4} combinations. The quantity L_{3} of each test reaches the noise level of the magnetogram. The optimal combination is denoted with a rhombus, where μ_{3} = 1.0 × 10^{−3} and μ_{4} = 8.9 × 10^{−4}. 

In the text 
Fig. 5. Panel a: potential field reconstruction. Panel b: NLFF reconstruction. Panel c: MHS reconstruction. The field lines originate from the same footpoints in the negative polarity of each panels. 

In the text 
Fig. 6. MHS field lines within the heights [600, 1400] km overplotted on the image observed in Ca II H core line with 1.1 Å wide filter. Footpoints of the field lines are uniformly selected on the photosphere. 

In the text 
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.