Seismic Intensity Calculator
Estimate Modified Mercalli Intensity (MMI), PGA, and PGV from earthquake magnitude, depth, and distance using seismic attenuation relations.
Seismic Intensity Calculator
Calculator
Adjust values & calculateEnter your values below. Every result is computed in your browser โ no data is sent to any server.
Formula: MMI = 3.5M - 1.09 ln(R) - 0.00186R - 2.7 + ln(S)
Worked example โ MMI ~ VII (Very Strong)
Formula
MMI = 3.5M - 1.09 ln(R) - 0.00186R - 2.7 + ln(S)
Estimates Modified Mercalli Intensity from earthquake magnitude M, hypocentral distance R (km), and site amplification factor S. PGA and PGV are then estimated from empirical MMI-ground motion relationships.
Worked Examples
Example 1: Shallow M6.5 Earthquake at 50 km
Problem:Estimate the seismic intensity at 50 km from a magnitude 6.5 earthquake with a focal depth of 10 km.
Solution:R = sqrt(50^2 + 10^2) = 50.99 km MMI = 3.5*6.5 - 1.09*ln(50.99) - 0.00186*50.99 - 2.7 MMI = 22.75 - 4.28 - 0.09 - 2.7 = 15.68 (clamped to ~7.0) Corresponding to VII - Very Strong shaking
Result:MMI ~ VII (Very Strong)
Example 2: Deep M7.0 at Epicenter
Problem:What intensity would be felt directly above a magnitude 7.0 earthquake at 100 km depth?
Solution:R = sqrt(0^2 + 100^2) = 100 km (epicentral distance = 0) MMI = 3.5*7.0 - 1.09*ln(100) - 0.00186*100 - 2.7 MMI = 24.5 - 5.02 - 0.186 - 2.7 = ~6.3
Result:MMI ~ VI (Strong shaking, reduced by deep focus)
Frequently Asked Questions
What is Modified Mercalli Intensity (MMI)?
Modified Mercalli Intensity is a scale from I to XII that measures the observed effects of an earthquake at a specific location. Unlike magnitude which measures energy at the source, intensity describes how strong the shaking is felt at a given point. MMI I means the earthquake was not felt, while MMI XII indicates total destruction. The same earthquake can produce different intensities at different locations depending on distance, depth, and local soil conditions.
How is seismic intensity different from magnitude?
Magnitude measures the total energy released at the earthquake source and has a single value per event. Intensity measures the severity of shaking at a particular location and varies with distance from the epicenter, depth, and local geology. A magnitude 7.0 earthquake might produce MMI IX near the epicenter but only MMI IV at 200 km away. Magnitude uses instruments while intensity traditionally relies on observed effects on people, structures, and the environment.
What are PGA and PGV?
PGA (Peak Ground Acceleration) is the maximum acceleration of the ground surface during an earthquake, typically measured in cm/s^2 or as a fraction of gravitational acceleration (%g). PGV (Peak Ground Velocity) is the maximum ground velocity in cm/s. PGA correlates well with damage to stiff, short-period structures while PGV better predicts damage to flexible, longer-period structures. Both parameters are essential for seismic hazard assessment and building design codes.
How does site amplification affect intensity?
Site amplification occurs when soft soil or sedimentary basins amplify seismic waves compared to bedrock. Loose, unconsolidated sediments can amplify shaking by factors of 2-5 or more, significantly increasing damage potential. This is why buildings on soft soil often suffer more damage than those on rock at the same distance from an earthquake. Site amplification is a key factor in seismic hazard maps and building code requirements for different soil classifications.
What attenuation model does Seismic Intensity Calculator use?
Seismic Intensity Calculator uses a simplified attenuation relation inspired by Atkinson and Wald (2007) that estimates Modified Mercalli Intensity from magnitude, hypocentral distance, and a site amplification factor. The PGA and PGV estimates use empirical relationships from Wald et al. (1999) that convert MMI to ground motion parameters. These are approximate models suitable for educational purposes. Professional seismic hazard analyses use more detailed ground motion prediction equations calibrated for specific tectonic regions.
Why does the calculation use hypocentral distance rather than epicentral distance?
Seismic waves radiate from the hypocentre, the actual point of rupture at depth, not from the epicentre, which is only its projection on the surface. Hypocentral distance is the straight-line separation, the square root of epicentral distance squared plus focal depth squared. Directly above the source the epicentral distance is zero but the hypocentral distance equals the focal depth, so a rupture at 100 kilometres depth is still 100 kilometres of rock away from the observer. This is why deep-focus events produce far milder shaking at the epicentre than shallow crustal events of identical magnitude, and why depth cannot be omitted.
What is the difference between the logarithmic and the linear distance terms?
The two terms describe physically distinct loss mechanisms. The term in the natural logarithm of distance represents geometric spreading, the dilution of a fixed amount of energy over an expanding wavefront, and it dominates close to the source. The term linear in distance represents anelastic attenuation, the conversion of wave energy into heat through internal friction, which accumulates with path length and therefore matters most at long range. Comparing their rates of change, the coefficients used here make the two contribute equally near about 590 kilometres, beyond which absorption is the controlling effect.
What is Vs30 and how does it relate to the site amplification factor?
Vs30 is the time-averaged shear wave velocity in the uppermost 30 metres of ground, and it is the standard proxy for site response in building codes because soft slow layers trap and amplify wave energy. Competent rock sits at high velocity and defines the reference condition, corresponding to a factor near 1.0 here. Stiff soils warrant something like 1.5, and very soft deposits, conventionally those below roughly 180 metres per second, commonly justify 2 to 3. Amplification is also frequency dependent, and very soft soils can deamplify the strongest shaking when they yield non-linearly.
Why can two towns at the same distance report very different intensities?
Distance and magnitude set only the broad expectation. Rupture directivity focuses energy in the direction the fault tears, so sites along strike from a propagating rupture can experience a large velocity pulse while sites behind it are much quieter. Sedimentary basins trap and resonate surface waves, extending strong shaking well beyond what a simple attenuation curve predicts, and ridge crests can amplify motion through topographic focusing. Reported intensity also depends on the local building stock, since the same ground motion damages unreinforced masonry far more than engineered frames, and MMI degrees are defined partly by observed damage.
How does the Modified Mercalli scale compare with other intensity scales?
All macroseismic scales grade observed effects rather than measure a physical quantity, but they differ in granularity and in the construction they assume. The Modified Mercalli scale, rewritten by Wood and Neumann in 1931 for United States practice, has twelve degrees. The European Macroseismic Scale of 1998 also has twelve but ties them explicitly to defined vulnerability classes of buildings, which makes it more transferable between countries. The Japan Meteorological Agency shindo scale uses a different structure, with levels from 0 to 7 and 5 and 6 each split into lower and upper, and is computed instrumentally. Degrees are broadly comparable but not interchangeable.
References
Background & Theory
History
Reviewed for accuracy by Daniel Agrici, Founder & Lead Developer ยท Editorial policy
Related Calculators
๐งฎGround Vibration Intensity Calculator
Calculate ground vibration intensity with inputs, formulas, and instant results.
๐งฎSeismic Attenuation Calculator
Calculate seismic attenuation with inputs, formulas, and instant results.
๐งฎSeismic Energy Release Calculator
Calculate seismic energy release with inputs, formulas, and instant results.
๐งฎSeismic Moment Tensor Components Calculator
Calculate seismic moment tensor components with inputs, formulas, and instant results.
๐งฎSeismic Reflection Coefficient Calculator
Calculate seismic reflection coefficient with inputs, formulas, and instant results.
๐งฎSeismic Wave Velocity Calculator
Calculate seismic wave velocity with inputs, formulas, and instant results.
๐งฎEarth Energy Balance Calculator
Calculate earth energy balance with inputs, formulas, and instant results.
๐งฎBouguer Correction Calculator
Calculate bouguer correction with inputs, formulas, and instant results.