Earthquakes and Engineering Seismology

Introduction

Earthquake engineering begins with an understanding of how earthquakes originate, how seismic energy travels through the Earth, and how shaking is measured and interpreted. Structural design is not based on the abstract existence of an earthquake alone; it depends on the characteristics of a particular seismic source, the path between source and site, local ground conditions, and the vulnerability of the built environment.

Engineering seismology applies geophysical knowledge to engineering decisions. It provides the basis for identifying seismic hazards, selecting design ground motions, interpreting instrumental records, estimating recurrence, and developing appropriate structural and non-structural protection. Because earthquake occurrence is uncertain, earthquake engineering is fundamentally a field of decision-making under uncertainty rather than a discipline of exact prediction.

Learning Objectives

After studying this chapter, the learner should be able to:

1. Explain plate tectonics and the principal geological causes of earthquakes.

2. Distinguish faults, hypocentre, epicentre, rupture dimensions, and seismic source mechanisms.

3. Describe body waves and surface waves and relate them to engineering damage.

4. Distinguish magnitude, intensity, energy, and ground-motion measures.

5. Explain recurrence relationships and the limitations of earthquake prediction.

6. Interpret basic seismological information for preliminary engineering assessment.

7. Identify the implications of uncertainty, site effects, and data quality for engineering decisions.

1. The Earth as a Seismic System

1.1 Internal structure and tectonic plates

The Earth consists broadly of the crust, mantle, outer core, and inner core. The mechanically strong outer shell, called the lithosphere, is divided into tectonic plates. These plates move relative to one another over the weaker, ductile asthenosphere. Plate motion is driven by processes including mantle convection, slab pull, ridge push, and gravitational adjustment.

Most damaging earthquakes occur in the shallow crust, generally within the upper tens of kilometres, although earthquakes can occur much deeper in subducting slabs. The major tectonic environments are:

• Divergent boundaries, where plates move apart and new crust forms. Earthquakes are commonly shallow and are associated with normal faulting.
• Convergent boundaries, where plates move together. Subduction zones can generate very large earthquakes and tsunamis; continental collision zones produce intense crustal deformation.
• Transform boundaries, where plates slide laterally past one another. Strike-slip earthquakes are characteristic.
• Intraplate regions, located within plates, where earthquakes may occur on ancient or reactivated weaknesses. Their recurrence is often more difficult to estimate because instrumental and historical records are sparse.

Plate boundaries explain the global distribution of seismicity, but they do not identify every active fault. Local geological investigations remain essential.

1.2 Stress, strain, and elastic rebound

Stress is the force per unit area acting within a material, while strain is deformation expressed relative to the original geometry. Tectonic plate motion gradually loads rocks surrounding a fault. At low stress levels, rocks may deform approximately elastically and recover much of their deformation when the load is removed. If frictional resistance and rock strength are exceeded, sudden slip occurs and stored elastic strain energy is released as seismic energy, heat, and permanent deformation.

The elastic rebound theory, developed from observations of the 1906 San Francisco earthquake, describes this cycle. It does not imply that all earthquakes are identical or that the entire fault slips simultaneously. Modern source models recognize heterogeneous rupture, variable friction, barriers, asperities, and dynamic stress changes.

2. Faults and Earthquake Sources

2.1 Fault terminology

A fault is a fracture or fracture zone across which displacement has occurred. The block above an inclined fault surface is the hanging wall; the block below it is the footwall.

A normal fault occurs when the hanging wall moves downward relative to the footwall, generally reflecting extension. A reverse fault occurs when the hanging wall moves upward, generally reflecting compression. A low-angle reverse fault associated with subduction is commonly called a thrust fault. A strike-slip fault involves predominantly horizontal displacement. It is right-lateral or left-lateral according to the apparent direction of motion observed from the opposite block.

Real faults may have oblique-slip motion, combining vertical and horizontal components. A fault trace at the surface is not necessarily the full rupture surface at depth.

2.2 Hypocentre, epicentre, and rupture

The hypocentre or focus is the point at depth where rupture initiates. The epicentre is the point on the Earth’s surface vertically above the hypocentre. The distance between a site and the hypocentre is the hypocentral distance; the distance between a site and the epicentre is the epicentral distance.

An earthquake rupture usually propagates across a finite fault area. Important source parameters include rupture length, rupture width, average slip, rupture velocity, rise time, and fault orientation. The initial rupture point is not necessarily the point of greatest slip or strongest shaking. Directivity can concentrate seismic energy in the direction of rupture propagation, creating unusually strong and long-period motion at some sites.

2.3 Focal mechanism

