The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

In this article, we present the application of the weighted horizontal gradient of magnetic field (WGM) flare prediction method to 3-dimensional (3D) extrapolated magnetic configurations of 13 flaring solar active regions (ARs). The main aim is to identify an optimal height range, if any, in the interface region between the photosphere and lower corona, where the flare onset time prediction capability of WGM is best exploited. The optimal height is where flare prediction, by means of the WGM method, is achieved earlier than at the photospheric level. 3D magnetic structures, based on potential and non-linear force-free field extrapolations, are constructed to study a vertical range from the photosphere up to the low corona with a 45 km step size. The WGM method is applied as a function of height to all 13 flaring AR cases that are subject to certain selection criteria. We found that applying the WGM method between 1000 and 1800 km above the solar surface would improve the prediction of the flare onset time by around 2-8 hours.Certain caveats and an outlook for future work along these lines are also discussed.

A crucial challenge to successful flare prediction is forecasting periods that transition between "flare-quiet" and "flare-active". Building on earlier studies in this series (Barnes et al. 2016; Leka et al. 2019a,b) in which we describe methodology, details, and results of flare forecasting comparison efforts, we focus here on patterns of forecast outcomes (success and failure) over multi-day periods. A novel analysis is developed to evaluate forecasting success in the context of catching the first event of flare-active periods, and conversely, of correctly predicting declining flare activity. We demonstrate these evaluation methods graphically and quantitatively as they provide both quick comparative evaluations and options for detailed analysis. For the testing interval 2016-2017, we determine the relative frequency distribution of two-day dichotomous forecast outcomes for three different event histories (i.e., event/event, no-event/event and event/no-event), and use it to highlight performance differences between forecasting methods. A trend is identified across all forecasting methods that a high/low forecast probability on day-1 remains high/low on day-2 even though flaring activity is transitioning. For M-class and larger flares, we find that explicitly including persistence or prior flare history in computing forecasts helps to improve overall forecast performance. It is also found that using magnetic/modern data leads to improvement in catching the first-event/first-no-event transitions. Finally, 15% of major (i.e., M-class or above) flare days over the testing interval were effectively missed due to a lack of observations from instruments away from the Earth-Sun line.

In analyses of rare-events, regardless of the domain of application, class-imbalance issue is intrinsic. Although the challenges are known to data experts, their explicit impact on the analytic and the decisions made based on the findings are often overlooked. This is in particular prevalent in interdisciplinary research where the theoretical aspects are sometimes overshadowed by the challenges of the application. To show-case these undesirable impacts, we conduct a series of experiments on a recently created benchmark data, named Space Weather ANalytics for Solar Flares (SWAN-SF). This is a multivariate time series dataset of magnetic parameters of active regions. As a remedy for the imbalance issue, we study the impact of data manipulation (undersampling and oversampling) and model manipulation (using class weights). Furthermore, we bring to focus the auto-correlation of time series that is inherited from the use of sliding window for monitoring flares' history. Temporal coherence, as we call this phenomenon, invalidates the randomness assumption, thus impacting all sampling practices including different cross-validation techniques. We illustrate how failing to notice this concept could give an artificial boost in the forecast performance and result in misleading findings. Throughout this study we utilized Support Vector Machine as a classifier, and True Skill Statistics as a verification metric for comparison of experiments. We conclude our work by specifying the correct practice in each case, and we hope that this study could benefit researchers in other domains where time series of rare events are of interest.

We investigate the relation between characteristics of coronal mass ejections and parameterizations of the eruptive capability of solar active regions widely used in solar flare prediction schemes. These parameters, some of which are explored for the first time, are properties related to topological features, namely, magnetic polarity inversion lines (MPILs) that indicate large amounts of stored non-potential (i.e. free) magnetic energy. We utilize the Space Weather Database of Notifications, Knowledge, Information (DONKI) and the Large Angle and Spectrometric Coronograph (LASCO) databases to find flare-associated coronal mass ejections and their kinematic characteristics while properties of MPILs are extracted from Helioseismic and Magnetic Imager (HMI) vector magnetic-field observations of active regions to extract the properties of source-region MPILs. The correlation between all properties and the characteristics of CMEs ranges from moderate to very strong. More significant correlations hold particularly for fast CMEs, which are most important in terms of adverse space-weather manifestations. Non-neutralized currents and the length of the main MPIL exhibit significantly stronger correlations than the rest of the properties. This finding supports a causal relationship between coronal mass ejections and non-neutralized electric currents in highly sheared, conspicuous MPILs. In addition, non-neutralized currents and MPIL length carry distinct, independent information as to the eruptive potential of active regions. The combined total amount of non-neutralized electric currents and the length of the main polarity inversion line, therefore, reflect more efficiently than other parameters the eruptive capacity of solar active regions and the CME kinematic characteristics stemming from these regions.

Solar flares originate from magnetically active regions but not all solar active regions give rise to a flare. Therefore, the challenge of solar flare prediction benefits by an intelligent computational analysis of physics-based properties extracted from active region observables, most commonly line-of-sight or vector magnetograms of the active-region photosphere. For the purpose of flare forecasting, this study utilizes an unprecedented 171 flare-predictive active region properties, mainly inferred by the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) in the course of the European Union Horizon 2020 FLARECAST project. Using two different supervised machine learning methods that allow feature ranking as a function of predictive capability, we show that: i) an objective training and testing process is paramount for the performance of every supervised machine learning method; ii) most properties include overlapping information and are therefore highly redundant for flare prediction; iii) solar flare prediction is still - and will likely remain - a predominantly probabilistic challenge.

