Seismic Site Classification and Correlation between Standard Penetration Test N Value and Shear Wave Velocity for Lucknow City in Indo-Gangetic Basin

Subsurface lithology and seismic site classification of Lucknow urban center located in the central part of the Indo-Gangetic Basin (IGB) are presented based on detailed shallow subsurface investigations and borehole analysis. These are done by carrying out 47 seismic surface wave tests using multichannel analysis of surface waves (MASW) and 23 boreholes drilled up to 30 m with standard penetration test (SPT) N values. Subsurface lithology profiles drawn from the drilled boreholes show low- to medium-compressibility clay and silty to poorly graded sand available till depth of 30 m. In addition, deeper boreholes (depth >150 m) were collected from the Lucknow Jal Nigam (Water Corporation), Government of Uttar Pradesh to understand deeper subsoil stratification. Deeper boreholes in this paper refer to those with depth over 150 m. These reports show the presence of clay mix with sand and Kankar at some locations till a depth of 150 m, followed by layers of sand, clay, and Kankar up to 400 m. Based on the available details, shallow and deeper cross-sections through Lucknow are presented. Shear wave velocity (SWV) and N-SPT values were measured for the study area using MASW and SPT testing. Measured SWV and N-SPT values for the same locations were found to be comparable. These values were used to estimate 30 m average values of N-SPT (N ₃₀) and SWV (V _s ³⁰ ) for seismic site classification of the study area as per the National Earthquake Hazards Reduction Program (NEHRP) soil classification system. Based on the NEHRP classification, the entire study area is classified into site class C and D based on V _s ³⁰ and site class D and E based on N ₃₀. The issue of larger amplification during future seismic events is highlighted for a major part of the study area which comes under site class D and E. Also, the mismatch of site classes based on N ₃₀ and V _s ³⁰ raises the question of the suitability of the NEHRP classification system for the study region. Further, 17 sets of SPT and SWV data are used to develop a correlation between N-SPT and SWV. This represents a first attempt of seismic site classification and correlation between N-SPT and SWV in the Indo-Gangetic Basin.     

Investigation of mining-induced earthquakes in Iran within a time window of 2006–2013

Iran is one of the most seismically active areas of the world and frequently suffers from destructive earthquakes. Rare studies on anthropogenic-induced seismicity in Iran may be related to less attention to triggered events and more concern to natural origin of earthquakes. Hence, the present study as a frontier research aims to investigate the mining-induced earthquakes in Iran. For this purpose, distribution of ~?76,000 seismic events was investigated between the years 2006 and 2013. This study considered a correlation test to investigate the possible mining triggering of the seismic events based on a network of 194 geographical pixels (1°?×?1°) in ArcGIS. Results conveniently confirmed a positive meaningful relation between all earthquake events with magnitudes M?>?0.5 and mining activities in Iran (R?=?0.42). Detailed results confirmed that the most of earthquake swarms (at least ~?60%) had mining-induced origin, which were spatially located in same pixels of metallic mineral mining sites. The correlation test between earthquake swarms and mining activities indicated positive and meaningful relationships in four regions of Alborz, Kopet Dag, Kerman, and Zagros, respectively (R?=?0.61, 0.54, 0.51, and 0.50). Hence, aforementioned seismic regions exposed sensitive seismic responses toward mining triggering effects in Iran.     

Detection and Analysis of Near-Surface Explosions on the Kola Peninsula

