Source Specific Station Corrections for Regional Phases at Fennoscandian Stations

—?The IASPEI91 global travel-time curves are used as the default for event location at the Prototype International Data Center (PIDC). In order to improve event location, a 1-D Baltic travel-time model was implemented at the PIDC in 1997 for locating events using regional phases from Fennoscandian stations. Where a single model is insufficient for characterizing the regional geology, path-dependent corrections, or Source Specific Station Corrections (SSSCs), are more appropriate for event locations. We have developed SSSCs for regional phases at the Fennoscandian stations by interpolating travel times through different 1-D models. SSSCs for stations NRIS and SPITS are also derived, given the fact that paths from both stations to high latitude events are within the Fennoscandia regionalization as Baltic.¶Validation testing of the SSSCs demonstrates that using SSSCs in event location is superior to not using SSSCs, a nd, in most cases, to using the 1-D model directly when locating events. For a ground-truth data set which includes events in the Baltic Shield with location accuracy better than 2?km, the average improvement in location due to SSSCs is 9?km, and the median coverage ellipse is reduced by 2710?km² (from 3830 to 1120?km²). These results are similar to those obtained using the 1-D Baltic model. For a CEB (Calibration Event Bulletin) data set which includes events along the North Atlantic oceanic ridge and in central/southern Europe, using SSSCs the ridge events move closer to the ridge axis, and the European events move closer to CEB locations than 1-D Baltic locations. For a constrained JHD (Joint Hypocenter Determination) data set of events in the Novaya Zemlya region, when using SSSCs or the 1-D Baltic model, relative to the JHD locations mislocations are less or similar to those without SSSCs. All coverage ellipses are smaller but sti ll contain the JHD solutions.¶Our SSSCs are strongly dependent on the 1-D regional models and regionalization. Future development in 1-D velocity models and travel-time curves should improve such SSSCs, event locations, and uncertainties. It is hoped that the implementation and demonstration of SSSCs in the PIDC software will encourage these further developments. These SSSCs were implemented at the PIDC for Reviewed Event Bulletin (REB) location in April 1999.     

Improving Regional Seismic Event Location in China

—?In an effort to improve our ability to locate seismic events in China using only regional data, we have developed empirical propagation path corrections and applied such corrections using traditional location routines. Thus far, we have concentrated on corrections to observed P arrival times for crustal events using travel-time observations available from the USGS Earthquake Data Reports, the International Seismic Centre Bulletin, the preliminary International Data Center Reviewed Event Bulletin, and our own travel-time picks from regional data. Location ground truth for events used in this study ranges from 25?km for well-located teleseimic events, down to 2?km for nuclear explosions located using satellite imagery. We also use eight events for which depth is constrained using several waveform methods. We relocate events using the EvLoc algorithm from a region encompassing much of China (latitude 20°–55°N; longitude 65°–115°E). We observe that travel-time residuals exhibit a distance-dependent bias using IASPEI91 as our base model. To remedy this bias, we have developed a new 1-D model for China, which removes a significant portion of the distance bias. For individual stations having sufficient P-wave residual data, we produce a map of the regional travel-time residuals from all well-located teleseismic events. Residuals are used only if they are smaller than 10?s in absolute value and if the seismic event is located with accuracy better than 25?km. From the residual data, correction surfaces are constructed using modified Bayesian kriging. Modified Bayesian kriging offers us the advantage of providing well-behaved interpolants and their errors, but requires that we have adequate error estimates associated with the travel-time residuals from which they are constructed. For our P-wave residual error estimate, we use the sum of measurement and modeling errors, where measurement error is based on signal-to-noise ratios when available, and on the published catalog estimate otherwise. Our modeling error originates from the variance of travel-time residuals for our 1-D China model. We calculate propagation path correction surfaces for 74 stations in and around China, including six stations from the International Monitoring System. The statistical significance of each correction surface is evaluated using a cross-validation technique. We show relocation results for nuclear tests from the Balapan and Lop Nor test sites, and for earthquakes located using interferometric synthetic aperture radar. These examples show that the use of propagation path correction surfaces in regional relocations eliminates distance bias in the residual curves and significantly improves the accuracy and precision of seismic event locations.     

