首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
v—vRegional seismic discriminants for the May 11, 1998 Indian underground nuclear test(s) and earthquakes recorded at station NIL (Nilore, Pakistan) provide new data to test strategies that can be used to monitor the Comprehensive Nuclear-Test-Ban Treaty (CTBT). Three categories of regional discriminants (ratios of P- and/or S-wave energy) were measured on short-period (0.5–6uHz) seismograms: P/S amplitude ratios (phase ratios) measured in the same frequency band, P- and S-wave spectral ratios (i.e., low frequency to high frequency for the same phase) and P/S cross-spectral ratios (i.e., low frequency S-wave to high frequency P-wave). The P/S amplitude ratios show good separation of the Indian nuclear test and regional earthquakes for Pn/Lg and Pn/Sn, however Pg/Lg does not discriminate as well. Pn/Lg and Pn/Sn discriminate well at frequencies as low as 0.5–2uHz, especially after accounting for path effects. This observation differs from previous studies that report poor separation of earthquakes and explosions at lower frequencies. The P/S amplitude ratios do not show any magnitude dependence, suggesting that forming the ratios in a fixed frequency band cancels the effects of source size-corner frequency scaling. Spatial variability of the observed discriminants arises from variations in crustal waveguide and/or attenuation structure (path propagation effects). Grouping amplitude ratios for earthquakes with paths similar to the Indian test greatly improves discrimination. Removing distance trends does not generally improve discrimination. Accounting for path effects with Bayesian kriging significantly improves discrimination. Spectral ratios (e.g., Pn [0.5–1uHz]/Pn [4–6uHz]) and cross-spectral ratios (e.g., Lg [1–2uHz]/Pn [4–6uHz]) show distance and magnitude dependence. We developed a technique for simultaneously removing the effects of distance and source size-corner frequency scaling on amplitude and spectral ratios. The technique uses a gridsearch to find several parameters that characterize the observed distance and magnitude dependence. Discrimination of the Indian test improved dramatically after the distance and magnitude trends were removed from the spectral and cross-spectral ratio data.  相似文献   

2.
The dynamic calibration method (DCM), using natural seismicity data and initially elaborated in [Kedrov, 2001; Kedrov et al., 2001; Kedrov and Kedrov, 2003], is applied to International Monitoring System (IMS) stations in Central Asia. The algorithm of the method is refined and a program is designed for calibrating diagnostic parameters (discriminants) that characterize a seismic source on the source-station traces. The DCM calibration of stations in relation to the region under study is performed by the choice of attenuation coefficients that adapt the diagnostic parameters to the conditions in a reference region. In this method, the stable Eurasia region is used as the latter. The calibration used numerical data samples taken from the archive of the International Data Centre (IDC) for the IMS stations MKAR, BVAR, EIL, ASF, and CMAR. In this paper, we used discriminants in the spectral and time domains that have the form
$D_i = X_i - a_m m_b - b_\Delta \log \Delta $
and are independent of the magnitude m b and the epicentral distance Δ; these discriminants were elaborated in [Kedrov et al., 1990; Kedrov and Lyuke, 1999] on the basis of a method used for identification of events at regional distances in Eurasia. Prerequisites of the DCM are the assumptions that the coefficient a m is regionindependent and the coefficient b Δ depends only on the geotectonic characteristics of the medium and does not depend on the source type. Thus, b Δ can be evaluated only from a sample of earthquakes in the region studied; it is used for adapting the discriminants D(X i ) in the region studied to the reference region. The algorithm is constructed in such a way that corrected values of D(X i) are calculated from the found values of the calibration coefficients b Δ, after which natural events in the region under study are selected by filtering. Empirical estimates of the filtering efficiency as a function of a station vary in a range of 95–100%. The DCM was independently tested using records obtained at the IRIS (Incorporated Research Institutions for Seismology) stations BRVK and MAKZ from explosions detonated in India on May 11, 1998, and Pakistan on May 28, 1998; these stations are similar in location and recording instrumentation characteristics to the IMS stations BVAR and MKAR. This test resulted in correct recognition of the source type and thereby directly confirmed the validity of the proposed calibration method of stations with the use of natural seismicity data. It is shown that the calibration coefficients b Δ for traces similar in the conditions of signal propagation (e.g., the traces from Iran to the stations EIL and ASF) are comparable for nearly all diagnostic parameters. We arrive at the conclusion that the method of dynamic calibration of stations using natural seismicity data in a region where no explosions were detonated can be significant for a rapid and inexpensive calibration of IMS stations. The DCM can also be used for recognition of industrial chemical explosions that are sometimes erroneously classified in regional catalogs as earthquakes.
  相似文献   

