USGS Science⌕ Search

SEARCH · USGS Science

Results for “Shock Waves”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
USGS2003

Shock-wave-induced fracturing of calcareous nannofossils from the Chesapeake Bay impact crater

Fractured calcareous nannofossils of the genus Discoaster from synimpact sediments within the Chesapeake Bay impact crater demonstrate that other petrographic shock indicators exist for the cratering process in addition to quartz minerals. Evidence for shock-induced taphonomy includes marginal fracturing of rosette-shaped Discoaster species into pentagonal shapes and pressure- and temperature-induced dissolution of ray tips and edges of discoasters. Rotational deformation of individual crystallites may be the mechanism that produces the fracture pattern. Shock-wave-fractured calcareous nannofossils were recovered from synimpact matrix material representing tsunami or resurge sedimentation that followed impact. Samples taken from cohesive clasts within the crater rubble show no evidence of shock-induced fracturing. The data presented here support growing evidence that microfossils can be used to determine the intensity and timing of wet-impact cratering.

Geology↗
USGS1965

Stishovite: Synthesis by shock wave

Small amounts of stishovite were separated from specimens of explosively shocked sandstones, novaculite, and single-crystal quartz. Estimated peak pressures for the syntheses ranged from 150 to 280 kilobars, and shock temperatures were from 150° to 900°C. No coesite was detected in any sample. It is suggested that quartz can invert during shock to a short-range-order phase, with sixfold coordination. A small portion of this phase may develop the long-range order of stishovite, and, during the more protracted decrease of the pressure pulse through the stability field of coesite accompanying meteorite crater formation, a portion may invert to coesite.

Science↗
USGS1992

Seismic excitation by space shuttles

