Vector magnetogram data are often used as photospheric boundary conditions for force-free coronal magnetic field extrapolations. In general, however, vector magnetogram data are not consistent with the force-free assumption. In this article, we demonstrate a way to deal with inconsistent boundary data, by generalizing the "self-consistency procedure" of Wheatland & Regnier (2009). In that procedure, the inconsistency is resolved by an iterative process of constructing two solutions based on the values of the force-free parameter alpha on the two polarities of the field in the boundary (the P and N polarities), and taking uncertainty-weighted averages of the boundary alpha values in the P and N solutions. When the alpha values in the P and N regions are very different, the self-consistent solution may lose high alpha values from the boundary conditions. We show how, by altering the weighting of the uncertainties in the P or N boundary conditions, we can preserve high alpha values in the self-consistent solution. The weighted self-consistent extrapolation method is demonstrated on an analytic bipole field and applied to vector magnetogram data taken by the Helioseismic and Magnetic Imager (HMI) instrument for NOAA active region AR 12017 on 2014 March 29.

Magnetic helicity is conserved under ideal magnetohydrodynamics (MHD) and quasi-conserved even under a resistive process. The standard definition for magnetic helicity cannot be applied directly to an open magnetic field in a volume, because it is gauge-dependent. Instead, the relative magnetic helicity is widely used. We find that the energy of a potential magnetic field in a rectangular domain with periodic lateral boundary conditions is less than that of the field with a fixed normal component on all six boundaries. To make use of this lower energy potential field in the analysis of relative magnetic helicity, we introducing a new definition for magnetic helicity for the magnetic field, which involves the periodic potential field. We apply this definition to a sequence of analytic solutions and a numerical simulation. The results show that our new gauge-invariant helicity is very close to the current-carrying part of the relative magnetic helicity of the original magnetic field. We find also that the ratio between the current-carrying helicity and the relative magnetic helicity for the original and our defined relative helicity show different behavior. It seems that the new helicity is more sensitive to the component of the field due to the electric current in the volume, which is the source for instabilities and solar eruptive phenomena.

The distributions of solar flare energies and waiting times have not been described simultaneously by a single physical model, yet. In this research, we investigate whether recent avalanche models can describe the distributions for both the released energies and waiting times of flares in an active region. Flaring events are simulated using the modified Lu and Hamilton model (Lu and Hamilton (1991), ApJ, 380, 89) and also the optimized model (Farhang et al. (2018), ApJ, 859, 41). Applying a probability balance equation approach, we study the statistics of the simulated flaring events and investigate the origin of the observed power law in the flare frequency-size distribution. The results indicate that the power law originates in the distribution of transition rates (the distribution of the probabilities of transitions between different energies) rather than the distribution of the energy of the active region. It is also observed that the waiting-time distribution of simulated flaring events follows a q-exponential function which approximates a simple Poisson distribution.

During eruptive flares, vector magnetograms show increasing horizontal magnetic field and downward Lorentz force in the Sun's photosphere around the polarity-inversion line. Such behavior has often been associated with the implosion conjecture and interpreted as the result of either momentum conservation while the eruption moves upward, or of the contraction of flare loops. We characterize the physical origin of these observed behaviors by analyzing a generic 3D MHD simulation of an eruptive flare. Even though the simulation was undesigned to recover the magnetic field and Lorentz force properties, it is fully consistent with them, and it provides key additional informations to understand them. The area where the magnetic field increases gradually develops between current ribbons, which spread away from each other and are connected to the coronal region. This area is merely the footprint of the coronal post-flare loops, whose contraction increases their shear field component and the magnetic energy density in line with the ideal induction equation. For simulated data, we computed the Lorentz force density map by applying the method used in observations. We obtained increase of the downward component of the Lorentz force density around the PIL -consistent with observations. However, this significantly differs from the Lorentz force density maps obtained directly from the 3D magnetic field and current. These results altogether question previous interpretations based on the implosion conjecture and momentum conservation with the CME, and rather imply that the observed increases in photospheric horizontal magnetic fields result from the reconnection-driven contraction of sheared flare-loops.

Rapid and stepwise changes of the magnetic field are often observed during flares but cannot be explained by models yet. Using a 45 min sequence of SDO/HMI 135 s fast-cadence vector magnetograms of the X1 flare on 2014-03-29 we construct, at each timestep, nonlinear force-free models for the coronal magnetic field. Observed flare-related changes in the line-of-sight magnetic field $B_{\rm LOS}$ at the photosphere and chromosphere are compared with changes in the magnetic fields in the models. We find a moderate agreement at the photospheric layer (the basis for the models), but no agreement at chromospheric layers. The observed changes at the photosphere and chromosphere are surprisingly different, and are unlikely to be reproduced by a force-free model. The observed changes are likely to require a change in the magnitude of the field, not just in its direction.

