About Earthquake Fault — Stick-Slip Simulator
This simulation models earthquake fault mechanics using the Burridge-Knopoff spring-slider framework, where a grid of fault patches accumulates tectonic stress at a constant loading rate until each patch exceeds its static friction threshold, then suddenly slips and drops to kinetic friction — the classic stick-slip cycle. The fault plane visualises stress levels by colour (blue = low, red = near-failure), and ruptures cascade through neighbouring patches via stress transfer, producing a realistic distribution of earthquake sizes. Users can observe how seismic moment M0 is computed from slip area and mean displacement, how moment magnitude Mw is derived from M0, and how the resulting event catalogue follows the Gutenberg-Richter frequency-magnitude relationship.
Stick-slip friction on geological faults is the physical mechanism behind virtually every tectonic earthquake on Earth, from microseisms below magnitude 0 to great megathrust events above magnitude 9 such as the 2011 Tohoku earthquake. The Gutenberg-Richter law — an empirical power-law describing how many earthquakes of each magnitude occur — has been confirmed across tectonic settings worldwide and underpins modern probabilistic seismic hazard analysis.
Frequently Asked Questions
What is the stick-slip mechanism and why does it cause earthquakes?
Stick-slip is a frictional instability in which two surfaces lock together (stick) while stress builds, then suddenly slide past each other (slip) when the applied stress exceeds static friction. On a geological fault, tectonic plate motion continuously loads the fault with elastic stress over years to centuries. When stress surpasses the fault's static friction strength, rupture initiates and propagates along the fault plane, releasing stored elastic energy as seismic waves — what we feel as an earthquake.
How do I use the simulation controls to explore fault behaviour?
Use the Loading rate (vL) slider to control how fast tectonic stress accumulates — higher rates produce more frequent earthquakes. The Static friction (mu_s) slider sets the stress threshold required to start rupture, while Kinetic friction (mu_k) controls the stress level after slip; a larger gap between mu_s and mu_k produces larger, more energetic events. Increase the Fault area to scale up seismic moment M0 and magnitude Mw. Watch the colour map for stress hotspots, the time-series chart for stress cycles and magnitude spikes, and the Gutenberg-Richter plot to see the b-value emerge as events accumulate.
What does the Gutenberg-Richter b-value mean, and what value does this simulation produce?
The Gutenberg-Richter law states that log N = a - bM, where N is the number of earthquakes with magnitude greater than or equal to M, and b is the slope describing the relative abundance of small versus large events. A b-value near 1.0 is typical of natural tectonic seismicity globally. In this simulation the b-value is estimated by linear regression on the log-frequency versus magnitude plot; adjusting friction contrast (mu_s minus mu_k) shifts the b-value, with greater contrast tending to produce more large events and a lower b.
How is seismic moment M0 calculated, and how does it relate to magnitude Mw?
Seismic moment is defined as M0 = mu * A * D, where mu is the shear modulus of the crust (approximately 30 GPa for the upper crust), A is the rupture area, and D is the mean co-seismic slip. This simulation uses mu = 3 x 10^10 Pa, computes D from the summed slip across ruptured patches, and applies the Hanks-Kanamori formula Mw = (log10(M0) - 9.1) / 1.5 to convert to moment magnitude. Unlike the older Richter local magnitude scale, Mw does not saturate for large earthquakes and is the standard scale used by seismological agencies worldwide.
What real-world faults exhibit the clearest stick-slip behaviour?
The San Andreas Fault in California is the textbook example: the locked section between Parkfield and San Bernardino accumulates elastic strain at roughly 35 mm/yr but ruptures infrequently in major events such as the 1906 San Francisco earthquake (Mw 7.9). By contrast, the creeping section near Hollister slips aseismically at near-plate-rate with only small earthquakes — a behaviour controlled by velocity-strengthening friction minerals (serpentinite and talc) rather than the velocity-weakening friction modelled here. Subduction megathrusts such as the Cascadia, Nankai, and Chilean zones show the largest stick-slip events of all.
Is it a misconception that earthquakes always occur on visible surface faults?
Yes. Many damaging earthquakes occur on blind faults — faults that do not break the surface and may not even be mapped before the event. The 1994 Northridge earthquake (Mw 6.7) in California and the 2010 Canterbury earthquake (Mw 7.1) in New Zealand both ruptured previously unknown blind thrust faults. Additionally, intraplate earthquakes far from tectonic plate boundaries, such as the 1811-1812 New Madrid sequence in the central United States, demonstrate that significant seismicity can occur well away from obvious fault zones.
Who developed the spring-slider model used here, and when?
Robert Burridge and Leon Knopoff introduced the spring-slider (or slider-block) model of earthquake faults in a landmark 1967 paper in the Bulletin of the Seismological Society of America. Their model replaced a continuous fault with a one-dimensional chain of blocks connected by springs to a moving driver plate, capturing the essential stick-slip dynamics while remaining mathematically tractable. The model was later extended to two dimensions and used by Per Bak, Chao Tang, and Kurt Wiesenfeld in 1987-1989 to demonstrate self-organised criticality — the idea that the fault system naturally evolves to a critical state producing power-law earthquake statistics without fine-tuning of parameters.
What other phenomena and simulations are closely related to earthquake fault dynamics?
Earthquake faults share physics with several other systems. Avalanche dynamics (snow, sand, and granular piles) exhibit the same self-organised critical statistics as seismicity. Acoustic emission in materials science uses the same Gutenberg-Richter power law to characterise micro-crack populations ahead of fracture. Tidal triggering links fault mechanics to fluid dynamics and celestial mechanics. Related simulations worth exploring include wave interference (seismic waves propagating away from the rupture), climate tipping points (another threshold-crossing system with positive feedbacks), and the double pendulum (chaotic sensitivity to initial conditions, analogous to earthquake unpredictability).
How is earthquake science used in seismic hazard engineering?
Probabilistic Seismic Hazard Analysis (PSHA) quantifies the likelihood that ground shaking will exceed a specified level at a site within a given time window, typically 50 years for building codes. It integrates Gutenberg-Richter recurrence rates, rupture area scaling laws, ground-motion prediction equations, and site amplification factors. The outputs feed directly into building code design spectra (such as Eurocode 8 and ASCE 7), bridge and dam safety standards, nuclear facility licensing, and earthquake insurance pricing. Real-time seismic monitoring networks also use automated magnitude and slip estimates — quantities this simulation computes — to trigger tsunami warnings within minutes of a large offshore earthquake.
What are the current frontiers in earthquake fault research?
Several open questions drive active research: slow slip events and episodic tremor on subduction zones (observed on Cascadia and Nankai) behave like very slow earthquakes lasting days to weeks, suggesting a continuum between seismic and aseismic slip. Earthquake early warning systems (such as ShakeAlert in the western United States and Japan's nationwide system) aim to provide seconds of warning before strong shaking arrives, exploiting the fact that P-waves travel faster than the destructive S-waves. Machine learning is being applied to detect subtle seismic signals that precede some foreshock sequences, though reliable short-term earthquake prediction remains an unsolved problem. Additionally, induced seismicity from wastewater injection (linked to oil and gas operations) has sharply increased seismic rates in previously quiet regions such as Oklahoma, raising questions about anthropogenic fault reactivation thresholds.