In this Letter, we present a seismological detection of a rising motion of magnetic flux in the shallow convection zone of the Sun, and show estimates of the emerging speed and its decelerating nature. In order to evaluate the speed of subsurface flux that creates an active region, we apply six Fourier filters to the Doppler data of NOAA AR 10488, observed with the Solar and Heliospheric Observatory/Michelson Doppler Imager, to detect the reduction of acoustic power at six different depths from -15 to -2 Mm. All the filtered acoustic powers show reductions, up to 2 hr before the magnetic flux first appears at the visible surface. The start times of these reductions show a rising trend with a gradual deceleration. The obtained velocity is first several km s(-1) in a depth range of 15-10 Mm, then similar to 1.5 km s(-1) at 10-5 Mm, and finally similar to 0.5 km s(-1) at 5-2 Mm. If we assume that the power reduction is actually caused by the magnetic field, the velocity of the order of 1 km s(-1) is well in accordance with previous observations and numerical studies. Moreover, the gradual deceleration strongly supports the theoretical model that the emerging flux slows down in the uppermost convection zone before it expands into the atmosphere to build an active region.
K. Kusano, Y. Bamba, T. T. Yamamoto, Y. Iida, S. Toriumi, A. Asai
ASTROPHYSICAL JOURNAL 760(1) 31 2012年11月 査読有り
Solar flares and coronal mass ejections, 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 2006 December 13 and the M-class flare on 2011 February 13, 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 timescale of changes in the small magnetic structures.
It is widely accepted that solar active regions including sunspots are formed by the emerging magnetic flux from the deep convection zone. In previous numerical simulations, we found that the horizontal divergent flow (HDF) occurs before the flux emergence at the photospheric height. This paper reports the HDF detection prior to the flux emergence of NOAA AR 11081, which is located away from the disk center. We use SDO/HMI data to study the temporal changes of the Doppler and magnetic patterns from those of the reference quiet Sun. As a result, the HDF appearance is found to come before the flux emergence by about 100 minutes. Also, the horizontal speed of the HDF during this time gap is estimated to be 0.6-1.5 km s(-1), up to 2.3 km s(-1). The HDF is caused by the plasma escaping horizontally from the rising magnetic flux. And the interval between the HDF and the flux emergence may reflect the latency during which the magnetic flux beneath the solar surface is waiting for the instability onset to the further emergence. Moreover, SMART H alpha images show that the chromospheric plages appear about 14 minutes later, located cospatial with the photospheric pores. This indicates that the plages are caused by plasma flowing down along the magnetic fields that connect the pores at their footpoints. One important result of observing the HDF may be the possibility of predicting the sunspot appearances that occur in several hours.
We have performed a three-dimensional magnetohydrodynamic simulation to study the emergence of a twisted magnetic flux tube from - 20 000 km of the solar convection zone to the corona through the photosphere and the chromosphere. The middle part of the initial tube is endowed with a density deficit to instigate a buoyant emergence. As the tube approaches the surface, it extends horizontally and makes a flat magnetic structure due to the photosphere ahead of the tube. Further emergence to the corona breaks out via the interchange-mode instability of the photospheric fields, and eventually several magnetic domes build up above the surface. What is new in this three-dimensional experiment is multiple separation events of the vertical magnetic elements are observed in the photospheric magnetogram, and they reflect the interchange instability. Separated elements are found to gather at the edges of the active region. These gathered elements then show shearing motions. These characteristics are highly reminiscent of active region observations. On the basis of the simulation results above, we propose a theoretical picture of the flux emergence and the formation of an active region that explains the observational features, such as multiple separations of faculae and the shearing motion.
We present the new results of the two-dimensional numerical experiments on the cross-sectional evolution of a twisted magnetic flux tube rising from the deeper solar convection zone (-20,000 km) to the corona through the surface. The initial depth is 10 times deeper than most of the previous calculations focusing on the flux emergence from the uppermost convection zone. We find that the evolution is illustrated by the following two-step process. The initial tube rises due to its buoyancy, subject to aerodynamic drag due to the external flow. Because of the azimuthal component of the magnetic field, the tube maintains its coherency and does not deform to become a vortex roll pair. When the flux tube approaches the photosphere and expands sufficiently, the plasma on the rising tube accumulates to suppress the tube's emergence. Therefore, the flux decelerates and extends horizontally beneath the surface. This new finding owes to our large-scale simulation, which simultaneously calculates the dynamics within the interior as well as above the surface. As the magnetic pressure gradient increases around the surface, magnetic buoyancy instability is triggered locally and, as a result, the flux rises further into the solar corona. We also find that the deceleration occurs at a higher altitude than assumed in our previous experiment using magnetic flux sheets. By conducting parametric studies, we investigate the conditions for the two-step emergence of the rising flux tube: field strength greater than or similar to 1.5 x 10(4) G and the twist greater than or similar to 5.0 x 10(-4) km(-1) at -20,000 km depth.
PUBLICATIONS OF THE ASTRONOMICAL SOCIETY OF JAPAN 63(2) 407-415 2011年4月 査読有り筆頭著者責任著者
We present a series of numerical experiments that model the evolution of magnetic flux tubes with a different amount of initial twist. As a result of calculations, tightly twisted tubes reveal a rapid two-step emergence to the atmosphere with a slight slowdown at the surface, while weakly twisted tubes show a slow two-step emergence waiting longer the secondary instability to be triggered. This picture of the two-step emergence is highly consistent with recent observations. These tubes show multiple magnetic domes above the surface, indicating that the secondary emergence is caused by an interchange mode of magnetic buoyancy instability. In the case of the weakest twist, the tube exhibits an elongated photospheric structure, and never rises into the corona. The formation of the photospheric structure is due to an inward magnetic tension force of the azimuthal field component of the rising flux tube (i.e., tube's twist). When the twist is weak, the azimuthal field cannot hold the tube's coherency, and the tube extends laterally at the subadiabatic surface. In addition, we newly found that the total magnetic energy measured above the surface depends on the initial twist. Strong twist tubes follow the initial relation between the twist and the magnetic energy, while weak twist tubes deviate from this relation, because these tubes store their magnetic energy in the photospheric structure.
We perform two-dimensional magnetodydrodynamic simulations of the flux emergence from the solar convection zone to the corona. The flux sheet is initially located moderately deep in the adiabatically stratified convection zone (-20,000 km) and is perturbed to trigger the Parker instability. The flux rises through the solar interior due to the magnetic buoyancy, but suffers a gradual deceleration and a flattening in the middle of the way to the surface since the plasma piled on the emerging loop cannot pass through the convectively stable photosphere. As the magnetic pressure gradient enhances, the flux becomes locally unstable to the Parker instability so that the further evolution to the corona occurs. The second-step nonlinear emergence is well described by the expansion law by Shibata et al. To investigate the condition for this "two-step emergence" model, we vary the initial field strength and the total flux. When the initial field is too strong, the flux exhibits the emergence to the corona without a deceleration at the surface and reveals an unrealistically strong flux density at each footpoint of the coronal loop, while the flux either fragments within the convection zone or cannot pass through the surface when the initial field is too weak. The condition for the "two-step emergence" is found to be 10(21)-10(22) Mx with 10(4) G at z = -20,000 km. We present some discussions in connection with recent observations and the results of the thin-flux-tube model.
T. Shimizu, S. Imada, T. Kawate, K. Ichimoto, Y. Suematsu, H. Hara, Y. Katsukawa, M. Kubo, S. Toriumi, T. Watanabe, T. Yokoyama, C.M. Korendyke, H.P. Warren, T. Tarbell, B. De Pontieu, L. Teriaca, U.H. Schühle, S. Solanki, L.K. Harra, S. Matthews, A. Fludra, F. Auchère, V. Andretta, G. Naletto, A. Zhukov
S. Toriumi, T. Shimizu, S. Imada, T. Kawate, C. Quitero Noda, K. Ichimoto, H. Hara, T. Watanabe, Y. Suematsu, Y. Katsukawa, Solar-C WG, C. Korendyke, H. Warren, T. Tarbell, S. Solanki, L. Teriaca, L. Harra, A. Fludra, F. Auchere, A. Vincenzo, A. Zhukov