We developed a reliable probabilistic solar flare forecasting model using a deep neural network, named Deep Flare Net-Reliable (DeFN-R). The model can predict the maximum classes of flares that occur in the following 24 h after observing images, along with the event occurrence probability. We detected active regions from 3x10^5 solar images taken during 2010-2015 by Solar Dynamic Observatory and extracted 79 features for each region, which we annotated with flare occurrence labels of X-, M-, and C-classes. The extracted features are the same as used by Nishizuka et al. (2018); for example, line-of-sight/vector magnetograms in the photosphere, brightening in the corona, and the X-ray emissivity 1 and 2 h before an image. We adopted a chronological split of the database into two for training and testing in an operational setting: the dataset in 2010-2014 for training and the one in 2015 for testing. DeFN-R is composed of multilayer perceptrons formed by batch normalizations and skip connections. By tuning optimization methods, DeFN-R was trained to optimize the Brier skill score (BSS). As a result, we achieved BSS = 0.41 for >=C-class flare predictions and 0.30 for >=M-class flare predictions by improving the reliability diagram while keeping the relative operating characteristic curve almost the same. Note that DeFN is optimized for deterministic prediction, which is determined with a normalized threshold of 50%. On the other hand, DeFN-R is optimized for a probability forecast based on the observation event rate, whose probability threshold can be selected according to users' purposes.

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 investigated the relationship between the spectral structures of type II solar radio bursts in the hectometric and kilometric wavelength ranges and solar energetic particles (SEPs). To examine the statistical relationship between type II bursts and SEPs, we selected 26 coronal mass ejection (CME) events with similar characteristics (e.g., initial speed, angular width, and location) observed by the Large Angle and Spectrometric Coronagraph (LASCO), regardless of the characteristics of the corresponding type II bursts and the SEP flux. Then, we compared associated type II bursts observed by the Radio and Plasma Wave Experiment (WAVES) onboard the Wind spacecraft and the SEP flux observed by the Geostationary Operational Environmental Satellite (GOES) orbiting around the Earth. We found that the bandwidth of the hectometric type II bursts and the peak flux of the SEPs has a positive correlation (with a correlation coefficient of 0.64). This result supports the idea that the nonthermal electrons of type II bursts and the nonthermal ions of SEPs are generated by the same shock and suggests that more SEPs may be generated for a wider or stronger CME shock with a longer duration. Our result also suggests that considering the spectral structures of type II bursts can improve the forecasting accuracy for the peak flux of gradual SEPs.

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.

In this work, we investigate the reliability of the probabilistic binary forecast. We mathematically prove that a necessary, but not sufficient, condition for achieving a reliable probabilistic forecast is maximizing the Peirce skill score (PSS) at the threshold probability of the climatological base rate. The condition is confirmed by using artificially synthesized forecast-outcome pair data and previously published probabilistic solar flare forecast models. The condition gives a partial answer as to why some probabilistic forecast system lack reliability, because the system, which does not satisfy the proved condition, can never be reliable. Therefore, the proved condition is very important for the developers of a probabilistic forecast system. The result implies that those who want to develop a reliable probabilistic forecast system must adjust or train the system so as to maximize PSS near the threshold probability of the climatological base rate.

Coronal mass ejections (CMEs) cause disturbances in the environment of the Earth when they arrive at the Earth. However, the prediction of the arrival of CMEs still remains a challenge. We have developed an interplanetary scintillation (IPS) estimation system based on a global magnetohydrodynamic (MHD) simulation of the inner heliosphere to predict the arrival time of CMEs. In this system, the initial speed of a CME is roughly derived from white light coronagraph observations. Then, the propagation of the CME is calculated by a global MHD simulation. The IPS response is estimated by the three-dimensional density distribution of the inner heliosphere derived from the MHD simulation. The simulated IPS response is compared with the actual IPS observations made by the Institute for Space-Earth Environmental Research, Nagoya University, and shows good agreement with that observed. We demonstrated how the simulation system works using a halo CME event generated by a X9.3 flare observed on September 5, 2017. We find that the CME simulation that best estimates the IPS observation can more accurately predict the time of arrival of the CME at the Earth. These results suggest that the accuracy of the CME arrival time can be improved if our current MHD simulations include IPS data.

