Key Equations
$M_w = \frac{2}{3}(\log_{10}M_0 - 9.1)$$\log_{10}E = 1.5M + 4.8$
$\log_{10}\text{PGA}\approx 0.5M - 1.85 - 1.68\log_{10} R$
Adjust magnitude, focal depth, and epicentral distance to instantly compute PGA, JMA seismic intensity, and MMI scale. Visualize seismic wave propagation from hypocenter to surface.
The core relationship is between the Moment Magnitude ($M_w$) and the Seismic Moment ($M_0$). The seismic moment is calculated from the fault's physical properties: area, average slip, and the rigidity of the rock. $M_w$ is then derived from it, creating a scale that doesn't saturate for very large earthquakes.
$$M_w = \frac{2}{3}(\log_{10}M_0 - 9.1)$$Where $M_0$ is the seismic moment in Newton-meters (N·m). The constant 9.1 is a scaling factor to make the magnitude values roughly consistent with older scales.
To estimate how shaking weakens with distance, we use an attenuation relationship. A simplified form predicts the Peak Ground Acceleration (PGA) based on the earthquake's magnitude and your distance from it.
$$\log_{10}(\text{PGA}) \approx 0.5M_w - 1.85 - 1.68\log_{10}R$$Here, PGA is in %g, $M_w$ is the moment magnitude, and $R$ is the epicentral distance in kilometers. The term $-1.68\log_{10} R$ models the geometric spreading and damping of seismic waves as they travel through the Earth.
Earthquake-Resistant Building Design: Engineers use PGA values from attenuation models to define the design basis earthquake for a region. The seismic intensity scales (like JMA or MMI) help create building codes that specify different construction requirements for zones expecting moderate vs. severe shaking.
Emergency Response & Rapid Assessment: Immediately after a quake, rapid estimates of magnitude and predicted intensity distribution (ShakeMaps) are generated. This allows emergency services to prioritize areas likely to have suffered the worst damage and direct rescue resources effectively.
Public Communication & Education: Translating magnitude into TNT equivalent or comparing it to historical quakes (e.g., "This M6.5 quake released energy similar to 20 Hiroshima bombs") helps the public grasp the immense power involved, fostering better preparedness.
Insurance & Risk Modeling: Insurance companies use probabilistic seismic hazard analysis, which relies on magnitude-frequency relationships and attenuation laws, to model financial risk and set premiums for properties in earthquake-prone areas.
There are several key points you should be aware of when using this tool, especially if you're considering practical applications. First, "the calculation results are only a guideline for average ground conditions". The tool's formulas assume generic bedrock (basement rock). In reality, the shaking can often be amplified 2 to 3 times on softer soil layers (alluvium) deposited on top of this bedrock. For example, even at the same 20km distance from the epicenter, while a hard mountainous area might calculate as a seismic intensity of 5 Lower, a reclaimed land or valley plain could potentially experience an intensity of 5 Upper or higher. Next, note that the relationship between magnitude and seismic intensity is not a simple proportion. Just because an M7 earthquake has 10 times the energy of an M6, it doesn't mean the seismic intensity increases by one full grade. At locations sufficiently far from the hypocenter, the intensity may hardly change even if the magnitude increases by 1. Finally, the difference between "epicentral distance" and "hypocentral distance". The tool's input is "Epicentral Distance R", but the actual shaking is determined by the straight-line "Hypocentral Distance" from the earthquake's focus. If you input 0 km for the epicentral distance for an earthquake with a depth of 50 km, it will be calculated using a hypocentral distance of 50 km. This difference can be ignored for calculations at points far from the epicenter, but you need to be mindful of it when dealing with shallow earthquakes near the epicenter.
The calculation logic behind this simulator is actually applied as a foundation in various engineering fields. The first to mention are earthquake engineering and seismic design. The Peak Ground Acceleration (PGA) calculated here is the first step in creating response spectra for designing buildings and bridges. Next is geotechnical engineering. To quantify the aforementioned ground amplification effects, the shear wave velocity structure of the ground is investigated, and detailed site response analyses using wave propagation theory are conducted. Furthermore, in the fields of risk engineering and urban disaster prevention planning, such simulated intensity distribution results are combined with functions of population distribution and building collapse rates to develop estimates of human and economic damage (seismic risk assessment). Moreover, similar concepts are used in the field of mechanical engineering as foundational data for evaluating the expected frequency components of shaking (not just PGA) in the transport of precision equipment or the base-isolation design within critical facilities.
If you're interested in the formulas used by this tool, the next step is to delve into "why those equations?". The first step is understanding attenuation relations. Equations like $\log_{10}\text{PGA} \approx 0.5M - 1.85 - 1.68\log_{10}R$ used in the tool are empirical formulas derived from statistical processing of numerous seismic records. To learn more deeply, explore the so-called attenuation models; you'll find there are various formulas with different coefficients depending on the region and ground conditions. Second, the relationship between seismic intensity and physical parameters. The Japan Meteorological Agency (JMA) seismic intensity is actually determined by a complex algorithm that considers not just PGA, but also period and duration. Learning the details, such as how the calculation method changes at intensity 4, will give you a clearer understanding of what observed intensity means. Finally, let's return to the starting point of the tool: the seismic moment $M_0$. This is defined as the product of the fault area $S$, the average slip $D$, and the rigidity $\mu$: $M_0 = \mu S D$. Starting from this physical quantity, learning the series of theories for estimating seismic wave radiated energy and fault parameters becomes a genuine first step towards understanding the physics of earthquakes themselves.