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.

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.

We study the predictive capabilities of magnetic feature properties (MF) generated by Solar Monitor Active Region Tracker (SMART) for solar flare forecasting from two datasets: the full dataset of SMART detections from 1996 to 2010 that has been previously studied by Ahmed et al. (2011) and a subset of that dataset which only includes detections that are NOAA active regions (ARs). Main contributions: we use marginal relevance as a filter feature selection method to identify most useful SMART MF properties for separating flaring from non-flaring detections and logistic regression to derive classification rules to predict future observations. For comparison, we employ a Random Forest, Support Vector Machine and a set of Deep Neural Network models, as well as Lasso for feature selection. Using the linear model with three features we obtain significantly better results (TSS=0.84) to those reported by Ahmed et al.(2011) for the full dataset of SMART detections. The same model produced competitive results (TSS=0.67) for the dataset of SMART detections that are NOAA ARs which can be compared to a broader section of flare forecasting literature. We show that more complex models are not required for this data.

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.

Most solar flares originate in sunspot groups, where magnetic field changes lead to energy build-up and release. However, few flare-forecasting methods use information of sunspot-group evolution, instead focusing on static point-in-time observations. Here, a new forecast method is presented based upon the 24-hr evolution in McIntosh classification of sunspot groups. Evolution-dependent $\geqslant$C1.0 and $\geqslant$M1.0 flaring rates are found from NOAA-numbered sunspot groups over December 1988 to June 1996 (Solar Cycle 22; SC22) before converting to probabilities assuming Poisson statistics. These flaring probabilities are used to generate operational forecasts for sunspot groups over July 1996 to December 2008 (SC23), with performance studied by verification metrics. Major findings are: i) considering Brier skill score (BSS) for $\geqslant$C1.0 flares, the evolution-dependent McIntosh-Poisson method ($\text{BSS}_{\text{evolution}}=0.09$) performs better than the static McIntosh-Poisson method ($\text{BSS}_{\text{static}} = -0.09$); ii) low BSS values arise partly from both methods over-forecasting SC23 flares from the SC22 rates, symptomatic of $\geqslant$C1.0 rates in SC23 being on average $\approx$80% of those in SC22 (with $\geqslant$M1.0 being $\approx$50%); iii) applying a bias-correction factor to reduce the SC22 rates used in forecasting SC23 flares yields modest improvement in skill relative to climatology for both methods ($\mathrm{BSS}^{\mathrm{corr}}_{\mathrm{static}} = 0.09$ and $\mathrm{BSS}^{\mathrm{corr}}_{\mathrm{evolution}} = 0.20$) and improved forecast reliability diagrams.

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

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.

Determining the mechanisms responsible for the heating of the coronal plasma and maintaining and accelerating the solar wind are long standing goals in solar physics. There is a clear need to constrain the energy, mass and momentum flux through the solar corona and advance our knowledge of the physical process contributing to these fluxes. Furthermore, the accurate forecasting of Space Weather conditions at the near-Earth environment and, more generally, the plasma conditions of the solar wind throughout the heliosphere, require detailed knowledge of these fluxes in the near-Sun corona. Here we present a short case for a space-based imaging-spectrometer coronagraph, which will have the ability to provide synoptic information on the coronal environment and provide strict constraints on the mass, energy, and momentum flux through the corona. The instrument would ideally achieve cadences of $\sim10$~s, spatial resolution of 1" and observe the corona out to 2~$R_{\sun}$. Such an instrument will enable significant progress in our understanding of MHD waves throughout complex plasmas, as well as potentially providing routine data products to aid Space Weather forecasting.

For almost 20 years the physical nature of globally propagating waves in the solar corona (commonly called "EIT waves") has been controversial and subject to debate. Additional theories have been proposed over the years to explain observations that did not fit with the originally proposed fast-mode wave interpretation. However, the incompatibility of observations made using the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory with the fast-mode wave interpretation was challenged by differing viewpoints from the twin Solar Terrestrial Relations Observatory spacecraft and higher spatial/temporal resolution data from the Solar Dynamics Observatory. In this article, we reexamine the theories proposed to explain "EIT waves" to identify measurable properties and behaviours that can be compared to current and future observations. Most of us conclude that "EIT waves" are best described as fast-mode large-amplitude waves/shocks that are initially driven by the impulsive expansion of an erupting coronal mass ejection in the low corona.

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.