A workshop was recently held at Nagoya University (31 October - 02 November 2017), sponsored by the Center for International Collaborative Research, at the Institute for Space-Earth Environmental Research, Nagoya University, Japan, to quantitatively compare the performance of today's operational solar flare forecasting facilities. Building upon Paper I of this series (Barnes et al. 2016), in Paper II (Leka et al. 2019) we described the participating methods for this latest comparison effort, the evaluation methodology, and presented quantitative comparisons. In this paper we focus on the behavior and performance of the methods when evaluated in the context of broad implementation differences. Acknowledging the short testing interval available and the small number of methods available, we do find that forecast performance: 1) appears to improve by including persistence or prior flare activity, region evolution, and a human "forecaster in the loop"; 2) is hurt by restricting data to disk-center observations; 3) may benefit from long-term statistics, but mostly when then combined with modern data sources and statistical approaches. These trends are arguably weak and must be viewed with numerous caveats, as discussed both here and in Paper II. Following this present work, we present in Paper IV a novel analysis method to evaluate temporal patterns of forecasting errors of both types (i.e., misses and false alarms; Park et al. 2019). Hence, most importantly, with this series of papers we demonstrate the techniques for facilitating comparisons in the interest of establishing performance-positive methodologies.

Solar flares are extremely energetic phenomena in our Solar System. Their impulsive, often drastic radiative increases, in particular at short wavelengths, bring immediate impacts that motivate solar physics and space weather research to understand solar flares to the point of being able to forecast them. As data and algorithms improve dramatically, questions must be asked concerning how well the forecasting performs; crucially, we must ask how to rigorously measure performance in order to critically gauge any improvements. Building upon earlier-developed methodology (Barnes et al, 2016, Paper I), international representatives of regional warning centers and research facilities assembled in 2017 at the Institute for Space-Earth Environmental Research, Nagoya University, Japan to - for the first time - directly compare the performance of operational solar flare forecasting methods. Multiple quantitative evaluation metrics are employed, with focus and discussion on evaluation methodologies given the restrictions of operational forecasting. Numerous methods performed consistently above the "no skill" level, although which method scored top marks is decisively a function of flare event definition and the metric used; there was no single winner. Following in this paper series we ask why the performances differ by examining implementation details (Leka et al. 2019, Paper III), and then we present a novel analysis method to evaluate temporal patterns of forecasting errors in (Park et al. 2019, Paper IV). With these works, this team presents a well-defined and robust methodology for evaluating solar flare forecasting methods in both research and operational frameworks, and today's performance benchmarks against which improvements and new methods may be compared.

Solar active regions (ARs) that produce major flares typically exhibit strong plasma shear flows around photospheric magnetic polarity inversion lines (MPILs). It is therefore important to quantitatively measure such photospheric shear flows in ARs for a better understanding of their relation to flare occurrence. Photospheric flow fields were determined by applying the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM) method to a large data set of 2,548 co-aligned pairs of AR vector magnetograms with 12-min separation over the period 2012-2016. From each AR flow-field map, three shear-flow parameters were derived corresponding to the mean (<S>), maximum (S_max) and integral (S_sum) shear-flow speeds along strong-gradient, strong-field MPIL segments. We calculated flaring rates within 24 hr as a function of each shear-flow parameter, and also investigated the relation between the parameters and the waiting time ({\tau}) until the next major flare (class M1.0 or above) after the parameter observation. In general, it is found that the larger S_sum an AR has, the more likely it is for the AR to produce flares within 24 hr. It is also found that among ARs which produce major flares, if one has a larger value of S_sum then {\tau} generally gets shorter. These results suggest that large ARs with widespread and/or strong shear flows along MPILs tend to not only be more flare productive, but also produce major flares within 24 hr or less.

Efficient prediction of solar flares relies on parameters that quantify the eruptive capability of solar active regions. Several such quantitative predictors have been proposed in the literature, inferred mostly from photospheric magnetograms and/or white-light observations. Two of them are the Ising energy and the sum of the total horizontal magnetic field gradient. The former has been developed from line-of-sight magnetograms, while the latter uses sunspot detections and characteristics, based on continuum images. Aiming to include these parameters in an automated prediction scheme, we test their applicability on regular photospheric magnetic field observations provided by the Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO). We test their efficiency as predictors of flaring activity on a representative sample of active regions and investigate possible modifications of these quantities. The Ising energy appears to be an efficient predictor, and the efficiency is even improved if it is modified to describe interacting magnetic partitions or sunspot umbrae. The sum of the horizontal magnetic field gradient appears to be slightly more promising than the three variations of the Ising energy we implement in this article. The new predictors are also compared with two very promising predictors: the effective connected magnetic field strength and the total unsigned non-neutralized current. Our analysis shows that the efficiency of morphological predictors depends on projection effects in a nontrivial way. All four new predictors are found useful for inclusion in an automated flare forecasting facility, such as the Flare Likelihood and Region Eruption Forecasting (FLARECAST), but their utility, among others, will ultimately be determined by the validation effort underway in the framework of the FLARECAST project.

We propose a forecasting approach for solar flares based on data from Solar Cycle 24, taken by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) mission. In particular, we use the Space-weather HMI Active Region Patches (SHARP) product that facilitates cut-out magnetograms of solar active regions (AR) in the Sun in near-real-time (NRT), taken over a five-year interval (2012 - 2016). Our approach utilizes a set of thirteen predictors, which are not included in the SHARP metadata, extracted from line-of-sight and vector photospheric magnetograms. We exploit several Machine Learning (ML) and Conventional Statistics techniques to predict flares of peak magnitude >M1 and >C1, within a 24 h forecast window. The ML methods used are multi-layer perceptrons (MLP), support vector machines (SVM) and random forests (RF). We conclude that random forests could be the prediction technique of choice for our sample, with the second best method being multi-layer perceptrons, subject to an entropy objective function. A Monte Carlo simulation showed that the best performing method gives accuracy ACC=0.93(0.00), true skill statistic TSS=0.74(0.02) and Heidke skill score HSS=0.49(0.01) for >M1 flare prediction with probability threshold 15% and ACC=0.84(0.00), TSS=0.60(0.01) and HSS=0.59(0.01) for >C1 flare prediction with probability threshold 35%.

