Local site effect of a clay site in Shiraz based on seismic hazard of Shiraz Plain

Propagation of seismic waves through soil layers would drastically change the frequency content and amplitude-based features of ground motions at the surface. These alterations are known as seismic site effects. Computation of site effects of high-populated areas such as large cities is of great importance (e.g., it is used in development of seismic microzonation of a region). Shiraz is one of the most populous cities of Iran and is located in a high seismic hazardous region. A representative clay site in this city is selected to assess local site effects. The time series and random vibration theory procedure in the frequency domain are implemented to analyze the aforementioned site. Furthermore, the nonlinear dynamic soil behavior is simulated by the equivalent linear method and the nonlinear method via DEEPSOIL program. Three types of soil column uncertainties such as shear wave velocity, modulus reduction, and damping ratio of soil layers as well as depth of underlying rock half-space (D _bed) are considered herein. The mean amplification and standard deviation of natural logarithm of amplification factors are computed for a variety of analysis types. The results of the current study show that the computed mean and standard deviation of amplification factor in ln units by considering only V _S uncertainty are in good agreement with the corresponding ones by considering V _S and modulus reduction and damping ratio variabilities simultaneously for the studied site. Furthermore, it seems that the effect of bedrock depth in definition of spectral shapes of the Iranian seismic building code should be taken into account.     

探地雷达数据的S变换时频分析

S变换自问世以来,凭借其优越性已经被广泛地应用于数字信号处理中。对模拟的探地雷达数据和实测的探地雷达数据进行S变换时频分析。其中探地雷达数据的模拟采用时域有限差分法,地电模型为层状介质,通过S变换时频分析,可以更有效地区分薄层状介质。对实测的探地雷达数据进行S变换时频分析,去除了部分噪音干扰,突出有效信号,图像更直观,易于图像的解译。     

时频～炮检距域面波衰减方法及应用

面波以低频强能量为特征,数据处理中常用的面波衰减技术有:频率滤波、F-K滤波和速度滤波等.由于面波特殊的频散现象,使得上述方法的应用效果受到一定的影响.基于面波的特征随炮检距变化这一特点,提出一种时频～炮检距域面波衰减方法:①对单炮记录中的面波区域做Fourier变换或时频分析,了解面波频率随时间和炮检距的变化规律;②...     

Surface wave studies for shear wave velocity and bedrock depth estimation over basalts

Shear wave velocity (V _S) estimation is of paramount importance in earthquake hazard assessment and other geotechnical/geo engineering studies. In our study, the shear wave velocity was estimated from ground roll using multichannel analysis of surface wave (MASW) technique making use of dispersive characteristics of Rayleigh type surface waves followed by imaging the shallow subsurface basaltic layers in an earthquake-prone region near Jabalpur, India. The reliability of MASW depends on the accurate determination of phase velocities for horizontally traveling fundamental mode Rayleigh waves. Inversion of data from surface waves resulted in a shear wave velocity (V _S) in the range of 200–1,200 m/s covering the top soil to weathering and up to bedrock corresponding to a depth of 10–30 m. The P-wave velocity (V _P) obtained from refraction seismic studies at these locations found to be comparable with V _S at an assumed specific Poisson’s ratio. A pair of selected set of V _S profiles over basalt which did not result in a hazardous situation in an earthquake of moderate magnitude are presented here as a case study; in other words, the shear wave velocity range of more than 200 m/s indicate that the area is highly unlikely prone to liquefaction during a moderate or strong earthquake. The estimated depth to basalt is found to be 10–12 m in both the cases which is also supported by refraction studies.     

Source parameters of earthquakes in the reservoir-triggered seismic (RTS) zone of Koyna–Warna, Western India

New empirical relations are derived for source parameters of the Koyna–Warna reservoir-triggered seismic zone in Western India using spectral analysis of 38 local earthquakes in the magnitude range M _L 3.5–5.2. The data come from a seismic network operated by the CSIR-National Geophysical Research Institute, India, during March 2005 to April 2012 in this region. The source parameters viz. seismic moment, source radius, corner frequency and stress drop for the various events lie in the range of 10¹³–10¹⁶ Nm, 0.1–0.4 km, 2.9–9.4 Hz and 3–26 MPa, respectively. Linear relationships are obtained among the seismic moment (M ₀), local magnitude (M _L), moment magnitude (M _w), corner frequency (fc) and stress drop (?σ). The stress drops in the Koyna–Warna region are found to increase with magnitude as well as focal depths of earthquakes. Interestingly, accurate depths derived from moment tensor inversion of earthquake waveforms show a strong correlation with the stress drops, seemingly characteristic of the Koyna–Warna region.     