Flares produce sudden and permanent changes in the horizontal photospheric magnetic field. In particular flares generally produce increased magnetic shear in the photospheric field along the neutral line. Recent observations show also that flares can produce sudden photospheric motion. We present a model for the observed changes as the response of the photosphere to a large-amplitude shear Alfv\'{e}n wave propagating down from the corona on either side of the neutral line. The Alfv\'{e}nic front is assumed to impact the photosphere close to the neutral line first, and then successively further away with time, such that the line of impact coincides with the flare ribbon. The wave introduces magnetic shear and velocity shear. The magnetic shear introduced at the photosphere has the same sign on either side of the neutral line, while the velocity shear has the opposite sign. We discuss the possibility that this process is responsible for particle acceleration in flares.

Solar flares are an abrupt release of magnetic energy in the Sun's atmosphere due to reconnection of the coronal magnetic field. This occurs in response to turbulent flows at the photosphere which twist the coronal field. Similar to earthquakes, solar flares represent the behavior of a complex system, and expectedly their energy distribution follows a power law. We present a statistical model based on the principle of minimum energy in a coronal loop undergoing magnetic reconnection, which is described as an avalanche process. We show the distribution of peaks for the flaring events in this self-organized critical system is scale-free. The obtained power law index 1.84 $ \pm $ 0.02 for the peaks is in good agreement with satellite observations of soft X-ray flares. The principle of minimum energy can be applied for general avalanche models to describe many other phenomena.

Quantitative models of sunspot and starspot decay predict the timescale of magnetic diffusion and may yield important constraints in stellar dynamo models. Motivated by recent measurements of starspot lifetimes, we investigate the disintegration of a magnetic flux tube by nonlinear diffusion. Previous theoretical studies are extended by considering two physically motivated functional forms for the nonlinear diffusion coefficient $D$: an inverse power-law dependence $D \propto B^{-\nu}$ and a step-function dependence of $D$ on the magnetic field magnitude $B$. Analytical self-similar solutions are presented for the power-law case, including solutions exhibiting "superfast" diffusion. For the step-function case, the heat-balance integral method yields approximate solutions, valid for moderately suppressed diffusion in the spot. The accuracy of the resulting solutions is confirmed numerically, using a method which provides an accurate description of long-time evolution by imposing boundary conditions at infinite distance from the spot. The new models may allow insight into differences and similarities between sunspots and starspots.

Since the early 1980s, decimetric spike bursts have been attributed to electron cyclotron maser emission (ECME) by the electrons that produce hard X-ray bursts as they precipitate into the chromosphere in the impulsive phase of a solar flare. Spike bursts are regarded as analogous to the auroral kilometric radiation (AKR), which is associated with the precipitation of auroral electrons in a geomagnetic substorm. Originally, a loss-cone-driven version of ECME, developed for AKR, was applied to spike bursts, but it is now widely accepted that a different, horseshoe-driven, version of EMCE applies to AKR. We explore the implications of the assumption that horseshoe-driven ECME also applies to spike bursts. We develop a 1D model for the acceleration of the electrons by a parallel electric field, and show that under plausible assumptions it leads to a horseshoe distribution of electrons in a solar flare. A second requirement for horseshoe-driven ECME is an extremely low plasma density, referred to as a density cavity. We argue that a coronal density cavity should develop in association with a hard X-ray burst, and that such a density cavity can overcome a long-standing problem with the escape of ECME through the second-harmonic absorption layer. Both the horseshoe distribution and the associated coronal density cavity are highly localized, and could not be resolved in the statistically large number of local precipitation regions needed to explain a hard X-ray burst. The model highlights the "number problem" in the supply of the electrons needed to explain a hard X-ray burst.

Prediction of solar flares is an important task in solar physics. The occurrence of solar flares is highly dependent on the structure and the topology of solar magnetic fields. A new method for predicting large (M and X class) flares is presented, which uses machine learning methods applied to the Zernike moments of magnetograms observed by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) for a period of six years from 2 June 2010 to 1 August 2016. Magnetic field images consisting of the radial component of the magnetic field are converted to finite sets of Zernike moments and fed to the Support Vector Machine (SVM) classifier. Zernike moments have the capability to elicit unique features from any 2-D image, which may allow more accurate classification. The results indicate whether an arbitrary active region has the potential to produce at least one large flare. We show that the majority of large flares can be predicted within 48 hours before their occurrence, with only 10 false negatives out of 385 flaring active region magnetograms, and 21 false positives out of 179 non-flaring active region magnetograms. Our method may provide a useful tool for prediction of solar flares which can be employed alongside other forecasting methods.

Solar flares produce radiation which can have an almost immediate effect on the near-Earth environment, making it crucial to forecast flares in order to mitigate their negative effects. The number of published approaches to flare forecasting using photospheric magnetic field observations has proliferated, with varying claims about how well each works. Because of the different analysis techniques and data sets used, it is essentially impossible to compare the results from the literature. This problem is exacerbated by the low event rates of large solar flares. The challenges of forecasting rare events have long been recognized in the meteorology community, but have yet to be fully acknowledged by the space weather community. During the interagency workshop on "all clear" forecasts held in Boulder, CO in 2009, the performance of a number of existing algorithms was compared on common data sets, specifically line-of-sight magnetic field and continuum intensity images from MDI, with consistent definitions of what constitutes an event. We demonstrate the importance of making such systematic comparisons, and of using standard verification statistics to determine what constitutes a good prediction scheme. When a comparison was made in this fashion, no one method clearly outperformed all others, which may in part be due to the strong correlations among the parameters used by different methods to characterize an active region. For M-class flares and above, the set of methods tends towards a weakly positive skill score (as measured with several distinct metrics), with no participating method proving substantially better than climatological forecasts.

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.