We developed a solar flare prediction model using a deep neural network (DNN), named Deep Flare Net (DeFN). The model can calculate the probability of flares occurring in the following 24 h in each active region, which is used to determine the most likely maximum classes of flares via a binary classification (e.g., >=M class versus <M class or >=C class versus <C class). From 3x10^5 observation images taken during 2010-2015 by Solar Dynamic Observatory, we automatically detected sunspots and calculated 79 features for each region, to which flare occurrence labels of X-, M-, and C-class were attached. We adopted the features used in Nishizuka et al. (2017) and added some features for operational prediction: coronal hot brightening at 131 A (T>=10^7 K) and the X-ray and 131 A intensity data 1 and 2 h before an image. For operational evaluation, we divided the database into two for training and testing: the dataset in 2010-2014 for training and the one in 2015 for testing. The DeFN model consists of deep multilayer neural networks, formed by adapting skip connections and batch normalizations. To statistically predict flares, the DeFN model was trained to optimize the skill score, i.e., the true skill statistic (TSS). As a result, we succeeded in predicting flares with TSS=0.80 for >=M-class flares and TSS=0.63 for >=C-class flares. Note that in usual DNN models, the prediction process is a black box. However, in the DeFN model, the features are manually selected, and it is possible to analyze which features are effective for prediction after evaluation.

In this article, we discuss a verification study of an operational solar flare forecast in the Regional Warning Center (RWC) Japan. The RWC Japan has been issuing four-categorical deterministic solar flare forecasts for a long time. In this forecast verification study, we used solar flare forecast data accumulated over 16 years (from 2000 to 2015). We compiled the forecast data together with solar flare data obtained with the Geostationary Operational Environmental Satellites (GOES). Using the compiled data sets, we estimated some conventional scalar verification measures with 95% confidence intervals. We also estimated a multi-categorical scalar verification measure. These scalar verification measures were compared with those obtained by the persistence method and recurrence method. As solar activity varied during the 16 years, we also applied verification analyses to four subsets of forecast-observation pair data with different solar activity levels. We cannot conclude definitely that there are significant performance difference between the forecasts of RWC Japan and the persistence method, although a slightly significant difference is found for some event definitions. We propose to use a scalar verification measure to assess the judgment skill of the operational solar flare forecast. Finally, we propose a verification strategy for deterministic operational solar flare forecasting.

We have developed a digital fast Fourier transform (FFT) spectrometer made of an analog-to-digital converter (ADC) and a field-programmable gate array (FPGA). The base instrument has independent ADC and FPGA modules, which allow us to implement different spectrometers in a relatively easy manner. Two types of spectrometers have been instrumented, one with 4.096 GS/s sampling speed and 2048 frequency channels and the other with 2.048 GS/s sampling speed and 32768 frequency channels. The signal processing in these spectrometers has no dead time and the accumulated spectra are recorded in external media every 8 ms. A direct sampling spectroscopy up to 8 GHz is achieved by a microwave track-and-hold circuit, which can reduce the analog receiver in front of the spectrometer. Highly stable spectroscopy with a wide dynamic range was demonstrated in a series of laboratory experiments and test observations of solar radio bursts.

We developed a flare prediction model using machine learning, which is optimized to predict the maximum class of flares occurring in the following 24 h. Machine learning is used to devise algorithms that can learn from and make decisions on a huge amount of data. We used solar observation data during the period 2010-2015, such as vector magnetogram, ultraviolet (UV) emission, and soft X-ray emission taken by the Solar Dynamics Observatory and the Geostationary Operational Environmental Satellite. We detected active regions from the full-disk magnetogram, from which 60 features were extracted with their time differentials, including magnetic neutral lines, the current helicity, the UV brightening, and the flare history. After standardizing the feature database, we fully shuffled and randomly separated it into two for training and testing. To investigate which algorithm is best for flare prediction, we compared three machine learning algorithms: the support vector machine (SVM), k-nearest neighbors (k-NN), and extremely randomized trees (ERT). The prediction score, the true skill statistic (TSS), was higher than 0.9 with a fully shuffled dataset, which is higher than that for human forecasts. It was found that k-NN has the highest performance among the three algorithms. The ranking of the feature importance showed that the previous flare activity is most effective, followed by the length of magnetic neutral lines, the unsigned magnetic flux, the area of UV brightening, and the time differentials of features over 24 h, all of which are strongly correlated with the flux emergence dynamics in an active region.

