7.1.3 Magnitudes

The phase readings for magnitude determination is described on page [*].

In SEISAN version 8.3, there are substantial changes in the way amplitudes are read and two new magnitude scales have been added (broad band body and surface wave magnitudes). Furthermore, the Richter attenuation curve is now used be default for the body wave magnitude. The phase names used for amplitudes have also changed. These changes are due to the new standards for magnitude calculation approved by the IASPEI. For more on the application of the different magnitude scales, see Havskov and Ottemöller [2010].

Magnitudes are calculated using coda, amplitude and spectral level. Parameters are given in the station file using the RESET TEST variables. For magnitude based on amplitude, the amplitude must be given in nanometers or nm/s in the input file (SEISAN standard).

Local magnitude Ml

The formula used to calculate local magnitude is

\bgroup\color{black}$\displaystyle Ml = a * log_{10} (amp) + b*log_{10}(dist) + c*dist + d
$\egroup

where a,b,c,d are constants, \bgroup\color{black}$ log_{10}$\egroup is logarithm to the base 10, amp is maximum ground amplitude (zero \bgroup\color{black}$ -$\egrouppeak) in nm and dist is hypocentral distance in km (RESET TEST 75-78). The default constants used, if not set in STATION0.HYP, are for California [Hutton and Boore, 1987] which gives the following relation

\bgroup\color{black}$\displaystyle Ml = log_{10}(amp) + 1.11 log_{10}(dist) + 0.00189 dist - 2.09
$\egroup

It is here assumed that the gain of the Wood-Anderson instrument is \bgroup\color{black}$ 2080$\egroup. An amplitude of \bgroup\color{black}$ 1 mm$\egroup of the Wood Anderson seismogram is then \bgroup\color{black}$ 10^6nm/2080$\egroup and inserting this amplitude above together with a distance of 100 km gives magnitude 3 as originally defined by Richter. If using the default STATION0.HYP, distributed with SEISAN, the Ml parameters are for Norway. It is assumed that the maximum amplitude is picked on a seismogram simulating the original Wood-Anderson seismogram, see program MULPLT. SEISAN uses hypocentral distance, while the original Ml scale used epicentral distance (no deep earthquakes in California). We use hypocentral distance so Ml also can be used for deep earthquakes, but the user should be aware that the Ml relation for deep earthquakes might be different from the relation for shallow earthquakes.

Local magnitudes are only calculated for events with epicentral distance LESS THAN TEST(57) (default 1500 km) and if the period is less than 5.0 secs. All amplitudes for the phases `L', `S ', Sg, SG, AMP, and AML, AMPL, IAML or blank are used. This means that if an amplitude is picked on both Lg and Sg, both will be used. The period is not used. The many possible phase names is a result of changes over time and thus to ensure that Ml is calculated correctly with older data. From version 8.3, MULPLT produces the standard IASPEI name IAML.

At very short distance, the traditional Ml scale may not be valid and it is possible to include additional terms in the Ml calculation [Luckett et al., 2018].

\bgroup\color{black}$\displaystyle Ml = a * log_{10} (amp) + b*log_{10}(dist) + c*dist + e*exp(-f*dist) + d
$\egroup

The values for e and f are given by TEST(xx) and TEST(yy), respectively.

Local magnitude Mw;

The mb(Pn/Sn) scale has recently been defined for the North Atlantic region [Kim and Ottemöller, 2017; Kim et al., 2020]. This scale is based on Pn and Sn amplitudes measured, just like ML amplitudes, on the Wood-Anderson simulated trace. The reason for introducing this scale was that in many regions ML is defined to correct for attenuation in continental crust where Sg/Lg are typically the phases with largest amplitudes. However, these are not propagating in Oceanic crust and amplitudes need to be measured from Pn or Sn phases. As the scale is calibrated against GCMT MW, we label it Mw in SEISAN

The Mw scale in SEISAN is implemented as described in Kim and Ottemöller [2017], given by

\bgroup\color{black}$\displaystyle MN = log_{10}(amp) - a log_{10}(100/dist) + S + E + b
$\egroup

where the amplitude is in units of nm and measured on Pn or Sn on the WA simulated trace, dist is in km, S and E are the station and event (or source region) correction terms. For the North Atlantic, a = 1.86 and b = 1.62.

The scale is activated with RESET TEST(117).

The scale is defined in mbn.par, an example based on Kim and Ottemöller [2017]; Kim et al. [2020] is:

# 
# definition of scale for Pn and Sn
# MN = log10(A) - a log10(ref/dist) + S + E + b
#                                       a         ref       b
SCALE MBPN                              +1.861    100.      1.618
SCALE MBSN                              +2.176    100.      1.426

#
# group velocity windows for Pn and Sn
#
GROUP VEL PN                            5.0       9.0
GROUP VEL SN                            2.0       5.0

