The great British earthquake of 1984 seen from afar

Summary. The Lleyn Peninsular earthquake which occurred at 0656 on 1984 July 19 was recorded at three medium aperture seismological array stations located at teleseismic distances. From these recordings an estimate is made of the hypocentre, origin time, magnitude and fault-plane solution of the earthquake. The fault-plane solution was determined using the Pearce algorithm and indicates that the focal mechanism is predominantly strike-slip. The fault-plane solution was used to generate synthetic seismograms for comparison with the observed to confirm the nature of the source and in particular the depth of focus which was estimated to be 20.5 km. It is concluded that the determination of the earthquake parameters using only three teleseismic seismograms is in good agreement with the results obtained from an analysis of 45 local and regional seismological stations.     

Applying the relative hypocentre location approach: where was the 1980 November 23 Irpinia earthquake?

The purpose of this work is to evaluate under what conditions it is feasible and with what accuracy it is possible to locate the nucleation point of a large earthquake, given the availability of aftershocks located with high precision by the deployment of a local network. We experiment with several approaches and apply them lo the location of the epicentre of the 1980 November 23 Irpinia earthquake ( M ^w= 6.9).
First we use local P ^g phases selected to optimize the azimuthal coverage, obtaining a well-constrained location with a small statistical error, which typically underestimates the true hypocentre uncertainty.
We then exploit the relative location technique, obtaining stable, almost coincident solutions under three conditions: (1) using multiple independent master events to derive an average epicentre; (2) fitting simultaneously the larger data set for all available master events, using a forward approach; (3) conducting an a priori evaluation of the statistics of station and master events to separate model uncertainties and improve the statistical accuracy of the relative locations. Moreover, only by introducing station statistics can we achieve the desired accuracy of ≅ 1 km in constraining the rupture nucleation point of this large earthquake, and we show that the application of the relative location technique to uncleaned, unweighted data for a single master event provides only a crude epicentre with a confidence ellipse deceivingly smaller than the true hypocentre uncertainty.
The revised epicentre for the 1980 November 23 Irpinia earthquake (48.803 °N-15.302°E) validates the class of multidisciplinary reconstructions of the source process such as the model of Valensise et al. (1989), based on the hypocentre of Westaway & Jackson (1987), and is shifted by almost 13 km to the NW of the epicentre recently proposed by Westaway (1992).     

The PWaves From the Amchitka Island Explosions

Summary. There is evidence that the equivalent seismic sources of the Amchitka Island explosions — Longshot, Milrow and Cannikin — depart significantly from the simple model of a point compressional-source in a layered elastic-medium. Consequently modelling the observed seismograms using standard source-models may not be the most efficient method of determining source properties. Here an alternative to modelling is used to obtain information on the seismic sources due to the explosions. Broad-band (BB) estimates of the P signals are obtained from the short-period (SP) seismograms, corrected for attenuation, and interpreted in terms of P, pP and radiation from secondary sources. the main conclusions are:
(i) BB estimates of the radiated displacement from the explosions can be obtained with only a small reduction in the signal-to-noise ratio seen on SP seismograms;
(ii) observations of differences in pulse amplitudes and spectra are not necessarily due to differences in anelastic attenuation;
(iii) P and pP at a given station may differ in shape so that notches in the signal spectrum may not be related to source depth;
(iv) there is evidence of arrivals that others have identified as due to slap-down but which could be interpreted as an overshoot to pP;
(v) direct interpretation of the estimated ground displacement is a better procedure for determining the seismic source properties of explosions than modelling SP seismograms using idealised models as a starting point.     

Inverting seismograms of weak events for empirical Green's tensor derivatives

Given a station and a suite of clustered weak event records whose focal mechanisms are known, empirical Green's tensor spatial derivatives (EGTD) can be calculated, almost free of the source effects. The inversion method solves an overdetermined system of linear algebraic equations in the frequency domain, and a FFT provides the time histories. They represent a least-squares averaged representation of the path between the station and the entire focal zone, including multiple P and S phases, coda waves, etc. The method can be used for modelling earthquakes whose focal mechanisms are quite different from the weak events employed. This is proved here using both synthetic seismograms and earthquake swarm records with variable focal mechanisms. It is expected that the EGTD method can improve strong motion syntheses and source inversions.     

Determination of earthquake early warning parameters, τc and Pd, for southern California

