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.

We report multi-wavelength monitoring observations of an M-dwarf flare star AD Leonis with Seimei Telescope (6150--7930 {\AA}), SCAT (Spectroscopic Chuo-university Astronomical Telescope; 3700--7500 {\AA}), 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$\alpha$, 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$\alpha$ emission line full width at 1/8 maximum dramatically increases to 14 {\AA} from 8 {\AA} in the low-resolution spectra (R$\sim$ 2000) accompanied with the large white-light flares, (2) some weak H$\alpha$/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$\alpha$ 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$\alpha$ 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$\alpha$ 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.

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.

We report results of astrometric VLBI observations toward a Mira variable star BX~Cam using the VLBI array "VERA". The observations were performed from February 2012 to November 2014. Obtained parallax is 1.73$\pm$0.03 mas corresponding to a distance of 0.58$\pm$0.01 kpc. Parallax of this source was also reported in Gaia DR2 as 4.13$\pm$0.25 mas, and there is a 240 \% difference between these two measurements. Astrometric results from our VLBI observations show that we exactly traced angular motions of the seven maser spots in BX~Cam. We calculated stellar luminosities using both parallaxes, and obtained luminosities of $L_{\ast}^{\mathrm{VERA}} = 4950\pm170 L_{\odot}$ and $L_{\ast}^{\mathrm{Gaia}} = 870\pm110 L_{\odot}$. Deduced luminosities also support a validity of the parallax that we determined with VERA. Evaluating the two parallaxes, we concluded that the parallax of 1.73$\pm$0.03 mas from the VERA observations is correct for BX~Cam. We obtained a systemic motion of BX~Cam as (${\mu_{\alpha}}\cos{\delta}^{\mathrm{sys}}$, ${\mu}_{\delta}^{\mathrm{sys}}$) $=$ (13.48$\pm$0.14, $-$34.30$\pm$0.18) mas\,yr$^{-1}$. A total of 73 H$_2$O maser spots detected from our VLBI observations show a spatial distribution of 30 au $\times$ 80 au with a strong elongation along north-south direction. They show outflows with a three-dimensional velocity of 14.79$\pm$1.40 km\,s$^{-1}$. From a comparison between time variations of $V$-band magnitudes and H$_2$O maser, we found that variation of the H$_2$O maser is relevant to that seen in $V$-band even though the H$_2$O maser does not recover its maximum flux in each cycle.

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/.

