We report on results of imaging and spectral studies of X-ray emission from Jupiter observed by Suzaku. In 2006 Suzaku had found diffuse X-ray emission in $1\unicode{x2013}5$ keV associated with Jovian inner radiation belts. It has been suggested that the emission is caused by the inverse-Compton scattering by ultra-relativistic electrons ($ \sim 50 $ MeV) in Jupiter's magnetosphere. To confirm the existence of this emission and to understand its relation to the solar activity, we conducted an additional Suzaku observation in 2014 around the maximum of the 24th solar cycle. As a result, we successfully found again the diffuse emission around Jupiter in $1\unicode{x2013}5$ keV and also point-like emission in $0.4\unicode{x2013}1$ keV. The luminosity of the point-like emission which was probably composed of solar X-ray scattering, charge exchange, or auroral bremsstrahlung emission increased by a factor of $ \sim 5$ with respect to 2006, most likely due to an increase of the solar activity. The diffuse emission spectrum in the $1\unicode{x2013}5$ keV band was well-fitted with a flat power-law function ($ \Gamma = 1.4 \pm 0.1 $) as in the past observation, which supported the inverse-Compton scattering hypothesis. However, its spatial distribution changed from $ \sim 12 \times 4 $ Jovian radius (Rj) to $ \sim 20 \times 7 $ Rj. The luminosity of the diffuse emission increased by a smaller factor of $ \sim 3 $. This indicates that the diffuse emission is not simply responding to the solar activity, which is also known to cause little effect on the distribution of high-energy electrons around Jupiter. Further sensitive study of the spatial and spectral distributions of the diffuse hard X-ray emission is important to understand how high-energy particles are accelerated in Jupiter's magnetosphere.

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.

On SOL2017-09-06 solar active region 12673 produced an X9.3 flare which is regarded as largest to occur in solar cycle 24. In this work we have preformed a magnetohydrodynamic (MHD) simulation in order to reveal the three-dimensional (3D) dynamics of the magnetic fields associated with the X9.3 solar flare. We first performed an extrapolation of the 3D magnetic field based on the observed photospheric magnetic field prior to the flare and then used it as the initial condition for an MHD simulation. Consequently, the simulation showed a dramatic eruption. In particular, we found that a large coherent flux rope composed of highly twisted magnetic field lines is formed during the eruption. A series of small flux ropes are found to lie along a magnetic polarity inversion line prior to the flare. Reconnection occurring between each small flux rope during the early stages of the eruption forms the large and highly twisted flux rope.Furthermore, we found a writhing motion of the erupting flux rope. The understanding of these dynamics is important in increasing the accuracy of space weather forecasting. We report on the detailed dynamics of the 3D eruptive flux rope and discuss the possible mechanisms of the writhing motion.

A flare kernel associated with a C4 class flare was observed in a spectral window including the He I triplet 10830 A and Si I 10827 A with a spectropolarimeter on the Domeless Solar Telescope at Hida Observatory on August 9th, 2015. Observed Stokes profiles of the He I triplet in the flare kernel in its post-maximum phase are well reproduced through inversions considering the Zeeman and the Paschen-Back effects with a three slab model of the flare kernel, in which two slabs having up and downward velocities produce emissions and one slab produces an absorption. The magnetic field strength inferred from the emission components of the He I line is 1400 G, which is significantly stronger than 690 G that is observed at the same location in the same line 6.5 hours before the flare. In addition, photospheric magnetic field vector derived from the Si I10827 A is similar to that of the flare kernel. To explain this results, we suggest that the emission in the He I triplet during the flare is produced in the deep layer, around which bombardment of non-thermal electrons leads to the formation of a coronal temperature plasma. Assuming a hydrogen column density at the location where the He I emissions are formed, and a power-law index of non-thermal electron energy distribution, we derived the low-energy cutoff of the non-thermal electron as 20 - 30 keV, which is consistent with that inferred from hard X-ray data obtained by RHESSI.

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

We estimated the accuracy of coronal magnetic fields derived from radio observations by comparing them to potential field calculations and the DEM measurements using EUV observations. We derived line of sight component of the coronal magnetic field from polarization observations of the thermal bremsstrahlung in the NOAA active region 11150, observed around 3:00 UT on February 3, 2011 using the Nobeyama Radioheliograph at 17 GHz. Because the thermal bremsstrahlung intensity at 17 GHz includes both chromospheric and coronal components, we extracted only the coronal component by measuring the coronal emission measure in EUV observations. In addition, we derived only the radio polarization component of the corona by selecting the region of coronal loops and weak magnetic field strength in the chromosphere along the line of sight. The upper limit of the coronal longitudinal magnetic fields were determined as 100 - 210 G. We also calculated the coronal longitudinal magnetic fields from the potential field extrapolation using the photospheric magnetic field obtained from the Helioseismic and Magnetic Imager (HMI). However, the calculated potential fields were certainly smaller than the observed coronal longitudinal magnetic field. This discrepancy between the potential and the observed magnetic field strengths can be explained consistently by two reasons; (a) the underestimation of the coronal emission measure resulting from the limitation of the temperature range of the EUV observations, (b) the underestimation of the coronal magnetic field resulting from the potential field assumption.