Shock waves generated by the space shuttles Columbia (August 13, 1989), Atlantis (April 11, 1991) and Discovery (September 18, 1991) on their return to Edwards Air Force Base, California, were recorded by TERRAscope (Caltech's broadband seismic network), the Caltech-U.S.G.S Southern California Seismic Network (SCSN), and the University of Southern California (USC) Los Angeles Basin Seismic Network. The spatial pattern of the arrival times exhibits hyperbolic shock fronts from which the path, velocity and altitude of the space shuttle could be determined. The shock wave was acoustically coupled to the ground, converted to a seismic wave, and recorded clearly at the broadband TERRAscope stations. The acoustic coupling occurred very differently depending on the conditions of the Earth's surface surrounding the station. For a seismic station located on hard bedrock, the shock wave (N wave) was clearly recorded with little distortion. Aside from the N wave, very little acoustic coupling of the shock wave energy to the ground occurred at these sites. The observed N wave record was used to estimate the overpressure of the shock wave accurately; a pressure change of 0.5 to 2.2 mbars was obtained. For a seismic station located close to the ocean or soft sedimentary basins, a significant amount of shock wave energy was transferred to the ground through acoustic coupling of the shock wave and the oceanic Rayleigh wave. A distinct topography such as a mountain range was found effective to couple the shock wave energy to the ground. Shock wave energy was also coupled to the ground very effectively through large man made structures such as high rise buildings and offshore oil drilling platforms. For the space shuttle Columbia, in particular, a distinct pulse having a period of about 2 to 3 seconds was observed, 12.5 s before the shock wave, with a broadband seismograph in Pasadena. This pulse was probably excited by the high rise buildings in downtown Los Angeles which were simultaneously hit by the space shuttle shock waves. The proximity of the natural periods of the high rise buildings and the modal periods of the Los Angeles basin enabled efficient energy transfer from shock wave to seismic wave. ?? 1992 Springer-Verlag.

Shock Waves↗
USGS1991

Seismic excitation by the space shuttle Columbia

Seismic stations in southern California recorded the atmospheric shock waves generated by the space shuttle Columbia on its return to the Edwards Air Force base on 13 August 1989 (Fig. 1). In addition to the shock wave, the broad-band IRIS–TERRAscope station at Pasadena recorded a distinct pulse with a period of ~2–3 seconds, which arrived 12.5 seconds before the shock wave (Fig. 2). This pulse was also recorded at the University of Southern California, near downtown Los Angeles, where it arrived 3 seconds after the shock wave. The origin of this pulse could not be readily identified. We show here that it was a seismic P wave excited by the motion of high-rise buildings in downtown Los Angeles, which were hit by the shock wave. The proximity of the natural period of the high-rise buildings to that of the Los Angeles basin enabled efficient energy transfer from shock wave to seismic wave.

California↗
USGS1997

A numerical investigation of choked flow dynamics and its application to the triggering mechanism of long-period events at Redoubt Volcano, Alaska

We use numerical simulations of transonic flow through a crack to study the dynamics of the formation of shock waves downstream from a nozzle-like constriction inside the crack. The model solves the full set of Navier-Stokes equations in two dimensions via an explicit multifield finite difference representation. The crack walls are assumed to be perfectly rigid, and elastic coupling to the solid is not considered. The simulations demonstrate how the behavior of unsteady shock waves near the walls can produce recurring step-like pressure transients in the flow, which in turn induce resonance of the fluid-filled crack. The motion of the shock waves is governed primarily by smooth, low-amplitude pressure fluctuations at the outlet of the crack. The force induced on the walls scales with the amplitude of the shock, which is a function of the magnitude of the inlet pressure, aperture of the constriction, and thickness of the boundary layer. The applied force also scales in proportion to the spatial extent of the shock excursion, which depends on the fluctuation rate of outlet pressure. Using the source parameters of long-period (LP) events at Redoubt Volcano, Alaska, as a guide for our simulations, we infer that coupling of the shock to the walls occurs for crack inlet to outlet pressure ratios p i /p o >2.31 and that the position of the shock front becomes most sensitive to outlet pressure fluctuations for flow regimes with p i /p o >2.48. For such regimes, fluctuations of outlet pressure of up to ±0.5 MPa at rates up to 3 MPa/s are sufficient to induce pressure transients with magnitudes up to 12.5 MPa over 0.1–2.5 m of the walls within ∼0.5 s. These flow parameters may be adequate for triggering the LP events in the precursory swarm to the December 14, 1989, eruption of Redoubt. According to the flow model the recurrence rate and amplitudes of L.P events are inferred to be a manifestation of the response of a shallow hydrothermal reservoir to the sustained injection of superheated steam from a magma column roofing below this reservoir.

Journal of Geophysical Research B: Solid Earth↗
USGS1979

Probable hydrologic effects of a hypothetical failure of Mackay Dam on the Big Lost River Valley from Mackay, Idaho to the Idaho National Engineering Laboratory

Mackay Dam is an irrigation reservoir on the Big Lost River, Idaho, approximately 7.2 kilometers northwest of Mackay, Idaho. Consequences of possible rupture of the dam have long concerned the residents of the river valley. The presence of reactors and of a management complex for nuclear wastes on the reservation of the Idaho National Engineering Laboratory (INEL), near the river , give additional cause for concern over the consequences of a rupture of Mackay Dam. The objective of this report is to calculate and route the flood wave resulting from the hypothetical failure of Mackay Dam downstream to the INEL. Both a full and a 50 percent partial breach of this dam are investigated. Two techniques are used to develop the dam-break model. The method of characteristics is used to propagate the shock wave after the dam fails. The linear implicit finite-difference solution is used to route the flood wave after the shock wave has dissipated. The time of travel of the flood wave, duration of flooding, and magnitude of the flood are determined for eight selected sites from Mackay Dam, Idaho, through the INEL diversion. At 4.2 kilometers above the INEL diversion, peak discharges of 1,550.2 and 1,275 cubic meters per second and peak flood elevations of 1,550.3 and 1,550.2 meters were calculated for the full and partial breach, respectively. Flood discharges and flood peaks were not compared for the area downstream of the diversion because of the lack of detailed flood plain geometry. (Kosco-USGS)

Idaho↗
USGS1979

Elevations and discharges produced by a simulated flood wave on the lower Sabine River, Louisiana and Texas, caused by a theoretical dam failure

The Toledo Bend Reservoir is located on the lower Sabine River between Louisiana and Texas. The objective of this study was to calculate the flood wave that would result from the theoretical failure of 25 percent of Toledo Bend Dam and route the wave downstream to Orange, Tex. Computations assumed failure (1) at the peak of the 100-year flood when discharge of the Sabine River would be 102,000 cubic feet per second and (2) when the average discharge would be 10,000 cubic feet per second. Two techniques were used in the dam-break model. The method of characteristics was used to propagate the shock wave after the dam fails. The linear implicit finite-difference solution was used to route the flood wave after the shock wave has dissipated. The magnitude of the flood was determined for sites near Burkeville, Bon Wier, Ruliff, and Orange, Tex., along the lower Sabine River. For these sites, respectively, the following peak elevations were calculated: 119, 82, 31, and 13 feet for the 100-year flood and 110, 75, 27, and 9 feet for the average discharge.

Louisiana, Texas↗
USGS1980

Geophysical investigations in deep horizontal holes drilled ahead of tunnelling

Deep horizontal drill holes have been used since 1967 by the Defense Nuclear Agency as a primary exploration tool for siting nuclear events in tunnels at the Nevada Test Site. The U.S. Geological Survey had developed geophysical logging techniques for obtaining resistivity and velocity in these holes, and to date 33 horizontal drill holes in excess of 300 m in depth have been successfully logged. The deepest hole was drilled to a horizontal depth of 1125 m. The purposes of the logging measurements are to define clay zones, because of the unstable ground conditions such zones can present to tunnelling, and to define zones of partially saturated rock, because of the attenuating effects such zones have on the shock wave generated by the nuclear detonation. Excessive attenuation is undesirable because the shock wave is used as a tunnel closure mechanism to contain debris and other undesirable explosion products. Measurements are made by pumping resistivity, sonic and geophone probes down the drill string and out of the bit into the open hole. Clay zones are defined by the electrical resistivity technique based on empirical data relating the magnitude of the resistivity measurement to qualitative clay content. Rock exhibiting resistivity of less than 20 ??-m is considered potentially unstable, and resistivities less than 10 ??-m indicate appreciable amounts of clay are present in the rock. Partially saturated rock zones are defined by the measurement of the rock sound speed. Zones in the rock which exhibit velocities less than 2450 m/sec are considered of potential concern. ?? 1980.

International Journal of Rock Mechanics and Mining↗
USGS1965

Kink-bands: Shock deformation of biotite resulting from a nuclear explosion

Microscopic examination of granodiorite samples from the shock region around a nuclear explosion reveals sharply folded lens-shaped zones (kink-bands) in the mineral biotite. Fifty percent of these zones are oriented approximately 90 degrees to the direction of shock-wave propagation, but other zones are symmetrically concentrated at shear angles of 50 degrees and 70 degrees to the direction of shock-wave propagation.

Science↗
USGS1987

Computer simulations of large asteroid impacts into oceanic and continental sites--Preliminary results on atmospheric, cratering and ejecta dynamics

Computer simulations have been completed that describe passage of a 10-km-diameter asteroid through the Earth's atmosphere and the subsequent cratering and ejecta dynamics caused by impact of the asteroid into both oceanic and continental sites. The asteroid was modeled as a spherical body moving vertically at 20 km/s with a kinetic energy of 2.6 × 10 30 ergs (6.2 × 10 7 Mt ). Detailed material modeling of the asteroid, ocean, crustal units, sedimentary unit, and mantle included effects of strength and fracturing, generic asteroid and rock properties, porosity, saturation, lithostatic stresses, and geothermal contributions, each selected to simulate impact and geologic conditions that were as realistic as possible. Calculation of the passage of the asteroid through a U.S. Standard Atmosphere showed development of a strong bow shock wave followed by a highly shock compressed and heated air mass. Rapid expansion of this shocked air created a large low-density region that also expanded away from the impact area. Shock temperatures in air reached ∼20, 000K near the surface of the uplifting crater rim and were as high as ∼2000K at more than 30 km range and 10 km altitude. Calculations to 30 s showed that the shock fronts in the air and in most of the expanding shocked air mass preceded the formation of the crater, ejecta, and rim uplift and did not interact with them. As cratering developed, uplifted rim and target material were ejected into the very low density, shock-heated air immediately above the forming crater, and complex interactions could be expected. Calculations of the impact events showed equally dramatic effects on the oceanic and continental targets through an interval of 120 s. Despite geologic differences in the targets, both cratering events developed comparable dynamic flow fields and by ∼29s had formed similar-sized transient craters ∼39km deep and ∼62km across. Transient-rim uplift of ocean and crust reached a maximum altitude of nearly 40 km at ∼30s and began to decay at velocities of 500 m/s to develop large-tsunami conditions. After ∼30s, strong gravitational rebound drove both craters toward broad flat-floored shapes. At 120 s, transient crater diameters were ∼80km (continental) and ∼105km (oceanic) and transient depths were ∼27km; crater floors consisting of melted and fragmented hot rock were rebounding rapidly upward. By 60 s, the continental crater had ejected ∼2 × 10 14 t, about twice the mass ejected from the oceanic crater. By 120 s, ∼70, 000km 3 (continental) and ∼90, 000km 3 (oceanic) target material were excavated (no mantle) and massive ejecta blankets were formed around the craters. We estimate that in excess of ∼70% of the ejecta would finally lie within ∼3 crater diameters of the impact, and the remaining ejecta (∼10 13 t), including the vaporized asteroid, would be ejected into the atmosphere to altitudes as high as the ionosphere. Effects of secondary volcanism and return of the ocean over hot oceanic crater floor could also be expected to contribute substantial material to the atmosphere.

International Journal of Impact Engineering↗
USGS1996

Distinguishing shocked from tectonically deformed quartz by the use of the SEM and chemical etching

Multiple sets of crystallographically-oriented planar deformation features (PDFs) are generated by high-strain-rate shock waves at pressures of > 12 GPa in naturally shocked quartz samples. On surfaces, PDFs appear as narrow (50-500 nm) lamellae filled with amorphosed quartz (diaplectic glass) which can be etched with hydrofluoric acid or with hydrothermal alkaline solutions. In contrast, slow-strain-rate tectonic deformation pressure produces wider, semi-linear and widely spaced arrays of dislocation loops that are not glass filled. Etching samples with HF before examination in a scanning electron microscope (SEM) allows for unambiguous visual distinction between glass-filled PDFs and glass-free tectonic deformation arrays in quartz. This etching also reveals the internal 'pillaring' often characteristic of shock-induced PDFs. This technique is useful for easily distinguishing between shock and tectonic deformation in quartz, but does not replace optical techniques for characterizing the shock features.

Earth and Planetary Science Letters↗
USGS2007

Directed blasts and blast-generated pyroclastic density currents: a comparison of the Bezymianny 1956, Mount St Helens 1980, and Soufrière Hills, Montserrat 1997 eruptions and deposits

We compare eruptive dynamics, effects and deposits of the Bezymianny 1956 (BZ), Mount St Helens 1980 (MSH), and Soufrière Hills volcano, Montserrat 1997 (SHV) eruptions, the key events of which included powerful directed blasts. Each blast subsequently generated a high-energy stratified pyroclastic density current (PDC) with a high speed at onset. The blasts were triggered by rapid unloading of an extruding or intruding shallow magma body (lava dome and/or cryptodome) of andesitic or dacitic composition. The unloading was caused by sector failures of the volcanic edifices, with respective volumes for BZ, MSH, and SHV c. 0.5, 2.5, and 0.05 km 3 . The blasts devastated approximately elliptical areas, axial directions of which coincided with the directions of sector failures. We separate the transient directed blast phenomenon into three main parts, the burst phase, the collapse phase, and the PDC phase. In the burst phase the pressurized mixture is driven by initial kinetic energy and expands rapidly into the atmosphere, with much of the expansion having an initially lateral component. The erupted material fails to mix with sufficient air to form a buoyant column, but in the collapse phase, falls beyond the source as an inclined fountain, and thereafter generates a PDC moving parallel to the ground surface. It is possible for the burst phase to comprise an overpressured jet, which requires injection of momentum from an orifice; however some exploding sources may have different geometry and a jet is not necessarily formed. A major unresolved question is whether the preponderance of strong damage observed in the volcanic blasts should be attributed to shock waves within an overpressured jet, or alternatively to dynamic pressures and shocks within the energetic collapse and PDC phases. Internal shock structures related to unsteady flow and compressibility effects can occur in each phase. We withhold judgment about published shock models as a primary explanation for the damage sustained at MSH until modern 3D numerical modeling is accomplished, but argue that much of the damage observed in directed blasts can be reasonably interpreted to have been caused by high dynamic pressures and clast impact loading by an inclined collapsing fountain and stratified PDC. This view is reinforced by recent modeling cited for SHV. In distal and peripheral regions, solids concentration, maximum particle size, current speed, and dynamic pressure are diminished, resulting in lesser damage and enhanced influence by local topography on the PDC. Despite the different scales of the blasts (devastated areas were respectively 500, 600, and >10 km 2 for BZ, MSH, and SHV), and some complexity involving retrogressive slide blocks and clusters of explosions, their pyroclastic deposits demonstrate strong similarity. Juvenile material composes >50% of the deposits, implying for the blasts a dominantly magmatic mechanism although hydrothermal explosions also occurred. The character of the magma fragmented by explosions (highly viscous, phenocryst-rich, variable microlite content) determined the bimodal distributions of juvenile clast density and vesicularity. Thickness of the deposits fluctuates in proximal areas but in general decreases with distance from the crater, and laterally from the axial region. The proximal stratigraphy of the blast deposits comprises four layers named A, B, C, D from bottom to top. Layer A is represented by very poorly sorted debris with admixtures of vegetation and soil, with a strongly erosive ground contact; its appearance varies at different sites due to different ground conditions at the time of the blasts. The layer reflects intense turbulent boundary shear between the basal part of the energetic head of the PDC and the substrate. Layer B exhibits relatively well-sorted fines depleted debris with some charred plant fragments; its deposition occurred by rapid suspension sedimentation in rapidly waning, high-concentration conditions. Layer C is mainly a poorly sorted massive layer enriched by fines with its uppermost part laminated, created by rapid sedimentation under moderate-concentration, weakly tractive conditions, with the uppermost laminated part reflecting a dilute depositional regime with grain-by-grain traction deposition. By analogy to laboratory experiments, mixing at the flow head of the PDC created a turbulent dilute wake above the body of a gravity current, with layer B deposited by the flow body and layer C by the wake. The uppermost layer D of fines and accretionary lapilli is an ash fallout deposit of the finest particles from the high-rising buoyant thermal plume derived from the sediment-depleted pyroclastic density current. The strong similarity among these eruptions and their deposits suggests that these cases represent similar source, transport and depositional phenomena.

Bezymianny, Mount St Helens, Soufriere Hills↗
USGS2022

Diverse tsunamigenesis triggered by the Hunga Tonga-Hunga Ha’apai eruption

On the evening of 15 January 2022, the Hunga Tonga-Hunga Ha’apai volcano 1 unleashed a violent underwater eruption, blanketing the surrounding land masses in ash and debris . The eruption generated tsunamis observed around the world. An event of this type last occurred in 1883 during the eruption of Krakatau , and thus we have the first observations of a tsunami from a large emergent volcanic eruption captured with modern instrumentation. Here we show that the explosive eruption generated waves through multiple mechanisms, including: (1) air–sea coupling with the initial and powerful shock wave radiating out from the explosion in the immediate vicinity of the eruption; (2) collapse of the water cavity created by the underwater explosion; and (3) air–sea coupling with the air-pressure pulse that circled the Earth several times, leading to a global tsunami. In the near field, tsunami impacts are strongly controlled by the water-cavity source whereas the far-field tsunami, which was unusually persistent, can be largely described by the air-pressure pulse mechanism. Catastrophic damage in some harbours in the far field was averted by just tens of centimetres, implying that a modest sea level rise combined with a future, similar event would lead to a step-function increase in impacts on infrastructure. Piecing together the complexity of this event has broad implications for coastal hazards in similar geophysical settings, suggesting a currently neglected source of global tsunamis.

Nature↗
USGS1967

Shock effects in certain rock-forming minerals

Shock effects in quartz, plagioclase, biotite, amphibole, and some accessory minerals have been observed in rocks subjected to various degrees of meta morphism by meteoritic impact. The shock features described are unique; they are never observed in rocks from normal geologic environments. Such features are described: 1) Multiple sets of closely spaced planar microstructures occur in quartz, plagioclase, and other rock-forming minerals. Those characteristic of shock consist of alternating platelets, with a range of reduced mean index of re fraction and birefringence; they con sist of platelets that have been partially or completely transformed to an amor phous phase. 2) Quartz and plagioclase are selec tively and completely transformed to silica glass and plagioclase glass in the solid state, whereas the associated mafic minerals remained crystalline. There is no reaction between adjacent minerals. 3) High-pressure polymorphs occur, such as coesite or stishovite. Coesite oc Curs exclusively within silica glass; it has not been observed as a reaction or breakdown product. 4) Nickel-iron spherules occur in the fused glass or impactites. 5) The occurrence of droplets of ilmenite, rutile, pseudobrookite, and baddaleyite in impactites indicates a temperature of formation exceeding 150°C. 6) Dense glass occurs, similar in composition to bulk rock, in which iron oxide, such as fine particles of mag netite, is completely dissolved. All these features are characteristic of a process involving the rapid rise and fall of extremiiely high pressures and temperatures. Minerals and mineral as semblages experiencing such high strain rates and sudden changes of pressures and temperatures react and change in dependently to the bulk chemical com position, under nonequilibrium condi tions. Many aspects of shock features re quire careful study. Kink bands in biotite and deformation lamellae in quartz occur in tectonically deformed rocks. These features should be studied with great care in order to determine whether reduction in mean index of refraction and total birefringence along the planar structures have resulted from vitrification or phase transition; their presence is additional evidence in favor of a shock mechanism. Vitreous phases or glasses formed by shock also have many unique prop erties; they have not been studied by such methods as thermoluminescence, electron spin resonance, low-angle x ray diffraction, or infrared spectroscopy. Shock-fused glass of high density needs to be studied in detail in carefully con trolled laboratory conditions. Experimental shock-wave studies of the equation-of-state of single minerals and mineral assemblages, under care fully controlled conditions, must pre cede estimates of peak pressures and peak and residual temperatures of shocked natural mineral assemblages. Detailed petrographic and mineralogic studies, however, have provided useful and definitive criteria for characteriza tion of impact events. Such data should be of paramount importance in the study of samples brought back from Moon.

