Source author record

Markus J. Aschwanden

Markus J. Aschwanden appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

50works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

50 published item(s)

preprint2026arXiv

The Energy-Duration Relationship in Astrophysical Self-Organized Criticality Systems

Scaling laws in astrophysical systems that involve the energy, the geometry, and the spatio-temporal evolution, provide the theoretical framework for physical models of energy dissipation processes. A leading model is the standard fractal-diffusive self-organized criticality (FD-SOC) model, which is built on four fundamental assumptions: (i) the dimensionality $d=3$, (ii) the fractal dimension $D_V=d-1/2=2.5$, (iii) classical diffusion $L \propto T^{(1/2)}$, and (iv) the proportionality of the dissipated energy to the fractal volume $E \propto V$. Based on these assumptions, the FD-SOC model predicts a scaling law of $T \propto E^k \propto E^{(4/5)} = E^{0.8}$. On the observational side, we find empirical scaling laws of $T \propto E^{0.81\pm0.03}$ by Peng et al.~(2023) and $T \propto E^{0.86\pm0.03}$ by Araujo \& Valio (2021) that are self-consistent with the theoretical prediction of the FD-SOC model. However, cases with a small time range $q_T = \log{(T_{max}/T_{min})} \lapprox 2$ have large statistical uncertainties and systematic errors, which produces smaller scaling law exponents ($k \approx 0.3, ..., 0.6$) as a consequence. The close correlation of the scaling exponent $k$ with the truncation bias $q_T$ implies that the dispersion of k-values is an observational effect, rather than a physical property.

preprint2022arXiv

Interface Region Imaging Spectrograph (IRIS) Observations of the Fractal Dimension in the Solar Atmosphere

While previous work explored the fractality and self-organized criticality (SOC) of flares and nanoflares in wavelengths emitted in the solar corona (such as in hard X-rays, soft X-rays, and EUV wavelenghts), we focus here on impulsive phenomena in the photosphere and transition region, as observed with the {\sl Interface Region Imaging Spectrograph (IRIS)} in the temperature range of $T_e \approx 10^4-10^6$ K. We find the following fractal dimensions (in increasing order): $D_A=1.21 \pm 0.07$ for photospheric granulation, $D_A=1.29 \pm 0.15$ for plages in the transition region, $D_A=1.54 \pm 0.16$ for sunspots in the transition region, $D_A=1.59 \pm 0.08$ for magnetograms in active regions, $D_A=1.56 \pm 0.08$ for EUV nanoflares, $D_A=1.76 \pm 0.14$ for large solar flares, and up to $D_A=1.89 \pm 0.05$ for the largest X-class flares. We interpret low values of the fractal dimension ($1.0 \lapprox D_A \lapprox 1.5$) in terms of sparse curvi-linear flow patterns, while high values of the fractal dimension ($1.5 \lapprox D_A \lapprox 2.0$) indicate near space-filling transport processes, such as chromospheric evaporation. Phenomena in the solar transition region appear to be consistent with SOC models, based on their size distributions of fractal areas $A$ and (radiative) energies $E$, which show power law slopes of $α_A^{obs}=2.51 \pm 0.21$ (with $α_A^{theo}=2.33$ predicted), and $α_E^{obs}=2.03 \pm 0.18$ (with $α_E^{theo}=1.80$ predicted).

preprint2022arXiv

Reconciling Power Law Slopes in Solar Flare and Nanoflare Size Distributions

We unify the power laws of size distributions of solar flare and nanoflare energies. We present three models that predict the power law slopes $α_E$ of flare energies defined in terms of the 2-D and 3-D fractal dimensions ($D_A, D_V$): (i) The spatio-temporal standard SOC model, defined by the power law slope $α_{E1}=1+2/(D_V+2)=(13/9)\approx 1.44$; (ii) the 2-D thermal energy model, $α_{E2}=1+2/D_A=(7/3)\approx 2.33$, and (iii) the 3-D thermal energy model, $α_{E3}=1+2/D_V=(9/5)\approx 1.80$. The theoretical predictions of energies are consistent with the observational values of these three groups, i.e., $α_{E1}=1.47 \pm 0.07$; $α_{E2}=2.38 \pm 0.09$, and $α_{E3}=1.80 \pm 0.18$. These results corroborate that the energy of nanoflares does not diverge at small energies, since $(α_{E1}<2$) and $(α_{E3}<2)$, except for the unphyiscal 2-D model $(α_{E2}>2)$. This conclusion adds an additional argument against the scenario of coronal heating by nanoflares.

preprint2021arXiv

Self-Organized Criticality in Stellar Flares

Power law size distributions are the hallmarks of nonlinear energy dissipation processes governed by self-organized criticality. Here we analyze 75 data sets of stellar flare size distributions, mostly obtained from the {\sl Extreme Ultra-Violet Explorer (EUVE)} and the {\sl Kepler} mission. We aim to answer the following questions for size distributions of stellar flares: (i) What are the values and uncertainties of power law slopes? (ii) Do power law slopes vary with time ? (iii) Do power law slopes depend on the stellar spectral type? (iv) Are they compatible with solar flares? (v) Are they consistent with self-organized criticality (SOC) models? We find that the observed size distributions of stellar flare fluences (or energies) exhibit power law slopes of $α_E=2.09\pm0.24$ for optical data sets observed with Kepler. The observed power law slopes do not show much time variability and do not depend on the stellar spectral type (M, K, G, F, A, Giants). In solar flares we find that background subtraction lowers the uncorrected value of $α_E=2.20\pm0.22$ to $α_E=1.57\pm0.19$. Furthermore, most of the stellar flares are temporally not resolved in low-cadence (30 min) Kepler data, which causes an additional bias. Taking these two biases into account, the stellar flare data sets are consistent with the theoretical prediction $N(x) \propto x^{-α_x}$ of self-organized criticality models, i.e., $α_E=1.5$. Thus, accurate power law fits require automated detection of the inertial range and background subtraction, which can be modeled with the generalized Pareto distribution, finite-system size effects, and extreme event outliers.

preprint2020arXiv

Global Energetics of Solar Flares: X. Petschek Reconnection Rate and Alfven Mach Number of Magnetic Reconnection Outflows

We investigate physical scaling laws for magnetic energy dissipation in solar flares, in the framework of the Sweet-Parker model and the Petschek model. We find that the total dissipated magnetic energy $E_{diss}$ in a flare depends on the mean magnetic field component $B_f$ associated with the free energy $E_f$, the length scale $L$ of the magnetic area, the hydrostatic density scale height $λ$ of the solar corona, the Alfvén Mach number $M_A=v_1/v_A$ (the ratio of the inflow speed $v_1$ to the Alfvénic outflow speed $v_A$), and the flare duration $τ_f$, i.e., $E_{diss} = (1/4π) B_f^2\ L\ λ v_A\ M_A\ τ_f$, where the Alfvén speed depends on the nonpotential field strength $B_{np}$ and the mean electron density $n_e$ in the reconnection outflow. Using MDI/SDO and AIA/SDO observations and 3-D magnetic field solutions obtained with the vertical-current approximation nonlinear force-free field code (VCA-NLFFF) we measure all physical parameters necessary to test scaling laws, which represents a new method to measure Alfvén Mach numbers $M_A$, the reconnection rate, and the total free energy dissipated in solar flares.

preprint2020arXiv

Global Energetics of Solar Flares. XI. Flare Magnitude Predictions of the GOES-Class