Probabilistic seismic hazard analysis and spectral accelerations for United Arab Emirates

In recent years, the United Arab Emirates (UAE) has experienced an unprecedented growth which is coupled with the increase in seismic activity in the surroundings. Previous studies presents significant variations in their results whereas some recent studies although very detailed focus on only few cities. This study reviews the results of previous studies and presents new findings for the whole of UAE based on the improved source model and use of next generation attenuation (NGA) equations. The peak ground accelerations, spectral accelerations and deaggregation of hazard for major cities are presented. Moreover, the breakdown of the range of mapped spectral accelerations (S _0.2 and S ₁) is proposed to form the basis for the development of site amplification factors in subsequent studies. The results of this study indicate almost similar values of ground motion compared to some recently published studies and smaller values compared to some earlier studies.     

Site classification of Pondicherry using shear-wave velocity and horizontal-to-vertical spectral ratio

Site classification studies play a vital role in earthquake hazard assessment since in situ ground conditions substantially affect the characteristics of incoming seismic waves during earthquakes. Flat areas along the coast and rivers generally consist of thick layers of soft clay and sand. Such deposits amplify certain frequencies of ground motion, thereby attributing to an increase in the damage due to an earthquake. Hence, site classification studies have been carried out using shear-wave velocity, ground response, and corresponding amplification at 83 locations in Pondicherry, a coastal city in India. The present study is aimed at estimating the shear-wave velocity through multichannel analysis of surface waves and to compute the average shear-wave velocity (V _S ³⁰ ), stiffness, and N values using empirical relations. Further, site-response studies (horizontal-to-vertical spectral ratio) were conducted to estimate the ground-response frequencies and corresponding amplifications through Nakamura technique. From the results, the study area was classified into three types, i.e., C-class: with V _S ³⁰ in the range of 360–760 m/s, D-class: with V _S ³⁰ in the range of 180–360 m/s, and E-class: with V _S ³⁰  < 180 m/s following the National Earthquake Hazard Reduction Programme norms (BSSC in NEHRP recommended provisions for seismic regulations for new buildings and other structures (FEMA 450), part 1: provisions. Building Seismic Safety Council for the Federal Emergency Management Agency, Washington, 2003). Finally, a site classification map for Pondicherry region has been prepared, which can be used in urban planning and strengthening of existing structures against future earthquakes.     

GIS-based multi-criteria earthquake hazards evaluation using analytic hierarchy process for a nuclear power plant site,west Alexandria,Egypt

Nuclear power plants are designed to prevent the hazardous effects of the earthquakes and any external events to keep the safety of the plant. Ninety-one shallow seismic refraction profiles were performed to determine shear wave velocity of the engineering layers at the site of El Dabaa area that is situated to the northern coastline of Egypt for seismic hazard microzonation evaluation according to hazard index values. A microzonation is a procedure of delineating an area into individual zones having different ranks of numerous seismic hazards. This will aid in classifying areas of high seismic risk which is vigorous for industrial design of nuclear structures. The site response analysis requires the characterization of subsurface materials considering local subsurface profiles of the site. Site classification of the area under investigation was undertaken using P- and S-waves and available borehole data. The studied nuclear power plant site has been characterized as per NEHRP site classification using an average velocity of transverse wave (V _s ³⁰ ) of depth 30 m which acquired from seismic survey. This site was categorized into two site classes: the major one is “site class B,” and the minor one is “site class A.” The attenuation coefficient, the damping ratio and the liquefaction potential are geotechnical parameters which were derived from P- and S-waves, and have their major effects on the seismic hazard contribution. 1D ground response analysis was carried out in the places of seismic profiles inside the site for estimating the amount of ground quaking using peak ground acceleration (PGA), site amplification, predominant frequency and spectral accelerations on the surface of ground by the DEEPSOIL software package. Seven factors (criteria) deliberated to assess the earthquake hazard index map are: (1) the peak ground acceleration at the bedrock, (2) the amplification of the site, (3) the liquefaction potential, (4) the main frequency of the earthquake signal, (5) the average V _s of the first 30 m from the ground surface, (6) the depth to the groundwater and (7) the depth to the bedrock. These features were exemplified in normalized maps after uniting them to 0–1 scores according to some criteria by the minimum and maximum values as linear scaling points. Multi-criteria evaluation is an application of multi-criteria decision analysis theory that used for developing a seismic hazard index map for a nuclear power plant site at El Dabaa area in ArcGIS 10.1 software. Two models of decision making were used in this work for seismic hazard microzonation. The analytic hierarchy process model was applied to conduct the relative weights of the criteria by pairwise comparison using Expert Choice Software. An earthquake hazard index map was combined using Weighted Linear Combination model of the raster weighted overlay tool of ArcGIS 10.1. The results indicated that most of the study site of the nuclear power plant is a region of low to moderate hazard; its values are ranging between 0.2 and 0.4.     

