Development and application of a GIS-based tool for earthquake-induced hazard prediction

An integrated GIS-based tool (GTIS) was constructed to estimate site effects related to the earthquake hazards in the Gyeongju area of Korea. To build the GTIS for the study area, intensive site investigations and geotechnical data collections were performed and a walk-over site survey was additionally carried out to acquire surface geo-knowledge data in accordance with the procedure developed to build the GTIS. For practical applications of the GTIS used to estimate the site effects associated with the amplification of ground motion, seismic microzoning maps of the characteristic site period and the mean shear wave velocity to a depth of 30 m were created and presented as a regional synthetic strategy addressing earthquake-induced hazards. Additionally, based on one-dimensional site response analyses, various seismic microzoning maps for short- and mid-period amplification potentials were created for the study area. Case studies of seismic microzonations in the Gyeongju area verified the usefulness of the GTIS for predicting seismic hazards in the region.     

Assessing shear wave velocity profiles using multiple passive techniques of Shillong region of northeast India

Shillong region of northeast India falls under seismic zone V. Being a commercial hub, urbanization in this study region is at its peak. In order to qualitatively assess the subsurface velocity profiling of this area, we have blended two robust techniques, namely spatial autocorrelation (SPAC) and frequency wavenumber (FK) method. Corresponding to array noise data collected at five sites, situated in the close proximity of boreholes, we have evaluated V_S and V_P as well. The shear wave velocity estimates yielded by these techniques are found to be consistent with each other. The computed V_s values up to depth of 30 m are observed to be in the range of 275–375 m/s, in most of the sites which implies prevalence of low-velocity zone at some pocket areas. The estimates so found are systematically analyzed and implications are outlined. The results were corroborated by substantial evidence of site geology as well as geotechnical information.     

Liquefaction Hazard Assessment of Earth Quake Prone Area: a Study Based on Shear Wave Velocity by Multichannel Analysis of Surface Waves (MASW)

The shear wave velocity (V_S) profile based on the dispersive characteristics of fundamental mode of Rayleigh type surface waves indicate underground stiffness change with depth as well as near surface stiffness. The most important utility of shear wave velocity (V_S) is to estimate the liquefaction hazard potential of an area particularly in seismically active region. Rayleigh type surface waves were utilized to estimate the velocity (V_S) of shallow subsurface covering a depth range of 30–50 m employing multichannel analysis of surface waves. The liquefaction hazard map predicts an approximate percentage of an area that will have surface manifestation of liquefaction during an earth quake. The surface wave data acquired in an earth quake prone region of Jabalpur (Seismic zone III), India, yields a velocity (V_S) range of 200–750 m/s corresponding to the subsurface depth of 30–35 m. The results were analyzed for possible liquefaction hazard in the study area and presented here besides the N values.     

Spatial distribution of soil shear-wave velocity and the fundamental period of vibration – a case study of the Saguenay region,Canada

The spatial distribution of soil shear-wave velocity and the fundamental period of vibration were selected as input parameters for the determination of potential seismic site effects in the Saguenay region, Canada. The methodology used in this study involved three clear steps. First, a 3D geological model of the surficial deposits was built taking into consideration the type, spatial distribution and thickness of the deposits. Second, representative average V_s values were determined for each of the major soil units. Finally, the average shear-wave velocity from the ground surface to bedrock (V_sav), the shear-wave velocity of the upper 30 m (V_s30) and the fundamental site resonance period (T₀) were calculated over a regular grid for the study area. The results include the spatial distribution of the fundamental site resonance period, the average shear-wave velocity in the first 30 m of the ground and the spatial distribution of National Building Code of Canada seismic soil classes for the Saguenay region.     

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.     

Seismic microzonation for Penang using geospatial contour mapping

