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.

     

Simultaneous inversion of the aftershock data of the 1993 Killari earthquake in Peninsular India and its seismotectonic implications

The aftershock sequence of the September 30th, 1993 Killari earthquake in the Latur district of Maharashtra state, India, recorded by 41 temporary seismograph stations are used for estimating 3-D velocity structure in the epicentral area. The local earthquake tomography (LET) method of Thurber (1983) is used. About 1500P and 1200S wave travel-times are inverted. TheP andS wave velocities as well asV _P/V_Sratio vary more rapidly in the vertical as well as in the horizontal directions in the source region compared to the adjacent areas. The main shock hypocentre is located at the junction of a high velocity and a low velocity zone, representing a fault zone at 6–7 km depth. The estimated average errors ofP velocity andV _P/V_Sratio are ±0.07 km/s and ±0.016, respectively. The best resolution ofP and S-wave velocities is obtained in the aftershock zone. The 3-D velocity structure and precise locations of the aftershocks suggest a ‘stationary concept’ of the Killari earthquake sequence.     

Deterministic seismic microzonation of Kolkata city

This paper presents the deterministic seismic microzonation of densely populated Kolkata city situated on the world’s largest delta island with very soft and thick soil deposit in the surficial layers. A fourth-order accurate staggered-grid finite-difference algorithm for SH-wave propagation simulation in visco-elastic medium is used for the linear computation of ground motion amplifications in sedimentary deposit. Different maps such as for fundamental frequency (F ₀), peak ground acceleration (PGA), peak ground velocity, and peak ground displacement are developed for variety of end-users communities, including structural and geotechnical engineers for performance-based designs, building officials, emergency managers, land-use planners, private businesses, and the general public. The scenario of simulated amplification factors in the different frequency bands revealed that the Kolkata city is very much prone to severe damage even during a moderate earthquake and very selective damage may occur at some of the localities during local and distant earthquakes. The deterministically predicted PGA at bedrock level is 0.0844 g and the maximum PGA predicted at the free surface is 0.6 g in Kolkata city due to maximum credible earthquake (M _w = 5.4) associated with Eocene Hinge Zone at a depth of 36 km. The seismic microzonation of Kolkata city reveals that the Nager Bazar and Nimtala areas are the safest regions with earthquake point of view.     

An Empirical Correlation of Index Geomechanical Parameters with the Compressional Wave Velocity

The geomechanical strength of rockmass plays a key role in planning and design of mining and civil construction projects. Determination of geomechanical properties in the field as well as laboratory is time consuming, tedious and a costly affair. In this study, density, slake durability index, uniaxial compressive strength (UCS) and P-wave velocity tests were conducted on four igneous, six sedimentary and three metamorphic rock varieties. These properties are crucial and used extensively in geotechnical engineering to understand the stability of the structures. The main aim of this study is to determine the various mechanical properties of 13 different rock types in the laboratory and establish a possible and acceptable correlation with P-wave velocity which can be determined in the field as well as laboratory with ease and accuracy. Empirical equations were developed to calculate the density, slake durability index and UCS from P-wave velocities. Strong correlations among P-wave velocity with the physical properties of different rock were established. The relations mainly follow a linear trend. Student’s ‘t’ test and ‘F’ test were performed to ensure proper analysis and validation of the proposed correlations. These correlations can save time and reduce cost during design and planning process as they represent a reliable engineering tool.     

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.     

Soil void ratio correlation with shear wave velocities and SPT N values for Indo-Gangetic basin

In this study attempt has been made to understand in-situ void ratio in Indo-Gangetic basin (IGB) and to form empirical relations between void ratio and shear wave velocity (V_s), N values considering subsoil investigation data. Multichannel analysis of surface wave (MASW) test and standard penetration test was carried out along with soil property measured at 25 locations. The general soil profile varied from silty sand to clay of low compressibility, ground water level fluctuated between 1-27 m, depth of borehole varied from 20-40 m. Regression analysis was conducted on 202 data sets of void ratio and shear wave velocity, 293 data sets of void ratio and SPT- N value, which resulted in inverse correlations between void ratio and V_s, SPT N value. The datas were segregated into fine, coarse grained data based on engineering classification and relations were developed separately. Until now, no studies have related in-situ void ratio to V_s and SPT N. These correlations will be useful to predict void ratio for sites having measured values of V_s and N value. These void ratios can be further used to assess liquefaction susceptibility.     