Electric currents in solar active regions are thought to provide the energy released via magnetic reconnection in solar flares. Vertical electric current densities $J_z$ at the photosphere may be estimated from vector magnetogram data, subject to substantial uncertainties. The values provide boundary conditions for nonlinear force- free modelling of active region magnetic fields. A method is presented for estimating values of $J_z$ taking into account uncertainties in vector magnetogram field values, and minimizing $J_z^2$ across the active region. The method is demonstrated using the boundary values of the field for a force-free twisted bipole, with the addition of noise at randomly chosen locations.

Disintegration of sunspots (and starspots) by fluxtube erosion, originally proposed by Simon and Leighton, is considered. A moving boundary problem is formulated for a nonlinear diffusion equation that describes the sunspot magnetic field profile. Explicit expressions for the sunspot decay rate and lifetime by turbulent erosion are derived analytically and verified numerically. A parabolic decay law for the sunspot area is obtained. For moderate sunspot magnetic field strengths, the predicted decay rate agrees with the results obtained by Petrovay and Moreno-Insertis. The new analytical and numerical solutions significantly improve the quantitative description of sunspot and starspot decay by turbulent erosion.

The shapes of solar coronal loops are sensitive to the presence of electrical currents that are the carriers of the nonpotential energy available for impulsive activity. We use this information in a new method for modeling the coronal magnetic field of AR 11158 as a nonlinear force-free field (NLFFF). The observations used are coronal images around time of major flare activity on 2011/02/15, together with the surface line-of-sight magnetic field measurements. The data are from the Helioseismic and Magnetic Imager and Atmospheric Imaging Assembly (HMI and AIA, respectively) onboard the Solar Dynamics Observatory (SDO). The model fields are constrained to approximate the coronal loop configurations as closely as possible, while also subject to the force-free constraints. The method does not use transverse photospheric magnetic field components as input, and is thereby distinct from methods for modeling NLFFFs based on photospheric vector magnetograms. We validate the method using observations of AR 11158 at a time well before major flaring, and subsequently review the field evolution just prior to and following an X2.2 flare and associated eruption. The models indicate that the energy released during the instability is about $1\times10^{32}$ erg, consistent with what is needed to power such a large eruptive flare. Immediately prior to the eruption the model field contains a compact sigmoid bundle of twisted flux that is not present in the post-eruption models, which is consistent with the observations. The core of that model structure is twisted by $\approx0.9$ full turns about its axis.

We present a code for solving the nonlinear force-free equations in spherical polar geometry, with the motivation of modeling the magnetic field in the corona. The code is an implementation of the Grad-Rubin method. Our method is applicable to a spherical domain of arbitrary angular size. The implementation is based on a global spectral representation for the magnetic field which makes no explicit assumptions about the form of the magnetic field at the transverse boundaries of the domain. We apply the code to a bipolar test case with analytic boundary conditions, and we demonstrate the convergence of the Grad-Rubin method, and the self-consistency of the resulting numerical solution.

Bulk energization of electrons to $10\,-\,20\,$keV in solar flares is attributed to dissipation of Alfv\'en waves that transport energy and potential downward to an acceleration region near the chromosphere. The acceleration involves the parallel electric field that develops in the limit of inertial Alfv\'en waves (IAWs). A two-potential model for IAWs is used to relate the parallel potential to the cross-field potential transported by the waves. We identify a maximum parallel potential in terms of a maximum current density that corresponds to the threshold for the onset of anomalous resistivity. This maximum is of order $10\,$kV when the threshold is that for the Buneman instability. We argue that this restricts the cross-field potential in an Alfv\'en wave to about $10\,$kV. Effective dissipation requires a large number of up- and down-current paths associated with multiple Alfv\'en waves. The electron acceleration occurs in localized, transient, anomalously-conducting regions (LTACRs) and is associated with the parallel electric field determined by Ohm's law with an anomalous resistivity. We introduce an idealized model in which the LTACRs are (upward-)current sheets, a few skin depths in thickness, separated by much-larger regions of weaker return current. We show that this model can account semi-quantitatively for bulk energization.

Alfv\'en waves play three related roles in the impulsive phase of a solar flare: they transport energy from a generator region to an acceleration region; they map the cross-field potential (associated with the driven energy release) from the generator region onto the acceleration region; and within the acceleration region they damp by setting up a parallel electric field that accelerates electrons and transfers the wave energy to them. The Alfv\'en waves may also be regarded as setting up new closed current loops, with field-aligned currents that close across field lines at boundaries. A model is developed for large-amplitude Alfv\'en waves that shows how Alfv\'en waves play these roles in solar flares. A picket-fence structure for the current flow is incorporated into the model to account for the "number problem" and the energy of the accelerated electrons.

