Moment magnitude scale

(Redirected from Moment Magnitude Scale)

The moment magnitude scale (MMS; denoted explicitly with M or Mw or Mwg, and generally implied with use of a single M for magnitude[1]) is a measure of an earthquake's magnitude ("size" or strength) based on its seismic moment. Mw was defined in a 1979 paper by Thomas C. Hanks and Hiroo Kanamori. Similar to the local magnitude/Richter scale (ML ) defined by Charles Francis Richter in 1935, it uses a logarithmic scale; small earthquakes have approximately the same magnitudes on both scales. Despite the difference, news media often use the term "Richter scale" when referring to the moment magnitude scale.

Moment magnitude (Mw ) is considered the authoritative magnitude scale for ranking earthquakes by size.[2] It is more directly related to the energy of an earthquake than other scales, and does not saturate – that is, it does not underestimate magnitudes as other scales do in certain conditions.[3] It has become the standard scale used by seismological authorities like the U.S. Geological Survey[4] for reporting large earthquakes (typically M > 4), replacing the local magnitude (ML ) and surface wave magnitude (Ms ) scales. Subtypes of the moment magnitude scale (Mww , etc.) reflect different ways of estimating the seismic moment.

History

edit

Richter scale: the original measure of earthquake magnitude

edit

At the beginning of the twentieth century, very little was known about how earthquakes happen, how seismic waves are generated and propagate through the Earth's crust, and what information they carry about the earthquake rupture process; the first magnitude scales were therefore empirical.[5] The initial step in determining earthquake magnitudes empirically came in 1931 when the Japanese seismologist Kiyoo Wadati showed that the maximum amplitude of an earthquake's seismic waves diminished with distance at a certain rate.[6] Charles F. Richter then worked out how to adjust for epicentral distance (and some other factors) so that the logarithm of the amplitude of the seismograph trace could be used as a measure of "magnitude" that was internally consistent and corresponded roughly with estimates of an earthquake's energy.[7] He established a reference point and the ten-fold (exponential) scaling of each degree of magnitude, and in 1935 published what he called the "magnitude scale", now called the local magnitude scale, labeled ML .[8] (This scale is also known as the Richter scale, but news media sometimes use that term indiscriminately to refer to other similar scales.)

The local magnitude scale was developed on the basis of shallow (~15 km (9 mi) deep), moderate-sized earthquakes at a distance of approximately 100 to 600 km (62 to 373 mi), conditions where the surface waves are predominant. At greater depths, distances, or magnitudes the surface waves are greatly reduced, and the local magnitude scale underestimates the magnitude, a problem called saturation. Additional scales were developed[9] – a surface-wave magnitude scale (Ms) by Beno Gutenberg in 1945,[10] a body-wave magnitude scale (mB) by Gutenberg and Richter in 1956,[11] and a number of variants[12] – to overcome the deficiencies of the ML  scale, but all are subject to saturation. A particular problem was that the Ms  scale (which in the 1970s was the preferred magnitude scale) saturates around Ms 8.0 and therefore underestimates the energy release of "great" earthquakes[13] such as the 1960 Chilean and 1964 Alaskan earthquakes. These had Ms  magnitudes of 8.5 and 8.4 respectively but were notably more powerful than other M 8 earthquakes; their moment magnitudes were closer to 9.6 and 9.3, respectively. [14]

Single couple or double couple

edit

The study of earthquakes is challenging as the source events cannot be observed directly, and it took many years to develop the mathematics for understanding what the seismic waves from an earthquake can tell about the source event. An early step was to determine how different systems of forces might generate seismic waves equivalent to those observed from earthquakes.[15]

The simplest force system is a single force acting on an object. If it has sufficient strength to overcome any resistance it will cause the object to move ("translate"). A pair of forces, acting on the same "line of action" but in opposite directions, will cancel; if they cancel (balance) exactly there will be no net translation, though the object will experience stress, either tension or compression. If the pair of forces are offset, acting along parallel but separate lines of action, the object experiences a rotational force, or torque. In mechanics (the branch of physics concerned with the interactions of forces) this model is called a couple, also simple couple or single couple. If a second couple of equal and opposite magnitude is applied their torques cancel; this is called a double couple.[16] A double couple can be viewed as "equivalent to a pressure and tension acting simultaneously at right angles".[17]

