Context. High resolution magnetic field measurements are routinely done only in the solar photosphere. Higher layers like the chromosphere and corona can be modeled by extrapolating the photospheric magnetic field upward. In the solar corona, plasma forces can be neglected and the Lorentz force vanishes. This is not the case in the upper photosphere and chromosphere where magnetic and non-magnetic forces are equally important. One way to deal with this problem is to compute the plasma and magnetic field self-consistently with a magnetohydrostatic (MHS) model. Aims. We aim to derive the magnetic field, plasma pressure and density of AR11768 by applying the newly developed extrapolation technique to the SUNRISE/IMaX data. Methods. An optimization method is used for the MHS modeling. The initial conditions consist of a nonlinear force-free field (NLFFF) and a gravity-stratified atmosphere. Results. In the non-force-free layer, which is spatially resolved by the new code, Lorentz forces are effectively balanced by the gas pressure gradient force and the gravity force. The pressure and density are depleted in strong field regions, which is consistent with observations. Denser plasma, however, is also observed at some parts of the active region edges. In the chromosphere, the fibril-like plasma structures trace the magnetic field nicely. Bright points in SUNRISE/SuFI 3000 {$\AA$} images are often accompanied by the plasma pressure and electric current concentrations. In addition, the average of angle between MHS field lines and the selected chromospheric fibrils is $11.8^\circ$, which is smaller than those computed from the NLFFF model ($15.7^\circ$) and linear MHS model ($20.9^\circ$). This indicates that the MHS solution provides a better representation of the magnetic field in the chromosphere.

For the extrapolation of magnetic fields into the solar corona from measurements taken in the photosphere (or chromosphere) force-free magnetic fields are typically used. This does not take into account that the lower layers of the solar atmosphere are not force-free. While some numerical extrapolation methods using magnetohydrostatic magnetic fields have been suggested, a complementary and numerically comparatively cheap method is to use analytical magnetohydrostatic equilibria to extrapolate the magnetic field. In this paper, we present a new family of solutions for a special class of analytical three-dimensional magnetohydrostatic equilibria, which can be of use for such magnetic field extrapolation. The new solutions allow for the more flexible modelling of a transition from non-force-free to (linear) force-free magnetic fields. In particular, the height and width of the region where this transition takes place can be specified by choosing appropriate model parameters.

Context. On the sun, the magnetic field vector is measured routinely only in the photosphere. By using these photospheric measurements as boundary condition, we developed the magnetohydrostatic (MHS) extrapolation to model the solar atmosphere. The model makes assumption about the relative importance of magnetic and non-magnetic forces. While the solar corona is force-free, this is not the case in photosphere and chromosphere. Aim. The model has been tested with an exact equilibria in \cite{zw18}. Here we present a more challenging and realistic test of our model with radiative MHD simulation of a solar flare. Methods. By using the optimization method, the MHS model computes self-consistently the magnetic field, plasma pressure and density. The nonlinear force-free field (NLFFF) and gravity stratified atmosphere along the field line are assumed as the initial condition of the optimization. Results. Compared with NLFFF, the MHS model gives an improved magnetic field not only in magnitude and direction, but also in the magnetic connectivity. Besides, the MHS model is able to recover the main structure of the plasma in the photosphere and chromosphere.

In a thorough study, we investigate the origin of a remarkable plasma and magnetic field configuration observed in situ on June 22, 2011 near L1, which appears to be a magnetic ejecta (ME) and a shock signature engulfed by a solar wind high-speed stream (HSS). We identify the signatures as an Earth-directed coronal mass ejection (CME), associated with a C7.7 flare on June 21, 2011, and its interaction with a HSS, which emanates from a coronal hole (CH) close to the launch site of the CME. The results indicate that the major interaction between the CME and the HSS starts at a height of 1.3 Rsun up to 3 Rsun. Over that distance range, the CME undergoes a strong north-eastward deflection of at least 30 degrees due to the open magnetic field configuration of the CH. We perform a comprehensive analysis for the CME-HSS event using multi-viewpoint data (from the Solar TErrestrial RElations Observatories, the Solar and Heliospheric Observatory and the Solar Dynamics Observatory), and combined modeling efforts (nonlinear force-free field modeling, Graduated Cylindrical Shell CME modeling, and the Forecasting a CMEs Altered Trajectory ForeCAT model). We aim at better understanding its early evolution and interaction process as well as its interplanetary propagation and related in situ signatures, and finally the resulting impact on the Earth's magnetosphere.

Understanding the relationship between the magnetic field and coronal heating is one of the central problems of solar physics. However, studies of the magnetic properties of impulsively heated loops have been rare. We present results from a study of 34 evolving coronal loops observed in the Fe XVIII line component of AIA/SDO 94 A filter images from three active regions with different magnetic conditions. We show that the peak intensity per unit cross-section of the loops depends on their individual magnetic and geometric properties. The intensity scales proportionally to the average field strength along the loop ($B_{avg}$) and inversely with the loop length ($L$) for a combined dependence of $(B_{avg}/L)^{0.52\pm0.13}$. These loop properties are inferred from magnetic extrapolations of the photospheric HMI/SDO line-of-sight and vector magnetic field in three approximations: potential and two Non Linear Force-Free (NLFF) methods. Through hydrodynamic modeling (EBTEL model) we show that this behavior is compatible with impulsively heated loops with a volumetric heating rate that scales as $\epsilon_H\sim B_{avg}^{0.3\pm0.2}/L^{0.2\pm^{0.2}_{0.1}}$.

Modeling the interface region between solar photosphere and corona is challenging, because the relative importance of magnetic and plasma forces change by several orders of magnitude. While the solar corona can be modeled by the force-free assumption, we need to take care about plasma forces (pressure gradient and gravity) in photosphere and chromosphere, here within the magneto-hydro-static (MHS) model. We solve the MHS equations with the help of an optimization principle and use vector magnetogram as boundary condition. Positive pressure and density are ensured by replacing them with two new basic variables. The Lorentz force during optimization is used to update the plasma pressure on the bottom boundary, which makes the new extrapolation works even without pressure measurement on the photosphere. Our code is tested by using a linear MHS model as reference. From the detailed analyses, we find that the newly developed MHS extrapolation recovers the reference model at high accuracy. The MHS extrapolation is, however, numerically more expensive than the nonlinear force-free field (NLFFF) extrapolation and consequently one should limit their application to regions where plasma forces become important, e.g. in a layer of about 2 Mm above the photosphere.

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

A wide variety of phenomena such as gentle but persistent brightening, dynamic slender features (~100 km), and compact (~1'') ultraviolet (UV) bursts are associated with the heating of the solar chromosphere. High spatio-temporal resolution is required to capture the finer details of the likely magnetic reconnection-driven, rapidly evolving bursts. Such observations are also needed to reveal their similarities to large-scale flares, which are also thought to be reconnection driven, and more generally their role in chromospheric heating. Here we report observations of chromospheric heating in the form of a UV burst obtained with the balloon-borne observatory, SUNRISE. The observed burst displayed a spatial morphology similar to that of a large-scale solar flare with circular ribbon. While the co-temporal UV observations at 1.5'' spatial resolution and 24s cadence from the Solar Dynamics Observatory showed a compact brightening, the SUNRISE observations at diffraction-limited spatial resolution of 0.1'' at 7s cadence revealed a dynamic sub-structure of the burst that it is composed of extended ribbon-like features and a rapidly evolving arcade of thin (~0.1'' wide) magnetic loop-like features, similar to post-flare loops. Such a dynamic sub-structure reveals the small-scale nature of chromospheric heating in these bursts. Furthermore, based on magnetic field extrapolations, this heating event is associated with a complex fan-spine magnetic topology. Our observations strongly hint at a unified picture of magnetic heating in the solar atmosphere from some large-scale flares to small-scale bursts, all being associated with such a magnetic topology.

It is widely believed that loops observed in the solar atmosphere trace out magnetic field lines. However, the degree to which magnetic field extrapolations yield field lines that actually do follow loops has yet to be studied systematically. In this paper we apply three different extrapolation techniques - a simple potential model, a NLFF model based on photospheric vector data, and a NLFF model based on forward fitting magnetic sources with vertical currents - to 15 active regions that span a wide range of magnetic conditions. We use a distance metric to assess how well each of these models is able to match field lines to the 12,202 loops traced in coronal images. These distances are typically 1-2". We also compute the misalignment angle between each traced loop and the local magnetic field vector, and find values of 5-12$^\circ$. We find that the NLFF models generally outperform the potential extrapolation on these metrics, although the differences between the different extrapolations are relatively small. The methodology that we employ for this study suggests a number of ways that both the extrapolations and loop identification can be improved.

We report on high-resolution imaging and spectral observations of eruptions of a spiral structure in the transition region, which were taken with the Interface Region Imaging Spectrometer (IRIS), the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI). The eruption coincided with the appearance of two series of jets, with velocities comparable to the Alfv\'en speeds in their footpoints. Several pieces of evidence of magnetic braiding in the eruption are revealed, including localized bright knots, multiple well-separated jet threads, transition region explosive events and the fact that all these three are falling into the same locations within the eruptive structures. Through analysis of the extrapolated three-dimensional magnetic field in the region, we found that the eruptive spiral structure corresponded well to locations of twisted magnetic flux tubes with varying curl values along their lengths. The eruption occurred where strong parallel currents, high squashing factors, and large twist numbers were obtained. The electron number density of the eruptive structure is found to be $\sim3\times10^{12}$ cm$^{-3}$, indicating that significant amount of mass could be pumped into the corona by the jets. Following the eruption, the extrapolations revealed a set of seemingly relaxed loops, which were visible in the AIA 94 \AA\ channel indicating temperatures of around 6.3 MK. With these observations, we suggest that magnetic braiding could be part of the mechanisms explaining the formation of solar eruption and the mass and energy supplement to the corona.

Getting insights into the 3D structure of the solar coronal magnetic field have been done in the past by two completely different approaches: (1.) Nonlinear force-free field (NLFFF) extrapolations, which use photospheric vector magnetograms as boundary condition. (2.) Stereoscopy of coronal magnetic loops observed in EUV coronal images from different vantage points. Both approaches have their strength and weaknesses. Extrapolation methods are sensitive to noise and inconsistencies in the boundary data and the accuracy of stereoscopy is affected by the ability of identifying the same structure in different images and by the separation angle between the view directions. As a consequence, for the same observational data, the computed 3D coronal magnetic field with the two methods do not necessarily coincide. In an earlier work (Paper I) we extended our NLFFF optimization code by the inclusion of stereoscopic constrains. The method was successfully tested with synthetic data and within this work we apply the newly developed code to a combined data-set from SDO/HMI, SDO/AIA and the two STEREO spacecraft. The extended method (called S-NLFFF) contains an additional term that monitors and minimizes the angle between the local magnetic field direction and the orientation of the 3D coronal loops reconstructed by stereoscopy. We find that prescribing the shape of the 3D stereoscopically reconstructed loops the S-NLFFF method leads to a much better agreement between the modeled field and the stereoscopically reconstructed loops. We also find an appreciable decrease by a factor of two in the angle between the current and the magnetic field which indicates the improved quality of the force-free solution obtained by S-NLFFF.

{We report on the first of a series of studies on coronal bright points investigating the physical mechanism that generates these phenomena.} {The aim of this paper is to understand the magnetic-field structure that hosts the bright points.} {We use longitudinal magnetograms taken by the Solar Optical Telescope with the Narrowband Filter Imager. For a single case, magnetograms from the Helioseismic and Magnetic Imager were added to the analysis. The longitudinal magnetic field component is used to derive the potential magnetic fields of the large regions around the bright points. A magneto-static field extrapolation method is tested to verify the accuracy of the potential field modelling. The three dimensional magnetic fields are investigated for the presence of magnetic null points and their influence on the local magnetic domain.} {In 9 out of 10 cases the bright point resides in areas where the coronal magnetic field contains an opposite polarity intrusion defining a magnetic null point above it. It is found that X-ray bright points reside, in these 9 cases, in a limited part of the projected fan dome area, either fully inside the dome or expanding over a limited area below which typically a dominant flux concentration resides. The 10th bright point is located in a bipolar loop system without an overlying null point.} {All bright points in coronal holes and two out of tree bright points in quiet Sun regions are seen to reside in regions containing a magnetic null point. An yet unidentified process(es) generates the BPs in specific regions of the fan-dome structure. }