We explore a practical approach to earthquake early warning in southern California by determining a ground-motion period parameter  τ _c   and a high-pass filtered displacement amplitude parameter Pd from the initial 3 s of the P waveforms recorded at the Southern California Seismic Network stations for earthquakes with M > 4.0. At a given site, we estimate the magnitude of an event from  τ _c   and the peak ground-motion velocity ( PGV ) from Pd . The incoming three-component signals are recursively converted to ground acceleration, velocity and displacement. The displacements are recursively filtered with a one-way Butterworth high-pass filter with a cut-off frequency of 0.075 Hz, and a P -wave trigger is constantly monitored. When a trigger occurs,  τ _c   and Pd are computed. We found the relationship between  τ _c   and magnitude ( M ) for southern California, and between Pd and PGV for both southern California and Taiwan. These two relationships can be used to detect the occurrence of a major earthquake and provide onsite warning in the area around the station where onset of strong ground motion is expected within seconds after the arrival of the P wave. When the station density is high, the methods can be applied to multistation data to increase the robustness of onsite early warning and to add the regional warning approach. In an ideal situation, such warnings would be available within 10 s of the origin time of a large earthquake whose subsequent ground motion may last for tens of seconds.     

Complex P waveforms from a Gulf of Aden earthquake

Summary. It is shown that complex teleseismic P waveforms from a shallow earthquake in a tectonic area can be interpreted using a simple source model embedded in a plane layer velocity structure (with sea layer) whose details are based upon independent evidence. This gives hope that structural complexity in tectonic regions may not always make distant P -wave seismograms impossible to interpret, and that, instead, source complexity may be responsible for some of the many complex waveforms observed, even for earthquakes below magnitude m _b 5.5.     

Source investigation of a small event using empirical Green's functions and simulated annealing

We propose a two-step inversion of three-component seismograms that (1) recovers the far-field source time function at each station and (2) estimates the distribution of co-seismic slip on the fault plane for small earthquakes (magnitude 3 to 4). The empirical Green's function (EGF) method consists of finding a small earthquake located near the one we wish to study and then performing a deconvolution to remove the path, site, and instrumental effects from the main-event signal.
The deconvolution between the two earthquakes is an unstable procedure: we have therefore developed a simulated annealing technique to recover a stable and positive source time function (STF) in the time domain at each station with an estimation of uncertainties. Given a good azimuthal coverage, we can obtain information on the directivity effect as well as on the rupture process. We propose an inversion method by simulated annealing using the STF to recover the distribution of slip on the fault plane with a constant rupture-velocity model. This method permits estimation of physical quantities on the fault plane, as well as possible identification of the real fault plane.
We apply this two-step procedure for an event of magnitude 3 recorded in the Gulf of Corinth in August 1991. A nearby event of magnitude 2 provides us with empirical Green's functions for each station. We estimate an active fault area of 0.02 to 0.15 km² and deduce a stress-drop value of 1 to 30 bar and an average slip of 0.1 to 1.6 cm. The selected fault of the main event is in good agreement with the existence of a detachment surface inferred from the tectonics of this half-graben.     

Time-domain computation of scattering from a heterogeneous layer

Adopting Born and ray approximations, time-domain synthetic seismograms for P-P and P-S scattering from a plane wave incident on a thin, laterally heterogeneous layer are presented in this paper. The time-domain P coda is a convolution between a structure function and the second-order derivative of the time function of the incident P wave. Examples of synthetic seismograms are given using a time function from a computed short-period seismogram for a point explosive source in a half-space. These show that it is impossible, with realistic values of the parameters involved, to generate significant codas when only single scattering is involved.     

A fault plane solution using theoretical P seismograms

We use the method of Hudson and Douglas, Hudson & Blarney to compute seismograms which simulate the codas of 10 short period P -wave seismograms from a shallow earthquake. The polarities and relative amplitudes of P and pP measured from seven of the observed seismograms are used to compute a fault plane solution with confidence limits, assuming that the source radiates as a double couple. This solution is in approximate agreement with that given for the same earthquake by Sykes & Sbar, who used only the onset polarities of short-period P waves. The small difference between the two solutions can be explained by interference between the true first motion of P and microseismic noise at two stations.
The results show that, for some shallow earthquakes, the relative amplitude method has the following advantages over the first motions method. First, a P/pP amplitude ratio (with appropriate confidence limits) can always be measured, even in seismograms which are so noisy that the first motion of P is uncertain. Second, the fault plane solutions obtained from relative amplitudes have known confidence limits. Finally, by using more information from each seismogram, the relative amplitude method requires considerably fewer seismograms than the first motions method.     

Seismic-wave propagation through a cracked solid: polarization as a possible dilatancy diagnostic

summary . A new technique is presented for modelling the elastic constants of cracked structures with application to systems with weak concentrations of parallel cracks, and of simple biplanar and triplanar cracks. The velocities and V_p/V_s ratios of these anisotropic structures are used to provide quantitative models for some earthquake precursors. These results indicate the great importance of crack geometry to the behaviour of precursors. The velocities of saturated cracks appear to favour the dilatancy-diffusion model of precursory phenomena. Synthetic seismograms are calculated for propagation through possible dilatancy zones. The seismograms show some characteristic features which may be useful for the investigation of earthquake dilatancy.     