In this study we determine scaling relationships of observed solar flares that can be used to predict upper limits of the GOES-class magnitude of solar flares. The flare prediction scheme is based on the scaling of the slowly-varying potential energy $E_p(t)$, which is extrapolated in time over an interval of $Δt \le$ 24 hrs. The observed scaling of the dissipated energy $E_{diss}$ scales with the potential field energy as $E_{diss} \propto E_p^{1.32}$. In addition, the observed scaling relationship of the flare volume, $V \propto E_{diss}^{1.17}$, the multi-thermal energy, $E_{th} \propto V^{0.76}$, the flare emission measure $EM \propto E_{th}^{0.79}$, the EM-weighted temperature $T_{w}$, and the GOES flux, $F_8(t) \propto E_p(t)^{0.92}$, allows us then to predict an upper limit of the GOES-class flare magnitude in the extrapolated time window. We find a good correlation (CCC$\approx 0.7$) between the observed and predicted GOES-class flare magnitudes (in 172 X and M-class events). This is the first algorithm that employs observed scaling laws of physical flare parameters to predict GOES flux upper limits, an important capability that complements previous flare prediction methods based on machine-learning algorithms used in space weather forecasting.

preprint2019arXiv

Global Energetics of Solar Flares. IX. Refined Magnetic Modeling

A more accurate analytical solution of the {\sl vertical-current approximation nonlinear force-free field (VCA3-NLFFF)} model is presented that includes besides the radial $(B_r)$ and the azimuthal $(B_φ)$ magnetic field components a poloidal component $(B_θ \neq 0)$ also. This new analytical solution is of second-order accuracy in the divergence-freeness condition, and of third-order accuracy in the force-freeness condition. We re-analyze the sample of 173 GOES M- and X-class flares observed with the {\sl Atmospheric Imaging Assembly (AIA)} and {\sl Helioseismic and Magnetic Imager (HMI)} onboard the {\sl Solar Dynamics Observatory (SDO)}. The new code reproduces helically twisted loops with a low winding number below the kink instability consistently, avoiding unstable, highly-twisted structures of the Gold-Hoyle flux rope type. The magnetic energies agree within $E_{VCA3}/E_W=0.99\pm0.21$ with the Wiegelmann (W-NLFFF) code. The time evolution of the magnetic field reveals multiple, intermittent energy build-up and releases in most flares, contradicting both the Rosner-Vaiana model (with gradual energy storage in the corona) and the principle of time scale separation ($τ_{flare} \ll τ_{storage}$) postulated in self-organized criticality models. The mean dissipated flare energy is found to amount to $7\%\pm3\%$ of the potential energy, or $60\%\pm26\%$ of the free energy, a result that can be used for predicting flare magnitudes based on the potential field of active regions.

preprint2019arXiv

Non-Stationary Fast-Driven Self-Organized Criticality in Solar Flares

The original concept of self-organized criticality (Bak et al.~1987), applied to solar flare statistics (Lu and Hamilton 1991), assumed a slow-driven and stationary flaring rate, which warrants time scale separation (between flare durations and inter-flare waiting times), it reproduces power-law distributions for flare peak fluxes and durations, but predicts an exponential waiting time distribution. In contrast to these classical assumptions we observe: (i) multiple energy dissipation episodes during most flares, (ii) violation of the principle of time scale separation, (iii) a fast-driven and non-stationary flaring rate, (iv) a power law distribution for waiting times $Δt$, with a slope of $α_{Δt} \approx 2.0$, as predicted from the universal reciprocality between mean flaring rates and mean waiting times; and (v) pulses with rise times and decay times of the dissipated magnetic free energy on time scales of $12\pm6$ min, up to 13 times in long-duration ($\lapprox 4$ hrs) flares. These results are inconsistent with coronal long-term energy storage (Rosner and Vaiana 1978), but require photospheric-chromospheric current injections into the corona.

preprint2016arXiv

Global Energetics of Solar Flares: IV. Coronal Mass Ejection Energetics

This study entails the fourth part of a global flare energetics project, in which the mass $m_{\mathrm{cme}}$, kinetic energy $E_{\mathrm{kin}}$, and the gravitational potential energy $E_{\mathrm{grav}}$ of coronal mass ejections (CMEs) is measured in 399 M and X-class flare events observed during the first 3.5 yrs of the Solar Dynamics Observatory (SDO) mission, using a new method based on the EUV dimming effect. The EUV dimming is modeled in terms of a radial adiabatic expansion process, which is fitted to the observed evolution of the total emission measure of the CME source region. The model derives the evolution of the mean electron density, the emission measure, the bulk plasma expansion velocity, the mass, and the energy in the CME source region. The EUV dimming method is truly complementary to the Thomson scattering method in white light, which probes the CME evolution in the heliosphere at $r > 2 R_{\odot}$, while the EUV dimming method tracks the CME launch in the corona. We compare the CME parameters obtained in white light with the LASCO/C2 coronagraph with those obtained from EUV dimming with the Atmospheric Imaging Assembly (AIA) onboard SDO for all identical events in both data sets. We investigate correlations between CME parameters, the relative timing with flare parameters, frequency occurrence distributions, and the energy partition between magnetic, thermal, nonthermal, and CME energies. CME energies are found to be systematically lower than the dissipated magnetic energies, which is consistent with a magnetic origin of CMEs.

preprint2016arXiv

The Vertical Current Approximation Nonlinear Force-Free Field Code - Description, Performance Tests, and Measurements of Magnetic Energies Dissipated in Solar Flares

In this work we provide an updated description of the Vertical Current Approximation Nonlinear Force-Free Field (VCA-NLFFF) code, which is designed to measure the evolution of the potential, nonpotential, free energies, and the dissipated magnetic energies during solar flares. This code provides a complementary and alternative method to existing traditional NLFFF codes. The chief advantages of the VCA-NLFFF code over traditional NLFFF codes are the circumvention of the unrealistic assumption of a force-free photosphere in the magnetic field extrapolation method, the capability to minimize the misalignment angles between observed coronal loops (or chromospheric fibril structures) and theoretical model field lines, as well as computational speed. In performance tests of the VCA-NLFFF code, by comparing with the NLFFF code of Wiegelmann (2004), we find agreement in the potential, nonpotential, and free energy within a factor of about 1.3, but the Wiegelmann code yields in the average a factor of 2 lower flare energies. The VCA-NLFFF code is found to detect decreases in flare energies in most X, M, and C-class flares. The successful detection of energy decreases during a variety of flares with the VCA-NLFFF code indicates that current-driven twisting and untwisting of the magnetic field is an adequate model to quantify the storage of magnetic energies in active regions and their dissipation during flares. - The VCA-NLFFF code is also publicly available in the Solar SoftWare (SSW).

preprint2015arXiv

25 Years of Self-Organized Criticality: Numerical Detection Methods

The detection and characterization of self-organized criticality (SOC), in both real and simulated data, has undergone many significant revisions over the past 25 years. The explosive advances in the many numerical methods available for detecting, discriminating, and ultimately testing, SOC have played a critical role in developing our understanding of how systems experience and exhibit SOC. In this article, methods of detecting SOC are reviewed; from correlations to complexity to critical quantities. A description of the basic autocorrelation method leads into a detailed analysis of application-oriented methods developed in the last 25 years. In the second half of this manuscript space-based, time-based and spatial-temporal methods are reviewed and the prevalence of power laws in nature is described, with an emphasis on event detection and characterization. The search for numerical methods to clearly and unambiguously detect SOC in data often leads us outside the comfort zone of our own disciplines - the answers to these questions are often obtained by studying the advances made in other fields of study. In addition, numerical detection methods often provide the optimum link between simulations and experiments in scientific research. We seek to explore this boundary where the rubber meets the road, to review this expanding field of research of numerical detection of SOC systems over the past 25 years, and to iterate forwards so as to provide some foresight and guidance into developing breakthroughs in this subject over the next quarter of a century.

preprint2015arXiv

Blind Stereoscopy of the Coronal Magnetic Field

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$.

preprint2015arXiv

Magnetic Energy Dissipation during the 2014 March 29 Solar Flares