The single couple and double couple models are important in seismology because each can be used to derive how the seismic waves generated by an earthquake event should appear in the "far field" (that is, at distance). Once that relation is understood it can be inverted to use the earthquake's observed seismic waves to determine its other characteristics, including fault geometry and seismic moment.[citation needed]

In 1923 Hiroshi Nakano showed that certain aspects of seismic waves could be explained in terms of a double couple model.[18] This led to a three-decade-long controversy over the best way to model the seismic source: as a single couple, or a double couple.[16] While Japanese seismologists favored the double couple, most seismologists favored the single couple.[19] Although the single couple model had some shortcomings, it seemed more intuitive, and there was a belief – mistaken, as it turned out – that the elastic rebound theory for explaining why earthquakes happen required a single couple model.[20] In principle these models could be distinguished by differences in the radiation patterns of their S-waves, but the quality of the observational data was inadequate for that.[21]

The debate ended when Maruyama (1963), Haskell (1964), and Burridge and Knopoff (1964) showed that if earthquake ruptures are modeled as dislocations the pattern of seismic radiation can always be matched with an equivalent pattern derived from a double couple,[citation needed] but not from a single couple.[22] This was confirmed as better and more plentiful data coming from the World-Wide Standard Seismograph Network (WWSSN) permitted closer analysis of seismic waves. Notably, in 1966 Keiiti Aki showed that the seismic moment of the 1964 Niigata earthquake as calculated from the seismic waves on the basis of a double couple was in reasonable agreement with the seismic moment calculated from the observed physical dislocation.[23]

Dislocation theory

edit

A double couple model suffices to explain an earthquake's far-field pattern of seismic radiation, but tells us very little about the nature of an earthquake's source mechanism or its physical features.[24] While slippage along a fault was theorized as the cause of earthquakes (other theories included movement of magma, or sudden changes of volume due to phase changes[25]), observing this at depth was not possible, and understanding what could be learned about the source mechanism from the seismic waves requires an understanding of the source mechanism.[5]

Modeling the physical process by which an earthquake generates seismic waves required much theoretical development of dislocation theory, first formulated by the Italian Vito Volterra in 1907, with further developments by E. H. Love in 1927.[26] More generally applied to problems of stress in materials,[27] an extension by F. Nabarro in 1951 was recognized by the Russian geophysicist A. V. Vvedenskaya as applicable to earthquake faulting.[28] In a series of papers starting in 1956 she and other colleagues used dislocation theory to determine part of an earthquake's focal mechanism, and to show that a dislocation – a rupture accompanied by slipping – was indeed equivalent to a double couple.[29]

In a pair of papers in 1958, J. A. Steketee worked out how to relate dislocation theory to geophysical features.[30] Numerous other researchers worked out other details,[31] culminating in a general solution in 1964 by Burridge and Knopoff, which established the relationship between double couples and the theory of elastic rebound, and provided the basis for relating an earthquake's physical features to seismic moment.[32]

Seismic moment

edit

Seismic moment – symbol M0  – is a measure of the fault slip and area involved in the earthquake. Its value is the torque of each of the two force couples that form the earthquake's equivalent double-couple.[33] (More precisely, it is the scalar magnitude of the second-order moment tensor that describes the force components of the double-couple.[34]) Seismic moment is measured in units of Newton meters (N·m) or Joules, or (in the older CGS system) dyne-centimeters (dyn-cm).[35]

The first calculation of an earthquake's seismic moment from its seismic waves was by Keiiti Aki for the 1964 Niigata earthquake.[36] He did this two ways. First, he used data from distant stations of the WWSSN to analyze long-period (200 second) seismic waves (wavelength of about 1,000 kilometers) to determine the magnitude of the earthquake's equivalent double couple.[37] Second, he drew upon the work of Burridge and Knopoff on dislocation to determine the amount of slip, the energy released, and the stress drop (essentially how much of the potential energy was released).[38] In particular, he derived an equation that relates an earthquake's seismic moment to its physical parameters:

M0 = μūS

