Search

Carbonaceous chondrites (CCs) are windows into the early Solar system and the histories of their parent bodies. Their infrared spectral signatures are powerful proxies for deciphering their composition and evolution history, but still present formidable challenges. In our study, we delved into the infrared spectra spanning 1-25 micron of 17 CCs, with distinct petrological characteristics and varying degrees of alteration. As aqueous alteration intensifies, the 3 micron-region absorption feature associated with OH-bearing minerals and water, and the 6 micron band indicative of water molecules, both grow in intensity. Simultaneously, their band centers shift towards shorter wavelengths. Moreover, as alteration progresses, a distinctive absorption feature emerges near 2.72 micron, resembling the OH absorption feature found in serpentine and saponite minerals. Comparison of aqueous alteration to laboratory-heated CCs suggests that the 3 micron region OH/H2O absorption feature differs between CC heated to less than or more than ~300C. The 12.4 micron/11.4 micron reflectance ratio diminishes, and the reflectance peak in the 9-14 micron range shifts towards shorter wavelengths. These changes are attributed to the transformation of anhydrous silicates into phyllosilicates. In the 15-25 micron region, the influence of thermal metamorphism becomes evident and results in the appearance of more spectral features, the single reflectance peak at 22.1 micron undergoes a transformation into two distinct peaks at 19 micron and 25 micron, which is primarily attributed to the increased presence of anhydrous silicates and olivine recrystallization. These findings offer novel insights into the volatile-rich compositions of C-complex asteroids and the thermal evolution histories of their parent bodies., Comment: Accepted by Icarus, 7028 total words, 11 figures, 1 table