3.
朝鲜2009年和2013年两次核爆的地震学特征对比研究   总被引:4,自引:1,他引:3       下载免费PDF全文
利用布设在长白山地区临时地震台站接收到的朝鲜核爆的波形资料,对2009年5月25日和2013年2月12日两次朝鲜核试验的地震学特征进行比较.震中距范围从145km到420km.采用P/S型谱比值方法识别朝鲜核爆,通过与2009年3月20日长春地震和2013年1月23日沈阳地震事件的比较,表明在频率大于3 Hz时P/S型谱比值能够有效识别发生在中朝边境地区的地下核试验.选定参考台站,利用区域震相Pg波的振幅谱比值计算朝鲜核爆至各台站路径上的相对衰减.结合介质速度模型,在一定程度上反映了长白山地区衰减情况,为进一步研究长白山地区衰减层析成像提供初始模型.  相似文献   

4.
v--vS/P amplitude ratios have proven to be a valuable discriminant in support of monitoring a Comprehensive Nuclear Test Ban Treaty. Regional S and P phases attenuate at different rates and the attenuation can vary geographically. Therefore, calibration is needed to apply the S/P discriminant in new regions. Calibration includes application of frequency-dependent source and distance corrections for regional Pn, Pg, Sn, and Lg phases.¶Jenkins et al. (1998) developed Pn, Pg, Sn, and Lg amplitude models for nine geographic regions and two global composite models, stable and tectonic. They determined frequency-dependent source and attenuation corrections from a large data set obtained from the Prototype International Data Center (PIDC). We use their corrections to evaluate calibrated S/P discriminants.¶Our discrimination data set includes >1000 amplitude ratios from earthquakes, industrial explosions, chemical explosions, and nuclear explosions from Lop Nor, India and Pakistan. We find that the calibrated S/P ratio is largest for earthquakes and smallest for the nuclear explosions, as expected. However, the discriminant is not universally valid. In particular, the S/P ratio for the Pakistan nuclear explosion fell within the normal range for the earthquakes. This event was recorded by only a few stations at far-regional distances and appears to have an anomalously high Sn amplitude. The industrial explosions overlap with the earthquake population, however the buried chemical explosions generally register lower S/P ratio than earthquakes.  相似文献   

5.
P-wave travel-time residuals at the Warramunga Seismic Array (WRA) in the Northern Territory, Australia, have been studied from 49 earthquakes with epicenters south of 19°S in the Fiji-Tonga region. Focal depths are between 42 and 679 km as determined from pP-P. Using the Jeffreys-Bullen and the Herrin travel-time tables the epicentral parameters have been redetermined by considering only “normal” seismic stations in the location procedure. These are those stations where P-wave travel times are probably not affected by lateral heterogeneities caused by the lithosphere descending beneath the Tonga trench. Epicenters of deep earthquakes below 300 km have been relocated by using stations at Δ > 25° only. Epicenters from shallower-depth earthquakes have been recalculated without using stations between 35 < Δ < 75° epicentral distance. In both cases focal depths were determined from pP-P times. The resulting pattern of P-residuals at WRA does not show any significant change with depth below 350 km. The residuals become more negative for shallower earthquakes above about 250 km. P-waves to WRA are advanced by approximately 2 s compared with those from deep earthquakes. The results do not essentially differ for the two different travel-time tables used. The observations can be interpreted by P-wave velocities that are higher in the sinking slab down to 350–400 km by 5±2% than in both the Jeffreys-Bullen and Herrin models. Without considering possible elevations of phase boundaries this estimate yields a temperature contrast of 1000±450°C between slab and normal mantle material in this depth range.  相似文献   

