The pebble accretion model has the potential to explain the formation of various types of planets. A growing planet embedded in a disk induces three-dimensional (3D) gas flow, which may influence pebble accretion. In this study, we investigate the influence of the 3D planet-induced gas flow on pebble accretion. Assuming a non-isothermal, inviscid gas disk, we perform 3D hydrodynamical simulations on the spherical polar grid. Then we numerically integrate the equation of motion of pebbles in 3D using hydrodynamical simulations data. We find that the trajectories of pebbles in the planet-induced gas flow differ significantly from those in the unperturbed shear flow for a wide range of pebble sizes investigated (${\rm St}=10^{-3}$--$10^{0}$, where ${\rm St}$ is the Stokes number). The horseshoe flow and outflow of the gas alter the motion of the pebbles, which leads to the reduction of the width of the accretion window, $w_{\rm acc}$, and the accretion cross section, $A_{\rm acc}$. On the other hand, the changes in trajectories also cause an increase in relative velocity of pebbles to the planet, which offsets the reduction of $w_{\rm acc}$ and $A_{\rm acc}$. As a consequence, in the Stokes regime, the accretion probability of pebbles, $P_{\rm acc}$, in the planet-induced gas flow is comparable to that in the unperturbed shear flow except when the Stokes number is small, ${\rm St}\sim10^{-3}$, in 2D accretion, or when the thermal mass of the planet is small, $m=0.03$ in 3D accretion. In contrast, in the Epstein regime, $P_{\rm acc}$ in the planet-induced gas flow becomes smaller than that in the shear flow in the Stokes regime in both 2D and 3D accretion, regardless of assumed ${\rm St}$ and $m$. Our results suggest that the 3D planet-induced gas flow may be helpful to explain the distribution of exoplanets as well as the architecture of the solar system.

The three-dimensional structure of the gas flow around a planet is thought to influence the accretion of both gas and solid materials. In particular, the outflow in the mid-plane region may prevent the accretion of the solid materials and delay the formation of super-Earths' cores. However, it is not yet understood how the nature of the flow field and outflow speed change as a function of the planetary mass. In this study, we investigate the dependence of gas flow around a planet embedded in a protoplanetary disc on the planetary mass. Assuming an isothermal, inviscid gas disc, we perform three-dimensional hydrodynamical simulations on the spherical polar grid, which has a planet located at its centre. We find that gas enters the Bondi or Hill sphere at high latitudes and exits through the mid-plane region of the disc regardless of the assumed dimensionless planetary mass $m=R_{\rm Bondi}/H$, where $R_{\rm Bondi}$ and $H$ are the Bondi radius of the planet and disc scale height, respectively. The altitude from where gas predominantly enters the envelope varies with the planetary mass. The outflow speed can be expressed as $|u_{\rm out}|=\sqrt{3/2}mc_{\rm s}$ $(R_{\rm Bondi}\leq R_{\rm Hill})$ or $|u_{\rm out}|=\sqrt{3/2}(m/3)^{1/3} c_{\rm s}$ ($R_{\rm Bondi}\geq R_{\rm Hill}$), where $c_{\rm s}$ is the isothermal sound speed and $R_{\rm Hill}$ is the Hill radius. The outflow around a planet may reduce the accretion of dust and pebbles onto the planet when $m\gtrsim\sqrt{\rm St}$, where St is the Stokes number. Our results suggest that the flow around proto-cores of super-Earths may delay their growth and, consequently, help them to avoid runaway gas accretion within the lifetime of the gas disc.