Application of 3-D Crustal and Upper Mantle Velocity Model of North America for Location of Regional Seismic Events

—?Seismic event locations based on regional 1-D velocity-depth sections can have bias errors caused by travel-time variations within different tectonic provinces and due to ray-paths crossing boundaries between tectonic provinces with different crustal and upper mantle velocity structures. Seismic event locations based on 3-D velocity models have the potential to overcome these limitations. This paper summarizes preliminary results for calibration of IMS for North America using 3-D velocity model. A 3-D modeling software was used to compute Source-Station Specific Corrections (SSSCs(3-D)) for Pn travel times utilizing 3-D crustal and upper mantle velocity model for the region. This research was performed within the framework of the United States/Russian Federation Joint Program of Seismic Calibration of the International Monitoring System (IMS) in Northern Eurasia and North America.¶An initial 3-D velocity model for North America was derived by combining and interpolating 1-D velocity-depth sections for different tectonic units. In areas where no information on 1-D velocity-depth sections was available, tectonic regionalization was used to extrapolate or interpolate. A Moho depth map was integrated. This approach combines the information obtained from refraction profiles with information derived from local and regional network data. The initial 3-D velocity model was tested against maps of Pn travel-time residuals for eight calibration explosions; corrections to the 3-D model were made to fit the observed residuals. Our goal was to find a 3-D crustal and upper mantle velocity model capable predicting Pn travel times with an accuracy of 1.0–1.5 seconds (r.m.s.).¶The 3-D velocity model for North America that gave the best fit to the observed travel times, was used to produce maps of SSSCs(3-D) for seismic stations. The computed SSSCs(3-D) vary approximately from +5 seconds to ?5 seconds for the western USA and the Pre-Cambrian platform, respectively. These SSSCs(3-D) along with estimated modeling and measurement errors were used to relocate, using regional data, an independent set of large chemical explosions (with known locations and origin times) detonated within various tectonic provinces of North America. Utilization of the 3-D velocity model through application of the computed SSSCs(3-D) resulted in a substantial improvement in seismic event location accuracy and in a significant decrease of error ellipse area for all events analyzed in comparison both with locations based on the IASPEI91 travel times and locations based on 1-D regional velocity models.     

A Study of Small Magnitude Seismic Events During 1961–1989 on and near the Semipalatinsk Test Site, Kazakhstan

—?Official Russian sources in 1996 and 1997 have stated that 340 underground nuclear tests (UNTs) were conducted during 1961–1989 at the Semipalatinsk Test Site (STS) in Eastern Kazakhstan. Only 271 of these nuclear tests appear to have been described with well-determined origin time, coordinates and magnitudes in the openly available technical literature. Thus, good open documentation has been lacking for 69 UNTs at STS.¶The main goal of our study was to provide detections, estimates of origin time and location, and magnitudes, for as many of these previously undocumented events as possible. We used data from temporary and permanent seismographic stations in the former USSR at distances from 500?km to about 1500?km from STS. As a result, we have been able to assign magnitude for eight previously located UNTs whose magnitude was not previously known. For 31 UNTs, we have estimated origin time an d assigned magnitude — and for 19 of these 31 we have obtained locations based on seismic signals. Of the remaining 30 poorly documented UNTs, 15 had announced yields that were less than one ton, and 13 occurred simultaneously with another test which was detected. There are only two UNTs, for which the announced yield exceeds one ton and we have been unable to find seismic signals.¶Most of the newly detected and located events were sub-kiloton. Their magnitudes range from 2.7 up to 5.1 (a multi-kiloton event on 1965 Feb. 4 that was often obscured at teleseismic stations by signals from an earthquake swarm in the Aleutians).¶For 17 small UNTs at STS, we compare the locations (with their uncertainties) that we had earlier determined in 1994 from analysis of regional seismic waves, with ground-truth information obtained in 1998. The average error of the seismically-determined locations is only about 5?km. The ground-truth location is almost alw ays within the predicted small uncertainty of the seismically-determined location.¶Seismically-determined yield estimates are in good agreement with the announced total annual yield of nuclear tests, for each year from 1964 to 1989 at Semipalatinsk.¶We also report the origin time, location, and seismic magnitude of 29 chemical explosions and a few earthquakes on or near STS during the years 1961–1989.¶Our new documentation of STS explosions is important for evaluating the detection, location, and identification capabilities of teleseismic and regional arrays and stations; and how these capabilities have changed with time.     