We apply a novel pre-flare tracking of sunspot groups towards improving the estimation of flare onset time by focusing on the evolution of the 3D magnetic field construction of AR 11429. The 3D magnetic structure is based on potential field extrapolation encompassing a vertical range from the photosphere through the chromosphere and transition region into the low corona. The basis of our proxy measure of activity prediction is the so-called weighted horizontal gradient of magnetic field (WG_M) defined between spots of opposite polarities close to the polarity inversion line of an active region. The temporal variation of the distance of the barycenter of the opposite polarities is also found to possess potentially important diagnostic information about the flare onset time estimation as function of height similar to its counterpart introduced initially in an application at the photosphere only in Korsos et al. (2015). We apply the photospheric pre-flare behavioural patterns of sunspot groups to the evolution of their associated 3D-constructed AR 11429 as function of height. We found that at a certain height in the lower solar atmosphere the onset time may be estimated much earlier than at the photosphere or at any other heights. Therefore, we present a tool and recipe that may potentially identify the optimum height for flare prognostic in the solar atmosphere allowing to improve our flare prediction capability and capacity.

The effect of using two representations of the normal-to-surface magnetic field to calculate photospheric measures that are related to active region (AR) potential for flaring is presented. Several AR properties were computed using line-of-sight ($B_{\rm los}$) and spherical-radial ($B_{r}$) magnetograms from the Space-weather HMI Active Region Patch (SHARP) products of the \emph{Solar Dynamics Observatory}, characterizing the presence and features of magnetic polarity inversion lines, fractality, and magnetic connectivity of the AR photospheric field. The data analyzed corresponds to $\approx$ 4,000 AR observations, achieved by randomly selecting 25 % of days between September 2012 and May 2016 for analysis at 6-hr cadence. Results from this statistical study include: i) the $B_{r}$ component results in a slight upwards shift of property values in a manner consistent with a field-strength underestimation by the $B_{\rm los}$ component; ii) using the $B_{r}$ component results in significantly lower inter-property correlation in one-third of the cases, implying more independent information about the state of the AR photospheric magnetic field; iii) flaring rates for each property vary between the field components in a manner consistent with the differences in property-value ranges resulting from the components; iv) flaring rates generally increase for higher values of properties, except Fourier spectral power index that has flare rates peaking around a value of $5/3$. These findings indicate that there may be advantages in using $B_{r}$ rather than $B_{\rm los}$ in calculating flare-related AR magnetic properties, especially for regions located far from central meridian.

We explore the association of non-neutralized currents with solar flare occurrence in a sizable sample of observations, aiming to show the potential of such currents in solar flare prediction. We use the regularly produced high quality vector magnetograms by the Helioseismic Magnetic Imager and more specifically the Space weather HMI Active Region Patches (SHARP). Through a newly established method, that incorporates detailed error analysis, we calculate the non-neutralized currents contained in active regions (AR). Two predictors are produced, namely the total and the maximum unsigned non-neutralized current. Both are tested in AR time-series and a representative sample of point-in-time observations during the interval 2012-2016. The average values of non-neutralized currents in flaring active regions are by more than an order of magnitude higher than in non-flaring regions and correlate very well with the corresponding flare index. The temporal evolution of these parameters appears to be connected to physical processes, such as flux emergence and/or magnetic polarity inversion line formation that are associated with increased solar flare activity. Using Bayesian inference of flaring probabilities, it is shown that the total unsigned non-neutralized current outperforms significantly the total unsigned magnetic flux and other well established current-related predictors. Thus, it shows good prospects for inclusion in an operational flare forecasting service. We plan to use the new predictor in the framework of the FLARECAST project along with other highly performing predictors.

Patsourakos et al. (Astrophys. J. 817, 14, 2016) and Patsourakos and Georgoulis (Astron. Astrophys. 595, A121, 2016) introduced a method to infer the axial magnetic field in flux-rope coronal mass ejections (CMEs) in the solar corona and farther away in the interplanetary medium. The method, based on the conservation principle of magnetic helicity, uses the relative magnetic helicity of the solar source region as input estimates, along with the radius and length of the corresponding CME flux rope. The method was initially applied to cylindrical force-free flux ropes, with encouraging results. We hereby extend our framework along two distinct lines. First, we generalize our formalism to several possible flux-rope configurations (linear and nonlinear force-free, non-force-free, spheromak, and torus) to investigate the dependence of the resulting CME axial magnetic field on input parameters and the employed flux-rope configuration. Second, we generalize our framework to both Sun-like and active M-dwarf stars hosting superflares. In a qualitative sense, we find that Earth may not experience severe atmosphere-eroding magnetospheric compression even for eruptive solar superflares with energies ~ 10^4 times higher than those of the largest Geostationary Operational Environmental Satellite (GOES) X-class flares currently observed. In addition, the two recently discovered exoplanets with the highest Earth-similarity index, Kepler 438b and Proxima b, seem to lie in the prohibitive zone of atmospheric erosion due to interplanetary CMEs (ICMEs), except when they possess planetary magnetic fields that are much higher than that of Earth.

We study the writhe, twist and magnetic helicity of different magnetic flux ropes, based on models of the solar coronal magnetic field structure. These include an analytical force-free Titov--D\'emoulin equilibrium solution, non force-free magnetohydrodynamic simulations, and nonlinear force-free magnetic field models. The geometrical boundary of the magnetic flux rope is determined by the quasi-separatrix layer and the bottom surface, and the axis curve of the flux rope is determined by its overall orientation. The twist is computed by the Berger--Prior formula that is suitable for arbitrary geometry and both force-free and non-force-free models. The magnetic helicity is estimated by the twist multiplied by the square of the axial magnetic flux. We compare the obtained values with those derived by a finite volume helicity estimation method. We find that the magnetic helicity obtained with the twist method agrees with the helicity carried by the purely current-carrying part of the field within uncertainties for most test cases. It is also found that the current-carrying part of the model field is relatively significant at the very location of the magnetic flux rope. This qualitatively explains the agreement between the magnetic helicity computed by the twist method and the helicity contributed purely by the current-carrying magnetic field.

