We studied a flare onset process in terms of stability of a three-dimensional (3D) magnetic field in active region 12371 producing an eruptive M6.5 flare in 2015 June 22. In order to reveal the 3D magnetic structure, we first extrapolated the 3D coronal magnetic fields based on time series of the photospheric vector magnetic fields under a nonlinear force-free field (NLFFF) approximation. The NLFFFs nicely reproduced the observed sigmoidal structure which is widely considered to be preeruptive magnetic configuration. In particular, we found that the sigmoid is composed of two branches of sheared arcade loops. On the basis of the NLFFFs, we investigated the sheared arcade loops to explore the onset process of the eruptive flare using three representative magnetohydrodynamic instabilities: the kink, torus, and double arc instabilities (DAI). The DAI, recently proposed by Ishiguro & Kusano, is a double arc loop that can be more easily destabilized than a torus loop. Consequently, the NLFFFs are found to be quite stable against the kink and torus instabilities. However, the sheared arcade loops formed prior to the flare possibly become unstable against the DAI. As a possible scenario for the onset process of the M6.5 flare, we suggest a three-step process: (1) double arc loops are formed by the sheared arcade loops through the tether-cutting reconnection during an early phase of the flare, (2) the DAI contributes to the expansion of destabilized double arc loops, and (3) finally, the torus instability makes the full eruption.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

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 nonlinear force-free field (NLFFF) extrapolation is widely used to reconstruct the three-dimensional magnetic field in the solar corona from the observed photospheric magnetic field. However, the pressure gradient and gravitational forces are ignored in the NLFFF model, even though the photospheric and chromospheric magnetic fields are not in general force-free. Here we develop a magnetohydrodynamic (MHD) relaxation method that reconstructs the solar atmospheric (chromospheric and coronal) magnetic field as a non-force-free magnetic field (NFFF) in magnetohydrostatic equilibrium where the Lorentz, pressure gradient, and gravitational forces are balanced. The system of basic equations for the MHD relaxation method is derived, and mathematical properties of the system are investigated. A robust numerical solver for the system is constructed based on the modern high-order shock capturing scheme. Two-dimensional numerical experiments that include the pressure gradient and gravitational forces are also demonstrated.

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.

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.

Coronal magnetic fields are responsible for the onset of solar flares and solar eruptions. However, the type of magnetic field parameters that can be used to measure the critical condition for a solar eruption is still unclear. As an effort to understand the possible condition for a solar flare, we have examined the non-dimensional parameter $\kappa$ introduced by Ishiguro & Kusano (2017), which contains information about magnetic twist distribution and magnetic flux in an active region (AR). We introduce a new parameter $\kappa^\ast$, as a proxy for $\kappa$, and we have analyzed the evolution of $\kappa^\ast$ during the flaring period of an AR using the nonlinear force-free field (NLFFF) extrapolated from the photospheric vector magnetic field data. Using data from the Solar Dynamics Observatory (SDO)/Helioseismic and Magnetic Imager (HMI), we have calculated $\kappa^\ast$ for the AR NOAA 11158 during its three-day flaring period. We found that $\kappa^\ast$ increased to a certain level before two large flares and decreased significantly after their onset. The results suggest that $\kappa^\ast$ may be used as an indicator of the necessary condition for the onset of a solar eruption in the AR. Based on this result, we propose a new method to assess the possibility of a large solar eruption from an AR by combining the parameter $\kappa^\ast$ and information about the magnetic energy of the AR.

The X1 flare and associated filament eruption occurring in NOAA Active Region 12017 on SOL2014-03-29 has been the source of intense study. In this work, we analyse the results of a series of non linear force free field extrapolations of the pre and post flare period of the flare. In combination with observational data provided by the IRIS, Hinode and SDO missions, we have confirmed the existence of two flux ropes present within the active region prior to flaring. Of these two flux ropes, we find that intriguingly only one erupts during the X1 flare. We propose that the reason for this is due to tether cutting reconnection allowing one of the flux ropes to rise to a torus unstable region prior to flaring, thus allowing it to erupt during the subsequent flare.