Shear wave velocity (V _s) and the fundamental site period of the subsurface condition are the primary parameters that affect seismic soil amplification in particular sites. Within the topmost layer of the soil, which measures 30 m, the average shear wave velocity V _s30 is commonly used to build codes for site classification for the design of earthquake-resistant structures and to conduct microzonation studies. In this study, the development of a microzonation map for V _s30 distribution, National Earthquake Hazard Reduction Program V _s30 site classification, and a fundamental site period for Penang are presented. The multichannel analysis of surface wave (MASW) test was conducted for more than 50 sites with available borehole data to develop the microzonation maps. The ten selected V _s profiles measured by MASW show a good correlation with the data obtained using empirical correlations in a previous study. The highest V _s values were identified at the northeastern and southeastern parts of Penang Island, corresponding to the shallow bedrock and the outcrop zone. Conversely, the lowest V _s values were found in the northwestern and southwestern parts of the Penang mainland owing to the thick layer of soft clay and silt deposits. The site period map shows the variation in site periods, with the highest value of 1.03 s at the western part of the Penang mainland and the lowest value of 0.02 s at the eastern part of the Penang Island. The microzonation maps developed in this study are vital to studies on seismic hazard and earthquake mitigation programs in Malaysia.     

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.     

Preliminary results for a semi-automated quantification of site effects using geomorphometry and ASTER satellite data for Mozambique,Pakistan and Turkey

Estimation of the degree of local seismic wave amplification (site effects) requires precise information about the local site conditions. In many regions of the world, local geologic information is either sparse or is not readily available. Because of this, seismic hazard maps for countries such as Mozambique, Pakistan and Turkey are developed without consideration of site factors and, therefore, do not provide a complete assessment of future hazards. Where local geologic information is available, details on the traditional maps often lack the precision (better than 1:10,000 scale) or the level of information required for modern seismic microzonation requirements. We use high-resolution (1:50,000) satellite imagery and newly developed image analysis methods to begin addressing this problem. Our imagery, consisting of optical data and digital elevation models (DEMs), is recorded from the ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) sensor system. We apply a semi-automated, object-oriented, multi-resolution feature segmentation method to identify and extract local terrain features. Then we classify the terrain types into mountain, piedmont and basin units using geomorphometry (topographic slope) as our parameter. Next, on the basis of the site classification schemes from the Wills and Silva (1998) study and the Wills et al (2000) and Wills and Clahan (2006) maps of California, we assign the local terrain units with V _s 30 (the average seismic shear-wave velocity through the upper 30m of the subsurface) ranges for selected regions in Mozambique, Pakistan and Turkey. We find that the applicability of our site class assignments in each region is a good first-approximation for quantifying local site conditions and that additional work, such as the verification of the terrain’s compositional rigidity, is needed.     

Multichannel analysis of surface waves (MASW) for seismic site characterization using 2D genetic algorithm at Bahrah area,Wadi Fatima,Saudi Arabia

The present study deals with the seismic site classification of Bahrah area, Wadi Fatima, to characterize the local site conditions. The dynamic behavior of sediments was studied by the application of surface wave inversion. The multichannel analysis of surface waves (MASW) shallow geophysical technique was utilized for site classification. MASW survey was carried out at 66 sites along with 13 seismic refraction profiles at suitable localities. MASW and seismic refraction profiles were processed and compared with the available borehole data. The integration of MASW techniques with seismic refraction and borehole data progressively enhanced the subsurface visualization and reliability of the shear wave velocity estimation in the subsurface in the study area. The subsurface shear-wave velocity model was achieved by the solution of an inverse problem-based dispersion of surface waves and propagation in a vertically heterogeneous medium. The 2D genetic algorithm was employed for the inversion of dispersion curves to obtain velocity and thickness of subsurface layers. The depth to engineering bedrock and velocity of shear waves in the first 30 m was deciphered and mapped. The depth of bedrock in study area varies from 4 to 30 m, and V _{S 30} ranges from 320 to 800 m/s. The most of study area falls in B and C class categories in addition to few sites of D class according to the NEHRP guidelines.     

Variations in shear wave velocity and soil site class in Kolkata city using regression and sensitivity analysis