We calculated the time evolution of the free magnetic energy during the 2014-Mar-29 flare (SOL2014-03-29T17:48), the first X-class flare detected by IRIS. The free energy was calculated from the difference between the nonpotential field, constrained by the geometry of observed loop structures, and the potential field. We use AIA/SDO and IRIS images to delineate the geometry of coronal loops in EUV wavelengths, as well as to trace magnetic field directions in UV wavelengths in the chromosphere and transition region. We find an identical evolution of the free energy for both the coronal and chromospheric tracers, as well as agreement between AIA and IRIS results, with a peak free energy of $E_{free}(t_{peak}) \approx (45 \pm 2) \times 10^{30}$ erg, which decreases by an amount of $ΔE_{free} \approx (29 \pm 3) \times 10^{30}$ erg during the flare decay phase. The consistency of free energies measured from different EUV and UV wavelengths for the first time here, demonstrates that vertical electric currents (manifested in form of helically twisted loops) can be detected and measured from both chromospheric and coronal tracers.

preprint2015arXiv

Thresholded Power Law Size Distributions of Instabilities in Astrophysics

Power law-like size distributions are ubiquitous in astrophysical instabilities. There are at least four natural effects that cause deviations from ideal power law size distributions, which we model here in a generalized way: (1) a physical threshold of an instability; (2) incomplete sampling of the smallest events below a threshold $x_0$; (3) contamination by an event-unrelated background $x_b$; and (4) truncation effects at the largest events due to a finite system size. These effects can be modeled in simplest terms with a "thresholded power law" distribution function (also called generalized Pareto [type II] or Lomax distribution), $N(x) dx \propto (x+x_0)^{-a} dx$, where $x_0 > 0$ is positive for a threshold effect, while $x_0 < 0$ is negative for background contamination. We analytically derive the functional shape of this thresholded power law distribution function from an exponential-growth evolution model, which produces avalanches only when a disturbance exceeds a critical threshold $x_0$. We apply the thresholded power law distribution function to terrestrial, solar (HXRBS, BATSE, RHESSI), and stellar flare (Kepler) data sets. We find that the thresholded power law model provides an adequate fit to most of the observed data. Major advantages of this model are the automated choice of the power law fitting range, diagnostics of background contamination, physical inastability thresholds, instrumental detection thresholds, and finite system size limits. When testing self-organized criticality models, which predict ideal power laws, we suggest to include these natural truncation effects.

preprint2014arXiv

25 Years of Self-Organized Criticality: Solar and Astrophysics

Shortly after the seminal paper {\sl "Self-Organized Criticality: An explanation of 1/f noise"} by Bak, Tang, and Wiesenfeld (1987), the idea has been applied to solar physics, in {\sl "Avalanches and the Distribution of Solar Flares"} by Lu and Hamilton (1991). In the following years, an inspiring cross-fertilization from complexity theory to solar and astrophysics took place, where the SOC concept was initially applied to solar flares, stellar flares, and magnetospheric substorms, and later extended to the radiation belt, the heliosphere, lunar craters, the asteroid belt, the Saturn ring, pulsar glitches, soft X-ray repeaters, blazars, black-hole objects, cosmic rays, and boson clouds. The application of SOC concepts has been performed by numerical cellular automaton simulations, by analytical calculations of statistical (powerlaw-like) distributions based on physical scaling laws, and by observational tests of theoretically predicted size distributions and waiting time distributions. Attempts have been undertaken to import physical models into the numerical SOC toy models, such as the discretization of magneto-hydrodynamics (MHD) processes. The novel applications stimulated also vigorous debates about the discrimination between SOC models, SOC-like, and non-SOC processes, such as phase transitions, turbulence, random-walk diffusion, percolation, branching processes, network theory, chaos theory, fractality, multi-scale, and other complexity phenomena. We review SOC studies from the last 25 years and highlight new trends, open questions, and future challenges, as discussed during two recent ISSI workshops on this theme.

preprint2014arXiv

Global Energetics of Solar Flares: I. Magnetic Energies

We present the first part of a project on the global energetics of solar flares and coronal mass ejections (CMEs) that includes about 400 M- and X-class flares observed with AIA and HMI onboard SDO. We calculate the potential energy, free energy, and the flare-dissipated magnetic energy. We calculate these magnetic parameters using two different NLFFF codes: The COR-NLFFF code uses the line-of-sight magnetic field component $B_z$ from HMI to define the potential field, and the 2D coordinates of automatically detected coronal loops in 6 coronal wavelengths from AIA to measure the helical twist of coronal loops caused by vertical currents, while the PHOT-NLFFF code extrapolates the photospheric 3D vector fields. We find agreement between the two codes in the measurement of free energies and dissipated energies within a factor of $ \approx 3$. The size distributions of magnetic parameters exhibit powerlaw slopes that are approximately consistent with the fractal-diffusive self-organized criticality model. The magnetic parameters exhibit scaling laws for the nonpotential energy, $E_{np} \propto E_p^{1.02}$, for the free energy, $E_{free} \propto E_p^{1.7}$ and $E_{free} \propto B_φ^{1.0} L^{1.5}$, for the dissipated energy, $E_{diss} \propto E_p^{1.6}$ and $E_{diss} \propto E_{free}^{0.9}$, and the energy dissipation volume, $V \propto E_{diss}^{1.2}$. The potential energies vary in the range of $E_p = 1 \times 10^{31} - 4 \times 10^{33}$ erg, while the free energy has a ratio of $E_{free}/E_p \approx 1%-25%$. The Poynting flux amounts to $F_{flare} \approx 5 \times 10^{8} - 10^{10}$ erg cm$^{-2}$ s$^{-1}$ during flares, which averages to $F_{AR} \approx 6 \times 10^6$ erg cm$^{-2}$ s$^{-1}$ during the entire observation period and is comparable with the coronal heating rate requirement in active regions.

preprint2014arXiv

The Compatibility of Flare Temperatures Observed with AIA, GOES, and RHESSI

We test the compatibility and biases of multi-thermal flare DEM (differential emission measure) peak temperatures determined with AIA with those determined by GOES and RHESSI using the isothermal assumption. In a set of 149 M- and X-class flares observed during the first two years of the SDO mission, AIA finds DEM peak temperatures at the time of the peak GOES 1-8 A flux to have an average of Tp = 12.0+/-2.9 MK and Gaussian DEM widths of log10(sigma_T) = 0.50+/-0.13. From GOES observations of the same 149 events, a mean temperature of Tp = 15.6+/-2.4 MK is inferred, which is systematically higher by a factor of TGOES/TAIA = 1.4+/-0.4. We demonstrate that this discrepancy results from the isothermal assumption in the inversion of the GOES filter ratio. From isothermal fits to photon spectra at energies of E ~ 6-12 keV of 61 of these events, RHESSI finds the temperature to be higher still by a factor of TRHESSI/TAIA = 1.9+/-1.0. We find that this is partly a consequence of the isothermal assumption. However, RHESSI is not sensitive to the low-temperature range of the DEM peak, and thus RHESSI samples only the high-temperature tail of the DEM function. This can also contribute to the discrepancy between AIA and RHESSI temperatures. The higher flare temperatures found by GOES and RHESSI imply correspondingly lower emission measures. We conclude that self-consistent flare DEM temperatures and emission measures require simultaneous fitting of EUV (AIA) and soft X-ray (GOES and RHESSI) fluxes.

preprint2014arXiv

The Magnetic Field of Active Region 11158 During the 2011 February 12-17 Flares : Differences between Photospheric Extrapolation and Coronal Forward-Fitting Methods