Seismic and infrasonic observations of signals from a sequence of near-surface explosions at a site on the Kola Peninsula have been analyzed. NORSAR’s automatic network processing of these events shows a significant scatter in the location estimates and, to improve the automatic classification of the events, we have performed full waveform cross-correlation on the data set. Although the signals from the different events share many characteristics, the waveforms do not exhibit a ripple-for-ripple correspondence and cross-correlation does not result in the classic delta-function indicative of repeating signals. Using recordings from the ARCES seismic array (250 km W of the events), we find that a correlation detector on a single channel or three-component station would not be able to detect subsequent events from this source without an unacceptable false alarm rate. However, performing the correlation on each channel of the full ARCES array, and stacking the resulting traces, generates a correlation detection statistic with a suppressed background level which is exceeded by many times its standard deviation on only very few occasions. Performing f-k analysis on the individual correlation coefficient traces, and rejecting detections indicating a non-zero slowness vector, results in a detection list with essentially no false alarms. Applying the algorithm to 8 years of continuous ARCES data identified over 350 events which we confidently assign to this sequence. The large event population provides additional confidence in relative travel-time estimates and this, together with the occurrence of many events between 2002 and 2004 when a temporary network was deployed in the region, reduces the variability in location estimates. The best seismic location estimate, incorporating phase information for many hundreds of events, is consistent with backazimuth measurements for infrasound arrivals at several stations at regional distances. At Lycksele, 800 km SW of the events, as well as at ARCES, infrasound is detected for most of the events in the summer and for few in the winter. At Apatity, some 230 km S of the estimated source location, infrasound is detected for most events. As a first step to providing a Ground Truth database for this useful source of infrasound, we provide the times of explosions for over 50 events spanning 1 year.     

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.     

四川、甘肃地区VS30经验估计研究

目前我国建筑工程抗震设计规范中对于工程场地条件的判断依据主要是地表以下20m深度范围内土层的等效剪切波速,简称VS20。相比之下,国外应用较广的是地表以下30m深度范围内的等效剪切波速,简称VS30。这种差别导致国内科研工作者在应用国外的地震工程、工程抗震模型时经常遇到对场地条件描述不准确的困难。为了解决这个问题,本文根据147个四川、甘肃地区国家强震动台站20m左右深度的钻孔剪切波速数据,利用延拓方法、场地分类统计方法以及基于地形特征的VS30估计方法研究各台站VS30与VS20的经验关系,对比发现基于速度梯度延拓的结果最为可取。参考国际上通用的Geomatrix Classification场地分类标准,最终得到四川、甘肃地区各类场地的平均VS30,此结果可以为缺乏钻孔数据的工程场地的VS30估计提供参考。     

“Great East Siberian” earthquake of February 1, 1725 (M = 8.2): Revision of the main parameters in the light of additional data and modern requirements

The earthquake of February 1 (January 21 in the Julian Calendar), 1725 was for a long time considered not only the earliest precisely dated historical seismic event in eastern Siberia, but also the most powerful earthquake for the entire period of recorded seismic events in the region: M = 8.2 (Novyi katalog…, 1977). The epicenter location (Stanovoy Upland) and the magnitude of the event were assessed on the basis of very scanty historical data, as well as using paleoseismogeological information. The February 1, 1725 event received the name “The Great East Siberian” earthquake and served for decades as decisive evidence for the assessment of the seismic hazard and seismic zoning of the northeastern flank of the Baikal rift zone. However, the solution of the focal parameters in the (Novyi katalog…, 1977) has caused serious doubts. In this paper a newly elaborated version is proposed that is based on a detailed reevaluation of the initially known macroseismic information, as well as additional historical data that previously had not come to the attention of seismologists. As the result, a different solution of the focal parameters (51.8° N; 113.0° E, eastern Transbaikalia) and a significantly lower magnitude (M = 6.0) compared with the parameters given in (Novyi katalog…, 1977) were obtained. The presented solution makes us more attentive to the estimates of seismic hazard in east Transbaikalia based on historical data, as well as to the historical data themselves.     

Seismotectonic study of Kishtwar region of Jammu Province using local broadband seismic data

Two moderate earthquakes of Mw 5.7 on the first of May and Mw 5.2 on the second of August occurred in the Kishtwar region in the year 2013. Our broadband seismic observatories located in the region recorded these events and the aftershocks. We analyzed these data to understand the seismotectonics of this region. Most of the events were located between 33.03° to 33.29° N latitude and 75.40° to 76.07° E longitude. Focal depths of these shallow earthquakes range from 7 to 12 km and are confined between Panjal Thrust (PT) and Kishtwar Window (KW). Spectral analysis of these events reveals that stress drop, source radius, corner frequency, and moment magnitude varied between 3.3 and 70.1 bars, 0.121 and 3.55 km, 0.397 and 6.06 Hz, and Mw 2.2 and Mw 5.7, respectively. The low stress drop of small-magnitude earthquakes reveals the brittle nature of the upper crust which is coincident with the field observations. The variation of stress drop with magnitude shows positive correlation whereas no such relation was observed between stress drop and depth of focus. The b value calculated (0.83) for the area reveals high stress accumulation within the incompetent rock zones in the area.     