6.
Improving Regional Seismic Event Location in China   总被引:1,自引:0,他引:1  
—?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.  相似文献   

7.
A method for analysing travel times measured at a large array or a network of seismographs from many earthquakes within a specific region has been developed. Approximate relative station corrections are calculated from the residuals on a least-squares line or least-squares quadratic form fitted through the times for each earthquake, and may be improved by iteration after a preliminary travel-time curve has been derived. Accurate relative baseline corrections for each earthquake are also calculated iteratively, and an optimum slowness-distance curve is determined from the combined corrected travel times from all earthquakes using a trade-off procedure. Calculations using synthetic travel-time data suggest that abrupt changes in slowness of ~ 0.4 s deg?1 due to the presence of triplications are generally resolvable, provided that the effects of lateral variations are small, even with random epicentre mislocations in the range ± 0.5°. Slowness measurements at a network of temporary stations deployed across Australia do not show any discontinuities in slowness greater than 0.2 s deg?1 in the distance range 45–54°. Similar measurements at the Warramunga array from the same source regions, however, suggest the presence of complexity in the slowness curve at distances close to 50°. Relative arrival times at the temporary network generally have standard deviations less than 0.25 s, thus suggesting that details of structure finer than those derived from conventional travel-time studies can be resolved.  相似文献   

8.
—?Identification of seismic events is a major scientific issue in the framework of verification of the Comprehensive Nuclear-Test-Ban Treaty (CTBT). Of special interest in this context is the identification of the numerous low-yield mining or blasting events, especially those occurring in the same area as earthquakes, such as the Vogtland area in the border region of Germany and the Czech Republic. Seismic events in this area were investigated by WÜSTER (1993, 1995), who achieved complete discrimination using measures of spectral decay and spectral variance at the GERESS array and a quadratic discrimination function.¶A subset of these events, for which ground-truth information is available, has been analyzed in this study using multivariate statistical analysis. Various parameters based on measurements from seismic waveforms of the broadband stations of the German Regional Seismic Network (GRSN) and short-period elements of the GERESS array are tested for statistical significance in a linear regression analysis, in particular spectral amplitude ratios for the L g phase and P g /S g amplitude ratios. The subset includes a total of 35 explosions and 24 earthquakes. The results of our study argue that identification based on spectral L g and high-frequency P g /S g ratios is promising. However, discrimination success is strongly varying from station to station; thus, weighting according to station success rates could improve the overall identification capability.  相似文献   

9.
We use the dense Israel Seismic Network (ISN) to discriminate between low magnitude earthquakes and explosions in the Middle East region. This issue is important for CTBT monitoring, especially when considering small nuclear tests which may be conducted under evasive conditions. We explore the performance of efficient discriminants based on spectral features of seismograms using waveforms of 50 earthquakes and 114 quarry and underwater blasts with magnitudes 1.0–2.8, recorded by ISN short-period stations at distances up to 200 km. The single-station spectral ratio of the low and high-frequency seismic energy shows an overlap between explosions and earthquakes. After averaging over a subnet of stations, the resolving power is enhanced and the two classes of events are separated. Different frequency bands were tested; the (1–3 Hz)/(6–8 Hz) ratio provided the best discriminant performance. We also estimated normalized r.m.s. spectral amplitudes in several sequential equal frequency windows within the 1–12 Hz band and applied multiparametric automatic classification procedures (Linear Discrimination Function and Artificial Neural Network) to the amplitudes averaged over a subnetwork. A leave-one-out test showed a low rate of error for the multiparametric procedures. An innovative multi-station discriminant is proposed, based on spectral modulation associated with ripple-firing in quarry blasts and with the bubbling effect in underwater explosions. It utilizes a distinct azimuth-invariant coherency of spectral shapes for different stations in the frequency range (1–12 Hz). The coherency is measured by semblance statistics commonly used in seismic prospecting for phase correlation in the time domain. After modification, the statistics applied to the network spectra provided event separation. A new feature of all the above mentioned procedures is that they are based on smoothed (0.5 Hz window), instrument-corrected FFT spectra of the whole signal; they are robust to the accuracy of onset time estimation and, thus well suited to automatic event identification.  相似文献   

