Showing 1–50 of 54 results for author: Genda, H

Large planets may not form fractionally large moons

Authors: Miki Nakajima, Hidenori Genda, Erik Asphaug, Shigeru Ida

Abstract: One of the unique aspects of Earth is that it has a fractionally large Moon, which is thought to have formed from a Moon-forming disk generated by a giant impact. The Moon stabilizes the Earth's spin axis at least by several degrees and contributes to Earth's stable climate. Given that impacts are common during planet formation, exomoons, which are moons around planets in extrasolar systems, shoul… ▽ More One of the unique aspects of Earth is that it has a fractionally large Moon, which is thought to have formed from a Moon-forming disk generated by a giant impact. The Moon stabilizes the Earth's spin axis at least by several degrees and contributes to Earth's stable climate. Given that impacts are common during planet formation, exomoons, which are moons around planets in extrasolar systems, should be common as well, but no exomoon has been confirmed. Here we propose that an initially vapor-rich moon-forming disk is not capable of forming a large moon because growing moonlets, which are building blocks of a moon, experience strong gas drag and quickly fall toward the planet. Our impact simulations show that terrestrial and icy planets that are larger than $\sim 1.3-1.6 R_\oplus$ produce entirely vapor disks, which fail to form a large moon. This indicates that (1) our model supports the Moon-formation models that produce vapor-poor disks and (2) rocky and icy exoplanets whose radii are smaller than $\sim 1.6 R_\oplus$ are ideal candidates for hosting fractionally large exomoons. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Comments: This paper was published in 2022

Journal ref: Nature Communications 13, 568 (2022)

Mixing model of Phobos' bulk elemental composition for the determination of its origin: Multivariate analysis of MMX/MEGANE data

Authors: Kaori Hirata, Tomohiro Usui, Ryuki Hyodo, Hidenori Genda, Ryota Fukai, David J. Lawrence, Nancy L. Chabot, Patrick N. Peplowski, Hiroki Kusano

Abstract: The formation process of the two Martian moons, Phobos and Deimos, is still debated with two main competing hypotheses: the capture of an asteroid or a giant impact onto Mars. In order to reveal their origin, the Martian Moons eXploration (MMX) mission by Japan Aerospace Exploration Agency (JAXA) plans to measure Phobos' elemental composition by a gamma-ray and neutron spectrometer called MEGANE.… ▽ More The formation process of the two Martian moons, Phobos and Deimos, is still debated with two main competing hypotheses: the capture of an asteroid or a giant impact onto Mars. In order to reveal their origin, the Martian Moons eXploration (MMX) mission by Japan Aerospace Exploration Agency (JAXA) plans to measure Phobos' elemental composition by a gamma-ray and neutron spectrometer called MEGANE. This study provides a model of Phobos' bulk elemental composition, assuming the two formation hypotheses. Using the mixing model, we established a MEGANE data analysis flow to discriminate between the formation hypotheses by multivariate analysis. The mixing model expresses the composition of Phobos in 6 key lithophile elements that will be measured by MEGANE (Fe, Si, O, Ca, Mg, and Th) as a linear mixing of two mixing components: material from Mars and material from an asteroid as represented by primitive meteorite compositions. The inversion calculation includes consideration of MEGANE's measurement errors ($E_P$) and derives the mixing ratio for a given Phobos composition, based on which the formation hypotheses are judged. For at least 65\% of the modeled compositions, MEGANE measurements will determine the origin uniquely ($E_P$ = 30\%), and this increases from 74 to 87\% as $E_P$ decreases from 20 to 10\%. Although the discrimination performance depends on $E_P$, the current operation plan for MEGANE predicts an instrument performance for $E_P$ of 20--30\%, resulting in ~70\% discrimination between the original hypotheses. MEGANE observations can also enable the determination of the asteroid type of the captured body or the impactor. The addition of other measurements, such as MEGANE's measurements of the volatile element K, as well as observations by other MMX remote sensing instruments, will also contribute to the MMX mission's goal to constrain the origin of Phobos. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 34 pages, 7 figures, accepted for publication in Icarus

Exploring the recycling model of Phobos formation: rubble-pile satellites

Authors: Gustavo Madeira, Sebastien Charnoz, Yun Zhang, Ryuki Hyodo, Patrick Michel, Hidenori Genda, Silvia Giuliatti Winter