The early solar system contained a short-lived radionuclide, $^{26}$Al (its half-life time $t_{1/2} = 0.7$ Myr). The decay energy $^{26}$Al is thought to have controlled the thermal evolution of planetesimals and, possibly, the water contents of planets. Many hypotheses have been proposed for the origin of $^{26}$Al in the solar system. One of the possible hypotheses is the `disk injection scenario'; when the protoplanetary disk of the solar system had already formed, a nearby $(<1 \,\mathrm{pc})$ supernova injected radioactive material directly into the disk. Such a $^{26}$Al injection hypothesis has been tested so far with limited setups for disk structure and supernova distance, and treated disk disruption and $^{26}$Al injection separately. Here, we revisit this problem to investigate whether there are self-consistent conditions under which the surviving disk radius can receive enough $^{26}$Al which can account for the abundance in the early solar system. We also consider a range of disk mass and structure, $^{26}$Al yields from supernova, and a large dust mass fraction $\eta_\mathrm{d}$. We find that $^{26}$Al yields of supernova are required as $\gtrsim 2.1\times10^{-3}M_\odot(\eta_\mathrm{d}/0.2)^{-1}$, challenging to achieve with known possible $^{26}$Al ejection and dust mass fraction ranges. Furthermore, we find that even if the above conditions are met, the supernova flow changes the disk temperature, which may not be consistent with the solar-system record. Our results place a strong constraint on the disk injection scenario. Rather, we suggest that the fresh $^{26}$Al of the early solar system must have been synthesized/injected in other ways., Comment: 12 pages, 5 figures, 2tables, published in The Astrophysical Journal

A growing planet embedded in a protoplanetary disk induces three-dimensional gas flow, which exhibits a midplane outflow that can suppress dust accretion onto the planet and form global dust substructures (rings and gaps). Because analytic formulae for the planet-induced outflow are useful for modeling its influences on local and global dust surface densities and planet accretion, we derive the analytic formulae that describe the morphology and velocity of the planet-induced outflow. We first perform three-dimensional, nonisothermal hydrodynamical simulations of the gas flow past a planet, which enables us to introduce a fitting formula describing the morphology of the outflow. We then derive an analytic formula for the outflow speed using Bernoulli's theorem. We successfully derived a fitting formula for the midplane outflow morphology (the shape of the streamline), which is valid when the dimensionless thermal mass falls below $m\lesssim0.6$. The obtained analytic formulae for the outflow, such as the maximum outflow speed and the velocity distributions of the outflow in the radial and vertical directions to the disk, show good agreement with the numerical results. We find the following trends: (1) the maximum outflow speed increases with the planetary mass and has a peak of $\sim$30--40$\%$ of the sound speed when the dimensionless thermal mass is $m\sim0.3$, corresponding to a super-Earth mass planet at 1 au for the typical steady accretion disk model, and (2) the presence of the headwind (namely, the global pressure force acting in the positive radial direction of the disk) enhances (reduces) the outflow toward the outside (inside) of the planetary orbit. The planet-induced outflow of the gas affects the dust motion when the dimensionless stopping time of dust falls below ${\rm St}\lesssim\min(10m^2,0.1)$., Comment: 25 pages, 29 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

GREX-PLUS (Galaxy Reionization EXplorer and PLanetary Universe Spectrometer) is a mission candidate for a JAXA's strategic L-class mission to be launched in the 2030s. Its primary sciences are two-fold: galaxy formation and evolution and planetary system formation and evolution. The GREX-PLUS spacecraft will carry a 1.2 m primary mirror aperture telescope cooled down to 50 K. The two science instruments will be onboard: a wide-field camera in the 2-8 $\mu$m wavelength band and a high resolution spectrometer with a wavelength resolution of 30,000 in the 10-18 $\mu$m band. The GREX-PLUS wide-field camera aims to detect the first generation of galaxies at redshift $z>15$. The GREX-PLUS high resolution spectrometer aims to identify the location of the water ``snow line'' in proto-planetary disks. Both instruments will provide unique data sets for a broad range of scientific topics including galaxy mass assembly, origin of supermassive blackholes, infrared background radiation, molecular spectroscopy in the interstellar medium, transit spectroscopy for exoplanet atmosphere, planetary atmosphere in the Solar system, and so on., Comment: This document is the first version of a collection of scientific themes which can be achieved with GREX-PLUS. Each section in Chapters 2 and 3 is based on the presentation at the GREX-PLUS Science Workshop held on 24-25 March, 2022 at Waseda University

We investigate the spin state of a protoplanet during the pebble accretion influenced by the gas flow in the gravitational potential of the protoplanet and how it depends on the planetary mass, the headwind speed, the distance from the host star, and the pebble size. We perform nonisothermal three-dimensional hydrodynamical simulations in a local frame to obtain the gas flow around the planet. We then numerically integrate three-dimensional orbits of pebbles under the obtained gas flow. Finally, assuming uniform spatial distribution of incoming pebbles, we calculate net spin by summing up specific angular momentum that individual pebbles transfer to the protoplanet at impacts. We find that a protoplanet with the envelope acquires prograde net spin rotation regardless of the planetary mass, the pebble size, and the headwind speed of the gas. This is because accreting pebbles are dragged by the envelope that commonly has prograde rotation. As the planetary mass or orbital radius increases, the envelope is thicker and the prograde rotation is faster, resulting in faster net prograde spin. When the dimensionless thermal mass of the planet, $m = R_{\mathrm{Bondi}} / H$, where $R_{\mathrm{Bondi}}$ and $H$ are the Bondi radius and the disk gas scale height, is larger than a certain critical mass ($m \gtrsim 0.3$ at $0.1 \, \mathrm{au}$ or $m \gtrsim 0.1$ at $1 \, \mathrm{au}$), the spin rotation exceeds the breakup one. The predicted spin frequency reaches the breakup one at the planetary mass $m_{\mathrm{iso,rot}} \sim 0.1 \, (a / 1 \, \mathrm{au})^{-1/2}$ (where $a$ is the orbital radius), suggesting that the protoplanet cannot grow beyond $m_{\mathrm{iso,rot}}$. It is consistent with the Earth's current mass and could help the formation of the Moon by a giant impact on fast-spinning proto-Earth., Comment: 14 pages, 10 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

Recent high-spatial-resolution observations have revealed dust substructures in protoplanetary disks such as rings and gaps, which do not always correlate with gas. Because radial gas flow induced by low-mass, non-gas-gap-opening planets could affect the radial drift of dust, it potentially forms these dust substructures in disks. We investigate the potential of gas flow induced by low-mass planets to sculpt the rings and gaps in the dust profiles. We first perform three-dimensional hydrodynamical simulations, which resolve the local gas flow past a planet. We then calculate the trajectories of dust influenced by the planet-induced gas flow. Finally, we compute the steady-state dust surface density by incorporating the influences of the planet-induced gas flow into a one-dimensional dust advection-diffusion model. The outflow of the gas toward the outside of the planetary orbit inhibits the radial drift of dust, leading to dust accumulation (the dust ring). The outflow toward the inside of the planetary orbit enhances the inward drift of dust, causing dust depletion around the planetary orbit (the dust gap). Under weak turbulence ($\alpha_{\rm diff}\lesssim10^{-4}$, where $\alpha_{\rm diff}$ is the turbulence strength parameter), the gas flow induced by the planet with $\gtrsim1\,M_{\oplus}$ (Earth mass) generates the dust ring and gap in the distribution of small dust grains ($\lesssim1$ cm) with the radial extent of $\sim1\text{--}10$ times gas scale height around the planetary orbit without creating a gas gap and pressure bump. The gas flow induced by low-mass, non-gas-gap-opening planets can be considered a possible origin of the observed dust substructures in disks. Our results may be helpful to explain the disks whose dust substructures were found not to correlate with those of the gas., Comment: 25 pages, 20 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

Earth's surface environment is largely influenced by its budget of major volatile elements: carbon (C), nitrogen (N), and hydrogen (H). Although the volatiles on Earth are thought to have been delivered by chondritic materials, the elemental composition of the bulk silicate Earth (BSE) shows depletion in the order of N, C, and H. Previous studies have concluded that non-chondritic materials are needed for this depletion pattern. Here, we model the evolution of the volatile abundances in the atmosphere, oceans, crust, mantle, and core through the accretion history by considering elemental partitioning and impact erosion. We show that the BSE depletion pattern can be reproduced from continuous accretion of chondritic bodies by the partitioning of C into the core and H storage in the magma ocean in the main accretion stage and atmospheric erosion of N in the late accretion stage. This scenario requires a relatively oxidized magma ocean ($\log_{10} f_{\rm O_2}$ $\gtrsim$ $\rm{IW}$$-2$, where $f_{\rm O_2}$ is the oxygen fugacity, ${\rm IW}$ is $\log_{10} f_{\rm O_2}^{\rm IW}$, and $f_{\rm O_2}^{\rm IW}$ is $f_{\rm O_2}$ at the iron-w\"{u}stite buffer), the dominance of small impactors in the late accretion, and the storage of H and C in oceanic water and carbonates in the late accretion stage, all of which are naturally expected from the formation of an Earth-sized planet in the habitable zone., Comment: 17 pages, 4 figures (comprising 11 panels in total), 1 table, Methods, and Supplementary Information (available online: https://www.nature.com/articles/s41598-021-99240-w)

Pebble accretion is one of the major theories in planet formation. Aerodynamically small particles, called pebbles, are highly affected by the gas flow. A growing planet embedded in a protoplanetary disk induces three-dimensional (3D) gas flow. In our previous study, Paper I, we focused on the shear regime of pebble accretion, and investigated the influence of planet-induced gas flow on pebble accretion. In Paper I, we found that pebble accretion is inefficient in the planet-induced gas flow compared to that in the unperturbed flow, in particular when ${\rm St}\lesssim10^{-3}$, where St is the Stokes number. Following Paper I, we investigate the influence of planet-induced gas flow on pebble accretion. We consider the headwind of the gas, which is not included in Paper I. We extend our study to the headwind regime of pebble accretion in this study. Assuming a nonisothermal, inviscid sub-Keplerian gas disk, we perform 3D hydrodynamical simulations on the spherical polar grid which has a planet with the dimensionless mass, $m=R_{\rm Bondi}/H$, located at its center, where $R_{\rm Bondi}$ and $H$ are the Bondi radius and the disk scale height. We then numerically integrate the equation of motion of pebbles in 3D using hydrodynamical simulation data. Combining our results with the spacial variety of turbulence strength and pebble size in a disk, we conclude that the planet-induced gas flow still allows for pebble accretion in the early stage of planet formation. Suppression of pebble accretion due to the planet-induced gas flow occurs only in the late stage of planet formation, in particular in the inner region of the disk. This may be helpful to explain the distribution of exoplanets and the architecture of the Solar System, both of which have small inner and large outer planets., Comment: 24 pages, 23 figures, Accepted for publication in Astronomy and Astrophysics (A&A). arXiv admin note: text overlap with arXiv:1911.12999

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., Comment: 19 pages, 13 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

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., Comment: 16 pages, 14 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

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., Comment: 35 pages, 9 figures, accepted for publication in EPSL

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., Comment: 15 pages, 11 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

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., Comment: 31 pages, 12 Figures, accepted for publication in Icarus

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., Comment: Accepted for Publication in Icarus (54 pages, 16 figures)

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., Comment: 31 pages, 9 figures, accepted for publication in ApJ

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., Comment: 31 pages, 8 figures, accepted for publication in Geochemical Journal

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., Comment: 29 pages, 4 figures, 2 tables, accepted for publication in Earth and Planetary Science Letters

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., Comment: 43 pages, 10 figures, accepted to be published in ApJ

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., Comment: 7 pages, 7 figures, accepted for publication in MNRAS