Source author record

Kazunari Shibata

Kazunari Shibata appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

70works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

70 published item(s)

preprint2026arXiv

Statistical Study of Appearance Timing of H$α$ Postflare Loops: Simple Scaling Law Based on Radiative Cooling

Recent Sun-as-a-star studies have shown that postflare loops can manifest as a secondary peak in the H$α$ light curve, suggesting that stellar postflare loops are detectable. To understand what determines the timing of such a secondary peak in the H$α$ light curve associated with postflare loops, we must quantitatively identify the key physical processes controlling the appearance of H$α$ postflare loops. Previous case studies have indicated that the appearance timing of H$α$ postflare loops is likely governed by radiative cooling. However, the statistical characteristics of the timing of H$α$ postflare loops appearance remain insufficiently investigated. In this study, we statistically investigated the appearance timing of H$α$ postflare loops to quantify their cooling processes. As a result, we found a negative correlation between the time difference between the soft X-ray peak and the appearance of the H$α$ postflare loops ($Δt$) and the soft X-ray peak flux ($F_\mathrm{X}$). This relationship is consistent with the theoretical scaling between radiative cooling timescale ($τ_{\mathrm{rad}}$) and $F_\mathrm{X}$, where $τ_{\mathrm{rad}} \propto~F_\mathrm{X}^{-1/2}$. This statistical result indicates that the appearance timing of H$α$ postflare loops relative to the soft X-ray peak is primarily controlled by radiative cooling. Furthermore, we examined the dependence of the scaling law on flare spatial scales ($L$). Consequently, we demonstrated that spatial scale of unresolved stellar flares can be estimated using the following scaling law: $L\propto F_\mathrm{X}^{1/3}Δt^{2/3}$. Our results are useful for interpreting secondary peaks in the H$α$ data of stellar flares and provide new method to estimate spatial scale of unresolved stellar flares.

preprint2023arXiv

Multiwavelength observation of an active M-dwarf star EV Lac and its stellar flare accompanied by a delayed prominence eruption

We conducted 4-night multiwavelength observations of an active M-dwarf star EV Lac on 2022 October 24$-$27 with simultaneous coverage of soft X-rays (NICER; 0.2$-$12 $\mathrm{keV}$, Swift XRT; 0.2$-$10 $\mathrm{keV}$), near-ultraviolet (Swift UVOT/UVW2; 1600$-$3500 Å), optical photometry (TESS; 6000$-$10000 Å), and optical spectroscopy (Nayuta/MALLS; 6350$-$6800 Å). During the campaign, we detected a flare starting at 12:28 UTC on October 25 with its white-light bolometric energy of $3.4 \times 10^{32}$ erg. At about 1 hour after this flare peak, our $\mathrm{Hα}$ spectrum showed a blue-shifted excess component at its corresponding velocity of $\sim 100 \: \mathrm{km \: s^{-1}}$. This may indicate that the prominence erupted with a 1-hour delay of the flare peak. Furthermore, the simultaneous 20-second cadence near-ultraviolet and white-light curves show gradual and rapid brightening behaviors during the rising phase at this flare. The ratio of flux in NUV to white light at the gradual brightening was $\sim 0.49$, which may suggest that the temperature of the blackbody is low ($< 9000 \: \mathrm{K}$) or the maximum energy flux of a nonthermal electron beam is less than $5\times10^{11} \: \mathrm{erg \: cm^{-2} \: s^{-1}}$. Our simultaneous observations of NUV and white-light flare raise the issue of a simple estimation of UV flux from optical continuum data by using a blackbody model.

preprint2022arXiv

Discovery of a Long-Duration Superflare on a Young Solar-Type Star EK Draconis with Nearly Similar Time Evolution for H$α$ and White-Light Emissions

Young solar-type stars are known to show frequent "superflares", which may severely influence the habitable worlds on young planets via intense radiations and coronal mass ejections. Here we report an optical spectroscopic and photometric observation of a long-duration superflare on the young solar-type star EK Draconis (50-120 Myr age) with the Seimei telescope and $Transiting$ $Exoplanet$ $Survey$ $Satellite$ ($TESS$). The flare energy 2.6$\times$10$^{34}$ erg and white-light flare duration 2.2 hr are much larger than those of the largest solar flares, and this is the largest superflare on a solar-type star ever detected by optical spectroscopy. The H$α$ emission profile shows no significant line asymmetry, meaning no signature of a filament eruption, unlike the only previous detection of a superflare on this star (Namekata et al. 2021, $Nat.Astron$). Also, it did not show significant line broadening, indicating that the non-thermal heating at the flare footpoints are not essential or that the footpoints are behind the limb. The time evolution and duration of the H$α$ flare are surprisingly almost the same as those of the white-light flare, which is different from general M-dwarf (super-)flares and solar flares. This unexpected time evolution may suggest that different radiation mechanisms than general solar flares are predominant, as follows: (1) radiation from (off-limb) flare loops, and (2) re-radiation via radiative backwarming, in both of which the cooling timescales of flare loops could determine the timescales of H$α$ and white light.

preprint2022arXiv

Extreme solar events

We trace the evolution of research on extreme solar and solar-terrestrial events from the 1859 Carrington event to the rapid development of the last twenty years. Our focus is on the largest observed/inferred/theoretical cases of sunspot groups, flares on the Sun and Sun-like stars, coronal mass ejections, solar proton events, and geomagnetic storms. The reviewed studies are based on modern observations, historical or long-term data including the auroral and cosmogenic radionuclide record, and Kepler observations of Sun-like stars. We compile a table of 100- and 1000-year events based on occurrence frequency distributions for the space weather phenomena listed above. Questions considered include the Sun-like nature of superflare stars and the existence of impactful but unpredictable solar "black swans" and extreme "dragon king" solar phenomena that can involve different physics from that operating in events which are merely large.

preprint2022arXiv

Sun-as-a-star Analysis of H$α$ Spectra of a Solar Flare Observed by SMART/SDDI: Time Evolution of Red Asymmetry and Line Broadening

Stellar flares sometimes show red/blue asymmetries of H$α$ line, which can indicate chromospheric dynamics and prominence activations. However, the origin of asymmetries is not completely understood. For a deeper understanding of stellar data, we performed a Sun-as-a-star analysis of H$α$ line profiles of an M4.2-class solar flare showing dominant emissions from flare ribbons by using the data of the Solar Dynamics Doppler Imager onboard the Solar Magnetic Activity Research Telescope at Hida Observatory. The Sun-as-a-star H$α$ spectra of the flare show red asymmetry of up to $\sim$95 km s$^{-1}$ and line broadening of up to $\sim$7.5 Å. The Sun-as-a-star H$α$ profiles are consistent with spectra from flare regions with weak intensity, but they take smaller redshift velocities and line widths by a factor of $\sim$2 than those with strong intensity. The redshift velocities, as well as line widths, peak out and decay more rapidly than the H$α$ equivalent widths, which is consistent with chromospheric condensation model and spatially-resolved flare spectra. This suggests that as a result of superposition, the nature of chromospheric condensation is observable even from stellar flare spectra. The time evolution of redshift velocities is found to be similar to that of luminosities of near-ultraviolet rays (1600 Å), while the time evolution of line broadening is similar to that of optical white lights. These H$α$ spectral behaviors in Sun-as-a-star view could be helpful to distinguish whether the origin of H$α$ red asymmetry of stellar flares is a flare ribbon or other phenomena.

preprint2021arXiv

A Three-Dimensional Velocity of an Erupting Prominence Prior to a Coronal Mass Ejection

We present a detailed three-dimensional (3D) view of a prominence eruption, coronal loop expansion, and coronal mass ejections (CMEs) associated with an M4.4 flare that occurred on 2011 March 8 in the active region NOAA 11165. Full-disk H$α$ images of the flare and filament ejection were successfully obtained by the Flare Monitoring Telescope (FMT) following its relocation to Ica University, Peru. Multiwavelength observation around the H$α$ line enabled us to derive the 3D velocity field of the H$α$ prominence eruption. Features in extreme ultraviolet were also obtained by the Atmospheric Imager Assembly onboard the {\it Solar Dynamic Observatory} and the Extreme Ultraviolet Imager on board the {\it Solar Terrestrial Relations Observatory - Ahead} satellite. We found that, following collision of the erupted filament with the coronal magnetic field, some coronal loops began to expand, leading to the growth of a clear CME. We also discuss the succeeding activities of CME driven by multiple interactions between the expanding loops and the surrounding coronal magnetic field.

preprint2020arXiv

Cost estimation for alternative aviation plans against potential radiation exposure associated with solar proton events for the airline industry

We present a systematic approach to effectively evaluate potential risk cost caused by exposure to solar proton events (SPEs) from solar flares for the airline industry. We also evaluate associated health risks from radiation, to provide relevant alternative ways to minimize economic loss and opportunity. The estimated radiation dose induced by each SPE for the passengers of each flight is calculated using ExoKyoto and PHITS. We determine a few scenarios for the estimated dose limit at 1 and 20mSv, corresponding to the effective dose limit for the general public and occupational exposure, respectively, as well as a higher dose induced an extreme superflare. We set a hypothetical airline shutdown scenario at 1mSv for a single flight per passenger, due to legal restrictions under the potential radiation dose. In such a scenario, we calculate the potential loss in direct and opportunity cost under the cancelation of the flight. At the same time, we considered that, even under such a scenario, if the airplane flies at a slightly lower altitude (from 12 to 9.5km: atmospheric depth from 234 to 365g/cm$^{2}$), the total loss becomes much smaller than flight cancelation, and the estimated total dose goes down from 1.2 to 0.45mSv, which is below the effective dose limit for the general public. In case of flying at an even lower altitude (7km: atmospheric depth 484g/cm$^{2}$), the estimated total dose becomes much smaller, 0.12 mSv. If we assume the increase of fuel cost is proportional to the increase in atmospheric depth, the increase in cost becomes 1.56 and 2.07 for the case of flying at 9.5 km and at 7 km, respectively. Lower altitude flights provide more safety for the potential risk of radiation doses induced by severe SPEs. At the same time, since there is total loss caused by flight cancelation, we propose that considering lower flight altitude is the best protection against solar flares.

preprint2020arXiv

Energy Transfer by Nonlinear Alfvén Waves in the Solar Chromosphere, and Its Effect on Spicule Dynamics, Coronal Heating, and Solar Wind Acceleration

