Search

The environmental effects of impact events differ with respect to time (seconds to decades) and spatial (local to global) scales. Short-term localized damage is produced by thermal radiation, blast-wave propagation in the atmosphere, crater excavation, earthquakes, and tsunami. Global and long-term effects are related to the ejection of dust and climate-active gases (carbon dioxide, sulfur oxides, water vapor, methane) into the atmosphere. At the end of the Cretaceous, the impact of a >10 km diameter asteroid led to a major mass extinction. Modern civilization is vulnerable to even relatively small impacts, which may occur in the near future, that is, tens to hundreds of years.

Modeling results of the water vapor plume produced by a comet impact on the Moon and of the resulting water ice deposits in the lunar cold traps are presented. The water vapor plume is simulated near the point of impact by the SOVA hydrocode and in the far field by the Direct Simulation Monte Carlo (DSMC) method using as input the SOVA hydrocode solution at a fixed hemispherical interface. The SOVA hydrocode models the physics of the impact event such as the surface deformation and material phase changes during the impact. The further transport and retention processes, including gravity, photodestruction processes, and variable surface temperature with local polar cold traps, are modeled by the DSMC method for months after impact. In order to follow the water from the near field of the impact to the full planetary induced atmosphere, the 3D parallel DSMC code used a collision limiting scheme and an unsteady multi-domain approach. 3D results for the 45° oblique impact of a 2 km in diameter comet on the surface of the Moon at 30 km/s are presented. Most of the cometary water is lost due to escape just after impact and only ∼3% of the cometary water is initially retained on the Moon. Early downrange focusing of the water vapor plume is observed but the later material that is moving more slowly takes on a more symmetric shape with time. Several locations for the point of impact were investigated and final retention rates of ∼0.1% of the comet mass were observed. Based on the surface area of the cold traps used in the present simulations, ∼1 mm of ice would have accumulated in the cold traps after such an impact. Estimates for the total mass of water accumulated in the polar cold traps over 1 byr are consistent with recent observations.

– Despite its centennial exploration history, there are still unresolved questions about Meteor Crater, the first recognized impact crater on Earth. This theoretical study addresses some of these questions by comparing model results with field and laboratory studies of Meteor Crater. Our results indicate that Meteor Crater was formed by a high-velocity impact of a fragmented projectile, ruling out a highly dispersed swarm as well as a very low impact velocity. Projectile fragmentation caused many fragments to fall separately from the main body of the impactor, making up the bulk of the Canyon Diablo meteorites; most of these fragments were engulfed in the expansion plume as they approached the surface without suffering high shock compression, and were redistributed randomly around the crater. Thus, the distribution of Canyon Diablo meteorites is not representative of projectile trajectory, as is usual for impactor fragments in smaller strewn fields. At least 50% of the main impactor was ejected from the crater during crater excavation and was dispersed mostly downrange of the crater as molten particles (spheroids) and highly shocked solid fragments (shrapnel). When compared with the known distribution, model results suggest an impactor from the SW. Overall, every model case produced much higher amounts of pure projectile material than observed. The projectile-target mixing was not considered in the models; however, this process could be the main sink of projectile melt, as all analyzed melt particles have high concentrations of projectile material. The fate of the solid projectile fragments is still not completely resolved. Model results suggest that the depth of melting in the target can reach the Coconino sandstone formation. However, most of the ejected melt originates from 30–40 m depth and, thus, is limited to Moenkopi and upper Kaibab material. Some melt remains in the target; based on the estimated volume of the breccia lens at Meteor Crater, our models suggest at most a 2% content of melt in the breccia. Finally, a high water table at the time of impact could have aided strong dispersion of target and projectile melt.