Location Calibration Data for CTBT Monitoring at the Prototype International Data Center

—?Ground-truth information is essential for location calibration of the International Monitoring System network being developed under the Comprehensive Nuclear-Test-Ban Treaty. The objective of the calibration effort is to improve the accuracy of seismic event locations and to reduce the size of the error ellipse, both in automatic and in human analyst-reviewed bulletins, in order to meet the On-Site Inspection requirement for the size of the inspection area. Several databases were compiled and are maintained at the Prototype International Data Center (PIDC) to support calibration efforts. The Nuclear Explosion Database contains most readily accessible information about all nuclear explosions worldwide. The events in the Calibration Event Bulletin (CEB) carefully selected well located events from the PIDC Reviewed Event Bulletin and relocated using additional arrivals from regional networks requested from various National Data Centers. The Ground-Truth Database contains carefully selected events with known or well estimated location accuracies from the Nuclear Explosion Database, CEB, as well as from bulletins of U.S. National Earthquake Information Center and International Seismic Centre. It also contains data on chemical explosions and quarry blasts when confirmed by local or national authorities. Ground-truth events are subdivided into various ground-truth categories according to their location accuracy. The databases have been used in various calibration studies to derive and test corrections to improve event locations. Several location calibration techniques are briefly described. The validation test for any proposed operational change requires that the results meet the location calibration metrics developed and implemented at the PIDC.     

Global Event Location with Full and Sparse Data Sets Using Three-dimensional Models of Mantle P-wave Velocity

—?In order to improve on the accuracy of event locations at teleseismic distances it is necessary to adequately correct for lateral variations in structure along the ray paths, either through deterministic model-based corrections, empirical path/station corrections, or a combination of both approaches. In this paper we investigate the ability of current three-dimensional models of mantle P-wave velocity to accurately locate teleseismic events. We test four recently published models; two are parameterized in terms of relatively long-wavelength spherical harmonic functions up to degree 12, and two are parameterized in terms of blocks of constant velocity which have a dimension of a few hundreds of km. These models, together with detailed crustal corrections, are used to locate a set of 112 global test events, consisting of both earthquakes and explosions with P-wave travel-time data compiled by the Internation al Seismological Centre (ISC). The results indicate that the supposedly higher resolution block models do not improve the accuracy of teleseismic event locations over the longer wavelength spherical harmonic models. For some source locations the block models do not predict the range of observed travel-time residuals as well as the longer wavelength models. The accuracy of the locations largely varies randomly with geographic position although events in central Asia are particularly well located. We also tested the effect of reduced data sets on the locations. Multiple location iterations using 30 P-wave travel times indicate that teleseismic events may be located within an area of 1000?km² of the true location 66% of the time with only the model-based corrections, and increasing to 75% if calibration information is available. If as few as 8 phases are available then this is possible only 50% of the time. Further refinement in models and/or procedure, such as the addition of P _n phases, azimuth data, and consideration of P-wave anisotropy may provide further improvement in the teleseismic location of small events.     

