Ground Acceleration PGA Calculator
Estimate peak ground acceleration from earthquake magnitude and distance using ground motion prediction equations.
Ground Acceleration PGA Calculator
Calculator
Adjust values & calculateEnter your values below. Every result is computed in your browser โ no data is sent to any server.
Formula: ln(PGA) = c1 + c2(M-6) + c3(M-6)^2 + c4 ln(Reff) + c5 R
Worked example โ PGA = 0.0182g (17.86 Gal) | MMI ~ III-IV | Light shaking, minimal damage
Formula
ln(PGA) = c1 + c2(M-6) + c3(M-6)^2 + c4 ln(Reff) + c5 R
Where M = moment magnitude, R = epicentral distance (km), Reff = sqrt(R^2 + h^2) is the effective distance accounting for source depth, and c1 through c5 are empirical coefficients. Site amplification factors are applied based on NEHRP site class to account for local soil conditions.
Worked Examples
Example 1: Moderate Earthquake at Medium Distance
Problem:Estimate the PGA for a magnitude 6.0 earthquake at 30 km from the epicenter on NEHRP Site Class D (stiff soil).
Solution:Using simplified attenuation: ln(PGA) = -1.715 + 0.5(6.0-6) + (-0.53)(6.0-6)^2 + (-0.778)ln(sqrt(30^2+7^2)) + (-0.0031)(30) = -1.715 + 0 + 0 + (-0.778)ln(30.81) - 0.093 = -1.715 - 2.667 - 0.093 = -4.475 PGA (rock) = e^(-4.475) = 0.0114g Site D amplification: 0.0114 x 1.6 = 0.0182g PGA = 0.0182g = 17.86 cm/sยฒ (Gal)
Result:PGA = 0.0182g (17.86 Gal) | MMI ~ III-IV | Light shaking, minimal damage
Example 2: Strong Earthquake Near Epicenter
Problem:Estimate PGA for a magnitude 7.5 earthquake at 10 km from the fault on Site Class C (dense soil/soft rock).
Solution:Using simplified attenuation: ln(PGA) = -1.715 + 0.5(7.5-6) + (-0.53)(7.5-6)^2 + (-0.778)ln(sqrt(10^2+7^2)) + (-0.0031)(10) = -1.715 + 0.75 - 1.191 + (-0.778)ln(12.21) - 0.031 = -1.715 + 0.75 - 1.191 - 1.949 - 0.031 = -4.136 PGA (rock) = e^(-4.136) = 0.0160g Site C amplification: 0.0160 x 1.2 = 0.0192g
Result:PGA = 0.0192g | These simplified estimates illustrate the method; real GMPEs give higher values for M7.5 near-source
Frequently Asked Questions
What is Peak Ground Acceleration (PGA)?
Peak Ground Acceleration is the maximum acceleration of the ground surface during an earthquake at a specific location. It is typically expressed as a fraction or percentage of the acceleration due to gravity (g = 9.81 m/s squared). PGA is one of the most commonly used parameters in earthquake engineering for characterizing the severity of ground shaking at a site. A PGA of 0.1g means the ground accelerated at one-tenth the acceleration of gravity, while a PGA of 1.0g equals the full force of gravity. Strong earthquakes near their epicenter can produce PGA values exceeding 1.0g, though values of 0.1g to 0.3g are more common for damaging earthquakes at moderate distances.
How do ground motion prediction equations (GMPEs) work?
Ground motion prediction equations, also called attenuation relationships, are empirical models that estimate ground shaking intensity as a function of earthquake magnitude, distance from the source, site conditions, and other parameters. They are developed by analyzing thousands of recorded ground motions from past earthquakes using regression analysis. Modern Next Generation Attenuation (NGA) models like those by Boore, Atkinson, Campbell, and Abrahamson include complex terms for magnitude scaling, geometric spreading, anelastic attenuation, hanging wall effects, and basin depth. These equations predict the median ground motion and its uncertainty, typically expressed as the natural logarithm of PGA or spectral acceleration, allowing engineers to estimate both expected values and probabilities of exceedance.
What is the NEHRP site classification system?
The National Earthquake Hazards Reduction Program site classification system categorizes soil and rock conditions at a site based on the average shear wave velocity in the upper 30 meters (Vs30). Class A represents hard rock with Vs30 greater than 1500 meters per second, typical of eastern US bedrock. Class B is rock with Vs30 between 760 and 1500 m/s. Class C is very dense soil or soft rock at 360 to 760 m/s. Class D is stiff soil at 180 to 360 m/s, which is the default assumption when site data is unavailable. Class E is soft soil below 180 m/s. Class F includes liquefiable soils and peat requiring site-specific analysis. Softer sites generally amplify ground motions, which is why site class significantly affects PGA estimates and building code design requirements.
How is PGA related to the Modified Mercalli Intensity scale?
PGA and Modified Mercalli Intensity (MMI) measure different aspects of earthquake effects but are correlated. PGA is an instrumental measurement of ground acceleration, while MMI is a qualitative scale based on observed damage and human perception, ranging from I (not felt) to XII (total destruction). Empirical relationships by Wald and others relate the two: MMI I to III corresponds to PGA below 0.039g with no damage, MMI IV to V corresponds to 0.039g to 0.092g with light shaking and minor damage, MMI VI to VII corresponds to 0.092g to 0.34g with strong shaking and moderate to considerable damage, and MMI VIII and above corresponds to PGA exceeding 0.34g with severe damage potential.
How is PGA used in seismic building design codes?
Building codes use PGA as a fundamental input for determining seismic design forces. In the International Building Code and ASCE 7, the mapped PGA with a 2 percent probability of exceedance in 50 years (approximately a 2,475-year return period) serves as the basis for calculating design spectral accelerations. The site-adjusted PGA is multiplied by amplification factors to create the design response spectrum, which defines the forces a building must resist at different vibration periods. Short-period design acceleration Sds equals two-thirds of the site-modified spectral acceleration at 0.2 seconds, while long-period design acceleration Sd1 covers longer period structures. These parameters determine the Seismic Design Category, which controls detailing requirements, system selection, and analysis procedures.
What is the difference between the geometric spreading and anelastic attenuation terms?
Geometric spreading is the c4 ln(Reff) term, and it describes energy being distributed over an expanding wavefront, an effect that is purely geometric and therefore identical for every frequency. Anelastic attenuation is the c5 R term, and it describes energy genuinely lost to internal friction and scattering along the travel path, quantified in seismology by the quality factor Q. Because one is logarithmic and the other linear in distance, spreading dominates the first few tens of kilometres while the anelastic term becomes decisive beyond roughly 100 km, where it strips high frequencies fastest. Regional differences in Q are why eastern North American motions carry much further than western ones at the same magnitude.
Why is a fictitious depth added inside the effective distance?
Without it, ln(R) would run to negative infinity as an observer approached the epicentre and the equation would predict unbounded acceleration. Adding a pseudo-depth of a few kilometres, so that Reff = sqrt(R^2 + h^2), floors the distance at h and reproduces the real observation that near-source ground motion saturates: a finite fault radiates from an extended surface rather than a point, and the closest patch of rupture is still kilometres away. The same term stands in for the finite depth of the hypocentre. Its practical consequence is that PGA changes little between a site directly above the source and one a few kilometres away, then falls steeply beyond that.
Which distance measure should be entered, epicentral or fault distance?
Every ground motion prediction equation is fitted to one specific metric and gives biased results with another. Epicentral distance R_epi is measured to the surface point above the hypocentre, hypocentral distance includes focal depth, Joyner-Boore distance R_jb is the shortest horizontal distance to the surface projection of the rupture plane, and rupture distance R_rup is the shortest distance to the rupture surface itself. For a small event at long range they nearly coincide, but for a large earthquake whose rupture is a hundred kilometres long they can differ by tens of kilometres, so a site sitting over the far end of the rupture has a small R_jb despite a large epicentral distance. Ground Acceleration PGA Calculator takes epicentral distance.
What does the scatter around a predicted median mean for design?
A ground motion model predicts a lognormal distribution, not a single value. The aleatory variability, expressed as a standard deviation of the natural logarithm and typically around 0.5 for PGA, means the 84th percentile motion is the median multiplied by e raised to 0.5, roughly 1.65 times, and the 16th percentile is the median divided by the same factor. Actual recordings routinely fall a factor of two either side of the prediction even when magnitude, distance and site class are known exactly. Probabilistic seismic hazard analysis integrates over that distribution together with all possible magnitudes and locations, which is why a code map at a 2 percent chance of exceedance in 50 years shows values well above any single median prediction.
Why does PGA correlate poorly with damage to some structures?
PGA is a single instantaneous peak dominated by high-frequency energy, often a brief spike that carries very little of the earthquake's total energy. Stiff, short-period structures such as low-rise masonry respond to it, but a tall building or a long-span bridge with a natural period of one to several seconds is driven by longer-period energy that PGA does not measure at all. As a rough guide the fundamental period of a regular frame building is on the order of a tenth of a second per storey. Peak ground velocity, spectral acceleration evaluated at the structure's own period, and the duration of strong shaking all correlate better with damage, which is why modern codes design against a response spectrum rather than a single acceleration.
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.
๐งฎDEM Ground Distance Calculator
Calculate dem resolution to ground distance 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.
๐งฎCrustal Density Calculator
Calculate crustal density with inputs, formulas, and instant results.
๐งฎCrustal Thickness Calculator (H-k)
Calculate crustal thickness from receiver functions with inputs, formulas, and instant results.
๐งฎEarthquake Magnitude to Energy Calculator
Calculate earthquake magnitude to energy with inputs, formulas, and instant results.
๐งฎGutenberg-Richter Recurrence Calculator
Calculate earthquake recurrence (gutenbergโrichter) with inputs, formulas, and instant results.