1. Introduction
Coronal holes (CHs) are extended dark patches on the solar surface observed in ultraviolet (UV) and X-ray radiation. The magnetic field in CHs is open, and the CHs are cooler and less dense than the corona. The CHs are supposed to have a low-
plasma (the ratio of gas to magnetic pressure) which is typically of the order of
. CHs have an open divergent magnetic structure and essentially lie above the unipolar magnetic regions of the photosphere. These structures have strong links with the solar wind which emanates from their base. See review [
1] about measurements of the plasma properties in CHs and how they are used to reveal details about the physical processes that heat the solar corona and accelerate the solar wind.
Physically different arrangements to CHs are active regions (ARs), composed of dense and hot plasma cores that are supposed to be distributed along closed magnetic field lines. The magnetic field in most of the ARs is bipolar, but more complex distributions can be formed due to the emergence of a new magnetic flux through the photosphere. UV-imaging spectroscopy during the early 1970s from Skylab already revealed that ARs are composed of filamentary structures, commonly called loops, rather than consisting of a simple diffuse plasma distribution. The AR corona has temperatures in excess of 3 MK and it is surrounded by a longer-lived halo with temperatures in the range 1.5 to 2.5 MK. The density typically changes from m−3 in the core to m−3 in the halo.
In spite of the dynamism revealed by the observations of coronal structures, it is possible to assume that they are in an approximate equilibrium if have long lifetimes. This leads to the construction of models based on the static assumption which under the presence of a magnetic field, is referred to as magnetohydrostatics (MHS). The MHS solutions are based on the force balance (between the magnetic, gas pressure and gravity) and also on the energy balance. In this regard, there is an apparent lack of literature on MHS equilibrium models of CHs, except the studies in Refs. [
2,
3,
4] of the 1980s. The MHS models can provide a better understanding of how CHs and ARs are kept in the two basic balance conditions. In addition, they can help to understand how CHs interact with other structures, such as active coronal regions. In this regard, the interaction between CHs and nearby ARs is a subject that was investigated mainly from the observational point of view, specially in relation to the solar wind generation (see, e.g. [
5]) and the reported outflows at the edge of ARs (see [
6]). In this context, Ref. [
7] studied, using magnetohydrodynamic (MHD) models, the interactions between open and closed magnetic field regions to understand the acceleration of particles. These 3D results indicate that an interchange reconnection allows accelerated particles to escape from deep within the coronal mass ejection (CME) flux rope.
The MHS models of CHs and ARs are required for other purposes. For example, to carry out investigations about the interaction of global MHD waves with these structures because, so far, relatively simple geometries were addressed, mostly based on a purely vertical magnetic field [
8,
9,
10,
11] in the case of CHs. Global 3D MHD simulations were also used to investigate this problem, e.g., [
12], but the results about the interaction with CHs are limited and a deeper analysis is needed, specially using elementary models like the ones proposed in the present paper.
Ref. [
13] explored different CH models that include a cold and underdense region (the CH) that connects with an atmosphere at typically 1 MK (corona) through a smooth CH boundary. These authors also investigated how simple MHS models can also reproduce the main features of ARs, mostly focusing on the high pressure and diffuse background of these structures instead of the single-loop structures often reported and analyzed in the observations; see, e.g., [
14,
15]. In Ref. [
13], CHs and ARs were considered as separate structures. The main goal of the present study is to build MHS models that include both CHs and ARs. Here, a common formalism is proposed to describe these beforehand antagonistic structures. The main characteristics of the obtained MHS models are analyzed in detail, paying particular attention to the parameters that can describe the interaction or interplay between CHs and ARs. Three basic types of MHS solutions are found and compared with the observations of the CHs with nearby ARs.
2. The Problem of Magnetohydrostatic Equilibrium in 2D
Let us start by introducing the basic equations that describe an MHS equilibrium. The force balance reads:
where
is the magnetic field,
p is the gas pressure,
the plasma density,
g the gravity acceleration on the solar surface, and
the magnetic permeability of free space. The magnetic field from Maxwell’s equations has to satisfy that
Let us suppose that the plasma is composed of a fully ionized hydrogen that satisfies the ideal gas law,
where
T is the temperature,
the gas constant and
the mean atomic weight (
for a fully ionized hydrogen plasma).
The aim here is to obtain solutions to Equations (
1)–(
3); however, these equations represent a system of five equations with six unknowns:
(three components),
p,
, and the temperature,
T. An energy equation is required to have a closed system, but here the approach of Ref. [
16] is adopted, where the energy equation is not solved directly (see also [
17,
18]). Instead, the temperature profile is chosen according to some observational constraints. Once a solution is obtained, one can calculate the corresponding energy balance that the system has to satisfy in order to keep a thermal equilibrium, but this is not the main goal of the present study; see Ref. [
13] for the calculation of the energy balance in some particular magnetic configurations.
The analysis here is restricted to two dimensions (2D) because the presence of gravity (pointing in the minus
z-direction) significantly complicates the calculation of MHS solutions in three dimensions. The magnetic field is written in terms of a magnetic flux function,
, ensuring Equation (
2), but that needs to be determined. The flux function is the
y-component of the vector potential. For simplicity, it is supposed here that there is no component of the magnetic field in the
-direction and write
It can be shown that the force balance condition in two dimensions given by Equation (
1) is equivalent to the following nonlinear elliptical partial differential equation in terms of the flux function
A (see, for example, [
16,
19,
20])
along with
Equation (
8) is straightforwardly solved using the ideal gas law,
where
is the temperature profile that can depend on the
z coordinate as well. Equation (
8) imposes a balance between the gas pressure gradient and the gravity force along the magnetic fields lines, while Equation (
7) represents the condition of force balance perpendicular to the magnetic field. The function
determines the profile of the gas pressure at
. Equation (
7) is of the Grad–Shafranov type but includes the effect of gravity. This equation must be solved under some boundary conditions. Here, a rectangular domain with
and
is considered, and
is set for the case with a symmetry axis. In this study, the pressure scale height, defined as
, where
is the reference coronal temperature, is used as the normalization spatial scale. Note that
h and
H are not necessarily equal, and
is taken in the calculations of
Section 5. The magnetic field and the thermodynamic variables are given in terms of
A,
and
. If the gas pressure is ignored, Equation (
7) reduces to a Laplace equation that leads to the potential solution.
The reference level used,
, is located at the base of the corona. If the magnetic field is known at this reference level, i.e.,
, then one can calculate the flux function at
by directly integrating Equation (
6):
which is used to find the gas pressure at the reference level:
One needs to solve Equation (
7) subject to the boundary value of
calculated above and in terms of some given
that has to be prescribed. Values at the rest of the boundaries of the rectangular box, i.e., at the left, right and top edges, are also needed. From the numerical point of view, the most straight condition to impose at these edges is that the magnetic field lines intersect perpendicularly with the boundaries; see, e.g., [
21]. This is straightforwardly achieved by imposing that the derivative of
A with respect to the perpendicular direction to the boundary is zero, according to Equations (
4) or (
6). With these boundary conditions (BCs), one has to make sure that the solution near the origin, or the region of interest, is not significantly affected by the values of
L and
H. In some special cases, it is possible to replace the last boundary condition by forcing the solution to be evanescent for large
z. Then, no finite height of the system is required. Normally, this condition can be applied when one solves the problem analytically; see [
13].
5. Numerical Results
The analytical approach considered in
Section 3 and
Section 4 is based on the values of the flux function, pressure and temperature at the lower boundary (
). However, one still needs to find the full solution in the plane
x-
z. It is crucial to understand how the topology of the magnetic field changes according to the values of the parameters and how the connection between the CH and the nearby bipolar region is established. In general, finding solutions to Equation (
7) requires numerical methods unless the profile for
and the boundary conditions are quite simple. Only in the case that Equation (
7) is a linear equation, general analytical methods are available. However, this is quite a complicated task at this point, and for this reason, Equation (
7) was solved numerically. For the details about the numerical method used, see [
13]. Here, it is sufficient to mention that the numerical code PDE2D [
22] was used which solves, via the method of finite elements, partial differential equations of the type found in the problem under consideration.
The numerical solutions for individual CHs and ARs were described in detail in Ref. [
13] and will be not repeated here. Let us focus on the most simple case of a configuration that combines these two structures and investigate the characteristics of the solutions in 2D. For simplicity reasons, a focus is on the symmetric case introduced earlier in
Section 4.3.1, and not on the results for the non-symmetric case discussed in
Section 4.3.2. Moreover, the analysis here is restricted to the basic model, although equivalent results were found for the extended model. It is convenient to remind that when the parameter
a is small, the extended model reduces to the basic model. There is also a technical reason as well to focus on the basic model; the convergence of the numerical method for the extended model is much slower than that for the basic model.
Figure 8 shows a numerical solution that contains both a CH (located at
) and an AR located on the right part of the spatial domain (
and
). In this example, the left foot of the AR has the same polarity as the CH (
) (case 1 in
Table 1). The temperature and density are low inside the CH, while these magnitudes are high in the AR in comparison with the coronal environment. Using the same functional form for the functions
and
, one obtains a rather smooth solution in the
x-
z plane. The maximum temperature and density values at the core of the AR agree with the derived analytical values given by Equations (
32) and (
35). For the present choice of parameters, this is the case pointing to the right of the reference value located at
in
Figure 2 and
Figure 3.
Figure 9 is similar to
Figure 8, but the left foot of the AR has the opposite polarity to that of the CH (
). It corresponds to a point to the left of the value
in
Figure 2 and
Figure 3 (case 2 in
Table 1). Again, the temperature and density are low inside the CH, but the differences arise essentially in the AR and near the corresponding footpoints. The temperature and density are high inside the AR, but next to the two footpoints (located at
and
), one finds a depletion in the temperature and density that remind of the behavior inside the CH. This is not surprising as soon as one realizes that the magnetic field lines are essentially open and expanding with a height around
and
in a way similar to that in the CH located at
. Another relevant feature of the MHS solution of the present example is that the temperature above the AR (
) does not tend to the coronal temperature as in
Figure 8 and it is closer to the CH temperature instead; therefore, the AR is essentially surrounded and embedded by a cool and low-density plasma.
The last example of the type of solutions allowed in the model here is shown in
Figure 10. The configuration is specifically chosen to have a magnetic field in the AR that is between the values given by Equations (
33) and (
34) (Case 3 in
Table 1), i.e., the difference with respect to the results shown in
Figure 8 and
Figure 9 is just the value of the magnetic field at the AR. For this case, one obtains a CH with basically the same properties as in
Figure 8 and
Figure 9, but the peculiarity relies on the bipole. The AR in this model is characterized by a low temperature and a low density in comparison with the coronal environment. Therefore, the closeness of the magnetic field lines does not necessarily indicate that the bipolar region must have a hot and dense core according to the solutions found. In this region, a void in the density and temperature is possible according to the model here. This configuration is the most peculiar, and most likely unrealistic, of the three examples shown in
Figure 8,
Figure 9 and
Figure 10.
6. Conclusions and Discussion
In this paper, MHS models that represent a compound of a CH and an AR are constructed using the formalism of Ref [
13]. 2D magnetostatic equilibria were generated by solving the Grad–Shafranov equation. The magnetic field arrangement, chosen to represent both open and closed magnetic field lines, is incorporated through the boundary conditions needed to solve the partial differential equation. A relatively simple functional form was used for the plasma pressure and temperature in terms of the flux function. The use of a common functional dependence for the thermal structure allows us to study the properties of an MHS equilibrium that contains a CH, taken as the reference structure, and an AR situated at a certain distance from the CH. This is the key assumption in this paper. Interestingly, the same functional form describes correctly the general features of both CHs (low pressure, temperature and density) and ARs (high pressure, temperature and density) at least in two of the three types of the models presented in this paper. The exact form of the pressure and temperature functions considered in the present study is not relevant. Other choices can also lead to physically acceptable solutions. The contribution of the model considered relies on the finding of solutions that can represent concurrently two different magnetic structures with properties similar to those reported by the observations. It is worth mentioning that the construction of models with different spatial dependencies according to the spatial location in the domain is still possible but much more complex than the situation explored in the present case. Hence, a different functional form for the CH and the AR could be considered. This needs to be addressed in future studies because it can lead to models that do not show unphysical or unobserved configurations like the ones reported here.
The study analyzed how the main parameters of one structure depend on the parameters of the other. In this regard, the novelty of the solutions obtained in the present paper is that they provide information about the coupling between the two structures, an effect that was not yet addressed in the literature, as to the author’s knowledge. This interplay between the two configurations is due to a) the common functional form assumed for the pressure and temperature and b) the topology of the magnetic field, i.e., how the global magnetic structure combines the joint effect of the presence of the CH and AR.
When the bipolar magnetic configuration is perfectly symmetric, i.e., has a net zero magnetic flux, keeping the parameters of the CH constant and using as a reference the mean coronal values for the temperature and pressure, three essentially different types of bipolar solutions or ARs are found. These solutions basically depend on the magnitude and especially on the sign of the magnetic field at the footpoints,
, see also
Table 1:
When the closest foot of the bipolar configuration to the CH and the CH itself have the same magnetic polarity, a solution resembling a typical AR with closed field lines embedding a hot and dense core is obtained.
When the magnetic field of the bipolar foot is below and therefore the CH and the closest AR foot have an opposite polarity, one again obtains a case of an AR with a hot and dense core. However, the magnetic field lines near the footpoints of the AR are essentially open and the plasma has low values of the pressure, temperature and density there. The AR is surrounded by cold and light plasma, instead. Hence, the configuration may resemble that of multiple CHs located near the AR or even a dark halo.
When the polarity of the CH and the closest foot of the bipolar region are opposite, but the magnetic field of the foot and is in the range from 0 to , one finds closed magnetic field lines filled with plasma with the pressure, temperature and density below the coronal reference values and reaching a minimum at CH values. The plasma conditions are therefore that of typical CHs but with closed magnetic field lines in nearby locations.
The results presented in this paper shed light on the possible relationship between CHs and nearby ARs and provide new information that should be carefully contrasted with observations. The models considered here predict the existence of three different types of bipolar regions interacting with CHs. It is apparent that the first type of solution, i.e., the standard AR configuration, is observed routinely. This is the most common case. The second type of configuration could be related to the existence of several CHs reported on the disk, and under some conditions, it could be related to the presence of outflows at the edge of the ARs (see [
5]) just because the magnetic field can be locally open there, as shown in the solutions found here. Actually, there are observations where the CH presumably extends to the edge of an AR or even surrounds it, but since the emission of the AR is too strong compared to the CH, these dark regions are difficult to be observed. The third type of solution is the most difficult to interpret because in general the observations do not show dark ARs. Nevertheless, there are cases where an AR disappears and a CH appears afterward precisely at the same location.
In such a configuration, the third type of solution could be possible but for a highly limited time before the AR transforms into a pure CH. Long observations are required to understand this case because the appearance of a CH at the same location of an AR is a process that takes several days/weeks, longer than the usual observing periods.
Simple analytical expressions are obtained here for the pressure, temperature and density at the core of the AR and their dependence were studied upon several major physical parameters, such as the intensity of the magnetic field and the width of the magnetic spatial distribution at the reference level . Two different models were investigated, the substantial difference among them being that in the extended model the temperature and pressure do not grow unbounded with the magnetic field of the AR, while in the basic model they do. In the extended model, it was imposed that the temperature and pressure can not be larger than a certain value and this leads to the redefinition of the dependence of the temperature and pressure on the flux function.
The situation for an AR with a certain imbalance in the magnetic flux was also considered, and this corresponds to the asymmetric case. This is a more realistic setting according to photospheric magnetograms. An interesting result under such asymmetry conditions is that the relevant magnitudes tend to asymptotic values when is significantly high. In this limit, the two models proposed here have restrained maximum pressure, temperature and density values at the core of the AR, and the exact magnetic field value of the AR in this situation is indeed not of much importance. Hence, an additional constrain anticipated by the model here is the existence of a maximum density at the core of the AR. This quantity is essentially independent of the magnetic field of the AR as long as is much larger than that at the center of the CH, . This maximum critical density depends only on the coronal reference values and on the plasma pressure and temperature at the center of the CH.
The family of solutions calculated numerically can be used, among other findings to study the process of the emergence of a magnetic flux adjacent to a CH. The emergence can involve a wide variety of time scales, but under some circumstances, the premise can be used that the magnetic field, which is coupled to the plasma in the model here considered, evolves through a series of equilibrium states (see, e.g., [
23,
24]). These equilibrium states are precisely the MHS solutions that were obtained from the calculations. This method, that can be considered complementary to the full MHD modeling, provides useful information about the topology of the magnetic field without the requirement of an accurate physical description of the reconnection process that needs to be implemented, often using various approximations, in the time-dependent MHD simulations. By smoothly increasing the magnetic field of the bipolar region (
) from zero and up to a certain value, the connectivity of the system has to change, and some field lines that were open become closed, and vice versa. In particular, open lines belonging to the right part of the CH become closed while initially closed lines rooted on the right edge of the bipolar region may transform into open. This is the essence of the process of an interchange reconnection (see [
7,
25,
26]). Interestingly, there is observational evidence of an emerging magnetic flux interacting with CHs supporting this mechanism ([
27,
28,
29]), and the theoretical models studied in the present paper can be useful.
It is imperative to mention that in the models proposed in this paper, several important processes that can have a relevant effect on the results were ignored. No presence of flows, ubiquitous in CHs as well as in ARs were considered. From the theoretical point of view, the inclusion of flows can be performed using the approaches of Refs. [
2,
3] (see also [
30]). In this case, the Grad–Shafranov equation is coupled to the Bernoulli equation. Nevertheless, fast flows of the order of the Alfvén velocity, which are typically below the observed ones (see [
1]), are needed to produce significant changes in the MHS equilibrium.
Another critical point of the approach used in this study is that the energy equation is not treated self-consistently. Among other difficulties, the exact form of the heating function in the corona is still unknown. For this reason, Ref. [
13] have tried to determine this heating function based on the constrains of the force balance and thermal balance. These authors find regions where the energy balance is not possible, and this depends on the scale of the temperature variation with a height that eventually modifies the conduction and the radiative losses in the system. In the present paper, it was preferred not to address this problem, while leaving it for future studies, and focussed on the force balance only. Including self-consistency and an energy equation can be performed using the approach of Ref. [
31].
The 2D models studied here might seem be too idealized to represent a real 3D structure, such as a CH or an AR. The 2D models are applicable when the structure is invariant in one direction (the
y-direction in this case). The use of the flux function allows us to simplify the problem. In 3D, the equilibrium of the forces require more complicated calculations and a numerical treatment is unavoidable. For example, using two Euler potentials, the equivalent counterparts of the flux function in 2D, one obtains two equations that provide the force balance across the magnetic field (see, e.g., [
32,
33]). This was recently investigated in Ref. [
34], where open and closed magnetic field lines in 3D are able to coexist under the presence of gravity. The procedure in Ref. [
34] can be used to investigate the features of the combination of a CH and an AR in 3D in an equivalent form, as was performed in 2D in the present paper.
Concluding, let us stress that an elementary understanding of the interplay between CHs and ARs using MHS models is essential before additional effects are added to the problem.