Magnetic helicity is a conserved quantity of ideal magneto-hydrodynamics characterized by an inverse turbulent cascade. Accordingly, it is often invoked as one of the basic physical quantities driving the generation and structuring of magnetic fields in a variety of astrophysical and laboratory plasmas. We provide here the first systematic comparison of six existing methods for the estimation of the helicity of magnetic fields known in a finite volume. All such methods are reviewed, benchmarked, and compared with each other, and specifically tested for accuracy and sensitivity to errors. To that purpose, we consider four groups of numerical tests, ranging from solutions of the three-dimensional, force-free equilibrium, to magneto-hydrodynamical numerical simulations. Almost all methods are found to produce the same value of magnetic helicity within few percent in all tests. In the more solar-relevant and realistic of the tests employed here, the simulation of an eruptive flux rope, the spread in the computed values obtained by all but one method is only 3%, indicating the reliability and mutual consistency of such methods in appropriate parameter ranges. However, methods show differences in the sensitivity to numerical resolution and to errors in the solenoidal property of the input fields. In addition to finite volume methods, we also briefly discuss a method that estimates helicity from the field lines' twist, and one that exploits the field's value at one boundary and a coronal minimal connectivity instead of a pre-defined three-dimensional magnetic-field solution.

Coronal heating through the explosive release of magnetic energy remains an open problem in solar physics. Recent hydrodynamical models attempt an investigation by placing swarms of 'nanoflares' at random sites and times in modeled one-dimensional coronal loops. We investigate the problem in three dimensions, using extrapolated coronal magnetic fields of observed solar active regions. We apply a nonlinear force-free field extrapolation above an observed photospheric magnetogram of NOAA active region (AR) 11158. We then determine the locations, energy contents, and volumes of 'unstable' areas, namely areas prone to releasing magnetic energy due to locally accumulated electric current density. Statistical distributions of these volumes and their fractal dimension are inferred, investigating also their dependence on spatial resolution. Further adopting a simple resistivity model, we infer the properties of the fractally distributed electric fields in these volumes. Next, we monitor the evolution of 10^5 particles (electrons and ions) obeying an initial Maxwellian distribution with a temperature of 10 eV, by following their trajectories and energization when subjected to the resulting electric fields. For computational convenience, the length element of the magnetic-field extrapolation is 1 arcsec, much coarser than the particles collisional mean free path in the low corona. The presence of collisions traps the bulk of the plasma around the unstable volumes, or current sheets (UCS), with only a tail of the distribution gaining substantial energy. Assuming that the distance between UCS is similar to the collisional mean free path we find that the low active-region corona is heated to 100-200 eV, corresponding to temperatures exceeding 2 MK, within tens of seconds for electrons and thousands of seconds for ions. Fractally distributed, nanoflare-triggening fragmented UCS ...

Aims. The magnetic field of coronal mass ejections (CMEs) determines their structure, evolution, and energetics, as well as their geoeffectiveness. However, we currently lack routine diagnostics of the near-Sun CME magnetic field, which is crucial for determining the subsequent evolution of CMEs. Methods. We recently presented a method to infer the near-Sun magnetic field magnitude of CMEs and then extrapolate it to 1 AU. This method uses relatively easy to deduce observational estimates of the magnetic helicity in CME-source regions along with geometrical CME fits enabled by coronagraph observations. We hereby perform a parametric study of this method aiming to assess its robustness. We use statistics of active region (AR) helicities and CME geometrical parameters to determine a matrix of plausible near-Sun CME magnetic field magnitudes. In addition, we extrapolate this matrix to 1 AU and determine the anticipated range of CME magnetic fields at 1 AU representing the radial falloff of the magnetic field in the CME out to interplanetary (IP) space by a power law with index aB. Results. The resulting distribution of the near-Sun (at 10 Rs ) CME magnetic fields varies in the range [0.004, 0.02] G, comparable to, or higher than, a few existing observational inferences of the magnetic field in the quiescent corona at the same distance. We also find that a theoretically and observationally motivated range exists around aB = -1.6 +-0.2, thereby leading to a ballpark agreement between our estimates and observationally inferred field magnitudes of magnetic clouds (MCs) at L1. Conclusions. In a statistical sense, our method provides results that are consistent with observations.

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.

We employ different shapes of apodizing windows in the local correlation tracking (LCT) routine to retrieve horizontal velocities using numerical simulations of convection. LCT was applied on a time sequence of temperature maps generated by the Nirvana code with four different apodizing windows, namely--Gaussian, Lorentzian, trapezoidal and triangular, with varying widths. In terms of correlations (between the LCT-retrieved and simulated flow field), the triangular and the trapezoidal perform the best and worst, respectively. On segregating the intrinsic velocities in the simulations on the basis of their magnitudes, we find that for all windows, a significantly higher correlation is obtained for the intermediate and high-velocity bins and only modest or weak values in the low-velocity bins. The differences between the LCT-retrieved and simulated flow fields were determined spatially which show large residuals at or close to the boundary of granules. The extent to which the horizontal flow vectors retrieved by LCT compare with the simulated values, depends entirely on the width of the central peak of the apodizing window for a given $\sigma$. Even though LCT suffers from a lack of spatial content as seen in simulations, its simplicity and speed could serve as a viable first-order tool to probe horizontal flows--one that is ideal for large data sets.

