Modeling coronal magnetic field using spherical geometry: cases with several active regions

The magnetic fields in the solar atmosphere structure the plasma, store free magnetic energy and produce a wide variety of active solar phenomena, like flare and coronal mass ejections (CMEs). The distribution and strength of magnetic fields are routinely measured in the solar surface (photosphere). Therefore, there is considerable interest in accurately modeling the 3D structure of the coronal magnetic field using photospheric vector magnetograms. Knowledge of the 3D structure of magnetic field lines also help us to interpret other coronal observations, e.g., EUV images of the radiating coronal plasma. Nonlinear force-free field (NLFFF) models are thought to be viable tools for those task. Usually those models use Cartesian geometry. However, the spherical nature of the solar surface cannot be neglected when the field of view is large. In this work, we model the coronal magnetic field above multiple active regions using NLFFF extrapolation code using vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS)/Vector Spectromagnetograph (VSM) as a boundary conditions. We compare projections of the resulting magnetic field lines solutions with their respective coronal EUV-images from the Atmospheric Imaging Assembly (SDO/AIA) observed on October 15, 2011 and November 13, 2012. This study has found that the NLFFF model in spherical geometry reconstructs the magnetic configurations for several active regions which agrees to some extent with observations. During October 15, 2011 observation, there are substantial number of trans-equatorial loops carrying electric current.     

How To Use Magnetic Field Information For Coronal Loop Identification

The structure of the solar corona is dominated by the magnetic field because the magnetic pressure is about four orders of magnitude higher than the plasma pressure. Due to the high conductivity the emitting coronal plasma (visible, e.g., in SOHO/EIT) outlines the magnetic field lines. The gradient of the emitting plasma structures is significantly lower parallel to the magnetic field lines than in the perpendicular direction. Consequently information regarding the coronal magnetic field can be used for the interpretation of coronal plasma structures. We extrapolate the coronal magnetic field from photospheric magnetic field measurements into the corona. The extrapolation method depends on assumptions regarding coronal currents, e.g., potential fields (current-free) or force-free fields (current parallel to magnetic field). As a next step we project the reconstructed 3D magnetic field lines on an EIT-image and compare with the emitting plasma structures. Coronal loops are identified as closed magnetic field lines with a high emissivity in EIT and a small gradient of the emissivity along the magnetic field.     

On the formation of coronal cavities

We present a theoretical study of the formation of a coronal cavity and its relation to a quiescent prominence. We argue that the formation of a coronal cavity is initiated by the condensation of plasma which is trapped by the coronal magnetic field in a closed streamer and which then flows down to the chromosphere along the field lines due to lack of stable magnetic support against gravity. The existence of a coronal cavity depends on the coronal magnetic field strength; with low strength, the plasma density is not high enough for condensation to occur. Furthermore, we suggest that prominence and cavity material is supplied from the chromospheric level. Whether a coronal cavity and a prominence coexist depends on the magnetic field configuration; a prominence requires stable magnetic support.We initiate the study by considering the stability of condensation modes of a plasma in the coronal streamer model obtained by Steinolfson et al. (1982) using a 2-D, time dependent, ideal MHD computer simulation; they calculated the dynamic interaction between outward flowing solar wind plasma and a global coronal magnetic field. In the final steady state, they found a density enhancement in the closed field region with the enhancement increasing with increasing strength of the magnetic field. Our stability calculation shows that if the density enhancement is higher than a critical value, the plasma is unstable to condensation modes. We describe how, depending on the magnetic field configuration, the condensation may produce a coronal cavity and/or initiate the formation of a prominence.NRC Research Associate.     

Magnetic Discontinuities in the Near-Earth Solar Wind: Evidence of In-Transit Turbulence or Remnants of Coronal Structure?

Fluctuations in the solar wind plasma and magnetic field are well described by the sum of two power law distributions. It has been postulated that these distributions are the result of two independent processes: turbulence, which contributes mainly to the smaller fluctuations, and crossing the boundaries of flux tubes of coronal origin, which dominates the larger variations. In this study we explore the correspondence between changes in the magnetic field with changes in other solar wind properties. Changes in density and temperature may result from either turbulence or coronal structures, whereas changes in composition, such as the alpha-to-proton ratio are unlikely to arise from in-transit effects. Observations spanning the entire ACE dataset are compared with a null hypothesis of no correlation between magnetic field discontinuities and changes in other solar wind parameters. Evidence for coronal structuring is weaker than for in-transit turbulence, with only ∼ 25% of large magnetic field discontinuities associated with a significant change in the alpha-to-proton ratio, compared to ∼ 40% for significant density and temperature changes. However, note that a lack of detectable alpha-to-proton signature is not sufficient to discount a structure as having a solar origin.     

