Quantification of Impact-Induced Melt Production in Numerical Modeling Revisited

Melting and vaporization of rocks in impact cratering is mostly attributed to be a consequence of shock compression. However, other mechanism such as plastic work and decompression by structural uplift also contribute to melt production. In this study we expand the commonly used method to determine shock-induced melting in numerical models from the peak shock pressure by a new approach to account for additional heating due plastic work and internal friction. We compare our new approach with the straight-forward method to simply quantify melting from the temperature relative to the solidus temperature at any arbitrary point in time in the course of crater formation. This much simpler method does account for plastic work but suffers from reduced accuracy due to numerical diffusion inherent to ongoing advection in impact crater formation models. We demonstrate that our new approach is more accurate than previous methods in particular for quantitative determination of impact melt distribution in final crater structures. In addition, we assess the contribution of plastic work to the overall melt volume and find, that melting is dominated by plastic work for impacts at velocities smaller than 7.5–12.5 km/s in rocks, depending on the material strength. At higher impact velocities shock compression is the dominating mechanism for melting. Here, the conventional peak shock pressure method provides similar results compared with our new model. Our method serves as a powerful tool to accurately determine impact-induced heating in particular at relatively low-velocity impacts

3D-simulation of lunar megaregolith evolution: Quantitative constraints on spatial variation and size of fragment

The early impact bombardment extensively fractured the lunar crust resulting in the formation of the so-called megaregolith. Previous estimates of megaregolith distribution vary significantly with respect to the vertical extent and the size-frequency distribution of fragments was rarely studied. We built a spatially resolved numerical model to simulate the process of cumulative impact fragmentation, aiming to backtrack the megaregolith evolution history and to constrain its fragment distribution. The results highlight the pivotal role of basin-forming events on the megaregolith formation. Especially the South-Pole Aitken (SPA) impact established the initial megaregolith structure which remained distinct after 0.5 Ga subsequent fragmentation. At 3.8 Ga, the megaregolith displays substantial lateral variation and layering: the highly fractured upper layer of ∼2.5 km is dominated by meter-scale fragments; the disturbed lower layer deeper than tens of kilometers is mainly consisting of kilometer-scale fragments; the transition zone >5 km contains fragments of various size scales

The Timeline of Early Lunar Bombardment Constrained by the Evolving Distributions of Differently Aged Melt

The timeline of the early lunar bombardment remains unclear. The bombardment rate as a function of time is commonly modeled by three types of shapes: tail-end, sawtooth, and terminal cataclysm. Differently aged melt records the occurrence time of impact events and thus is crucial for constraining the timeline of the early lunar bombardment. Based on a spatially resolved numerical model, we simulate the evolving distribution of differently aged melt with a long-term impact mixing, where different shapes of impact rate function are considered. We compare the outcome of melt age distribution from different scenarios with the actual data from the lunar meteorites and the returned samples. The results suggest that, if the present data are representative of the melt age distribution on the Moon, the shape of the impact rate function is more likely comparable to the tail-end over the sawtooth and the terminal cataclysm, with the terminal cataclysm being least likely. In addition, using state-of-the-art U–Pb dating techniques, more abundant ancient basin melt is likely to be found in returned samples

Melting efficiency of troilite-iron assemblages in shock-darkening : insight from numerical modeling

We studied shock-darkening in ordinary chondrites by observing the propagation of shock waves and melting through mixtures of metals and iron sulfides. We used the shock physics code iSALE at the mesoscale to simulate shock compression of modeled ordinary chondrites (using olivine, iron and troilite). We introduced FeS-FeNi eutectic properties and partial melting in a series of chosen configurations of iron and troilite grains mixtures in a sample plate. We observed, at a nominal pressure of 45 GPa, partial melting of troilite in all models. Only few of the models showed partial melting of iron (a phase difficult to melt in shock heating) due to the eutectic properties of the mixtures. Iron melting only occurred in models presenting either strong shock wave concentration effects or effects of heating by pore crushing, for which we provided more details. Further effects are discussed such as the frictional heating between iron and troilite and the heat diffusion in scenarios with strongly heated troilite. We also characterized troilite melting in the 32-60 GPa nominal pressure range. We concluded that specific dispositions of the iron and troilite grains in mixtures exist that lead to melting of iron and explain why it is possible to find a mix of metals and iron sulfides in shock-darkened ordinary chondrites.Peer reviewe