We present the first astrometry catalog from the Japanese VLBI (very long baseline interferometer) project VERA (VLBI Exploration of Radio Astrometry). We have compiled all the astrometry results from VERA, providing accurate trigonometric annual parallax and proper motion measurements. In total, 99 maser sources are listed in the VERA catalog. Among them, 21 maser sources are newly reported while the rest of 78 sources are referred to previously published results or those in preparation for forthcoming papers. The accuracy in the VERA astrometry are revisited and compared with those from the other VLBI astrometry projects such as BeSSeL (The Bar and Spiral Structure Legacy) Survey and GOBELINS (the Gould's Belt Distances Survey) with the VLBA (Very Long Baseline Array). We have confirmed that most of the astrometry results are consistent with each other, and the largest error sources are due to source structure of the maser features and their rapid variation, along with the systematic calibration errors and different analysis methods. Combined with the BeSSeL results, we estimate the up-to-date fundamental Galactic parameter of $R_{0}=7.92\pm0.16_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~kpc and $\Omega_{\odot}=30.17\pm0.27_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~km~s$^{-1}$~kpc$^{-1}$, where $R_{0}$ and $\Omega_{\odot}$ are the distance from the Sun to the Galactic center and the Sun's angular velocity of the Galactic circular rotation, respectively.

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.

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\AA\ 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.

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.

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($\Gamma\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

We report measurements of parallax and proper motion for four 22 GHz water maser sources as part of VERA Outer Rotation Curve project. All sources show Galactic latitudes of $>$ 2$^{\circ}$ and Galactocentric distances of $>$ 11 kpc at the Galactic longitude range of 95$^{\circ}$ $< l <$ 126$^{\circ}$. The sources trace the Galactic warp reaching to 200$\sim$400 pc, and indicate the signature of the warp to 600 pc toward the north Galactic pole. The new results along with previous results in the literature show the maximum height of the Galactic warp is increased with Galactocentric distance. Also, we examined velocities perpendicular to the disk for the sample, and found an oscillatory behavior between the vertical velocities and Galactic heights. This behavior suggests the existence of the bending (vertical density) waves, possibly induced by a perturbing satellite (e.g. passage of the Sagittarius dwarf galaxy).

As part our investigation into the Galactic rotation curve, we carried out Very Long Baseline Interferometry (VLBI) observations towards the star-forming region IRAS 01123+6430 using VLBI Exploration of Radio Astrometry (VERA) to measure its annual parallax and proper motion. The annual parallax was measured to be 0.151+/-0.042 mas, which corresponds to a distance of D=6.61^{+2.55}_{-1.44} kpc, and the obtained proper motion components were ({\mu}_{\alpha}cos{\delta}, {\mu}_{\delta})=(-1.44+/-0.15, -0.27+/-0.16) mas yr^{-1} in equatorial coordinates. Assuming Galactic constants of (R_0, {\Theta}_0)=(8.05+/-0.45 kpc, 238+/-14 km s^{-1}), the Galactocentric distance and rotation velocity were measured to be (R, {\Theta})=(13.04+/-2.24 kpc, 239+/-22 km s^{-1}), which are consistent with a flat Galactic rotation curve. The newly estimated distance provides a more accurate bolometric luminosity of the central young stellar object, L_Bol=(3.11+/-2.86)\times 10^3 L_solar, which corresponds to a spectral type of B1-B2. The analysis of 12CO(J=1-0) survey data obtained with the Five College Radio Astronomical Observatory (FCRAO) 14 m telescope shows that the molecular cloud associated with IRAS 01123+6430 consists of arc-like and linear components, which well matches a structure predicted by numerical simulation of the cloud-cloud collision (CCC) phenomenon. The coexistence of arc-like and linear components implies that the relative velocity of initial two clouds was as slow as 3-5 km s^{-1}, which meets the expected criteria of massive star formation where the core mass is effectively increased in the presence of low relative velocity (~3-5 km s^{-1}), as suggested by Takahira et al.(2014).

On 2014 March 29, an intense solar flare classified as X1.0 occurred in the active region 12017. Several associated phenomena accompanied this event, among them a fast-filament eruption, large-scale propagating disturbances in the corona and the chromosphere including a Moreton wave, and a coronal mass ejection. This flare was successfully detected in multiwavelength imaging in H-alpha line by the Flare Monitoring Telescope (FMT) at Ica University, Peru. We present a detailed study of the Moreton wave associated with the flare in question. Special attention is paid to the Doppler characteristics inferred from the FMT wing (H-alpha$\pm0.8$~{\AA}) observations, which are used to examine the downward/upward motion of the plasma in the chromosphere. Our findings reveal that the downward motion of the chromospheric material at the front of the Moreton wave attains a maximum velocity of 4 km/s, whereas the propagation speed ranges between 640 and 859 km/s. Furthermore, utilizing the weak shock approximation in conjunction with the velocity amplitude of the chromospheric motion induced by the Moreton wave, we derive the Mach number of the incident shock in the corona. We also performed the temperature-emission measure analysis of the coronal wave based on the Atmospheric Imaging Assembly (AIA) observations, which allowed us to derive the compression ratio, and to estimate the Alfv\'en and fast-mode Mach numbers of the order of 1.06-1.28 and 1.05-1.27. Considering these results and the MHD linear theory we discuss the characteristics of the shock front and the interaction with the chromospheric plasma.

High-energy radiation caused by exoplanetary space weather events from planet-hosting stars can play a crucial role in conditions promoting or destroying habitability in addition to the conventional factors. In this paper, we present the first quantitative impact evaluation system of stellar flares on the habitability factors with an emphasis on the impact of Stellar Proton Events. We derive the maximum flare energy from stellar starspot sizes and examine the impacts of flare associated ionizing radiation on CO$_2$, H$_2$, N$_2$+O$_2$ --rich atmospheres of a number of well-characterized terrestrial type exoplanets. Our simulations based on the Particle and Heavy Ion Transport code System [PHITS] suggest that the estimated ground level dose for each planet in the case of terrestrial-level atmospheric pressure (1 bar) for each exoplanet does not exceed the critical dose for complex (multi-cellular) life to persist, even for the planetary surface of Proxima Centauri b, Ross-128 b and TRAPPIST-1 e. However, when we take into account the effects of the possible maximum flares from those host stars, the estimated dose reaches fatal levels at the terrestrial lowest atmospheric depth on TRAPPIST-1 e and Ross-128 b. Large fluxes of coronal XUV radiation from active stars induces high atmospheric escape rates from close-in exoplanets suggesting that the atmospheric depth can be substantially smaller than that on the Earth. In a scenario with the atmospheric thickness of 1/10 of Earth's, the radiation dose from close-in planets including Proxima Centauri b and TRAPPIST-1 e reach near fatal dose levels with annual frequency of flare occurrence from their hoststars.

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

We report the latest view of Kepler solar-type (G-type main-sequence) superflare stars, including recent updates with Apache Point Observatory (APO) 3.5m telescope spectroscopic observations and Gaia-DR2 data. First, we newly conducted APO3.5m spectroscopic observations of 18 superflare stars found from Kepler 1-min time cadence data. More than half (43 stars) are confirmed to be "single" stars, among 64 superflare stars in total that have been spectroscopically investigated so far in this APO3.5m and our previous Subaru/HDS observations. The measurements of $v\sin i$ (projected rotational velocity) and chromospheric lines (Ca II H\&K and Ca II 8542\AA) support the brightness variation of superflare stars is caused by the rotation of a star with large starspots. We then investigated the statistical properties of Kepler solar-type superflare stars by incorporating Gaia-DR2 stellar radius estimates. As a result, the maximum superflare energy continuously decreases as the rotation period $P_{\mathrm{rot}}$ increases. Superflares with energies $\lesssim 5\times10^{34}$ erg occur on old, slowly-rotating Sun-like stars ($P_{\mathrm{rot}}\sim$25 days) approximately once every 2000--3000 years, while young rapidly-rotating stars with $P_{\mathrm{rot}}\sim$ a few days have superflares up to $10^{36}$ erg. The maximum starspot area does not depend on the rotation period when the star is young, but as the rotation slows down, it starts to steeply decrease at $P_{\mathrm{rot}}\gtrsim$12 days for Sun-like stars. These two decreasing trends are consistent since the magnetic energy stored around starspots explains the flare energy, but other factors like spot magnetic structure should also be considered.

The KVN(Korean VLBI Network)-style simultaneous multi-frequency receiving mode is demonstrated to be promising for mm-VLBI observations. Recently, other Very long baseline interferometry (VLBI) facilities all over the globe start to implement compatible optics systems. Simultaneous dual/multi-frequency VLBI observations at mm wavelengths with international baselines are thus possible. In this paper, we present the results from the first successful simultaneous 22/43 GHz dual-frequency observation with KaVA(KVN and VERA array), including images and astrometric results. Our analysis shows that the newly implemented simultaneous receiving system has brought a significant extension of the coherence time of the 43 GHz visibility phases along the international baselines. The astrometric results obtained with KaVA are consistent with those obtained with the independent analysis of the KVN data. Our results thus confirm the good performance of the simultaneous receiving systems for the non-KVN stations. Future simultaneous observations with more global stations bring even higher sensitivity and micro-arcsecond level astrometric measurements of the targets.

Magnetic reconnection may be the fundamental process allowing energy stored in magnetic fields to be released abruptly, solar flares and coronal mass ejection (CME) being archetypal natural plasma examples. Magnetic reconnection is much too slow a process to be efficient on the large scales, but accelerates once small enough scales are formed in the system. For this reason, the fractal reconnection scenario was introduced (Shibata and Tanuma 2001) to explain explosive events in the solar atmosphere: it was based on the recursive triggering and collapse via tearing instability of a current sheet originally thinned during the rise of a filament in the solar corona. Here we compare the different fractal reconnection scenarios that have been proposed, and derive generalized scaling relations for the recursive triggering of fast, `ideal' - i.e. Lundquist number independent - tearing in collapsing current sheet configurations with arbitrary current profile shapes. An important result is that the Sweet-Parker scaling with Lundquist number, if interpreted as the aspect ratio of the singular layer in an ideally unstable sheet, is universal and does not depend on the details of the current profile in the sheet. Such a scaling however must not be interpreted in terms of stationary reconnection, rather it defines a step in the accelerating sequence of events of the ideal tearing mediated fractal cascade. We calculate scalings for the expected number of plasmoids for such generic profiles and realistic Lundquist numbers.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

Filaments, the dense cooler plasma floating in the solar corona supported by magnetic fields, generally exhibit certain activations before they erupt. In our previous study (Seki et al. 2017 ), we observed that the standard deviation of the line-of-sight (LOS) velocities of the small-scale motions in a filament increased prior to its eruption. However, because that study only analyzed one event, it is unclear whether such an increase in the standard deviation of LOS velocities is common in filament eruptions. In this study, 12 filaments that vanished in H{\alpha} line center images were analyzed in a manner similar to the one in our previous work; these included two quiescent filaments, four active region filaments, and six intermediate filaments. We verified that in all the 12 events, the standard deviation of the LOS velocities increased before the filaments vanished. Moreover, we observed that the quiescent filaments had approximately 10 times longer duration of an increase in the standard deviation than the other types of filaments. We concluded that the standard deviation of the LOS velocities of the small-scale motions in a filament can potentially be used as the precursor of a filament eruption.

We report the first unambiguous quasi-periodic large-scale extreme-ultraviolet (EUV) wave or shock that was detected by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. During the whip-like unwinding eruption of a small filament on 2012 April 24, multiple consecutive large-scale wavefronts emanating from AR11467 were observed simultaneously along the solar surface and a closed transequatorial loop system. In the meantime, an upward propagating dome-shaped wavefront was also observed, whose initial speed and deceleration are about 1392 km/s and 1.78 km/s^2, respectively. Along the solar surface, the quasi-peridoic wavefronts had a period of about 163 +/- 21 seconds and propagated at a nearly constant speed of 747 +/- 26 km/s; they interacted with active region AR11469 and launched a sympathetic upward propagating secondary EUV wave. The wavefronts along the loop system propagated at a speed of 897 km/s, and they were reflected back at the southern end of the loop system at a similar speed. In addition to the propagating waves, a standing kink wave was also present in the loop system simultaneously. Periodicity analysis reveals that the period of the wavefronts was consistent with that of the unwinding helical structures of the erupting filament. Based on these observational facts, we propose that the observed quasi-periodic EUV wavefronts were most likely excited by the periodic unwinding motion of the filament helical structures. In addition, two different seismological methods are applied to derive the magnetic field strength of the loop system, and for the first time the reliability of these inversion techniques are tested with the same magnetic structure.

Active solar-type stars show large quasi-periodic brightness variations caused by stellar rotations with star spots, and the amplitude changes as the spots emerge and decay. The Kepler data are suitable for investigations on the emergence and decay processes of star spots, which are important to understand underlying stellar dynamo and stellar flares. In this study, we measured temporal evolutions of star spot area with Kepler data by tracing local minima of the light curves. In this analysis, we extracted temporal evolutions of star spots showing clear emergence and decay without being disturbed by stellar differential rotations. We applied this method to 5356 active solar-type stars observed by Kepler and obtained temporal evolutions of 56 individual star spots. We calculated lifetimes, emergence and decay rates of the star spots from the obtained temporal evolutions of spot area. As a result, we found that lifetimes ($T$) of star spots are ranging from 10 to 350 days when spot areas ($A$) are 0.1-2.3 percent of the solar hemisphere. We also compared them with sunspot lifetimes, and found that the lifetimes of star spots are much shorter than those extrapolated from an empirical relation of sunspots ($T\propto A$), while being consistent with other researches on star spot lifetimes. The emerging and decay rates of star spots are typically $5 \times 10^{20}$ $\rm Mx\cdot h^{-1}$ ($8$ $\rm MSH\cdot h^{-1}$) with the area of 0.1-2.3 percent of the solar hemisphere and are mostly consistent with those expected from sunspots, which may indicate the same underlying processes.

In this article, we insist on the importance and the challenges of the prediction of solar eruptive phenomena including flares, coronal mass ejections (CME), and filament eruptions fully based on the ground-based telescopes. It is true that satellites' data are indispensable for the space weather prediction, but they are vulnerable to the space weather effects. Therefore, the ground-based telescopes can be complementary to them from the viewpoint of space weather prediction. From this view point, one possible new flare prediction method that makes use of H-alpha, red wings, and blue wings images obtained by the SDDI/SMART, the ground-based telescope at Hida Observatory, is presented. And in order to show the possibility for the actual operation based on that method, the recent progress of CHAIN project, the international observation network, is mentioned in terms of their outcomes and capacity buildings.

The study of historical great geomagnetic storms is crucial for assessing the possible risks to the technological infrastructure of a modern society, caused by extreme space-weather events. The normal benchmark has been the great geomagnetic storm of September 1859, the so-called "Carrington Event". However, there are numerous records of another great geomagnetic storm in February 1872. This storm, about 12 years after the Carrington Event, resulted in comparable magnetic disturbances and auroral displays over large areas of the Earth. We have revisited this great geomagnetic storm in terms of the auroral and sunspot records in the historical documents from East Asia. In particular, we have surveyed the auroral records from East Asia and estimated the equatorward boundary of the auroral oval to be near 24.3 deg invariant latitude (ILAT), on the basis that the aurora was seen near the zenith at Shanghai (20 deg magnetic latitude, MLAT). These results confirm that this geomagnetic storm of February 1872 was as extreme as the Carrington Event, at least in terms of the equatorward motion of the auroral oval. Indeed, our results support the interpretation of the simultaneous auroral observations made at Bombay (10 deg MLAT). The East Asian auroral records have indicated extreme brightness, suggesting unusual precipitation of high-intensity, low-energy electrons during this geomagnetic storm. We have compared the duration of the East Asian auroral displays with magnetic observations in Bombay and found that the auroral displays occurred in the initial phase, main phase, and early recovery phase of the magnetic storm.

We have performed 5 night spectroscopic observation of the Halpha line of EV Lac with a medium wavelength resolution (R~ 10,000) using the 2m Nayuta telescope at the Nishi-Harima Astronomical Observatory. EV Lac always possesses the Halpha emission line; however, its intensity was stronger on August 15, 2015 than during other four-night periods. On this night, we observed a rapid rise (~ 20min) and a subsequent slow decrease (~ 1.5h) of the emission-line intensity of Halpha, which was probably caused by a flare. We also found an asymmetrical change in the Halpha line on the same night. The enhancement has been observed in the blue wing of the Ha line during each phase of this flare (from the flare start to the flare end), and absorption components were present in its red wing during the early and later phases of the flare. Such blue enhancement (blue asymmetry) of the Halpha line is sometimes seen during solar flares, but only during the early phases. Even for solar flares, little is known about the origin of the blue asymmetry. Compared with solar-flare models, the presented results can lead to the understanding of the dynamics of stellar flares.

Since the discovery of stellar superflares by Kepler satellite, these extremely energetic events have been studied in analogy to solar flares. Their white-light (WL) continuum emission has been interpreted as being produced by heated ribbons. In this paper we compute the WL emission from overlying flare loops depending on their density and temperature and show that, under conditions expected during superflares, the continuum brightening due to extended loop arcades can significantly contribute to stellar flux detected by Kepler. This requires electron densities in the loops $10^{12} - 10^{13}$ cm$^{-3}$ or higher. We show that such densities, exceeding those typically present in solar flare loops, can be reached on M-dwarf and solar-type superflare stars with large starspots and much stronger magnetic fields. Quite importantly, the WL radiation of loops is not very sensitive to their temperature and thus both cool as well as hot loops may contribute. We show that the WL intensity emergent from optically-thin loops is lower than the blackbody radiation from flare ribbons, but the contribution of loops to total stellar flux can be quite important due to their significant emitting areas. This new scenario for interpreting superflare emission suggests that the observed WL flux is due to a mixture of the ribbon and loop radiation and can be even loop-dominated during the gradual phase of superflares.

Sunspot observations with telescopes in 18th century were carried out in Japan as well. One of these sunspot observations is recorded in an account called Sansaizusetsu narabini Kansei irai Jissoku Zusetsu (Charts of Three Worlds and Diagrams of Actual Observations since Kansei Era). We analyze manuscripts of this account to show in total 15 sunspot drawings in 1749-1750. These observations were carried out by contemporary official astronomers in Japan, with telescopes covered by zongurasus (< zonglas in Dutch, corresponding to "sunglass" in English). We count their group number of sunspots to locate them in long-term solar activity and show that their observations were situated around the solar maximum in 1749 or 1750. We also computed their locations and areas, while we have to admit the difference of variant manuscripts with one another. These observational records show the spread of sunspot observations not only in Europe but also in Japan and hence may contribute to crosscheck or possibly improve the known sunspot indices.

Recently, many superflares on solar-type stars were discovered as white-light flares (WLFs). A correlation between the energies (E) and durations (t) of superflares is derived as $t\propto E^{0.39}$, and this can be theoretically explained by magnetic reconnection ($t\propto E^{1/3}$). In this study, we carried out a statistical research on 50 solar WLFs with SDO/HMI to examine the t-E relation. As a result, the t-E relation on solar WLFs ($t\propto E^{0.38}$) is quite similar stellar superflares, but the durations of stellar superflares are much shorter than those extrapolated from solar WLFs. We present the following two interpretations; (1) in solar flares, the cooling timescale of WL emission may be longer than the reconnection one, and the decay time can be determined by the cooling timescale; (2) the distribution can be understood by applying a scaling law $t\propto E^{1/3}B^{-5/3}$ derived from the magnetic reconnection theory.

We performed coordinated observations of AR 12205, which produced a C-class flare on 2014 November 11, with the Interface Region Imaging Spectrograph (IRIS) and the Domeless Solar Telescope (DST) at Hida Observatory. Using spectral data in the Si IV 1403 \AA, C II 1335 \AA, and Mg II h and k lines from IRIS and the Ca II K, Ca II 8542 \AA, and H$\alpha$ lines from DST, we investigated a moving flare kernel during the flare. In the Mg II h line, the leading edge of the flare kernel showed the intensity enhancement in the blue wing, and the smaller intensity of the blue-side peak (h2v) than that of the red-side one (h2r). The blueshift lasted for 9-48 s with a typical speed of 10.1 $\pm$ 2.6 km s$^{-1}$ and it was followed by the high intensity and the large redshift with a speed of up to 51 km s$^{-1}$ detected in the Mg II h line. The large redshift was a common property for all six lines but the blueshift prior to it was found only in the Mg II lines. A cloud modeling of the Mg II h line suggests that the blue wing enhancement with such peak difference can be caused by a chromospheric-temperature (cool) upflow. We discuss a scenario in which an upflow of cool plasma is lifted up by expanding hot plasma owing to the deep penetration of non-thermal electrons into the chromosphere. Furthermore, we found that the blueshift persisted without any subsequent redshift in the leading edge of the flare kernel during its decaying phase. The cause of such long-lasting blueshift is also discussed.

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

Three Japanese sunspot drawings associated with Iwahashi Zenbei (1756-1811) are shown here from contemporary manuscripts and woodprint documents with the relevant texts. We revealed the observational date of one of the drawings to be 26 August 1793, and the overall observations lasted for over a year. Moreover, we identified the observational site for the dated drawing at Fushimi in Japan. We then compared his observations with group sunspot number and raw group count from Sunspot Index and Long-term Solar Observations (SILSO) to reveal its data context, and concluded that these drawings filled the gaps in understanding due to the fragmental sunspots observations around 1793. These drawings are important as a clue to evaluate astronomical knowledge of contemporary Japan in the late 19 th century and are valuable as a non-European observation, considering that most sunspot observations up to mid-19 th century are from Europe.

Dim red aurora at low magnetic latitudes is a visual and recognized manifestation of geomagnetic storms. The great low-latitude auroral displays seen throughout East Asia on 16-18 September 1770 are considered to manifest one of the greatest storms. Recently found 111 historical documents in East Asia attest that these low-latitude auroral displays were succeeding for almost 9 nights during 10-19 September 1770 in the lowest magnetic latitude areas (< 30{\deg}). This suggests that the duration of the great magnetic storm is much longer than usual. Sunspot drawings from 1770 reveals the fact that sunspots area was twice as large as those observed in another great storm of 1859, which substantiates this unusual storm activities in 1770. These spots likely ejected several huge, sequential magnetic structures in short duration into interplanetary space, resulting in spectacular world-wide aurorae in mid-September 1770. These findings provide new insights about the history, duration, and effects of extreme magnetic storms that may be valuable for those who need to mitigate against extreme events.

Recently, many superflares on solar-type stars have been discovered as white-light flares (WLFs). The statistical study found a correlation between their energies ($E$) and durations ($\tau$): $\tau \propto E^{0.39}$ (Maehara et al. 2017 $EP\& S$, 67, 59), similar to those of solar hard/soft X-ray flares: $\tau \propto E^{0.2-0.33}$. This indicates a universal mechanism of energy release on solar and stellar flares, i.e., magnetic reconnection. We here carried out a statistical research on 50 solar WLFs observed with \textit{SDO}/HMI and examined the correlation between the energies and durations. As a result, the $E$--$\tau$ relation on solar WLFs ($\tau \propto E^{0.38}$) is quite similar to that on stellar superflares ($\tau \propto E^{0.39}$). However, the durations of stellar superflares are one order of magnitude shorter than those expected from solar WLFs. We present the following two interpretations for the discrepancy. (1) In solar flares, the cooling timescale of WLFs may be longer than the reconnection one, and the decay time of solar WLFs can be elongated by the cooling effect. (2) The distribution can be understood by applying a scaling law ($\tau \propto E^{1/3}B^{-5/3}$) derived from the magnetic reconnection theory. In this case, the observed superflares are expected to have 2-4 times stronger magnetic field strength than solar flares.

We present the results of comparative study of amplitude calibrations for East-Asia VLBI Network (EAVN) at 22 and 43 GHz using two different methods of an "a-priori" and a "template spectrum", particularly on lower declination sources. Using observational data sets of early EAVN observations, we investigated the elevation-dependence of the gain values at seven stations of the KaVA (KVN and VERA Array) and three additional telescopes in Japan (Takahagi 32m, Yamaguchi 32m and Nobeyama 45m). By comparing the independently obtained gain values based on these two methods, we found that the gain values from each method were consistent within 10% at elevations higher than 10 degree. We also found that the total flux densities of two images produced from the different amplitude calibrations were in agreement within 10% at both 22 and 43 GHz. By using the template spectrum method, furthermore, the additional radio telescopes can participate in the KaVA (i.e. EAVN) so that it can give a notable sensitivity increase. Therefore, our results will constrain the detailed conditions to reliably measure the VLBI amplitude using EAVN and give a potential to extend possible telescopes comprising EAVN.

We report a solar jet phenomenon associated with the C5.4 class flare on 2014 November 11. The data of jet was provided by Solar Dynamics Observatory (SDO), X-Ray Telescope (XRT) aboard Hinode, Interface Region Imaging Spectrograph (IRIS) and Domeless Solar Telescope (DST) at Hida Observatory, Kyoto University. These plentiful data enabled us to present this series of papers to discuss the entire processes of the observed phenomena including the energy storage, event trigger, and energy release. In this paper, we focus on the energy release process of the observed jet, and mainly describe our spectral analysis on the H-alpha data of DST to investigate the internal structure of the H-alpha jet and its temporal evolution. This analysis reveals that in the physical quantity distributions of the H-alpha jet, such as line-of-sight velocity and optical thickness, there is a significant gradient in the direction crossing the jet. We interpret this internal structure as the consequence of the migration of energy release site, based on the idea of ubiquitous reconnection. Moreover, by measuring the horizontal flow of the fine structures in jet, we succeeded in deriving the three-dimensional velocity field and the line-of-sight acceleration field of the H-alpha jet. The analysis result indicates a part of ejecta in the H-alpha jet experienced the additional acceleration after it had been ejected from the lower atmosphere. This secondary acceleration was specified to occur in the vicinity of the intersection between the trajectories of the H-alpha jet and X-ray jet observed by Hinode/XRT. We propose that a fundamental cause of this phenomenon is the magnetic reconnection involving the plasmoid in the observed jet.

We propose a mechanism for quasi-periodic oscillations of both coronal mass ejections (CMEs) and flare loops as related to magnetic reconnection in eruptive solar flares. We perform two-dimensional numerical MHD simulations of magnetic flux rope eruption, with three different values of the global Lundquist number. In the low Lundquist number run, no oscillatory behavior is found. In the moderate Lunquist number run, on the other hand, quasi-periodic oscillations are excited both at the bottom of the flux rope and at the flare loop-top. In the high Lundquist number run, quasi-periodic oscillations are also excited; in the meanwhile, the dynamics become turbulent due to the formation of multiple plasmoids in the reconnection current sheet. In high and moderate Lundquist number runs, thin reconnection jet collide with the flux rope bottom or flare loop-top and dig them deeply. Steep oblique shocks are formed as termination shocks where reconnection jet is bent (rather than decelerated) in horizontal direction, resulting in supersonic back-flows. The structure becomes unstable, and quasi-periodic oscillation of supersonic back-flows appear at locally confined high-beta region at both the flux rope bottom and flare loop-top. We compare the observational characteristics of quasi-periodic oscillations in erupting flux ropes, post-CME current sheets, flare ribbons and light curves, with corresponding dynamical structures found in our simulation.

There is a subclass of the X-ray jets from young stellar objects which are heated very close to the footpoint of the jets, particularly DG Tau jets. Previous models attribute the strong heating to shocks in the jets. However, the mechanism that localizes the heating at the footpoint remains puzzling. We presented a different model of such X-ray jets, in which the disk atmosphere is magnetically heated. Our disk corona model is based on the so-called nanoflare model for the solar corona. We show that the magnetic heating near the disks can result in the formation of a hot corona with a temperature of > 10^6 K even if the average field strength in the disk is moderately weak, > 1 G. We determine the density and the temperature at the jet base by considering the energy balance between the heating and cooling. We derive the scaling relations of the mass loss rate and terminal velocity of jets. Our model is applied to the DG Tau jets. The observed temperature and estimated mass loss rate are consistent with the prediction of our model in the case of the disk magnetic field strength of ~20 G and the heating region of < 0.1 au. The derived scaling relation of the temperature of X-ray jets could be a useful tool to estimate the magnetic field strength. We also found that the jet X-ray can have a significant impact on the ionization degree near the disk surface and the dead-zone size.

In this article, we present the results of the surveys on sunspots and auroral candidates in Rikkokushi, Japanese Official Histories from the early 7th century to 887 to review the solar and auroral activities. In total, we found one sunspot record and 13 auroral candidates in Rikkokushi. We then examine the records of the sunspots and auroral candidates, compare the auroral candidates with the lunar phase to estimate the reliability of the auroral candidates, and compare the records of the sunspots and auroral candidates with the contemporary total solar irradiance reconstructed from radioisotope data. We also identify the locations of the observational sites to review possible equatorward expansion of auroral oval. These discussions suggest a major gap of auroral candidates from the late 7th to early 9th century, which includes the minimum number of candidates reconstructed from the radioisotope data, a similar tendency as the distributions of sunspot records in contemporary China, and a relatively high magnetic latitude with a higher potential for observing aurorae more frequently than at present.

We report a solar jet phenomenon associated with successive flares on November 10th 2014. These explosive events were involved with the satellite spots' emergence around a delta-type sunspot in the decaying active region NOAA 12205. The data of this jet was provided by Solar Dynamic Observatory (SDO), X-Ray Telescope (XRT) aboard Hinode, Interface Region Imaging Spectrograph (IRIS) and Domeless Solar Telescope (DST) at Hida Observatory, Kyoto University. These plentiful data enabled us to present this series of papers to discuss the entire processes of the observed phenomena including the energy storage, event trigger, and energy release. In this paper, we focus on the energy build-up and trigger phases, by analyzing the photospheric horizontal flow field around the active region with an optical flow method. The analysis reveals the following three. (i) The observed explosive phenomena involved three satellite spots, the magnetic fluxes of which successively reconnected with their pre-existing ambient fields. (ii) All of these satellite spots emerged in the moat region of a pivotal delta-type sunspot, especially near its convergent boundary with the neighboring supergranules or moat regions of adjacent sunspots. (iii) Around the jet ejection site, the positive polarities of satellite spot and adjacent emerging flux encountered the global magnetic field with negative polarity in the moat region of the pivotal delta-type sunspot, and thus the polarity inversion line was formed along the convergent boundary of the photospheric horizontal flow channels.

We report the initial results of our high-cadence monitoring program on the radio jet in the active galaxy M87, obtained by the KVN and VERA Array (KaVA) at 22 GHz. This is a pilot study that preceded a larger KaVA-M87 monitoring program, which is currently ongoing. The pilot monitoring was mostly performed every two to three weeks from December 2013 to June 2014, at a recording rate of 1 Gbps, obtaining the data for a total of 10 epochs. We successfully obtained a sequence of good quality radio maps that revealed the rich structure of this jet from <~1 mas to 20 mas, corresponding to physical scales (projected) of ~0.1-2 pc (or ~140-2800 Schwarzschild radii). We detected superluminal motions at these scales, together with a trend of gradual acceleration. The first evidence for such fast motions and acceleration near the jet base were obtained from recent VLBA studies at 43 GHz, and the fact that very similar kinematics are seen at a different frequency and time with a different instrument suggests these properties are fundamental characteristics of this jet. This pilot program demonstrates that KaVA is a powerful VLBI array for studying the detailed structural evolution of the M87 jet and also other relativistic jets.

We analyze the correlation between starspots and superflares on solar-type stars using observations from the Kepler mission. The analysis shows that the observed fraction of stars with superflares decreases as the rotation period increases and as the amplitude of photometric variability associated with rotation decreases. We found that the fraction of stars with superflares among the stars showing large-amplitude rotational variations, which are thought to be the signature of the large starspots, also decreases as the rotation period increases. The small fraction of superflare stars among the stars with large starspots in the longer-period regime suggests that some of the stars with large starspots show a much lower flare activity than the superflare stars with the same spot area. Assuming simple relations between spot area and lifetime and between spot temperature and photospheric temperature, we compared the size distribution of large starspot groups on slowly-rotating solar-type stars with that of sunspot groups. The size distribution of starspots shows the power-law distribution and the size distribution of larger sunspots lies on this power-law line. We also found that frequency-energy distributions for flares originating from spots with different sizes are the same for solar-type stars with superflares and the Sun. These results suggest that the magnetic activity we observe on solar-type stars with superflares and that on the Sun is caused by the same physical processes.

Fast interplanetary coronal mass ejections (interplanetary CMEs, or ICMEs) are the drivers of strongest space weather storms such as solar energetic particle events and geomagnetic storms. The connection between space weather impacting solar wind disturbances associated with fast ICMEs at Earth and the characteristics of causative energetic CMEs observed near the Sun is a key question in the study of space weather storms as well as in the development of practical space weather prediction. Such shock-driving fast ICMEs usually expand at supersonic speed during the propagation, resulting in the continuous accumulation of shocked sheath plasma ahead. In this paper, we propose the "sheath-accumulating propagation" (SAP) model that describe the coevolution of the interplanetary sheath and decelerating ICME ejecta by taking into account the process of upstream solar wind plasma accumulation within the sheath region. Based on the SAP model, we discussed (1) ICME deceleration characteristics, (2) the fundamental condition for fast ICME at Earth, (3) thickness of interplanetary sheath, (4) arrival time prediction and (5) the super-intense geomagnetic storms associated with huge solar flares. We quantitatively show that not only speed but also mass of the CME are crucial in discussing the above five points. The similarities and differences among the SAP model, the drag-based model and the`snow-plough' model proposed by \citet{tappin2006} are also discussed.

Coronal disturbances associated with solar flares, such as H$\alpha$ Moreton waves, X-ray waves, and extreme ultraviolet (EUV) coronal waves are discussed herein in relation to magnetohydrodynamics fast-mode waves or shocks in the corona. To understand the mechanism of coronal disturbances, full-disk solar observations with high spatial and temporal resolution over multiple wavelengths are of crucial importance. We observed a filament eruption, whose shape is like a "dandelion", associated with the M1.6 flare that occurred on 2011 February 16 in the H$\alpha$ images taken by the Flare Monitoring Telescope at Ica University, Peru. We derive the three-dimensional velocity field of the erupting filament. We also identify winking filaments that are located far from the flare site in the H$\alpha$ images, whereas no Moreton wave is observed. By comparing the temporal evolution of the winking filaments with those of the coronal wave seen in the extreme ultraviolet images data taken by the Atmospheric Imaging Assembly on board the {\it Solar Dynamics Observatory} and by the Extreme Ultraviolet Imager on board the {\it Solar Terrestrial Relations Observatory-Ahead}, we confirm that the winking filaments were activated by the EUV coronal wave.

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\AA 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}^\alpha$, $F_p\propto F_{SXR}^\beta$ and $F_p\propto V_{CME}^\gamma$ with $\alpha = 0.30\pm 0.04$, $\beta = 1.19 \pm 0.08$ and $\gamma = 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.

Recent studies of solar-type superflare stars have suggested that even old slowly rotating stars similar to the Sun can have large starspots and superflares. We conducted high dispersion spectroscopy of 49 nearby solar-analog stars (G-type main sequence stars with $T_{\rm{eff}}\approx5,600\sim6,000$ K) identified as ROSAT soft X-ray sources, which are not binary stars from the previous studies. We expected that these stars can be used as a proxy of bright solar-analog superflare stars, since superflare stars are expected to show strong X-ray luminosity. More than half (37) of the 49 target stars show no evidence of binarity, and their atmospheric parameters ($T_{\rm{eff}}$, $\log g$, and [Fe/H]) are within the range of ordinary solar-analog stars. We measured Ca II 8542 and H$\alpha$ lines, which are good indicators of the chromospheric activity. The intensity of these lines indicates that all the target stars have large starspots. We also measured $v\sin i$ (projected rotational velocity) and Lithium abundance for the target stars. Li abundance is a key to understanding the evolution of the stellar convection zone, which reflects the stellar age, mass and rotational history. We confirmed that many of the target stars rapidly rotate and have high Li abundance, compared with the Sun, as suggested by many previous studies. There are, however, also some target stars that rotate slowly ($v\sin i$=2$\sim$3 km s$^{-1}$) and have low Li abundance like the Sun. These results support that old and slowly rotating stars similar to the Sun could have high activity level and large starspots. This is consistent with the results of our previous studies of solar-type superflare stars. In the future, it is important to conduct long-term monitoring observations of these active solar-analog stars in order to investigate detailed properties of large starspots from the viewpoint of stellar dynamo theory.

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{\AA} 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.

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.

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.

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.

We have measured an annual parallax of the Mira variable R~Ursae~Majoris (R~UMa) with the VLBI exploration for Radio Astronomy (VERA). From the monitoring VLBI observations spanning about two years, we detected H$_2$O maser spots in the LSR velocities ranges from 37 to 42 km\,s$^{-1}$. We derived an annual parallax of 1.97$\pm$0.05\,mas, and it gives a corresponding distance of 508$\pm$13\,pc. The VLBI maps revealed 72 maser spots distributed in $\sim$110 au area around an expected stellar position. Circumstellar kinematics of the maser spots were also revealed by subtracting a systemic motion in the Hipparcos catalog from proper motions of each maser spots derived from our VLBI observations. Infrared photometry is also conducted to measure a $K$ band apparent magnitude, and we obtained a mean magnitude of $m_K$ = 1.19$\pm$0.02\,mag. Using the trigonometric distance, the $m_K$ is converted to a $K$ band absolute magnitude of $M_K = -$7.34$\pm$0.06\,mag. This result gives a much more accurate absolute magnitude of R~UMa than previously provided. We solved a zero-point of $M_K - \log P$ relation for the Galactic Mira variables and obtained a relation of $M_K = -$3.52 $\log P$ + (1.09 $\pm$ 0.14). Other long period variables including red supergiants, whose distances were determined from astrometric VLBI, were also compiled to explore the different sequences of $M_K - \log P$ relation.

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.

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.

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.

The recent observations of Singh et al. (2012) have shown multiple plasma ejections and the intermittent nature of magnetic reconnection in the solar chromosphere, highlighting the need for fast reconnection to occur in highly collisional plasma. However, the physical process through which fast magnetic reconnection occurs in partially ionized plasma, like the solar chromosphere, is still poorly understood. It has been shown that for sufficiently high magnetic Reynolds numbers, Sweet-Parker current sheets can become unstable leading to tearing mode instability and plasmoid formation, but when dealing with a partially ionized plasma the strength of coupling between the ions and neutrals plays a fundamental role in determining the dynamics of the system. We propose that as the reconnecting current sheet thins and the tearing instability develops, plasmoid formation passes through strongly, intermediately, and weakly coupled (or decoupled) regimes, with the time scale for the tearing mode instability depending on the frictional coupling between ions and neutrals. We present calculations for the relevant time scales for fractal tearing in all three regimes. We show that as a result of the tearing mode instability and the subsequent non-linear instability due to the plasmoid-dominated reconnection, the Sweet-Parker current sheet tends to have a fractal-like structure, and when the chromospheric magnetic field is sufficiently strong the tearing instability can reach down to kinetic scales, which are hypothesized to be necessary for fast reconnection.

We searched for superflares on solar-type stars using the Kepler short-cadence (1-min sampling) data 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. 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 ^{-\alpha}$) with $\alpha=1.5$ for $10^{33}<E<10^{36}$ erg. 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 amplitude of brightness variations. We also found that the duration of superflares ($\tau$) increases with the flare energy ($E$) as $\tau \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.

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 $\delta$-spot regions that have a strong magnetic shear and emerge with polarity orientations not following Hale-Nicholson and Joy Laws.

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$\sigma$, 2.1$\sigma$, 3.8$\sigma$ confidence levels, respectively. The best three TSS values were $0.75\pm0.07$, $0.48\pm0.02$, and $0.56\pm0.04$, respectively.