The hydrogen isotopic (D/H) ratio reflects the global cycling and evolution of water on Earth as it fractionates through planetary processes. We model the water cycle taking seafloor hydrothermal alteration, chemical alteration of continental crust, slab subduction, hydrogen escape from the early Earth, and degassing at mid-ocean ridges, hot spots, and arcs into account. The differences in D/H ratios between present-day oceans, oceanic and continental crust, and mantle are thought to reflect isotopic fractionation through seafloor alteration, chemical alteration, and slab dehydration. However, if the speed of plate tectonics has been nearly constant through out Earth's history, the degassing and regassing rates are too small to reach the present-day D/H ratios. We show that (a) hydrogen escape from reduced early atmosphere, (b) secular net regassing, or (c) faster plate tectonics on early Earth is needed to reproduce the present-day D/H ratios of the water reservoirs. The low D/H ratio of Archean seawater at 3.8 Ga has previously been interpreted as a signature of (a) hydrogen escape, but we find it can also be explained either by (b) secular net regassing or by (c) faster plate tectonics on early Earth. The rates of hydrogen escape from early Earth and secular regassing on present-day Earth are constrained to be lower than 2.1 x 10^11 kg/yr and 3.9 x 10^11 kg/yr. Consequently, the volume of water in the present-day mantle could result entirely from the regassing through Earth's history. In that case, the volume of initial oceans could be 2 to 3 times larger than that of current Earth. We suggest that, in addition to the D/H ratio of Archean seawater, identifying the D/H ratios of both seawater and mantle throughout Earth's history would allow to distinguish these evolutionary scenarios.

The ubiquity of super-Earths poses a problem for planet formation theory to explain how they avoided becoming gas giants. Rapid recycling of the envelope gas of planets embedded in a protoplanetary disc has been proposed to delay the cooling and following accretion of disc gas. We compare isothermal and non-isothermal 3D hydrodynamical simulations of the gas flow past a planet to investigate the influence on the feasibility of the recycling mechanism. Radiative cooling is implemented by using the beta cooling model. We find that, in either case, gas enters the Bondi sphere at high latitudes and leaves through the midplane regions, or vice versa when disc gas rotates sub-Keplerian. However, in contrast to the isothermal case where the recycling flow reaches the deeper part of the envelope, the inflow is inhibited from reaching the deep envelope in the non-isothermal case. Once the atmosphere starts cooling, buoyant force prevents the high-entropy disc gas from intruding the low-entropy atmosphere. We suggest that the buoyancy barrier isolates the lower envelope from the recycling and allows further cooling, which may lead runaway gas accretion onto the core.

The atmospheres of the terrestrial planets are known to have been modified as a consequence of the impact degassing and atmospheric erosion during the late accretion. Despite the commonality of these processes, there are distinct gaps -- roughly two orders of magnitude -- between the abundances of noble gases and nitrogen in the present-day atmospheres on Venus, Earth, and Mars. The element partitioning on planetary surfaces is thought to be significantly different between the three planets ~4 Ga: the runaway greenhouse on Venus, the carbon-silicate cycle and ocean formation on Earth, and the CO2-ice and H2O-ice formation on Mars. Consequences of element partitioning for the atmospheric evolution during the late accretion onto Venus, Earth, and Mars are investigated with a numerical model. We set upper limits to the partial pressures of CO2 and H2O on Earth and Mars, which corresponds to the state of phase equilibrium and carbon-silicate cycle. The final N2 mass shrinks by ~40% and ~15% for Earth and Mars, respectively. The effect of element partitioning is found to be insufficient to reproduce the gaps. For Venus, the survival of the primordial atmosphere through the late accretion may partially account for the present-day atmosphere. Whereas on Mars, the atmospheric escape due to solar extreme UV and wind may have also influenced the atmospheric evolution.

