Earthquake Energy: How Much Power Does a Magnitude 9 Release?
On May 22, 1960, the ground beneath southern Chile ruptured along a fault zone 1,000 kilometers long. The rupture lasted approximately ten minutes. In that time, the Valdivia earthquake released more seismic energy than every other earthquake recorded in the twentieth century combined. The moment magnitude was 9.5 β the largest instrumentally recorded earthquake in human history. Tsunamis generated by the rupture crossed the Pacific, killing people in Hawaii, Japan, and the Philippines. The ground oscillated for days afterward, ringing like a struck bell at periods of hundreds of seconds as Earth's free oscillations β the planet's own natural resonant frequencies β slowly damped away.
The number 9.5 on the moment magnitude scale carries almost no intuitive weight for most people. It is a single digit larger than 8.5, which sounds manageable enough. But the magnitude scale is logarithmic in a very specific way, and the energy difference between those two numbers is not small. The 1960 Chile earthquake released approximately 32 times more seismic energy than an M8.5, roughly 1,000 times more than an M8.0, about 32,000 times more than an M7.5, and somewhere in the neighborhood of one million times more energy than the atomic bomb dropped on Hiroshima. These numbers are not metaphors or approximations β they follow directly from the mathematics of the magnitude scale and the physics of fault rupture, and understanding them is fundamental to understanding what earthquakes actually are.
This post is about those numbers: what they mean, where they come from, how seismologists actually measure earthquake energy, and why the raw energy figure β impressive as it is β tells you less about destruction than you might expect.
The Magnitude Scale: A Brief History of a Misunderstood Number
Charles Richter introduced his magnitude scale in 1935 not as a measure of energy but as a practical way to compare the sizes of Southern California earthquakes using the amplitude of seismic waves recorded on a specific instrument β the Wood-Anderson torsion seismograph β at a standard distance of 100 kilometers. The Richter scale was local, instrumental, and explicitly designed to be a relative ranking rather than an absolute physical measurement. Richter himself was emphatic that magnitude was not a direct measure of energy, a point that has been consistently lost in popular coverage of earthquakes ever since.
The original Richter scale worked well for moderate California earthquakes but broke down at large magnitudes and for earthquakes far from the recording stations. Over the following decades, seismologists developed a family of magnitude scales β body wave magnitude (mb), surface wave magnitude (Ms), and others β each optimized for different earthquake sizes and distance ranges, each with its own saturation problem at the top end of the scale where the physical size of the rupture outgrows the sensitivity of the seismic waves being measured.
The Moment Magnitude Scale: The Modern Standard
The modern standard β moment magnitude, Mw β was introduced by Hiroo Kanamori in 1977 and is grounded directly in the physics of fault rupture rather than in the amplitude of any particular seismic wave type. Moment magnitude is derived from the seismic moment M0, a quantity with units of force times distance (Newton-meters) that captures three fundamental properties of the rupture: the area of the fault that slipped, the average amount of slip across that area, and the shear modulus of the rock β its resistance to the shearing deformation that fault slip represents.
βοΈ The Seismic Moment Formula
Seismic moment is defined as M0 = ΞΌ Γ A Γ D, where ΞΌ is the shear modulus of the crust (typically ~30 GPa for crustal rock), A is the rupture area in square meters, and D is the average slip displacement in meters. For the 2011 Tohoku earthquake, A was approximately 400 km Γ 200 km = 8 Γ 1010 mΒ², D was roughly 10β60 meters of slip (averaging perhaps 20 m), and ΞΌ β 30 Γ 109 Pa. That gives M0 β 3.0β5.0 Γ 1022 NΒ·m β consistent with the independently determined value of ~5.3 Γ 1022 NΒ·m that corresponds to Mw 9.1. The moment magnitude formula is: Mw = (2/3) Γ log10(M0) β 6.07.
The moment magnitude scale does not saturate β it remains physically meaningful for earthquakes of any size, from the smallest detectable microseismic events at Mβ2 to the largest conceivable megathrust ruptures at M9.5 or beyond. It is the scale used in all modern seismological reporting, and when news outlets report an earthquake's magnitude today, the number given is almost always Mw β even when they call it "the Richter scale," a habit that persists despite Richter's own scale having been retired from scientific use decades ago.
The Logarithmic Scale: What Each Unit Actually Means
The moment magnitude scale is logarithmic, but in a way that is slightly more complex than the simple "each unit is 10 times larger" rule that often appears in popular science writing. The relationship between magnitude and seismic moment β the fundamental physical quantity β is logarithmic with base 10 and a factor of 2/3 in the exponent. The consequence is that each unit increase in magnitude corresponds to a factor of 103/2 β 31.6 increase in seismic moment.
Energy scales similarly, though the relationship between seismic moment and radiated seismic energy introduces an additional complication. Not all of the seismic moment is converted to radiated seismic waves β a significant fraction goes into heat generated by friction on the fault surface, permanent deformation of the surrounding rock, and other inelastic processes. The fraction that becomes radiated seismic energy (what actually travels outward as seismic waves and shakes buildings) is expressed by the apparent stress or the radiation efficiency of the earthquake, which varies between events.
π The Energy-Magnitude Relationship
The empirical relationship between radiated seismic energy Es (in joules) and moment magnitude Mw is approximately: log10(Es) = 1.5 Γ Mw + 4.8. This means each unit increase in Mw multiplies radiated energy by 101.5 β 31.6. A two-unit increase multiplies energy by 31.6Β² β 1,000. A two-unit increase in magnitude β say from M7 to M9 β therefore corresponds to about 1,000 times more radiated seismic energy. The total seismic moment scales even more steeply: each magnitude unit corresponds to a factor of ~31.6 in moment, or roughly 1,000 per two units, the same ratio that governs the energy relationship.
Working through the numbers concretely: an M5.0 earthquake radiates approximately 2 Γ 1012 joules of seismic energy β roughly equivalent to the energy released by a small tactical nuclear weapon, or about 500 tons of TNT. An M7.0 radiates about 2 Γ 1015 joules β 500,000 tons of TNT, comparable to a large thermonuclear weapon. An M9.0 radiates approximately 2 Γ 1018 joules β roughly 500 million tons of TNT, or about 25,000 times the yield of the largest nuclear weapon ever detonated (the Soviet Tsar Bomba, at approximately 50 megatons). These are radiated energy figures only β they do not include the energy dissipated as heat on the fault, which is typically several times larger.
Seismic Moment vs. Radiated Energy: An Important Distinction
One of the most common sources of confusion in earthquake energy discussions is the conflation of seismic moment β the physical measure of fault slip β with radiated seismic energy β the energy that actually travels as seismic waves. They are related but distinct quantities, and the ratio between them carries important physical information about the rupture process.
A fault that slips very slowly β like the slow earthquakes discussed in seismological research β can have an enormous seismic moment (reflecting large slip on a large fault area) while radiating almost no seismic energy, because the slow slip velocity produces no dynamic stress drop and no high-frequency seismic radiation. Conversely, a small, sharp earthquake on a highly stressed fault segment can radiate a relatively high fraction of its stored elastic energy as seismic waves, producing strong shaking disproportionate to its moment.
The ratio of radiated seismic energy to seismic moment β the apparent stress divided by the shear modulus, or the radiation efficiency β varies by roughly two orders of magnitude across the population of earthquakes. Typical crustal earthquakes radiate roughly 5 Γ 10β5 of their seismic moment as energy. Subduction megathrust earthquakes, which occur on shallow, warm, fluid-saturated fault interfaces, tend to radiate a somewhat smaller fraction β they are mechanically "slower" in character than continental strike-slip earthquakes of equivalent moment, producing somewhat less high-frequency shaking per unit of seismic moment.
Where Does the Rest of the Energy Go?
For a typical crustal earthquake, the energy budget breaks down roughly as follows. The total elastic strain energy released by the rupture is partitioned between radiated seismic energy (the waves that travel outward and cause shaking, typically 5β10% of the total), frictional heat generated on the fault surface during slip (the largest single component, typically 50β80%), and fracture energy consumed in creating new crack surfaces and deforming the rock immediately adjacent to the fault (5β20%). A small fraction goes into permanent displacement of rock and topographic change. The heat generated on the fault during a major earthquake is enormous β the Tohoku rupture generated enough frictional heat to raise the temperature of the fault zone by several hundred degrees Celsius, detectable in borehole temperature measurements made years later.
The Great Earthquake Catalog: Energy in Numbers
Putting the largest earthquakes in the instrumental record side by side against familiar reference points makes the energy scale viscerally comprehensible in a way that abstract joule counts cannot.
| Event | Mw | Radiated Energy (J) | TNT Equivalent |
|---|---|---|---|
| 1960 Valdivia, Chile | 9.5 | ~1.1 Γ 1019 | ~2.6 billion tons |
| 1964 Alaska (Good Friday) | 9.2 | ~1.4 Γ 1018 | ~330 million tons |
| 2004 Sumatra-Andaman | 9.1 | ~1.1 Γ 1018 | ~260 million tons |
| 2011 Tohoku, Japan | 9.1 | ~1.9 Γ 1018 | ~450 million tons |
| 2010 Maule, Chile | 8.8 | ~2.5 Γ 1017 | ~60 million tons |
| 1906 San Francisco | 7.9 | ~2.0 Γ 1015 | ~480,000 tons |
| Soviet Tsar Bomba (1961) | β | ~2.1 Γ 1017 | 50 million tons |
| Hiroshima bomb (1945) | β | ~6.3 Γ 1013 | ~15,000 tons |
Several things stand out from this table. The 1960 Chile earthquake released roughly five times more energy than the 2011 Tohoku event, despite both being listed as the largest earthquakes of their respective centuries. The Tsar Bomba β the most powerful nuclear weapon ever detonated β released about as much energy as the 2010 Chile earthquake and roughly one-fifth the energy of Tohoku. The 1906 San Francisco earthquake, which destroyed a major city, falls several hundred times below the yield of the Tsar Bomba on the energy scale. And Hiroshima β the cultural reference point most people use to anchor their sense of catastrophic energy release β is overwhelmed by even a modest M7 continental earthquake.
Annual Global Seismic Energy Budget
The total seismic energy released by all earthquakes worldwide in a given year is dominated, overwhelmingly, by the largest events. In a typical year without a great earthquake, the annual global seismic energy release is on the order of 1017 to 1018 joules β comparable to a handful of M8 events and their aftershocks. In years that contain a great earthquake, the budget is transformed: the 2011 Tohoku event alone accounted for a significant fraction of the total seismic energy released globally in the entire twentieth century up to that point.
This concentration of energy in rare, extreme events is a direct consequence of the Gutenberg-Richter frequency-magnitude relationship. The number of earthquakes decreases by a factor of roughly ten for each unit increase in magnitude β there are about ten times as many M6s as M7s, ten times as many M7s as M8s, and so on. But the energy released per event increases by a factor of 31.6 per magnitude unit. The net effect is that the energy budget skews enormously toward the largest earthquakes: approximately 75β90% of all seismic energy released in any given decade comes from events M8.5 and larger, which may occur only a handful of times per decade.
β‘ Global Seismic vs. Solar Energy
The total seismic energy released by all earthquakes worldwide in a typical year β roughly 1017β1018 joules β sounds enormous until compared to other planetary energy flows. Solar energy reaching Earth's surface totals approximately 5.5 Γ 1024 joules per year. Global annual seismic energy release is therefore roughly 10 million to 100 million times smaller than the solar energy intercepted by Earth in the same period. Earthquakes are catastrophic on a human scale but energetically minor perturbations in Earth's overall thermodynamic budget β most of the planet's internal heat escapes through slow conduction and volcanic activity, not through seismic rupture.
Earth's Free Oscillations: The Planet Ringing Like a Bell
The most spectacular demonstration of the raw energy in a great earthquake is not the immediate ground shaking but what happens afterward. Large earthquakes excite the free oscillations of Earth β the natural resonant modes at which the entire planet vibrates, with periods ranging from a few minutes to nearly an hour. These normal modes of Earth were theoretically predicted decades before they were observed, but the 1960 Chile earthquake provided the first clear instrumental detection: seismographs worldwide recorded the planet ringing at its characteristic frequencies for weeks after the mainshock.
Earth's free oscillations fall into two families. Spheroidal modes involve radial motion β the planet breathing, or expanding and contracting at its resonant frequencies. Toroidal modes involve twisting motion, with the solid Earth rotating back and forth slightly around its axis. The fundamental spheroidal mode, designated βSβ, has a period of approximately 54 minutes β meaning the entire planet subtly expands and contracts 54 minutes apart with each oscillation. The fundamental toroidal mode, βTβ, has a period of approximately 44 minutes. These are extraordinarily long-period oscillations that can only be excited by earthquakes large enough to deform the entire planet.
For the 2004 Sumatra-Andaman earthquake, free oscillations were detected on superconducting gravimeters worldwide and were still measurable two weeks after the mainshock. The amplitude of the oscillations provided an independent check on the earthquake's seismic moment, and the pattern of which modes were most strongly excited constrained the spatial distribution of slip along the 1,300-kilometer rupture. Free oscillation analysis has since become a standard tool for characterizing the largest earthquakes and for probing Earth's internal structure β the frequencies and decay rates of the normal modes depend on the density, elasticity, and attenuation of the materials through which they propagate.
Why Energy Alone Does Not Predict Destruction
If raw seismic energy determined earthquake damage, the 1960 Chile earthquake β by far the most energetic in the instrumental record β would have been by far the most destructive. It was not. The 2010 Haiti earthquake, with a moment magnitude of 7.0, released roughly 500,000 times less seismic energy than the 1960 Chile event but killed approximately 160,000 people compared to Chile's estimated 1,655 fatalities. The 1976 Tangshan earthquake in China, M7.8, killed approximately 240,000 people. Energy and destruction are related, but several intervening factors β each as important as the raw energy budget β determine what a given earthquake actually does to people and infrastructure.
Depth
Seismic energy spreads outward from the hypocenter as it travels through the Earth, and its intensity at the surface decreases with the square of the distance from the source. A shallow earthquake with its hypocenter at 5β10 km depth delivers far more energy per unit area to the surface than a deep earthquake of identical magnitude with its hypocenter at 50 or 100 km. The Haiti earthquake occurred at approximately 13 km depth directly beneath a densely populated city. The 2013 M8.0 Santa Cruz Islands earthquake occurred at 10 km depth in a remote region and caused minimal casualties despite releasing roughly 1,000 times more energy than Haiti.
Distance from Population Centers
The 1960 Chile earthquake struck a sparsely populated section of southern Chile. The energy was enormous, but the population density in the immediate rupture zone was low. The 2010 Haiti earthquake struck 25 km southwest of Port-au-Prince, a capital city of more than two million people with virtually no earthquake-resistant construction. The proximity of the hypocenter to a densely populated, seismically vulnerable urban area is the dominant factor in casualty counts, overwhelmingly outweighing differences in total energy release for events in the M7βM8 range.
Local Geology and Site Amplification
Seismic waves are amplified dramatically when they pass from hard rock into soft sediments. Mexico City, built on the drained bed of an ancient lake, experienced catastrophic amplification during the 1985 MichoacΓ‘n earthquake β ground motion in the lake bed sediments was amplified by a factor of 20β50 relative to surrounding bedrock, and the resonant period of the sediment column matched the dominant period of the seismic waves, producing sustained, large-amplitude oscillations that collapsed buildings designed without this resonance in mind. The source earthquake was more than 350 km away and M8.0 β but soft basin sediments transformed distant shaking into a local disaster.
Building Stock and Construction Quality
The single most powerful predictor of earthquake casualties in the modern era is the quality and type of building construction. Unreinforced masonry β adobe, brick, and stone construction without steel reinforcement or engineered connections β performs catastrophically in earthquakes. It is widespread in precisely the parts of the world most exposed to seismic hazard: the Middle East, South Asia, Central Asia, and Latin America. Japan, with comparably high seismic exposure, has invested massively in engineered construction, retrofitting, and building codes, and routinely experiences major earthquakes with far lower casualty rates than countries with equivalent shaking but weaker building stock.
Duration, Frequency Content, and Structural Resonance
Beyond total energy, two seismological properties of the shaking matter enormously for structural damage: the duration of strong ground motion and the frequency content of the seismic waves relative to the natural frequencies of the structures being shaken.
Duration matters because structures that survive the first few cycles of shaking may be progressively weakened by continued motion. Low-cycle fatigue β the degradation of structural connections and materials under repeated large-amplitude loading β is responsible for many collapses that occur in the later stages of a long rupture rather than in the initial seconds. Great earthquakes have rupture durations of two to ten minutes, compared to 10β30 seconds for M7 events, and this extended shaking duration elevates the damage potential of large events beyond what peak ground acceleration alone would predict.
Frequency content determines which structures are most at risk. Short, stiff structures β one to three stories β resonate at high frequencies (1β10 Hz) and are most vulnerable to the high-frequency content of near-field seismic waves. Tall, flexible structures β ten stories and above β resonate at low frequencies (0.1β1 Hz) and are most vulnerable to the long-period energy that large, distant earthquakes generate in abundance. This is why tall buildings in Los Angeles are specifically engineered against the long-period energy expected from a great Cascadia subduction zone earthquake several hundred kilometers away, while shorter masonry buildings face their primary risk from nearby, moderate events.
Measuring Energy in Real Time: How Seismologists Do It
Determining the energy of an earthquake in real time β within minutes of the event β is a problem with significant operational implications. Earthquake early warning systems, tsunami warning centers, and emergency management agencies all need rapid, accurate magnitude estimates to issue timely alerts. The challenge is that the most accurate magnitude measurement β seismic moment from long-period waveform analysis β requires waiting for the slowest seismic waves to arrive at distant stations, which takes ten to twenty minutes for a global network. Faster measurements are possible but less accurate.
P-wave Amplitude Methods
The first energy arriving at a seismic station from any earthquake is the P-wave β the compressional wave that travels through the body of the Earth at 6β8 km/s. Earthquake early warning systems exploit the observation that the amplitude and frequency content of the initial P-wave arrival correlate with the final earthquake magnitude, allowing rapid magnitude estimates within 3β10 seconds of rupture initiation. These estimates are used to issue alerts before the destructive S-waves and surface waves arrive β providing seconds to tens of seconds of warning, enough to automatically halt trains, open fire station doors, and alert the public to take cover.
W-phase and Long-Period Methods
For great earthquakes, the most reliable moment magnitude determination comes from analysis of the W-phase β a long-period wave that arrives between the P-wave and the surface wave and whose amplitude at periods of 100β1,000 seconds is directly proportional to seismic moment. The W-phase method was developed specifically for rapid characterization of large subduction zone earthquakes and is now used operationally by tsunami warning centers worldwide. A reliable W-phase magnitude estimate for the Tohoku earthquake was available within approximately 20 minutes of the mainshock β fast enough to inform tsunami warning decisions.
The Upper Bound: How Large Can an Earthquake Get?
A natural question when contemplating the energy of the largest known earthquakes is whether there is a physical upper limit β and if so, where it lies. The answer involves the geometry of the fault systems capable of producing great earthquakes and the mechanical constraints on how much elastic strain energy can accumulate before rupture.
The largest earthquakes occur on subduction zone megathrusts β the interfaces between converging tectonic plates β because these are the longest and broadest fault surfaces on Earth. The Valdivia rupture was approximately 1,000 km long and 200 km wide. Longer ruptures are possible in principle β the 2004 Sumatra-Andaman earthquake ruptured 1,300 km of the Sunda megathrust β but the total length of continuous, locked subduction interface limits the maximum rupture length, and therefore the maximum seismic moment, for any given subduction zone.
Theoretical estimates place the upper bound for subduction megathrust earthquakes at approximately M9.5βM10.0, corresponding to ruptures of the longest continuous megathrust segments on Earth (the Japan Trench, the Cascadia subduction zone, the Tonga-Kermadec system) with maximum slip. Whether M10 is physically achievable is uncertain β it would require simultaneous rupture of an extraordinarily long fault segment with unusually large slip, and no event approaching that scale is known in the geological or instrumental record. The 1960 Chile M9.5 remains the upper benchmark, and extrapolating beyond it requires geological evidence rather than empirical scaling.
π Could All Faults Rupture Simultaneously?
A purely theoretical exercise: if every major fault on Earth ruptured simultaneously, releasing the total accumulated elastic strain energy stored in the global seismogenic crust, the total energy release would be on the order of 1022β1023 joules β equivalent to roughly 10,000 great earthquakes occurring at once. While physically impossible as a single event, this thought experiment illustrates how the elastic energy stored in Earth's crust dwarfs the output of any individual earthquake and underscores why fault locking and episodic rupture β rather than continuous aseismic creep β represents the fundamental mechanism by which tectonic strain energy is episodically catastrophically released.
The Gutenberg-Richter Relationship and Energy Statistics
The Gutenberg-Richter law β the observation that earthquake frequency decreases by a factor of ten for each unit increase in magnitude β applies remarkably consistently across tectonic settings, magnitude ranges, and time periods. It describes a power-law distribution of earthquake sizes, implying that there is no characteristic earthquake size and that the fault system is in a state of self-organized criticality β operating near a threshold where a small initial rupture can cascade, with some probability, to any larger size.
The energy implications of the Gutenberg-Richter relationship are profound. Because each magnitude unit represents a factor of 31.6 in energy but only a factor of 10 in frequency, the energy per unit time released by earthquakes of a given magnitude band increases with magnitude β there is no magnitude range that dominates the energy budget except the very largest events. The extreme right tail of the earthquake size distribution carries the vast majority of the total energy, which is why a single great earthquake can release more seismic energy than all smaller earthquakes combined over decades.
Conclusion
The energy of a magnitude 9 earthquake defies easy intuition. A billion tons of TNT is not a concept that scales from human experience. The abstractions β joules, Newton-meters, seismic moments β are precise and physically meaningful, but they require translation to land with appropriate force. The most useful translation may simply be this: the 1960 Chile earthquake released more energy than all other earthquakes in the twentieth century combined, set the entire planet ringing at its natural resonant frequencies for weeks, generated tsunamis that crossed an ocean and killed people 17,000 kilometers away, and was caused by ten minutes of slip on a fault β the same physical process, governed by the same rate-and-state friction laws, that produces the microseismic background noise of thousands of imperceptible M1 events every day.
The magnitude scale compresses this 10-billion-fold range of energy into a single-digit number, which is both its greatest convenience and its greatest liability. An M9 and an M5 look almost identical on a news ticker. The physics separating them is the difference between a tremor you might feel in your office chair and an event that deforms the planet's rotational axis, shortens the length of the day, and permanently rearranges the gravity field of an entire ocean basin. Understanding the logarithm β really understanding it, not just accepting it intellectually β is the key to reading earthquake news with an accurate sense of what is actually being reported.
Support Earthquake Radar
Earthquake Radar provides free, real-time earthquake monitoring and in-depth educational resources to help communities understand seismic risk. If this guide was useful, please consider supporting our mission:
Your support helps maintain free earthquake monitoring services and fund more educational resources for communities worldwide.
Twitter/X