Correlations between electrical and elastic properties of solid–liquid composites with interfacial energy-controlled equilibrium microstructures

Electrical conductivity and seismic velocity are studied for plausible pore geometries in the Earth's interior for reliable quantitative analysis of experimental data such as seismic tomography and magnetotelluric explorations. Electrical conductivity of a two-phase system with equilibrium, interfacial energy-controlled phase geometry is calculated for the dihedral angles θ = 40°–100° that are typical for rock–aqueous fluid and θ = 20°–60° for rock–melt systems of lower crust and upper mantle for the case of tetrakaidecahedral grains. Electrical conductivity vs. seismic velocity correlations are acquired by combining of the simulated electrical conductivities with the seismic velocity calculated with the help of equilibrium geometry model Takei [Takei, Y., Effect of pore geometry on V_P/V_S: From equilibrium geometry to crack. J. Geophys. Res. 107 (2002): 10.1029/2001JB000522.] for the same pore geometries. The results show that electrical conductivity gradually decreases reaching zero when seismic velocities reach seismic velocities of intact rock for rock–melt systems, while for rock–aqueous fluid systems with θ ≥ 60° conductivity drops to zero at velocities up to 10% smaller. This can explain the seeming discrepancy of the low seismic velocity region, attributed to the high fluid fraction, and the low electrical conductivity of the same region, which is sometimes faced at collocated electromagnetic and seismic experiments.     

Site response and source spectra of S waves in the Zagros region, Iran

S wave amplitude spectra from shallow earthquakes with magnitudes ranging between 4.2 and 6.2 in the Zagros region of Iran that occurred between 1998 and 2008 are used to examine source parameters and site response of S waves. A generalized inversion scheme has been used to separate the source, propagation path, and local site effects from S wave spectra. For removing the trade-off between source and site terms and propagation effects (including geometric and anelastic attenuation), the spectral amplitudes of the records used were corrected for attenuation and geometrical spreading function using a path model proposed by Zafarani and Soghrat (Bull Seism Soc Am 102:2031–2045, 2012) for the region. We assume a Brune’s point source model to retrieve source parameters like corner frequency, moment magnitude, and high-frequency fall off coefficient, for each event. When the source spectra are interpreted in terms of Brune’s model, the average stress drops obtained are about 7.1 and 5.9 MPa (71 and 59 bars), respectively for the eastern and western Zagros regions. Stress drops range from 1.4 to 35.0 MPa (14 to 350 bars), with no clear dependence on magnitude. The results in terms of stress drop and S wave seismic energy indicate that the Zagros events are more similar to interplate earthquakes of western North America than to intraplate events of eastern North America. The method also provides us with site responses for all 40 stations individually and is an interesting alternative to other methods, such as the H/V method. A new empirical relationship between body-wave magnitudes and moment magnitude has been proposed for the Iranian plateau using derived seismic moment from the inversion.     

Earthquake relocation and 3-dimensional crustal structure of P-wave velocity in cen-tral-western China

IntroductionUnderstandingthemechanismofcontinentalearthquakesisveryimportantforseismichaz-ardpreventionandearthquakeprediction.Themodernseismotectonictheoryandtheideaofearthquakepredictionaredevelopedmainlyfromthestudiesoninterplateearthquakes,whicharedifficulttoexplainthephenomenaofintraplateearthquakes,suchasthecontinentalearthquakesoccurredinChinesemainland.Whiletheinterplateearthquakesoccurredalongtheplatebounda-ries,theintraplateearthquakesdistributediffuselyintheinterioroftheplates.Thus…     

Detection capability of seismic network based on noise analysis and magnitude of completeness

