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

Metal pollution in white dwarf photospheres originates from the accretion of some combination of planets, moons, asteroids, comets, boulders, pebbles and dust. When large bodies reside in dynamically stagnant locations -- unable themselves to pollute nor even closely approach the white dwarf -- then smaller reservoirs of impact debris may become a complementary or the primary source of metal pollutants. Here, we take a first step towards exploring this possibility by computing limits on the recoil mass that escapes the gravitational pull of the target object following a single impact onto an atmosphere-less surface. By considering vertical impacts only with the full-chain analytical prescription from Kurosawa & Takada (2019), we provide lower bounds for the ejected mass for basalt, granite, iron and water-rich target objects across the radii range 10^{0-3} km. Our use of the full-chain prescription as opposed to physical experiments or hydrocode simulations allows us to quickly sample a wide range of parameter space appropriate to white dwarf planetary systems. Our numerical results could be used in future studies to constrain freshly-generated small debris reservoirs around white dwarfs given a particular planetary system architecture, bombardment history, and impact geometries.

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.

We investigate the patterns observed in ejecta curtain induced by hypervelocity impact (2-6 km/s) with a variety of the size and shape of target particles. We characterize the patterns by an angle, defined as the ratio of the characteristic length of the pattern obtained by Fourier transformation to the distance from the impact point. This angle is found to be almost the same as that obtained by the reanalysis of the patterns in the previous study at lower impact velocities (Kadono et al., 2015, Icarus 250, 215-221), which are consistent with lunar crater-ray systems. Assuming that the pattern is formed by mutual collision of particles with fluctuation velocity in excavation flow, we evaluate an angle at which the pattern growth stops and show that this angle is the same in the order of magnitude as the ratio of the fluctuation velocity and the radial velocity. This relation is confirmed in the results of experiments and numerical simulations. Finally, we discuss the dependence of the patterns on impact conditions. The experiments show no dependence of the angle on impact velocity. This indicates that the ratio between the fluctuation and radial velocity components in excavation flow does not depend on impact velocity. Moreover, the independences on particle size and particle shape suggest that the angle characterizing the structure of the patterns does not depend on cohesive force. Since cohesive forces should be related with elastic properties of particles, the structure does not depend on elastic properties, though inelastic collisions are important for the persistence and contrast of the patterns.

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.

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

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.

Dry lakebeds might constitute large volatile reservoirs on Mars. Hypervelocity impacts onto ancient dry lakebeds would have affected the volatile distribution on Mars. We developed a new experimental method to investigate the response of evaporitic minerals (halite and gypsum) to impact shocks in an open system. This technique does not result in chemical contamination from the operation of the gas gun. The technique is termed the two-valve method and the gun system is located in the Planetary Exploration Research Center, Chiba Institute of Technology, Japan. We detected the vaporization of halite at 31 GPa and devolatilization from gypsum at 11 GPa, suggesting that impact-induced volatile release from dry lakebeds has periodically occurred throughout Martian history. The vaporization of halite deposits might have enhanced the production of perchlorates, which are found globally on Mars. The water loss from gypsum possibly explains the coexisting types of Ca-sulfates found in Gale Crater.

Impact craters are among the most prominent topographic features on planetary bodies. Crater scaling laws allow us to extract information about the impact histories on the host bodies. The pi-group scaling laws have been constructed based on the point-source approximation, dimensional analysis, and the results from impact experiments. Recent impact experiments, however, demonstrated that the scaling parameters themselves exhibits complex behavior against the change in the impact conditions and target properties. Here, we propose an alternative, fully analytical method to predict impact outcomes based on impact cratering mechanics. This approach is based on the Maxwell Z-model and the residual velocity. We present analytical expressions of (1) the proportionality relation between the ejection velocity and the ejection position, (2) the radius of a growing crater as a function of time, and (3) the transient crater radii. Since we focused on obtaining analytical solutions, a number of simplifications are employed. Due to the simplifications in the strength model, the accuracy of the prediction in the strength-dominated regime is relatively low. Our model reproduces the power-law behavior of the ejecta velocity distribution and the approximate time variation of a growing crater. The predicted radii under typical impact conditions mostly converge to a region between the two typical scaling lines by the pi-group scaling laws, strongly supporting the notion that the new method is one of the simplest ways to predict impact outcomes, as it provides analytical solutions. Our model could serve as a quick-look tool to estimate the impact outcome under a given set of conditions, and it might provide new insights into the nature of impact excavation processes.

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.