We investigated the role of photospheric plasma motions in the formation and evolution of polar magnetic patches using time-sequence observations with high spatial resolution. The observations were obtained with the spectropolarimeter on board the Hinode satellite. From the statistical analysis using 75 magnetic patches, we found that they are surrounded by strong converging, supergranulation associated flows during their apparent life time and that the converging flow around the patch boundary is better observed in the Doppler velocity profile in the deeper photosphere. Based on our analysis we suggest that the like-polarity magnetic fragments in the polar region are advected and clustered by photospheric converging flows thereby resulting in the formation of polar magnetic patches. Our observations show that, in addition to direct cancellation magnetic patches decay by fragmentation followed by unipolar disappearance or unipolar disappearance without fragmentation. It is possible that the magnetic patches of existing polarity fragment or diffuse away into smaller elements and eventually cancel out with opposite polarity fragments that reach the polar region around solar cycle maximum. This could be one of the possible mechanisms by which the existing polarity decay during the reversal of the polar magnetic field.

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.

The magnetic properties of the Sun's polar faculae are investigated with spectropolarimetric observations of the north polar region obtained by the Hinode satellite in 2007 September. Polar faculae are embedded in nearly all magnetic patches with fluxes greater than $10^{18}$ Mx, while magnetic patches without polar faculae dominate in the flux range below $10^{18}$ Mx. The faculae are considerably smaller than their parent patches, and single magnetic patches contain single or multiple faculae. The faculae in general have higher intrinsic magnetic field strengths than the surrounding regions within their parent patches. Less than 20% of the total magnetic flux contributed by the large (${\ge}10^{18}$ Mx) concentrations, which are known to be modulated by the solar cycle, is accounted for by the associated polar faculae.

We study the three-dimensional magnetic structure of solar active region 11158, which produced one X-class and several M-class flares on 2011 February 13$-$16. We focus on the magnetic twist in four flare events, M6.6, X2.2, M1.0, and M1.1. The magnetic twist is estimated from the nonlinear force-free field extrapolated from the vector fields obtained from the Helioseismic and Magnetic Imager on board the Solar Dynamic Observatory using magnetohydrodynamic relaxation method developed by \cite{2011ApJ...738..161I}. We found that strongly twisted lines ranging from half-turn to one-turn twist were built up just before the M6.6- and X2.2 flares and disappeared after that. Because most of the twist remaining after these flares was less than half-turn twist, this result suggests that the buildup of magnetic twist over the half-turn twist is a key process in the production of large flares. On the other hand, even though these strong twists were also built up just before the M1.0 and M1.1 flares, most of them remained afterwords. Careful topological analysis before the M1.0 and M1.1 flares shows that the strongly twisted lines were surrounded mostly by the weakly twisted lines formed in accordance with the clockwise motion of the positive sunspot, whose footpoints are rooted in strong magnetic flux regions. These results imply that these weakly twisted lines might suppress the activity of the strongly twisted lines in the last two M-class flares.

We analyze the temporal evolution of the three-dimensional (3D) magnetic structure of the flaring active region (AR) NOAA 10930 by using the nonlinear force-free fields extrapolated from the photospheric vector magnetic fields observed by the Solar Optical Telescope on board {\it Hinode}. This AR consisted mainly of two types of twisted magnetic field lines: One has a strong negative (left-handed) twist due to the counterclockwise motion of the positive sunspot and is rooted in the regions of both polarities in the sunspot at a considerable distance from the polarity inversion line (PIL). In the flare phase, dramatic magnetic reconnection occurs in those negatively twisted lines in which the absolute value of the twist is greater than a half-turn. The other type consists of both positively and negatively twisted field lines formed relatively close to the PIL between two sunspots. A strong CaII image began to brighten in this region of mixed polarity, in which the positively twisted field lines were found to be injected within one day across the pre-existing negatively twisted region, along which strong currents were embedded. Consequently, the central region near the PIL distributed with a mix of differently twisted field lines and the strong currents may play a prominent role in flare onset.

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 have been monitoring yearly variation in the Sun's polar magnetic fields with the Solar Optical Telescope aboard {\it Hinode} to record their evolution and expected reversal near the solar maximum. All magnetic patches in the magnetic flux maps are automatically identified to obtain the number density and magnetic flux density as a function of th total magnetic flux per patch. The detected magnetic flux per patch ranges over four orders of magnitude ($10^{15}$ -- $10^{20}$ Mx). The higher end of the magnetic flux in the polar regions is about one order of magnitude larger than that of the quiet Sun, and nearly that of pores. Almost all large patches ($ \geq 10^{18}$ Mx) have the same polarity, while smaller patches have a fair balance of both polarities. The polarity of the polar region as a whole is consequently determined only by the large magnetic concentrations. A clear decrease in the net flux of the polar region is detected in the slow rising phase of the current solar cycle. The decrease is more rapid in the north polar region than in the south. The decrease in the net flux is caused by a decrease in the number and size of the large flux concentrations as well as the appearance of patches with opposite polarity at lower latitudes. In contrast, we do not see temporal change in the magnetic flux associated with the smaller patches ($ < 10^{18}$ Mx) and that of the horizontal magnetic fields during the years 2008--2012.