Alfvén waves are responsible for the transfer of magnetic energy in the magnetized plasma. They are involved in heating solar atmosphere and driving solar wind through various nonlinear processes. Since the magnetic field configurations directly affect the nonlinearity of Alfvén waves, it is important to investigate how they relate to the solar atmosphere and wind structure through the nonlinear propagation of Alfvén waves. In this study, we carried out the one-dimensional magnetohydrodynamic simulations to realize the above relation. The results show that when the nonlinearity of Alfvén waves in the chromosphere exceeds a critical value, the dynamics of the solar chromosphere (e.g., spicule) and the mass loss rate of solar wind tend to be independent of the energy input from the photosphere. In a situation where the Alfvén waves are highly nonlinear, the strong shear torsional flow generated in the chromosphere ``fractures'' the magnetic flux tube. This corresponds to the formation of chromospheric intermediate shocks, which limit the transmission of the Poynting flux into the corona by Alfvén waves and also inhibits the propagation of chromospheric slow shock.

preprint2020arXiv

Nonlinear Alfvén Wave Model of Stellar Coronae and Winds from the Sun to M dwarfs

M dwarf's atmosphere and wind is expected to be highly magnetized. The nonlinear propagation of Alfvén wave could play a key role in both heating the stellar atmosphere and driving the stellar wind. Along this Alfvén wave scenario, we carried out the one-dimensional compressive magnetohydrodynamic (MHD) simulation about the nonlinear propagation of Alfvén wave from the M dwarf's photosphere, chromosphere to the corona and interplanetary space. Based on the simulation results, we develop the semi-empirical method describing the solar and M dwarf's coronal temperature, stellar wind velocity, and wind's mass loss rate. We find that M dwarfs' coronae tend to be cooler than solar corona, and that M dwarfs' stellar winds would be characterized with faster velocity and much smaller mass loss rate compared to those of the solar wind.

preprint2020arXiv

Numerical simulation of solar photospheric jet-like phenomena caused by magnetic reconnection

Jet phenomena with a bright loop in their footpoint, called anemone jets, have been observed in the solar corona and chromosphere. These jets are formed as a consequence of magnetic reconnection, and from the scale universality of magnetohydrodynamics (MHD), it can be expected that anemone jets exist even in the solar photosphere. However, it is not necessarily apparent that jets can be generated as a result of magnetic reconnection in the photosphere, where the magnetic energy is not dominant. Furthermore, MHD waves generated from the photospheric jets could contribute to chromospheric heating and spicule formation; however, this hypothesis has not yet been thoroughly investigated. In this study, we perform 3D MHD simulation including gravity with the solar photospheric parameter to investigate anemone jets in the solar photosphere. In the simulation, jet-like structures were induced by magnetic reconnection in the solar photosphere. We determined that these jet-like structures were caused by slow shocks formed by the reconnection and were propagated approximately in the direction of the background magnetic field. We also suggested that MHD waves from the jet-like structures could influence local atmospheric heating and spicule formation.

preprint2020arXiv

Optical and X-ray observations of stellar flares on an active M dwarf AD Leonis with Seimei Telescope, SCAT, NICER and OISTER

We report multi-wavelength monitoring observations of an M-dwarf flare star AD Leonis with Seimei Telescope (6150--7930 Å), SCAT (Spectroscopic Chuo-university Astronomical Telescope; 3700--7500 Å), NICER (Neutron Star Interior Composition Explorer; 0.2--12.0 keV), and collaborations of OISTER (Optical and Infrared Synergetic Telescopes for Education and Research) program. Twelve flares are detected in total which include ten H$α$, four X-ray, and four optical-continuum flares; one of them is a superflare with the total energy of $\sim$ 2.0$\times$10$^{33}$ erg. We found that (1) during the superflare, the H$α$ emission line full width at 1/8 maximum dramatically increases to 14 Å from 8 Å in the low-resolution spectra (R$\sim$ 2000) accompanied with the large white-light flares, (2) some weak H$α$/X-ray flares are not accompanied with white-light emissions, and (3) the non-flaring emissions show clear rotational modulations in X-ray and H$α$ intensity in the same phase. To understand these observational features, one-dimensional hydrodynamic flare simulations are performed by using the RADYN code. As a result of simulations, we found the simulated H$α$ line profiles with hard and high-energy non-thermal electron beams are consistent with that of the initial phase line profiles of the superflares, while those with more soft- and/or weak-energy beam are consistent with those in decay phases, indicating the changes in the energy fluxes injected to the lower atmosphere. Also, we found that the relation between optical continuum and H$α$ intensity is nonlinear, which can be one cause of the non-white-light flares. The flare energy budget exhibits diversity in the observations and models, and more observations of stellar flares are necessary for constraining the occurrence of various emission line phenomena in stellar flares.

preprint2020arXiv

SMART/SDDI Filament Disappearance Catalogue

This paper describes a new SMART/SDDI Filament Disappearance Catalogue, in which we listed almost all the filament disappearance events that the Solar Dynamics Doppler Imager (SDDI) has observed since its installation on the Solar Magnetic Activity Research Telescope (SMART) in May 2016. Our aim is to build a database that can help predict the occurrence and severity of coronal mass ejections (CMEs). The catalogue contains miscellaneous information associated with filament disappearance such as flare, CME, active region, three-dimensional trajectory of erupting filaments, detection in Interplanetary Scintillation (IPS), occurrence of interplanetary CME (ICME) and Dst index. We also provide statistical information on the catalogue data. The catalogue is available from the following website: https://www.kwasan.kyoto-u.ac.jp/observation/event/sddi-catalogue/.

preprint2020arXiv

Statistical Properties of Superflares on Solar-type Stars: Results Using All of the Kepler Primary Mission Data

We report the latest statistical analyses of superflares on solar-type (G-type main-sequence; effective temperature is 5100 - 6000 K) stars using all of the $Kepler$ primary mission data, and $Gaia$-DR2 (Data Release 2) catalog. We updated the flare detection method from our previous studies by using high-pass filter to remove rotational variations caused by starspots. We also examined the sample biases on the frequency of superflares, taking into account gyrochronology and flare detection completeness. The sample size of solar-type stars and Sun-like stars (effective temperature is 5600 - 6000 K and rotation period is over 20 days in solar-type stars) are $\sim$4 and $\sim$12 times, respectively, compared with Notsu et al. (2019, ApJ, 876, 58). As a result, we found 2341 superflares on 265 solar-type stars, and 26 superflares on 15 Sun-like stars: the former increased from 527 to 2341 and the latter from 3 to 26 events compared with our previous study. This enabled us to have a more well-established view on the statistical properties of superflares. The observed upper limit of the flare energy decreases as the rotation period increases in solar-type stars. The frequency of superflares decreases as the stellar rotation period increases. The maximum energy we found on Sun-like stars is $4 \times 10^{34}$ erg. Our analysis of Sun-like stars suggest that the Sun can cause superflares with energies of $\sim 7 \times 10^{33}$ erg ($\sim$X700-class flares) and $\sim 1 \times 10^{34}$ erg ($\sim$X1000-class flares) once every $\sim$3,000 years and $\sim$6,000 years, respectively.

preprint2020arXiv

Temporal Evolution of Spatially-Resolved Individual Star Spots on a Planet-Hosting Solar-type Star: Kepler 17

Star spot evolution is visible evidence of the emergence/decay of the magnetic field on stellar surface, and it is therefore important for the understanding of the underlying stellar dynamo and consequential stellar flares. In this paper, we report the temporal evolution of individual star spot area on the hot-Jupiter-hosting active solar-type star Kepler 17 whose transits occur every 1.5 days. The spot longitude and area evolution are estimated (1) from the stellar rotational modulations of Kepler data and (2) from the brightness enhancements during the exoplanet transits caused by existence of large star spots. As a result of the comparison, number of spots, spot locations, and the temporal evolution derived from the rotational modulations is largely different from those of in-transit spots. We confirm that although only two light curve minima appear per rotation, there are clearly many spots present on the star. We find that the observed differential intensity changes are sometimes consistent with the spot pattern detected by transits, but they sometimes do not match with each other. Although the temporal evolution derived from the rotational modulation differs from those of in-transit spots to a certain degree, the emergence/decay rates of in-transit spots are within an order of magnitude of those derived for sunspots as well as our previous research based only on rotational modulations. This supports a hypothesis that the emergence/decay of sunspots and extremely-large star spots on solar-type stars occur through the same underlying processes.

preprint2020arXiv

Time-resolved spectroscopy and photometry of an M dwarf flare star YZ Canis Minoris with OISTER and TESS: Blue asymmetry in H$α$ line during the non-white light flare

In this paper, we present the results from spectroscopic and photometric observations of the M-type flare star YZ CMi in the framework of the Optical and Infrared Synergetic Telescopes for Education and Research (OISTER) collaborations during the Transiting Exoplanet Survey Satellite (TESS) observation period. We detected 145 white-light flares from the TESS light curve and 4 H$α$ flares from the OISTER observations performed between 2019-01-16 and 2019-01-18. Among them, 3 H$α$ flares were associated with white-light flares. However, one of them did not show clear brightening in continuum; during this flare, the H$α$ line exhibited blue-asymmetry which has lasted for $\sim 60$ min. The line of sight velocity of the blue-shifted component is $-80$ - $-100$ km s$^{-1}$. This suggests that there can be upward flows of chromospheric cool plasma even without detectable red/NIR continuum brightening. By assuming that the blue-asymmetry in H$α$ line was caused by a prominence eruption on YZ CMi, we estimated the mass and kinetic energy of the upward-moving material to be $10^{16}$ - $10^{18}$ g and $10^{29.5}$ - $10^{31.5}$ erg, respectively. The estimated mass is comparable to expectations from the empirical relation between the flare X-ray energy and mass of upward-moving material for stellar flares and solar CMEs. In contrast, the estimated kinetic energy for the non-white-light flare on YZ CMi is roughly $2$ orders of magnitude smaller than that expected from the relation between flare X-ray energy and kinetic energy for solar CMEs. This could be understood by the difference in the velocity between CMEs and prominence eruptions.

preprint2019arXiv

Estimating the temperature and density of a spicule from 100 GHz data obtained with ALMA

We succeeded in observing two large spicules simultaneously with the Atacama Large Millimeter/submillimeter Array (ALMA), the Interface Region Imaging Spectrograph (IRIS), and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory. One is a spicule seen in the IRIS Mg II slit-jaw images and AIA 304Å images (MgII/304A spicule). The other one is a spicule seen in the 100GHz images obtained with ALMA (100GHz spicule). Although the 100GHz spicule overlapped with the MgII/304A spicule in the early phase, it did not show any corresponding structures in the IRIS Mg II and AIA 304A images after the early phase. It suggests that the spicules are individual events and do not have a physical relationship. To obtain the physical parameters of the 100GHz spicule, we estimate the optical depths as a function of temperature and density using two different methods. One is using the observed brightness temperature by assuming a filling factor, and the other is using an emission model for the optical depth. As a result of comparing them, the kinetic temperature of the plasma and the number density of ionized hydrogens in the 100GHz spicule are ~6800 K and 2.2 x 10^10 cm^-3. The estimated values can explain the absorbing structure in the 193A image, which appear as a counterpart of the 100GHz spicule. These results suggest that the 100GHz spicule presented in this paper is classified to a macrospicule without a hot sheath in former terminology.