Assessing the detection threshold of seismic networks becomes of increased importance namely in the context of monitoring induced seismicity due to underground operations. Achieving the maximum possible sensitivity of industrial seismic monitoring is a precondition for successful control of technological procedures. Similarly, the lowest detection threshold is desirable when monitoring the natural seismic activity aimed to imaging the fault structures in 3D and to understanding the ongoing processes in the crust. We compare the application of two different methods to the data of the seismic network WEBNET that monitors the earthquake swarm activity of the West-Bohemia/Vogtland region. First, we evaluate the absolute noise level and its possible non-stationary character that results in hampering the detectability of the seismic network by producing false alarms. This is realized by the statistical analysis of the noise amplitudes using the ratio of 99 and 95 percentiles. Second, the magnitude of completeness is determined for each of the nine stations by analysing the automatic detections of an intensive swarm period from August 2011. The magnitude–frequency distributions of all detected events and events detected at individual stations are compared to determine the magnitude of completeness at a selected completeness level. The resulting magnitude of completeness M _c of most of the stations varies between ?0.9 and ?0.5; an anomalous high M _c of 0.0 is found at the most distant station, which is probably due to inadequate correction for attenuation. We find that while the absolute noise level has no significant influence to the station sensitivity, the noise stationarity correlates with station sensitivity expressed in low magnitude of completeness and vice versa. This qualifies the method of analysing the stationary character of seismic noise as an effective tool for site surveying during the seismic station deployment.     

A surface-wave investigation of the rupture mechanism of the Gobi-Altai (December 4, 1957) earthquake

Long-period records of multiple Love waves from the 1957 earthquake in Mongolia (M_S = 8.0) at Pasadena are analysed and compared to synthetic seismograms, generated by the method of Kanamori. A fit in the time domain shows that the records are not consistent with the previous solution, achieved through a frequency-domain analysis of directivity by Ben-Menahem and Toksöz. The solution asks for a shorter rupture of 270 km at a velocity of 3.5 km/s. The focal parameters are constrained by updating all the reported first motion and are found to be: strike = 103°, dip = 53°, slip = 32°. A seismic moment of 1.8 · 10²⁸ dyn · cm is obtained. These figures are also consistent with a time-domain analysis of Love waves at Palisades and Strasbourg, and of Rayleigh waves at Pasadena, with a directivity study of Love waves at Pasadena, and with static deformation and isoseismal data. A discussion is given of the relation between moment, magnitude and rupture area, and a comparison is made with other events in the same region: it is concluded that this earthquake does not exhibit an “intra-plate” behavior, but rather compares better with “inter-plate” events, such as the great Assam earthquake.     

Seismic activity and deep conductivity structure of the Eastern Carpathians

We present results of a study of the seismicity and the geoelectric structure of the Eastern Carpathians. After the evaluation of the seismicity, new methods of processing and analyzing seismic data are developed, which allow constructing an averaged horizontal-layered velocity model of the crust in the Carpathian region of Ukraine, tracing the seismic active faults and localizing the seismic events both in horizontal and in vertical direction with a better precision. For the study of the conductivity structure beneath the Eastern Carpathians, the collected magnetovariation and magnetotelluric data are used. The depth of electrical conductivity anomalies are estimated and the resulting quasi-3D model of the conductivity structure beneath the Carpathians is compared with the seismicity in the depth range of 10 ± 2.5 km. The comparison suggests possible geological mechanisms: the seismic events occur mainly in resistive solid rock domains which surround aseismic high-conductivity zones, consisting of at least partially melted material. Aqueous fluids, or a joint effect of several mechanisms, may also play an active role in this distribution.     

Earthquake Source Zones in Northeast India: Seismic Tomography, Fractal Dimension and b Value Mapping

We have imaged earthquake source zones beneath the northeast India region by seismic tomography, fractal dimension and b value mapping. 3D P-wave velocity (Vp) structure is imaged by the Local Earthquake Tomography (LET) method. High precision P-wave (3,494) and S-wave (3,064) travel times of 980 selected earthquakes, m _d ≥ 2.5, are used. The events were recorded by 77 temporary/permanent seismic stations in the region during 1993–1999. By the LET method simultaneous inversion is made for precise location of the events as well as for 3D seismic imaging of the velocity structure. Fractal dimension and seismic b value has been estimated using the 980 LET relocated epicenters. A prominent northwest–southeast low Vp structure is imaged between the Shillong Plateau and Mikir hills; that reflects the Kopili fault. At the fault end, a high-Vp structure is imaged at a depth of 40 km; this is inferred to be the source zone for high seismic activity along this fault. A similar high Vp seismic source zone is imaged beneath the Shillong Plateau at 30 km depth. Both of the source zones have high fractal dimension, from 1.80 to 1.90, indicating that most of the earthquake associated fractures are approaching a 2D space. The spatial fractal dimension variation map has revealed the seismogenic structures and the crustal heterogeneities in the region. The seismic b value in northeast India is found to vary from 0.6 to 1.0. Higher b value contours are obtained along the Kopili fault (~1.0), and in the Shillong Plateau (~0.9) The correlation coefficient between the fractal dimension and b value is found to be 0.79, indicating that the correlation is positive and significant. To the south of Shillong Plateau, a low Vp structure is interpreted as thick (~20 km) sediments in the Bengal basin, with almost no seismic activity in the basin.     