Solar energetic particles (SEPs) are one of the extreme space weather phenomena. A huge SEP event increases the radiation dose received by aircrews, who should be warned of such events as early as possible. We developed a warning system for aviation exposure to SEPs. This article describes one component of the system, which calculates the temporal evolution of the SEP intensity and the spectrum immediately outside the terrestrial magnetosphere. To achieve this, we performed numerical simulations of SEP transport in interplanetary space, in which interplanetary SEP transport is described by the focused transport equation. We developed a new simulation code to solve the equation using a set of stochastic differential equations. In the code, the focused transport equation is expressed in a magnetic field line coordinate system, which is a non-orthogonal curvilinear coordinate system. An inverse Gaussian distribution is employed as the injection profile of SEPs at an inner boundary located near the Sun. We applied the simulation to observed SEP events as a validation test. The results show that our simulation can closely reproduce observational data for the temporal evolution of particle intensity. By employing the code, we developed the WArning System for AVIation Exposure to Solar energetic particles (WASAVIES).

The variation in the figure axis in the Earth owing to the geophysical causes directly reflects the physical state of the Earth and therefore, it is important for understanding the Earth. Of the variations in the figure axis arising from different causes, the variation owing to the Earth's geophysical causes is determined from polar motion data. Equations that give the relation between the figure axis and rotational axis are derived and the motion of the figure pole with respect to the Earth's surface or the figure axis in the Earth is determined by using these equations. Estimation of error suggests that the accuracy of the determination for the position of the figure pole at each moment is sufficiently good and that even a variation that occurs during a short duration of time, such as one or two days, is significant. The determined position of the figure pole exhibits a rapid and complicated motion in a short duration. On the other hand, a long-range average of the pole shows stable annual and semi-annual variations, but it is quite independent of the Chandler wobble.

The polar motion data is analyzed to obtain accurate position of the figure axis referred to the Earth-fixed frame. The variation of the figure axis should be the basic object to which the geophysical events are linked. By the method of rigid dynamics, the relation between the rotational and the figure axes is derived. The polar motion is the motion of the rotational axis on the Earth's surface and therefore the exact position of the figure axis on the surface at each moment is obtained from the polar motion. Since the accuracy of the recent data for the polar motion is very high, the obtained position of the figure axis is considered to keep good quality. As an average for a long duration of time, the obtained figure axis exhibits stable annual and semi-annual variations while it has no component with the period of the Chandler wobble. Besides this general feature, it shows different patterns from year to year as well as small irregularities with shorter periods, most of which are considered to be significant but not attributed to the observational errors. Further, a simple model is introduced for the cause of the seasonal variation of the figure axis. The model explains both the variation of the figure axis obtained above and that of the rotational speed which is so far known.

We show examples of excitation of coronal waves by flare-related abrupt eruptions of magnetic rope structures. The waves presumably rapidly steepened into shocks and freely propagated afterwards like decelerating blast waves that showed up as Moreton waves and EUV waves. We propose a simple quantitative description for such shock waves to reconcile their observed propagation with drift rates of metric type II bursts and kinematics of leading edges of coronal mass ejections (CMEs). Taking account of different plasma density falloffs for propagation of a wave up and along the solar surface, we demonstrate a close correspondence between drift rates of type II bursts and speeds of EUV waves, Moreton waves, and CMEs observed in a few known events.