We developed a {\sl coronal non-linear force-free field (COR-NLFFF)} forward-fitting code that fits an approximate {\sl non-linear force-free field (NLFFF)} solution to the observed geometry of automatically traced coronal loops. In contrast to photospheric NLFFF codes, which calculate a magnetic field solution from the constraints of the transverse photospheric field, this new code uses coronal constraints instead, and this way provides important information on systematic errors of each magnetic field calculation method, as well as on the non-forcefreeness in the lower chromosphere. In this study we applied the COR-NLFFF code to active region NOAA 11158, during the time interval of 2011 Feb 12 to 17, which includes an X2.2 GOES-class flare plus 35 M and C-class flares. We calcuated the free magnetic energy with a 6-minute cadence over 5 days. We find good agreement between the two types of codes for the total nonpotential $E_N$ and potential energy $E_P$, but find up to a factor of 4 discrepancy in the free energy $E_{free}=E_N-E_P$, and up to a factor of 10 discrepancy in the decrease of the free energy $ΔE_{free}$ during flares. The coronal NLFFF code exhibits a larger time variability, and yields a decrease of free energy during the flare that is sufficient to satisfy the flare energy budget, while the photospheric NLFFF code shows much less time variability and an order of magnitude less free energy decrease during flares. The discrepancy may partly be due to the pre-processing of photospheric vector data, but more likely due to the non-forcefreeness in the lower chromosphere. We conclude that the coronal field cannot be correctly calculated based on photospheric data alone, but requires additional information on coronal loop geometries.

preprint2013arXiv

A Macroscopic Description of a Generalized Self-Organized Criticality System: Astrophysical Applications

We suggest a generalized definition of self-organized criticality (SOC) systems: SOC is a critical state of a nonlinear energy dissipation system that is slowly and continuously driven towards a critical value of a system-wide instability threshold, producing scale-free, fractal-diffusive, and intermittent avalanches with powerlaw-like size distributions. We develop here a macroscopic description of SOC systems that provides an equivalent description of the complex microscopic fine structure, in terms of fractal-diffusive transport (FD-SOC). Quantitative values for the size distributions of SOC parameters (length scales $L$, time scales $T$, waiting times $Δt$, fluxes $F$, and energies $E$) are derived from first principles, using the scale-free probability conjecture, $N(L) dL \propto L^{-d}$, for Euclidean space dimension $d$. We apply this model to astrophysical SOC systems, such as lunar craters, the asteroid belt, Saturn ring particles, magnetospheric substorms, radiation belt electrons, solar flares, stellar flares, pulsar glitches, soft gamma-ray repeaters, black-hole objects, blazars, and cosmic rays. The FD-SOC model predicts correctly the size distributions of 8 out of these 12 astrophysical phenomena, and indicates non-standard scaling laws and measurement biases for the others.

preprint2013arXiv

Multi-Wavelength Observations of the Spatio-Temporal Evolution of Solar Flares with AIA/SDO: I. Universal Scaling Laws of Space and Time Parameters

We extend a previous statistical solar flare study of 155 GOES M- and X-class flares observed with AIA/SDO (Aschwanden 2012) to all 7 coronal wavelengths (94, 131, 171, 193, 211, 304, 335 \ang) to test the wavelength-dependence of scaling laws and statistical distributions. Except for the 171 and 193 \ang\ wavelengths, which are affected by EUV dimming caused by coronal mass ejections (CMEs), we find near-identical size distributions of geometric (lengths $L$, flare areas $A$, volumes $V$, fractal dimension $D_2$), temporal (flare durations $T$), and spatio-temporal parameters (diffusion coefficient $κ$, spreading exponent $β$, and maximum expansion velocities $v_{max}$) in different wavelengths, which are consistent with the universal predictions of the fractal-diffusive avalanche model of a slowly-driven self-organized criticality (FD-SOC) system, i.e., $N(L) \propto L^{-3}$, $N(A) \propto A^{-2}$, $N(V) \propto V^{-5/3}$, $N(T) \propto T^{-2}$, $D_2=3/2$, for a Euclidean dimension $d=3$. Empirically we find also a new strong correlation $κ\propto L^{0.94\pm0.01}$ and the 3-parameter scaling law $L \propto κ T^{0.1}$, which is more consistent with the logistic-growth model than with classical diffusion. The findings suggest long-range correlation lengths in the FD-SOC system that operate in the vicinity of a critical state, which could be used for predictions of individual extreme events. We find also that eruptive flares (with accompanying CMEs), have larger volumes $V$, longer flare durations $T$, higher EUV and soft X-ray fluxes, and somewhat larger diffusion coefficients $κ$ than confined flares (without CMEs).

preprint2013arXiv

Multi-Wavelength Observations of the Spatio-Temporal Evolution of Solar Flares with AIA/SDO: II. Hydrodynamic Scaling Laws and Thermal Energies

In this study we measure physical parameters of the same set of 155 M and X-class solar flares observed with AIA/SDO as analyzed in Paper I, by performing a {\sl differential emission measure (DEM)} analysis to determine the flare peak emission measure $EM_p$, peak temperature $T_p$, electron density $n_p$, and thermal energy $E_{th}$, in addition to the spatial scales $L$, areas $A$, and volumes $V$ measured in Paper I. The parameter ranges for M and X-class flares are: $\log(EM_p)=47.0-50.5$, $T_p=5.0-17.8$ MK, $n_p=4 \times 10^9-9 \times 10^{11}$ cm$^{-3}$, and thermal energies of $E_{th}=1.6 \times 10^{28}-1.1 \times 10^{32}$ erg. We find that these parameters obey the Rosner-Tucker-Vaiana (RTV) scaling law $T_p^2 \propto n_p L$ and $H \propto T^{7/2} L^{-2}$ during the peak time $t_p$ of the flare density $n_p$, when energy balance between the heating rate $H$ and the conductive and radiative loss rates is achieved for a short instant, and thus enables the applicability of the RTV scaling law. The application of the RTV scaling law predicts powerlaw distributions for all physical parameters, which we demonstrate with numerical Monte-Carlo simulations as well as with analytical calculations. A consequence of the RTV law is also that we can retrieve the size distribution of heating rates, for which we find $N(H) \propto H^{-1.8}$, which is consistent with the magnetic flux distribution $N(Φ) \propto Φ^{-1.85}$ observed by Parnell et al.(2009) and the heating flux scaling law $F_H \propto H L \propto B/L$ of Schrijver et al.(2004). The fractal-diffusive self-organized criticality model in conjunction with the RTV scaling law reproduces the observed powerlaw distributions and their slopes for all geometrical and physical parameters and can be used to predict the size distributions for other flare datasets, instruments, and detection algorithms.

preprint2013arXiv

Optimization of Curvi-Linear Tracing Applied to Solar Physics and Biophysics

We developed an automated pattern recognition code that is particularly well suited to extract one-dimensional curvi-linear features from two-dimensional digital images. A former version of this {\sl Oriented Coronal CUrved Loop Tracing (OCCULT)} code was applied to spacecraft images of magnetic loops in the solar corona, recorded with the NASA spacecraft {\sl Transition Region And Coronal Explorer (TRACE)} in extreme ultra-violet wavelengths. Here we apply an advanced version of this code ({\sl OCCULT-2}) also to similar images from the {\sl Solar Dynamics Observatory (SDO)}, to chromospheric H-$α$ images obtained with the {\sl Swedish Solar Telescope (SST)}, and to microscopy images of microtubule filaments in live cells in biophysics. We provide a full analytical description of the code, optimize the control parameters, and compare the automated tracing with visual/manual methods. The traced structures differ by up to 16 orders of magnitude in size, which demonstrates the universality of the tracing algorithm.

preprint2012arXiv

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. I. Theory

We derive an analytical approximation of nonlinear force-free magnetic field solutions (NLFFF) that can efficiently be used for fast forward-fitting to solar magnetic data, constrained either by observed line-of-sight magnetograms and stereoscopically triangulated coronal loops, or by 3D vector-magnetograph data. The derived NLFFF solutions provide the magnetic field components $B_x({\bf x})$, $B_y({\bf x})$, $B_z({\bf x})$, the force-free parameter $α({\bf x})$, the electric current density ${\bf j}({\bf x})$, and are accurate to second-order (of the nonlinear force-free $α$-parameter). The explicit expressions of a force-free field can easily be applied to modeling or forward-fitting of many coronal phenomena.