We carried out spectroscopic observations with Subaru/HDS of 50 solar-type superflare stars found from Kepler data. More than half (34 stars) of the target stars show no evidence of the binary system, and we confirmed atmospheric parameters of these stars are roughly in the range of solar-type stars. We then conducted the detailed analyses for these 34 stars. First, the value of the "$v\sin i$" (projected rotational velocity) measured from spectroscopic results is consistent with the rotational velocity estimated from the brightness variation. Second, there is a correlation between the amplitude of the brightness variation and the intensity of Ca II IR triplet line. All the targets expected to have large starspots because of their large amplitude of the brightness variation show high chromospheric activities compared with the Sun. These results support that the brightness variation of superflare stars is explained by the rotation of a star with large starspots.

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.

We report the trigonometric parallax of IRAS 07427-2400 with VERA to be 0.185 $\pm$ 0.027 mas, corresponding to a distance of 5.41$^{+0.92}_{-0.69}$ kpc. The result is consistent with the previous result of 5.32$^{+0.49}_{-0.42}$ kpc obtained by Choi et al. (2014) within error. To remove the effect of internal maser motions (e.g., random motions), we observed six maser features associated with IRAS 07427-2400 and determined systematic proper motions of the source by averaging proper motions of the six maser features. The obtained proper motions are ($\mu_{\alpha}$cos$\delta$, $\mu_{\delta}$) = ($-$1.79 $\pm$ 0.32, 2.60 $\pm$ 0.17) mas yr$^{-1}$ in equatorial coordinates, while Choi et al. (2014) showed ($\mu_{\alpha}$cos$\delta$, $\mu_{\delta}$) = ($-$2.43 $\pm$ 0.02, 2.49 $\pm$ 0.09) mas yr$^{-1}$ with one maser feature. Our astrometry results place the source in the Perseus arm, the nearest main arm in the Milky Way. Using our result with previous astrometry results obtained from observations of the Perseus arm, we conducted direct (quantitative) comparisons between 27 astrometry results and an analytic gas dynamics model based on the density-wave theory and obtained two results. First is the pitch angle of the Perseus arm determined by VLBI astrometry, 11.1 $\pm$ 1.4 deg, differing from what is determined by the spiral potential model (probably traced by stars), $\sim$ 20 deg. The second is an offset between a dense gas region and the bottom of the spiral potential model. The dense gas region traced by VLBI astrometry is located downstream of the spiral potential model, which was previously confirmed in the nearby grand-design spiral galaxy M51 in Egusa et al. (2011).

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-{\beta} 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.

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.

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^{-\alpha}$) with $\alpha \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 ($\tau$) increases with the flare energy ($E$) as $\tau \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.