We assess the validity of the free magnetic energy - relative magnetic helicity diagram for solar magnetic structures. We used two different methods of calculating the free magnetic energy and the relative magnetic helicity budgets: a classical, volume-calculation nonlinear force-free (NLFF) method applied to finite coronal magnetic structures and a surface-calculation NLFF derivation that relies on a single photospheric or chromospheric vector magnetogram. Both methods were applied to two different data sets, namely synthetic active-region cases obtained by three-dimensional magneto-hydrodynamic (MHD) simulations and observed active-region cases, which include both eruptive and noneruptive magnetic structures. The derived energy--helicity diagram shows a consistent monotonic scaling between relative helicity and free energy with a scaling index 0.84$\pm$0.05 for both data sets and calculation methods. It also confirms the segregation between noneruptive and eruptive active regions and the existence of thresholds in both free energy and relative helicity for active regions to enter eruptive territory. We consider the previously reported energy-helicity diagram of solar magnetic structures as adequately validated and envision a significant role of the uncovered scaling in future studies of solar magnetism.

In earlier works we introduced and tested a nonlinear force-free (NLFF) method designed to self-consistently calculate the free magnetic energy and the relative magnetic helicity budgets of the corona of observed solar magnetic structures. The method requires, in principle, only a single, photospheric or low-chromospheric, vector magnetogram of a quiet-Sun patch or an active region and performs calculations in the absence of three-dimensional magnetic and velocity-field information. In this work we strictly validate this method using three-dimensional coronal magnetic fields. Benchmarking employs both synthetic, three-dimensional magnetohydrodynamic simulations and nonlinear force-free field extrapolations of the active-region solar corona. We find that our time-efficient NLFF method provides budgets that differ from those of more demanding semi-analytical methods by a factor of ~3, at most. This difference is expected from the physical concept and the construction of the method. Temporal correlations show more discrepancies that, however, are soundly improved for more complex, massive active regions, reaching correlation coefficients of the order of, or exceeding, 0.9. In conclusion, we argue that our NLFF method can be reliably used for a routine and fast calculation of free magnetic energy and relative magnetic helicity budgets in targeted parts of the solar magnetized corona. As explained here and in previous works, this is an asset that can lead to valuable insight into the physics and the triggering of solar eruptions.

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 investigate the free magnetic energy and relative magnetic helicity budgets of solar quiet regions. Using a novel non-linear force-free method requiring single solar vector magnetograms we calculate the instantaneous free magnetic energy and relative magnetic helicity budgets in 55 quiet-Sun vector magnetograms. As in a previous work on active regions, we construct here for the first time the (free) energy-(relative) helicity diagram of quiet-Sun regions. We find that quiet-Sun regions have no dominant sense of helicity and show monotonic correlations a) between free magnetic energy/relative helicity and magnetic network area and, consequently, b) between free magnetic energy and helicity. Free magnetic energy budgets of quiet-Sun regions represent a rather continuous extension of respective active-region budgets towards lower values, but the corresponding helicity transition is discontinuous due to the incoherence of the helicity sense contrary to active regions. We further estimate the instantaneous free magnetic-energy and relative magnetic-helicity budgets of the entire quiet Sun, as well as the respective budgets over an entire solar cycle. Derived instantaneous free magnetic energy budgets and, to a lesser extent, relative magnetic helicity budgets over the entire quiet Sun are comparable to the respective budgets of a sizeable active region, while total budgets within a solar cycle are found higher than previously reported. Free-energy budgets are comparable to the energy needed to power fine-scale structures residing at the network, such as mottles and spicules.

In previous works we introduced a nonlinear force-free method that self-consistently calculates the instantaneous budgets of free magnetic energy and relative magnetic helicity in solar active regions (ARs). Calculation is expedient and practical, using only a single vector magnetogram per computation. We apply this method to a timeseries of 600 high-cadence vector magnetograms of the eruptive NOAA AR 11158 acquired by the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory over a five-day observing interval. Besides testing our method extensively, we use it to interpret the dynamical evolution in the AR, including eruptions. We find that the AR builds large budgets of both free magnetic energy and relative magnetic helicity, sufficient to power many more eruptions than the ones it gave within the interval of interest. For each of these major eruptions, we find eruption-related decreases and subsequent free-energy and helicity budgets that are consistent with the observed eruption (flare and coronal-mass-ejection [CME]) sizes. In addition, we find that (1) evolution in the AR is consistent with the recently proposed (free) energy - (relative) helicity diagram of solar ARs, (2) eruption-related decreases occur {\it before} the flare and the projected CME-launch times, suggesting that CME progenitors precede flares, and (3) self-terms of free energy and relative helicity most likely originate from respective mutual-terms, following a progressive mutual-to-self conversion pattern that most likely stems from magnetic reconnection. This results in the non-ideal formation of increasingly helical pre-eruption structures and instigates further research on the triggering of solar eruptions with magnetic helicity firmly placed in the eruption cadre.

We model nanoflare heating of extrapolated active-region coronal loops via the acceleration of electrons and protons in Harris-type current sheets. The kinetic energy of the accelerated particles is estimated using semi-analytical and test-particle-tracing approaches. Vector magnetograms and photospheric Doppler velocity maps of NOAA active region 09114, recorded by the Imaging Vector Magnetograph (IVM), were used for this analysis. A current-free field extrapolation of the active-region corona was first constructed. The corresponding Poynting fluxes at the footpoints of 5000 extrapolated coronal loops were then calculated. Assuming that reconnecting current sheets develop along these loops, we utilized previous results to estimate the kinetic-energy gain of the accelerated particles and we related this energy to nanoflare heating and macroscopic loop characteristics. Kinetic energies of 0.1 to 8 keV (for electrons) and 0.3 to 470 keV (for protons) were found to cause heating rates ranging from $10^{-6}$ to 1 $\mathrm{erg\, s^{-1} cm^{-3}}$. Hydrodynamic simulations show that such heating rates can sustain plasma in coronal conditions inside the loops and generate plasma thermal distributions which are consistent with active region observations. We concluded the analysis by computing the form of X-ray spectra generated by the accelerated electrons using the thick target approach that were found to be in agreement with observed X-ray spectra, thus supporting the plausibility of our nanoflare-heating scenario.