preprint2019arXiv

IRIS Mg II Observations and Non-LTE Modeling of Off-limb Spicules in a Solar Polar Coronal Hole

We investigated the off-limb spicules observed in the Mg II h and k lines by IRIS in a solar polar coronal hole. We analyzed the large dataset of obtained spectra to extract quantitative information about the line intensities, shifts, and widths. The observed Mg II line profiles are broad and double-peaked at lower altitudes, broad but flat-topped at middle altitudes, and narrow and single-peaked with the largest Doppler shifts at higher altitudes. We use 1D non-LTE vertical slab models (i.e. models which consider departures from Local Thermodynamic Equilibrium) in single-slab and multi-slab configurations to interpret the observations and to investigate how a superposition of spicules along the line of sight (LOS) affects the synthetic Mg II line profiles. The used multi-slab models are either static, i.e. without any LOS velocities, or assume randomly assigned LOS velocities of individual slabs, representing the spicule dynamics. We conducted such single-slab and multi-slab modeling for a broad set of model input parameters and showed the dependence of the Mg II line profiles on these parameters. We demonstrated that the observed line widths of the h and k line profiles are strongly affected by the presence of multiple spicules along the LOS. We later showed that the profiles obtained at higher altitudes can be reproduced by single-slab models representing individual spicules. We found that the multi-slab model with a random distribution of the LOS velocities ranging from -25 to 25 km s$^{-1}$ can well reproduce the width and the shape of Mg II profiles observed at middle altitudes.

preprint2016arXiv

A Deep-Learning Approach for Operation of an Automated Realtime Flare Forecast

Automated forecasts serve important role in space weather science, by providing statistical insights to flare-trigger mechanisms, and by enabling tailor-made forecasts and high-frequency forecasts. Only by realtime forecast we can experimentally measure the performance of flare-forecasting methods while confidently avoiding overlearning. We have been operating unmanned flare forecast service since August, 2015 that provides 24-hour-ahead forecast of solar flares, every 12 minutes. We report the method and prediction results of the system.

preprint2016arXiv

Above-the-loop-top Oscillation and Quasi-periodic Coronal Wave Generation in Solar Flares

Observations revealed that various kinds of oscillations are excited in solar flare regions. Quasi-periodic pulsations (QPPs) in the flare emissions are commonly observed in a wide range of wavelengths. Recent observations have found that fast-mode magnetohydrodynamic (MHD) waves are quasi-periodically emitted from some flaring sites (quasi-periodic propagating fast-mode magnetoacoustic waves; QPFs). Both of QPPs and QPFs imply a cyclic disturbance originating from the flaring sites. However, the physical mechanisms remain puzzling. By performing a set of two-dimensional MHD simulations of a solar flare, we discovered the local oscillation above the loops filled with evaporated plasma (above-the-loop-top region) and the generation of QPFs from such oscillating regions. Unlike all previous models for QPFs, our model includes essential physics for solar flares, such as magnetic reconnection, heat conduction, and chromospheric evaporation. We revealed that QPFs can be spontaneously excited by the above-the-loop-top oscillation. It was found that this oscillation is controlled by the backflow of the reconnection outflow. The new model revealed that flare loops and the above-the-loop-top region are full of shocks and waves, which is different from the previous expectations based on a standard flare model and previous simulations. In this paper, we will show the QPF generation process based on our new picture of flare loops and will briefly discuss a possible relationship between QPFs and QPPs. Our findings will change the current view of solar flares to a new view in which they are a very dynamic phenomenon with full of shocks and waves.

preprint2016arXiv

East Asian Observations of Low Latitude Aurora during the Carrington Magnetic Storm

The magnetic storm around 1859 September 2, caused by so-called Carrington flare, was the most intense in the history of modern scientific observations, and hence is considered to be the benchmark event for space weather. The magnetic storm caused worldwide observations of auroras even at very low latitudes such as Hawaii, Panama, or Santiago, and the available magnetic field measurement at Bombay, India, showed two peaks: the main was the Carrington event which occurred in day time in East Asia, and a second storm after the Carrington event which occurred at night in East Asia. In this paper, we present a result from surveys of aurora records in East Asia, which provides new information of the aurora activity of this important event. We found some new East Asian records of low latitude aurora observations caused by the storm which occurred after the Carrington event. The size of the aurora belt of the second peak of the Carrington magnetic storm was even wider than usual low-latitude aurora events.

preprint2016arXiv

Fractal Reconnection in Solar and Stellar Environments

Recent space based observations of the Sun revealed that magnetic reconnection is ubiquitous in the solar atmosphere, ranging from small scale reconnection (observed as nanoflares) to large scale one (observed as long duration flares or giant arcades). Often the magnetic reconnection events are associated with mass ejections or jets, which seem to be closely related to multiple plasmoid ejections from fractal current sheet. The bursty radio and hard X-ray emissions from flares also suggest the fractal reconnection and associated particle acceleration. We shall discuss recent observations and theories related to the plasmoid-induced-reconnection and the fractal reconnection in solar flares, and their implication to reconnection physics and particle acceleration. Recent findings of many superflares on solar type stars that has extended the applicability of the fractal reconnection model of solar flares to much a wider parameter space suitable for stellar flares are also discussed.

preprint2016arXiv

Observational Evidence of Particle Acceleration Associated with Plasmoid Motions

We report a strong association between the particle acceleration and plasma motions found in the 2010 August 18 solar flare. The plasma motions are tracked in the extreme-ultraviolet (EUV) images taken by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory and the Extreme UltraViolet Imager (EUVI) on the Solar Terrestrial Relation Observatory spacecraft Ahead, and the signature of particle acceleration was investigated by using Nobeyama Radioheliograph data. In our previous paper, we reported that in EUV images many plasma blobs appeared in the current sheet above the flare arcade. They were ejected bidirectionally along the current sheet, and the blobs that were ejected sunward collided with the flare arcade. Some of them collided or merged with each other before they were ejected from the current sheet. We discovered impulsive radio bursts associated with such plasma motions (ejection, coalescence, and collision with the post flare loops). The radio bursts are considered to be the gyrosynchrotron radiation by nonthermal high energy electrons. In addition, the stereoscopic observation by AIA and EUVI suggests that plasma blobs had a three-dimensionally elongated structure. We consider that the plasma blobs were three-dimensional plasmoids (i.e. flux ropes) moving in a current sheet. We believe that our observation provides clear evidence of particle acceleration associated with the plasmoid motions. We discuss possible acceleration mechanisms on the basis of our results.

preprint2016arXiv

Scaling Relations in Coronal Mass Ejections and Energetic Proton Events associated with Solar Superflares

In order to discuss the potential impact of solar 'superflares' on space weather, we investigated statistical relations among energetic proton peak flux with energy higher than $ 10 \rm MeV$ ($F_p$), CME speed near the Sun ($V_{CME}$) obtained by {\it SOHO}/LASCO coronagraph and flare soft X-ray peak flux in 1-8Åband ($F_{SXR}$) during 110 major solar proton events (SPEs) recorded from 1996 to 2014. The linear regression fit results in the scaling relations $V_{CME} \propto F_{SXR}^α$, $F_p\propto F_{SXR}^β$ and $F_p\propto V_{CME}^γ$ with $α= 0.30\pm 0.04$, $β= 1.19 \pm 0.08$ and $γ= 4.35 \pm 0.50$, respectively. On the basis of simple physical assumptions, on the other hand, we derive scaling relations expressing CME mass ($M_{CME}$), CME speed and energetic proton flux in terms of total flare energy ($E_{flare}$) as, $M_{CME}\propto E_{flare}^{2/3}$, $V_{CME}\propto E_{flare}^{1/6}$ and $F_{p}\propto E_{flare}^{5/6}\propto V_{CME}^5$, respectively. We then combine the derived scaling relations with observation, and estimated the upper limit of $V_{CME}$ and $F_p$ to be associated with possible solar superflares.

preprint2016arXiv

Statistical properties of superflares on solar-type stars with Kepler data

Superflares are flares that release total energy 10$\sim$10$^{4}$ times greater than that of the biggest solar flares with energy of $\sim$10$^{32}$ erg. We searched superflares on solar-type stars (G-type main sequence stars) using the Kepler 30-min (long) and 1-min (short) cadence data. We found more than 1500 superflares on 279 stars from 30-min cadence data (Q0-6) and 187 superflares on 23 stars from 1-min cadence data (Q0-17). The bolometric energy of detected superflares ranges from the order of 10$^{32}$ erg to 10$^{36}$ erg. Using these data, we found that the occurrence frequency ($dN/dE$) of superflares is expressed as a power-law function of flare energy ($E$) with the index of -1.5 for $10^{33}<E<10^{36}$erg. Most of the superflare stars show quasi-periodic light variations with the amplitude of a few percent, which can be explained by the rotation of the star with large starspots. The bolometric energy released by flares is consistent with the magnetic energy stored around such large starspots. Furthermore, our analyses indicate that the occurrence frequency of superflares depends on the rotation period, and that the flare frequency increases as the rotation period decreases. However, the energy of the largest flares observed in a given period bin does not show any clear correlation with the rotation period. We also found that the duration of superflares increases with the flare energy as $E^{0.39+/-0.03}$. This can be explained if we assume the time-scale of flares is determined by the Alfv$\acute{\rm{e}}$n time.

preprint2016arXiv

Validation of a Scaling Law for the Coronal Magnetic Field Strengths and Loop Lengths of Solar and Stellar Flares