How the solar corona is heated to high temperatures remains an unsolved mystery in solar physics. In the present study we analyse observations of 50 whole active-region loops taken with the Extreme-ultraviolet Imaging Spectrometer (EIS) on board the Hinode satellite. Eleven loops were classified as cool (<1 MK) and 39 as warm (1-2 MK) loops. We study their plasma parameters such as densities, temperatures, filling factors, non-thermal velocities and Doppler velocities. We combine spectroscopic analysis with linear force-free magnetic-field extrapolation to derive the three-dimensional structure and positioning of the loops, their lengths and heights as well as the magnetic field strength along the loops. We use density-sensitive line pairs from Fe XII, Fe XIII, Si X and Mg VII ions to obtain electron densities by taking special care of intensity background-subtraction. The emission-measure loci method is used to obtain the loop temperatures. We find that the loops are nearly isothermal along the line-of-sight. Their filling factors are between 8% and 89%. We also compare the observed parameters with the theoretical RTV scaling law. We find that most of the loops are in an overpressure state relative to the RTV predictions. In a followup study, we will report a heating model of a parallel-cascade-based mechanism and will compare the model parameters with the loop plasma and structural parameters derived here.

A dense forest of slender bright fibrils near a small solar active region is seen in high-quality narrowband Ca II H images from the SuFI instrument onboard the Sunrise balloon-borne solar observatory. The orientation of these slender Ca II H fibrils (SCF) overlaps with the magnetic field configuration in the low solar chromosphere derived by magnetostatic extrapolation of the photospheric field observed with Sunrise/IMaX and SDO/HMI. In addition, many observed SCFs are qualitatively aligned with small-scale loops computed from a novel inversion approach based on best-fit numerical MHD simulation. Such loops are organized in canopy-like arches over quiet areas that differ in height depending on the field strength near their roots.

Many magnetic structures in the solar atmosphere evolve rather slowly so that they can be assumed as (quasi-)static or (quasi-)stationary and represented via magneto-hydrostatic (MHS) or stationary magneto-hydrodynamic (MHD) equilibria, respectively. While exact 3D solutions would be desired, they are extremely difficult to find in stationary MHD. We construct solutions with magnetic and flow vector fields that have three components depending on all three coordinates. We show that the non-canonical transformation method produces quasi-3D solutions of stationary MHD by mapping 2D or 2.5D MHS equilibria to corresponding stationary MHD states, i.e., states that display the same field line structure as the original MHS equilibria. These stationary MHD states exist on magnetic flux surfaces of the original 2D MHS states. Although the flux surfaces and therefore also the equilibria have a 2D character, these stationary MHD states depend on all three coordinates and display highly complex currents. The existence of geometrically complex 3D currents within symmetric field-line structures provide the base for efficient dissipation of the magnetic energy in the solar corona by Ohmic heating. We also discuss the possibility of maintaining an important subset of non-linear MHS states, namely force-free fields, by stationary flows. We find that force-free fields with non-linear flows only arise under severe restrictions of the field-line geometry and of the magnetic flux density distribution.

Magneto-static models may overcome some of the issues facing force-free magnetic field extrapolations. So far they have seen limited use and have faced problems when applied to quiet-Sun data. Here we present a first application to an active region. We use solar vector magnetic field measurements gathered by the IMaX polarimeter during the flight of the \sunrise{} balloon-borne solar observatory in June 2013 as boundary condition for a magneto-static model of the higher solar atmosphere above an active region. The IMaX data are embedded in active region vector magnetograms observed with SDO/HMI. This work continues our magneto-static extrapolation approach, which has been applied earlier ({\it Paper I}) to a quiet Sun region observed with \sunrise{} I. In an active region the signal-to-noise-ratio in the measured Stokes parameters is considerably higher than in the quiet Sun and consequently the IMaX measurements of the horizontal photospheric magnetic field allow us to specify the free parameters of the model in a special class of linear magneto-static equilibria. The high spatial resolution of IMaX (110-130 km, pixel size 40 km) enables us to model the non-force-free layer between the photosphere and the mid chromosphere vertically by about 50 grid points. In our approach we can incorporate some aspects of the mixed beta layer of photosphere and chromosphere, e.g., taking a finite Lorentz force into account, which was not possible with lower resolution photospheric measurements in the past. The linear model does not, however, permit to model intrinsic nonlinear structures like strongly localized electric currents.

The prediction of the arrival time for fast coronal mass ejections (CMEs) and their associated shocks is highly desirable in space weather studies. In this paper, we use two shock propagation models, i.e. Data Guided Shock Time Of Arrival (DGSTOA) and Data Guided Shock Propagation Model (DGSPM), to predict the kinematical evolution of interplanetary shocks associated with fast CMEs. DGSTOA is based on the similarity theory of shock waves in the solar wind reference frame, and DGSPM on the non-similarity theory in the stationary reference frame. The inputs are the kinematics of the CME front at the maximum speed moment obtained from the geometric triangulation method applied to STEREO imaging observations together with the Harmonic Mean approximation. The outputs provide the subsequent propagation of the associated shock. We apply these models to the CMEs on 2012 January 19, January 23, and March 7. We find that the shock models predict reasonably well the shock's propagation after the impulsive acceleration. The shock's arrival time and local propagation speed at Earth predicted by these models are consistent with in situ measurements of WIND. We also employ the Drag-Based Model (DBM) as a comparison, and find that it predicts a steeper deceleration than the shock models after the rapid deceleration phase. The predictions of DBM at 1 AU agree with the following ICME or sheath structure, not the preceding shock. These results demonstrate the applicability of the shock models used here for future arrival time prediction of interplanetary shocks associated with fast CMEs.

We report the evolution of magnetic field and its energy in NOAA active region 11158 over 5 days based on a vector magnetogram series from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamic Observatory (SDO). Fast flux emergence and strong shearing motion led to a quadrupolar sunspot complex that produced several major eruptions, including the first X-class flare of Solar Cycle 24. Extrapolated non-linear force-free coronal fields show substantial electric current and free energy increase during early flux emergence near a low-lying sigmoidal filament with sheared kilogauss field in the filament channel. The computed magnetic free energy reaches a maximum of ~2.6e32 erg, about 50% of which is stored below 6 Mm. It decreases by ~0.3e32 erg within 1 hour of the X-class flare, which is likely an underestimation of the actual energy loss. During the flare, the photospheric field changed rapidly: horizontal field was enhanced by 28% in the core region, becoming more inclined and more parallel to the polarity inversion line. Such change is consistent with the conjectured coronal field "implosion", and is supported by the coronal loop retraction observed by the Atmospheric Imaging Assembly (AIA). The extrapolated field becomes more "compact" after the flare, with shorter loops in the core region, probably because of reconnection. The coronal field becomes slightly more sheared in the lowest layer, relaxes faster with height, and is overall less energetic.

We investigate the characteristics and the sources of the slow (< 450 km/s) solar wind during the four years (2006-2009) of low solar activity between Solar Cycles 23 and 24. We use a comprehensive set of in-situ observations in the near-Earth solar wind (Wind and ACE) and remove the periods when large-scale interplanetary coronal mass ejections were present. The investigated period features significant variations in the global coronal structure, including the frequent presence of low-latitude active regions in 2006-2007, long-lived low- and mid-latitude coronal holes in 2006 - mid-2008 and mostly the quiet Sun in 2009. We examine both Carrington Rotation averages of selected solar plasma, charge state and compositional parameters and distributions of these parameters related to Quiet Sun, Active Region Sun and the Coronal Hole Sun. While some of the investigated parameters (e.g., speed, the C^{+6}/C^{+4} and He/H ratio) show clear variations over our study period and with solar wind source type, some (Fe/O) exhibit very little changes. Our results highlight the difficulty in distinguishing between the slow solar wind sources based on the inspection of the solar wind conditions.

We investigate the evolution of NOAA Active Region 11817 during 2013 August 10--12, when it developed a complex field configuration and produced four confined, followed by two eruptive, flares. These C-and-above flares are all associated with a magnetic flux rope (MFR) located along the major polarity inversion line, where shearing and converging photospheric flows are present. Aided by the nonlinear force-free field modeling, we identify the MFR through mapping magnetic connectivities and computing the twist number $\mathcal{T}_w$ for each individual field line. The MFR is moderately twisted ($|\mathcal{T}_w| < 2$) and has a well-defined boundary of high squashing factor $Q$. We found that the field line with the extremum $|\mathcal{T}_w|$ is a reliable proxy of the rope axis, and that the MFR's peak $|\mathcal{T}_w|$ temporarily increases within half an hour before each flare while it decreases after the flare peak for both confined and eruptive flares. This pre-flare increase in $|\mathcal{T}_w|$ has little effect on the active region's free magnetic energy or any other parameters derived for the whole region, due to its moderate amount and the MFR's relatively small volume, while its decrease after flares is clearly associated with the stepwise decrease in free magnetic energy due to the flare. We suggest that $\mathcal{T}_w$ may serve as a useful parameter in forewarning the onset of eruption, and therefore, the consequent space weather effects. The helical kink instability is identified as the prime candidate onset mechanism for the considered flares.

Our aim is to model the 3D magnetic field structure of the upper solar atmosphere, including regions of non-negligible plasma beta. We use high-resolution photospheric magnetic field measurements from SUNRISE/IMaX as boundary condition for a magneto-static magnetic field model. The high resolution of IMaX allows us to resolve the interface region between photosphere and corona, but modelling this region is challenging for the following reasons. While the coronal magnetic field is thought to be force-free (the Lorentz-force vanishes), this is not the case in the mixed plasma $\beta$ environment in the photosphere and lower chromosphere. In our model, pressure gradients and gravity forces are taken self-consistently into account and compensate the non-vanishing Lorentz-force. Above a certain height (about 2 Mm) the non-magnetic forces become very weak and consequently the magnetic field becomes almost force-free. Here we apply a linear approach, where the electric current density consists of a superposition of a field-line parallel current and a current perpendicular to the Sun's gravity field. We illustrate the prospects and limitations of this approach and give an outlook for an extension towards a non-linear model.

The dynamic properties of flare ribbons and the often associated filament eruptions can provide crucial information on the flaring coronal magnetic field. This Letter analyzes the GOES-class X1.0 flare on 2014 March 29 (SOL2014-03-29T17:48), in which we found an asymmetric eruption of a sigmoidal filament and an ensuing circular flare ribbon. Initially both EUV images and a preflare nonlinear force-free field model show that the filament is embedded in magnetic fields with a fan-spine-like structure. In the first phase, which is defined by a weak but still increasing X-ray emission, the western portion of the sigmoidal filament arches upward and then remains quasi-static for about five minutes. The western fan-like and the outer spine-like fields display an ascending motion, and several associated ribbons begin to brighten. Also found is a bright EUV flow that streams down along the eastern fan-like field. In the second phase that includes the main peak of hard X-ray (HXR) emission, the filament erupts, leaving behind two major HXR sources formed around its central dip portion and a circular ribbon brightened sequentially. The expanding western fan-like field interacts intensively with the outer spine-like field, as clearly seen in running difference EUV images. We discuss these observations in favor of a scenario where the asymmetric eruption of the sigmoidal filament is initiated due to an MHD instability and further facilitated by reconnection at a quasi-null in corona; the latter is in turn enhanced by the filament eruption and subsequently produces the circular flare ribbon.

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.