preprint2012arXiv

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. II. Numeric Code and Tests

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.

preprint2012arXiv

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. III. The Free Energy

An analytical approximation of a nonlinear force-free magnetic field (NLFFF) solution was developed in Paper I, while a numerical code that performs fast forward-fitting of this NLFFF approximation to a line-of-sight magnetogram and coronal 3D loops has been described and tested in Paper II. Here we calculate the free magnetic energy $E_{\rm free}=E_{\rm N}-E_{\rm P}$, i.e. the difference of the magnetic energies between the nonpotential field and the potential field. A second method to estimate the free energy is obtained from the mean misalignment angle change $Δμ=μ_{\rm P}-μ_{\rm N}$ between the potential and nonpotential field, which scales as $E_{\rm free}/E_{\rm P} \approx \tan^2{(Δμ)}$. For four active regions observed with STEREO in 2007 we find free energies in the range of $q_{\rm free}=(E_{\rm free}/E_{\rm P}) \approx 1%-10%$, with an uncertainty of less than $\pm 2%$ between the two methods, while the free energies obtained from 11 other NLFFF codes exhibit a larger scatter of order $\approx\pm10%$. We find also a correlation between the free magnetic energy and the GOES flux of the largest flare that occurred during the observing period, which can be quantified by an exponential relationship, $F_{\rm GOES} \propto \exp{(q_{\rm free}/0.015)}$, implying an exponentiation of the dissipated currents.

preprint2012arXiv

Automated Solar Flare Statistics in Soft X-rays over 37 Years of GOES Observations - The Invariance of Self-Organized Criticality during Three Solar Cycles

We analyzed the soft X-ray light curves from the {\sl Geostationary Operational Environmental Satellites (GOES)} over the last 37 years (1975-2011) and measured with an automated flare detection algorithm over 300,000 solar flare events (amounting to $\approx 5$ times higher sensitivity than the NOAA flare catalog). We find a powerlaw slope of $α_F=1.98\pm0.11$ for the (background-subtracted) soft X-ray peak fluxes that is invariant through three solar cycles and agrees with the theoretical prediction $α_F=2.0$ of the {\sl fractal-diffusive self-organized criticality (FD-SOC)} model. For the soft X-ray flare rise times we find a powerlaw slope of $α_T =2.02\pm0.04$ during solar cycle minima years, which is also consistent with the prediction $α_T=2.0$ of the FD-SOC model. During solar cycle maxima years, the powerlaw slope is steeper in the range of $α_T \approx 2.0-5.0$, which can be modeled by a solar cycle-dependent flare pile-up bias effect. These results corroborate the FD-SOC model, which predicts a powerlaw slope of $α_E=1.5$ for flare energies and thus rules out significant nanoflare heating. While the FD-SOC model predicts the probability distribution functions of spatio-temporal scaling laws of nonlinear energy dissipation processes, additional physical models are needed to derive the scaling laws between the geometric SOC parameters and the observed emissivity in different wavelength regimes, as we derive here for soft X-ray emission. The FD-SOC model yields also statistical probabilities for solar flare forecasting.

preprint2012arXiv

First 3D Reconstructions of Coronal Loops with the STEREO A+B Spacecraft: IV. Magnetic Modeling with Twisted Force-Free Fields

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.

preprint2012arXiv

Flares in the Crab Nebula Driven by Untwisting Magnetic Fields

The recent discovery of PeV electrons from the Crab nebula, produced on rapid time scales of one day or less with a sharply peaked gamma-ray spectrum without hard X-rays, challenges traditional models of diffusive shock acceleration followed by synchrotron radiation. Here we outline an accleration model involving a DC electric field parallel to the magnetic field in a twisted toroidal field around the pulsar. Sudden developments of resistivity in localized regions of the twisted field are thought to drive the particle acceleration, up to PeV energies, resulting in flares. This model can reproduce the observed time scales of $T \approx 1$ day, the peak photon energies of $U_{Φ,rr} \approx 1$ MeV, maximum electron energies of $U_{e,rr} \approx 1$ PeV, and luminosities of $L \approx 10^{36}$ erg s$^{-1}$.

preprint2012arXiv

Nonlinear Force-Free Magnetic Field Fitting to Coronal Loops with and without Stereoscopy

We developed a new nonlinear force-free magnetic field (NLFFF) forward-fitting algorithm based on an analytical approximation of force-free and divergence-free NLFFF solutions, which requires as input a line-of-sight magnetogram and traced 2D loop coordinates of coronal loops only, in contrast to stereoscopically triangulated 3D loop coordinates used in previous studies. Test results of simulated magnetic configurations and from four active regions observed with STEREO demonstrate that NLFFF solutions can be fitted with equal accuracy with or without stereoscopy, which relinquishes the necessity of STEREO data for magnetic modeling of active regions (on the solar disk). The 2D loop tracing method achieves a 2D misalignment of $μ_2=2.7^\circ\pm 1.3^\circ$ between the model field lines and observed loops, and an accuracy of $\approx 1.0%$ for the magnetic energy or free magnetic energy ratio. The three times higher spatial resolution of TRACE or SDO/AIA (compared with STEREO) yields also a proportionally smaller misalignment angle between model fit and observations. Visual/manual loop tracings are found to produce more accurate magnetic model fits than automated tracing algorithms. The computation time of the new forward-fitting code amounts to a few minutes per active region.

preprint2012arXiv

Quasi-periodic Fast-mode Wave Trains Within a Global EUV Wave and Sequential Transverse Oscillations Detected by SDO/AIA

We present the first unambiguous detection of quasi-periodic wave trains within the broad pulse of a global EUV wave (so-called "EIT wave") occurring on the limb. These wave trains, running ahead of the lateral CME front of 2-4 times slower, coherently travel to distances $>R_{sun}/2$ along the solar surface, with initial velocities up to 1400 km/s decelerating to ~650 km/s. The rapid expansion of the CME initiated at an elevated height of 110 Mm produces a strong downward and lateral compression, which may play an important role in driving the primary EUV wave and shaping its front forwardly inclined toward the solar surface. The waves have a dominant 2 min periodicity that matches the X-ray flare pulsations, suggesting a causal connection. The arrival of the leading EUV wave front at increasing distances produces an uninterrupted chain sequence of deflections and/or transverse (likely fast kink mode) oscillations of local structures, including a flux-rope coronal cavity and its embedded filament with delayed onsets consistent with the wave travel time at an elevated (by ~50%) velocity within it. This suggests that the EUV wave penetrates through a topological separatrix surface into the cavity, unexpected from CME caused magnetic reconfiguration. These observations, when taken together, provide compelling evidence of the fast-mode MHD wave nature of the {\it primary (outer) fast component} of a global EUV wave, running ahead of the {\it secondary (inner) slow} component of CME-caused restructuring.

preprint2012arXiv

Self-Organized Criticality Systems in Astrophysics (Chapter 13)

Chapter 13: SOC Systems in Astrophysics --- Content list: 13.1 Theory -- 13.1.1 The Sacle-Free Probability Theorem - 13.1.2 The Fractal-Diffusive Spatio-Temporal Relationship - 13.1.3 Size Distributions of Astrophysical Observables - 13.1.4 Scaling Laws for Thermal Emission of Astrophysical Plasmas - 13.1.5 Scaling Laws for Astrophysical Acceleration Mechanisms - 13.2 Observations -- 13.2.1 Lunar Craters - 13.2.2 Asteroid Belt - 13.2.3 Saturn Ring - 13.2.4 Magnetospheric Substorms and Auroras - 13.2.5 Solar Flares - 13.2.6 Stellar Flares - 13.2.7 Pulsars - 13.2.8 Soft Gamma-Ray Repeaters - 13.2.9 Black-Hole Objects - 13.2.10 Blazars - 13.2.11 Cosmic Rays - 13.3 Conclusions

