Source author record

Kanya Kusano

Kanya Kusano 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

13works
4topics
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

13 published item(s)

preprint2022arXiv

Super-Droplet Method for the Numerical Simulation of Clouds and Precipitation: a Particle-Based Microphysics Model Coupled with Non-hydrostatic Model

A novel, particle based, probabilistic approach for the simulation of cloud microphysics is proposed, which is named the Super-Droplet Method (SDM). This method enables accurate simulation of cloud microphysics with less demanding cost in computation. SDM is applied to a warm-cloud system, which incorporates sedimentation, condensation/evaporation, and stochastic coalescence. The methodology to couple super-droplets and a non-hydrostatic model is also developed. It is confirmed that the result of our Monte Carlo scheme for the stochastic coalescence of super-droplets agrees fairly well with the solutions of the stochastic coalescence equation. The behavior of the model is evaluated using a simple test problem, that of a shallow maritime cumulus formation initiated by a warm bubble. Possible extensions of SDM are briefly discussed. A theoretical analysis suggests that the computational cost of SDM becomes lower than the spectral (bin) method when the number of attributes - the variables that identify the state of each super-droplet - becomes larger than some critical value, which we estimate to be in the range $2\sim4$.

preprint2022arXiv

The Effects of Oscillations & Collisions of Emerging Bipolar Regions on the Triggering of Solar Flares

The ability to predict the occurrence of solar flares in advance is important to humankind due to the potential damage they can cause to Earth's environment and infrastructure. It has been shown in Kusano et al. (2012) that a small-scale bipolar region (BR), with its flux reversed relative to the potential component of the overlying field, appearing near the polarity inversion line (PIL) is sufficient to effectively trigger a solar flare. In this study we perform further 3D magnetohydrodynamic simulations to study the effect that the motion of these small-scale BRs has on the effectiveness of flare triggering. The effect of two small-scale BRs colliding is also simulated. The results indicate that the strength of the triggered flare is dependent on how much of the overlying field is disrupted by the BR. Simulations of linear oscillations of the BR showed that oscillations along the PIL increase the flare strength whilst oscillations across the PIL detract from the flare strength. The flare strength is affected more by larger amplitude oscillations but is relatively insensitive to the frequency of oscillations. In the most extreme case the peak kinetic energy of the flare increased more than threefold compared to a non-oscillating BR. Simulations of torsional oscillations of the BR showed a very small effect on the flare strength. Finally, simulations of colliding BRs showed the generation of much stronger flares as the flares triggered by each individual BR coalesce. These results show that significantly stronger flares can result from motion of the BR along the PIL of a sheared field or from the presence of multiple BRs in the same region.

preprint2020arXiv

A Comparison of Flare Forecasting Methods. IV. Evaluating Consecutive-Day Forecasting Patterns

A crucial challenge to successful flare prediction is forecasting periods that transition between "flare-quiet" and "flare-active". Building on earlier studies in this series (Barnes et al. 2016; Leka et al. 2019a,b) in which we describe methodology, details, and results of flare forecasting comparison efforts, we focus here on patterns of forecast outcomes (success and failure) over multi-day periods. A novel analysis is developed to evaluate forecasting success in the context of catching the first event of flare-active periods, and conversely, of correctly predicting declining flare activity. We demonstrate these evaluation methods graphically and quantitatively as they provide both quick comparative evaluations and options for detailed analysis. For the testing interval 2016-2017, we determine the relative frequency distribution of two-day dichotomous forecast outcomes for three different event histories (i.e., event/event, no-event/event and event/no-event), and use it to highlight performance differences between forecasting methods. A trend is identified across all forecasting methods that a high/low forecast probability on day-1 remains high/low on day-2 even though flaring activity is transitioning. For M-class and larger flares, we find that explicitly including persistence or prior flare history in computing forecasts helps to improve overall forecast performance. It is also found that using magnetic/modern data leads to improvement in catching the first-event/first-no-event transitions. Finally, 15% of major (i.e., M-class or above) flare days over the testing interval were effectively missed due to a lack of observations from instruments away from the Earth-Sun line.