Operational Processing of Hydroacoustics at the Prototype International Data Center

— The Prototype International Data Center (PIDC) has designed and implemented a system to process data from the International Monitoring System's hydroacoustic network. The automatic system detects and measures various signal characteristics that are then used to classify the signal into one of three categories. The detected signals are combined with the seismic and infrasonic detections to automatically form event hypotheses. The automatic results are reviewed by human analysts to form the Reviewed Event Bulletin (REB). Continuous processing of hydroacoustic data has been in place since May 1997 and during that time a large database of hydroacoustic signals has been accumulated. For a two-year period, the REB contains 13,582 T phases that are associated to 8,437 events. This is roughly 25% of REB events after taking station downtime into account. Predicted travel times used in locations are based on the arrival time of the peak a mplitude mode calculated from a normal mode propagation model. Global sound velocity and bathymetry databases are used to obtain reliable 2-D, seasonally dependent, travel-time tables for each hydroacoustic station in the PIDC. A limited number of ground-truth observations indicate that the predicted travel times are good to within 5 seconds for paths extending to over 7,000?km – corresponding to a relative error of less than 0.1%. The ground truth indicates that the random errors in measuring arrival times for impulsive signals are between 1 and 6 seconds. This paper describes and evaluates the automatic hydroacoustic processing compared to the analyst reviewed results. In addition, special studies help characterize the overall performance of the hydroacoustic network.     

Best Practice in Earthquake Location Using Broadband Three-component Seismic Waveform Data

—?We present an earthquake location algorithm, the Broadband Waveform Regional Earthquake Location Program (BW_RELP), which utilizes phase onset times and wave azimuths recorded by three-component broadband seismic stations and an adaptive migrating grid search algorithm to find the global minimum in an arbitrary normed misfit parameter. The performance of BW_RELP is demonstrated using regional (300–800?km distant) broadband recordings to locate events in the 1995 Ridgecrest, California earthquake sequence. The purpose of this study is to introduce the BW_RELP algorithm in detail and to expand on the previous paper by Deger et?al. (BSSA, 88, 1353–1362, 1998), using one Berkeley Digital Seismic Network (BDSN) station (YBH) and two USNSN stations (ELK and MNV) which span 300–800?km in distance and 55 degrees in azimuth, to further investigate the capability of a sparse broadband network of three-component stations at monitoring a region located outside of the network, as will be the case in the monitoring of the Comprehensive Test-Ban-Treaty (CTBT) for low magnitude seismic events. We assess the capability of this sparse three-station broadband network and we compare locations estimated from phase onset time and wave azimuth measurements to a ground-truth catalog of high-quality earthquake locations derived from data recorded by the Southern California Seismic Network (SCSN). The results indicate that in the regional distance range it is possible, when an appropriate calibration event is available, to obtain absolute event locations to within 18?km as is prescribed by the CTBT.     

Analysis of the IMS Location Accuracy in Northern Eurasia and North America Using Regional and Global Pn Travel-time Tables

—?Joint Research Program of Seismic Calibration of the International Monitoring System (IMS) in Northern Eurasia and North America has been signed by the Nuclear Treaty Programs Office (NTPO), Department of Defense USA, and the Special Monitoring Service (SMS) of the Ministry of Defense, Russian Federation (RF). Under the Program historical data from nuclear and large chemical explosions of known location and shot time, together with appropriate geological and geophysical data, has been used to derive regional Pn/P travel-time tables for seismic event location within the lower 48 States of the USA and the European part of the RF. These travel-time tables are up to 5?seconds faster in shields than the IASPEI91 tables, and up to 5?seconds slower in the Western USA. Relocation experiments using the regional Pn travel-time curves and surrogate networks for the IMS network generally improved locations for regional seismic events. The distance between true and estimated location (mislocation) was decreased from an average of 18.8?km for the IASPEI91 tables to 10.1?km for the regional Pn travel-time tables. However, the regional travel-time table approach has limitations caused by travel-time variations inside major tectonic provinces and paths crossing several tectonic provinces with substantially different crustal and upper mantle velocity structure.¶The RF members of the Calibration Working Group (WG): Colonel Vyacheslav Gordon (chairman); Dr. Prof. Marat Mamsurov, and Dr. Nikolai Vasiliev. The US members of the WG: Dr. Anton Dainty (chairman), Dr. Douglas Baumgardt, Mr. John Murphy, Dr. Robert North, and Dr. Vladislav Ryaboy.     