with μ being the rigidity (or resistance to moving) of a fault with a surface area of S over an average dislocation (distance) of . (Modern formulations replace ūS with the equivalent D̄A, known as the "geometric moment" or "potency".[39]) By this equation the moment determined from the double couple of the seismic waves can be related to the moment calculated from knowledge of the surface area of fault slippage and the amount of slip. In the case of the Niigata earthquake the dislocation estimated from the seismic moment reasonably approximated the observed dislocation.[40]

Seismic moment is a measure of the work (more precisely, the torque) that results in inelastic (permanent) displacement or distortion of the Earth's crust.[41] It is related to the total energy released by an earthquake. However, the power or potential destructiveness of an earthquake depends (among other factors) on how much of the total energy is converted into seismic waves.[42] This is typically 10% or less of the total energy, the rest being expended in fracturing rock or overcoming friction (generating heat).[43]

Nonetheless, seismic moment is regarded as the fundamental measure of earthquake size,[44] representing more directly than other parameters the physical size of an earthquake.[45] As early as 1975 it was considered "one of the most reliably determined instrumental earthquake source parameters".[46]

Introduction of an energy-motivated magnitude Mw

edit

Most earthquake magnitude scales suffered from the fact that they only provided a comparison of the amplitude of waves produced at a standard distance and frequency band; it was difficult to relate these magnitudes to a physical property of the earthquake. Gutenberg and Richter suggested that radiated energy Es could be estimated as

 

(in Joules). Unfortunately, the duration of many very large earthquakes was longer than 20 seconds, the period of the surface waves used in the measurement of Ms . This meant that giant earthquakes such as the 1960 Chilean earthquake (M 9.5) were only assigned an Ms 8.2. Caltech seismologist Hiroo Kanamori[47] recognized this deficiency and took the simple but important step of defining a magnitude based on estimates of radiated energy, Mw , where the "w" stood for work (energy):

 

Kanamori recognized that measurement of radiated energy is technically difficult since it involves the integration of wave energy over the entire frequency band. To simplify this calculation, he noted that the lowest frequency parts of the spectrum can often be used to estimate the rest of the spectrum. The lowest frequency asymptote of a seismic spectrum is characterized by the seismic moment, M0 . Using an approximate relation between radiated energy and seismic moment (which assumes stress drop is complete and ignores fracture energy),

 

(where E is in Joules and M0  is in N m), Kanamori approximated Mw  by

 

Moment magnitude scale

edit

The formula above made it much easier to estimate the energy-based magnitude Mw , but it changed the fundamental nature of the scale into a moment magnitude scale. USGS seismologist Thomas C. Hanks noted that Kanamori's Mw  scale was very similar to a relationship between ML  and M0  that was reported by Thatcher & Hanks (1973)

 

Hanks & Kanamori (1979) combined their work to define a new magnitude scale based on estimates of seismic moment

 

where   is defined in newton meters (N·m).

Current use

edit

Moment magnitude is now the most common measure of earthquake size for medium to large earthquake magnitudes,[48][scientific citation needed] but in practice, seismic moment (M0 ), the seismological parameter it is based on, is not measured routinely for smaller quakes. For example, the United States Geological Survey does not use this scale for earthquakes with a magnitude of less than 3.5,[citation needed] which includes the great majority of quakes.

Popular press reports most often deal with significant earthquakes larger than M~ 4. For these events, the preferred magnitude is the moment magnitude Mw , not Richter's local magnitude ML .[49][4]

Definition

edit

The symbol for the moment magnitude scale is Mw , with the subscript "w" meaning mechanical work accomplished. The moment magnitude Mw  is a dimensionless value defined by Hiroo Kanamori[50] as

 

where M0  is the seismic moment in dyne⋅cm (10−7 N⋅m).[51] The constant values in the equation are chosen to achieve consistency with the magnitude values produced by earlier scales, such as the local magnitude and the surface wave magnitude. Thus, a magnitude zero microearthquake has a seismic moment of approximately 1.1×109 N⋅m, while the Great Chilean earthquake of 1960, with an estimated moment magnitude of 9.4–9.6, had a seismic moment between 1.4×1023 N⋅m and 2.8×1023 N⋅m.

Seismic moment magnitude (M wg or Das Magnitude Scale ) and moment magnitude (M w) scales

To understand the magnitude scales based on Mo detailed background of Mwg and Mw scales is given below.

M w scale

