Atmospheric erosion and replenishment induced by impacts of cosmic bodies upon the Earth and Mars

The atmospheric erosion induced by impacts of cosmic bodies with sizes from ～100 m to 10 km is calculated for the Earth with its present atmosphere and for Mars with a dense carbon dioxide atmosphere that could be at the early stages of planetary evolution. Numerical results are compared to simple analytic models and calculations performed by other authors; approximate formulas are suggested. The evolutions of early atmospheres, which could exist at the late stage of the planetary accumulation, are numerically simulated using an integral model of impact-induced atmospheric erosion and replenishment in the approximation of a one-component atmosphere with a composition determined by the basic atmosphile component of the bodies falling onto the planet.     

Atmospheric erosion and replenishment induced by impacts upon the Earth and Mars during a heavy bombardment

Consequences of a heavy bombardment for the atmospheres of Earth and Mars are investigated with a stochastic model. The main result is the dominance of the accumulation. The atmospheric pressure is strongly increasing both for Earth and Mars in the course of an enhanced bombardment. The effect of atmospheric erosion is found to be minor, regarding escape during meteorite entry, in the expanding vapor plume, and ejection due to free-surface motion. The initial atmospheric surface pressure if comparable to the modern value turns out as a less important additive constant of the final pressure. Impactor retention and atmospheric erosion are parametrized in terms of scaling laws, compatible with recent numerical simulations. The dependence on impactor size, atmospheric and planetary parameters is analyzed among alternative models and numerical results. The stochastic model is fed with the net replenishment originating from impactor material and the loss of preexisting atmospheric gas. Major input parameters are the total cumulative impactor mass and the relative mass of atmophile molecules in comets and asteroids. Input size distributions of the impactor ensemble correspond to presently observed main belt asteroids and KBOs. Velocity distributions are taken from dynamical simulations for the Nice model. Depending on the composition of large cometary impactors, the Earth could acquire a more massive atmosphere, a few bars in terms of surface pressure, mostly as CO and CO₂. For Mars accumulation of 1–4 bars of CO and CO₂ requires an asteroidal ‘late veneer’ of the order of 10²⁴ g containing 2% atmophiles.     

Projectile preservation during oblique hypervelocity impacts

Impact angle plays a significant role in determining the fate of the projectile. In this study, we use a suite of hypervelocity impact experiments to reveal how impact angle affects the preservation, distribution, and physical state of projectile residues in impact craters. Diverse types of projectiles, including amorphous silicates, crystalline silicates, and aluminum, in two sizes (6.35 and 12.7 mm), were launched into blocks of copper or 6061 aluminum at speeds between 1.9 and 5.7 km s⁻¹. Crater interiors preserve projectile residues in all cases, including conditions relevant to the asteroid belt. These residues consist of projectile fragments or projectile-rich glasses, depending on impact conditions. During oblique impacts at 30° and 45°, the uprange crater wall preserves crystalline fragments of the projectile. The fragments of water-rich projectiles such as antigorite remain hydrated. Several factors contribute to enhanced preservation on the uprange wall, including a weaker shock uprange, uprange acceleration as the shock reflects off the back of the projectile, and rapid quenching of melts along the projectile–target interface. These findings have two broader implications. First, the results suggest a new collection strategy for flyby sample return missions. Second, these results predict that the M-type asteroid Psyche should bear exogenic, impactor-derived debris.     

Launch of martian meteorites in oblique impacts

A high-velocity oblique impact into the martian surface accelerates solid target material to escape velocity. A fraction of that material eventually falls as meteorites on Earth. For a long time they were called the SNC meteorites (Shergotty, Nakhla, and Chassigny). We study production of potential martian meteorites numerically within the frame of 3D hydrodynamic modeling. The ratio of the volume of escaping solid ejecta to projectile volume depends on the impact angle, impact velocity and the volatile content in the projectile and in the target. The size distribution of ejected fragments appears to be of crucial importance for the atmosphere-ejecta interaction in the case of a relatively small impact (with final crater size <3 km): 10-cm-sized particles are decelerated efficiently, while 30-50% of larger fragments could escape Mars. The results of numerical modeling are compared with shock metamorphic features in martian meteorites, their burial depth, and preatmospheric mass. Although it is impossible to accelerate ejected fragments to escape velocity without substantial compression (above 10 GPa), the maximum temperature increase in dunite (Chassigny) or ortopyroxenite (ALH84001) may be lower than 200 degree. This result is consistent with the observed chaotic magnetization of ALH84001. The probability of microbes' survival may be rather high even for the extreme conditions during the ejection process.     