preprint2020arXiv

Onset Mechanism of M6.5 Solar Flare Observed in Active Region 12371

We studied a flare onset process in terms of stability of a three-dimensional (3D) magnetic field in active region 12371 producing an eruptive M6.5 flare in 2015 June 22. In order to reveal the 3D magnetic structure, we first extrapolated the 3D coronal magnetic fields based on time series of the photospheric vector magnetic fields under a nonlinear force-free field (NLFFF) approximation. The NLFFFs nicely reproduced the observed sigmoidal structure which is widely considered to be preeruptive magnetic configuration. In particular, we found that the sigmoid is composed of two branches of sheared arcade loops. On the basis of the NLFFFs, we investigated the sheared arcade loops to explore the onset process of the eruptive flare using three representative magnetohydrodynamic instabilities: the kink, torus, and double arc instabilities (DAI). The DAI, recently proposed by Ishiguro & Kusano, is a double arc loop that can be more easily destabilized than a torus loop. Consequently, the NLFFFs are found to be quite stable against the kink and torus instabilities. However, the sheared arcade loops formed prior to the flare possibly become unstable against the DAI. As a possible scenario for the onset process of the M6.5 flare, we suggest a three-step process: (1) double arc loops are formed by the sheared arcade loops through the tether-cutting reconnection during an early phase of the flare, (2) the DAI contributes to the expansion of destabilized double arc loops, and (3) finally, the torus instability makes the full eruption.

preprint2019arXiv

A Magnetohydrodynamic Relaxation Method for Non-Force-Free Magnetic Field in Magnetohydrostatic Equilibrium

A nonlinear force-free field (NLFFF) extrapolation is widely used to reconstruct the three-dimensional magnetic field in the solar corona from the observed photospheric magnetic field. However, the pressure gradient and gravitational forces are ignored in the NLFFF model, even though the photospheric and chromospheric magnetic fields are not in general force-free. Here we develop a magnetohydrodynamic (MHD) relaxation method that reconstructs the solar atmospheric (chromospheric and coronal) magnetic field as a non-force-free magnetic field (NFFF) in magnetohydrostatic equilibrium where the Lorentz, pressure gradient, and gravitational forces are balanced. The system of basic equations for the MHD relaxation method is derived, and mathematical properties of the system are investigated. A robust numerical solver for the system is constructed based on the modern high-order shock capturing scheme. Two-dimensional numerical experiments that include the pressure gradient and gravitational forces are also demonstrated.

preprint2014arXiv

Comparison between Hinode/SOT and SDO/HMI, AIA Data for the Study of the Solar Flare Trigger Process

To elucidate the flare trigger mechanism, we have analyzed several flare events which were observed by Hinode/Solar Optical Telescope (SOT), in our previous study. Because of the limitation of SOT field of view, however, only four events in the Hinode data sets have been utilizable. Therefore, increasing the number of events is required for evaluating the flare trigger models. We investigated the applicability of data obtained by the Solar Dynamics Observatory (SDO) to increase the data sample for a statistical analysis of the flare trigger process. SDO regularly observes the full disk of the sun and all flares although its spatial resolution is lower than that of Hinode. We investigated the M6.6 flare which occurred on 13 February 2011 and compared the analyzed data of SDO with the results of our previous study using Hinode/SOT data. Filter and vector magnetograms obtained by the Helioseismic and Magnetic Imager (HMI) and filtergrams from the Atmospheric Imaging Assembly (AIA) 1600A were employed. From the comparison of small-scale magnetic configurations and chromospheric emission prior to the flare onset, we confirmed that the trigger region is detectable with the SDO data. We also measured the magnetic shear angles of the active region and the azimuth and strength of the flare-trigger field. The results were consistent with our previous study. We concluded that statistical studies of the flare trigger process are feasible with SDO as well as Hinode data. We also investigated the temporal evolution of the magnetic field before the flare onset with SDO.

preprint2014arXiv

Coronal Behavior Before the Large Flare Onset