#
# station corrections (S)
#                                       Pn        Sn
STATION CORRECTION  ARE0                 0.168     0.168
STATION CORRECTION  BJO1                -0.217    -0.088
STATION CORRECTION  DAG                 -0.285    -0.166
STATION CORRECTION  DBG                  0.101    -0.052
STATION CORRECTION  HAMF                 0.143     0.094
STATION CORRECTION  HOPEN               -0.348    -0.332
STATION CORRECTION  HSPB                 0.047     0.086
STATION CORRECTION  JMIC                -0.308     0.143
STATION CORRECTION  KBS                 -0.072     0.001
STATION CORRECTION  KEV                  0.191     0.133
STATION CORRECTION  KTK1                 0.072     0.155
STATION CORRECTION  LOF                  0.153    -0.140
STATION CORRECTION  MOR8                 0.065    -0.137
STATION CORRECTION  NOR                  0.120     0.145
STATION CORRECTION  NSS                 -0.046    -0.298
STATION CORRECTION  SCO                  0.331     0.343
STATION CORRECTION  SPA0                -0.259    -0.086
STATION CORRECTION  STEI                 0.023     0.002
STATION CORRECTION  TRO                  0.120     0.028

#
# event corrections E
#                   name of area         minlon maxlin minlat maxlat Pn    Sn
SOURCE CORRECTION   Gakkel Ridge         -10.0   10.0   82.0   85.0 -0.090-0.082
SOURCE CORRECTION   Spitsbergen  TF      -10.0   10.0   79.7   82.0 -0.194-0.352
SOURCE CORRECTION   Molloy FZ              0.0   10.0   77.0   79.7 -0.047-0.118
SOURCE CORRECTION   Knipovich Ridge        0.0   11.0   73.0   77.0 +0.215+0.215
SOURCE CORRECTION   Mohns Ridge           -4.2    8.0   71.0   73.0 +0.043+0.277
SOURCE CORRECTION   Jan Mayen TF         -14.0   -4.2   70.0   72.0 -0.260-0.077
SOURCE CORRECTION   Kolbeinsey Ridge     -20.0  -14.0   69.3   71.0 +0.274+0.379

Coda magnitude Mc

The coda magnitude is calculated using

\bgroup\color{black}$\displaystyle Mc = a * log_{10}(coda) + b * dist + c
$\egroup