Hiroo Kanamori[50] defined a magnitude scale (Log W0 = 1.5 Mw + 11.8, where W0 is the minimum strain energy) for great earthquakes using Gutenberg Richter Eq. (1).

Log Es = 1.5 Ms + 11.8                                                                                     (A)

Hiroo Kanamori[50] used W0 in place of Es (dyn.cm) and consider a constant term (W0/Mo = 5 × 10−5) in Eq. (A) and estimated Ms and denoted as Mw (dyn.cm). The energy Eq. (A) is derived by substituting m = 2.5 + 0.63 M in the energy equation Log E = 5.8 + 2.4 m (Richter 1958), where m is the Gutenberg unified magnitude and M is a least squares approximation to the magnitude determined from surface wave magnitudes. After replacing the ratio of seismic Energy (E) and Seismic Moment (Mo), i.e., E/Mo = 5 × 10−5, into the Gutenberg–Richter energy magnitude Eq. (A), Hanks and Kanamori [51] provided Eq. (B):

Log M0 = 1.5 Ms + 16.1                                                                                   (B)

Note that Eq. (B) was already derived by Hiroo Kanamori[50] and termed it as Mw. Eq. (B) was based on large earthquakes; hence, in order to validate Eq. (B) for intermediate and smaller earthquakes, Hanks and Kanamori (1979) compared this Eq. (B) with Eq. (1) of Percaru and Berckhemer (1978) for the magnitude 5.0 ≤ Ms ≤ 7.5 (Hanks and Kanamori 1979). Note that Eq. (1) of Percaru and Berckhemer (1978) for the magnitude range 5.0 ≤ Ms ≤ 7.5 is not reliable due to the inconsistency of defined magnitude range (moderate to large earthquakes defined as Ms ≤ 7.0 and Ms = 7–7.5) and scarce data in lower magnitude range (≤ 7.0) which rarely represents the global seismicity (e.g., see Figs. 1A, B, 4 and Table 2 of Percaru and Berckhemer 1978). Furthermore, Equation (1) of Percaru and Berckhemer 1978) is only valid for (≤ 7.0).[52]

Relations between seismic moment, potential energy released and radiated energy

edit

Seismic moment is not a direct measure of energy changes during an earthquake. The relations between seismic moment and the energies involved in an earthquake depend on parameters that have large uncertainties and that may vary between earthquakes. Potential energy is stored in the crust in the form of elastic energy due to built-up stress and gravitational energy.[53] During an earthquake, a portion   of this stored energy is transformed into

  • energy dissipated   in frictional weakening and inelastic deformation in rocks by processes such as the creation of cracks
  • heat  
  • radiated seismic energy  

The potential energy drop caused by an earthquake is related approximately to its seismic moment by

 

where   is the average of the absolute shear stresses on the fault before and after the earthquake (e.g., equation 3 of Venkataraman & Kanamori 2004) and   is the average of the shear moduli of the rocks that constitute the fault. Currently, there is no technology to measure absolute stresses at all depths of interest, nor method to estimate it accurately, and   is thus poorly known. It could vary highly from one earthquake to another. Two earthquakes with identical   but different   would have released different  .

The radiated energy caused by an earthquake is approximately related to seismic moment by

 

where   is radiated efficiency and   is the static stress drop, i.e., the difference between shear stresses on the fault before and after the earthquake (e.g., from equation 1 of Venkataraman & Kanamori 2004). These two quantities are far from being constants. For instance,   depends on rupture speed; it is close to 1 for regular earthquakes but much smaller for slower earthquakes such as tsunami earthquakes and slow earthquakes. Two earthquakes with identical   but different   or   would have radiated different  .

Because   and   are fundamentally independent properties of an earthquake source, and since   can now be computed more directly and robustly than in the 1970s, introducing a separate magnitude associated to radiated energy was warranted. Choy and Boatwright defined in 1995 the energy magnitude[54]

 

where   is in J (N·m).

Comparative energy released by two earthquakes

edit

Assuming the values of σ̄/μ are the same for all earthquakes, one can consider Mw  as a measure of the potential energy change ΔW caused by earthquakes. Similarly, if one assumes   is the same for all earthquakes, one can consider Mw  as a measure of the energy Es radiated by earthquakes.