We examine the history of the loss and replenishment of the Martian atmosphere using elemental and isotopic compositions of nitrogen and noble gases. The evolution of the atmosphere is calculated by taking into consideration various processes: impact erosion and replenishment by asteroids and comets, atmospheric escape induced by solar radiation and wind, volcanic degassing, and gas deposition by interplanetary dust particles. Our model reproduces the elemental and isotopic compositions of N and noble gases (except for Xe) in the Martian atmosphere, as inferred from exploration missions and analyses of Martian meteorites. Other processes such as ionization-induced fractionation, which are not included in our model, are likely to make a large contribution in producing the current Xe isotope composition. Since intense impacts during the heavy bombardment period greatly affect the atmospheric mass, the atmospheric pressure evolves stochastically. Whereas a dense atmosphere preserves primitive isotopic compositions, a thin atmosphere on early Mars is severely influenced by stochastic impact events and following escape-induced fractionation. The onset of fractionation following the decrease in atmospheric pressure is explained by shorter timescales of isotopic fractionation under a lower atmospheric pressure. The comparison of our numerical results with the less fractionated N ($^15$N/$^14$N) and Ar ($^38$Ar/$^36$Ar) isotope compositions of the ancient atmosphere recorded in the Martian meteorite Allan Hills 84001 provides a lower limit of the atmospheric pressure in 4 Ga to preserve the primitive isotopic compositions. We conclude that the atmospheric pressure was higher than approximately 0.5 bar at 4 Ga.

Observations have revealed that a significant number of hot Jupiters have anomalously large radii. Layered convection induced by compositional inhomogeneity has been proposed to account for the radius anomaly of hot Jupiters. To reexamine the impact of the compositional inhomogeneity, we perform an evolutionary calculation by determining convection regime at each evolutionary time step according to the criteria from linear analyses. It is shown that the impact is limited in the case of the monotonic gradient of heavy element abundance. The layered convection is absent for the first 1 Gyr from the formation of hot Jupiters and instead overturning convection develops. The super-adiabaticity of the temperature gradient is limited by the neutrally stable state for the Ledoux stability criterion. The effect of the increased mass of heavy elements essentially compensates the effect of the delayed contraction on the planetary radius caused by compositional inhomogeneity. In addition, even in the case where the layered convection is artificially imposed, this mechanism requires extremely thin layers (~ 10^1-10^3 cm) to account for the observed radius anomaly. The long-term stability of such thin layers remains to be studied. Therefore, if the criteria adopted in this paper are adequate, it might be difficult to explain the inflated radii of hot Jupiters by monotonic gradient of heavy element abundance alone.

Remote sensing data from orbiter missions have proposed that ground ice may currently exist on Mars, although the volume is still uncertain. Recent analyses of Martian meteorites have suggested that the water reservoirs have at least three distinct hydrogen isotope compositions (D/H ratios): primordial and high D/H ratios, which are approximately the same and six times that of ocean water on Earth, respectively, and a newly identified intermediate D/H ratio, which is approximately two to three times higher than that in ocean water on Earth. We calculate the evolution of the D/H ratios and the volumes of the water reservoirs on Mars by modeling the exchange of hydrogen isotopes between multiple water reservoirs and the atmospheric escape. The D/H ratio is slightly higher in the topmost thin surface-ice layer than that in the atmosphere because of isotopic fractionation by sublimation, whereas the water-ice reservoir just below the exchangeable topmost surface layer retains the intermediate D/H signature found in Martian meteorites. We propose two possible models for constraining the volume of the ground ice considering the observed D/H ratios and geomorphological estimates of Paleo-oceans. The first assumes that the atmospheric loss is dominated by the Jeans escape. In this case, the volume of ground ice should be larger than the total volume of the observable surface ice that mainly occurs as polar layered ice deposits. The other model assumes diffusion-limited atmospheric loss in which the interactive evolution of the multiple water reservoirs naturally accounts for the observed D/H ratios. In this scenario, a large volume of ground ice does not necessarily exist currently on Mars as opposed to the perspective view proposed on the basis of recent orbiter missions.