Using solar vector magnetograms of the highest available spatial resolution and signal-to-noise ratio we perform a detailed study of electric current patterns in two solar active regions: a flaring/eruptive, and a flare-quiet one. We aim to determine whether active regions inject non-neutralized (net) electric currents in the solar atmosphere, responding to a debate initiated nearly two decades ago that remains inconclusive. We find that well-formed, intense magnetic polarity inversion lines (PILs) within active regions are the only photospheric magnetic structures that support significant net current. More intense PILs seem to imply stronger non-neutralized current patterns per polarity. This finding revises previous works that claim frequent injections of intense non-neutralized currents by most active regions appearing in the solar disk but also works that altogether rule out injection of non-neutralized currents. In agreement with previous studies, we also find that magnetically isolated active regions remain globally current-balanced. In addition, we confirm and quantify the preference of a given magnetic polarity to follow a given sense of electric currents, indicating a dominant sense of twist in active regions. This coherence effect is more pronounced in more compact active regions with stronger PILs and must be of sub-photospheric origin. Our results yield a natural explanation of the Lorentz force, invariably generating velocity and magnetic shear along strong PILs, thus setting a physical context for the observed pre-eruption evolution in solar active regions.

Using a recently proposed nonlinear force-free method designed for single vector magnetograms of solar active regions we calculate the instantaneous free magnetic energy and relative magnetic helicity budgets in 162 vector magnetograms corresponding to 42 different active regions. We find a statistically robust, monotonic correlation between the free magnetic energy and the relative magnetic helicity in the studied regions. This correlation implies that magnetic helicity, besides free magnetic energy, may be an essential ingredient for major solar eruptions. Eruptive active regions appear well segregated from non-eruptive ones in both free energy and relative helicity with major (at least M-class) flares occurring in active regions with free energy and relative helicity exceeding 4x10^{31} erg and 2x10^{42} Mx^2, respectively. The helicity threshold agrees well with estimates of helicity contents of typical coronal mass ejections.

Expanding on an earlier work that relied on linear force-free magnetic fields, we self-consistently derive the instantaneous free magnetic energy and relative magnetic helicity budgets of an unknown three-dimensional nonlinear force-free magnetic structure extending above a single, known lower-boundary magnetic field vector. The proposed method does not rely on the detailed knowledge of the three-dimensional field configuration but is general enough to employ only a magnetic connectivity matrix on the lower boundary. The calculation yields a minimum free magnetic energy and a relative magnetic helicity consistent with this free magnetic energy. The method is directly applicable to photospheric or chromospheric vector magnetograms of solar active regions. Upon validation, it basically reproduces magnetic energies and helicities obtained by well-known, but computationally more intensive and non-unique, methods relying on the extrapolated three-dimensional magnetic field vector. We apply the method to three active regions, calculating the photospheric connectivity matrices by means of simulated annealing, rather than a model-dependent nonlinear force-free extrapolation. For two of these regions we correct for the inherent linear force-free overestimation in free energy and relative helicity that is larger for larger, more eruptive, active regions. In the third studied region, our calculation can lead to a physical interpretation of observed eruptive manifestations. We conclude that the proposed method, including the proposed inference of the magnetic connectivity matrix, is practical enough to contribute to a physical interpretation of the dynamical evolution of solar active regions.

We study plasma flows along selected coronal loops in NOAA Active Region 10926, observed on 3 December 2006 with Hinode's EUV Imaging Spectrograph (EIS). From the shape of the loops traced on intensity images and the Doppler shifts measured along their length we compute their three-dimensional (3D) shape and plasma flow velocity using a simple geometrical model. This calculation was performed for loops visible in the Fe VIII 185 Ang., Fe X 184 Ang., Fe XII 195 Ang., Fe XIII 202 Ang., and Fe XV 284 Ang. spectral lines. In most cases the flow is unidirectional from one footpoint to the other but there are also cases of draining motions from the top of the loops to their footpoints. Our results indicate that the same loop may show different flow patterns when observed in different spectral lines, suggesting a dynamically complex rather than a monolithic structure. We have also carried out magnetic extrapolations in the linear force-free field approximation using SOHO/MDI magnetograms, aiming toward a first-order identification of extrapolated magnetic field lines corresponding to the reconstructed loops. In all cases, the best-fit extrapolated lines exhibit left-handed twist (alpha < 0), in agreement with the dominant twist of the region.

In a recent paper, Leka at al. (Solar Phys. 260, 83, 2009)constructed a synthetic vector magnetogram representing a three-dimensional magnetic structure defined only within a fraction of an arcsec in height. They rebinned the magnetogram to simulate conditions of limited spatial resolution and then compared the results of various azimuth disambiguation methods on the resampled data. Methods relying on the physical calculation of potential and/or non-potential magnetic fields failed in nearly the same, extended parts of the field of view and Leka et al. (2009) attributed these failures to the limited spatial resolution. This study shows that the failure of these methods is not due to the limited spatial resolution but due to the narrowly defined test data. Such narrow magnetic structures are not realistic in the real Sun. Physics-based disambiguation methods, adapted for solar magnetic fields extending to infinity, are not designed to handle such data; hence, they could only fail this test. I demonstrate how an appropriate limited-resolution disambiguation test can be performed by constructing a synthetic vector magnetogram very similar to that of Leka et al. (2009) but representing a structure defined in the semi-infinite space above the solar photosphere. For this magnetogram I find that even a simple potential-field disambiguation method manages to resolve the ambiguity very successfully, regardless of limited spatial resolution. Therefore, despite the conclusions of Leka et al. (2009), a proper limited-spatial-resolution test of azimuth disambiguation methods is yet to be performed in order to identify the best ideas and algorithms.