preprint2012arXiv

Solar Stereoscopy with STEREO/EUVI A and B spacecraft from small (6 deg) to large (170 deg) spacecraft separation angles

We performed for the first time stereoscopic triangulation of coronal loops in active regions over the entire range of spacecraft separation angles ($α_{sep}\approx 6^\circ, 43^\circ, 89^\circ, 127^\circ$, and $170^\circ$). The accuracy of stereoscopic correlation depends mostly on the viewing angle with respect to the solar surface for each spacecraft, which affects the stereoscopic correspondence identification of loops in image pairs. From a simple theoretical model we predict an optimum range of $α_{sep} \approx 22^\circ-125^\circ$, which is also experimentally confirmed. The best accuracy is generally obtained when an active region passes the central meridian (viewed from Earth), which yields a symmetric view for both STEREO spacecraft and causes minimum horizontal foreshortening. For the extended angular range of $α_{sep}\approx 6^\circ-127^{\circ}$ we find a mean 3D misalignment angle of $μ_{PF} \approx 21^\circ-39^\circ$ of stereoscopically triangulated loops with magnetic potential field models, and $μ_{FFF} \approx 15^\circ-21^\circ$ for a force-free field model, which is partly caused by stereoscopic uncertainties $μ_{SE} \approx 9^\circ$. We predict optimum conditions for solar stereoscopy during the time intervals of 2012--2014, 2016--2017, and 2021--2023.

preprint2012arXiv

The Spatio-Temporal Evolution of Solar Flares Observed with AIA/SDO: Fractal Diffusion, Sub-Diffusion, or Logistic Growth ?

We explore the spatio-temporal evolution of solar flares by fitting a radial expansion model $r(t)$ that consists of an exponentially growing acceleration phase, followed by a deceleration phase that is parameterized by the generalized diffusion function $r(t) \propto κ(t-t_1)^{β/2}$, which includes the logistic growth limit ($β=0$), sub-diffusion ($β= 0-1$), classical diffusion ($β=1$), super-diffusion ($β= 1-2$), and the linear expansion limit ($β=2$). We analyze all M and X-class flares observed with GOES and AIA/SDO during the first two years of the SDO mission, amounting to 155 events. We find that most flares operate in the sub-diffusive regime ($β=0.53\pm0.27$), which we interpret in terms of anisotropic chain reactions of intermittent magnetic reconnection episodes in a low plasma-$β$ corona. We find a mean propagation speed of $v=15\pm12$ km s$^{-1}$, with maximum speeds of $v_{max}=80 \pm 85$ km s$^{-1}$ per flare, which is substantially slower than the sonic speeds expected for thermal diffusion of flare plasmas. The diffusive characteristics established here (for the first time for solar flares) is consistent with the fractal-diffusive self-organized criticality (FD-SOC) model, which predicted diffusive transport merely based on cellular automaton simulations.

preprint2012arXiv

Theoretical Models of Self-Organized Criticality (SOC) Systems

In this chapter 2 of the e-book "Self-Organized Criticality Systems" we summarize the classical cellular automaton models, which consist of a statistical aspect that is universal to all SOC systems, and a physical aspect that depends on the physical definition of the observable. Then we derive some general analytical formulations of SOC processes, such as the exponential-growth SOC model and the fractal-diffusive SOC model, which also have universal validity for SOC processes, while specific applications to observations require additional physical scaling laws (e.g., for astrophysical or geophysical observations). Finally we discuss alternative SOC processes, SOC-related, or non-SOC processes, such as: self-organization (without criticality), forced SOC model, Brownian motion or classical diffusion, hyper-diffusion and Levy flight, nonextensive Tsallis entropy, turbulence, percolation, phase transitions, network systems, and chaotic systems. We synthesize a metrics that specifies which observational SOC properties are shared by these processes.

preprint2011arXiv

A Statistical Fractal-Diffusive Avalanche Model of a Slowly-Driven Self-Organized Criticality System

We develop a statistical analytical model that predicts the occurrence frequency distributions and parameter correlations of avalanches in nonlinear dissipative systems in the state of a slowly-driven self-organized criticality (SOC) system. This model, called the fractal-diffusive SOC model, is based on the following four assumptions: (i) The avalanche size $L$ grows as a diffusive random walk with time $T$, following $L \propto T^{1/2}$; (ii) The instantaneous energy dissipation rate $f(t)$ occupies a fractal volume with dimension $D_S$, which predicts the relationships $F = f(t=T) \propto L^{D_S} \propto T^{D_S/2}$, $P \propto L^{S} \propto T^{S/2}$ for the peak energy dissipation rate, and $E \propto F T \propto T^{1+D_S/2}$ for the total dissipated energy; (iii) The mean fractal dimension of avalanches in Euclidean space $S=1,2,3$ is $D_S \approx (1+S)/2$; and (iv) The occurrence frequency distributions $N(x) \propto x^{-α_x}$ based on spatially uniform probabilities in a SOC system are given by $N(L) \propto L^{-S}$, which predicts powerlaw distributions for all parameters, with the slopes $α_T=(1+S)/2$, $α_F=1+(S-1)/D_S$, $α_P=2-1/S$, and $α_E=1+(S-1)/(D_S+2)$. We test the predicted fractal dimensions, occurrence frequency distributions, and correlations with numerical simulations of cellular automaton models in three dimensions $S=1,2,3$ and find satisfactory agreement within $\approx 10%$. One profound prediction of this universal SOC model is that the energy distribution has a powerlaw slope in the range of $α_E=1.40-1.67$, and the peak energy distribution has a slope of $α_P=1.67$ (for any fractal dimension $D_S=1,...,3$ in Euclidean space S=3), and thus predicts that the bulk energy is always contained in the largest events, which rules out significant nanoflare heating in the case of solar flares.

preprint2011arXiv

Coronal Loop Oscillations Observed with AIA - Kink-Mode with Cross-Sectional and Density Oscillations

A detailed analysis of a coronal loop oscillation event is presented, using data from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) for the first time. The loop oscillation event occurred on 2010 Oct 16, 19:05-19:35 UT, was triggered by an M2.9 GOES-class flare, located inside a highly inclined cone of a narrow-angle CME. This oscillation event had a number of unusual features: (i) Excitation of kink-mode oscillations in vertical polarization (in the loop plane); (ii) Coupled cross-sectional and density oscillations with identical periods; (iii) no detectable kink amplitude damping over the observed duration of four kink-mode periods ($P=6.3$ min); (iv) multi-loop oscillations with slightly ($\approx 10%$) different periods; and (v) a relatively cool loop temperature of $T\approx 0.5$ MK. We employ a novel method of deriving the electron density ratio external and internal to the oscillating loop from the ratio of Alfvénic speeds deduced from the flare trigger delay and the kink-mode period, i.e., $n_e/n_i=(v_A/v_{Ae})^2=0.08\pm0.01$. The coupling of the kink mode and cross-sectional oscillations can be explained as a consequence of the loop length variation in the vertical polarization mode. We determine the exact footpoint locations and loop length with stereoscopic triangulation using STEREO/EUVI-A data. We model the magnetic field in the oscillating loop using HMI/SDO magnetogram data and a potential field model and find agreement with the seismological value of the magnetic field, $B_{kink}=4.0\pm0.7$ G, within a factor of two.

preprint2011arXiv

Direct Imaging by SDO AIA of Quasi-periodic Fast Propagating Waves of ~2000 km/s in the Low Solar Corona