Impact-cratering processes on small bodies are thought to be mainly controlled by the local material strength because of their low surface gravity, and craters that are as large as the parent bodies should be affected by the target curvature. Although cratering processes on planar surfaces in the strength-controlled regime have been studied extensively, the mechanism by which target curvature affects the cratering processes remains unclear. Herein, we report on a series of impact experiments that used spherical targets with various diameters. The resultant craters consisted of a deep circular pit and an irregular-shaped spall region around the pit, which is consistent with the features reported in a number of previous cratering experiments on planar surfaces. However, the volume and radius of the craters increased with the normalized curvature. The results indicate that the increase in the spall-region volume and radius mainly contributes to the increase in the whole crater volume and radius, although the volume, depth, and radius of pits remain constant with curvature. The results of our model indicate that the geometric effect due to curvature (i.e., whereby the distance from the equivalent center to the target free surface is shorter for higher curvature values) contributes to increases in the cratering efficiency. Our results suggest that the impactors that produce the largest craters (basins) on some asteroids are thus smaller than what is estimated by current scaling laws, which do not take into account the curvature effects.

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.

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

We carried out laser shock experiments and wholly recovered shocked olivine and quartz samples. We investigated the petrographic features based on optical micrographs of sliced samples and found that each recovered sample comprises three regions, I (optically dark), II (opaque) and III (transparent). Scanning electron microscopy combined with electron back-scattered diffraction shows that there are no crystal features in the region I; the materials in the region I have once melted. Moreover, numerical calculations performed with the iSALE shock physics code suggest that the boundary between regions II and III corresponds to Hugoniot Elastic Limit (HEL). Thus, we succeeded in the recovery of the entire shocked samples experienced over a wide range of pressures from HEL (~ 10 GPa) to melting pressure (~ 100 GPa) in a hierarchical order.

The fate of surface water on Venus is one of the most important outstanding problems in comparative planetology. Here a new concept is proposed to explain water removal on a steam-covered proto Venus, referred to as impact-driven planetary desiccation. Since a steam atmosphere is photochemically unstable, water vapor dissociates into hydrogen and oxygen. Then, hydrogen escapes easily into space through hydrodynamic escape driven by strong extreme ultraviolet radiation from the young Sun. The focus is on the intense impact bombardment during the terminal stage of planetary accretion as generators of a significant amount of reducing agent. The fine-grained ejecta remove the residual oxygen, the counter part of escaped hydrogen, via the oxidation of iron-bearing rocks in a hot atmosphere. Thus, hypervelocity impacts cause net desiccation of the planetary surface. I constructed a stochastic cratering model using a Monte Carlo approach to investigate the cumulative mass of nonoxidized, ejected rocks due to the intense impact bombardment. It is shown that a thick steam atmosphere with a mass equivalent to that of the terrestrial oceans would be removed. The cumulative mass of rocky ejecta released into the atmosphere reaches 1 wt% of the host planet, which is 10000 times of the current mass of the Earths atmosphere. These results strongly suggest that chemical reactions between such large amounts of ejecta and planetary atmospheres are among the key factors required to understand atmospheric mass and its composition, not only in the Solar System but also in extrasolar systems.

The bulk chemical compositions of planets are uncertain, even for major elements such as Mg and Si. This is due to the fact that the samples available for study all originate from relatively shallow depths. Comparison of the stable isotope compositions of planets and meteorites can help overcome this limitation. Specifically, the non-chondritic Si isotope composition of the Earth's mantle was interpreted to reflect the presence of Si in the core, which can also explain its low density relative to pure Fe-Ni alloy. However, we have found that angrite meteorites display a heavy Si isotope composition similar to the lunar and terrestrial mantles. Because core formation in the angrite parent-body (APB) occurred under oxidizing conditions at relatively low pressure and temperature, significant incorporation of Si in the core is ruled out as an explanation for this heavy Si isotope signature. Instead, we show that equilibrium isotopic fractionation between gaseous SiO and solid forsterite at 1370 K in the solar nebula could have produced the observed Si isotope variations. Nebular fractionation of forsterite should be accompanied by correlated variations between the Si isotopic composition and Mg/Si ratio following a slope of 1, which is observed in meteorites. Consideration of this nebular process leads to a revised Si concentration in the Earth's core of 3.6 (+6.0/-3.6) wt% and provides estimates of Mg/Si ratios of bulk planetary bodies.