In this paper, we report our multiwavelength observations of a partial filament eruption event in NOAA active region 11283 on 2011 September 8. A magnetic null point and the corresponding spine and separatrix surface are found in the active region. Beneath the null point, a sheared arcade supports the filament along the highly complex and fragmented polarity inversion line. After being activated, the sigmoidal filament erupted and split into two parts. The major part rose at the speeds of 90$-$150 km s$^{-1}$ before reaching the maximum apparent height of $\sim$115 Mm. Afterwards, it returned to the solar surface in a bumpy way at the speeds of 20$-$80 km s$^{-1}$. The rising and falling motions were clearly observed in the extreme-ultravoilet (EUV), UV, and H$\alpha$ wavelengths. The failed eruption of the main part was associated with an M6.7 flare with a single hard X-ray source. The runaway part of the filament, however, separated from and rotated around the major part for $\sim$1 turn at the eastern leg before escaping from the corona, probably along large-scale open magnetic field lines. The ejection of the runaway part resulted in a very faint coronal mass ejection (CME) that propagated at an apparent speed of 214 km s$^{-1}$ in the outer corona. The filament eruption also triggered transverse kink-mode oscillation of the adjacent coronal loops in the same AR. The amplitude and period of the oscillation were 1.6 Mm and 225 s. Our results are important for understanding the mechanisms of partial filament eruptions and provide new constraints to theoretical models. The multiwavelength observations also shed light on space weather prediction.

The geometrical shapes and the physical properties of stellar wind -- interstellar medium interaction regions form an important stage for studying stellar winds and their embedded magnetic fields as well as cosmic ray modulation. Our goal is to provide a proper representation and classification of counter-flow configurations and counter-flow interfaces in the frame of fluid theory. In addition we calculate flows and large-scale electromagnetic fields based on which the large-scale dynamics and its role as possible background for particle acceleration, e.g. in the form of anomalous cosmic rays, can be studied. We find that for the definition of the boundaries, which are determining the astropause shape, the number and location of magnetic null points and stagnation points is essential. Multiple separatrices can exist, forming a highly complex environment for the interstellar and stellar plasma. Furthermore, the formation of extended tail structures occur naturally, and their stretched field and streamlines provide surroundings and mechanisms for the acceleration of particles by field-aligned electric fields.

This publication provides an overview of magnetic fields in the solar atmosphere with the focus lying on the corona. The solar magnetic field couples the solar interior with the visible surface of the Sun and with its atmosphere. It is also responsible for all solar activity in its numerous manifestations. Thus, dynamic phenomena such as coronal mass ejections and flares are magnetically driven. In addition, the field also plays a crucial role in heating the solar chromosphere and corona as well as in accelerating the solar wind. Our main emphasis is the magnetic field in the upper solar atmosphere so that photospheric and chromospheric magnetic structures are mainly discussed where relevant for higher solar layers. Also, the discussion of the solar atmosphere and activity is limited to those topics of direct relevance to the magnetic field. After giving a brief overview about the solar magnetic field in general and its global structure, we discuss in more detail the magnetic field in active regions, the quiet Sun and coronal holes.

We carry out a comprehensive investigation comparing the three-dimensional magnetic field restructuring, flare energy release, and the helioseismic response, of two homologous flares, the 2011 September 6 X2.1 (FL1) and September 7 X1.8 (FL2) flares in NOAA AR 11283. In our analysis, (1) a twisted flux rope (FR) collapses onto the surface at a speed of 1.5 km/s after a partial eruption in FL1. The FR then gradually grows to reach a higher altitude and collapses again at 3 km/s after a fuller eruption in FL2. Also, FL2 shows a larger decrease of the flux-weighted centroid separation of opposite magnetic polarities and a greater change of the horizontal field on the surface. These imply a more violent coronal implosion with corresponding more intense surface signatures in FL2. (2) The FR is inclined northward, and together with the ambient fields, it undergoes a southward turning after both events. This agrees with the asymmetric decay of the penumbra observed in the peripheral regions. (3) The amounts of free magnetic energy and nonthermal electron energy released during FL1 are comparable to those of FL2 within the uncertainties of the measurements. (4) No sunquake was detected in FL1; in contrast, FL2 produced two seismic emission sources S1 and S2 both lying in the penumbral regions. Interestingly, S1 and S2 are connected by magnetic loops, and the stronger source S2 has weaker vertical magnetic field. We discuss these results in relation to the implosion process in the low corona and the sunquake generation.

The solar coronal magnetic field produces solar activity, including extremely energetic solar flares and coronal mass ejections (CMEs). Knowledge of the structure and evolution of the magnetic field of the solar corona is important for investigating and understanding the origins of space weather. Although the coronal field remains difficult to measure directly, there is considerable interest in accurate modeling of magnetic fields in and around sunspot regions on the Sun using photospheric vector magnetograms as boundary data. In this work, we investigate effects of the size of the domain chosen for coronal magnetic field modeling on resulting model solution. We apply spherical Optimization procedure to vector magnetogram data of Helioseismic and Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO) with four Active Region observed on 09 March 2012 at 20:55UT. The results imply that quantities like magnetic flux density, electric current density and free magnetic energy density of ARs of interest are significantly different from the corresponding quantities obtained in the same region within the wider field of view. The difference is even more pronounced in the regions where there are connections to outside the domain.

Magnetic fields and flows in coronal structures, for example, in gradual phases in flares, can be described by 2D and 3D magnetohydrostatic (MHS) and steady magnetohydrodynamic (MHD) equilibria. Within a physically simplified, but exact mathematical model, we study the electric currents and corresponding electric fields generated by shear flows. Starting from exact and analytically calculated magnetic potential fields, we solveid the nonlinear MHD equations self-consistently. By applying a magnetic shear flow and assuming a nonideal MHD environment, we calculated an electric field via Faraday's law. The formal solution for the electromagnetic field allowed us to compute an expression of an effective resistivity similar to the collisionless Speiser resistivity. We find that the electric field can be highly spatially structured, or in other words, filamented. The electric field component parallel to the magnetic field is the dominant component and is high where the resistivity has a maximum. The electric field is a potential field, therefore, the highest energy gain of the particles can be directly derived from the corresponding voltage. In our example of a coronal post-flare scenario we obtain electron energies of tens of keV, which are on the same order of magnitude as found observationally. This energy serves as a source for heating and acceleration of particles.

We study the quasi-static evolution of coronal magnetic fields constructed from the Non Linear Force Free Field (NLFFF) approximation aiming to understand the relation between the magnetic field topology and ribbon emission during an X1.5 flare in active region (AR) NOAA 11166. The flare with a quasi-elliptical, and two remote ribbons occurred on March 9, 2011 at 23:13UT over a positive flux region surrounded by negative flux at the center of the bipolar AR. Our analysis of the coronal magnetic structure with potential and NLFFF solutions unveiled the existence of a single magnetic null point associated with a fan-spine topology and is co-spatial with the hard X-ray source. The footpoints of the fan separatrix surface agree with the inner edge of the quasi-elliptical ribbon and the outer spine is linked to one of the remote ribbons. During the evolution, the slow footpoint motions stressed the fieldlines along the polarity inversion line and caused electric current layers in the corona around the fan separatrix surface. These current layers trigger magnetic reconnection as a consequence of dissipating currents, which are visible as cusped shape structures at lower heights. The reconnection process reorganised the magnetic field topology whose signatures are observed at the separatrices/QSL structure both in the photosphere and corona during the pre-to-post flare evolution. In agreement with previous numerical studies, our results suggest that the line-tied footpoint motions perturb the fan-spine system and cause null point reconnection, which eventually causes the flare emission at the footpoints of the fieldlines.

Photospheric vector magnetograms from Helioseismic and Magnetic Imager on board the Solar Dynamic Observatory are utilized as the boundary conditions to extrapolate both non-linear force-free and potential magnetic fields in solar corona. Based on the extrapolations, we are able to determine the free magnetic energy (FME) stored in active regions (ARs). Over 3000 vector magnetograms in 61 ARs were analyzed. We compare FME with ARs' flare index (FI) and find that there is a weak correlation ($<60\%$) between FME and FI. FME shows slightly improved flare predictability relative to total unsigned magnetic flux of ARs in the following two aspects: (1) the flare productivity predicted by FME is higher than that predicted by magnetic flux and (2) the correlation between FI and FME is higher than that between FI and magnetic flux. However, this improvement is not significant enough to make a substantial difference in time-accumulated FI, rather than individual flare, predictions.

We studied the temporal evolution of the magnetic topology of the active region (AR) 11158 based on the reconstructed three-dimensional magnetic fields in the corona. The \nlfff\ extrapolation method was applied to the 12 minutes cadence data obtained with the \hmi\ (HMI) onboard the \sdo\ (SDO) during five days. By calculating the squashing degree factor Q in the volume, the derived quasi-separatrix layers (QSLs) show that this AR has an overall topology, resulting from a magnetic quadrupole, including an hyperbolic flux tube (HFT) configuration which is relatively stable at the time scale of the flare ($\sim 1-2$ hours). A strong QSL, which corresponds to some highly sheared arcades that might be related to the formation of a flux rope, is prominent just before the M6.6 and X2.2 flares, respectively. These facts indicate the close relationship between the strong QSL and the high flare productivity of AR 11158. In addition, with a close inspection of the topology, we found a small-scale HFT which has an inverse tear-drop structure above the aforementioned QSL before the X2.2 flare. It indicates the existence of magnetic flux rope at this place. Even though a global configuration (HFT) is recognized in this AR, it turns out that the large-scale HFT only plays a secondary role during the eruption. In final, we dismiss a trigger based on the breakout model and highlight the central role of the flux rope in the related eruption.

The theoretical concept that braided magnetic field lines in the solar corona may dissipate a sufficient amount of energy to account for the brightening observed in the active-region corona, has been substantiated by high-resolution observations only recently. From the analysis of coronal images obtained with the High Resolution Coronal Imager, first observational evidence of the braiding of magnetic field lines was reported by Cirtain et al. 2013 (hereafter CG13). We present nonlinear force-free reconstructions of the associated coronal magnetic field based on vector SDO/HMI magnetograms. We deliver estimates of the free magnetic energy associated to a braided coronal structure. Our model results suggest (~100 times) more free energy at the braiding site than analytically estimated by CG13, strengthening the possibility of the active-region corona being heated by field line braiding. We were able to assess the coronal free energy appropriately by using vector field measurements and attribute the lower energy estimate of CG13 to the underestimated (by a factor of 10) azimuthal field strength. We also quantify the increase of the overall twist of a flare-related flux rope which had been claimed by CG13. From our models we find that the overall twist of the flux rope increased by about half a turn within twelve minutes. Unlike another method, to which we compare our results to, we evaluate the winding of the flux rope's constituent field lines around each other purely based on their modeled coronal 3D field line geometry -- to our knowledge for the first time.

Magnetic fields in the solar corona are responsible for a wide range of phenomena. However, any direct measurements of the coronal magnetic fields are very difficult due to lack of suitable spectral lines, weak magnetic fields, and high temperatures. Therefore, one extrapolates photospheric field measurements into the corona. Owing to low coronal plasma $\beta$, we can apply a force-free model in lowest order to study the slow evolution of active region (AR) magnetic fields. On applying these models to AR 11158 and compared with coronal plasma tracers, we found that (1) the approximation of potential field to coronal structures over large length scales is a reasonable one, 2) linear force-free (LFF) assumption to AR coronal fields may not be applicable model as it assumes uniform twist over the entire AR, and 3) for modeling fields at sheared, stressed locations where energy release in the form of flares are usually observed, non-linear force free fields (NLFFF) seem to provide a good approximation. The maximum available free-energy profile shows step-wise decrease that is sufficient to power an M-class flare as observed.

