We have developed an efficient method to determine high-resolution hypocenter locations over large distances. The location method incorporates ordinary absolute travel-time measurements and/or cross-correlation P-and S-wave differential travel-time measurements. Residuals between observed and theoretical travel-time differences (or double-differences) are minimized for pairs of earthquakes at each station while linking together all observed event-station pairs. A least-squares solu- tion is found by iteratively adjusting the vector difference between hypocentral pairs. The double-difference algorithm minimizes errors due to unmodeled velocity struc- ture without the use of station corrections. Because catalog and cross-correlation data are combined into one system of equations, interevent distances within multiplets are determined to the accuracy of the cross-correlation data, while the relative lo- cations between multiplets and uncorrelated events are simultaneously determined to the accuracy of the absolute travel-time data. Statistical resampling methods are used to estimate data accuracy and location errors. Uncertainties in double-difference locations are improved by more than an order of magnitude compared to catalog locations. The algorithm is tested, and its performance is demonstrated on two clus- ters of earthquakes located on the northern Hayward fault, California. There it col- lapses the diffuse catalog locations into sharp images of seismicity and reveals hor- izontal lineations of hypocenters that define the narrow regions on the fault where stress is released by brittle failure.

/pdf/a-double-difference-earthquake-location-algorithm-method-and-13ttz2fahg.pdf

A Double-Difference Earthquake Location Algorithm: Method and Application to the Northern Hayward Fault, California

The inverse problem of three-dimensional (3-D) local earthquake tomography is formulated as a linear approximation to a nonlinear function. Thus the solutions obtained and the reliability estimates depend on the initial reference model. Inappropriate models may result in artifacts of significant amplitude. Here, we advocate the application of the same inversion formalism to determine hypocenters and one-dimensional (1-D) velocity model parameters, including station corrections, as the first step in the 3-D modeling process. We call the resulting velocity model the minimum 1-D model. For test purposes, a synthetic data set based on the velocity structure of the San Andreas fault zone in central California was constructed. Two sets of 3-D tomographic P velocity results were calculated with identical travel time data and identical inversion parameters. One used an initial 1-D model selected from a priori knowledge of average crustal velocities, and the other used the minimum 1-D model. Where the data well resolve the structure, the 3-D image obtained with the minimum 1-D model is much closer to the true model than the one obtained with the a priori reference model. In zones of poor resolution, there are fewer artifacts in the 3-D image based on the minimum 1-D model. Although major characteristics of the 3-D velocity structure are present in both images, proper interpretation of the results obtained with the a priori 1-D model is seriously compromised by artifacts that distort the image and that go undetected by either resolution or covariance diagnostics.

/pdf/initial-reference-models-in-local-earthquake-tomography-1uwqgwr2xs.pdf

Initial reference models in local earthquake tomography

[1] Since their development in the 1960s, seismic arrays have given a new impulse to seismology. Recordings from many uniform seismometers in a well-defined, closely spaced configuration produce high-quality and homogeneous data sets, which can be used to study the Earth's structure in great detail. Apart from an improvement of the signal-to-noise ratio due to the simple summation of the individual array recordings, seismological arrays can be used in many different ways to study the fine-scale structure of the Earth's interior. They have helped to study such different structures as the interior of volcanos, continental crust and lithosphere, global variations of seismic velocities in the mantle, the core-mantle boundary and the structure of the inner core. For this purpose many different, specialized array techniques have been developed and applied to an increasing number of high-quality array data sets. Most array methods use the ability of seismic arrays to measure the vector velocity of an incident wave front, i.e., slowness and back azimuth. This information can be used to distinguish between different seismic phases, separate waves from different seismic events and improve the signal-to-noise ratio by stacking with respect to the varying slowness of different phases. The vector velocity information of scattered or reflected phases can be used to determine the region of the Earth from whence the seismic energy comes and with what structures it interacted. Therefore seismic arrays are perfectly suited to study the small-scale structure and variations of the material properties of the Earth. In this review we will give an introduction to various array techniques which have been developed since the 1960s. For each of these array techniques we give the basic mathematical equations and show examples of applications. The advantages and disadvantages and the appropriate applications and restrictions of the techniques will also be discussed. The main methods discussed are the beam-forming method, which forms the basis for several other methods, different slant stacking techniques, and frequency–wave number analysis. Finally, some methods used in exploration geophysics that have been adopted for global seismology are introduced. This is followed by a description of temporary and permanent arrays installed in the past, as well as existing arrays and seismic networks. We highlight their purposes and discuss briefly the advantages and disadvantages of different array configurations.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2000RG000100

Array seismology: Methods and applications

We have developed a double-difference (DD) seismic tomography method that makes use of both absolute and relative arrival times. By reducing systematic errors using the more accurate relative arrival times, the method produces an improved velocity model. Simultaneously, it yields event locations of a quality equivalent to those of the DD location method, while avoiding simplifying assumptions of that method. We test this method on a synthetic dataset and find that it produces a more accurate velocity model and event locations than standard tomography. We also test this method on a Hayward fault, California, earthquake dataset spanning 1984–1998. The earthquakes relocated by this method collapse to a thin line along the fault trace, consistent with previous results. The DD velocity model has sharper velocity contrasts near the source region than the standard tomography model.

Double-Difference Tomography: The Method and Its Application to the Hayward Fault, California