We interpret solar flares as events originating from active regions that have reached the Self Organized Critical state, by using a refined Cellular Automaton model with initial conditions derived from observations. Aims: We investigate whether the system, with its imposed physical elements,reaches a Self Organized Critical state and whether well-known statistical properties of flares, such as scaling laws observed in the distribution functions of characteristic parameters, are reproduced after this state has been reached. Results: Our results show that Self Organized Criticality is indeed reached when applying specific loading and relaxation rules. Power law indices obtained from the distribution functions of the modeled flaring events are in good agreement with observations. Single power laws (peak and total flare energy) as well as power laws with exponential cutoff and double power laws (flare duration) are obtained. The results are also compared with observational X-ray data from GOES satellite for our active-region sample. Conclusions: We conclude that well-known statistical properties of flares are reproduced after the system has reached Self Organized Criticality. A significant enhancement of our refined Cellular Automaton model is that it commences the simulation from observed vector magnetograms, thus facilitating energy calculation in physical units. The model described in this study remains consistent with fundamental physical requirements, and imposes physically meaningful driving and redistribution rules.

Multiple recent investigations of solar magnetic field measurements have raised claims that the scale-free (fractal) or multiscale (multifractal) parameters inferred from the studied magnetograms may help assess the eruptive potential of solar active regions, or may even help predict major flaring activity stemming from these regions. We investigate these claims here, by testing three widely used scale-free and multiscale parameters, namely, the fractal dimension, the multifractal structure function and its inertial-range exponent, and the turbulent power spectrum and its power-law index, on a comprehensive data set of 370 timeseries of active-region magnetograms (17,733 magnetograms in total) observed by SOHO's Michelson Doppler Imager (MDI) over the entire Solar Cycle 23. We find that both flaring and non-flaring active regions exhibit significant fractality, multifractality, and non-Kolmogorov turbulence but none of the three tested parameters manages to distinguish active regions with major flares from flare-quiet ones. We also find that the multiscale parameters, but not the scale-free fractal dimension, depend sensitively on the spatial resolution and perhaps the observational characteristics of the studied magnetograms. Extending previous works, we attribute the flare-forecasting inability of fractal and multifractal parameters to i) a widespread multiscale complexity caused by a possible underlying self-organization in turbulent solar magnetic structures, flaring and non-flaring alike, and ii) a lack of correlation between the fractal properties of the photosphere and overlying layers, where solar eruptions occur. However useful for understanding solar magnetism, therefore, scale-free and multiscale measures may not be optimal tools for active-region characterization in terms of eruptive ability or, ultimately,for major solar-flare prediction.

We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday's equation, give rise to a respective normal field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modelled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition - the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.

Observations from the Hinode X-ray telescope (XRT) are used to study the structure of X-ray bright points (XBPs), sources of coronal jets. Several jet events are found to erupt from S-shaped bright points, suggesting that coronal micro-sigmoids are progenitors of the jets. The observations may help to explain numerous characteristics of coronal jets, such as helical structures and shapes. They also suggest that solar activity may be self-similar within a wide range of scales in terms of both properties and evolution of the observed coronal structures.

Observations from the X-ray telescope (XRT) on Hinode are used to study the nature of X-ray bright points, sources of coronal jets. Several jet events in the coronal holes are found to erupt from small-scale, S-shaped bright regions. This finding suggests that coronal micro-sigmoids may well be progenitors of coronal jets. Moreover, the presence of these structures may explain numerous observed characteristics of jets such as helical structures, apparent transverse motions, and shapes. In analogy to large-scale sigmoids giving rise to coronal mass ejections (CMEs), a promising future task would perhaps be to investigate whether solar eruptive activity, from coronal jets to CMEs, is self-similar in terms of properties and instability mechanisms.

This work examines the relation between the fractal properties of the photospheric magnetic patterns and those of the coronal magnetic fields in solar active regions.Our results show no correlation between the fractal dimensions of the observed 2d photospheric structures and the extrapolated unstable volumes in the corona, when nonlinear force-free extrapolation is used. This result is independent of efforts to (1) bring the photospheric magnetic fields closer to a nonlinear force-free equilibrium and (2) omit the lower part of the modeled magnetic field volume that is almost completely filled by unstable volumes. A significant correlation between the fractal dimensions of the photospheric and coronal magnetic features is only observed at the zero level (lower limit) of approximation of a current-free (potential) magnetic field extrapolation. We conclude that the complicated transition from photospheric non-force-free fields to coronal force-free ones hampers any direct correlation between the fractal dimensions of the 2d photospheric patterns and their 3d counterparts in the corona at the nonlinear force-free limit.

The Vector Spectromagnetograph (VSM) instrument has been recording photospheric and chromospheric magnetograms daily since August 2003. Full-disk photospheric vector magnetograms are observed at least weekly and, since November 2006, area-scans of active regions daily. Quick-look vector magnetic images, plus X3D and FITS formated files, are now publicly available daily. In the near future, Milne-Eddington inversion parameter data will also be available and a typical observing day will include three full-disk photospheric vector magnetograms. Besides full-disk observations, the VSM is capable of high temporal cadence area-scans of both the photosphere and chromosphere. Carrington rotation and daily synoptic maps are also available from the photospheric magnetograms and coronal hole estimate images.

Using an efficient magnetic complexity index in the active-region solar photosphere, we quantify the preflare strength of the photospheric magnetic polarity inversion lines in 23 eruptive active regions with flare/CME/ICME events tracked all the way from the Sun to the Earth. We find that active regions with more intense polarity inversion lines host statistically stronger flares and faster, more impulsively accelerated, CMEs. No significant correlation is found between the strength of the inversion lines and the flare soft X-ray rise times, the ICME transit times, and the peak $Dst indices of the induced geomagnetic storms. Corroborating these and previous results, we speculate on a possible interpretation for the connection between source active regions, flares, and CMEs. Further work is needed to validate this concept and uncover its physical details.

