Planetary Geology · Impact Physics
July 2026 · 14 min read · Shock Waves · Scaling Laws · Crater Morphology · Last updated: 3 July 2026

Impact Craters: The Physics of Hypervelocity Collisions

Written by MySimulator Team · Reviewed by MySimulator Editorial Review

Every rocky surface in the Solar System keeps a scoreboard of violence written in circles. A crater is not simply a hole punched by a falling rock — it is the frozen record of a shock wave that swept through solid rock at kilometres per second, vaporising, melting, fracturing, and hurling material across a planet in a matter of seconds. This article develops the physics of that process stage by stage: the contact and compression that launches the shock, the excavation flow that carves the bowl, the collapse that turns a simple pit into a mountain-ringed basin, the scaling laws that let geologists work backwards from crater size to impactor energy, the mineralogical fingerprints shock leaves behind, and how crater counting becomes a stopwatch for reading the age of a planet's surface.

1. Contact and Compression

Impact cratering unfolds in three overlapping stages: contact and compression, excavation, and modification. The first stage is astonishingly brief — for a kilometre-scale asteroid it lasts a fraction of a second — yet it sets the energy budget for everything that follows.

The moment an impactor touches a planetary surface, its leading face decelerates almost instantly while the rest of the body is still travelling at full speed. This mismatch launches two shock waves: one propagating downward and outward into the target rock, and a second, the reflected shock, propagating backward up into the impactor itself. Because natural impact velocities (typically 11–72 km/s for asteroids and comets hitting Earth) vastly exceed the speed of sound in rock (~5–6 km/s), the collision is hypervelocity: the impactor cannot simply push material out of the way, it must shock-compress it.

Impactor kinetic energy delivered to the target: E = (1/2) · m · v² For a 1 km diameter, 2600 kg/m³ stony asteroid at v = 20 km/s: m ≈ 1.36 × 10¹² kg E ≈ 2.7 × 10²⁰ J (roughly 6.5 million megatons of TNT)

Contact and compression ends when the reflected shock reaches the back of the impactor and releases as a rarefaction (relief) wave, which unloads the compressed impactor material and typically vaporises or melts it outright. By this point a hemispherical shock front is already racing into the target rock, and the excavation stage begins.

2. Shock Waves and Peak Pressure

A shock wave is not an ordinary sound wave — it is a discontinuity across which pressure, density, and temperature jump almost instantaneously, governed by the Rankine-Hugoniot jump conditions that express conservation of mass, momentum, and energy across the shock front.

Rankine-Hugoniot relations across a shock front: Mass: ρ₀·Us = ρ₁·(Us − Up) Momentum: P₁ − P₀ = ρ₀·Us·Up Energy: E₁ − E₀ = (1/2)·(P₁ + P₀)·(1/ρ₀ − 1/ρ₁) where ρ₀, P₀, E₀ = density, pressure, energy ahead of the shock (undisturbed) ρ₁, P₁, E₁ = density, pressure, energy behind the shock (shocked state) Us = shock front velocity, Up = particle velocity behind the shock

Peak shock pressures near the impact point can reach hundreds of gigapascals — far beyond the strength of any rock, which is why, at the moment of impact, both the target and the impactor briefly behave like fluids. As the shock front expands away from the point of impact its energy is spread over an ever-larger hemispherical surface, so peak pressure falls off steeply with distance:

Approximate peak pressure decay with distance r from the impact point: P(r) ≈ P₀ · (r₀ / r)ⁿ with n ≈ 2 to 3 (material dependent) Rough shock-pressure regimes and their effects on silicate rock: > 60 GPa → whole-rock melting, vaporisation near the point of impact 10–60 GPa → high-pressure mineral phases (e.g. coesite, stishovite), diaplectic glass 2–10 GPa → shatter cones, planar deformation features in quartz < 2 GPa → simple brittle fracturing, familiar "impact breccia"

This steep falloff explains why a crater has concentric zones: totally melted rock at the centre, intensely shocked and mixed breccia around it, and merely fractured bedrock at the rim — a shock-pressure gradient frozen into the geology.

3. The Excavation Stage

As the shock wave expands, it leaves in its wake a rarefaction wave that accelerates rock away from the impact point, setting up the well-known excavation flow — a symmetric pattern first described in the Maxwell Z-model, in which material moves outward and upward along curved streamlines centred on a point below the surface.

Rock closest to the impact point is thrown out ballistically as ejecta, landing beyond the crater rim and building the characteristic ejecta blanket, with the largest, latest-launched fragments forming rays that can stretch for hundreds of kilometres (as seen around lunar craters like Tycho). Rock farther from the centre is displaced outward and downward without ever leaving the ground, forming the crater walls — this deeper material is not excavated so much as folded and overturned, producing an overturned flap at the rim where the original stratigraphy is flipped upside down.

Transient crater diameter (before any wall collapse), Maxwell Z-model streamline: z = Z / (r^(Z−1)) [schematic streamline exponent, Z ≈ 2.7–3.5] The transient crater's depth-to-diameter ratio is roughly constant: d_t / D_t ≈ 1/3 (for both simple and complex craters, before modification)

Excavation continues until the growing crater has converted essentially all of the impactor's kinetic energy into kinetic energy of ejected and displaced rock, heat, and the energy locked into permanent deformation. At this point the crater reaches its maximum transient size — a bowl about three times deeper (relative to its width) than the crater that will ultimately be preserved, because what happens next is collapse.

4. Modification: Collapse and Rebound

The transient crater is gravitationally unstable: its walls are far steeper than the angle of repose of shattered rock, and gravity immediately begins to pull it back toward equilibrium. What happens during this modification stage depends strongly on the crater's size.