Soil liquefaction potential induced by the andalusian earthquake of 25 December 1884

On 25 December 1884, an earthquake of epicentral intensityI ₀ = IX in the MSK scale caused great damage in a large area in the provinces of Granada and Málaga, in the south of Spain. The reports of the Spanish, Italian and French Commissions that studied the earthquake described ground phenomena in seven different sites which can be identified as soil liquefaction.By means of dynamic penetration tests carried out in the above sites, the corresponding soil profiles (based on SPT data and water table depth) were established, and the occurrence of liquefaction was proved in five out of seven of these sites. Also, the intensities at such locations and the magnitude of the earthquake were estimated.From the geotechnical data and the cyclic stress ratio induced by the earthquake, liquefaction conditions were confirmed in all the five sites which presumably liquefied. Then, possible values of the minimum ground surface accelerations necessary for the onset of liquefaction at each location were calculated. The results obtained were completed with data reported in six liquefaction case studies from Japan and the United States, from which design charts relating soil acceleration with normalized SPT values for different intensity levels were drawn.Finally, by using standard attenuation curves, the above data were translated into epicentral distances, and good agreement with the known epicentral area was found. As a result, a consistent approach for liquefaction hazard and source location problems has been developed. The proposed method combines in its formulation historical evidence and earthquake engineering techniques.     

Directional properties of small strain shear stiffness in soils

A compiled database of shear wave velocity measurements in a variety of clays, silts and sands shows directional hierarchies between downhole (V_sVH), standard crosshole (V_sHV), and rotary crosshole (V_sHH) tests. The special in situ database has been collected from 33 well-documented geotechnical test sites. Expressions relating the small-strain shear modulus in terms of effective confining stress level, stress history and void ratio are explored for each of these three modes of directional shear wave velocity. The relationships are separated initially into soil groups (intact clays, fissured clays, sands and silts), and then generalised to consider all soil types together.     

Seismic reliability-based analysis and GIS mapping of cyclic mobility of clayey soils of Mumbai city,India

Cyclic mobility is a mechanism of ground failure due to lateral spreading of soils during an earthquake that usually occurs in soft or medium stiff saturated soils. The simplified procedures developed by the researchers give a factor of safety for judging the cyclic mobility potential. However, the simplified procedures do not take into account the uncertainty in the parameters required to estimate the cyclic stresses in the soil. In this study, a reliability framework based on the simplified procedure, considering the parameter uncertainty, has been proposed for computing the probability of cyclic mobility of clay deposits for a metro city of India, i.e., Mumbai city (latitudes 18°53′N–19°19′N and longitudes 72°47′E–72°58′E). Extensive geotechnical borehole data from 1028 boreholes across 50 locations in the city area of 390 km² and laboratory test data are collected and analyzed thoroughly. A correlation between undrained shear strength (Su) and other parameters such as natural water content (w), SPT N value, liquid limit (LL) and plasticity index (PI) has been established for Mumbai city and has been used in the proposed approach. The sensitivity analysis of the proposed approach predicts that Su has significant influence in the evaluation of the cyclic mobility. Cyclic mobility hazard maps are prepared using the geo-statistical analysis tool in GIS, and it shows that the clayey soils at few locations have a 60–90 % probability of cyclic mobility for a moment magnitude (M _w) of an earthquake of 7.5. These hazard maps can be used by the geotechnical engineers for the cyclic mobility hazard assessment of Mumbai city.     

Visualization of geotechnical data by means of geographic information system: a case study in Eskisehir city (NW Turkey)

Recently, geographic information system Software’s are widely used in engineering geology applications. Such a study was performed in the Eskisehir urban area of about 30 km² via geo-engineering maps, cross-sections and panel-block diagrams, which are formed to predict the properties of foundation soil. In order to determine the engineering geology properties of the region, physical and mechanical characteristics of the soil profile were tested on 1,011 disturbed and 383 undisturbed samples from a total of 170 boreholes. Topographic, lithological and geotechnical data of the region were compiled as a database in computer. Five geological cross-sections were drawn to decipher specifically the alluvial deposits in detail in densely populated areas. Hence, geological characteristics of the area were evaluated in a three-dimensional mode. Five standard penetration test (SPT) cross-sections at the same directions and SPT, uniaxial compressive strength (qu) and soil class Zonation maps were produced at different depths and panel-block diagrams were also created to define the geo-engineering properties of the investigation area. In this respect, the sub-surface of the study area is fully visualized and useful predictor maps, cross-sections and panel-block diagrams are performed to be used in further studies.     