Recently Pop ({\it Solar Phys.} {\bf 276}, 351, 2012) identified a Laplace (or double exponential) distribution in the number of days with a given absolute value in the change over a day, in sunspot number, for days on which the sunspot number does change. We show this phenomenological rule has a physical origin attributable to sunspot formation, evolution, and decay, rather than being due to the changes in sunspot number caused by groups rotating onto and off the visible disc. We also demonstrate a simple method to simulate daily sunspot numbers over a solar cycle using the \cite{2012SoPh..276..351P} result, together with a model for the cycle variation in the mean sunspot number. The procedure is applied to three recent solar cycles. We check that the simulated sunspot numbers reproduce the observed distribution of daily changes over those cycles.

The coronal magnetic field cannot be directly observed, but in principle it can be reconstructed from the comparatively well observed photospheric magnetic field. A popular approach uses a nonlinear force-free model. Non-magnetic forces at the photosphere are significant meaning the photospheric data are inconsistent with the force-free model, and this causes problems with the modeling (De Rosa et al., Astrophys. J. 696, 1780, 2009). In this paper we present a numerical implementation of the Grad-Rubin method for reconstructing the coronal magnetic field using a magnetostatic model. This model includes a pressure force and a non-zero magnetic Lorentz force. We demonstrate our implementation on a simple analytic test case and obtain the speed and numerical error scaling as a function of the grid size.

Presently, many models of the coronal magnetic field rely on photospheric vector magnetograms but these data have been shown to be problematic as the sole boundary information for nonlinear force-free field (NLFFF) extrapolations. Magnetic fields in the corona manifest themselves in high-energy images (X-rays and EUV) in the shapes of coronal loops, providing an additional constraint that at present is not used due to the mathematical complications of incorporating such input into numerical models. Projection effects and the limited number of usable loops further complicate the use of coronal information. We develop and test an algorithm to use images showing coronal loops in the modeling of the solar coronal magnetic field. We first fit projected field lines with field lines of constant-\als force-free fields to approximate the three-dimensional distribution of currents in the corona along a sparse set of trajectories. We then apply a Grad-Rubin-like iterative technique to obtain a volume-filling nonlinear force-free model of the magnetic field, modifying method presented in \citet{Wheatland2007}. We thoroughly test the technique on known analytical and solar-like model magnetic fields previously used for comparing different extrapolation techniques \citep{Schrijver2006, Schrijver2008} and compare the results with those obtained by presently available methods that rely only on the photospheric data. We conclude that we have developed a functioning method of modeling the coronal magnetic field by combining the line-of-sight component of photospheric magnetic field with information from coronal images. Vector magnetograms over the full or partial photospheric boundary of the numerical domain could optionally be used.

A Bayesian method for forecasting solar cycles is presented. The approach combines a Fokker--Planck description of short--timescale (daily) fluctuations in sunspot number (\citeauthor{NobleEtAl2011}, 2011, \apj{} \textbf{732}, 5) with information from other sources, such as precursor and/or dynamo models. The forecasting is illustrated in application to two historical cycles (cycles 19 and 20), and then to the current solar cycle (cycle 24). The new method allows the prediction of quantiles, i.e. the probability that the sunspot number falls outside large or small bounds at a given future time. It also permits Monte Carlo simulations to identify the expected size and timing of the peak daily sunspot number, as well as the smoothed sunspot number for a cycle. These simulations show how the large variance in daily sunspot number determines the actual reliability of any forecast of the smoothed maximum of a cycle. For cycle 24 we forecast a maximum daily sunspot number of $166\pm 24$, to occur in March 2013, and a maximum value of the smoothed sunspot number of $66\pm5$, indicating a very small solar cycle.

