Datasheet

DS 3.1

Topic

Magnitude calibration functions and complementary data

Compiled
by
Version

Peter Bormann (formerly GeoforschungsZentrum Potsdam, Telegrafenberg,

D-14473 Potsdam, Germany); E-mail: pb65@gmx.net
May 2001

1 Local magnitude Ml
Table 1 Calibration function L() = - log Ao for local magnitudes Ml according to Richter
(1958). Ao are the trace amplitudes in mm recorded by a Wood-Anderson Standard Torsion
Seismometer from an earthquake of Ml = 0. - epicentral distance in km.
(km)
0
10
20
30
40
50
60
70
80

L()
1.4
1.5
1.7
2.1
2.4
2.6
2.8
2.8
2.9

(km)
90
100
120
140
160
180
200
220
240

L()
3.0
3.0
3.1
3.2
3.3
3.4
3.5
3.65
3.7

(km)
260
280
300
320
340
360
380
400
420

L()
3.8
3.9
4.0
4.1
4.2
4.3
4.4
4.5
4.5

(km)
440
460
480
500
520
540
560
580
600

L()
4.6
4.6
4.7
4.7
4.8
4.8
4.9
4.9
4.9

Table 2 Regional calibration functions L() = - log Ao for Ml determinations. - epicentral

distance and R - hypocentral ("slant") distance with R = (2 + h2), both in km ; h
hypocentral depth in km, T - period in s; Com. - recording component.
Region

L() = - log Ao

Com. Range (km)

Reference

Southern California 1.110 log (R/100) + 0.00189(R - 100) + 3.0 horiz.

10 R 700

Hutton&Boore (1987)

Central California

1.000 log (R/100) + 0.00301(R - 100) + 3.0 horiz.

0 400

Bakun&Joyner (1984)

Great Basin,
Western USA

1.00.log (R/100) + 0.0069(R - 100) + 3.0

0.83 log (R/100) + 0.0026(R - 100) + 3.0

horiz.
horiz.

0 90
90 600

Chvez&Priestley
(1985)

Eastern N-America

1.55 log - 0.22

1.45 log + 0.11

horiz.
vertic.

100 800 Kim (1998)

100 800

Greece

1.58 log (R/100) + 3.0; for ML 3.7

2.00 log (R/100) + 3.0; for ML > 3.7

horiz.

100 800

Kiratzi&Papazachos
(1984)

Albania

1.6627 log + 0.0008 - 0.433

horiz

10 600

Muco&Minga (1991)

Central Europe

0.83 log R + (0.0017/T) (R - 100) + 1.41

vertic.

CentralEurope

1.11 lg R + 0.95 R/1000 + 0.69

vertic.

100 650 Wahlstrm&Strauch

(1984)
10< R < 1000 Stange (2001)

Norway/Fennoskan. 0.91 log R + 0.00087 R + 1.010

vertic.

0 < R 1500

Alsaker et al. (1991)

Tanzania

0.776 log(R/17) + 0.000902 (R - 17) + 2.0

horiz.

0 < R 1000

Langston et al. (1998)

South Australia

1.10 log + 0.0013 + 0.7

vertic.

40 < < 700

Greenhalgh&Singh
(1986)

Datasheet

DS 3.1

2 Teleseismic surface wave magnitude Ms

Table 3 Tabulated magnitude calibration values S () as published in Richter (1958) for MS
determinations according to equation Ms = log AHmax () + S (). AHmax is the (vectorially
combined) maximum horizontal surface-waves displacement amplitude in m for periods
around 20 2 s. Between 20 and 120 the values correspond (rounded to the nearest tenth
magnitude unit), to the values calculated according to the Gutenberg (1945a) relation for the
surface-wave magnitude: Ms = log A + 1.656 log + 1.818. However, for larger distances
this formula yields between 0.05 to 0.55 m.u. larger magnitudes than the tabulated values.
(degrees)
20
25
30
40
45
50

S ()
4.0
4.1
4.3
4.5
4.6
4.6

(degrees)
60
70
80
90
100
110

S ()
4.8
4.9
5.0
5.05
5.1
5.2

(degrees)
120
140
160
170
180

S ()
5.3
5.3
5.35
5.3
5.0

