Gas convection is observed in the solar photosphere as the granulation, i.e., having highly time-dependent cellular patterns, consisting of numerous bright cells called granules and dark surrounding-channels called intergranular lanes. Many efforts have been made to characterize the granulation, which may be used as an energy source for various types of dynamical phenomena. Although the horizontal gas flow dynamics in intergranular lanes may play a vital role, but they are poorly understood. This is because the Doppler signals can be obtained only at the solar limb, where the signals are severely degraded by a foreshortening effect. To reduce such a degradation, we use Hinode's spectroscopic data, which are free from a seeing-induced image degradation, and improve its image quality by correcting for straylight in the instruments. The dataset continuously covers from the solar disk to the limb, providing a multidirectional line-of-sight (LOS) diagnosis against the granulation. The obtained LOS flow-field variation across the disk indicates a horizontal flow speed of 1.8-2.4 km/s. We also derive the spatial distribution of the horizontal flow speed, which is 1.6 km/s in granules and 1.8 km/s in intergranular lanes, and where the maximum speed is inside intergranular lanes. This result newly suggests the following sequence of horizontal flow: A hot rising gas parcel is strongly accelerated from the granular center, even beyond the transition from the granules to the intergranular lanes, resulting in the fastest speed inside the intergranular lanes, and the gas may also experience decelerations in the intergranular lane.

Statistical dependencies among features of coronal mass ejections (CMEs), solar flares, and sigmoidal structures in soft-X-ray images were investigated. We applied analysis methods to all the features in the same way in order to investigate the reproducibility of the correlations among them, which may be found from the combination of previous statistical studies. The samples of 211 M-class and X-class flares, which were observed between 2006 and 2015 by Hinode/X-ray telescope, Solar and Heliospheric Observatory/Large Angle and Spectrometric Coronagraph, and GOES, were examined statistically. Five kinds of analysis were performed: Occurrence rate analysis, linear-correlation analysis, association analysis, the Kolmogorov--Smirnov test, and Anderson-Darling test. The analyses show three main results. First, the sigmoidal structure and long duration events (LDEs) has stronger dependency on the CME occurrence than large X-ray class events in on-disk events. Second, for the limb events, the significant dependency exists between LDEs and CME occurrence, and between X-ray class and CME occurrence. Third, there existed 32.4% of on-disk flare events, which had sigmoidal structure and were not accompanied by CMEs. However, the occurrence probability of CMEs without sigmoidal structures is very small, 8.8 %, in on-disk events. While the first and second results are consistent with previous studies, we newly provided the difference between the on-disk events and limb events. The third result that non-sigmoidal regions produce less eruptive events is also different from previous results. We suggest that sigmoidal structures in soft X-ray images will be a helpful feature for CME prediction regarding on-disk flare events.

Solar granules are bright patterns surrounded by dark channels called intergranular lanes in the solar photosphere and are a manifestation of overshooting convection. Observational studies generally find stronger upflows in granules and weaker downflows in intergranular lanes. This trend is, however, inconsistent with the results of numerical simulations in which downflows are stronger than upflows through the joint action of gravitational acceleration/deceleration and pressure gradients. One cause of this discrepancy is the image degradation caused by optical distortion and light diffraction and scattering that takes place in an imaging instrument. We apply a deconvolution technique to Hinode/SP data in an attempt to recover the original solar scene. Our results show a significant enhancement in both, the convective upflows and downflows, but particularly for the latter. After deconvolution, the up- and downflows reach maximum amplitudes of -3.0 km/s and +3.0 km/s at an average geometrical height of roughly 50 km, respectively. We found that the velocity distributions after deconvolution match those derived from numerical simulations. After deconvolution the net LOS velocity averaged over the whole FOV lies close to zero as expected in a rough sense from mass balance.