The NOAA active region AR 11029 was a small but highly active sunspot region which produced 73 GOES soft X-ray flares. The flares appear to show a departure from the well known power-law frequency-size distribution. Specifically, too few GOES C-class and no M-class flares were observed by comparison with a power-law distribution (Wheatland in Astrophys. J. 710, 1324, 2010). This was conjectured to be due to the region having insufficient magnetic energy to power large events. We construct nonlinear force-free extrapolations of the coronal magnetic field of active region AR 11029 using data taken on 24 October by the SOLIS Vector-SpectroMagnetograph (SOLIS/VSM), and data taken on 27 October by the Hinode Solar Optical Telescope SpectroPolarimeter (Hinode/SP). Force-free modeling with photospheric magnetogram data encounters problems because the magnetogram data are inconsistent with a force-free model, and we employ a recently developed `self-consistency' procedure which addresses this and accommodates uncertainties in the boundary data (Wheatland and Regnier in Astrophys. J. 700, L88, 2009). We calculate the total energy and free energy of the self-consistent solution and find that the free energy was 4x10^29 erg on 24 October, and 7x10^31 erg on 27 October. An order of magnitude scaling between RHESSI non-thermal energy and GOES peak X-ray flux is established from a sample of flares from the literature and is used to estimate flare energies from observed GOES peak X-ray flux. Based on the scaling, we conclude that the estimated free energy of AR 11029 on 27 October when the flaring rate peaked is sufficient to power M-class or X-class flares, and hence the modeling does not appear to support the hypothesis that the absence of large flares is due to the region having limited energy.

Sunspot numbers exhibit large short-timescale (daily-monthly) variation in addition to longer timescale variation due to solar cycles. A formal statistical framework is presented for estimating and forecasting randomness in sunspot numbers on top of deterministic (including chaotic) models for solar cycles. The Fokker-Planck approach formulated assumes a specified long-term or secular variation in sunspot number over an underlying solar cycle via a driver function. The model then describes the observed randomness in sunspot number on top of this driver function. We consider a simple harmonic choice for the driver function, but the approach is general and can easily be extended to include other drivers which account for underlying physical processes and/or empirical features of the sunspot numbers. The framework is consistent during both solar maximum and minimum, and requires no parameter restrictions to ensure non-negative sunspot numbers. Model parameters are estimated using statistically optimal techniques. The model agrees both qualitatively and quantitatively with monthly sunspot data even with the simplistic representation of the periodic solar cycle. This framework should be particularly useful for solar cycle forecasters and is complementary to existing modeling techniques. An analytic approximation for the Fokker-Planck equation is presented, which is analogous to the Euler approximation, which which allows for efficient maximum likelihood estimation of large data sets and/or when using difficult to evaluate driver functions.

A nonlinear force-free solution is constructed for the coronal magnetic field in NOAA solar active region AR 10953 based on a photospheric vector magnetogram derived from Hinode satellite observations on 30 April 2007, taking into account uncertainties in the boundary data and using improved methods for merging multiple-instrument data. The solution demonstrates the "self-consistency" procedure of Wheatland & Regnier (2009), for the first time including uncertainties. The self-consistency procedure addresses the problem that photospheric vector magnetogram data are inconsistent with the force-free model, and in particular that the boundary conditions on vertical electric current density are over-specified and permit the construction of two different nonlinear force-free solutions. The procedure modifies the boundary conditions on current density during a sequence of cycles until the two nonlinear force-free solutions agree. It hence constructs an accurate single solution to the force-free model, with boundary values close, but not matched exactly, to the vector magnetogram data. The inclusion of uncertainties preserves the boundary conditions more closely at points with smaller uncertainties. The self-consistent solution obtained for active region AR 10953 is significantly non-potential, with magnetic energy E/E_0 = 1.08, where E_0 is the energy of the reference potential (current-free) magnetic field. The self-consistent solution is shown to be robust against changes in the details of the construction of the two force-free models at each cycle. This suggests that reliable nonlinear force-free modeling of active regions is possible if uncertainties in vector magnetogram boundary data are included.

A time-dependent model for the energy of a flaring solar active region is presented based on a stochastic jump-transition model (Wheatland and Glukhov 1998; Wheatland 2008; Wheatland 2009). The magnetic free energy of the model active region varies in time due to a prescribed (deterministic) rate of energy input and prescribed (random) flare jumps downwards in energy. The model has been shown to reproduce observed flare statistics, for specific time-independent choices for the energy input and flare transition rates. However, many solar active regions exhibit time variation in flare productivity, as exemplified by NOAA active region AR 11029 (Wheatland 2010). In this case a time-dependent model is needed. Time variation is incorporated for two cases: 1. a step change in the rates of flare jumps; and 2. a step change in the rate of energy supply to the system. Analytic arguments are presented describing the qualitative behavior of the system in the two cases. In each case the system adjusts by shifting to a new stationary state over a relaxation time which is estimated analytically. The new model retains flare-like event statistics. In each case the frequency-energy distribution is a power law for flare energies less than a time-dependent rollover set by the largest energy the system is likely to attain at a given time. For Case 1, the model exhibits a double exponential waiting-time distribution, corresponding to flaring at a constant mean rate during two intervals (before and after the step change), if the average energy of the system is large. For Case 2 the waiting-time distribution is a simple exponential, again provided the average energy of the system is large. Monte Carlo simulations of Case~1 are presented which confirm the analytic estimates. The simulation results provide a qualitative model for observed flare statistics in active region AR 11029.

Active region 11029 was a small, highly flare-productive solar active region observed at a time of extremely low solar activity. The region produced only small flares: the largest of the $>70$ Geostationary Observational Environmental Satellite (GOES) events for the region has a peak 1--$8{\AA}$ flux of $2.2\times 10^{-6} {\rm W} {\rm m}^{-2}$ (GOES C2.2). The background-subtracted GOES peak-flux distribution suggests departure from power-law behavior above $10^{-6} {\rm W} {\rm m}^{-2}$, and a Bayesian model comparison strongly favors a power-law plus rollover model for the distribution over a simple power-law model. The departure from the power law is attributed to this small active region having a finite amount of energy. The rate of flaring in the region varies with time, becoming very high for two days coinciding with the onset of an increase in complexity of the photospheric magnetic field. The observed waiting-time distribution for events is consistent with a piecewise-constant Poisson model. These results present challenges for models of flare statistics and of energy balance in solar active regions.

Type III solar radio storms, observed at frequencies below approximately 16 MHz by space borne radio experiments, correspond to the quasi-continuous, bursty emission of electron beams onto open field lines above active regions. The mechanisms by which a storm can persist in some cases for more than a solar rotation whilst exhibiting considerable radio activity are poorly understood. To address this issue, the statistical properties of a type III storm observed by the STEREO/WAVES radio experiment are presented, examining both the brightness distribution and (for the first time) the waiting-time distribution. Single power law behavior is observed in the number distribution as a function of brightness; the power law index is approximately 2.1 and is largely independent of frequency. The waiting-time distribution is found to be consistent with a piecewise-constant Poisson process. This indicates that during the storm individual type III bursts occur independently and suggests that the storm dynamics are consistent with avalanche type behavior in the underlying active region.

Nonlinear force-free solutions for the magnetic field in the solar corona constructed using photospheric vector magnetic field boundary data suffer from a basic problem: the observed boundary data are inconsistent with the nonlinear force-free model. Specifically, there are two possible choices of boundary conditions on vertical current provided by the data, and the two choices lead to different force-free solutions. A novel solution to this problem is described. Bayesian probability is used to modify the boundary values on current density, using field-line connectivity information from the two force-free solutions and taking into account uncertainties, so that the boundary data are more consistent with the two nonlinear force-free solutions. This procedure may be iterated until a set of self-consistent boundary data (the solutions for the two choices of boundary conditions are the same) is achieved. The approach is demonstrated to work in application to Hinode/SOT observations of NOAA active region 10953.

Nonlinear force-free field (NLFFF) models are thought to be viable tools for investigating the structure, dynamics and evolution of the coronae of solar active regions. In a series of NLFFF modeling studies, we have found that NLFFF models are successful in application to analytic test cases, and relatively successful when applied to numerically constructed Sun-like test cases, but they are less successful in application to real solar data. Different NLFFF models have been found to have markedly different field line configurations and to provide widely varying estimates of the magnetic free energy in the coronal volume, when applied to solar data. NLFFF models require consistent, force-free vector magnetic boundary data. However, vector magnetogram observations sampling the photosphere, which is dynamic and contains significant Lorentz and buoyancy forces, do not satisfy this requirement, thus creating several major problems for force-free coronal modeling efforts. In this article, we discuss NLFFF modeling of NOAA Active Region 10953 using Hinode/SOT-SP, Hinode/XRT, STEREO/SECCHI-EUVI, and SOHO/MDI observations, and in the process illustrate the three such issues we judge to be critical to the success of NLFFF modeling: (1) vector magnetic field data covering larger areas are needed so that more electric currents associated with the full active regions of interest are measured, (2) the modeling algorithms need a way to accommodate the various uncertainties in the boundary data, and (3) a more realistic physical model is needed to approximate the photosphere-to-corona interface in order to better transform the forced photospheric magnetograms into adequate approximations of nearly force-free fields at the base of the corona. We make recommendations for future modeling efforts to overcome these as yet unsolved problems.

A Monte-Carlo approach to solving a stochastic jump transition model for active-region energy (Wheatland and Glukhov, Astrophys. J. 494, 1998; Wheatland, Astrophys. J. 679, 2008) is described. The new method numerically solves the stochastic differential equation describing the model, rather than the equivalent master equation. This has the advantages of allowing more efficient numerical solution, the modelling of time-dependent situations, and investigation of details of event statistics. The Monte-Carlo approach is illustrated by application to a Gaussian test case, and to the class of flare-like models presented in Wheatland (2008), which are steady-state models with constant rates of energy supply, and power-law distributed jump transition rates. These models have two free parameters: an index ($\delta $), which defines the dependence of the jump transition rates on active-region energy, and a non-dimensional ratio ($\overline{r})$ of total flaring rate to rate of energy supply. For $\overline{r}\ll 1$ the non-dimensional mean energy $<\overline{E}>$ of the active-region satisfies $<\overline{E}> \gg 1$, resulting in a power-law distribution of flare events over many decades in energy. The Monte-Carlo method is used to explore the behavior of the waiting-time distributions for the flare-like models. The models with $\delta\neq 0$ are found to have waiting times which depart significantly from simple Poisson behavior when $<\overline{E}> \gg 1$. The original model from Wheatland and Glukhov (1998), with $\delta=0$ (no dependence of transition rates on active-region energy), is identified as being most consistent with observed flare statistics.

A stochastic model for the energy of a flaring solar active region is presented, generalising and extending the approach of Wheatland & Glukhov (1998). The probability distribution for the free energy of an active region is described by the solution to a master equation involving deterministic energy input and random jump transitions downwards in energy (solar flares). It is shown how two observable distributions, the flare frequency-energy distribution and the flare waiting-time distribution, may be derived from the steady-state solution to the master equation, for given choices for the energy input and for the rates of flare transitions. An efficient method of numerical solution of the steady-state master equation is presented. Solutions appropriate for flaring, involving a constant rate of energy input and power-law distributed jump transition rates, are numerically investigated. The flare-like solutions exhibit power-law flare frequency-energy distributions below a high energy rollover, set by the largest energy the active region is likely to have. The solutions also exhibit approximately exponential (i.e. Poisson) waiting-time distributions, despite the rate of flaring depending on the free energy of the system.

Solar flares and coronal mass ejections are associated with rapid changes in field connectivity and powered by the partial dissipation of electrical currents in the solar atmosphere. A critical unanswered question is whether the currents involved are induced by the motion of pre-existing atmospheric magnetic flux subject to surface plasma flows, or whether these currents are associated with the emergence of flux from within the solar convective zone. We address this problem by applying state-of-the-art nonlinear force-free field (NLFFF) modeling to the highest resolution and quality vector-magnetographic data observed by the recently launched Hinode satellite on NOAA Active Region 10930 around the time of a powerful X3.4 flare. We compute 14 NLFFF models with 4 different codes and a variety of boundary conditions. We find that the model fields differ markedly in geometry, energy content, and force-freeness. We discuss the relative merits of these models in a general critique of present abilities to model the coronal magnetic field based on surface vector field measurements. For our application in particular, we find a fair agreement of the best-fit model field with the observed coronal configuration, and argue (1) that strong electrical currents emerge together with magnetic flux preceding the flare, (2) that these currents are carried in an ensemble of thin strands, (3) that the global pattern of these currents and of field lines are compatible with a large-scale twisted flux rope topology, and (4) that the ~10^32 erg change in energy associated with the coronal electrical currents suffices to power the flare and its associated coronal mass ejection.

According to current practice, the results of each run of a radiochemical solar neutrino experiment comprise an estimate of the flux and upper and lower error estimates. These estimates are derived by a maximum-likelihood procedure from the times of decay events in the analysis chamber. This procedure has the following shortcomings: (a) Published results sometimes include negative flux estimates. (b) Even if the flux estimate is non-negative, the probability distribution function implied by the flux and error estimates will extend into negative territory; and (c) The overall flux estimate derived from the results of a sequence of runs may differ substantially from an estimate made by a global analysis of all of the timing data taken together. These defects indicate that the usual packaging of data in radiochemical solar neutrino experiments provides an inadequate summary of the data, which implies a loss of information. This article reviews this problem from a Bayesian perspective, and suggests an alternative scheme for the packaging of radiochemical solar neutrino data, which is we believe free from the above objections.

The MHD virial theorem may be used to estimate the magnetic energy of active regions based on vector magnetic fields measured at the photosphere or chromosphere. However, the virial estimate depends on the measured vector magnetic field being force-free. Departure from force-freeness leads to an unknown systematic error in the virial energy estimate, and an origin dependence of the result. We present a method for estimating the systematic error by assuming that magnetic forces are confined to a thin layer near the photosphere. If vector magnetic field measurements are available at two levels in the low atmosphere (e.g. the photosphere and the chromosphere), the systematic error may be directly calculated using the observed horizontal and vertical field gradients, resulting in an energy estimate which is independent of the choice of origin. If (as is generally the case) measurements are available at only one level, the systematic error may be approximated using the observed horizontal field gradients together with a simple linear force-free model for the vertical field gradients. The resulting `improved' virial energy estimate is independent of the choice of origin, but depends on the choice of the model for the vertical field gradients, i.e. the value of the linear force-free parameter $\alpha$. This procedure is demonstrated for five vector magnetograms, including a chromospheric magnetogram.