Shibata & Yokoyama (1999, 2002) proposed a method of estimating the coronal magnetic field strengths ($B$) and magnetic loop lengths ($L$) of solar and stellar flares, on the basis of magnetohydrodynamic simulations of the magnetic reconnection model. Using the scaling law provided by Shibata & Yokoyama (1999, 2002), $B$ and $L$ are obtained as functions of the emission measure ($EM=n^2L^3$) and temperature ($T$) at the flare peak. Here, $n$ is the coronal electron density of the flares. This scaling law enables the estimation of $B$ and $L$ for unresolved stellar flares from the observable physical quantities $EM$ and $T$, which is helpful for studying stellar surface activities. To apply this scaling law to stellar flares, we discuss its validity for spatially resolved solar flares. $EM$ and $T$ were calculated from GOES soft X-ray flux data, and $B$ and $L$ are theoretically estimated using the scaling law. For the same flare events, $B$ and $L$ were also observationally estimated with images taken by Solar Dynamics Observatory (SDO)/ Helioseismic and Magnetic Imager (HMI) Magnetogram and Atmospheric Imaging Assembly (AIA) 94Å pass band. As expected, a positive correlation was found between the theoretically and observationally estimated values. We interpret this result as indirect evidence that flares are caused by magnetic reconnection. Moreover, this analysis makes us confident in the validity of applying this scaling law to stellar flares as well as solar flares.

preprint2015arXiv

A Theoretical Model of Pinching Current Sheet in Low-beta Plasmas

Magnetic reconnection is an important physical process in various explosive phenomena in the universe. In the previous studies, it was found that fast re- connection takes place when the thickness of a current sheet becomes on the order of a microscopic length such as the ion larmor radius or the ion inertial length. In this study, we investigated the pinching process of a current sheet by the Lorentz force in a low-β plasma using one-dimensional magnetohydrodynam- ics (MHD) simulations. It is known that there is an exact self-similar solution for this problem that neglects gas pressure. We compared the non-linear MHD dynamics with the analytic self-similar solution. From the MHD simulations, we found that with the gas pressure included the implosion process deviates from the analytic self-similar solution as t {\rightarrow} t 0, where t 0 is the explosion time when the thickness of a current sheet of the analytic solution becomes 0. We also found a pair of MHD fast-mode shocks are generated and propagate after the formation of the pinched current sheet as t {\rightarrow} t 0 . On the basis of the Rankine-Hugoniot relations, we derived the scaling law of the physical quantities with respect to the initial plasma beta in the pinched current sheet. Our study could help us to estimate the physical quantities in the pinched current sheet formed in a low-\b{eta} plasma.

preprint2015arXiv

High Dispersion Spectroscopy of Solar-type Superflare Stars. I. Temperature, Surface Gravity, Metallicity, and $v \sin i$

We conducted high dispersion spectroscopic observations of 50 superflare stars with Subaru/HDS, and measured the stellar parameters of them. These 50 targets were selected from the solar-type (G-type main sequence) superflare stars that we had discovered from the Kepler photometric data. As a result of these spectroscopic observations, we found that more than half (34 stars) of our 50 targets have no evidence of binary system. We then estimated effective temperature ($T_{\rm{eff}}$), surface gravity ($\log g$), metallicity ([Fe/H]), and projected rotational velocity ($v\sin i$) of these 34 superflare stars on the basis of our spectroscopic data. The accuracy of our estimations is higher than that of Kepler Input Catalog (KIC) values, and the differences between our values and KIC values ($(ΔT_{\rm{eff}})_{\rm{rms}} \sim 219$K, $(Δ\log g)_{\rm{rms}} \sim 0.37$ dex, and $(Δ\rm{[Fe/H]})_{\rm{rms}} \sim 0.46$ dex) are comparable to the large uncertainties and systematic differences of KIC values reported by the previous researches. We confirmed that the estimated $T_{\rm{eff}}$ and $\log g$ values of the 34 superflare stars are roughly in the range of solar-type stars. In particular, these parameters and the brightness variation period ($P_{0}$) of 9 stars are in the range of "Sun-like" stars ($5600\leq T_{\rm{eff}}\leq 6000$K, $\log g\geq$4.0, and $P_{0}>$10 days). Five of the 34 target stars are fast rotators ($v \sin i \geq 10$km s$^{-1}$), while 22 stars have relatively low $v \sin i$ values ($v \sin i<5$km s$^{-1}$). These results suggest that stars whose spectroscopic properties similar to the Sun can have superflares, and this supports the hypothesis that the Sun might cause a superflare.

preprint2015arXiv

High Dispersion Spectroscopy of Solar-type Superflare Stars. III. Lithium Abundances

We report on the abundance analysis of Li in solar-type (G-type main sequence) superflare stars which were found by the analysis of Kepler photometric data. Li is a key element to understand the evolution of the stellar convection zone which reflects the age of solar-type stars. We performed the high dispersion spectroscopy of solar-type superflare stars with Subaru/HDS, and confirmed that 34 stars show no evidence of binarity in our previous study. In this study, we derived the Li abundances of these 34 objects. We investigate correlations of Li abundance with stellar atmospheric parameters, rotational velocity, and superflare activities to understand the nature of superflare stars and the possibility of the nucleosynthesis of Li by superflares. We confirm the large dispersion in the Li abundance, and the correlation with stellar parameters is not seen. As compared with the Li abundance in Hyades cluster which is younger than the Sun, it is suggested that half of the observed stars are younger than Hyades cluster. The measured value of $v \sin i$ (projected rotational velocity) supports those objects are younger than the Sun. However, there are some objects which show the low Li abundance and slowly rotate on the basis of the estimated $v \sin i$ and $P$ (period of brightness variation). This result indicates that the superflare stars are not only young stars but also old stars like our Sun. In our observations, we could not find any evidence of Li productions by superflares. Further research on Li isotope abundances of superflare stars would clarify the Li production by stellar flares.

preprint2015arXiv

Magnetohydrodynamic shocks in and above post-flare loops: two-dimensional simulation and a simplified model

Solar flares are an explosive phenomenon, where super-sonic flows and shocks are expected in and above the post-flare loops. To understand the dynamics of post-flare loops, a two-dimensional magnetohydrodynamic (2D MHD) simulation of a solar flare has been carried out. We found new shock structures in and above the post-flare loops, which were not resolved in the previous work by Yokoyama and Shibata 2001. To study the dynamics of flows along the reconnected magnetic field, kinematics and energetics of the plasma are investigated along selected field lines. It is found that shocks are crucial to determine the thermal and flow structures in the post-flare loops. On the basis of the 2D MHD simulation, we have developed a new post-flare loop model which we call the pseudo-2D MHD model. The model is based on the 1D MHD equations, where all the variables depend on one space dimension and all the three components of the magnetic and velocity fields are considered. Our pseudo-2D model includes many features of the multi-dimensional MHD processes related to magnetic reconnection (particularly MHD shocks), which the previous 1D hydrodynamic models are not able to include. We compare the shock formation and energetics of a specific field line in the 2D calculation with those in our pseudo-2D MHD model, and we found that they give similar results. This model will allow us to study the evolution of the post-flare loops in a wide parameter space without expensive computational cost and without neglecting important physics associated with magnetic reconnection.

preprint2015arXiv

Numerical Study on Emergence of Kinked Flux Tube for Understanding of Possible Origin of Delta-spot Regions

We carried out a magnetohydrodynamics simulation where a subsurface twisted kink-unstable flux tube emerges from the solar interior to the corona. Unlike the previous expectations based on the bodily emergence of a knotted tube, we found that the kinked tube can spontaneously form a complex quadrupole structure at the photosphere. Due to the development of the kink instability before the emergence, the magnetic twist at the kinked apex of the tube is greatly reduced, although the other parts of the tube is still strongly twisted. This leads to the formation of a complex quadrupole structure: a pair of the coherent, strongly twisted spots and a narrow complex bipolar pair between it. The quadrupole is formed by the submergence of a portion of emerged magnetic fields. This result is relevant for understanding of the origin of the complex multipolar $δ$-spot regions that have a strong magnetic shear and emerge with polarity orientations not following Hale-Nicholson and Joy Laws.

preprint2015arXiv

Prominence Activation by Coronal Fast Mode Shock

An X5.4 class flare occurred in active region (AR) NOAA11429 on 2012 March 7. The flare was associated with very fast coronal mass ejection (CME) with its velocity of over 2500 km/s. In the images taken with STEREO-B/COR1, a dome-like disturbance was seen to detach from expanding CME bubble and propagated further. A Type-II radio burst was also observed at the same time. On the other hand, in EUV images obtained by SDO/AIA, expanding dome-like structure and its foot print propagating to the north were observed. The foot print propagated with its average speed of about 670 km/s and hit a prominence located at the north pole and activated it. While the activation, the prominence was strongly brightened. On the basis of some observational evidence, we concluded that the foot print in AIA images and the ones in COR1 images are the same, that is MHD fast mode shock front. With the help of a linear theory, the fast mode mach number of the coronal shock is estimated to be between 1.11 and 1.29 using the initial velocity of the activated prominence. Also, the plasma compression ratio of the shock is enhanced to be between 1.18 and 2.11 in the prominence material, which we consider to be the reason of the strong brightening of the activated prominence. The applicability of linear theory to the shock problem is tested with nonlinear MHD simulation.

preprint2015arXiv

Statistical properties of superflares on solar-type stars based on 1-min cadence data

We searched for superflares on solar-type stars using Kepler data with 1 min sampling in order to detect superflares with short duration. We found 187 superflares on 23 solar-type stars whose bolometric energy ranges from the order of $10^{32}$ erg to $10^{36}$ erg. Some superflares show multiple peaks with the peak separation of the order of $100$-$1000$ seconds which is comparable to the periods of quasi-periodic pulsations in solar and stellar flares. Using these new data combined with the results from the data with 30 min sampling, we found the occurrence frequency ($dN/dE$) of superflares as a function of flare energy ($E$) shows the power-law distribution ($dN/dE \propto E^{-α}$) with $α\sim -1.5$ for $10^{33}<E<10^{36}$ erg which is consistent with the previous results. The average occurrence rate of superflares with the energy of $10^{33}$ erg which is equivalent to X100 solar flares is about once in 500-600 years. The upper limit of energy released by superflares is basically comparable to a fraction of the magnetic energy stored near starspots which is estimated from the photometry. We also found that the duration of superflares ($τ$) increases with the flare energy ($E$) as $τ\propto E^{0.39\pm 0.03}$. This can be explained if we assume the time-scale of flares is determined by the Alfv$\acute{\rm e}$n time.

preprint2015arXiv

UFCORIN: A Fully Automated Predictor of Solar Flares in GOES X-Ray Flux