S-wave velocity structure in Tangshan earthquake region and its adjacent areas from joint inversion of receiver functions and surface wave dispersion*

Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas (39°N–41.5°N, 115.5°E–119.5°E), we conducted a nonlinear joint inversion of receiver functions and surface wave dispersion. We obtained some detailed information about the Tangshan earthquake region and its adjacent areas, including sedimentary thickness, Moho depth, and crustal and upper mantle S-wave velocity. Meanwhile, we also obtained the v_P/v_S structure along two sections across the Tangshan region. The results show that: (1) the Moho depth ranges from 30 km to 38 km, and it becomes shallower from Yanshan uplift area to North China basin; (2) the thickness of sedimentary layer ranges from 0 km to 3 km, and it thickens from Yanshan uplift region to North China basin; (3) the S-wave velocity structure shows that the velocity distribution of the upper crust has obvious correlation with the surface geological structure, while the velocity characteristics of the middle and lower crust are opposite to that of the upper crust. Compared with the upper crust, the heterogeneity of the middle and lower crust is more obvious; (4) the discontinuity of Moho on the two sides of Tangshan fault suggests that Tangshan fault cut the whole crust, and the low v_S and high v_P/v_S beneath the Tangshan earthquake region may reflect the invasion of mantle thermal material through Tangshan fault.     

Earthquakes Along the Passive Margin of Greenland: Evidence for Postglacial Rebound Control

—?An intriguing observation in Greenland is a clear spatial correlation between seismicity and deglaciated areas along passive continental margins, a piece of evidence for earthquake triggering due to postglacial rebound. Another piece of evidence for induced seismicity due to deglaciation derives from earthquake source mechanisms. Sparse, low magnitude seismicity has made it difficult to determine focal mechanisms from Greenland earthquakes. On the basis of two normal faulting events along deglaciated margins and from the spatial distribution of epicenters, earlier investigators suggested that the earthquakes of Greenland are due to postglacial rebound. This interpretation is tested here by using more recent data. Broadband waveforms of teleseismic P waves from the August 10, 1993 (m _b = 5.4) and October 14, 1998 (m _b = 5.1) earthquakes have been inverted for moment tensors and source parameters. Both mechanisms indicate normal faulting with small strike-slip components: the 1993 event, strike = 348.9°, dip = 41.0°, rake =?56.3°, focal depth = 11?km, seismic moment = 1.03?×?10²⁴ dyne-cm, and M _w = 5.3; the 1998 event, strike = 61.6°, dip = 58.0°, rake =?95.5°, focal depth = 5?km, seismic moment = 5.72?×?10²³ dyne-cm, and M _w = 5.1. These and the two prior events support the theory that the shallow part of the lithosphere beneath the deglaciated margins is under horizontal extension. The observed stress field can be explained as flexural stresses due to removal of ice loads and surface loads by glacial erosion. These local extensional stresses are further enhanced by the spreading stress of continental crust and reactivate preexisting faults. Earthquake characteristics observed from Greenland suggest that the dominant seismogenic stresses are from postglacial rebound and spreading of the continental lithosphere.     

Local seismic network for monitoring of a potential nuclear power plant area

