sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Inhester, B.

Normalized to: Inhester, B.

45 article(s) in total. 114 co-authors, from 1 to 24 common article(s). Median position in authors list is 3,0.

[1]  oai:arXiv.org:2004.03790  [pdf] - 2097397
Initiation and Early Kinematic Evolution of Solar Eruptions
Comments: Accepted for publication in ApJ; 24 pages, 12 figures, 3 tables
Submitted: 2020-04-07
We investigate the initiation and early evolution of 12 solar eruptions, including six active region hot channel and six quiescent filament eruptions, which were well observed by the \textsl{Solar Dynamics Observatory}, as well as by the \textsl{Solar TErrestrial RElations Observatory} for the latter. The sample includes one failed eruption and 11 coronal mass ejections, with velocities ranging from 493 to 2140~km~s$^{-1}$. A detailed analysis of the eruption kinematics yields the following main results. (1) The early evolution of all events consists of a slow-rise phase followed by a main-acceleration phase, the height-time profiles of which differ markedly and can be best fit, respectively, by a linear and an exponential function. This indicates that different physical processes dominate in these phases, which is at variance with models that involve a single process. (2) The kinematic evolution of the eruptions tends to be synchronized with the flare light curve in both phases. The synchronization is often but not always close. A delayed onset of the impulsive flare phase is found in the majority of the filament eruptions (5 out of 6). This delay, and its trend to be larger for slower eruptions, favor ideal MHD instability models. (3) The average decay index at the onset heights of the main acceleration is close to the threshold of the torus instability for both groups of events (although based on a tentative coronal field model for the hot channels), suggesting that this instability initiates and possibly drives the main acceleration.
[2]  oai:arXiv.org:1902.06953  [pdf] - 2124495
Electron Density Reconstruction of Solar Coronal Mass Ejections Based on Genetic Algorithm: Method and Application
Comments: 30 pages, 23 figures, 1 table. Submitted to ApJ, making the second revision
Submitted: 2019-02-19, last modified: 2019-10-09
We present a new method to reconstruct three dimensional electron density of Coronal Mass Ejection (CME) based on genetic algorithm, namely genetic reconstruction method (GRM). GRM is firstly applied to the model CMEs with different orientations and shapes. A set of analytic GL98 model CMEs are employed to produce synthetic CME images for GRM reconstruction. Model CMEs with longitude of 0, 45, 90, 135, 180 degrees and latitude of 0, 15, 30, 45 degrees are used to test performance of GRM. We also obscure the CME with simulated occulter of coronagraph to find out the influence of incompleteness of CME brightness. GRM reconstruction is carried out using data from STEREO A and B with separation angle of 90 degrees and from STEREO A and SOHO with separation angle of 73 degrees, respectively. Pearson correlation coefficient (pcc) is calculated to analyze similarity between model and reconstructed CMEs for brightness and electron density. Comparisons based on the similarity analysis under various conditions stated above give us valuable insights of advantages and limitations of GRM reconstruction. Then the method is applied to real coronagraph data from STEREO-A, B and SOHO on September 30th, 2013.
[3]  oai:arXiv.org:1810.13003  [pdf] - 1868002
On the error analyses of polarization measurements of the white-light coronagraph aboard ASO-S
Comments: 16 pages, 5 figure, Research in Astronomy and Astrophysics, accepted
Submitted: 2018-10-29
The Advanced Space-based Solar Observatory (ASO-S) mission aims to explore two most spectacular eruptions in the Sun: solar flares and coronal mass ejections (CMEs), and their magnetism. For the studies of CMEs, the payload Lyman-alpha Solar Telescope (LST) has been proposed. It includes a traditional white-light coronagraph and a Lyman-alpha coronagraph which opens a new window to CME observations. Polarization measurements taken by white-light coronagraphs are crucial to derive fundamental physical parameters of CMEs. To make such measurements, there are two options of Stokes polarimeter which have been used by existing white-light coronagraphs for space missions. One uses a single or triple linear polarizers, the other involves both a half-wave plate and a linear polarizer. We find that the former option subjects to less uncertainty in the derived Stokes vector propagated from detector noise. The latter option involves two plates which are prone to internal reflections and may have a reduced transmission factor. Therefore, the former option is adopted as our Stokes polarimeter scheme for LST. Based on the parameters of the intended linear polarizer(s) colorPol provided by CODIXX and the half-wave plate 2-APW-L2-012C by Altechna, it is further shown that the imperfect maximum transmittance of the polarizer significantly increases the variance amplification of Stokes vector by at least about 50% when compared with the ideal case. The relative errors of Stokes vector caused by the imperfection of colorPol polarizer and the uncertainty due to the polarizer assembling in the telescope are estimated to be about 5%. Among the considered parameters, we find that the dominant error comes from the uncertainty in the maximum transmittance of the polarizer.
[4]  oai:arXiv.org:1709.04177  [pdf] - 1588314
Nonlinear force-free coronal magnetic stereoscopy
Comments: 9 pages, 7 figures
Submitted: 2017-09-13
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.
[5]  oai:arXiv.org:1612.09360  [pdf] - 1533951
Measure the Propagation of a halo CME and Its Driven Shock with the Observations from a Single Perspective at Earth
Comments: 13 pages, 12 figures, 2 tables, accepted by APJ
Submitted: 2016-12-29
We present a detailed study of an earth-directed coronal mass ejection (Full halo CME) event happened on 2011 February 15 making use of white light observations by three coronagraphs and radio observations by Wind/WAVES. We applied three different methods to reconstruct the propagation direction and traveling distance of the CME and its driven shock. We measured the kinematics of the CME leading edge from white light images observed by STEREO A and B, tracked the CME-driven shock using the frequency drift observed by Wind/WAVES together with an interplanetary density model, and obtained the equivalent scattering centers of the CME by Polarization Ratio(PR) method. For the first time, we applied PR method to different features distinguished from LASCO/C2 polarimetric observations and calculated their projections onto white light images observed by STEREO A and B. By combining the graduated cylindrical shell (GCS) forward modeling with the PR method, we proposed a new GCS-PR method to derive 3D parameters of a CME observed from a single perspective at Earth. Comparisons between different methods show a good degree of consistence in the derived 3D results.
[6]  oai:arXiv.org:1607.05533  [pdf] - 1498353
Comparison of CME/shock propagation models with heliospheric imaging and in situ observations
Comments: accepted for publication in ApJS
Submitted: 2016-07-19
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.
[7]  oai:arXiv.org:1512.00651  [pdf] - 1385425
Thomson Scattering in the Solar Corona
Comments: 104 pages, 31 figures
Submitted: 2015-12-02, last modified: 2016-04-05
The basis for the application of Thomson scattering to the analysis of coronagraph images has been laid decades ago.Even though the basic formulation is undebated, a discussion has grown in recent years about the spatial distribution of Thomson scatter sensitivity in the corona and the inner heliosphere. These notes are an attempt to clarify the understanding of this topic. We reformulate the classical scattering calculations in a more transparent way using modern SI-compatible quantities extended to field correlation matrices. The resulting concise formulation is easily extended to the case of relativistic electrons. For relativistic electrons we calculate the Stokes parameters of the scattered radiation and determine changes in degree and orientation of its polarisation, blue-shift and radiant intensities depending on the electron velocity magnitude and direction. We discuss the probability to see these relativistic effects in white-light coronagraph observations of the solar corona. Many mathematical and some basic physical ingredients are made explicit in several chapters of the appendix.
[8]  oai:arXiv.org:1603.06379  [pdf] - 1439777
Analysis of UV and EUV emission from impacts on the Sun after 2011 June 7 eruptive flare
Comments: 10 pages, 14 figures To be published in A&A
Submitted: 2016-03-21
On 2011 June 7 debris from a large filament eruption fell back to the Sun causing bright ultraviolet (UV) and extreme ultraviolet (EUV) splashes across the surface. These impacts may give clues on the process of stellar accretion. The aim is to investigate how the impact emission is influenced by structures in the falling ejecta and at the solar surface. We determine the UV and EUV light curves of a sample of impacts. The ballistic impact velocity is estimated from the ejection and landing times and, where possible, compared with the velocity derived by tracking the downflows in SDO/AIA and STEREO/EUVI images. Estimates of the column density before impact are made from the darkness of the falling plasma in the 193 A channel. The impact velocities were between 230 and 450 km/s. All impacts produced bright EUV emission at the impact site but bright UV was only observed when the impacting fragments reached the chromosphere. There was no clear relation between EUV intensity and kinetic energy. Low UV to EUV intensity ratios (I{UV}/I{EUV}) were seen (i) from impacts of low column-density fragments, (ii) when splashes, produced by some impacts, prevented subsequent fragments from reaching the chromosphere, and (iii) from an impact in an active region. The earliest impacts with the lowest velocity (~250 km/s) had the highest I{UV}/I{EUV}. The I{UV}/I{EUV} decreases with impact velocity, magnetic field at the impact site, and EUV ionising flux. Many of the infalling fragments dissipate above the chromosphere either due to ionisation and trapping in magnetic structures, or to them encountering a splash from an earlier impact. If the same happens in accreting stars then the reduced X-ray compared to optical emission that has been observed is more likely due to absorption by the trailing stream than locally at the impact site.
[9]  oai:arXiv.org:1509.02246  [pdf] - 1292340
Why does the apparent mass of a coronal mass ejection increase?
Comments: 27 pages, 2 tables, 9 figures, accepted by ApJ
Submitted: 2015-09-07
Mass is one of the most fundamental parameters characterizing the dynamics of a coronal mass ejection (CME). It has been found that CME apparent mass measured from the brightness enhancement in coronagraph images shows an increasing trend during its evolution in the corona. However, the physics behind it is not clear. Does the apparent mass gain come from the mass outflow from the dimming regions in the low corona, or from the pileup of the solar wind plasma around the CME when it propagates outwards from the Sun? We analyzed the mass evolution of six CME events. Their mass can increase by a factor of 1.6 to 3.2 from 4 to 15 Rs in the field of view (FOV) of the coronagraph on board the Solar Terrestrial Relations Observatory (STEREO). Over the distance about 7 to 15 Rs, where the coronagraph occulting effect can be negligible, the mass can increase by a factor of 1.3 to 1.7. We adopted the `snow-plough' model to calculate the mass contribution of the piled-up solar wind in the height range from about 7 to 15 Rs. For 2/3 of the events, the solar wind pileup is not sufficient to explain the measured mass increase. In the height range from about 7 to 15 Rs, the ratio of the modeled to the measured mass increase is roughly larger than 0.55. Although the ratios are believed to be overestimated, the result gives evidence that the solar wind pileup probably makes a non-negligible contribution to the mass increase. It is not clear yet whether the solar wind pileup is a major contributor to the final mass derived from coronagraph observations. However, our study suggests that the solar wind pileup plays increasingly important role in the mass increase as a CME moves further away from the Sun.
[10]  oai:arXiv.org:1503.08502  [pdf] - 972677
Radial Flow Pattern of a Slow CME
Comments: ApJ, accepted
Submitted: 2015-03-29
Height-time plots of the leading edge of coronal mass ejections (CME) have often been used to study CME kinematics. We propose a new method to analyze the CME kinematics in more detail by determining the radial mass transport process throughout the entire CME. Thus our method is able to estimate not only the speed of the CME front but also the radial flow speed inside the CME. We have applied the method to a slow CME with an average leading edge speed about 480 km s$^{-1}$. In the Lagrangian frame, the speed of the individual CME mass elements stay almost constant within 2 and 15 R$_S$, the range over which we analyzed the CME. Hence we have no evidence of net radial forces acting on parts of the CME in this range nor of a pile-up of mass ahead of the CME. We find evidence that the leading edge trajectory obtained by tie-pointing may gradually lag behind the Lagrangian front-side trajectories derived from our analysis. Our results also allow a much more precise estimate of the CME energy. Compared with conventional estimates using the CME total mass and leading-edge motion, we find that the latter may overestimate the kinetic energy and the gravitational potential energy.
[11]  oai:arXiv.org:1409.3896  [pdf] - 867191
SDO/AIA Observations of a Reflecting Longitudinal Wave in a Coronal Loop
Comments: 8 pages, 5 figures
Submitted: 2014-09-12
We report high resolution observations from the SDO/AIA of intensity oscillations in a hot, T~8-10 MK, loop. The AIA images show a large coronal loop that was rapidly heated following plasma ejection from one of the loop's footpoints. A wave-like intensity enhancement, seen very clearly in the 131 and 94 \AA\ channel images, propagated ahead of the ejecta along the loop, and was reflected at the opposite footpoint. The wave reflected four times before fading. It was only seen in the hot, 131 and 94 \AA\ channels. The characteristic period and the decay time of the oscillation was $\sim$630 and $\sim$440 s, respectively. The phase speed was about 460-510 km/s which roughly matches the sound speed of the loop (430-480 km/s). The observed properties of the oscillation are consistent with the observations of Doppler shift oscillations discovered by the SOHO/SUMER and with their interpretation as slow magnetoacoustic waves. We suggest that the impulsive injection of plasma, following reconnection at one of the loop footpoints, led to rapid heating and the propagation of a longitudinal compressive wave along the loop. The wave bounces back and forth a couple of times before fading.
[12]  oai:arXiv.org:1309.5916  [pdf] - 1179429
Polar plumes' orientation and the Sun's global magnetic field
Comments:
Submitted: 2013-09-20
We characterize the orientation of polar plumes as a tracer of the large-scale coronal magnetic field configuration. We monitor in particular the north and south magnetic pole locations and the magnetic opening during 2007-2008 and provide some understanding of the variations in these quantities. The polar plume orientation is determined by applying the Hough-wavelet transform to a series of EUV images and extracting the key Hough space parameters of the resulting maps. The same procedure is applied to the polar cap field inclination derived from extrapolating magnetograms generated by a surface flux transport model. We observe that the position where the magnetic field is radial (the Sun's magnetic poles) reflects the global organization of magnetic field on the solar surface, and we suggest that this opens the possibility of both detecting flux emergence anywhere on the solar surface (including the far side) and better constraining the reorganization of the corona after flux emergence.
[13]  oai:arXiv.org:1307.5120  [pdf] - 1172846
Kelvin-Helmholtz instability of a coronal streamer
Comments: ApJ, accepted
Submitted: 2013-07-18
The shear-flow-driven instability can play an important role in energy transfer processes in coronal plasma. We present for the first time the observation of a kink-like oscillation of a streamer probably caused by the streaming kinkmode Kelvin-Helmholtz instability. The wave-like behavior of the streamer was observed by Large Angle and Spectrometric Coronagraph Experiment (LASCO) C2 and C3 aboard SOlar and Heliospheric Observatory (SOHO). The observed wave had a period of about 70 to 80 minutes, and its wavelength increased from 2 Rsun to 3 Rsun in about 1.5 hours. The phase speeds of its crests and troughs decreased from 406 \pm 20 to 356 \pm 31kms^{-1} during the event. Within the same heliocentric range, the wave amplitude also appeared to increase with time. We attribute the phenomena to the MHD Kelvin-Helmholtz instability which occur at a neutral sheet in a fluid wake. The free energy driving the instability is supplied by the sheared flow and sheared magnetic field across the streamer plane. The plasma properties of the local environment of the streamer were estimated from the phase speed and instability threshold criteria.
[14]  oai:arXiv.org:1212.5639  [pdf] - 1158673
A Comparison Between Nonlinear Force-Free Field and Potential Field Models Using Full-Disk SDO/HMI Magnetogram
Comments: 5 figures, Submitted to Solar Physics Journal
Submitted: 2012-12-21, last modified: 2012-12-28
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.
[15]  oai:arXiv.org:1210.3668  [pdf] - 1152097
Full-disk nonlinear force-free field extrapolation of SDO/HMI and SOLIS/VSM magnetograms
Comments: 5 figures, Accepted by Astronomy and Astrophysics Journal
Submitted: 2012-10-12
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.
[16]  oai:arXiv.org:1207.2656  [pdf] - 1214918
Comparisons of CME morphological characteristics derived from five 3D reconstruction methods
Comments: to appear in Solar Physics
Submitted: 2012-07-11
We compare different methods to reconstruct the three-dimensional (3D) CME morphology. The explored methods include geometric localisation, mask fitting, forward modeling, polarisation ratio and local correlation tracking plus triangulation. The five methods are applied to the same CME event, which occurred on August 7 2010. Their corresponding results are presented and compared, especially in their propagation direction and spatial extent in 3D. We find that mask fitting and geometric localisation method produce consistent results. Reconstructions including three-view observations are more precise than reconstructions done with only two views. Compared to the forward modeling method, in which a-priori shape of the CME geometry is assumed, mask fitting has more flexibility. Polarisation ratio method makes use of the Thomson scattering geometry. We find spatially the 3D CME derived from mask fitting lies mostly in the overlap region obtained with the polarisation method from COR2 A and B. In addition, mask fitting can help resolve the front/back ambiguity inherent in the polarisation ratio method. However, local correlation tracking plus triangulation did not show a consistent result with the other four methods. For reconstructions of a diffuse CME, when the separation angle between STEREO A and B is large, finding two corresponding points in a STEREO image pair becomes very difficult. Excluding the local correlation tracking method, the latitude of the CME's centre of gravity derived from the other methods deviates within one degree and longitude differs within 19 degrees.
[17]  oai:arXiv.org:1203.6732  [pdf] - 513418
Study of a Prominence Eruption using PROBA2/SWAP and STEREO/EUVI Data
Comments: 15 Pages, Accepted for publication in Solar Physics, PROBA2 Topical Issue
Submitted: 2012-03-30, last modified: 2012-05-16
Observations of the early rise and propagation phases of solar eruptive prominences can provide clues about the forces acting on them through the behavior of their acceleration with height. We have analyzed such an event, observed on 13 April 2010 by SWAP on PROBA2 and EUVI on STEREO. A feature at the top of the erupting prominence was identified and tracked in images from the three spacecraft. The triangulation technique was used to derive the true direction of propagation of this feature. The reconstructed points were fitted with two mathematical models: i) a power-law polynomial function and ii) a cubic smoothing spline, in order to derive the accelerations. The first model is characterized by five degrees of freedom while the second one is characterized by ten degrees of freedom. The results show that the acceleration increases smoothly and it is continuously increasing with height. We conclude that the prominence is not accelerated immediately by local reconnection but rather is swept away as part of a large-scale relaxation of the coronal magnetic field.
[18]  oai:arXiv.org:1203.3261  [pdf] - 488628
Morphological evolution of a 3D CME cloud reconstructed from three viewpoints
Comments: ApJ, accepted
Submitted: 2012-03-14
The propagation properties of coronal mass ejections (CMEs) are crucial to predict its geomagnetic effect. A newly developed three dimensional (3D) mask fitting reconstruction method using coronagraph images from three viewpoints has been described and applied to the CME ejected on August 7, 2010. The CME's 3D localisation, real shape and morphological evolution are presented. Due to its interaction with the ambient solar wind, the morphology of this CME changed significantly in the early phase of evolution. Two hours after its initiation, it was expanding almost self-similarly. CME's 3D localisation is quite helpful to link remote sensing observations to in situ measurements. The investigated CME was propagating to Venus with its flank just touching STEREO B. Its corresponding ICME in the interplanetary space shows a possible signature of a magnetic cloud with a preceding shock in VEX observations, while from STEREO B only a shock is observed. We have calculated three principle axes for the reconstructed 3D CME cloud. The orientation of the major axis is in general consistent with the orientation of a filament (polarity inversion line) observed by SDO/AIA and SDO/HMI. The flux rope axis derived by the MVA analysis from VEX indicates a radial-directed axis orientation. It might be that locally only the leg of the flux rope passed through VEX. The height and speed profiles from the Sun to Venus are obtained. We find that the CME speed possibly had been adjusted to the speed of the ambient solar wind flow after leaving COR2 field of view and before arriving Venus. A southward deflection of the CME from the source region is found from the trajectory of the CME geometric center. We attribute it to the influence of the coronal hole where the fast solar wind emanated from.
[19]  oai:arXiv.org:1109.4889  [pdf] - 1084318
Coronal Magnetic Field Structure and Evolution for Flaring AR 11117 and its Surroundings
Comments: 14 pages, 14 figures, Accepted by Solar Physics Journal
Submitted: 2011-09-22, last modified: 2012-02-28
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.
[20]  oai:arXiv.org:1202.4981  [pdf] - 1116832
Break up of returning plasma after the 7 June 2011 filament eruption by Rayleigh-Taylor instabilities
Comments: 6 pages, 7 figures, submitted A&A. Movies are at http://www.mps.mpg.de/data/outgoing/innes/arcs_movie.avi and http://www.mps.mpg.de/data/outgoing/innes/horns_movie.avi
Submitted: 2012-02-22
A prominence eruption on 7 June 2011 produced spectacular curtains of plasma falling through the lower corona. At the solar surface they created an incredible display of extreme ultraviolet brightenings. The aim is to identify and analyze some of the local instabilities which produce structure in the falling plasma. The structures were investigated using SDO/AIA 171A and 193A images in which the falling plasma appeared dark against the bright coronal emission. Several instances of the Rayleigh-Taylor instability were investigated. In two cases the Alfven velocity associated with the dense plasma could be estimated from the separation of the Rayleigh-Taylor fingers. A second type of feature, which has the appearance of self-similar branching horns, is also discussed.
[21]  oai:arXiv.org:1202.3601  [pdf] - 1116692
How to optimize nonlinear force-free coronal magnetic field extrapolations from SDO/HMI vector magnetograms?
Comments: 17 Pages, 6 Figures, Sol. Phys., accepted
Submitted: 2012-02-16
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.
[22]  oai:arXiv.org:1112.2255  [pdf] - 450193
Particle kinetic analysis of a polar jet from SECCHI COR data
Comments: A&A, in press
Submitted: 2011-12-10
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.
[23]  oai:arXiv.org:1111.4123  [pdf] - 1091735
2D and 3D Polar Plume Analysis from the Three Vantage Positions of STEREO/EUVI A, B, and SOHO/EIT
Comments: 22 pages, 10 figures, Solar Physics article
Submitted: 2011-11-17
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.
[24]  oai:arXiv.org:1108.5304  [pdf] - 1083655
Solar magnetism eXplorer (SolmeX)
Comments: Submitted to Experimental Astronomy. 30 pages, 19 figures
Submitted: 2011-08-26
The magnetic field plays a pivotal role in many fields of Astrophysics. This is especially true for the physics of the solar atmosphere. Measuring the magnetic field in the upper solar atmosphere is crucial to understand the nature of the underlying physical processes that drive the violent dynamics of the solar corona -- that can also affect life on Earth. SolmeX, a fully equipped solar space observatory for remote-sensing observations, will provide the first comprehensive measurements of the strength and direction of the magnetic field in the upper solar atmosphere. The mission consists of two spacecraft, one carrying the instruments, and another one in formation flight at a distance of about 200m carrying the occulter to provide an artificial total solar eclipse. This will ensure high-quality coronagraphic observations above the solar limb. Solmex integrates two spectro-polarimetric coronagraphs for off-limb observations, one in the EUV and one in the IR, and three instruments for observations on the disk. The latter comprises one imaging polarimeter in the EUV for coronal studies, a spectro-polarimeter in the EUV to investigate the low corona, and an imaging spectro-polarimeter in the UV for chromospheric studies. SOHO and other existing missions have investigated the emission of the upper atmosphere in detail (not considering polarization), and as this will be the case also for missions planned for the near future. Therefore it is timely that SolmeX provides the final piece of the observational quest by measuring the magnetic field in the upper atmosphere through polarimetric observations.
[25]  oai:arXiv.org:1105.3391  [pdf] - 359976
Low polarized emission from the core of coronal mass ejections
Comments: 5 pages, 3 figures
Submitted: 2011-05-17
In white-light coronagraph images, cool prominence material is sometimes observed as bright patches in the core of coronal mass ejections (CMEs). If, as generally assumed, this emission is caused by Thomson-scattered light from the solar surface, it should be strongly polarised tangentially to the solar limb. However, the observations of a CME made with the SECCHI/STEREO coronagraphs on 31 August 2007 show that the emission from these bright core patches is exceptionally low polarised. We used the polarisation ratio method of Moran and Davila (2004) to localise the barycentre of the CME cloud. By analysing the data from both STEREO spacecraft we could resolve the plane-of-the-sky ambiguity this method usually suffers from. Stereoscopic triangulation was used to independently localise the low-polarisation patch relative to the cloud. We demonstrated for the first time that the bright core material is located close to the centre of the CME cloud. We show that the major part of the CME core emission, more than 85% in our case, is H$\alpha$ radiation and only a small fraction is Thomson-scattered light. Recent calculations also imply that the plasma density in the patch is 8 10$^8$ cm$^{-3}$ or more compared to 2.6 10$^6$ cm$^{-3}$ for the Thomson-scattering CME environment surrounding the core material.
[26]  oai:arXiv.org:1104.2246  [pdf] - 342958
Magnetic Connectivity between Active Regions 10987, 10988, and 10989 by Means of Nonlinear Force-Free Field Extrapolation
Comments: Solar Physics
Submitted: 2011-04-12
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.
[27]  oai:arXiv.org:1011.6285  [pdf] - 1042249
Nonlinear force-free field extrapolation in spherical geometry: improved boundary data treatment applied to a SOLIS/VSM vector magnetogram
Comments: Accepted on Astronomy and Astrophysics
Submitted: 2010-11-26
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.
[28]  oai:arXiv.org:1008.0771  [pdf] - 210790
The solar magnetic field
Comments: 83 pages, 42 figures, about 2.5 MB in size
Submitted: 2010-08-04
The magnetic field of the Sun is the underlying cause of the many diverse phenomena combined under the heading of solar activity. Here we describe the magnetic field as it threads its way from the bottom of the convection zone, where it is built up by the solar dynamo, to the solar surface, where it manifests itself in the form of sunspots and faculae, and beyond into the outer solar atmosphere and, finally, into the heliosphere. On the way it, transports energy from the surface and the subsurface layers into the solar corona, where it heats the gas and accelerates the solar wind.
[29]  oai:arXiv.org:1001.1924  [pdf] - 32690
Observations of a rotating macrospicule associated with an X-ray jet
Comments: 4 pages, 3 figures, accepted for publication in A&A
Submitted: 2010-01-12
We attempt to understand the driving mechanism of a macrospicule and its relationship with a coronal jet. We study the dynamics of a macrospicule and an associated coronal jet captured by multi-spacecraft observations. Doppler velocities both in the macrospicule and the coronal jet are determined by EIS and SUMER spectra. Their temporal evolution is studied using X-ray and He II 304 images. A blueshift of -120+/-15 km/s is detected on one side of the macrospicule, while a redshift of 50+/-6 km/s is found at the base of the other side. The inclination angle of the macrospicule inferred from a stereoscopic analysis with STEREO suggests that the measured Doppler velocities can be attributed to a rotating motion of the macrospicule rather than a radial flow or an expansion. The macrospicule is driven by the unfolding motion of a twisted magnetic flux rope, while the associated X-ray jet is a radial outflow.
[30]  oai:arXiv.org:0912.1514  [pdf] - 902675
Nonlinear force-free coronal magnetic field modelling and preprocessing of vector magnetograms in spherical geometry
Comments: 12 pages, 5 figures
Submitted: 2009-12-08
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.
[31]  oai:arXiv.org:0912.1267  [pdf] - 1018657
Solar stereoscopy - where are we and what developments do we require to progress?
Comments: 12 Pages, 9 Figures, a Review article
Submitted: 2009-12-07
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.
[32]  oai:arXiv.org:0908.2365  [pdf] - 27361
Stereoscopic Polar Plume Reconstructions from Stereo/Secchi Images
Comments: 25 pages, 13 figures
Submitted: 2009-08-17
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.
[33]  oai:arXiv.org:0902.1007  [pdf] - 315587
A Critical Assessment of Nonlinear Force-Free Field Modeling of the Solar Corona for Active Region 10953
Comments: Accepted to ApJ
Submitted: 2009-02-05
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.
[34]  oai:arXiv.org:0802.0773  [pdf] - 9808
First Stereoscopic Coronal Loop Reconstructions from Stereo Secchi Images
Comments: 6 pages, 5 figures
Submitted: 2008-02-06
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.
[35]  oai:arXiv.org:0802.0687  [pdf] - 9790
Magnetic modelling and tomography: First steps towards a consistent reconstruction of the solar corona
Comments: 29 pages, 6 figures
Submitted: 2008-02-05
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.
[36]  oai:arXiv.org:0801.4573  [pdf] - 9543
How to use magnetic field information for coronal loop identification?
Comments: 14 pages, 3 figures
Submitted: 2008-01-30
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.
[37]  oai:arXiv.org:0801.4519  [pdf] - 9521
Comparing magnetic field extrapolations with measurements of magnetic loops
Comments: 5 pages, 3 figures
Submitted: 2008-01-29
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.
[38]  oai:arXiv.org:0801.3237  [pdf] - 9231
Nonlinear force-free magnetic field extrapolations: comparison of the Grad-Rubin and Wheatland-Sturrock-Roumeliotis algorithm
Comments: 21 pages, 13 figures
Submitted: 2008-01-21, last modified: 2008-01-23
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.
[39]  oai:arXiv.org:0801.3240  [pdf] - 9232
Segmentation of Loops from Coronal EUV Images
Comments: 13 pages, 9 figures, Solar Physics, online first
Submitted: 2008-01-21
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.
[40]  oai:arXiv.org:0801.2916  [pdf] - 9164
Optimization approach for the computation of magnetohydrostatic coronal equilibria in spherical geometry
Comments: 6 pages, 4 figures
Submitted: 2008-01-18
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.
[41]  oai:arXiv.org:astro-ph/0612636  [pdf] - 87974
Magnetic Stereoscopy
Comments: 18 pages, 5 figures
Submitted: 2006-12-21
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.
[42]  oai:arXiv.org:astro-ph/0612639  [pdf] - 87977
Inversion of coronal Zeeman and Hanle Observations to reconstruct the coronal magnetic field
Comments: 6 pages, 6 figures
Submitted: 2006-12-21
Hanle-effect observations of forbidden coronal line transitions and recently also longitudinal Zeeman-effect measurements of coronal lines show quantitative signatures of the weak coronal magnetic field. The interpretation of these observations is, however, complicated by the fact that they are the result of line-of-sight integrations through the optically thin corona. We study by means of simulated observations the possibility of applying tomographic techniques in order to reconstruct the 3D magnetic field configuration in the solar corona from these observations. The reconstruction problem relates to a family of similar problems termed vector tomography. It is shown that Zeeman data and Hanle data alone obtained from vantage points in the ecliptic plane alone are sensitive only to certain magnetic field structures. For a full reconstruction it is necessary to combine the longitudinal Zeeman and Hanle effect data.
[43]  oai:arXiv.org:astro-ph/0612641  [pdf] - 87979
Preprocessing of vector magnetograph data for a non-linear force-free magnetic field reconstruction
Comments: 20 pages, 5 figures
Submitted: 2006-12-21
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.
[44]  oai:arXiv.org:astro-ph/0612649  [pdf] - 87987
Stereoscopy basics for the STEREO mission
Comments: 28 pages, 43 figures
Submitted: 2006-12-21
We discuss some basic principles of stereoscopy and their relevance to the reconstruction of coronal loops. The aim of the paper is to make the solar physicist familiar with basic stereoscopy principles and to give hints how they may apply to the analysis of data from the forthcoming STEREO mission. We disucss the geometry of the solar coronal stereo problem, give the basic principles of a tie-point reconstruction algorithm and consider ambiguities and resolution errors. Finally we mention extensions to plain stereoscopy such as a third view, a tomography-like approach and how magnetic field information can be used to improve the reconstruction.
[45]  oai:arXiv.org:astro-ph/0612650  [pdf] - 87988
Testing non-linear force-free coronal magnetic field extrapolations with the Titov-Demoulin equilibrium
Comments: 6 pages, 2 figures, Research Note
Submitted: 2006-12-21
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.