Probabilistic seismic hazard of Pakistan, Azad-Jammu and Kashmir

The seismic hazard study for Pakistan and Azad Jammu and Kashmir has been conducted by using probabilistic approach in terms of peak ground acceleration (PGA) in m/s² and also seismic hazard response spectra for different cities. A new version of Ambraseys et al. (Bull Earthq Eng 3:1–53, 2005) ground acceleration model is used, and parameterization is based on most recent updated earthquake catalogs that consisted of 14,000 events. The threshold magnitude was fixed at M _w 4.8, but seismic zones like northern Pakistan–Tajikistan, Hindukush and northern Afghanistan–Tajikistan border had M _w 5.2. The average normalized ‘a’ and ‘b’ values for all zones are 6.15 and 0.95, respectively. Seismicity of study area was modeled, and ground motion was computed for eight frequencies (0.025, 0.1, 0.2, 0.5, 1.0, 1.5, 2.0, 2.5 s) for different annual exceedance rates of 0.02, 0.01, 0.005, 0.002 and 0.001 (return periods 50, 100, 200, 500 and 1,000 years) for stiff rocks at the gridding of 0.1° × 0.1°. Seismic hazard maps based on computed PGA for 0.02, 0.01 and 0.002 annual exceedance are prepared. These maps indicate the earthquake hazard of Pakistan and surrounding areas in the form of acceleration contour lines, which are in agreement with geological and seismotectonic characteristics of the study area. The maximum seismic hazard values are found at Muzaffarabad, Gilgit and Quetta areas.     

A genetic algorithm approach for assessing soil liquefaction potential based on reliability method

Deterministic approaches are unable to account for the variations in soil’s strength properties, earthquake loads, as well as source of errors in evaluations of liquefaction potential in sandy soils which make them questionable against other reliability concepts. Furthermore, deterministic approaches are incapable of precisely relating the probability of liquefaction and the factor of safety (FS). Therefore, the use of probabilistic approaches and especially, reliability analysis is considered since a complementary solution is needed to reach better engineering decisions. In this study, Advanced First-Order Second-Moment (AFOSM) technique associated with genetic algorithm (GA) and its corresponding sophisticated optimization techniques have been used to calculate the reliability index and the probability of liquefaction. The use of GA provides a reliable mechanism suitable for computer programming and fast convergence. A new relation is developed here, by which the liquefaction potential can be directly calculated based on the estimated probability of liquefaction (P _L ), cyclic stress ratio (CSR) and normalized standard penetration test (SPT) blow counts while containing a mean error of less than 10% from the observational data. The validity of the proposed concept is examined through comparison of the results obtained by the new relation and those predicted by other investigators. A further advantage of the proposed relation is that it relates P _L and FS and hence it provides possibility of decision making based on the liquefaction risk and the use of deterministic approaches. This could be beneficial to geotechnical engineers who use the common methods of FS for evaluation of liquefaction. As an application, the city of Babolsar which is located on the southern coasts of Caspian Sea is investigated for liquefaction potential. The investigation is based primarily on in situ tests in which the results of SPT are analysed.     

Soil liquefaction studies at Mumbai city

Mumbai city is the economical capital of India and is situated about midway on the western coast of stable continental region of Peninsular India. Major part of the city being of reclaimed land, the soil type is of alluvium, sand, and recent conglomerate. There are some bigger water bodies within the city range. In this study, an attempt has been made to study the susceptibility of soil liquefaction using simplified empirical procedure based on number of blow counts (N values) of the soil layers from standard penetration test. The liquefaction susceptibility is quantified in terms of factor of safety along the borehole depths at available borehole locations using earthquake-induced cyclic stress on the soil and the cyclic resistance of the soil to withstand the load. The factor of safety against liquefaction is evaluated at different sites for two peak ground acceleration (PGA) levels pertaining to 10 and 2?% probability of exceedance in 50?years corresponding to uniform hazard response spectra for Mumbai city with 475- and 2,475-year return period, respectively. Contour maps are prepared that display the factor of safety at different depths for earthquake magnitude of M _w 6.5. These contour maps show the liquefaction vulnerability at different sites in the city.     

A characterisation of the physical properties of soil and the implications for landslide occurrence on the slopes of Mount Elgon,Eastern Uganda

