The foregoing section has described the appearance of the records at various distances from the focus of an earthquake. The seismologist who first looks at a record, however, may have no knowledge of the source of the disturbance and must therefore use the record itself to find out what he can.
If the station has a full set of long-period and short-period instruments, look for the first observable motion on all the records. The normal situation is that P or PKP is strongest on the short-period vertical instrument. If the same phase is readable on the horizontal records, the approximate azimuth can be determined from the fact that the P motion is along the ray, either towards or away from the source. Hence note the following:
The apparent azimuth of the horizontal motion is estimated by measuring the amplitudes of corresponding longitudinal waves on the N-S and E-W components, and determining the azimuth from the equation
The resulting value of tan alpha corresponds to two possible azimuths differing by 180°. The decision between these depends on the direction of vertical motion, as summarised in the following table:
The polarisation of surface waves on long-period instruments may also be used to estimate the azimuth of an earthquake. If the first clearly-recorded surface waves are on the horizontal instruments only, they will be Love waves, and their polarisation will be perpendicular to the direction of approach. Thus, if corresponding peaks are to the north and west, the Love waves will be approaching from the northeast, or southwest. This ambiguity is resolved by the Rayleigh waves, which are also recorded on the vertical component. The horizontal components will be in phase along the line of approach of the waves, but will be 90° out of phase relative to the vertical components. The direction from which the waves approach is the direction of the peak on the horizontal record which follows one quarter of a cycle after the upward peak on the vertical. Thus, if the order of successive peaks is U, N, D, S and U, E, D, W, the earthquake is in the north-eastern quadrant from the station.
The angle of incidence can also, in principle, be determined by measuring the vertical component of amplitude, and using the following formula:
In practice, these estimates may be severely disturbed by local crustal irregularities, and should therefore be interpreted with caution.
If the first P phase is weak or absent on the short-period instruments, the source is probably in the shadow zone, and the first motion seen on the long-period instruments is a diffracted P.
When the P arrival has been identified and read, look for S, which is usually the first phase which is more conspicuous on the long-period horizontal instruments than on the short-period vertical. Usually the general character of the seismic signal will suggest whether we have a normal-focus or a deep-focus event. Large shallow earthquakes have extended complex body-wave signals and excite large-amplitude surface waves. This is due to near-source multiple reflections of the seismic signal and perhaps to multiple ruptures within the main shock. Large deep earthquakes have compact body-wave signals and are poor exciters of surface waves.
Now read the times of all other phases, make a list of them and subtract the time of the first arrival. Then take a piece of paper, and put it down against the time axis of the travel-time curve. Make a mark at the origin to indicate the time of the first onset, and other marks to represent the times of the other phases on the scale of the graph. Move the piece of paper over the graph, keeping the first onset mark on the P or PKP curve, until as many of the other phases as possible fit on to the lines. If no fit can be found on the normal-focus curves, try the deep-focus ones, and interpolate mentally to find the depth as accurately as possible. Then enter the phase identifications, and check the diagnostic characteristics which have been given in the section on Record Content. Epicentral distance and time of origin can be read immediately from the curve. The process is illustrated by Fig. 1.1.2.
A very practical device for preliminary identification of phases is a set of traveltime curves plotted on transparent paper or on plexiglass plates; the time scale must be selected according to the recording speed of the seismograph producing the most complete records, usually classes B, C, D, E. Again by moving the graph over the seismogram we attempt to find the best coincidence between recorded phases and travel-time curves.
Although the general appearance of the record gives some indication of focal depth, a more precise estimate is often desirable. Shallow earthquakes are far more numerous than deep earthquakes. Shallow events are more likely to produce damage or to generate tsunamis. Also, good estimates of depth of focus are of considerable geophysical importance.
Certain phase arrivals may correspond to the near-source surface reflections, pP and sP, which are the most important 'depth phases'. Figs. 1.1.3a and 1.1.3b show the time differences pP-P and sP-P as a function of observer distance and source depth. This figure is based on the tables of Jeffreys and Shimshoni (1964) who used the same data as in the preparation of the J-B Tables (Jeffreys and Bullen, 1940). Note that the time difference pP-P and sP-P are not simply related as a function of observer distance and focal depth.
The pP phase, having been reflected from the free surface, will have opposite polarity from that of the P arrival, provided that both emerge from the same focal mechanism quadrant. Depth can be unambiguously determined when both the pP and sP arrivals are present, provided that multiple rupture arrivals and suchlike are not prominent. However in the case of shallow earthquakes, radiation pattern factors often lead to the strong excitation of only one of these depth phases. Experience at a given seismic station will often allow a good guess at this phase identification, depending on the source region. For example it is nearly always true that Chilean earthquakes will show a strong sP phase and a very weak pP phase at eastern United States stations. In the distance range 70° to 80°, PcP-P has values similar to pP-P and sP-P for shallow earthquakes and care should be taken to avoid confusion in phase identification. Generally, even if positive identification of a potential depth phase cannot be made, it is important to report such a reading to central seismological agencies, where additional data may permit positive phase identification.
If the pP reflection point is oceanic, the strongest reflector will be the ocean-air interface. The subsequent pP-P times will thus reflect significant water travel path and lead to estimates of depth that can be too great by 10 to 20 km. The corresponding sP phase will be reflected at the ocean bottom. Often weak reflections of the pP type will occur at acoustical boundaries in the source-region crust. These phases will look like pP but have much lower amplitude.
Less important depth phases include PcP, pPcP, sPcP, pPKP, sPKP, (SKP-PKP), and (PKS-SKP). These can be used with standard seismological tables (Jeffreys and Bullen, 1940; Herrin and others, 1968) for estimates of focal depth.
If the coordinates of the earthquake focus have been determined by an epicentral agency, the station operator will wish to determine azimuth and distance as a basis for final interpretation. The formulae used for this purpose are based on the equations of spherical trigonometry. In a sphere, the radius vector at any point on the surface coincides with the normal to the surface at the same point, so that the directions of either of these vectors may be used in the calculations. The real Earth, however, has a noticeable degree of ellipticity. 'Geographic' latitudes are determined by measuring the angle between the vertical direction of any point and the equatorial plane, and this can differ significantly from the 'Geocentric' latitude which is the angle between the radius vector and the equatorial plane.
The use of geocentric latitude substantially simplifies the task of calculating distances over the surface of an ellipsoidal earth, or of comparing travel times along ray paths which penetrate into it. Conversions from geographic to geocentric latitude may be calculated from the formula
Or may be read from Fig. 1.2.1. The geocentric latitudes so obtained will be appropriate for all subsequent charts and formulae.
The angular distance between a station (phi0, lambda0) and a known epicentre (phi, lambda) can be determined numerically by means of the following basic relation of spherical trigonometry:
In determining the azimuth from a station to an epicentre measured from north to east we use
To obtain azimuth from epicentre to station, interchange phi0, lambda0 with phi, lambda.
The formulae can be simplified for a numerical procedure by introducing the direction cosines. If we take rectangular axes at the centre of the Earth, the axes x and y being directed to the points on the equator with longitudes 0° and 90° respectively, and the axis z being directed positively to the north pole, the direction cosines of the radius to a place with geographic co-ordinates phi0, lambda0 are given by
where phi' is the geocentric latitude. If A, B, and C are analogous quantities for the epicentre (phi, lambda) we can write
Alternative formulae are
The first formula is recommended for Delta = 60°-120°, the second for Delta = 0°-60°, and the third one for Delta = 120°-180°. The values of a, b, c were tabulated by Comrie and Jeffreys (1938) and periodically by the International Seismological Centre (from 1967).
An approximate formula
can be used when Delta<6.5°. Angular distances can be converted to kilometres by converting degrees to radians, and multiplying by the local radius r of the earth, where
and h is the value from Table PAR 2 at the latitude corresponding to the centre of the ray path. As an alternative, N-S and E-W distances may be taken directly from geodetic tables, and the classical approximation (1° = 111.1 km) is good to within 0.1% over most of the surface of the Earth.
The chart which is most generally applicable for approximate determinations of azimuth and distance is the 'stereographic net' which is widely used for the solution of spherical triangles in crystallography.
For use in seismology, the meridians should be numbered in both directions across the equator. In referring to points in the eastern hemisphere, the zero meridian is taken to be the left-hand half of the bounding circle of the projection, and the longitudes are the positive numbers reading from left to right. For the western hemisphere, the zero meridian is the right-hand boundary, and the longitudes are the negative numbers reading from right to left. The latitude corresponding to each parallel is entered on the bounding circle, and the angular distances of the parallels from the south pole are marked along the vertical axis. When in use, the chart is fixed to a smooth board, and covered by a sheet of tracing paper held in place by a single pin through the centre.
To find the distance and azimuth of a point P (whose latitude and longitude are phi0, lambda0) from Q (phi, lambda) first subtract the angle lambda from the longitudes of both points. Plot P on the tracing paper at phi0, (lambda0 - lambda) and Q at phi, 0, using the sign convention for longitude outlined in the preceding paragraph. (Fig. 1.2.3). Now rotate the tracing paper until Q is carried to the south pole of the projection, and let P1 and Q1 be the displaced positions of P and Q. Then the distance PQ being equal to P1Q1 can be read immediately on the vertical scale of the projection. The azimuth of P is equal to the longitude of P1, and is read on the equatorial scale of the projection.
If it is sufficient to refer all azimuths and distances to a single station, the stereographic net may be used more conveniently by drawing the azimuth and distance circles on the base chart (Willmore, 1957). Assuming that the complete circle of Fig 1.2.4 represents the bounding circle of a stereographic net, the construction of the distance circles proceeds as follows:
In carrying out the construction, difficulty may be encountered when CG or CH becomes almost parallel with AB, for under these conditions the intersection point moves off towards infinity. To avoid this difficulty we use the relations
where theta = angle OCR. These relations can easily be proved by considering the magnitudes of the angles OCH and OCI. By yielding a value of the diameter OH-OI, the equations enable the larger circles to be constructed by using a bent spline, and therefore eliminate the requirement for a long beam compass.
The azimuth circles, being great circles through R, are orthogonal to all the small circles about R. Hence the small circle which projects as a straight line is the common diameter of the azimuth circles. Let this line cut AB at H'. As the line can be regarded as a circle of infinite diameter, we find theta from equation  by setting (theta + Delta/2) = 90°, and hence find OH from [ . The common diameter of the azimuth circles is the normal to AB through H' .
The N-S azimuth circle is the meridian through R, and its centre O'' is the point at which the diameter through H' cuts the equator. As the stereographic projection conserves angles, the great circle which corresponds to an azimuth Z will make an angle Z with the meridian NRS. Its tangent at R is therefore constructed by laying off an angle Z from the tangent of NRS, and its diameter is perpendicular to this tangent. The centre of the E-W azimuth circle is therefore the point O'' at which the diameter through R intersects the diameter through H'. For interpreting seismograms it is often sufficient to draw only the E-W azimuth circle in addition to the meridian. This circle is shown as WRE in Fig. 1.2.4. Base nets, 20 inches in diameter, may be obtained from the Dominion Observatory, Ottawa, Canada. A typical finished chart is shown in Fig. 1.2.4a.
To find the azimuth and distance of an earthquake epicentre in the same hemisphere as the recording station, the coordinates of the epicentre are located on the base net and the required quantities are then read off from the circles. For an epicentre in the opposite hemisphere, one works from the antipode of the epicentre instead of from the epicentre itself and obtains a value of 180°-Delta instead of Delta The necessary reversal may be obtained by reading the longitudes from the opposite end of the equatorial scale when entering the coordinates of the epicentre, by reading latitude South in place of North, and by numbering the distance circles so that the smallest one corresponds to the largest value of Delta. The azimuth of the antipode of the epicentre will be in the opposite quadrant from the azimuth of the epicentre itself.
A general method for determining azimuth is given in Fig. 1.2.5. Here, one enters the difference between the epicentral and station longitude on the vertical scale, and draws a line down to the epicentral distance on the horizontal scale of the inner part of axes. A parallel line from the epicentral latitude on the horizontal scale of the outer pair of axes cuts the vertical scale at a point corresponding to the azimuth alpha. Uncertainties relating to the sign or quadrant of the azimuth may be resolved by inspection of a sphere or stereogram.
Additional nomograms and charts are reviewed in the papers of Karnik (1955) and Tsuboi (1951).
To measure the approximate azimuth and distance between the epicentre and the station, a large terrestrial globe can be used. If a transparent plastic cover fitting over the globe with latitude and longitude grid lines is prepared, the measurement will be much easier, putting the pole of the cover at the station or the epicentre. The globe should be held in an open-topped frame, to avoid the obstruction of the conventional fixed axis through the North and South poles. The method of making commercially available globes is, however, liable to introduce errors of one or two degrees, and the method is, therefore, less accurate than the best of the other graphical methods.
The major advantage of using geocentric coordinates instead of geographic ones is that the geometry of ray paths through an ellipsoidal earth is greatly simplified. The order of magnitude of the effects with which we are concerned can be seen from the fact that the time of propagation of P waves along an equatorial diameter of the Earth is about 4 seconds longer than that required for propagation along the polar axis. The 'Jeffreys-Bullen' travel-time tables are worked out for a fictitious spherical earth, defined so that the volume within any constant-velocity sphere is equal to that of the corresponding ellipsoid in the real earth.
The exact tabulation of the ellipticity corrections requires a triple-entry table, in which the corrections are found in terms of epicentral latitude and the azimuth and distance of the observing station (Bullen, 1937a). A much more convenient approximation (Bullen, 1937b) is obtained by putting the correction in the form
where h0 and h1 are the heights of the epicentre and recording station above the standard sphere. The values of h as a function of latitude, and of f(Delta) are given in Tables 2 and 2a below (in Bullen's original paper, the times given in the table for h should be reversed in sign). Tables for other phases are also available (Bullen, 1938a and b)).
If the latitudes used for calculating and for entering the table are geocentric, the corrections are accurate to within ±0.2 s for distances out to 90°, increasing to ±0.3 s for certain paths over greater distances. In theory, they may be reduced to ±0.08 s by the use of a 'seismological latitude' phis and slightly modified values of f(Delta). phis is defined by the equations
where phic and phig are the geocentric and geographic latitudes respectively. In fact, the theoretical improvement achieved by the use of this formula is barely significant in view of the actual departure of the earth's interior from the assumed ellipsoidal form.
For deep-focus earthquakes (Bullen, 1938c) the ellipticity corrections are almost the same as those for a surface focus at the same epicentre. The altitude of the station theoretically requires special treatment in so far as the high points of the earth's crust generally reflect the presence of great thicknesses of low velocity material. For this purpose, Bullen suggests an additional correction of 0.14 s per kilometre of station altitude for P waves, and 0.25 s per kilometre for S waves. In fact, recent evaluation of station corrections (see, for example, Herrin and Taggart, 1968) have shown that altitude is not the only factor, or even the dominant factor, in determining the station correction, and the progressive refinement of corrections based directly on observations seems to be the most appropriate course of action.
The concept of magnitude was originally introduced to provide an instrumental measure of the size of earthquakes, using measurements of earth motion adjusted to take account of epicentral distance and focal depth. Later development led to the utilisation of many different earthquake phases, and to efforts either to 'unify' the results obtained by different methods into a common measure of earthquake energy, or to utilise differences between individual estimates for a given earthquake to build up a picture of the character of the source. Magnitude statistics have thereby come to provide a basis for studying and comparing the seismic activity of different regions of the earth, or for studying the variation of activity with time in a given region.
The first approach to this problem was the 'local magnitude' ML designed by Richter (1935) for the classification of local shocks in Southern California, and defined as
where a is the maximum trace amplitude recorded by the 'standard' Wood-Anderson torsion seismometer (T0 = 0.8 s, h = 0.8, V0 = 2800) at a given distance, and a0 is that for the shock of zero magnitude at the same distance.
The amplitude function a0(Delta) was determined empirically for distances Delta = 25-600 km (Richter, 1935) and later supplemented for 0 < Delta < 25 km (Gutenberg and Richter, 1942). Richter (1935) also gave log a = 3.37-log Delta (a in millimetres and Delta in kilometres) for 200 < Delta < 600 km with the hope of extending the formula to larger distances, although this extension revealed itself in later years to be unapplicable over regional distances. Outside California, the Richter scaling function of distance, and the relationship with other magnitude scales, is of provisional utility only.
It should be noted that the focal depth of California earthquakes lies in a very limited range, say 10-20 km, while in many other seismic regions the depth range is much wider than in California. Thus a considerable change in the amplitude function is expected for the deeper shocks in other regions.
The content of the amplitude tables is included in a nomogram (Fig. 3.1.1) in which the original format has been modified by Eiby and Muir (1961). To use the figure set a straight-edge between the distance or P-S interval on the left-hand scale and the maximum trace amplitude on the right-hand scale, and read ML where it crosses the centre scale. If a two-component Wood-Anderson system is available, the arithmetic mean of the two values of amplitude should be used, irrespective of the arrival times.
If a Wood-Anderson instrument is not available, records from other short-period seismographs can be used if the difference of magnification is properly allowed for. For this purpose, a magnification curve for the standard Wood-Anderson instrument is provided (Fig. 3.1.1a). If the magnification of an available seismograph is plotted on this chart, the separation of the curves for earth motion of any given period will be a measure of the amplitude ratio expected on records taken by the two instruments.
It is convenient to use the logarithmic scale of the figure like a slide rule. After measuring the amplitude and period of an appropriate phase on the seismogram, make two marks on a strip of paper to indicate the sensitivity ratio, and move the paper until the mark which represents the sensitivity of the available seismograph comes against a grid line which corresponds to the amplitude of the observed signal. The second mark on the paper will then indicate the trace amplitude which would have been recorded by a Wood-Anderson seismograph on the same pier.
In performing this operation, note that any difference between the shapes of the two response curves will affect the relative amplitudes of earthquake phases of different period. It may therefore happen that the largest amplitude recorded by the available seismograph may not correspond to the same phase as that which would give the largest record on the Wood-Anderson. For this reason, several phases should be read, and the one which gives the largest converted amplitude should be used.
More direct methods, based on the determination of true ground motion, are now available for assigning magnitudes to near earthquakes. Calibration functions for a number of areas have been published (see, for example, Karnik, 1963; Aranovich and others, 1966, etc.). Some examples which are tied to teleseismic magnitudes are set out in Sections 3.2 and 3.3 below.
Methods of determining magnitude from intensities experienced in the vicinity of the epicentre are summarised in Section MAC 3.2.3.
Further development led to the extension of the magnitude scale to all epicentral distances, to shocks deeper than normal and to the use of an estimate of particle velocity of the ground instead of trace amplitude.
Body and surface-wave magnitudes are defined as
where A is a ground amplitude in microns* (in the SI now called micrometers, µm), T is a period in seconds, corresponding to the particle velocity maximum (A/T)max, in a wave of a particular type; the empirical amplitude-distance function sigma (Delta, h), a so-called calibrating function variant for different types of waves, expresses a change in the (A/T)max value with epicentral distance Delta, and focal depth h, and corresponds to zero magnitude. Sigma delta mb and Sigma delta MS are station corrections or regional ones allowing for effects of a recording site, wave path, variation of depth, focal mechanism, etc. The procedure for measuring amplitude and period is given in OP 3.4.
In routine practice, the last term is often assumed to be zero, but it can be useful for individual stations to study their corrections for mb and MS in order to apply them to improve preliminary estimates of magnitudes, before obtaining information from regional and/or international centres. It is especially useful for a local tsunami warning service. Regional characteristics of the corrections for different epicentral regions may also be worthy of note.
*NEIS and ISC use A for mb in millimicrons (now called nanometres (nm) in the SI).
Magnitudes which can be determined by body-waves are mPV and mPH by P waves, mPPV and mPPH by PP waves, and mSH by S waves. The classical study is that of Gutenberg and Richter (1956) in which the calibrating functions, called Q(Delta) are given.
At the IASPEI General Assembly in Zurich (1967) the Committee on Magnitudes recommended stations to report the magnitude for all waves for which calibration functions are available, as well as publishing amplitude and period values separately.
Apart from the multiplicity of possible phases, the following problems arise in practice:
None of the recent studies has received world-wide endorsement for routine use, and the major international agencies are therefore continuing to apply the tables and nomograms of Gutenberg and Richter (1956) as recommended in 1967 by the Committee on Magnitudes. In the case of NEIS and ISC, the data stream is dominated by the input from short-period instruments, and the instructions for reading required operators, until recently, to concentrate attention on the first few cycles of the record. The statistics of the resulting magnitude estimates therefore differ markedly from those based on older data, and from those which depend on the instruments and practice of the USSR and affiliated networks.
It is hoped that the revised practices now recommended in sections OP and OUT of this Manual will bring about a progressive improvement of the situation. In the meantime, it is imperative to ensure that all studies based on the comparison or merging of sets of magnitude data should include adequate specification for class of instrument, reading procedure and method of data reduction.
For shallow earthquakes at distances less than 20°, the amplitude distance functions vary substantially from place to place, and regional studies must be invoked. Some examples are Vanek and Stelzner (1960) for Europe, Jordan and others (1965) for North America, Antonova and others (1967) for Central Asia and Soloviev and Solovieva (1970) for the Far East.
Miyamura (1974) proposed a different approach by using PKP, and constructed provisional calibrating functions for the New Zealand and Macquarie Islands regions.
Calibrating functions for depths greater than 30 km are given by Gutenberg and Richter (1956). Their figures 3-5 are repeated as Fig. 184.108.40.206a, Fig. 220.127.116.11b, and Fig. 18.104.22.168c of this Manual.
Historically the surface-wave magnitude was first calculated by Gutenberg and Richter (1936) as an extension of local magnitude at teleseismic distances using the trace maximum of the Wood-Anderson seismographs. Gutenberg (1945a) gave a formula based on maximum horizontal ground amplitudes, and from 1949 to 1959 various authors calculated similar empirical magnitude formulae for different stations (cf. Bath, 1966).
Soloviev (1955) proposed the use of the maximum ground particle velocity (A/T)max instead of the maximum ground displacement Amax and the corresponding calibrating functions were obtained by Soloviev and Shebalin (1957), Vanek and Stelzner (1959), Christoskov (1965) and others.
Collaboration between research teams in Prague and Moscow (Karnik and others, 1962, and Vanek and others, 1962) yielded the standardised calibrating function
for the epicentral distances 2° < Delta < 160°, as a conventional synthesis of the then published calibrating functions for Amax and (A/T)max. The IASPEI Committee on Magnitudes, Zurich 1967, recommended the use of this standardised formula for sigma(Delta) and the determination of the station and epicentre correction Sigma deltaM The mean periods corresponding to the maximum amplitudes of surface waves to be used for the MS determination were also given, and are reproduced in the following table:
When the period differs considerably from the values given in Table 22.214.171.124, it may be advisable not to use the data for the MS determination.
In order to find (A/T)max on the seismogram, calculate A/T for several trace maxima and select the largest value of A/T among them. Combine the two horizontal components vectorially in the form
irrespective of the types of wave (Rayleigh or Love) and the arrival times. If one of the two components is absent, use sqrt(AN) or sqrt(2AE).
Recent studies, e.g. Evernden (1971), Marshall and Basham (1973) and others, have shown that the condition of amplitude decrease with distance expressed by the term 1.66 log Delta of the standard calibrating function is not met for distances shorter than 25°.
The magnitude determination from the horizontal components arose historically when seismological stations were principally equipped only with horizontal instruments. During the last 20 years the stable operation of medium and long-period vertical seismometers has become easier and amplitude and period data of the vertical components have increased considerably. In addition, it is theoretically more desirable to use the vertical components than to use the horizontal ones, because the former contain exclusively Rayleigh type waves, while the latter contain superposed Love and Rayleigh waves which introduce instability in MS values to hamper the determination procedure.
Many seismological stations and centres have published MS or MLH and MLV routinely using different calibrating functions, but current practice is to use the standard formula for all components. The NEIS and ISC collect the amplitude and period data of surface waves and calculate MLH and MLV using the standard calibrating function for (A/T)max with periods T = 20±2 s at 20 < Delta < 160°. NEIS has published estimates of MS by averaging values of MLH at individual stations, but adopted the vertical component as its standard in May 1975. The ISC takes the average of available A/T values for all components at individual stations, but does not adopt a mean value for the complete earthquake.
Most of the routine determinations of surface-wave magnitudes have been made by using (A/T)max for the period range 18-22 s, without specifying the type of waves, although the standard calibrating function is assumed valid for smaller and larger period waves. Recently there has been a tendency to determine the surface-wave magnitude specifying the type of the waves used, e.g. MLRH or MLRV from Rayleigh waves and MLQH from Love waves, over a wider range of period.
Brune and King (1967) showed that the amplitudes of the surface waves at a 100-second period may be a better measure of the energy rating for larger earthquakes than MS at 20 s. Recent studies indicate that comparing MS at 30-40 s with mb is more effective than using MS at 20 s for the discrimination between underground nuclear explosions and natural shallow earthquakes (Marshall and Basham, 1972). Considering these results, we may see more use of surface-wave magnitudes at different periods, e.g. MS20, MS40, MS100 and so on, in conjunction with multiperiod estimates of mb to describe the earthquakes more completely.
Although historically the local magnitude ML, the surface-wave magnitude MS, and the body-wave magnitude mb were intended to coincide with each other, separate determinations revealed inconsistencies, especially between MS and mb. Gutenberg and Richter (1956a, b) finally obtained the relation
where a = 1/4, b = 7 and a = 0.37, b = 6.76 were possible pairs of values.
Although Richter had thought of magnitude as only a simple measure of size, Gutenberg was insistent on the need to define a magnitude representing the earthquake energy and Gutenberg and Richter (1956b) eventually defined the unified magnitude m as 'the weighted mean of mb and mS, where mb is the body-wave magnitude as found directly from body waves and mS, the corresponding value derived from MS by applying the relation
MS being the surface-wave magnitude as found directly from surface wave with period around 20 s. They also defined Mb as the value of surface-wave magnitude derived from mb and M without subscript as weighted mean between MS and mb. Using the unified magnitude m, Gutenberg and Richter (1956b) defined the earthquake energy E by the relation
which yields the often-used relation
after applying the relation
and putting mb = m and MS = M (Richter, 1958).
The basic thesis which underlies the concept of unified magnitude is that spectral content (which determines the partition of energy between body waves and surface waves) is a function of magnitude only. Modern seismologists, as we have seen in Section 3.2.2 above, are increasingly conscious of a multiplicity of constituents in earthquake character. Unification of magnitude data, therefore involves the loss of some of the information which was contained in the original magnitudes, and current work on the transference of magnitude estimates from one scale to another is tending to be aimed at the provision of an infill of approximate data for sets of events for which direct determinations do not provide complete cover.
Many agencies are faced with the problem of making magnitude estimates from the data of their own stations, independently of those listed by regional or international agencies. The general procedure is to adopt a measure of particle motion derived from available seismograms, and to fit a calibration function to magnitude data from an external agency for a set of events recorded at both levels. When this has been done, the derived formula becomes available for subsequent internal use. Some examples are as follows:
The Japan Meteorological Agency (JMA) has been publishing the magnitudes of Japanese shallow earthquakes using the formula
where Delta is the epicentral distance in kilometres, and AN and AE are the maximum ground amplitudes in microns measured on the N-S and E-W components of horizontal Wiechert seismographs in JMA stations (Tsuboi, 1954). M is the magnitude equivalent to that given Gutenberg and Richter (1949, 1954), which appears to correspond with MS as defined in Section 3.2.2. above.
For deeper earthquakes (h „ 60 km), JMA uses the formula
where K is a depth-distance factor given by Katsumata (1964).
A quantity called mb* has been defined for distances as short as 200 km by the formula
where V is the maximum ground velocity in microns/second in the P-wave train, and Delta is measured in kilometres. The equation was originally produced for use of surface explosions in the United States (Navarro and Brockman, 1970), and has been checked in Europe against underwater explosions (see, for example, Jacob and Willmore, 1972) and a few earthquakes. The available data show good agreement with teleseismic estimates of mb.
Bistricsany (1958) proposed the use of the duration of the surface-wave train for the determination of magnitude, and obtained the following formulae for Wiechert seismographs at the Budapest station:
where F and eL are the end and commencement times (in minutes) of the recorded surface waves, Delta is the epicentral distance in degrees and M is the magnitude equivalent to the surface-wave magnitude determined at Praha station.
Tsumura (1967) obtained
for local earthquakes in Kii Peninsula, Central Japan, where F - P is the total duration time in seconds, Delta is the epicentral distance in kilometres and M is the magnitude equivalent to MJMA.
With small epicentral distances the term containing Delta may be omitted in all of the formulae.
Owing to the narrow dynamic range recorded by a single seismometer, it is often impossible to measure the maximum amplitude of strong seismic movements, and the magnitude determination by duration is a useful and simple way of overcoming this problem. Thus, the method has been adopted recently at many seismological stations with their own definition of the duration time and empirical magnitude formula. The method is especially convenient for local shocks in determining the magnitude without requiring an exact knowledge of epicentral distance.
We have seen above that a single determination of magnitude is not enough to specify completely the spectral content of an earthquake source. ML for instance depends on the spectral level at a period of 0.8 s. Similarly, seismic moment, defined as M0 =µAu can be thought of as representing the zero-frequency spectral level. (Here A is the area of the ruptured fault, u the displacement of the fault movement, and µ the shear modulus of the source rock.) For a simple source of the type proposed by Brune, combined measurements of M0 and ML determine the amplitude spectrum completely, and enable the associated parameters source area, fault displacement, stress drop and seismic energy to be found. Seismic moment may be determined directly from the zero-frequency level of the radiation spectrum obtained from spectral analysis of earthquake records, or from the sum of the envelopes of the earthquake records on the three components of a standard long-period seismograph (Brune and others, 1963).
In the course of the joint General Assembly of the IASPEI the Subcommission on Magnitudes of the Commission on Practice held a number of discussions on the current status of magnitude determinations, and on the nomenclature to be used in future publications. The principal conclusions are summarised below.
The increasing use of telegraphic communication and computer processing has made it necessary to express magnitude information within the constraints of the upper-case alphabetic character set, without using subscripts. The following system is therefore recommended
When a magnitude estimate is derived from a specific instrument, it is desirable to indicate the class of instrument by adding the class letter as defined in INST 1.1, giving complete descriptions such as MPVA, MPVB, MSHC, etc. Note that this procedure is not applicable to the unified magnitude M.
Studies by Gorbunova and Kondorskaya (1977) have established the following relationship between MPVC and MPVA:
The importance of this distinction stems from the fact that the early work of Gutenberg and the current determinations by Eastern countries are based largely on Class C instruments, whereas the Western agencies have relied strongly on data from Class A.
BULLEN, K. E. 1963. Introduction to the theory of seismology. (Cambridge: University Press, third edition.) 381 pp.
COMRIE, L. J., JEFFREYS, H. 1938. The geocentric direction cosines of Seismological Observatories. (London: British Association for the Advancement of Science.)
JEFFREYS, H., and BULLEN, K. E. 1940, 1958. Seismological tables. (London: British Association for the Advancement of Science; reprinted by the ISC, Newbury.)
and SHIMSHONI, M. 1964. The times of pP, sS, S, sP and pS. Geophys. J., R. Astr. Soc., Vol. 8, pp. 324-337.
KARNIK, V. 1955. Zur graphischen Distanz- und Azimuthbestimmung eines Epizentrums. [Graphical determination of epicentral distance and azimuth.] [Trav. Inst. Geophys. Ac. Tche. Cosl. Sc. (Geofisikalni Sbornik) No. 22, pp. 45-56.
MACELWANE, J. B., and SOHON, F. W. 1949. Introduction to theoretical seismology. (St Louis: University).
RICHTER, C. F. 1958. Appendix II, VIII, in Elementary seismology. (San Francisco: W. H. Freeman.)
TSUBOI, C. 1951. A simple instrument useful for finding the angular distance on a point from a fixed point (phi, lambda). Geophys. Notes Geophysical Inst., Tokyo Univ., Vol. 4, p. 1.
WILLMORE, P. L. 1957. The use of the stereographic net for the rapid approximate determination of distance and azimuth. Earthquake Notes Vol. 28, pp. 4-6.
and HODGSON, J. H. 1955. Charts for measuring azimuth and distance and for tracing seismic rays through the earth. Pub. Dom. Obs., Vol. 14, No. 14, pp. 405-414.
BULLEN, K. E. 1937a. A suggested new seismological latitude. Mon. Not. RAS (Geophys. Supp.), Vol. 4, pp. 158-164.
1937b. The ellipticity correction to travel-times of P and S earthquake waves. Mon. Not. RAS (Geophys. Supp.), Vol. 4, pp. 143-157.
1938a. Ellipticity corrections to waves through the earth's central core. Mon. Not. RAS (Geophys. Supp.), Vol. 4, pp. 317-331.
1938b. Ellipticity corrections to earthquake waves reflected at the central core. Mon. Not. RAS (Geophys. Supp.), Vol. 4, pp. 332-335.
1938c. The ellipticity correction and deep-focus earthquakes. Mon. Not. RAS (Geophys. Supp.), Vol. 4, pp. 469-471.
HERRIN, E., and TAGGART, J. N. 1968. Regional variations in P travel times. Bull. Seism. Soc. Am., Vol. 58, pp. 1325-1337.
JEFFREYS, H. 1935. On the ellipticity correction in seismology. Mon.. Not. RAS (Geophys. Supp.), Vol. 3, pp. 271-274.
KOLOSENKO, M. N. 1958. Allowance for the ellipticity of the earth for epicentral determinations. Izv. Akad. Nauk. SSSR Ser. Geofiz., No. 1.
ADAMS, R. D. 1977. Survey of practice in determining magnitudes of near earthquakes. WDCA Report SE-8. (Boulder: World Data Center A for solid-earth geophysics.) 67 pp.
ANTONOVA, L., and HALTURIN, V. 1967. The periods of seismic waves during earthquakes. In Papers presented at the Ninth Assembly of the European Seismological Commission, pp. 125-135. (Copenhagen: Akademisk Forlag.)
ARANOVICH, Z. I., VVEDENSKAYA, N. A., GUBIN, I. E., KONDORSKAYA, N. V., FEDOROV, S. A., and SHEBALIN, N. V. 1966. Instruktsiya o poryyadke proizvodstva i obrabotik nablyudeniy na seismichekikh stantsiakh yedinoy systemy seismicheskikh nablyudeniy SSSR. [Instructions on operation and interpretation for the seismological stations of the unified seismological network system of the USSR.] (Moscow: Institut Fiziki Zemli Akad. Nauk SSSR.)
BATH, M. 1966. Earthquake energy and magnitude. Phys. & Chem. of the earth, Vol. 7, pp. 115-165.
1967. Handbook on earthquake magnitude determinations. (Uppsala: Seismol. Inst.)
BISZTRICSANY, E. 1958. On the problem of magnitude determination. Zeitch. f. Geophys., Bd. 24, pp. 153-160.
BRUNE, J. N., ESPINOSA, A., and OLIVER, J. 1963. Relative excitation of surface waves by earthquake and underground explosions in the California-Nevada Region. J. Geophys. Res., Vol. 68, pp. 3501-3513.
and KING, C. Y. 1967. Excitation of mantle Rayleigh waves of period 100 seconds as a function of magnitude. Bull. Seism. Soc. Am., Vol. 57, pp. 1355-1365.
BUNK, V. I., VVEDENSKAYA, N. A., GORBONOVA, I. V., KONDORSKAYA, N. V., LANDYREVA, N. S., and GEDOROVA, I. V. 1970. Correlation of MLM and mpv by data of the network of seismic stations of the USSR. Geophys. J. R. Astr. Soc., Vol. 19, pp. 533-542.
CARPENTER, E. W., MARSHALL, P. D., and DOUGLAS, A. 1967. The amplitude-distance curve for short-period teleseismic P-waves. Geophys. ]. R. Astr. Soc., Vol. 13, pp. 61-70.
CHRISTOSKOV, L. 1965. Magnitude-dependent calibrating functions of surface waves for Sofia. Studia Geoph. et Geod., Vol. 9, pp. 331-340.
EIBY, G. A., and MUIR, M. G. 1961. Tables to facilitate the study of near earthquakes. Seismological Observatory Bulletin S-109. (Wellington: Geophysics Division, New Zealand Department of Scientific and Industrial Research.)
EVERNDEN, J. F. 1971. Variation of Rayleigh-Wave amplitude with distance. Bull. Seism. Soc. Am., Vol. 61, pp. 231-240.
GORBUNOVA, I. V., and KONDORSKAYA, N. V. 1977. Magnitudes in the seismological practice of the USSR. Izv. an. SSSR, ser Fiziki Zemli, No. 2, Moscow. [In Russian.]
GUTENBERG, B. 1945a. Amplitudes of surface waves and magnitudes of shallow earthquakes. Bull. Seism. Soc. Am., Vol. 35, pp. 3-12.
1945b. Amplitudes of P, PP and S and magnitudes of shallow earthquakes. Bull. Seism. Soc. Am., Vol. 35, pp. 57-69.
1945c. Magnitude determination of deep focus earthquakes. Bull. Seism. Soc. Am., Vol. 35, pp. 117-130.
and RICHTER, C. F. 1936. On seismic waves. Gerlands Beitrage zur Geophysik, Vol. 47, pp. 73-131.
1942. Earthquake magnitude, intensity, energy and acceleration. Bull. Seism. Soc. Am., Vol. 32, pp. 163-191.
1st edition 1949, 2nd edition 1954. Seismicity of the earth and associated phenomena. Second edition. (Princeton University Press.) Reprinted 1965. (New York: Hafner.) 310 pp.
1956a. Magnitude and energy of earthquakes, Annali di geofisica, Vol. IX, pp. 115.
1956b. Earthquake magnitude, intensity, energy, and acceleration. Bull. Seism. Soc. Am., Vol. 46, pp. 105-143.
JACOB, A. W. B., and WILLMORE, P. L. 1972. Teleseismic P waves from a 10-ton explosion. Nature, Lond. Vol. 236, pp. 305-306.
JORDAN, J., BLACK, R., and BATES, CH. 1965. Patterns of maximum amplitudes of Pn and P waves over regional and continental areas. Bull. Seism. Soc. Am., Vol. 55, pp. 693-720.
KARNIK, V. 1963. Seismicity of Europe III. IUGG Monograph, No. 23, pp. 1-23.
1968. Seismicity of the European Area, Part I (Prague: Academia.) 3644 pp.
1969. Seismicity of the European Area, Part II (Prague: Academia.) 218 pp + maps.
KONDORSKAYA, N. V., RIZNICHENKO, Yu. V., SAVARENSKY, YE. F., SOLOVEV, S. L., SHEBALIN, N. V., VANEK,J., and ZATOPEK, A. 1962. Standardisation of the earthquake magnitude scales. Studia Geophys. et Geod., Vol. 6, pp. 41-48. Also as VANEK et al. 1962. Akad Nauk USSR Izv. Ser. Geofix., No. 2, pp. 153-158. (In Russian. English translation of Russian version, pp. 108-111, Am. Geophys. U., 1962.)
KATSUMATA, M. 1964. A method to determine the magnitude of deep focus earthquakes in and near Japan. Zisin (J. Seism. Soc. Japan) Ser. 11, Vol. 17, pp. 152-165.
MARSHALL, P. D., and BASHAM, P. W. 1972. Discrimination between earthquakes and underground explosions using an improved Ms scale. Geophys. J. R. Astr. Soc., Vol. 28, pp. 431-458.
1973. Rayleigh wave magnitude scale Ms. Pure and Appl. Geophys., Vol. 103, pp. 237-434.
MIYAMURA, S. 1974. Determination of body-wave magnitudes for shallow earthquakes in New Zealand and Macquaries Loop regions using PKP data. Phys. Earth planet. Interiors. Vol. 8, pp. 167-176.
NAVARRO, R., and BROCKMAN, F. R. 1970. Seismic activity in September 1969 near Ruliston Test Site. Publication USCGS, No. 746-5.
RICHTER, C. F. 1935. An instrumental earthquake magnitude scale. Bull. Seism. Soc. Am., Vol. 25, pp. 1-32.
1953. Elementary seismology (San Francisco: W. H. Freeman.) 768 pp.
SOLOVIEV, S. L. 1955. O klassifikatsiy zemletrayaseniy po velichine ikh energii. (Classification of earthquakes in order of energy.) Trudy Geofiz. Inst. AN SSSR, No. 30 (157), pp. 3-31.
and SHEBALIN, N. V. 1957. Opredelenie intensivnosti zemletryaseniya po smeshcheniyu pochvy v poverkhnostynkh. (Determination of intensity of earthquakes according to ground displacements in the surface waves.) Izv. AN SSSR, ser. geopfiz., No. 7, pp. 926-930.
and SOLOVIEVA, O. N. 1970. Amplitudnye krivye voln PV, PH, SH negluboko-forkusnykh tikho-okeanskikh zemletryaseniy na srasstoyaniyakh 2°-40°. [Amplitude curves of PV, PH, SH waves of shallow-focus pacific-ocean earthquakes from distance of 2°-40°.] Vortrage des soprones Symposiums der 4. Subkomission von KAPG.
TSUBOI, C. 1954. Determination of the Gutenberg-Richter's magnitude of earthquakes occurring in and near Japan. Zisin. (Japan: J. Seism. Soc.) Ser. II, Vol. 7, pp. 185-193.
TSUMARA, K. 1967. Determination of earthquake magnitude from total duration of oscillation. Bull. Earthq. Res. Inst., Tokyo, Vol. 45, pp. 7-18.
VANEK, J., and STELZNER, J. 1959. Bestimmung der Magnitudengleichen fur Jena. (Determination of the magnitude function for Jena. ) Gerl. Beitr. zu Geophys., Bd 77, pp. 105-119.
and STELZNER, J. 1960. Magnitude calibrating functions for body waves. Ann. di Geofis., Vol. 13, pp. 393-407.
et al. 1962. Standardisation of magnitude scales. See KARNIK et al. 1962.
KONDORSKAYA, N. V., and CHRISTOSKOV, L.1977. Magnitude in Seismological practice, Vol. 1, PV and PVs waves. (Sofia: Bulgarian Academy of Science.) [In Russian, with extensive English abstract.]
VEITH, K. F., and CLAWSON, G. E. 1972. Magnitude from short-period P-wave data. Bull. Seism. Soc. Am., Vol. 62, pp. 435-452.
VESANAN, E. 1942. Uber die typenanalytische auswertung der seismogramme. (On the character interpretation of seismograms.) Ann. Acad. Scient. Fennicae, Ser. A. III, No. 5.
Date created: 8/10/96 Last modified: 9/9/97 Copyright © 1996, 1997, Global Seismological Services Maintained by: Eric Bergman email@example.com