Calibration of the Regional Crustal Waveguide and the Retrieval of Source Parameters Using Waveform Modeling

v--vRegional crustal waveguide calibration is essential to the retrieval of source parameters and the location of smaller (M < 4.8) seismic events. This path calibration of regional seismic phases is strongly dependent on the accuracy of hypocentral locations of calibration (or master) events. This information can be difficult to obtain, especially for smaller events. Generally, explosion or quarry blast generated travel-time data with known locations and origin times are useful for developing the path calibration parameters, but in many regions such data sets are scanty or do not exist. We present a method which is useful for regional path calibration independent of such data, i.e. with earthquakes, which is applicable for events down to M_w = 4 and which has successfully been applied in India, central Asia, western Mediterranean, North Africa, Tibet and the former Soviet Union. These studies suggest that reliably determining depth is essential to establishing accurate epicentral location and origin time for events. We find that the error in source depth does not necessarily trade-off only with the origin time for events with poor azimuthal coverage, but with the horizontal location as well, thus resulting in poor epicentral locations. For example, hypocenters for some events in central Asia were found to move from their fixed-depth locations by about 20 km. Such errors in location and depth will propagate into path calibration parameters, particularly with respect to travel times. The modeling of teleseismic depth phases (pP, sP) yields accurate depths for earthquakes down to magnitude M_w = 4.7. This M_w threshold can be lowered to four if regional seismograms are used in conjunction with a calibrated velocity structure model to determine depth, with the relative amplitude of the P_nl waves to the surface waves and the interaction of regional sPmP and pPmP phases being good indicators of event depths. We also found that for deep events a seismic phase which follows an S-wave path to the surface and becomes critical, developing a head wave by S to P conversion is also indicative of depth. The detailed characteristic of this phase is controlled by the crustal waveguide. The key to calibrating regionalized crustal velocity structure is to determine depths for a set of master events by applying the above methods and then by modeling characteristic features that are recorded on the regional waveforms. The regionalization scheme can also incorporate mixed-path crustal waveguide models for cases in which seismic waves traverse two or more distinctly different crustal structures. We also demonstrate that once depths are established, we need only two-stations travel-time data to obtain reliable epicentral locations using a new adaptive grid-search technique which yields locations similar to those determined using travel-time data from local seismic networks with better azimuthal coverage.     

基于克里金方法的区域震相走时标定技术研究div.

在台网比较稀疏的情况下,台站走时标定是提高低震级事件定位、识别能力的重要手段．为了提高稀疏台网的定位精度,首先利用标定事件和IASPEI1991走时表以及中国大陆走时表,计算了初至P波和Lg波到AAK,MAKZ,NIL,TLG以及WMQ等5个台站的走时残差;接着,采用非平稳贝叶斯克里金方法和走时残差数据构造上述台站的走时校正曲面;最后,通过加入和不加入走时校正定位一系列位置准确的发生在新疆地区的地震事件评估克里金走时校正的有效性．结果表明,克里金走时校正能够较大地提高稀疏台网的定位精度,同时有效地缩小误差椭圆的面积．      

Signal Processing for Indian and Pakistan Nuclear Tests Recorded at IMS Stations Located in Israel

