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.     

RENDERING THREE-DIMENSIONAL SOLAR CORONAL STRUCTURES

An X-ray or EUV image of the corona or chromosphere is a 2D representation of an extended 3D complex for which a general inversion process is impossible. A specific model must be incorporated in order to understand the full 3D structure. We approach this problem by modeling a set of optically-thin 3D plasma flux tubes which we render these as synthetic images. The resulting images allow the interpretation of the X-ray/EUV observations to obtain information on (1) the 3D structure of X-ray images, i.e., the geometric structure of the flux tubes, and on (2) the internal structure using specific plasma characteristics, i.e., the physical structure of the flux tubes. The data-analysis technique uses magnetograms to characterize photospheric magnetic fields and extrapolation techniques to form the field lines. Using a new set of software tools, we have generated 3D flux tube structures around these field lines and integrated the plasma emission along the line of sight to obtain a rendered image. A set of individual flux-tube images is selected by a non-negative least-squares technique to provide a match with an observed X-ray image. The scheme minimizes the squares of the differences between the synthesized image and the observed image with a non-negative constraint on the coefficients of the brightness of the individual flux-tube loops. The derived images are used to determine the specific photospheric foot points and physical data, i.e., scaling laws for densities and loop lengths. The development has led to computer efficient integration and display software that is compatible for comparison with observations (e.g., Yohkoh SXT data, NIXT, or EIT). This analysis is important in determining directly the magnetic field configuration, which provides the structure of coronal loops, and indirectly the electric currents or waves, which provide the energy for the heating of the plasma. We have used very simple assumptions (i.e., potential magnetic fields and isothermal corona) to provide an initial test of the techniques before complex models are introduced. We have separated the physical and geometric contributions of the emission for a set of flux tubes and concentrated, in this initial study, on the geometric contributions by making approximations to the physical contributions. The initial results are consistent with the scaling laws derived from the Yohkoh SXT data.     

On Analysis of Dual Spacecraft Stereoscopic Observations to Determine the Three-Dimensional Morphology and Plasma Properties of Solar Coronal Flux Tubes

By using two spacecraft equipped with multi-bandpass X-ray telescopes, it is possible to obtain direct 3-dimensional morphology of coronal structures which is essential for understanding the energetics and dynamics of the solar atmosphere. X-ray observations taken only in orbit about the Earth are inadequate to fully resolve the 3-dimensional nature of the solar corona. These Earth-orbit observations produce 2-dimensional images and an appropriate model must be included to derive the 3-dimensional structures from the line-of-sight information. Stereoscopic observations from space will remove this limitation and are needed if we are to improve our knowledge of the 3-dimensional morphology of the corona.Several important points regarding a stereoscopic mission are investigated and illustrated using model coronal flux tubes and image-rendering techniques. Synthesized images are formed by integrating the emission from volume elements along the line-of-sight path through a 3-dimensional volume in which a set of model flux tubes are located. The flux tubes are defined by (1) a plasma model defining the emissivity for a specific density, temperature, and pressure distribution, and (2) a magnetic field model from which a set of field lines are selected to define the geometry of the flux tubes. The field lines are used to define the flux-tube volume by assuming an initial base radius and conservation of flux. An effective instrumental spectral-response function is folded into the integration. Analysis of pairs of these synthesized images with various angular perspectives are used to investigate the effect of angular separation on mission objectives. The resulting images and analysis provide guidelines for developing a stereoscopic mission.Our study produced four important results, namely: (1) An angular separation of 30 degrees maximizes the scientific return by direct triangulation analysis because of the tradeoff between increased line-of-sight resolution of position and decreased recognition of individual loop structures arising from the overlapping of multiple loops with increasing angular separation. (2) The analysis benefits from the use of time-differential images to select flux tubes from the collection of numerous overlapping systems by selecting only recently heated or cooled flux tubes. (3) An analysis needs to be developed for algebraic reconstruction techniques applying a priori information, specific to the solar coronal structures, i.e., flux-tube continuity, maximum emission strength, non-negative emission, previous history, and maximum gradients of emission. (4) An analysis strategy combining triangulation, modeling techniques, and algebraic restoration is necessary to derive a complete understanding of the 3-dimensional morphology of the magnetic field. In the same way that helioseismology is classical viewing of the Sun with a tailored set of analysis tools for probing the interior of the Sun, heliostereoscopy is classical viewing of the X-ray emitting corona and requires a tailored set of analysis tools to deduce the true 3-dimensional structure of the corona.     