A global study of S-to-P and P-to-S conversions from the upper mantle transition zone

b
Long-period data of the Global Digital Seismograph Network (GDSN) recorded over the three-year period from 1984 to 1986 were studied for the occurrence of S-P and P-S conversions from the upper mantle transition zone that appear as precursors to teleseismic S arrivals. Conversions of this type were identified on a large number of single-station records. Simple stacking of many records enhanced the appearance of converted phases and demonstrated that no major lateral variations in the nature of the transition zone exist between various tectonic regions. S-P and P-S conversions from the 400 km discontinuity were best observed at distances between 70 and 85 while conversions from the 670 km discontinuity showed up best at distances beyond 87. The analysis of published source mechanisms and comparison with synthetic seismograms suggests that the appearance of converted phases is primarily governed by the earthquake radiation pattern. Phases that have undergone S-P conversions beneath the receiver are best observed from dip-slip events that radiate strong SV - and weak P -waves towards the station. P-S conversions beneath the source area, on the other hand, are frequently observed from events that radiate strong P and little SV energy towards the station, and also from some strike-slip events. Comparison of observed with synthetic seismograms suggests that the PREM model of Dziewonski & Anderson (1981) explains most of the observations. Observed S-P and P-S conversions from the 670 km discontinuity, however, often have larger amplitudes than in the synthetics. Constructive interference of converted waves with the P -wave coda, source radiation effects and a velocity contrast across the 670 km discontinuity which is higher than in PREM may all contribute to the discrepancy.     

Long-period seismic source inversions using global tomographic models

We investigate the effect of laterally varying earth structure on centroid moment tensor inversions using fundamental mode mantle waves. Theoretical seismograms are calculated using a full formulation of surface wave ray theory. Calculations are made using a variety of global tomographic earth models. Results are compared with those obtained using the so-called great-circle approximation, which assumes that phase corrections are given in terms of mean phase slowness along the great circle, and which neglects amplitude effects of heterogeneity. Synthetic tests suggest that even source parameters which fit the data very well may have large errors due to incomplete knowledge of lateral heterogeneity. The method is applied to 31 shallow, large earthquakes. For a given earthquake, the focal mechanisms calculated using different earth models and different forward modelling techniques can significantly vary. We provide a range of selected solutions based on the fit to the data, rather than one single solution. Difficulties in constraining the dip-slip components of the seismic moment tensor often produce overestimates of seismic moment, leading to near vertical dip-slip mechanisms. This happens more commonly for earth models not fitting the data well, confirming that more accurate modelling of lateral heterogeneity can help to constrain the dip-slip components of the seismic moment tensor.     

Earthquake seismograms that show Doppler effects due to crack propagation

Summary. Examples are presented of earthquake P -wave pulses seen on broadband seismograms, to show that on such recordings the pulse shapes are more clearly seen than on conventional short-period and long-period seismograms. Most of the broadband seismograms have been chosen because they show marked differences between the pulse lengths of P and those of the surface reflections. In addition some of the pulses appear to have smooth onsets and abrupt trailing edges so that the onset of the pulse is difficult to observe and the largest amplitude arrivals seen on the seismogram coincide not with the onset of motion but with the termination of motion: that is the large arrivals mark stopping phases of motion.
We assume that the differences in pulse length are due to the effects of a moving source – that is a Doppler effect – and that the pulses with smooth onsets and abrupt trailing edges can be modelled simply by a source propagating on a line with low radiation amplitude at the start of motion. A trial and error method guided by a published fault plane solution is then used to obtain a fit between observed and computed seismograms for one of the earthquakes. This process leads to an estimate of the crack speed of about 1.4 times the 5-wave speed.
The errors that may arise in estimating source depths and orientation, if stopping phases are not recognized as such, is discussed.     

Kinematic multistation discriminator between local quarry blasts and earthquakes

It has long been known that S waves on seismograms of local explosions are often accompanied by strong low-frequency, low-velocity, Rg surface wave trains, often significantly diminished for earthquakes. We utilize this fact to construct a new formal discriminator between earthquakes and explosions by measuring the S -surface-wave group velocity. The method is based on analysing the velogram; that is, the display of the envelope of ground motion versus group velocity V = R/T , where R is the epicentral distance and T  the traveltime. We examine the distribution of seismic energy in time and space using envelopes of records from the Israel Seismic Network (ISN), from which we compute the velograms and observe differences in the velograms of quarry blasts and earthquakes. The data include 143 seismic events occurring in three areas (Galilee, Dead Sea, and Gilad) monitored by the ISN; the magnitude range is M _L = 1.0–2.8 at distances of 15–310 km. From the velograms we measure the group velocity, V _{m s} , within the 1–4 km s⁻¹ range at which the velogram reaches its maximum for each available station. The resulting V _{m s} (R) function is closely fitted by the empirical relationship a + b  ln  R , with a and b coefficients varying from event to event. A simple linear function c = b + 0.33a at a threshold C = 0.69 completely separates ( a,b ) pairs for the 67 Galilee events, and, for the 76 remaining events, one earthquake and four explosions are wrongly classified. After data validation and application of the Fisher linear discriminator, adapted to the events from Galilee, only two misidentified events remain for the whole data set.     