We have developed UFCORIN, a platform for studying and automating space weather prediction. Using our system we have tested 6,160 different combinations of SDO/HMI data as input data, and simulated the prediction of GOES X-ray flux for 2 years (2011-2012) with one-hour cadence. We have found that direct comparison of the true skill statistics (TSS) from small cross-validation sets is ill-posed, and used the standard scores ($z$) of the TSS to compare the performance of the various prediction strategies. The $z$ of a strategy is a stochastic variable of the stochastically-chosen cross-validation dataset, and the $z$ for the three strategies best at predicting X, $\geq$M and $\geq$C class flares are better than the average $z$ of the 6,160 strategies by 2.3$σ$, 2.1$σ$, 3.8$σ$ confidence levels, respectively. The best three TSS values were $0.75\pm0.07$, $0.48\pm0.02$, and $0.56\pm0.04$, respectively.

preprint2014arXiv

A Chain of Winking (Oscillating) Filaments Triggered by an Invisible Extreme-Ultraviolet Wave

Winking (oscillating) filaments have been observed for many years. However, observations of successive winking filaments in one event have not been reported yet. In this paper, we present the observations of a chain of winking filaments and a subsequent jet that are observed right after the X2.1 are in AR11283. The event also produced an Extreme-ultraviolet (EUV) wave that has two components: upward dome-like wave (850 km/s) and lateral surface wave (554 km/s) which was very weak (or invisible) in imaging observations. By analyzing the temporal and spatial relationships between the oscillating filaments and the EUV waves, we propose that all the winking filaments and the jet were triggered by the weak (or invisible) lateral surface EUV wave. The oscillation of the filaments last for two or three cycles, and their periods, Doppler velocity amplitudes, and damping times are 11 - 22 minutes, 6 - 14 km/s, and 25 - 60 minutes, respectively. We further estimate the radial component magnetic field and the maximum kinetic energy of the filaments, and they are 5 - 10 Gauss and ~ 10^19 J, respectively. The estimated maximum kinetic energy is comparable to the minimum energy of ordinary EUV waves, suggesting that EUV waves can efficiently launch filament oscillations on their path. Based on our analysis results, we conclude that the EUV wave is a good agent for triggering and connecting successive but separated solar activities in the solar atmosphere, and it is also important for producing solar sympathetic eruptions.

preprint2014arXiv

High Dispersion Spectroscopy of Solar-type Superflare Stars. II. Stellar Rotation, Starspots, and Chromospheric Activities

We conducted high dispersion spectroscopic observations of 50 superflare stars with Subaru/HDS. These 50 stars were selected from the solar-type superflare stars that we had discovered from the Kepler data. More than half (34 stars) of these 50 target superflare stars show no evidence of binarity, and we estimated stellar parameters of these 34 stars in our previous study (Notsu et al. 2015, hereafter called Paper I). According to our previous studies using Kepler data, superflare stars show quasi-periodic brightness variations whose amplitude (0.1-10\%) is much larger than that of the solar brightness variations (0.01-0.1\%) caused by the existence of sunspots on the rotating solar surface. In this study, we investigated whether these quasi-periodic brightness variations of superflare stars are explained by the rotation of a star with fairly large starspots, by using stellar parameters derived in Paper I. First, we confirmed that the value of the projected rotational velocity $v \sin i$ is consistent with the rotational velocity estimated from the period of the brightness variation. Next, we measured the intensity of Ca II infrared triplet lines and H$α$ line, good indicators of the stellar chromospheric activity, and compared them with other stellar properties. The intensity of Ca II infrared triplet lines indicates that the mean magnetic field strength ($\langle fB\rangle$) of the target superflare stars can be higher than that of the Sun. A correlation between the amplitude of the brightness variation and the intensity of Ca II triplet line was found. All the targets expected to have large starspots because of their large amplitude of the brightness variation show high chromospheric activities compared to the Sun. These results support that the brightness variation of superflare stars is due to the rotation with large starspots.

preprint2014arXiv

Superflare occurrence and energies on G-, K- and M-type dwarfs

Kepler data from G-, K- and M-type stars are used to study conditions that lead to superflares with energies above $10^{34} {\rm erg}$. From the 117,661 stars included, 380 show superflares with a total of 1690 such events. We study whether parameters, like effective temperature or the rotation rate, have any effect on the superflare occurrence rate or energy. With increasing effective temperature we observe a decrease in the superflare rate, which is analogous to the previous findings of a decrease in dynamo activity with increasing effective temperature. For slowly rotating stars, we find a quadratic increase of the mean occurrence rate with the rotation rate up to a critical point, after which the rate decreases linearly. Motivated by standard dynamo theory, we study the behavior of the relative starspot coverage, approximated as the relative brightness variation. For faster rotating stars, an increased fraction of stars shows higher spot coverage, which leads to higher superflare rates. A turbulent dynamo is used to study the dependence of the Ohmic dissipation as a proxy of the flare energy on the differential rotation or shear rate. The resulting statistics of the dissipation energy as a function of dynamo number is similar to the observed flare statistics as a function of the inverse Rossby number and shows similarly strong fluctuations. This supports the idea that superflares might well be possible for solar-type G stars.

preprint2014arXiv

Two Sun-like Superflare Stars Rotating as Slow as the Sun