We present results of an investigation aimed at characterizing the effects of oceanic impacts of 500 m and 1 km diameter asteroids on the lower and middle atmosphere, estimating ozone loss and potential danger from UV radiation at the Earth's surface. Little work has been done so far to assess the atmospheric perturbation from the impact of objects in this size range, even though their sizes are close to the threshold for causing global environmental effects. In particular, at the Earth's surface oceanic impacts are twice more likely to occur than land impacts. This work represents the first attempt at combining impact simulations with a three-dimensional shock physics code (SOVA), and atmospheric simulations using the general circulation model with interactive chemistry WACCM. SOVA simulations provided an estimate of the amount and state of material ejected into the atmosphere by the impacts. Estimated water vapor in the upper atmosphere was then introduced in the initial conditions for the WACCM simulations that modeled the subsequent perturbation of atmospheric chemistry Final estimates of the change over time in UV flux at the surface due to the impact-induced ozone change are then carried out using the TUV radiative transfer model. The results suggest that mid-latitude oceanic impacts of 1 km asteroids can produce a significant, global perturbation of upper atmospheric chemistry, including multi-year global ozone depletion comparable to ozone hole records registered in the mid-1990 s. Asteroids 500 m in diameter cause limited perturbations of upper atmospheric chemistry with significant ozone depletion confined to the hemisphere in which the impact occurred. Impact-induced ozone depletion affects UV irradiance at the Earth's surface, resulting in levels of UV-B irradiance that can be dangerous for living organisms, and, in the tropics and initial midlatitude summers, far exceed levels currently experienced anywhere on Earth.