Magnetic Stereoscopy

The space mission Solar TErrestrial RElations Observatory (STEREO) will provide images from two viewpoints. An important aim of the STEREO mission is to get a 3D view of the solar corona. We develop a program for the stereoscopic reconstruction of 3D coronal loops from images taken with the two STEREO spacecraft. A pure geometric triangulation of coronal features leads to ambiguities because the dilute plasma emissions complicates the association of features in image 1 with features in image 2. As a consequence of these problems, the stereoscopic reconstruction is not unique and multiple solutions occur. We demonstrate how these ambiguities can be resolved with the help of different coronal magnetic field models (potential, linear, and non-linear force-free fields). The idea is that, due to the high conductivity in the coronal plasma, the emitting plasma outlines the magnetic field lines. Consequently, the 3D coronal magnetic field provides a proxy for the stereoscopy, which allows to eliminate inconsistent configurations. The combination of stereoscopy and magnetic modelling is more powerful than one of these tools alone. We test our method with the help of a model active region and plan to apply it to the solar case as soon as STEREO data become available.     

Connections Between the White-Light Eclipse Corona and Magnetic Fields over the Solar Cycle

Observations of ten solar eclipses (1973–1999) enabled us to reveal and describe mutual relations between the white-light corona structures (e.g., global coronal forms and most conspicuous coronal features, such as helmet streamers and coronal holes) and the coronal magnetic field strength and topology. The magnetic field strength and topology were extrapolated from the photospheric data under the current-free assumption. In spite of this simplification the found correspondence between the white-light corona structure and magnetic field organization strongly suggests a governing role of the field in the appearance and evolution of local and global structures. Our analysis shows that the study of white-light corona structures over a long period of time can provide valuable information on the magnetic field cyclic variations. This is particularly important for the epoch when the corresponding measurements of the photospheric magnetic field are absent.     

The Dependence of the EIT Wave Velocity on the Magnetic Field Strength

“EIT waves” are a wavelike phenomenon propagating in the corona, which was initially observed in the extreme ultraviolet (EUV) wavelength by the EUV Imaging Telescope (EIT). Their nature is still elusive, with the debate on-going between fast-mode wave model and non-wave model. In order to distinguish between these models, we investigate the relation between the EIT wave velocity and the local magnetic field in the corona. It is found that the two parameters show significant negative correlation in most of the EIT wave fronts, i.e., the EIT wave propagates more slowly in the regions of stronger magnetic field. Such a result poses a big challenge to the fast-mode wave model, which would predict a strong positive correlation between the two parameters. However, it is demonstrated that such a result can be explained by the fieldline stretching model, i.e., that “EIT waves” are the propagation of apparent brightenings, which are generated by successive stretching of closed magnetic field lines pushed by the erupting flux rope during coronal mass ejections (CMEs).     

Energy released by the interaction of coronal magnetic fields

Comparisons between coronal spectroheliograms and photospheric magnetograms are presented to support the idea that as coronal magnetic fields interact, a process of field line reconnection usually takes place as a natural way of preventing magnetic stresses from building up in the lower corona. This suggests that the energy which would have been stored in stressed fields is continuously released as kinetic energy of material being driven aside to make way for the reconnecting fields. However, this kinetic energy is negligible compared to the thermal energy of the coronal plasma. Therefore, it appears that these slow adjustments of coronal magnetic fields cannot account for even the normal heating of the corona, much less the energetic events associated with solar flares.Visiting Scientist, Kitt Peak National Observatory, Tucson Arizona.Operated by the Association of Universities for Research in Astronomy, Inc., under contract with 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.     

Coronal general magnetic field evolution as a new parameter of the solar cycle

