In this study we test 30 variants of 5 physical scaling laws that describe different aspects of solar flares. We express scaling laws in terms of the magnetic potential field energy $E_p$, the mean potential field strength $B_p$, the free energy $E_{free}$, the dissipated magnetic flare energy $E_{diss}$, the mean loop length scale $L$, the mean helically twisted flux tube radius $R$, the sunspot radius $r$, the emission measure-weighted flare temperature $T_w$, the electron density $n_e$, and the total emission measure $EM$, measured from a data set of $\lapprox 400$ GOES M- and X-class flare events. The 5 categories of physical scaling laws include (i) a scaling law of the potential-field energy, (ii) a scaling law for helical twisting, (iii) a scaling law for Petschek-type magnetic reconnection, (iv) the Rosner-Tucker-Vaiana scaling law, and (v) the Shibata-Yokoyama scaling law. We test the self-consistency of these theoretical scaling laws with observed parameters by requiring two conditions: a cross-corrleation coefficient of CCC$>$0.5 between the observed and theoretically predicted scaling laws, and a linear regression fit with a slope of $\alpha \approx 1$. With these two criteria we find that 10 out of the 30 tested scaling law variants are consistent with the observed data, which strongly corroborates the existence and validity of the tested flare scaling laws.

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 $\lambda$ of the solar corona, the Alfv\'en Mach number $M_A=v_1/v_A$ (the ratio of the inflow speed $v_1$ to the Alfv\'enic outflow speed $v_A$), and the flare duration $\tau_f$, i.e., $E_{diss} = (1/4\pi) B_f^2\ L\ \lambda\ v_A\ M_A\ \tau_f$, where the Alfv\'en 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\'en Mach numbers $M_A$, the reconnection rate, and the total free energy dissipated in solar flares.

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

