Helmet Streamers with Triple Structures: Simulations of resistive dynamics

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.     

Nonlinear Force-Free Magnetic Field Extrapolations: Comparison of the Grad Rubin and Wheatland Sturrock Roumeliotis Algorithm

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.     

Estimating the Relative Helicity of Coronal Magnetic Fields

To quantify changes of the solar coronal field connectivity during eruptive events, one can use magnetic helicity, which is a measure of the shear or twist of a current-carrying (non-potential) field. To find a physically meaningful quantity, a relative measure, giving the helicity of a current-carrying field with respect to a reference (potential) field, is often evaluated. This requires a knowledge of the three-dimensional vector potential. We present a method to calculate the vector potential for a solenoidal magnetic field as the sum of a Laplacian part and a current-carrying part. The only requirements are the divergence freeness of the Laplacian and current-carrying magnetic field and the sameness of their normal field component on the bounding surface of the considered volume.     

Computer Vision for the Solar Dynamics Observatory (SDO)

In Fall 2008 NASA selected a large international consortium to produce a comprehensive automated feature-recognition system for the Solar Dynamics Observatory (SDO). The SDO data that we consider are all of the Atmospheric Imaging Assembly (AIA) images plus surface magnetic-field images from the Helioseismic and Magnetic Imager (HMI). We produce robust, very efficient, professionally coded software modules that can keep up with the SDO data stream and detect, trace, and analyze numerous phenomena, including flares, sigmoids, filaments, coronal dimmings, polarity inversion lines, sunspots, X-ray bright points, active regions, coronal holes, EIT waves, coronal mass ejections (CMEs), coronal oscillations, and jets. We also track the emergence and evolution of magnetic elements down to the smallest detectable features and will provide at least four full-disk, nonlinear, force-free magnetic field extrapolations per day. The detection of CMEs and filaments is accomplished with Solar and Heliospheric Observatory (SOHO)/Large Angle and Spectrometric Coronagraph (LASCO) and ground-based Hα data, respectively. A?completely new software element is a trainable feature-detection module based on a generalized image-classification algorithm. Such a trainable module can be used to find features that have not yet been discovered (as, for example, sigmoids were in the pre-Yohkoh era). Our codes will produce entries in the Heliophysics Events Knowledgebase (HEK) as well as produce complete catalogs for results that are too numerous for inclusion in the HEK, such as the X-ray bright-point metadata. This will permit users to locate data on individual events as well as carry out statistical studies on large numbers of events, using the interface provided by the Virtual Solar Observatory. The operations concept for our computer vision system is that the data will be analyzed in near real time as soon as they arrive at the SDO Joint Science Operations Center and have undergone basic processing. This will allow the system to produce timely space-weather alerts and to guide the selection and production of quicklook images and movies, in addition to its prime mission of enabling solar science. We briefly describe the complex and unique data-processing pipeline, consisting of the hardware and control software required to handle the SDO data stream and accommodate the computer-vision modules, which has been set up at the Lockheed-Martin Space Astrophysics Laboratory (LMSAL), with an identical copy at the Smithsonian Astrophysical Observatory (SAO).     

Global Solar Free Magnetic Energy and Electric Current Density Distribution of Carrington Rotation 2124

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 to only 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 nonlinear force-free field (NLFFF) extrapolation codes implemented to a synoptic maps of photospheric vector magnetic field synthesized from the 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 magnetic energy throughout the solar atmosphere for Carrington rotation 2124 using our global extrapolation code. We found that spatially, the low-lying, current-carrying core field demonstrates a strong concentration of free energy in the active-region core, from the photosphere to the lower corona (about 70 Mm). The free energy density appears largely co-spatial with the electric current distribution.     

How Should One Optimize Nonlinear Force-Free Coronal Magnetic Field Extrapolations from SDO/HMI Vector Magnetograms?

The Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) provides 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 semianalytic and numeric equilibria and applied to vector magnetograms from Hinode and ground-based observations. Recently we implemented a new version which takes into account measurement errors in photospheric vector magnetograms. Photospheric field measurements are often affected by measurement errors and finite nonmagnetic forces inconsistent for use as a boundary for a force-free field in the corona. To deal with these uncertainties, we developed two improvements: i) preprocessing of the surface measurements to make them compatible with a force-free field, and ii) new code which keeps a balance between the force-free constraint and deviation from the photospheric field measurements. Both methods contain free parameters, which must be optimized for use with data from SDO/HMI. In this work we describe the corresponding analysis method and evaluate the force-free equilibria by how well force-freeness and solenoidal conditions are fulfilled, by the angle between magnetic field and electric current, and by comparing projections of magnetic field lines with coronal images from the Atmospheric Imaging Assembly (SDO/AIA). We also compute the available free magnetic energy and discuss the potential influence of control parameters.     

Segmentation of Loops from Coronal EUV Images

We present a procedure to extract 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; these “ridgels” 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 that fit the observed loops. Finally, a half-automated user interface allows one to merge or split curves or eliminate or select loop fits obtained from this 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 near-simultaneous magnetograms measured by the SOHO/MDI Doppler imager. The field lines were calculated by 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.     

Evolution of a Coronal Loop System

The temporal variation of a loop system that appears to be changing rapidly is examined. The analyzed data were obtained on 15 May 1999, with the Transition Region and Coronal Explorer (TRACE) during an observing campaign and consist of observations in the Fe ix/Fe x 171?Å and Fe xii 195?Å passbands taken at a cadence of ～10 min. The special interest in this loop system is that it looks like one expanding loop; however, careful examination reveals that the loop consists of several strands and that new loop strands become visible successively at higher altitudes and lower loop strands fade out during the one hour of our observations. These strands have different widths, densities, and temperatures and are most probably consisting of, at least, a few unresolved thinner threads. Several geometric and physical parameters are derived for two of the strands and an effort is made to determine their 3D structure based on the extrapolation of the magnetic field lines. Electron density estimates allow us to derive radiative and conductive cooling times and to conclude that these loop strands are cooling by radiation.