Showing 1–26 of 26 results for author: Shiota, D

The Solar and Geomagnetic Storms in May 2024: A Flash Data Report

Authors: Hisashi Hayakawa, Yusuke Ebihara, Alexander Mishev, Sergey Koldobskiy, Kanya Kusano, Sabrina Bechet, Seiji Yashiro, Kazumasa Iwai, Atsuki Shinbori, Kalevi Mursula, Fusa Miyake, Daikou Shiota, Marcos V. D. Silveira, Robert Stuart, Denny M. Oliveira, Sachiko Akiyama, Kouji Ohnishi, Yoshizumi Miyoshi

Abstract: In May 2024, the scientific community observed intense solar eruptions that resulted in an extreme geomagnetic storm and auroral extension, highlighting the need to document and quantify these events. This study mainly focuses on their quantification. The source active region (AR 13664) evolved from 113 to 2761 millionths of the solar hemisphere between 4 May and 14 May 2024. AR 13664's magnetic f… ▽ More In May 2024, the scientific community observed intense solar eruptions that resulted in an extreme geomagnetic storm and auroral extension, highlighting the need to document and quantify these events. This study mainly focuses on their quantification. The source active region (AR 13664) evolved from 113 to 2761 millionths of the solar hemisphere between 4 May and 14 May 2024. AR 13664's magnetic free energy surpassed 1033 erg on 7 May 2024, triggering 12 X-class flares. Multiple interplanetary coronal mass ejections (ICMEs) came out from this AR, accelerating solar energetic particles toward Earth. At least four ICMEs seemingly piled up to disturb the interplanetary space, according to the satellite data and interplanetary scintillation data. The shock arrival at 17:05 UT on 10 May 2024 significantly compressed the magnetosphere down to ~ 5.04 R_E, and triggered a deep Forbush decrease. GOES satellite data and ground-based neutron monitors confirmed a ground-level enhancement from 2 UT to 10 UT on 11 May 2024. The ICMEs induced extreme geomagnetic storms, peaking at a Dst index of -412 nT at 2 UT on 11 May 2024, marking the sixth-largest storm since 1957. The AE and AL indices showed extreme auroral extensions that located the AE/AL stations into the polar cap. We gathered auroral records at that time and reconstructed the equatorward boundary of the visual auroral oval to 29.8° invariant latitude. We compared naked-eye and camera auroral visibility, providing critical caveats on their difference. We also confirmed global enhancements of storm-enhanced density of the ionosphere. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 39 pages, 20 figures, and 2 tables. Figures 15 and 17 will appear on the record version. Under Review

Magnetohydrodynamic simulation of coronal mass ejections using interplanetary scintillation data observed from radio sites ISEE and LOFAR

Authors: Kazumasa Iwai, Richard A. Fallows, Mario M. Bisi, Daikou Shiota, Bernard V. Jackson, Munetoshi Tokumaru, Ken'ichi Fujiki

Abstract: Interplanetary scintillation (IPS) is a useful tool for detecting coronal mass ejections (CMEs) throughout interplanetary space. Global magnetohydrodynamic (MHD) simulations of the heliosphere, which are usually used to predict the arrival and geo-effectiveness of CMEs, can be improved using IPS data. In this study, we demonstrate an MHD simulation that includes IPS data from multiple stations to… ▽ More Interplanetary scintillation (IPS) is a useful tool for detecting coronal mass ejections (CMEs) throughout interplanetary space. Global magnetohydrodynamic (MHD) simulations of the heliosphere, which are usually used to predict the arrival and geo-effectiveness of CMEs, can be improved using IPS data. In this study, we demonstrate an MHD simulation that includes IPS data from multiple stations to improve CME modelling. The CMEs, which occurred on 09-10 September 2017, were observed over the period 10-12 September 2017 using the Low-Frequency Array (LOFAR) and IPS array of the Institute for Space-Earth Environmental Research (ISEE), Nagoya University, as they tracked through the inner heliosphere. We simulated CME propagation using a global MHD simulation, SUSANOO-CME, in which CMEs were modeled as spheromaks, and the IPS data were synthesised from the simulation results. The MHD simulation suggests that the CMEs merged in interplanetary space, forming complicated IPS g-level distributions in the sky map. We found that the MHD simulation that best fits both LOFAR and ISEE data provided a better reconstruction of the CMEs and a better forecast of their arrival at Earth than from measurements when these simulations were fit from the ISEE site alone. More IPS data observed from multiple stations at different local times in this study can help reconstruct the global structure of the CME, thus improving and evaluating the CME modelling. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 17 pages, 9 figures, and 2 tables, accepted for publication in Advances in Space Research