Under these assumptions, the following formula, obtained by solving for M0  the equation defining Mw , allows one to assess the ratio   of energy release (potential or radiated) between two earthquakes of different moment magnitudes,   and  :

 

As with the Richter scale, an increase of one step on the logarithmic scale of moment magnitude corresponds to a 101.5 ≈ 32 times increase in the amount of energy released, and an increase of two steps corresponds to a 103 = 1000 times increase in energy. Thus, an earthquake of Mw  of 7.0 contains 1000 times as much energy as one of 5.0 and about 32 times that of 6.0.

Comparison with TNT equivalents

edit

To make the significance of the magnitude value plausible, the seismic energy released during the earthquake is sometimes compared to the effect of the conventional chemical explosive TNT. The seismic energy   results from the above-mentioned formula according to Gutenberg and Richter to

 

or converted into Hiroshima bombs:

 

For comparison of seismic energy (in joules) with the corresponding explosion energy, a value of 4.2 x 109 joules per ton of TNT applies. The table[55] illustrates the relationship between seismic energy and moment magnitude.

Mw ES
(Joules)
TNT-
equivalency
(tons)
equivalence
Hiroshima-
bomb
(12.5 kT TNT)
3 2.0 · 109 - -
4 6.3 · 1010 000.000.015 00.000.0012
5 2.0 · 1012 000.000.475 00.000.0380
6 6.3 · 1013 000.015,000 00.001.2000
7 2.0 · 1015 000.475,000 00.038,0000
8 6.3 · 1016 015,000,000 01,200,0000
9 2.0 · 1018 475,000,000 38,000,0000
10 6.3 · 1019 15,000,000,000 1,200,000,0000

The end of the scale is at the value 10.6, corresponding to the assumption that at this value the Earth's crust would have to break apart completely.[56]

Subtypes of Mw

edit

Various ways of determining moment magnitude have been developed, and several subtypes of the Mw  scale can be used to indicate the basis used.[57]

  • Mwb – Based on moment tensor inversion of long-period (~10 – 100 s) body-waves.
  • Mwr – From a moment tensor inversion of complete waveforms at regional distances (~ 1,000 miles). Sometimes called RMT.
  • Mwc – Derived from a centroid moment tensor inversion of intermediate- and long-period body- and surface-waves.
  • Mww – Derived from a centroid moment tensor inversion of the W-phase.
  • Mwp (Mi) – Developed by Seiji Tsuboi[58] for quick estimation of the tsunami potential of large near-coastal earthquakes from measurements of the P-waves, and later extended to teleseismic earthquakes in general.[59]
  • Mwpd – A duration-amplitude procedure which takes into account the duration of the rupture, providing a fuller picture of the energy released by longer lasting ("slow") ruptures than seen with Mw .[60]
  • Mdt –Rapidly estimates earthquake magnitude by combining maximum displacements of teleseismic P-wave and source durations.[61]

See also

edit

Notes