Cagniard de-Hoop method for a Haskell type vertical fault

Summary. There are two ways to apply the Cagniard de-Hoop method when generating synthetic seismograms due to a source in a three-dimensional medium. One is the Hankel transform method (hereafter called 'the cylindrical wave representation method'), which utilizes the property of the modified Bessel function. The other is 'the plane wave representation method', which replaces the Bessel function by a superposition of plane waves. In extending a point source to a finite dimensional source, the latter method is extremely useful, because it enables the integrations on the fault surface to be performed analytically.
Using the latter method, expressions for displacements due to a Haskell type vertical fault in a uniform half-space are obtained. A solution is given as a sum of four quadrantal source contributions, which is similar to Madariaga's solution for a source in the whole space. Each contribution consists of a single finite range integration or a single integration plus a pole contribution, depending on the location of the observation point with respect to the source. The procedure can be extended to other fault types.     

Modelling of ground motion in the Higashinada (Kobe) area for an aftershock of the 1995 January 17 Hyogo-ken Nanbu, Japan, earthquake

We model the ground motion from an aftershock of the 1995 January 17 Hyogo-ken Nanbu (Kobe) earthquake to investigate basin edge effects on wave propagation in Higashinada ward, downtown Kobe. Point-source finite-difference seismograms calculated using a double-couple solution and a 2-D basin structure are compared with the ground-motion velocity seismograms recorded in a small station array deployed at sites within and outside the heavily damage zone in Higashinada ward. The comparison suggests that in the frequency range 0.1-2 Hz that was analysed, the observed spatial amplitude variation of the aftershock ground motion is attributable mainly to the basin edge effects. We found that the basin edge effect, caused by the superposition of the direct S wave and the basin-edge-diffracted waves, amplified the ground motion in a narrow zone that is offset by about 0.7 km from the basin edge.     

Far-field simulation of the 1946 Aleutian tsunami

We present hydrodynamic far-field simulations of the Aleutian tsunami of 1946 April 1, using both a dislocation source representing a slow earthquake and a dipolar one modelling a large landslide. The earthquake source is derived from the recent seismological study by López and Okal, while the landslide source was previously used to explain the exceptional run-up at Scotch Cap in the near field. The simulations are compared to a field data set previously compiled from testimonies of elderly witnesses at 27 far-field locations principally in the Austral and Marquesas Islands, with additional sites at Pitcairn, Easter and Juan Fernández. We find that the data set is modelled satisfactorily by the dislocation source, while the landslide fails to match the measured amplitudes, and to give a proper rendition of the physical interaction of the wavefield with the shore, in particular at Nuku Hiva, Marquesas. The emerging picture is that the event involved both a very slow earthquake, responsible for the far-field tsunami, and a major landslide explaining the near-field run-up, but with a negligible contribution in the far field.     

Shear-wave polarizations on a curved wavefront at an isotropic free surface

Summary. We present polarization diagrams of the particle motions at the free surface of an isotropic half-space generated by incident shear waves from a local buried point source. The reflectivity technique is used to calculate synthetic seismograms from which the particle motions are plotted. The particle motions are examined over a range of epicentral distances in a uniform isotropic half-space for different source frequencies and polarization angles, and for different Poisson's ratios. The particle motions due to a curved wavefront possess different characteristics from those generated by plane wavefronts at corresponding incidence angles. A curved wavefront generates a local SP -phase: a P -headwave which propagates along the free surface, and arrives shortly before the direct S -wave. These two arrivals give rise to cruciform particle motions in the sagittal and horizontal planes, which could be misinterpreted as anisotropy-induced shear-wave splitting. An examination of the particle motion in the transverse plane, mutually orthogonal to the sagittal and horizontal planes, can be used to discriminate between isotropic and anisotropic interpretations. The amplitude of the SP -phase is enhanced when it propagates in a low-velocity surface layer overlying the source layer, and may then become the dominant phase on radial-component seismograms. The presence of even a single surface layer may introduce considerable complexity into the seismogram, and we examine the effects of layer thickness, velocity contrast, and source depth on the corresponding polarization diagrams. Reliable information on the source and propagation path characteristics of shear waves from a buried local point source can only be obtained from free-surface records if they are recorded within a very limited epicentral distance range.