v—vIn compliance with the Comprehensive Nuclear-Test-Ban-Treaty (CTBT) the International Monitoring System (IMS) was designed for detection and location of the clandestine Nuclear Tests (NT). Two auxiliary IMS seismic stations MRNI and EIL, deployed recently, were subjected to detectability, travel-time calibration and discrimination analysis. The study is based on the three recent 1998 underground nuclear explosions: one of India and two of Pakistan, which provided a ground-truth test of the existing IMS. These events, attaining magnitudes of 5.2, 4.8 and 4.6 correspondingly, were registered by many IMS and other seismic stations.¶The MRNI and EIL broadband (BB) stations are located in Israel at teleseismic distances (from the explosions) of 3600, 2800 and 2700ukm, respectively, where the signals from the tests are already weak. The Indian and the second Pakistan NT were not detected by the short-period Israel Seismic Network (ISN), using standard STA/LTA triggering. Therefore, for the chosen IMS stations we compare the STA/LTA response to the results of the more sensitive Murdock-Hutt (MH) and the Adaptive Statistically Optimal Detector (OD) that showed triggering for these three events. The second Pakistan NT signal arrived at the ISN and the IMS stations in the coda of a strong Afghanistan earthquake and was further disturbed by a preceding signal from a local earthquake. However, the NT signal was successfully extracted at EIL and MRNI stations using MH and OD procedures. For comparison we provide the signal analysis of the cooperating BB station JER, with considerably worse noise conditions than EIL and MRNI, and show that OD can detect events when the other algorithms fail. Using the most quiet EIL station, the most sensitive OD and different bandpass filters we tried in addition to detect the small Kazakh chemical 100-ton calibration explosion of 1998, with magnitude 3.7 at a distance approaching 4000ukm. The detector response curve showed uprising in the expected signal time interval, but yet was low for a reliable decision.¶After an NT is detected it should be recognized. Spectra were calculated in a 15-sec window including P and P-coda waves. The spectra for the first Pakistan NT showed a pronounced spectral null at 1.7uHz for all three components of the EIL station. The effect was confirmed by observation of the same spectral null at the vertical component of the ISN stations. For this ground-truth explosion with a reported shallow source depth, the phenomenon can be explained in terms of the interference of P and pP phases. However, the spectral null feature, considered separately, cannot serve as a reliable identification characteristic of nuclear explosions, because not all the tests provide the nulls, whereas some earthquakes show this feature. Therefore, the multi-channel spectral discrimination analysis, based on a spectral ratio of low-to-high frequency energy (in the 0.6–1uHz and 1–3uHz bands), and a semblance of spectral curves (in the 0.6–2uHz band), was conducted. Both statistics were calculated for the vertical component of the ISN stations as well for the three components of the EIL station. The statistics provided a reliable discrimination between the recent NT and several nearby earthquakes, and showed compliance with the former analysis of Soviet and Chinese NT, where nuclear tests demonstrated lower values of energy ratio and spectral semblance than earthquakes. ¶Accurate location of NT requires calibration of travel time for IMS stations. Using known source locations, IASPEI91 travel-time tables and NEIC origin times we calculated expected arrival time for the P waves to the EIL and MRNI stations and showed that the measured arrival time has a delay of about 4 sec. Similar results were obtained for the nearby Pakistan earthquakes. The analysis was complimented by the P travel-time measurements for the set of Semipalatinsk NT, which showed delays of about 3.7usec to the short-period MBH station which is a surrogate station for EIL. Similar delays at different stations evidence a path- rather than site-effect. The results can be used for calibration of the IMS stations EIL and MRNI regarding Asian seismic events.     

The Kirovskiy Explosion of September 29, 1996: Example of a CTB Event Notification for a Routine Mining Blast