Flares are a major explosive event in our solar system. They are often followed by coronal mass ejection that has a potential to trigger the geomagnetic storms. There are various studies aiming to predict when and where the flares are likely to occur. Most of these studies mainly discuss the photospheric and chromospheric activity before the flare onset. In this paper we study the coronal features before the famous large flare occurrence on December 13th, 2006. Using the data from Hinode/EUV Imaging Spectrometer (EIS), X-Ray Telescope (XRT), and Solar and Heliospheric Observatory (SOHO) /Extreme ultraviolet Imaging Telescope (EIT), we discuss the coronal features in the large scale (~ a few 100 arcsec) before the flare onset. Our findings are as follows: 1) The upflows in and around active region start growing from ~10 to 30 km /s a day before the flare. 2) The expanding coronal loops are clearly observed a few hours before the flare. 3) Soft X-ray and EUV intensity are gradually reduced. 4) The upflows are further enhanced after the flare. From these observed signatures, we conclude that the outer part of active region loops with low density were expanding a day before the flare onset, and the inner part with high density were expanding a few hours before the onset.

preprint2013arXiv

Magnetic Systems Triggering the M6.6-class Solar Flare in NOAA Active Region 11158

We report a detailed event analysis on the M6.6-class flare in the active region (AR) NOAA 11158 on 2011 February 13. AR 11158, which consisted of two major emerging bipoles, showed prominent activities including one X- and several M-class flares. In order to investigate the magnetic structures related to the M6.6 event, particularly the formation process of a flare-triggering magnetic region, we analyzed multiple spacecraft observations and numerical results of a flare simulation. We observed that, in the center of this quadrupolar AR, a highly sheared polarity inversion line (PIL) was formed through proper motions of the major magnetic elements, which built a sheared coronal arcade lying over the PIL. The observations lend support to the interpretation that the target flare was triggered by a localized magnetic region that had an intrusive structure, namely a positive polarity penetrating into a negative counterpart. The geometrical relationship between the sheared coronal arcade and the triggering region was consistent with the theoretical flare model based on the previous numerical study. We found that the formation of the trigger region was due to a continuous accumulation of the small-scale magnetic patches. A few hours before the flare occurrence, the series of emerged/advected patches reconnected with a preexisting fields. Finally, the abrupt flare eruption of the M6.6 event started around 17:30 UT. Our analysis suggests that, in a triggering process of a flare activity, all magnetic systems of multiple scales, not only the entire AR evolution but also the fine magnetic elements, are altogether involved.

preprint2013arXiv

Study of Magnetic Helicity Injection in the Active Regions NOAA 9236 Producing Multiple Flare-associated CME Events

To better understand a preferred magnetic field configuration and its evolution during Coronal Mass Ejection events, we investigated the spatial and temporal evolution of photospheric magnetic fields in the active region NOAA 9236 that produced eight flare-associated CMEs during the time period of 2000 November 23-26. The time variations of the total magnetic helicity injection rate and the total unsigned magnetic flux are determined and examined not only in the entire active region but also in some local regions such as the main sunspots and the CME-associated flaring regions using SOHO/MDI magnetogram data. As a result, we found that: (1) in the sunspots, a large amount of postive (right-handed) magnetic helicity was injected during most of the examined time period, (2) in the flare region, there was a continuous injection of negative (left-handed) magnetic helicity during the entire period, accompanied by a large increase of the unsigned magnetic flux, and (3) the flaring regions were mainly composed of emerging bipoles of magnetic fragments in which magnetic field lines have substantially favorable conditions for making reconnection with large-scale, overlying, and oppositely directed magnetic field lines connecting the main sunspots. These observational findings can also be well explained by some MHD numerical simulations for CME initiation (e.g., reconnection-favored emerging flux models). We therefore conclude that reconnection-favored magnetic fields in the flaring emerging flux regions play a crucial role in producing the multiple flare-associated CMEs in NOAA 9236.

preprint2012arXiv

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