In small craters, the walls simply slump inward under gravity, widening the crater slightly and burying the floor under a layer of slumped breccia — the result is a simple crater, a clean bowl shape.

In larger craters, the floor itself is weak enough — briefly acting almost like a fluid due to acoustic fluidisation and the sheer scale of the shattered, shock-heated rock — that it rebounds elastically-plastically upward from the point of maximum compression, like the recoil of a liquid after a stone is dropped in it. This central rebound overshoots and produces a central peak, and the outer walls collapse inward along listric normal faults to form concentric terraces. In the very largest impacts the central uplift itself becomes unstable and collapses outward into a ring, producing a peak ring or, in the biggest basins, multiple concentric rings.

Approximate final crater diameter after collapse (for the simple-to-complex transition): D_final ≈ 1.3 to 1.7 × D_transient (widening from wall slumping/collapse) Central peak height scales with final crater diameter D (km), gravity g, and target/impactor properties captured in an empirical constant k: h_peak ≈ k · D^n (n ≈ 1, empirically fit per planetary body)

Modification typically finishes within minutes for even the largest impacts — geologically instantaneous compared to the millions of years of erosion that follow and slowly erase the crater's fine structure.

5. Crater Scaling Laws

Turning an impactor's size, speed, and angle into a final crater diameter is the central engineering problem of impact cratering, and it is solved with Pi-group scaling (dimensional analysis using the Buckingham Pi theorem), calibrated against explosion craters, laboratory gun experiments, and numerical hydrocode simulations.

Pi-group scaling (gravity-dominated regime, simplified form): π_D = D · (g / v²)^μ dimensionless crater-size group π_2 = 1.61 · g · a / v² dimensionless gravity-scaled size where D = crater diameter g = surface gravity v = impact velocity a = impactor radius μ ≈ 0.41–0.55 (coupling exponent, depends on target: sand, rock, ice) A commonly used gravity-regime scaling relation (Holsapple & Housen): D ∝ (m/ρ_t)^(1/3) · (g·a/v²)^(−μ/(2+μ)) · (ρ_i/ρ_t)^ν where m = impactor mass, ρ_i, ρ_t = impactor and target density, ν ≈ 0.4

The key insight of scaling laws is that crater growth is not simply proportional to impactor energy. For most impacts above a few tens of metres, crater growth is gravity-dominated — the final size is set by the competition between the outward-moving excavation flow and the inward pull of gravity, so the same energy delivered on a low-gravity body (the Moon, an asteroid) produces a much larger crater than on Earth. Below roughly ten metres, cratering is instead strength-dominated, controlled by the target rock's tensile and shear strength rather than its weight.

A useful rule of thumb from these relations: a stony impactor traveling at typical asteroid encounter speeds excavates a crater roughly 10–20 times its own diameter on a rocky planet with Earth-like gravity — a 1 km asteroid can carve a crater 15–20 km across.

6. Simple vs Complex Craters

Whether a crater ends up simple or complex depends almost entirely on its final diameter relative to a threshold set by the target's gravity and material strength — the simple-to-complex transition diameter.

Approximate simple-to-complex transition diameters (rocky targets): Earth (g = 9.8 m/s²): D_c ≈ 2–4 km (lower for sedimentary, higher for crystalline rock) Moon (g = 1.6 m/s²): D_c ≈ 15–20 km Mars (g = 3.7 m/s²): D_c ≈ 5–10 km The transition scales approximately as: D_c ∝ 1/g (weaker gravity → larger craters can stay "simple")

7. Shock Metamorphism

Because impacts are the only common geological process that generates pressures above about 2 GPa in an instant (as opposed to the slow burial pressures of ordinary metamorphism), the minerals left behind are diagnostic — they are how geologists distinguish a genuine impact structure from a volcanic caldera or collapsed sinkhole that merely looks similar from orbit.

These signatures, combined with the presence of a geochemical anomaly such as elevated iridium (rare in Earth's crust but common in asteroids), form the standard evidence trail used to confirm and date impact structures — famously, the globally distributed iridium-rich clay layer at the Cretaceous-Paleogene boundary that led Walter and Luis Alvarez to propose the Chicxulub impact in 1980.

8. Crater Chronology: Dating a Surface

On worlds without plate tectonics or running water to erase old terrain — the Moon, Mars, Mercury — impact craters accumulate steadily over billions of years, and their density becomes a natural stopwatch. This is crater chronology: the older a surface, the more craters per unit area it has collected, following an approximately power-law size-frequency distribution punctuated by a period of intense early bombardment.

Crater production function (cumulative number N of craters ≥ diameter D per unit area): N(≥D) ∝ D^(−b) with b ≈ 2 to 3 for most planetary surfaces Relative age from crater density: older surface → more time to accumulate craters → higher N(≥D) at fixed D

Absolute ages are pinned down by calibrating this relative crater-counting method against radiometric dating of Apollo and Luna sample-return rocks from the Moon — the only body where scientists can directly compare a surface's crater density to the measured age of the rock beneath it. This calibration is then extrapolated (with real uncertainty) to crater counts on Mars, Mercury, and icy moons, where no samples yet exist. The technique revealed the Late Heavy Bombardment, an apparent spike in impact rates across the inner Solar System roughly 3.9–4.1 billion years ago, still debated as either a real cataclysm or partly an artefact of how the earliest lunar surface has been sampled.

9. Famous Craters

A handful of impact structures illustrate the full range of the physics above:

☄️
Asteroid Deflection Simulator
Model impactor mass, velocity, and trajectory to see how much deflection it takes to miss a planet
🌋
Volcanic Eruption Simulator
Compare an impact's instantaneous shock energy against the sustained power of an eruption
🌕
Moon Phases Simulator
See the cratered lunar surface that this article's chronology techniques are used to date