Meteor Crater is one of the first impact structures systematically studied on Earth. Its location in arid northern Arizona has been ideal for the preservation of the structure and the surviving meteoric material. The recovery of a large amount of meteoritic material in and around the crater has allowed a rough reconstruction of the impact event: an iron object 50 m in diameter impacted the Earths surface after breaking up in the atmosphere. The details of the disruption, however, are still debated. The final crater morphology (deep, bowl-shaped crater) rules out the formation of the crater by an open or dispersed swarm of fragments, in which the ratio of swarm radius to initial projectile radius Cd is larger than 3 (the final crater results from the sum of the craters formed by individual fragments). On the other hand, the lack of significant impact melt in the crater has been used to suggest that the impactor was slowed down to 12 km/s by the atmosphere, implying significant fragmentation and fragments separation up to 4 initial radii. This paper focuses on the problem of entry and motion through the atmosphere for a possible Canyon Diablo impactor as a first but necessary step for constraining the initial conditions of the impact event which created Meteor Crater. After evaluating typical models used to investigate meteoroid disruption, such as the pancake and separated fragment models, we have carried out a series of hydrodynamic simulations using the 3D code SOVA to model the impactor flight through the atmosphere, both as a continuum object and a disrupted swarm. Our results indicate that the most probable pre-atmospheric mass of the Meteor Crater projectile was in the range of 4x10^8 to 1.2x10^9 kg (equivalent to a sphere 4666 m in diameter). During the entry process the projectile lost probably 30% to 70% of its mass, mainly because of mechanical ablation and gross fragmentation. Even in the case of a tight swarm of particles (Cd

Large meteorite impacts, such as the one that created the Vredefort structure in South Africa ~2 Ga ago, result in significant heating of the target. The temperatures achieved in these events have important implications for post-impact metamorphism as well as for the development of hydrothermal systems. To investigate the post-impact thermal evolution and the size of the Vredefort structure, we have analyzed impact-induced shock heating in numerical simulations of terrestrial impacts by projectiles of a range of sizes thought to be appropriate for creating the Vredefort structure. When compared with the extent of estimated thermal shock metamorphism observed at different locations around Vredefort, our model results support our earlier estimates that the original crater was 120­-160 km in diameter, based on comparison of predicted to observed locations of shock features. The simulations demonstrate that only limited shock heating of the target occurs outside the final crater and that the cooling time was at least 0.3 Myr but no more than 30 Myr.

Using detailed geological, petrographic, geochemical, and geographical constraints we have performed numerical modeling studies that relate the Steinheim crater (apparent diameter Da = 3.8 km), the Ries crater (Da = 24 km) in southern Germany, and the moldavite (tektite) strewn field in Bohemia and Moravia (Czech Republic), Lusatia (East Germany), and Lower Austria. The moldavite strewn field extends from ~200 to 450 km from the center of the Ries to the east-northeast forming a fan with an angle of ~57°. An oblique impact of a binary asteroid from a west-southwest direction appears to explain the locations of the craters and the formation and distribution of the moldavites. The impactor must have been a binary asteroid with two widely separated components (some 1.5 and 0.15 km in diameter, respectively). We carried out a series of three-dimensional hydrocode simulations of a Ries-type impact. The results confirm previous results suggesting that impacts around 30-50° (from the horizontal) are the most favorable angles for near-surface melting, and, consequently for the formation of tektites. Finally, modeling of the motion of impact-produced tektite particles through the atmosphere produces, in the downrange direction, a narrow-angle distribution of the moldavites tektites in a fan like field with an angle of ~75°. An additional result of modeling the motion of melt inside and outside the crater is the preferred flow of melt from the main melt zone of the crystalline basement downrange towards the east-northeast rim. This explains perfectly the occurrence of coherent impact melt bodies (some tens of meters in size) in a restricted zone of the downrange rim of the Ries crater. The origin of these melt bodies, which represent chemically a mixture of crystalline basement rocks similar to the main melt mass contained (as melt particles

Jupiter's moon Europa may harbor an ocean beneath its ice cover, but the composition of that ocean and the overlying ice is nearly entirely unknown. Regardless of uncertainties in models for Europa's formation, we estimate lower limits for Europa's inventory of biogenic elements (such as C, N, O, and P) by investigating the contribution to the inventory of impact events over Europa's geologic history. A series of high-resolution hydrocode simulations were carried out over a range of comet densities (1.1, 0.8, and 0.6 g/cm3, corresponding to porosities between 0 and 45%) and impact velocities (16, 21.5, 26.5, and 30.5 km/s). We found that at typical impact velocities on Europa most impactor material reaches escape velocity, and it is assumed to be lost from Europa. For a nonporous comet, some fraction (20% or higher) of the projectile is retained by Europa even at the highest impact velocity modeled, 30.5 km/s. For porous comets, however, a significant fraction of the projectile (above 25%) is retained only for the lowest impact velocity modeled, 16 km/s. Integrated over solar system history, this suggests that 1 to 10 Gt of carbon could have been successfully delivered to Europa's surface by impacts of large comets (around 1 km in diameter). This is a few times more carbon than is contained in the procaryotic biomass of the upper 200 meters of the Earth's oceans, but about 2 orders of magnitude less if the whole depth of the oceans is considered. Therefore, regardless of its initial formation conditions, Europa should have a substantial inventory of “biogenic” elements, with implications for the chemistry of its oceans, ice cover, and the possibility of life.

Hydrocode modeling is a fundamental tool for the study of melt production in planetary impact events. Until recently, however, numerical modeling of impacts for melt production studies has been limited to vertical impacts. We present the first results of the investigation of melt production in oblique impacts. Simulations were carried out using Sandia's three-dimensional hydrocode CTH, coupled to the SESAME equation of state. While keeping other impact parameters constant, the calculations span impact angles (measured from the surface) from 90° (vertical impact) to 15°. The results show that impact angle affects the strength and distribution of the shock wave generated in the impact. As a result, both the isobaric core and the regions of melting in the target appear asymmetric and concentrated in the downrange, shallower portion of the target. The use of a pressure-decay power law (which describes pressure as function of linear distance from the impact point) to reconstruct the region of melting and vaporization is therefore complicated by the asymmetry of the shock wave. As an analog to the pressure decay versus distance from the impact point, we used a “volumetric pressure decay,” where the pressure decay is modeled as a function of volume of target material shocked at or above the given shock pressure. We find that the volumetric pressure decay exponent is almost constant for impact angles from 90° to 30°, dropping by about a factor of two for a 15° impact. In the range of shock pressures at which most materials of geologic interest melt or begin to vaporize, we find that the volume of impact melt decreases by at most 20% for impacts from 90° down to 45°. Below 45°, however, the amount of melt in the target decreases rapidly with impact angle. Compared to the vertical case, the reduction in volume of melt is about 50% for impacts at 30° and more than 90% for a 15° impact. These estimates do not include possible melting due to shear heating, which can contribute to the amount of melt production especially in very oblique impacts. Studies of melt production in vertical impacts suggest an energy scaling law in agreement with the point source limit. An energy scaling law, however, does not seem to hold for oblique impacts, even when the impact velocity is substituted by its vertical component. However, we find that for impact angles between about 30° and 90° (a range that includes 75% of impact events on planetary surfaces) the volume of melt is directly proportional to the volume of the transient crater generated by the impact.

Natural impacts in which the projectile strikes the target vertically are virtually nonexistent. Nevertheless, our inherent drive to simplify nature often causes us to suppose most impacts are nearly vertical. Recent theoretical, observational, and experimental work is improving this situation, but even with the current wealth of studies on impact cratering, the effect of impact angle on the final crater is not well understood. Although craters’ rims may appear circular down to low impact angles, the distribution of ejecta around the crater is more sensitive to the angle of impact and currently serves as the best guide to obliquity of impacts. Experimental studies established that crater dimensions depend only on the vertical component of the impact velocity. The shock wave generated by the impact weakens with decreasing impact angle. As a result, melting and vaporization depend on impact angle; however, these processes do not seem to depend on the vertical component of the velocity alone. Finally, obliquity influences the fate of the projectile: in particular, the amount and velocity of ricochet are a strong function of impact angle.

— 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 sin2θ, 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 sin3/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.

— A significant fraction of the Earth's prebiotic volatile inventory may have been delivered by asteroidal and cometary impacts during the period of heavy bombardment. The realization that comets are particularly rich in organic material seemed to strengthen this suggestion. Previous modeling studies, however, indicated that most organics would be entirely destroyed in large comet and asteroid impacts. The availability of new kinetic parameters for the thermal degradation of amino acids in the solid phase made it possible to readdress this question. We present the results of new high-resolution hydrocode simulations of asteroid and comet impact coupled with recent experimental data for amino acid pyrolysis in the solid phase. Differences due to impact velocity as well as projectile material have been investigated. Effects of angle of impacts were also addressed. The results suggest that some amino acids would survive the shock heating of large (kilometer-radius) cometary impacts. At the time of the origins of life on Earth, the steady-state oceanic concentration of certain amino acids (like aspartic and glutamic acid) delivered by comets could have equaled or substantially exceeded concentrations due to Miller-Urey synthesis in a CO2-rich atmosphere. Furthermore, in the unlikely case of a grazing impact (impact angle ∼5° from the horizontal), an amount of some amino acids comparable to that due to the background steady-state production or delivery would be delivered to the early Earth.

Two main types of material survive from the Canyon Diablo impactor, which produced Meteor Crater in Arizona: iron meteorites, which did not melt during the impact; and spheroids, which did. Ultrasensitive measurements using accelerator mass spectrometry show that the meteorites contain about seven times as much nickel-59 as the spheroids. Lower average nickel-59 contents in the spheroids indicate that they typically came from 0.5 to 1 meter deeper in the impactor than did the meteorites. Numerical modeling for an impact velocity of 20 kilometers per second shows that a shell 1.5 to 2 meters thick, corresponding to 16 percent of the projectile volume, remained solid on the rear surface; that most of the projectile melted; and that little, if any, vaporized.

Since the confirmation that the buried Chicxulub structure is the long-sought K/T boundary crater, numerous efforts have been devoted to modeling the impact event and estimating the amount of target material that underwent melting and vaporization. Previous hydrocode simulations modeled the Chicxulub event as a vertical impact. We carried out a series of three-dimensional (3D) hydrocode simulations of the Chicxulub impact event to study how the impact angle affects the results of impact events. The simulations model an asteroid, 10 km in diameter, impacting at 20 km/s on a target resembling the lithology of the Chicxulub site. The angles of impact modeled are 90° (vertical), 60°, 45°, 30°, and 15°. We find that the amount of sediments (surface layer) vaporized in the impact reaches a maximum for an impact angle of 30° from the surface, corresponding to less than two times the amount of vaporization for the vertical case. The degassing of the sedimentary layer, however, drops abruptly for a 15° impact angle. The amount of continental crust melted in the impact decreases monotonically (a consequence of the decrease in the maximum depth of melting) from the vertical impact case to the 15° impact. Melting and vaporization occur primarily in the downrange direction for oblique impacts, due to asymmetries in the strength of the shock wave with respect to the point of impact. The results can be used to scale the information from the available vertical simulations to correct for the angle of impact. A comparison of a 3D vertical impact simulation with a similar two-dimensional (2D) simulation shows good agreement between vertical 3D and 2D simulations.

We constructed a numerical model of the Chicxulub impact event using the Chart-D Squared (CSQ) code coupled with the ANalytic Equation Of State (ANEOS) package. In the simulations we utilized a target stratigraphy based on borehole data and employed newly developed equations of state for the materials that are believed to play a crucial role in the impact-related extinction hypothesis: carbonates (calcite) and evaporites (anhydrite). Simulations explored the effects of different projectile sizes (10 to 30 km in diameter) and porosity (0 to 50%). The effect of impact speed is addressed by doing simulations of asteroid impacts (vi = 20 km/s) and comet impacts (vi = 50 km/s). The masses of climatically important species injected into the upper atmosphere by the impact increase with the energy of the impact event, ranging from 350 to 3500 Gt for CO2, from 40 to 560 Gt for S, and from 200 to 1400 Gt for water vapor. While our results are in good agreement with those of Ivanov et al. [1996], our estimated CO2 production is 1 to 2 orders of magnitude lower than the results of Takata and Ahrens [1994], indicating that the impact event enhanced the end-Cretaceous atmospheric CO2 inventory by, at most, 40%. Consequently, sulfur may have been the most important climatically active gas injected into the stratosphere. The amount of S released by the impact is several orders of magnitude higher than any known volcanic eruption and, with H2O, is high enough to produce a sudden and significant perturbation of Earth's climate.

— The Vredefort structure in South Africa was created by a meteorite impact about two billion years ago. Since that time, the crater has been deeply eroded; so to estimate its original size, researchers have had to rely heavily upon comparison to other terrestrial impact structures. Recent estimates of the original crater diameter range from 160 km to as much as 400 km. In this study, we combined the capabilities of both hydrocode and finite-element modeling, using the former to predict where the pressure of an impact-generated shock wave would have been high enough to form planar deformation features (PDFs) and shatter cones and the latter to follow the subsequent displacement of these shock isobars during the collapse of the crater. We established constraints on the sizes of the projectile and the transient crater (and, thus, on the size of the final crater) by comparing the observed locations of PDFs around Vredefort to the results of our simulations of impacts by projectiles of various sizes. These simulations indicate that a rocky projectile with a diameter of ∼10 km, impacting vertically at a velocity of 20 km/s generates shock pressures that are consistent with the distribution of PDFs around Vredefort. These projectile parameters correspond to a transient crater ∼80 km in diameter or a final crater ∼120–160 km in diameter. Allowing for uncertainties in our modeling procedures, we consider final craters 120 to 200 km in diameter to be consistent with the observed locations of PDFs at Vredefort. The shock pressure contour corresponding to the formation of shatter cones is almost horizontal near the surface, making the locations of these features less useful constraints on the crater size. However, they may provide a constraint on the amount of erosion that has occurred since the impact.

The production of melt and vapor is an important process in impact cratering events. Because significant melting and vaporization do not occur in impacts at velocities currently achievable in the laboratory, a detailed study of the production of melt and vapor in planetary impact events is carried out with hydrocode simulations. Sandia's two-dimensional axisymmetric hydrocode CSQ was used to estimate the amount of melt and vapor produced for widely varying initial conditions: 10 to 80 km/sec for impact velocity, 0.2 to 10 km for the projectile radius. Runs with different materials demonstrate the material dependency of the final result. These results should apply to any size projectile (for given impact velocity and material), since the results can be dynamically scaled so long as gravity is unimportant in affecting the early-time flow. In contrast with the assumptions of previous analytical models, a clear difference in shape, impact-size dependence, and depth of burial has been found between the melt regions and the isobaric core. In particular, the depth of the isobaric core is not a good representation of the depth of the melt regions, which form deeper in the target. While near-surface effects cause the computed melt region shapes to look like “squashed spheres” the spherical shape is still a good analytical analog. One of the goals of melt production studies is to find proper scaling laws to infer melt production for any impact event of interest. We tested the point source limit scaling law for melt volumes (μ = 0.55–0.6) proposed by M. D. Bjorkman and K. A. Holsapple (1987,Int. J. Impact Eng.5, 155–163). Our results indicate that the point source limit concept does not apply to melt and vapor production. Rather, melt and vapor production follows an energy scaling law (μ = 0.67), in good agreement with previous results of T. J. Ahrens and J. D. O'Keefe [1977, inImpact and Explosion Cratering(D. J. Roddy, R. O. Pepin, and R. B. Merrill, Eds.), pp. 639–656, Pergamon Press, Elmsford, NY]. Finally we tested the accuracy of our melt production calculation against a terrestrial dataset compiled by R. A. F. Grieve and M. J. Cintala (1992,Meteorities27, 526–538). The hydrocode melt volumes are in good agreement with the estimated volumes of that set of terrestrial craters on crystalline basements. At present there is no good model for melt production from impact craters on sedimentary targets.

The Fall of the Dinosaurs According to the fossil record, the rule of dinosaurs came to an abrupt end ∼65 million years ago, when all nonavian dinosaurs and flying reptiles disappeared. Several possible mechanisms have been suggested for this mass extinction, including a large asteroid impact and major flood volcanism. Schulte et al. (p. 1214 ) review how the occurrence and global distribution of a global iridium-rich deposit and impact ejecta support the hypothesis that a single asteroid impact at Chicxulub, Mexico, triggered the extinction event. Such an impact would have instantly caused devastating shock waves, a large heat pulse, and tsunamis around the globe. Moreover, the release of high quantities of dust, debris, and gases would have resulted in a prolonged cooling of Earth's surface, low light levels, and ocean acidification that would have decimated primary producers including phytoplankton and algae, as well as those species reliant upon them.

The Letters by Archibald et al. , Keller et al. , and Courtillot and Fluteau question our conclusion that the Cretaceous-Paleogene mass extinction was caused by the asteroid impact at Chicxulub. All three Letters stress that Deccan flood basalt volcanism played a major role in the extinction. Keller

THE LETTERS BY ARCHIBALD ET AL., KELLER ET al., and Courtillot and Fluteau question our conclusion that the Cretaceous-Paleogene mass extinction was caused by the asteroid impact at Chicxulub. All three Letters stress that Deccan flood basalt volcanism played a major role in the extinction. Keller et al. and Archibald et al. also mention that climate change was a factor, and Archibald et al. We disagree with the hypothesis that vol- canic activity can explain the extinction.

A systematic comparison of a variable cosmogenic-production rate model (proposed in a previous paper) with the conventional constant-production rate model is carried out. Attention is focussed on the time-integrated concentrations in meteorites. A graphical method for the estimate of the exposure and terrestrial ages is applied with reference to the ten pairs of the five long-lived cosmogenic nuclides36Cl,26Al,10Be,53Mn and129I of interest in the Accelerator Mass Spectrometry.

The cratering event produced by the Deep Impact mission is a unique experimental opportunity, beyond the capability of Earth-based laboratories with regard to the impacting energy, target material, space environment, and extremely low-gravity field. Consequently, impact cratering theory and modeling play an important role in this mission, from initial inception to final data analysis. Experimentally derived impact cratering scaling laws provide us with our best estimates for the crater diameter, depth, and formation time: critical in the mission planning stage for producing the flight plan and instrument specifications. Cratering theory has strongly influenced the impactor design, producing a probe that should produce the largest possible crater on the surface of Tempel 1 under a wide range of scenarios. Numerical hydrocode modeling allows us to estimate the volume and thermodynamic characteristics of the material vaporized in the early stages of the impact. Hydrocode modeling will also aid us in understanding the observed crater excavation process, especially in the area of impacts into porous materials. Finally, experimentally derived ejecta scaling laws and modeling provide us with a means to predict and analyze the observed behavior of the material launched from the comet during crater excavation, and may provide us with a unique means of estimating the magnitude of the comet’s gravity field and by extension the mass and density of comet Tempel 1.

Mars is the most Earth-like of the Solar System s planets, and the first place to look for any sign of present or past extraterrestrial life. Its surface shows many features indicative of the presence of surface and sub-surface water, while impact cratering and volcanism have provided temporary and local surface heat sources throughout Mars geologic history. Impact craters are widely used ubiquitous indicators for the presence of sub-surface water or ice on Mars. In particular, the presence of significant amounts of ground ice or water would cause impact-induced hydrothermal alteration at Martian impact sites. The realization that hydrothermal systems are possible sites for the origin and early evolution of life on Earth has given rise to the hypothesis that hydrothermal systems may have had the same role on Mars. Rough estimates of the heat generated in impact events have been based on scaling relations, or thermal data based on terrestrial impacts on crystalline basements. Preliminary studies also suggest that melt sheets and target uplift are equally important heat sources for the development of a hydrothermal system, while its lifetime depends on the volume and cooling rate of the heat source, as well as the permeability of the host rocks. We present initial results of two-dimensional (2D) and three-dimensional (3D) simulations of impacts on Mars aimed at constraining the initial conditions for modeling the onset and evolution of a hydrothermal system on the red planet. Simulations of the early stages of impact cratering provide an estimate of the amount of shock melting and the pressure-temperature distribution in the target caused by various impacts on the Martian surface. Modeling of the late stage of crater collapse is necessary to characterize the final thermal state of the target, including crater uplift, and distribution of the heated target material (including the melt pool) and hot ejecta around the crater.

Numerical modeling is a fundamental tool for understanding the dynamics of impact cratering, especially at planetary scales. In particular, processes like melting/vaporization and crater collapse, typical of planetary-scale impacts, are not reproduced in the laboratory, and can only be investigated by numerical modeling. The continuum dynamics of impact cratering events is fairly well understood and implemented in numerical codes; however, the response of materials to shocks is governed by specific material properties. Accurate material models are thus crucial for realistic simulation of impact cratering, and still represent one of the major problems associated with numerical modeling of impacts.

[1] Crustal magnetic fields on Mars are inhomogeneously distributed with the strongest fields occurring over the southern highlands in a longitude sector between approximately 130°E and 240°E. Using analytic approximations and empirical scaling laws, it is estimated that much of the weakly magnetized southern highlands (i.e., that between 110°W and 130°E) was shocked to pressures exceeding 1–2 GPa during the Hellas and Argyre impacts. Possible primary remanence carriers in the martian crust include iron oxides and iron sulfides (pyrrhotite). If pyrrhotite is the main remanence carrier, extensive demagnetization of crustal regions (∼90%) may occur at shock pressures of 2 GPa or more. Thus, at least for this remanence carrier, impact shock demagnetization can potentially explain the distribution of crustal fields in the southern highlands.

Impact glass associated with 11 elongate depressions in the Pampean Plain of Argentina, north of the city of Rı́o Cuarto, was suggested to be proximal ejecta related to a highly oblique impact event. We have identified about 400 additional elongate features in the area that indicate an aeolian, rather than an impact, origin. We have also dated fragments of glass found at the Rı́o Cuarto depressions; the age is similar to that of glass recovered 800 kilometers to the southeast. This material may be tektite glass from an impact event around 0.48 million years ago, representing a new tektite strewn field.

Several impact craters on Jupiter's satellite Europa exhibit central peaks. On the terrestrial planets, central peaks consist of fractured but competent rock uplifted during cratering. Therefore, the observation of central peaks on Europa indicates that an ice layer must be sufficiently thick that the impact events did not completely penetrate it. We conducted numerical simulations of vapor and melt production during cratering of water ice layers overlying liquid water to estimate the thickness of Europa's icy crust. Because impacts disrupt material well beyond the zone of partial melting, our simulations put a lower limit on ice thickness at the locations and times of impact. We conclude that the ice must be more than 3 to 4 kilometers thick.

Proceedings of the workshop, Impact cratering: Bridging thegap between modeling and observations

Given that impact cratering is the most common geological process in the solar system, study of this topic provides an excellent means of interplanetary comparison. One can gain insight into the crustal properties of various solar system bodies, the population of asteroids and comets that potentially could cause impacts, and the varied ways in which impact cratering has influenced the geologic histories of the planets. The First International Conference on Impact Cratering in the Solar System recently was held at the 40th ESLAB (European Space Laboratory) Symposium at the European Space Agency's European Space Research and Technology Centre, in Noordwijk, the Netherlands. About 150 participants from around the world came to discuss the latest understanding of and issues related to impact cratering in the solar system. The uniqueness of this conference lay in its broad scope, covering everything from the late heavy bombardment controversy to impact-induced mass extinctions. Thus, the conference brought together many scientists of varying expertise who may otherwise have never interacted.

Formation of a planetary-scale impact crater hundreds of meters to hundreds of kilometers in diameter involves the interplay of multiphase processes occurring over size- and time-scales that range over many orders of magnitude. To further complicate matters, a hypervelocity impact into geologic materials has never been observed at greater than laboratory scales. Understanding the formation process clearly is a very complicated problem that requires a highly interdisciplinary approach. Significant work has been done recently in several key areas of impact studies, but in many respects, there is a “disconnect” among groups employing different approaches. This pertains, in particular, to modeling versus observations.

One of the most fundamental geologic processes that have shaped the surfaces of the solar system's planetary bodies is the impact of solid bodies. It affects the surface of a planetary body not only by creating a visible crater, but also by modifying, locally and temporally, the crater subsurface and a wide area around the excavated crater. In particular, terrestrial impact structures are in the unique position of providing the only ground truth available for investigating the impact cratering process. The establishment of shock metamorphic effects on rocks as reliable criteria for assigning impact origins to terrestrial structures has led to the recognition of over 160 impact structures on Earth over the past 40 years. Investigation of these structures is fundamental for validating theoretical and experimental impact studies.

The Δ 14C series of yearly sampled cosmogenic 14C in wines (1909–1952) was analyzed with the Fourier Vector Amplitude (FVA) method, using cyclograms as a graphic tool, to find information on periodicities imprinted by the sun. Because the high sensitivity of the FVA algorithm in detecting periodicities and their variations is emphasized by immediate visualization of its cyclograms, a suggestion has been found for a modulation event. Data are compared with a frequency modulation model, the extension of which to the long 9000–a Δ 14C series suggests a first approach to interpret the Suess wiggles.

We use one-dimensional (1D) atmospheric models coupled to a sulfate aerosol model to investigate climate forcing and short-term response to stratospheric sulfate aerosols produced by the reaction of S-bearing gases and water vapor released in the Chicxulub impact event. A 1D radiation model is used to assess the climate forcing due to the impact-related loading of S-bearing gases. The model suggests that a climate forcing 100 times larger than that from the Pinatubo volcanic eruption is associated with the Chicxulub impact event for at least 2 years after the impact. In particular, we find a saturation effect in the forcing, that is, there is no significant difference in the maximum forcing between the highest (approximately 300 Gt) and lowest (approximately 30 Gt) estimated stratospheric S-loading from the Chicxulub impact. However, higher S-loads increase the overall duration of the forcing by several months. We use a single column model for a preliminary investigation of the short-term climate response to the impact-related production of sulfate aerosols (the lack of horizontal feedbacks limits the usefulness of the single column model to the first few days after the impact). Compared with the present steady-state climate, the introduction of large amounts of sulfate aerosols in the stratosphere results in a significant cooling of the Earth's surface. A long-term climate response can only be investigated with the use of a three-dimensional atmospheric model, which allows for the atmospheric circulation to adjust to the perturbation. Overall, although the climate perturbation to the forcing appears to be relatively large, the geologic record shows no sign of a significant long-term climatic shift across the K/T boundary, which is indicative of a fast post-impact climatic recovery.

Over the last few decades, rapid improvement of computer capabilities has allowed impact cratering to be modeled with increasing complexity and realism, and has paved the way for a new era of numerical modeling of the impact process, including full, three-dimensional (3D) simulations. When properly benchmarked and validated against observation, computer models offer a powerful tool for understanding the mechanics of impact crater formation. This work presents results from the first phase of a project to benchmark and validate shock codes. A variety of 2D and 3D codes were used in this study, from commercial products like AUTODYN, to codes developed within the scientific community like SOVA, SPH, ZEUS-MP, iSALE, and codes developed at U.S. National Laboratories like CTH, SAGE/RAGE, and ALE3D. Benchmark calculations of shock wave propagation in aluminum-on-aluminum impacts were performed to examine the agreement between codes for simple idealized problems. The benchmark simulations show that variability in code results is to be expected due to differences in the underlying solution algorithm of each code, artificial stability parameters, spatial and temporal resolution, and material models. Overall, the inter-code variability in peak shock pressure as a function of distance is around 10 to 20%. In general, if the impactor is resolved by at least 20 cells across its radius, the underestimation of peak shock pressure due to spatial resolution is less than 10%. In addition to the benchmark tests, three validation tests were performed to examine the ability of the codes to reproduce the time evolution of crater radius and depth observed in vertical laboratory impacts in water and two well-characterized aluminum alloys. Results from these calculations are in good agreement with experiments. There appears to be a general tendency of shock physics codes to underestimate the radius of the forming crater. Overall, the discrepancy between the model and experiment results is between 10 and 20%, similar to the inter-code variability.