10.
The Method of Dynamic Calibration (MDC) of stations of the International Monitoring System (IMS) was developed for calibrating regions where no underground nuclear explosions were carried out, with the purpose of providing conditions for implementation of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) in nontrivial cases. Initially, the MDC had been presented in [Kedrov, 2001; Kedrov et al., 2001; Kedrov O.K. and Kedrov E.O., 2003] and then considered in detail in [Kedrov et al., 2008]. The core of MDC relates to adapting diagnostic parameters for the identification of underground nuclear explosions (UNE) and earthquakes elaborated for the region of Eurasia, taken as a basic region (BR), for other researched regions that differ from BR in the character of the attenuation of seismic waves. The unique characteristic of this method lies in the fact that calibration of diagnostic parameters with the help of attenuation coefficients b Δ at varied source-station traces is implemented using only natural seismicity data within the limits of an explored region and does not require special underground chemical explosions. The MDC algorithm is implemented in the research program ”Kalibr”, which was tested by using the experimental data from Eurasia region. It is shown in this work that MDC can be used for calibration of regions where a very low level of natural seismicity is observed. According to the results of the calibration of diagnostic parameters at IMS stations in several regions of North America, Africa, and Asia, the approximate classification of propagation conditions for seismic signals at source-station traces in platform and tectonically active regions is made. The results for the development of two research programs, “Spektr” and “Signal”, are presented; this software is intended for automation of calculation procedures for spectral diagnostic parameters of UNEs’ and earthquakes’ identification by amplitude spectra of P waves and by the maximal amplitudes of P, S, and LR signals. The application of these programs allowed us to accelerate the whole calibration procedure for a particular source-station trace using the ”Kalibr” program.  相似文献   

11.
v—vAt the German NDC initial work on seismic event identification has focused on the application of spectral amplitude ratios for Lg in order to discriminate naturally occurring seismic events from other events associated with mining and quarry activities. Only about 10% of all seismic events occurring in Germany and adjacent areas are due to natural seismicity and are mostly constrained to the Alpine regions and areas along the Rhinegraben, Rhenish massif, Swabian Jura, and the Bohemian massif (Vogtland region). Using data from the broadband GRSN network, spectral amplitude ratios are calculated from maximum trace amplitudes in the 1–2uHz and 6–8uHz frequency bands, which are within the passbands of the deployed STS-2 instruments and the recorded 20uHz data streams. These amplitude ratios then must be corrected with an appropriate attenuation model in order to remove propagation paths effects. For event identification, a scoring scheme is applied across the GRSN network, based on station-dependent scoring thresholds. In a case study aimed at testing the identification scheme, events are investigated from a quarry in southern Germany that provided ground-truth information for six events in 1997 to demonstrate the suitability of this identification approach. Except for one event with a rather strong earthquake signature, i.e., a low spectral Lg ratio, these events could be screened out from the earthquake population by their large Lg ratios. In a second step, aimed at applying the identification scheme, all events in Germany and neighboring areas that occurred in 1995 were processed, with approximately 800 out of more than 1200 events showing explosion-type Lg ratios, while only 10% remain in the earthquake population. However, specific mining areas appear to consistently produce earthquake-type spectral ratios indicative of particular blasting practices.  相似文献   

12.
—?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.  相似文献   

13.
Unified Local Magnitude Scale for Earthquakes of South Balkan Area   总被引:1,自引:0,他引:1  
—A homogeneous earthquake catalog spanning 1964–1995 for the southern Balkan area is presented, by expressing the size of the earthquakes in a unified local magnitude scale. The strategy followed to produce this catalog is also presented. Local magnitudes calculated by six Balkan seismological centers (ATH, THE, ISK, TIR, TTG, SKO) have been used in order to obtain relations between the estimated ML values from the Greek seismological centers (ATH and THE) and the remaining four Balkan stations. Since it was found that local magnitudes estimated by ATH and THE are almost identical, they have been used as one data set so as to correlate with the data of each one of the remaining four seismological stations. Based on the proposed relations, a unified local magnitude, MLGR, is given for each earthquake of the regional catalog. A published scaling relation between the ML values from ATH and THE networks and the corresponding seismic moment magnitude, Mw, was used in combination with the above relations, in order to enable the conversion of any ML value from any station into Mw. The catalog completeness has been checked and the b-value has been calculated for the complete data sample.  相似文献   