Solar eruptive phenomena, like flares and coronal mass ejections(CMEs) are governed by magnetic fields. To describe the structure of these phenomena one needs information on the magnetic flux density and the electric current density vector components in three dimensions throughout the atmosphere. However, current spectro-polarimetric measurements typically limit the determination of the vector magnetic field only to the photosphere. Therefore, there is considerable interest in accurate modeling of the solar coronal magnetic field using photospheric vector magnetograms as boundary data. In this work, we model the coronal magnetic field for global solar atmosphere using a nonlinear force-free field(NLFFF) extrapolation codes implemented to a synoptic maps of photospheric vector magnetic field synthesized from Vector Spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) as boundary condition. Using the resulting three dimensional magnetic field, we calculate the three dimensional electric current density and free magnetic energy throughout the solar atmosphere for Carrington rotation 2124. We found that spatially, the low-lying, current-carrying core field demonstrates strong concentration of free energy in the AR core, from the photosphere to the lower corona; the coronal field becomes slightly more sheared in the lowest layer and it relaxes to the potential field configuration with height. The free energy density appears largely co-spatial with the electric current distribution.

Magnetic reconnection is one of the primary mechanisms for triggering solar eruptive events, but direct observation of its rapid process has been of challenge. In this Letter we present, using a nonlinear force-free field (NLFFF) extrapolation technique, a visualization of field line connectivity changes resulting from tether-cutting reconnection over about 30 minutes during the 2011 February 13 M6.6 flare in NOAA AR 11158. Evidence for the tether-cutting reconnection was first collected through multiwavelength observations and then by the analysis of the field lines traced from positions of four conspicuous flare 1700 A footpoints observed at the event onset. Right before the flare, the four footpoints are located very close to the regions of local maxima of magnetic twist index. Especially, the field lines from the inner two footpoints form two strongly twisted flux bundles (up to ~1.2 turns), which shear past each other and reach out close to the outer two footpoints, respectively. Immediately after the flare, the twist index of regions around the footpoints greatly diminish and the above field lines become low lying and less twisted (~0.6 turns), overarched by loops linking the later formed two flare ribbons. About 10% of the flux (~3x10^19 Mx) from the inner footpoints has undergone a footpoint exchange. This portion of flux originates from the edge regions of the inner footpoints that are brightened first. These rapid changes of magnetic field connectivity inferred from the NLFFF extrapolation are consistent with the tether-cutting magnetic reconnection model.

The magnetic field permeating the solar atmosphere is generally thought to provide the energy for much of the activity seen in the solar corona, such as flares, coronal mass ejections (CMEs), etc. To overcome the unavailability of coronal magnetic field measurements, photospheric magnetic field vector data can be used to reconstruct the coronal field. Currently there are several modelling techniques being used to calculate three-dimension of the field lines into the solar atmosphere. For the first time, synoptic maps of photospheric vector magnetic field synthesized from Vector Spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) are used to model the coronal magnetic field and estimate free magnetic energy in the global scale. The free energy (i.e., the energy in excess of the potential field energy) is one of the main indicators used in space weather forecasts to predict the eruptivity of active regions. We solve the nonlinear force-free field equations using optimization principle in spherical geometry. The resulting three-dimensional magnetic fields are used to estimate the magnetic free energy content E_{free}=E_{nlfff}-E_{pot}, i.e., the difference of the magnetic energies between the nonpotential field and the potential field in the global solar corona. For comparison, we overlay the extrapolated magnetic field lines with the extreme ultraviolet (EUV) observations by the Atmospheric Imaging Assembly on board SDO. For a single Carrington rotation 2121, we find that the global NLFFF magnetic energy density is 10.3% higher than the potential one. Most of this free energy is located in active regions.

We consider a magnetic configuration consisting of an arcade structure and a detached plasmoid, resulting from a magnetic reconnection process, as is typically found in connection with solar flares. We study spontaneous current fragmentation caused by shear and vortex plasma flows. An exact analytical transformation method was applied to calculate self-consistent solutions of the nonlinear stationary MHD equations. The assumption of incompressible field-aligned flows implies that both the Alfven Mach number and the mass density are constant on field lines. We first calculated nonlinear MHS equilibria with the help of the Liouville method, emulating the scenario of a solar eruptive flare configuration with plasmoids and flare arcade. Then a Mach number profile was constructed that describes the upflow along the open magnetic field lines and implements a vortex flow inside the plasmoid. This Mach number profile was used to map the MHS equilibrium to the stationary one. We find that current fragmentation takes place at different locations within our configuration. Steep gradients of the Alfven Mach number are required, implying the strong influence of shear flows on current amplification and filamentation of the MHS current sheets. Crescent- or ring-like structures appear along the outer separatrix, butterfly structures between the upper and lower plasmoids, and strong current peaks close the lower boundary. Impressing an intrinsic small-scale structure on the upper plasmoid results in strong fragmentation of the plasmoid. Hence fragmentation of current sheets and plasmoids is an inherent property of MHD theory. Transformations from MHS into MHD steady-states deliver fine-structures needed for plasma heating and acceleration of particles and bulk plasma flows in dissipative events that are typically connected to magnetic reconnection processes in flares and coronal mass ejections.

Photospheric magnetic vector maps from two different instruments are used to model the nonlinear force-free coronal magnetic field above an active region. We use vector maps inferred from polarization measurements of the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (HMI) and the Solar Optical Telescope Spectropolarimeter (SP) aboard Hinode. Besides basing our model calculations on HMI data, we use both, SP data of original resolution and scaled down to the resolution of HMI. This allows us to compare the model results based on data from different instruments and to investigate how a binning of high-resolution data effects the model outcome. The resulting 3D magnetic fields are compared in terms of magnetic energy content and magnetic topology. We find stronger magnetic fields in the SP data, translating into a higher total magnetic energy of the SP models. The net Lorentz forces of the HMI and SP lower boundaries verify their force-free compatibility. We find substantial differences in the absolute estimates of the magnetic field energy but similar relative estimates, e.g., the fraction of excess energy and of the flux shared by distinct areas. The location and extension of neighboring connectivity domains differs and the SP model fields tend to be higher and more vertical. Hence, conclusions about the magnetic connectivity based on force-free field models are to be drawn with caution. We find that the deviations of the model solution when based on the lower-resolution SP data are small compared to the differences of the solutions based on data from different instruments.

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 11, 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 with observations. During October 11, 2011 observation, there are substantial number of trans-equatorial loops carrying electric current.

Measurements of magnetic fields and electric currents in the pre-eruptive corona are crucial to study solar eruptive phenomena, like flare and coronal mass ejections(CMEs). However, spectro-polarimetric measurements of certain photospheric lines permit a determination of the vector magnetic field at the photosphere. Thus, substantial collection of magnetograms relate to the photospheric surface field only. Numerical modeling is carried out by applying state-of-the-art nonlinear force-free field (NLFFF) reconstruction. Cartesian nonlinear force-free field (NLFFF) codes are not well suited for larger domains, since the spherical nature of the solar surface cannot be neglected when the field of view is large. One of the most significant results of Solar Dynamic Observatory (SDO) mission to date has been repeated observations of large, almost global scale events in which large scale connection between active regions may play fundamental role. Therefore, it appears prudent to implement a NLFFF procedure in spherical geometry for use when large scale boundary data are available, such as from the Helioseismic and Magnetic Imager (HMI) on board SDO. In this work, we model the coronal magnetic field above multiple active regions with the help of a potential field and a NLFFF extrapolation codes in a full-disk using HMI data as a boundary conditions. We compare projections of the resulting magnetic field lines solutions with full-disk coronal images from the Atmospheric Imaging Assembly (SDO/AIA) for both models. This study has found that the NLFFF model reconstructs the magnetic configuration better than the potential field model. We have concluded that much of trans-equatorial loops connecting the two solar hemispheres are current-free.

The moss is the area at the footpoint of the hot (3 to 5 MK) loops forming the core of the active region where emission is believed to result from the heat flux conducted down to the transition region from the hot loops. Studying the variation of Doppler shift as a function of line formation temperatures over the moss area can give clues on the heating mechanism in the hot loops in the core of the active regions. We investigate the absolute Doppler shift of lines formed at temperatures between 1 MK and 2 MK in a moss area within active region NOAA 11243 using a novel technique that allows determining the absolute Doppler shift of EUV lines by combining observations from the SUMER and EIS spectrometers. The inner (brighter and denser) part of the moss area shows roughly constant blue shift (upward motions) of 5 km/s in the temperature range of 1 MK to 1.6 MK. For hotter lines the blue shift decreases and reaches 1 km/s for Fe xv 284 {\AA} (~2 MK). The measurements are discussed in relation to models of the heating of hot loops. The results for the hot coronal lines seem to support the quasi-steady heating models for non-symmetric hot loops in the core of active regions.

Extrapolation codes in Cartesian geometry for modelling the magnetic field in the corona do not take the curvature of the Sun's surface into account and can only be applied to relatively small areas, e.g., a single active region. We compare the analysis of the photospheric magnetic field and subsequent force-free modeling based on full-disk vector maps from Helioseismic and Magnetic Imager (HMI) on board solar dynamics observatory (SDO) and Vector Spectromagnetograph (VSM) of the Synoptic Optical Long-term Investigations of the Sun (SOLIS). We use Helioseismic and Magnetic Imager and Vector Spectromagnetograph photospheric magnetic field measurements to model the force-free coronal field above multiple solar active regions, assuming magnetic forces to dominate. We solve the nonlinear force-free field equations by minimizing a functional in spherical coordinates over a full disk excluding the poles. After searching for the optimum modeling parameters for the particular data sets, we compare the resulting nonlinear force-free model fields. We compare quantities like the total magnetic energy content and free magnetic energy, the longitudinal distribution of the magnetic pressure and surface electric current density using our spherical geometry extrapolation code. The magnetic field lines obtained from nonlinear force-free extrapolation based on Helioseismic and Magnetic Imager and Vector Spectromagnetograph data have good agreement. However, the nonlinear force-free extrapolation based on Helioseismic and Magnetic Imager data have more contents of total magnetic energy, free magnetic energy, the longitudinal distribution of the magnetic pressure and surface electric current density compared to the one from Vector Spectromagnetograph data.

We test a nonlinear force-free field (NLFFF) optimization code in spherical geometry using an analytical solution from Low and Lou. Several tests are run, ranging from idealized cases where exact vector field data are provided on all boundaries, to cases where noisy vector data are provided on only the lower boundary (approximating the solar problem). Analytical tests also show that the NLFFF code in the spherical geometry performs better than that in the Cartesian one when the field of view of the bottom boundary is large, say, $20^\circ \times 20^\circ$. Additionally, We apply the NLFFF model to an active region observed by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) both before and after an M8.7 flare. For each observation time, we initialize the models using potential field source surface (PFSS) extrapolations based on either a synoptic chart or a flux-dispersal model, and compare the resulting NLFFF models. The results show that NLFFF extrapolations using the flux-dispersal model as the boundary condition have slightly lower, therefore better, force-free and divergence-free metrics, and contain larger free magnetic energy. By comparing the extrapolated magnetic field lines with the extreme ultraviolet (EUV) observations by the Atmospheric Imaging Assembly (AIA) on board SDO, we find that the NLFFF performs better than the PFSS not only for the core field of the flare productive region, but also for large EUV loops higher than 50 Mm.

The structure and dynamics of the solar corona is dominated by the magnetic field. In most areas in the corona magnetic forces are so dominant that all non-magnetic forces like plasma pressure gradient and gravity can be neglected in the lowest order. This model assumption is called the force-free field assumption, as the Lorentz force vanishes. This can be obtained by either vanishing electric currents (leading to potential fields) or the currents are co-aligned with the magnetic field lines. First we discuss a mathematically simpler approach that the magnetic field and currents are proportional with one global constant, the so-called linear force-free field approximation. In the generic case, however, the relation between magnetic fields and electric currents is nonlinear and analytic solutions have been only found for special cases, like 1D or 2D configurations. For constructing realistic nonlinear force-free coronal magnetic field models in 3D, sophisticated numerical computations are required and boundary conditions must be obtained from measurements of the magnetic field vector in the solar photosphere. This approach is currently of large interests, as accurate measurements of the photospheric field become available from ground-based (for example SOLIS) and space-born (for example Hinode and SDO) instruments. If we can obtain accurate force-free coronal magnetic field models we can calculate the free magnetic energy in the corona, a quantity which is important for the prediction of flares and coronal mass ejections. 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.