The triggering mechanism(s) and critical condition(s) of solar flares are still not completely clarified, although various studies have attempted to elucidate them. We have also proposed a theoretical flare-trigger model based on MHD simulations Kusano et al. 2012, in which two types of small-scale bipole field, the so-called Opposite Polarity (OP) and Reversed Shear (RS) types of field, can trigger flares. In this study, we evaluated the applicability of our flare-trigger model to observation of 32 flares that were observed by the Solar Dynamics Observatory (SDO), by focusing on geometrical structures. We classified the events into six types, including the OP and RS types, based on photospheric magnetic field configuration, presence of precursor brightenings, and shape of the initial flare ribbons. As a result, we found that approximately 30% of the flares were consistent with our flare-trigger model, and the number of RS type triggered flares is larger than that of the OP type. We found none of the sampled events contradicts our flare model, although we cannot clearly determine the trigger mechanism of 70% of the flares in this study. We carefully investigated the applicability of our flare-trigger model and the possibility that other models can explain the other 70% of the events. Consequently, we concluded that our flare-trigger model has certainly proposed important conditions for flare-triggering.

Aims. We report the small temporal variation of the axial dipole moment near the solar minimum and its application to the solar cycle prediction by the surface flux transport (SFT) model. Methods. We measure the axial dipole moment using the photospheric synoptic magnetogram observed by the Wilcox Solar Observatory (WSO), the ESA/NASA Solar and Heliospheric Observatory Michelson Doppler Imager (MDI), and the NASA Solar Dynamics Observatory Helioseismic and Magnetic Imager (HMI). We also use the surface flux transport model for the interpretation and prediction of the observed axial dipole moment. Results. We find that the observed axial dipole moment becomes approximately constant during the period of several years before each cycle minimum, which we call the axial dipole moment plateau. The cross-equatorial magnetic flux transport is found to be small during the period, although the significant number of sunspots are still emerging. The results indicates that the newly emerged magnetic flux does not contributes to the build up of the axial dipole moment near the end of each cycle. This is confirmed by showing that the time variation of the observed axial dipole moment agrees well with that predicted by the SFT model without introducing new emergence of magnetic flux. These results allows us to predict the axial dipole moment in Cycle 24/25 minimum using the SFT model without introducing new flux emergence. The predicted axial dipole moment of Cycle 24/25 minimum is 60--80 percent of Cycle 23/24 minimum, which suggests the amplitude of Cycle 25 even weaker than the current Cycle 24. Conclusions. The plateau of the solar axial dipole moment is an important feature for the longer prediction of the solar cycle based on the SFT model.

In order to understand the flare trigger mechanism, we conducted three-dimensional magnetohydrodynamic simulations using a coronal magnetic field model derived from data observed by the Hinode satellite. Several types of magnetic bipoles were imposed into the photospheric boundary of the Non-linear Force-Free Field (NLFFF) model of Active Region NOAA 10930 on 2006 December 13 to investigate what kind of magnetic disturbance may trigger the flare. As a result, we confirm that certain small bipole fields, which emerge into the highly sheared global magnetic field of an active region, can effectively trigger a flare. These bipole fields can be classified into two groups based on their orientation relative to the polarity inversion line: the so called opposite polarity (OP) and reversed shear (RS) structures as it was suggested by Kusano et al. (2012). We also investigated the structure of the footpoints of reconnected field lines. By comparing the distribution of reconstructed field lines and the observed flare ribbons, the trigger structure of the flare can be inferred. Our simulation suggests that the data-constrained simulation taking into account both the large-scale magnetic structure and the small-scale magnetic disturbance such as emerging fluxes is a good way to find out a flare productive active region for space weather prediction.