Differential melt scaling for oblique impacts on terrestrial planets

Analytical estimates of melt volumes produced by a given projectile and contained in a given impact crater are derived as a function of impact velocity, impact angle, planetary gravity, target and projectile densities, and specific internal energy of melting. Applications to impact events and impact craters on the Earth, Moon, and Mars are demonstrated and discussed. The most probable oblique impact (45°) produces ～1.6 times less melt volume than a vertical impact, and ～1.6 and 3.7 times more melt volume than impacts with 30° and 15° trajectories, respectively. The melt volume for a particular crater diameter increases with planetary gravity, so a crater on Earth should have more melt than similar-size craters on Mars and the Moon. The melt volume for a particular projectile diameter does not depend on gravity, but has a strong dependence on impact velocity, so the melt generated by a given projectile on the Moon is significantly larger than on Mars. Higher surface temperatures and geothermal gradients increase melt production, as do lower energies of melting. Collectively, the results imply thinner central melt sheets and a smaller proportion of melt particles in impact breccias on the Moon and Mars than on Earth. These effects are illustrated in a comparison of the Chicxulub crater on Earth, linked to the Cretaceous–Tertiary mass extinction, Gusev crater on Mars, where the Mars Exploration Rover Spirit landed, and Tsiolkovsky crater on the Moon. The results are comparable to those obtained from field and spacecraft observations, other analytical expressions, and hydrocode simulations.     

Numerical modeling of oblique hypervelocity impacts on strong ductile targets

Abstract– The majority of meteorite impacts occur at oblique incidence angles. However, many of the effects of obliquity on impact crater size and morphology are poorly understood. Laboratory experiments and numerical models have shown that crater size decreases with impact angle, the along‐range crater profile becomes asymmetric at low incidence angles, and below a certain threshold angle the crater planform becomes elliptical. Experimental results at approximately constant impact velocity suggest that the elliptical threshold angle depends on target material properties. Herein, we test the hypothesis that the threshold for oblique crater asymmetry depends on target material strength. Three‐dimensional numerical modeling offers a unique opportunity to study the individual effects of both impact angle and target strength; however, a systematic study of these two parameters has not previously been performed. In this work, the three‐dimensional shock physics code iSALE‐3D is validated against laboratory experiments of impacts into a strong, ductile target material. Digital elevation models of craters formed in laboratory experiments were created from stereo pairs of scanning electron microscope images, allowing the size and morphology to be directly compared with the iSALE‐3D craters. The simulated craters show excellent agreement with both the crater size and morphology of the laboratory experiments. iSALE‐3D is also used to investigate the effect of target strength on oblique incidence impact cratering. We find that the elliptical threshold angle decreases with decreasing target strength, and hence with increasing cratering efficiency. Our simulations of impacts on ductile targets also support the prediction from Chapman and McKinnon (1986) that cratering efficiency depends on only the vertical component of the velocity vector.     

Hydrocode modeling of oblique impacts: The fate of the projectile