We use SDO/HMI and SOLIS/VSM photospheric magnetic field measurements to model the force-free coronal field above a solar active region, assuming magnetic forces to dominate. We take measurement uncertainties caused by, e.g., noise and the particular inversion technique into account. After searching for the optimum modeling parameters for the particular data sets, we compare the resulting nonlinear force-free model fields. We show the degree of agreement of the coronal field reconstructions from the different data sources by comparing the relative free energy content, the vertical distribution of the magnetic pressure and the vertically integrated current density. Though the longitudinal and transverse magnetic flux measured by the VSM and HMI is clearly different, we find considerable similarities in the modeled fields. This indicates the robustness of the algorithm we use to calculate the nonlinear force-free fields against differences and deficiencies of the photospheric vector maps used as an input. We also depict how much the absolute values of the total force-free, virial and the free magnetic energy differ and how the orientation of the longitudinal and transverse components of the HMI- and VSM-based model volumes compares to each other.

Magnetohydrodynamic configurations with strong localized current concentrations and vortices play an important role for the dissipation of energy in space and astrophysical plasma. Within this work we investigate the relation between current sheets and vortex sheets in incompressible, stationary equilibria. For this approach it is helpful that the similar mathematical structure of magnetohydrostatics and stationary incompressible hydrodynamics allows us to transform static equilibria into stationary ones. The main control function for such a transformation is the profile of the Alfven-Mach number M_A, which is always constant along magnetic field lines, but can change from one field line to another. In the case of a global constant M_A, vortices and electric current concentrations are parallel. More interesting is the nonlinear case, where M_A varies perpendicular to the field lines. This is a typical situation at boundary layers like the magnetopause, heliopause, the solar wind flowing around helmet streamers and at the boundary of solar coronal holes. The corresponding current and vortex sheets show in some cases also an alignment, but not in every case. For special density distributions in 2D it is possible to have current but no vortex sheets. In 2D vortex sheets of field aligned-flows can also exist without strong current sheets, taking the limit of small Alfven Mach numbers into account. The current sheet can vanish if the Alfven Mach number is (almost) constant and the density gradient is large across some boundary layer. It should be emphasized that the used theory is not only valid for small Alfven Mach numbers M_A<<1, but also for M_A~1. Connection to other theoretical approaches and observations and physical effects in space plasmas are presented. Differences in the various aspects of theoretical investigations of current sheets and vortex sheets are given.

In this study, photospheric vector magnetograms obtained with the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS), are used as boundary conditions to model the three-dimensional nonlinear force-free (NLFF) coronal magnetic fields as a sequence of nonlinear force-free equilibria in spherical geometry. We study the coronal magnetic field structure inside active regions and its temporal evolution. We compare the magnetic field configuration obtained from NLFF extrapolation before and after flaring event in active region (AR) 11117 and its surroundings observed on 27 October 2010. We compare the magnetic field topologies and the magnetic energy densities and study the connectivities between AR 11117 and its surroundings. During the investigated time period, we estimate the change in free magnetic energy from before to after the flare to be 1.74x10^{32}erg which represents about 13.5% of nonlinear force-free magnetic energy before the flare. In this study, we find that electric currents from AR 11117 to its surroundings were disrupted after the flare.

The SDO/HMI instruments provide photospheric vector magnetograms with a high spatial and temporal resolution. Our intention is to model the coronal magnetic field above active regions with the help of a nonlinear force-free extrapolation code. Our code is based on an optimization principle and has been tested extensively with semi-analytic and numeric equilibria and been applied before to vector magnetograms from Hinode and ground based observations. Recently we implemented a new version which takes measurement errors in photospheric vector magnetograms into account. Photospheric field measurements are often due to measurement errors and finite nonmagnetic forces inconsistent as a boundary for a force-free field in the corona. In order to deal with these uncertainties, we developed two improvements: 1.) Preprocessing of the surface measurements in order to make them compatible with a force-free field 2.) The new code keeps a balance between the force-free constraint and deviation from the photospheric field measurements. Both methods contain free parameters, which have to be optimized for use with data from SDO/HMI. Within this work we describe the corresponding analysis method and evaluate the force-free equilibria by means of how well force-freeness and solenoidal conditions are fulfilled, the angle between magnetic field and electric current and by comparing projections of magnetic field lines with coronal images from SDO/AIA. We also compute the available free magnetic energy and discuss the potential influence of control parameters.

The rapid, irreversible change of the photospheric magnetic field has been recognized as an important element of the solar flare process. This Letter reports such a rapid change of magnetic fields during the 2011 February 13 M6.6 flare in NOAA AR 11158 that we found from the vector magnetograms of the Helioseismic and Magnetic Imager with 12-min cadence. High-resolution magnetograms of Hinode that are available at ~-5.5, -1.5, 1.5, and 4 hrs relative to the flare maximum are used to reconstruct three-dimensional coronal magnetic field under the nonlinear force-free field (NLFFF) assumption. UV and hard X-ray images are also used to illuminate the magnetic field evolution and energy release. The rapid change is mainly detected by HMI in a compact region lying in the center of the magnetic sigmoid, where the mean horizontal field strength exhibited a significant increase by 28%. The region lies between the initial strong UV and hard X-ray sources in the chromosphere, which are cospatial with the central feet of the sigmoid according to the NLFFF model. The NLFFF model further shows that strong coronal currents are concentrated immediately above the region, and that more intriguingly, the coronal current system underwent an apparent downward collapse after the sigmoid eruption. These results are discussed in favor of both the tether-cutting reconnection producing the flare and the ensuing implosion of the coronal field resulting from the energy release.

Aims. We analyze coronagraph observations of a polar jet observed by the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) instrument suite onboard the Solar TErrestrial RElations Observatory (STEREO) spacecraft. Methods. In our analysis we compare the brightness distribution of the jet in white-light coronagraph images with a dedicated kinetic particle model. We obtain a consistent estimate of the time that the jet was launched from the solar surface and an approximate initial velocity distribution in the jet source. The method also allows us to check the consistency of the kinetic model. In this first application, we consider only gravity as the dominant force on the jet particles along the magnetic field. Results. We find that the kinetic model explains the observed brightness evolution well. The derived initiation time is consistent with the jet observations by the EUVI telescope at various wavelengths. The initial particle velocity distribution is fitted by Maxwellian distributions and we find deviations of the high energy tail from the Maxwellian distributions. We estimate the jet's total electron content to have a mass between 3.2 \times 1014 and 1.8 \times 1015 g. Mapping the integrated particle number along the jet trajectory to its source region and assuming a typical source region size, we obtain an initial electron density between 8 \times 109 and 5 \times 1010 cm-3 that is characteristic for the lower corona or the upper chromosphere. The total kinetic energy of all particles in the jet source region amounts from 2.1 \times 1028 to 2.4 \times 1029 erg.

Polar plumes are seen as elongated objects starting at the solar polar regions. Here, we analyze these objects from a sequence of images taken simultaneously by the three spacecraft telescopes STEREO/EUVI A and B, and SOHO/EIT. We establish a method capable of automatically identifying plumes in solar EUV images close to the limb at 1.01 - 1.39 R in order to study their temporal evolution. This plume-identification method is based on a multiscale Hough-wavelet analysis. Then two methods to determined their 3D localization and structure are discussed: First, tomography using the filtered back-projection and including the differential rotation of the Sun and, secondly, conventional stereoscopic triangulation. We show that tomography and stereoscopy are complementary to study polar plumes. We also show that this systematic 2D identification and the proposed methods of 3D reconstruction are well suited, on one hand, to identify plumes individually and on the other hand, to analyze the distribution of plumes and inter-plume regions. Finally, the results are discussed focusing on the plume position with their cross-section area.

Understanding the solar outer atmosphere requires concerted, simultaneous solar observations from the visible to the vacuum ultraviolet (VUV) and soft X-rays, at high spatial resolution (between 0.1" and 0.3"), at high temporal resolution (on the order of 10 s, i.e., the time scale of chromospheric dynamics), with a wide temperature coverage (0.01 MK to 20 MK, from the chromosphere to the flaring corona), and the capability of measuring magnetic fields through spectropolarimetry at visible and near-infrared wavelengths. Simultaneous spectroscopic measurements sampling the entire temperature range are particularly important. These requirements are fulfilled by the Japanese Solar-C mission (Plan B), composed of a spacecraft in a geosynchronous orbit with a payload providing a significant improvement of imaging and spectropolarimetric capabilities in the UV, visible, and near-infrared with respect to what is available today and foreseen in the near future. The Large European Module for solar Ultraviolet Research (LEMUR), described in this paper, is a large VUV telescope feeding a scientific payload of high-resolution imaging spectrographs and cameras. LEMUR consists of two major components: a VUV solar telescope with a 30 cm diameter mirror and a focal length of 3.6 m, and a focal-plane package composed of VUV spectrometers covering six carefully chosen wavelength ranges between 17 and 127 nm. The LEMUR slit covers 280" on the Sun with 0.14" per pixel sampling. In addition, LEMUR is capable of measuring mass flows velocities (line shifts) down to 2 km/s or better. LEMUR has been proposed to ESA as the European contribution to the Solar C mission.

Extrapolation codes for modelling the magnetic field in the corona in cartesian geometry do not take the curvature of the Sun's surface into account and can only be applied to relatively small areas, \textit{e.g.}, a single active region. We apply a method for nonlinear force-free coronal magnetic field modelling of photospheric vector magnetograms in spherical geometry which allows us to study the connectivity between multi-active regions. We use vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS)/Vector Spectromagnetograph(VSM) to model the coronal magnetic field, where we study three neighbouring magnetically connected active regions (ARs: 10987, 10988, 10989) observed on 28, 29, and 30 March 2008, respectively. We compare the magnetic field topologies and the magnetic energy densities and study the connectivities between the active regions(ARs). We have studied the time evolution of magnetic field over the period of three days and found no major changes in topologies as there was no major eruption event. From this study we have concluded that active regions are much more connected magnetically than the electric current.

The relation of the solar surface magnetic field with mesogranular cells is studied using high spatial (~ 100 km) and temporal (~ 30 sec) resolution data obtained with the IMaX instrument aboard SUNRISE. First, mesogranular cells are identified using Lagrange tracers (corks) based on horizontal velocity fields obtained through Local Correlation Tracking. After ~ 20 min of integration, the tracers delineate a sharp mesogranular network with lanes of width below about 280 km. The preferential location of magnetic elements in mesogranular cells is tested quantitatively. Roughly 85% of pixels with magnetic field higher than 100 G are located in the near neighborhood of mesogranular lanes. Magnetic flux is therefore concentrated in mesogranular lanes rather than intergranular ones. Secondly, magnetic field extrapolations are performed to obtain field lines anchored in the observed flux elements. This analysis, therefore, is independent of the horizontal flows determined in the first part. A probability density function (PDF) is calculated for the distribution of distances between the footpoints of individual magnetic field lines. The PDF has an exponential shape at scales between 1 and 10 Mm, with a constant characteristic decay distance, indicating the absence of preferred convection scales in the mesogranular range. Our results support the view that mesogranulation is not an intrinsic convective scale (in the sense that it is not a primary energy-injection scale of solar convection), but also give quantitative confirmation that, nevertheless, the magnetic elements are preferentially found along mesogranular lanes.

