Researcher profile

Kazunari Shibata

Kazunari Shibata contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 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 &#34;superflares&#34;, 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 &#34;black swans&#34; and extreme &#34;dragon king&#34; 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&#39;&#39; 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&#39;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&#39;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&#39;s coronal temperature, stellar wind velocity, and wind&#39;s mass loss rate. We find that M dwarfs&#39; coronae tend to be cooler than solar corona, and that M dwarfs&#39; 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.