Sunspot groups are the main source of solar flares, with the energy to power them being supplied by magnetic-field evolution (e.g. flux emergence or twisting/shearing). To date, few studies have investigated the statistical relation between sunspot-group evolution and flaring, with none considering evolution in the McIntosh classification scheme. Here we present a statistical analysis of sunspot groups from Solar Cycle 22, focusing on 24-hour changes in the three McIntosh classification components. Evolution-dependent >C1.0, >M1.0, and >X1.0 flaring rates are calculated, leading to the following results: (i) flaring rates become increasingly higher for greater degrees of upward evolution through the McIntosh classes, with the opposite found for downward evolution; (ii) the highest flaring rates are found for upward evolution from larger, more complex, classes (e.g. Zurich D- and E-classes evolving upward to F-class produce >C1.0 rates of 2.66 +/- 0.28 and 2.31 +/- 0.09 flares per 24 hours, respectively); (iii) increasingly complex classes give higher rates for all flare magnitudes, even when sunspot groups do not evolve over 24 hours. These results support the hypothesis that injection of magnetic energy by flux emergence (i.e. increasing in Zurich or compactness classes) leads to a higher frequency and magnitude of flaring.

Context. The Sun is an active source of radio emission ranging from long duration radio bursts associated with solar flares and coronal mass ejections to more complex, short duration radio bursts such as solar S bursts, radio spikes and fibre bursts. While plasma emission is thought to be the dominant emission mechanism for most radio bursts, the electron-cyclotron maser (ECM) mechanism may be responsible for more complex, short-duration bursts as well as fine structures associated with long-duration bursts. Aims. We investigate the conditions for ECM in the solar corona by considering the ratio of the electron plasma frequency {\omega}p to the electron-cyclotron frequency {\Omega}e. The ECM is theoretically possible when {\omega}p/{\Omega}e < 1. Methods. Two-dimensional electron density, magnetic field, plasma frequency, and electron cyclotron frequency maps of the off- limb corona were created using observations from SDO/AIA and SOHO/LASCO, together with potential field extrapolations of the magnetic field. These maps were then used to calculate {\omega}p/{\Omega}e and Alfven velocity maps of the off-limb corona. Results. We found that the condition for ECM emission ({\omega}p/{\Omega}e < 1) is possible at heights < 1.07 R_sun in an active region near the limb; that is, where magnetic field strengths are > 40 G and electron densities are greater than 3x10^8 cm-3. In addition, we found comparatively high Alfv\'en velocities (> 0.02 c or > 6000 km s-1) at heights < 1.07 R_sun within the active region. Conclusions. This demonstrates that the condition for ECM emission is satisfied within areas of the corona containing large magnetic fields, such as the core of a large active region. Therefore, ECM could be a possible emission mechanism for high-frequency radio and microwave bursts.

The physical processes that trigger solar flares are not well understood and significant debate remains around processes governing particle acceleration, energy partition, and particle and energy transport. Observations at high resolution in energy, time, and space are required in multiple energy ranges over the whole course of many flares in order to build an understanding of these processes. Obtaining high-quality, co-temporal data from ground- and space- based instruments is crucial to achieving this goal and was the primary motivation for starting the Max Millennium program and Major Flare Watch (MFW) alerts, aimed at coordinating observations of all flares >X1 GOES X-ray classification (including those partially occulted by the limb). We present a review of the performance of MFWs from 1 February 2001 to 31 May 2010, inclusive, that finds: (1) 220 MFWs were issued in 3,407 days considered (6.5% duty cycle), with these occurring in 32 uninterrupted periods that typically last 2-8 days; (2) 56% of flares >X1 were caught, occurring in 19% of MFW days; (3) MFW periods ended at suitable times, but substantial gain could have been achieved in percentage of flares caught if periods had started 24 h earlier; (4) MFWs successfully forecast X-class flares with a true skill statistic (TSS) verification metric score of 0.500, that is comparable to a categorical flare/no-flare interpretation of the NOAA Space Weather Prediction Centre probabilistic forecasts (TSS = 0.488).