Probabilistic seismic hazard analysis of Tripura and Mizoram states

A probabilistic seismic hazard analysis for the states of Tripura and Mizoram in North East India is presented in this paper to evaluate the ground motion at bedrock level. Analyses were performed considering the available earthquake catalogs collected from different sources since 1731–2010 within a distance of 500 km from the political boundaries of the states. Earthquake data were declustered to remove the foreshocks and aftershocks in time and space window and then statistical analysis was carried out for data completeness. Based on seismicity, tectonic features and fault rupture mechanism, this region was divided into six major seismogenic zones and subsequently seismicity parameters (a and b) were calculated using Gutenberg–Richter (G–R) relationship. Faults data were extracted from SEISAT (Seismotectonic atlas of India, Geological Survey of India, New Delhi, 2000) published by Geological Survey of India and also from satellite images. The study area was divided into small grids of size 0.05° × 0.05° (approximately 5 km × 5 km), and the hazard parameters (rock level peak horizontal acceleration and spectral accelerations) were calculated at the center of each of these grid cells considering all the seismic sources within a radius of 500 km. Probabilistic seismic hazard analyses were carried out for Tripura and Mizoram states using the predictive ground motion equations given by Atkinson and Boore (Bull Seismol Soc Am 93:1703–1729, 2003) and Gupta (Soil Dyn Earthq Eng 30:368–377, 2010) for subduction belt. Attenuation relations were validated with the observed PGA values. Results are presented in the form of hazard curve, peak ground acceleration (PGA) and uniform hazard spectra for Agartala and Aizawl city (respective capital cities of Tripura and Mizoram states). Spatial variation of PGA at bedrock level with 2 and 10 % probability of exceedance in 50 years has been presented in the paper.     

Decomposing of mixed pattern of arsenic using fractal model in Gangdese belt,Tibet, China

Decomposing mixed geochemical patterns is a challenge in geochemical exploration and environmental assessment. In this paper, the spectrum–area technique (S–A) is used to decompose a mixed pattern of arsenic in Gangdese belt based on stream sediment data. S–A is a multifractal model based on power–law relationships between area of the set consisting of wave numbers with spectral energy density above S[A(>S)] on the 2D frequency domain. The original spatial distribution map of arsenic obtained by inverse distance weighted (IDW) shows a mixed pattern due to superposition of different geological processes or events and is converted into the frequency domain by means of Fourier transformation. Two components, including power spectrum density and phases, are obtained. The spectrum energy density (S) and the area (A) enclosed by the above-threshold spectrum energy density is plotted on a log–log scale. Two cutoff values determined by three straight lines define three filters which decompose the original map of arsenic into background, anomalous, and high frequency (noise) components. The areas with high anomaly of arsenic mainly are located surrounding known Cu deposits, indicating that arsenic anomalies may be related to Cu mineralization.     

Soil height randomness influence on seismic response: Case of an Algiers site

The paper deals with the effect of the depth to bedrock randomness on the seismic response of an Algiers site in time and frequency domains. Stochastic soil profile seismic analysis is carried out via Monte Carlo simulations coupled to the stiffness matrix method. The soil height is assumed to be a random variable with a log-normal distribution. The stochastic behavior of extreme ground acceleration and its response spectrum, transfer function, fundamental frequency, maximum amplitude, short- and mid-period amplification factors are derived from 1000 samples, as a parametric study is performed accounting for the influence of the coefficient of variation of the depth to bedrock. As the study herein considers a SH wave propagation pattern (SH stands for shear horizontal), the accelerations under study correspond to the E–W and N–S horizontal directions. The seismic acceleration corresponds to the Boumerdes earthquake (Algeria, May 21, 2003, M_w = 6.5). The probability density functions of the output parameters are derived using the maximum entropy principle, and compared to the log-normal distribution.The soil height heterogeneity causes an increase in the fundamental frequency of the soil profile, as well as a wider frequency content, in such a way that a larger number of structures are concerned by the resonance phenomenon.     