Science↗
USGS2004

The Rock Elm meteorite impact structure, Wisconsin: Geology and shock-metamorphic effects in quartz

The Rock Elm structure in southwest Wisconsin is an anomalous circular area of highly deformed rocks, ∼6.5 km in diameter, located in a region of virtually horizontal undeformed sedimentary rocks. Shock-produced planar microstructures (PMs) have been identified in quartz grains in several lithologies associated with the structure: sandstones, quartzite pebbles, and breccia. Two distinct types of PMs are present: P1 features, which appear identical to planar fractures (PFs or cleavage), and P2 features, which are interpreted as possible incipient planar deformation features (PDFs). The latter are uniquely produced by the shock waves associated with meteorite impact events. Both types of PMs are oriented parallel to specific crystallographic planes in the quartz, most commonly to c(0001), xi{11 (2) over bar2}, and r/k{10 (1) over bar1} . The association of unusual, structurally deformed strata with distinct shock-produced microdeformation features in their quartz-bearing rocks establishes Rock Elm as a meteorite impact structure and supports the view that the presence of multiple parallel cleavages in quartz may be used independently as a criterion for meteorite impact. Preliminary paleontological studies indicate a minimum age of Middle Ordovician for the Rock Elm structure. A similar age estimate (450-400 Ma) is obtained independently by combining the results of studies of the general morphology of complex impact structures with estimated rates of sedimentation for the region. Such methods may be applicable to dating other old and deeply eroded impact structures formed in sedimentary target rocks.