Large and magnetically complex sunspot groups are known to be associated with flares. To date, the Mount Wilson scheme has been used to classify sunspot groups based on their morphological and magnetic properties. The most flare prolific class, the delta sunspot-group, is characterised by opposite polarity umbrae within a common penumbra, separated by less than 2 degrees. In this article, we present a new system, called the Solar Monitor Active Region Tracker - Delta Finder (SMART-DF), that can be used to automatically detect and classify magnetic deltas in near-realtime. Using continuum images and magnetograms from the Helioseismic and Magnetic Imager (HMI) onboard NASA's Solar Dynamics Observatory (SDO), we first estimate distances between opposite polarity umbrae. Opposite polarity pairs having distances of less that 2 degrees are then identified, and if these pairs are found to share a common penumbra, they are identified as a magnetic delta configuration. The algorithm was compared to manual delta detections reported by the Space Weather Prediction Center (SWPC), operated by the National Oceanic and Atmospheric Administration (NOAA). SMART-DF detected 21 out of 23 active regions (ARs) that were marked as delta spots by NOAA during 2011 - 2012 (within +/- 60 degrees longitude). SMART-DF in addition detected five ARs which were not announced as delta spots by NOAA. The near-relatime operation of SMART-DF resulted in many deltas being identified in advance of NOAA's daily notification. SMART-DF will be integrated with SolarMonitor (www.solarmonitor.org) and the near-realtime information will be available to the public.