Martian surface morphology implies that Mars was once warm enough to maintain persistent liquid water on its surface. While the high D/H ratios (~6 times the Earth's ocean water) of the current martian atmosphere suggest that significant water has been lost from the surface during martian history, the timing, processes, and the amount of the water loss have been poorly constrained. Recent technical developments of ion-microprobe analysis of martian meteorites have provided accurate estimation of hydrogen isotope compositions (D/H) of martian water reservoirs at the time when the meteorites formed. Based on the D/H data from the meteorites, this study demonstrates that the water loss during the pre-Noachian (>41-99 m global equivalent layers, GEL) was more significant than in the rest of martian history (>10-53 m GEL). Combining our results with geological and geomorphological evidence for ancient oceans, we propose that undetected subsurface water/ice (~100-1000 m GEL) should have existed, and it exceeds the observable present water inventory (~20-30 m GEL) on Mars.

During their evolution, short-period exoplanets may lose envelope mass through atmospheric escape owing to intense XUV (X-ray and extreme ultraviolet) radiation from their host stars. Roche-lobe overflow induced by orbital evolution or intense atmospheric escape can also contribute to mass loss. To study the effects of mass loss on inner planet populations, we calculate the evolution of hot Jupiters considering mass loss of their envelopes and thermal contraction. Mass loss is assumed to occur through XUV-driven atmospheric escape and the following Roche-lobe overflow. The runaway effect of mass loss results in a dichotomy of populations: hot Jupiters that retain their envelopes and super Earths whose envelopes are completely lost. Evolution primarily depends on the core masses of planets and only slightly on migration history. In hot Jupiters with small cores (~10 Earth masses), runaway atmospheric escape followed by Roche-lobe overflow may create sub-Jupiter deserts, as observed in both mass and radius distributions of planetary populations. Comparing our results with formation scenarios and observed exoplanets populations, we propose that populations of closely orbiting exoplanets are formed by capturing of planets at/inside the inner edges of protoplanetary disks and subsequent evaporation of sub-Jupiters.

Short-period exoplanets potentially lose envelope masses during their evolution because of atmospheric escape caused by the intense XUV radiation from their host stars. We develop a combined model of atmospheric mass loss calculation and thermal evolution calculation of a planet to simulate its evolution and explore the dependences on the formation history of the planet. Thermal atmospheric escape as well as the Roche-lobe overflow contributes to mass loss. The maximum initial planetary model mass depends primarily on the assumed evolution model of the stellar XUV luminosity. We adapt the model to CoRoT-7b and Kepler-10b to explore the evolution of both planets and the maximum initial mass of these planets. We take the recent X-ray observation of CoRoT-7 into account and exploring the effect of different XUV evolution models on the planetary initial mass. Our calculations indicate that both hot super Earths could be remnants of Jupiter mass gas planets.

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.

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 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 present conspicuous activities of plasma ejections along a light bridge of a stable and mature sunspot in NOAA 8971 on 2000 May 02. We found the ejections both in the H-alpha (10^4 K) images obtained with the Domeless Solar Telescope (DST) at Hida Observatory, and in the 171 A (Fe IX/X, ~10^6 K) images obtained with the Transition Region and Coronal Explorer (TRACE). Main characteristics of the ejections are as follows: (1)Ejections occur intermittently and recurrently, (2)The velocities and the timings of the 171 A-ejections are the same as those of H-alpha ejections, (3)The appearance of ejections are different from one another, i.e. the H-alpha ejections have jet-like appearance, while that of 171 A-ejections is like a loop.

Sunspot proper motions and flares of a super active region NOAA 5395, which was the biggest and the most flare-active region in the 22nd sunspot cycle, were analyzed in details. We measured sunspot proper motions by using the H-alpha - 5.0 A images obtained with the 60-cm Domeless Solar Telescope (DST) at Hida Observatory, Kyoto University and found some peculiar vortex-like motions of small satellite spots successively emerged from the leading edge of this sunspot group. To explain these motions of small sunspots, we proposed a schematic model of the successive emergence of twisted and winding magnetic flux ropes coiling around a trunk of magnetic flux tube. The location of the strongest flare activity was found to coincide with very the site of the vortex-like motions of sunspots. We conclude that the flare-productive magnetic shear is produced by the emergence of the twisted magnetic flux bundle. Magnetic energy is stored in the twisted flux tubes which are originally formed in the convection zone and released as flares in the course of the emergence of the twisted flux tubes above the photosphere.