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.

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.

Accurate forecasting of the arrival time and subsequent geomagnetic impacts of Coronal Mass Ejections (CMEs) at Earth is an important objective for space weather forecasting agencies. Recently, the CME Arrival and Impact working team has made significant progress towards defining community-agreed metrics and validation methods to assess the current state of CME modeling capabilities. This will allow the community to quantify our current capabilities and track progress in models over time. Firstly, it is crucial that the community focuses on the collection of the necessary metadata for transparency and reproducibility of results. Concerning CME arrival and impact we have identified 6 different metadata types: 3D CME measurement, model description, model input, CME (non-)arrival observation, model output data and metrics and validation methods. Secondly, the working team has also identified a validation time period, where all events within the following two periods will be considered: 1 January 2011-31 December 2012 and January 2015-31 December 2015. Those two periods amount to a total of about 100 hit events at Earth and a large amount of misses. Considering a time period will remove any bias in selecting events and the event set will represent a sample set that will not be biased by user selection. Lastly, we have defined the basic metrics and skill scores that the CME Arrival and Impact working team will focus on.

Magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. Alfven waves and magneto-sonic waves are particular classes of MHD waves. These wave modes are clearly different and have pure properties in uniform plasmas of infinite extent only. Due to plasma non-uniformity MHD waves have mixed properties and cannot be classified as pure Alfven or magneto-sonic waves. However, vorticity is a quantity unequivocally related to Alfven waves as compression is for magneto-sonic waves. Here, we investigate MHD waves superimposed on a one-dimensional non-uniform straight cylinder with constant magnetic field. For a piecewise constant density profile we find that the fundamental radial modes of the non-axisymmetric waves have the same properties as surface Alfven waves at a true discontinuity in density. Contrary to the classic Alfveen waves in a uniform plasma of infinite extent, vorticity is zero everywhere except at the cylinder boundary. If the discontinuity in density is replaced with a continuous variation of density, vorticity is spread out over the whole interval with non-uniform density. The fundamental radial modes of the non-axisymmetric waves do not need compression to exist unlike the radial overtones. In thin magnetic cylinders the fundamental radial modes of the non-axisymmetric waves with phase velocities between the internal and the external Alfven velocities can be considered as surface Alfven waves. On the contrary, the radial overtones can be related to fast-like magneto-sonic modes.

Context. Magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. In magnetic waveguides resonant absorption due to plasma inhomogeneity naturally transfers wave energy from large-scale motions to small-scale motions. In the cooler parts of the solar atmosphere as, e.g., the chromosphere, effects due to partial ionization may be relevant for wave dynamics and heating. Aims. We study resonant Alfven waves in partially ionized plasmas. Methods. We use the multifluid equations in the cold plasma approximation. We investigate propagating resonant MHD waves in partially ionized flux tubes. We use approximate analytical theory based on normal modes in the thin tube and thin boundary approximations along with numerical eigenvalue computations. Results. We find that the jumps of the wave perturbations across the resonant layer are the same as in fully ionized plasmas. The damping length due to resonant absorption is inversely proportional to the frequency, while that due to ion-neutral collisions is inversely proportional to the square of the frequency. For observed frequencies in the solar atmosphere, the amplitude of MHD kink waves is more efficiently damped by resonant absorption than by ion-neutral collisions. Conclusions. Most of the energy carried by chromospheric kink waves is converted into localized azimuthal Alfven waves that can deposit energy in the coronal medium. The dissipation of wave energy in the chromosphere due to ion-neutral collisions is only effective for high-frequency waves. The chromosphere acts as a filter for kink waves with periods shorter than 10 s.

Oscillations in coronal loops are usually interpreted in terms of uncoupled magnetohydrodynamic (MHD) waves. Examples of these waves are standing transverse motions, interpreted as the kink MHD modes, and propagating slow modes, commonly reported at the loop footpoints. Here we study a simple system in which fast and slow MHD waves are coupled. The goal is to understand the fingerprints of the coupling when boundary conditions are imposed in the model. The reflection problem of a fast and slow MHD wave interacting with a rigid boundary, representing the line-tying effect of the photosphere, is analytically investigated. Both propagating and standing waves are analysed and the time-dependent problem of the excitation of these waves is considered. An obliquely incident fast MHD wave on the photosphere inevitably generates a slow mode. The frequency of the generated slow mode at the photosphere is exactly the same as the frequency of the incident fast MHD mode, but its wavelength is much smaller, assuming that the sound speed is smaller than the Alfv\'en speed. The main signatures of the generated slow wave are density fluctuations at the loop footpoints. We have derived a simple formula that relates the velocity amplitude of the transverse standing mode with the density enhancements at the footpoints due to the driven slow modes. Using these results it is shown that there are possible evidences in the observations of the coupling between these two modes.