The solar corona as a minimum energy system?

This paper is an exploration of the possibility that the large-scale equilibrium of plasma and magnetic fields in the solar corona is a minimum energy state. Support for this conjecture is sought by considering the simplest form of that equilibrium in a dipole solar field, as suggested by the observed structure of the corona at times of minimum solar activity. Approximate, axisymmetric solutions to the MHD equations are constructed to include both a magnetically closed, hydrostatic region and a magnetically open region where plasma flows along field lines in the form of a transonic, thermally-driven wind. Sequences of such solutions are obtained for various degrees of magnetic field opening, and the total energy of each solution is computed, including contributions from both the plasma and magnetic field. It is shown that along a sequence of increasingly closed coronal magnetic field, the total energy curve is a non-monotonic function of the parameter measuring the degree of magnetic field opening, with a minimum occurring at moderate field opening.For reasonable choices of model parameters (coronal temperature, base density, base magnetic field strength, etc.), the morphology of the minimum energy solution resembles the observed quiet, solar minimum corona. The exact location energy minimum along a given sequence depends rather sensitively on some of the adopted parameter values. It is nevertheless argued that the existence of an energy minimum along the sequences of solutions should remain a robust property of more realistic coronal wind models that incorporate the basic characteristics of the equilibrium corona- the presence of both open and closed magnetic regions.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

The coronal and interplanetary magnetic fields at the time of the solar eclipse of 7 March, 1970

The magnetic field in the outer corona and in interplanetary space has been calculated from the photospheric magnetic fields measured around the time of the 7 March, 1970 eclipse. The field-line maps are compared with eclipse photographs showing coronal structures out to about 12 r . The projected field lines as well as the observed streamers appear straight. This is caused by the rapid expansion of the outer corona and is not an indication of corotation. The calculations show that the angular velocity of the coronal plasma decreases rapidly with distance.The relation between magnetic fields and density enhancements is discussed. The field strength in the photosphere seems to determine the amount of mechanical heating of the lower corona. The density structure higher up in the corona will, however, depend decisively on the topology of the field, particularly on whether we are on open or closed field lines, and not simply on field strength.The calculations show a sector structure of the interplanetary field, which agrees well with spacecraft observations. Also the magnitudes of the observed and calculated interplanetary field agree after the Mt. Wilson magnetograph data have been corrected to account for the temperature and saturation effects in the spectral line Fei 5250 Å.On leave from the Astronomical Observatory, Lund, Sweden.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Thermal instability in slab geometry in the presence of anisotropical thermal conduction

Prominences and filaments are thought to arise as a consequence of a magnetized plasma undergoing thermal instability. Therefore, the thermal stability of a magnetized plasma is investigated under coronal conditions. The equilibrium structure of the plasma is approximated by a 1-D slab configuration. This is investigated in thermal instability taking into account optically thin plasma radiation and anisotropic thermal conduction. The thermal conduction perpendicular to the magnetic field is taken to be small but non-zero.The classic rigid wall boundary conditions which are often applied in the literature, either directly on the plasma or indirectly through some other medium, are replaced by a more physical situation in which the plasma column is placed in a low-density background stretching towards infinity. Results for a uniform equilibrium structure indicate the major effect of this change is on the eigenfunctions rather than on the growth rate. Essentially, perpendicular thermal conduction introduces field-aligned fine structure. It is also shown that in the presence of perpendicular thermal conduction, thermal instability in a slab model is only possible if the inner plasma has the shortest thermal instability time scale.Research Assistant of the National Fund for Scientific Research (Belgium).     

Bidirectional Type III Solar Radio Bursts