We report on the results of high dispersion spectroscopy of two `superflare stars', KIC 9766237, and KIC 9944137 with Subaru/HDS. Superflare stars are G-type main sequence stars, but show gigantic flares compared to the Sun, which have been recently discovered in the data obtained with the Kepler spacecraft. Though most of these stars are thought to have a rotation period shorter than 10 days on the basis of photometric variabilities, the two targets of the present paper are estimated to have a rotation period of 21.8 d, and 25.3 d. Our spectroscopic results clarified that these stars have stellar parameters similar to those of the Sun in terms of the effective temperature, surface gravity, and metallicity. The projected rotational velocities derived by us are consistent with the photometric rotation period, indicating a fairy high inclination angle. The average strength of the magnetic field on the surface of these stars are estimated to be 1-20 G, by using the absorption line of Ca II 8542. We could not detect any hint of binarity in our spectra, although more data are needed to firmly rule out the presence of an unseen low-mass companion. These results claim that the spectroscopic properties of these superflare stars are very close to those of the Sun, and support the hypothesis that the Sun might cause a superflare.

preprint2013arXiv

Diffraction, Refraction, and Reflection of An Extreme-Ultraviolet Wave Observed during Its Interactions with Remote Active Regions

We present observations of the diffraction, refraction, and reflection of a global extreme-ultraviolet (EUV) wave propagating in the solar corona. These intriguing phenomena are observed when the wave interacts with two remote active regions, and they together exhibit the wave property of this EUV wave. When the wave approached AR11465, it became weaker and finally disappeared in the active region, but a few minutes latter a new wavefront appeared behind the active region, and it was not concentric with the incoming wave. In addition, a reflected wave was also observed simultaneously on the wave incoming side. When the wave approached AR11459, it transmitted through the active region directly and without reflection. The formation of the new wavefront and the transmission could be explained with diffraction and refraction effects, respectively. We propose that the different behaviors observed during the interactions may caused by different speed gradients at the boundaries of the two active regions. For the origin of the EUV wave, we find that it formed ahead of a group of expanding loops a few minutes after the start of the loops' expansion, which represents the initiation of the associated coronal mass ejection (CME). Based on these results, we conclude that the EUV wave should be a nonlinear magnetosonic wave or shock driven by the associated CME, which propagated faster than the ambient fast-mode speed and gradually slowed down to an ordinary linear wave. Our observations support the hybrid model that includes both fast wave and slow non-wave components.

preprint2013arXiv

High Dispersion Spectroscopy of the Superflare Star KIC6934317

We conducted the high-resolution spectroscopic observation with Subaru/HDS for a G-type star (KIC6934317). We selected this star from the data of the Kepler spacecraft. This star produces a lot of superflares, and the total energy of the largest superflare on this star is ~10^3 times larger (~2.2*10^35 erg) than that of the most energetic flare on the Sun (~10^32 erg). The core depth and emission flux of Ca II infrared triplet lines and Ha line show high chromospheric activity in this star, in spite of its low lithium abundance and the small amplitude of the rotational modulation. Using the empirical relations between emission flux of chromospheric lines and X-ray flux, this star is considered to show much higher coronal activity than that of the Sun. It probably has large starspots which can store a large amount of magnetic energy enough to give rise to superflares. We also estimated the stellar parameters, such as effective temperature, surface gravity, metallicity, projected rotational velocity (v sin_i), and radial velocity. KIC6934317 is then confirmed to be an early G-type main sequence star. The value of v sin_i is estimated to be ~1.91 km s^-1. In contrast, the rotational velocity is calculated to be ~20 km s^-1 by using the period of the brightness variation as the rotation period. This difference can be explained by its small inclination angle (nearly pole-on). The small inclination angle is also supported by the contrast between the large superflare amplitude and the small stellar brightness variation amplitude. The lithium abundance and isochrones implies that the age of this star is more than about a few Gyr, though a problem why this star with such an age has a strong activity remains unsolved.

preprint2013arXiv

Numerical Simulations of Solar Chromospheric Jets Associated with Emerging Flux

We studied the acceleration mechanisms of chromospheric jets associated with emerging flux using a two dimensional magnetohydrodynamic (MHD) simulation. We found that slow mode shock waves generated by magnetic reconnection in the chromosphere and the photosphere play key roles in the acceleration mechanisms of chromospheric jets. An important parameter is the height of magnetic reconnection. When magnetic reconnection takes place near the photosphere, the reconnection out- flow collides with the region where the plasma beta is much larger than unity. Then the plasma moves along a magnetic field. This motion generates a slow mode wave. The slow mode wave develops to a strong slow shock as it propagates upward. When the slow shock crosses the transition region, the transition region is lifted up. As a result, we obtain a chromospheric jet as the lifted transition region. When magnetic reconnection takes place in the upper chromosphere, the chromospheric plasma is accelerated due to the combination of the Lorentz force and the whip-like motion of magnetic field. We found that the chromospheric plasma is further accelerated through the interaction between the transition region (steep density gradient) and a slow shock emanating from the reconnection point. This is an MHD effect which has not been discussed before.

preprint2013arXiv

Superflares on Solar Type Stars Observed with Kepler I. Statistical Properties of Superflares

By extending our previous study by Maehara et al. (2012), we searched for superflares on G-type dwarfs (solar type stars) using Kepler data for a longer period (500 days) than that (120 days) in our previous study. As a result, we found 1547 superflares on 279 G-type dwarfs, which are much more than previous 365 superflares on 148 stars. Using these new data, we studied the statistical properties of occurrence frequency of superflares, and basically confirmed the previous results, i.e., the occurrence frequency (dN/dE) of superflares vs flare energy (E) shows power-law distribution with dN/dE \propto E^{-α}, where α~ 2. It is interesting that this distribution is roughly on the same line as that for solar flares. In the case of the Sun-like stars (with surface temperature 5600-6000K and slowly rotating with a period longer than 10 days), the occurrence frequency of superflares with energy of 10^34 -10^35 erg is once in 800-5000 years. We also studied long term (500 days) stellar brightness variation of these superflare stars, and found that in some G-type dwarfs the occurrence frequency of superflares was extremely high, ~ 57 superflares in 500 days (i.e., once in 10 days). In the case of Sun-like stars, the most active stars show the frequency of one superflares (with 10^34 erg) in 100 days. There is evidence that these superflares have extremely large starspots with a size about 10 times larger than that of the largest sunspot. We argue that the physical origin of extremely high occurrence frequency of superflares in these stars may be attributed to the existence of extremely large starspots.

preprint2013arXiv

Superflares on Solar-Type Stars Observed with Kepler II. Photometric Variability of Superflare-Generating Stars : A Signature of Stellar Rotation and Starspots

We performed simple spot-model calculations for quasi-periodic brightness variations of solar-type stars showing superflares, by using Kepler photometric data. Most of superflare stars show quasi-periodic brightness modulations with the typical period of one to a few tens of days. Our results indicate that these brightness variations of superflare stars can be explained by the rotation of the star with fairly large starspots. Using the result of the period analysis, we investigated the relation between the energy and frequency of superflares and the rotation period. Stars with relatively slower rotation rates can still produce flares that are as energetic as those of more rapidly rotating stars, although the average flare frequency is lower for more slowly rotating stars. We found that the energy of superflares are related to the total coverage of starspots. The correlation between the spot coverage and the flare energy in superflares is similar to that in solar flares. These results suggest that the energy of superflares can be explained by the magnetic energy stored around starspots.

preprint2013arXiv

The Role of a Flux Rope Ejection in Three-dimensional Magnetohydrodynamic Simulation of a Solar Flare

We investigated the dynamic evolution of a 3-dimensional (3D) flux rope eruption and magnetic reconnection process in a solar flare, by simply extending 2-dimensional (2D) resistive magnetohydrodynamic simulation model of solar flares with low $β$ plasma to 3D model. We succeeded in reproducing a current sheet and bi-directional reconnection outflows just below the flux rope during the eruption in our 3D simulations. We calculated four cases of a strongly twisted flux rope and a weakly twisted flux rope in 2D and 3D simulations. The time evolution of a weakly twisted flux rope in 3D simulation shows similar behaviors to 2D simulation, while a strongly twisted flux rope in 3D simulation shows clearly different time evolution from 2D simulation except for the initial phase evolution. The ejection speeds of both strongly and weakly twisted flux ropes in 3D simulations are larger than 2D simulations, and the reconnection rates in 3D cases are also larger than 2D cases. This indicates a positive feedback between the ejection speed of a flux rope and the reconnection rate even in the 3D simulation, and we conclude that the plasmoid-induced reconnection model can be applied to 3D. We also found that small scale plasmoids are formed inside a current sheet and make it turbulent. These small scale plasmoid ejections has role in locally increasing reconnection rate intermittently as observed in solar flares, coupled with a global eruption of a flux rope.

preprint2012arXiv

Effect of Interacting Rarefaction Waves on Relativistically Hot Jets

The effect of rarefaction acceleration on the propagation dynamics and structure of relativistically hot jets is studied through relativistic hydrodynamic simulations. We emphasize the nonlinear interaction of rarefaction waves excited at the interface between a cylindrical jet and the surrounding medium. From simplified one-dimensional models with radial jet structure, we find that a decrease in the relativistic pressure due to the interacting rarefaction waves in the central zone of the jet transiently yields a more powerful boost of the bulk jet than that expected from single rarefaction acceleration. This leads to a cyclic in-situ energy conversion between thermal and bulk kinetic energies which induces radial oscillating motion of the jet. The oscillation timescale is characterized by the initial pressure ratio of the jet to the ambient medium, and follows a simple scaling relation $τ_{\rm oscillation} \propto (P_{\rm jet,0}/P_{\rm amb,0})^{1/2}$. It is confirmed from extended two-dimensional simulations that this radial oscillating motion in the one-dimensional system manifests as modulation of the structure of the jet in a more realistic situation where a relativistically hot jet propagates through an ambient medium. It is found that when the ambient medium has a power law pressure distribution, the size of the reconfinement region along the propagation direction of the jet in the modulation structure $λ$ evolves according to a self-similar relation $λ\propto t^{α/2}$ where $α$ is the power-law index of the pressure distribution.

preprint2011arXiv

A Study on Red Asymmetry of H-alpha Flare Ribbons Using Narrowband Filtergram in the 2001 April 10 Solar Flare

We report a detailed examination of the "red asymmetry" of H-alpha emission line seen during the 2001 April 10 solar flare by using a narrowband filtergram. We investigated the temporal evolution and the spatial distribution of the red asymmetry by using the H-alpha data taken with the 60cm Domeless Solar Telescope at Hida Observatory, Kyoto University. We confirmed that the red asymmetry clearly appeared all over the flare ribbons, and the strong red asymmetry is located on the outer narrow edges of the flare ribbons, with the width of about 1.5" - 3.0" (1000 - 2000 km), where the strong energy releases occur. Moreover, we found that the red asymmetry, which also gives a measure of the Doppler shift of the H-alpha emission line concentrates on a certain value, not depending on the intensity of the H-alpha kernels. This implies not only that the temporal evolutions of the red asymmetry and those of the intensity are not in synchronous in each flare kernel, but also that the peak asymmetry (or velocity of the chromospheric condensation) of individual kernel is not a strong function of their peak intensity.

preprint2011arXiv

Dependence of the Magnetic Energy of Solar Active Regions on the Twist Intensity of the Initial Flux Tubes

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 interchange mode of magnetic buoyancy instability. As for the weakest twist case, the tube exhibits an elongated photospheric structure and never rises into the corona. The formation of the photospheric structure is due to inward magnetic tension force of the azimuthal field component of the rising flux tube (i.e., tube's twist). When the twist is weak, azimuthal field cannot hold the tube's coherency, and the tube extends laterally at the subadiabatic surface. In addition, we newly find 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 deviates from this relation, because these tubes store their magnetic energy in the photospheric structures.

preprint2011arXiv

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

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.

preprint2011arXiv

Numerical Simulations of the Magnetic Rayleigh-Taylor Instability in the Kippenhahn-Schlüter Prominence Model

The launch of the Hinode satellite has allowed unprecedented high resolution, stable images of solar quiescent prominences to be taken over extended periods of time. These new images led to the discovery of dark upflows that propagated from the base of prominences developing highly turbulent profiles. As yet, how these flows are driven is not fully understood. To study the physics behind this phenomena we use 3-D magnetohydrodynamic (MHD) simulations to investigate the nonlinear stability of the Kippenhahn-Shlüter prominence model to the magnetic Rayleigh-Taylor instability. The model simulates the rise of a buoyant tube inside a quiescent prominence, where the upper boundary between the tube and prominence model is perturbed to excite the interchange of magnetic field lines. We found upflows of constant velocity (maximum found 6\,km\,s$^{-1}$) and a maximum plume width $\approx 1500$\,km which propagate through a height of approximately 6\,Mm in the no guide field case. The case with the strong guide field (initially $B_y=2B_x$) results in a large plume that rises though the prominence model at $\sim 5$\,km\,s$^{-1}$ with width $\sim 900$\,km (resulting in width of 2400\,km when viewed along the axis of the prominence) reaching a height of $\sim 3.1$\,Mm. In both cases nonlinear processes were important for determining plume dynamics.

preprint2011arXiv

Simultaneous Observation of Reconnection Inflow and Outflow Associated with the 2010 August 18 Solar Flare

We report the simultaneous extreme ultraviolet observation of magnetic reconnection inflow and outflow in a flare on 2010 August 18 observed by the Atmospheric Imaging Assembly on board the {\it Solar Dynamics Observatory}. We found that during the rise phase of the flare, some plasma blobs appeared in the sheet structure above the hot loops. The plasma blobs were ejected bidirectionally along the sheet structure (outflow), at the same time as the threads visible in extreme ultraviolet images moved toward the sheet structure (inflow). The upward and downward ejection velocities are 220-460 km s$^{-1}$ and 250-280 km s$^{-1}$, respectively. The inflow speed changed from 90 km s$^{-1}$ to 12 km s$^{-1}$ in 5 minutes. By using these velocities, we estimated the nondimensional reconnection rate, which we found to vary during this period from 0.20 to 0.055. We also found that the plasma blobs in the sheet structure collided or merged with each other before they were ejected from the sheet structure. We hypothesize that the sheet structure is the current sheet and that these plasma blobs are plasmoids or magnetic islands, which could be important for understanding the dynamics of the reconnection region.

preprint2011arXiv

Special Relativistic Magnetohydrodynamic Simulation of Two-Component Outflow Powered by Magnetic Explosion on Compact Stars

The nonlinear dynamics of outflows driven by magnetic explosion on the surface of a compact star is investigated through special relativistic magnetohydrodynamic simulations. We adopt, as the initial equilibrium state, a spherical stellar object embedded in hydrostatic plasma which has a density $ρ(r) \propto r^{- α}$ and is threaded by a dipole magnetic field. The injection of magnetic energy at the surface of compact star breaks the equilibrium and triggers a two-component outflow. At the early evolutionary stage, the magnetic pressure increases rapidly around the stellar surface, initiating a magnetically driven outflow. A strong forward shock driven outflow is then excited. The expansion velocity of the magnetically driven outflow is characterized by the Alfvén velocity on the stellar surface, and follows a simple scaling relation $v_{\rm mag} \propto {v_{\rm A}}^{1/2}$. When the initial density profile declines steeply with radius, the strong shock is accelerated self-similarly to relativistic velocity ahead of the magnetically driven component. We find that it evolves according to a self-similar relation $Γ_{\rm sh} \propto r_{\rm sh}$, where $Γ_{\rm sh}$ is the Lorentz factor of the plasma measured at the shock surface $r_{\rm sh}$. Purely hydrodynamic process would be responsible for the acceleration mechanism of the shock driven outflow. Our two-component outflow model, which is the natural outcome of the magnetic explosion, can provide a better understanding of the magnetic active phenomena on various magnetized compact stars.

preprint2011arXiv

Statistical Study on the Nature of Solar Flux Emergence

We studied 101 flux emergence events ranging from small ephemeral regions to large emerging flux regions which were observed with Hinode Solar Optical Telescope filtergram. We investigated how the total magnetic flux of the emergence event controls the nature of emergence. To determine the modes of emergences, horizontal velocity fields of global motion of the magnetic patches in the flux emerging sites were measured by the local correlation tracking. Between two main polarities of the large emerging flux regions with more than around 2 \times 10^19 Mx, there were the converging flows of anti-polarity magnetic patches. On the other hand, small ephemeral regions showed no converging flow but simple diverging pattern. When we looked into the detailed features in the emerging sites, irrespective of the total flux and the spatial size, all the emergence events were observed to consist of single or multiple elementary emergence unit(s). The typical size of unitary emergence is 4 Mm and consistent with the simulation results. From the statistical study of the flux emergence events, the maximum spatial distance between two main polarities, the magnetic flux growth rate and the mean separation speed were found to follow the power-law functions of the total magnetic flux with the indices of 0.27, 0.57, and -0.16, respectively. From the discussion on the observed power-law relations, we got a physical view of solar flux emergence that emerging magnetic fields float and evolve balancing to the surrounding turbulent atmosphere. Key words: Sun: magnetic fields - Sun: emerging flux - Sun: photosphere - Sun: chromosphere

preprint2010arXiv

Ca II K Spectral Study of an Emerging Flux Region using Domeless Solar Telescope in Hida Observatory

A cooperative observation with Hida observatory and Hinode satellite was performed on an emerging flux region. The successive Ca II K spectro-heliograms of the emerging flux region were taken by the Domeless Solar Telescope of Hida observatory. Hinode observed the emerging flux region with Ca II H and Fe I Stokes IQUV filtergrams. In this study, detailed dynamics and temporal evolution of the magnetic flux emergence was studied observationally. The event was first detected in the photospheric magnetic field signals. 3 minutes later, the horizontal expansion of the dark area was detected. And then, 7 minutes later than the horizontal expansion, the emerging loops were detected with the maximal rise speed of 2.1 km/s at chromospheric heights. The observed dynamics of emerging magnetic flux from the photosphere to the upper chromosphere is well consistent with the results of previous simulation works. The gradual rising phase of flux tubes with a weak magnetic strength was confirmed by our observation.

preprint2010arXiv

Fan-shaped jets in three dimensional reconnection simulation as a model of ubiquitous solar jets

Magnetic reconnection is a fundamental process in space and astrophysical plasmas in which oppositely directed magnetic fields changes its connectivity and eventually converts its energy into kinetic and thermal energy of the plasma. Recently, ubiquitous jets (for example, chromospheric anemone jets, penumbral microjets, umbral light bridge jets) have been observed by Solar Optical Telescope on board the satellite Hinode. These tiny and frequently occurring jets are considered to be a possible evidence of small-scale ubiquitous reconnection in the solar atmosphere. However, the details of three dimensional magnetic configuration are still not very clear. Here we propose a new model based on three dimensional simulations of magnetic reconnection using a typical current sheet magnetic configuration with a strong guide field. The most interesting feature is that the jets produced by the reconnection eventually move along the guide field lines. This model provides a fresh understanding of newly discovered ubiquitous jets and moreover a new observational basis for the theory of astrophysical magnetic reconnection.

preprint2010arXiv

Generation of Alfven Waves by Magnetic Reconnection

In this paper, results of 2.5-dimensional magnetohydrodynamical simulations are reported for the magnetic reconnection of non-perfectly antiparallel magnetic fields. The magnetic field has a component perpendicular to the computational plane, that is, guide field. The angle theta between magnetic field lines in two half regions is a key parameter in our simulations whereas the initial distribution of the plasma is assumed to be simple; density and pressure are uniform except for the current sheet region. Alfven waves are generated at the reconnection point and propagate along the reconnected field line. The energy fluxes of the Alfven waves and magneto-acoustic waves (slow mode and fast mode) generated by the magnetic reconnection are measured. Each flux shows the similar time evolution independent of theta. The percentage of the energies (time integral of energy fluxes) carried by the Alfven waves and magneto-acoustic waves to the released magnetic energy are calculated. The Alfven waves carry 38.9%, 36.0%, and 29.5% of the released magnetic energy at the maximum (theta=80^\circ) in the case of beta=0.1, 1, and 20 respectively, where beta is the plasma beta (the ratio of gas pressure to magnetic pressure). The magneto-acoustic waves carry 16.2% (theta=70^\circ), 25.9% (theta=60^\circ), and 75.0% (theta=180^\circ) of the energy at the maximum. Implications of these results for solar coronal heating and acceleration of high-speed solar wind are discussed.

preprint2010arXiv

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

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.

preprint2010arXiv

Nonlinear propagation of Alfven waves driven by observed photospheric motions: Application to the coronal heating and spicule formation

We have performed MHD simulations of Alfven wave propagation along an open flux tube in the solar atmosphere. In our numerical model, Alfven waves are generated by the photospheric granular motion. As the wave generator, we used a derived temporal spectrum of the photospheric granular motion from G-band movies of Hinode/SOT. It is shown that the total energy flux at the corona becomes larger and the transition region height becomes higher in the case when we use the observed spectrum rather than white/pink noise spectrum as the wave generator. This difference can be explained by the Alfven wave resonance between the photosphere and the transition region. After performing Fourier analysis on our numerical results, we have found that the region between the photosphere and the transition region becomes an Alfven wave resonant cavity. We have confirmed that there are at least three resonant frequencies, 1, 3 and 5 mHz, in our numerical model. Alfven wave resonance is one of the most effective mechanisms to explain the dynamics of the spicules and the sufficient energy flux to heat the corona.

preprint2010arXiv

Observations of Chromospheric Anemone Jets with Hinode SOT and Hida Ca II Spectroheliogram

We present the first simultaneous observations of chromospheric "anemone" jets in solar active regions with Hinode SOT Ca II H broadband filetergram and Ca II K spetroheliogram on the Domeless Solar Telescope (DST) at Hida Observatory. During the coordinated observation, 9 chromospheric anemone jets were simultaneously observed with the two instruments. These observations revealed three important features, i.e.: (1) the jets are generated in the lower chromosphere, (2) the length and lifetime of the jets are 0.4-5 Mm and 40-320 sec, (3) the apparent velocity of the jets with Hinode SOT are 3-24 km/s, while Ca II K3 component at the jets show blueshifts (in 5 events) in the range of 2- 6 km/s. The chromospheric anemone jets are associated with mixed polarity regions which are either small emerging flux regions or moving magnetic features. It is found that the Ca II K line often show red or blue asymmetry in K2/K1 component: the footpoint of the jets associated with emerging flux regions often show redshift (2-16 km/s), while the one with moving magnetic features show blueshift (around 5 km/s). Detailed analysis of magnetic evolution of the jet foaming regions revealed that the reconnection rate (or canceling rate) of the total magnetic flux at the footpoint of the jets are of order of 10^{16} Mx/s, and the resulting magnetic energy release rate (1.1-10) x 10^{24} erg/s, with the total energy release (1-13) x 10^{26} erg for the duration of the magnetic cancellations, 130s. These are comparable to the estimated total energy, 10^{26} erg, in a single chromospheric anemone jet. An observation-based physical model of the jet is presented. The relation between chromospheric anemone jets and Ellerman bombs is discussed.

preprint2010arXiv

Solar-type Magnetic Reconnection Model for Magnetar Giant Flare

We present a theoretical model describing magnetar giant flares on the basis of solar flare/coronal mass ejection theory. In our model, a preflare activity plays a crucial role in driving evaporating flows and supplying baryonic matters into the magnetosphere. The loaded baryonic matter, that is called "prominence", is then gradually uplifted via crustal cracking with maintaining a quasi-force-free equilibrium of the magnetosphere. Finally the prominence is erupted by the magnetic pressure force due to the loss of equilibrium triggered by the explosive magnetic reconnection. The giant flare should be induced as a final outcome of the prominence eruption accompanied by large-scale field reconfigurations. An essential difference between solar and magnetar flares is the control process of their evolutionary dynamics. The flaring activity on magnetars is mainly controlled by the radiative process unlike the solar flare governed by the electron conduction. It is highly suggestive that our model is accountable for the physical properties of the extraordinary giant flare observed on 2004 December 27 from SGR1806-20, including the source of baryonic matters loaded in the expanding ejecta observed after the giant burst.

preprint2010arXiv

Spicule Dynamics over Plage Region

We studied spicular jets over a plage area and derived their dynamic characteristics using Hinode Solar Optical Telescope (SOT) high-resolution images. The target plage region was near the west limb of the solar disk. This location permitted us to study the dynamics of spicular jets without the overlapping effect of spicular structures along the line of sight. In this work, to increase the ease with which we can identify spicules on the disk, we applied the image processing method `MadMax' developed by Koutchmy et al. (1989). It enhances fine, slender structures (like jets), over a diffuse background. We identified 169 spicules over the target plage. This sample permits us to derive statistically reliable results regarding spicular dynamics. The properties of plage spicules can be summarized as follows: (1) In a plage area, we clearly identified spicular jet features. (2) They were shorter in length than the quiet region limb spicules, and followed ballistic motion under constant deceleration. (3) The majority (80%) of the plage spicules showed the cycle of rise and retreat, while 10% of them faded out without a complete retreat phase. (4) The deceleration of the spicule was proportional to the velocity of ejection (i.e. the initial velocity).