Quasi-periodic, propagating fast mode magnetosonic waves in the corona were difficult to observe in the past due to relatively low instrument cadences. We report here evidence of such waves directly imaged in EUV by the new SDO AIA instrument. In the 2010 August 1 C3.2 flare/CME event, we find arc-shaped wave trains of 1-5% intensity variations (lifetime ~200 s) that emanate near the flare kernel and propagate outward up to ~400 Mm along a funnel of coronal loops. Sinusoidal fits to a typical wave train indicate a phase velocity of 2200 +/- 130 km/s. Similar waves propagating in opposite directions are observed in closed loops between two flare ribbons. In the k-$ω$ diagram of the Fourier wave power, we find a bright ridge that represents the dispersion relation and can be well fitted with a straight line passing through the origin. This k-$ω$ ridge shows a broad frequency distribution with indicative power at 5.5, 14.5, and 25.1 mHz. The strongest signal at 5.5 mHz (period 181 s) temporally coincides with quasi-periodic pulsations of the flare, suggesting a common origin. The instantaneous wave energy flux of $(0.1-2.6) \times 10^7 ergs/cm^2/s$ estimated at the coronal base is comparable to the steady-state heating requirement of active region loops.

preprint2011arXiv

Energy-Dependent Timing of Thermal Emission in Solar Flares

We report solar flare plasma to be multi-thermal in nature based on the theoretical model and study of the energy-dependent timing of thermal emission in ten M-class flares. We employ high-resolution X-ray spectra observed by the Si detector of the "Solar X-ray Spectrometer" (SOXS). The SOXS onboard the Indian GSAT-2 spacecraft was launched by the GSLV-D2 rocket on 8 May 2003. Firstly we model the spectral evolution of the X-ray line and continuum emission flux F(ε) from the flare by integrating a series of isothermal plasma flux. We find that multi-temperature integrated flux F(ε) is a power-law function of εwith a spectral index (γ) \approx -4.65. Next, based on spectral-temporal evolution of the flares we find that the emission in the energy range E= 4 - 15 keV is dominated by temperatures of T= 12 - 50 MK, while the multi-thermal power-law DEM index (γ) varies in the range of -4.4 and -5.7. The temporal evolution of the X-ray flux F(ε,t) assuming a multi-temperature plasma governed by thermal conduction cooling reveals that the temperature-dependent cooling time varies between 296 and 4640 s and the electron density (n_e) varies in the range of n_e= (1.77-29.3)*10^10 cm-3. Employing temporal evolution technique in the current study as an alternative method for separating thermal from non-thermal components in the energy spectra, we measure the break-energy point ranging between 14 and 21\pm1.0 keV.

preprint2011arXiv

GeV Particle Acceleration in Solar Flares and Ground Level Enhancement (GLE) Events

{\sl Ground Level Enhancement (GLE)} events represent the largest class of {\sl solar energetic particle (SEP)} events that require acceleration processes to produce $\gapprox 1$ GeV ions in order to produce showers of secondary particles in the Earth's atmosphere with sufficient intensity to be detected by ground-level neutron monitors, above the background of cosmics rays. Although the association of GLE events with both solar flares and coronal mass ejections (CMEs) is undisputed, the question arises about the location of the responsible acceleration site: coronal flare sites or heliospheric CME-associated shocks? To investigate the first possibility we explore the timing of GLE events with respect to hard X-ray production in solar flares, considering the height and magnetic topology of flares, the role of extended acceleration, and particle trapping. We find that 50% (6 out of 12) of recent (non-occulted) GLE events are accelerated during the impulsive flare phase, while the remaining half are accelerated significantly later. It appears that the prompt GLE component, which is observed in virtually all GLE events, is caused by flare-accelerated particles in the lower corona, while the delayed gradual GLE component can be produced by both, either by extended acceleration and/or trapping in flare sites, or by particles accelerated in CME-associated shocks during their propagation through the heliosphere.

preprint2011arXiv

Implications of X-ray Observations for Electron Acceleration and Propagation in Solar Flares

High-energy X-rays and gamma-rays from solar flares were discovered just over fifty years ago. Since that time, the standard for the interpretation of spatially integrated flare X-ray spectra at energies above several tens of keV has been the collisional thick-target model. After the launch of the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) in early 2002, X-ray spectra and images have been of sufficient quality to allow a greater focus on the energetic electrons responsible for the X-ray emission, including their origin and their interactions with the flare plasma and magnetic field. The result has been new insights into the flaring process, as well as more quantitative models for both electron acceleration and propagation, and for the flare environment with which the electrons interact. In this article we review our current understanding of electron acceleration, energy loss, and propagation in flares. Implications of these new results for the collisional thick-target model, for general flare models, and for future flare studies are discussed.

preprint2011arXiv

Solar Corona Loop Studies with AIA: I. Cross-Sectional Temperature Structure

We present a first systematic study on the cross-sectional temperature structure of coronal loops using the six coronal temperature filters of the Atmospheric Imaging Assembly (AIA) instrument on the Solar Dynamics Observatory (SDO). We analyze a sample of 100 loop snapshots measured at 10 different locations and 10 different times in active region NOAA 11089 on 2010 July 24, 21:00-22:00 UT. The cross-sectional flux profiles are measured and a cospatial background is subtracted in 6 filters in a temperature range of $T \approx 0.5-16$ MK, and 4 different parameterizations of differential emission measure (DEM) distributions are fitted. We find that the reconstructed DEMs consist predominantly of narrowband peak temperature components with a thermal width of $σ_{log(T)} \le 0.11\pm0.02$, close to the temperature resolution limit of the instrument, consistent with earlier triple-filter analysis from TRACE by Aschwanden and Nightingale (2005) and from EIS/Hinode by Warren et al. (2008) or Tripathi et al. (2009). We find that 66% of the loops could be fitted with a narrowband single-Gaussian DEM model, and 19% with a DEM consisting of two narrowband Gaussians (which mostly result from pairs of intersecting loops along the same line-of-sight). The mostly isothermal loop DEMs allow us also to derive an improved empirical response function of the AIA 94 \ang\ filter, which needs to be boosted by a factor of $q_{94} = 6.7\pm 1.7$ for temperatures at $log(T) \lapprox 6.3$. The main result of near-isothermal loop cross-sections is not consistent with the predictions of standard nanoflare scenarios, but can be explained by flare-like heating mechanisms that drive chromospheric evaporation and upflows of heated plasma coherently over loop cross-sections of $w \approx 2-4$ Mm.

preprint2010arXiv

A Universal Scaling Law for the Fractal Energy Dissipation Domain in Self-Organized Criticality Systems

Nonlinear dissipative systems in the state of self-organized criticality release energy sporadically in avalanches of all sizes, such as in earthquakes, auroral substorms, solar and stellar flares, soft gamma-ray repeaters, and pulsar glitches. The statistical occurrence frequency distributions of event energies $E$ generally exhibit a powerlaw-like function $N(E)\propto E^{-α_E}$ with a powerlaw slope of $α_E \approx 1.5$. The powerlaw slope $α_E$ of energies can be related to the fractal dimension $D$ of the spatial energy dissipation domain by $D=3/α_E$, which predicts a powerlaw slope $α_E=1.5$ for area-rupturing or area-spreading processes with $D=2$. For solar and stellar flares, 2-D area-spreading dissipation domains are naturally provided in current sheets or separatrix surfaces in a magnetic reconnection region. Thus, this universal scaling law provides a useful new diagnostic on the topology of the spatial energy dissipation domain in geophysical and astrophysical observations.

preprint2010arXiv

Bootstrapping the Coronal Magnetic Field with STEREO: I. Unipolar Potential Field Modeling