Routine measurements of the solar magnetic field are mainly carried out in the photosphere. Therefore, one has to infer the field strength in the higher layers of the solar atmosphere from the measured photospheric field based on the assumption that the corona is force-free. Meanwhile, those measured data are inconsistent with the above force-free assumption. Therefore, one has to apply some transformations to these data before nonlinear force-free extrapolation codes can be applied. Extrapolation codes in cartesian geometry for modelling the magnetic field in the corona do not take the curvature of the Sun's surface into account and can only be applied to relatively small areas, e.g., a single active region. Here we apply a method for nonlinear force-free coronal magnetic field modelling and preprocessing of photospheric vector magnetograms in spherical geometry using the optimization procedure.We solve the nonlinear force-free field equations by minimizing a functional in spherical coordinates over a restricted area of the Sun. We extend the functional by an additional term, which allows to incorporate measurement error and treat regions with lacking observational data. We use vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS) to model the coronal magnetic field. We study two neighbouring magnetically connected active regions observed on May 15 2009. For vector magnetograms with variable measurement precision and randomly scattered data gaps (e.g., SOLIS/VSM) the new code yields field models which satisfy the solenoidal and force-free condition significantly better as it allows deviations between the extrapolated boundary field and observed boundary data within measurement errors. Data gaps are assigned to an infinite error. We extend this new scheme to spherical geometry and apply it for the first time to real data.

Extrapolations of solar photospheric vector magnetograms into three-dimensional magnetic fields in the chromosphere and corona are usually done under the assumption that the fields are force-free. The field calculations can be improved by preprocessing the photospheric magnetograms. We compare two preprocessing methods presently in use, namely the methods of Wiegelmann et al. (2006) and Fuhrmann et al. (2007). The two preprocessing methods were applied to a recently observed vector magnetogram. We examine the changes in the magnetogram effected by the two preprocessing algorithms. Furthermore, the original magnetogram and the two preprocessed magnetograms were each used as input data for nonlinear force-free field extrapolations by means of two different methods, and we analyze the resulting fields. Both preprocessing methods managed to significantly decrease the magnetic forces and magnetic torques that act through the magnetogram area and that can cause incompatibilities with the assumption of force-freeness in the solution domain. Both methods also reduced the amount of small-scale irregularities in the observed photospheric field, which can sharply worsen the quality of the solutions. For the chosen parameter set, the Wiegelmann et al. method led to greater changes in strong-field areas, leaving weak-field areas mostly unchanged, and thus providing an approximation of the magnetic field vector in the chromosphere, while the Fuhrmann et al. method weakly changed the whole magnetogram, thereby better preserving patterns present in the original magnetogram. Both preprocessing methods raised the magnetic energy content of the extrapolated fields to values above the minimum energy, corresponding to the potential field. Also, the fields calculated from the preprocessed magnetograms fulfill the solenoidal condition better than those calculated without preprocessing.

We investigate the fine structure of magnetic fields in the atmosphere of the quiet Sun. We use photospheric magnetic field measurements from {\sc Sunrise}/IMaX with unprecedented spatial resolution to extrapolate the photospheric magnetic field into higher layers of the solar atmosphere with the help of potential and force-free extrapolation techniques. We find that most magnetic loops which reach into the chromosphere or higher have one foot point in relatively strong magnetic field regions in the photosphere. $91%$ of the magnetic energy in the mid chromosphere (at a height of 1 Mm) is in field lines, whose stronger foot point has a strength of more than 300 G, i.e. above the equipartition field strength with convection. The loops reaching into the chromosphere and corona are also found to be asymmetric in the sense that the weaker foot point has a strength $B < 300$ G and is located in the internetwork. Such loops are expected to be strongly dynamic and have short lifetimes, as dictated by the properties of the internetwork fields.

Strong gradients in plasma flows play a major role in space and astrophysical plasmas. A typical situation is that a static plasma equilibrium is surrounded by a plasma flow, which can lead to strong plasma flow gradients at the separatrices between field lines with different magnetic topologies, e.g., planetary magnetospheres, helmet streamers in the solar corona, or at the boundary between the heliosphere and interstellar medium. Within this work we make a first step to understand the influence of these flows towards the occurrence of current sheets in a stationary state situation. We concentrate here on incompressible plasma flows and 2D equilibria, which allow us to find analytic solutions of the stationary magnetohydrodynamics equations (SMHD). First we solve the magnetohydrostatic (MHS) equations with the help of a Grad-Shafranov equation and then we transform these static equilibria into a stationary state with plasma flow. We are in particular interested to study SMHD-equilibria with strong plasma flow gradients perpendicular to separatrices. We find that induced thin current sheets occur naturally in such situations. The strength of the induced currents depend on the Alfv\'en Mach number and its gradient, and on the magnetic field.

A magnetic channel - a series of polarity reversals separating elongated flux threads with opposite polarities - may be a manifestation of a highly non-potential magnetic configuration in active regions. To understand its formation we have carried out a detailed analysis of the magnetic channel in AR 10930 using data taken by the Solar Optical Telescope/Hinode. As a result, we found upflows (-0.5 to -1.0 km/s) and downflows (+1.5 to +2.0 km/s) inside and at both tips of the thread respectively, and a pair of strong vertical currents of opposite polarity along the channel. Moreover, our analysis of the nonlinear force-free fields constructed from the photospheric magnetic field indicates that the current density in the lower corona may have gradually increased as a result of the continuous emergence of the highly sheared flux along the channel. With these results, we suggest that the magnetic channel originates from the emergence of the twisted flux tube that has formed below the surface before the emergence.

In this letter, we study in detail the evolution of the post-flare loops on 2005 January 15 that occurred between two consecutive solar eruption events, both of which generated a fast halo CME and a major flare. The post-flare loop system, formed after the first CME/flare eruption, evolved rapidly, as manifested by the unusual accelerating rise motion of the loops. Through nonlinear force-free field (NLFFF) models, we obtain the magnetic structure over the active region. It clearly shows that the flux rope below the loops also kept rising accompanied with increasing twist and length. Finally, the post-flare magnetic configuration evolved to a state that resulted in the second CME/flare eruption. This is an event in which the post-flare loops can re-flare in a short period of $\sim$16 hr following the first CME/flare eruption. The observed re-flaring at the same location is likely driven by the rapid evolution of the flux rope caused by the magnetic flux emergence and the rotation of the sunspot. This observation provides valuable information on CME/flare models and their prediction.

Context: Solar magnetic fields are regularly extrapolated into the corona starting from photospheric magnetic measurements that can suffer from significant uncertainties. Aims: Here we study how inaccuracies introduced into the maps of the photospheric magnetic vector from the inversion of ideal and noisy Stokes parameters influence the extrapolation of nonlinear force-free magnetic fields. Methods: We compute nonlinear force-free magnetic fields based on simulated vector magnetograms, which have been produced by the inversion of Stokes profiles, computed froma 3-D radiation MHD simulation snapshot. These extrapolations are compared with extrapolations starting directly from the field in the MHD simulations, which is our reference. We investigate how line formation and instrumental effects such as noise, limited spatial resolution and the effect of employing a filter instrument influence the resulting magnetic field structure. The comparison is done qualitatively by visual inspection of the magnetic field distribution and quantitatively by different metrics. Results: The reconstructed field is most accurate if ideal Stokes data are inverted and becomes less accurate if instrumental effects and noise are included. The results demonstrate that the non-linear force-free field extrapolation method tested here is relatively insensitive to the effects of noise in measured polarization spectra at levels consistent with present-day instruments. Conclusions heading: Our results show that we can reconstruct the coronal magnetic field as a nonlinear force-free field from realistic photospheric measurements with an accuracy of a few percent, at least in the absence of sunspots.

Context: Knowledge about the coronal magnetic field is important to the understanding the structure of the solar corona. We compute the field in the higher layers of the solar atmosphere from the measured photospheric field under the assumption that the corona is force-free. Aims: Here we develop a method for nonlinear force-free coronal magnetic field medelling and preprocessing of photospheric vector magnetograms in spherical geometry using the optimization procedure. Methods: We describe a newly developed code for the extrapolation of nonlinear force-free coronal magnetic fields in spherical coordinates over a restricted area of the Sun. The program uses measured vector magnetograms on the solar photosphere as input and solves the force-free equations in the solar corona. We develop a preprocessing procedure in spherical geometry to drive the observed non-force-free data towards suitable boundary conditions for a force-free extrapolation. Results: We test the code with the help of a semi-analytic solution and assess the quality of our reconstruction qualitatively by magnetic field line plots and quantitatively with a number of comparison metrics for different boundary conditions. The reconstructed fields from the lower boundary data with the weighting function are in good agreement with the original reference fields. We added artificial noise to the boundary conditions and tested the code with and without preprocessing. The preprocessing recovered all main structures of the magnetogram and removed small-scale noise. The main test was to extrapolate from the noisy photospheric vector magnetogram with and without preprocessing. The preprocessing was found to significantly improve the agreement between the extrapolated and the exact field.

Observations from the two STEREO-spacecraft give us for the first time the possibility to use stereoscopic methods to reconstruct the 3D solar corona. Classical stereoscopy works best for solid objects with clear edges. Consequently an application of classical stereoscopic methods to the faint structures visible in the optically thin coronal plasma is by no means straight forward and several problems have to be treated adequately: 1.)First there is the problem of identifying one dimensional structures -e.g. active region coronal loops or polar plumes- from the two individual EUV-images observed with STEREO/EUVI. 2.) As a next step one has the association problem to find corresponding structures in both images. 3.) Within the reconstruction problem stereoscopic methods are used to compute the 3D-geometry of the identified structures. Without any prior assumptions, e.g., regarding the footpoints of coronal loops, the reconstruction problem has not one unique solution. 4.) One has to estimate the reconstruction error or accuracy of the reconstructed 3D-structure, which depends on the accuracy of the identified structures in 2D, the separation angle between the spacecraft, but also on the location, e.g., for east-west directed coronal loops the reconstruction error is highest close to the loop top. 5.) Eventually we are not only interested in the 3D-geometry of loops or plumes, but also in physical parameters like density, temperature, plasma flow, magnetic field strength etc. Helpful for treating some of these problems are coronal magnetic field models extrapolated from photospheric measurements, because observed EUV-loops outline the magnetic field. This feature has been used for a new method dubbed 'magnetic stereoscopy'. As examples we show recent application to active region loops.

We present stereoscopic reconstructions of the location and inclination of polar plumes of two data sets based on the two simultaneously recorded images taken by the EUVI telescopes in the SECCHI instrument package onboard the \emph{STEREO (Solar TErrestrial RElations Observatory)} spacecraft. The ten plumes investigated show a superradial expansion in the coronal hole in 3D which is consistent with the 2D results. Their deviations from the local meridian planes are rather small with an average of $6.47^{\circ}$. By comparing the reconstructed plumes with a dipole field with its axis along the solar rotation axis, it is found that plumes are inclined more horizontally than the dipole field. The lower the latitude is, the larger is the deviation from the dipole field. The relationship between plumes and bright points has been investigated and they are not always associated. For the first data set, based on the 3D height of plumes and the electron density derived from SUMER/\emph{SOHO} Si {\sc viii} line pair, we found that electron densities along the plumes decrease with height above the solar surface. The temperature obtained from the density scale height is 1.6 to 1.8 times larger than the temperature obtained from Mg {\sc ix} line ratios. We attribute this discrepancy to a deviation of the electron and the ion temperatures. Finally, we have found that the outflow speeds studied in the O {\sc vi} line in the plumes corrected by the angle between the line of sight and the plume orientation are quite small with a maximum of 10 $\mathrm{km s^{-1}}$. It is unlikely that plumes are a dominant contributor to the fast solar wind.