Wisconsin↗
USGS1999

Air blasts generated by rockfall impacts: Analysis of the 1996 Happy Isles event in Yosemite National Park

The July 10, 1996, Happy Isles rockfall in Yosemite National Park, California, released 23,000 to 38,000 m 3 of granite in four separate events. The impacts of the first two events which involved a 550-m free fall, generated seismic waves and atmospheric pressure waves (air blasts). We focus on the dynamic behavior of the second air blast that downed over 1000 trees, destroyed a bridge, demolished a snack bar, and caused one fatality and several injuries. Calculated velocities for the air blast from a two-phase, finite difference model are compared to velocities estimated from tree damage. From tornadic studies of tree damage, the air blast is estimated to have traveled <108–120 m/s within 50 m from the impact and decreased to <10–20 m/s within 500 m from the impact. The numerical model simulates the two-dimensional propagation of an air blast through a dusty atmosphere with initial conditions defined by the impact velocity and pressure. The impact velocity (105–107 m/s) is estimated from the Colorado Rockfall Simulation Program that simulates rockfall trajectories. The impact pressure (0.5 MPa) is constrained by the kinetic energy of the impact (10 10 –10 12 J) estimated from the seismic energy generated by the impact. Results from the air blast simulations indicate that the second Happy Isles air blast (weak shock wave) traveled with an initial velocity above the local sound speed. The size and location of the first impact are thought to have injected <50 wt% dust into the atmosphere. This amount of dust lowered the local atmospheric sound speed to ∼220 m/s. The discrepancy between calculated velocity data and field estimated velocity data (∼220 m/s versus ∼110 m/s) is attributed to energy dissipated by the downing of trees and additional entrainment of debris into the atmosphere not included in the calculations.

Journal of Geophysical Research B: Solid Earth↗