14.
This study analyzes and compares the P- and S-wave displacement spectra from local earthquakes and explosions of similar magnitudes. We propose a new approach to discrimination between low-magnitude shallow earthquakes and explosions by using ratios of P- to S-wave corner frequencies as a criterion. We have explored 2430 digital records of the Israeli Seismic Network (ISN) from 456 local events (226 earthquakes, 230 quarry blasts, and a few underwater explosions) of magnitudes Md?=?1.4–3.4, which occurred at distances up to 250 km during 2001–2013 years. P-wave and S-wave displacement spectra were computed for all events following Brune’s source model of earthquakes (1970, 1971) and applying the distance correction coefficients (Shapira and Hofstetter, Teconophysics 217:217–226, 1993; Ataeva G, Shapira A, Hofstetter A, J Seismol 19:389-401, 2015), The corner frequencies and moment magnitudes were determined using multiple stations for each event, and then the comparative analysis was performed.The analysis showed that both P-wave and especially S-wave displacement spectra of quarry blasts demonstrate the corner frequencies lower than those obtained from earthquakes of similar magnitudes. A clear separation between earthquake and explosion populations was obtained for ratios of P- to S-wave corner frequency f 0(P)/f 0(S). The ratios were computed for each event with corner frequencies f 0 of P- and S-wave, which were obtained from the measured f 0 I at individual stations, then corrected for distance and finally averaged. We obtained empirically the average estimation of f 0(P)/f 0(S)?=?1.23 for all used earthquakes, and 1.86 for all explosions. We found that the difference in the ratios can be an effective discrimination parameter which does not depend on estimated moment magnitude M w .The new multi-station Corner Frequency Discriminant (CFD) for earthquakes and explosions in Israel was developed based on ratios P- to S-wave corner frequencies f 0(P)/f 0(S), with the empirical threshold value of the ratio for Israel as 1.48.  相似文献   

15.
A 12-year period experimental data, from 1 January 1995, to 20 August 2007, have been examined for possible correlations between Seismic Electric Signals (SES) of the VAN method and source parameters of the corresponding earthquakes in Western Greece. During that period 13 earthquakes of magnitude M W ≥ 5 with epicenters in the area 19–24°E, 36–41°N and available the CMT solutions (Centroid Moment Tensor focal mechanism solutions) have been found to be preceded by a SES each time at one of the three VAN stations of IOA, PIR or PAT. The results of IOA and PIR stations are compared to those reported by Uyeda et al. 1999, for the previous period 1983–1994. The IOA station seems to be sensitive to earthquakes with thrust type mechanisms being mainly consistent to its past behavior. The PIR station, which is moved from its previous position by few kilometers, exhibits some changes. It detects now both strike slip and thrust type earthquakes and becomes sensitive to new areas indicating probably the strong dependence of the SES station detection ability upon its location.  相似文献   

16.
利用重复地震研究山东乳山地区地壳介质波速变化   总被引:3,自引:0,他引:3       下载免费PDF全文
王鹏  郑建常  谭毅培 《地震学报》2016,38(5):728-738
利用2013年10月—2015年6月山东乳山震群的地震波形资料,基于距离乳山震群最近的3个台站所记录的地震波形互相关系数c≥0.9的地震即为重复地震的原则,识别了乳山震群的重复地震.利用波形互相关时延法,计算分析了文登台和招远台这两个固定台站的走时差随时间的变化特征.结果表明:乳山震群的地震相似性很好,且震中位置较为集中,重复地震的时间跨度较大,有利于提高观测结果的时间分辨率;乳山震群中3次MS>4.0地震前均出现了短期的走时差低值异常现象,反映了乳山震群震源区在中强震前会出现短期地壳介质速度明显升高的过程,可以为该区域的中强震预测提供一定依据.   相似文献   