Table 4 Ms magnitude calibration values S () for vectorially combined maximum

horizontal surface-waves displacement amplitude (in m) from shallow earthquakes (h 60
km) as derived by the Prague-Moscow group (Vank et al. (1962)). The IASPEI
recommended standard magnitude formula Ms = log (AH/T)max + 1.66 log + 3.3 fits
these tabulated values between 1 and 140 with deviations < 0.05 magnitude units. For larger
distances the formula overestimates the magnitude between 0.05 and 0.55 (at 180) m.u..
Note: For periods T = 20s the IASPEI standard formula for Ms yields magnitude values that
are 0.18 m.u. larger than those derived from the Gutenberg formula for Ms (see Tab. 3). Since
- in average maximum horizontal and vertical component surface-wave amplitudes agree
well, Tab. 4 and (up to = 160) also the standard Ms calibration function are nowadays
commonly used for Ms determination from vertical component surface-wave records.

0
10
20
30
40
50
60
70
80
90
100
110
120
130
140
150
160
170
180

0
4.96
5.46
5.75
5.96
6.12
6.25
6.36
6.46
6.55
6.62
6.69
6.75
6.79
6.82
6.84
6.84
6.81
6.49

1
3.30
5.03
5.50
5.78
5.98
6.13
6.26
6.37
6.47
6.55
6.63
6.70
6.76
6.79
6.82
6.84
6.84
6.81

2
3.80
5.09
5.53
5.80
5.99
6.14
6.27
6.38
6.48
6.56
6.64
6.70
6.76
6.80
6.82
6.84
6.83
6.80

3
4.09
5.15
5.56
5.82
6.01
6.16
6.28
6.39
6.49
6.57
6.64
6.71
6.76
6.80
6.82
6.84
6.83
6.79

4.30
5.20
5.59
5.84
6.03
6.17
6.30
6.40
6.49
6.58
6.65
6.72
6.77
6.80
6.83
6.84
6.83
6.77

4.46
5.25
5.62
5.86
6.04
6.18
6.31
6.41
6.50
6.58
6.66
6.72
6.77
6.81
6.83
6.84
6.82
6.74

6
4.59
5.29
5.65
5.88
6.06
6.20
6.32
6.42
6.51
6.59
6.66
6.73
6.78
6.81
6.83
6.84
6.82
6.71

7
4.70
5.34
5.68
5.90
6.07
6.21
6.33
6.43
6.52
6.60
6.67
6.74
6.78
6.81
6.83
6.84
6.82
6.69

8
4.80
5.38
5.71
5.92
6.09
6.22
6.34
6.44
6.53
6.61
6.68
6.74
6.78
6.81
6.83
6.84
6.82
6.64

9
4.88
5. 42
5. 73
5. 94
6. 10
6. 24
6. 35
6. 45
6. 54
6. 61
6 .69
6 .75
6 .79
6 .81
6 .83
6 .84
6 .82
6 .59

Surface-wave magnitudes are determined from the maximum amplitude or A/T ratio
measured in the surface-wave train. This is usually the Airy phase of Rayleigh waves (Rmax,
see 2.3 in Chapter 2). It is well developed for shallow earthquakes (depth h < 70 km). Table 5
gives the time difference between the Rmax and the P wave as a function of distance.
2

Datasheet

DS 3.1

Table 5 Time interval (tRmax - tP) between the arrival of the maximum phase of the Rayleigh
wave and the first onset of P waves as a function of according to Archangelskaya (1959)
and Gorbunova and Kondorskaya (1977) (From Willmore, 1979).

10
15
20
25
30
35
40
45
50

tRmax - tP
(min)
4-5
6-8
9-10
10-12
13-14
15-16
18-19
21
24

tRmax - tP
(min)
26
28-29
31
33
35
37
39-40
42
43

55
60
65
70
75
80
85
90
95

100
105
110
115
120
125
130
140
150

tRmax - tP
(min)
45-46
47-48
48-50
53
55
57
60
64
70

3 Teleseismic body-wave magnitudes mB