A focal mechanism represents the orientation and sense of motion of the causative fault. It is commonly inferred from the first motions of P waves recorded at several stations, waveform inversion, or the distribution of seismic radiation. The familiar “beach-ball” diagram contains two nodal planes: the actual fault plane and an auxiliary plane. Without geological or aftershock evidence, seismological data alone may not uniquely identify which plane ruptured.

The principal mechanism types are normal, reverse or thrust, and strike-slip. Focal mechanisms are valuable for tectonic interpretation, seismic hazard models, and selecting realistic scenario earthquakes.

3. Seismic Waves

Earthquake energy propagates through the Earth as elastic waves. Their arrival times, amplitudes, frequencies, and polarization provide information about the source and the intervening medium.

3.1 Body waves

P waves, or primary waves, are compressional waves. Particle motion is approximately parallel to the direction of propagation. They can travel through solids, liquids, and gases and are the fastest seismic waves. They generally arrive first at a recording station and are usually less damaging than S waves for ordinary structures, although large P-wave effects can be important near the source and for sensitive equipment.

S waves, or secondary waves, are shear waves. Particle motion is perpendicular to the direction of propagation. They travel through solids but not through fluids. S waves usually have greater damaging potential than P waves because they generate shear deformation and commonly have larger amplitudes.

For a homogeneous, isotropic elastic medium, approximate wave velocities are

vP = sqrt((K + 4G/3) / rho)

and

vS = sqrt(G / rho),

where K is the bulk modulus, G is the shear modulus, and rho is mass density. These expressions assume linear elastic behaviour and a homogeneous medium. In real geological formations, velocity varies with depth, direction, fractures, saturation, and frequency. Since K + 4G/3 is greater than G, vP is greater than vS for ordinary positive elastic moduli.

3.2 Surface waves

Surface waves travel near the ground surface and generally attenuate with depth. Rayleigh waves involve retrograde or elliptical particle motion in a vertical plane. Love waves involve horizontal shear motion and require suitable layered conditions for their development. Surface waves often have relatively large amplitudes and long durations and can be particularly important for flexible structures and long-period infrastructure.

3.3 Engineering significance of wave characteristics

The damaging effect of a wave is not determined by amplitude alone. Frequency content, duration, phase characteristics, direction, and the relationship between ground-motion periods and structural natural periods are important. A short, high-frequency motion may affect stiff, low-rise structures and non-structural components, whereas long-period motion may be critical for tall buildings, bridges, tanks, and long-span structures.

Seismic waves are modified by geometric spreading, material damping, scattering, reflection, refraction, and conversion between wave types. Local soil layers can amplify selected frequencies and lengthen shaking duration. Thus, an earthquake record at bedrock cannot automatically be used as the expected motion at the foundation level of every structure.

4. Measuring Earthquakes

4.1 Seismometers and accelerometers

A seismometer measures ground motion, commonly velocity or displacement, over a broad frequency range. An accelerometer records ground acceleration and is especially useful for strong-motion engineering applications. A complete record requires accurate timing, calibrated sensors, suitable sampling, dynamic range, and careful processing.

The three principal components are two horizontal components and one vertical component. Instrument orientation, baseline correction, filtering, clipping, and noise must be considered before using a record for analysis.

4.2 Magnitude

Magnitude is a quantitative measure related to the size of an earthquake source. It is logarithmic: an increase of one magnitude unit corresponds to a tenfold increase in the measured wave amplitude used by the relevant scale, although energy increases by approximately a factor of 32 for commonly used magnitude-energy approximations.

The original local magnitude, ML, was developed by Richter for a specified region, instrument, and reference distance. It is not universally reliable for very large earthquakes because of saturation. Other scales include body-wave magnitude, mb, and surface-wave magnitude, Ms, each with frequency and source limitations.

For large earthquakes, moment magnitude, Mw, is preferred because it is tied to physical source parameters. Seismic moment is

M0 = mu A D,

where M0 is seismic moment in newton-metres, mu is the fault-zone shear modulus in pascals, A is the rupture area in square metres, and D is average slip in metres. The commonly used SI relation is

Mw = (2/3)(log10 M0 - 9.1),

where log10 denotes the base-10 logarithm and M0 is expressed in N m. The numerical constant is convention-dependent; alternative unit systems use different constants. The equation assumes the standard seismic-moment definition and an empirically calibrated magnitude conversion.

Magnitude is not a direct measure of shaking at a particular site. Two earthquakes with similar Mw can produce very different motions because of depth, distance, mechanism, rupture directivity, attenuation, site conditions, and path effects.

4.3 Intensity