The stability of the magnetic field in the solar corona is important for understanding the causes of solar eruptions. Although various scenarios have been suggested to date, the tether-cutting reconnection scenario proposed by Moore et al.(2001) is one of the widely accepted models to explain the onset process of solar eruptions. Although the tether-cutting reconnection scenario proposed that sigmoidal field formed by the internal reconnection is the magnetic field in pre-eruptive state, the stability of the sigmoidal field has not yet been investigated quantitatively. In this paper, in order to elucidate the stability problem of pre-eruptive state, we developed a simple numerical analysis, in which the sigmoidal field is modeled by a double arc electric current loop and its stability is analyzed. As a result, we found that the double arc loop is more easily destabilized than the axisymmetric torus, and it becomes unstable even if the external field does not decay with altitude, which is in contrast to the axisymmetric torus instability. This suggests that the tether-cutting reconnection may well work as the onset mechanism of solar eruptions, and if so the critical condition for eruption under certain geometry may be determined by a new type of instability rather than the torus instability. Based on them, we propose a new type of instability called double arc instability (DAI). We discuss the critical conditions for DAI and derive a new parameter $\kappa$ defined as the product of the magnetic twist and the normalized flux of tether-cutting reconnection.

The physical properties and its contribution to the onset of solar flare are still unclear although chromospheric brightening is considered a precursor phenomenon of flare. Many studies suggested that photospheric magnetic field changes cause destabilization of large-scale coronal structure. We aim to understand how a small photospheric change contributes to a flare and to reveal how the intermediary chromosphere behaves in the precursor phase. We analyzed the precursor brightening of the X1.6 flare on 2014 October 22 in the AR 12192 using the Interface Region Imaging Spectrograph (IRIS) and Hinode/EUV Imaging Spectrometer (EIS) data. We investigated a localized jet with the strong precursor brightening, and compared the intensity, Doppler velocity, and line width in C II, Mg II k, Si IV lines by IRIS and He II, Fe XII, Fe XV lines by Hinode/EIS. We also analyzed photospheric magnetic field and chromospheric/coronal structures using Solar Dynamics Observatory (SDO)/Helioseismic and Magnetic Imager (HMI) and Atmospheric Imaging Assembly (AIA). We found a significant blueshift (~ 100 km/s), which is related to the strong precursor brightening over a characteristic magnetic field structure, and the blueshift was observed at all the temperature. This might indicate that the flow is accelerated by Lorentz force. Moreover, the large-scale coronal loop that connects the foot-points of the flare ribbons was destabilized just after the precursor brightening with the blueshift. It suggests that magnetic reconnection locally occurred in the lower chromosphere and it triggered magnetic reconnection of the X1.6 flare in the corona.

The solar magnetic field in a flare-producing active region (AR) is much more complicated than theoretical models, which assume a very simple magnetic field structure. The X1.0 flare, which occurred in AR 12192 on 2014 October 25, showed a complicated three-ribbon structure. To clarify the trigger process of the flare and to evaluate the applicability of a simple theoretical model, we analyzed the data from Hinode/Solar Optical Telescope and the Solar Dynamics Observatory/Helioseismic and Magnetic Imager, Atmospheric Imaging Assembly. We investigated the spatio-temporal correlation between the magnetic field structures, especially the non-potentiality of the horizontal field, and the bright structures in the solar atmosphere. As a result, we determined that the western side of the positive polarity, which is intruding on a negative polarity region, is the location where the flare was triggered. This is due to the fact that the sign of the magnetic shear in that region was opposite that of the major shear of the AR, and the significant brightenings were observed over the polarity inversion line (PIL) in that region before flare onset. These features are consistent with the recently proposed flare-trigger model that suggests that small reversed shear (RS) magnetic disturbances can trigger solar flares. Moreover, we found that the RS field was located slightly off the flaring PIL, contrary to the theoretical prediction. We discuss the possibility of an extension of the RS model based on an extra numerical simulation. Our result suggests that the RS field has a certain flexibility for displacement from a highly sheared PIL, and that the RS field triggers more flares than we expected.