Abstract— All impacts are oblique to some degree. Only rarely do projectiles strike a planetary surface (near) vertically. The effects of an oblique impact event on the target are well known, producing craters that appear circular even for low impact angles (>15° with respect to the surface). However, we still have much to learn about the fate of the projectile, especially in oblique impact events. This work investigates the effect of angle of impact on the projectile. Sandia National Laboratories' three‐dimensional hydrocode CTH was used for a series of high‐resolution simulations (50 cells per projectile radius) with varying angle of impact. Simulations were carried out for impacts at 90, 60, 45, 30, and 15° from the horizontal, while keeping projectile size (5 km in radius), type (dunite), and impact velocity (20 km/s) constant. The three‐dimensional hydrocode simulations presented here show that in oblique impacts the distribution of shock pressure inside the projectile (and in the target as well) is highly complex, possessing only bilateral symmetry, even for a spherical projectile. Available experimental data suggest that only the vertical component of the impact velocity plays a role in an impact. If this were correct, simple theoretical considerations indicate that shock pressure, temperature, and energy would depend on sin²θ, where θ is the angle of impact (measured from the horizontal). However, our numerical simulations show that the mean shock pressure in the projectile is better fit by a sin θ dependence, whereas shock temperature and energy depend on sin^3/2 θ. This demonstrates that in impact events the shock wave is the result of complex processes that cannot be described by simple empirical rules. The mass of shock melt or vapor in the projectile decreases drastically for low impact angles as a result of the weakening of the shock for decreasing impact angles. In particular, for asteroidal impacts the amount of projectile vaporized is always limited to a small fraction of the projectile mass. In cometary impacts, however, most of the projectile is vaporized even at low impact angles. In the oblique impact simulations a large fraction of the projectile material retains a net downrange motion. In agreement with experimental work, the simulations show that for low impact angles (30 and 15°), a downrange focusing of projectile material occurs, and a significant amount of it travels at velocities larger than the escape velocity of Earth.     

Ejecta deposition after oblique impacts: An influence of impact scale

Abstract– Simple estimates suggest that ejecta blankets around larger craters should be more asymmetric than around smaller craters for the same oblique impact angle. Numerical simulations presented in the paper confirm that an increase in the scale of gravity‐dominated craters (and in the size of the corresponding projectiles) increases the asymmetry of both impact craters and ejecta blankets around them.     

Laboratory simulation of volatilisation from melts induced by micrometeoroid impacts

The laboratory simulation experiments on volatilization from the melts formed due to micrometeoroid impacts on the lunar surface were carried out. The simulation was performed using modulated laser pulses incident on rocks and minerals in vacuum; in so doing microcraters on the targets, glass particles, condensates were formed and gases solved in the bulk of the rock and mineral grains were released. It is shown that volatilization from only the crater glass layers is negligible, which fully confirms the theoretical predictions of Nussinov and Chernyak.The volatilizations from the drops formed by the micrometeoric impacts can be noticeable. For micron-sized drops, predominant among the others after the impact, the volatilization of Na, K and other volatiles can be up to 10 to 20%. For larger-sized (d10² m) drops the volatilization should lead to the appearance of the significant gradients of the element concentrations over the particle cross-section and as a result of their normalization it leads to the surface enrichment by some elements (Si and others).The mechanism of dust formation due to the surface rocks attack by volatilizing alkaline metals (Naughtonet al.) on the Moon probably is not effective. It is the consequence of such a fact that the condensate of the target materials evaporated due to other simultaneous micrometeoric impacts and had the same (as the target has) elemental composition is the very strong background for volatilizing and condensating alkaline elements.Preliminary conclusions about a possible correlation between the composition and the sizes of dust particles in the Solar system and in the Universe (at similar initial composition) have been drawn.     

Atmospheric impacts, fragmentation, and small craters on Venus