—?On September 29, 1996, a routine mining blast of about 390 metric tons was detonated underground at the Kirovskiy mine in the central Kola Peninsula. The United States was notified two weeks in advance that the blast was to take place and was given the date, approximate time, location and total charge. The explosion was detected and located by the prototype International Data Center (pIDC) and published in the Reviewed Event Bulletin (REB). Detailed information about the blast, including the type and depth of mining operation, the underground charge configuration, and the blasting delay pattern, is reviewed and combined with a seismological analysis of the event. The seismic analysis points to a possible associated tectonic component to the blast, consisting of a small rock burst or induced tremor, spall, or some combination of these mechanisms, that may have enhanced the shear waves, produced large Rg waves at low frequency, and small Pn/Sn and Pn/Lg amplitude ratios at high frequency. While these discriminants might identify the event as an earthquake, the spectral/cepstral analysis of the event clearly shows the ripple-fire delays. This event provides important confidence-building measures for both location calibration, in the form of travel-time corrections for location of mine events in this region, and for improved understanding of seismic discriminants expected for large mine blasts that may have an associated induced tectonic component (e.g., spall, mine tremor or rock burst).     

Instrumental Seismicity of the Western Alps: A Revised Catalogue

—The western Alpine regions have been instrumented since the beginning of the century, and the number of seismological stations largely increased since 1980. This dense network has allowed an important improvement in the hypocentral determination, even for low magnitude events. This condition was a good opportunity to perform a synthesis of 32 years of instrumental seismicity in the Western Alps and southeast of France (1962–1993) and to attempt an improvement of the older event location with the assistance of the more recent locations.¶The revised catalogue of seismicity is built using station corrections and regional crustal models. After the elimination of non-natural events, the catalogue is composed of 6697 events. Another improvement corresponds to the revision of magnitudes. We performed several tests to evaluate the reliability of our results location of quarry events and rock bursts, epicentral correlation with geological features, coherence in depth with interpreted seismic profile (ECORS line), Moho isobaths. A first use of this catalogue is presented for the Haute-Ubaye region in the southwestern Alps.     

Ground truth seismic events and location capability at Degelen mountain, Kazakhstan

We utilized nuclear explosions from the Degelen Mountain sub-region of the Semipalatinsk Test Site (STS), Kazakhstan, to assess seismic location capability directly. Excellent ground truth information for these events was either known or was estimated from maps of the Degelen Mountain adit complex. Origin times were refined for events for which absolute origin time information was unknown using catalog arrival times, our ground truth location estimates, and a time baseline provided by fixing known origin times during a joint hypocenter determination (JHD). Precise arrival time picks were determined using a waveform cross-correlation process applied to the available digital data. These data were used in a JHD analysis. We found that very accurate locations were possible when high precision, waveform cross-correlation arrival times were combined with JHD. Relocation with our full digital data set resulted in a mean mislocation of 2 km and a mean 95% confidence ellipse (CE) area of 6.6 km² (90% CE: 5.1 km²), however, only 5 of the 18 computed error ellipses actually covered the associated ground truth location estimate. To test a more realistic nuclear test monitoring scenario, we applied our JHD analysis to a set of seven events (one fixed) using data only from seismic stations within 40° epicentral distance. Relocation with these data resulted in a mean mislocation of 7.4 km, with four of the 95% error ellipses covering less than 570 km² (90% CE: 438 km²), and the other two covering 1730 and 8869 km² (90% CE: 1331 and 6822 km²). Location uncertainties calculated using JHD often underestimated the true error, but a circular region with a radius equal to the mislocation covered less than 1000 km² for all events having more than three observations.     

Locating Seismic Events in the CTBT Context

—?The verification of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) requires the determination of accurate location of seismic events from a fixed network of seismic stations across the globe. The requirements of possible on-site inspections mean that the goal is to place the location estimate in a zone smaller than 1000 km² that includes the true location. Because a defined set of stations will be used, corrections can be refined to represent the influence of departures from the global reference model IASPEI91. The primary stations in the International Monitoring Scheme (IMS) are mostly seismic arrays and therefore the present location scheme is based on minimisation of a misfit function built from arrival time, azimuth and array slowness residuals. The effective network will change markedly with the magnitude of the event and as a result regional information has to be integrated into the location process.     