Solar flares are generally believed to be powered by free magnetic energy stored in the corona, but the build up of coronal energy alone may be insufficient for the imminent flare occurrence. The flare onset mechanism is a critical but less understood problem, insights into which could be gained from small-scale energy releases known as precursors, which are observed as small pre-flare brightenings in various wavelengths, and also from certain small-scale magnetic configurations such as the opposite polarity fluxes, where magnetic orientation of small bipoles is opposite to that of the ambient main polarities. However, high-resolution observations of flare precursors together with the associated photospheric magnetic field dynamics are lacking. Here we study precursors of a flare using unprecedented spatiotemporal resolution of the 1.6 m New Solar Telescope, complemented by novel microwave data. Two episodes of precursor brightenings are initiated at a small-scale magnetic channel (a form of opposite polarity fluxes) with multiple polarity inversions and enhanced magnetic fluxes and currents, lying near the footpoints of sheared magnetic loops. The low-atmospheric origin of these precursor emissions is corroborated by microwave spectra. We propose that the emerging magnetic channel field interacts with the sheared arcades to cause precursor brightenings at the main flare core region. These high-resolution results provide evidence of low-atmospheric small-scale energy release and possible relationship to the onset of the main flare.

Observations of the sun suggest that solar activities systematically create north-south hemispheric asymmetries. For instance, the hemisphere in which the sunspot activity is more active tends to switch after the early half of each solar cycle. Svalgaard & Kamide (2013) recently pointed out that the time gaps of polar field reversal between the north and south hemispheres are simply consequences of the asymmetry of sunspot activity. However, the mechanism underlying the asymmetric feature in solar cycle activities is not yet well understood. In this paper, in order to explain the cause of the asymmetry from the theoretical point of view, we investigate the relationship between the dipole- and quadrupole-type components of the magnetic field in the solar cycle using the mean-field theory based on the flux transport dynamo model. As a result, we found that there are two different attractors of the solar cycle, in which either the north or the south polar field is first reversed, and that the flux transport dynamo model well explains the phase-asymmetry of sunspot activity and the polar field reversal without any ad hoc source of asymmetry.

We analyze a three-dimensional (3D) magnetic structure and its stability in large solar active region(AR) 12192, using the 3D coronal magnetic field constructed under a nonlinear force-free field (NLFFF) approximation. In particular, we focus on the magnetic structure that produced an X3.1-class flare which is one of the X-class flares observed in AR 12192. According to our analysis, the AR contains multiple-flux-tube system, {\it e.g.}, a large flux tube, both of whose footpoints are anchored to the large bipole field, under which other tubes exist close to a polarity inversion line (PIL). These various flux tubes of different sizes and shapes coexist there. In particular, the later are embedded along the PIL, which produces a favorable shape for the tether-cutting reconnection and is related to the X-class solar flare. We further found that most of magnetic twists are not released even after the flare, which is consistent with the fact that no observational evidence for major eruptions was found. On the other hand, the upper part of the flux tube is beyond a critical decay index, essential for the excitation of torus instability before the flare, even though no coronal mass ejections (CMEs) were observed. We discuss the stability of the complicated flux tube system and suggest the reason for the existence of the stable flux tube. In addition, we further point out a possibility for tracing the shape of flare ribbons, on the basis of a detailed structural analysis of the NLFFF before a flare.