The magnetic field lines of the corona associated with the solar-cycle surface general magnetic field are calculated by a potential-field approximation to study the solar-cycle evolution of the geometry of the coronal field. The surface field evolution used here is the radial field evolution, predicted by a model of the solar cycle driven by the dynamo action of the global convection, and justified observationally using Mount Wilson magnetic synoptic chart data. The evolution of the calculated coronal general field is now good for comparison with observations and shows the following. (i) The field of the polar and high-altitude corona has dipolar structure in almost all phases of the solar cycle except in a short time interval around maximum phase despite the quadrupolar structure of the general magnetic field at the surface; quadrupolar field forms loop-like structure in the lower corona. The almost-dipolar structure of the polar and high-latitude corona and the loop formation of the equatorial lower corona explain the appearance of the undisturbed minimum corona observed at eclipses. (ii) The polar field lines are directed almost radially at the minimum phase, which should be responsible for polar plumes. The field lines slowly open up to participate in the loop-like structure of the equatorial lower corona, and rapidly change their structure and polarity at the maximum phase, to resume the almost radial configuration slowly, (iii) During the rapidly changing maximum phase, the field lines do not penetrate deep into the interplanetary space resulting in the absence of polar plumes and the appearance of the circular corona- the maximum corona. In this phase, the coronal field should not be approximated by a dipole field. The surface field evolution which can explain such behaviors of the corona is characteristic of the solar-cycle process dominated by the latitudinal gradient of the differential rotation. If the radial gradient dominated in the subsurface process, the coronal evolution would look quite different and would show latitudinal propagation of enhancement of activity. Although nonaxisymmetric features should be superposed on the axisymmetric general field to express the real corona, the general field can be a basic coronal field in studying long-term interaction between the convection zone and the interstellar space especially in studying the magnetic braking of the solar rotation.     

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 .     

2D Feature Recognition And 3d Reconstruction In Solar Euv Images

EUV images show the solar corona in a typical temperature range of T >rsim 1 MK, which encompasses the most common coronal structures: loops, filaments, and other magnetic structures in active regions, the quiet Sun, and coronal holes. Quantitative analysis increasingly demands automated 2D feature recognition and 3D reconstruction, in order to localize, track, and monitor the evolution of such coronal structures. We discuss numerical tools that “fingerprint” curvi-linear 1D features (e.g., loops and filaments). We discuss existing finger-printing algorithms, such as the brightness-gradient method, the oriented-connectivity method, stereoscopic methods, time-differencing, and space–time feature recognition. We discuss improved 2D feature recognition and 3D reconstruction techniques that make use of additional a priori constraints, using guidance from magnetic field extrapolations, curvature radii constraints, and acceleration and velocity constraints in time-dependent image sequences. Applications of these algorithms aid the analysis of SOHO/EIT, TRACE, and STEREO/SECCHI data, such as disentangling, 3D reconstruction, and hydrodynamic modeling of coronal loops, postflare loops, filaments, prominences, and 3D reconstruction of the coronal magnetic field in general.     

MHD Modeling of Coronal Large-Amplitude Waves Related to CME Lift-off

We have employed a two-dimensional magnetohydrodynamic simulation code to study mass motions and large-amplitude coronal waves related to the lift-off of a coronal mass ejection (CME). The eruption of the filament is achieved by an artificial force acting on the plasma inside the flux rope. By varying the magnitude of this force, the reaction of the ambient corona to CMEs with different acceleration profiles can be studied. Our model of the ambient corona is gravitationally stratified with a quadrupolar magnetic field, resulting in an ambient Alfvén speed that increases as a function of height, as typically deduced for the low corona. The results of the simulations show that the erupting flux rope is surrounded by a shock front, which is strongest near the leading edge of the erupting mass, but also shows compression near the solar surface. For rapidly accelerating filaments, the shock front forms already in the low corona. Although the speed of the driver is less than the Alfvén speed near the top of the atmosphere, the shock survives in this region as well, but as a freely propagating wave. The leading edge of the shock becomes strong early enough to drive a metric type II burst in the corona. The speed of the weaker part of the shock front near the surface is lower, corresponding to the magnetosonic speed there. We analyze the (line-of-sight) emission measure of the corona during the simulation and recognize a wave receding from the eruption site, which strongly resembles EIT waves in the low corona. Behind the EIT wave, we clearly recognize a coronal dimming, also observed during CME lift-off. We point out that the morphology of the hot downstream region of the shock would be that of a hot erupting loop, so care has to be taken not to misinterpret soft X-ray imaging observations in this respect. Finally, the geometry of the magnetic field around the erupting mass is analyzed in terms of precipitation of particles accelerated in the eruption complex. Field lines connected to the shock are further away from the photospheric neutral line below the filament than the field lines connected to the current sheet below the flux rope. Thus, if the DC fields in the current sheet accelerate predominantly electrons and the shock accelerates ions, the geometry is consistent with recent observations of gamma rays being emitted further out from the neutral line than hard X-rays.     

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.     

Impulses of activity and the Solar cycle