The solar photosphere is the visible surface of the Sun, where many bright granules, surrounded by narrow dark intergranular lanes, are observed everywhere. The granular pattern is a manifestation of convective motion at the photospheric level, but its velocity structure in the height direction is poorly understood observationally. Applying bisector analysis to a photospheric spectral line recorded by the Hinode Solar Optical Telescope, we derived the velocity structure of the convective motion in granular regions and intergranular lanes separately. The amplitude of motion of the convective material decreases from 0.65 to 0.40 km/s as the material rises in granules, whereas the amplitude of motion increases from 0.30 to 0.50 km/s as it descends in intergranular lanes. These values are significantly larger than those obtained in previous studies using bisector analysis. The acceleration of descending materials with depth is not predicted from the convectively stable condition in a stratified atmosphere. Such convective instability can be developed more efficiently by radiative cooling and/or a gas pressure gradient, which can control the dynamical behavior of convective material in intergranular lanes. Our analysis demonstrated that bisector analysis is a useful method for investigating the long-term dynamic behavior of convective material when a large number of pixels is available. In addition, one example is the temporal evolution of granular fragmentation, in which downflowing material develops gradually from a higher layer downward.

The solar dynamo problem is the question of how the cyclic variation in the solar magnetic field is maintained. One of the important processes is the transport of magnetic flux by surface convection. To reveal this process, the dependence of the squared displacement of magnetic flux concentrations upon the elapsed time is investigated in this paper via a feature-recognition technique and a continual five-day magnetogram. This represents the longest time scale over which a satellite observation has ever been performed for this problem. The dependence is found to follow a power-law and differ significantly from that of diffusion transport. Furthermore there is a change in the behavior at a spatial scale of 10^{3.8} km. A super-diffusion behavior with an index of 1.4 is found on smaller scales, while changing to a sub-diffusion behavior with an index of 0.6 on larger ones. I interpret this difference in the transport regime as coming from the network-flow pattern.

In this second paper in the series, we investigate occurrence frequencies of apparent unipolar processes, cancellation, and emergence of patch structures in quiet regions. Apparent unipolar events are considerably more frequent than cancellation and emergence as per our definition, which is consistent with Lamb et al. (2013). Furthermore, we investigate the frequency distributions of changes in flux during apparent unipolar processes are and found that they concentrate around the detection limit of the analysis. Combining these findings with the results of our previous paper, Iida et al. (2012), that merging and splitting are more dominant than emergence and cancellation, these results support the understanding that apparent unipolar processes are actually interactions with and among patches below the detection limit and that there still are numerous flux interactions between the flux range in this analysis and below the detection limit. We also investigate occurrence frequency distributions of flux decrease during cancellation. We found a relatively strong dependence, 2.48$\pm$0:26 as a power-law index. This strong dependence on flux is consistent with the model, which is suggested in the previous paper.

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.

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

We investigate surface processes of magnetic patches, namely merging, splitting, emergence, and cancellation, by using an auto-detection technique. We find that merging and splitting are locally predominant in the surface level, while the frequencies of the other two are less by one or two orders of magnitude. The frequency dependences on flux con- tent of surface processes are further investigated. Based on these observations, we discuss a possible whole picture of the maintenance. Our conclusion is that the photospheric magnetic field structure, especially its power-law nature, is maintained by the processes locally in the surface not by the interactions between different altitudes. We suggest a scenario of the flux maintenance as follows: The splitting and merging play a crucial role for the generation of the power-law distribution, not the emergence nor cancellation do. This power-law distribution results in another power-law one of the cancellation with an idea of the random convective transport. The cancellation and emergence have a common value for the power-law indices in their frequency distributions, which may suggest a "recycle of fluxes by submergence and re-emergence".

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.