Locating Small Aftershocks Using a Small-Aperture Temporary Seismic Array

In this paper, we developed a specialized method to locate small aftershocks using a small-aperture temporary seismic array. The array location technique uses the first P arrival times to determine the horizontal slowness vector of the incoming P wave, then combines it with S–P times to determine the event location. In order to reduce the influence of lateral velocity variation on the location determinations, we generated slowness corrections using events well-located by the permanent broadband network as calibration events, then we applied the corrections to the estimated slownesses. Applications of slowness corrections significantly improved event locations. This method can be a useful tool to locate events recorded by temporary fault-zone arrays in the near field but unlocated by the regional permanent seismic network. As a test, we first applied this method to 64 well-located aftershocks of the 1992 Landers, California, earthquake, recorded by both the Caltech/USGS Southern California Seismic Network and a small-aperture, temporary seismic array. The average horizontal and vertical separations between our locations and the well-determined catalogue locations are 1.35 and 1.75 km, respectively. We then applied this method to 132 unlocated aftershocks recorded only by the temporary seismic array. The locations show a clear tendency to follow the surface traces of the mainshock rupture.     

Scoring IDC SEL1 events for identifying misformed events

The International Data Center (IDC) produces several automatic seismic event lists followed by a reviewed event list as part of the verification regime for the Comprehensive Nuclear-Test-Ban Treaty. For the time period from 2005 to 2011, on average, approximately 160 events per day were included in the first IDC automatic Standard Event List (SEL1). For the same time period, “only” 86 events per day were included in the IDC Reviewed Event Bulletin (REB). Several possible reasons may explain why an event included in SEL1 is excluded from the REB; for example: misassociation, namely, the association process generated an artificial event by associating information from stations that recorded several different events. Reducing the number of SEL1 events that are ultimately excluded from the REB would improve the automatic process relying on the SEL1 Bulletin and would reduce the workload of the analysts at the IDC. In this work, we present a scoring method that indicates the probability that an SEL1 event will “survive” the analysis process and be included in the REB. The method tries to imitate “good” analyst practice. Namely, based on the experience gained in analyzing information from a specific region with a seismic network, the analyst can decide if an event is a “real” event or not. The score was created for the International Monitoring System primary seismic stations based on the REB for the time period of 2005 to the beginning of 2011. The score is designed to indicate if the set of stations contributing/not contributing to the event is consistent with what is expected based on the analysis of past events that occurred in the same region. A byproduct of the analysis needed for the score is a statistical measure that assesses the typical contribution of a specific station to a specific region. This information can be useful in its own right, for example, to design subnetworks that will increase automatic processing efficiency.     

Development and Applications of Double-difference Seismic Tomography

Double-difference (DD) tomography is a generalization of DD location; it simultaneously solves for the three-dimensional velocity structure and seismic event locations. DD tomography uses a combination of absolute and more accurate differential arrival times and hierarchically determines the velocity structure from larger scale to smaller scale. This method is able to produce more accurate event locations and velocity structure near the source region than standard tomography, which uses only absolute arrival times. We conduct a stability and uncertainty analysis of DD tomography based on a synthetic data set. Currently three versions of the DD tomography algorithms exist: tomoDD, tomoFDD and tomoADD. TomoDD assumes a flat earth model and uses a pseudo-bending ray-tracing algorithm to find rays between events and stations while tomoFDD uses a finite-difference travel-time algorithm and the curvature of the Earth is considered. Both codes are based on a regularly distributed inversion grid, with the former for a local scale and the latter for a regional scale. In contrast, tomoADD adapts the inversion mesh to match with the data distribution based on tetrahedral and Voronoi diagrams. We discuss examples of applying DD tomography to characterize fault zone structure, image high-resolution structure of subduction zones, and determine the velocity structure of volcanoes.