Intensity describes the effects of an earthquake at a particular location. The Modified Mercalli Intensity (MMI) scale uses observed effects on people, objects, buildings, and the natural environment. Intensity is therefore spatially variable: one earthquake has one magnitude but many intensity values.

Instrumental intensity measures, such as macroseismic intensity derived from recordings, may complement traditional descriptions. Intensity data are useful for historical earthquakes, damage mapping, emergency management, and estimating shaking where instruments were absent. However, observed damage also depends on construction quality, occupancy, soil conditions, and exposure, so intensity is not a pure physical measure.

4.4 Ground-motion parameters

Important engineering measures include peak ground acceleration (PGA), peak ground velocity (PGV), peak ground displacement (PGD), and response spectral acceleration. PGA is the largest absolute acceleration in a record; it is not the same as the acceleration experienced by every structure. A response spectrum shows the maximum response of idealized single-degree-of-freedom oscillators with specified damping over a range of natural periods.

A ground-motion record should be selected or scaled with attention to magnitude, distance, fault mechanism, site class, duration, and period range relevant to the structure. A single parameter such as PGA cannot fully characterize seismic demand.

5. Earthquake Occurrence and Recurrence

5.1 Foreshocks, mainshocks, and aftershocks

An aftershock is an earthquake occurring after a larger event in the same general source region, usually as the crust adjusts to changed stress. A foreshock is an event later identified as preceding a larger earthquake. The designation is retrospective; an earthquake cannot be known with certainty to be a foreshock when it occurs.

Aftershock activity commonly decreases with time. The modified Omori relation is often written as

n(t) = K/(t+c)^p,

where n(t) is the aftershock rate at time t after the mainshock, K is a productivity constant, c accounts for early-time incompleteness or source effects, and p is commonly near 1 but must be estimated for the sequence. The relation is empirical, not a deterministic prediction, and parameters vary among regions and sequences.

5.2 Frequency–magnitude distribution

The Gutenberg–Richter relationship is commonly expressed as

log10 N(M) = a - bM,

where N(M) is the cumulative number of earthquakes with magnitude equal to or greater than M during a specified time and region, and a and b are estimated parameters. The b-value describes the relative proportion of small to large earthquakes. The relationship is valid only over a magnitude range in which the catalogue is complete and the seismicity is reasonably stationary. Apparent changes in b may reflect catalogue incompleteness, changing detection capability, or genuine physical processes.

5.3 Recurrence and probability

A recurrence interval is an average statistical interval, not a timetable. If an annual occurrence rate lambda is assumed constant and events are modelled as a Poisson process, the probability of at least one event in time T is

P(N >= 1) = 1 - exp(-lambda T).

This model assumes independent events and a time-independent rate. It does not account for stress renewal, fault interaction, clustering, elapsed time since the last event, or changing detection thresholds. Renewal and time-dependent models may be more appropriate where fault histories and recurrence intervals are available.

For engineering, hazard is commonly expressed probabilistically as the likelihood that a ground-motion parameter will exceed a specified level during a period. Probabilistic seismic hazard analysis integrates uncertainties in earthquake size, location, recurrence, attenuation, and site response. Deterministic scenarios instead examine specified earthquakes and are useful for understanding credible worst cases. Both approaches have legitimate applications.

6. Earthquake Prediction, Forecasting, and Uncertainty

Reliable short-term prediction requires specifying the time, location, and magnitude of an earthquake. No consistently validated method currently provides such predictions for general operational use. Earthquake forecasting is more defensible: it estimates probabilities over a defined region and time period using seismicity, geological, geodetic, and statistical information.

Global Navigation Satellite System measurements, InSAR, geological mapping, paleoseismology, dense seismic arrays, and machine-learning methods are improving observation and interpretation. AI can assist with phase picking, event detection, catalogue association, waveform classification, and rapid damage mapping. It does not remove physical uncertainty, data bias, or the need for expert validation. Models must be tested on independent data and should report uncertainty rather than only a categorical prediction.

7. Applied Case: The 2001 Bhuj Earthquake

Background

The 26 January 2001 Bhuj earthquake in Gujarat, India, had a moment magnitude commonly reported as approximately Mw 7.7. It occurred in an intraplate region associated with the Kachchh rift system rather than at an active plate boundary. The event caused extensive casualties, building damage, infrastructure disruption, and economic loss across a broad area.

Problem

The earthquake demonstrated that serious hazard exists in regions with lower instrumental seismicity than major plate boundaries. Many buildings had inadequate seismic resistance, poor construction details, heavy roofs or upper stories, weak connections, and irregular configurations. Local ground conditions also contributed to variations in damage.