Melt Production and Ejection From Lunar Intermediate-Sized Impact Craters: Where Is the Molten Material Deposited?

Differently aged impact melt in lunar samples is key to unveiling the early bombardment history of the Moon. Due to the mixing of melt products ejected from distant craters, the interpretations of the origin of lunar samples are difficult. We use numerical modeling for a better quantitative understanding of the production of impact-induced melt and in particular its distribution in ejecta blankets for lunar craters with sizes ranging from 1.5 to 50 km. We approximate the lunar stratigraphy with a porosity gradient, which represents the gradual transition from upper regolith via megaregolith to the solid crustal material. For this lunar setting, we quantify the melt production relative to crater volume and derive parameters describing its increasing trend with increasing transient crater size. We found that about 30%–40% of the produced melt is ejected from the crater. The melt concentration in the ejecta blanket increases almost linearly with distance from the crater center, while the thickness of the ejecta blanket decreases following a power law. Our study demonstrates that if in lunar samples the concentration of a melt with a certain age is interpreted to be of a nonlocal origin, these melts could be the impact products of a large crater (>10 km) located hundreds of kilometers away

Mineralogy and Provenance of Quaternary Sediments in the Gaxun Nur Basin, Northwestern Inner Mongolia, China

Abstract HKT-ISTP 2013 B

Shock physics mesoscale modeling of shock stage 5 and 6 in ordinary and enstatite chondrites

Shock-darkening, the melting of metals and iron sulfides into a network of veins within silicate grains, altering reflectance spectra of meteorites, was previously studied using shock physics mesoscale modeling. Melting of iron sulfides embedded in olivine was observed at pressures of 40-50 GPa. This pressure range is at the transition between shock stage 5 (C-S5) and 6 (C-S6) of the shock metamorphism classification in ordinary and enstatite chondrites. To better characterize C-S5 and C-S6 with a mesoscale modeling approach and assess post-shock heating and melting, we used multi-phase (i.e. olivine/enstatite, troilite, iron, pores, and plagioclase) meshes with realistic configurations of grains. We carried out a systematic study of shock compression in ordinary and enstatite chondrites at pressures between 30 and 70 GPa. To setup mesoscale sample meshes with realistic silicate, metal, iron sulfide, and open pore shapes, we converted backscattered electron microscope images of three chondrites. The resolved macroporosity in meshes was 3-6%. Transition from shock C-S5 to C-S6 was observed through (1) the melting of troilite above 40 GPa with melt fractions of similar to 0.7-0.9 at 70 GPa, (2) the melting of olivine and iron above 50 GPa with melt fraction of similar to 0.001 and 0.012, respectively, at 70 GPa, and (3) the melting of plagioclase above 30 GPa (melt fraction of 1, at 55 GPa). Post-shock temperatures varied from similar to 540 K at 30 GPa to similar to 1300 K at 70 GPa. We also constructed models with increased porosity up to 15% porosity, producing higher post-shock temperatures (similar to 800 K increase) and melt fractions (similar to 0.12 increase) in olivine. Additionally we constructed a pre-heated model to observe post-shock heating and melting during thermal metamorphism. This model presented similar results (melting) at pressures 10-15 GPa lower compared to the room temperature models. Finally, we demonstrated dependence of post-shock heating and melting on the orientation of open cracks relative to the shock wave front. In conclusion, the modeled melting and post-shock heating of individual phases were mostly consistent with the current shock classification scheme (Stoffler et al., 1991, 2018).Peer reviewe

Insight into the Distribution of High-pressure Shock Metamorphism in Rubble-pile Asteroids