Resonant absorption and mode conversion are both extensively studied mechanisms for wave "absorption" in solar magnetohydrodynamics (MHD). But are they really distinct? We re-examine a well-known simple resonant absorption model in a cold MHD plasma that places the resonance inside an evanescent region. The normal mode solutions display the standard singular resonant features. However, these same normal modes may be used to construct a ray bundle which very clearly undergoes mode conversion to an Alfv\'en wave with no singularities. We therefore conclude that resonant absorption and mode conversion are in fact the same thing, at least for this model problem. The prime distinguishing characteristic that determines which of the two descriptions is most natural in a given circumstance is whether the converted wave can provide a net escape of energy from the conversion/absorption region of physical space. If it cannot, it is forced to run away in wavenumber space instead, thereby generating the arbitrarily small scales in situ that we recognize as fundamental to resonant absorption and phase mixing. On the other hand, if the converted wave takes net energy way, singularities do not develop, though phase mixing may still develop with distance as the wave recedes.

Magnetohydrodynamic (MHD) waves are often reported in the solar atmosphere and usually classified as slow, fast, or Alfv\'en. The possibility that these waves have mixed properties is often ignored. The goal of this work is to study and determine the nature of MHD kink waves. This is done by calculating the frequency, the damping rate and the eigenfunctions of MHD kink waves for three widely different MHD waves cases: a compressible pressure-less plasma, an incompressible plasma and a compressible plasma with non-zero plasma pressure which allows for MHD radiation. In all three cases the frequency and the damping rate are for practical purposes the same as they differ at most by terms proportional to $(k_z R)^2$. In the magnetic flux tube the kink waves are in all three cases, to a high degree of accuracy incompressible waves with negligible pressure perturbations and with mainly horizontal motions. The main restoring force of kink waves in the magnetised flux tube is the magnetic tension force. The total pressure gradient force cannot be neglected except when the frequency of the kink wave is equal or slightly differs from the local Alfv\'{e}n frequency, i.e. in the resonant layer. Kink waves are very robust and do not care about the details of the MHD wave environment. The adjective fast is not the correct adjective to characterise kink waves. If an adjective is to be used it should be Alfv\'{e}nic. However, it is better to realize that kink waves have mixed properties and cannot be put in one single box.

First results from a high-resolution three-dimensional nonlinear numerical study of the kink oscillation are presented. We show in detail the development of a shear instability in an untwisted line-tied magnetic flux tube. The instability produces significant deformations of the tube boundary. An extended transition layer may naturally evolve as a result of the shear instability at a sharp transition between the flux tube and the external medium. We also discuss the possible effects of the instability on the process of resonant absorption when an inhomogeneous layer is included in the model. One of the implications of these results is that the azimuthal component of the magnetic field of a stable flux tube in the solar corona, needed to prevent the shear instability, is probably constrained to be in a very specific range.

We study the excitation and damping of transverse oscillations in a multi-stranded model of a straight line-tied coronal loop. The transverse geometry of our equilibrium configuration is quite irregular and more realistic than the usual cylindrical loop model. By numerically solving the time-dependent ideal magnetohydrodynamic equations in two dimensions we show how the global motion of the whole bundle of strands, excited by an external disturbance, is converted into localized Alfv\'enic motions due to the process of resonant absorption. This process produces the attenuation of the transverse oscillations. At any location in the structure two dominant frequencies are found, the frequency of the global mode, or quasi-mode, and the local Alfv\'en frequency. We find that the mechanism of mode conversion, due to the coupling between fast and Alfv\'en waves, is not compromised by the complicated geometry of the model. We also show that it is possible to have energy conversion not only at the external edge of the composite loop but also inside the structure. The implications of these results and their relationship with the observations are discussed.

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