We report the discovery of torsional Alfv\'enic oscillations in solar flares, which modulate the time evolution of the magnetic free energy $E_f(t)$, while the magnetic potential energy $E_p(t)$ is uncorrelated, and the nonpotential energy varies as $E_{np}(t) = E_p + E_f(t)$. The mean observed time period of the torsional oscillations is $P_{obs}=15.1 \pm 3.9$ min, the mean field line length is $L=135\pm35$ Mm, and the mean phase speed is $v_{phase} =315 \pm 120$ km s$^{-1}$, which we interpret as torsional Alfv\'enic waves in flare loops with enhanced electron densities. Most of the torsional oscillations are found to be decay-less, but exhibit a positive or negative trend in the evolution of the free energy, indicating new emerging flux (if positive), magnetic cancellation, or flare energy dissipation (if negative). The time evolution of the free energy has been calculated in this study with the {\sl Vertical-Current Approximation (Version 4) Nonlinear Force-Free Field (VCA4-NLFFF)} code, which incorporates automatically detected coronal loops in the solution and bypasses the non-forcefreeness of the photospheric boundary condition, in contrast to traditional NLFFF codes.

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 $\Delta t$, with a slope of $\alpha_{\Delta 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.

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_\varphi)$ magnetic field components a poloidal component $(B_{\theta} \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 ($\tau_{flare} \ll \tau_{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.

The free energy that is dissipated in a magnetic reconnection process of a solar flare, generally accompanied by a coronal mass ejection (CME), has been considered as the ultimate energy source of the global energy budget of solar flares in previous statistical studies. Here we explore the effects of the aerodynamic drag force on CMEs, which supplies additional energy from the slow solar wind to a CME event, besides the magnetic energy supply. For this purpose we fit the analytical aerodynamic drag model of Cargill (2004) and Vrsnak et al.~{\bf (2013)} to the height-time profiles $r(t)$ of LASCO/SOHO data in 14,316 CME events observed during the first 8 years (2010-2017) of the SDO era {\bf (ensuring EUV coverage with AIA)}. Our main findings are: (i) a mean solar wind speed of $w=472 \pm 414$ km s$^{-1}$, (ii) a maximum drag-accelerated CME energy of $E_{drag} \lapprox 2 \times 10^{32}$ erg, (iii) a maximum flare-accelerated CME energy of $E_{flare} \lapprox 1.5 \times 10^{33}$ erg; (iv) the ratio of the summed kinetic energies of all flare-accelerated CMEs to the drag-accelerated CMEs amounts to a factor of 4; (v) the inclusion of the drag force slightly lowers the overall energy budget of CME kinetic energies in flares from $\approx 7\%$ to $\approx 4\%$; and (vi) the arrival times of CMEs at Earth can be predicted with an accuracy of $\approx 23\%$.

One of the key problems in solar flare physics is the determination of the low-energy cut-off; the value that determines the energy of nonthermal electrons and hence flare energetics. We discuss different approaches to determine the low-energy cut-off in the spectrum of accelerated electrons: (i) the total electron number model, (ii) the time-of-flight model (based on the equivalence of the time-of-flight and the collisional deflection time); (iii) the warm target model of Kontar et al.~(2015), and (iv) the model of the spectral cross-over between thermal and nonthermal components. We find that the first three models are consistent with a low-energy cutoff with a mean value of $\approx 10$ keV, while the cross-over model provides an upper limit for the low-energy cutoff with a mean value of $ \approx 21$ keV. Combining the first three models we find that the ratio of the nonthermal energy to the dissipated magnetic energy in solar flares has a mean value of $q_E=0.57\pm0.08$, which is consistent with an earlier study based on the simplified approximation of the warm target model alone ($q_E=0.51\pm0.17$). This study corroborates the self-consistency between three different low-energy cutoff models in the calculation of nonthermal flare energies.

We search for outliers in extreme events of statistical size distributions of astrophysical data sets, motivated by the {\sl Dragon-King hypothesis} of Sornette (2009), which suggests that the most extreme events in a statistical distribution may belong to a different population, and thus may be generated by a different phyiscal mechanism, in contrast to the strict power law behavior of {\sl self-organized criticality (SOC)} models. Identifying such disparate outliers is important for space weather predictions. Possible physical mechanisms to produce such outliers could be generated by sympathetic flaring. However, we find that Dragon-King events are not common in solar and stellar flares, identified in 4 out of 25 solar and stellar flare data sets only. Consequently, small, large, and extreme flares are essentially scale-free and can be modeled with a single physical mechanism. In very large data sets ($N \gapprox 10^4$) we find significant deviations from ideal power laws in almost all data sets. Neverthess, the fitted power law slopes constrain physcial scaling laws in terms of flare areas and volumes, which have the highest nonlinearity in their scaling laws.

Coronal loops in active regions are often characterized by quasi-circular and helically twisted (sigmoidal) geometries, which are consistent with dipolar potential field models in the former case, and with nonlinear force-free field models with vertical currents in the latter case. Alternatively, Parker-type nanoflare models of the solar corona hypothesize that a braiding mechanism operates between unresolved loop strands, which is a more complex topological model. In this study we use the vertical-current approximation of a nonpotential magnetic field solution (that fulfills the divergence-free and force-free conditions) to characterize the number of helical turns $N_{twist}$ in twisted coronal loops. We measure the helical twist in 15 active regions observed with AIA and HMI/SDO and find a mean nonpotentiality angle (between the potential and nonpotential field directions) of $\mu_{NP} = 15^\circ \pm 3^\circ$. The resulting mean rotational twist angle is $\varphi = 49^\circ \pm 11^\circ$, which corresponds to $N_{twist}=\varphi/360^\circ = 0.14\pm0.03$ turns with respect to the untwisted potential field, with an absolute upper limit of $N_{twist} \lapprox 0.5$, which is far below the kink instability limit of $|N_{twist}| \gapprox 1$. The number of twist turns $N_{twist}$ corresponds to the Gauss linkage number $N_{link}$ in braiding topologies. We conclude that any braided topology (with $|N_{link}| \ge 1$) cannot explain the observed stability of loops in a force-free corona, nor the observed low twist number. Parker-type nanoflaring can thus occur in non-forcefree environments only, such as in the chromosphere and transition region.

Magneto-hydrodynamic (MHD) Alfv\'en waves have been a focus of laboratory plasma physics and astrophysics for over half a century. Their unique nature makes them ideal energy transporters, and while the solar atmosphere provides preferential conditions for their existence, direct detection has proved difficult as a result of their evolving and dynamic observational signatures. The viability of Alfv\'en waves as a heating mechanism relies upon the efficient dissipation and thermalization of the wave energy, with direct evidence remaining elusive until now. Here we provide the first observational evidence of Alfv\'en waves heating chromospheric plasma in a sunspot umbra through the formation of shock fronts. The magnetic field configuration of the shock environment, alongside the tangential velocity signatures, distinguish them from conventional umbral flashes. Observed local temperature enhancements of 5% are consistent with the dissipation of mode-converted Alfv\'en waves driven by upwardly propagating magneto-acoustic oscillations, providing an unprecedented insight into the behaviour of Alfv\'en waves in the solar atmosphere and beyond.

Parker's coronal braiding and nanoflaring scenario predicts the development of tangential discontinuities and highly misaligned magnetic field lines, as a consequence of random buffeting of their footpoints due to the action of sub-photospheric convection. The increased stressing of magnetic field lines is thought to become unstable above some critical misalignment angle and to result into local magnetic reconnection events, which is generally referred to as Parker's `nanoflaring scenario'. In this study we show that the {\sl minimum (magnetic) energy principle} leads to a bifurcation of force-free field solutions for helical twist angles at $|\varphi(t)| = \pi$, which prevents the build-up of arbitrary large free energies and misalignment angles. The minimum energy principle predicts that neighbored magnetic field lines are almost parallel (with misalignment angles of $\Delta \mu \approx 1.6^\circ-1.8^\circ$), and do not reach a critical misalignment angle prone to nanoflaring. Consequently, no nanoflares are expected in the divergence-free and force-free parts of the solar corona, while they are more likely to occur in the chromosphere and transition region.

It is widely believed that loops observed in the solar atmosphere trace out magnetic field lines. However, the degree to which magnetic field extrapolations yield field lines that actually do follow loops has yet to be studied systematically. In this paper we apply three different extrapolation techniques - a simple potential model, a NLFF model based on photospheric vector data, and a NLFF model based on forward fitting magnetic sources with vertical currents - to 15 active regions that span a wide range of magnetic conditions. We use a distance metric to assess how well each of these models is able to match field lines to the 12,202 loops traced in coronal images. These distances are typically 1-2". We also compute the misalignment angle between each traced loop and the local magnetic field vector, and find values of 5-12$^\circ$. We find that the NLFF models generally outperform the potential extrapolation on these metrics, although the differences between the different extrapolations are relatively small. The methodology that we employ for this study suggests a number of ways that both the extrapolations and loop identification can be improved.

Self-organization is a property of dissipative nonlinear processes that are governed by an internal driver and a positive feedback mechanism, which creates regular geometric and/or temporal patterns and decreases the entropy, in contrast to random processes. Here we investigate for the first time a comprehensive number of 16 self-organization processes that operate in planetary physics, solar physics, stellar physics, galactic physics, and cosmology. Self-organizing systems create spontaneous {\sl order out of chaos}, during the evolution from an initially disordered system to an ordered stationary system, via quasi-periodic limit-cycle dynamics, harmonic mechanical resonances, or gyromagnetic resonances. The internal driver can be gravity, rotation, thermal pressure, or acceleration of nonthermal particles, while the positive feedback mechanism is often an instability, such as the magneto-rotational instability, the Rayleigh-B\'enard convection instability, turbulence, vortex attraction, magnetic reconnection, plasma condensation, or loss-cone instability. Physical models of astrophysical self-organization processes involve hydrodynamic, MHD, and N-body formulations of Lotka-Volterra equation systems.

The geometric arrangement of planet and moon orbits into a regularly spaced pattern of distances is the result of a self-organizing system. The positive feedback mechanism that operates a self-organizing system is accomplished by harmonic orbit resonances, leading to long-term stable planet and moon orbits in solar or stellar systems. The distance pattern of planets was originally described by the empirical Titius-Bode law, and by a generalized version with a constant geometric progression factor (corresponding to logarithmic spacing). We find that the orbital periods $T_i$ and planet distances $R_i$ from the Sun are not consistent with logarithmic spacing, but rather follow the quantized scaling $(R_{i+1}/R_i) = (T_{i+1}/T_i)^{2/3} = (H_{i+1}/H_i)^{2/3}$, where the harmonic ratios are given by five dominant resonances, namely $(H_{i+1} : H_{i}) = (3:2), (5:3), (2:1), (5:2), (3:1)$. We find that the orbital period ratios tend to follow the quantized harmonic ratios in increasing order. We apply this harmonic orbit resonance model to the planets and moons in our solar system, and to the exo-planets of 55 Cnc and HD 10180 planetary systems. The model allows us a prediction of missing planets in each planetary system, based on the quasi-regular self-organizing pattern of harmonic orbit resonance zones. We predict 7 (and 4) missing exo-planets around the star 55 Cnc (and HD 10180). The accuracy of the predicted planet and moon distances amounts to a few percents. All analyzed systems are found to have $\approx 10$ resonant zones that can be occupied with planets (or moons) in long-term stable orbits.

Magnetic elements of the solar surface are studied in magnetograms recorded with the high-resolution Solar Dynamics Observatory / Helioseismic and Magnetic Imager . To extract some statistical and physical properties of these elements (e.g., filling factors, magnetic flux, size, lifetimes), the Yet Another Feature Tracking Algorithm (YAFTA), a region-based method, is employed. An area with 400$^{\prime\prime}\times$400$^{\prime\prime}$ was selected to investigate the magnetic characteristics during the year 2011. The correlation coefficient between filling factors of negative and positive polarities is 0.51. A broken power law fit was applied to the frequency distribution of size and flux. Exponents of the power-law distributions for sizes smaller and greater than 16 arcsec$^2$ were found to be -2.24 and -4.04, respectively. The exponents of power$-$law distributions for fluxes smaller and greater than 2.63$\times$10$^{19}$ Mx were found to be -2.11 and -2.51, respectively. The relationship between the size ($S$) and flux ($F$) of elements can be expressed by a power-law behavior in the form of $S\propto F~^{0.69}$. The lifetime and its relationship with the flux and size of quiet-Sun (QS) elements are studied during three days. The code detected patches with lifetimes of about 15 hours, which we call long-duration events. It is found that more than 95\% of the magnetic elements have lifetimes of less than 100 minutes. About 0.05\% of the elements were found with lifetimes of more than 6 hours. The relationships between the size (S), lifetime (T), and the flux (F) for patches in the QS, indicate the power$-$law relationships $S\propto T~^{0.25}$ and $F\propto T~^{0.38}$, respectively. Executing a detrended fluctuation analysis of the time series of new emerged magnetic elements, we find a Hurst exponent of 0.82, which implies long-range temporal correlation in the system.

The current exoplanet database includes 5454 confirmed planets and candidate planets observed with the KEPLER mission. We find 932 planet pairs from which we extract distance and orbital period ratios. While earlier studies used the Titius-Bode law or a generalized version with logarithmic spacing, which both lack a physical model, we employ here the theory of harmonic orbit resonances, which contains quantized ratios instead, to explain the observed planet distance ratios and to predict undetected exoplanets. We find that the most prevailing harmonic ratios are (2:1), (3:2), and (5:3), in 73\% of the cases, while alternative harmonic ratios of (5:4), (4:3), (5:2), (3:1) occur in 27\% of the other cases. Our orbital predictions includes 171 exoplanets, 2 Jupiter moons, one Saturn moon, 3 Uranus moons, and 4 Neptune moons. The accuracy of the predicted planet distances amounts to a few percent, which fits the data significantly better than the Titius-Bode law or a logarithmic spacing. This information may be useful for targeted exoplanet searches with Kepler data and to estimate the number of live-carrying planets in habitable zones.

In this study we refine a CME model presented in an earlier study on the global energetics of solar flares and associated CMEs, and apply it to all (860) GOES M- and X-class flare events observed during the first 7 years (2010-2016) of the Solar Dynamics Observatory (SDO) mission, which doubles the statistics of the earlier study. The model refinements include: (1) the CME geometry in terms of a 3D sphere undergoing self-similar adiabatic expansion; (2) the inclusion of solar gravitational deceleration during the acceleration and propagation of the CME, which discriminates eruptive and confined CMEs; (4) a self-consistent relationship between the CME center-of-mass motion detected during EUV dimming and the leading-edge motion observed in white-light coronagraphs; (5) the equi-partition of the CME kinetic and thermal energy; and (6) the Rosner-Tucker-Vaiana (RTV) scaling law. The refined CME model is entirely based on EUV dimming observations (using AIA/SDO data) and complements the traditional white-light scattering model (using LASCO/SOHO data), and both models are independently capable to determine fundamental CME parameters such as the CME mass, speed, and energy. Comparing the two methods we find that: (1) LASCO is less sensitive than AIA in detecting CMEs (in 24$\%$ of the cases); (2) CME masses below $m_{cme} \sim 10^{14}$ g are under-estimated by LASCO; (3) AIA and LASCO masses, speeds, and energy agree closely in the statistical mean after elimination of outliers; (4) the CMEs parameters of the speed $v$, emission measure-weighted flare peak temperature $T_e$, and length scale $L$ are consistent with the following scaling laws (derived from first principles): $v \propto T_e^{1/2}$, $v \propto (m_{cme})^{1/4}$, and $m_{cme} \propto L^2$.

In this study we analyze {\sl Atmospheric Imaging Assembly (AIA)} and Hi-C images in order to investigate absolute limits for the finest loop strands. We develop a model of the occurrence-size distribution function of coronal loop widths, characterized by a lower limit of widths $w_{min}$, a peak width $w_p$, a peak occurrence number $n_p$, and a power law slope $a$. Our data analysis includes automated tracing of curvi-linear features with the OCCULT-2 code, automated sampling of the cross-sectional widths of coronal loops, and fitting of the theoretical size distribution to the observed distribution. With Monte-Carlo simulations and variable pixel sizes $\Delta x$ we derive a first diagnostic criterion to discriminate whether the loop widths are unresolved $(w_p/\Delta x \approx 2.5\pm0.2)$, or fully resolved (if $w_p/\Delta x > 2.7$). For images with resolved loop widths we can apply a second diagnostic criterion that predicts the lower limit of loop widths, $w_{min} \approx 3 (\Delta x_{crit}-0.37")$ as a function of the critical resolution $\Delta x_{crit}$. We find that the loop widths are marginally resolved in AIA images, but are fully resolved in Hi-C images, where our model predicts a lower limit of loop widths at $w_{min} \approx 100$ km and a most frequent (peak) value at $w_p \approx 300$ km, in agreement with recent results of Brooks et al. This result agrees with the statistics of photospheric granulation sizes and thus supports coronal heating mechanisms operating on the macroscopic scale of photospheric magneto-convection, rather than nanoflare heating models with unresolved microscopic scales.

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles ($E_{\mathrm{nt}}$), the energy of direct heating ($E_{\mathrm{dir}}$), and the energy in coronal mass ejections ($E_{\mathrm{CME}}$), which are the primary energy dissipation processes in a flare, is found to have a ratio of $(E_{\mathrm{nt}}+E_{\mathrm{dir}}+ E_{\mathrm{CME}})/E_{\mathrm{mag}} = 0.87 \pm 0.18$, compared with the dissipated magnetic free energy $E_{\mathrm{mag}}$, which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: $0.51\pm0.17$ in nonthermal energy of $\ge 6$ keV electrons, $0.17\pm0.17$ in nonthermal $\ge 1$ MeV ions, $0.07\pm0.14$ in CMEs, and $0.07\pm0.17$ in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction $\approx 0.03$ of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at $e_c \approx 6$ keV, based on the mean differential emission measure (DEM) peak temperature of $T_e=8.6$ MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

This study entails the third part of a global flare energetics project, in which Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) data of 191 M and X-class flare events from the first 3.5 yrs of the Solar Dynamics Observatory (SDO) mission are analyzed. We fit a thermal and a nonthermal component to RHESSI spectra, yielding the temperature of the differential emission measure (DEM) tail, the nonthermal power law slope and flux, and the thermal/nonthermal cross-over energy $e_{\mathrm{co}}$. From these parameters we calculate the total nonthermal energy $E_{\mathrm{nt}}$ in electrons with two different methods: (i) using the observed cross-over energy $e_{\mathrm{co}}$ as low-energy cutoff, and (ii) using the low-energy cutoff $e_{\mathrm{wt}}$ predicted by the warm thick-target bremsstrahlung model of Kontar et al. {\bf Based on a mean temperature of $T_e=8.6$ MK in active regions we find low-energy cutoff energies of $e_{\mathrm{wt}} =6.2\pm 1.6$ keV for the warm-target model, which is significantly lower than the cross-over energies $e_{\mathrm{co}}=21 \pm 6$ keV. Comparing with the statistics of magnetically dissipated energies $E_{\mathrm{mag}}$ and thermal energies $E_{\mathrm{th}}$ from the two previous studies, we find the following mean (logarithmic) energy ratios with the warm-target model: $E_{\mathrm{nt}} = 0.41 \ E_{\mathrm{mag}}$, $E_{\mathrm{th}} = 0.08 \ E_{\mathrm{mag}}$, and $E_{\mathrm{th}} = 0.15 \ E_{\mathrm{nt}}$. The total dissipated magnetic energy exceeds the thermal energy in 95% and the nonthermal energy in 71% of the flare events, which confirms that magnetic reconnection processes are sufficient to explain flare energies. The nonthermal energy exceeds the thermal energy in 85\% of the events, which largely confirms the warm thick-target model.

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.

The aim of this study is to explore the suitability of chromospheric images for magnetic modeling of active regions. We use high-resolution images (0.1") from the Interferometric Bidimensional Spectrometer (IBIS) in the Ca II 8542 A line, the Rapid Oscillations in the Solar Atmosphere (ROSA) instrument in the H-alpha 6563 A line, the Interface Region Imaging Spectrograph (IRIS) in the 2796 A line, and compare non-potential magnetic field models obtained from those chromospheric images with those obtained from images of the Atmospheric Imaging Assembly (AIA) in coronal (171 A, etc.) and in chromospheric (304 A) wavelengths. Curvi-linear structures are automatically traced in those images with the OCCULT-2 code, to which we forward-fitted magnetic field lines computed with the Vertical-Current Approximation Non-Linear Force Free Field (VCA-NLFFF) code. We find that the chromospheric images: (1) reveal crisp curvi-linear structures (fibrils, loop segments, spicules) that are extremely well-suited for constraining magnetic modeling; (2) that these curvi-linear structures are field-aligned with the best-fit solution by a median misalignment angle of ~4-7 deg; (3) the free energy computed from coronal data may underestimate that obtained from cromospheric data by a factor of ~ 2-4, (4) the height range of chromospheric features is confined to h ~ 4000$ km, while coronal features are detected up to h ~ 35,000$ km; and (5) the plasma-beta parameter is beta ~ 10^(-5)-10^(-1) for all traced features. We conclude that chromospheric images reveal important magnetic structures that are complementary to coronal images and need to be included in comprehensive magnetic field models, a quest that is not accomodated in standard NLFFF codes.

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

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.

We compare the ability of 11 Differential Emission Measure (DEM) forward-fitting and inversion methods to constrain the properties of active regions and solar flares by simulating synthetic data using the instrumental response functions of SDO/AIA, SDO/EVE, RHESSI, and GOES/XRS. The codes include the single-Gaussian DEM, a bi-Gaussian DEM, a fixed-Gaussian DEM, a linear spline DEM, the spatial synthesis DEM, the Monte-Carlo Markov chain DEM, the regularized DEM inversion, the Hinode/XRT method, a polynomial spline DEM, an EVE+GOES, and an EVE+RHESSI method. Averaging the results from all 11 DEM methods, we find the following accuracies in the inversion of physical parameters: the EM-weighted temperature $T_w^{fit}/T_w^{sim}=0.9\pm0.1$, the peak emission measure $EM_p^{fit}/EM_p^{sim}=0.6\pm0.2$, the total emission measure $EM_t^{fit}/EM_t^{sim}=0.8\pm0.3$, and the multi-thermal energies $E_{th}^{fit}/EM_{th}^{sim}=1.2\pm0.4$. We find that the AIA spatial synthesis, the EVE+GOES, and the EVE+RHESSI method yield the most accurate results.

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 $\alpha_s \approx 5^\circ$.

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

We present the second part of a project on the global energetics of solar flares and CMEs that includes about 400 M- and X-class flares observed with AIA/SDO during the first 3.5 years of its mission. In this Paper II we compute the differential emission measure (DEM) distribution functions and associated multi-thermal energies, using a spatially-synthesized Gaussian DEM forward-fitting method. The multi-thermal DEM function yields a significantly higher (by an average factor of $\approx 14$), but more comprehensive (multi-)thermal energy than an isothermal energy estimate from the same AIA data. We find a statistical energy ratio of $E_{th}/E_{diss} \approx 2\%-40\%$ between the multi-thermal energy $E_{th}$ and the magnetically dissipated energy $E_{diss}$, which is an order of magnitude higher than the estimates of Emslie et al.~2012. For the analyzed set of M and X-class flares we find the following physical parameter ranges: $L=10^{8.2}-10^{9.7}$ cm for the length scale of the flare areas, $T_p=10^{5.7}-10^{7.4}$ K for the DEM peak temperature, $T_w=10^{6.8}-10^{7.6}$ K for the emission measure-weighted temperature, $n_p=10^{10.3}-10^{11.8}$ cm$^{-3}$ for the average electron density, $EM_p=10^{47.3}-10^{50.3}$ cm$^{-3}$ for the DEM peak emission measure, and $E_{th}=10^{26.8}-10^{32.0}$ erg for the multi-thermal energies. The deduced multi-thermal energies are consistent with the RTV scaling law $E_{th,RTV} = 7.3 \times 10^{-10} \ T_p^3 L_p^2$, which predicts extremal values of $E_{th,max} \approx 1.5 \times 10^{33}$ erg for the largest flare and $E_{th,min} \approx 1 \times 10^{24}$ erg for the smallest coronal nanoflare. The size distributions of the spatial parameters exhibit powerlaw tails that are consistent with the predictions of the fractal-diffusive self-organized criticality model combined with the RTV scaling law.

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_{\varphi}^{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.

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.

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

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.

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

We study the temporal, spatial and spectral evolution of the M1.8 flare, which occurred in NOAA AR 11195 (S17E31) on 22 April 2011, and explore the underlying physical processes during the precursors and their relation to the main phase. The study of the source morphology using the composite images in 131 {\deg}A wavelength observed by the SDO/AIA and 6-14 keV revealed a multiloop system that destabilized systematically during the precursor and main phases. In contrast, HXR emission (20-50 keV) was absent during the precursor phase, appearing only from the onset of the impulsive phase in the form of foot-points of emitting loop/s. This study has also revealed the heated loop-top prior to the loop emission, although no accompanying foot-point sources were observed during the precursor phase. We estimate the flare plasma parameters viz. T, EM, power-law index, and photon turn-over energy by forward fitting RHESSI spectral observations. The energy released in the precursor phase was thermal and constituted ~1 per cent of the total energy released during the flare. The study of morphological evolution of the filament in conjunction with synthesized T and EM maps has been carried out which reveals (a) Partial filament eruption prior to the onset of the precursor emission, (b) Heated dense plasma over the polarity inversion line and in the vicinity of the slowly rising filament during the precursor phase. Based on the implications from multi-wavelength observations, we propose a scheme to unify the energy release during the precursor and main phase emissions in which, the precursor phase emission has been originated via conduction front formed due to the partial filament eruption. Next, the heated leftover S-shaped filament has undergone slow rise and heating due to magnetic reconnection and finally erupted to produce emission during the impulsive and gradual phases.

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(\Phi) \propto \Phi^{-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.

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 $\kappa$, spreading exponent $\beta$, 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 $\kappa \propto L^{0.94\pm0.01}$ and the 3-parameter scaling law $L \propto \kappa\ 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 $\kappa$ than confined flares (without CMEs).

We study the association of solar flares with coronal mass ejections (CMEs) during the deep, extended solar minimum of 2007-2009, using extreme-ultraviolet (EUV) and white-light (coronagraph) images from the {\it Solar Terrestrial Relations Observatory} (STEREO). Although all of the fast (v $>$ 900 km s$^{-1}$) {\it and} wide ($\theta >$ 100$\arcdeg$) CMEs are associated with a flare that is at least identified in GOES soft X-ray light curves, a majority of flares with relatively high X-ray intensity for the deep solar minimum (e.g. $\gtrsim$1 \times 10^{-6}$ W m$^{-2}$ or C1) are not associated with CMEs. Intense flares tend to occur in active regions with strong and complex photospheric magnetic field, but the active regions that produce CME-associated flares tend to be small, including those that have no sunspots and therefore no NOAA active-region numbers. Other factors on scales comparable to and larger than active regions seem to exist that contribute to the association of flares with CMEs. We find the possible low coronal signatures of CMEs, namely eruptions, dimmings, EUV waves, and Type III bursts, in 91%, 74%, 57%, and 74%, respectively, of the 35 flares that we associate with CMEs. None of these observables can fully replace direct observations of CMEs by coronagraphs.

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

With increasing solar activity since 2010, many flares from the backside of the Sun have been observed by the Extreme Ultraviolet Imager (EUVI) on either of the twin STEREO spacecraft. Our objective is to estimate their X-ray peak fluxes from EUVI data by finding a relation of the EUVI with GOES X-ray fluxes. Because of the presence of the Fe xxiv line at 192 A, the response of the EUVI 195 A channel has a secondary broad peak around 15 MK, and its fluxes closely trace X-ray fluxes during the rise phase of flares. If the flare plasma is isothermal, the EUVI flux should be directly proportional to the GOES flux. In reality, the multithermal nature of the flare and other factors complicate the estimation of the X-ray fluxes from EUVI observations. We discuss the uncer- tainties, by comparing GOES fluxes with the high cadence EUV data from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). We conclude that the EUVI 195 A data can provide estimates of the X-ray peak fluxes of intense flares (e.g., above M4 in the GOES scale) with uncertainties of a factor of a few. Lastly we show examples of intense flares from regions far behind the limb, some of which show eruptive signatures in AIA images.

We use multiwavelength data from space and ground based instruments to study the solar flares and coronal mass ejections (CMEs) on January 23, 2012 that were responsible for one of the largest solar energetic particle (SEP) events of solar cycle 24. The eruptions consisting of two fast CMEs (1400 km/s and 2000 km/s) and M-class flares that occurred in active region 11402 located at N28 W36. The two CMEs occurred in quick successions, so they interacted very close to the Sun. The second CME caught up with the first one at a distance of 11-12 Rsun. The CME interaction may be responsible for the elevated SEP flux and significant changes in the intensity profile of the SEP event. The compound CME resulted in a double-dip moderate geomagnetic storm (Dst = -73 nT). The two dips are due to the southward component of the interplanetary magnetic field in the shock sheath and the ICME intervals. One possible reason for the lack of a stronger geomagnetic storm may be that the ICME delivered a glancing blow to Earth.

Employing coronagraphic and EUV observations close to the solar surface made by the Solar Terrestrial Relations Observatory (STEREO) mission, we determined the heliocentric distance of coronal mass ejections (CMEs) at the starting time of associated metric type II bursts. We used the wave diameter and leading edge methods and measured the CME heights for a set of 32 metric type II bursts from solar cycle 24. We minimized the projection effects by making the measurements from a view that is roughly orthogonal to the direction of the ejection. We also chose image frames close to the onset times of the type II bursts, so no extrapolation was necessary. We found that the CMEs were located in the heliocentric distance range from 1.20 to 1.93 solar radii (Rs), with mean and median values of 1.43 and 1.38 Rs, respectively. We conclusively find that the shock formation can occur at heights substantially below 1.5 Rs. In a few cases, the CME height at type II onset was close to 2 Rs. In these cases, the starting frequency of the type II bursts was very low, in the range 25 to 40 MHz, which confirms that the shock can also form at larger heights. The starting frequencies of metric type II bursts have a weak correlation with the measured CME/shock heights and are consistent with the rapid decline of density with height in the inner corona.

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

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 $\Delta\mu=\mu_{\rm P}-\mu_{\rm N}$ between the potential and nonpotential field, which scales as $E_{\rm free}/E_{\rm P} \approx \tan^2{(\Delta\mu)}$. 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.

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 \kappa (t-t_1)^{\beta/2}$, which includes the logistic growth limit ($\beta=0$), sub-diffusion ($\beta = 0-1$), classical diffusion ($\beta=1$), super-diffusion ($\beta = 1-2$), and the linear expansion limit ($\beta=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 ($\beta=0.53\pm0.27$), which we interpret in terms of anisotropic chain reactions of intermittent magnetic reconnection episodes in a low plasma-$\beta$ 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.

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

We performed for the first time stereoscopic triangulation of coronal loops in active regions over the entire range of spacecraft separation angles ($\alpha_{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 $\alpha_{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 $\alpha_{sep}\approx 6^\circ-127^{\circ}$ we find a mean 3D misalignment angle of $\mu_{PF} \approx 21^\circ-39^\circ$ of stereoscopically triangulated loops with magnetic potential field models, and $\mu_{FFF} \approx 15^\circ-21^\circ$ for a force-free field model, which is partly caused by stereoscopic uncertainties $\mu_{SE} \approx 9^\circ$. We predict optimum conditions for solar stereoscopy during the time intervals of 2012--2014, 2016--2017, and 2021--2023.

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 $\mu=19^\circ-46^\circ$ for a potential field model, which is reduced to $\mu=14^\circ-19^\circ$ for a non-potential field model parameterized by twist parameters. The residual error is commensurable with stereoscopic measurement errors ($\mu_{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 $\alpha$-parameter (median $\alpha \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 $\alpha$), represents the first code that enables fast forward-fitting to photospheric magnetograms and stereoscopically triangulated loops in the solar corona.

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 $\alpha({\bf x})$, the electric current density ${\bf j}({\bf x})$, and are accurate to second-order (of the nonlinear force-free $\alpha$-parameter). The explicit expressions of a force-free field can easily be applied to modeling or forward-fitting of many coronal phenomena.

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 $\mu \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 $\alpha$-parameters ($\alpha_{\rm fit}/\alpha_{\rm model} \approx 0.9-1.0$ for simulations and $\alpha_{\rm fit}/\alpha_{\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.

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 $\alpha_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 $\alpha_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 $\alpha_T =2.02\pm0.04$ during solar cycle minima years, which is also consistent with the prediction $\alpha_T=2.0$ of the FD-SOC model. During solar cycle maxima years, the powerlaw slope is steeper in the range of $\alpha_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 $\alpha_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.

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_{\Phi,rr} \approx 1$ MeV, maximum electron energies of $U_{e,rr} \approx 1$ PeV, and luminosities of $L \approx 10^{36}$ erg s$^{-1}$.

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.

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.

{\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.

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^{-\alpha_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 $\alpha_T=(1+S)/2$, $\alpha_F=1+(S-1)/D_S$, $\alpha_P=2-1/S$, and $\alpha_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 $\alpha_E=1.40-1.67$, and the peak energy distribution has a slope of $\alpha_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.

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.

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

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\'enic 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.

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(\epsilon) from the flare by integrating a series of isothermal plasma flux. We find that multi-temperature integrated flux F(\epsilon) is a power-law function of \epsilon with a spectral index (\gamma) \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 (\gamma) varies in the range of -4.4 and -5.7. The temporal evolution of the X-ray flux F(\epsilon,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.

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 $\sigma_{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.

We present the science case for a broadband X-ray imager with high-resolution spectroscopy, including simulations of X-ray spectral diagnostics of both active regions and solar flares. This is part of a trilogy of white papers discussing science, instrument (Bandler et al. 2010), and missions (Bookbinder et al. 2010) to exploit major advances recently made in transition-edge sensor (TES) detector technology that enable resolution better than 2 eV in an array that can handle high count rates. Combined with a modest X-ray mirror, this instrument would combine arcsecondscale imaging with high-resolution spectra over a field of view sufficiently large for the study of active regions and flares, enabling a wide range of studies such as the detection of microheating in active regions, ion-resolved velocity flows, and the presence of non-thermal electrons in hot plasmas. It would also enable more direct comparisons between solar and stellar soft X-ray spectra, a waveband in which (unusually) we currently have much better stellar data than we do of the Sun.

The coming decade will see the routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness. To capitalize on the new (or soon to be available) facilities such as SDO, ATST and FASR, and the challenges they present in the visualization and synthesis of multi-wavelength datasets, we propose that realistic, sophisticated, 3D active region and flare modeling is timely and critical, and will be a forefront of coronal studies over the coming decade. To make such modeling a reality, a broad, concerted effort is needed to capture the wealth of information resulting from the data, develop a synergistic modeling effort, and generate the necessary visualization, interpretation and model-data comparison tools to accurately extract the key physics.

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.

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 $\alpha_P=1.72\pm0.08$ for the peak flux, $\alpha_E=1.60\pm0.14$ for the total flux, and $\alpha_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 $\alpha_P(t)=\alpha_0+\Delta \alpha \cos{[2\pi (t-t_0)/T_{cycle}]}$, with a mean of $\alpha_0=1.73$, a variation of $\Delta \alpha =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.

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^{-\alpha_E}$ with a powerlaw slope of $\alpha_E \approx 1.5$. The powerlaw slope $\alpha_E$ of energies can be related to the fractal dimension $D$ of the spatial energy dissipation domain by $D=3/\alpha_E$, which predicts a powerlaw slope $\alpha_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.

We investigate the recently quantified misalignment of $\alpha_{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 $\alpha_{PF} \approx 13^\circ-20^\circ$. We evaluate also stereoscopic measurement errors and find a contribution of $\alpha_{SE}\approx 7^\circ-12^\circ$, which constrains the residual misalignment to $\alpha_{NP}=\alpha_{PF}-\alpha_{SE}\approx 5^\circ -9^\circ$, which is likely due to the nonpotentiality of the active regions. The residual misalignment angle $\alpha_{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.

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.

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 ($\Delta t \approx 10^{-3}- 10^3$ hrs), can be reconciled with a single distribution function, $N(\Delta t) \propto \lambda_0 (1 + \lambda_0 \Delta t)^{-2}$, which has a powerlaw slope of $p \approx 2.0$ at large waiting times ($\Delta t \approx 1-1000$ hrs) and flattens out at short waiting times $\Delta t \lapprox \Delta t_0 = 1/\lambda_0$. We find a consistent breakpoint at $\Delta t_0 = 1/\lambda_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, $\Delta t_0 \propto 1/n_{det}$. This waiting time distribution can be modeled with a nonstationary Poisson process with a flare rate $\lambda=1/\Delta t$ that varies as $f(\lambda) \propto \lambda^{-1} \exp{-(\lambda/\lambda_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.

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.

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.

Our understanding of fundamental processes in the solar corona has been greatly progressed based on the space observations of SMM, Yohkoh, Compton GRO, SOHO, TRACE, RHESSI, and STEREO. We observe now acoustic waves, MHD oscillations, turbulence-related line broadening, magnetic configurations related to reconnection processes, and radiation from high-energy particles on a routine basis. We review a number of key observations in EUV, soft X-rays, and hard X-rays that innovated our physical understanding of the solar corona, in terms of hydrodynamics, MHD, plasma heating, and particle acceleration processes.

In this study we compile for the first time comprehensive data sets of solar and stellar flare parameters, including flare peak temperatures T_p, flare peak volume emission measures EM_p, and flare durations t_f from both solar and stellar data, as well as flare length scales L from solar data. Key results are that both the solar and stellar data are consistent with a common scaling law of EM_p ~ T_p^4.7, but the stellar flares exhibit ~250 times higher emission measures (at the same flare peak temperature). For solar flares we observe also systematic trends for the flare length scale L(T_p) ~ T_p^0.9 and the flare duration t_F(T_p) ~ T_p^0.9 as a function of the flare peak temperature. Using the theoretical RTV scaling law and the fractal volume scaling observed for solar flares, i.e., V(L) ~ L^2.4, we predict a scaling law of EM_p ~ T_p^4.3, which is consistent with observations, and a scaling law for electron densities in flare loops, n_p ~ T_p^2/L ~ T_p^1.1. The RTV-predicted electron densities were also found to be consistent with densities inferred from total emission measures, n_p=(EM_p/q_V*V)^1/2, using volume filling factors of q_V=0.03-0.08 constrained by fractal dimensions measured in solar flares. Our results affect also the determination of radiative and conductive cooling times, thermal energies, and frequency distributions of solar and stellar flare energies.

The fastest pulsar and the slowest nova; the oldest galaxies and the youngest stars; the weirdest life forms and the commonest dwarfs; the highest energy particles and the lowest energy photons. These were some of the extremes of Astrophysics 2006. We attempt also to bring you updates on things of which there is currently only one (habitable planets, the Sun, and the universe) and others of which there are always many, like meteors and molecules, black holes and binaries.

We bring you, as usual, the Sun and Moon and stars, plus some galaxies and a new section on astrobiology. Some highlights are short (the newly identified class of gamma-ray bursts, and the Deep Impact on Comet 9P/ Tempel 1), some long (the age of the universe, which will be found to have the Earth at its center), and a few metonymic, for instance the term "down-sizing" to describe the evolution of star formation rates with redshift.

In the recent work of Nakariakov et al. (2004), it has been shown that the time dependences of density and velocity in a flaring loop contain pronounced quasi-harmonic oscillations associated with the 2nd harmonic of a standing slow magnetoacoustic wave. That model used a symmetric heating function (heat deposition was strictly at the apex). This left outstanding questions: A) is the generation of the 2nd harmonic a consequence of the fact that the heating function was symmetric? B) Would the generation of these oscillations occur if we break symmetry? C) What is the spectrum of these oscillations? Is it consistent with a 2nd spatial harmonic? The present work (and partly Tsiklauri et al. (2004b)) attempts to answer these important outstanding questions. Namely, we investigate the physical nature of these oscillations in greater detail: we study their spectrum (using periodogram technique) and how heat positioning affects the mode excitation. We found that excitation of such oscillations is practically independent of location of the heat deposition in the loop. Because of the change of the background temperature and density, the phase shift between the density and velocity perturbations is not exactly a quarter of the period, it varies along the loop and is time dependent, especially in the case of one footpoint (asymmetric) heating. We also were able to model successfully SUMER oscillations observed in hot coronal loops.

Low-frequency longitudinal oscillations of a flaring coronal loop are studied numerically. In the recent work of Nakariakov et al., Astron. and Astrophys., 414, L25-L28 (2004) it has been shown that the time dependences of density and velocity in a flaring loop contain well-pronounced quasi-harmonic oscillations associated with a 2nd harmonic of a standing slow magnetoacoustic wave. In this work we investigate physical nature of these oscillations in greater detail, namely, we study their spectrum (using periodogram technique) and how does heat positioning affects the mode excitation. We found that excitation of such oscillations practically independent of the positioning of the heat deposition in the loop. Because of the change of the background temperature and density, the phase shift between the density and velocity perturbations is not exactly equal to the quarter of the period, it varies along the loop and is time dependent, especially in the case of one footpoint (asymmetric) heating. KEYWORDS -- Sun: flares -- Sun: oscillations -- Sun: Corona -- Stars: flare -- Stars: oscillations -- Stars: coronae

A 1-D loop radiative hydrodynamic model that incorporates the effects of gravitational stratification, heat conduction, radiative losses, external heat input, presence of helium, and Braginskii viscosity is used to simulate elementary flare loops. The physical parameters for the input are taken from observations of the Bastille-Day flare of 2000 July 14. The present analysis shows that: (a) The obtained maximum values of the electron density can be considerably higher (4.2 X 10^{11} cm^{-3} or more) in the case of footpoint heating than in the case of apex heating (2.5 X 10^{11} cm^{-3}). (b) The average cooling time after the flare peak takes less time in the case of footpoint heating than in the case of apex heating. (c) The peak apex temperatures are significantly lower (by about 10 MK) for the case of footpoint heating than for apex heating (for the same average loop temperature of about 30 MK). This characteristic allows us to discriminate between different heating positioning. (d) In both cases (of apex and footpoint heating), the maximum obtained apex temperature T^{max} is practically independent of the heating duration sigma_{t}, but scales directly with the heating rate E_{H0}. (e) The maximum obtained densities at the loop apex, n_e^{max}, increase with the heating rate E_{H0} and heating duration sigma_{t} for both footpoint and apex heating. In Paper II we will use the outputs of these hydrodynamic simulations, which cover a wide range of the parameter space of heating rates and durations, as an input for forward-fitting of the multi-loop arcade of the Bastille-day flare. KEYWORDS: Sun: Flares -- Sun: Activity -- Sun: Corona

Evolution of a coronal loop in response to an impulsive energy release is numerically modelled. It is shown that the loop density evolution curves exhibit quasi-periodic perturbations with the periods given approximately by the ratio of the loop length to the average sound speed, associated with the second standing harmonics of an acoustic wave. The density perturbations have a maximum near the loop apex. The corresponding field-aligned flows have a node near the apex. We suggest that the quasi-periodic pulsations with periods in the range 10-300 s, frequently observed in flaring coronal loops in the radio, visible light and X-ray bands, may be produced by the second standing harmonic of the acoustic mode.

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.

Recent observations with TRACE reveal that the time delay between the appearance of a cooling loop in different EUV temperature filters is proportional to the loop length, dt_12 ~ L. We model this cooling delay in terms of radiative loss and confirm this linear relationship theoretically. We derive an expression that can be used to constrain the coronal iron enhancement alpha_Fe=A_Fe^cor/A_Fe^Ph relative to the photospheric value as function of the cooling delay dt_12, flux F_2, loop width w, and filling factor q_w < 1. With this relation we find upper limits on the iron abundance enhancement of alpha_Fe < 4.8+/-1.7 for 10 small-scale nanoflare loops, and alpha_Fe < 1.4+/-0.4 for 5 large-scale loops, in the temperature range of T~1.0-1.4 MK. This result supports the previous finding that low-FIP elements, including Fe, are enhanced in the corona. The same relation constitutes also a lower limit for the filling factor, which is q_w > 0.2+/-0.1 and q_w > 0.8+/-0.2 for the two groups of coronal loops.

Recent observations show a variety of oscillation modes in the corona. Early non-imaging observations in radio wavelengths showed a number of fast-period oscillations in the order of seconds, which have been interpreted as fast sausage mode oscillations. TRACE observations from 1998 have for the first time revealed the lateral displacements of fast kink mode oscillations, with periods of ~3-5 minutes, apparently triggered by nearby flares and destabilizing filaments. Recently, SUMER discovered with Doppler shift measurements loop oscillations with longer periods (10-30 minutes) and relatively short damping times in hot (7 MK) loops, which seem to correspond to longitudinal slow magnetoacoustic waves. In addition, propagating longitudinal waves have also been detected with EIT and TRACE in the lowest density scale height of loops near sunspots. All these new observations seem to confirm the theoretically predicted oscillation modes and can now be used as a powerful tool for ``coronal seismology'' diagnostic.

We compare the photometric accuracy of spectra and images in flares observed with the Ramaty High Energy Solar Spectroscopic Imager (RHESSI)}spacecraft. We test the accuracy of the photometry by comparing the photon fluxes obtained in different energy ranges from the spectral-fitting software SPEX with those fluxes contained in the images reconstructed with the Clean, MEM, MEM-Vis, Pixon, and Forward-fit algorithms. We quantify also the background fluxes, the fidelity of source geometries, and spatial spectra reconstructed with the five image reconstruction algorithms. We investigate the effects of grid selection, pixel size, field-of-view, and time intervals on the quality of image reconstruction. The detailed parameters and statistics are provided in an accompanying CD-ROM and web page. We find that Forward-fit, Pixon, and Clean have a robust convergence behavior and a photometric accuracy in the order of a few percents, while MEM does not converge optimally for large degrees of freedom (for large field-of-views and/or small pixel sizes), and MEM-Vis suffers in the case of time-variable sources. This comparative study documents the current status of the RHESSI spectral and imaging software, one year after launch.

A recent study by Nakariakov et al. pointed out that the dispersion relation of MHD sausage mode oscillations has been incorrectly applied to coronal loops, neglecting the highly dispersive nature of the phase speed and the long-wavelength cutoff of the wave number. In the light of these new insights we revisit previous observations that have been interpreted in terms of MHD sausage mode oscillations in coronal loops and come to the following conclusions: (1) Fast sausage MHD mode oscillations require such a high electron density imposed by the wave number cutoff that they can only occur in flare loops; (2) In the previously reported radio observations ($\nu \approx 100$ MHz to 1 GHz) with periods of $P\approx 0.5-5$ s, the fast sausage MHD mode oscillation is likely to be confined to a small segment (corresponding to a high harmonic node) near the apex of the loop, rather than involving a global oscillation over the entire loop length. The recent microwave and soft X-ray observations of fast periods ($P\approx 6-17$ s) by Asai et al. and Melnikov et al., however, are consistent with fast sausage MHD oscillations at the fundamental harmonic.

One of the proposed damping mechanisms of coronal (transverse) loop oscillations in the kink-mode is resonant absorption as a result of the Alfven speed variation at the outer boundary of coronal loops. Analytical expressions for the period and damping time exist for loop models with thin non-uniform boundaries. Here we measure the thickness of the non-uniform layer in oscillating loops for 11 events, by forward-fitting of the cross-sectional density profile and line-of-sight integration to the cross-sectional fluxes observed with TRACE 171 A. This way we model the internal and external electron density of the coronal plasma in oscillating loops. This allows us to test the theoretically predicted damping rates for thin boundaries as function of the density ratio. We find that the density ratio predicted by the damping time is higher than the density ratio estimated from the background fluxes. The lower densities modeled from the background fluxes are likely to be a consequence of the neglected hotter plasma that is not detected with the TRACE 171 A filter. Taking these correction into account, resonant absorption predicts damping times of kink-mode oscillations that are commensurable with the observed ones and provides a new diagnostic of the density contrast of oscillating loops.

An analytical MHD model of coronal loops with compressible flows and including heating is compared to observational data. The model is constructed via a systematic nonlinear separation of the variables method used to calculate several classes of exact MHD equilibria in Cartesian geometry and uniform gravity. By choosing a particularly versatile solution class with a large parameter space we are able to calculate models whose loop length, shape, plasma density, temperature and velocity profiles are fitted to loops observed with TRACE, SoHO/CDS and SoHO/SUMER. Synthetic emission profiles are also calculated and fitted to the observed emission patterns. An analytical discussion is given of the two-dimenional balance of the Lorentz force and the gas pressure gradient, gravity and inertial forces acting along and across the loop. These models are the first to include a fully consistent description of the magnetic field, 2D geometry, plasma density and temperature, flow velocity and thermodynamics of loops. The consistently calculated heating profiles which are largely dominated by radiative losses are influenced by the flow and are asymmetrical being concentrated at the upflow footpoint.

We investigate the effect of hydrostatic scale heights ${\lambda}(T)$ in coronal loops on the determination of the vertical temperature structure $T(h)$ of the solar corona. Every method that determines an average temperature at a particular line-of-sight from optically thin emission (e.g. in EUV or soft X-ray wavelengths) of a mutli-temperature plasma, is subject to the emission measure-weighted contributions $dEM(T)/dT$ from different temperatures. Because most of the coronal structures (along open or closed field lines) are close to hydrostatic equilibrium, the hydrostatic temperature scale height introduces a height-dependent weighting function that causes a systematic bias in the determination of the temperature structure $T(h)$ as function of altitude $h$. The net effect is that the averaged temperature seems to increase with altitude, $dT(h)/dh > 0$, even if every coronal loop (of a multi-temperature ensemble) is isothermal in itself. We simulate this effect with differential emission measure distributions observed by {\sl SERTS} for an instrument with a broadband temperature filter such as {\sl Yohkoh/SXT} and find that the apparent temperature increase due to hydrostatic weighting is of order $\Delta T \approx T_0 \times h/r_{\sun}$. We suggest that this effect largely explains the systematic temperature increase in the upper corona reported in recent studies (e.g. by Sturrock et al., Wheatland et al., or Priest et al.), rather than being an intrinsic signature of a coronal heating mechanism.

Both electrons and ions are accelerated in solar flares and carry nonthermal energy from the acceleration site to the chromospheric energy loss site, but the relative amount of energy carried by electrons versus ions is subject of debate. In this {\sl Letter} we test whether the observed energy-dependent timing delays of 20-200 keV HXR emission can be explained in terms of propagating electrons versus protons. For a typical flare, we show that the timing delays of fast ($\lapprox 1$ s) {\sl HXR pulses} is consistent with time-of-flight differences of directly precipitating electrons, while the timing delays of the {\sl smooth HXR} flux is consistent with collisional deflection times of trapped electrons. We show that these HXR timing delays cannot be explained either by $\le 1$ MeV protons (as proposed in a model by Simnett \& Haines 1990), because of their longer propagation and trapping times, or by $\approx 40$ MeV protons (which have the same velocity as $\approx 20$ keV electrons), because of their longer trapping times and the excessive fluxes required to generate the HXRs. Thus, the HXR timing results clearly rule out protons as the primary generators of $\ge 20$ keV HXR emission.