A Bayesian approach to solar flare prediction has been developed, which uses only the event statistics of flares already observed. The method is simple, objective, and makes few ad hoc assumptions. It is argued that this approach should be used to provide a baseline prediction for certain space weather purposes, upon which other methods, incorporating additional information, can improve. A practical implementation of the method for whole-Sun prediction of Geostationary Observational Environment Satellite (GOES) events is described in detail, and is demonstrated for 4 November 2003, the day of the largest recorded GOES flare. A test of the method is described based on the historical record of GOES events (1975-2003), and a detailed comparison is made with US National Oceanic and Atmospheric Administration (NOAA) predictions for 1987-2003. Although the NOAA forecasts incorporate a variety of other information, the present method out-performs the NOAA method in predicting mean numbers of event days, for both M-X and X events. Skill scores and other measures show that the present method is slightly less accurate at predicting M-X events than the NOAA method, but substantially more accurate at predicting X events, which are important contributors to space weather.

In solar physics, especially in exploratory stages of research, it is often necessary to compare the power spectra of two or more time series. One may, for instance, wish to estimate what the power spectrum of the combined data sets might have been, or one may wish to estimate the significance of a particular peak that shows up in two or more power spectra. One may also on occasion need to search for a complex of peaks in a single power spectrum, such as a fundamental and one or more harmonics, or a fundamental plus sidebands, etc. Visual inspection can be revealing, but it can also be misleading. This leads one to look for one or more ways of forming statistics, which readily lend themselves to significance estimation, from two or more power spectra. We derive formulas for statistics formed from the sum, the minimum, and the product of two or more power spectra. A distinguishing feature of our formulae is that, if each power spectrum has an exponential distribution, each statistic also has an exponential distribution.