Flares are a major explosive event in our solar system. They are often followed by coronal mass ejection that has a potential to trigger the geomagnetic storms. There are various studies aiming to predict when and where the flares are likely to occur. Most of these studies mainly discuss the photospheric and chromospheric activity before the flare onset. In this paper we study the coronal features before the famous large flare occurrence on December 13th, 2006. Using the data from Hinode/EUV Imaging Spectrometer (EIS), X-Ray Telescope (XRT), and Solar and Heliospheric Observatory (SOHO) /Extreme ultraviolet Imaging Telescope (EIT), we discuss the coronal features in the large scale (~ a few 100 arcsec) before the flare onset. Our findings are as follows: 1) The upflows in and around active region start growing from ~10 to 30 km /s a day before the flare. 2) The expanding coronal loops are clearly observed a few hours before the flare. 3) Soft X-ray and EUV intensity are gradually reduced. 4) The upflows are further enhanced after the flare. From these observed signatures, we conclude that the outer part of active region loops with low density were expanding a day before the flare onset, and the inner part with high density were expanding a few hours before the onset.

To elucidate the flare trigger mechanism, we have analyzed several flare events which were observed by Hinode/Solar Optical Telescope (SOT), in our previous study. Because of the limitation of SOT field of view, however, only four events in the Hinode data sets have been utilizable. Therefore, increasing the number of events is required for evaluating the flare trigger models. We investigated the applicability of data obtained by the Solar Dynamics Observatory (SDO) to increase the data sample for a statistical analysis of the flare trigger process. SDO regularly observes the full disk of the sun and all flares although its spatial resolution is lower than that of Hinode. We investigated the M6.6 flare which occurred on 13 February 2011 and compared the analyzed data of SDO with the results of our previous study using Hinode/SOT data. Filter and vector magnetograms obtained by the Helioseismic and Magnetic Imager (HMI) and filtergrams from the Atmospheric Imaging Assembly (AIA) 1600A were employed. From the comparison of small-scale magnetic configurations and chromospheric emission prior to the flare onset, we confirmed that the trigger region is detectable with the SDO data. We also measured the magnetic shear angles of the active region and the azimuth and strength of the flare-trigger field. The results were consistent with our previous study. We concluded that statistical studies of the flare trigger process are feasible with SDO as well as Hinode data. We also investigated the temporal evolution of the magnetic field before the flare onset with SDO.

We present a comparison of the Solar Dynamics Observatory (SDO) analysis of NOAA Active Region (AR) 11158 and numerical simulations of flux-tube emergence, aiming to investigate the formation process of this flare-productive AR. First, we use SDO/Helioseismic and Magnetic Imager (HMI) magnetograms to investigate the photospheric evolution and Atmospheric Imaging Assembly (AIA) data to analyze the relevant coronal structures. Key features of this quadrupolar region are a long sheared polarity inversion line (PIL) in the central delta-sunspots and a coronal arcade above the PIL. We find that these features are responsible for the production of intense flares, including an X2.2-class event. Based on the observations, we then propose two possible models for the creation of AR 11158 and conduct flux-emergence simulations of the two cases to reproduce this AR. Case 1 is the emergence of a single flux tube, which is split into two in the convection zone and emerges at two locations, while Case 2 is the emergence of two isolated but neighboring tubes. We find that, in Case 1, a sheared PIL and a coronal arcade are created in the middle of the region, which agrees with the AR 11158 observation. However, Case 2 never builds a clear PIL, which deviates from the observation. Therefore, we conclude that the flare-productive AR 11158 is, between the two cases, more likely to be created from a single split emerging flux than from two independent flux bundles.