We report results of the performance evaluation of a new hardware correlator in Korea, the Daejeon correlator, developed by the Korea Astronomy and Space Science Institute (KASI) and the National Astronomical Observatory of Japan (NAOJ). We conducted Very Long Baseline Interferometry (VLBI) observations at 22~GHz with the Korean VLBI Network (KVN) in Korea and the VLBI Exploration of Radio Astrometry (VERA) in Japan, and correlated the aquired data with the Daejeon correlator. For evaluating the performance of the new hardware correlator, we compared the correlation outputs from the Daejeon correlator for KVN observations with those from a software correlator, the Distributed FX (DiFX). We investigated the correlated flux densities and brightness distributions of extragalactic compact radio sources. The comparison of the two correlator outputs show that they are consistent with each other within $<8\%$, which is comparable with the amplitude calibration uncertainties of KVN observations at 22~GHz. We also found that the 8\% difference in flux density is caused mainly by (a) the difference in the way of fringe phase tracking between the DiFX software correlator and the Daejeon hardware correlator, and (b) an unusual pattern (a double-layer pattern) of the amplitude correlation output from the Daejeon correlator. The visibility amplitude loss by the double-layer pattern is as small as 3\%. We conclude that the new hardware correlator produces reasonable correlation outputs for continuum observations, which are consistent with the outputs from the DiFX software correlator.