A test of a new Bayesian approach to solar flare prediction (Wheatland 2004a) is presented. The approach uses the past history of flaring together with phenomenological rules of flare statistics to make a prediction for the probability of occurrence of a large flare within an interval of time, or to refine an initial prediction (which may incorporate other information). The test of the method is based on data from the Geostationary Observational Environmental Satellites (GOES), and involves whole-Sun prediction of soft X-ray flares for 1976-2003. The results show that the method somewhat over-predicts the probability of all events above a moderate size, but performs well in predicting large events.

A number of methods of flare prediction rely on classification of physical characteristics of an active region, in particular optical classification of sunspots, and historical rates of flaring for a given classification. However these methods largely ignore the number of flares the active region has already produced, in particular the number of small events. The past history of occurrence of flares (of all sizes) is an important indicator to future flare production. We present a Bayesian approach to flare prediction, which uses the flaring record of an active region together with phenomenological rules of flare statistics to refine an initial prediction for the occurrence of a big flare during a subsequent period of time. The initial prediction is assumed to come from one of the extant methods of flare prediction. The theory of the method is outlined, and simulations are presented to show how the refinement step of the method works in practice.

Circuit models involving bulk currents and inductances are often used to estimate the energies of coronal magnetic field configurations, in particular configurations associated with solar flares. The accuracy of circuit models is tested by comparing calculated energies of linear force-free fields with specified boundary conditions with corresponding circuit estimates. The circuit models are found to provide reasonable (order of magnitude) estimates for the energies of the non-potential components of the fields, and to reproduce observed functional dependences of the energies. However, substantial departure from the circuit estimates is observed for large values of the force-free parameter, and this is attributed to the influence of the non-potential component of the field on the path taken by the current.