Gutenberg (1945) developed a magnitude relationship for teleseismic body waves such as P,
PP and S in the period range 0.5 s to 12 s (i.e., mostly based on medium-period readings):
mB = log (A/T)max + Q(, h).
Gutenberg and Richter (1956a) published a table with Q() values for P, PP and S waves in
vertical (V=Z) and horizontal (H) components for shallow events (Table 6) as well as
diagrams giving for all these waves Q values as a function of and source depth h (Figs. 1ac). These Q values are valid only when A is given in m.

Figure 1a

Datasheet

DS 3.1

Figure 1b

Figure 1c

Datasheet

DS 3.1

Table 6 Values of Q() for P, PP and S waves for shallow shocks (h < 70 km) according to
Gutenberg and Richter (1956a) if the ground amplitude is given in m.
PV PH PPV PPH SH PV PH PPV PPH SH

PV PH PPV PPH SH

16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55

7.3
7.4
7.5
7.5
7.4
7.3
7.4
7.5
7.6
7.7
7.8
7.9
7.9
8.0
8.1
8.2
8.6
8.8
9.0

5.9
5.9
5.9
6.0
6.0
6.1
6.2
6.3
6.3
6.5
6.4
6.5
6.6
6.6
6.6
6.7
6.7
6.7
6.7
6.7
6.6
6.5
6.5
6.4
6.4
6.5
6.5
6.5
6.5
6.7
6.8
6.9
6.9
6.8
6.7
6.7
6.7
6.7
6.8
6.8

6.0
6.0
6.0
6.1
6.1
6.2
6.3
6.4
6.5
6.6
6.6
6.7
6.7
6.7
6.8
6.9
6.9
6.9
6.9
6.9
6.9
6.7
6.7
6.6
6.6
6.7
6.7
6.7
6.7
6.9
7.1
7.2
7.2
7.1
7.0
7.0
7.0
7.0
7.1
7.1

6.7
6.7
6.8
6.8
6.8
6.8
6.7
6.7
6.7
6.6
6.6
6.5
6.5
6.6
6.7
6.7
6.7
6.7
6.7
6.7
6.7
6.7
6.7
6.7
6.8
6.9

6.8
6.8
6.9
6.9
6.9
6.9
6.8
6.8
6.8
6.7
6.7
6.6
6.6
6.7
6.8
6.8
6.8
6.8
6.8
6.8
6.8
6.8
6.8
6.8
6.9
7.0

7.2
6.8
6.2
5.8
5.8
6.0
6.2
6.2
6.2
6.2
6.2
6.3
6.3
6.3
6.3
6.3
6.4
6.4
6.5
6.6
6.6
6.6
6.6
6.7
6.7
6.6
6.5
6.5
6.5
6.5
6.6
6.6
6.7
6.7
6.6
6.5
6.5
6.6
6.6
6.6

56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95

6.8
6.8
6.8
6.8
6.8
6.9
7.0
6.9
7.0
7.0
7.0
7.0
7.0
7.0
6.9
6.9
6.9
6.9
6.8
6.8
6.9
6.9
6.9
6.8
6.7
6.8
6.9
7.0
7.0
7.0
6.9
7.0
7.1
7.0
7.0
7.1
7.1
7.2
7.1
7.2

7.1
7.1
7.1
7.1
7.1
7.2
7.3
7.3
7.3
7.4
7.4
7.4
7.4
7.4
7.3
7.3
7.3
7.2
7.1
7.1
7.2
7.2
7.3
7.2
7.1
7.2
7.2
7.4
7.4
7.4
7.3
7.3
7.5
7.4
7.3
7.5
7.4
7.5
7.4
7.6

6.9
6.9
7.0
7.0
7.1
7.2
7.3
7.3
7.3
7.3
7.3
7.2
7.1
7.0
7.0
7.1
7.1
7.1
7.0
6.9
6.9
6.9
6.9
6.9
6.9
7.0
7.1
7.2
7.3
7.3
7.3
7.2
7.2
7.2
7.2
7.2
7.2
7.2
7.2
7.2

7.0
7.0
7.1
7.2
7.3
7.4
7.4
7.4
7.5
7.5
7.4
7.4
7.3
7.2
7.2
7.3
7.3
7.3
7.2
7.1
7.1
7.1
7.1
7.1
7.1
7.2
7.3
7.4
7.5
7.5
7.5
7.4
7.4
7.4
7.4
7.4
7.4
7.4
7.4
7.4