Extreme-ultraviolet data from EIT/SOHO (1996–2002), soft X-ray data from Yohkoh (1991–2001), and magnetic field data from MDI/SOHO (1996–2002) and Kitt Peak Observatory, NSO/NOAO (1991–2002) are analyzed together in the form of synoptic maps for the investigation of solar cycle variations of the corona and their relation to the magnetic field. These results show new interesting relations between the evolution of the topological structure of the corona, coronal heating and the large-scale magnetic field. The long-lived coronal structures are related to complexes of solar activity and display quasi-periodic behavior (in the form of impulses of coronal activity) with periods of 1.0–1.5 year, in the axisymmetric distribution of EUV and X-ray fluxes during the current solar cycle 23. In particular, during the second maximum of this cycle the solar corona became somewhat hotter than it was in the period of the first maximum.     

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.     

Differential rotation and the structure and energy content of coronal magnetic fields

It is argued that differential rotation of the photospheric magnetic fields will induce currents in the corona. The work done against surface magnetic stresses will increase the energy content of the coronal magnetic field. The electrical conductivities are high and the foot points of field lines move with the differential rotation. The force-free field equations are solved with this constraint to obtain a minimum estimate of the energy increase for a quadrupole field. During a solar rotation the magnetic energy increases by 25%. Local release of this energy in the corona would have a significant effect. The expansion of field lines as a result of the differential rotation should increase the amount of flux and the field strength in the solar wind region.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Solar coronal heating by magnetosonic waves

Solar coronal heating by magnetohydrodynamic (MHD) waves is investigated. ultraviolet (UV) and X-ray emission lines of the corona show non-thermal broadenings. The wave rms velocities inferred from these observations are of the order of 25–60 km s⁻¹ . Assuming that these values are not negligible, we solved MHD equations in a quasi-linear approximation, by retaining the lowest order non-linear term in rms velocity. Plasma density distribution in the solar corona is assumed to be inhomogeneous. This plasma is also assumed to be permeated by dipole-like magnetic loops. Wave propagation is considered along the magnetic field lines. As dissipative processes, only the viscosity and parallel (to the local magnetic field lines) heat conduction are assumed to be important. Two wave modes emerged from the solution of the dispersion relation. The fast mode magneto-acoustic wave, if originated from the coronal base can propagate upwards into the corona and dissipate its mechanical energy as heat. The damping length-scale of the fast mode is of the order of 500 km. The wave energy flux associated with these waves turned out to be of the order of 2.5×10⁵ ergs cm⁻² s⁻¹ which is high enough to replace the energy lost by thermal conduction to the transition region and by optically thin coronal emission. The fast magneto-acoustic waves prove to be a likely candidate to heat the solar corona. The slow mode is absent, in other words cannot propagate in the solar corona.     

Solar-interplanetary modeling: 3-D solar wind solutions in prescribed non-radial magnetic field geometries

A model is presented which describes the 3-dimensional non-radial solar wind expansion between the Sun and the Earth in a specified magnetic field configuration subject to synoptically observed plasma properties at the coronal base. In this paper, the field is taken to be potential in the inner corona based upon the Mt. Wilson magnetograph observations and radial beyond a certain chosen surface. For plasma boundary conditions at the Sun, we use deconvoluted density profiles obtained from synopticK-coronameter brightness observations. The temperature is taken to be 2 × 10⁶ K at the base of closed field lines and 1.6 x 10⁶K at the base of open field lines. For a sample calculation, we employ data taken during the period of the 12 November 1966 eclipse. Although qualitative agreement with observations at 1 AU is obtained, important discrepancies emerge which are not apparent from spherically symmetric models or those models which do not incorporate actual observations in the lower corona. These discrepancies appear to be due to two primary difficulties - the rapid geometric divergence of the open field lines in the inner corona as well as the breakdown in the validity of the Spitzer heat conduction formula even closer to the Sun than predicted by radial flow models. These two effects combine to produce conductively dominated solutions and lower velocities, densities, and field strengths at the Earth than those observed. The traditional difficulty in solar wind theory in that unrealistically small densities must be assumed at the coronal base in order to obtain observed densities at 1 AU is more than compensated for here by the rapid divergence of field lines in the inner corona. For these base conditions, the value ofβ(ratio of gas pressure to magnetic pressure) is shown to be significantly greater than one over most of the lower corona - suggesting that, for the coronal boundary conditions used here, the use of a potential or force-free magnetic field configuration may not be justified. The calculations of this paper point to the directions where future research on solar-interplanetary modelling should receive priority:

1. better models for the coronal magnetic field structure
2. improved understanding of the thermal conductivity relevant for the solar wind plasma.