Testing the Alfvén-wave model of the solar wind with interplanetary scintillation

Authors: Munehito Shoda, Kazumasa Iwai, Daikou Shiota

Abstract: Understanding the mechanism(s) of the solar wind acceleration is important in astrophysics and geophysics. A promising model of the solar wind acceleration is known as the wave/turbulence-driven (WTD) model, in which Alfvén waves feed energy to the solar wind. In this study, we tested the WTD model with global measurement of wind speed from interplanetary scintillation (IPS) observations. For Carr… ▽ More Understanding the mechanism(s) of the solar wind acceleration is important in astrophysics and geophysics. A promising model of the solar wind acceleration is known as the wave/turbulence-driven (WTD) model, in which Alfvén waves feed energy to the solar wind. In this study, we tested the WTD model with global measurement of wind speed from interplanetary scintillation (IPS) observations. For Carrington rotations in minimal and maximal activity phases, we selected field lines calculated by the potential-field source-surface method in high- and mid-latitudes and compared the simulated and observed wind velocities. The simulation was performed in a self-consistent manner by solving the magnetohydrodynamic equations from the photosphere to the solar wind. In high-latitude regions, the simulated solar wind velocity agrees better with the IPS observation than with the classical Wang--Sheeley empirical estimation, both in maximal and minimal activity phases. In mid-latitude regions, the agreement worsens, possibly because of the inaccuracy of the WTD model and/or the magnetic-field extrapolation. Our results indicate that the high-latitude solar wind is likely to be driven by waves and turbulence, and that the physics-based prediction of the solar wind velocity is highly feasible with an improved magnetic-field extrapolation. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted for publication in the Astrophysical Journal

Validation of Coronal Mass Ejection Arrival-Time Forecasts by Magnetohydrodynamic Simulations based on Interplanetary Scintillation Observations

Authors: Kazumasa Iwai, Daikou Shiota, Munetoshi Tokumaru, Kenichi Fujiki, Mitsue Den, Yûki Kubo

Abstract: Coronal mass ejections (CMEs) cause various disturbances of the space environment; therefore, forecasting their arrival time is very important. However, forecasting accuracy is hindered by limited CME observations in interplanetary space. This study investigates the accuracy of CME arrival times at the Earth forecasted by three-dimensional (3D) magnetohydrodynamic (MHD) simulations based on interp… ▽ More Coronal mass ejections (CMEs) cause various disturbances of the space environment; therefore, forecasting their arrival time is very important. However, forecasting accuracy is hindered by limited CME observations in interplanetary space. This study investigates the accuracy of CME arrival times at the Earth forecasted by three-dimensional (3D) magnetohydrodynamic (MHD) simulations based on interplanetary scintillation (IPS) observations. In this system, CMEs are approximated as spheromaks with various initial speeds. Ten MHD simulations with different CME initial speed are tested, and the density distributions derived from each simulation run are compared with IPS data observed by the Institute for Space-Earth Environmental Research (ISEE), Nagoya University. The CME arrival time of the simulation run that most closely agrees with the IPS data is selected as the forecasted time. We then validate the accuracy of this forecast using 12 halo CME events. The average absolute arrival-time error of the IPS-based MHD forecast is approximately 5.0 h, which is one of the most accurate predictions that ever been validated, whereas that of MHD simulations without IPS data, in which the initial CME speed is derived from white-light coronagraph images, is approximately 6.7 h. This suggests that the assimilation of IPS data into MHD simulations can improve the accuracy of CME arrival-time forecasts. The average predicted arrival times are earlier than the actual arrival times. These early predictions may be due to overestimation of the magnetic field included in the spheromak and/or underestimation of the drag force from the background solar wind, the latter of which could be related to underestimation of CME size or background solar wind density. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: 17 pages, 4 figures, and 1 table, accepted for publication in Earth, Planets and Space

Suzaku observation of Jupiter's X-rays around solar maximum

Authors: Masaki Numazawa, Yuichiro Ezoe, Takaya Ohashi, Kumi Ishikawa, Yoshizumi Miyoshi, Tomoki Kimura, Yasunobu Uchiyama, Daikou Shiota, Graziella Branduardi-Raymont

