Automatic Location of Swarm Earthquakes from Local Network Data

An advanced method of automated seismic phase picking and exact location and magnitude determination of swarm micro-earthquakes from local network data is presented. The phase picker is applied in two steps: first, S-wave groups are identified using a polarisation detector, and then corresponding P-wave groups are searched for. The times of maximum P- and S-amplitudes are then used as starting points for the determination of accurate P- and S-arrival times. The maximum S-wave amplitudes are utilised for determining local magnitudes. The whole procedure is checked by simultaneous preliminary hypocentre location providing estimates of local magnitudes and a compatibility check of the candidate P- and S-phases. The closest station to the earthquake cluster is used as a master, and the phase search at the remaining stations is governed by the P- and S-phases identified at the master station. Thanks to the use of apriori information on the approximate position of hypocentres, the procedure is also capable of picking the individual P- and S-phases of sequences of overlapping swarm events. The performance of the procedure was tested by comparison of the automatically and interactively created catalogues of the January 1997 NW-Bohemia micro-earthquake swarm. With stations located at epicentral distances between 0 and 20 km, the difference between hypocentre coordinates obtained by automatic and interactive processing did not exceed 80 m for 86% events. All events above magnitude 0.5 were identified, and the automatically determined polarity of first P-wave motion proved to be correct in 89% of them.     

The role of constraints on joint hypocentre determination in seismotectonic studies

In this paper we show how the performance of the joint hypocentre determination (JHD) method can be improved, leading to reduced instability in cases close to singularity. The method has been extended by imposing a number of constraints introduced by other authors, and adding a new one. We tested the stability of the method and the relative advantages of the various constraints by simulating a geometrical space distribution of hypocentres recorded by a regional seismic network. We then applied this method to deep earthquakes that occurred in the Southern Tyrrhenian Sea subduction zone and to the seismicity of the Northern Apennines, which is moderately deeper than the typical shallow seismicity of the Apennines. The results obtained from the analysis of synthetic data and actual earthquakes confirm that the JHD method produces less scatter in the hypocentral determinations with respect to the standard locations. The main conclusion of our study is that we can significantly reduce the systematic mislocations that result from applying JHD to very clustered seismicity if we introduce the appropriate set of constraints.     

New techniques for the analysis of earthquake sources from local array data with an application to the 1993 Scotts Mills, Oregon, aftershock sequence

We analysed aftershocks recorded by a temporary digital seismic network following the moderate M _w = 5.5 1993, Scotts Mills, Oregon, earthquake. A technique to retrieve source moment tensors from local waveforms was developed, tested, and applied to 41 small earthquakes ( M _w ranging from 1.6 to 3.2). The derived focal mechanisms, although well resolved, are highly variable and do not share a common nodal plane. In contrast, the majority of the events, relocated with a joint hypocentre determination algorithm, collapse to a well-focused plane. The incompatibility of the nodal planes of most events with the plane defined by their locations suggests that the aftershocks did not occur on the fault plane, but tightly around it, outlining the rupture area rather than defining it. Furthermore, the moment tensors reveal stable P -axes, whereas T  -axes plunges are highly dispersed. We detect a rotation of average T  -axis plunge with depth, indicating a change from shallower, predominantly dip-slip mechanisms to deeper strike-slip mechanisms. These characteristics are difficult to explain by remnant stress concentrations on the main-shock rupture plane or asperity- and barrier-type models. We suggest that the aftershocks occurred under the ambient regional stress, triggered by a sudden weakening of the region surrounding the main-shock slip, rather than from a shear stress increase due to the main shock.     

Analysis of the seismicity in the central part of the Pyrenees (France), and tectonic implications

Two temporary seismological networks have been set up in 2000 and 2002 in the central part of the Pyrenees, in a region, which appears as a transition between two domains where both the seismic activity and the tectonic regime are different. Together with the permanent networks, they allowed us to obtain precise hypocenter locations for more than 400 events with local magnitudes ranging from 1.5 to 4.6, as well as 30 new focal mechanisms. The seismicity is distributed in several clusters, which are not located along the North Pyrenean Fault, considered as the major tectonic accident resulting from the suture of the Iberian and Eurasian convergent plates when the range formed. Several small fault segments dipping to the north are identified. The maximum focal depth varies from 10 to 20 km, with variations which are roughly parallel to those of the Moho, indicating a thickening of the seismogenic layer to the east of the studied area. The obtained focal solutions reveal a predominance of normal faulting to the West and reverse faulting to the East, with strike-slip motions in between. The largest fault segment to the East, with a length of about 20 km, could possibly be related to a large historical event which occurred in 1660, with intensity IX, close to cities which have become since then important touristic centres.     

Statistical discrimination of foreshocks from other earthquake clusters

When earthquake activity begins, it may be a foreshock sequence to a larger earthquake, a swarm, or a simple main-shock-aftershock sequence. This paper is concerned with the conditional probability that it will be foreshock activity of a later larger earthquake, depending on the occurrence pattern of some early events in the sequence. The earthquake catalogue of the Japan Meteorological Agency (1926-1993, M_J≥4) is decomposed into a large number of clusters in time and space in order to compare statistical features of foreshocks with those of swarms and aftershocks. Using such a data set, Ogata, Utsu & Katsura (1995) revealed some discriminating features of foreshocks relative to the other types of clusters, for example the events' closer proximity in time and space, and a tendency towards chronologically increasing magnitudes, which encouraged us to construct models which forecast the probability of the earthquakes being foreshocks. Specifically, the probability is a function of the history of magnitude differences, spans between origin times and distances between epicentres within a cluster. For purposes of illustration, the models were fitted to the early part of the data (1926-1975) and the validity of the forecasting procedure was checked on data from the later period (1976-1993). Two procedures for evaluating the performance of the probability forecast are suggested. Furthermore, for the case where only a single event is available (i.e. either it is the first event in a cluster or an isolated event), we also forecast the probability of the event being a foreshock as a function of its geographic location. Then, the validity of the forecast is demonstrated in a similar manner. Finally, making use of the multi-element prediction formula, we show that the forecasting performance is enhanced by the joint use of the information in the location of the first event, and that in the subsequent interevent history in the cluster.     

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

A Temporary Network For Seismological Monitoring In West-Saxony: First Results

Since December 1997 a network of 10 autonomously recording digital seismic stations operates in W-Saxony to monitor microseismic activity generated by tectonic movements in the Gera-Jáchymov fault zone and induced by the flooding of the uranium mining area near Aue. Further goals to be followed with these data are the in situ verification of engineering seismological parameters and the study of the transmission properties of seismic waves across the Gera-Jáchymov fault zone to improve seismic risk evaluation. So far 19 tectonic and more than 50 induced events have been located and source parameters have been determined. Event clusters near Zwickau and within the Aue mining area could be spatially resolved with a master event technique. The cluster near Zwickau occured on a common fault plane striking parallel to the Gera-Jáchymov fault zone and dipping at an angle of about 60° towards ENE. The induced events near Aue occured within a volume of approximately 5 km ³ . The majority of these events were located in the granite body underneath the mining area. Fault plane solutions show a variety of focal mechanisms with no clear relation to the regional stress field.