17.
Forensic seismology revisited   总被引:1,自引:0,他引:1  
The first technical discussions, held in 1958, on methods of verifying compliance with a treaty banning nuclear explosions, concluded that a monitoring system could be set up to detect and identify such explosions anywhere except underground: the difficulty with underground explosions was that there would be some earthquakes that could not be distinguished from an explosion. The development of adequate ways of discriminating between earthquakes and underground explosions proved to be difficult so that only in 1996 was a Comprehensive Nuclear Test Ban Treaty (CTBT) finally negotiated. Some of the important improvements in the detection and identification of underground tests—that is in forensic seismology—have been made by the UK through a research group at the Atomic Weapons Establishment (AWE). The paper describes some of the advances made in identification since 1958, particularly by the AWE Group, and the main features of the International Monitoring System (IMS), being set up to verify the Test Ban. Once the Treaty enters into force, then should a suspicious disturbance be detected the State under suspicion of testing will have to demonstrate that the disturbance was not a test. If this cannot be done satisfactorily the Treaty has provisions for on-site inspections (OSIs): for a suspicious seismic disturbance for example, an international team of inspectors will search the area around the estimated epicentre of the disturbance for evidence that a nuclear test really took place. Early observations made at epicentral distances out to 2,000 km from the Nevada Test Site showed that there is little to distinguish explosion seismograms from those of nearby earthquakes: for both source types the short-period (SP: ∼1 Hz) seismograms are complex showing multiple arrivals. At long range, say 3,000–10,000 km, loosely called teleseismic distances, the AWE Group noted that SP P waves—the most widely and well-recorded waves from underground explosions—were in contrast simple, comprising one or two cycles of large amplitude followed by a low-amplitude coda. Earthquake signals on the other hand were often complex with numerous arrivals of similar amplitude spread over 35 s or more. It therefore appeared that earthquakes could be recognised on complexity. Later however, complex explosion signals were observed which reduced the apparent effectiveness of complexity as a criterion for identifying earthquakes. Nevertheless, the AWE Group concluded that for many paths to teleseismic distances, Earth is transparent for P signals and this provides a window through which source differences will be most clearly seen. Much of the research by the Group has focused on understanding the influence of source type on P seismograms recorded at teleseismic distances. Consequently the paper concentrates on teleseismic methods of distinguishing between explosions and earthquakes. One of the most robust criteria for discriminating between earthquakes and explosions is the m b : M s criterion which compares the amplitudes of the SP P waves as measured by the body-wave magnitude m b, and the long-period (LP: ∼0.05 Hz) Rayleigh-wave amplitude as measured by the surface-wave magnitude M s; the P and Rayleigh waves being the main wave types used in forensic seismology. For a given M s, the m b for explosions is larger than for most earthquakes. The criterion is difficult to apply however, at low magnitude (say m b < 4.5) and there are exceptions—earthquakes that look like explosions. A difficulty with identification criteria developed in the early days of forensic seismology was that they were in the main empirical—it was not known why they appeared to work and if there were test sites or earthquakes where they would fail. Consequently the AWE Group in cooperation with the University of Cambridge used seismogram modelling to try and understand what controls complexity of SP P seismograms, and to put the m b : M s criterion on a theoretical basis. The results of this work show that the m b : M s criterion is robust because several factors contribute to the separation of earthquakes and explosions. The principal reason for the separation however, is that for many orientations of the earthquake source there is at least one P nodal plane in the teleseismic window and this biases m b low. Only for earthquakes with near 45° dip-slip mechanisms where the antinode of P is in the source window is the m b:M s criterion predicted to fail. The results from modelling are consistent with observation—in particular there are earthquakes, “anomalous events”, which look explosion-like on the m b:M s criterion, that turn out to have mechanisms close to 45° dip-slip. Fortunately the P seismograms from such earthquakes usually show pP and sP, the reflections from the free surface of P and S waves radiated upwards. From the pP–P and sP–P times the focal depth can be estimated. So far the estimated depth of the anomalous events have turned out to be ∼20 km, too deep to be explosions. Studies show that the observation that P seismograms are more complex than predicted by simple models can be explained on the weak-signal hypothesis: the standard phases, direct P and the surface reflections, are weak because of amongst other things, the effects of the radiation pattern or obstacles on the source-to-receiver path; other non-standard arrivals then appear relatively large on the seismograms. What has come out of the modelling of P seismograms is a criterion for recognising suspicious disturbances based on simplicity rather than complexity. Simple P seismograms for earthquakes at depths of more than a few kilometres are likely to be radiated only to stations that lie in a confined range of azimuths and distances. If then, simple seismograms are recorded over a wide range of distances and particularly azimuths, it is unlikely the source is an earthquake at depth. It is possible to test this using the relative amplitudes of direct P and later arrivals that might be surface reflections. The procedure is to use only the simple P seismograms on the assumption that whereas the propagation through Earth may make a signal more complex it is unlikely to make it simpler. From the amplitude of the coda of these seismograms, bounds can be placed on the size of possible pP and sP. The relative-amplitude method is then used to search for orientations of the earthquake source that are compatible with the observations. If no such orientations are found the source must be shallow so that any surface reflections merge with direct P, and hence could be an explosion. The IMS when completed will be a global network of 321 monitoring stations, including 170 seismological stations principally to detect the seismic waves from earthquakes and underground explosions. The IMS will also have stations with hydrophones, microbarographs and radionuclide detectors to detect explosions in the oceans and the atmosphere and any isotopes in the air characteristic of a nuclear test. The Global Communications Infrastructure provides communications between the IMS stations and the International Data Centre (IDC), Vienna, where the recordings from the monitoring stations is collected, collated, and analysed. The IDC issues bulletins listing geophysical disturbances, to States Signatories to the CTBT. The assessment of the disturbances to decide whether any are possible explosions, is a task for State Signatories. For each Signatory to do a detailed analysis of all disturbances would be expensive and time consuming. Fortunately many disturbances can be readily identified as earthquakes and removed from consideration—a process referred to as “event screening”. For example, many earthquakes with epicentres over the oceans can be distinguished from underwater explosions, because an explosion signal is of much higher frequency than that of earthquakes that occur below the ocean bed. Further, many earthquakes could clearly be identified at the IDC on the m b : M s criterion, but there is a difficulty—how to set the decision line. The possibility has to be very small that an explosion will be classed by mistake, as an earthquake. The decision line has therefore to be set conservatively, consequently with routine application of current screening criteria, only about 50% of earthquakes can be positively identified as such. Various methods have been proposed whereby a “determined violator” could avoid the provisions of a CTBT and carry out a test that would be either undetected or detected but not identified as an explosion. The increase in complexity and cost of such a test should discourage any State from attempting it. In addition, there is always the possibility of some stations detecting the test, the test being identified as suspicious, and so subject to an OSI. With time as the IMS becomes more efficient and effective it will act increasingly to deter anyone contemplating a clandestine test, from going ahead. What has emerged is several robust criteria. The criteria include: location, which when combined with hydro-acoustic data can identify earthquakes under the sea; m b : M s; and depth of focus. More detailed study is required of any remaining seismic disturbance that is regarded as suspicious: for example, is close to a site where nuclear tests have been carried out in the past. Any disturbance that is shown to be explosion-like, may be the subject of an OSI. One surprise is how little plate tectonics has contributed to resolving problems in forensic seismology. Much of the evidence for plate tectonics comes from seismological studies so it would be expected that the implications for Earth structure arising from forensic seismology would be consistent with plate-tectonic models. So far the AWE Group have found little synergy between plate tectonics and forensic seismology. It is to be hoped that the large volume of seismological data of high quality now being collected by the IMS and the increasing number of digital stations, will result in a revised Earth model that is consistent with the findings of forensic seismology, so that a future review of progress will show that the forensic seismologist can draw on this model in attempting to interpret apparently anomalous seismograms.
A. DouglasEmail:
  相似文献   

