10.1 Collision Tectonics
When Continents Collide
Continental collision is the terminal stage of the Wilson Cycle. It begins when subduction consumes all oceanic lithosphere separating two continental masses, forcing continent–continent convergence. Unlike oceanic lithosphere, continental crust is too buoyant ($\rho_c \approx 2,750$ kg/m³ vs. $\rho_m \approx 3,300$ kg/m³) to be subducted deeply into the mantle. The result is not consumption but crustal thickening, producing the highest mountain belts and most extensive plateaus on Earth.
The India–Eurasia collision, which began approximately 50 Ma when the Tethyan oceanic lithosphere was fully consumed, is the archetype. It has produced the Himalaya–Tibet orogenic system, the most dramatic expression of continental collision on the modern Earth: an orogen spanning over 2,500 km with crustal thicknesses approaching 70 km and elevations exceeding 8 km.
Suture Zones: Scars of Closed Oceans
A suture zone marks the boundary where two continental plates have been welded together following the closure of an intervening ocean basin. Suture zones preserve fragments of the destroyed oceanic lithosphere and the sediments that once lay upon it, providing the geological evidence for former ocean basins.
Ophiolites
Fragments of oceanic crust and upper mantle emplaced tectonically onto continental margins during collision. A complete ophiolite sequence (Penrose definition, 1972) includes, from bottom to top: ultramafic tectonite (harzburgite), layered gabbro, sheeted dike complex, and pillow basalts. Their presence in a mountain belt is a diagnostic indicator of a suture zone.
Mélanges
Chaotic mixtures of blocks of diverse lithology (basalt, chert, limestone, serpentinite, blueschist) set in a sheared sedimentary or serpentinite matrix. Mélanges form in the subduction channel and accretionary prism prior to collision, and are subsequently incorporated into the suture zone during the final stages of ocean closure.
High-Pressure Metamorphic Rocks
Blueschists and eclogites found along suture zones record the extreme P–T conditions of subduction (pressures of 1.5–3 GPa at temperatures of 300–600°C). Their exhumation to the surface, often via buoyancy-driven return flow in the subduction channel, preserves a record of the subduction process that preceded collision.
The Indus–Tsangpo Suture Zone (ITSZ) is the most prominent suture on Earth, extending over 2,000 km along the southern margin of the Tibetan Plateau. It marks where the Tethyan ocean closed as India collided with Eurasia. The ITSZ contains ophiolitic fragments (e.g., Xigaze ophiolite), mélanges, and remnants of the Tethyan sedimentary sequence, providing unequivocal evidence that a major ocean basin once separated India from Asia.
Crustal Thickening & Shortening
The primary mechanical response to continental collision is crustal thickening through thrust stacking. Normal continental crust has a thickness of approximately 35 km. In the Himalaya–Tibet system, the crust has been thickened to approximately 70 km — a doubling achieved through imbricate thrust sheets that stack slices of crust atop one another.
The total amount of crustal shortening can be estimated from balanced cross-sections, which reconstruct the pre-deformation geometry by unfolding and unfaulting the observed structures. The shortening $\Delta L$ is the difference between the restored length and the present-day length of the section:
Crustal shortening from balanced cross-section:
\[ \Delta L = L_0 - L_f \]
where $L_0$ is the original (restored) length and $L_f$ is the final (deformed) length of the section. Shortening strain $e = \Delta L / L_0$.
For the India–Eurasia collision, the total convergence since ~50 Ma is approximately 2,500 km (based on plate kinematic reconstructions from seafloor magnetic anomalies). This enormous convergence has been accommodated through a combination of mechanisms, often summarized as the strain budget:
Strain budget for continental collision:
\[ \Delta L_{\text{total}} = \Delta L_{\text{thickening}} + \Delta L_{\text{extrusion}} + \Delta L_{\text{erosion}} + \Delta L_{\text{subduction}} \]
Total convergence is partitioned among crustal thickening (thrust stacking), lateral extrusion, erosional removal, and possible subduction of continental lithosphere.
~2,500 km
Total Convergence Since 50 Ma
~70 km
Tibetan Crustal Thickness
~35 km
Normal Continental Crust
Lateral Extrusion & Escape Tectonics
Not all convergence is accommodated by crustal thickening in the collision zone itself. A significant fraction is absorbed by lateral extrusion — the tectonic escape of large crustal blocks sideways along major strike-slip fault systems. This concept, pioneered by Tapponnier & Molnar (1976), was inspired by analog experiments using a rigid indenter (India) pushed into a deformable medium (Asia).
In the India–Eurasia system, the Indochina block has been extruded southeastward along two major left-lateral fault systems: the Ailao Shan–Red River fault and the Sagaing fault. GPS geodesy shows that present-day eastward motion of eastern Tibet and Yunnan occurs at rates of 10–15 mm/yr. Paleomagnetic data and structural reconstructions suggest that Indochina may have rotated clockwise by 15–25° and translated ~500–1,000 km to the southeast since collision began.
Other major strike-slip faults associated with lateral extrusion include the Altyn Tagh fault (left-lateral, ~10 mm/yr) bounding the northern margin of Tibet, and the Karakorum fault(right-lateral, ~10 mm/yr) on the western margin. Together, these faults partition the broad India–Eurasia deformation zone into discrete blocks that move laterally, distributing the effects of collision over an area of approximately 5 × 10&sup6; km².
Gravitational Potential Energy & Orogenic Collapse
As the crust thickens and the orogen grows, the elevated mass of the mountain belt represents an enormous reservoir of gravitational potential energy (GPE). When the GPE of the orogen significantly exceeds that of the surrounding lowlands, the excess drives extensional deformation and lateral spreading — a process known as orogenic collapse or gravitational collapse.
Gravitational potential energy per unit area:
\[ \text{GPE} = \int_{-h_c}^{h_{\text{elev}}} \rho(z)\, g\, z \; dz \]
where the integration extends from the base of the crust ($-h_c$) to the surface elevation ($h_{\text{elev}}$), $\rho(z)$ is the density profile, and $g$ is gravitational acceleration. The excess GPE ($\Delta$GPE) relative to a reference column drives horizontal deviatoric stresses.
For the Tibetan Plateau (elevation ~5 km, crustal thickness ~70 km), the excess GPE relative to normal continental lithosphere generates horizontal deviatoric stresses of approximately 50–100 MPa. These stresses are sufficient to drive the normal faulting and east–west extension observed across Tibet today, as well as the lateral extrusion of eastern Tibet.
The balance between continued tectonic convergence (building the orogen) and gravitational collapse (tearing it down) determines whether the orogen continues to grow in elevation, reaches a steady state, or collapses. The Tibetan Plateau appears to be in a quasi-steady state where the convergence-driven compressive stress is approximately balanced by the excess GPE.
Thin Viscous Sheet Model
The thin viscous sheet model, developed by England & McKenzie (1982), provides a continuum mechanics framework for understanding distributed deformation in continental collision zones. The model treats the lithosphere as a thin layer of viscous fluid deforming under the combined effects of boundary forces (indentation) and internal body forces (gravity acting on topography).
The key assumption is that the lithosphere is thin compared to the horizontal scale of deformation, so vertical variations in stress can be averaged out. The vertically averaged stress balance is:
Vertically averaged stress equilibrium:
\[ \frac{\partial \bar{\sigma}_{xx}}{\partial x} + \frac{\partial \bar{\tau}_{xy}}{\partial y} = \frac{\partial}{ \partial x}\left(\text{GPE}\right) \]
where $\bar{\sigma}_{xx}$ and $\bar{\tau}_{xy}$ are vertically averaged stress components, and the right-hand side represents the gradient in gravitational potential energy that drives flow from thick to thin crust.
The behavior of the thin viscous sheet is characterized by the Argand number, which measures the ratio of gravitational (buoyancy) stresses to viscous stresses arising from the imposed convergence:
Argand number:
\[ \text{Ar} = \frac{\Delta\rho \, g \, L}{\eta \, v} \]
where $\Delta\rho$ is the density contrast between crust and mantle,$L$ is a characteristic crustal thickness, $\eta$ is the effective viscosity, and $v$ is the convergence velocity.
Low Argand Number (Ar << 1)
Gravitational stresses are small compared to viscous stresses. The lithosphere is strong and deformation is concentrated along narrow fault zones. This regime characterizes rigid-block behavior with discrete, localized faults separating relatively undeformed blocks — similar to the indenter-escape model of Tapponnier & Molnar.
High Argand Number (Ar >> 1)
Gravitational stresses dominate. The crust cannot sustain significant topographic gradients and flows laterally to maintain a nearly uniform thickness. This produces broad, flat plateaus with diffuse deformation — consistent with the appearance of the Tibetan Plateau, where elevations are remarkably uniform at ~5 km despite active deformation throughout.
Key Numbers
~50 Ma
India–Eurasia Collision Onset
~5 cm/yr
Current Convergence Rate
2,500 km
Total Shortening
5 × 10&sup6; km²
Deformation Zone Area