preprint2009arXiv

Prominence Formation Associated with an Emerging Helical Flux Rope

The formation and evolution process and magnetic configuration of solar prominences remain unclear. In order to study the formation process of prominences, we examine continuous observations of a prominence in NOAA AR 10953 with the Solar Optical Telescope on the \emph{Hinode} satellite. As reported in our previous Letter, we find a signature suggesting that a helical flux rope emerges from below the photosphere under a pre-existing prominence. Here we investigate more detailed properties and photospheric indications of the emerging helical flux rope, and discuss their relationship to the formation of the prominence. Our main conclusions are: (1) A dark region with absence of strong vertical magnetic fields broadens and then narrows in Ca \textsc{ii} H-line filtergrams. This phenomenon is consistent with the emergence of the helical flux rope as photospheric counterparts. The size of the flux rope is roughly 30,000 km long and 10,000 km wide. The width is larger than that of the prominence. (2) No shear motion or converging flows are detected, but we find diverging flows such as mesogranules along the polarity inversion line. The presence of mesogranules may be related to the emergence of the helical flux rope. (3) The emerging helical flux rope reconnects with magnetic fields of the pre-existing prominence to stabilize the prominence for the next several days. We thus conjecture that prominence coronal magnetic fields emerge in the form of helical flux ropes that contribute to the formation and maintenance of the prominence.

preprint2008arXiv

Characteristics of Anemone Active Regions Appearing in Coronal Holes Observed with {\it Yohkoh} Soft X-ray Telescope