6.6
6.6
6.6
6.6
6.6
6.7
6.7
6.7
6.8
6.9
6.9
6.9
6.9
6.9
6.9
7.0
7.0
6.9
6.8
6.8
6.8
6.8
6.9
6.8
6.7
6.8
6.9
6.9
6.9
6.8
6.7
6.8
6.8
6.8
6.8
6.9
6.9
6.9
7.0
7.0

96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
112
114
116
118
120
122
124
126
128
130
132
134
136
138
140
142
144
146
148
150
152
154
156
158
160
170

7.6
7.8
7.8
7.8
7.7
7.6
7.7
7.9
7.9
8.1
8.2
8.3
8.3
8.4
8.5
8.6
9.0

7.2
7.2
7.2
7.2
7.2
7.2
7.2
7.2
7.3
7.3
7.4
7.4
7.4
7.4
7.4
7.4
7.5
7.5
7.5
7.5
7.4
7.3
7.2
7.1
7.0
7.0
6.9
6.9
7.0
7.1
7.1
7.0
6.9
6.9
6.9
6.9
6.9
6.9
6.9
6.9
6.9

7.4
7.4
7.4
7.4
7.4
7.4
7.4
7.4
7.5
7.5
7.6
7.6
7.6
7.6
7.6
7.6
7.7
7.7
7.7
7.7
7.6
7.5
7.4
7.4
7.3
7.3
7.2
7.2
7.3
7.4
7.4
7.3
7.2
7.2
7.2
7.2
7.2
7.2
7.2
7.2
7.2

7.1
7.2
7.3
7.3
7.4
7.4
7.4
7.3
7.3
7.2
7.2
7.2
7.2
7.2
7.2

At the IASPEI General Assembly in Zrich (1967) the Committee on Magnitudes

recommended stations to report the magnitude for all waves for which calibration functions
are available, as well as to publish amplitude and period values separately. Q(
, h)PZ is now
the accepted standard calibration function for mb magnitude determinations at
international data centers based on short-period vertical component P-wave readings. This is
not fully correct because the Q values have been derived mainly from intermediate-period
seismic recordings.

Datasheet

DS 3.1

4 Complementary short-period body-wave magnitude scales

Another calibration function P(, h) for mb determination has been elaborated by Veith and
Clawson (1972). It is based on large sets of short-period vertical-component P-wave
amplitudes from large explosions at 19 different sites. Although specifically derived from
short-period data it is not yet accepted as IASPEI standard for mb. It looks much smoother
than Q(, h)PZ and resembles better an inverse A- relationship for short-period P as shown in
Fig. 3.13. It is currently used by the preliminary International Data Center established for the
monitoring of the Comprehensive Test-Ban Treaty (CTBT) however with a non-standard
instrument response.

Figure 2 Calibration functions P(, h) for mb determination from narrow-band short-period

vertical-component records with peak displacement magnification around 1 Hz (WWSSN-SP
characteristic) according to Veith and Clawson (1972). Note: P values have to be used in
conjunction with maximum P-wave peak-to-trough (2A!) amplitudes in units of nanometers
(1 nm = 10-9m) (modified from Veith and Clawson, Magnitude from short-period P-wave
data, BSSA, 62, 2, p. 446, Seismological Society of America).
An experimental calibration function for magnitude determinations based on short-period
vertical-component readings of various PKP phases in the distance range 145 to 164 has
been developed by Wendt (Bormann and Wendt, 1999). The following relationship is used:
mb(PKP) = log10 (A/T) + Q(, h)PKPab,bc,df
with amplitude A in m (10-6 m) (see Figure 3). Extensive use of this relationship at station
CLL proved that mb determinations from core phases are possible with a standard deviation
of less than 0.2 magnitude units as compared to P-wave mb determinations by NEIC and
ISC. If more than one PKP phase can be identified and A and T been measured then the
average value from all individual magnitude determinations provides a more stable estimate.
The applicability of these calibration functions should be tested with data from other stations
of the world-wide network.
6

Datasheet

DS 3.1

Figure 3 Calibration functions according to S. Wendt for the determination of mb(PKP) for
PKPdf, PKPbc and PKPab (see Bormann and Wendt, 1999).
7

Datasheet

DS 3.1

References (see References under Miscellaneous in Volume 2)