We conducted astrometric VLBI observations of water-vapor maser emission in the massive star forming region IRAS 21379+5106 to measure the annual parallax and proper motion, using VERA. The annual parallax was measured to be $0.262 \pm 0.031$ mas corresponding to a trigonometric distance of $3.82^{+0.51}_{-0.41}$ kpc. The proper motion was $(\mu_\alpha\cos{\delta}, \mu_\delta)=(-2.74 \pm 0.08, -2.87 \pm 0.18)$ mas yr$^{-1}$. Using this result, the Galactic rotational velocity was estimated to be $V_\theta=218\pm 19$ km s$^{-1}$ at the Galactocentric distance $R=9.22\pm0.43$ kpc, when we adopted the Galactic constants $R_0=8.05\pm 0.45$ kpc and $V_0=238\pm 14$ km s$^{-1}$. With newly determined distance, {the bolometric luminosity of the central young stellar object was re-evaluated to $(2.15\pm 0.54)\times 10^3 L_\odot$, which corresponds to spectral type of} B2--B3. Maser features were found to be distributed along a straight line from south-west to north-east. In addition, a vector map of the internal motions constructed from the residual proper motions implies that maser features trace a bipolar flow and that it cannot be explained by simple ballistic motion.

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 ($(\Delta T_{\rm{eff}})_{\rm{rms}} \sim 219$K, $(\Delta \log g)_{\rm{rms}} \sim 0.37$ dex, and $(\Delta\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.

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.

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$\alpha$ 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.

We have developed a software-based polarization spectrometer, PolariS, to acquire full-Stokes spectra with a very high spectral resolution of 61 Hz. The primary aim of PolariS is to measure the magnetic fields in dense star-forming cores by detecting the Zeeman splitting of molecular emission lines. The spectrometer consists of a commercially available digital sampler and a Linux computer. The computer is equipped with a graphics processing unit (GPU) to process FFT and cross-correlation using the CUDA (Compute Unified Device Architecture) library developed by NVIDIA. Thanks to a high degree of precision in quantization of the analog-to-digital converter and arithmetic in the GPU, PolariS offers excellent performances in linearity, dynamic range, sensitivity, bandpass flatness and stability. The software has been released under the MIT License and is available to the public. In this paper, we report the design of PolariS and its performance verified through engineering tests and commissioning observations.

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.

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

We conducted phase referencing VLBI observations of the Mira variable T~Lepus (T~Lep) using VERA, from 2003 to 2006. The distance to the source was determined from its annual parallax which was measured to be 3.06$\pm$0.04 mas, corresponding to a distance of 327$\pm$4\,pc. Our observations revealed the distribution and internal kinematics of H$_2$O masers in T~Lep, and we derived a source systemic motions of 14.60$\pm$0.50 mas yr$^{-1}$ and $-$35.43$\pm$0.79 mas yr$^{-1}$ in right ascension and declination, respectively. We also determined a LSR velocity of $V_\mathrm{LSR}^{\ast} = -$27.63 km s$^{-1}$. Comparison of our result with an image recently obtained from the VLTI infrared interferometer reveals a linear scale picture of the circumstellar structure of T~Lep. Analysis of the source systemic motion in the Galacto-centric coordinate frame indicates a large peculiar motion, which is consistent with the general characteristics of AGB stars. This source makes a contribution to the calibration of the period-luminosity relation of Galactic Mira variables. From the compilation of data for nearby Mira variables found in the literature, whose distances were derived from astrometric VLBI observations, we have calibrated the Galactic Mira period-luminosity relation to a high degree of accuracy.

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.

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.

Multi-epoch phase-referencing VLBI (Very Long Baseline Interferometry) observations with VERA (VLBI Exploration of Radio Astrometry) were performed for the symbiotic star R Aquarii (R Aqr) from September 2005 to Oct 2006. Tracing one of the $v=2$, $J=1-0$ SiO maser spots, we measured an annual parallax of $\pi = 4.59\pm0.24$ mas, corresponding to a distance of $218_{-11}^{+12}$ pc. Our result is consistent with earlier distance measurements, but yields the highest accuracy of about $5\%$ level. Applying our distance, we derived an absolute K-band magnitude of $M_{\mathrm{K}} = -7.71 \pm 0.11$, which is consistent with the recent Period-Luminosity relation by VLBI parallax measurements for 5 OH-Mira variables. In addition, the expansion age of an inner nebulae around R Aqr is found to be about 240 years, corresponds to about the year 1773.

Very long baseline interferometry (VLBI) monitoring of the 6.7 GHz methanol maser allows us to measure the internal proper motions of the maser spots and therefore study the gas motion around high-mass young stellar objects. To this end, we have begun monitoring observations with the East-Asian VLBI Network. In this paper we present the results of the first epoch observation for 36 sources, including 35 VLBI images of the methanol maser. Since two independent sources were found in three images, respectively, images of 38 sources were obtained. In 34 sources, more than or equal to 10 spots were detected. The observed spatial scale of the maser distribution was from 9 to 4900 astronomical units, and the following morphological categories were observed: elliptical, arched, linear, paired, and complex. The position of the maser spot was determined to an accuracy of approximately 0.1 mas, sufficiently high to measure the internal proper motion from two years of monitoring observations. The VLBI observation, however, detected only approximately 20% of all maser emission, suggesting that the remaining 80% of the total flux was spread into an undetectable extended distribution. Therefore, in addition to high-resolution observations, it is important to observe the whole structure of the maser emission including extended low-brightness structures, to reveal the associated site of the maser and gas motion.

We present the first measurement of the absolute proper motions of IRAS 00259+5625 (CB3, LBN594) associated with the HI loop called the "NGC281 superbubble" that extends from the Galactic plane over ~300 pc toward decreasing galactic latitude. The proper motion components measured with VERA are (mu_alpha cos(delta), mu_delta) = (-2.48 +/- 0.32, -2.85 +/- 0.65) mas yr^{-1}, converted into (mu_l cos(b), mu_b) = (-2.72 +/- 0.32, -2.62 +/- 0.65) mas yr^{-1} in the Galactic coordinates. The measured proper motion perpendicular to the Galactic plane (mu_b) shows vertical motion away from the Galactic plane with a significance of about ~4-sigma. As for the source distance, the distance measured with VERA is marginal, 2.4^{+1.0}_{-0.6} kpc. Using the distance, an absolute vertical motion (v_{b}) of -17.9 +/- 12.2 km s^{-1} is determined with ~1.5-sigma significance. The tendency of the large vertical motion is consistent with previous VLBI results for NGC 281 associated with the same superbubble. Thus, our VLBI results indicate the superbubble expansion motion whose origin is believed to be sequential supernova explosions.

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.

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^{-\alpha}, where \alpha ~ 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.

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 $\beta$ 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.

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.

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.

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.

We propose the particle acceleration model coupled with multiple plasmoid ejections in a solar flare. Unsteady reconnection produces plasmoids in a current sheet and ejects them out to the fast shocks, where particles in a plasmoid are reflected upstream the shock front by magnetic mirror effect. As the plasmoid passes through the shock front, the reflection distance becomes shorter and shorter driving Fermi acceleration, until it becomes proton Larmor radius. The fractal distribution of plasmoids may also have a role in naturally explaining the power-law spectrum in nonthermal emissions.

We report a detailed examination of the fine structure inside flare ribbons and the temporal evolution of this fine structure during the X2.5 solar flare that occurred on 2004 November 10. We examine elementary bursts of the C IV (1550{\AA}) emission lines seen as local transient brightenings inside the flare ribbons in the ultraviolet (1600{\AA}) images taken with Transition Region and Coronal Explorer, and we call them C IV kernels. This flare was also observed in Ha with the Sartorius 18 cm Refractor telescope at Kwasan observatory, Kyoto University, and in hard X-rays (HXR) with Reuven Ramaty High Energy Solar Spectroscopic Imager. Many C IV kernels, whose sizes were comparable to or less than 2", were found to brighten successively during the evolution of the flare ribbon. The majority of them were well correlated with the Ha kernels in both space and time, while some of them were associated with the HXR emission. These kernels were thought to be caused by the precipitation of nonthermal particles at the footpoints of the reconnecting flare loops. The time profiles of the C IV kernels showed intermittent bursts, whose peak intensity, duration, and time interval were well described by power-law distribution functions. This result is interpreted as evidence for "self-organized criticality" in avalanching behavior in a single flare event, or for fractal current sheets in the impulsive reconnection region.

The Soft X-ray Telescope (SXT) on board Yohkoh revealed that the ejection of X-ray emitting plasmoid is sometimes observed in a solar flare. It was found that the ejected plasmoid is strongly accelerated during a peak in the hard X-ray emission of the flare. In this paper we present an examination of the GOES X 2.3 class flare that occurred at 14.51 UT on 2000 November 24. In the SXT images we found multiple plasmoid ejections with velocities in the range of 250-1500 km/s, which showed blob-like or loop-like structures. Furthermore, we also found that each plasmoid ejection is associated with an impulsive burst of hard X-ray emission. Although some correlation between plasmoid ejection and hard X-ray emission has been discussed previously, our observation shows similar behavior for multiple plasmoid ejection such that each plasmoid ejection occurs during the strong energy release of the solar flare. As a result of temperature-emission measure analysis of such plasmoids, it was revealed that the apparent velocities and kinetic energies of the plasmoid ejections show a correlation with the peak intensities in the hard X-ray emissions.

This paper reports very-long-baseline interferometry observations of the radio-loud broad absorption line (BAL) quasar J1020+4320 at 1.7, 2.3, 6.7, and 8.4 GHz using the Japanese VLBI network (JVN) and European VLBI network (EVN). The radio morphology is compact with a size of ~10 pc. The convex radio spectrum is stable over the last decade; an observed peak frequency of 3.2 GHz is equivalent to 9.5 GHz in the rest frame, suggesting an age of the order of ~100 years as a radio source, according to an observed correlation between linear size and peak frequency of compact steep spectrum (CSS) and giga-hertz peaked spectrum (GPS) radio sources. A low-frequency radio excess suggests relic of past jet activity. J1020+4320 may be one of the quasars with recurrent and short-lived jet activity during a BAL-outflowing phase.

Recent observations of solar type stars with the Kepler satellite by Maehara et al. have revealed the existence of superflares (with energy of 10^33 - 10^35 erg) on Sun-like stars, which are similar to our Sun in their surface temperature (5600 K - 6000 K) and slow rotation (rotational period > 10 days). From the statistical analysis of these superflares, it was found that superflares with energy 10^34 erg occur once in 800 years and superflares with 10^35 erg occur once in 5000 years on Sun-like stars. In this paper, we examine whether superflares with energy of 10^33 - 10^35 erg could occur on the present Sun through the use of simple order-of-magnitude estimates based on current ideas relating to the mechanisms of the solar dynamo.

We present analyses to determine the fundamental parameters of the Galaxy based on VLBI astrometry of 52 Galactic maser sources obtained with VERA, VLBA and EVN. We model the Galaxy's structure with a set of parameters including the Galaxy center distance R_0, the angular rotation velocity at the LSR Omega_0, mean peculiar motion of the sources with respect to Galactic rotation (U_src, V_src, W_src), rotation-curve shape index, and the V component of the Solar peculiar motions V_sun. Based on a Markov chain Monte Carlo method, we find that the Galaxy center distance is constrained at a 5% level to be R_0 = 8.05 +/- 0.45 kpc, where the error bar includes both statistical and systematic errors. We also find that the two components of the source peculiar motion U_src and W_src are fairly small compared to the Galactic rotation velocity, being U_src = 1.0 +/- 1.5 km/s and W_src = -1.4 +/- 1.2 km/s. Also, the rotation curve shape is found to be basically flat between Galacto-centric radii of 4 and 13 kpc. On the other hand, we find a linear relation between V_src and V_sun as V_src = V_sun -19 (+/- 2) km/s, suggesting that the value of V_src is fully dependent on the adopted value of V_sun. Regarding the rotation speed in the vicinity of the Sun, we also find a strong correlation between Omega_0 and V_sun. We find that the angular velocity of the Sun, Omega_sun, which is defined as Omega_sun = Omega_0 + V_sun/R_0, can be well constrained with the best estimate of Omega_sun = 31.09 +/- 0.78 km/s/kpc. This corresponds to Theta_0 = 238 +/- 14 km/s if one adopts the above value of R_0 and recent determination of V_sun ~ 12 km/s.

We introduce the GENJI program (Gamma-ray Emitting Notable AGN Monitoring by Japanese VLBI), which is a monitoring program of gamma-ray bright AGNs with the VERA array (VLBI Exploration of Radio Astrometry). The GENJI programme aims a dense monitoring at 22 GHz towards the $\gamma$-ray emitting active galactic nuclei (AGNs) to investigate the radio time variation of the core and possible ejection of new radio component, motion of jets, and their relation with the emission at other wavelengths especially in $\gamma$-rays. Currently we are monitoring 8 $\gamma$-ray-emitting notable AGNs (DA 55, 3C 84, M 87, PKS 1510-089, DA 406, NRAO 530, BL Lac, 3C 454.3) about once every two weeks. This programme is promising to trace the trend of radio time variation on shorter timescale than conventional VLBI monitoring programme and to provide complimentary data with them (e.g., MOJAVE, Boston University Blazar Project). In particular, we successfully coordinated quick follow-up observations after the GeV $\gamma$-ray flare in NRAO 530 and 3C 454.3 reported by the Fermi Gamma-ray Space Telescope. Here we present the initial results of morphology and light curves for the first 7-month operation.

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 $\tau_{\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 $\lambda$ evolves according to a self-similar relation $\lambda \propto t^{\alpha/2}$ where $\alpha$ is the power-law index of the pressure distribution.

We report measurement of trigonometric parallax of IRAS 05168+3634 with VERA. The parallax is 0.532 +/- 0.053 mas, corresponding to a distance of 1.88+0.21/-0.17 kpc. This result is significantly smaller than the previous distance estimate of 6 kpc based on kinematic distance. This drastic change in the source distance revises not only physical parameters of IRAS 05168+3634, but also its location of the source, placing it in the Perseus arm rather than the Outer arm. We also measure proper motions of the source. A combination of the distance and the proper motions with systemic velocity yields rotation velocity ({\Theta}) of 227+9/-11 km s-1 at the source, assuming {\Theta}0 = 240 km s-1. Our result combined with previous VLBI results for six sources in the Perseus arm indicates that the sources rotate systematically slower than the Galactic rotation velocity at the LSR. In fact, we show observed disk peculiar motions averaged over the seven sources in the Perseus arm as (Umean, Vmean) = (11 +/- 3, -17 +/- 3) km s-1, indicating that these seven sources are systematically moving toward the Galactic center, and lag behind the Galactic rotation.

Solar spicules are the fundamental magnetic structures in the chromosphere and considered to play a key role in channelling the chromosphere and corona. Recently, it was suggested by De Pontieu et al. that there were two types of spicules with very different dynamic properties, which were detected by space- time plot technique in the Ca ii H line (3968 A) wavelength from Hinode/SOT observations. 'Type I' spicule, with a 3-7 minute lifetime, undergoes a cycle of upward and downward motion; in contrast, 'Type II' spicule fades away within dozens of seconds, without descending phase. We are motivated by the fact that for a spicule with complicated 3D motion, the space-time plot, which is made through a slit on a fixed position, could not match the spicule behavior all the time and might lose its real life story. By revisiting the same data sets, we identify and trace 105 and 102 spicules in quiet sun (QS) and coronal hole (CH), respectively, and obtain their statistical dynamic properties. First, we have not found a single convincing example of 'Type II' spicules. Secondly, more than 60% of the identified spicules in each region show a complete cycle, i.e., majority spicules are 'Type I'. Thirdly, the lifetime of spicules in QS and CH are 148 s and 112 s, respectively, but there is no fundamental lifetime difference between the spicules in QS and CH reported earlier. Therefore, the suggestion of coronal heating by 'Type II' spicules should be taken with cautions. Subject headings: Sun: chromosphere Sun:transition region Sun:corona

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.

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.

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.

We observed at 22 GHz with the VLBI array VERA a sample of 1536 sources with correlated flux densities brighter than 200 mJy at 8 GHz. One half of target sources has been detected. The detection limit was around 200 mJy. We derived the correlated flux densities of 877 detected sources in three ranges of projected baseline lengths. The objective of these observations was to determine the suitability of given sources as phase calibrators for dual-beam and phase-referencing observations at high frequencies. Preliminary results indicate that the number of compact extragalactic sources at 22 GHz brighter than a given correlated flux density level is twice less than at 8 GHz.

We present a numerical investigation of the coronal evolution of a coronal mass ejection (CME) on 2005 August 22 using a 3-D thermodynamics magnetohydrodynamic model, the SWMF. The source region of the eruption was anemone active region (AR) 10798, which emerged inside a coronal hole. We validate our modeled corona by producing synthetic extreme ultraviolet (EUV) images, which we compare to EIT images. We initiate the CME with an out-of-equilibrium flux rope with an orientation and chirality chosen in agreement with observations of a H-alpha filament. During the eruption, one footpoint of the flux rope reconnects with streamer magnetic field lines and with open field lines from the adjacent coronal hole. It yields an eruption which has a mix of closed and open twisted field lines due to interchange reconnection and only one footpoint line-tied to the source region. Even with the large-scale reconnection, we find no evidence of strong rotation of the CME as it propagates. We study the CME deflection and find that the effect of the Lorentz force is a deflection of the CME by about 3 deg/Rsun towards the East during the first 30 minutes of the propagation. We also produce coronagraphic and EUV images of the CME, which we compare with real images, identifying a dimming region associated with the reconnection process. We discuss the implication of our results for the arrival at Earth of CMEs originating from the limb and for models to explain the presence of open field lines in magnetic clouds.

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\"{u}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.

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

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 $\rho(r) \propto r^{- \alpha}$ 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\'en 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 $\Gamma_{\rm sh} \propto r_{\rm sh}$, where $\Gamma_{\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.

We performed three dimensional resistive magnetohydrodynamic simulations to study the magnetic reconnection using an initially shearing magnetic field configuration (force free field with a current sheet in the middle of the computational box). It is shown that there are two types of reconnection jets: the ordinary reconnection jets and fan-shaped jets, which are formed along the guide magnetic field. The fan-shaped jets are much different from the ordinary reconnection jets which are ejected by magnetic tension force. There are two driving forces for accelerating the fan-shaped jets. The one is the Lorentz force which dominates the motion of fluid elements at first and then the gas pressure gradient force accelerates the fluid elements in the later stage. The dependence on magnetic reconnection angle and resistivity value has also been studied. The formation and evolution of these jets provide a new understanding of dynamic magnetohydrodynamic jets.

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.

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.

We present the results from 1.5D diffusion simulations of the Kippenhahn-Schlueter prominence model magnetic field evolution under the influence of the ambipolar terms of Cowling resistivity. We show that initially the evolution is determined by the ratio of the horizontal and vertical magnetic fields, which gives current sheet thinning (thickening) when this ratio is large (small) and a marginal case where a new characteristic current sheet length scale is formed. After a timespan greater than the Cowling resistivity time, the current sheet thickens as a power law of $t$ independent of the ratio of the field strengths. These results imply that when Cowling resistivity is included in the model, the tearing instability time scale is reduced by more than one order of magnitude when the ratio of the horizontal field to the vertical field is 20\% or less. These results imply that, over the course of its lifetime, the structure of the prominence can be significantly altered by Cowling resistivity, and in some cases will allow the tearing instability to occur.

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.

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.

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.

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.

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).

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.

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.

We report on trigonometric parallax measurements for the Galactic star forming region G14.33-0.64 toward the Sagittarius spiral arm. We conducted multi-epoch phase-referencing observations of an H2O maser source in G14.33-0.64 with the Japanese VLBI array VERA. We successfully detected a parallax of 0.893+/-0.101 mas, corresponding to a source distance of 1.12+/-0.13 kpc, which is less than half of the kinematic distance for G14.33-0.64. Our new distance measurement demonstrates that the Sagittarius arm lies at a closer distance of ~1 kpc, instead of previously assumed ~2-3 kpc from kinematic distances. The previously suggested deviation of the Sagittarius arm toward the Galactic center from the symmetrically fitted model (Taylor & Cordes 1993) is likely due to large errors of kinematic distances at low galactic longitudes. G14.33-0.64 most likely traces the near side of the Sagittarius arm. We attempted fitting the pitch angle of the arm with other parallax measurements along the arm, which yielded two possible pitch angles of i=34.7+/-2.7 degrees and i=11.2+/-10.5 degrees. Our proper motion measurements suggest G14.33-0.64 has no significant peculiar motion relative to the differential rotation of the Galaxy (assumed to be in a circular orbit), indicating that the source motion is in good agreement with the Galactic rotation.

Reported observations in H-alpha, Ca II H and K or or other chromospheric lines of coronal rain trace back to the days of the Skylab mission. Offering a high contrast in intensity with respect to the background (either bright in emission if observed at the limb, or dark in absorption if observed on disk) these cool blobs are often observed falling down from high coronal heights above active regions. A physical explanation for this spectacular phenomenon has been put forward thanks to numerical simulations of loops with footpoint concentrated heating, a heating scenario in which cool condensations naturally form in the corona. This effect has been termed 'catastrophic cooling' and is the predominant explanation for coronal rain. In this work we further investigate the link between this phenomenon and the heating mechanisms acting in the corona. We start by analyzing observations of coronal rain at the limb in the Ca II H line performed by the SOT instrument on board of the Hinode satellite. We then compare the observations with 1.5-dimensional MHD simulations of loops being heated by small-scale discrete events concentrated towards the footpoints (that could come, for instance, from magnetic reconnection events), and by Alfven waves generated at the photosphere. It is found that if a loop is heated predominantly from Alfven waves coronal rain is inhibited due to the characteristic uniform heating they produce. Hence coronal rain may not only point to the spatial distribution of the heating in coronal loops but also to the agent of the heating itself. We thus propose coronal rain as a marker for coronal heating mechanisms.

In the context of coronal heating, among the zoo of MHD waves that exist in the solar atmosphere, Alfven waves receive special attention. Indeed, these waves constitute an attractive heating agent due to their ability to carry over the many different layers of the solar atmosphere sufficient energy to heat and maintain a corona. However, due to their incompressible nature these waves need a mechanism such as mode conversion (leading to shock heating), phase mixing, resonant absorption or turbulent cascade in order to heat the plasma. New observations with polarimetric, spectroscopic and imaging instruments such as those on board of the japanese satellite Hinode, or the SST or CoMP, are bringing strong evidence for the existence of energetic Alfven waves in the solar corona. In order to assess the role of Alfven waves in coronal heating, in this work we model a magnetic flux tube being subject to Alfven wave heating through the mode conversion mechanism. Using a 1.5-dimensional MHD code we carry out a parameter survey varying the magnetic flux tube geometry (length and expansion), the photospheric magnetic field, the photospheric velocity amplitudes and the nature of the waves (monochromatic or white noise spectrum). It is found that independently of the photospheric wave amplitude and magnetic field a corona can be produced and maintained only for long (> 80 Mm) and thick (area ratio between photosphere and corona > 500) loops. Above a critical value of the photospheric velocity amplitude (generally a few km/s) the corona can no longer be maintained over extended periods of time and collapses due to the large momentum of the waves. These results establish several constraints on Alfven wave heating as a coronal heating mechanism, especially for active region loops.

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.

Alfven waves created by sub-photospheric motions or by magnetic reconnection in the low solar atmosphere seem good candidates for coronal heating. However, the corona is also likely to be heated more directly by magnetic reconnection, with dissipation taking place in current sheets. Distinguishing observationally between these two heating mechanisms is an extremely difficult task. We perform 1.5-dimensional MHD simulations of a coronal loop subject to each type of heating and derive observational quantities that may allow these to be differentiated.

We present multi-frequency simultaneous VLBA observations at 15, 22 and 43 GHz towards the nucleus of the nearby radio galaxy NGC 1052. These three continuum images reveal a double-sided jet structure, whose relative intensity ratios imply that the jet axis is oriented close to the sky plane. The steeply rising spectra at 15-43 GHz at the inner edges of the jets strongly suggest that synchrotron emission is absorbed by foreground thermal plasma. We detected H_2O maser emission in the velocity range of 1550-1850 km/s, which is redshifted by 50-350 km/s with respect to the systemic velocity of NGC 1052. The redshifted maser gas appears projected against both sides of the jet, similar to the HI seen in absorption. The H_2O maser gas is located where the free-free absorption opacity is large. This probably implies that the masers in NGC 1052 are associated with a circumnuclear torus or disk as in the nucleus of NGC 4258. Such circumnuclear structure can be the sense of accreting onto the central engine.

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$\alpha$ 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.

The heating of the solar chromosphere and corona is a long-standing puzzle in solar physics. Hinode observations show the ubiquitous presence of chromospheric anemone jets outside sunspots in active regions. They are typically 3 to 7 arc seconds = 2000 to 5000 kilometers long and 0.2 to 0.4 arc second = 150 to 300 kilometers wide, and their velocity is 10 to 20 kilometers per second. These small jets have an inverted Y-shape, similar to the shape of x-ray anemone jets in the corona. These features imply that magnetic reconnection similar to that in the corona is occurring at a much smaller spatial scale throughout the chromosphere and suggest that the heating of the solar chromosphere and corona may be related to small-scale ubiquitous reconnection.

Hinode discovered a beautiful giant jet with both cool and hot components at the solar limb on 2007 February 9. Simultaneous observations by the Hinode SOT, XRT, and TRACE 195 satellites revealed that hot (5x10^6 K) and cool (10^4 K) jets were located side by side and that the hot jet preceded the associated cool jet (1-2 minutes). A current-sheet-like structure was seen in optical (Ca IIH), EUV (195A), and soft X-ray emissions, suggesting that magnetic reconnection is occurring in the transition region or upper chromosphere. Alfven waves were also observed with Hinode SOT. These propagated along the jet at velocities of 200 km/s with amplitudes (transverse velocity) of 5-15 km/s and a period of 200 s. We performed two-dimensional MHD simulation of the jets on the basis of the emerging flux-reconnection model, by extending Yokoyama and Shibata's model. We extended the model with a more realistic initial condition (10^6 K corona) and compared our model with multiwavelength observations. The improvement of the coronal temperature and density in the simulation model allowed for the first time the reproduction of the structure and evolution of both the cool and hot jets quantitatively, supporting the magnetic reconnection model. The generation and the propagation of Alfven waves are also reproduced self-consistently in the simulation model.

The plasmoid-induced-reconnection model explaining solar flares based on bursty reconnection produced by an ejecting plasmoid suggests a possible relation between the ejection velocity of a plasmoid and the rate of magnetic reconnection. In this study, we focus on the quantitative description of this relation. We performed magnetohydrodynamic (MHD) simulations of solar flares by changing the values of resistivity and the plasmoid velocity. The plasmoid velocity has been changed by applying an additional force to the plasmoid to see how the plasmoid velocity affects the reconnection rate. An important result is that the reconnection rate has a positive correlation with the plasmoid velocity, which is consistent with the plasmoid-induced-reconnection model for solar flares. We also discuss an observational result supporting this positive correlation.

We have made multi-epoch VLBI observations of H2O maser emission in the massive star forming region IRAS 06061+2151 with the Japanese VLBI network (JVN) from 2005 May to 2007 October. The detected maser features are distributed within an 1\arcsec$\times$1\arcsec (2000 au$\times$2000 au at the source position) around the ultra-compact H {\small\bf II} region seen in radio continuum emission. Their bipolar morphology and expanding motion traced through their relative proper motions indicate that they are excited by an energetic bipolar outflow. Our three-dimensional model fitting has shown that the maser kinematical structure in IRAS 06061+2151 is able to be explained by a biconical outflow with a large opening angle ($>$ 50\degr). The position angle of the flow major axis coincides very well with that of the large scale jet seen in 2.1$\:\mu\rmn{m}$ hydrogen emission. This maser geometry indicates the existence of dual structures composed of a collimated jet and a less collimated massive molecular flow. We have also detected a large velocity gradient in the southern maser group. This can be explained by a very small (on a scale of several tens of au) and clumpy (the density contrast by an order of magnitude or more) structure of the parental cloud. Such a structure may be formed by strong instability of shock front or splitting of high density core.

We present a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred to the quasar J1147-0724 located at an angular separation of 1.23$^{\circ}$. In observations spanning nearly two years, we have detected the maser spot at the LSR velocity of 34.7 km s$^{-1}$, for which we measured the annual parallax of 2.33$\pm$0.13 mas corresponding to a distance of 430$^{+25}_{-23}$ pc. This measurement has an accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate the stellar radius of S Crt and compare it with those of Mira variables.

We report astrometric observations of H2O masers around the red supergiant VY Canis Majoris (VY CMa) carried out with VLBI Exploration of Radio Astrometry (VERA). Based on astrometric monitoring for 13 months, we successfully measured a trigonometric parallax of 0.88 +/- 0.08 mas, corresponding to a distance of 1.14 +0.11/-0.09 kpc. This is the most accurate distance to VY CMa and the first one based on an annual parallax measurement. The luminosity of VY CMa has been overestimated due to a previously accepted distance. With our result, we re-estimate the luminosity of VY CMa to be (3 +/- 0.5) x 10^5 L_sun using the bolometric flux integrated over optical and IR wavelengths. This improved luminosity value makes location of VY CMa on the Hertzsprung-Russel (HR) diagram much closer to the theoretically allowable zone (i.e. the left side of the Hayashi track) than previous ones, though uncertainty in the effective temperature of the stellar surface still does not permit us to make a final conclusion.

We present results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H$_{2}$O masers associated with a young stellar object (YSO) IRAS 22198+6336 in a dark cloud L1204G. Based on the phase-referencing VLBI astrometry, we derive an annual parallax of IRAS 22198+6336 to be 1.309$\pm$0.047 mas, corresponding to the distance of 764$\pm$27 pc from the Sun. Although the most principal error source of our astrometry is attributed to the internal structure of the maser spots, we successfully reduce the errors in the derived annual parallax by employing the position measurements for all of the 26 detected maser spots. Based on this result, we reanalyze the spectral energy distribution (SED) of IRAS 22198+6336 and find that the bolometric luminosity and total mass of IRAS 22198+6336 are 450$L_{\odot}$ and 7$M_{\odot}$, respectively. These values are consistent with an intermediate-mass YSO deeply embedded in the dense dust core, which has been proposed to be an intermediate-mass counterpart of a low-mass Class 0 source. In addition, we obtain absolute proper motions of the H$_{2}$O masers for the most blue-shifted components. We propose that the collimated jets aligned along the east-west direction are the most plausible explanation for the origin of the detected maser features.

We have used the Japanese VLBI array VERA to perform high-precision astrometry of an H2O maser source in the Galactic star-forming region NGC 281 West, which has been considered to be part of a 300-pc superbubble. We successfully detected a trigonometric parallax of 0.355+/-0.030 mas, corresponding to a source distance of 2.82+/-0.24 kpc. Our direct distance determination of NGC 281 has resolved the large distance discrepancy between previous photometric and kinematic studies; likely NGC 281 is in the far side of the Perseus spiral arm. The source distance as well as the absolute proper motions were used to demonstrate the 3D structure and expansion of the NGC 281 superbubble, ~650 pc in size parallel to the Galactic disk and with a shape slightly elongated along the disk or spherical, but not vertically elongated, indicating the superbubble expansion may be confined to the disk. We estimate the expansion velocity of the superbubble as ~20 km/s both perpendicular to and parallel to the Galactic disk with a consistent timescale of ~20 Myr.

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.

High spatial resolution observation of an emerging flux region (EFR) was done using a vector magnetograph and a H-alpha Lyot filtergraph with Domeless Solar Telescope at Hida Observatory on October 22, 2006. In H-alpha wing images, we could see many Ellerman bombs (EBs) in the EFR. Two observation modes, slit scan and slit fixed, were performed with the vector magnetograph, along with H-alpha filtergraph. Using H-alpha wing images, we detected 12 EBs during slit scan observation period and 9 EBs during slit fixed observation period. With slit scan observation, we found that all the EBs were distributed in the area where the spatial gradient of vertical field intensity was large, which indicates the possibility of rapid topological change in magnetic field in the area of EBs. With slit fixed observation, we found that EBs were distributed in the areas of undulatory magnetic fields, both in vertical and horizontal components. This paper is the first to report the undulatory pattern in horizontal components of magnetic field, which is also evidence for emerging magnetic flux by Parker instability. These results allow us to confirm the association between EBs and emerging flux tubes. Three triggering mechanism of EBs is discussed with respect to emerging flux tubes: 9 out of 21 EBs occurred at the footpoints of emerging flux tubes, 8 out of 21 EBs occurred at the top of emerging flux tubes, and 4 out of 21 EBs occurred at unipolar region. Each case can be explained by magnetic reconnection in the low chromosphere.

We present multi-frequency simultaneous VLBA observations at 15, 22 and 43 GHz towards the nucleus of the nearby radio galaxy NGC 1052. These three continuum images reveal a double-sided jet structure, whose relative intensity ratios imply that the jet axis is oriented close to the sky plane. The steeply rising spectra at 15-43 GHz at the inner edges of the jets strongly suggest that synchrotron emission is absorbed by foreground thermal plasma. We detected H2O maser emission in the velocity range of 1550-1850 km/s, which is redshifted by 50-350 km/s with respect to the systemic velocity of NGC 1052. The redshifted maser gas appears projected against both sides of the jet, in the same manner as the HI seen in absorption. The H2O maser gas are located where the free-free absorption opacity is large. This probably imply that the masers in NGC 1052 are associated with a circumnuclear torus or disk as in the nucleus of NGC 4258. Such circumnuclear structure can be the sence of accreting onto the central engine.

Continuous observations were obtained of active region 10953 with the Solar Optical Telescope (SOT) on board the \emph{Hinode} satellite during 2007 April 28 to May 9. A prominence was located over the polarity inversion line (PIL) in the south-east of the main sunspot. These observations provided us with a time series of vector magnetic fields on the photosphere under the prominence. We found four features: (1) The abutting opposite-polarity regions on the two sides along the PIL first grew laterally in size and then narrowed. (2) These abutting regions contained vertically-weak, but horizontally-strong magnetic fields. (3) The orientations of the horizontal magnetic fields along the PIL on the photosphere gradually changed with time from a normal-polarity configuration to a inverse-polarity one. (4) The horizontal-magnetic field region was blueshifted. These indicate that helical flux rope was emerging from below the photosphere into the corona along the PIL under the pre-existing prominence. We suggest that this supply of a helical magnetic flux into the corona is associated with evolution and maintenance of active-region prominences.

Solar prominences are cool 10$^4$ Kelvin plasma clouds supported in the surrounding 10$^6$ Kelvin coronal plasma by as-yet undetermined mechanisms. Observations from \emph{Hinode} show fine-scale threadlike structures oscillating in the plane of the sky with periods of several minutes. We suggest these transverse magnetohydrodynamic waves may represent Alfv\'en waves propagating on coronal magnetic field lines and these may play a role in heating the corona.

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.

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.

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.

We report on the results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio Astrometry) of the 22 GHz H2O masers associated with the young stellar object SVS 13 in the NGC 1333 region. We have carried out phase-referencing VLBI astrometry and measured an annual parallax of the maser features in SVS 13 of 4.25+/-0.32 mas, corresponding to the distance of 235+/-18 pc from the Sun. Our result is consistent with a photometric distance of 220 pc previously reported. Even though the maser features were detectable only for 6 months, the present results provide the distance to NGC 1333 with much higher accuracy than photometric methods. The absolute positions and proper motions have been derived, revealing that the H2O masers with the LSR (local standard of rest) velocities of 7-8 km s-1 are most likely associated with VLA4A, which is a radio counterpart of SVS 13. The origin of the observed proper motions of the maser features are currently difficult to attribute to either the jet or the rotating circumstellar disk associated with VLA4A, which should be investigated through future high-resolution astrometric observations of VLA4A and other radio sources in NGC 1333.

We have performed high-precision astrometry of H2O maser sources in Galactic star forming region Sharpless 269 (S269) with VERA. We have successfully detected a trigonometric parallax of 189+/-8 micro-arcsec, corresponding to the source distance of 5.28 +0.24/-0.22 kpc. This is the smallest parallax ever measured, and the first one detected beyond 5 kpc. The source distance as well as proper motions are used to constrain the outer rotation curve of the Galaxy, demonstrating that the difference of rotation velocities at the Sun and at S269 (which is 13.1 kpc away from the Galaxy's center) is less than 3%. This gives the strongest constraint on the flatness of the outer rotation curve and provides a direct confirmation on the existence of large amount of dark matter in the Galaxy's outer disk.

On 2005 September 13 a filament eruption accompanied by a halo CME occurred in the most flare-productive active region NOAA 10808 in Solar Cycle 23. Using multi-wavelength observations before the filament eruption on Sep. 13th, we investigate the processes leading to the catastrophic eruption. We find that the filament slowly ascended at a speed of 0.1km/s over two days before the eruption. During slow ascending, many small flares were observed close to the footpoints of the filament, where new magnetic elements were emerging. On the basis of the observational facts we discuss the triggering mechanism leading to the filament eruption. We suggest the process toward the eruption as follows: First, a series of small flares played a role in changing the topology of the loops overlying the filament. Second, the small flares gradually changed the equilibrium state of the filament and caused the filament to ascend slowly over two days. Finally, a C2.9 flare that occurred when the filament was close to the critical point for loss of equilibrium directly led to the catastrophic filament eruption right after itself.

We present the initial results of multi-epoch VLBI observations of the 22 GHz H2O masers in the Orion KL region with VERA (VLBI Exploration of Radio Astrometry). With the VERA dual-beam receiving system, we have carried out phase-referencing VLBI astrometry and successfully detected an annual parallax of Orion KL to be 2.29+/-0.10 mas, corresponding to the distance of 437+/-19 pc from the Sun. The distance to Orion KL is determined for the first time with the annual parallax method in these observations. Although this value is consistent with that of the previously reported, 480+/-80 pc, which is estimated from the statistical parallax method using proper motions and radial velocities of the H2O maser features, our new results provide the much more accurate value with an uncertainty of only 4%. In addition to the annual parallax, we have detected an absolute proper motion of the maser feature, suggesting an outflow motion powered by the radio source I along with the systematic motion of source I itself.

We report on absolute proper-motion measurements of an H2O maser source in the NGC 281 West molecular cloud, which is located ~320 pc above the Galactic plane and is associated with an HI loop extending from the Galactic plane. We have conducted multi-epoch phase-referencing observations of the maser source with VERA (VLBI Exploration of Radio Astrometry) over a monitoring period of 6 months since May 2006. We find that the H2O maser features in NGC 281 West are systematically moving toward the southwest and further away from the Galactic plane with a vertical velocity of ~20-30 km/s at its estimated distance of 2.2-3.5 kpc. Our new results provide the most direct evidence that the gas in the NGC 281 region on the HI loop was blown out from the Galactic plane, most likely in a superbubble driven by multiple or sequential supernova explosions in the Galactic plane.

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.

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.

Neutrino radiation takes a major role in the momentum, heat, and lepton transports in proto-neutron stars (PNSs). These diffusive processes affect the growth of magnetorotational instability (MRI) in PNSs. We perform a local linear analysis for the axisymmetric and nonaxisymmetric MRI including the effects of neutrino transports and ohmic dissipation. We find that the MRI can grow even in the multi-diffusive situations that are realized in neutrino loaded PNSs. When the toroidal magnetic component dominates over the poloidal one, nonaxisymmetric MRI modes grow much faster than axisymmetric modes. These results suggest the importance of the nonaxisymmetric MRI in PNSs. Thus the understandings of three-dimensional nonlinear evolutions of the MRI are necessary to reveal the explosion mechanism of core-collapse supernovae.

This paper presents results of a survey search for bright compact radio sources at 22 GHz with the VERA radio-interferometer. Each source from a list of 2494 objects was observed in one scan for 2 minutes. The purpose of this survey was to find compact extragalactic sources bright enough at 22 GHz to be useful as phase calibrators. Observed sources were either a) within 6 degrees of the Galactic plane, or b) within 11 degrees from the Galactic center; or c) within 2 degrees from known water masers. Among the observed sources, 549 were detected, including 180 extragalactic objects which were not previously observed with the very long baseline interferometry technique. Estimates of the correlated flux densities of the detected sources are presented.

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity nearly 0.3c) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3r_S. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfven waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star.