Bidirectional coronal type III bursts are modeled by combining a model of coronal electron heating and beam generation via time-of-flight effects with semiquantitative estimates of quasilinear relaxation. Electromagnetic emissivities are estimated by extending the recently developed theory of interplanetary type III bursts to coronal emissions, including its features of stochastic Langmuir-wave growth and three-wave interactions. The results are investigated for heating on open and closed coronal field lines and are compared with observations of normal, reverse-slope, bidirectional, and inverted-J and -U coronal type III radio bursts. Harmonic emission is predicted to dominate at plasma frequencies above roughly 100 MHz where the efficiency of fundamental emission falls off steeply, while its free-free reabsorption rises. The model also explains the observed trends in the likelihood of occurrence of normal, reverse-slope, and bidirectional coronal type III bursts.     

Constructing the Coronal Magnetic Field By Correlating Parameterized Magnetic Field Lines With Observed Coronal Plasma Structures

A method is presented for constructing the coronal magnetic field from photospheric magnetograms and observed coronal loops. A set of magnetic field lines generated from magnetogram data is parameterized and then deformed by varying the parameterized values. The coronal flux tubes associated with this field are adjusted until the correlation between the field lines and the observed coronal loops is maximized. A mathematical formulation is described which ensures that (i) the normal component of the photospheric field remains unchanged, (ii) the field is given in the entire corona over an active region, (iii) the field remains divergence-free, and (iv) electric currents are introduced into the field. It is demonstrated that a parameterization of a potential field, comprising a radial stretching of the field, can provide a match for a simple bipolar active region, AR 7999, which crossed the central meridian on 1996 November 26. The result is a non-force-free magnetic field with the Lorentz force being of the order of 10^–5.5 g cm s^–2 resulting from an electric current density of 0.079 A m^–2. Calculations show that the plasma beta becomes larger than unity at a relatively low height of 0.25 r supporting the non-force-free conclusion. The presence of such strong non-radial currents requires large transverse pressure gradients to maintain a magnetostatic atmosphere, required by the relatively persistent nature of the coronal structures observed in AR 7999. This scheme is an important tool in generating a magnetic field solution consistent with the coronal flux tube observations and the observed photospheric magnetic field.     

Evolving coronal holes and interplanetary erupting stream disturbances

Coronal holes and interplanetary disturbances are important aspects of the physics of the Sun and heliosphere. Interplanetary disturbances are identified as an increase in the density turbulence compared with the ambient solar wind. Erupting stream disturbances are transient large-scale structures of enhanced density turbulence in the interplanetary medium driven by the high-speed flows of low-density plasma trailing behind for several days. Here, an attempt has been made to investigate the solar cause of erupting stream disturbances, mapped by Hewish & Bravo (1986) from interplanetary scintillation (IPS) measurements made between August 1978 and August 1979 at 81.5 MHz. The position of the sources of 68 erupting stream disturbances on the solar disk has been compared with the locations of newborn coronal holes and/or the areas that have been coronal holes previously. It is found that the occurrence of erupting stream disturbances is linked to the emergence of new coronal holes at the eruption site on the solar disk. A coronal hole is indicative of a radial magnetic field of a predominant magnetic polarity. The newborn coronal hole emerges on the Sun, owing to the changes in magnetic field configuration leading to the opening of closed magnetic structure into the corona. The fundamental activity for the onset of an erupting stream seems to be a transient opening of pre-existing closed magnetic structures into a new coronal hole, which can support highspeed flow trailing behind the compression zone of the erupting stream for several days.     

The structure of coronal arcades and the formation of solar prominences

The temperature and density are obtained for coronal plasma in thermal and hydrostatic equilibrium and located in a force-free magnetic arcade. The isotherms are found to be inclined to the magnetic field lines and so care should be taken in inferring the magnetic structure from observed emission.When the coronal pressure becomes too great, the equilibrium ceases to exist and the material cools to form a quiescent prominence. The same process can be initiated at low heating rates when the width or shear of the arcade exceeds a critical value.We suggest that the prominence should be modelled as a dynamic structure with plasma always draining downwards. Material is continually sucked up along field lines of the ambient arcade and into the region lacking a hot equilibrium, where it cools to form new prominence material.     

Cme Associated with Transequatorial Loops and a Bald Patch Flare

