•Impact-induced seismicity plays an important role on planetary surfaces.•Seismic signals are recorded using numerical modeling.•Seismic parameters are significantly affected by different target ...properties.•Seismic magnitudes resulting from impacts can be predicted using numerical modeling.
We quantify the seismicity of impact events using a combined numerical and experimental approach. The objectives of this work are (1) the calibration of the numerical model by utilizing real-time measurements of the elastic wave velocity and pressure amplitudes in laboratory impact experiments; (2) the determination of seismic parameters, such as quality factor Q and seismic efficiency k, for materials of different porosity and water saturation by a systematic parameter study employing the calibrated numerical model. By means of “numerical experiments” we found that the seismic efficiency k decreases slightly with porosity from k=3.4×10−3 for nonporous quartzite to k=2.6×10−3 for 25% porous sandstone. If pores are completely or partly filled with water, we determined a seismic efficiency of k=8.2×10−5, which is approximately two orders of magnitude lower than in the nonporous case. By measuring the attenuation of the seismic wave with distance in our numerical experiments we determined the seismic quality factor Q to range between ∼35 for the solid quartzite and 80 for the porous dry targets. For water saturated target materials, Q is much lower, <10. The obtained values are in the range of literature values. Translating the seismic efficiency into seismic magnitudes we show that the seismic magnitude of an impact event is about one order of magnitude smaller considering a water saturated target in comparison to a solid or porous target. Obtained seismic magnitudes decrease linearly with distance to the point of impact and are consistent with empirical data for distances closer to the point of impact. The seismic magnitude decreases more rapidly with distance for a water saturated material compared to a dry material.
–
Generation and propagation of shock waves by meteorite impact is significantly affected by material properties such as porosity, water content, and strength. The objective of this work was to ...quantify processes related to the shock‐induced compaction of pore space by numerical modeling, and compare the results with data obtained in the framework of the Multidisciplinary Experimental and Modeling Impact Research Network (MEMIN) impact experiments. We use mesoscale models resolving the collapse of individual pores to validate macroscopic (homogenized) approaches describing the bulk behavior of porous and water‐saturated materials in large‐scale models of crater formation, and to quantify localized shock amplification as a result of pore space crushing. We carried out a suite of numerical models of planar shock wave propagation through a well‐defined area (the “sample”) of porous and/or water‐saturated material. The porous sample is either represented by a homogeneous unit where porosity is treated as a state variable (macroscale model) and water content by an equation of state for mixed material (ANEOS) or by a defined number of individually resolved pores (mesoscale model). We varied porosity and water content and measured thermodynamic parameters such as shock wave velocity and particle velocity on meso‐ and macroscales in separate simulations. The mesoscale models provide additional data on the heterogeneous distribution of peak shock pressures as a consequence of the complex superposition of reflecting rarefaction waves and shock waves originating from the crushing of pores. We quantify the bulk effect of porosity, the reduction in shock pressure, in terms of Hugoniot data as a function of porosity, water content, and strength of a quartzite matrix. We find a good agreement between meso‐, macroscale models and Hugoniot data from shock experiments. We also propose a combination of a porosity compaction model (ε–α model) that was previously only used for porous materials and the ANEOS for water‐saturated quartzite (all pore space is filled with water) to describe the behavior of partially water‐saturated material during shock compression. Localized amplification of shock pressures results from pore collapse and can reach as much as four times the average shock pressure in the porous sample. This may explain the often observed localized high shock pressure phases next to more or less unshocked grains in impactites and meteorites.
A combination of shock recovery experiments and numerical modeling of shock deformation in the low‐shock pressure range from 2.5 to 20 GPa for two dry sandstone types of different porosity, a ...completely water‐saturated sandstone, and a well‐indurated quartzite provides new insights into strongly heterogeneous distribution of different shock features. (1) For nonporous quartzo‐feldspathic rocks, the traditional classification scheme (Stöffler ) is suitable with slight changes in pressure calibration. (2) For water‐saturated quartzose rocks, a cataclastic texture (microbreccia) seems to be typical for the shock pressure range up to 20 GPa. This microbreccia does not show formation of PDFs but diaplectic quartz glass/SiO2 melt is formed at 20 GPa (~1 vol%). (3) For porous quartzose rocks, the following sequence of shock features is observed with progressive increase in shock pressure (1) crushing of pores, (2) intense fracturing of quartz grains, and (3) increasing formation of diaplectic quartz glass/SiO2 melt replacing fracturing. The formation of diaplectic quartz glass/SiO2 melt, together with SiO2 high‐pressure phases, is a continuous process that strongly depends on porosity. This experimental observation is confirmed by our concomitant numerical modeling. Recalibration of the shock classification scheme results in a porosity versus shock pressure diagram illustrating distinct boundaries for the different shock stages.
The Double Asteroid Redirection Test (DART) is a NASA-sponsored mission that will be the first direct test of the kinetic impactor technique for planetary defense. The DART spacecraft will impact ...into Didymos-B, the moon of the binary system 65803 Didymos, and the resulting period change will be measured from Earth. Impact simulations will be used to predict the crater size and momentum enhancement expected from the DART impact. Because the specific material properties (strength, porosity, internal structure) of the Didymos-B target are unknown, a wide variety of numerical simulations must be performed to better understand possible impact outcomes. This simulation campaign will involve a large parameter space being simulated using multiple different shock physics hydrocodes. In order to understand better the behaviors and properties of numerical simulation codes applicable to the DART impact, a benchmarking and validation program using different numerical codes to solve a set of standard problems was designed and implemented. The problems were designed to test the effects of material strength, porosity, damage models, and target geometry on the ejecta following an impact and thus the momentum transfer efficiency. Several important results were identified from comparing simulations across codes, including the effects of model resolution and porosity and strength model choice: 1) momentum transfer predictions almost uniformly exhibit a larger variation than predictions of crater size; 2) the choice of strength model, and the values used for material strength, are significantly more important in the prediction of crater size and momentum enhancement than variation between codes; 3) predictions for crater size and momentum enhancement tend to be similar (within 15‐20%) when similar strength models are used in different codes. These results will be used to better design a modeling plan for the DART mission as well as to better understand the potential results that may be expected due to unknown target properties. The DART impact simulation team will determine a specific desired material parameter set appropriate for the Didymos system that will be standardized (to the extent possible) across the different codes when making predictions for the DART mission. Some variation in predictions will still be expected, but that variation can be bracketed by the results shown in this study.
•DART is the first space-based test of a kinetic impactor for momentum enhancement.•Porosity has large effects on crater size and significantly affects β predictions.•Predictions of momentum enhancement strongly depend on resolution.•Momentum enhancement predictions show larger variation than crater size predictions.•Strength model choice is more important than inherent inter-code variation.
This paper reviews major findings of the Multidisciplinary Experimental and Modeling Impact Crater Research Network (MEMIN). MEMIN is a consortium, funded from 2009 till 2017 by the German Research ...Foundation, and is aimed at investigating impact cratering processes by experimental and modeling approaches. The vision of this network has been to comprehensively quantify impact processes by conducting a strictly controlled experimental campaign at the laboratory scale, together with a multidisciplinary analytical approach. Central to MEMIN has been the use of powerful two‐stage light‐gas accelerators capable of producing impact craters in the decimeter size range in solid rocks that allowed detailed spatial analyses of petrophysical, structural, and geochemical changes in target rocks and ejecta. In addition, explosive setups, membrane‐driven diamond anvil cells, as well as laser irradiation and split Hopkinson pressure bar technologies have been used to study the response of minerals and rocks to shock and dynamic loading as well as high‐temperature conditions. We used Seeberger sandstone, Taunus quartzite, Carrara marble, and Weibern tuff as major target rock types. In concert with the experiments we conducted mesoscale numerical simulations of shock wave propagation in heterogeneous rocks resolving the complex response of grains and pores to compressive, shear, and tensile loading and macroscale modeling of crater formation and fracturing. Major results comprise (1) projectile–target interaction, (2) various aspects of shock metamorphism with special focus on low shock pressures and effects of target porosity and water saturation, (3) crater morphologies and cratering efficiencies in various nonporous and porous lithologies, (4) in situ target damage, (5) ejecta dynamics, and (6) geophysical survey of experimental craters.
Abstract
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%).
We present the outcomes of simulations of the formation of the Vista Alegre impact structure, Paraná Basin, Brazil. The target comprised a thick sequence of volcanic rocks of predominantly basaltic ...composition of the Serra Geral Formation that had been deposited on top of sedimentary rocks (sandstones) of the Pirambóia/Botucatu formations. The cratering process was modeled using the iSALE shock physics code. Our best‐fit model suggests that (1) the crater was originally ~10 km in size; (2) it was formed in ~115 s by a stony projectile of 1000 m in diameter, for an assumed impact velocity of 12 km s−1; (3) target rocks underwent a peak pressure of ~20 GPa, in agreement with previous petrographic investigations of shock deformation. Furthermore, the model points out that the sedimentary strata below the layer of volcanic rocks were raised by ~650 meters at the central part of the crater, which resulted in the current partial exposure of the sandstones at the surface. The outcomes of our modeling suggest that parameters like cohesion and strength of the target rocks, after shock compression, determined the final morphology of the crater, especially the absence of a topographically prominent central peak. Finally, the results of the numerical modeling are roughly in agreement with gravity data over the structure, in particular with respect to the presence of the uplifted sedimentary strata, which are responsible for a low gravity signature at the center of the structure.
The Double Asteroid Redirection Test (DART) is a NASA-sponsored mission that will be the first direct test of the kinetic impactor technique for planetary defense. Here, the DART spacecraft will ...impact into Didymos-B, the moon of the binary system 65803 Didymos and the resulting period change will be measured from Earth. Impact simulations will be used to predict the crater size and momentum enhancement expected from the DART impact. Because the specific material properties (strength, porosity, internal structure) of the Didymos-B target are unknown, a wide variety of numerical simulations must be performed to better understand possible impact outcomes. This simulation campaign will involve a large parameter space being simulated using multiple different shock physics hydrocodes. In order to understand better the behaviors and properties of numerical simulation codes applicable to the DART impact, a benchmarking and validation program using different numerical codes to solve a set of standard problems was designed and implemented.
In this paper, the behavior of natural sandstone under hypervelocity impact is investigated numerically. Sandstone essentially consists of cemented quartz grains, it can be generally considered as a ...polygranular structure with air or water-filled pores. In order to study the influence of pores and water filling on the behavior of sandstone under shock loading, a numerical investigation at the mesoscale is performed. We explicitly resolve pores and the closure of pores in a homogeneous and isotropic quartz matrix. The quartz matrix is described with an equation of state and a strength model; the pore geometry is simplified. Two-dimensional planar shock simulations are conducted in order to analyze effects of quartz strength and water-saturation on pore crushing. Averaging is used to perform a macroscale analysis of the dependencies of pressure and shock velocity on particle velocity. The simulation results shall enhance the basis for the development of accurate macroscale sandstone models that take the water content into account.
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 (\sim 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. 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 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%).