18.
选取重庆数字地震台网记录的武隆地区107次ML≥1.5地震波形和观测报告资料,采用波形互相关技术进行相关计算,识别出12组同时被3个地震台站记录且各台波形互相关系数不小于0.8的重复地震。利用射线追踪方法,对武隆地区重复地震位置做归一化处理,得到所选台站记录地震P波走时差变化,结果显示:在2017年武隆MS5.0地震前4年,武隆、仙女山、涪陵地震台记录到P波走时差均存在长期的负异常变化,持续时间长达339天,共出现18次负值,表明该时段地壳介质速度有明显升高现象。   相似文献   

19.
It is a common opinion that only crustal earthquakes can occur in the Crimea–Black Sea region. Since the existence of deep earthquakes in the Crimea–Black Sea region is extremely important for the construction of a geodynamic model for this region, an attempt is made to verify the validity of this widespread view. To do this, the coordinates of all earthquakes recorded by the stations of the Crimean seismological network are reinterpreted with an algorithm developed by one of the authors. The data published in the seismological catalogs and bulletins of the Crimea–Black Sea region for 1970–2012 are used for the analysis. To refine the coordinates of hypocenters of earthquakes in the Crimea–Black Sea region, in addition to the data from stations of the Crimean seismological network, information from seismic stations located around the Black Sea coast are used. In total, the data from 61 seismic stations were used to determine the hypocenter coordinates. The used earthquake catalogs for 1970–2012 contain information on ~2140 events with magnitudes from–1.5 to 5.5. The bulletins provide information on the arrival times of P- and S-waves at seismic stations for 1988 events recorded by three or more stations. The principal innovation of this study is the use of the original author’s hypocenter determination algorithm, which minimizes the functional of distances between the points (X, Y, H) and (x, y, h) corresponding to the theoretical and observed seismic wave travel times from the earthquake source to the recording stations. The determination of the coordinates of earthquake hypocenters is much more stable in this case than the usual minimization of the residual functional for the arrival time of an earthquake wave at a station (the difference between the theoretical and observed values). Since determination of the hypocenter coordinates can be influenced by the chosen velocity column beneath each station, special attention is focused on collecting information on velocity profiles. To evaluate the influence of the upper mantle on the results of calculating the velocity model, two different low-velocity and high-velocity models are used; the results are compared with each other. Both velocity models are set to a depth of 640 km, which is fundamentally important in determining hypocenters for deep earthquakes. Studies of the Crimea–Black Sea region have revealed more than 70 earthquakes with a source depth of more than 60 km. The adequacy of the obtained depth values is confirmed by the results of comparing the initial experimental data from the bulletins with the theoretical travel-time curves for earthquake sources with depths of 50 and 200 km. The sources of deep earthquakes found in the Crimea–Black Sea region significantly change our understanding of the structure and geotectonics of this region.  相似文献   