This study presents a plan for seismic monitoring of a region around a potential nuclear power plant. Seismic monitoring is needed to evaluate seismic risk. The International Atomic Energy Agency has set guidelines on seismic hazard evaluation and monitoring of such areas. According to these guidelines, we have made a plan for a local network of seismic stations to collect data for seismic source characterization and seismotectonic interpretations, as well as to monitor seismic activity and natural hazards. The detection and location capability of the network were simulated using different station configurations by computing spatial azimuthal coverages and detection threshold magnitudes. Background noise conditions around Pyhäjoki were analyzed by comparing data from different stations. The annual number of microearthquakes that should be detected with a dense local network centered around Pyhäjoki was estimated. The network should be dense enough to fulfill the requirements of azimuthal coverage better than 180° and automatic event location capability down to ML?～?0 within a distance of 25 km from the site. A network of 10 stations should be enough to reach these goals. With this setup, the detection threshold magnitudes are estimated to be ML?=??0.1 and ML?=?0.1 within a radius of 25 and 50 km from Pyhäjoki, respectively. The annual number of earthquakes detected by the network is estimated to be 2 (ML?≥?～ ?0.1) within 25 km radius and 5 (ML?≥?～?0.1 to ～0.1) within 50 km radius. The location accuracy within 25 km radius is estimated to be 1–2 and 4 km for horizontal coordinates and depth, respectively. Thus, the network is dense enough to map out capable faults with horizontal accuracy of 1–2 km within 25 km radius of the site. The estimation is based on the location accuracies of five existing networks in northern Europe. Local factors, such as seismic noise sources, geology and infrastructure might limit the station configuration and detection and location capability of the network.     

丽江地震台新台址勘选

介绍了丽江地震台新台址勘选过程。使用GL-PS2型一体化地震计进行测点勘选工作,计算分析了2021年12月28日13时至30日12时连续48 h内的地脉动速度均方根值(RMS)和噪声功率谱密度,并对该时段内记录到的天然地震事件进行分析。结果表明,记录到的4次天然地震事件震相清晰易识,拟选台址背景噪声达到Ⅰ类台基噪声水平。综合分析认为,该拟选台址观测环境较好,符合地震台建设要求。     

Evaluating Tsunami Hazard in the Northwestern Indian Ocean

We evaluate here the tsunami hazard in the northwestern Indian Ocean. The maximum regional earthquake calculated from seismic hazard analysis, was used as the characteristic earthquake for our tsunami hazard assessment. This earthquake, with a moment magnitude of M _w 8.3 and a return period of about 1000 years, was moved along the Makran subduction zone (MSZ) and its possible tsunami wave height along various coasts was calculated via numerical simulation. Both seismic hazard analysis and numerical modeling of the tsunami were validated using historical observations of the Makran earthquake and tsunami of the 1945. Results showed that the possible tsunami may reach a maximum height of 9.6 m in the region. The distribution of tsunami wave height along various coasts is presented. We recommend the development of a tsunami warning system in the region, and emphasize the value of education as a measure to mitigate the death toll of a possible tsunami in this region.     

Discrimination of earthquakes and quarry blasts in the eastern Black Sea region of Turkey

In recent years, a large number of quarry blasts have been detonated in the eastern Black Sea region. When these blasts are recorded by seismic stations, they contaminate the regional earthquake catalog. It is necessary to discriminate quarry blast records from the earthquake catalogs in order to determine the real seismicity of the region. Earthquakes and quarry blasts can be separated through different methods. These methods should be applied concurrently in order to safely distinguish these events. In this study, we discriminated quarry blasts from earthquakes in the eastern Black Sea region of Turkey. We used 186 seismic events recorded by the Karadeniz Technical University and Bogaziçi University Kandilli Observatory Earthquake Research Institute stations which are Trabzon, Espiye, Pazar, Borçka, Ayd?ntepe, and Gümü?hane between years of 2002 and 2010. For the discrimination of quarry blasts from earthquakes, we used both, statistical methods (calculation of the maximum ratio of S to P waves (S/P), complexity (C)) and spectral methods (spectrogram calculation). These methods included measuring the maximum amplitude S/P, C, spectral ratio, and time-frequency analysis. We especially relied on two-dimensional time-frequency analysis methods to discriminate quarry blasts from earthquakes in Turkey. As a result of this study, 68 % of the examined seismic events were determined to be quarry blasts and 32 % to be earthquakes. The earthquakes occurring on land are related to small faults and the blasts are concentrated in large quarries. Nearly 40 % of the earthquakes occurred in the Black Sea, most of them are related to the Black Sea thrust belt, where the largest earthquake was observed in the time period studied. The areas with the largest earthquake potential in the eastern Black Sea region are in the sea.