We performed magnetohydrodynamic simulation of coronal mass ejections (CMEs) and associated giant arcade formations, and the results suggested new interpretations of observations of CMEs. We performed two cases of the simulation: with and without heat conduction. Comparing between the results of the two cases, we found that reconnection rate in the conductive case is a little higher than that in the adiabatic case and the temperature of the loop top is consistent with the theoretical value predicted by the Yokoyama-Shibata scaling law. The dynamical properties such as velocity and magnetic fields are similar in the two cases, whereas thermal properties such as temperature and density are very different.In both cases, slow shocks associated with magnetic reconnectionpropagate from the reconnection region along the magnetic field lines around the flux rope, and the shock fronts form spiral patterns. Just outside the slow shocks, the plasma density decreased a great deal. The soft X-ray images synthesized from the numerical results are compared with the soft X-ray images of a giant arcade observed with the Soft X-ray Telescope aboard {\it Yohkoh}, it is confirmed that the effect of heat conduction is significant for the detailed comparison between simulation and observation. The comparison between synthesized and observed soft X-ray images provides new interpretations of various features associated with CMEs and giant arcades.

We report the results of 3-dimensional magnetohydrodynamic (MHD) simulations of a jet formation by the interaction between an accretion disk and a large scale magnetic field. The disk is not treated as a boundary condition but is solved self-consistently. To investigate the stability of MHD jet, the accretion disk is perturbed with a non-axisymmetric sinusoidal or random fluctuation of the rotational velocity. The dependences of the jet velocity $(v_z)$, mass outflow rate $(\dot{M}_w)$, and mass accretion rate $(\dot{M}_a)$ on the initial magnetic field strength in both non-axisymmetric cases are similar to those in the axisymmetric case. That is, $v_z \propto B_0^{1/3}$, $\dot{M}_w \propto B_0$ and $\dot{M}_a \propto B_0^{1.4}$ where $B_0$ is the initial magnetic field strength. The former two relations are consistent with the Michel's steady solution, $v_z \propto (B_0^2/\dot{M}_w)^{1/3}$, although the jet and accretion do not reach the steady state. In both perturbation cases, a non-axisymmetric structure with $m=2$ appears in the jet, where $m$ means the azimuthal wave number. This structure can not be explained by Kelvin-Helmholtz instability and seems to originate in the accretion disk. Non-axisymmetric modes in the jet reach almost constant levels after about 1.5 orbital periods of the accretion disk, while all modes in the accretion disk grow with oscillation. As for the angular momentum transport by Maxwell stress, the vertical component, $<B_\phi B_z/4\pi>$, is comparable to the radial component, $<B_\phi B_r/4\pi>$, in the wide range of initial magnetic field strength.

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.