Frequencies of magnetic patch processes on supergranule boundary, namely flux emergence, splitting, merging, and cancellation, are investigated through an automatic detection. We use a set of line of sight magnetograms taken by the Solar Optical Telescope (SOT) on board Hinode satellite. We found 1636 positive patches and 1637 negative patches in the data set, whose time duration is 3.5 hours and field of view is 112" \times 112". Total numbers of magnetic processes are followed: 493 positive and 482 negative splittings, 536 positive and 535 negative mergings, 86 cancellations, and 3 emergences. Total numbers of emergence and cancellation are significantly smaller than those of splitting and merging. Further, frequency dependences of merging and splitting processes on flux content are investigated. Merging has a weak dependence on flux content only with a power- law index of 0.28. Timescale for splitting is found to be independent of parent flux content before splitting, which corresponds to \sim 33 minutes. It is also found that patches split into any flux contents with a same probability. This splitting has a power-law distribution of flux content with an index of -2 as a time independent solution. These results support that the frequency distribution of flux content in the analyzed flux range is rapidly maintained by merging and splitting, namely surface processes. We suggest a model for frequency distributions of cancellation and emergence based on this idea.

We investigate the value of horizontal turbulent diffusivity {\eta} by numerical calculation of thermal convection. In this study, we introduce a new method whereby the turbulent diffusivity is estimated by monitoring the time devel- opment of the passive scalar, which is initially distributed in a given Gaussian function with a spatial scale d0. Our conclusions are as follows: (1) Assuming the relation {\eta} = Lcvrms/3 where vrms is the RMS velocity, the characteristic length Lc is restricted by the shortest one among the pressure (density) scale height and the region depth. (2) The value of turbulent diffusivity becomes greater with the larger initial distribution scale d0. (3) The approximation of turbulent diffusion holds better when the ratio of the initial distribution scale d0 to the characteristic length Lc is larger.

We investigate the frequency of magnetic activities, namely flux emergence, splitting, merging, and cancellation, through an automatic detection in order to understand the generation of the power-law distribution of magnetic flux reported by Parnell et al. (2009). The quiet Sun magnetograms observed in the Na I line by Hinode/ SOT is used in this study. Investigated patches range from 10^{17}Mx to 10^{19}Mx. Emergence and cancellation are much less frequent than merging and splitting. The time scale for splitting is found to be ~33 minutes and is independent of the flux contained in the splitting patch. Moreover magnetic patches split into any flux contents with an equal probability. It is shown that such a fragmentation process leads to a distribution with a power-law index -2. Merging has a very weak dependence on flux content only with a power-law index -0.33. These results suggest that 1) magnetic patches are fragmented by splitting, merging, and tiny cancellation; 2) flux is removed from the photosphere through tiny cancellations after these fragmentations.

Context. The anelastic approximation is often adopted in numerical calculation with low Mach number, such as stellar internal convection. This approximation requires frequent global communication, because of an elliptic partial differential equation. Frequent global communication is negative factor for the parallel computing with a large number of CPUs. Aims. The main purpose of this paper is to test the validity of a method that artificially reduces the speed of sound for the compressible fluid equations in the context of stellar internal convection. The reduction of speed of sound allows for larger time steps in spite of low Mach number, while the numerical scheme remains fully explicit and the mathematical system is hyperbolic and thus does not require frequent global communication. Methods. Two and three dimensional compressible hydrodynamic equations are solved numerically. Some statistical quantities of solutions computed with different effective Mach numbers (due to reduction of speed of sound) are compared to test the validity of our approach. Results. Numerical simulations with artificially reduced speed of sound are a valid approach as long as the effective Mach number (based on the reduced speed of sound) remains less than 0.7.

A cancellation is thought to be a basic process of the photospheric magnetic field and plays an important role in magnetic flux budget and in various solar activities. There are two major theoretical scenarios for this phenomena, i.e. the "U-loop emergence" and the Omega-loop submergence models. It is important to clarify which is the dominant process during the cancellation for the estimation of the solar magnetic flux transport through the surface. We study the vector magnetic field and velocity structures around a quiet Sun cancellation by using the Solar Optical Telescope on board Hinode satellite. Transverse magnetic field connecting the canceling magnetic features and strong long-lasting Doppler red-shift signal are found. The transverse field is observed in the first spectropolarimetric observation after the occurrence of the cancellation while the red-shift clearly delayed to the cancellation by 20 minutes. These results indicate that the observed cancellation is an Omega-loop submergence.