时频域高分辨地震层序识别

地震信号频率随着时间的不稳定变化包含了地下反射层序特征的重要信息,时频分析是检测地震信号局部谱特征的重要工具。为适应实际地震信号处理,对S变换改进后得到一种新的广义S变换,用于分析和补偿地震信号的高频成分,得到了精细的地震层序识别剖面。实际资料处理表明,它对砂岩油气藏的成层特征的分析具有分辨率高和高信噪比的优点。地震信号的广义S变换时频谱分解可作为地震层序识别和解释的重要补充。     

Correlation of seismic intensity with Fourier acceleration spectra

We present a method for estimating the seismic intensity in terms of MMI or MSK scale using Fourier amplitude spectra of ground acceleration. The method implies that severity of earthquake ground motion is determined by spectral amplitudes in relatively narrow frequency band: so-called “representative frequencies”, at decreasing frequencies (from 7–8 Hz for small intensities to 0.7 – 1.0 Hz for MMI(MSK) = VIII–IX) with increasing intensity level. It is examined through estimation of probable intensity at a site using recordings of recent earthquakes in several seismic regions and prediction of intensity distribution patterns for some earthquakes. Seismic hazard maps, in terms of intensity levels based upon the proposed approach, should describe regional features of seismic waves excitation and propagation, as well as local ground conditions.     

广义S变换多频解释技术及其在薄层评价中的应用

为了利用时频转换后地震资料完成高速围岩屏蔽薄层定量评价,应用广义S变换对目标地震开展解释性处理,落实薄层时频域地震可识别程度。结合实钻砂体发育数据、频谱分析、薄层时频特征等,开展时频转换数据优选,提取频率梯度属性定性分析目标区沉积储层平面展布特征。构建实钻数据与频率梯度数学关系,完成薄层定量评价。对中东地区实际资料处理解释表明,基于广义S变换的多频解释技术对于高速围岩屏蔽的薄层评价具有较好的应用效果,提高了储层预测精度,为利用地震资料开展储层定量评价提供了一套技术思路,具有广阔的应用前景。     

Groundwater Systems and Eco-hydrological Features in the Main Karst Regions of China

1 Introduction Karst topography is widespread in China (Lu et al., 1973; Lu, 1986a), but its basic features are quite different in the main karstified areas in southern and northern China (Lu, 1999, 2003). The karst development appears to differ not only in the karst geomorphological types but also in the forms of the different karst water systems. The eco-hydrological features have obviously been influenced by human activities. The eco-hydrogeology is thus of great importance to study. Th…     

Geological variation in S-wave velocity structures in Northern Taiwan and implications for seismic hazards based on ambient noise analysis

Ambient noise analysis in Northern Taiwan revealed obvious lateral variations related to major geological units. The empirical Green’s functions extracted from interstation ambient noise were regarded as Rayleigh waves, from which we analyzed the group velocities for period from 3 to 6 s. According to geological features, we divided Northern Taiwan into seven subregions, for which regionalized group velocities were derived by using the pure-path method. On average, the group velocities in mountain areas were higher than those in the plain areas. We subsequently inverted the S-wave velocity structure for each subregion down to 6 km in depth. Following the analysis, we proposed the first models of geology-dependent shallow S-wave structures in Northern Taiwan. Overall, the velocity increased substantially from west to east; specifically, the mountain areas, composed of metamorphic rocks, exhibited higher velocities than did the coastal plain and basin, which consist of soft sediment. At a shallow depth, the Western Coastal Plain, Taipei Basin, and Ilan Plain displayed a larger velocity gradient than did other regions. At the top 3 km of the model, the average velocity gradient was 0.39 km/s per km for the Western Coastal Plain and 0.15 km/s per km for the Central Range. These S-wave velocity models with large velocity gradients caused the seismic waves to become trapped easily in strata and, thus, the ground motion was amplified. The regionalized S-wave velocity models derived from ambient noises can provide useful information regarding seismic wave propagation and for assessing seismic hazards in Northern Taiwan.     

Locating, monitoring, and characterizing typhoon-linduced landslides with real-time seismic signals