Analysis and response

Post-earthquake investigations combined seismological records, geological studies, field damage surveys, remote sensing, and structural reconnaissance. The event prompted strengthening of seismic codes and enforcement practices, greater attention to ductile detailing, improved public awareness, and expansion of seismic monitoring and hazard assessment in India.

Outcome and learning points

The principal lesson is that low apparent frequency does not mean negligible risk. Intraplate faults may have long recurrence intervals, and the instrumental catalogue may cover only a small fraction of the relevant geological time. Engineering decisions should therefore combine historical evidence, paleoseismology, active-fault studies, regional tectonics, and conservative treatment of uncertainty. Code compliance, quality control, construction supervision, and community preparedness are as important as theoretical hazard estimates.

8. Professional Practice, Governance, and Limitations

Seismic hazard assessment is multidisciplinary. Geologists investigate faults and Quaternary activity; seismologists analyse sources and waves; geotechnical engineers assess soil response, liquefaction, and slope stability; structural engineers translate hazard into design demand; planners and emergency managers address exposure and recovery.

Important limitations include incomplete earthquake catalogues, uncertain fault geometry, sparse strong-motion records, non-stationary seismicity, model disagreement, and imperfect knowledge of site response. Results should therefore communicate epistemic uncertainty, arising from limited knowledge, separately from aleatory variability, arising from natural randomness.

Professional practice should use current national codes, documented assumptions, peer review for critical facilities, traceable data, and transparent communication of uncertainty. Ethical responsibility includes avoiding false precision, not suppressing inconvenient evidence, and recognizing that vulnerable communities may face disproportionate risk. Regulatory compliance is a minimum requirement, not a guarantee of zero damage.

Best-Practice Workflow for Preliminary Seismic Assessment

1. Define the engineering decision, performance objective, and period of interest.

2. Establish regional tectonic setting and identify potential seismic sources.

3. Compile and assess historical and instrumental catalogues.

4. Investigate active faults, paleoseismic evidence, and geodetic deformation where relevant.

5. Characterize source, path, and site effects.

6. Select appropriate deterministic scenarios or probabilistic methods.

7. Use ground-motion measures relevant to the structure, not PGA alone.

8. Quantify model and parameter uncertainty.

9. Validate assumptions using independent geological, geophysical, and engineering evidence.

10. Update the assessment when new records, maps, codes, or scientific findings become available.

Chapter Summary

Earthquakes result primarily from sudden slip on faults caused by the release of accumulated tectonic strain, although volcanic, induced, and collapse-related earthquakes also occur. Plate tectonics explains broad seismicity patterns, while local faults and geological conditions control many engineering effects. The hypocentre is the rupture-initiation point, the epicentre is its surface projection, and the finite rupture has dimensions, slip, duration, and directivity that influence shaking.

P waves are compressional body waves, S waves are shear body waves, and Rayleigh and Love waves are surface waves. Magnitude describes earthquake size, whereas intensity describes local effects. Moment magnitude is related to seismic moment through source area, average slip, and shear modulus, but magnitude alone does not determine site damage. Ground-motion parameters such as PGA, PGV, PGD, and response spectra provide complementary engineering information.

Earthquake occurrence is described statistically using aftershock relations, Gutenberg–Richter frequency–magnitude distributions, recurrence models, and probability calculations. These models are conditional on assumptions and data quality; they do not provide exact prediction. Modern seismic hazard assessment integrates tectonics, fault studies, monitoring, geotechnical conditions, structural vulnerability, and uncertainty. The Bhuj case demonstrates why intraplate earthquakes and construction quality must be taken seriously. Sound engineering practice combines scientific evidence, appropriate regulation, quality assurance, ethical communication, and continual updating.

Key Terms

Fault: A fracture or fracture zone across which displacement has occurred.

Hypocentre: The point within the Earth where earthquake rupture initiates.

Epicentre: The point on the surface vertically above the hypocentre.

Hanging wall: The block above an inclined fault plane.

Footwall: The block below an inclined fault plane.

Normal fault: A fault in which the hanging wall moves downward relative to the footwall.

Reverse fault: A fault in which the hanging wall moves upward relative to the footwall.

Strike-slip fault: A fault with predominantly horizontal displacement.

Seismic moment: A physical measure of source size equal to shear modulus times rupture area times average slip.

Magnitude: A logarithmic measure related to earthquake source size.

Intensity: A measure of earthquake effects at a particular location.

P wave: A compressional body wave and usually the fastest seismic wave.

S wave: A shear body wave that cannot propagate through fluids.