We report on the first simultaneous observation of an H-alpha Moreton wave, the corresponding EUV fast coronal waves, and a slow and bright EUV wave (typical EIT wave). Associated with an X6.9 flare that occurred on 2011 August 9 at the active region NOAA 11263, we observed a Moreton wave in the H-alpha images taken by the Solar Magnetic Activity Research Telescope (SMART) at Hida Observatory of Kyoto University. In the EUV images obtained by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) we found not only the corresponding EUV fast "bright" coronal wave, but also the EUV fast "faint" wave that is not associated with the H-alpha Moreton wave. We also found a slow EUV wave, which corresponds to a typical EIT wave. Furthermore, we observed, for the first time, the oscillations of a prominence and a filament, simultaneously, both in the H-alpha and EUV images. To trigger the oscillations by the flare-associated coronal disturbance, we expect a coronal wave as fast as the fast-mode MHD wave with the velocity of about 570 - 800 km/s. These velocities are consistent with those of the observed Moreton wave and the EUV fast coronal wave.

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.

Observations of the polar region of the Sun are critically important for understanding the solar dynamo and the acceleration of solar wind. We carried out precise magnetic observations on both the North polar region and the quiet Sun at the East limb with the Spectro-Polarimeter of the Solar Optical Telescope aboard Hinode to characterize the polar region with respect to the quiet Sun. The average area and the total magnetic flux of the kG magnetic concentrations in the polar region appear to be larger than those of the quiet Sun. The magnetic field vectors classified as vertical in the quiet Sun have symmetric histograms around zero in the strengths, showing balanced positive and negative flux, while the histogram in the North polar region is clearly asymmetric, showing a predominance of the negative polarity. The total magnetic flux of the polar region is larger than that of the quiet Sun. In contrast, the histogram of the horizontal magnetic fields is exactly the same between the polar region and the quiet Sun. This is consistent with the idea that a local dynamo process is responsible for the horizontal magnetic fields. A high-resolution potential field extrapolation shows that the majority of magnetic field lines from the kG-patches in the polar region are open with a fanning-out structure very low in the atmosphere, while in the quiet Sun, almost all the field lines are closed.

Alfven waves created by sub-photospheric motions or by magnetic reconnection in the low solar atmosphere seem good candidates for coronal heating. However, the corona is also likely to be heated more directly by magnetic reconnection, with dissipation taking place in current sheets. Distinguishing observationally between these two heating mechanisms is an extremely difficult task. We perform 1.5-dimensional MHD simulations of a coronal loop subject to each type of heating and derive observational quantities that may allow these to be differentiated.

The plasmoid-induced-reconnection model explaining solar flares based on bursty reconnection produced by an ejecting plasmoid suggests a possible relation between the ejection velocity of a plasmoid and the rate of magnetic reconnection. In this study, we focus on the quantitative description of this relation. We performed magnetohydrodynamic (MHD) simulations of solar flares by changing the values of resistivity and the plasmoid velocity. The plasmoid velocity has been changed by applying an additional force to the plasmoid to see how the plasmoid velocity affects the reconnection rate. An important result is that the reconnection rate has a positive correlation with the plasmoid velocity, which is consistent with the plasmoid-induced-reconnection model for solar flares. We also discuss an observational result supporting this positive correlation.

We performed magnetohydrodynamic simulation of coronal mass ejections (CMEs) and associated giant arcade formations, and the results suggested new interpretations of observations of CMEs. We performed two cases of the simulation: with and without heat conduction. Comparing between the results of the two cases, we found that reconnection rate in the conductive case is a little higher than that in the adiabatic case and the temperature of the loop top is consistent with the theoretical value predicted by the Yokoyama-Shibata scaling law. The dynamical properties such as velocity and magnetic fields are similar in the two cases, whereas thermal properties such as temperature and density are very different.In both cases, slow shocks associated with magnetic reconnectionpropagate from the reconnection region along the magnetic field lines around the flux rope, and the shock fronts form spiral patterns. Just outside the slow shocks, the plasma density decreased a great deal. The soft X-ray images synthesized from the numerical results are compared with the soft X-ray images of a giant arcade observed with the Soft X-ray Telescope aboard {\it Yohkoh}, it is confirmed that the effect of heat conduction is significant for the detailed comparison between simulation and observation. The comparison between synthesized and observed soft X-ray images provides new interpretations of various features associated with CMEs and giant arcades.