By using the method presented by Isobe et al. (2002), the non-dimensional reconnection rate Vin/Va has been determined for the impulsive phase of three two-ribbon flares, where Vin is the velocity of the reconnection inflow and Va is the Alfven velocity. The non-dimensional reconnection rate is important to make a constraint on the theoretical models of magnetic reconnection. In order to reduce the uncertainty of the reconnection rate, it is important to determine the energy release rate of the flares from observational data as accurately as possible. To this end, we have carried out one dimensional hydrodynamic simulations of a flare loop and synthesized the count rate detected by the soft X-ray telescope (SXT) aboard Yohkoh satellite. We found that the time derivative of the thermal energy contents in a flare arcade derived from SXT data is smaller than the real energy release rate by a factor of 0.3 - 0.8, depending on the loop length and the energy release rate. The result of simulation is presented in the paper and used to calculate the reconnection rate. We found that reconnection rate is 0.047 for the X2.3 flare on 2000 November 24, 0.015 for the M3.7 flare on 2000 July 14, and 0.071 for the C8.9 flare on 2000 November 16. These values are similar to that derived from the direct observation of the reconnection inflow by Yokoyama et al. (2001), and consistent with the fast reconnection models such as that of Petschek (1964).

We report on the multi-epoch observations of H2O maser emission in star forming region OH 43.8-0.1 carried out with VLBI Exploration of Radio Astrometry (VERA). The large-scale maser distributions obtained by single-beam VLBI mapping reveal new maser spots scattered in area of 0.7 x 1.0 arcsec, in addition to a `shell-like' structure with a scale of 0.3 x 0.5 arcsec which was previously mapped by Downes et al.(1979). Proper motions are also obtained for 43 spots based on 5-epoch monitoring with a time span of 281 days. The distributions of proper motions show a systematic outflow in the north-south direction with an expansion velocity of ~8 km/s, and overall distributions of maser spots as well as proper motions are better represented by a bipolar flow plus a central maser cluster with a complex structure, rather than a shell with a uniform expansion such as those found in Cep A R5 and W75N VLA2. The distance to OH 43.8-0.1 is also estimated based on the statistical parallax, yielding D = 2.8 +/- 0.5 kpc. This distance is consistent with a near kinematic distance and rules out a far kinematics distance (~9 kpc), and the LSR velocity of OH 43.8-0.1 combined with the distance provides a constraint on the flatness of the galactic rotation curve, that there is no systematic difference in rotation speeds at the Sun and at the position of OH 43.8-0.1, which is located at the galacto-centric radius of ~6.3 kpc.

Space solar missions such as {\it Yohkoh} and {\it RHESSI} observe the hard X- and gamma-ray emission from energetic electrons in impulsive solar flares. Their energization mechanism, however, is unknown. In this paper, we suggest that the internal shocks are created in the reconnection jet and that they are possible sites of particle acceleration. We examine how magnetic reconnection creates the multiple shocks by performing two-dimensional resistive magnetohydrodynamic simulations. In this paper, we use a very small grid to resolve the diffusion region. As a result, we find that the current sheet becomes thin due to the tearing instability, and it collapses to a Sweet-Parker sheet. The thin sheet becomes unstable to the secondary tearing instability. Fast reconnection starts by the onset of anomalous resistivity immediately after the secondary tearing instability. During the bursty, time-dependent magnetic reconnection, the secondary tearing instability continues in the diffusion region where the anomalous resistivity is enhanced. As a result, many weak shocks are created in the reconnection jet. This situation produces turbulent reconnection. We suggest that multiple fast shocks are created in the jet and that the energetic electrons can be accelerated by these shocks.

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity $\sim 0.3c$) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units ($\tau_{\rm S} = r_{\rm S}/c$ where $r_{\rm S} \equiv 2GM/c^2$) which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around $r = 3 r_{\rm S}$. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfv\'{e}n waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star.

We analyzed the results of non-linear resistive magnetohydrodynamical (MHD) simulations of jet formation to study the acceleration mechanism of axisymmetric, resistive MHD jets. The initial state is a constant angular momentum, polytropic torus threaded by weak uniform vertical magnetic fields. The time evolution of the torus is simulated by applying the CIP-MOCCT scheme extended for resistive MHD equations. We carried out simulations up to 50 rotation period at the innermost radius of the disk created by accretion from the torus. The acceleration forces and the characteristics of resistive jets were studied by computing forces acting on Lagrangian test particles. Since the angle between the rotation axis of the disk and magnetic field lines is smaller in resistive models than in ideal MHD models, magnetocentrifugal acceleration is smaller. The effective potential along a magnetic field line has maximum around $z \sim 0.5r_0$ in resistive models, where $r_0$ is the radius where the density of the initial torus is maximum. Jets are launched after the disk material is lifted to this height by pressure gradient force. Even in this case, the main acceleration force around the slow magnetosonic point is the magnetocentrifugal force. The power of the resistive MHD jet is comparable to the mechanical energy liberated in the disk by mass accretion. Joule heating is not essential for the formation of jets.

We report on a strong outburst of the W49N H2O maser observed with VERA. Single-dish monitoring with VERA 20 m telescopes detected a strong outburst of the maser feature at V_LSR = -30.7 km/s in 2003 October. The outburst had a duration of ~100 days and a peak intensity of 7.9 x 10^4 Jy, being one of the strongest outbursts in W49N observed so far. VLBI observations with the VERA array were also carried out near to the maximum phase of the outburst, and the outburst spot was identified in the VLBI map. While the map was in good agreement with previous studies, showing three major concentrations of maser spots, we found a newly formed arc-like structure in the central maser concentration, which may be a shock front powered by a forming star or a star cluster. The outburst spot was found to be located on the arc-like structure, indicating a possible connection of the present outburst to a shock phenomenon.

The Kelvin--Helmholtz (KH) and tearing instabilities are likely to be important for the process of fast magnetic reconnection that is believed to explain the observed explosive energy release in solar flares. Theoretical studies of the instabilities, however, typically invoke simplified initial magnetic and velocity fields that are not solutions of the governing magnetohydrodynamic (MHD) equations. In the present study, the stability of a reconnecting current sheet is examined using a class of exact global MHD solutions for steady state incompressible magnetic reconnection, discovered by Craig & Henton. Numerical simulation indicates that the outflow solutions where the current sheet is formed by strong shearing flows are subject to the KH instability. The inflow solutions where the current sheet is formed by a fast and weakly sheared inflow are shown to be tearing unstable. Although the observed instability of the solutions can be interpreted qualitatively by applying standard linear results for the KH and tearing instabilities, the magnetic field and plasma flow, specified by the Craig--Henton solution, lead to the stabilization of the current sheet in some cases. The sensitivity of the instability growth rate to the global geometry of magnetic reconnection may help in solving the trigger problem in solar flare research.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of collapsars including a rotating black hole. This paper is an extension of our previous paper (Mizuno et al. 2004). The current calculation focuses on the effect of black hole rotation using general relativistic MHD with simplified microphysics. Initially, we assume that the core collapse is failed in this star. A few solar mass rotating black hole is inserted by hand into the calculation. We consider two cases, the co-rotating case and counter-rotating case with respect to the black hole rotation. Although the counter-rotating case may be unrealistic for collapsar, we perform as a maximally dragging case of magnetic field. The simulation results show the formation of a disk-like structure and the generation of a jet-like outflow near the central black hole. The jet-like outflow propagates outwardly with the twisted magnetic field and becomes collimated. We have found that the jets are generated and accelerated mainly by the magnetic field. The total jet velocity in the rotating black hole case is comparable to that of the non-rotating black hole case (Mizuno et al. 2004), 0.3c. When the rotation of the black hole is faster, the magnetic field is twisted strongly owing to the frame-dragging effect. The magnetic energy stored by the twisting magnetic field is converted to kinetic energy of the jet directly rather than propagating as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the poloidal velocity of the jet becomes faster. In the rapidly rotating black hole case the jet-like outflow can be produced by the frame dragging effect only through the twisting of magnetic field even if there is no stellar rotation.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of the gravitational collapse of a magnetized rotating massive star as a model of gamma ray bursts (GRBs). The current calculation focuses on general relativistic MHD with simplified microphysics (we ignore neutrino cooling, physical equation of state, and photodisintegration). Initially, we assume that the core collapse has failed in this star. A few $M_{\odot}$ black hole is inserted by hand into the calculation. The simulations presented in the paper follow the accretion of gas into a black hole that is assumed to have formed before the calculation begins.The simulation results show the formation of a disk-like structure and the generation of a jetlike outflow inside the shock wave launched at the core bounce. We have found that the jet is accelerated by the magnetic pressure and the centrifugal force and is collimated by the pinching force of the toroidal magnetic field amplified by the rotation and the effect of geometry of the poloidal magnetic field. The maximum velocity of the jet is mildly relativistic (0.3c). The velocity of the jet becomes larger as the initial rotational velocity of stellar matter gets faster. On the other hand, the dependence on the initial magnetic field strength is a bit more complicated: the velocity of the jet increases with the initial field strength in the weak field regime, then is saturated at some intermediate field strength, and decreases beyond the critical field strength. These results are related to the stored magnetic energy determined by the balance between the propagation time of the Alfven wave and the rotation time of the disk (or twisting time).

We investigated the statistical distributions of physical quantities of solar flares and associated coronal mass ejections (CMEs). We found that the distributions of the X-ray peak fluxes of CME-related flares, their time intervals, and speeds of associated CMEs are in good agreement with log-normal distributions. One possible interpretation of this is that only large energetic mass ejections can escape from the solar corona, which become CMEs in the interplanetary space. This ``filtering effect'' may be the reason why statistical distributions of some physical quantities are similar to log-normal distribution. It is known that the distributions of several physical quantities of gamma-ray bursts (GRBs) are also log-normal distributions, and that the time variations of gamma-ray intensity in GRBs are very similar to those of hard X-rays of solar flares. On the basis of these similarities and physical consideration of magnetohydrodynamic properties of an accretion disk around a black hole, which is supposed to be situated in the central engine of a GRB, we propose a new model of the central engine of GRBs, the {\it flare/CME model}, in which GRBs are formed by nonsteady jets consisting of intermittent ejections of mass (like solar CMEs) associated with reconnection events (like solar flares) in the accretion disk corona. Such nonsteady MHD jets would produce many shocks as well as local reconnection events far from the central engine. In this model, only large energetic mass ejections can escape from the accretion disk corona, so that statistical distributions of some physical quantities are similar to log-normal distributions.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of the gravitational collapse of a magnetized rotating massive star as a model of gamma ray bursts (GRBs). This simulation showed the formation of a disk-like structure and the generation of a jet-like outflow inside the shock wave launched at the core bounce. We have found the jet is accelerated by the magnetic pressure and the centrifugal force and is collimated by the pinching force of the toroidal magnetic field amplified by the rotation and the effect of geometry of the poloidal magnetic field. The maximum velocity of the jet is mildly relativistic ($\sim$ 0.3 c).

We present three-dimensional (3-D) magnetohydrodynamical (MHD) simulations of radiatively inefficient accretion flow around black holes. General relativistic effects are simulated by using the pseudo-Newtonian potential. We start calculations with a rotating torus threaded by localized poloidal magnetic fields with plasma beta, a ratio of the gas pressure to the magnetic pressure, $\beta =10$ and 100. When the bulk of torus material reaches the innermost region close to a central black hole, a magnetically driven jet emerges. This magnetic jet is derived by vertically inflating toroidal fields (`magnetic tower') and has a two-component structure: low-$\beta$ ($\lsim 1$) plasmas threaded with poloidal (vertical) fields are surrounded by that with toroidal fields. The collimation width of the jet depends on external pressure, pressure of ambient medium; the weaker the external pressure is, the wider and the longer-lasting becomes the jet. Unless the external pressure is negligible, the bipolar jet phase ceases after several dynamical timescales at the original torus position and a subsequent quasi-steady state starts. The black hole is surrounded by quasi-spherical zone with highly inhomogeneous structure in which toroidal fields are dominant except near the rotation axis. Mass accretion takes place mainly along the equatorial plane. Comparisons with other MHD simulation results and observational implications are discussed.

We present the results of first dual-beam observations with VERA (VLBI Exploration of Radio Astrometry). The first dual-beam observations of a pair of H2O maser sources W49N and OH43.8-0.1 have been carried out on 2002 May 29 and July 23, and fringes of H2O maser lines at 22 GHz have been successfully detected. While the residual fringe phases of both sources show rapid variations over 360 degree due to the atmospheric fluctuation, the differential phase between the two sources remains constant for 1 hour with r.m.s. of 8 degree, demonstrating that the atmospheric phase fluctuation is removed effectively by the dual-beam phase referencing. The analysis based on Allan standard deviation reveals that the differential phase is mostly dominated by white phase noise, and the coherence function calculated from the differential phase shows that after phase referencing the fringe visibility can be integrated for arbitrarily long time. These results demonstrate VERA's high capability of phase referencing, indicating that VERA is a promising tool for phase referencing VLBI astrometry at 10 microarcsec level accuracy.