The Vector Spectromagnetograph (VSM) of the NSO's Synoptic Optical Long-Term Investigations of the Sun (SOLIS) facility is now operational and obtains the first-ever vector magnetic field measurements of the entire visible solar hemisphere. To fully exploit the unprecedented SOLIS/VSM data, however, one must first address two critical problems: first, the study of solar active regions requires an automatic, physically intuitive, technique for active-region identification in the solar disk. Second, use of active-region vector magnetograms requires removal of the azimuthal $180^o$-ambiguity in the orientation of the transverse magnetic field component. Here we report on an effort to address both problems simultaneously and efficiently. To identify solar active regions we apply an algorithm designed to locate complex, flux-balanced, magnetic structures with a dominant E-W orientation on the disk. Each of the disk portions corresponding to active regions is thereafter extracted and subjected to the Nonpotential Magnetic Field Calculation (NPFC) method that provides a physically-intuitive solution of the 180-degree ambiguity. Both algorithms have been integrated into the VSM data pipeline and operate in real time, without human intervention. We conclude that this combined approach can contribute meaningfully to our emerging capability for full-disk vector magnetography as pioneered by SOLIS today and will be carried out by ground-based and space-borne magnetographs in the future.

We self-consistently derive the magnetic energy and relative magnetic helicity budgets of a three-dimensional linear force-free magnetic structure rooted in a lower boundary plane. For the potential magnetic energy we derive a general expression that gives results practically equivalent to those of the magnetic Virial theorem. All magnetic energy and helicity budgets are formulated in terms of surface integrals applied to the lower boundary, thus avoiding computationally intensive three-dimensional magnetic field extrapolations. We analytically and numerically connect our derivations with classical expressions for the magnetic energy and helicity, thus presenting a so-far lacking unified treatment of the energy/helicity budgets in the constant-alpha approximation. Applying our derivations to photospheric vector magnetograms of an eruptive and a noneruptive solar active regions, we find that the most profound quantitative difference between these regions lies in the estimated free magnetic energy and relative magnetic helicity budgets. If this result is verified with a large number of active regions, it will advance our understanding of solar eruptive phenomena. We also find that the constant-alpha approximation gives rise to large uncertainties in the calculation of the free magnetic energy and the relative magnetic helicity. Therefore, care must be exercised when this approximation is applied to photospheric magnetic field observations. Despite its shortcomings, the constant-alpha approximation is adopted here because this study will form the basis of a comprehensive nonlinear force-free description of the energetics and helicity in the active-region solar corona, which is our ultimate objective.

We outline a general methodology to infer the inductive velocity field vector in solar active regions. For the first time, both the field-aligned and the cross-field velocity components are reconstructed. The cross-field velocity solution accounts for the changes of the vertical magnetic field seen between a pair of successive active-region vector magnetograms via the ideal induction equation. The field-aligned velocity is obtained using the Doppler velocity and the calculated cross-field velocity. Solving the ideal induction equation in vector magnetograms measured at a given altitude in the solar atmosphere is an under-determined problem. In response, our general formalism allows the use of any additional constraint for the inductive cross-field velocity to enforce a unique solution in the induction equation. As a result, our methodology can give rise to new velocity solutions besides the one presented here. To constrain the induction equation, we use a special case of the minimum structure approximation that was introduced in previous studies and is already employed here to resolve the $180^o$-ambiguity in the input vector magnetograms. We reconstruct the inductive velocity for three active regions, including NOAA AR 8210 for which previous results exist. Our solution believably reproduces the horizontal flow patterns in the studied active regions but breaks down in cases of localized rapid magnetic flux emergence or submergence. Alternative approximations and constraints are possible and can be accommodated into our general formalism.

Intermittent magnetohydrodynamical turbulence is most likely at work in the magnetized solar atmosphere. As a result, an array of scaling and multi-scaling image-processing techniques can be used to measure the expected self-organization of solar magnetic fields. While these techniques advance our understanding of the physical system at work, it is unclear whether they can be used to predict solar eruptions, thus obtaining a practical significance for space weather. We address part of this problem by focusing on solar active regions and by investigating the usefulness of scaling and multi-scaling image-processing techniques in solar flare prediction. Since solar flares exhibit spatial and temporal intermittency, we suggest that they are the products of instabilities subject to a critical threshold in a turbulent magnetic configuration. The identification of this threshold in scaling and multi-scaling spectra would then contribute meaningfully to the prediction of solar flares. We find that the fractal dimension of solar magnetic fields and their multi-fractal spectrum of generalized correlation dimensions do not have significant predictive ability. The respective multi-fractal structure functions and their inertial-range scaling exponents, however, probably provide some statistical distinguishing features between flaring and non-flaring active regions. More importantly, the temporal evolution of the above scaling exponents in flaring active regions probably shows a distinct behavior starting a few hours prior to a flare and therefore this temporal behavior may be practically useful in flare prediction. The results of this study need to be validated by more comprehensive works over a large number of solar active regions.

The formation and evolution of active regions is an inherently complex phenomenon. Magnetic fields generated at the base of the convection zone follow a chaotic evolution before reaching the solar surface. In this article, we use a 2-D probabilistic Cellular Automaton (CA) to model the statistical properties of the magnetic patterns formed on the solar surface and to estimate the magnetic energy released in the interaction of opposite polarities. We assume that newly emerged magnetic flux tubes stimulate the emergence of new magnetic flux in their neighborhood. The flux-tubes move randomly on the surface of the sun, and they cancel and release their magnetic energy when they collide with magnetic flux of opposite polarity, or diffuse into the ``empty'' photosphere. We assume that cancellation of magnetic flux in collisions causes "flares" and determine the released energy as the difference in the square of the magnetic field flux. The statistics of the simulated "flares" follow a power-law distribution in energy, with power-law index a=2.2 +/- 0.1. The size distribution function of the simulated active regions exhibits a power law behavior with index k = 1.93 +/- 0.08, and the fractal dimension of the magnetized areas on the simulated solar surface is close to D = 1.42 +/- 0.12. Both quantities, D and k, are inside the range of the observed values.