The detrimental effects of an earthquake are strongly influenced by the response of soils subjected to dynamic loading. The behavior of soils under dynamic loading is governed by the dynamic soil properties such as shear wave velocity, damping characteristics and shear modulus. Worldwide, it is a common practice to obtain shear wave velocity (V _s in m/s) using the correlation with field standard penetration test (SPT) N values in the absence of sophisticated dynamic field test data. In this paper, a similar but modified advanced approach has been proposed for a major metro city of eastern India, i.e., Kolkata city (latitudes 22°20′N–23°00′N and longitudes 88°04′E–88°33′E), to obtain shear wave velocity profile and soil site classification using regression and sensitivity analyses. Extensive geotechnical borehole data from 434 boreholes located across 75 sites in the city area of 185 km² and laboratory test data providing information on the thickness of subsoil strata, SPT N values, consistency indices and percentage of fines are collected and analyzed thoroughly. A correlation between shear wave velocity (V _s) and SPT N value for various soil profiles of Kolkata city has been established by using power model of nonlinear regression analysis and compared with existing correlations for other Indian cities. The present correlations, having regression coefficients (R ²) in excess of 0.96, indicated good prediction capability. Sensitivity analysis predicts that significant influence of soil type exists in determining V _s values, for example, typical silty sand shows 30.4 % increase in magnitude of V _s as compared to silt of Kolkata city. Moreover, the soil site classification shows Class D and Class E category of soil that exists typically in Kolkata city as per NEHRP (Recommended provisions for seismic regulations for new buildings and other structures—Part 1: Provisions. Prepared by the Building Seismic Safety Council for the Federal Emergency Management Agency (Report FEMA 450), Washington, DC, 2003) guidelines and thereby highlighting the seismic vulnerability of the city. The results presented in this study can be utilized for seismic microzonation, ground response analysis and hazard assessment for Kolkata city.

     

Seismic microzonation of Bangalore,India

In the present study, an attempt has been made to evaluate the seismic hazard considering local site effects by carrying out detailed geotechnical and geophysical site characterization in Bangalore, India to develop microzonation maps. An area of 220 km², encompassing Bangalore Mahanagara Palike (BMP) has been chosen as the study area. Seismic hazard analysis and microzonation of Bangalore are addressed in three parts: in the first part, estimation of seismic hazard is done using seismotectonic and geological information. Second part deals with site characterization using geotechnical and shallow geophysical techniques. In the last part, local site effects are assessed by carrying out one-dimensional (1-D) ground response analysis (using the program SHAKE2000) using both standard penetration test (SPT) data and shear wave velocity data from multichannel analysis of surface wave (MASW) survey. Further, field experiments using microtremor studies have also been carried out for evaluation of predominant frequency of the soil columns. The same has been assessed using 1-D ground response analysis and compared with microtremor results. Further, the Seed and Idriss simplified approach has been adopted to evaluate the soil liquefaction susceptibility and liquefaction resistance assessment. Microzonation maps have been prepared with a scale of 1:20,000. The detailed methodology, along with experimental details, collated data, results and maps are presented in this paper.     

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.     

Geostatistical modelling of spatial and depth variability of SPT data for Bangalore

The purpose of this study is to develop a geostatistical model to evaluate the spatial and depth variability of Standard Penetration Test (SPT) data from Bangalore, India. The database consists of 766 boreholes spread over a 220 km² area, with several SPT values (N) in each of them. The geostatistical analysis is done for N corrected (N _corrected) values. The N _corrected value has been corrected for different parameters such as overburden stress, size of the bore hole, type of sampler, hammer energy and length of the connecting rod. The knowledge of the semivariogram of the SPT data is used with kriging theory to estimate the values at points in the subsurface of Bangalore where field measurements are not available. The model is used to compute the variance of estimated data. The model predicts reasonably well the SPT data. The geostatistical model provides valuable results that can be used for seismic hazard analysis, site response and liquefaction studies for the development of microzonation maps. The predicted N values from the developed model can also be used to estimate the subsurface information, allowable bearing pressure of soils and elastic modulus of soils.     

S-waves profiles from noise cross correlation at small scale