It is often necessary to compare the power spectra of two or more time series: one may, for instance, wish to estimate what the power spectrum of the combined data sets might have been, or one may wish to estimate the significance of a particular peak that shows up in two or more power spectra. Also, one may occasionally need to search for a complex of peaks in a single power spectrum, such as a fundamental and one or more harmonics, or a fundamental plus sidebands, etc. Visual inspection can be revealing, but it can also be misleading. This leads one to look for one or more ways of forming statistics, which readily lend themselves to significance estimation, from two or more power spectra. The familiar chi-square statistic provides a convenient mechanism for combining variables drawn from normal distributions, and one may generalize the chi-square statistic to be any function of any number of variables with arbitrary distributions. In dealing with power spectra, we are interested mainly in exponential distributions. One well-known statistic, formed from the sum of two or more variables with exponential distributions, satisfies the gamma distribution. We show that a transformation of this statistic has the convenient property that it has an exponential distribution. We introduce two additional statistics formed from two or more variables with exponential distributions. For certain investigations, we may wish to study the minimum power (as a function of frequency) drawn from two or more power spectra. In other investigations, it may be helpful to study the product of the powers. We give numerical examples and an example drawn from our solar-neutrino research.

The distribution of times $\Delta t$ between coronal mass ejections (CMEs) in the Large Angle and Spectrometric Coronagraph (LASCO) CME catalog for the years 1996-2001 is examined. The distribution exhibits a power-law tail $\propto (\Delta t)^{\gamma}$ with an index $\gamma\approx -2.36\pm 0.11$ for large waiting times ($\Delta t>10 {\rm hours}$). The power-law index of the waiting-time distribution varies with the solar cycle: for the years 1996-1998 (a period of low activity), the power-law index is $\gamma\approx-1.86\pm 0.14$, and for the years 1999-2001 (a period of higher activity), the index is $\gamma\approx-2.98\pm 0.20$. The observed CME waiting-time distribution, and its variation with the cycle, may be understood in terms of CMEs occurring as a time-dependent Poisson process. The CME waiting-time distribution is compared with that for greater than C1 class solar flares in the Geostationary Operational Environmental Satellite (GOES) catalog for the same years. The flare and CME waiting-time distributions exhibit power-law tails with very similar indices and time variation.

The question of whether flares occur as a Poisson process has important consequences for flare physics. Recently Lepreti et al. presented evidence for local departure from Poisson statistics in the Geostationary Operational Environmental Satellite (GOES) X-ray flare catalog. Here it is argued that this effect arises from a selection effect inherent in the soft X-ray observations; namely that the slow decay of enhanced flux following a large flare makes detection of subsequent flares less likely. It is also shown that the power-law tail of the GOES waiting-time distribution varies with the solar cycle. This counts against any intrinsic significance to the appearance of a power law, or to the value of its index.

It was recently pointed out that the distribution of times between solar flares (the flare waiting-time distribution) follows a power law, for long waiting times. Based on 25 years of soft X-ray flares observed by Geostationary Operational Environmental Satellite (GOES) instruments it is shown that 1. the waiting-time distribution of flares is consistent with a time-dependent Poisson process, and 2. the fraction of time the Sun spends with different flaring rates approximately follows an exponential distribution. The second result is a new phenomenological law for flares. It is shown analytically how the observed power-law behavior of the waiting times originates in the exponential distribution of flaring rates. These results are argued to be consistent with a non-stationary avalanche model for flares.

Recent analysis of the Homestake data indicates that the solar neutrino flux contains a periodic variation that may be attributed to rotational modulation occurring deep in the solar interior, either in the tachocline or in the radiative zone. This paper presents an analysis of GALLEX data that yields supporting evidence of this rotational modulation at the 0.1% significance level. The depth of modulation inferred from the rotational signature is large enough to explain the neutrino deficit. The Rieger 157-day periodicity, first discovered in solar gamma-ray flares, is present also in Homestake data. A related oscillation with period 52 days is found in the GALLEX data. The relationship of these periods to the rotational period inferred from neutrino data suggests that they are due to r-mode oscillations.