We study a flare which occurred on 3 November 1997 at 10:31 UT in the vicinity of a parasitic polarity of AR 8100. Using SOHO/EIT 195 Å observations, we identify the brightening of thin transequatorial loops connecting AR 8100 and AR 8102, and dimmings located between the two active regions. Difference images highlight the presence of a loop-like structure rooted near the flare location usually called an EIT wave. The coronal magnetic field derived from potential extrapolations from a SOHO/MDI magnetogram shows that the topology is complex near the parasitic polarity. There, a `bald patch' (defined as the locations where the magnetic field is tangent to the photosphere) is present. We conclude that the flare was a `bald patch flare'. Moreover, the extrapolation confirms that there is a large coronal volume filled with transequatorial field lines interconnecting AR 8100 and AR 8102, and overlaying the bald patch. We show that the dimmings are located at the footpoints of these large field lines, which can be also related to the thin bright loops observed during the flare. As this event was related to a coronal mass ejection (CME) observed by SOHO/LASCO, we propose that the observed dimmings are due to a decrease in plasma density during the opening of the transequatorial loops connecting both ARs. We propose a scenario where these large field lines are in fact pushed up by the opening of low-lying sheared field lines forming the bald patch. We finally discuss how the fast opening of these field lines can produce the brightening near the footpoints of the separatrix, observed as an `EIT wave'.     

Formation of neutral current sheet and loop coronal transients

An eruption of opposite magnetic flux into a bipolar background field is likely to lead to the formation of a natural current sheet between the new emerging field and the background. A numerical study is made on this process, based on the ideal MMD equations, taking into account the interaction between the magnetic field and the coronal plasma. The result shows that a subsonic eruption will give rise to a four region structure; 1) a cool and dense prominence made of the erupting material in the innermost region; 2) a cool and tenuous region further out; 3) a hot and dense loop formed by the concentration of both the erupting material and the coronal material in the neutral current sheet; and 4) a forerunner region outside the loop with density slightly above the background, due to fast magneto-acoustic waves. This structure agrees with the observed features of typical loop coronal transients. Therefore the eruption of opposite magnetic flux into a bipolar background is probably an important mechanism for triggering off such transients.     

Including stereoscopic information in the reconstruction of coronal magnetic fields

We present a method to include stereoscopic information about the three-dimensional structure of flux tubes into the reconstruction of the coronal magnetic field. Due to the low plasma beta in the corona we can assume a force-free magnetic field, with the current density parallel to the magnetic field lines. Here we use linear force-free fields for simplicity. The method uses the line-of-sight magnetic field on the photosphere as observational input. The value of is determined iteratively by comparing the reconstructed magnetic field with the observed structures. The final configuration is the optimal linear force-free solution constrained by both the photospheric magnetogram and the observed plasma structures. As an example we apply our method to SOHO MDI/EIT data of an active region. In the future it is planned to apply the method to analyse data from the SECCHI instrument aboard the STEREO mission.     

Flares and plasma flow caused by interacting coronal loops

Active region NOAA 7360 was observed in 1992 December with various instruments including the Yohkoh satellite. In this region, a small loop emerged near one of the footpoints of a pre-existing large coronal loop. These loops show evidence that interactions between coronal loops cause flares, microflares, and plasma flow. All of the four flares observed in this region show that brightenings in the small loop occurred first, and then the large loop flared up. The brightenings in the large loop can not occur by themselves, but must be triggered by the brightenings in the small loop. There must be interactions between the loops to cause these flares. As well as the flares, many microflares occurred in the small loop. More than half of them are accompanied by plasma ejection phenomena from the small loop into the large loop. The large loop is filled with ejected plasma with velocities of about 1000 km s^–1. These ejection phenomena are considered as X-ray jets. The associated occurrences of the microflares and the jets suggest that they are also caused by interactions between the loops. The recurrent occurrences of the homologous flares and microflares mean that the magnetic field structure in this region inevitably causes the activity due to loop-loop interactions; the flares and jets occur under a common magnetic field structure.     

Relating White-Light Coronal Images to Magnetic Fields and Plasma Flow