Rayleigh wave: A surface wave involving elliptical particle motion.

Love wave: A horizontally polarized surface shear wave.

PGA: Peak ground acceleration in a recorded ground-motion component.

Response spectrum: Maximum response of idealized oscillators over a range of periods and specified damping.

Aftershock: An earthquake occurring after a larger event in the same general source region.

Recurrence interval: A statistical average interval between events meeting specified criteria.

Seismic hazard: The probability or expected level of earthquake-induced ground motion at a location.

Review Questions

Short Answer Questions

1. Distinguish the hypocentre from the epicentre.

2. What is the principal difference between magnitude and intensity?

3. Why do S waves generally have greater engineering significance than P waves?

4. Define a strike-slip fault and distinguish right-lateral from left-lateral motion.

5. What is meant by an aftershock?

Long Answer Questions

1. Explain plate tectonics, elastic rebound, and the generation of earthquakes.

2. Describe P waves, S waves, Rayleigh waves, and Love waves, including their engineering significance.

3. Explain seismic moment and moment magnitude, including the meaning and units of each variable.

4. Discuss the Gutenberg–Richter relationship and the limitations of using it for recurrence assessment.

5. Compare deterministic and probabilistic seismic hazard assessment.

Analytical / Application Questions

1. A catalogue reports 100 earthquakes of magnitude 4 or greater and 10 earthquakes of magnitude 5 or greater in the same complete observation period. Estimate the b-value using the Gutenberg–Richter relationship and explain the assumptions involved.

2. An engineering site is near an active strike-slip fault but has few instrumental records. Propose a defensible workflow for estimating seismic hazard, identifying the principal uncertainties.

3. Two sites experience the same Mw 6.5 earthquake, but one records substantially greater damage. Analyse at least five physical and engineering reasons why the observed effects may differ.

Multiple Choice Questions with Answers

Question 1

Which statement best distinguishes magnitude from intensity?

A. Magnitude varies from site to site, while intensity is constant

B. Magnitude describes source size, while intensity describes effects at a location

C. Magnitude applies only to surface waves, while intensity applies only to body waves

D. Magnitude measures building damage directly, while intensity measures fault slip

Correct Answer: B. Magnitude describes source size, while intensity describes effects at a location

Question 2

In a normal fault, the hanging wall generally moves:

A. Upward relative to the footwall

B. Downward relative to the footwall

C. Horizontally rightward only

D. Without any displacement

Correct Answer: B. Downward relative to the footwall

Question 3

Which seismic wave cannot propagate through a fluid?

A. P wave

B. Rayleigh wave

C. S wave

D. Atmospheric wave

Correct Answer: C. S wave

Question 4

The epicentre is:

A. The point of maximum slip on a fault

B. The deepest point of an earthquake

C. The surface point vertically above the hypocentre

D. The end point of rupture propagation

Correct Answer: C. The surface point vertically above the hypocentre

Question 5

Seismic moment is calculated as the product of:

A. Magnitude, distance, and duration

B. Shear modulus, rupture area, and average slip

C. P-wave velocity, density, and frequency

D. PGA, PGV, and PGD

Correct Answer: B. Shear modulus, rupture area, and average slip

Question 6

Which measure is most directly associated with the physical size of a large earthquake source?

A. Modified Mercalli intensity

B. Moment magnitude

C. Epicentral intensity only

D. Building damage grade

Correct Answer: B. Moment magnitude

Question 7

The Gutenberg–Richter b-value primarily describes:

A. The ratio of P-wave to S-wave velocity

B. The relative frequency of small and large earthquakes

C. The duration of aftershock shaking

D. The orientation of a fault plane

Correct Answer: B. The relative frequency of small and large earthquakes

Question 8

A recurrence interval should generally be interpreted as:

A. An exact date for the next earthquake

B. A guaranteed maximum time between earthquakes

C. A statistical average interval under specified assumptions

D. The duration of fault rupture

Correct Answer: C. A statistical average interval under specified assumptions

Question 9

Which factor can amplify or modify bedrock motion at a site?

A. Local soil layering and resonance

B. The alphabetical name of the earthquake

C. The colour of the seismometer

D. The number of buildings outside the region

Correct Answer: A. Local soil layering and resonance

Question 10

Which statement about earthquake prediction is most accurate?

A. Short-term prediction is routinely exact worldwide

B. Forecasting probabilities over defined regions and times is more defensible

C. Magnitude alone identifies the time and location of the next event

D. Aftershocks can always be predicted individually

Correct Answer: B. Forecasting probabilities over defined regions and times is more defensible