We perform four numerical magnetohydrodynamic simulations in 2.5 dimensions (2.5D) of fast Coronal Mass Ejections (CMEs) and their associated shock fronts between 10Rs and 300Rs. We investigate the relative change in the shock standoff distance, Sd, as a fraction of the CME radial half-width, Dob (i.e. Sd/Dob). Previous hydrodynamic studies have related the shock standoff distance for Earths magnetosphere to the density compression ratio (DR,Ru/Rd) measured across the bow shock (Spreiter, Summers and Alksne 1966). The DR coefficient, kdr, which is the proportionality constant between the relative standoff distance (Sd/Dob) and the compression ratio, was semi-empirically estimated as 1.1. For CMEs, we show that this value varies linearly as a function of heliocentric distance and changes significantly for different radii of curvature of the CMEs leading edge. We find that a value of 0.8+-0.1 is more appropriate for small heliocentric distances (<30Rs) which corresponds to the spherical geometry of a magnetosphere presented by Seiff (1962). As the CME propagates its cross section becomes more oblate and the kdr value increases linearly with heliocentric distance, such that kdr= 1.1 is most appropriate at a heliocentric distance of about 80Rs. For terrestrial distances (215Rs) we estimate kdr= 1.8+-0.3, which also indicates that the CME cross-sectional structure is generally more oblate than that of Earths magnetosphere. These alterations to the proportionality coefficients may serve to improve investigations into the estimates of the magnetic field in the corona upstream of a CME as well as the aspect ratio of CMEs as measured in situ.

preprint2012arXiv

Decorrelation Times of Photospheric Fields and Flows

We use autocorrelation to investigate evolution in flow fields inferred by applying Fourier Local Correlation Tracking (FLCT) to a sequence of high-resolution (0.3 \arcsec), high-cadence ($\simeq 2$ min) line-of-sight magnetograms of NOAA active region (AR) 10930 recorded by the Narrowband Filter Imager (NFI) of the Solar Optical Telescope (SOT) aboard the {\em Hinode} satellite over 12--13 December 2006. To baseline the timescales of flow evolution, we also autocorrelated the magnetograms, at several spatial binnings, to characterize the lifetimes of active region magnetic structures versus spatial scale. Autocorrelation of flow maps can be used to optimize tracking parameters, to understand tracking algorithms' susceptibility to noise, and to estimate flow lifetimes. Tracking parameters varied include: time interval $Δt$ between magnetogram pairs tracked, spatial binning applied to the magnetograms, and windowing parameter $σ$ used in FLCT. Flow structures vary over a range of spatial and temporal scales (including unresolved scales), so tracked flows represent a local average of the flow over a particular range of space and time. We define flow lifetime to be the flow decorrelation time, $τ$. For $Δt > τ$, tracking results represent the average velocity over one or more flow lifetimes. We analyze lifetimes of flow components, divergences, and curls as functions of magnetic field strength and spatial scale. We find a significant trend of increasing lifetimes of flow components, divergences, and curls with field strength, consistent with Lorentz forces partially governing flows in the active photosphere, as well as strong trends of increasing flow lifetime and decreasing magnitudes with increases in both spatial scale and $Δt$.

preprint2011arXiv

Coronal Electron Distribution in Solar Flares: Drift-Kinetic Model

Using a model of particle acceleration and transport in solar flares, we investigate the height distribution of coronal electrons by focusing on the energy-dependent pitch-angle scattering. When pitch-angle scattering is not included, the peak heights of loop-top electrons are constant, regardless of their energy, owing to the continuous acceleration and compression of the electrons via shrinkage of magnetic loops. On the other hand, under pitch-angle scattering, the electron heights are energy dependent; intermediate energy electrons are at a higher altitude, whereas lower and higher energy electrons are at lower altitudes. This implies that the intermediate energy electrons are inhibited to follow the shrinking field lines to lower altitudes because pitch-angle scattering causes efficient precipitation of these electrons into the footpoint and their subsequent loss from the loop. This result is qualitatively consistent with the position of the above-the-loop-top hard X-ray (HXR) source that is located above coronal HXR loops emitted by lower energy electrons and microwaves emitted by higher energy electrons. Quantitative agreement with observations might be achieved by considering primary acceleration before the onset of loop shrinkage and additional pitch-angle scattering via wave-particle interactions.

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.