Researcher profile

Markus J. Aschwanden

Markus J. Aschwanden contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
24works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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

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

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 &#34;EIT wave&#34;) 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

Theoretical Models of Self-Organized Criticality (SOC) Systems

In this chapter 2 of the e-book &#34;Self-Organized Criticality Systems&#34; 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

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 &#34;Solar X-ray Spectrometer&#34; (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&#39;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.

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

Self-Organized Criticality in Solar Physics and Astrophysics

The concept of &#34;self-organized criticality&#34; (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

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.