Soil properties of major landslides that occurred recently on the mid-altitude slopes of Mount Elgon, eastern Uganda were analysed. A mudflow, located at the Kitati protected forest site, and two deep debris flows on the Nametsi and Buwabwala deforested steep slopes (36°–58°) were surveyed. In order to test the hypothesis that ‘soils at the landslide sites are particularly ‘problem soils’ and thus prone to landslides’, the following analyses were undertaken: particle size distribution, Atterberg limits, shear strength and factor of safety (F_s). Soils at the Kitati and Buwabwala sites exhibited expansive potential, owing to clay contents well above 20%. A clay content exceeding 32% was identified at the Nametsi debris flow site implying an extremely high expansive potential of the soil. High liquid limits (LLs) at Kitati (59%) and Buwabwala (53%) meant that the soils qualified as vertisols susceptible to landslides. High plasticity indices (PIs) (averaging 33%) also confirmed the vertic nature of soils at the Nametsi debris flow site. Whereas the value of F _s  < 1 for the Kitati site signifies an inherently unstable slope, Nametsi and Buwabwala are supposedly stable slopes (F _s  > 1). Despite this finding, the stable sites could be described as only conditionally stable because of the interplay of various physical, pedological and anthropogenic factors. The results point to the fact that soils at the landslide sites are inherently ‘problem soils’ where slope failure can occur even without human intervention. Therefore, the hypothesis that soils at three landslide sites are inherently ‘problem soils’ and prone to landslides, is accepted.     

Advanced performance in geotechnical engineering using tomography analysis

This paper presents the technique to demonstrate reliable two-dimensional (2-D) tomography of near-surface soil through multichannel analysis of surface wave (MASW) method. The MASW method is used as seismic method for determining the shear wave velocity profile of near-surface soil with better performance. Normally, this method represents one-dimensional (1-D) soil layer profile to delineate anomalous subsurface materials and detect soil characteristics in geotechnical investigations. Conventionally, the MASW method is able to represent depth corresponding shear wave velocity through 1-D representation. In this paper, this method is improved through representation of depth and distance versus shear wave velocity profile using 2-D tomography analysis in geotechnical investigations. The outcomes of the newly developed method are more reliable and informative in comparison to the outcome of the conventional MASW method. The significance of this research is incorporating advanced tomography technique with MASW method to obtain 2-D tomography of geotechnical characteristics with consistency.     

Use of Surface Waves in Statistical Correlations of Shear Wave Velocity and Penetration Resistance of Chennai Soils

Shear wave velocity (V _s) is one of the most important input parameter to represent the stiffness of the soil layers. It is preferable to measure V _s by in situ wave propagation tests, however it is often not economically feasible to perform the tests at all locations. Hence, a reliable correlation between V _s and standard penetration test blow counts (SPT-N) would be a considerable advantage. This paper presents the development of empirical correlations between V _s and SPT-N value for different categories of soil in Chennai city characterized by complex variation of soil conditions. The extensive shear wave velocity measurement was carried out using Multichannel Analysis of Surface Waves (MASW) technique at the sites where the SPT-N values are available. The bender element test is performed to compare the field MASW test results for clayey soils. The correlations between shear wave velocity and SPT-N with and without energy corrections were developed for three categories of soil: all soils, sand and clay. The proposed correlations between uncorrected and energy corrected SPT-N were compared with regression equations proposed by various other investigators and found that the developed correlations exhibit good prediction performance. The proposed uncorrected and energy corrected SPT-N relationships show a slight variation in the statistical analysis indicating that both the uncorrected and energy corrected correlations can predict shear wave velocity with equal accuracy. It is also found that the soil type has a little effect on these correlations below SPT-N value of about 10.     

A Simple Limit Equilibrium Approach for Calculation of Ultimate Bearing Capacity of Shallow Foundations on Two-Layered Granular Soils

The bearing capacity of shallow foundations in a non-homogeneous soil profile has been a challenging task in geotechnical engineering. In this paper, a limit equilibrium method is used for calculating bearing capacity factors of shallow foundations constructed on a two-layered granular soil profile. The main objective has been to determine the ultimate bearing capacity computed from equivalent bearing capacity factors N_q and N_γ and comparing that with numerical analysis using finite element methods. It will be shown that the data obtained form the developed method are well comparable with those obtained from FE approach, specially when the difference between shear strength parameters of layers is low which is a practical case for sedimentary soil profiles and also for artificially compacted soils. A computer program has been developed to investigate the influence of various parameters on bearing capacity factors.     