We use high-resolution three-dimensional numerical models of aerodynamically disrupted asteroids to predict the characteristic properties of small impact craters on Venus. We map the mass and kinetic energy of the impactor passing though a plane near the surface for each simulation, and find that the typical result is that mass and energy sort themselves into one to several strongly peaked regions, which we interpret as more-or-less discrete fragments. The fragments are sufficiently well separated as to imply the formation of irregular or multiple craters that are quite similar to those found on Venus. We estimate the diameters of the resulting craters using a scaling law derived from the experiments of Schultz and Gault (1985, J. Geophys. Res. 90 (B5), 3701-3732) of dispersed impactors into targets. We compare the spacings and sizes of our estimated craters with measured diameters tabulated in a Venus crater database (Herrick and Phillips, 1994a, Icarus 111, 387-416; Herrick et al., 1997, in: Venus II, Univ. of Arizona Press, Tucson, AZ, pp. 1015-1046; Herrick, 2003, http://www.lpi.usra.edu/research/vc/vchome.html) and find quite satisfactory agreement, despite the uncertainty in our crater diameter estimates. The comparison of the observed crater characteristics with the numerical results is an after-the-fact test of our model, namely the fluid-dynamical treatment of large impacts, which the model appears to pass successfully.     

Investigating pressure magnitudes at depth for oblique impacts into layered targets: Applications to terrestrial impacts in sedimentary targets

The positive identification of the Rock Elm impact structure (Wisconsin, USA) and the Upheaval Dome (Utah, USA) as impact craters was complicated by a lack of distinctive shock features in the record. Low‐impedance surface layers over high‐impedance bedrock affect energy coupling and shock effects in the substrate; in both cases, removal of surface sediments erased most of the original impact structures, thereby making identification of the impact origin difficult. In this study, a combination of laboratory and 3‐D numerical experiments reveals the underlying processes controlling subsurface deformation and demonstrates that a low‐impedance layer can reduce expression of peak shock pressures left in the rock record, as at the Rock Elm and the Upheaval Dome impact sites. 3‐D CTH models of the Rock Elm impact structure predict that peak shock pressures should fall below the hugoniot elastic limit of quartz in the basement rocks, yet still induce permanent deformation. The model predicts peak pressures around 5–10 GPa, levels consistent with field observations of shocked quartz from both Rock Elm and the Upheaval Dome. Consequently, other impact sites exhibiting minimal shock features might be explained.     

Experimental ejection angles for oblique impacts: Implications for the subsurface flow‐field

Abstract— A simple analytical solution for subsurface particle motions during impact cratering is useful for tracking the evolution of the transient crater shape at late times. A specific example of such an analytical solution is Maxwell's Z‐Model, which is based on a point‐source assumption. Here, the parameters for this model are constrained using measured ejection angles from both vertical and oblique experimental impacts at the NASA Ames Vertical Gun Range. Data from experiments reveal that impacts at angles as high as 45° to the target's surface generate subsurface flow‐fields that are significantly different from those created by vertical impacts. The initial momentum of the projectile induces a subsurface momentum‐driven flow‐field that evolves in three dimensions of space and in time to an excavation flow‐field during both vertical and oblique impacts. A single, stationary point‐source model (specifically Maxwell's Z‐Model), however, is found inadequate to explain this detailed evolution of the subsurface flow‐field during oblique impacts. Because 45° is the most likely impact angle on planetary surfaces, a new analytical model based on a migrating point‐source could prove quite useful. Such a model must address the effects of the subsurface flow‐field evolution on crater excavation, ejecta deposition, and transient crater morphometry.     

Scaling of oblique impacts in frictional targets: Implications for crater size and formation mechanisms

Almost every meteorite impact occurs at an oblique angle of incidence, yet the effect of impact angle on crater size or formation mechanism is only poorly understood. This is, in large part, due to the difficulty of inferring impactor properties, such as size, velocity and trajectory, from observations of natural craters, and the expense and complexity of simulating oblique impacts using numerical models. Laboratory oblique impact experiments and previous numerical models have shown that the portion of the projectile’s kinetic energy that is involved in crater excavation decreases significantly with impact angle. However, a thorough quantification of planetary-scale oblique impact cratering does not exist and the effect of impact angle on crater size is not considered by current scaling laws. To address this gap in understanding, we developed iSALE-3D, a three-dimensional multi-rheology hydrocode, which is efficient enough to perform a large number of well-resolved oblique impact simulations within a reasonable time. Here we present the results of a comprehensive numerical study containing more than 200 three-dimensional hydrocode-simulations covering a broad range of projectile sizes, impact angles and friction coefficients. We show that existing scaling laws in principle describe oblique planetary-scale impact events at angles greater than 30° measured from horizontal. The displaced mass of a crater decreases with impact angle in a sinusoidal manner. However, our results indicate that the assumption that crater size scales with the vertical component of the impact velocity does not hold for materials with a friction coefficient significantly lower than 0.7 (sand). We found that increasing coefficients of friction result in smaller craters and a formation process more controlled by impactor momentum than by energy.     

Release of neutral sodium atoms from the surface of Mercury induced by meteoroid impacts

Meteoroid impact has been shown to be a source of sodium, and most likely of other elements, on the Moon. The same process could be also relevant for Mercury. In this work we calculate the vapor and neutral Na production rates on Mercury due to the impacts of meteoroids in the radius range of 10⁻⁸-10⁻¹ m. We limit our calculations to this size range, because meteoroids with radius larger than 10⁻¹ m have not to be found important for the daily production of the exosphere. This work is based on a new dynamical model of the meteoroid flux at the heliocentric distance of Mercury, regarding objects in the size range 10⁻²-10⁻¹ m. This size range, never investigated before, is not affected by nongravitational forces, such as the Poynting-Robertson effect, which is dominant for particles smaller than 10⁻² m. In order to evaluate the release of neutral sodium atoms also for smaller meteoroids we have used the distribution reported by M.J. Cintala [1992. Impact-induced thermal effects in the lunar and mercurian regoliths. J. Geophys. Res. 97, 947-973] calculated for particle size range 10⁻⁸-10⁻³ m. We have extrapolated this distribution up to 10⁻² m and we have based the impact calculations on a new surface composition assuming 90% plagioclase and 10% pyroxene. The results of our model are that (i) the total mass of vapor produced by the impact of meteoroids in the size range 10⁻⁸-10⁻¹ m is 4.752×¹⁰⁸ g per year, and (ii) the production rate of neutral sodium atoms is 1.5×¹⁰²² s⁻¹.     

Neutral sodium atoms release from the surface of the Moon induced by meteoroid impacts

In this work, we calculate the vapour and neutral Na production rates on the Moon, as due to the impacts of meteoroids in the radius range of 10⁻⁸–0.15 m. We limit our calculations to this size range, since meteoroids with radius larger than 0.15 m have not been found to be important for the production of the exosphere in a time interval comparable with that of the observations.
We have considered a new dynamical model of the flux of meteoroids at the heliocentric distance of the Moon, regarding objects in the radius range of 10⁻²–0.15 m. Instead, the flux of smaller meteoroids (radius range 10⁻⁸–10⁻² m) has been calculated using the two distributions adopted by Cintala and Love & Brownlee.
The results of our model are that (i) the neutral Na production rate is  ∼3–4.9 × 10⁴ atoms cm⁻² s⁻¹  , slightly larger than the previous estimates  (∼2–3 × 10⁴ atoms cm⁻² s⁻¹)  , and (ii) only about 6 per cent of neutral Na is produced by the impacts of meteoroids in the size range 10⁻³–0.15 m, whereas about 94 per cent of the Na comes from the  10⁻⁵–10⁻³ m  size range.     

Neutral sodium atoms release from the surfaces of the Moon and Mercury induced by meteoroid impacts

In this work, we calculate the neutral Na production rates on the Moon and Mercury, as due to the impacts of meteoroids having an impact probability on the surface that can influence the daily observations of the exosphere: the meteoroids radius range considered for the Moon and Mercury are 10⁻⁸–0.15 and 10⁻⁸–0.10 m, respectively. We also estimate the mass of meteoroids that has impacted the surfaces of the Moon and Mercury in the last 3.8 Gy (after the end of the Late Heavy Bombardment).The results of our model are that (i) the Na production rates are ∼(3–4.9)×10⁴ and ∼(1.8–2.3)×10⁶ atoms cm⁻² s⁻¹, for Moon and Mercury, respectively, and (ii) in the last 3.8 Gy, the mass of meteoroids that has impacted the whole surface of the Moon and Mercury has been 8.86×10¹⁸ and 2.66×10¹⁹ g, respectively.     

Dynamical erosion of the asteroid belt and implications for large impacts in the inner Solar System

The cumulative effects of weak resonant and secular perturbations by the major planets produce chaotic behavior of asteroids on long timescales. Dynamical chaos is the dominant loss mechanism for asteroids with diameters in the current asteroid belt. In a numerical analysis of the long-term evolution of test particles in the main asteroid belt region, we find that the dynamical loss history of test particles from this region is well described with a logarithmic decay law. In our simulations the loss rate function that is established at persists with little deviation to at least . Our study indicates that the asteroid belt region has experienced a significant amount of depletion due to this dynamical erosion—having lost as much as ∼50% of the large asteroids—since 1 Myr after the establishment of the current dynamical structure of the asteroid belt. Because the dynamical depletion of asteroids from the main belt is approximately logarithmic, an equal amount of depletion occurred in the time interval 10-200 Myr as in 0.2-4 Gyr, roughly ∼30% of the current number of large asteroids in the main belt over each interval. We find that asteroids escaping from the main belt due to dynamical chaos have an Earth-impact probability of ∼0.3%. Our model suggests that the rate of impacts from large asteroids has declined by a factor of 3 over the last 3 Gyr, and that the present-day impact flux of objects on the terrestrial planets is roughly an order of magnitude less than estimates currently in use in crater chronologies and impact hazard risk assessments.     

Atmospheric shock waves after impacts of cosmic bodies up to 1000 m in diameter

The results of numerical modeling of meteoroids' interaction with Earth's atmosphere are presented. We model the entry in two dimensions and then interpolate the results into a 3‐D model to calculate interaction of shock waves with the surface. Maximum shock pressures, wind speeds, and areas subjected to substantial overpressure are calculated for oblique impacts of asteroids and comets. We show that vertical impacts produce a smaller damage zone on the surface than oblique ones. Damage caused by shock waves covers an order of magnitude larger area than any other hazardous effects. The function of energy release in the atmosphere, which is traditionally used in meteoritics, has a limited application if cosmic bodies are larger than tens of meters in diameter: at each time moment energy is smoothed along a substantial length of the trajectory; both emitted radiation (routinely used for calibration of semi‐analytical models) and shock wave amplitude are complex functions of temperature–density distributions in atmosphere.     

Atmospheric mass deposition by captured planetesimals

We examine the deposition of heavy elements in the envelope of a protoplanet growing according to the core accretion scenario of Pollack et al. [Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y., 1996. Icarus 124, 62-85]. We use their atmospheric models and the deposition rates that they computed, and we calculate the amount of heavy elements that can be dissolved in the envelope. For planetesimals composed of a mixture of water, CHON, and rock, we find that almost all of the water is dissolved in the atmosphere. A substantial amount of CHON is also dissolved but it remains sequestered at high temperatures near the core. Some fraction of the rock is also dissolved in the very high temperature region near the core envelope boundary. If this dissolved material can be mixed upward later in the planet's evolution, the resulting structure would be much closer to that determined by matching the moments of Jupiter's gravitational field.     

Laboratory simulation of cometary erosion by space plasma

A laboratory experiment has been made where a plasma stream collides with targets made of different materials of cosmic interest. The experiment can be viewed as a process simulation of the solar wind particle interaction with solid surfaces in space — e.g., cometary dust. Special interest is given to sputtering of OH and Na.It is also shown that the erosion of solid particles in interplanetary space at large heliocentric distances is most likely dominated by sputtering and by sublimation near the Sun. The heliocentric distance of the limit between the two regions is determined mainly by the material properties of the eroded surface, e.g., heat of sublimation and sputtering yield, a typical distance being 0.5 AU.It is concluded that the observations of Na in comets at large solar distances, in some cases also near the Sun, is most likely to be explained by solar wind sputtering. OH emission in space could be of importance also from dry, water-free, matter by means of molecule sputtering. The observed OH production rates in comets are however too large to be explained in this way and are certainly the results of sublimation and dissociation of H₂O from an icy nucleus.