-
MHD modeling of a geoeffective interplanetary CME with the magnetic topology informed by in-situ observations
Authors:
E. Provornikova,
V. G. Merkin,
A. Vourlidas,
A. Malanushenko,
S. E. Gibson,
E. Winter,
N. Arge
Abstract:
Variations of the magnetic field within solar coronal mass ejections (CMEs) in the heliosphere depend on the CME`s magnetic structure as it leaves the solar corona and its subsequent evolution through interplanetary space. To account for this evolution, we developed a new numerical model of the inner heliosphere that simulates the propagation of a CME through a realistic background solar wind and…
▽ More
Variations of the magnetic field within solar coronal mass ejections (CMEs) in the heliosphere depend on the CME`s magnetic structure as it leaves the solar corona and its subsequent evolution through interplanetary space. To account for this evolution, we developed a new numerical model of the inner heliosphere that simulates the propagation of a CME through a realistic background solar wind and allows various CME magnetic topologies. To this end, we incorporate the Gibson-Low CME model within our global MHD model of the inner heliosphere, GAMERA-Helio. We apply the model to study the propagation of the geoeffective CME that erupted on 3 April, 2010 with the aim to reproduce the temporal variations of the magnetic field vector during the CME passage by Earth. Parameters of the Gibson-Low CME are informed by STEREO white-light observations near the Sun. The magnetic topology for this CME - the tethered flux rope - is informed by in-situ magnetic field observations near Earth. We performed two simulations testing different CME propagation directions. For an in-ecliptic direction, the simulation shows a rotation of all three magnetic field components within the CME, as seen at Earth, similar to that observed. With a southward propagation direction, suggested by coronal imaging observations, the modeled By and Bz components are consistent with the ACE data, but the Bx component lacks the observed change from negative to positive. In both cases, the model favors the East-West orientation of the CME flux rope, consistent with the orientation previously inferred from the STEREO/HI heliospheric images.
△ Less
Submitted 20 May, 2024;
originally announced May 2024.
-
Towards data-driven modeling and real-time prediction of solar flares and coronal mass ejections
Authors:
M. Rempel,
Y. Fan,
M. Dikpati,
A. Malanushenko,
M. D. Kazachenko,
M. C. M. Cheung,
G. Chintzoglou,
X. Sun,
G. H. Fisher,
T. Y. Chen
Abstract:
Modeling of transient events in the solar atmosphere requires the confluence of 3 critical elements: (1) model sophistication, (2) data availability, and (3) data assimilation. This white paper describes required advances that will enable statistical flare and CME forecasting (e.g. eruption probability and timing, estimation of strength, and CME details, such as speed and magnetic field orientatio…
▽ More
Modeling of transient events in the solar atmosphere requires the confluence of 3 critical elements: (1) model sophistication, (2) data availability, and (3) data assimilation. This white paper describes required advances that will enable statistical flare and CME forecasting (e.g. eruption probability and timing, estimation of strength, and CME details, such as speed and magnetic field orientation) similar to weather prediction on Earth.
△ Less
Submitted 29 December, 2022;
originally announced December 2022.
-
The Coronal Veil
Authors:
A. Malanushenko,
M. C. M. Cheung,
C. E. DeForest,
J. A. Klimchuk,
M. Rempel
Abstract:
Coronal loops, seen in solar coronal images, are believed to represent emission from magnetic flux tubes with compact cross-sections. We examine the 3D structure of plasma above an active region in a radiative magnetohydrodynamic simulation to locate volume counterparts for coronal loops. In many cases, a loop cannot be linked to an individual thin strand in the volume. While many thin loops are p…
▽ More
Coronal loops, seen in solar coronal images, are believed to represent emission from magnetic flux tubes with compact cross-sections. We examine the 3D structure of plasma above an active region in a radiative magnetohydrodynamic simulation to locate volume counterparts for coronal loops. In many cases, a loop cannot be linked to an individual thin strand in the volume. While many thin loops are present in the synthetic images, the bright structures in the volume are fewer, and of complex shape. We demonstrate that this complexity can form impressions of thin bright loops, even in the absence of thin bright plasma strands. We demonstrate the difficulty of discerning from observations whether a particular loop corresponds to a strand in the volume, or a projection artifact. We demonstrate how apparently isolated loops could deceive observers, even when observations from multiple viewing angles are available.
While we base our analysis on a simulation, the main findings are independent from a particular simulation setup and illustrate the intrinsic complexity involved in interpreting observations resulting from line-of-sight integration in an optically thin plasma.
We propose alternative interpretation for strands seen in EUV images of the corona. The "coronal veil" hypothesis is mathematically more generic, and naturally explains properties of loops that are difficult to address otherwise -- such as their constant cross section and anomalously high density scale height. We challenge the paradigm of coronal loops as thin magnetic flux tubes, offering new understanding of solar corona and, by extension, of other magnetically confined bright, hot plasmas.
△ Less
Submitted 29 November, 2021; v1 submitted 28 June, 2021;
originally announced June 2021.
-
Multi-spacecraft observations of coronal loops to verify a force-free field reconstruction and infer loop cross sections
Authors:
Marika I. McCarthy,
Dana W. Longcope,
Anna Malanushenko
Abstract:
Active region EUV loops are believed to trace a subset of magnetic field lines through the corona. Malanushenko et al. (2009) proposed a method, using loop images and line-of-sight photospheric magnetograms, to infer the three-dimensional shape and field strength along each loop. McCarthy et al. (2019) used this novel method to compute the total magnetic flux interconnecting a pair of active regio…
▽ More
Active region EUV loops are believed to trace a subset of magnetic field lines through the corona. Malanushenko et al. (2009) proposed a method, using loop images and line-of-sight photospheric magnetograms, to infer the three-dimensional shape and field strength along each loop. McCarthy et al. (2019) used this novel method to compute the total magnetic flux interconnecting a pair of active regions observed by SDO/AIA. They adopted the common assumption that each loop had a circular cross section. The accuracy of inferred shape and circularity of cross sections can both be tested using observations of the same loops from additional vantage points as provided by STEREO/EUVI. Here, we use multiple viewing angles to confirm the three-dimensional structure of loops. Of 151 viable cases, 105 (69.5%) matched some form of visible coronal structure when viewed approximately in quadrature. A loop with a circular cross-section should appear of a similar width in different perspectives. In contradiction to this, we find a puzzling lack of correlation between loop diameters seen from different perspectives, even an anti-correlation in some cases. Features identified as monolithic loops in AIA may, in fact, be more complex density enhancements. The 30.5% of reconstructions from AIA which did not match any feature in EUVI might be such enhancements. Others may be genuine loop structures, but with elliptical cross sections. We observe an anti-correlation between diameter and brightness, lending support to the latter hypothesis. Of 13 suitable for width analysis, four loops are consistent with non-circular cross sections, where we find anti-correlation in both comparisons.
△ Less
Submitted 6 April, 2021;
originally announced April 2021.
-
Untangling the global coronal magnetic field with multiwavelength observations
Authors:
S. E. Gibson,
A. Malanushenko,
G. de Toma,
S. Tomczyk,
K. Reeves,
H. Tian,
Z. Yang,
B. Chen,
G. Fleishman,
D. Gary,
G. Nita,
V. M. Pillet,
S. White,
U. Bąk-Stęślicka,
K. Dalmasse,
T. Kucera,
L. A. Rachmeler,
N. E. Raouafi,
J. Zhao
Abstract:
Magnetism defines the complex and dynamic solar corona. Coronal mass ejections (CMEs) are thought to be caused by stresses, twists, and tangles in coronal magnetic fields that build up energy and ultimately erupt, hurling plasma into interplanetary space. Even the ever-present solar wind possesses a three-dimensional morphology shaped by the global coronal magnetic field, forming geoeffective coro…
▽ More
Magnetism defines the complex and dynamic solar corona. Coronal mass ejections (CMEs) are thought to be caused by stresses, twists, and tangles in coronal magnetic fields that build up energy and ultimately erupt, hurling plasma into interplanetary space. Even the ever-present solar wind possesses a three-dimensional morphology shaped by the global coronal magnetic field, forming geoeffective corotating interaction regions. CME evolution and the structure of the solar wind depend intimately on the coronal magnetic field, so comprehensive observations of the global magnetothermal atmosphere are crucial both for scientific progress and space weather predictions. Although some advances have been made in measuring coronal magnetic fields locally, synoptic measurements of the global coronal magnetic field are not yet available.
We conclude that a key goal for 2050 should be comprehensive, ongoing 3D synoptic maps of the global coronal magnetic field. This will require the construction of new telescopes, ground and space-based, to obtain complementary, multiwavelength observations sensitive to the coronal magnetic field. It will also require development of inversion frameworks capable of incorporating multi-wavelength data, and forward analysis tools and simulation testbeds to prioritize and establish observational requirements on the proposed telescopes.
△ Less
Submitted 17 December, 2020;
originally announced December 2020.
-
Measuring and modeling the rate of separator reconnection between an emerging and existing active region
Authors:
Marika I. McCarthy,
Dana W. Longcope,
Anna Malanushenko,
David E. McKenzie
Abstract:
Magnetic reconnection occurs when new flux emerges into the corona and becomes incorporated into the existing coronal field. A new active region (AR) emerging in the vicinity of an existing AR provides a convenient laboratory in which reconnection of this kind can be quantified. We use high time-cadence 171 $Å$ data from SDO/AIA focused on new/old active region pair 11147/11149, to quantify reconn…
▽ More
Magnetic reconnection occurs when new flux emerges into the corona and becomes incorporated into the existing coronal field. A new active region (AR) emerging in the vicinity of an existing AR provides a convenient laboratory in which reconnection of this kind can be quantified. We use high time-cadence 171 $Å$ data from SDO/AIA focused on new/old active region pair 11147/11149, to quantify reconnection. We identify new loops as brightenings within a strip of pixels between the regions. This strategy is premised on the assumption that the energy brightening a loop originates in magnetic reconnection. We catalog 301 loops observed in the 48-hour time period beginning with the emergence of AR 11149. The rate at which these loops appear between the two ARs is used to calculate the reconnection rate between them. We then fit these loops with magnetic field, solving for each loop's field strength, geometry, and twist (via its proxy, coronal $α$). We find the rate of newly-brightened flux overestimates the flux which could be undergoing reconnection. This excess can be explained by our finding that the interconnecting region is not at its lowest energy (constant-$α$) state; the extrapolations exhibit loop-to-loop variation in $α$. This flux overestimate may result from the slow emergence of AR 11149, allowing time for Taylor relaxation internal to the domain of the reconnected flux to bring the $α$ distribution towards a single value which provides another mechanism for brightening loops after they are first created.
△ Less
Submitted 14 November, 2019;
originally announced November 2019.
-
The Influence of Spatial Resolution on Nonlinear Force-Free Modeling
Authors:
M. L. DeRosa,
M. S. Wheatland,
K. D. Leka,
G. Barnes,
T. Amari,
A. Canou,
S. A. Gilchrist,
J. K. Thalmann,
G. Valori,
T. Wiegelmann,
C. J. Schrijver,
A. Malanushenko,
X. Sun,
S. Régnier
Abstract:
The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLF…
▽ More
The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.
△ Less
Submitted 21 August, 2015;
originally announced August 2015.
-
Blind Stereoscopy of the Coronal Magnetic Field
Authors:
Markus J. Aschwanden,
Carolus J. Schrijver,
Anna Malanushenko
Abstract:
We test the feasibility of 3D coronal-loop tracing in stereoscopic EUV image pairs, with the ultimate goal of enabling efficient 3D reconstruction of the coronal magnetic field that drives flares and coronal mass ejections (CMEs). We developed an automated code designed to perform triangulation of coronal loops in pairs (or triplets) of EUV images recorded from different perspectives. The automate…
▽ More
We test the feasibility of 3D coronal-loop tracing in stereoscopic EUV image pairs, with the ultimate goal of enabling efficient 3D reconstruction of the coronal magnetic field that drives flares and coronal mass ejections (CMEs). We developed an automated code designed to perform triangulation of coronal loops in pairs (or triplets) of EUV images recorded from different perspectives. The automated (or blind) stereoscopy code includes three major tasks: (i) automated pattern recognition of coronal loops in EUV images, (ii) automated pairing of corresponding loop patterns from two different aspect angles, and (iii) stereoscopic triangulation of 3D loop coordinates. We perform tests with simulated stereoscopic EUV images and quantify the accuracy of all three procedures. In addition we test the performance of the blind stereoscopy code as a function of the spacecraft-separation angle and as a function of the spatial resolution. We also test the sensitivity to magnetic non-potentiality. The automated code developed here can be used for analysis of existing {\sl Solar TErrestrial RElationship Observatory (STEREO)} data, but primarily serves for a design study of a future mission with dedicated diagnostics of non-potential magnetic fields. For a pixel size of 0.6\arcsec (corresponding to the {\sl Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA)} spatial resolution of 1.4\arcsec), we find an optimum spacecraft-separation angle of $α_s \approx 5^\circ$.
△ Less
Submitted 15 June, 2015;
originally announced June 2015.
-
Thermal Diagnostics with the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory: A Validated Method for Differential Emission Measure Inversions
Authors:
Mark C. M. Cheung,
P. Boerner,
C. J. Schrijver,
P. Testa,
F. Chen,
H. Peter,
A. Malanushenko
Abstract:
We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include…
▽ More
We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a non-linear force-free field, and (3) thermodynamic models from a fully-compressible, 3D MHD simulation of AR corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and XRT data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.
△ Less
Submitted 13 April, 2015;
originally announced April 2015.
-
Deciphering Solar Magnetic Activity I: On The Relationship Between The Sunspot Cycle And The Evolution Of Small Magnetic Features
Authors:
Scott W. McIntosh,
Xin Wang,
Robert J. Leamon,
Rachel Howe,
Larisza D. Krista,
Anna V. Malanushenko,
Jonathan W. Cirtain,
Joseph B. Gurman,
Michael J. Thompson
Abstract:
Sunspots are a canonical marker of the Sun's internal magnetic field which flips polarity every ~22-years. The principal variation of sunspots, an ~11-year variation in number, modulates the amount of magnetic field that pierces the solar surface and drives significant variations in our Star's radiative, particulate and eruptive output over that period. This paper presents observations from the So…
▽ More
Sunspots are a canonical marker of the Sun's internal magnetic field which flips polarity every ~22-years. The principal variation of sunspots, an ~11-year variation in number, modulates the amount of magnetic field that pierces the solar surface and drives significant variations in our Star's radiative, particulate and eruptive output over that period. This paper presents observations from the Solar and Heliospheric Observatory and Solar Dynamics Observatory indicating that the 11-year sunspot variation is intrinsically tied it to the spatio-temporal overlap of the activity bands belonging to the 22-year magnetic activity cycle. Using a systematic analysis of ubiquitous coronal brightpoints, and the magnetic scale on which they appear to form, we show that the landmarks of sunspot cycle 23 can be explained by considering the evolution and interaction of the overlapping activity bands of the longer scale variability.
△ Less
Submitted 1 April, 2014; v1 submitted 12 March, 2014;
originally announced March 2014.
-
Using Coronal Loops to Reconstruct the Magnetic Field of an Active Region Before and After a Major Flare
Authors:
A. Malanushenko,
C. Schrijver,
M. L. DeRosa,
M. S. Wheatland
Abstract:
The shapes of solar coronal loops are sensitive to the presence of electrical currents that are the carriers of the nonpotential energy available for impulsive activity. We use this information in a new method for modeling the coronal magnetic field of AR 11158 as a nonlinear force-free field (NLFFF). The observations used are coronal images around time of major flare activity on 2011/02/15, toget…
▽ More
The shapes of solar coronal loops are sensitive to the presence of electrical currents that are the carriers of the nonpotential energy available for impulsive activity. We use this information in a new method for modeling the coronal magnetic field of AR 11158 as a nonlinear force-free field (NLFFF). The observations used are coronal images around time of major flare activity on 2011/02/15, together with the surface line-of-sight magnetic field measurements. The data are from the Helioseismic and Magnetic Imager and Atmospheric Imaging Assembly (HMI and AIA, respectively) onboard the Solar Dynamics Observatory (SDO). The model fields are constrained to approximate the coronal loop configurations as closely as possible, while also subject to the force-free constraints. The method does not use transverse photospheric magnetic field components as input, and is thereby distinct from methods for modeling NLFFFs based on photospheric vector magnetograms. We validate the method using observations of AR 11158 at a time well before major flaring, and subsequently review the field evolution just prior to and following an X2.2 flare and associated eruption. The models indicate that the energy released during the instability is about $1\times10^{32}$ erg, consistent with what is needed to power such a large eruptive flare. Immediately prior to the eruption the model field contains a compact sigmoid bundle of twisted flux that is not present in the post-eruption models, which is consistent with the observations. The core of that model structure is twisted by $\approx0.9$ full turns about its axis.
△ Less
Submitted 18 December, 2013;
originally announced December 2013.
-
On Anisotropy in Expansion of Magnetic Flux Tubes in the Solar Corona
Authors:
A. Malanushenko,
C. J. Schrijver
Abstract:
Most 1d hydrodynamic models of plasma confined to magnetic flux tubes assume circular cross-section of these tubes. We use potential field models to show that flux tubes in circumstances relevant to the solar corona do not in general maintain the same cross-sectional shape through their length and therefore the assumption of a circular cross-section is rarely true. We support our hypothesis with m…
▽ More
Most 1d hydrodynamic models of plasma confined to magnetic flux tubes assume circular cross-section of these tubes. We use potential field models to show that flux tubes in circumstances relevant to the solar corona do not in general maintain the same cross-sectional shape through their length and therefore the assumption of a circular cross-section is rarely true. We support our hypothesis with mathematical reasoning and numeric experiments. We demonstrate that lifting this assumption in realistic non-circular loops make apparent expansion of magnetic flux tubes consistent with that of observed coronal loops. We propose that in a bundle of ribbon-like loops those that are viewed along the wide direction would stand out against those that are viewed across the wide direction, due to the difference in their column depths. That would impose a bias towards selecting loops that appear not to be expanding seen projected in the plane of sky. An implication of this selection bias is that the preferentially selected non-circular loops would appear to have increased pressure scale height even if they are resolved by current instruments.
△ Less
Submitted 28 June, 2013;
originally announced July 2013.
-
First 3D Reconstructions of Coronal Loops with the STEREO A+B Spacecraft: IV. Magnetic Modeling with Twisted Force-Free Fields
Authors:
Markus J. Aschwanden,
Jean-Pierre Wuelser,
Nariaki V. Nitta,
James R. Lemen,
Marc L. DeRosa,
Anna Malanushenko
Abstract:
The three-dimensional (3D) coordinates of stereoscopically triangulated loops provide strong constraints for magnetic field models of active regions in the solar corona. Here we use STEREO/A and B data from some 500 stereoscopically triangulated loops observed in four active regions (2007 Apr 30, May 9, May 19, Dec 11), together with SOHO/MDI line-of-sight magnetograms. We measure the average misa…
▽ More
The three-dimensional (3D) coordinates of stereoscopically triangulated loops provide strong constraints for magnetic field models of active regions in the solar corona. Here we use STEREO/A and B data from some 500 stereoscopically triangulated loops observed in four active regions (2007 Apr 30, May 9, May 19, Dec 11), together with SOHO/MDI line-of-sight magnetograms. We measure the average misalignment angle between the stereoscopic loops and theoretical magnetic field models, finding a mismatch of $μ=19^\circ-46^\circ$ for a potential field model, which is reduced to $μ=14^\circ-19^\circ$ for a non-potential field model parameterized by twist parameters. The residual error is commensurable with stereoscopic measurement errors ($μ_{SE} \approx 8^\circ-12^\circ$). We developed a potential field code that deconvolves a line-of-sight magnetogram into three magnetic field components $(B_x, B_y, B_z)$, as well as a non-potential field forward-fitting code that determines the full length of twisted loops ($L \approx 50-300$ Mm), the number of twist turns (median $N_{twist}=0.06$), the nonlinear force-free $α$-parameter (median $α\approx 4 \times 10^{-11}$ cm$^{-1}$), and the current density (median $j_z \approx 1500$ Mx cm$^{-2}$ s$^{-1}$). All twisted loops are found to be far below the critical value for kink instability, and Joule dissipation of their currents is found be be far below the coronal heating requirement. The algorithm developed here, based on an analytical solution of nonlinear force-free fields that is accurate to second order (in the force-free parameter $α$), represents the first code that enables fast forward-fitting to photospheric magnetograms and stereoscopically triangulated loops in the solar corona.
△ Less
Submitted 11 July, 2012;
originally announced July 2012.
-
A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. II. Numeric Code and Tests
Authors:
Markus J. Aschwanden,
Anna Malanushenko
Abstract:
Based on a second-order approximation of nonlinear force-free magnetic field solutions in terms of uniformly twisted field lines derived in Paper I, we develop here a numeric code that is capable to forward-fit such analytical solutions to arbitrary magnetogram (or vector magnetograph) data combined with (stereoscopically triangulated) coronal loop 3D coordinates. We test the code here by forward-…
▽ More
Based on a second-order approximation of nonlinear force-free magnetic field solutions in terms of uniformly twisted field lines derived in Paper I, we develop here a numeric code that is capable to forward-fit such analytical solutions to arbitrary magnetogram (or vector magnetograph) data combined with (stereoscopically triangulated) coronal loop 3D coordinates. We test the code here by forward-fitting to six potential field and six nonpotential field cases simulated with our analytical model, as well as by forward-fitting to an exactly force-free solution of the Low and Lou (1990) model. The forward-fitting tests demonstrate: (i) a satisfactory convergence behavior (with typical misalignment angles of $μ\approx 1^\circ-10^\circ$), (ii) relatively fast computation times (from seconds to a few minutes), and (iii) the high fidelity of retrieved force-free $α$-parameters ($α_{\rm fit}/α_{\rm model} \approx 0.9-1.0$ for simulations and $α_{\rm fit}/α_{\rm model} \approx 0.7\pm0.3$ for the Low and Lou model). The salient feature of this numeric code is the relatively fast computation of a quasi-forcefree magnetic field, which closely matches the geometry of coronal loops in active regions, and complements the existing {\sl nonlinear force-free field (NLFFF)} codes based on photospheric magnetograms without coronal constraints.
△ Less
Submitted 11 July, 2012;
originally announced July 2012.
-
Direct Measurements of Magnetic Twist in the Solar Corona
Authors:
A. Malanushenko,
M. H. Yusuf,
D. W. Longcope
Abstract:
In the present work we study evolution of magnetic helicity in the solar corona. We compare the rate of change of a quantity related to the magnetic helicity in the corona to the flux of magnetic helicity through the photosphere and find that the two rates are similar. This gives observational evidence that helicity flux across the photosphere is indeed what drives helicity changes in solar corona…
▽ More
In the present work we study evolution of magnetic helicity in the solar corona. We compare the rate of change of a quantity related to the magnetic helicity in the corona to the flux of magnetic helicity through the photosphere and find that the two rates are similar. This gives observational evidence that helicity flux across the photosphere is indeed what drives helicity changes in solar corona during emergence.
For the purposes of estimating coronal helicity we neither assume a strictly linear force-free field, nor attempt to construct a non-linear force-free field. For each coronal loop evident in Extreme Ultraviolet (EUV) we find a best-matching line of a linear force-free field and allow the twist parameter alpha to be different for each line. This method was introduced and its applicability was discussed in Malanushenko et. al. (2009).
The object of the study is emerging and rapidly rotating AR 9004 over about 80 hours. As a proxy for coronal helicity we use the quantity <alpha_i*L_i/2> averaged over many reconstructed lines of magnetic field. We argue that it is approximately proportional to "flux-normalized" helicity H/Phi^2, where H is helicity and Phi is total enclosed magnetic flux of the active region. The time rate of change of such quantity in the corona is found to be about 0.021 rad/hr, which is compatible with the estimates for the same region obtained using other methods Longcope et. al. (2007), who estimated the flux of normalized helicity of about 0.016 rad/hr.
△ Less
Submitted 24 February, 2012;
originally announced February 2012.
-
Guiding Nonlinear Force-Free Modeling Using Coronal Observations: First Results Using a Quasi Grad-Rubin Scheme
Authors:
A. Malanushenko,
C. J. Schrijver,
M. L. DeRosa,
M. S. Wheatland,
S. A. Gilchrist
Abstract:
Presently, many models of the coronal magnetic field rely on photospheric vector magnetograms but these data have been shown to be problematic as the sole boundary information for nonlinear force-free field (NLFFF) extrapolations. Magnetic fields in the corona manifest themselves in high-energy images (X-rays and EUV) in the shapes of coronal loops, providing an additional constraint that at prese…
▽ More
Presently, many models of the coronal magnetic field rely on photospheric vector magnetograms but these data have been shown to be problematic as the sole boundary information for nonlinear force-free field (NLFFF) extrapolations. Magnetic fields in the corona manifest themselves in high-energy images (X-rays and EUV) in the shapes of coronal loops, providing an additional constraint that at present is not used due to the mathematical complications of incorporating such input into numerical models. Projection effects and the limited number of usable loops further complicate the use of coronal information. We develop and test an algorithm to use images showing coronal loops in the modeling of the solar coronal magnetic field. We first fit projected field lines with field lines of constant-\als force-free fields to approximate the three-dimensional distribution of currents in the corona along a sparse set of trajectories. We then apply a Grad-Rubin-like iterative technique to obtain a volume-filling nonlinear force-free model of the magnetic field, modifying method presented in \citet{Wheatland2007}. We thoroughly test the technique on known analytical and solar-like model magnetic fields previously used for comparing different extrapolation techniques \citep{Schrijver2006, Schrijver2008} and compare the results with those obtained by presently available methods that rely only on the photospheric data. We conclude that we have developed a functioning method of modeling the coronal magnetic field by combining the line-of-sight component of photospheric magnetic field with information from coronal images. Vector magnetograms over the full or partial photospheric boundary of the numerical domain could optionally be used.
△ Less
Submitted 24 February, 2012;
originally announced February 2012.
-
Reconstructing the Local Twist of Coronal Magnetic Fields and the Three-Dimensional Shape of the Field Lines from Coronal Loops in EUV and X-Ray Images
Authors:
A. Malanushenko,
D. W. Longcope,
D. E. McKenzie
Abstract:
Non-linear force-free fields are the most general case of force-free fields, but the hardest to model as well. There are numerous methods of computing such fields by extrapolating vector magnetograms from the photosphere, but very few attempts have so far made quantitative use of coronal morphology. We present a method to make such quantitative use of X-Ray and EUV images of coronal loops. Each…
▽ More
Non-linear force-free fields are the most general case of force-free fields, but the hardest to model as well. There are numerous methods of computing such fields by extrapolating vector magnetograms from the photosphere, but very few attempts have so far made quantitative use of coronal morphology. We present a method to make such quantitative use of X-Ray and EUV images of coronal loops. Each individual loop is fit to a field line of a linear force-free field, allowing the estimation of the field line's twist, three-dimensional geometry and the field strength along it.
We assess the validity of such a reconstruction since the actual corona is probably not a linear force-free field and that the superposition of linear force-free fields is generally not itself a force-free field. To do so, we perform a series of tests on non-linear force-free fields, described in Low & Lou (1990). For model loops we project field lines onto the photosphere. We compare several results of the method with the original field, in particular the three-dimensional loop shapes, local twist (coronal alpha), distribution of twist in the model photosphere and strength of the magnetic field. We find that, (i) for these trial fields, the method reconstructs twist with mean absolute deviation of at most 15% of the range of photospheric twist, (ii) that heights of the loops are reconstructed with mean absolute deviation of at most 5% of the range of trial heights and (iii) that the magnitude of non-potential contribution to photospheric field is reconstructed with mean absolute deviation of at most 10% of the maximal value.
△ Less
Submitted 28 September, 2009;
originally announced September 2009.
-
Additive Self Helicity as a Kink Mode Threshold
Authors:
A. Malanushenko,
D. W. Longcope,
Y. Fan,
S. E. Gibson
Abstract:
In this paper we propose that additive self helicity, introduced by Longcope and Malanushenko (2008), plays a role in the kink instability for complex equilibria, similar to twist helicity for thin flux tubes (Hood and Priest (1979), Berger and Field (1984)). We support this hypothesis by a calculation of additive self helicity of a twisted flux tube from the simulation of Fan and Gibson (2003).…
▽ More
In this paper we propose that additive self helicity, introduced by Longcope and Malanushenko (2008), plays a role in the kink instability for complex equilibria, similar to twist helicity for thin flux tubes (Hood and Priest (1979), Berger and Field (1984)). We support this hypothesis by a calculation of additive self helicity of a twisted flux tube from the simulation of Fan and Gibson (2003). As more twist gets introduced, the additive self helicity increases, and the kink instability of the tube coincides with the drop of additive self helicity, after the latter reaches the value of $H_A/Φ^2\approx 1.5$ (where $Φ$ is the flux of the tube and $H_A$ is additive self helicity).
We compare additive self helicity to twist for a thin sub-portion of the tube to illustrate that $H_A/Φ^2$ is equal to the twist number, studied by Berger and Field (1984), when the thin flux tube approximation is applicable. We suggest, that the quantity $H_A/Φ^2$ could be treated as a generalization of a twist number, when thin flux tube approximation is not applicable. A threshold on a generalized twist number might prove extremely useful studying complex equilibria, just as twist number itself has proven useful studying idealized thin flux tubes. We explicitly describe a numerical method for calculating additive self helicity, which includes an algorithm for identifying a domain occupied by a flux bundle and a method of calculating potential magnetic field confined to this domain. We also describe a numerical method to calculate twist of a thin flux tube, using a frame parallelly transported along the axis of the tube.
△ Less
Submitted 27 September, 2009;
originally announced September 2009.