Unconventional oil and gas production provides a rapidly growing energy source; however, high-production states in the United States, such as Oklahoma, face sharply rising numbers of earthquakes. Subsurface pressure data required to unequivocally link earthquakes to wastewater injection are rarely accessible. Here we use seismicity and hydrogeological models to show that fluid migration from high-rate disposal wells in Oklahoma is potentially responsible for the largest swarm. Earthquake hypocenters occur within disposal formations and upper basement, between 2- and 5-kilometer depth. The modeled fluid pressure perturbation propagates throughout the same depth range and tracks earthquakes to distances of 35 kilometers, with a triggering threshold of ~0.07 megapascals. Although thousands of disposal wells operate aseismically, four of the highest-rate wells are capable of inducing 20% of 2008 to 2013 central U.S. seismicity.

https://science.sciencemag.org/content/sci/345/6195/448.full.pdf

Sharp increase in central Oklahoma seismicity since 2008 induced by massive wastewater injection

Using the database provided by the Reviewed Event Bulletins (REBs) for the first 2.5 years of the Group of Scientific Experts Technical Test-3 (GSETT-3) experiment, we compiled mislocation vectors for both arrays and selected three-component stations of the primary network from the published slowness and azimuth information gained through f-k- and polarization analysis. Imposing constraints such as a minimum signal-to-noise ratio (SNR) and number of defining phases, we aim at eliminating location bias as the hypocentral parameters are taken from the REBs. Results from 14 arrays with apertures from about 1 km to more than 20 km are presented as well as from 18 three-component stations, which indicate that the mislocation vectors in many cases can improve location accuracy considerably. If these mislocation vectors are compiled to provide coverage of a sufficient portion of the slowness domain these empirical corrections can easily be applied prior to location processing. In the context of the Comprehensive Nuclear Test-Ban Treaty (CTBT), these mislocation patterns could be essential for providing accurate event location of suspicious low-magnitude events, as these location parameters will be used to pinpoint the area where to conduct an on-site inspection.

Determination of Mislocation Vectors to Evaluate bias at GSETT-3 primary stations

This report of the Swiss Seismological Service summarizes the seismic activity in Switzerland and surrounding regions during 2010. During this period, 407 earthquakes and 85 quarry blasts were detected and located in the region under consideration. With a total of only 19 events with ML ≥ 2.5, the seismic activity in the year 2010 was below the average over the previous 35 years. The two most noteworthy earthquakes were the ML 3.4 Barrhorn event near Sankt Niklaus (VS) and the ML 3.0 event of Feldkirch, both of which produced shaking of intensity IV.

/pdf/earthquakes-in-switzerland-and-surrounding-regions-during-1fzv5991pi.pdf

Earthquakes in Switzerland and surrounding regions during 2010

This report of the Swiss Seismological Service summarizes the seismic activity in Switzerland and surrounding regions during 2008. During this period, 451 earthquakes and 75 quarry blasts were detected and located in the region under consideration. The three strongest events occurred in the Valais, near Lac des Toules (M
 L
 3.6), and in Graubunden, near Ilanz (M
 L
 3.7) and Paspels (M
 L
 4.0). Although felt by the population, they were not reported to have caused any damage. However, with a total of only 15 events with M
 L
 ≥ 2.5, the seismic activity in the year 2008 was far below the average over the previous 33 years.

https://link.springer.com/content/pdf/10.1007%2Fs00015-009-1339-8.pdf

Earthquakes in Switzerland and surrounding regions during 2008

Regional Moment-Tensor Inversion in the European-Mediterranean

It's time to revisit one of the more successful and popular seismic data access systems available on the Internet. The Automatic Data Request Manager (AutoDRM) is approaching ten years old but is going strong, providing access to ever-increasing amounts of seismic data for more and more users. A fairly complete description of the development of the AutoDRM was one of the first guest columns of the Electronic Seismologist back in the summer of 1996. The system has undergone relatively few changes and enhancements since then, attesting to its original good design and thoughtful execution. If anything, its success is also one of its current problems; there are so many data available via the AutoDRM mechanism that it is hard to keep up with what is available when and from where. To help with this problem the father of AutoDRM is back at it, having now produced a system for summarizing data availability from the myriad of AutoDRM's and providing this summary via a simple set of Web pages and, yes, via an AutoDRM e-mail request method as well. Urs Kradolfer is again our guest columnist to explain this new service.

Seismograms are the key data for many seismological researchers, and since seismic waveforms are now widely recorded and archived in digital form, global communications facilities, such as the Internet, potentially allow the seismological community to access lots of data. An important step toward making seismic signals available for research was the foundation of international data centers such as those at IRIS and ORFEUS. These centers collect and provide waveforms for significant earthquakes originating from many stations worldwide.

Mainly in order to allow access to near-real-time waveform data, but also to facilitate international data exchange, many observatories have implemented an AutoDRM, which is an e-mail-based data exchange system (Kradolfer 1993, …

Urs Kradolfer

Papers

Determination of Mislocation Vectors to Evaluate bias at GSETT-3 primary stations

Earthquakes in Switzerland and surrounding regions during 2010

Earthquakes in Switzerland and surrounding regions during 2008

Regional Moment-Tensor Inversion in the European-Mediterranean

Waves4U: Waveform Availability Through AutoDRM's