Coronal structure of active regions appearing in coronal holes is studied by using the data obtained with the Soft X-Ray Telescope (SXT) aboard {\it Yohkoh} from 1991 November to 1993 March. The following characteristics are found; Many of active regions appearing in coronal holes show a structure that looks like a ``sea-anemone''. Such active regions are called {\it anemone ARs}. About one-forth of all active regions that were observed with SXT from their births showed the anemone structure. For almost all the anemone ARs, the order of magnetic polarities is consistent with the Hale-Nicholson's polarity law. These anemone ARs also showed more or less east-west asymmetry in X-ray intensity distribution, such that the following (eastern) part of the ARs is brighter than its preceding (western) part. This, as well as the anemone shape itself, is consistent with the magnetic polarity distribution around the anemone ARs. These observations also suggest that an active region appearing in coronal holes has simpler (less sheared) and more preceding-spot-dominant magnetic structure than those appearing in other regions.

preprint2008arXiv

Evolution of Anemone AR NOAA 10798 and the Related Geo-Effective Flares and CMEs

We present a detailed examination of the features of the Active Region (AR) NOAA 10798. This AR generated coronal mass ejections (CMEs) that caused a large geomagnetic storm on 24 August 2005 with the minimum Dst index of -216 nT. We examined the evolution of the AR and the features on/near the solar surface and in the interplanetary space. The AR emerged in the middle of a small coronal hole, and formed a {\it sea anemone} like configuration. H$α$ filaments were formed in the AR, which have southward axial field. Three M-class flares were generated, and the first two that occurred on 22 August 2005 were followed by Halo-type CMEs. The speeds of the CMEs were fast, and recorded about 1200 and 2400 km s$^{-1}$, respectively. The second CME was especially fast, and caught up and interacted with the first (slower) CME during their travelings toward Earth. These acted synergically to generate an interplanetary disturbance with strong southward magnetic field of about -50 nT, which was followed by the large geomagnetic storm.

preprint2007arXiv

Dead Zone Formation and Nonsteady Hyperaccretion in Collapsar Disks : A Possible Origin of Short-Term Variability in the Prompt Emission of Gamma-Ray Bursts

The central engine of gamma-ray bursts (GRBs) is believed to be a hot and dense disk with hyperaccretion onto a few solar-mass black hole. We investigate where the magnetorotational instability (MRI) actively operates in the hyperaccretion disk, which can cause angular momentum transport in the disk. The inner region of hyperaccretion disks can be neutrino opaque, and the energy- and momentum-transport by neutrinos could affect the growth of the MRI significantly. Assuming reasonable disk models and a weak magnetic field $B \lesssim 10^{14} \rm{G}$, it is found that the MRI is strongly suppressed by the neutrino viscosity in the inner region of hyperaccretion disks. On the other hand, the MRI can drive active MHD turbulence in the outer neutrino-transparent region regardless of the field strength. This suggests that the baryonic matter is accumulated into the inner dead zone where the MRI grows inactively and the angular momentum transport is inefficient. When the dead zone gains a large amount of mass and becomes gravitationally unstable, intense mass accretion onto the central black hole would occur episodically through the gravitational torque. This process can be a physical mechanism of the short-term variability in the prompt emission of GRBs. Finally, the origin of flaring activities in the X-ray afterglow is predicted in the context of our episodic accretion scenario.

preprint2007arXiv

Long Term Simulations Of Astrophysical Jets; Energy Structure and Quasi-Periodic Ejection

We have performed self-consistent 2.5-dimensional nonsteady MHD numerical simulations of jet formation as long as possible, including the dynamics of accretion disks. Although the previous nonsteady MHD simulations for astrophysical jets revealed that the characteristics of nonsteady jets are similar to those of steady jets, the calculation time of these simulations is very short compared with the time scale of observed jets. Thus we have investigated long term evolutions of mass accretion rate, mass outflow rate, jet velocity, and various energy flux. We found that the ejection of jet is quasi-periodic. The period of the ejection is related to the time needed for the initial magnetic filed to be twisted to generate toroidal filed. We compare our results with both the steady state theory and previous 2.5-dimensional nonsteady MHD simulations.

preprint2007arXiv

Small Scale Magnetic Flux Emergence Observed with Hinode/Solar Optical Telescope

We observed small scale magnetic flux emergence in a sunspot moat region by the Solar Optical Telescope (SOT) aboard the Hinode satellite. We analyzed filtergram images observed in the wavelengths of Fe 6302 angstrom, G-band and Ca II H. In Stokes I images of Fe 6302 angstrom, emerging magnetic flux were recognized as dark lanes. In G-band, they showed their shapes almost the same as in Stokes I images. These magnetic flux appeared as dark filaments in Ca II H images. Stokes V images of Fe 6302 angstrom showed pairs of opposite polarities at footpoints of each filament. These magnetic concentrations are identified to correspond to bright points in G-band/Ca II H images. From the analysis of time-sliced diagrams, we derived following properties of emerging flux, which are consistent with the previous works. (1) Two footpoints separate each other at a speed of 4.2 km/s during the initial phase of evolution and decreases to about 1 km/s in 10 minutes later. (2) Ca II H filaments appear almost simultaneously with the formation of dark lanes in Stokes I in the observational cadence of 2 minutes. (3) The lifetime of the dark lanes in Stokes I and G-band is 8 minutes, while that of Ca filament is 12 minutes. An interesting phenomena was observed that an emerging flux tube expands laterally in the photosphere with a speed of 3.8 km/s. Discussion on the horizontal expansion of flux tube will be given with refernce to previous simulation studies.

preprint2007arXiv

Umbral Fine Structures in Sunspots Observed with Hinode Solar Optical Telescope

High resolution imaging observation of a sunspot umbra was done with Hinode Solar Optical Telescope (SOT). Filtergrams in wavelengths of blue and green continuum were taken during three consecutive days. The umbra consisted of a dark core region, several diffuse components and numerous umbral dots. We derived basic properties of umbral dots (UDs), especially their temperatures, lifetimes, proper motions, spatial distribution and morphological evolution. Brightness of UDs is confirmed to depend on the brightness of their surrounding background. Several UDs show fission and fusion. Thanks to the stable condition of space observation, we could first follow the temporal behavior of these events. The derived properties of internal structure of the umbra are discussed in viewpoint of magnetoconvection in a strong magnetic field.

preprint2007arXiv

Velocity Structure of Jets in Coronal Hole

Velocity structures of jets in a coronal hole have been derived for the first time. Hinode observations revealed the existence of many bright points in coronal holes. They are loop-shaped and sometimes associated with coronal jets. Spectra obtained with the Extreme ultraviolet Imaging Spectrometer (EIS) on board Hinode are analyzed to infer Doppler velocity of bright loops and jets in a coronal hole of the north polar region. Elongated jets above bright loops are found to be blue-shifted by 30 km/s at maximum, while foot points of bright loops are red-shifted. Blue-shifts detected in coronal jets are interpreted as upflows produced by magnetic reconnection between emerging flux and the ambient field in the coronal hole.

preprint2005arXiv

Magnetic-Reconnection and Current-Sheet Model for the Radio Arc and Threads in the Galactic Center

We propose a new mechanism to explain the radio Arc and threads in the Galactic Center by current sheets produced by local magnetic shears due to interaction of a moving cloud and the vertical field based on three-dimensional magneto-hydrodynamical simulations. Magnetic reconnetion and acceleration of cosmic-ray electrons in the current sheet will result in high contrast of radio emissivity inside and outside the Arc and threads.

preprint2002arXiv

A HR-like Diagram for Solar/Stellar Flares and Corona -- Emission Measure vs Temperature Diagram

In our previous paper, we have presented a theory to explain the observed universal correlation between the emission measure ($EM=n^2 V$) and temperature (T) for solar/stellar flares on the basis of the magnetic reconnection model with heat conduction and chromospheric evaporation. Here n is the electron density and V is the volume. By extending our theory to general situations, we examined the EM-T diagram in detail, and found the following properties: 1) The universal correlation sequence (``main sequence flares'') with $EM \propto T^{17/2}$ corresponds to the case of constant heating flux or equivalently the case of constant magnetic field strength in the reconnection model. 2) The EM-T diagram has a forbidden region, where gas pressure of flares exceeds magnetic pressure. 3) There is a coronal branch with $EM \propto T^{15/2}$ for $T< 10^7$ K, and $EM \propto T^{13/2}$ for $T> 10^7$ K. This branch is situated left side of the main sequence flares in the EM-T diagram. 4) There is another forbidden region determined by the length of flare loop; a lower limit of flare loop is $10^7$ cm. Small flares near this limit correspond to nanoflares observed by SOHO/EIT. 5) We can plot flare evolution track on the EM-T diagram. A flare evolves from the coronal branch to main sequence flares, then returns to the coronal branch eventually. These properties of the EM-T diagram are similar to those of the HR diagram for stars, and thus we propose that the EM-T diagram is quite useful to estimate the physical quantities (loop length, heating flux, magnetic field strength, total energy and so on) of flares and corona when there is no spatially resolved imaging observations.

preprint2001arXiv

Plasmoid-Induced-Reconnection and Fractal Reconnection

As a key to undertanding the basic mechanism for fast reconnection in solar flares, plasmoid-induced-reconnection and fractal reconnection are proposed and examined. We first briefly summarize recent solar observations that give us hints on the role of plasmoid (flux rope) ejections in flare energy release. We then discuss the plasmoid-induced-reconnection model, which is an extention of the classical two-ribbon-flare model which we refer to as the CSHKP model. An essential ingredient of the new model is the formation and ejection of a plasmoid which play an essential role in the storage of magnetic energy (by inhibiting reconnection) and the induction of a strong inflow into reconnection region. Using a simple analytical model, we show that the plasmoid ejection and acceleration are closely coupled with the reconnection process, leading to a nonlinear instability for the whole dynamics that determines the macroscopic reconnection rate uniquely. Next we show that the current sheet tends to have a fractal structure via the following process path: tearing, sheet thinning, Sweet- Parker sheet, secondary tearing, further sheet thinning... These processes occur repeatedly at smaller scales until a microscopic plasma scale (either the ion Larmor radius or the ion inertial length) is reached where anomalous resistivity or collisionless reconnection can occur. The current sheet eventually has a fractal structure with many plasmoids (magnetic islands) of different sizes. When these plasmoids are ejected out of the current sheets, fast reconnection occurs at various different scales in a highly time dependent manner. Finally, a scenario is presented for fast reconnection in the solar corona on the basis of above plasmoid-induced-reconnection in a fractal current sheet.