Cross correlations of ambient seismic noise recordings are studied to infer shear-wave velocities (V_S) profiles versus depth. Experiments are performed, over inter-station distances ranging from 60 m to about 400 m, in the crowded and noisy historical centre of Napoli, where the conditions for the use of active (controlled source) seismic spreadings are prohibitive, even for just one receiver. From the noise cross correlation, group velocity dispersion curves are extracted, with FTAN method, and their non-linear inversion supplies V_S profiles versus depth. The information of nearby stratigraphies and the range of V_S variability for samples of Neapolitan soils and rocks confirms the validity and reliability of the results obtained with our expeditious procedure. If successfully tested in other geological settings, the proposed approach is a low cost methodology to get reliable V_S velocity profiles versus depth.     

A newly developed seismic microzonation model of Erbaa (Tokat, Turkey) located on seismically active eastern segment of the North Anatolian Fault Zone (NAFZ)

A methodology to model seismic microzonation maps is required in the hazard mitigation decision plans of the earthquake prone areas. The stage of disaster preparedness for new residential places is of great importance for detailed seismic microzonation models. The effects of local geological and geotechnical site conditions were considered in order to establish site characterization as the initial stage of the models in this study. Dynamic soil properties based on the empirical correlations between shear wave velocity (V _s) and standard penetration test blow counts were taken into account in order to define representative soil profiles extending down to the engineering bedrock. One-dimensional site response analyses were performed to analyze earthquake characteristics on the ground surface. The layers for soil classification, geology, depth to groundwater level, amplification, distance to fault, slope and aspect, and liquefaction-induced ground deformation potential of the study area were prepared in seismic microzonation models. The study area, Erbaa, is placed along the seismically active North Anatolian Fault Zone. Final seismic microzonation map of the study area was evaluated applying different GIS-based Multi-Criteria Decision Analysis (MCDA) techniques. Two of the MCDA techniques, simple additive weighting and analytical hierarchical process (AHP), are considered during the evaluation step of the final seismic microzonation map. The comparison is made in order to distinguish two different maps based on these MCDA techniques. Eventually, AHP-based seismic microzonation map is more preferable for the seismic design purposes in this study.     

Lower crustal fluid distribution in the northeastern Japan arc revealed by high-resolution 3D seismic tomography

The Ou Backbone Range strikes northwards through the central northeastern Japan arc and is bounded on both sides by the active reverse Uwandaira and Sen'ya faults. We have applied a traveltime inversion method (seismic tomography) with spatial velocity correlation to active and passive seismic data in order to investigate a three-dimensional (3-D) velocity structure. The data set contains 33,993 P- and 18,483 S-wave arrivals from 706 natural sources and 40 blasts, as well as 2803 P-wave traveltime data from 10 explosions detonated during the 1997 controlled source experiment. The traveltime inversion reveals a zone beneath the Ou Backbone Range in which P-wave velocities (V_P) are approximately 6–8% lower than the average velocity at equivalent depths. The low V_P and a low V_P to S-wave velocity (V_S) ratio (V_P/V_S) of about 1.65 suggest the presence of aqueous fluids in the middle crust.     

Estimation of shearwave velocity in drifts using multichannel analysis of surface wave (MASW) technique — A case study from Jammu & Kashmir,India

Multichannel analysis of surface waves (MASW) is a non-destructive seismic prospecting method utilizing Rayleigh waves for imaging and characterizing shallow sub-surface structure. Multichannel analysis of surface waves (MASW) studies were conducted in drift areas of two bridge sites in the hilly terrain of J&K for imaging and characterizing shallow sub-surface structure. The purpose of the present study is to estimate the shear wave velocity (V_S) and subsurface structure in four drifts made in a hilly terrain for construction of two bridges. Rayleigh waves are having dispersive properties, travelling along or near the ground surface and are usually characterized by relatively low velocity, low frequency, and high amplitude. The study area comprises of Tertiary group of rocks which are underlain by Siwalik group. The main rock type in the study area is dolomite which has undergone various geological processes like weathering, jointing, fracturing and shearing. MASW data was collected inside four drifts in the mountainous terrain of J&K state which are located on either sides of Chenab river. The data was analyzed by relevant processing software using dispersion and inversion technique. Shear wave velocities were estimated up to 30 m depth. Average shear wave velocity (V_S ³⁰) up to top 30m was also computed. It is observed that, V_S in the range 400–800 m/s upto 10–15 m corresponding to weathered rock, followed by compact dolomite rock up to the depth of about 30 m with V_S in the range 1200–1600 m/s. Some low velocity zones are also identified from these sections which represent shear zones.     