We present a 2.5D MHD simulation of a magnetic flux rope (FR) propagating in the heliosphere and investigate the cause of the observed sharp plasma beta transition. Specifically, we consider a strong internal magnetic field and an explosive fast start, such that the plasma beta is significantly lower in the FR than the sheath region that is formed ahead. This leads to an unusual FR morphology in the first stage of propagation, while the more traditional view (e.g. from space weather simulations like Enlil) of a `pancake' shaped FR is observed as it approaches 1 AU. We investigate how an equipartition line, defined by a magnetic Weber number, surrounding a core region of a propagating FR can demarcate a boundary layer where there is a sharp transition in the plasma beta. The substructure affects the distribution of toroidal flux, with the majority of the flux remaining in a small core region which maintains a quasi-cylindrical structure. Quantitatively, we investigate a locus of points where the kinetic energy density of the relative inflow field is equal to the energy density of the transverse magnetic field (i.e. effective tension force). The simulation provides compelling evidence that at all heliocentric distances the distribution of toroidal magnetic flux away from the FR axis is not linear; with 80% of the toroidal flux occurring within 40% of the distance from the FR axis. Thus our simulation displays evidence that the competing ideas of a pancaking structure observed remotely can coexist with a quasi-cylindrical magnetic structure seen in situ.

We develop a nonlinear force-free field (NLFFF) extrapolation code based on the magnetohydrodynamic (MHD) relaxation method. We extend the classical MHD relaxation method in two important ways. First, we introduce an algorithm initially proposed by cite{2002JCoPh.175..645D} to effectively clean the numerical errors associated with $nabla cdot vec{B}$. Second, the multi-grid type method is implemented in our NLFFF to perform direct analysis of the high-resolution magnetogram data. As a result of these two implementations, we successfully extrapolated the high resolution force-free field introduced by cite{1990ApJ...352..343L} with better accuracy in a drastically shorter time. We also applied our extrapolation method to the MHD solution obtained from the flux-emergence simulation by cite{2012ApJ...748...53M}. We found that NLFFF extrapolation may be less effective for reproducing areas higher than a half-domain, where some magnetic loops are found in a state of continuous upward expansion. However, an inverse S shaped structure consisting of the sheared and twisted loops formed in the lower region can be captured well through our NLFFF extrapolation method. We further discuss how well these sheared and twisted fields are reconstructed by estimating the magnetic topology and twist quantitatively.

We investigated four major solar flare events that occurred in active regions NOAA 10930 (December 13 and 14, 2006) and NOAA 11158 (February 13 and 15, 2011) by using data observed by the Solar Optical Telescope (SOT) onboard the Hinode satellite. To reveal the trigger mechanism of solar flares, we analyzed the spatio-temporal correlation between the detailed magnetic field structure and the emission image of the Ca H line at the central part of flaring regions for several hours prior to the onset of flares. We observed in all the flare events that the magnetic shear angle in the flaring regions exceeded 70 degrees, as well as that characteristic magnetic disturbances developed at the centers of flaring regions in the pre-flare phase. These magnetic disturbances can be classified into two groups depending on the structure of their magnetic polarity inversion lines; the so-called "Opposite-Polarity" and "Reversed-Shear" magnetic field recently proposed by our group, although the magnetic disturbance in one event of the four samples is too subtle to clearly recognize the detailed structure. The result suggests that some major solar flares are triggered by rather small magnetic disturbances. We also show that the critical size of the flare-trigger field varies among flare events and briefly discuss how the flare-trigger process depends on the evolution of active regions.

To better understand a preferred magnetic field configuration and its evolution during Coronal Mass Ejection events, we investigated the spatial and temporal evolution of photospheric magnetic fields in the active region NOAA 9236 that produced eight flare-associated CMEs during the time period of 2000 November 23-26. The time variations of the total magnetic helicity injection rate and the total unsigned magnetic flux are determined and examined not only in the entire active region but also in some local regions such as the main sunspots and the CME-associated flaring regions using SOHO/MDI magnetogram data. As a result, we found that: (1) in the sunspots, a large amount of postive (right-handed) magnetic helicity was injected during most of the examined time period, (2) in the flare region, there was a continuous injection of negative (left-handed) magnetic helicity during the entire period, accompanied by a large increase of the unsigned magnetic flux, and (3) the flaring regions were mainly composed of emerging bipoles of magnetic fragments in which magnetic field lines have substantially favorable conditions for making reconnection with large-scale, overlying, and oppositely directed magnetic field lines connecting the main sunspots. These observational findings can also be well explained by some MHD numerical simulations for CME initiation (e.g., reconnection-favored emerging flux models). We therefore conclude that reconnection-favored magnetic fields in the flaring emerging flux regions play a crucial role in producing the multiple flare-associated CMEs in NOAA 9236.

We report a detailed event analysis on the M6.6-class flare in the active region (AR) NOAA 11158 on 2011 February 13. AR 11158, which consisted of two major emerging bipoles, showed prominent activities including one X- and several M-class flares. In order to investigate the magnetic structures related to the M6.6 event, particularly the formation process of a flare-triggering magnetic region, we analyzed multiple spacecraft observations and numerical results of a flare simulation. We observed that, in the center of this quadrupolar AR, a highly sheared polarity inversion line (PIL) was formed through proper motions of the major magnetic elements, which built a sheared coronal arcade lying over the PIL. The observations lend support to the interpretation that the target flare was triggered by a localized magnetic region that had an intrusive structure, namely a positive polarity penetrating into a negative counterpart. The geometrical relationship between the sheared coronal arcade and the triggering region was consistent with the theoretical flare model based on the previous numerical study. We found that the formation of the trigger region was due to a continuous accumulation of the small-scale magnetic patches. A few hours before the flare occurrence, the series of emerged/advected patches reconnected with a preexisting fields. Finally, the abrupt flare eruption of the M6.6 event started around 17:30 UT. Our analysis suggests that, in a triggering process of a flare activity, all magnetic systems of multiple scales, not only the entire AR evolution but also the fine magnetic elements, are altogether involved.

Solar flares and coronal mass ejections (CMEs), the most catastrophic eruptions in our solar system, have been known to affect terrestrial environments and infrastructure. However, because their triggering mechanism is still not sufficiently understood, our capacity to predict the occurrence of solar eruptions and to forecast space weather is substantially hindered. Even though various models have been proposed to determine the onset of solar eruptions, the types of magnetic structures capable of triggering these eruptions are still unclear. In this study, we solved this problem by systematically surveying the nonlinear dynamics caused by a wide variety of magnetic structures in terms of three-dimensional magnetohydrodynamic simulations. As a result, we determined that two different types of small magnetic structures favor the onset of solar eruptions. These structures, which should appear near the magnetic polarity inversion line (PIL), include magnetic fluxes reversed to the potential component or the nonpotential component of major field on the PIL. In addition, we analyzed two large flares, the X-class flare on December 13, 2006 and the M-class flare on February 13, 2011, using imaging data provided by the Hinode satellite, and we demonstrated that they conform to the simulation predictions. These results suggest that forecasting of solar eruptions is possible with sophisticated observation of a solar magnetic field, although the lead time must be limited by the time scale of changes in the small magnetic structures.

We perform four numerical magnetohydrodynamic simulations in 2.5 dimensions (2.5D) of fast Coronal Mass Ejections (CMEs) and their associated shock fronts between 10Rs and 300Rs. We investigate the relative change in the shock standoff distance, Sd, as a fraction of the CME radial half-width, Dob (i.e. Sd/Dob). Previous hydrodynamic studies have related the shock standoff distance for Earths magnetosphere to the density compression ratio (DR,Ru/Rd) measured across the bow shock (Spreiter, Summers and Alksne 1966). The DR coefficient, kdr, which is the proportionality constant between the relative standoff distance (Sd/Dob) and the compression ratio, was semi-empirically estimated as 1.1. For CMEs, we show that this value varies linearly as a function of heliocentric distance and changes significantly for different radii of curvature of the CMEs leading edge. We find that a value of 0.8+-0.1 is more appropriate for small heliocentric distances (<30Rs) which corresponds to the spherical geometry of a magnetosphere presented by Seiff (1962). As the CME propagates its cross section becomes more oblate and the kdr value increases linearly with heliocentric distance, such that kdr= 1.1 is most appropriate at a heliocentric distance of about 80Rs. For terrestrial distances (215Rs) we estimate kdr= 1.8+-0.3, which also indicates that the CME cross-sectional structure is generally more oblate than that of Earths magnetosphere. These alterations to the proportionality coefficients may serve to improve investigations into the estimates of the magnetic field in the corona upstream of a CME as well as the aspect ratio of CMEs as measured in situ.

We use autocorrelation to investigate evolution in flow fields inferred by applying Fourier Local Correlation Tracking (FLCT) to a sequence of high-resolution (0.3 \arcsec), high-cadence ($\simeq 2$ min) line-of-sight magnetograms of NOAA active region (AR) 10930 recorded by the Narrowband Filter Imager (NFI) of the Solar Optical Telescope (SOT) aboard the {\em Hinode} satellite over 12--13 December 2006. To baseline the timescales of flow evolution, we also autocorrelated the magnetograms, at several spatial binnings, to characterize the lifetimes of active region magnetic structures versus spatial scale. Autocorrelation of flow maps can be used to optimize tracking parameters, to understand tracking algorithms' susceptibility to noise, and to estimate flow lifetimes. Tracking parameters varied include: time interval $\Delta t$ between magnetogram pairs tracked, spatial binning applied to the magnetograms, and windowing parameter $\sigma$ used in FLCT. Flow structures vary over a range of spatial and temporal scales (including unresolved scales), so tracked flows represent a local average of the flow over a particular range of space and time. We define flow lifetime to be the flow decorrelation time, $\tau$. For $\Delta t > \tau$, tracking results represent the average velocity over one or more flow lifetimes. We analyze lifetimes of flow components, divergences, and curls as functions of magnetic field strength and spatial scale. We find a significant trend of increasing lifetimes of flow components, divergences, and curls with field strength, consistent with Lorentz forces partially governing flows in the active photosphere, as well as strong trends of increasing flow lifetime and decreasing magnitudes with increases in both spatial scale and $\Delta t$.

Using a model of particle acceleration and transport in solar flares, we investigate the height distribution of coronal electrons by focusing on the energy-dependent pitch-angle scattering. When pitch-angle scattering is not included, the peak heights of loop-top electrons are constant, regardless of their energy, owing to the continuous acceleration and compression of the electrons via shrinkage of magnetic loops. On the other hand, under pitch-angle scattering, the electron heights are energy dependent; intermediate energy electrons are at a higher altitude, whereas lower and higher energy electrons are at lower altitudes. This implies that the intermediate energy electrons are inhibited to follow the shrinking field lines to lower altitudes because pitch-angle scattering causes efficient precipitation of these electrons into the footpoint and their subsequent loss from the loop. This result is qualitatively consistent with the position of the above-the-loop-top hard X-ray (HXR) source that is located above coronal HXR loops emitted by lower energy electrons and microwaves emitted by higher energy electrons. Quantitative agreement with observations might be achieved by considering primary acceleration before the onset of loop shrinkage and additional pitch-angle scattering via wave-particle interactions.

We performed magnetohydrodynamic simulation of a formation process of coronal mass ejections (CMEs), focusing on interaction (reconnection) between an ejecting flux rope and its ambient field. We examined three cases with different ambient fields: no ambient field, and cases with dipole field of two opposite directions which are parallel and anti-parallel to that of the flux rope surface. As a result, while the flux rope disappears in the anti-parallel case, in other cases the flux ropes can evolve to CMEs and show different amounts of rotation of the flux rope. The results imply that the interaction between an ejecting flux rope and its ambient field is an important process for determining CME formation and CME orientation, and also show that the amount and direction of magnetic flux within the flux rope and the ambient field are key parameters for CME formation. Especially, the interaction (reconnection) plays a significant role to the rotation of the flux rope, with a process similar to "tilting instability" in a spheromak-type experiment of laboratory plasma.

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