Funding Information: This work was supported by the Academy of Finland, project Nos. 293975 and 335595, the European Regional Development Fund, the Mobilitas Pluss programme (grant No. MOBJD639), and the NASA Solar System Exploration Research Virtual Institute Center for Lunar and Asteroid Surface Science, and it was conducted within institutional support RVO 67985831 of the Institute of Geology of the Czech Academy of Sciences. R.L. appreciates funding from the European Union’s Horizon 2020 research and innovation program, NEO-MAPP, grant agreement No. 870377. Publisher Copyright: © 2022. The Author(s). Published by the American Astronomical Society.Shock metamorphism in ordinary chondrites allows for reconstructing impact events between asteroids in the main asteroid belt. Shock-darkening of ordinary chondrites occurs at the onset of complete shock melting of the rock (>70 GPa) or injection of sulfide and metal melt into the cracks within solid silicates (∼50 GPa). Darkening of ordinary chondrites masks diagnostic silicate features observed in the reflectance spectrum of S-complex asteroids so they appear similar to C/X-complex asteroids. In this work, we investigate the shock pressure and associated metamorphism pattern in rubble-pile asteroids at impact velocities of 4–10 km s−1. We use the iSALE shock physics code and implement two-dimensional models with simplified properties in order to quantify the influence of the following parameters on shock-darkening efficiency: impact velocity, porosity within the asteroid, impactor size, and ejection efficiency. We observe that, in rubble-pile asteroids, the velocity and size of the impactor are the constraining parameters in recording high-grade shock metamorphism. Yet, the recorded fraction of higher shock stages remains low (<0.2). Varying the porosity of the boulders from 10% to 30% does not significantly affect the distribution of pressure and fraction of shock-darkened material. The pressure distribution in rubble-pile asteroids is very similar to that of monolithic asteroids with the same porosity. Thus, producing significant volumes of high-degree shocked ordinary chondrites requires strong collision events (impact velocities above 8 km s−1 and/or large sizes of impactors). A large amount of asteroid material escapes during an impact event (up to 90%); however, only a small portion of the escaping material is shock-darkened (6%).Peer reviewe

Earth-like Habitats in Planetary Systems

Understanding the concept of habitability is related to an evolutionary knowledge of the particular planet-in-question. Additional indications so-called "systemic aspects" of the planetary system as a whole governs a particular planet's claim on habitability. Here we focus on such systemic aspects and discuss their relevance to the formation of an 'Earth-like' habitable planet. We summarize our results obtained by lunar sample work and numerical models within the framework of the Research Alliance "Planetary Evolution and Life". We consider various scenarios which simulate the dynamical evolution of the Solar System and discuss the likelihood of forming an Earth-like world orbiting another star. Our model approach is constrained by observations of the modern Solar System and the knowledge of its history. Results suggest that the long-term presence of terrestrial planets is jeopardized due to gravitational interactions if giant planets are present. But habitability of inner rocky planets may be supported in those planetary systems hosting giant planets. Gravitational interactions within a complex multiple-body structure including giant planets may supply terrestrial planets with materials which formed in the colder region of the proto-planetary disk. During these processes, water, the prime requisite for habitability, is delivered to the inner system. This may occur either during the main accretion phase of terrestrial planets or via impacts during a post-accretion bombardment. Results for both processes are summarized and discussed with reference to the lunar crater record. Starting from a scenario involving migration of the giant planets this contribution discusses the delivery of water to Earth, the modification of atmospheres by impacts in a planetary system context and the likelihood of the existence of extrasolar Earth-like habitable worlds.Comment: 36 Pages, 6 figures, 2014, Special Issue in Planetary and Space Science on the Helmholtz Research Alliance on Planetary Evolution and Lif

Constraining the characteristics of tsunami waves from deformable submarine slides

As a marine hazard, submarine slope failures have the potential to directly destroy offshore infrastructure, and, if a tsunami is generated, it also endangers the life of those who live and work at the coastline. The hazard and risk from tsunamis generated by submarine mass failure is difficult to quantify and evaluate due to the problems to constrain the characteristics of the triggered submarine landslide, which introduces unquantifiable uncertainty to hazard assessments based on numerical modelling. To lower the uncertainty, we present a method that determines material parameters for the slide body to constrain the generated tsunami waves. Our method employs the distribution of landslide run-out masses and their comparison with simulations. It assumes that the slide material can be approximated by bulk values during the slide motion. To demonstrate our method, we make use of Valdes slide run-out masses off the Chilean coas