Structural Geology

2.1 The Stress Tensor

The state of stress at a point in a deforming rock is described by the Cauchy stress tensor, a second-rank tensor with nine components. In a Cartesian coordinate system, the stress tensor $\sigma_{ij}$ relates the traction vector $\mathbf{T}$ acting on a surface with outward unit normal $\hat{n}$:

The full stress tensor in 3D is:

By conservation of angular momentum, the stress tensor is symmetric: $\sigma_{ij} = \sigma_{ji}$, reducing the independent components from nine to six. The diagonal elements are normal stresses and the off-diagonal elements are shear stresses.

Derivation 1: Principal Stresses from Eigenvalue Problem

There exist special orientations where the traction vector is purely normal to the surface (no shear component). These are the principal stress directions. The principal stresses are found by solving the eigenvalue problem:

For a non-trivial solution, the determinant must vanish:

Expanding this determinant yields the characteristic equation, a cubic polynomial in $\sigma$:

where the three stress invariants are:

The three roots $\sigma_1 \geq \sigma_2 \geq \sigma_3$ are the principal stresses. By convention in geology, compressive stress is positive. The maximum compressive stress $\sigma_1$ controls fault orientation and type.

Historical Context

Augustin-Louis Cauchy (1789-1857) formalized the stress tensor concept in the 1820s, building on the earlier work of Euler and Navier on continuum mechanics. The application to geological structures was pioneered by M. King Hubbert (1951) and John Ramsay (1967), who systematically related stress fields to fault and fold geometries observed in the field. Anderson's theory of faulting (1905, 1951) connected principal stress orientations to the three fundamental fault types, providing a mechanical framework for structural geology.

2.2 Mohr Circle Analysis

The Mohr circle provides a graphical representation of the stress state on planes of all orientations through a point. For the 2D case, consider principal stresses $\sigma_1$ and $\sigma_3$.

Derivation 2: Normal and Shear Stress on an Arbitrary Plane

Consider a plane whose normal makes angle $\theta$ with the $\sigma_1$ direction. The traction vector has components:

The normal stress on this plane is $\sigma_n = T_1 \cos\theta + T_3 \sin\theta$:

The shear stress is $\tau = -T_1 \sin\theta + T_3 \cos\theta$:

These two equations define a circle in $(\sigma_n, \tau)$ space centered at$\left(\frac{\sigma_1 + \sigma_3}{2}, 0\right)$ with radius $R = \frac{\sigma_1 - \sigma_3}{2}$. This is the Mohr circle. Maximum shear stress occurs at $\theta = 45°$ and equals $R$.

The Coulomb Failure Criterion

Rock failure occurs when the stress state reaches the failure envelope. The Mohr-Coulomb criterion states:

where $c$ is cohesion, $\mu = \tan\phi$ is the coefficient of internal friction, and $\phi$ is the angle of internal friction. For most rocks, $\phi \approx 30°$, giving $\mu \approx 0.6$ (Byerlee's law). The failure plane makes an angle $\theta = 45° - \phi/2$ with $\sigma_1$, which is approximately $30°$.

Otto Mohr (1835-1918) introduced the graphical circle construction in 1882. Charles-Augustin de Coulomb had earlier (1773) proposed the linear failure criterion for soils. The combined Mohr-Coulomb criterion remains the most widely used failure criterion in structural geology and geotechnical engineering.

2.3 Strain Analysis

Derivation 3: The Strain Tensor from Displacement Gradients

When a body deforms, a point at position $\mathbf{x}$ displaces to $\mathbf{x} + \mathbf{u}$, where $\mathbf{u}$ is the displacement vector. The displacement gradient tensor is:

This tensor can be decomposed into symmetric and antisymmetric parts:

The symmetric part is the infinitesimal strain tensor:

The antisymmetric part represents rigid-body rotation:

For finite strains common in structural geology, we use the Green-Lagrange strain tensor:

The stretch ratio $\lambda = l/l_0$ gives the ratio of deformed to undeformed length. Natural (logarithmic) strain is $\varepsilon_{\text{nat}} = \ln\lambda$, which is additive for sequential deformations, making it preferred for large strains.

The Flinn Diagram

Derek Flinn (1962) introduced a diagram to classify strain ellipsoids using the ratios of principal stretches. The parameter $K$ is defined as:

where $a = \lambda_1/\lambda_2$ and $b = \lambda_2/\lambda_3$. When $K = 1$, the strain is plane strain; $K > 1$ is prolate (constrictional); $K < 1$ is oblate (flattening).

2.4 Anderson's Theory of Faulting

Derivation 4: Fault Type from Principal Stress Orientation

E.M. Anderson (1905, 1951) recognized that Earth's surface is a free surface with zero shear traction, so one principal stress must be vertical (equal to the overburden pressure $\sigma_v = \rho g z$). The fault type depends on which principal stress is vertical:

The fault plane dips at $\theta = 45° - \phi/2$ from $\sigma_1$. For $\phi = 30°$:

• Normal faults dip at ~60° from horizontal
• Reverse faults dip at ~30° from horizontal (thrust faults dip ← 30°)
• Strike-slip faults are vertical, striking at ~30° from $\sigma_1$

Derivation 5: Effective Stress and Pore Pressure

Karl Terzaghi (1923) introduced the concept of effective stress. In rocks with fluid-filled pores at pressure $P_f$, the effective stress is:

The Mohr-Coulomb criterion becomes:

High pore fluid pressure reduces the effective normal stress, shifting the Mohr circle to the left and bringing it closer to the failure envelope. This explains how thrust faults can slip at low apparent shear stresses (the Hubbert-Rubey paradox, 1959). The pore pressure ratio is defined as:

Hydrostatic pore pressure gives $\lambda_p \approx 0.4$. Lithostatic (overpressured) conditions approach $\lambda_p = 1$, where the effective stress vanishes and even small differential stresses can cause failure. This is critical for understanding fault reactivation, induced seismicity from wastewater injection, and the mechanics of accretionary wedges at subduction zones.

2.5 Fold Geometry and Mechanics

Folds form when layered rocks undergo lateral compression. The geometry of folds is described by wavelength $\lambda$, amplitude $A$, interlimb angle, and axial plane orientation. Ramsay (1967) classified folds using dip isogons into Classes 1A, 1B (parallel), 1C, 2 (similar), and 3.

Biot's Folding Theory

Maurice Biot (1961) analyzed the buckling instability of a competent layer of viscosity $\eta_1$ and thickness $h$ embedded in a weaker medium of viscosity $\eta_2$. The dominant wavelength is:

This predicts that strong layers in weak matrix produce long-wavelength folds, while layers with small viscosity contrast produce short-wavelength or no folds. The wavelength-to-thickness ratio $\lambda_d/h$ is a key observable in the field and directly constrains the viscosity ratio.

Applications in Structural Analysis

Structural geology has direct applications in resource exploration and hazard assessment. Fold-and-thrust belts host major petroleum reserves (e.g., the Zagros, Appalachian fold belt). Understanding fault geometry is essential for seismic hazard assessment. Modern techniques include balanced cross-sections, kinematic modeling, and integration with seismic reflection data.

2.6 Rock Rheology and Deformation Mechanisms

The brittle-ductile transition marks a fundamental change in deformation behavior with depth. In the upper crust, rocks deform by fracture (brittle behavior governed by the Mohr-Coulomb criterion). Below ~10-15 km (depending on geothermal gradient and rock type), rocks deform by crystal-plastic mechanisms (ductile flow).

Strength Envelopes

The lithospheric strength profile (yield stress envelope or "Christmas tree" diagram) combines Byerlee's law for brittle failure with power-law creep for ductile flow. The brittle strength increases linearly with depth:

where $\beta$ depends on faulting regime (extension vs compression) and $\lambda_p$is the pore pressure ratio. The ductile strength decreases exponentially with temperature through the power-law creep equation:

where $\dot{\varepsilon}$ is strain rate, $n \approx 3$-4 for olivine and quartz, and $E^*$ is the activation energy (500 kJ/mol for olivine, 135 kJ/mol for wet quartz). The intersection of the brittle and ductile curves defines the brittle-ductile transition depth and the peak lithospheric strength. Continental lithosphere is weakest at the Moho (due to the quartz-olivine transition), while oceanic lithosphere is strongest in the upper mantle.

Shear Zones

At depth, deformation localizes into ductile shear zones rather than discrete faults. These zones show progressive grain size reduction (mylonites) and development of crystallographic preferred orientation (CPO). The shear zone thickness depends on the strain rate, temperature, and grain size. Pseudotachylites (friction-melted rock) record rare seismic slip events at depths normally considered aseismic, bridging the gap between structural geology and seismology.

2.7 Methods: Stereographic Projection and Cross-Sections

Stereographic projection maps 3D orientation data (planes and lines) onto a 2D circle, enabling rapid analysis of structural relationships. A plane with strike $\alpha$ and dip $\delta$ projects as a great circle; a line with trend and plunge projects as a point.

The equal-area (Schmidt net) projection preserves areal relationships and is used for statistical analysis of fabric data. Pi ($\pi$) diagrams plot poles to bedding to determine fold axis orientations. Beta ($\beta$) diagrams plot intersections of great circles. Contoured stereonets reveal preferred orientations in large datasets.

Balanced cross-sections enforce conservation of area (or line length) to ensure geometric validity. The technique, developed by Dahlstrom (1969), requires that deformed and restored sections have the same area, providing a test of proposed structural interpretations. This is fundamental to petroleum exploration in fold-and-thrust belts and to understanding mountain building processes.

Strain Measurement Techniques

Strain in rocks can be measured using deformed objects of known initial shape. Common strain markers include:

• Ooids and reduction spots: Initially spherical objects deformed to ellipsoids. The R$_f$/$\phi$ method (Ramsay, 1967) uses the final axial ratio and orientation to determine strain.
• Fossils: Brachiopods, trilobites, and other fossils of known undeformed geometry. The Wellman (1962) method uses three non-parallel lines to construct the strain ellipse.
• Pressure shadows around rigid objects: The geometry of fibrous overgrowths on pyrite, magnetite, or garnet records the incremental strain history.
• Deformed veins and dikes: Boudinaged layers indicate extension; folded veins indicate shortening. The boudin aspect ratio constrains the viscosity contrast.

The strain ratio $R_s = \lambda_1/\lambda_3$ (the ratio of maximum to minimum stretch) ranges from 1 (undeformed) to →10 in highly strained rocks. In mylonite zones, strain ratios can exceed 50, with principal axes oriented at small angles to the shear zone boundary. The relationship between finite strain and shear zone width constrains the displacement across the zone.

Computational Lab: Structural Geology