20.
Automated Detection, Extraction, and Measurement of Regional Surface Waves   总被引:3,自引:0,他引:3  
v--vOur goal is to develop and test an effective method to detect, identify, extract, and quantify surface wave signals for weak events observed at regional stations. We describe an automated surface wave detector and extractor designed to work on weak surface wave signals across Eurasia at intermediate periods (8 s-40 s). The method is based on phase-matched filters defined by the Rayleigh wave group travel-time predictions from the broadband group velocity maps presented by Ritzwoller and Levshin (1998) and Ritzwoller et al. (1998) and proceeds in three steps: Signal compression, signal extraction or cleaning, and measurement. First, the dispersed surface wave signals are compressed in time by applying an anti-dispersion or phase-matched filter defined from the group velocity maps. We refer to this as the `compressed signal.' Second, the surface wave is then extracted by filtering `noise' temporally isolated from the time-compressed signal. This filtered signal is then redispersed by applying the inverse of the phase-matched filter. Finally, we adaptively estimate spectral amplitude as well as group and phase velocity on the filtered signal. The method is naturally used as a detector by allowing origin time to slide along the time axis. We describe preliminary results of the application of this method to a set of nuclear explosions and earthquakes that occurred on or near the Chinese Lop Nor test site from 1992 through 1996 and one explosion on the Indian Rajasthan test site that occurred in May of 1998.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号