We investigate the recently quantified misalignment of $α_{mis} \approx 20^\circ-40^\circ$ between the 3-D geometry of stereoscopically triangulated coronal loops observed with STEREO/EUVI (in four active regions) and theoretical (potential or nonlinear force-free) magnetic field models extrapolated from photospheric magnetograms. We develop an efficient method of bootstrapping the coronal magnetic field by forward-fitting a parameterized potential field model to the STEREO-observed loops. The potential field model consists of a number of unipolar magnetic charges that are parameterized by decomposing a photospheric magnetogram from MDI. The forward-fitting method yields a best-fit magnetic field model with a reduced misalignment of $α_{PF} \approx 13^\circ-20^\circ$. We evaluate also stereoscopic measurement errors and find a contribution of $α_{SE}\approx 7^\circ-12^\circ$, which constrains the residual misalignment to $α_{NP}=α_{PF}-α_{SE}\approx 5^\circ -9^\circ$, which is likely due to the nonpotentiality of the active regions. The residual misalignment angle $α_{NP}$ of the potential field due to nonpotentiality is found to correlate with the soft X-ray flux of the active region, which implies a relationship between electric currents and plasma heating.

preprint2010arXiv

Reconciliation of Waiting Time Statistics of Solar Flares Observed in Hard X-Rays

We study the waiting time distributions of solar flares observed in hard X-rays with ISEE-3/ICE, HXRBS/SMM, WATCH/GRANAT, BATSE/CGRO, and RHESSI. Although discordant results and interpretations have been published earlier, based on relatively small ranges ($< 2$ decades) of waiting times, we find that all observed distributions, spanning over 6 decades of waiting times ($Δt \approx 10^{-3}- 10^3$ hrs), can be reconciled with a single distribution function, $N(Δt) \propto λ_0 (1 + λ_0 Δt)^{-2}$, which has a powerlaw slope of $p \approx 2.0$ at large waiting times ($Δt \approx 1-1000$ hrs) and flattens out at short waiting times $Δt \lapprox Δt_0 = 1/λ_0$. We find a consistent breakpoint at $Δt_0 = 1/λ_0 = 0.80\pm0.14$ hours from the WATCH, HXRBS, BATSE, and RHESSI data. The distribution of waiting times is invariant for sampling with different flux thresholds, while the mean waiting time scales reciprocically with the number of detected events, $Δt_0 \propto 1/n_{det}$. This waiting time distribution can be modeled with a nonstationary Poisson process with a flare rate $λ=1/Δt$ that varies as $f(λ) \propto λ^{-1} \exp{-(λ/λ_0)}$. This flare rate distribution represents a highly intermittent flaring productivity in short clusters with high flare rates, separated by quiescent intervals with very low flare rates.

preprint2010arXiv

Self-Organized Criticality in Solar Physics and Astrophysics

The concept of "self-organized criticality" (SOC) has been introduced by Bak, Tang, and Wiesenfeld (1987) to describe the statistics of avalanches on the surface of a sandpile with a critical slope, which produces a scale-free powerlaw size distribution of avalanches. In the meantime, SOC behavior has been identified in many nonlinear dissipative systems that are driven to a critical state. On a most general level, SOC is the statistics of coherent nonlinear processes, in contrast to the Poisson statistics of incoherent random processes. The SOC concept has been applied to laboratory experiments (of rice or sand piles), to human activities (population growth, language, economy, traffic jams, wars), to biophysics, geophysics (earthquakes, landslides, forest fires), magnetospheric physics, solar physics (flares), stellar physics (flares, cataclysmic variables, accretion disks, black holes, pulsar glitches, gamma ray bursts), and to galactic physics and cosmology.

preprint2010arXiv

The State of Self-Organized Criticality of the Sun During the Last 3 Solar Cycles. I. Observations

We analyze the occurrence frequency distributions of peak fluxes $P$, total fluxes $E$, and durations $T$ of solar flares over the last three solar cycles (during 1980--2010) from hard X-ray data of HXRBS/SMM, BATSE/CGRO, and RHESSI. From the synthesized data we find powerlaw slopes with mean values of $α_P=1.72\pm0.08$ for the peak flux, $α_E=1.60\pm0.14$ for the total flux, and $α_T=1.98\pm0.35$ for flare durations. We find a systematic anti-correlation of the powerlaw slope of peak fluxes as a function of the solar cycle, varying with an approximate sinusoidal variation $α_P(t)=α_0+Δα\cos{[2π(t-t_0)/T_{cycle}]}$, with a mean of $α_0=1.73$, a variation of $Δα=0.14$, a solar cycle period $T_{cycle}=12.6$ yrs, and a cycle minimum time $t_0=1984.1$. The powerlaw slope is flattest during the maximum of a solar cycle, which indicates a higher magnetic complexity of the solar corona that leads to an overproportional rate of powerful flares.

preprint2010arXiv

The State of Self-Organized Criticality of the Sun During the Last Three Solar Cycles. II. Theoretical Model

The observed powerlaw distributions of solar flare parameters can be interpreted in terms of a nonlinear dissipative system in the state of self-organized criticality (SOC). We present a universal analytical model of a SOC process that is governed by three conditions: (i) a multiplicative or exponential growth phase, (ii) a randomly interrupted termination of the growth phase, and (iii) a linear decay phase. This basic concept approximately reproduces the observed frequency distributions. We generalize it to a randomized exponential-growth model, which includes also a (log-normal) distribution of threshold energies before the instability starts, as well as randomized decay times, which can reproduce both the observed occurrence frequency distributions and the scatter of correlated parametyers more realistically. With this analytical model we can efficiently perform Monte-Carlo simulations of frequency distributions and parameter correlations of SOC processes, which are simpler and faster than the iterative simulations of cellular automaton models. Solar cycle modulations of the powerlaw slopes of flare frequency distributions can be used to diagnose the thresholds and growth rates of magnetic instabilities responsible for solar flares.

preprint2009arXiv

4D Modeling of CME expansion and EUV dimming observed with STEREO/EUVI

This is the first attempt to model the kinematics of a CME launch and the resulting EUV dimming quantitatively with a self-consistent model. Our 4D-model assumes self-similar expansion of a spherical CME geometry that consists of a CME front with density compression and a cavity with density rarefaction, satisfying mass conservation of the total CME and swept-up corona. The model contains 14 free parameters and is fitted to the 2008 March 25 CME event observed with STEREO/A and B. Our model is able to reproduce the observed CME expansion and related EUV dimming during the initial phase from 18:30 UT to 19:00 UT. The CME kinematics can be characterized by a constant acceleration (i.e., a constant magnetic driving force). While the observations of EUVI/A are consistent with a spherical bubble geometry, we detect significant asymmetries and density inhomogeneities with EUVI/B. This new forward-modeling method demonstrates how the observed EUV dimming can be used to model physical parameters of the CME source region, the CME geometry, and CME kinematics.

preprint2009arXiv

A Critical Assessment of Nonlinear Force-Free Field Modeling of the Solar Corona for Active Region 10953

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.

preprint2003arXiv

Tomographic 3D-Modeling of the Solar Corona with FASR

The "Frequency-Agile Solar Radiotelescope" (FASR) litteraly opens up a new dimension in addition to the 3D Euclidian geometry: the frequency dimension. The 3D geometry is degenerated to 2D in all images from astronomical telescopes, but the additional frequency dimension allows us to retrieve the missing third dimension by means of physical modeling. We call this type of 3D reconstruction "Frequency Tomography". In this study we simulate a realistic 3D model of an active region, composed of 500 coronal loops with the 3D geometry [x(s),y(s),z(s)] constrained by magnetic field extrapolations and the physical parameters of the density n_e(s) and temperature T_e(s) given by hydrostatic solutions. We simulate a series of 20 radio images in a frequency range of f=0.1-10 GHz, anticipating the capabilities of FASR, and investigate what physical information can be retrieved from such a dataset. We discuss also forward-modeling of the chromospheric and Quiet Sun density and temperature structure, another primary goal of future FASR science.