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

where a,b,c,d are constants, is logarithm to the base 10, amp is maximum ground amplitude (zero peak) 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

It is here assumed that the gain of the Wood-Anderson instrument is . An amplitude of of the Wood Anderson seismogram is then 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].

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

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

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

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).

Ms is calculated using the standard

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 (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

or

where 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 and 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 can be used.

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

or

where 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.

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 : Fri Nov 12 10:33:10 UTC 2021