An Experimental Characterization of Shear Wave Velocity (V<Subscript>s</Subscript>) in Clean and Hydrocarbon-Contaminated Sand

The characteristics of hydrocarbon-contaminated soils have been among major concerns of geotechnical engineers due to its significant frequency of event and also its influential consequences on our surroundings from various environmental and engineering viewpoints. Heretofore, the effects of diverse kinds of hydrocarbon contaminants on majority of geotechnical properties of fine- and coarse-grained soils such as grain size, hydraulic conductivity, plasticity, compressibility, internal friction, cohesion, and shear strength have been investigated. However, there has not been a concentrated research study examining shear wave velocity (\({\text{V}}_{\text{s}}\)) of hydrocarbon-contaminated soils as an important geotechnical property of soil due to this fact that, in small/very small strain levels, the maximum shear modulus of soils (\({\text{G}}_{ \hbox{max} }\)) can be determined using shear wave velocity (\({\text{G}}_{ \hbox{max} } =\uprho{\text{V}}_{\text{s}}^{2}\)). This paper aims to investigate effects of hydrocarbon contamination on shear wave velocity of sandy soils by comparing shear wave velocities in identically prepared clean and contaminated samples. To this aim, an Iranian light crude oil, a standard type of silica sand (Ottawa sand), and a bender element apparatus were used to minutely measure shear wave velocity of clean and crude oil contaminated sand samples. Moreover, dry and quasi-moist tamping methods were employed in order to provide comparable clean and contaminated specimens (containing 4, 6, 8, 10, and 12 wt% of crude oil), respectively. Firstly, a comprehensive bender element (BE) and resonant column tests were conducted on the identically prepared clean sand samples at various amounts of frequency (2–20 kHz) and under various confining pressure (50–500 kPa) to find the best methods of accurately determining shear wave travel time in BE tests. Thereafter, BE tests were conducted to examine shear wave velocity in contaminated specimens. Based on the results, it was found that there was a critical value for crude oil content with the maximum shear wave velocity so that shear wave velocity of 4 wt% contaminated sand (V_s-4 wt%) was about 1.2 times higher than clean one (V_s-clean), and contrastingly adding further crude oil up to 6 wt% made a significant reduction in value of shear wave velocity to some extent that V_s-6 wt% was slightly lower than V_s-clean (V_s-6 wt% = 0.95–0.97V_s-clean). Moreover, adding more contaminant (8–12 wt%) into sand had negligible influences on shear wave velocity. In this paper, the effects of crude oil contamination on sand microstructure were also evaluated using scanning electron microscopy.     

基于互相关函数的单孔法波速测试优化算法

土的剪切波速是岩土工程领域的一个重要参数。建筑场地类别划分、地基饱和砂土液化判别、地震反应分析等都依赖于剪切波速,由于其重要性,对它的测定技术和精确程度应重点研究。针对工程中广泛采用的单孔法剪切波速测试,提出了一种基于互相关函数计算土层剪切波速的优化算法。克服了传统数据处理方法中经常遇到的负波速、波速数值明显过大或过小和同一土层各测点波速差别较大等4种不合理因素,提高了数据分析的精度和效率。     

Influence of magnetic fabric anisotropy on seismic wave velocity in paramagnetic granites from NW Himalaya: Results from preliminary investigations

Anisotropy of Magnetic Susceptibility (AMS) and seismic wave velocity studies of some paramagnetic Himalayan granitoids show good correlation between magnetic fabric anisotropy and P wave velocity (Vp). Vp shows strong positive correlation with magnetic lineation (L) and degree of magnetic anisotropy (P′) having correlation coefficient (r) values of 0.93 and 0.89 respectively. Both Vp and Vs show positive correlation with the SiO₂ content of Proterozoic and Paleozoic granitoids. Velocity of S wave (Vs) shows negative correlation with mean magnetic susceptibility (K_m) having ‘r’ value of 0.86. The correlation between Vs-K_m, V_p-P′, V_p-L also shows >95% probability in Spearman’s rank correlation. Based on the results from the present sample size it is suggested that, in paramagnetic granites, V_p is proportional to intensity of deformation and preferred orientation of minerals as well as the mineralogy. On the other hand, Vs is more dependent on the mineralogy alone.