We aim at studying the small-scale evolution at the boundaries of an equatorial coronal hole connected with a channel of open magnetic flux with the polar region and an `isolated' one in the extreme-ultraviolet spectral range. We intend to determine the spatial and temporal scale of these changes. Imager data from TRACE in the Fe IX/X 171 A passband and EIT on-board Solar and Heliospheric Observatory in the Fe XII 195 A passband were analysed. We found that small-scale loops known as bright points play an essential role in coronal holes boundaries evolution at small scales. Their emergence and disappearance continuously expand or contract coronal holes. The changes appear to be random on a time scale comparable with the lifetime of the loops seen at these temperatures. No signature was found for a major energy release during the evolution of the loops. Although coronal holes seem to maintain their general shape during a few solar rotations, a closer look at their day-by-day and even hour-by-hour evolution demonstrates a significant dynamics. The small-scale loops (10" - 40" and smaller) which are abundant along coronal hole boundaries have a contribution to the small-scale evolution of coronal holes. Continuous magnetic reconnection of the open magnetic field lines of the coronal hole and the closed field lines of the loops in the quiet Sun is more likely to take place.

We study the correlation heights, which indicate the formation height of Extreme Ultraviolet (EUV) lines in an active region using observations from the EUV Imaging Spectrometer (EIS) and Solar Optical Telescope (SOT) on board \emph{Hinode}. The nonlinear force-free field (NLFFF) optimization method is adopted to extrapolate the 3D magnetic fields to higher layers. Three subregions with different characteristics are selected in the active region for this study. The results show that the formation heights in different subregions vary with their different magnetic fields or velocity patterns. After solving the line blending problem between the He {\sc \romannumeral 2} 256.32 \AA and Si {\sc \romannumeral 10} 256.37 \AA lines by the double Gaussian curve fitting, we find that the transition region lies higher in a strong magnetic area. In a pre-flare heating area there possibly exist multithermal loops as implied by comparing the Doppler velocity and the magnetic field on the solar disk.

Nonlinear force-free field (NLFFF) models are thought to be viable tools for investigating the structure, dynamics and evolution of the coronae of solar active regions. In a series of NLFFF modeling studies, we have found that NLFFF models are successful in application to analytic test cases, and relatively successful when applied to numerically constructed Sun-like test cases, but they are less successful in application to real solar data. Different NLFFF models have been found to have markedly different field line configurations and to provide widely varying estimates of the magnetic free energy in the coronal volume, when applied to solar data. NLFFF models require consistent, force-free vector magnetic boundary data. However, vector magnetogram observations sampling the photosphere, which is dynamic and contains significant Lorentz and buoyancy forces, do not satisfy this requirement, thus creating several major problems for force-free coronal modeling efforts. In this article, we discuss NLFFF modeling of NOAA Active Region 10953 using Hinode/SOT-SP, Hinode/XRT, STEREO/SECCHI-EUVI, and SOHO/MDI observations, and in the process illustrate the three such issues we judge to be critical to the success of NLFFF modeling: (1) vector magnetic field data covering larger areas are needed so that more electric currents associated with the full active regions of interest are measured, (2) the modeling algorithms need a way to accommodate the various uncertainties in the boundary data, and (3) a more realistic physical model is needed to approximate the photosphere-to-corona interface in order to better transform the forced photospheric magnetograms into adequate approximations of nearly force-free fields at the base of the corona. We make recommendations for future modeling efforts to overcome these as yet unsolved problems.

The super-storm of November 20, 2003 was associated with a high speed coronal mass ejection which originated in the NOAA AR 10501 on November 18. This coronal mass ejection had severe terrestrial consequences leading to a geomagnetic storm with DST index of -472 nT, the strongest of the current solar cycle. In this paper, we attempt to understand the factors that led to the coronal mass ejection on November 18. We have also studied the evolution of the photospheric magnetic field of NOAA AR 10501, the source region of this coronal mass ejection. For this purpose, the MDI line-of-sight magnetograms and vector magnetograms from Solar Flare Telescope, Mitaka, obtained during November, 17-19, 2003 were analysed. In particular, quantitative estimates of the temporal variation in magnetic flux, energy and magnetic field gradient were estimated for the source active region. The evolution of these quantities was studied for the 3-day period with an objective to understand the pre-flare configuration leading up to the moderate flare which was associated with the geo-effective coronal mass ejection. We also examined the chromospheric images recorded in H-alpha from Udaipur Solar Observatory to compare the flare location with regions of different magnetic field and energy. Our observations provide evidence that the flare associated with the CME occurred at a location marked by high magnetic field gradient which led to release of free energy stored in the active region.

We study the coronal magnetic field structure inside active regions and its temporal evolution. We attempt to compare the magnetic configuration of an active region in a very quiet period with that for the same region during a flare. Probably for the first time, we use vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS) to model the coronal magnetic field as a sequence of nonlinear force-free equilibria. We study the active region NOAA 10960 observed on 2007 June 7 with three snapshots taken during a small C1.0 flare of time cadence 10 minutes and six snapshots during a quiet period. The total magnetic energy in the active region was approximately $3 \times 10^{25}$ J. Before the flare the free magnetic energy was about 5% of the potential field energy. A part of this excess energy was released during the flare, producing almost a potential configuration at the beginning of the quiet period. During the investigated period, the coronal magnetic energy was only a few percent higher than that of the potential field and consequently only a small C1.0 flare occurred. This was compared with an earlier investigated active region 10540, where the free magnetic energy was about 60% higher than that of the potential field producing two M-class flares. However, the free magnetic energy accumulates before and is released during the flare which appears to be the case for both large and small flares.

The solar magnetic field is key to understanding the physical processes in the solar atmosphere. Nonlinear force-free codes have been shown to be useful in extrapolating the coronal field upward from underlying vector boundary data. However, we can only measure the magnetic field vector routinely with high accuracy in the photosphere, and unfortunately these data do not fulfill the force-free condition. We must therefore apply some transformations to these data before nonlinear force-free extrapolation codes can be self-consistently applied. To this end, we have developed a minimization procedure that yields a more chromosphere-like field, using the measured photospheric field vectors as input. The procedure includes force-free consistency integrals, spatial smoothing, and -- newly included in the version presented here -- an improved match to the field direction as inferred from fibrils as can be observed in, e.g., chromospheric H$\alpha$ images. We test the procedure using a model active-region field that included buoyancy forces at the photospheric level. The proposed preprocessing method allows us to approximate the chromospheric vector field to within a few degrees and the free energy in the coronal field to within one percent.

We present the first reconstruction of the three-dimensional shape of magnetic loops in an active region from two different vantage points based on simultaneously recorded images. The images were taken by the two EUVI telescopes of the SECCHI instrument onboard the recently launched STEREO spacecraft when the heliocentric separation of the two space probes was 12 degrees. We demostrate that these data allow to obtain a reliable three-dimensional reconstruction of sufficiently bright loops. The result is compared with field lines derived from a coronal magnetic field model extrapolated from a photospheric magnetogram recorded nearly simultaneously by SOHO/MDI. We attribute discrepancies between reconstructed loops and extrapolated field lines to the inadequacy of the linear force-free field model used for the extrapolation.

We undertake a first attempt towards a consistent reconstruction of the coronal magnetic field and the coronal density structure. We consider a stationary solar corona which has to obey the equations of magnetohydrostatics. We solve these equations with help of a newly developed optimization scheme. As a first step we illustrate how tomographic information can be included into the reconstruction of coronal magnetic fields. In a second step we use coronal magnetic field information to improve the tomographic inversion process. As input the scheme requires magnetic field measurements on the photosphere from vector-magnetographs and the line-of-sight integrated density distribution from coronagraphs. We test our codes with well known analytic magnetohydrostatic equilibria and models. The program is planed for use within the STEREO mission.

Coronal holes (CH) emit significantly less at coronal temperatures than quiet-Sun regions (QS), but can hardly be distinguished in most chromospheric and lower transition region lines. A key quantity for the understanding of this phenomenon is the magnetic field. We use data from SOHO/MDI to reconstruct the magnetic field in coronal holes and the quiet Sun with the help of a potential magnetic model. Starting from a regular grid on the solar surface we then trace field lines, which provide the overall geometry of the 3D magnetic field structure. We distinguish between open and closed field lines, with the closed field lines being assumed to represent magnetic loops. We then try to compute some properties of coronal loops. The loops in the coronal holes (CH) are found to be on average flatter than in the QS. High and long closed loops are extremely rare, whereas short and low-lying loops are almost as abundant in coronal holes as in the quiet Sun. When interpreted in the light of loop scaling laws this result suggests an explanation for the relatively strong chromospheric and transition region emission (many low-lying, short loops), but the weak coronal emission (few high and long loops) in coronal holes. In spite of this contrast our calculations also suggest that a significant fraction of the cool emission in CHs comes from the open flux regions. Despite these insights provided by the magnetic field line statistics further work is needed to obtain a definite answer to the question if loop statistics explain the differences between coronal holes and the quiet Sun.

We developed a code for the reconstruction of nonlinear force-free and non-force-free coronal magnetic fields. The 3D magnetic field is computed numerically with the help of an optimization principle. The force-free and non-force-free codes are compiled in one program. The force-free approach needs photospheric vector magnetograms as input. The non-force-free code additionally requires the line-of-sight integrated coronal density distribution in combination with a tomographic inversion code. Previously the optimization approach has been used to compute magnetic fields using all six boundaries of a computational box. Here we extend this method and show how the coronal magnetic field can be reconstructed only from the bottom boundary, where the boundary conditions are measured with vector magnetographs. The program is planed for use within the Stereo mission.

We compare the small-scale features visible in the Ne viii Doppler-shift map of an equatorial coronal hole (CH) as observed by SUMER with the small-scale structures of the magnetic field as constructed from a simultaneous photospheric magnetogram by a potential magnetic-field extrapolation. The combined data set is analysed with respect to the small-scale flows of coronal matter, which means that the Ne viii Doppler-shift used as tracer of the plasma flow is investigated in close connection with the ambient magnetic field. Some small closed-field regions in this largely open CH are also found in the coronal volume considered. The Doppler-shift patterns are found to be clearly linked with the field topology.

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.

We compare magnetic field extrapolations from a photospheric magnetogram with the observationally inferred structure of magnetic loops in a newly developed active region. This is the first time that the reconstructed 3D-topology of the magnetic field is available to test the extrapolations. We compare the observations with potential fields, linear force-free fields and non-linear force-free fields. This comparison reveals that a potential field extrapolation is not suitable for a reconstruction of the magnetic field in this young, developing active region. The inclusion of field-line-parallel electric currents, the so called force-free approach, gives much better results. Furthermore, a non-linear force-free computation reproduces the observations better than the linear force-free approximation, although no free parameters are available in the former case.

We compare the performance of two alternative algorithms which aim to construct a force-free magnetic field given suitable boundary conditions. For this comparison, we have implemented both algorithms on the same finite element grid which uses Whitney forms to describe the fields within the grid cells. The additional use of conjugate gradient and multigrid iterations result in quite effective codes. The Grad-Rubin and Wheatland-Sturrock-Roumeliotis algorithms both perform well for the reconstruction of a known analytic force-free field. For more arbitrary boundary conditions the Wheatland-Sturrock-Roumeliotis approach has some difficulties because it requires overdetermined boundary information which may include inconsistencies. The Grad-Rubin code on the other hand loses convergence for strong current densities. For the example we have investigated, however, the maximum possible current density seems to be not far from the limit beyond which a force free field cannot exist anymore for a given normal magnetic field intensity on the boundary.

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 $\alpha$ is determined iteratively by comparing the reconstructed magnetic field with the observed structures. The final configuration is the optimal linear force 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.

We present a procedure which extracts bright loop features from solar EUV images. In terms of image intensities, these features are elongated ridge-like intensity maxima. To discriminate the maxima, we need information about the spatial derivatives of the image intensity. Commonly, the derivative estimates are strongly affected by image noise. We therefore use a regularized estimation of the derivative which is then used to interpolate a discrete vector field of ridge points ``ridgels'' which are positioned on the ridge center and have the intrinsic orientation of the local ridge direction. A scheme is proposed to connect ridgels to smooth, spline-represented curves which fit the observed loops. Finally, a half-automated user interface allows one to merge or split, eliminate or select loop fits obtained form the above procedure. In this paper we apply our tool to one of the first EUV images observed by the SECCHI instrument onboard the recently launched STEREO spacecraft. We compare the extracted loops with projected field lines computed from almost-simultaneously-taken magnetograms measured by the SOHO/MDI Doppler imager. The field lines were calculated using a linear force-free field model. This comparison allows one to verify faint and spurious loop connections produced by our segmentation tool and it also helps to prove the quality of the magnetic-field model where well-identified loop structures comply with field-line projections. We also discuss further potential applications of our tool such as loop oscillations and stereoscopy.

Knowledge of the structure of the coronal magnetic field is important for our understanding of many solar activity phenomena, e.g. flares and CMEs. However, the direct measurement of coronal magnetic fields is not possible with present methods, and therefore the coronal field has to be extrapolated from photospheric measurements. Due to the low plasma beta the coronal magnetic field can usually be assumed to be approximately force free, with electric currents flowing along the magnetic field lines. There are both observational and theoretical reasons which suggest that at least prior to an eruption the coronal magnetic field is in a nonlinear force free state. Unfortunately the computation of nonlinear force free fields is way more difficult than potential or linear force free fields and analytic solutions are not generally available. We discuss several methods which have been proposed to compute nonlinear force free fields and focus particularly on an optimization method which has been suggested recently. We compare the numerical performance of a newly developed numerical code based on the optimization method with the performance of another code based on an MHD relaxation method if both codes are applied to the reconstruction of a semi-analytic nonlinear force-free solution. The optimization method has also been tested for cases where we add random noise to the perfect boundary conditions of the analytic solution, in this way mimicking the more realistic case where the boundary conditions are given by vector magnetogram data. We find that the convergence properties of the optimization method are affected by adding noise to the boundary data and we discuss possibilities to overcome this difficulty.

Recent observations of the solar corona with the LASCO coronagraph on board of the SOHO spacecraft have revealed the occurrence of triple helmet streamers even during solar minimum, which occasionally go unstable and give rise to large coronal mass ejections. There are also indications that the slow solar wind is either a combination of a quasi-stationary flow and a highly fluctuating component or may even be caused completely by many small eruptions or instabilities. As a first step we recently presented an analytical method to calculate simple two-dimensional stationary models of triple helmet streamer configurations. In the present contribution we use the equations of time- dependent resistive magnetohydrodynamics to investigate the stability and the dynamical behaviour of these configurations. We particularly focus on the possible differences between the dynamics of single isolated streamers and triple streamers and on the way in which magnetic reconnection initiates both small scale and large scale dynamical behaviour of the streamers. Our results indicate that small eruptions at the helmet streamer cusp may incessantly accelerate small amounts of plasma without significant changes of the equilibrium configuration and might thus contribute to the non-stationary slow solar wind. On larger time and length scales, large coronal eruptions can occur as a consequence of large scale magnetic reconnection events inside the streamer configuration. Our results also show that triple streamers are usually more stable than a single streamer.

Context: This paper presents a method which can be used to calculate models of the global solar corona from observational data. Aims: We present an optimization method for computing nonlinear magnetohydrostatic equilibria in spherical geometry with the aim to obtain self-consistent solutions for the coronal magnetic field, the coronal plasma density and plasma pressure using observational data as input. Methods: Our code for the self-consistent computation of the coronal magnetic fields and the coronal plasma solves the non-force-free magnetohydrostatic equilibria using an optimization method. Previous versions of the code have been used to compute non-linear force-free coronal magnetic fields from photospheric measurements in Cartesian and spherical geometry, and magnetostatic-equilibria in Cartesian geometry. We test our code with the help of a known analytic 3D equilibrium solution of the magnetohydrostatic equations. The detailed comparison between the numerical calculations and the exact equilibrium solutions is made by using magnetic field line plots, plots of density and pressure and some of the usual quantitative numerical comparison measures. Results: We find that the method reconstructs the equilibrium accurately, with residual forces of the order of the discretisation error of the analytic solution. The correlation with the reference solution is better than 99.9% and the magnetic energy is computed accurately with an error of <0.1%. Conclusions: We applied the method so far to an analytic test case. We are planning to use this method with real observational data as input as soon as possible.

The coronal magnetic field is an important quantity because the magnetic field dominates the structure of the solar corona. Unfortunately direct measurements of coronal magnetic fields are usually not available. The photospheric magnetic field is measured routinely with vector magnetographs. These photospheric measurements are extrapolated into the solar corona. The extrapolated coronal magnetic field depends on assumptions regarding the coronal plasma, e.g. force-freeness. Force-free means that all non-magnetic forces like pressure gradients and gravity are neglected. This approach is well justified in the solar corona due to the low plasma beta. One has to take care, however, about ambiguities, noise and non-magnetic forces in the photosphere, where the magnetic field vector is measured. Here we review different numerical methods for a nonlinear force-free coronal magnetic field extrapolation: Grad-Rubin codes, upward integration method, MHD-relaxation, optimization and the boundary element approach. We briefly discuss the main features of the different methods and concentrate mainly on recently developed new codes.

The solar magnetic field is key to understanding the physical processes in the solar atmosphere. Nonlinear force-free codes have been shown to be useful in extrapolating the coronal field from underlying vector boundary data [see Schrijver et al. 2006 for an overview]. However, we can only measure the magnetic field vector routinely with high accuracy in the photosphere with, e.g., Hinode/SOT, and unfortunately these data do not fulfill the force-free consistency condition as defined by Aly (1989). We must therefore apply some transformations to these data before nonlinear force-free extrapolation codes can be legitimately applied. To this end, we have developed a minimization procedure that uses the measured photospheric field vectors as input to approximate a more chromospheric like field The method was dubbed preprocessing. See Wiegelmann et al. 2006 for details]. The procedure includes force-free consistency integrals and spatial smoothing. The method has been intensively tested with model active regions [see Metcalf et al. 2008] and been applied to several ground based vector magnetogram data before. Here we apply the preprocessing program to photospheric magnetic field measurements with the Hinode/SOT instrument.

Our understanding of coronal phenomena, such as coronal plasma thermodynamics, faces a major handicap caused by missing coronal magnetic field measurements. Several lines in the UV wavelength range present suitable sensitivity to determine the coronal magnetic field via the Hanle effect. The latter is a largely unexplored diagnostic of coronal magnetic fields with a very high potential. Here we study the magnitude of the Hanle-effect signal to be expected outside the solar limb due to the Hanle effect in polarized radiation from the H {\sc{i}} Ly$\alpha$ and $\beta$ lines, which are among the brightest lines in the off-limb coronal FUV spectrum. For this purpose we use a magnetic field structure obtained by extrapolating the magnetic field starting from photospheric magnetograms. The diagnostic potential of these lines for determining the coronal magnetic field, as well as their limitations are studied. We show that these lines, in particular H {\sc{i}} Ly$\beta$, are useful for such measurements.

Solar flares and coronal mass ejections are associated with rapid changes in field connectivity and powered by the partial dissipation of electrical currents in the solar atmosphere. A critical unanswered question is whether the currents involved are induced by the motion of pre-existing atmospheric magnetic flux subject to surface plasma flows, or whether these currents are associated with the emergence of flux from within the solar convective zone. We address this problem by applying state-of-the-art nonlinear force-free field (NLFFF) modeling to the highest resolution and quality vector-magnetographic data observed by the recently launched Hinode satellite on NOAA Active Region 10930 around the time of a powerful X3.4 flare. We compute 14 NLFFF models with 4 different codes and a variety of boundary conditions. We find that the model fields differ markedly in geometry, energy content, and force-freeness. We discuss the relative merits of these models in a general critique of present abilities to model the coronal magnetic field based on surface vector field measurements. For our application in particular, we find a fair agreement of the best-fit model field with the observed coronal configuration, and argue (1) that strong electrical currents emerge together with magnetic flux preceding the flare, (2) that these currents are carried in an ensemble of thin strands, (3) that the global pattern of these currents and of field lines are compatible with a large-scale twisted flux rope topology, and (4) that the ~10^32 erg change in energy associated with the coronal electrical currents suffices to power the flare and its associated coronal mass ejection.

The space mission 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.

AIMS: We develop an optimization principle for computing stationary MHD equilibria. METHODS: Our code for the self-consistent computation of the coronal magnetic fields and the coronal plasma uses non-force-free MHD equilibria. Previous versions of the code have been used to compute non-linear force-free coronal magnetic fields from photospheric measurements. The program uses photospheric vector magnetograms and coronal EUV images as input. We tested our reconstruction code with the help of a semi-analytic MHD-equilibrium. The quality of the reconstruction was judged by comparing the exact and reconstructed solution qualitatively by magnetic field-line plots and EUV-images and quantitatively by several different numerical criteria. RESULTS: Our code is able to reconstruct the semi-analytic test equilibrium with high accuracy. The stationary MHD optimization code developed here has about the same accuracy as its predecessor, a non-linear force-free optimization code. The computing time for MHD-equilibria is, however, longer than for force-free magnetic fields. We also extended a well-known class of nonlinear force-free equilibria to the non-force-free regime for purposes of testing the code. CONCLUSIONS: We demonstrate that the code works in principle using tests with analytical equilibria, but it still needs to be applied to real data.

Knowledge regarding the coronal magnetic field is important for the understanding of many phenomena, like flares and coronal mass ejections. Because of the low plasma beta in the solar corona the coronal magnetic field is often assumed to be force-free and we use photospheric vector magnetograph data to extrapolate the magnetic field into the corona with the help of a non-linear force-free optimization code. Unfortunately the measurements of the photospheric magnetic field contain inconsistencies and noise. In particular the transversal components (say Bx and By) of current vector magnetographs have their uncertainties. Furthermore the magnetic field in the photosphere is not necessary force-free and often not consistent with the assumption of a force-free field above. We develop a preprocessing procedure to drive the observed non force-free data towards suitable boundary conditions for a force-free extrapolation. As a result we get a data set which is as close as possible to the measured data and consistent with the force-free assumption.

CONTEXT: As the coronal magnetic field can usually not be measured directly, it has to be extrapolated from photospheric measurements into the corona. AIMS: We test the quality of a non-linear force-free coronal magnetic field extrapolation code with the help of a known analytical solution. METHODS: The non-linear force-free equations are numerically solved with the help of an optimization principle. The method minimizes an integral over the force-free and solenoidal condition. As boundary condition we use either the magnetic field components on all six sides of the computational box in Case I or only on the bottom boundary in Case II. We check the quality of the reconstruction by computing how well force-freeness and divergence-freeness are fulfilled and by comparing the numerical solution with the analytical solution. The comparison is done with magnetic field line plots and several quantitative measures, like the vector correlation, Cauchy Schwarz, normalized vector error, mean vector error and magnetic energy. RESULTS: For Case I the reconstructed magnetic field shows good agreement with the original magnetic field topology, whereas in Case II there are considerable deviations from the exact solution. This is corroborated by the quantitative measures, which are significantly better for Case I. CONCLUSIONS: Despite the strong nonlinearity of the considered force-free equilibrium, the optimization method of extrapolation is able to reconstruct it; however, the quality of reconstruction depends significantly on the consistency of the input data, which is given only if the known solution is provided also at the lateral and top boundaries, and on the presence or absence of flux concentrations near the boundaries of the magnetogram.

We describe a newly developed code for the extrapolation of nonlinear force-free coronal magnetic fields in spherical coordinates. The program uses measured vector magnetograms on the solar photosphere as input and solves the force-free equations in the solar corona. The method is based on an optimization principle and the heritage of the newly developed code is a corresponding method in Cartesian geometry. We test the newly developed code with the help of a semi-analytic solution and rate the quality of our reconstruction qualitatively by magnetic field line plots and quantitatively with a number of comparison metrics. We find that we can reconstruct the original test field with high accuracy. The method is fast if the computation is limited to low co-latitudes (say $30^\circ \leq \theta \leq 150^\circ$), but becomes significantly slower if the polar regions are included.

Recent observations of the solar corona with the LASCO coronagraph on board of the SOHO spacecraft have revealed the occurrence of triple helmet streamers even during solar minimum, which occasionally go unstable and give rise to particularly huge coronal mass ejections. We present a method to calculate (semi-)analytically self-consistent stationary configurations of triple helmet streamers which can serve as input for stability considerations and dynamical calculations. The method is based on an asymptotic expansion procedure using the elongated structure of the streamers. The method is very flexible and can be used in both Cartesian and spherical geometry. We discuss the effects of magnetic shear, gravity and field-aligned flow on open field lines. Example solutions illustrating the influence of each of these features on the solution are presented.