-
A New Approach to Solar Flare Prediction
Authors:
Michael L. Goodman,
Chiman Kwan,
Bulent Ayhan,
Eric L. Shang
Abstract:
All three components of the current density are required to compute the heating rate due to free magnetic energy dissipation. Here we present a first test of a new model developed to determine if the times of increases in the resistive heating rate in active region (AR) photospheres are correlated with the subsequent occurrence of M and X flares in the corona. A data driven, 3 D, non-force-free ma…
▽ More
All three components of the current density are required to compute the heating rate due to free magnetic energy dissipation. Here we present a first test of a new model developed to determine if the times of increases in the resistive heating rate in active region (AR) photospheres are correlated with the subsequent occurrence of M and X flares in the corona. A data driven, 3 D, non-force-free magnetohydrodynamic model restricted to the near-photospheric region is used to compute time series of the complete current density and the resistive heating rate per unit volume $(Q(t))$ in each pixel in neutral line regions (NLRs) of 14 ARs. The model is driven by time series of the magnetic field ${\bf B}$ measured by the Helioseismic & Magnetic Imager on the Solar Dynamics Observatory (SDO) satellite. Spurious Doppler periods due to SDO orbital motion are filtered out of the time series for ${\bf B}$ in every AR pixel. For each AR, the cumulative distribution function (CDF) of the values of the NLR area integral $Q_i(t)$ of $Q(t)$ is found to be a scale invariant power law distribution essentially identical to the observed CDF for the total energy released in coronal flares. This suggests that coronal flares and the photospheric $Q_i$ are correlated, and powered by the same process. The model predicts spikes in $Q_i$ with values orders of magnitude above background values. These spikes are driven by spikes in the non-force free component of the current density. The times of these spikes are plausibly correlated with times of subsequent M or X flares a few hours to a few days later. The spikes occur on granulation scales, and may be signatures of heating in horizontal current sheets. It is also found that the times of relatively large values of the rate of change of the NLR unsigned magnetic flux are also plausibly correlated with the times of subsequent M and X flares, and spikes in $Q_i$.
△ Less
Submitted 27 April, 2020; v1 submitted 3 March, 2020;
originally announced March 2020.
-
Conditions for Photospherically Driven Alfvenic Oscillations to Heat the Solar Chromosphere by Pedersen Current Dissipation
Authors:
Michael L. Goodman
Abstract:
A magnetohydrodynamic model that includes a complete electrical conductivity tensor is used to estimate conditions for photospherically driven, linear, non-plane Alfvenic oscillations extending from the photosphere to the lower corona to drive a chromospheric heating rate due to Pedersen current dissipation that is comparable to the net chromospheric net radiative loss of $\sim 10^7$ ergs-cm…
▽ More
A magnetohydrodynamic model that includes a complete electrical conductivity tensor is used to estimate conditions for photospherically driven, linear, non-plane Alfvenic oscillations extending from the photosphere to the lower corona to drive a chromospheric heating rate due to Pedersen current dissipation that is comparable to the net chromospheric net radiative loss of $\sim 10^7$ ergs-cm$^{-2}$-sec$^{-1}$. The heating rates due to electron current dissipation in the photosphere and corona are also computed. The wave amplitudes are computed self-consistently as functions of an inhomogeneous background (BG) atmosphere. The effects of the conductivity tensor are resolved numerically using a resolution of 3.33 m. The oscillations drive a chromospheric heating flux $F_{Ch} \sim 10^7 - 10^8$ ergs-cm$^{-2}$-sec$^{-1}$ at frequencies $ν\sim 10^2 - 10^3$ mHz for BG magnetic field strengths $B \gtrsim 700$ G and magnetic field perturbation amplitudes $\sim 0.01 - 0.1$ $B$. The total resistive heating flux increases with $ν$. Most heating occurs in the photosphere. Thermalization of Poynting flux in the photosphere due to electron current dissipation regulates the Poynting flux into the chromosphere, limiting $F_{Ch}$. $F_{Ch}$ initially increases with $ν$, reaches a maximum, and then decreases with increasing $ν$ due to increasing electron current dissipation in the photosphere. The resolution needed to resolve the oscillations increases from $\sim 10$ m in the photosphere to $\sim 10$ km in the upper chromosphere, and is proportional to $ν^{-1/2}$. Estimates suggest that these oscillations are normal modes of photospheric flux tubes with diameters $\sim 10-20$ km, excited by magnetic reconnection in current sheets with thicknesses $\sim 0.1$ km.
△ Less
Submitted 30 October, 2014;
originally announced October 2014.
-
Radiating Current Sheets in the Solar Chromosphere
Authors:
Michael L. Goodman,
Philip G. Judge
Abstract:
An MHD model of a Hydrogen plasma with flow, an energy equation, NLTE ionization and radiative cooling, and an Ohm's law with anisotropic electrical conduction and thermoelectric effects is used to self-consistently generate atmospheric layers over a $50$ km height range. A subset of these solutions contain current sheets, and have properties similar to those of the lower and middle chromosphere.…
▽ More
An MHD model of a Hydrogen plasma with flow, an energy equation, NLTE ionization and radiative cooling, and an Ohm's law with anisotropic electrical conduction and thermoelectric effects is used to self-consistently generate atmospheric layers over a $50$ km height range. A subset of these solutions contain current sheets, and have properties similar to those of the lower and middle chromosphere. The magnetic field profiles are found to be close to Harris sheet profiles, with maximum field strengths $\sim 25-150$ G. The radiative flux $F_R$ emitted by individual sheets is $\sim 4.9 \times 10^5 - 4.5 \times 10^6$ ergs-cm$^{-2}$-s$^{-1}$, to be compared with the observed chromospheric emission rate of $\sim 10^7$ ergs-cm$^{-2}$-s$^{-1}$. Essentially all emission is from regions with thicknesses $\sim 0.5 - 13$ km containing the neutral sheet. About half of $F_R$ comes from sub-regions with thicknesses 10 times smaller. A resolution $\lesssim 5-130$ m is needed to resolve the properties of the sheets. The sheets have total H densities $\sim 10^{13}-10^{15}$ cm$^{-3}$. The ionization fraction in the sheets is $\sim 2-20$ times larger, and the temperature is $\sim 2000-3000$ K higher than in the surrounding plasma. The Joule heating flux $F_J$ exceeds $F_R$ by $\sim 4-34 \%$, the difference being balanced in the energy equation mainly by a negative compressive heating flux. Proton Pedersen current dissipation generates $\sim 62-77\%$ of the positive contribution to $F_J$. The remainder of this contribution is due to electron current dissipation near the neutral sheet where the plasma is weakly magnetized.
△ Less
Submitted 4 June, 2014;
originally announced June 2014.
-
Acceleration of Type 2 Spicules in the Solar Chromosphere - 2: Viscous Braking and Upper Bounds on Coronal Energy Input
Authors:
Michael L. Goodman
Abstract:
A magnetohydrodynamic model is used to determine conditions under which the Lorentz force accelerates plasma to type 2 spicule speeds in the chromosphere. The model generalizes a previous model to include a more realistic pre-spicule state, and the vertical viscous force. Two cases of acceleration under upper chromospheric conditions are considered. The magnetic field strength for these cases is <…
▽ More
A magnetohydrodynamic model is used to determine conditions under which the Lorentz force accelerates plasma to type 2 spicule speeds in the chromosphere. The model generalizes a previous model to include a more realistic pre-spicule state, and the vertical viscous force. Two cases of acceleration under upper chromospheric conditions are considered. The magnetic field strength for these cases is <= 12.5 and 25 G. Plasma is accelerated to terminal vertical speeds of 66 and 78 km/s in 100 s, compared with 124 and 397 km/s for the case of zero viscosity. The flows are localized within horizontal diameters ~ 80 and 50 km. The total thermal energy generated by viscous dissipation is ~ 10 times larger than that due to Joule dissipation, but the magnitude of the total cooling due to rarefaction is >~ this energy. Compressive heating dominates during the early phase of acceleration. The maximum energy injected into the corona by type 2 spicules, defined as the energy flux in the upper chromosphere, may largely balance total coronal energy losses in quiet regions, possibly also in coronal holes, but not in active regions. It is proposed that magnetic flux emergence in inter-granular regions drives type 2 spicules.
△ Less
Submitted 11 March, 2014;
originally announced March 2014.