On 17 January 2010, STEREO-B observed in extreme ultraviolet (EUV) and white light a large-scale dome-shaped expanding coronal transient with perfectly connected off-limb and on-disk signatures. Veronig et al. (2010, ApJL 716, 57) concluded that the dome was formed by a weak shock wave. We have revealed two EUV components, one of which corresponded to this transient. All of its properties found from EUV, white light, and a metric type II burst match expectations for a freely expanding coronal shock wave including correspondence to the fast-mode speed distribution, while the transient sweeping over the solar surface had a speed typical of EUV waves. The shock wave was presumably excited by an abrupt filament eruption. Both a weak shock approximation and a power-law fit match kinematics of the transient near the Sun. Moreover, the power-law fit matches expansion of the CME leading edge up to 24 solar radii. The second, quasi-stationary EUV component near the dimming was presumably associated with a stretched CME structure; no indications of opening magnetic fields have been detected far from the eruption region.

We discuss a mechanism for cosmic ray penetration into an interplanetary magnetic flux rope, particularly the effect of the finite Larmor radius and magnetic field irregularities. First, we derive analytical solutions for cosmic ray behavior inside a magnetic flux rope, on the basis of the Newton-Lorentz equation of a particle, to investigate how cosmic rays penetrate magnetic flux ropes under an assumption of there being no scattering by small-scale magnetic field irregularities. Next, we perform a numerical simulation of a cosmic ray penetration into an interplanetary magnetic flux rope by adding small-scale magnetic field irregularities. This simulation shows that a cosmic ray density distribution is greatly different from that deduced from a guiding center approximation because of the effect of the finite Larmor radius and magnetic field irregularities for the case of a moderate to large Larmor radius compared to the flux rope radius.

This article discusses statistical models for solar flare interval distribution in individual active regions. We analyzed solar flare data in 55 active regions that are listed in the GOES soft X-ray flare catalog. We discuss some problems with a conventional procedure to derive probability density functions from any data set and propose a new procedure, which uses the maximum likelihood method and Akaike Information Criterion (AIC) to objectively compare some competing probability density functions. We found that lognormal and inverse Gaussian models are more likely models than the exponential model for solar flare interval distribution in individual active regions. The results suggest that solar flares do not occur randomly in time; rather, solar flare intervals appear to be regulated by solar flare mechanisms. We briefly mention a probabilistic solar flare forecasting method as an application of a solar flare interval distribution analysis.

The spacetime structure of the spatially uniformly expanding universe is described in terms of a kind of global space and global time instead of the space and time we usually recognize. The global space at some instant is a space in which the global time is equal and the global time is equal to the proper time of any point moving together with the expansion of the universe which has elapsed since the big bang. The universe descibed in terms of the global space and time thus defined is open spatially and the Einstein's equations give a solution that the space expands eternally regardless of the mass density of the universe indicating that any expanding universe is open dynamically as well, as it is the case only when the density is less than the critical value in the standard model. In fact the critical density has not any particular meaning in the new frame. Finally the equations for the light path in the new frame are studied and are found similar to those in the standard model.

Among the major uncertainties involved in the Chandrasekhar mass models for Type Ia supernovae are the companion star of the accreting white dwarf (or the accretion rate that determines the carbon ignition density) and the flame speed after ignition. We present nucleosynthesis results from relatively slow deflagration (1.5 - 3 % of the sound speed) to constrain the rate of accretion from the companion star. Because of electron capture, a significant amount of neutron-rich species such as ^{54}Cr, ^{50}Ti, ^{58}Fe, ^{62}Ni, etc. are synthesized in the central region. To avoid the too large ratios of ^{54}Cr/^{56}Fe and ^{50}Ti/^{56}Fe, the central density of the white dwarf at thermonuclear runaway must be as low as \ltsim 2 \e9 \gmc. Such a low central density can be realized by the accretion as fast as $\dot M \gtsim 1 \times 10^{-7} M_\odot yr^{-1}$. These rapidly accreting white dwarfs might correspond to the super-soft X-ray sources.

Presupernova evolution and explosive nucleosynthesis in massive stars for main-sequence masses from 13 $M_\odot$ to 70 $M_\odot$ are calculated. We examine the dependence of the supernova yields on the stellar mass, $^{12}C(\alpha, \gamma) ^{16}O}$ rate, and explosion energy. The supernova yields integrated over the initial mass function are compared with the solar abundances.