Recent numerical simulations of MHD jets from accretion disks are briefly reviewed with emphasis on the scaling law for jet speed and the role of magnetic reconnection in relation to time variability in accretion disks, jets, and flares. On the basis of these studies, possible interpretation is given on why statistical properties of peak intensity, peak interval, and peak duration of gamma ray bursts (log-normal distribution) are different from those in solar flares and black hole accretion disks (power-law distribution). From these considerations, a new model, ``magnetized plasmoid model'', is proposed for a central engine of gamma ray bursts. (see http://www.kwasan.kyoto-u.ac.jp/~shibata/grb_final/grb.pdf for pdf files including both text and figures)

We report observations of water masers around the semiregular variable RT Virginis (RT Vir), which have been made with the Very Long Baseline Array (VLBA) of the National Radio Astronomy Observatory (NRAO) at five epochs, each separated by three weeks of time. We detected about 60 maser features at each epoch. Overall, 61 features, detected at least twice, were tracked by their radial velocities and proper motions. The 3-D maser kinematics exhibited a circumstellar envelope that is expanding roughly spherically with a velocity of about 8 km/s. Asymmetries in both the spatial and velocity distributions of the maser features were found in the envelope, but less significant than that found in other semiregular variables. Systematic radial-velocity drifts of individual maser features were found with amplitudes of <= 2 km/s/yr. For one maser feature, we found a quadratic position shift with time along a straight line on the sky. This apparent motion indicates an acceleration with an amplitude of 33 km/s/yr, implying the passage of a shock wave driven by the stellar pulsation of RT Vir. The acceleration motion is likely seen only on the sky plane because of a large velocity gradient formed in the accelerating maser region. We estimated the distance to RT Vir to be about 220 pc on the basis of both the statistical parallax and model-fitting methods for the maser kinematics.

We present the relation of the spatial distribution of H-alpha kernels with the distribution of hard X-ray (HXR) sources seen during the 2001 April 10 solar flare. This flare was observed in H-alpha with the (Sartorius) telescope at Kwasan Observatory, Kyoto University, and in hard X-rays (HXRs) with the Hard X-ray Telescope (HXT) onboard Yohkoh. We compared the spatial distribution of the HXR sources with that of the H-alpha kernels. While many H-alpha kernels are found to brighten successively during the evolution of the flare ribbons, only a few radiation sources are seen in the HXR images. We measured the photospheric magnetic field strengths at each radiation source in the H-alpha images, and found that the H-alpha kernels accompanied by HXR radiation have magnetic strengths about 3 times larger than those without HXR radiation. We also estimated the energy release rates based on the magnetic reconnection model. The release rates at the H-alpha kernels with accompanying HXR sources are 16-27 times larger than those without HXR sources. These values are sufficiently larger than the dynamic range of HXT, which is about 10, so that the difference between the spatial distributions of the H$\alpha$ kernels and the HXR sources can be explained.

We propose the Galactic flare model for the origin of the X-ray gas in the Galactic halo. For this purpose, we examine the magnetic reconnection triggered by Parker instability (magnetic buoyancy instability), by performing the two-dimensional resistive numerical magnetohydrodynamic simulations. As a result of numerical simulations, the system evolves as following phases: Parker instability occurs in the Galactic disk. In the nonlinear phase of Parker instability, the magnetic loop inflates from the Galactic disk into the Galactic halo, and collides with the anti-parallel magnetic field, so that the current sheets are created in the Galactic halo. The tearing instability occurs, and creates the plasmoids (magnetic islands). Just after the plasmoid ejection, further current-sheet thinning occurs in the sheet, and the anomalous resistivity sets in. Petschek reconnection starts, and heats the gas quickly in the Galactic halo. It also creates the slow and fast shock regions in the Galactic halo. The magnetic field ($B\sim 3 \mu$G), for example, can heat the gas ($n\sim 10^{-3}$ cm$^{-3}$) to temperature of $\sim 10^6$ K via the reconnection in the Galactic halo. The gas is accelerated to Alfv\'en velocity ($\sim 300$ km s$^{-1}$). Such high velocity jets are the evidence of the Galactic flare model we present in this paper, if the Doppler shift of the bipolar jet is detected in the Galactic halo. Full size figures are available at http://www.kwasan.kyoto-u.ac.jp/~tanuma/study/ApJ2002/ApJ2002.html

We have performed a first fully 3-D GRMHD simulation with Schwarzschild black hole with a free falling corona. The initial simulation results show that a jet is created as in the previous simulations using the axisymmetric geometry with the mirror symmetry at the equator. However, the time to generate the jet is longer than in the 2-D axisymmetric simulations. We expect that due to the additional freedom in the azimuthal dimension without axisymmetry with respect to the z axis and reflection symmetry with respect to the equatorial plane, the dynamics of jet formation can be modified. Further simulations are required for the study of instabilities along the azimuthal direction such as accretion-eject instability

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.

We construct a simple model for a magnetic reconnection heated corona above a thin accretion disk in AGNs and Galactic black hole candidates(GBHCs). The basic assumptions are that (1) the magnetic reconnection heat is cooled down overwhelmingly by Compton scattering in the corona and that (2) thermal conduction is dominantly cooled by evaporation of the chromospheric plasma in the disk-corona interface before Compton cooling sets in. With these two basic equations as well as equi-partition of magnetic energy with gas energy in the disk, we can consistently determine the fraction of accretion energy dissipated in the corona without free parameters, and thus determine the temperature and all other quantities in both the corona and disk for given black hole mass and accretion rate. Then, we calculate the luminosity contributed from the disk and corona and the coronal flux weighted Compton $y$ parameter. It is found that, at a low luminosity (less than 0.1 Eddington luminosity) the spectrum is hard with energy spectral index of around 1, while at a high luminosity (larger than 0.1 Eddington luminosity) the spectrum can be either soft or hard. We also find that the situation is almost same for supermassive and stellar-mass black holes. These features are consistent with observations of AGNs and GBHCs.

We have performed a first fully 3-D GRMHD simulation with Schwarzschild black hole with a free falling corona. The initial simulation results show that a jet is created as in previous axisymmetric simulations. However, the time to generate the jet is longer than in the 2-D simulations. We expect that due to the additional azimuthal dimension the dynamics of jet formation can be modified.

We report on the first space-VLBI observations of the OH masers in two main-line OH transitions at 1665 and 1667 MHz. The observations involved the space radio telescope on board the Japanese satellite HALCA and an array of ground radio telescopes. The maps of the maser region and images of individual maser spots were produced with an angular resolution of 1 mas, which is several times higher than the angular resolution available on the ground. The maser spots were only partly resolved and a lower limit to the brightness temperature 6x10^12 K was obtained. The masers seem to be located in the direction of low interstellar scattering.

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.

We present a new model for time evolution of fast magnetic reconnection in free space, which is characterized by self-similarity. Reconnection triggered by locally enhanced resistivity assumed at the center of the current sheet can self-similarly and unlimitedly evolve until external factors affect the evolution. The possibility and stability of this type of evolution are verified by numerical simulations in a very wide spatial dynamic range. Actual astrophysical reconnection in solar flares and geomagnetospheric substorms can be treated as an evolutionary process in free space, because the resultant scale is much larger than the initial scale. In spite of this fact, most of the previous numerical works focused on the evolutionary characters strongly affected by artificial boundary conditions on the simulation boundary. Our new model clarifies a realistic evolution for such cases. The characteristic structure around the diffusion region is quite similar to the Petschek model which is characterized by a pair of slow-mode shocks and the fast-mode rarefaction-dominated inflow. However, in the outer region, a vortex-like return flow driven by the fast-mode compression caused by the piston effect of the plasmoid takes place. The entire reconnection system expands self-similarly.

We carried out 2.5-dimensional resistive magnetohydrodynamic simulations to study the effects of magnetic diffusivity on magnetically driven mass accretion and jet formation. We found that (1) when the normalized magnetic diffusivity, is small, mass accretion and jet formation take place intermittently; (2) when diffusivity is middle, the system evolves toward a quasi-steady state; and the system evolves toward a quasi-steady state; and (3) when diffusivity is large, the accretion/mass outflow rate decreases with diffusivity and approaches 0. The results of these simulations indicate magnetic braking provide a mass accretion rate which is sufficient to explain the activity of AGNs.

We examine the magnetic reconnection triggered by a supernova (or a point explosion) in interstellar medium, by performing two-dimensional resistive magnetohydrodynamic (MHD) numerical simulations with high spatial resolution. We found that the magnetic reconnection starts long after a supernova shock (fast-mode MHD shock) passes a current sheet. The current sheet evolves as follows: (i) Tearing-mode instability is excited by the supernova shock, and the current sheet becomes thin in its nonlinear stage. (ii) The current-sheet thinning is saturated when the current-sheet thickness becomes comparable to that of Sweet-Parker current sheet. After that, Sweet-Parker type reconnection starts, and the current-sheet length increases. (iii) ``Secondary tearing-mode instability'' occurs in the thin Sweet-Parker current sheet. (iv) As a result, further current-sheet thinning occurs and anomalous resistivity sets in, because gas density decreases in the current sheet. Petschek type reconnection starts and heats interstellar gas. Magnetic energy is released quickly while magnetic islands are moving in the current sheet during Petschek type reconnection. The released magnetic energy is determined by the interstellar magnetic field strength, not energy of initial explosion nor distance to explosion. We suggest that magnetic reconnection is a possible mechanism to generate X-ray gas in Galaxy.

We report on the first space-VLBI observations of the OH34.26+0.15 maser in two main line OH transitions at 1665 and 1667 MHz. The observations involved the space radiotelescope on board the Japanese satellite HALCA and an array of ground radio telescopes. The map of the maser region and images of individual maser spots were produced with an angular resolution of 1 milliarcsec which is several times higher than the angular resolution available on the ground. The maser spots were only partly resolved and a lower limit to the brightness temperature 6x10^{12} K was obtained. The maser seems to be located in the direction of low interstellar scattering, an order of magnitude lower than the scattering of a nearby extragalactic source and pulsar.

Rapid fluctuation with a frequency dependence of $1/f^{\alpha}$ (with $\alpha \simeq 1 - 2$) is characteristic of radiation from black-hole objects. Its origin remains poorly understood. We examine the three-dimensional magnetohydrodynamical (MHD) simulation data, finding that a magnetized accretion disk exhibits both $1/f^\alpha$ fluctuation (with $\alpha \simeq 2$) and a fractal magnetic structure (with the fractal dimension of $D \sim 1.9$). The fractal field configuration leads reconnection events with a variety of released energy and of duration, thereby producing $1/f^\alpha$ fluctuations.

Recent observations of Galactic Black Hole Candidates (BHCs) suggest that those that are superluminal jet sources have more rapid black hole spin rates than otherwise normal BHCs. This provides observational support for models of relativistic jet formation that extract rotational energy of the central black hole. To investigate this mechanism, we have developed a new general relativistic magnetohydrodynamic code in Kerr geometry. Here we report on the first numerical simulation of the formation of a relativistic jet in a rapidly-rotating (a=0.95) Kerr black hole magnetosphere. We assume that the initial velocity of the disk is zero. We find that the maximum velocity of the jet reaches 0.93c (Lorentz factor, 2.7) and the terminal velocity of the jet is 0.85c (Lorentz factor, 1.9). On the other hand, for a non-rotating (a=0) Schwarzschild black hole, the maximum outflow velocity is less than 0.6c for initial magnetospheric conditions similar to to those of the Kerr black hole case. These numerical results show the importance of the rapidly rotating black hole for the relativistic jet formation.

To investigate the formation mechanism of relativistic jets in active galactic nuclei and micro-quasars, we have developed a new general relativistic magnetohydrodynamic code in Kerr geometry. Here we report on the first numerical simulation of jet formation in a rapidly-rotating (a=0.95) Kerr black hole magnetosphere. We study cases in which the Keplerian accretion disk is both co-rotating and counter-rotating with respect to the black hole rotation. In the co-rotating disk case, our results are almost the same as those in Schwarzschild black hole cases: a gas pressure-driven jet is formed by a shock in the disk, and a weaker magnetically-driven jet is also generated outside the gas pressure-driven jet. On the other hand, in the counter-rotating disk case, a new powerful magnetically-driven jet is formed inside the gas pressure-driven jet. The newly found magnetically-driven jet in the latter case is accelerated by a strong magnetic field created by frame dragging in the ergosphere. Through this process, the magnetic field extracts the energy of the black hole rotation.

We propose a model of hard X-ray flares in protostars observed by ASCA satellite. Assuming that the dipole magnetic field of the protostar threads the protostellar disk, we carried out 2.5-dimensional magnetohydrodynamic (MHD) simulations of the disk-star interaction. The closed magnetic loops connecting the central star and the disk are twisted by the rotation of the disk. As the twist accumulates, magnetic loops expand and finally approach to the open field configuration. A current sheet is formed inside the expanding loops. In the presence of resistivity, magnetic reconnection takes place in the current sheet. Outgoing magnetic island and post flare loops are formed as a result of the reconnection. The time scale of this `flare' is the order of the rotation period of the disk. The released magnetic energy partly goes into the thermal energy and heats up the flaring plasma up to $10^8$ K. The length of the flaring loop is several times of the radius of the central star, consistent with observations. The speed of the hot plasmoid ejected by the reconnection is $200-400 $ km s$^{-1}$ when the footpoint of the loop is at 0.03 AU from 1 M$_\odot$ protostar. The hot plasma outflow can explain the speed and mass flow rate of optical jets. Dense, cold, magnetically accelerated wind ($v \sim 150-250$ km s$^{-1}$) emanates from the surface of the disk along the partially open magnetic field lines threading the disk. This dense, cold wind may correspond to high velocity neutral winds.

The effect of the magnetic skew on the Parker instability is investigated by means of the linear stability analysis for a gravitationally stratified gas layer permeated by a horizontal magnetic field. When the magnetic field is skewed (i.e., the field line direction is a function of the height), the wavelength of the most unstable mode is $ \lambda \; \sim \; 10 H $ where $ H $ is the pressure scale height. The growth rate of the short wavelength modes is greatly reduced when the gradient in the magnetic field direction exceeds 0.5 radian per scale height. Our results indicate that the Parker instability in a skewed magnetic field preferentially forms large scale structures like giant molecular clouds.