Geospatial contour mapping of shear wave velocity for Mumbai city

Shear wave velocity is one of the most important input parameter in the analysis of geotechnical earthquake engineering problems, particularly to estimate site-specific amplification factor and ground response study. Dynamic in situ tests such as spectral analysis of surface waves (SASW) or multichannel analysis of surface waves (MASW) are very expensive. Also due to lack of specialized personnel, these tests are generally avoided in many soil investigation programs. Worldwide, several researchers have developed correlations between the SPT ‘N’ value and shear wave velocity ‘V _s’, which are useful for determining the dynamic soil properties. In the present study, more than 400 numbers of soil borehole data were collected from various geotechnical investigation agencies, government engineering institutes and geotechnical laboratories from different parts of Mumbai city, which is financial capital of India with highest population density. In this paper, an attempt has been made to develop the correlation between the SPT ‘N’ value and shear wave velocity ‘V _s’ for various soil profile of Mumbai city and compared with other existing correlations for different cities in India. Using Geographical Information System (GIS), a geospatial contour map of shear wave velocity profile for Mumbai city is prepared with contour intervals of 25 and 50 m/s. The scarcity of database or maps of shear wave velocity profile for Mumbai city will make the present geospatial contour maps extremely useful and beneficial to the designer, practitioners for seismic hazard study involved in geotechnical earthquake engineering.     

Lateral lithospheric conductivity patterns,from European synchronous observatory data

Much information on the regional lithospheric structure may come from MTS data acquired by synchronous 2D arrays and processed with regard to the nonlocal response of a laterally inhomogeneous subsurface. We suggest to invert the nonlocal MT responses applying correlation of all surface horizontal and vertical components of the geomagnetic field recorded simultaneously at all stations. The inversion algorithm has been applied to 2004–2005 European observatory data of diurnal S_q variations for first five harmonics and yielded lateral conductivity patterns for different periods. The maps show spatial correlation between conductivity maxima and lithospheric thickness minima and, specifically, highlight the contours of the Pannonian basin, where lithosphere is as thin as ～50 km, from seismic data.     

Comparison of P- and S-wave velocity profiles obtained from surface seismic refraction/reflection and downhole data

High-resolution seismic-reflection/refraction data were acquired on the ground surface at six locations to compare with near-surface seismic-velocity downhole measurements. Measurement sites were in Seattle, WA, the San Francisco Bay Area, CA, and the San Fernando Valley, CA. We quantitatively compared the data in terms of the average shear-wave velocity to 30-m depth (V_s30), and by the ratio of the relative site amplification produced by the velocity profiles of each data type over a specified set of quarter-wavelength frequencies. In terms of V_s30, similar values were determined from the two methods. There is <15% difference at four of the six sites. The V_s30 values at the other two sites differ by 21% and 48%. The relative site amplification factors differ generally by less than 10% for both P- and S-wave velocities. We also found that S-wave reflections and first-arrival phase delays are essential for identifying velocity inversions. The results suggest that seismic reflection/refraction data are a fast, non-invasive, and less expensive alternative to downhole data for determining V_s30. In addition, we emphasize that some P- and S-wave reflection travel times can directly indicate the frequencies of potentially damaging earthquake site resonances. A strong correlation between the simple S-wave first-arrival travel time/apparent velocity on the ground surface at 100 m offset from the seismic source and the V_s30 value for that site is an additional unique feature of the reflection/refraction data that could greatly simplify V_s30 determinations.