Normalized to: Inhester, B.
[1]
oai:arXiv.org:2004.03790 [pdf] - 2097397
Initiation and Early Kinematic Evolution of Solar Eruptions
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
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
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
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
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
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
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
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?
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
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
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
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
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
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
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
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
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
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
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
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?
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
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
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)
Peter, H.;
Abbo, L.;
Andretta, V.;
Auchere, F.;
Bemporad, A.;
Berrilli, F.;
Bommier, V.;
Braukhane, A.;
Casini, R.;
Curdt, W.;
Davila, J.;
Dittus, H.;
Fineschi, S.;
Fludra, A.;
Gandorfer, A.;
Griffin, D.;
Inhester, B.;
Lagg, A.;
Degl'Innocenti, E. Landi;
Maiwald, V.;
Sainz, R. Manso;
Pillet, V. Martinez;
Matthews, S.;
Moses, D.;
Parenti, S.;
Pietarila, A.;
Quantius, D.;
Raouafi, N. -E.;
Raymond, J.;
Rochus, P.;
Romberg, O.;
Schlotterer, M.;
Schuehle, U.;
Solanki, S.;
Spadaro, D.;
Teriaca, L.;
Tomczyk, S.;
Bueno, J. Trujillo;
Vial, J. -C.
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
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
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
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
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
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
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?
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
Feng, L.;
Inhester, B.;
Solanki, S. K.;
Wilhelm, K.;
Wiegelmann, T.;
Podlipnik, B.;
Howard, R. A.;
Plunkett, S. P.;
Wuelser, J. P.;
Gan, W. Q.
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
DeRosa, Marc L.;
Schrijver, Carolus J.;
Barnes, Graham;
Leka, K. D.;
Lites, Bruce W.;
Aschwanden, Markus J.;
Amari, Tahar;
Canou, Aurelien;
McTiernan, James M.;
Regnier, Stephane;
Thalmann, Julia K.;
Valori, Gherardo;
Wheatland, Michael S.;
Wiegelmann, Thomas;
Cheung, Mark C. M.;
Conlon, Paul A.;
Fuhrmann, Marcel;
Inhester, Bernd;
Tadesse, Tilaye
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
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
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?
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
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
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
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
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
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
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
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
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
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.