where coda is coda length in secs and a,b and are constants (RESET TEST 7-9). If `a' is given as a negative number, the following formula will be used

\bgroup\color{black}$\displaystyle Mc = abs(a)*log_{10}(coda)*log_{10}(coda) + b * dist + c
$\egroup

If both Mc and Ml are calculated, Ml is written first on the header line.

Coda magnitude is only calculated if the epicentral distance is less than TEST(57).

Surface wave magnitude Ms

Ms is calculated using the standard

\bgroup\color{black}$\displaystyle Ms = log_{10}(amp/T)+1.66log_{10}(dist)+3.3
$\egroup

where T is period. Amplitude is in micrometer and distance in degrees, however in the Nordic format nm and km are used and the program converts. Ms is only calculated if the period is larger than 10.0 seconds in which case the program automatically assumes that Ms is the wanted magnitude. The phase used can be AMS, AMPS, AMP or blank. The current version of MULPLT produces the standard IASPEI name IAMs_20. The many possible phase names are a result of changes over time and thus to ensure that Ms is calculated correctly with older data. It is assumed that the amplitude has been picked on a WWSSN standard LP trace and that the period is in the range \bgroup\color{black}$ 18-22 s$\egroup (see program MULPLT). Ms will be calculated even if the period is outside this range, but it will not be correct according to the standard. The distance range is between TEST(114) and 100 deg. before version 10.5, there were no distance limits. TEST(114)=20 deg by default. Depth must be less than TEST(115). Before version 10.5, there were no depth limit. TEST(115) is 60 km by default.

Broadband surface wave magnitude MS (IASPEI code MS_BB, but SEISAN uses MS for simplicity, new from SEISAN version 8.3)

MS is calculated using the standard

\bgroup\color{black}$\displaystyle MS = log_{10}(amp/T)_{max}+1.66log_{10}(dist)+3.3
$\egroup

or

\bgroup\color{black}$\displaystyle MS = log_{10}(V_{max}/2 \pi )+1.66log_{10}(dist)+3.3
$\egroup

where \bgroup\color{black}$ V_{max}$\egroup is the maximum velocity. The IASPEI definition is to use velocity and the period is thus not needed but read for information. The velocity is in micrometer/s and distance in degrees, however in the Nordic format \bgroup\color{black}$ nm/s$\egroup and \bgroup\color{black}$ km$\egroup are used and the program converts when calculating magnitudes. MS is only calculated if the period is larger than 3 seconds and less then 60 seconds, distance must be larger than or equal to 222 km (2 degrees) and less or equal to 160 degrees. The depth must be less than TEST(115) (default 60 km). Before version 10.5, there was no check of depth. The phase used to report the amplitude and period must be called IVMs_BB which the current version of MULPLT produces. The biggest advantage using MS compared to Ms, is that any period in the range \bgroup\color{black}$ 2-60 s$\egroup can be used.

Body wave magnitude mb

mb is calculated using

\bgroup\color{black}$\displaystyle mb = log_{10}(amp/T) + Q(dist,depth)
$\egroup

where Q is a hardwired function of distance and depth and amp is the amplitude in nm. There are two possibilities: The default (set by REST TEST(108) is the standard Gutenberg and Richter (1956) curve while alternatively the Veith-Clawson curve can be used [Veith and Clawson, 1972]. Before SEISAN version 8.3, Veith-Clawson was always used. mb is only calculated if the epicentral distance is less than or equal to 100 degrees and larger than or equal to TEST(113), default 20 deg (IASPEI standard and SEISAN default is 20 degrees).Before version 10.5, the lower distance was TEST(57) (default 1500 km). The period must be smaller than 3s and larger then 0.2 s and the phase is P, AMP, AMb, AMB, AMPB, AMPb, blank character or IAmb. The current version of MULPLT produces the standard IASPEI name IAmb. The many possible phase names are a result of changes over time and thus to ensure that mb is calculated correctly with older data.

Broad band body wave magnitude mB (new from SEISAN version 8.3)

The broad band magnitude mB (official IASPEI name is mB_BB) is calculated using

\bgroup\color{black}$\displaystyle mB = log_{10}(amp/T)_{max} + Q(dist,depth)
$\egroup

or

\bgroup\color{black}$\displaystyle mB = log_{10}(V_{max}/2 \pi ) + Q(dist,depth)
$\egroup

where \bgroup\color{black}$ V_{max}$\egroup is the maximum velocity and Q is a hardwired function of distance and depth. The IASPEI standard is to use velocity and SEISAN store the velocity in nm/s. There are two possibilities for the atteneuation function: The default (set by RESET TEST(108) is the standard Gutenberg and Richter (1956) curve while alternatively the Veith-Clawson curve can be used [Veith and Clawson, 1972]. mB is only calculated if the epicentral distance is less than or equal to 100 degrees and larger than or equal to TEST(113) (IASPEI standard and SEISAN default 2 degrees. Before version 10.5, the lower limit was TEST(57), default 1500 km) and the period is larger than 0.2s and less than 30s and the phase name is IVmB_BB . The current version of MULPLT produces the standard IASPEI name IVmB_BB. The biggest advantage using mB compared to mb, is that the mB scale does not saturate before magnitude 8.

Moment magnitude Mw

Mw is calculated as

\bgroup\color{black}$\displaystyle Mw = 2/3 * \Big( log_{10}(moment) - 9.1 \Big)
$\egroup

where moment is in Nm (see also section 8.12). When an event is relocated, the moment is also recalculated according to revised hypocentral distance. Mw is only calculated from moment as given on the SPEC lines, see MULPLT.

NOTE: If an amplitude has a given period between 5 and 10 secs, it is not used for Ml and mb magnitude calculation, see above. If an event is not located, there will normally be no magnitude calculation and all magnitude and distance information is deleted from the output S-file (hyp.out) except, the magnitude in the 3rd position on the header line if it has an agency different from the default agency. The only exception is that if a coda is given, the epicentral distance is retained and coda magnitude will therefore be calculated. This means that for events, which cannot be located, it is still possible to calculate coda magnitudes by manually entering the epicentral distance on the line containing the coda length.

On the first header line, there is room for 3 magnitudes. If there is a magnitude in the 3rd position, it is not overwritten unless the default agency is overwritten, so there will often only be room for 2 calculated magnitudes on the first header line. If more magnitudes are calculated, they will be written on a subsequent hypocenter line, which is identified by having the same year, month, day and hypocenter agency as the first header line. This means that there is room for a total of 6 magnitudes, which can each, be updated when relocating. Hypocenter info and all 6 magnitudes can be printed out on one line with program REPORT.

All magnitudes can have a station dependent correction given in the station file. This correction does not affect the Mc in print.out file. Mb and mB use the same correction and Ms and MS use the same correction. Only calculate magnitude: If TEST(106) is set to 1.0, only magnitudes are calculated, provided a distance is given.

NOTE*********

Amplitude phases and coda length (END-phase) can be given weight 4 and the corresponding magnitude is not calcualted. Other weights cannot be used.

Peter Voss : Tue Jun 8 13:38:42 UTC 2021