Abstract: 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 th… ▽ More 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 ($ Γ= 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. △ Less

Submitted 14 June, 2019; originally announced June 2019.

Comments: 17 pages, 4 figures, 1 table

Development of a coronal mass ejection arrival time forecasting system using interplanetary scintillation observations

Authors: Kazumasa Iwai, Daikou Shiota, Munetoshi Tokumaru, Kenichi Fujiki, Mitsue Den, Yûki Kubo

Abstract: 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 ini… ▽ More 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. △ Less

Submitted 27 March, 2019; originally announced March 2019.

Comments: 39 pages, 6 figures, accepted for publication in Earth, Planets and Space

Magnetohydrodynamic Modeling of a Solar Eruption Associated with X9.3 Flare Observed in Active Region 12673

Authors: Satoshi Inoue, Daikou Shiota, Yumi Bamba, Sung-Hong Park

Abstract: 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 photosph… ▽ More 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. △ Less

Submitted 7 September, 2018; originally announced September 2018.

Comments: Accepted for publication in ApJ

Measurement of Vector Magnetic Field in a Flare kernel with a Spectropolarimetric Observation in He I 10830 A

Authors: Tetsu Anan, Takurou Yoneya, Kiyoshi Ichimoto, Satoru Ueno, Daikou Shiota, Satoshi Nozawa, Shinsuke Takasao, Tomoko Kawate

Abstract: 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 Pasche… ▽ More 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. △ Less

Submitted 21 August, 2018; originally announced August 2018.

Comments: 15 pages, 8 figures, 1 table, accepted for publication in PASJ (Hida topical issue)

Improvement of solar cycle prediction: Plateau of solar axial dipole moment

Authors: H. Iijima, H. Hotta, S. Imada, K. Kusano, D. Shiota

Abstract: 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 (MD… ▽ More 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. △ Less

Submitted 17 October, 2017; originally announced October 2017.

Comments: 5 pages, 3 figures, accepted for publication in A&A Letter

Journal ref: A&A 607, L2 (2017)

Magnetohydrodynamic Simulations for Studying Solar Flare Trigger Mechanism

Authors: Johan Muhamad, Kanya Kusano, Satoshi Inoue, Daikou Shiota

Abstract: 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 k… ▽ More 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. △ Less

Submitted 21 June, 2017; originally announced June 2017.

Comments: 28 pages, 10 figures

Triggering Process of the X1.0 Three-ribbon Flare in the Great Active Region NOAA 12192

Authors: Yumi Bamba, Satoshi Inoue, Kanya Kusano, Daikou Shiota

Abstract: 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… ▽ More 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. △ Less

Submitted 4 April, 2017; originally announced April 2017.

Journal ref: The Astrophysical Journal, 838:134 (16pp), 2017 April 1

Coronal Magnetic Fields Derived from Simultaneous Microwave and EUV Observations and Comparison with the Potential Field Model

Authors: Shun Miyawaki, Kazumasa Iwai, Kiyoto Shibasaki, Daikou Shiota, Satoshi Nozawa

Abstract: 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 Nob… ▽ More 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. △ Less

Submitted 14 December, 2015; originally announced December 2015.

Comments: accepted for publication in ApJ

Photospheric Flow Field Related to the Evolution of the Sun's Polar Magnetic Patches Observed by Hinode SOT

Authors: Anjali John Kaithakkal, Y. Suematsu, M. Kubo, Y. Iida, D. Shiota, S. Tsuneta

Abstract: 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 asso… ▽ More 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. △ Less

Submitted 3 January, 2015; v1 submitted 27 December, 2014; originally announced December 2014.

Comments: 13 pages, 21 figures, Accepted for Publication in ApJ

Nonlinear Force-Free Extrapolation of the Coronal Magnetic Field Based on the MHD Relaxation Method

Authors: S. Inoue, T. Magara, V. S. Pandey, D. Shiota., K. Kusano, G. S. Choe, K. S. Kim

Abstract: 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 implemen… ▽ More 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. △ Less

Submitted 14 November, 2013; originally announced November 2013.

Comments: 36 pages, 12 figures, accepted for publication in The Astrophysical Journal

The Association of Polar Faculae with Polar Magnetic Patches Examined with Hinode Observations

Authors: Anjali John Kaithakkal, Y. Suematsu, M. Kubo, D. Shiota, S. Tsuneta

Abstract: 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 considerabl… ▽ More 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. △ Less

Submitted 5 November, 2013; originally announced November 2013.

Comments: 14 pages, 10 figures

Journal ref: The Astrophysical Journal, 776:122 (6pp), 2013 October 20

A Plasma β Transition Within a Propagating Flux Rope

Authors: Neel P. Savani, A. Vourlidas, D. Shiota, M. G. Linton, K. Kusano, N. Lugaz, A. P. Rouillard

Abstract: 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 t… ▽ More 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. △ Less

Submitted 27 November, 2013; v1 submitted 17 October, 2013; originally announced October 2013.

Comments: 11 pages of text + 6 figures. Accepted to ApJ on 16 Oct 2013

Magnetic Structure Producing X- and M-Class Solar Flares in Solar Active Region 11158

Authors: S. Inoue, K. Hayashi, D. Shiota, T. Magara, G. S. Choe

Abstract: 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 b… ▽ More 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. △ Less

Submitted 29 April, 2013; originally announced April 2013.

Comments: 30 pages, 10 figures, The Astrophysical Journal, accepted

Buildup and Release of Magnetic Twist during the X3.4 Solar Flare of December 13, 2006

Authors: S. Inoue, D. Shiota, T. T. Yamamoto, V. S. Pandey, T. Magara, G. S. Choe

Abstract: 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)… ▽ More 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. △ Less

Submitted 27 September, 2012; originally announced September 2012.

Comments: 28 pages, 8 figures, The Astrophysical Journal, accepted

A study of the Heliocentric dependence of Shock Standoff Distance and Geometry using 2.5D MHD Simulations of CME-driven shocks

Authors: Neel P. Savani, Daikou Shiota, Kanya Kusano, Angelos Vourlidas, Noé Lugaz

Abstract: 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 magne… ▽ More 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. △ Less

Submitted 10 September, 2012; originally announced September 2012.

Comments: 28 Pages; 9 figures embedded in the article.pdf. Accepted in Astrophysical Journal

Polar Field Reversal Observations with Hinode

Authors: D. Shiota, S. Tsuneta, M. Shimojo, N. Sako, D. Orozco Suárez, R. Ishikawa

Abstract: 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 magn… ▽ More 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. △ Less

Submitted 10 May, 2012; originally announced May 2012.

Comments: 21 pages, 7 figures. Accepted for publication in ApJ

First Simultaneous Observation of H-alpha Moreton Wave, EUV Wave, and Filament/Prominence Oscillations

Authors: Ayumi Asai, Takako T. Ishii, Hiroaki Isobe, Reizaburo Kitai, Satoru UeNo, Shin'ichi Nagata, Satoshi Morita, Keisuke Nishida, Daikou Shiota, Akihito Oi, Maki Akioka, Kazunari Shibata

Abstract: 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 K… ▽ More 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. △ Less

Submitted 26 December, 2011; originally announced December 2011.

Comments: 15 pages, 5 figures, accepted for publication in ApJL

MHD Modeling for Formation Process of Coronal Mass Ejections: Interaction between Ejecting Flux Rope and Ambient Field

Authors: Daikou Shiota, Kanya Kusano, Takahiro Miyoshi, Kazunari Shibata

Abstract: 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 res… ▽ More 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. △ Less

Submitted 2 June, 2010; originally announced June 2010.

Comments: 24 pages, 5 figures. Accepted for publication in ApJ

Is the Polar Region Different from the Quiet Region of the Sun?

Authors: Hiroaki Ito, Saku Tsuneta, Daikou Shiota, Munetoshi Tokumaru, Ken'ichi Fujiki

Abstract: 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… ▽ More 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. △ Less

Submitted 20 May, 2010; originally announced May 2010.

Comments: Accepted for publication in APJ

Signatures of Coronal Heating Mechanisms

Authors: P. Antolin, K. Shibata, T. Kudoh, D. Shiota, D. Brooks

Abstract: 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… ▽ More 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. △ Less

Submitted 10 March, 2009; originally announced March 2009.

Comments: To appear in "Magnetic Coupling between the Interior and the Atmosphere of the Sun", eds. S.S. Hasan and R.J. Rutten, Astrophysics and Space Science Proceedings, Springer-Verlag, Heidelberg, Berlin, 2009

Numerical examination of plasmoid-induced reconnection model for solar flares: the relation between plasmoid velocity and reconnection rate

Authors: Keisuke Nishida, Masaki Shimizu, Daikou Shiota, Hiroyuki Takasaki, Tetsuya Magara, Kazunari Shibata

Abstract: 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 val… ▽ More 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. △ Less

Submitted 4 September, 2008; originally announced September 2008.

Comments: 27 pages, 12 figures, Accepted for publication in ApJ

Journal ref: Astrophys.J.690:748-757,2009

Self-Consistent MHD Modeling of a Coronal Mass Ejection, Coronal Dimming, and a Giant Cusp-Shaped Arcade Formation

Authors: Daikou Shiota, Hiroaki Isobe, P. F. Chen, Tetsuya T. Yamamoto, Takuma Sakajiri, Kazunari Shibata

Abstract: 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… ▽ More 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. △ Less

Submitted 23 August, 2005; originally announced August 2005.

Comments: 39 pages, 18 figures. Accepted for publication in the Astrophysical Journal. The PDF file with high resplution figures can be downloaded from http://www.kwasan.kyoto-u.ac.jp/~shiota/study/ApJ62426.preprint.pdf"

Journal ref: Astrophys.J.634:663-678,2005