Abstract: Phobos is the target of the return sample mission Martian Moons eXploration by JAXA that will analyze in great details the physical and compositional properties of the satellite from orbit, from the surface and in terrestrial laboratories, giving clues about its formation. Some models propose that Phobos and Deimos were formed after a giant impact giving rise to an extended debris disk. Assuming t… ▽ More Phobos is the target of the return sample mission Martian Moons eXploration by JAXA that will analyze in great details the physical and compositional properties of the satellite from orbit, from the surface and in terrestrial laboratories, giving clues about its formation. Some models propose that Phobos and Deimos were formed after a giant impact giving rise to an extended debris disk. Assuming that Phobos formed from a cascade of disruptions and re-accretions of several parent bodies in this disk, and that they are all characterized by a low material cohesion, Hesselbrock & Milton (2017) have showed that a recycling process may happen during the assembling of Phobos, by which Phobos' parents are destroyed into a Roche-interior ring and reaccreted several times. In the current paper we explore in details the recycling model, and pay particular attention to the characteristics of the disk using 1D models of disk/satellite interactions. In agreement with previous studies we confirm that, if Phobos' parents bodies are gravitational aggregates (rubble piles), then the recycling process does occur. However, Phobos should be accompanied today by a Roche-interior ring. Furthermore, the characteristics of the ring are not reconcilable with today`s observations of Mars' environment, which put stringent constraints on the existence of a ring around Mars. The recycling mechanism may or may not have occurred at the Roche limit for an old moon population, depending on their internal cohesion. However, the Phobos we see today cannot be the outcome of such a recycling process. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Comments: Accept in The Astronomical Journal

Early water delivery to terrestrial planet regions during the stages of Jupiter's formation and migration in the Grand Tack model

Authors: Masahiro Ogihara, Hidenori Genda, Yasuhito Sekine

Abstract: The formation and subsequent migration of gas giants could significantly affect the material mixing in the Solar System. In this study, we use N-body simulations to investigate how much water is transported into the region of the terrestrial planet formation during the growth and migration phases of Jupiter in the Grand Tack model. We found that Jupiter's growth was accompanied by significant mass… ▽ More The formation and subsequent migration of gas giants could significantly affect the material mixing in the Solar System. In this study, we use N-body simulations to investigate how much water is transported into the region of the terrestrial planet formation during the growth and migration phases of Jupiter in the Grand Tack model. We found that Jupiter's growth was accompanied by significant mass transport, and that a substantial amount of water (about 10 times Earth's ocean mass for the initial planetesimal distribution based on the minimum-mass solar nebula) was transported into the terrestrial planet region. The total amount delivered increased further during Jupiter's migration phase (totaling about 10-40 times Earth's ocean mass), which was less dependent on simulation parameters. In addition, at these stages, terrestrial planets were not fully grown. Therefore, water supplied during these early stages could interact with metallic iron during the core formation of protoplanets and/or growing Earth. Since hydrogen in water molecules can dissolve into their cores, this could explain the density deficit observed in the current Earth core. Notably, Jupiter could play an important role as a ''barrier'' in explaining the dichotomy of the isotopic compositions between noncarbonaceous (NC) and carbonaceous (CC) meteorites. This study's results show that Jupiter's growth necessitates some mixing of NC and CC materials. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 14 pages, 12 figures, accepted for publication in PSJ, a 6-minute video summary can be found here: https://exoplanet-talks.org/talk/405

Experimentally shock-induced melt veins in basalt: Improving the shock classification of eucrites

Authors: Haruka Ono, Kosuke Kurosawa, Takafumi Niihara, Takashi Mikouchi, Naotaka Tomioka, Junko Isa, Hiroyuki Kagi, Takuya Matsuzaki, Hiroshi Sakuma, Hidenori Genda, Tatsuhiro Sakaiya, Tadashi Kondo, Masahiro Kayama, Mizuho Koike, Yuji Sano, Masafumi Murayama, Wataru Satake, Takafumi Matsui

Abstract: Basaltic rocks occur widely on the terrestrial planets and differentiated asteroids, including the asteroid 4 Vesta. We conducted a shock recovery experiment with decaying compressive pulses on a terrestrial basalt at Chiba Institute of Technology, Japan. The sample recorded a range of pressures, and shock physics modeling was conducted to add a pressure scale to the observed shock features. The s… ▽ More Basaltic rocks occur widely on the terrestrial planets and differentiated asteroids, including the asteroid 4 Vesta. We conducted a shock recovery experiment with decaying compressive pulses on a terrestrial basalt at Chiba Institute of Technology, Japan. The sample recorded a range of pressures, and shock physics modeling was conducted to add a pressure scale to the observed shock features. The shocked sample was examined by optical and electron microscopy, electron back-scattered diffractometry, and Raman spectroscopy. We found that localized melting occurs at a lower pressure (~10 GPa) than previously thought (>20 GPa). The shocked basalt near the epicenter represents shock degree C of a recently proposed classification scheme for basaltic eucrites and, as such, our results provide a pressure scale for the classification scheme. Finally, we estimated the total fraction of the basaltic eucrites classified as shock degree C to be ~15% by assuming the impact velocity distribution onto Vesta. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 17 pages, 4 figures, 1 table, 1 supplementary Information, accepted for publication in Geophysical Research Letters

Effect of impact velocity and angle on deformational heating and post-impact temperature

Authors: Shigeru Wakita, Hidenori Genda, Kosuke Kurosawa, Thomas M. Davison, Brandon C. Johnson

Abstract: The record of impact induced shock-heating in meteorites is an important key for understanding the collisional history of the solar system. Material strength is important for impact heating, but the effect of impact angle and impact velocity on shear heating remains poorly understood. Here, we report three-dimensional oblique impact simulations, which confirm the enhanced heating due to material s… ▽ More The record of impact induced shock-heating in meteorites is an important key for understanding the collisional history of the solar system. Material strength is important for impact heating, but the effect of impact angle and impact velocity on shear heating remains poorly understood. Here, we report three-dimensional oblique impact simulations, which confirm the enhanced heating due to material strength and explore the effects of impact angle and impact velocity. We find that oblique impacts with an impact angle that is steeper than 45 degree produce a similar amount of heated mass as vertical impacts. On the other hand, grazing impacts produce less heated mass and smaller heated regions compared to impacts at steeper angles. We derive an empirical formula of the heated mass, as a function of the impact angle and velocity. This formula can be used to estimate the impact conditions (velocity and angle) that had occurred and caused Ar loss in the meteoritic parent bodies. Furthermore, our results indicate that grazing impacts at higher impact velocities could generate a similar amount of heated material as vertical impacts at lower velocities. As the heated material produced by grazing impacts has experienced lower pressure than the material heated by vertical impacts, our results imply that grazing impacts may produce weakly shock-heated meteorites. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 21 pages, 11 figures, accepted for publication in JGR: Planets

Challenges in forming Phobos and Deimos directly from a splitting of an ancestral single moon

Authors: Ryuki Hyodo, Hidenori Genda, Ryosuke Sekiguchi, Gustavo Madeira, Sébastien Charnoz

Abstract: The origin and evolution of Martian moons have been intensively debated in recent years. It is proposed that Phobos and Deimos may originate directly from a splitting of an ancestral moon orbiting at around the Martian synchronous orbit. At this hypothetical splitting, the apocenter of the inner moon (presumed as Phobos) and the pericenter of the outer moon (presumed as Deimos) are reported to coi… ▽ More The origin and evolution of Martian moons have been intensively debated in recent years. It is proposed that Phobos and Deimos may originate directly from a splitting of an ancestral moon orbiting at around the Martian synchronous orbit. At this hypothetical splitting, the apocenter of the inner moon (presumed as Phobos) and the pericenter of the outer moon (presumed as Deimos) are reported to coincide, in that, their semi-major axes reside inside and outside the Martian synchronous orbit with non-zero eccentricities, respectively. However, the successive orbital evolution of the two moons is not studied. Here, we perform direct $N$-body orbital integrations of the moons, including the Martian oblateness of the $J_2$ and $J_4$ terms. We show that the two moons, while they precess, likely collide within $\sim 10^4$ years with an impact velocity of $v_{\rm imp} \sim 100-300$ m s$^{-1}$ ($\sim 10-30$ times moons' escape velocity) and with an isotropic impact direction. The impact occurs around the apocenter and the pericenter of the inner and outer moons, respectively, where the timescale of this periodic orbital alignment is regulated by the precession. By performing additional impact simulations, we show that such a high-velocity impact likely results in a disruptive outcome, forming a debris ring at around the Martian synchronous orbit, from which several small moons would accrete. Such an evolutionary path would eventually form a different Martian moons system from the one we see today. Therefore, it seems unlikely that Phobos and Deimos are split directly from a single ancestral moon. △ Less

Submitted 22 August, 2022; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: 11 pages, 5 Figures, accepted for publication in Planetary Science Journal (PSJ)

Shock recovery with decaying compressive pulses: Shock effects in calcite (CaCO$_3$) around the Hugoniot elastic limit

Authors: Kosuke Kurosawa, Haruka Ono, Takafumi Niihara, Tatsuhiro Sakaiya, Tadashi Kondo, Naotaka Tomioka, Takashi Mikouchi, Hidenori Genda, Takuya Matsuzaki, Masahiro Kayama, Mizuho Koike, Yuji Sano, Masafumi Murayama, Wataru Satake, Takafumi Matsui

Abstract: Shock metamorphism of minerals in meteorites provides insights into the ancient Solar System. Calcite is an abundant aqueous alteration mineral in carbonaceous chondrites. Return samples from the asteroids Ryugu and Bennu are expected to contain calcite-group minerals. Although shock metamorphism in silicates has been well studied, such data for aqueous alteration minerals are limited. Here, we in… ▽ More Shock metamorphism of minerals in meteorites provides insights into the ancient Solar System. Calcite is an abundant aqueous alteration mineral in carbonaceous chondrites. Return samples from the asteroids Ryugu and Bennu are expected to contain calcite-group minerals. Although shock metamorphism in silicates has been well studied, such data for aqueous alteration minerals are limited. Here, we investigated the shock effects in calcite with marble using impact experiments at the Planetary Exploration Research Center of Chiba Institute of Technology. We produced decaying compressive pulses with a smaller projectile than the target. A metal container facilitates recovery of a sample that retains its pre-impact stratigraphy. We estimated the peak pressure distributions in the samples with the iSALE shock physics code. The capability of this method to produce shocked grains that have experienced different degrees of metamorphism from a single experiment is an advantage over conventional uniaxial shock recovery experiments. The shocked samples were investigated by polarizing microscopy and X-ray diffraction analysis. We found that more than half of calcite grains exhibit undulatory extinction when peak pressure exceeds 3 GPa. This shock pressure is one order of magnitude higher than the Hugoniot elastic limit (HEL) of marble, but it is close to the HEL of a calcite crystal, suggesting that the undulatory extinction records dislocation-induced plastic deformation in the crystal. Finally, we propose a strategy to re-construct the maximum depth of calcite grains in a meteorite parent body, if shocked calcite grains are identified in chondrites and/or return samples from Ryugu and Bennu. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Comments: 25 pages, 9 figures, 1 table, 1 supporting information, accepted for publication in Journal of Geophysical Research Planets

Giant Impact onto a Vesta-Like Asteroid and Formation of Mesosiderites through Mixing of Metallic Core and Surface Crust

Authors: Keisuke Sugiura, Makiko K. Haba, Hidenori Genda

Abstract: Mesosiderites are a type of stony-iron meteorites composed of a mixture of silicates and Fe-Ni metals. The mesosiderite silicates and metals are considered to have originated from the crust and metal core, respectively, of a differentiated asteroid. In contrast, mesosiderites rarely contain the olivine that is mainly included in a mantle. Although a giant impact onto a differentiated asteroid is c… ▽ More Mesosiderites are a type of stony-iron meteorites composed of a mixture of silicates and Fe-Ni metals. The mesosiderite silicates and metals are considered to have originated from the crust and metal core, respectively, of a differentiated asteroid. In contrast, mesosiderites rarely contain the olivine that is mainly included in a mantle. Although a giant impact onto a differentiated asteroid is considered to be a probable mechanism to mix crust and metal materials to form mesosiderites, it is not obvious how such a giant impact can form mesosiderite-like materials without including mantle materials. We conducted numerical simulations of giant impacts onto differentiated asteroids, using the smoothed particle hydrodynamics method, to investigate the detailed distribution of mixed materials on the resultant bodies. For the internal structure of a target body, we used a thin-crust model derived from the magma ocean crystallization model of the asteroid Vesta and a thick-crust and a large-core model suggested from the proximity observation of Vesta by the Dawn probe. In the simulations with the former model, excavation of the metal core requires nearly catastrophic impacts and mantle is exposed over large surface areas. Thus, stony-iron materials produced on its surface are likely to include mantle materials and it is difficult to produce mesosiderite-like materials. Conversely, in the simulations with the latter model, mantle materials are exposed only at impact sites, even when the impacts excavate the metal core, and the formation of a surface with little mantle material and the formation of mesosiderite-like materials are possible. Therefore, our simulations suggest that an internal structure with a thick crust and a large core is more likely as a mesosiderite parent body rather than the thin-crust internal structure inferred from the conventional magma ocean model. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 21 pages, 8 figures, accepted for publication in Icarus

The Onset of a Globally Ice-covered State for a Land Planet

Authors: T. Kodama, H. Genda, J. Leconte, A. Abe-Ouchi

Abstract: The climates of terrestrial planets with a small amount of water on their surface, called land planets, are significantly different from the climates of planets having a large amount of surface water. Land planets have a higher runaway greenhouse threshold than aqua planets, which extends the inner edge of the habitable zone inward. Land planets also have the advantage of avoiding global freezing… ▽ More The climates of terrestrial planets with a small amount of water on their surface, called land planets, are significantly different from the climates of planets having a large amount of surface water. Land planets have a higher runaway greenhouse threshold than aqua planets, which extends the inner edge of the habitable zone inward. Land planets also have the advantage of avoiding global freezing due to drier tropics, leading to a lower planetary albedo. In this study, we systematically investigate the complete freezing limit for various surface water distribution using a three-dimensional dynamic atmospheric model. As in a previous study, we found that a land planet climate has dry tropics that result in less snow and fewer clouds. The complete freezing limit decreases from that for aqua planets (92% S0, where S0 is Earth's present insolation) to that for land planets (77% S0) with an increasing dry area. Values for the complete freezing limit for zonally uniform surface water distributions are consistently lower than those for meridionally uniform surface water distribution. This is because the surface water distribution in the tropics in the meridionally uniform cases causes ice-albedo feedback until a planet lapses into the complete freezing state. For a surface water distribution using the topographies of the terrestrial planets, the complete freezing limit has values near those for the meridionally uniform cases. Our results indicate that the water distribution is important for the onset of a global ice-covered state for Earth-like exoplanets. △ Less

Submitted 25 November, 2021; originally announced November 2021.

Comments: 29 pages, 8 figures. Accepted for publication in Journal of Geophysical Research - Planets

Numerous chondritic impactors and oxidized magma ocean set Earth's volatile depletion

Authors: Haruka Sakuraba, Hiroyuki Kurokawa, Hidenori Genda, Kenji Ohta

Abstract: 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 ne… ▽ More 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ü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. △ Less

Submitted 23 October, 2021; originally announced October 2021.

Comments: 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)

Journal ref: Scientific Reports, volume 11, Article number: 20894 (2021)

Tidal evolution of the eccentric moon around dwarf planet (225088) Gonggong

Authors: Sota Arakawa, Ryuki Hyodo, Daigo Shoji, Hidenori Genda

Abstract: Recent astronomical observations revealed that (225088) Gonggong, a 1000-km-sized trans-Neptunian dwarf planet, hosts an eccentric satellite, Xiangliu, with an eccentricity of approximately 0.3. As the majority of known satellite systems around trans-Neptunian dwarf planets have circular orbits, the observed eccentricity of Gonggong--Xiangliu system may reflect the singular properties of the syste… ▽ More Recent astronomical observations revealed that (225088) Gonggong, a 1000-km-sized trans-Neptunian dwarf planet, hosts an eccentric satellite, Xiangliu, with an eccentricity of approximately 0.3. As the majority of known satellite systems around trans-Neptunian dwarf planets have circular orbits, the observed eccentricity of Gonggong--Xiangliu system may reflect the singular properties of the system. In this study, we assumed that Gonggong--Xiangliu system formed via a giant impact and investigated the following secular tidal evolution of Gonggong--Xiangliu system under the simplifying assumption of homogeneous bodies and of zero orbital inclination. We conducted coupled thermal--orbital evolution simulations using the Andrade viscoelastic model and included higher-order eccentricity functions. The distribution of the final eccentricity from a large number of simulations with different initial conditions revealed that the radius of Xiangliu is not larger than 100 km. We also derived the analytical solution of the semilatus rectum evolution, a function of the radius of Xiangliu. From the point of view of the final semilatus rectum, the radius of Xiangliu was estimated to be close to 100 km. Together with the results of the Hubble Space Telescope observations, our findings suggest Gonggong and Xiangliu have similar albedos. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: 36 pages, 22 figures. Accepted for publication in AJ

SPH Simulations for Shape Deformation of Rubble-Pile Asteroids Through Spinup: The Challenge for Making Top-Shaped Asteroids Ryugu and Bennu

Authors: Keisuke Sugiura, Hiroshi Kobayashi, Sei-ichiro Watanabe, Hidenori Genda, Ryuki Hyodo, Shu-ichiro Inutsuka

Abstract: Asteroid Ryugu and asteroid Bennu, which were recently visited by spacecraft Hayabusa2 and OSIRIS-REx, respectively, are spinning top-shaped rubble piles. Other axisymmetric top-shaped near-Earth asteroids have been observed with ground-based radar, most of which rotate near breakup rotation periods of ~ 3 hours. This suggests that rotation-induced deformation of asteroids through rotational spinu… ▽ More Asteroid Ryugu and asteroid Bennu, which were recently visited by spacecraft Hayabusa2 and OSIRIS-REx, respectively, are spinning top-shaped rubble piles. Other axisymmetric top-shaped near-Earth asteroids have been observed with ground-based radar, most of which rotate near breakup rotation periods of ~ 3 hours. This suggests that rotation-induced deformation of asteroids through rotational spinup produces top shapes. Although some previous simulations using the Discrete Element Method showed that spinup of rubble piles may produce oblate top shapes, it is still unclear what kinds of conditions such as friction angles of constituent materials and spinup timescales are required for top-shape formation. Here we show, through Smoothed Particle Hydrodynamics simulations of granular bodies spinning-up at different rates, that the rotation-induced deformation of spherical rubble piles before breakup can be classified into three modes according to the friction angle φ_{d}: quasi-static and internal deformation for φ_{d} < 40 degrees, dynamical and internal deformation for 50 degrees < φ_{d} < 60 degrees, and surface landslides for φ_{d} > 70 degrees. Note that these apparent large values of friction angle can be acceptable if we consider the effect of cohesion among blocks of a rubble pile under weak gravity. Bodies with φ_{d} < 60 degrees evolve into oblate spheroids through internal deformation, but never form pronounced equators defining a top shape. In contrast, bodies with φ_{d} > 70 degrees deform into axisymmetric top shapes through an axisymmetric surface landslides if spinup timescales are < a few days. In addition, through slow spinups with timescales > 1 month, bodies with φ_{d} > 70 degrees deform into non-axisymmetric shapes via localized landslides. We suggest that rapid spinup mechanisms are preferable for the formation of axisymmetric top shapes. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 29 pages, 11 figures, accepted for publication in Icarus

Erosion and accretion by cratering impacts on rocky and icy bodies

Authors: Ryuki Hyodo, Hidenori Genda

Abstract: During planet formation, numerous small impacting bodies result in cratering impacts on large target bodies. A fraction of the target surface is eroded, while a fraction of the impactor material accretes onto the surface. These fractions depend upon the impact velocities, the impact angles, and the escape velocities of the target. This study uses smoothed particle hydrodynamics simulations to mode… ▽ More During planet formation, numerous small impacting bodies result in cratering impacts on large target bodies. A fraction of the target surface is eroded, while a fraction of the impactor material accretes onto the surface. These fractions depend upon the impact velocities, the impact angles, and the escape velocities of the target. This study uses smoothed particle hydrodynamics simulations to model cratering impacts onto a planar icy target for which gravity is the dominant force and material strength is neglected. By evaluating numerical results, scaling laws are derived for the escape mass of the target material and the accretion mass of the impactor material onto the target surface. Together with recently derived results for rocky bodies in a companion study, a conclusion is formulated that typical cratering impacts on terrestrial planets, except for those on Mercury, led to a net accretion, while those on the moons of giant planets, e.g., Rhea and Europa, led to a net erosion. Our newly derived scaling laws would be useful for predicting the erosion of the target body and the accretion of the impactor for a variety of cratering impacts that would occur on large rocky and icy planetary bodies during planet formation and collisional evolution from ancient times to today. △ Less

Submitted 11 April, 2021; originally announced April 2021.

Comments: 12 pages, 8 figures. Accepted for publication in ApJ. arXiv admin note: substantial text overlap with arXiv:2006.00732

The role of post-shock heating by plastic deformation during impact devolatilization of calcite

Authors: Kosuke Kurosawa, Hidenori Genda, Shintaro Azuma, Keishi Okazaki

Abstract: An accurate understanding of the relationship between the impact conditions and the degree of shock-induced thermal metamorphism in meteorites allows the impact environment in the early Solar System to be understood. A recent hydrocode has revealed that impact heating is much higher than previously thought. This is because plastic deformation of the shocked rocks causes further heating during deco… ▽ More An accurate understanding of the relationship between the impact conditions and the degree of shock-induced thermal metamorphism in meteorites allows the impact environment in the early Solar System to be understood. A recent hydrocode has revealed that impact heating is much higher than previously thought. This is because plastic deformation of the shocked rocks causes further heating during decompression, which is termed post-shock heating. Here we compare impact simulations with laboratory experiments on the impact devolatilization of calcite to investigate whether the post-shock heating is also significant in natural samples. We calculated the mass of CO$_2$ produced from the calcite, based on thermodynamics. We found that iSALE can reproduce the devolatilization behavior for rocks with the strength of calcite. In contrast, the calculated masses of CO2 at lower rock strengths are systematically smaller than the experimental values. Our results require a reassessment of the interpretation of thermal metamorphism in meteorites. △ Less

Submitted 4 March, 2021; originally announced March 2021.

Comments: 30 pages, 4 figures, 1 Supporting Information, accepted for publication in Geophysical Research Letters

Modification of the composition and density of Mercury from late accretion

Authors: Ryuki Hyodo, Hidenori Genda, Ramon Brasser

Abstract: Late accretion is a process that strongly modulated surface geomorphic and geochemical features of Mercury. Yet, the fate of the impactors and their effects on Mercury's surface through the bombardment epoch are not clear. Using Monte-Carlo and analytical approaches of cratering impacts, we investigate the physical and thermodynamical outcomes of late accretion on Mercury. Considering the uncertai… ▽ More Late accretion is a process that strongly modulated surface geomorphic and geochemical features of Mercury. Yet, the fate of the impactors and their effects on Mercury's surface through the bombardment epoch are not clear. Using Monte-Carlo and analytical approaches of cratering impacts, we investigate the physical and thermodynamical outcomes of late accretion on Mercury. Considering the uncertainties in late accretion, we develop scaling laws for the following parameters as a function of impact velocity and total mass of late accretion: (1) depth of crustal erosion, (2) the degree of resurfacing, and (3) mass accreted from impactor material. Existing dynamical models indicate that Mercury experienced an intense impact bombardment (a total mass of $\sim 8 \times 10^{18} - 8 \times 10^{20}$ kg with a typical impact velocity of $30-40$ km s$^{-1}$) after $4.5$ Ga. With this, we find that late accretion could remove 50 m to 10 km of the early crust of Mercury, but the change to its core-to-mantle ratio is negligible. Alternatively, the mantles of putative differentiated planetesimals in the early solar system could be more easily removed and their respective core fraction increased, if Mercury ultimately accreted from such objects. Although the cratering is notable for erasing the older geological surface records on Mercury, we show that $\sim 40-50$ wt% of the impactor's exogenic materials, including the volatile-bearing materials, can be heterogeneously implanted on Mercury's surface as a late veneer ($>3\times 10^{18}-1.6 \times 10^{19}$ kg in total). About half of the accreted impactor's materials are vaporized, and the rest is completely melted upon the impact. We expect that the further interplay between our theoretical results and forthcoming surface observations of Mercury, including the BepiColombo mission, will lead us to a better understanding of Mercury's origin and evolution. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: 19 pages, 7 figures. Accepted for publication in Icarus

Escape and accretion by cratering impacts: Formulation of scaling relations for high-speed ejecta

Authors: Ryuki Hyodo, Hidenori Genda

Abstract: Numerous small bodies inevitably lead to cratering impacts on large planetary bodies during planet formation and evolution. As a consequence of these small impacts, a fraction of the target material escapes from the gravity of the large body, and a fraction of the impactor material accretes onto the target surface, depending on the impact velocities and angles. Here, we study the mass of the high-… ▽ More Numerous small bodies inevitably lead to cratering impacts on large planetary bodies during planet formation and evolution. As a consequence of these small impacts, a fraction of the target material escapes from the gravity of the large body, and a fraction of the impactor material accretes onto the target surface, depending on the impact velocities and angles. Here, we study the mass of the high-speed ejecta that escapes from the target gravity by cratering impacts when material strength is neglected. We perform a large number of cratering impact simulations onto a planar rocky target using the smoothed particle hydrodynamics method. We show that the escape mass of the target material obtained from our numerical simulations agrees with the prediction of a scaling law under a point-source assumption when $v_{\rm imp} \gtrsim 12 v_{\rm esc}$, where $v_{\rm imp}$ is the impact velocity and $v_{\rm esc}$ is the escape velocity of the target. However, we find that the point-source scaling law overestimates the escape mass up to a factor of $\sim 70$, depending on the impact angle, when $v_{\rm imp} \lesssim 12 v_{\rm esc}$. Using data obtained from numerical simulations, we derive a new scaling law for the escape mass of the target material for $v_{\rm imp} \lesssim 12 v_{\rm esc}$. We also derive a scaling law that predicts the accretion mass of the impactor material onto the target surface upon cratering impacts by numerically evaluating the escape mass of the impactor material. Our newly derived scaling laws are useful for predicting the escape mass of the target material and the accretion mass of the impactor material for a variety of cratering impacts that would occur on large planetary bodies during planet formation. △ Less

Submitted 1 June, 2020; originally announced June 2020.

Comments: 13 pages, 8 figures. Accepted for publication in ApJ

Impact Ejecta near the Impact Point Observed using Ultra-high-speed Imaging and SPH Simulations, and a Comparison of the Two Methods

Authors: Takaya Okamoto, Kosuke Kurosawa, Hidenori Genda, Takafumi Matsui

Abstract: High-speed impact ejecta at velocities comparable to the impact velocity are expected to contribute to material transport between planetary bodies and deposition of ejecta far from the impact crater. We investigated the behavior of high-speed ejecta produced at angles of 45 and 90 degrees, using both experimental and numerical methods. The experimental system developed at the Planetary Exploration… ▽ More High-speed impact ejecta at velocities comparable to the impact velocity are expected to contribute to material transport between planetary bodies and deposition of ejecta far from the impact crater. We investigated the behavior of high-speed ejecta produced at angles of 45 and 90 degrees, using both experimental and numerical methods. The experimental system developed at the Planetary Exploration Research Center of Chiba Institute of Technology (Japan) allowed us to observe the initial growth of the ejecta. We succeeded in imaging high-speed ejecta at 0.2 $\mathrm{μ}$s intervals for impacts of polycarbonate projectiles of 4.8 mm diameter onto a polycarbonate plate at an impact velocity of ~4 km s$^{-1}$. Smoothed particle hydrodynamics (SPH) simulations of various numerical resolutions were conducted for the same impact conditions as pertaining to the experiments. We compared the morphology and velocities of the ejecta for the experiments and simulations, and we confirmed a close match for high-resolution simulations (with $\geq10^6$ SPH particles representing the projectile). According to the ejecta velocity distributions obtained from our high-resolution simulations, the ejection velocities of the high-speed ejecta for oblique impacts are much greater than those for vertical impacts. The translational motion of penetrating projectiles parallel to the target surface in oblique impacts could cause long-term, sustained acceleration at the root of the ejecta. △ Less

Submitted 18 March, 2020; originally announced March 2020.

Comments: 33 pages, 16 figures, 1 supporting information, accepted for publication in Journal of Geophysical Research: Planets

Transport of impact ejecta from Mars to its moons as a means to reveal Martian history

Authors: Ryuki Hyodo, Kosuke Kurosawa, Hidenori Genda, Tomohiro Usui, Kazuhisa Fujita

Abstract: Throughout the history of the solar system, Mars has experienced continuous asteroidal impacts. These impacts have produced impact-generated Mars ejecta, and a fraction of this debris is delivered to Earth as Martian meteorites. Another fraction of the ejecta is delivered to the moons of Mars, Phobos and Deimos. Here, we studied the amount and condition of recent delivery of impact ejecta from Mar… ▽ More Throughout the history of the solar system, Mars has experienced continuous asteroidal impacts. These impacts have produced impact-generated Mars ejecta, and a fraction of this debris is delivered to Earth as Martian meteorites. Another fraction of the ejecta is delivered to the moons of Mars, Phobos and Deimos. Here, we studied the amount and condition of recent delivery of impact ejecta from Mars to its moons. Using state-of-the-art numerical approaches, we report, for the first time, that materials delivered from Mars to its moons are physically and chemically different from the Martian meteorites, which are all igneous rocks with a limited range of ages. We show that Mars ejecta mixed in the regolith of its moons potentially covers all its geological eras and consists of all types of rocks, from sedimentary to igneous. A Martian moons sample-return mission will bring such materials back to Earth, and the samples will provide a wealth of "time-resolved" geochemical information about the evolution of Martian surface environments. △ Less

Submitted 27 December, 2019; originally announced December 2019.

Comments: 9 pages, 3 figures, published in Scientific Reports

Journal ref: Sci Rep 9, 19833 (2019)

Enhancement of impact heating in pressure-strengthened rocks in oblique impacts

Authors: Shigeru Wakita, Hidenori Genda, Kosuke Kurosawa, Thomas M Davison

Abstract: Shock-induced metamorphism in meteorites informs us about the collisional environment and history of our solar system. Recently the importance of material strength in impact heating was reported from head-on impact simulations. Here, we perform three-dimensional oblique impact simulations, and confirm the additional heating due to material strength for oblique impacts. Despite a large difference i… ▽ More Shock-induced metamorphism in meteorites informs us about the collisional environment and history of our solar system. Recently the importance of material strength in impact heating was reported from head-on impact simulations. Here, we perform three-dimensional oblique impact simulations, and confirm the additional heating due to material strength for oblique impacts. Despite a large difference in the peak pressure at the impact point at a given impact velocity, we find that the heated mass for an oblique impact is nearly the same as that for a head-on impact. Thus, our results differ from the previous finding that the heated mass decreases as the impact becomes more oblique, and show that the additional shear heating is more effective for oblique impacts than for head-on impacts. This also indicates that material ejected during oblique impact tends to experience lower shock pressures but higher temperatures. △ Less

Submitted 1 December, 2019; originally announced December 2019.

Comments: 14 pages, 5 figures, accepted for publication in Geophysical Research Letters

The formation of the Martian moons

Authors: Pascal Rosenblatt, Ryuki Hyodo, Francesco C. Pignatale, Antony Trinh, Sébastien Charnoz, Kevin M. Dunseath, Mariko Terao-Dunseath, Hidenori Genda

Abstract: Almost all the planets of our solar system have moons. Each planetary system has however unique characteristics. The Martian system has not one single big moon like the Earth, not tens of moons of various sizes like for the giant planets, but two small moons: Phobos and Deimos. How did form such a system? This question is still being investigated on the basis of the Earth-based and space-borne obs… ▽ More Almost all the planets of our solar system have moons. Each planetary system has however unique characteristics. The Martian system has not one single big moon like the Earth, not tens of moons of various sizes like for the giant planets, but two small moons: Phobos and Deimos. How did form such a system? This question is still being investigated on the basis of the Earth-based and space-borne observations of the Martian moons and of the more modern theories proposed to account for the formation of other moon systems. The most recent scenario of formation of the Martian moons relies on a giant impact occurring at early Mars history and having also formed the so-called hemispheric crustal dichotomy. This scenario accounts for the current orbits of both moons unlike the scenario of capture of small size asteroids. It also predicts a composition of disk material as a mixture of Mars and impactor materials that is in agreement with remote sensing observations of both moon surfaces, which suggests a composition different from Mars. The composition of the Martian moons is however unclear, given the ambiguity on the interpretation of the remote sensing observations. The study of the formation of the Martian moon system has improved our understanding of moon formation of terrestrial planets: The giant collision scenario can have various outcomes and not only a big moon as for the Earth. This scenario finds a natural place in our current vision of the early solar system when conditions were favorable for giant collisions to occur. The next step in exploration of Martian moon is a sample return mission to test the giant collision scenario for their origin, and to provide tests of models of early solar system dynamics since Mars may retain material exchanged between the inner and outer solar system. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: 30 pages, 8 figures, 2 tables. Accepted for publication in the Oxford Research Encyclopedia of Planetary Science

Inner Edge of Habitable Zones for Earth-sized Planets with Various Surface Water Distributions

Authors: Takanori Kodama, Hidenori Genda, Ryouta O'ishi, Ayako Abe-Ouchi, Yutaka Abe

Abstract: When planets receive insolation above a certain critical value called the runaway threshold, liquid surface water vaporizes completely, which forms the inner edge of the habitable zone. Because land planets can emit a large amount of radiation from the dry tropics, they have a higher runaway threshold than aqua planets do. Here we systematically investigated the runaway threshold for various surfa… ▽ More When planets receive insolation above a certain critical value called the runaway threshold, liquid surface water vaporizes completely, which forms the inner edge of the habitable zone. Because land planets can emit a large amount of radiation from the dry tropics, they have a higher runaway threshold than aqua planets do. Here we systematically investigated the runaway threshold for various surface water distributions using a three-dimensional dynamic atmosphere model. The runaway threshold for the meridionally uniform surface water distribution increases from the typical value for the aqua-planet regime (~130% S0) to one for the land-planet regime (~155% S0) as the dry surface area increases, where S0 is the present Earth's insolation. Although this result is similar to the previous work considering zonally uniform surface water distributions, the runaway threshold for the land-planet regime is quite low compared to that of the previous work. This is because a part of the tropical atmosphere is always wet for the meridionally uniform case. We also considered the surface water distributions determined by the Earth's, Mars' and Venus' topographies. We found that their runaway thresholds are close to that for the meridionally uniform cases, and the amount of water at the boundary between an aqua- and land-planet regime is around 10% of the Earth's ocean. This clearly shows that the runaway threshold is not determined uniquely by the luminosity of the central star, but it has a wide range caused by the surface water distribution of the terrestrial water planet itself. △ Less

Submitted 16 August, 2019; originally announced August 2019.

Comments: 28 pages, 12 figures, accepted for publication in Journal of Geophysical Research - Planets

Assessment of the probability of microbial contamination for sample return from Martian moons II: The fate of microbes on Martian moons

Authors: Kosuke Kurosawa, Hidenori Genda, Ryuki Hyodo, Akihiko Yamagishi, Takashi Mikouchi, Takafumi Niihara, Shingo Matsuyama, Kazuhisa Fujita

Abstract: This paper presents a case study of microbe transportation in the Mars-satellites system. We examined the spatial distribution of potential impact-transported microbes on the Martian moons using impact physics by following a companion study (Fujita et al.). We used sterilization data from the precede studies. We considered that the microbes came mainly from the Zunil crater on Mars. We found that… ▽ More This paper presents a case study of microbe transportation in the Mars-satellites system. We examined the spatial distribution of potential impact-transported microbes on the Martian moons using impact physics by following a companion study (Fujita et al.). We used sterilization data from the precede studies. We considered that the microbes came mainly from the Zunil crater on Mars. We found that 70-80% of the microbes are likely to be dispersed all over the moon surface and are rapidly sterilized due to radiation except for those microbes within a thick ejecta deposit produced by meteoroids. The other 20-30% might be shielded from radiation by thick regolith layers that formed at collapsed layers in craters produced by Mars rock impacts. The total number of potentially surviving microbes at the thick ejecta deposits is estimated to be 3-4 orders of magnitude lower than at the Mars rock craters. The microbe concentration is irregular in the horizontal direction and is largely depth-dependent due to the radiation sterilization. The surviving fraction of transported microbes would be only 1 ppm on Phobos and 100 ppm on Deimos, suggesting that the transport processes and radiation severely affect microbe survival. The microbe sampling probability from the Martian moons was also investigated. We suggest that sample return missions from the Martian moons are classified into Unrestricted Earth-Return missions for 30 g samples and 10 cm depth sampling, even in our conservative scenario. We also conducted a full statistical analysis for sampling the regolith of Phobos to include the effects of uncertainties in input parameters on the sampling probability. The most likely probability of microbial contamination for return samples is estimated to be two orders of magnitude lower than the criterion defined by the planetary protection policy of the Committee on Space Research (COSPAR). △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: 50 pages, 16 figures, 2 tables, The abstract was shortened to meet the arXiv words limit, accepted for publication in Life Sciences in Space Research, special issues on planetary protection

Assessment of the probability of microbial contamination for sample return from Martian moons I: Departure of microbes from Martian surface

Authors: Kazuhisa Fujita, Kosuke Kurosawa, Hidenori Genda, Ryuki Hyodo, Shingo Matsuyama, Akihiko Yamagishi, Takashi Mikouchi, Takafumi Niihara

Abstract: Potential microbial contamination of Martian moons, Phobos and Deimos, which can be brought about by transportation of Mars ejecta produced by meteoroid impacts on the Martian surface, has been comprehensively assessed in a statistical approach, based on the most probable history of recent major gigantic meteoroid collisions on the Martian surface. This article is the first part of our study to as… ▽ More Potential microbial contamination of Martian moons, Phobos and Deimos, which can be brought about by transportation of Mars ejecta produced by meteoroid impacts on the Martian surface, has been comprehensively assessed in a statistical approach, based on the most probable history of recent major gigantic meteoroid collisions on the Martian surface. This article is the first part of our study to assess potential microbial density in Mars ejecta departing from the Martian atmosphere, as a source of the second part where statistical analysis of microbial contamination probability is conducted. Potential microbial density on the Martian surface as the source of microorganisms was estimated by analogy to the terrestrial areas having the similar arid and cold environments, from which a probabilistic function was deduced as the asymptotic limit. Microbial survival rate during hypervelocity meteoroid collisions was estimated by numerical analysis of impact phenomena with and without taking internal friction and plastic deformation of the colliding meteoroid and the target ground into consideration. Trajectory calculations of departing ejecta through the Martian atmosphere were conducted with taking account of aerodynamic deceleration and heating by the aid of computational fluid dynamic analysis. It is found that Mars ejecta smaller than 0.03 m in diameter hardly reach the Phobos orbit due to aerodynamic deceleration, or mostly sterilized due to significant aerodynamic heating even though they can reach the Phobos orbit and beyond. Finally, the baseline dataset of microbial density in Mars ejecta departing for Martian moons has been presented for the second part of our study. △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: 24 pages, 12 figures, and 5 tables, accepted for publication in Life Sciences in Space Research, special issues on planetary protection

Early formation of moons around large trans-Neptunian objects via giant impacts

Authors: Sota Arakawa, Ryuki Hyodo, Hidenori Genda

Abstract: Recent studies have revealed that all large (over 1000 km in diameter) trans-Neptunian objects (TNOs) form satellite systems. Although the largest Plutonian satellite, Charon, is thought to be an intact fragment of an impactor directly formed via a giant impact, whether giant impacts can explain the variations in secondary-to-primary mass ratios and spin/orbital periods among all large TNOs remain… ▽ More Recent studies have revealed that all large (over 1000 km in diameter) trans-Neptunian objects (TNOs) form satellite systems. Although the largest Plutonian satellite, Charon, is thought to be an intact fragment of an impactor directly formed via a giant impact, whether giant impacts can explain the variations in secondary-to-primary mass ratios and spin/orbital periods among all large TNOs remains to be determined. Here we systematically perform hydrodynamic simulations to investigate satellite formation via giant impacts. We find that the simulated secondary-to-primary mass ratio varies over a wide range, which overlaps with observed mass ratios. We also reveal that the satellite systems' current distribution of spin/orbital periods and small eccentricity can be explained only when their spins and orbits tidally evolve: initially as fluid-like bodies, but finally as rigid bodies. These results suggest that all satellites of large TNOs were formed via giant impacts in the early stage of solar system formation, before the outward migration of Neptune, and that they were fully or partially molten during the giant impact era. △ Less

Submitted 25 June, 2019; originally announced June 2019.

Comments: 24 pages, 8 figures. Accepted for publication in Nature Astronomy. The final authenticated version is available online at: https://rdcu.be/bHBXX

Mars in the aftermath of a colossal impact

Authors: Jason Man Yin Woo, Hidenori Genda, Ramon Brasser, Stephen J. Mojzsis

Abstract: The abundance of highly siderophile elements (HSEs) inferred for Mars' mantle from martian meteorites implies a Late Veneer (LV) mass addition of ~0.8 wt% with broadly chondritic composition. Late accretion to Mars by a differentiated Ceres-sized (~1000 km diameter) object can account for part of the requisite LV mass, and geochronological constraints suggests that this must have occurred no later… ▽ More The abundance of highly siderophile elements (HSEs) inferred for Mars' mantle from martian meteorites implies a Late Veneer (LV) mass addition of ~0.8 wt% with broadly chondritic composition. Late accretion to Mars by a differentiated Ceres-sized (~1000 km diameter) object can account for part of the requisite LV mass, and geochronological constraints suggests that this must have occurred no later than ca. 4480 Ma. Here, we analyze the outcome of the hypothetical LV giant impact to Mars with smoothed particle hydrodynamics simulations together with analytical theory. Results show that, in general about 50% of the impactor's metallic core shatters into ~10m fragments that subsequently fragment into sub-mm metallic hail at re-accretion. This returns a promising delivery of HSEs into martian mantle compared to either a head-on and hit-and-run collision; in both cases,<10% of impactor's core materials are fragmented and finally embedded in the martian mantle. Isotopic evidence from martian meteorites, and interpretations from atmospheric mapping data show that a global surface water reservoir could be present during the early Noachian (before ca. 4100 Ma). The millimeter-sized metal hail could thus react with a martian hydrosphere to generate ~3 bars of H2, which is adequate to act as a greenhouse and keep early Mars warm. Yet, we also find that this atmosphere is transient. It typically survives shorter than 3 Myr based on the expected extreme ultraviolet (EUV) flux of the early Sun; if the Sun was a slow rotator an accordingly weaker EUV flux could extend this lifetime to >10 Myr. A dense pre-Noachian CO2 atmosphere should lower the escape efficiency of hydrogen by IR emission. A more detailed hydrodynamic atmospheric model of this early hydrogen atmosphere is warranted to examine its effect on pre-Noachian Mars. △ Less

Submitted 20 June, 2019; originally announced June 2019.

Journal ref: Icarus 333 (2019) 87-95

Fates of hydrous materials during planetesimal collisions

Authors: Shigeru Wakita, Hidenori Genda

Abstract: Hydrous minerals are found on the surfaces of asteroids, but their origin is not clear. If their origin is endogenic, the hydrous minerals that were formed in the inner part of a planetesimal (or parent body) should come out on to the surface without dehydration. If their origin is exogenic, the source of hydrous minerals accreting onto asteroids is needed. Collisions in the asteroid belt would be… ▽ More Hydrous minerals are found on the surfaces of asteroids, but their origin is not clear. If their origin is endogenic, the hydrous minerals that were formed in the inner part of a planetesimal (or parent body) should come out on to the surface without dehydration. If their origin is exogenic, the source of hydrous minerals accreting onto asteroids is needed. Collisions in the asteroid belt would be related to both origins because collisions excavate the surface and eject the materials. However, the fate of hydrous minerals in large planetesimals during the collisional process has not been well investigated. Here, we explore planetesimal collisions by using the iSALE-2D code, and investigate the effect of an impact for the target planetesimal containing hydrous minerals. Our numerical results for the fiducial case (5 km/s of the impact velocity) show that hydrous minerals are slightly heated during the collisions. This moderate heating indicates that they can avoid the dehydration reaction and keep their original composition. Some hydrous minerals have larger velocity than the escape velocity of the collision system. This means that hydrous minerals can escape from the planetesimal and support the theory of exogenic origin for the hydrous minerals on asteroids. Meanwhile, the velocity of other hydrous minerals is smaller than the escape velocity of the system. This also indicates the possibility of an endogenic origin for the hydrous minerals on asteroids. Our results suggest that hydrous minerals on asteroids can be provided by planetesimal collisions. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 38 pages, 14 figures, accepted for publication in Icarus

Life-Detection Technologies for the Next Two Decades

Authors: Chaitanya Giri, Tony Z. Jia, H. James Cleaves II, Tomohiro Usui, Dhananjay Bodas, Christopher Carr, Huan Chen, Yuka Fujii, Yoshihiro Furukawa, Hidenori Genda, Richard J. Gillams, Keiko Hamano, Shingo Kameda, Ramanarayanan Krishnamurthy, Cornelia Meinert, Markus Meringer, Kishore Paknikar, Sudha Rajamani, Sisinthy Shivaji, Andrew Steele, Masateru Taniguchi, Hikaru Yabuta, Akihiko Yamagishi

Abstract: Since its inception six decades ago, astrobiology has diversified immensely to encompass several scientific questions including the origin and evolution of Terran life, the organic chemical composition of extraterrestrial objects, and the concept of habitability, among others. The detection of life beyond Earth forms the main goal of astrobiology, and a significant one for space exploration in gen… ▽ More Since its inception six decades ago, astrobiology has diversified immensely to encompass several scientific questions including the origin and evolution of Terran life, the organic chemical composition of extraterrestrial objects, and the concept of habitability, among others. The detection of life beyond Earth forms the main goal of astrobiology, and a significant one for space exploration in general. This goal has galvanized and connected with other critical areas of investigation such as the analysis of meteorites and early Earth geological and biological systems, materials gathered by sample-return space missions, laboratory and computer simulations of extraterrestrial and early Earth environmental chemistry, astronomical remote sensing, and in-situ space exploration missions. Lately, scattered efforts are being undertaken towards the R&D of the novel and as-yet-space-unproven life-detection technologies capable of obtaining unambiguous evidence of extraterrestrial life, even if it is significantly different from Terran life. As the suite of space-proven payloads improves in breadth and sensitivity, this is an apt time to examine the progress and future of life-detection technologies. △ Less

Submitted 14 October, 2018; originally announced October 2018.

Comments: 6 pages, the white paper was submitted to and cited by the National Academy of Sciences in support of the Astrobiology Science Strategy for the Search for Life in the Universe

Collisional Disruption of Planetesimals in the Gravity Regime with iSALE Code: Comparison with SPH code for Purely Hydrodynamic Bodies

Authors: Ryo Suetsugu, Hidekazu Tanaka, Hiroshi Kobayashi, Hidenori Genda

Abstract: In most of the previous studies related to collisional disruption of planetesimals in the gravity regime, Smoothed Particle Hydrodynamics (SPH) simulations have been used. On the other hand, impact simulations using grid-based hydrodynamic code have not been sufficiently performed. In the present study, we execute impact simulations in the gravity regime using the shock-physics code iSALE, which i… ▽ More In most of the previous studies related to collisional disruption of planetesimals in the gravity regime, Smoothed Particle Hydrodynamics (SPH) simulations have been used. On the other hand, impact simulations using grid-based hydrodynamic code have not been sufficiently performed. In the present study, we execute impact simulations in the gravity regime using the shock-physics code iSALE, which is a grid-based Eulerian hydrocode. We examine the dependence of the critical specific impact energy Q_RD* on impact conditions for a wide range of specific impact energy (Q_R) from disruptive collisions to erosive collisions, and compare our results with previous studies. We find collision outcomes of the iSALE simulation agree well with those of the SPH simulation. Detailed analysis mainly gives three results. (1) The value of Q_RD* depends on numerical resolution, and is close to convergence with increasing numerical resolution. The difference in converged value of Q_RD* between the iSALE code and the SPH code is within 30%. (2) Ejected mass normalized by total mass (M_ej/M_tot) generally depends on various impact conditions. However, when Q_R is normalized by Q_RD* that is calculated for each impact simulation, M_ej/M_tot can be scaled by Q_R/Q_RD*, and is independent of numerical resolution, impact velocity and target size. (3) This similarity law for Q_R/Q_RD* is confirmed for a wide range of specific impact energy. We also derive a semi-analytic formula for Q_RD* based on the similarity law and the crater scaling law. We find that the semi-analytic formula for the case with a non-porous object is consistent with numerical results. △ Less

Submitted 29 May, 2018; originally announced May 2018.

Comments: Accepted for publication in Icarus

Impact degassing and atmospheric erosion on Venus, Earth, and Mars during the late accretion

Authors: Haruka Sakuraba, Hiroyuki Kurokawa, Hidenori Genda

Abstract: 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 partit… ▽ More 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. △ Less

Submitted 18 May, 2018; originally announced May 2018.

Comments: 31 pages, 12 Figures, accepted for publication in Icarus

On the Impact Origin of Phobos and Deimos IV: Volatile Depletion

Authors: Ryuki Hyodo, Hidenori Genda, Sébastien Charnoz, Francesco C. Pignatale, Pascal Rosenblatt

Abstract: Recent works have shown that Martian moons Phobos and Deimos may have accreted within a giant impact-generated disk whose composition is about an equal mixture of Martian material and impactor material. Just after the giant impact, the Martian surface is heated up to $\sim3000-6000$ K and the building blocks of moons, including volatile-rich vapor, are heated up to $\sim2000$ K. In this paper, we… ▽ More Recent works have shown that Martian moons Phobos and Deimos may have accreted within a giant impact-generated disk whose composition is about an equal mixture of Martian material and impactor material. Just after the giant impact, the Martian surface is heated up to $\sim3000-6000$ K and the building blocks of moons, including volatile-rich vapor, are heated up to $\sim2000$ K. In this paper, we investigate the volatile loss from the building blocks of Phobos and Deimos by hydrodynamic escape of vapor and radiation pressure on condensed particles. We show that a non-negligible amount of volatiles ($> 10\%$ of the vapor with temperature $> 1000$ K via hydrodynamic escape, and moderately volatile dusts that condense at $\sim700-2000$ K via radiation pressure) could be removed just after the impact during their first signle orbit from their pericenters to apocenters. Our results indicate that bulk Phobos and Deimos are depleted in volatile elements. Together with future explorations such as JAXA's MMX (Martian Moons eXploration) mission, our results would be used to constrain the origin of Phobos and Deimos. △ Less

Submitted 19 June, 2018; v1 submitted 24 April, 2018; originally announced April 2018.

Comments: 15 pages, 8 figures. Accepted for publication in ApJ

Implantation of Martian materials in the inner solar system by a mega impact on Mars

Authors: Ryuki Hyodo, Hidenori Genda

Abstract: Observations and meteorites indicate that the Martian materials are enigmatically distributed within the inner solar system. A mega impact on Mars creating a Martian hemispheric dichotomy and the Martian moons can potentially eject Martian materials. A recent work has shown that the mega-impact-induced debris is potentially captured as the Martian Trojans and implanted in the asteroid belt. Howeve… ▽ More Observations and meteorites indicate that the Martian materials are enigmatically distributed within the inner solar system. A mega impact on Mars creating a Martian hemispheric dichotomy and the Martian moons can potentially eject Martian materials. A recent work has shown that the mega-impact-induced debris is potentially captured as the Martian Trojans and implanted in the asteroid belt. However, the amount, distribution, and composition of the debris has not been studied. Here, using hydrodynamic simulations, we report that a large amount of debris ($\sim 1\%$ of Mars' mass), including Martian crust/mantle and the impactor's materials ($\sim 20:80$), are ejected by a dichotomy-forming impact, and distributed between $\sim 0.5-3.0$ astronomical units. Our result indicates that unmelted Martian mantle debris ($\sim 0.02\%$ of Mars' mass) can be the source of Martian Trojans, olivine-rich asteroids in the Hungarian region and the main asteroid belt, and some of which even hit early Earth. A mega impact can naturally implant Martian mantle materials within the inner solar system. △ Less

Submitted 19 March, 2018; originally announced March 2018.

Comments: 8 pages, 5 figures. Accepted for publication in ApJL

Dependence of the onset of the runaway greenhouse effect on the latitudinal surface water distribution of Earth-like planets

Authors: T. Kodama, A. Nitta, H. Genda, Y. Takao, R. O'ishi, A. Abe-Ouchi, Y. Abe

Abstract: Liquid water is one of the most important materials affecting the climate and habitability of a terrestrial planet. Liquid water vaporizes entirely when planets receive insolation above a certain value, which is called the runaway greenhouse threshold. This threshold forms the inner most limit of the habitable zone. Here, we investigate the effect of the distribution of surface water on the runawa… ▽ More Liquid water is one of the most important materials affecting the climate and habitability of a terrestrial planet. Liquid water vaporizes entirely when planets receive insolation above a certain value, which is called the runaway greenhouse threshold. This threshold forms the inner most limit of the habitable zone. Here, we investigate the effect of the distribution of surface water on the runaway greenhouse threshold for Earth-sized planets using a three-dimensional dynamic atmosphere model. We considered a 1-bar atmosphere whose composition is similar to the current Earth's atmosphere with a zonally uniform distribution of surface water. As previous studies have already shown, we also recognized two climate regimes: the land planet regime, which has dry low latitude and wet high latitude region, and the aqua planet regime, which is globally wet. We showed that each regime is controlled by the width of the Hadley circulation, the amount of surface water, and the planetary topography. We found that the runaway greenhouse threshold varies continuously with the surface water distribution from about 130% (an aqua planet) to 180% (the extreme case of a land planet) of the present insolation at Earth's orbit. Our results indicate that the inner edge of the habitable zone is not a single sharp boundary, but a border whose location varies depending on planetary surface condition, such as the amount of surface water. Since land planets have wider habitable zones and less cloud cover, land planets would be good targets for future observations investigating planetary habitability. △ Less

Submitted 22 January, 2018; originally announced January 2018.

Comments: 25 pages, 9 figures. Accepted for publication in Journal of Geophysical Research - Planets

Effects of friction and plastic deformation in shock-comminuted damaged rocks on impact heating

Authors: Kosuke Kurosawa, Hidenori Genda

Abstract: Hypervelocity impacts cause significant heating of planetary bodies. Such events are recorded by a reset of 40Ar-36Ar ages and/or impact melts. Here, we investigate the influence of friction and plastic deformation in shock-generated comminuted rocks on the degree of impact heating using the iSALE shock-physics code. We demonstrate that conversion from kinetic to internal energy in the targets wit… ▽ More Hypervelocity impacts cause significant heating of planetary bodies. Such events are recorded by a reset of 40Ar-36Ar ages and/or impact melts. Here, we investigate the influence of friction and plastic deformation in shock-generated comminuted rocks on the degree of impact heating using the iSALE shock-physics code. We demonstrate that conversion from kinetic to internal energy in the targets with strength occurs during pressure release, and additional heating becomes significant for low-velocity impacts (<10 km/s). This additional heat reduces the impact-velocity thresholds required to heat the targets with the 0.1 projectile mass to temperatures for the onset of Ar loss and melting from 8 and 10 km/s, respectively, for strengthless rocks to 2 and 6 km/s for typical rocks. Our results suggest that the impact conditions required to produce the unique features caused by impact heating span a much wider range than previously thought. △ Less

Submitted 3 January, 2018; originally announced January 2018.

Comments: 16 pages, 4 Figures, 1 Supporting Information, Accepted for publication in Geophysical Research Letters

On the impact origin of Phobos and Deimos III: resulting composition from different impactors

Authors: Francesco C. Pignatale, Sébastien Charnoz, Pascal Rosenblatt, Ryuki Hyodo, Tomoki Nakamura, Hidenori Genda

Abstract: The origin of Phobos and Deimos in a giant impact generated disk is gaining larger attention. Although this scenario has been the subject of many studies, an evaluation of the chemical composition of the Mars' moons in this framework is missing. The chemical composition of Phobos and Deimos is unconstrained. The large uncertainty about the origin of the mid-infrared features, the lack of absorptio… ▽ More The origin of Phobos and Deimos in a giant impact generated disk is gaining larger attention. Although this scenario has been the subject of many studies, an evaluation of the chemical composition of the Mars' moons in this framework is missing. The chemical composition of Phobos and Deimos is unconstrained. The large uncertainty about the origin of the mid-infrared features, the lack of absorption bands in the visible and near-infrared spectra, and the effects of secondary processes on the moons' surface make the determination of their composition very difficult from remote sensing data. Simulations suggest a formation of a disk made of gas and melt with their composition linked to the nature of the impactor and Mars. Using thermodynamic equilibrium we investigate the composition of dust (condensates from gas) and solids (from a cooling melt) that result from different types of Mars impactors (Mars-, CI-, CV-, EH-, comet-like). Our calculations show a wide range of possible chemical compositions and noticeable differences between dust and solids depending on the considered impactors. Assuming Phobos and Deimos as result of the accretion and mixing of dust and solids, we find that the derived assemblage (dust rich in metallic-iron, sulphides and/or carbon, and quenched solids rich in silicates) can be compatible with the observations. The JAXA's MMX (Martian Moons eXploration) mission will investigate the physical and chemical properties of the Maroons, especially sampling from Phobos, before returning to Earth. Our results could be then used to disentangle the origin and chemical composition of the pristine body that hit Mars and suggest guidelines for helping in the analysis of the returned samples. △ Less

Submitted 15 December, 2017; v1 submitted 14 December, 2017; originally announced December 2017.

Comments: Accepted for publication in ApJ (13/12/2017)

On the Impact Origin of Phobos and Deimos II: True Polar Wander and Disk Evolution

Authors: Ryuki Hyodo, Pascal Rosenblatt, Hidenori Genda, Sébastien Charnoz

Abstract: Phobos and Deimos are the two small Martian moons, orbiting almost on the equatorial plane of Mars. Recent works have shown that they can accrete within an impact-generated inner dense and outer light disk, and that the same impact potentially forms the Borealis basin, a large northern hemisphere basin on the current Mars. However, there is no a priori reason for the impact to take place close to… ▽ More Phobos and Deimos are the two small Martian moons, orbiting almost on the equatorial plane of Mars. Recent works have shown that they can accrete within an impact-generated inner dense and outer light disk, and that the same impact potentially forms the Borealis basin, a large northern hemisphere basin on the current Mars. However, there is no a priori reason for the impact to take place close to the north pole (Borealis present location) nor to generate a debris disk in the equatorial plane of Mars (in which Phobos and Deimos orbit). In this paper, we investigate these remaining issues on the giant impact origin of the Martian moons. First, we show that the mass deficit created by the Borealis impact basin induces a global reorientation of the planet to realign its main moment of inertia with the rotation pole (True Polar Wander). This moves the location of the Borealis basin toward its current location. Next, using analytical arguments, we investigate the detailed dynamical evolution of the eccentric inclined disk from the equatorial plane of Mars that is formed by the Martian-moon-forming impact. We find that, as a result of precession of disk particles due to the Martian dynamical flattening $J_{2}$ term of its gravity field and particle-particle inelastic collisions, eccentricity and inclination are damped and an inner dense and outer light equatorial circular disk is eventually formed. Our results strengthen the giant impact origin of Phobos and Deimos that can finally be tested by a future sample return mission such as JAXA's Martian Moons eXploration (MMX) mission. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: 13 pages, 5 figures. Accepted for publication in ApJ

The terrestrial late veneer from core disruption of a lunar-sized impactor

Authors: H. Genda, R. Brasser, S. J. Mojzsis

Abstract: Overabundances in highly siderophile elements (HSEs) of Earth's mantle can be explained by conveyance from a singular, immense (3000 km in a diameter) "Late Veneer" impactor of chondritic composition, subsequent to lunar formation and terrestrial core-closure. Such rocky objects of approximately lunar mass (about 0.01 M_E) ought to be differentiated, such that nearly all of their HSE payload is se… ▽ More Overabundances in highly siderophile elements (HSEs) of Earth's mantle can be explained by conveyance from a singular, immense (3000 km in a diameter) "Late Veneer" impactor of chondritic composition, subsequent to lunar formation and terrestrial core-closure. Such rocky objects of approximately lunar mass (about 0.01 M_E) ought to be differentiated, such that nearly all of their HSE payload is sequestered into iron cores. Here, we analyze the mechanical and chemical fate of the core of such a Late Veneer impactor, and trace how its HSEs are suspended - and thus pollute - the mantle. For the statistically most-likely oblique collision (about 45degree), the impactor's core elongates and thereafter disintegrates into a metallic hail of small particles (about 10 m). Some strike the orbiting Moon as sesquinary impactors, but most re-accrete to Earth as secondaries with further fragmentation. We show that a single oblique impactor provides an adequate amount of HSEs to the primordial terrestrial silicate reservoirs via oxidation of (<m-sized) metal particles with a hydrous, pre-impact, early Hadean Earth. △ Less

Submitted 21 September, 2017; originally announced September 2017.

Comments: Accepted for publication in Earth and Planetary Science Letters, 17 pages, 4 figures

Hydrocode modeling of the spallation process during hypervelocity impacts: Implications for the ejection of Martian meteorites

Authors: Kosuke Kurosawa, Takaya Okamoto, Hidenori Genda

Abstract: Hypervelocity ejection of material by impact spallation is considered a plausible mechanism for material exchange between two planetary bodies. We have modeled the spallation process during vertical impacts over a range of impact velocities from 6 to 21 km/s using both grid- and particle-based hydrocode models. The Tillotson equations of state, which are able to treat the nonlinear dependence of d… ▽ More Hypervelocity ejection of material by impact spallation is considered a plausible mechanism for material exchange between two planetary bodies. We have modeled the spallation process during vertical impacts over a range of impact velocities from 6 to 21 km/s using both grid- and particle-based hydrocode models. The Tillotson equations of state, which are able to treat the nonlinear dependence of density on pressure and thermal pressure in the strongly shocked matter, were used to study the hydrodynamic and thermodynamic response after impacts. The effects of material strength and gravitational acceleration were not considered. A two-dimensional time-dependent pressure field within a 1.5-fold projectile radius from the impact point was investigated in cylindrical coordinates to address the generation of spalled material. A resolution test was also performed to reject ejected materials with peak pressures that were too low due to artificial viscosity. The relationship between ejection velocity veject and peak pressure Ppeak was also derived. Our approach shows that late stage acceleration in an ejecta curtain occurs due to the compressible nature of the ejecta, resulting in an ejection velocity that can be higher than the ideal maximum of the resultant particle velocity after passage of a shock wave. We also calculate the ejecta mass that can escape from a planet like Mars (i.e., veject higher than 5 km/s) that matches the petrographic constraints from Martian meteorites, and which occurs when Ppeak from 30-50 GPa. Although the mass of such ejecta is limited to from 0.1-1 percent of the projectile mass in vertical impacts, this is sufficient for spallation to have been a plausible mechanism for the ejection of Martian meteorites. Finally, we propose that impact spallation is a plausible mechanism for the generation of tektites. △ Less

Submitted 2 September, 2017; originally announced September 2017.

Comments: 67 pages, 28 figures, accepted for publication in Icarus

On the Impact Origin of Phobos and Deimos I: Thermodynamic and Physical Aspects

Authors: Ryuki Hyodo, Hidenori Genda, Sébastien Charnoz, Pascal Rosenblatt

Abstract: Phobos and Deimos are the two small moons of Mars. Recent works have shown that they can accrete within an impact-generated disk. However, the detailed structure and initial thermodynamic properties of the disk are poorly understood. In this paper, we perform high-resolution SPH simulations of the Martian moon-forming giant impact that can also form the Borealis basin. This giant impact heats up t… ▽ More Phobos and Deimos are the two small moons of Mars. Recent works have shown that they can accrete within an impact-generated disk. However, the detailed structure and initial thermodynamic properties of the disk are poorly understood. In this paper, we perform high-resolution SPH simulations of the Martian moon-forming giant impact that can also form the Borealis basin. This giant impact heats up the disk material (around $\sim 2000$ K in temperature) with an entropy increase of $\sim 1500$ J K$^{-1}$ kg$^{-1}$. Thus, the disk material should be mostly molten, though a tiny fraction of disk material ($< 5\%$) would even experience vaporization. Typically, a piece of molten disk material is estimated to be meter sized due to the fragmentation regulated by their shear velocity and surface tension during the impact process. The disk materials initially have highly eccentric orbits ($e \sim 0.6-0.9$) and successive collisions between meter-sized fragments at high impact velocity ($\sim 3-5$ km s$^{-1}$) can grind them down to $\sim100 μ$m-sized particles. On the other hand, a tiny amount of vaporized disk material condenses into $\sim 0.1 μ$m-sized grains. Thus, the building blocks of the Martian moons are expected to be a mixture of these different sized particles from meter-sized down to $\sim 100 μ$m-sized particles and $\sim 0.1 μ$m-sized grains. Our simulations also suggest that the building blocks of Phobos and Deimos contain both impactor and Martian materials (at least 35%), most of which come from the Martian mantle (50-150 km in depth; at least 50%). Our results will give useful information for planning a future sample return mission to Martian moons, such as JAXA's MMX (Martian Moons eXploration) mission. △ Less

Submitted 31 July, 2017; v1 submitted 19 July, 2017; originally announced July 2017.

Comments: 11 pages, 6 figures. Accepted for publication in ApJ

Ejection of iron-bearing giant-impact fragments and the dynamical and geochemical influence of the fragment re-accretion

Authors: Hidenori Genda, Tsuyoshi Iizuka, Takanori Sasaki, Yuichiro Ueno, Masahiro Ikoma

Abstract: The Earth was born in violence. Many giant collisions of protoplanets are thought to have occurred during the terrestrial planet formation. Here we investigated the giant impact stage by using a hybrid code that consistently deals with the orbital evolution of protoplanets around the Sun and the details of processes during giant impacts between two protoplanets. A significant amount of materials (… ▽ More The Earth was born in violence. Many giant collisions of protoplanets are thought to have occurred during the terrestrial planet formation. Here we investigated the giant impact stage by using a hybrid code that consistently deals with the orbital evolution of protoplanets around the Sun and the details of processes during giant impacts between two protoplanets. A significant amount of materials (up to several tens of percent of the total mass of the protoplanets) is ejected by giant impacts. We call these ejected fragments the giant-impact fragments (GIFs). In some of the erosive hit-and-run and high-velocity collisions, metallic iron is also ejected, which comes from the colliding protoplanets' cores. From ten numerical simulations for the giant impact stage, we found that the mass fraction of metallic iron in GIFs ranges from ~ 1wt% to ~ 25wt%. We also discussed the effects of the GIFs on the dynamical and geochemical characteristics of formed terrestrial planets. We found that the GIFs have the potential to solve the following dynamical and geochemical conflicts: (1) The Earth, currently in a near circular orbit, is likely to have had a highly eccentric orbit during the giant impact stage. The GIFs are large enough in total mass to lower the eccentricity of the Earth to its current value via their dynamical friction. (2) The concentrations of highly siderophile elements (HSEs) in the Earth's mantle are greater than what was predicted experimentally. Re-accretion of the iron-bearing GIFs onto the Earth can contribute to the excess of HSEs. In addition, the estimated amount of iron-bearing GIFs provides significant reducing agent that could transform primitive CO2-H2O atmosphere and ocean into more reducing H2-bearing atmosphere. Thus, GIFs are important for the origin of Earth's life and its early evolution. △ Less

Submitted 18 April, 2017; originally announced April 2017.

Comments: Accepted for publication in Earth and Planetary Science Letters, 33 pages, 7 figures, 1 table

Impact Erosion Model for Gravity-Dominated Planetesimals

Authors: H. Genda, T. Fujita, H. Kobayashi, H. Tanaka, R. Suetsugu, Y. Abe

Abstract: Disruptive collisions have been regarded as an important process for planet formation, while non-disruptive, small-scale collisions (hereafter called erosive collisions) have been underestimated or neglected by many studies. However, recent studies have suggested that erosive collisions are also important to the growth of planets, because they are much more frequent than disruptive collisions. Alt… ▽ More Disruptive collisions have been regarded as an important process for planet formation, while non-disruptive, small-scale collisions (hereafter called erosive collisions) have been underestimated or neglected by many studies. However, recent studies have suggested that erosive collisions are also important to the growth of planets, because they are much more frequent than disruptive collisions. Although the thresholds of the specific impact energy for disruptive collisions (Q_RD^*) have been investigated well, there is no reliable model for erosive collisions. In this study, we systematically carried out impact simulations of gravity-dominated planetesimals for a wide range of specific impact energy (Q_R) from disruptive collisions (Q_R ~ Q_RD^*) to erosive ones (Q_R << Q_RD^*) using the smoothed particle hydrodynamics method. We found that the ejected mass normalized by the total mass (M_ej/M_tot) depends on the numerical resolution, the target radius (R_tar) and the impact velocity (v_imp), as well as on Q_R, but that it can be nicely scaled by Q_RD^* for the parameter ranges investigated (R_tar = 30-300 km, v_imp = 2-5 km/s). This means that M_ej/M_tot depends only on Q_R/Q_RD^* in these parameter ranges. We confirmed that the collision outcomes for much less erosive collisions (Q_R < 0.01 Q_RD^*) converge to the results of an impact onto a planar target for various impact angles and that M_ej/M_tot = C * QR/QRD* holds. For disruptive collisions (Q_R ~ Q_RD^*), the curvature of the target has a significant effect on Mej/Mtot. We also examined the angle-averaged value of M_ej/M_tot and found that the numerically obtained relation between angle-averaged M_ej/M_tot and Q_R/Q_RD^* is very similar to the cases for 45-degree impacts. We proposed a new erosion model based on our numerical simulations for future research on planet formation with collisional erosion. △ Less

Submitted 8 March, 2017; originally announced March 2017.

Comments: Accepted for publication in Icarus, 41 pages, 16 figures

Ring Formation around Giant Planets by Tidal Disruption of a Single Passing Large Kuiper Belt Object

Authors: Ryuki Hyodo, Sébastien Charnoz, Keiji Ohtsuki, Hidenori Genda

Abstract: The origin of rings around giant planets remains elusive. Saturn's rings are massive and made of 90-95% of water ice. In contrast, the much less massive rings of Uranus and Neptune are dark and likely to have higher rock fraction. Here we investigate, for the first time, the tidal disruption of a passing object, including the subsequent formation of planetary rings. First, we perform SPH simulatio… ▽ More The origin of rings around giant planets remains elusive. Saturn's rings are massive and made of 90-95% of water ice. In contrast, the much less massive rings of Uranus and Neptune are dark and likely to have higher rock fraction. Here we investigate, for the first time, the tidal disruption of a passing object, including the subsequent formation of planetary rings. First, we perform SPH simulations of the tidal destruction of big differentiated objects ($M_{\rm body}=10^{21-23}$) that experience close encounters with Saturn or Uranus. We find that about $0.1-10$% of the mass of the passing body is gravitationally captured around the planet. However, these fragments are initially big chunks and have highly eccentric orbits around the planet. Then, we perform N-body simulations including the planet's oblateness, starting with data obtained from the SPH simulations. Our N-body simulations show that the chunks are tidally destroyed during their next several orbits. Their individual orbits then start to precess incoherently around the planet's equator, which enhances their encounter velocities on longer-term evolution, resulting in more destructive impacts. These collisions would damp their eccentricities resulting in a progressive collapse of the debris cloud into a thin equatorial and low-eccentricity ring. These high energy impacts are expected to be catastrophic enough to produce small particles. Our numerical results also show that the mass of formed rings is large enough to explain current rings including inner regular satellites around Saturn and Uranus. In the case of Uranus, a body can go deeper inside the planet's Roche limit resulting in a more efficient capture of rocky material compared to Saturn's case in which mostly ice is captured. Thus, our results can naturally explain the compositional difference between the rings of Saturn, Uranus and Neptune. △ Less

Submitted 8 September, 2016; originally announced September 2016.

Comments: 33 pages, 21 figures; accepted in Icarus

Formation of Centaurs' rings through their partial tidal disruption during planetary encounters

Authors: Ryuki Hyodo, Sébastien Charnoz, Hidenori Genda, Keiji Ohtsuki

Abstract: Centaurs are minor planets orbiting between Jupiter and Neptune that have or had crossing orbits with one or more giant planets. Recent observations and reinterpretation of previous observations have revealed the existence of ring systems around 10199 Chariklo and 2060 Chiron. However, the origin of the ring systems around such a minor planet is still an open question. Here, we propose that the ti… ▽ More Centaurs are minor planets orbiting between Jupiter and Neptune that have or had crossing orbits with one or more giant planets. Recent observations and reinterpretation of previous observations have revealed the existence of ring systems around 10199 Chariklo and 2060 Chiron. However, the origin of the ring systems around such a minor planet is still an open question. Here, we propose that the tidal disruption of a differentiated object that experiences a close encounter with a giant planet could naturally form diverse ring-satellite systems around the Centaurs. During the close encounter, the icy mantle of the passing object is preferentially ripped off by the planet's tidal force and the debris is distributed mostly within the Roche limit of the largest remnant body. Assuming the existence of $20-50$wt% silicate core below the icy mantle, a disk of particles is formed when the objects pass within $0.4-0.8$ of the planet's Roche limit with the relative velocity at infinity $3-6$km s$^{-1}$ and 8h initial spin period of the body. The resultant ring mass is $0.1-10$% of the central object's mass. Such particle disks are expected to spread radially, and materials spreading beyond the Roche limit would accrete into satellite(s). Our numerical results suggest that ring formation would be a natural outcome of such extreme close encounters and Centaurs can naturally have such ring systems because they cross the orbits of the giant planets. △ Less

Submitted 1 September, 2016; v1 submitted 11 August, 2016; originally announced August 2016.

Comments: 8 pages, 3 figures, published in ApJL

Journal ref: The Astrophysical Journal Letters 2016, Volume 828, Number 1

The Giant Impact Simulations with Density Independent Smoothed Particle Hydrodynamics

Authors: Natsuki Hosono, Takayuki R. Saitoh, Junichiro Makino, Hidenori Genda, Shigeru Ida

Abstract: At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free s… ▽ More At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free surface such as a planetary surface. This difficulty comes from the assumption of differentiability of density in SSPH. We have developed an alternative formulation of SPH, density independent SPH (DISPH), which is based on differentiability of pressure instead of density to solve the problem of a contact discontinuity. In this paper, we report the results of the GI simulations with DISPH and compare them with those obtained with SSPH. We found that the disk properties, such as mass and angular momentum produced by DISPH is different from that of SSPH. In general, the disks formed by DISPH are more compact: while formation of a smaller mass moon for low-oblique impacts is expected with DISPH, inhibition of ejection would promote formation of a larger mass moon for high-oblique impacts. Since only the improvement of core-mantle boundary significantly affects the properties of circumplanetary disks generated by GI and DISPH has not been significantly improved from SSPH for a free surface, we should be very careful when some conclusions are drawn from the numerical simulations for GI. And it is necessary to develop the numerical hydrodynamical scheme for GI that can properly treat the free surface as well as the contact discontinuity. △ Less

Submitted 2 February, 2016; originally announced February 2016.

Comments: Accepted for publication in Icarus

Rapid Water Loss can Extend the Lifetime of the Planetary Habitability

Authors: T. Kodama, H. Genda, Y. Abe, K. J. Zahnle

Abstract: Two habitable planetary states are proposed: an aqua planet like the Earth and a land planet that has a small amount of water. Land planets keep liquid water under larger solar radiation compared to aqua planets. Water loss may change an aqua planet into a land planet, and the planet can remain habitable for a longer time than had it stayed an aqua planet. We calculate planetary evolution with hyd… ▽ More Two habitable planetary states are proposed: an aqua planet like the Earth and a land planet that has a small amount of water. Land planets keep liquid water under larger solar radiation compared to aqua planets. Water loss may change an aqua planet into a land planet, and the planet can remain habitable for a longer time than had it stayed an aqua planet. We calculate planetary evolution with hydrogen escape for different initial water inventories and different distances from the central star. We find that there are two conditions to evolve an aqua planet into a land planet: the critical amount of water on the surface M_ml consistent with a planet being a land planet, and the critical amount of water vapor in the atmosphere M_cv that defines the onset of the runaway greenhouse state. We find that Earth-size aqua planets with initial oceans < 10 % of the Earth's can evolve into land planets if M_cv = 3 m in precipitable water and M_ml = 5 % of the Earth's ocean mass. Such planets can keep liquid water on their surface for another 2 Gyrs. The initial amount of water and M_cv are shown to be important dividing parameters of the planetary evolution path. Our results indicate that massive hydrogen escape could give a fresh start as another kind of habitable planet rather than the end of its habitability. △ Less

Submitted 12 September, 2015; originally announced September 2015.

Comments: Accepted for publication in The Astrophysical Journal, 31 pages, 11 figures

Resolution Dependence of Disruptive Collisions between Planetesimals in the Gravity Regime

Authors: H. Genda, T. Fujita, H. Kobayashi, H. Tanaka, Y. Abe

Abstract: Collisions are a fundamental process in planet formation. If colliding objects simply merge, a planetary object can grow. However, if the collision is disruptive, planetary growth is prevented. Therefore, the impact conditions under which collisions are destructive are important in understanding planet formation. So far, the critical specific impact energy for a disruptive collision Q_D^* has been… ▽ More Collisions are a fundamental process in planet formation. If colliding objects simply merge, a planetary object can grow. However, if the collision is disruptive, planetary growth is prevented. Therefore, the impact conditions under which collisions are destructive are important in understanding planet formation. So far, the critical specific impact energy for a disruptive collision Q_D^* has been investigated for various types of collisions between objects ranging in scale from centimeters to thousands of kilometers. Although the values of Q_D^* have been calculated numerically while taking into consideration various physical properties such as self-gravity, material strength, and porosity, the dependence of Q_D^* on numerical resolution has not been sufficiently investigated. In this paper, using the smoothed particle hydrodynamics (SPH) method, we performed numerical simulations of collisions between planetesimals at various numerical resolutions (from 2 x 10^4 to 5 x 10^6 SPH particles) and investigated the resulting variation in Q_D^*. The value of Q_D^* is shown to decrease as the number of SPH particles increases, and the difference between the Q_D^* values for the lowest and highest investigated resolutions is approximately a factor of two. Although the results for 5 x 10^6 SPH particles do not fully converge, higher-resolution simulations near the impact site show that the value of Q_D^* for the case with 5 x 10^6 SPH particles is close to the expected converged value. Although Q_D^* depends on impact parameters and material parameters, our results indicate that at least 5 x 10^6 SPH particles are required for numerical simulations in disruptive collisions to obtain the value of Q_D^* within 20% error. △ Less

Submitted 22 August, 2015; originally announced August 2015.

Comments: Accepted for publication in Icarus, 29 pages, 11 figures

Warm Debris Disks Produced by Giant Impacts During Terrestrial Planet Formation

Authors: H. Genda, H. Kobayashi, E. Kokubo

Abstract: In our solar system, Mars-sized protoplanets frequently collided with each other during the last stage of terrestrial planet formation called the giant impact stage. Giant impacts eject a large amount of material from the colliding protoplanets into the terrestrial planet region, which may form debris disks with observable infrared excesses. Indeed, tens of warm debris disks around young solar-typ… ▽ More In our solar system, Mars-sized protoplanets frequently collided with each other during the last stage of terrestrial planet formation called the giant impact stage. Giant impacts eject a large amount of material from the colliding protoplanets into the terrestrial planet region, which may form debris disks with observable infrared excesses. Indeed, tens of warm debris disks around young solar-type stars have been observed. Here, we quantitatively estimate the total mass of ejected materials during the giant impact stages. We found that $\sim$0.4 times the Earth's mass is ejected in total throughout the giant impact stage. Ejected materials are ground down by collisional cascade until micron-sized grains are blown out by radiation pressure. The depletion timescale of these ejected materials is determined primarily by the mass of the largest body among them. We conducted high-resolution simulations of giant impacts to accurately obtain the mass of the largest ejected body. We then calculated the evolution of the debris disks produced by a series of giant impacts and depleted by collisional cascades to obtain the infrared excess evolution of the debris disks. We found that the infrared excess is almost always higher than the stellar infrared flux throughout the giant impact stage ($\sim$100 Myr) and is sometimes $\sim$10 times higher immediately after a giant impact. Therefore, giant impact stages would explain the infrared excess from most observed warm debris disks. The observed fraction of stars with warm debris disks indicates that the formation probability of our solar system-like terrestrial planets is approximately 10%. △ Less

Submitted 5 August, 2015; originally announced August 2015.

Comments: Accepted to ApJ

Formation of Phobos and Deimos via a Giant Impact

Authors: Robert Citron, Hidenori Genda, Shigeru Ida

Abstract: Although the two moons of Mars, Phobos and Deimos, have long been thought to be captured asteroids, recent observations of their compositions and orbits suggest that they may have formed from debris generated by one or more giant impacts of bodies with ~ 0.01 x target mass. Recent studies have both analytically estimated debris produced by giant impacts on Mars and numerically examined the evoluti… ▽ More Although the two moons of Mars, Phobos and Deimos, have long been thought to be captured asteroids, recent observations of their compositions and orbits suggest that they may have formed from debris generated by one or more giant impacts of bodies with ~ 0.01 x target mass. Recent studies have both analytically estimated debris produced by giant impacts on Mars and numerically examined the evolution of circum-Mars debris disks. We perform a numerical study (Smoothed Particle Hydrodynamics simulation) of debris retention from giant impacts onto Mars, particularly in relation to a Borealis-scale giant impact (E ~ 3 x 10^29 J) capable of producing the Borealis basin. We find that a Borealis-scale impact is capable of producing a disk of mass ~ 5 x 10^20 kg (~ 1 - 4 % of the impactor mass), sufficient debris to form at least one of the martian moons according to recent numerical studies of martian debris disk evolution. While a Borealis-scale impact may generate sufficient debris to form both Phobos and Deimos, further studies of the debris disk evolution are necessary. Our results can serve as inputs for future studies of martian debris disk evolution. △ Less

Submitted 18 March, 2015; originally announced March 2015.

Comments: Accepted to Icarus

Journal ref: Icarus 252 (2015) 334-338

Merging Criteria for Giant Impacts of Protoplanets

Authors: H. Genda, E. Kokubo, S. Ida

Abstract: At the final stage of terrestrial planet formation, known as the giant impact stage, a few tens of Mars-sized protoplanets collide with one another to form terrestrial planets. Almost all previous studies on the orbital and accretional evolution of protoplanets in this stage have been based on the assumption of perfect accretion, where two colliding protoplanets always merge. However, recent impac… ▽ More At the final stage of terrestrial planet formation, known as the giant impact stage, a few tens of Mars-sized protoplanets collide with one another to form terrestrial planets. Almost all previous studies on the orbital and accretional evolution of protoplanets in this stage have been based on the assumption of perfect accretion, where two colliding protoplanets always merge. However, recent impact simulations have shown that collisions among protoplanets are not always merging events, that is, two colliding protoplanets sometimes move apart after the collision (hit-and-run collision). As a first step towards studying the effects of such imperfect accretion of protoplanets on terrestrial planet formation, we investigated the merging criteria for collisions of rocky protoplanets. Using the smoothed particle hydrodynamic (SPH) method, we performed more than 1000 simulations of giant impacts with various parameter sets, such as the mass ratio of protoplanets, $γ$, the total mass of two protoplanets, $M_{\rm T}$, the impact angle, $θ$, and the impact velocity, $v_{\rm imp}$. We investigated the critical impact velocity, $v_{\rm cr}$, at the transition between merging and hit-and-run collisions. We found that the normalized critical impact velocity, $v_{\rm cr}/v_{\rm esc}$, depends on $γ$ and $θ$, but does not depend on $M_{\rm T}$, where $v_{\rm esc}$ is the two-body escape velocity. We derived a simple formula for $v_{\rm cr}/v_{\rm esc}$ as a function of $γ$ and $θ$, and applied it to the giant impact events obtained by \textit{N}-body calculations in the previous studies. We found that 40% of these events should not be merging events. △ Less

Submitted 16 September, 2011; originally announced September 2011.

Comments: Accepted to ApJ

Giant impacts in the Saturnian System: a possible origin of diversity in the inner mid-sized satellites

Authors: Yasuhito Sekine, Hidenori Genda

Abstract: It is widely accepted that Titan and the mid-sized regular satellites around Saturn were formed in the circum-Saturn disk. Thus, if these mid-sized satellites were simply accreted by collisions of similar ice-rock satellitesimals in the disk, the observed wide diversity in density (i.e., the rock fraction) of the Saturnian mid-sized satellites is enigmatic. A recent circumplanetary disk model sugg… ▽ More It is widely accepted that Titan and the mid-sized regular satellites around Saturn were formed in the circum-Saturn disk. Thus, if these mid-sized satellites were simply accreted by collisions of similar ice-rock satellitesimals in the disk, the observed wide diversity in density (i.e., the rock fraction) of the Saturnian mid-sized satellites is enigmatic. A recent circumplanetary disk model suggests satellite growth in an actively supplied circumplanetary disk, in which Titan-sized satellites migrate inward by interaction with the gas and are eventually lost to the gas planet. Here we report numerical simulations of giant impacts between Titan-sized migrating satellites and smaller satellites in the inner region of the Saturnian disk. Our results suggest that in a giant impact with impact velocity > 1.4 times the escape velocity and impact angle of ~45 degree, a smaller satellite is destroyed, forming multiple mid-sized satellites with a very wide diversity in satellite density (the rock fraction = 0-92 wt%). Our results of the relationship between the mass and rock fraction of the satellites resulting from giant impacts reproduce the observations of the Saturnian mid-sized satellites. Giant impacts also lead to internal melting of the formed mid-sized satellites, which would initiate strong tidal dissipation and geological activity, such as those observed on Enceladus today and Tethys in the past. Our findings also imply that giant impacts might have affected the fundamental physical property of the Saturnian mid-sized satellites as well as those of the terrestrial planets in the solar system and beyond. △ Less

Submitted 20 June, 2011; originally announced June 2011.

Comments: 18 pages, 3 figures, Planetary and Space Science, in press