Landslides induced by typhoon Morakot during its passage across Taiwan on 7–9 Aug 2009 claimed more than 700 lives and caused heavy economic loss. Unlike earthquake monitoring, precise locations of landslides could not be determined in near-real time because their seismic phases are difficult to identify. Here, we show that large, damaging landslide events are characterized seismically by a distinct waveform pattern of frequent intermixes of P and S waves over a time window of several tens of seconds. The predominant frequency band during these time windows ranges from 0.5 to 5?Hz. The high-frequency content is clearly deficient relative to that of local earthquakes by about one to two orders. We also demonstrate that large landslide events can be located and monitored with algorithms specifically designed for real-time seismic applications. This near-real-time monitoring capability would be particularly useful for emergency responders and government organizations to coordinate effective relief-and-rescue operations.     

Estimation of velocity function parameters in unconsolidated sands using semblance velocity analysis

To properly understand seismic wave propagation in unconsolidated sand layers, it is important to estimate the parameters of their continuous velocity–depth functions. This study proposes a procedure to estimate the V ₀ and k parameters of a specific velocity function, where V ₀ is the direct P-wave velocity at the ground surface and k is the velocity gradient. The V ₀ and k parameters are generally independent of each other. However, it is possible to relate them numerically because both depend strongly on the porosity (?) and water saturation (S _w). The proposed procedure starts by tabulating V ₀ and k for 0.05?≤???≤?0.5 sampled at Δ??=?0.05 and S _w?=?0.6, so that only V ₀ is needed for fitting. Then, time–distance (T-X) type curves of the direct arrival are calculated for the corresponding values of V ₀ and k parameters values. The type curves are fitted then to the observed shot gather through a modification of the classic semblance velocity analysis method. Once the best-fit V ₀ value is found, the corresponding k, ?, and S _w values are picked from a V ₀–k–? lookup table. The procedure is applied on synthetic shot gathers with various amounts of additive Gaussian random noise. Results show that the method is robust and tolerant to low to moderate amounts of noise.     

Seismic velocity and Poisson’s ratio tomography of the crust beneath East Anatolia

Eastern Anatolia is a region in the early stages of continent–continent collision and so provides a unique opportunity to study the early development of continental plateau. Located within the Alpine–Himalayan fold-thrust fault belt, the Anatolian plateau is geologically very complex, with over half of the surface area covered with late Cenozoic volcanics of diverse composition. The plateau is also seismically active and is dissected by numerous seismogenic faults predominantly of strike-slip motion. In this study, we determine 3-D tomographic images of the crust beneath eastern Anatolia by inverting a large number of arrival time data of P- and S-waves. From the obtained P- and S-wave velocity models, we estimated the Poisson’s ratio structures for a more reliable interpretation of the obtained velocity anomalies. Our tomographic results are generally consistent with the major tectonic features of the region. High P- and S-wave velocity anomalies are recognized near the surface, while at deeper crustal layers, low seismic wave velocities are widely distributed. Poisson’s ratio exhibits significant structural heterogeneities compared to the imaged velocity structure. The seismic activity is intense along highly heterogeneous zones and is closely associated with pre-existing faults in the central and western parts of the study area. Results of the checkerboard resolution test indicate that the imaged anomalies are reliable features down to a depth of about 40 km. The low-velocity/high Poisson’s ratio zones in the middle to lower crust are consistent with many geophysical observations such as strong Sn attenuation, low Pn and Sn velocity, and the absence of mantle lid, implying the presence of partial melt in the uppermost mantle.     

Effect of soil layering on the ground response parameters: a parametric study

Soil layering has a pivotal role on the behavior and propagation of seismic waves, hence the ground response during seismic loading. Parametric study to estimate the effect of soil layering on ground response parameters is of prime importance considering the engineering significance of structure founded in seismic zone; nonetheless, it is yet to be well understood. The objective of this study was to investigate the effect of soil layering and soil properties on the ground response parameters. One-dimensional linear ground response analysis was conducted with variation in soil layer parameters including impedance ratio (?? _z ), layer thickness (d _r ), and damping ratio (D). The acceleration time history of the Bhuj Earthquake (M _w ?=?7.7), India, was used in the analysis. The results obtained from the analysis were presented as variation of ground response parameters such as spectral acceleration (SA), amplification ratio (M) with the soil layer parameters. Results showed higher values of SA at lower D and then decreased with increase in D, that in fact depict the resistance offered to the particle oscillation at comparatively higher values of D. Similarly, variation in SA and M was very less or negligible when the ?? _z was varied from 1 to 3 and the d _r equal to 0.2, while for d _r greater than 0.2 the variation increased with ?? _z and d _r . The outcome from the parametric study presented in this paper clearly demonstrates the significance of ?? _z , D, and d _r of the soil layers on the ground response parameters.