The solar magnetic field is key to a detailed understanding of the Sun's atmosphere and its transition to the solar wind. However, the lack of detailed magnetic field measurements everywhere except at the photosphere has made it challenging to determine its topology and to understand how it produces the observed plasma properties of the corona and solar wind. Recent progress based on the synthesis of diversified observations has shown that the corona is highly filamentary, that the coronal magnetic field is predominantly radial, and that the ability of closed fields to trap plasma at the base of the corona is a manifestation of how the solar field controls the solar wind. In this paper, we explain how these results are consistent with the relationship between density structure of white-light images and fields and flow. We point out that the ‘shape’ of the corona observed in white-light images is a consequence of the steep fall-off in density with radial distance, coupled with the inherent limitation in the sensitivity of the observing instrument. We discuss how the significant variation in radial density fall-off with latitude leads to a coronal shape that is more precisely revealed when a radial gradient filter is used, but which also gives a false impression of the tracing of highly non-radial fields. Instead, the coronal field is predominantly radial, and the two magnetic features that influence the shape of the corona are the closed fields at the base of the corona, and the polarity reversal forming the heliospheric current sheet in the outer corona. An erratum to this article is available at .     

Model for flare loops,fast motions,and opening of magnetic field in the corona

A model is presented for the penetration into the corona of a new magnetic field of a developing bipolar region and for its interaction with an old large-scale coronal field. An important feature of the model is a reconnection of the old and new fields inside the current sheet arising along the zero line of the total magnetic field calculated in the potential approximation. The magnetic reconnection and accumulation of plasma inside the current sheet can explain the appearance of dense coronal loops and the energy source at their tops. The plasma together with the magnetic lines is flowed into the sheet from both its sides. This fact explains the appearance of coronal cavities above the loops. If the large-scale field gradually decreases with the height, the loop motion is slowed down. The account of the dipolar structure of the magnetic field at large heights explains the possibility of a rapid break of the new field through the corona and the appearance of transients and open field regions - the coronal holes. In this case a fast rising current sheet can be a source of accelerated particles and of type II radio burst, instead of the shock wave considered usually.     

Attenuation of small-amplitude oscillations in a prominence–corona model with a transverse magnetic field

Observations show that small-amplitude prominence oscillations are usually damped after a few periods. This phenomenon has been theoretically investigated in terms of non-ideal magnetoacoustic waves, non-adiabatic effects being the best candidates to explain the damping in the case of slow modes. We study the attenuation of non-adiabatic magnetoacoustic waves in a slab prominence embedded in the coronal medium. We assume an equilibrium configuration with a transverse magnetic field to the slab axis and investigate wave damping by thermal conduction and radiative losses. The magnetohydrodynamic equations are considered in their linearised form and terms representing thermal conduction, radiation and heating are included in the energy equation. The differential equations that govern linear slow and fast modes are numerically solved to obtain the complex oscillatory frequency and the corresponding eigenfunctions. We find that coronal thermal conduction and radiative losses from the prominence plasma reveal as the most relevant damping mechanisms. Both mechanisms govern together the attenuation of hybrid modes, whereas prominence radiation is responsible for the damping of internal modes and coronal conduction essentially dominates the attenuation of external modes. In addition, the energy transfer between the prominence and the corona caused by thermal conduction has a noticeable effect on the wave stability, radiative losses from the prominence plasma being of paramount importance for the thermal stability of fast modes. We conclude that slow modes are efficiently damped, with damping times compatible with observations. On the contrary, fast modes are less attenuated by non-adiabatic effects and their damping times are several orders of magnitude larger than those observed. The presence of the corona causes a decrease of the damping times with respect to those of an isolated prominence slab, but its effect is still insufficient to obtain damping times of the order of the period in the case of fast modes.     

Waves and Oscillations in the Corona – (Invited Review)

It has long been suggested on theoretical grounds that MHD waves must occur in the solar corona, and have important implications for coronal physics. An unequivocal identification of such waves has however proved elusive, though a number of events were consistent with an interpretation in terms of MHD waves. Recent detailed observations of waves in events observed by SOHO and TRACE removes that uncertainty, and raises the importance of MHD waves in the corona to a higher level. Here we review theoretical aspects of how MHD waves and oscillations may occur in a coronal medium. Detailed observations of waves and oscillations in coronal loops, plumes and prominences make feasible the development of coronal seismology, whereby parameters of the coronal plasma (notably the Alfvén speed and through this the magnetic field strength) may be determined from properties of the oscillations. MHD fast waves are refracted by regions of low Alfvén speed and slow waves are closely field-guided, making regions of dense coronal plasma (such as coronal loops and plumes) natural wave guides for MHD waves. There are analogies with sound waves in ocean layers and with elastic waves in the Earth's crust. Recent observations also indicate that coronal oscillations are damped. We consider the various ways this may be brought about, and its implications for coronal heating.