Cosmic rays and solar energetic particles may be accelerated to relativistic energies by shock waves in astrophysical plasmas. On the Sun, shocks and particle acceleration are often associated with the eruption of magnetized plasmoids, called coronal mass ejections (CMEs). However, the physical relationship between CMEs and shock particle acceleration is not well understood. Here, we use extreme ultraviolet, radio and white-light imaging of a solar eruptive event on 22 September 2011 to show that a CME-induced shock (Alfv\'en Mach number 2.4$^{+0.7}_{-0.8}$) was coincident with a coronal wave and an intense metric radio burst generated by intermittent acceleration of electrons to kinetic energies of 2-46 keV (0.1-0.4 c). Our observations show that plasmoid-driven quasi-perpendicular shocks are capable of producing quasi-periodic acceleration of electrons, an effect consistent with a turbulent or rippled plasma shock surface.

The continuous stream of data available from the Atmospheric Imaging Assembly (AIA) telescopes onboard the Solar Dynamics Observatory (SDO) spacecraft has allowed a deeper understanding of the Sun. However, the sheer volume of data has necessitated the development of automated techniques to identify and analyse various phenomena. In this article, we describe the Coronal Pulse Identification and Tracking Algorithm (CorPITA) for the identification and analysis of coronal "EIT waves". CorPITA uses an intensity-profile technique to identify the propagating pulse, tracking it throughout its evolution before returning estimates of its kinematics. The algorithm is applied here to a data-set from February 2011, allowing its capabilities to be examined and critiqued. This algorithm forms part of the SDO Feature Finding Team initiative and will be implemented as part of the Heliophysics Event Knowledgebase (HEK). This is the first fully automated algorithm to identify and track the propagating "EIT wave" rather than any associated phenomena and will allow a deeper understanding of this controversial phenomenon.

Super-Alfv\'enic shock waves associated with coronal mass ejections (CMEs) can produce radio emission known as Type II bursts. In the absence of direct imaging, accurate estimates of coronal electron densities, magnetic field strengths and Alfv\'en speeds are required in order to calculate the kinematics of shocks. To date, 1D radial models have been used, but these are not appropriate for shocks propagating in non-radial directions. Here, we study a coronal shock wave associated with a CME and Type II radio burst using 2D electron density and Alfv\'en speed maps to determine the locations that shocks are excited as the CME expands through the corona. Coronal density maps were obtained from emission measures derived from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) and polarized brightness measurements from the Large Angle and Spectrometric Coronagraph (LASCO) on board the Solar and Heliospheric Observatory (SOHO). Alfv\'en speed maps were calculated using these density maps and magnetic field extrapolations from the Helioseismic and Magnetic Imager (SDO/HMI). The computed density and Alfv\'en speed maps were then used to calculate the shock kinematics in non-radial directions. Using the kinematics of the Type II burst and associated shock, we find our observations to be consistent with the formation of a shock located at the CME flanks where the Alfv\'en speed has a local minimum. 1D density models are not appropriate for shocks that propagate non-radially along the flanks of a CME. Rather, the 2D density, magnetic field and Alfv\'en speed maps described here give a more accurate method for determining the fundamental properties of shocks and their relation to CMEs.

The study of solar eruptive events and associated phenomena is of great importance in the context of solar and heliophysics. Coronal mass ejections (CMEs) and coronal waves are energetic manifestations of the restructuring of the solar magnetic field and mass motion of the plasma. Characterising this motion is vital for deriving the dynamics of these events and thus understanding the physics driving their initiation and propagation. The development and use of appropriate methods for measuring event kinematics is therefore imperative. Traditional approaches to the study of CME and coronal wave kinematics do not return wholly accurate nor robust estimates of the true event kinematics and associated uncertainties. We highlight the drawbacks of these approaches, and demonstrate improved methods for accurate and reliable determination of the kinematics. The Savitzky-Golay filter is demonstrated as a more appropriate fitting technique for CME and coronal wave studies, and a residual resampling bootstrap technique is demonstrated as a statistically rigorous method for the determination of kinematic error estimates and goodness-of-fit tests. It is shown that the scatter on distance-time measurements of small sample size can significantly limit the ability to derive accurate and reliable kinematics. This may be overcome by (i) increasing measurement precision and sampling cadence, and (ii) applying robust methods for deriving the kinematics and reliably determining their associated uncertainties. If a priori knowledge exists and a pre-determined model form for the kinematics is available (or indeed any justified fitting-form to be tested against the data), then its precision can be examined using a bootstrapping technique to determine the confidence interval associated with the model/fitting parameters.

Solar flares are powered by energy stored in the coronal magnetic field, a portion of which is released when the field reconfigures into a lower energy state. Investigation of sunspot magnetic field topology during flare activity is useful to improve our understanding of flaring processes. Here we investigate the deviation of the non-linear field configuration from that of the linear and potential configurations, and study the free energy available leading up to and after a flare. The evolution of the magnetic field in NOAA region 10953 was examined using data from Hinode/SOT-SP, over a period of 12 hours leading up to and after a GOES B1.0 flare. Previous work on this region found pre- and post-flare changes in photospheric vector magnetic field parameters of flux elements outside the primary sunspot. 3D geometry was thus investigated using potential, linear force-free, and non-linear force-free field extrapolations in order to fully understand the evolution of the field lines. Traced field line geometrical and footpoint orientation differences show that the field does not completely relax to a fully potential or linear force-free state after the flare. Magnetic and free magnetic energies increase significantly ~ 6.5-2.5 hours before the flare by ~ 10^31 erg. After the flare, the non-linear force-free magnetic energy and free magnetic energies decrease but do not return to pre-flare 'quiet' values. The post-flare non-linear force-free field configuration is closer (but not equal) to that of the linear force-free field configuration than a potential one. However, the small degree of similarity suggests that partial Taylor relaxation has occurred over a time scale of ~ 3-4 hours.

The Sun Watcher with Active Pixels and Image Processing (SWAP) is an EUV solar telescope on board ESA's Project for Onboard Autonomy 2 (PROBA2) mission launched on 2 November 2009. SWAP has a spectral bandpass centered on 17.4 nm and provides images of the low solar corona over a 54x54 arcmin field-of-view with 3.2 arcsec pixels and an imaging cadence of about two minutes. SWAP is designed to monitor all space-weather-relevant events and features in the low solar corona. Given the limited resources of the PROBA2 microsatellite, the SWAP telescope is designed with various innovative technologies, including an off-axis optical design and a CMOS-APS detector. This article provides reference documentation for users of the SWAP image data.

Solar flares occur in complex sunspot groups, but it remains unclear how the probability of producing a flare of a given magnitude relates to the characteristics of the sunspot group. Here, we use Geostationary Operational Environmental Satellite X-ray flares and McIntosh group classifications from solar cycles 21 and 22 to calculate average flare rates for each McIntosh class and use these to determine Poisson probabilities for different flare magnitudes. Forecast verification measures are studied to find optimum thresholds to convert Poisson flare probabilities into yes/no predictions of cycle 23 flares. A case is presented to adopt the true skill statistic (TSS) as a standard for forecast comparison over the commonly used Heidke skill score (HSS). In predicting flares over 24 hr, the maximum values of TSS achieved are 0.44 (C-class), 0.53 (M-class), 0.74 (X-class), 0.54 (>=M1.0), and 0.46 (>=C1.0). The maximum values of HSS are 0.38 (C-class), 0.27 (M-class), 0.14 (X-class), 0.28 (>=M1.0), and 0.41 (>=C1.0). These show that Poisson probabilities perform comparably to some more complex prediction systems, but the overall inaccuracy highlights the problem with using average values to represent flaring rate distributions.

The solar surface and atmosphere are highly dynamic plasma environments, which evolve over a wide range of temporal and spatial scales. Large-scale eruptions, such as coronal mass ejections, can be accelerated to millions of kilometres per hour in a matter of minutes, making their automated detection and characterisation challenging. Additionally, there are numerous faint solar features, such as coronal holes and coronal dimmings, which are important for space weather monitoring and forecasting, but their low intensity and sometimes transient nature makes them problematic to detect using traditional image processing techniques. These difficulties are compounded by advances in ground- and space- based instrumentation, which have increased the volume of data that solar physicists are confronted with on a minute-by-minute basis; NASA's Solar Dynamics Observatory for example is returning many thousands of images per hour (~1.5 TB/day). This chapter reviews recent advances in the application of images processing techniques to the automated detection of active regions, coronal holes, filaments, CMEs, and coronal dimmings for the purposes of space weather monitoring and prediction.

Solar flares occur due to the sudden release of energy stored in active-region magnetic fields. To date, the pre-cursors to flaring are still not fully understood, although there is evidence that flaring is related to changes in the topology or complexity of an active region's magnetic field. Here, the evolution of the magnetic field in active region NOAA 10953 was examined using Hinode/SOT-SP data, over a period of 12 hours leading up to and after a GOES B1.0 flare. A number of magnetic-field properties and low-order aspects of magnetic-field topology were extracted from two flux regions that exhibited increased Ca II H emission during the flare. Pre-flare increases in vertical field strength, vertical current density, and inclination angle of ~ 8degrees towards the vertical were observed in flux elements surrounding the primary sunspot. The vertical field strength and current density subsequently decreased in the post-flare state, with the inclination becoming more horizontal by ~7degrees. This behaviour of the field vector may provide a physical basis for future flare forecasting efforts.

One of the most dramatic manifestations of solar activity are large-scale coronal bright fronts (CBFs) observed in extreme ultraviolet (EUV) images of the solar atmosphere. To date, the energetics and kinematics of CBFs remain poorly understood, due to the low image cadence and sensitivity of previous EUV imagers and the limited methods used to extract the features. In this paper, the trajectory and morphology of CBFs was determined in order to investigate the varying properties of a sample of CBFs, including their kinematics and pulse shape, dispersion, and dissipation. We have developed a semi-automatic intensity profiling technique to extract the morphology and accurate positions of CBFs in 2.5-10 min cadence images from STEREO/EUVI. The technique was applied to sequences of 171A and 195A images from STEREO/EUVI in order to measure the wave properties of four separate CBF events. Following launch at velocities of ~240-450kms^{-1} each of the four events studied showed significant negative acceleration ranging from ~ -290 to -60ms^{-2}. The CBF spatial and temporal widths were found to increase from ~50 Mm to ~200 Mm and ~100 s to ~1500 s respectively, suggesting that they are dispersive in nature. The variation in position-angle averaged pulse-integrated intensity with propagation shows no clear trend across the four events studied. These results are most consistent with CBFs being dispersive magnetoacoustic waves.

The interaction of open and closed field lines at coronal hole boundaries is widely accepted to be due to interchange magnetic reconnection. To date, it is unclear how the boundaries vary on short timescales and at what velocity this occurs. Here, we describe an automated boundary tracking method used to determine coronal hole boundary displacements on short timescales. The bound- ary displacements were found to be isotropic and to have typical expansion/contraction speeds of \leq2 km s^-1, which indicate magnetic reconnection rates of \leq 3 \times 10^-3. The observed displacements were used in conjunction with the interchange reconnection model to derive typical diffusion coeffi- cients of \leq 3 \times 10^13 cm^2 s^-1. These results are consistent with an interchange reconnection process in the low corona driven by the random granular motion of open and closed fields in the photosphere.

We present an automated system for detecting, tracking, and cataloging emerging active regions throughout their evolution and decay using SOHO Michelson Doppler Interferometer (MDI) magnetograms. The SolarMonitor Active Region Tracking (SMART) algorithm relies on consecutive image differencing to remove both quiet-Sun and transient magnetic features, and region-growing techniques to group flux concentrations into classifiable features. We determine magnetic properties such as region size, total flux, flux imbalance, flux emergence rate, Schrijver's R-value, R* (a modified version of R), and Falconer's measurement of non-potentiality. A persistence algorithm is used to associate developed active regions with emerging flux regions in previous measurements, and to track regions beyond the limb through multiple solar rotations. We find that the total number and area of magnetic regions on disk vary with the sunspot cycle. While sunspot numbers are a proxy to the solar magnetic field, SMART offers a direct diagnostic of the surface magnetic field and its variation over timescale of hours to years. SMART will form the basis of the active region extraction and tracking algorithm for the Heliophysics Integrated Observatory (HELIO).

The kinematics of a globally propagating disturbance (also known as an ``EIT wave") is discussed using Extreme UltraViolet Imager (EUVI) data Solar Terrestrial Relations Observatory (STEREO). We show for the first time that an impulsively generated propagating disturbance has similar kinematics in all four EUVI passbands (304, 171, 195, and 284 A). In the 304 A passband the disturbance shows a velocity peak of 238+/-20 kms-1 within ~28 minutes of its launch, varying in acceleration from 76 ms-2 to -102 ms-2. This passband contains a strong contribution from a Si XI line (303.32 A) with a peak formation temperature of ~1.6 MK. The 304 A emission may therefore be coronal rather than chromospheric in origin. Comparable velocities and accelerations are found in the coronal 195 A passband, while lower values are found in the lower cadence 284 A passband. In the higher cadence 171 A passband the velocity varies significantly, peaking at 475+/-47 kms-1 within ~20 minutes of launch, with a variation in acceleration from 816 ms-2 to -413 ms-2. The high image cadence of the 171 A passband (2.5 minutes compared to 10 minutes for the similar temperature response 195 A passband) is found to have a major effect on the measured velocity and acceleration of the pulse, which increase by factors of ~2 and ~10, respectively. This implies that previously measured values (e.g., using EIT) may have been underestimated. We also note that the disturbance shows strong reflection from a coronal hole in both the 171 and 195 A passbands. The observations are consistent with an impulsively generated fast-mode magnetoacoustic wave.

High-cadence, synchronized, multiwavelength optical observations of a solar active region (NOAA 10794) are presented. The data were obtained with the Dunn Solar Telescope at the National Solar Observatory/Sacramento Peak using a newly developed camera system : the Rapid Dual Imager. Wavelet analysis is undertaken to search for intensity related oscillatory signatures, and periodicities ranging from 20 to 370 s are found with significance levels exceeding 95%. Observations in the H-alpha blue wing show more penumbral oscillatory phenomena when compared to simultaneous G-band observations. The H-alpha oscillations are interpreted as the signatures of plasma motions with a mean velocity of 20 km/s. The strong oscillatory power over H-alpha blue-wing and G-band penumbral bright grains is an indication of the Evershed flow with frequencies higher than previously reported.

We seek to clarify the nature of running penumbral (RP) waves: are they chromospheric trans-sunspot waves or a visual pattern of upward-propagating waves? Full Stokes spectropolarimetric time series of the photospheric Si I 10827 \AA line and the chromospheric He I 10830 \AA multiplet were inverted using a Milne-Eddington atmosphere. Spatial pixels were paired between the outer umbral/inner penumbral photosphere and the penumbral chromosphere using inclinations retrieved by the inversion and the dual-height pairings of line-of-sight velocity time series were studied for signatures of wave propagation using a Fourier phase difference analysis. The dispersion relation for radiatively cooling acoustic waves, modified to incorporate an inclined propagation direction, fits well the observed phase differences between the pairs of photospheric and chromospheric pixels. We have thus demonstrated that RP waves are in effect low-beta slow-mode waves propagating along the magnetic field.

High-cadence optical observations of an H-alpha blue-wing bright point near solar AR NOAA 10794 are presented. The data were obtained with the Dunn Solar Telescope at the National Solar Observatory/Sacramento Peak using a newly developed camera system, the Rapid Dual Imager. Wavelet analysis is undertaken to search for intensity-related oscillatory signatures, and periodicities ranging from 15 to 370 s are found with significance levels exceeding 95%. During two separate microflaring events, oscillation sites surrounding the bright point are observed to twist. We relate the twisting of the oscillation sites to the twisting of physical flux tubes, thus giving rise to reconnection phenomena. We derive an average twist velocity of 8.1 km/s and detect a peak in the emitted flux between twist angles of 180 and 230 degrees.

Aims: The primary objective of this study is to search for and identify wave modes within a sunspot penumbra. Methods: Infrared spectropolarimetric time series data are inverted using a model comprising two atmospheric components in each spatial pixel. Fourier phase difference analysis is performed on the line-of-sight velocities retrieved from both components to determine time delays between the velocity signals. In addition, the vertical separation between the signals in the two components is calculated from the Stokes velocity response functions. Results: The inversion yields two atmospheric components, one permeated by a nearly horizontal magnetic field, the other with a less-inclined magnetic field. Time delays between the oscillations in the two components in the frequency range 2.5-4.5 mHz are combined with speeds of atmospheric wave modes to determine wave travel distances. These are compared to expected path lengths obtained from response functions of the observed spectral lines in the different atmospheric components. Fast-mode (i.e., modified p-mode) waves exhibit the best agreement with the observations when propagating toward the sunspot at an angle ~50 degrees to the vertical.

Spectropolarimetric time series data of the primary spot of active region NOAA 9448 were obtained in the Si I 10827 \AA line and the He I 10830 \AA multiplet with the Tenerife Infrared Polarimeter. Throughout the time series the spectrograph slit was fixed over a region covering umbra, a light bridge, penumbra, and quiet sun. We present speeds of running penumbral waves in the chromosphere, their relation to both photospheric and chromospheric umbral oscillations, and their dependence on the magnetic field topology.

Two sequences of solar images obtained by the Transition Region and Coronal Explorer in three UV passbands are studied using wavelet and Fourier analysis and compared to the photospheric magnetic flux measured by the Michelson Doppler Interferometer on the Solar Heliospheric Observatory to study wave behaviour in differing magnetic environments. Wavelet periods show deviations from the theoretical cutoff value and are interpreted in terms of inclined fields. The variation of wave speeds indicates that a transition from dominant fast-magnetoacoustic waves to slow modes is observed when moving from network into plage and umbrae. This implies preferential transmission of slow modes into the upper atmosphere, where they may lead to heating or be detected in coronal loops and plumes.

High time resolution observations of a white--light flare on the active star EQ PegB show evidence of intensity variations with a period of approximately 10 s. The period drifts to longer values during the decay phase of the flare. If the oscillation is interpreted as an impulsively--excited, standing--acoustic wave in a flare loop, the period implies a loop length of 1.7 Mm and 3.4 Mm for the case of the fundamental mode and the second harmonic, respectively. However, the small loop lengths imply a very high modulation depth making the acoustic interpretation unlikely. A more realistic interpretation may be that of a fast--MHD wave, with the modulation of the emission being due to the magnetic field. Alternatively, the variations could be due to a series of reconnection events. The periodic signature may then arise as a result of the lateral separation of individual flare loops or current sheets with oscillatory dynamics (i.e. periodic reconnection).

We present FUV and UV spectroscopic observations of AD Leonis, with the aim of investigating opacity effects in the transition regions of late-type stars. The C III lines in FUSE spectra show significant opacity during both the quiescent and flaring states of AD Leonis, with up to 30% of the expected flux being lost during the latter. Other FUSE emission lines tested for opacity include those of O VI, while C IV, Si IV and N V transitions observed with STIS are also investigated. These lines only reveal modest amounts of opacity with losses during flaring of up to 20%. Optical depths have been calculated for homogeneous and inhomogeneous geometries, giving path lengths of ~20-60 km and \~10-30 km, respectively, under quiescent conditions. However path lengths derived during flaring are ~2-3 times larger. These values are in excellent agreement with both estimates of the small-scale structure observed in the solar transition region, and path lengths derived previously for several other active late-type stars.

Recent calculations of atomic data for Fe XV have been used to generate theoretical line ratios involving n = 3-4 transitions in the soft X-ray spectral region (52-83 A), for a wide range of electron temperatures and densities applicable to solar and stellar coronal plasmas. A comparison of these with solar flare observations from a rocket-borne spectrograph (XSST) reveals generally good agreement between theory and experiment. In particular, the 82.76 A emission line in the XSST spectrum is identified, for the first time to our knowledge in an astrophysical source. Most of the Fe XV transitions which are blended have had the species responsible clearly identified, although there remain a few instances where this has not been possible. The line ratio calculations are also compared with a co-added spectrum of Capella obtained with the Chandra satellite, which is probably the highest signal-to-noise observation achieved for a stellar source in the 25-175 A soft X-ray region. Good agreement is found between theory and experiment, indicating that the Fe XV lines are reliably detected in Chandra spectra, and hence may be employed as diagnostics to determine the temperature and/or density of the emitting plasma. However the line blending in the Chandra data is such that individual emission lines are difficult to measure accurately, and fluxes may only be reliably determined via detailed profile fitting of the observations. The co-added Capella spectrum is made available to hopefully encourage further exploration of the soft X-ray region in astronomical sources.

Previously, large discrepancies have been found between theory and observation for Fe XV emission line ratios in solar flare spectra covering the 224-327 A wavelength range, obtained by the Naval Research Laboratory's S082A instrument on board Skylab. These discrepancies have been attributed to either errors in the adopted atomic data or the presence of additional atomic processes not included in the modelling, such as fluorescence. However our analysis of these plus other S082A flare observations (the latter containing Fe XV transitions between 321-482 A), performed using the most recent Fe XV atomic physics calculations in conjunction with a CHIANTI synthetic flare spectrum, indicate that blending of the lines is primarily responsible for the discrepancies. As a result, most Fe XV lines cannot be employed as electron density diagnostics for solar flares, at least at the spectral resolution of S082A and similar instruments (i.e. ~ 0.1 A). An exception is the intensity ratio I(321.8 A)/I(327.0 A), which appears to provide good estimates of the electron density at this spectral resolution.

Simultaneous observations of explosive chromospheric evaporation are presented using data from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and the Coronal Diagnostic Spectrometer (CDS) onboard SOHO. For the first time, co-spatial imaging and spectroscopy have been used to observe explosive evaporation within a hard X-ray emitting region. RHESSI X-ray images and spectra were used to determine the flux of non-thermal electrons accelerated during the impulsive phase of an M2.2 flare. Assuming a thick-target model, the injected electron spectrum was found to have a spectral index of ~7.3, a low energy cut-off of ~20 keV, and a resulting flux of >4x10^10 ergs cm^-2 s^-1. The dynamic response of the atmosphere was determined using CDS spectra, finding a mean upflow velocity of 230+/-38 km s^-1 in Fe XIX (592.23A), and associated downflows of 36+/-16 km s^-1 and 43+/-22 km s^-1 at chromospheric and transition region temperatures, respectively, relative to an averaged quiet-Sun spectra. The errors represent a 1 sigma dispersion. The properties of the accelerated electron spectrum and the corresponding evaporative velocities were found to be consistent with the predictions of theory.

A detailed study is presented of the decaying solar active region NOAA 10103 observed with the Coronal Diagnostic Spectrometer (CDS), the Michelson Doppler Imager (MDI) and the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory (SOHO). Electron density maps formed using Si X (356.03A/347.41A) show that the density varies from ~10^10 cm^-3 in the active region core, to ~7x10^8 cm^-3 at the region boundaries. Over the five days of observations, the average electron density fell by ~30%. Temperature maps formed using Fe XVI(335.41A)/Fe XIV(334.18A) show electron temperatures of \~2.34x10^6 K in the active region core, and ~2.10x10^6 K at the region boundaries. Similarly to the electron density, there was a small decrease in the average electron temperature over the five day period. The radiative, conductive, and mass flow losses were calculated and used to determine the resultant heating rate (P_H). Radiative losses were found to dominate the active region cooling process. As the region decayed, the heating rate decreased by almost a factor of five between the first and last day of observations. The heating rate was then compared to the total unsigned magnetic flux (Phi_tot), yielding a power-law of the form P_H ~ Phi_tot^(0.81 +/- 0.32). This result suggests that waves rather than nanoflares may be the dominant heating mechanism in this active region.

New calculations of radiative rates and electron impact excitation cross sections for Fe XI are used to derive emission line intensity ratios involving 3s^23p^4 - 3s^23p^33d transitions in the 180-223 A wavelength range. These ratios are subsequently compared with observations of a solar active region, obtained during the 1995 flight Solar EUV Research Telescope and Spectrograph (SERTS). The version of SERTS flown in 1995 incorporated a multilayer grating that enhanced the instrumental sensitivity for features in the 170 - 225 A wavelength range, observed in second-order between 340 and 450 A. This enhancement led to the detection of many emission lines not seen on previous SERTS flights, which were measured with the highest spectral resolution (0.03 A) ever achieved for spatially resolved active region spectra in this wavelength range. However, even at this high spectral resolution, several of the Fe XI lines are found to be blended, although the sources of the blends are identified in the majority of cases. The most useful Fe XI electron density diagnostic line intensity ratio is I(184.80 A)/I(188.21 A). This ratio involves lines close in wavelength and free from blends, and which varies by a factor of 11.7 between N_e = 10^9 and 10^11 cm^-3, yet shows little temperature sensitivity. An unknown line in the SERTS spectrum at 189.00 A is found to be due to Fe XI, the first time (to our knowledge) this feature has been identified in the solar spectrum. Similarly, there are new identifications of the Fe XI 192.88, 198.56 and 202.42 A features, although the latter two are blended with S VIII/Fe XII and Fe XIII, respectively.