edit
  1. ^ These are normally not bolded. In the technical literature a single bolded "M" – with or without italicization – is used for several related concepts.[example needed]
  2. ^ Bormann, Wendt & Di Giacomo 2013, p. 86.
  3. ^ Bormann, Wendt & Di Giacomo 2013, p. 18.
  4. ^ a b The "USGS Earthquake Magnitude Policy" for reporting earthquake magnitudes to the public as formulated by the USGS Earthquake Magnitude Working Group was implemented January 18, 2002, and posted at https://earthquake.usgs.gov/aboutus/docs/020204mag_policy.php. That page was removed following a web redesign; a copy is archived at the Internet Archive.
  5. ^ a b Miyake 2017, p. 112.
  6. ^ Suzuki 2001, p. 121. See also Figure 2-22 in Richter 1958 (copy in Bormann, Wendt & Di Giacomo 2013, p. 60), which replicates Wadati's curves.
  7. ^ Gutenberg & Richter 1956a.
  8. ^ Richter 1935.
  9. ^ See Bormann & Saul 2009 for an overview.
  10. ^ Gutenberg 1945a.
  11. ^ Gutenberg 1945b, Gutenberg & Richter 1956b.
  12. ^ See Seismic magnitude scales.
  13. ^ Kanamori 1977, p. 2981.
  14. ^ ISC-EHB Event 879136 [IRIS]. ISC-EHB Event 869809 [IRIS].
  15. ^ Miyake 2017, pp. 112–113; Stauder 1962, p. 39.
  16. ^ a b Miyake 2017, p. 115.
  17. ^ Ben-Menahem 1995, p. 1210; Maruyama 1963, p. 484.
  18. ^ Ben-Menahem 1995, p. 1210.
  19. ^ Miyake 2017, p. 115. See Byerly 1960 for a contemporary account of why many seismologists favored a single couple model.
  20. ^ Miyake 2017, pp. 116, 117.
  21. ^ Pujol 2003b, p. 164.
  22. ^ Pujol 2003b, p. 165; Miyake 2017, pp. 117–118.
  23. ^ Aki 1966b, p. 84; Pujol 2003b, p. 167.
  24. ^ Julian, Miller & Foulger 1998, §2.2.1.
  25. ^ Miyake 2017, pp. 114, 117; Maruyama 1963, p. 483.
  26. ^ Miyake 2017, p. 117.
  27. ^ Steketee 1958b, pp. 1168–1169.
  28. ^ Stauder 1962, p. 42; Aki & Richards 2002, p. 48.
  29. ^ Honda 1962, pp. 32, 65, and see bibliography; Ben-Menahem 1995, p. 1212; Udías 1991, p. 90; Maruyama 1963, p. 467.
  30. ^ Miyake 2017, p. 467; Steketee 1958a, 1958b.
  31. ^ Udías 1991 provides a partial overview.
  32. ^ Pujol 2003b, pp. 165, 167; Miyake 2017, p. 118.
  33. ^ Bormann, Wendt & Di Giacomo 2013, p. 14.
  34. ^ Aki 1966b, p. 73; Kassaras & Kapetanidis 2018, p. 410.
  35. ^ Beroza & Kanamori 2015, p. 5.
  36. ^ Dziewonski, Chou & Woodhouse 1981, p. 2826; Aki 1966b.
  37. ^ Aki 1966a, pp. 24, 36.
  38. ^ Aki 1966a, p. 24.
  39. ^ Bormann, Wendt & Di Giacomo 2013, p. 12, equation 3.1.
  40. ^ Aki 1966b, p. 84.
  41. ^ Bormann, Wendt & Di Giacomo 2013, p. 14; Bormann & Di Giacomo 2011, p. 412.
  42. ^ Bormann, Wendt & Di Giacomo 2013, pp. 39–40.
  43. ^ Bormann, Wendt & Di Giacomo 2013, p. 7.
  44. ^ Deichmann 2006, p. 1268.
  45. ^ Abe 1982, p. 322.
  46. ^ Kanamori & Anderson 1975, p. 1076.
  47. ^ Kanamori 1977.
  48. ^ Boyle 2008.
  49. ^ Bormann, Wendt & Di Giacomo 2013, p. 86
  50. ^ a b c d Kanamori 1977.
  51. ^ a b Hanks & Kanamori 1979.
  52. ^ Das, Menesis & Urrutia 2023
  53. ^ Kostrov 1974; Dahlen 1977.
  54. ^ Choy & Boatwright 1995
  55. ^ FAQs – Measuring Earthquakes: How much energy is released in an earthquake? United States Geological Survey
  56. ^ Quarks & Co. "Erdbeben – wenn die Erde zurückschlägt" (PDF). Archived (PDF) from the original on 2014-08-26. Retrieved 2022-03-17.
  57. ^ USGS Technical Terms used on Event Pages.
  58. ^ Tsuboi et al. 1995.
  59. ^ Bormann, Wendt & Di Giacomo 2013, §3.2.8.2, p. 135.
  60. ^ Bormann, Wendt & Di Giacomo 2013, §3.2.8.3, pp. 137–128.
  61. ^ Wang, Dun; Kawakatsu, Hitoshi; Zhuang, Jiancang; Mori, Jim; Maeda, Takuto; Tsuruoka, Hiroshi; Zhao, Xu (2017-06-16). "Automated determination of magnitude and source length of large earthquakes using backprojection and P wave amplitudes". Geophysical Research Letters. 44 (11): 5447–5456. Bibcode:2017GeoRL..44.5447W. doi:10.1002/2017GL073801. ISSN 0094-8276.

Sources

edit
edit