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.     

Determination of shear wave velocity and depth to basement using multichannel analysis of surface wave technique

The dispersive characteristics of Rayleigh type surface waves were utilized to estimate shear wave velocity (Vs) profile followed by imaging the shallow subsurface granitic layers in the heart of Hyderabad. The reliability of Multichannel Analysis of Surface Waves (MASW) depends on the accurate determination of phase velocities for horizontally traveling fundamental mode Rayleigh waves. Multichannel recording leads to effective identification and isolation of various factors of noise. Calculating the 1-D shear wave velocity (Vs) field from surface waves ensures high degree of accuracy irrespective of cultural noise. The main advantage of mapping the bed rock surface with shear wave velocity (Vs) is the insensitivity of MASW to velocity inversion besides being free from many constraints such as contrast in physical properties etc. Modeling of surface waves data results a shear wave velocity (Vs) of 250?C750 m/s covering the top soil to weathering and up to bedrock corresponding to a depth range of 10?C30 m. Further, the computed N values (which is an indicator of site characteristic) based on the harmonic shear wave velocity up to a depth of 5 m is found to be quite high (> 25?C30) well above 5 indicating the site to be safe and strong enough and not prone to liquefaction. A pair of selected set of results over granites are presented here as a case study highlighting the salient features of MASW.     

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.     

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.     

Seismic investigation of near-surface geological structure in the Paducah, Kentucky, area: application to earthquake hazard evaluation

Seismics method were used to evaluate shallow geological conditions at 33 sites in the vicinity of Paducah, Kentucky. A combined set of P- and S-wave seismic refraction and reflection soundings were used, in addition to local borehole information, to produce structure maps of (1) a shallow (< 30 m deep) horizon believed to represent an unconformity surface at the top of the Eocene, and (2) the Paleozoic bedrock surface (< 85 to > 160 m deep). Shear-wave velocity contrasts across the shallow unconformity were generally 2-to-1 while the contrast at the top of the Paleozoic bedrock exceeds 5-to-1. These seismic boundaries have been determined to be very important in modelling and interpreting earthquake ground motion amplification in the Paducah area. The quality and accuracy of the data, and the cost effective nature of the methods, suggest that other communities in areas at risk to damage from seismic activity, with foundation conditions comparable to Paducah, might benefit from similar characterization in order (1) to identify seismically hazardous, near-surface, geological conditions, and (2) to develop geological models that could be used in computer simulations of site response.     

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 soil assessment in the Ubhur area,north of Jeddah,western Saudi Arabia,using a multichannel analysis of surface waves method

Ten boreholes drilled in Ubhur area up to the depth of bedrock indicted the shallow depth of bedrock where the average depth ranges between 10 and 15 m. The standard penetration test N-values of these boreholes were measured and averaged. Based on N-values to the depth of bedrock, Ubhur area can be classified as site class C and D. Multichannel analysis of surface waves technique has been applied along seventy six profiles using 24-channel geophone array and 4.5Hz vertical geophones with 1m geophone spacing and sledgehammer and/or weight drop as seismic energy sources. Values of shear wave velocity to 30 m are calculated and then averaged (Vs30) where it ranges between 310.08 m/s and 1139.8 m/s. Therefore, Ubhur area can be classified into site class B, C and D based on site classification of the national earthquake hazards reduction program (NEHRP) recommendations. Accordingly, the greatest part of the study area falls in site class C while class B and D covered limited areas in the western and the eastern parts respectively. Depending on the shallow depth of bedrock in the study area, the Vs30 parameter is not applicable in the study area so the average values of Vs for the soil thickness, excluding the bedrock, have been calculated and mapped for site class C and D only. So Vs30 approach is not applicable for areas with shallow depth of bedrock which gives higher classification.     

Seismic vulnerability assessment in the new urban area of Diriyah Governorate,Riyadh, Saudi Arabia

Nine seismic refraction profiles were conducted and processed to study the near-surface sediments in the new urban area of Diriyah. The 2D geoseismic models illustrate two layers: a surface layer of soft sediments and weathered to hard limestone bedrock. Moreover, microtremor measurements were performed at 38 sites for 40 min using three-component seismographs and processed to assess the peak spectral amplitude and the corresponding fundamental resonance frequency. The seismic vulnerability index at each measurement site was estimated. These results correlate well with the geotechnical borehole data. The north-western zone is highly vulnerable due to the great thickness of the soft sediments.     

Investigation of the relationship between ground and engineering bedrock at northern part of the Gulf of İzmir by borehole data supported geophysical works

Loss of life and property that may occur as a result of a possible earthquake can be reduced by earthquake resistant building designs. In order to investigate possible ground motion amplification in earthquake resistant building design, relationship between the ground and engineering bedrock must be ensured. In order to provide this relation, structure, basic characteristics, and thickness of the ground are investigated. In this context, calculating ground transfer function, obtaining horizontal earthquake acceleration changes, calculating Vs values and defining the engineering bedrock are necessary. In this study, Menemen plain, the nothern part of Izmir metropolitan located in active earthquake zone and its immediate vicinity have been examined to define the structure, ground, engineering and bedrock relation. In this context, Menemen plain has been investigated by geophysical methods, which are supported with borehole data (microtremor, MASW – multichannel analysis of surface waves, microgravity measurements, and vertical electrical sounding–VES). Microtremor method was conducted at 377 points in average in the investigation area to define fundamental period and empirical transfer function; after that in order to create basin model and to define the shallow subsurface geometry, microgravity measurements were carried out by using Scintrex CG-5. Also, MASW measurements were carried out in approximately 277 profiles and Schlumberger VES measurements were conducted at approximately 7 points in the investigation area. The existence of a linear relation between H/V peak period values obtained by microtremor measurements and ground thickness in the investigation area is also supported by geothermal drilling logs (depth of 600 m) with microgravity survey. Also, in some parts of the investigation area, it was observed that high S velocity (Vs) values affected H/V peak period values in sections of the ground close to the surface and there was an inversely correlated relation between this ground thickness and peak period values. This event occurring in shallow depths is supported by both VES sections and 2nd order vertical gravity derivative. As a result, depth of the engineering bedrock was obtained between 200 and 700 m and this unit was proposed as Bornova Melange for the investigation area in the scope of the works carried out. Also, it is observed that the area from ground to the engineering bedrock consists of four different layers which were defined by individual S velocities and densities. According to all results, characterictics of the shallow subsurface show that there is a high heterogeneity. Therefore, according to Eurocode8 (EC8 2004) regulations, soil characteristic of the Menemen plain and its vicinity are in the S1–S2 soil class.     

Application of MASW in road failure investigation

Multichannel analysis of surface waves (MASW) survey was conducted to measure shear wave velocities in order to ascertain the likely causes of road failure along LASU-IBA expressway in Alimosho local government area, Lagos, Nigeria. MASW data were acquired along the express road. The acquired dataset was processed and transformed into two-dimensional structure reflective of depth and surface wave velocity distribution within a depth of investigation using SurfSeis software. The MASW shear wave velocity data were compared to geophysical data that was acquired along the same profile. The comparison was also done with geotechnical data that had been acquired prior to the study some meters away from the study area. The correlation between N values to measured shear wave velocity using MASW was generated. The comparison illustrates the accuracy and consistency of MASW-derived shear wave velocity profiles. We concluded that (1) the low-velocity region that varies between 100 and 250 m/s at surface down to 4 m beneath the surface is characterized by loose/peat materials and may have been responsible for the road failure within the study area; this region depicts a very loose compaction area. (2) The MASW technique is a time–cost-effective tool for obtaining reliable shear wave velocity profiles, and (3) the MASW is particularly attractive in areas that cannot be readily assessed by other geophysical and geotechnical tools.     

场地类别分类方案研究

基于性能的抗震设计,要求工程师设计出具有预期抗震性能的结构,一个关键因素是地震作用的确定,这在很大程度上取决于局部场地条件。通过收集和分析北京、苏州和唐山城区956个钻孔资料,建立地表20 m和30 m深土层走时平均剪切波速VS20和VS30的关系式;现场钻探获取北京城区深105 m的典型钻孔原状土样,试验给出各类土体动剪切模量和阻尼比曲线;建立北京城区170个钻孔的场地反应计算模型,采用Nakamura提出的HVSR法和陈国兴等提出的弱震法估算场地基本周期TS值,结合国内外现行抗震规范的场地分类及一些学者对场地分类的研究成果,提出两种新的场地分类建议方案：基于等效剪切波速VSE和覆盖土层厚度H（地表至剪切波速VS ≥ 500 m/s的基岩深度）的双指标场地分类方案及基于VSE、H和TS的三指标场地分类方案。提出的场地分类方案对我国现行抗震规范场地分类方法的改进有参考价值。     

Investigation of groundwater potentiality in the eastern part of Qena,central Eastern Desert,Egypt, using 2D homogeneous function method

The present study represents an application of the 2D homogeneous function method to investigate the near-surface groundwater potentiality to the east of Qena town, Egypt. The homogeneous function automatically inverts first arrival refractions to derive a 2D velocity distribution, which involves seismic boundaries. A complex set of observed refracted traveltimes along four seismic profiles (12 spreads) is used to construct refraction velocity field section and structural section. Such sections allow viewing the complex layered structure and delineating its boundaries and faults. The interpreted geoseismic layers were traced, according to their own gradients and velocity ranges. The resultant models were successful in fulfilling the proposed objective, as they provide: (1) the required details on the delineation of the subsurface layers, where three geoseismic layers with different ranges of velocities were detected; (2) determination of the thickness and the inversion boundary of water-bearing layer; and (3) detection of normal and reverse faults with displacement of about 10–40 m and variable dipping directions. The obtained result was found to be completely compatible with the information gained from identified boundaries of the nearby borehole.     

Application of seismic velocity tomography in investigation of karst collapse hazards,Guangzhou, China

This paper describes seismic velocity tomography applied to the investigation and assessment of karst collapse hazards to facilitate accurate characterization of geological conditions of karst sinkhole formation. In the survey areas of Xiamao, Guangzhou, China, and Huangchi, Foshan, China, seismic velocity tomography was used to explore the structures of rock and soil associated with karst collapse. The results show that sand intercalated with clay or clay intercalated with soft soil dominates the cover of these two areas. The overburden is 20–33 m thick and underlain by Carboniferous limestone. In the limestone, there are well-developed karst caves and cracks as well as highly fluctuating bedrock surfaces. The seismic velocities are less than 2500 m/s in the cover, 2500–4500 m/s in the karst fracture zones and caves of Xiamao, and 1500–2000 m/s in the Huangchi collapse area. The karst fracture zones, relief of bedrock surfaces, and variations of soil thicknesses revealed by seismic velocity tomography are well constrained and in agreement with those in the drilling borehole profiles. This paper demonstrates that seismic velocity tomography can delineate anomalies of rock and soil with the advantages of speed, intuitive images, and high resolution.     

A mobile multi-depth borehole sensor set-up to study the surface-to-base seismic transfer functions

The best method to evaluate the seismic site response is by means of borehole vertical arrays that use earthquake records from different depths. In this paper we introduce the implementation of a single borehole sensor system (synchronized to a sensor on the surface) that is fixed at variable depths within a single well. This system is used for recording small amplitude earthquake signals at variable stiffness conditions in depth to compute empirical borehole transfer functions. The computed average empirical borehole transfer functions allow the estimation of an S-wave velocity model that is constrained using the frequency peak observed in the H/V ratio curve.Pairs of surface and borehole earthquake records were obtained with the borehole sensor placed at − 10, −20, −50, and − 100 m depth in a test site in Managua, Nicaragua. The average velocity of the final model down to − 100 m appeared to be in good agreement with the average velocity computed via cross-correlation using the surface and borehole signals. Likewise, an inverted MASW profile and H/V ratio at the same site agree with the S-wave velocity model obtained.     

High-Frequency Rayleigh-Wave Method

High-frequency (≥2 Hz) Rayleigh-wave data acquired with a multichannei recording sys-tem have been utilized to determine shear (S)-wave velocities in near-surface geophysics since the early 1980s. This overview article discusses the main research results of high-frequency surface-wave tech-niques achieved by research groups at the Kansas Geological Survey and China University of Geosciences in the last 15 years. The multichannel analysis of surface wave (MASW) method is a nou-iuvasive acoustic approach to estimate near-surface S-wave velocity. The differences between MASW results and direct borehole measurements are approximately 15% or less and random. Studies show that simultaneous inversion with higher modes and the fundamental mode can increase model resolution and an investigation depth. The other important seismic property, quality factor (Q), can also be estimated with the MASW method by inverting attenuation coefficients of Rayleigh waves. An inverted model (S-wave velocity or Q) obtained using a damped least-squares method can be assessed by an optimal damping vector in a vicinity of the inverted model determined by an objective function, which is the trace of a weighted sum of model-resolution and model-covariance matrices. Current developments include modeling high-frequency Rayleigh-waves in near-surface media, which builds a foundation for shallow seismic or Rayleigh-wave inversion in the time-offset domain; imaging dispersive energy with high resolution in the frequency-velocity domain and possibly with data in an arbitrary acquisition geometry, which opens a door for 3D surface-wave techniques; and successfully separating surface-wave modes, which provides a valuable tool to perform S-wave velocity profiling with high-horizontal resolution.     

Efficiency of seismic attributes in detecting near-surface cavities

The purpose of this study is to evaluate the efficacy of using seismic attributes to detect near-surface cavities. The methods used in this study include interpreting dispersion curves and amplitude mapping of the multichannel analysis of surface wave (MASW) technique and interpreting the delay in first arrivals of compressional waves. To test these methods, a seismic survey was conducted above a known near-surface cavity in Al-Suman Area, Saudi Arabia. The cause of the cavity is carbonization in the area; there are many cavities similar to this one. The seismic data were collected using a seismograph system with 48 vertical geophones. Both techniques show a tangible result for detecting the cavity. The 2D section of shear wave velocity, which was obtained by inverting the dispersion curves from the MASW technique, leads us to determine the shape of the cavity, as described by a low-velocity zone. Frequency against relative offset is plotted and shows a significant frequency drop in the presence of the cavity, which also provides an indication to the presence of cavity underneath. This interpretation is matched by the interpretation of observed delays in first arrivals of compressional waves. The integration of both P wave seismic refraction and MASW gives confidence in the result and matches observations of the existing cavity closely.     

Identification of cross-valley faults in the Maynardville Limestone,Oak Ridge Reservation,Tennessee, using seismic refraction tomography

First arrival times from P-wave refraction and reflection seismic surveys along Bear Creek Valley on the Oak Ridge Reservation, Tennessee, were inverted to produce refraction tomographic velocity images showing seismic velocity variations within thinly mantled karstic bedrock to a depth of approximately 20 m. Inverted velocities are consistent with two distinct bedrock groups: the Nolichucky Shale (2,730–5,150 m/s) and Maynardville Limestone (3,940–7,575 m/s). Low-velocity zones (2,700–4,000 m/s) in the tomographic images correspond to previously inferred cross-valley strike-slip faults; in places, these faults create permeability barriers that offset or block groundwater flowing along Bear Creek Valley. These faults may also force groundwater contaminants, such as dense non-aqueous phase liquids, to migrate laterally or downward, spreading contamination throughout the groundwater system. Other, previously unmapped cross-valley faults may also be visible in the tomographic images. Borehole logs suggest the low-velocity values are caused by low rigidity fractured and vuggy rock, water zones, cavities and collapse features. Surface streams, including Bear Creek, tend to lie directly above these low-velocity zones, suggesting fault and fracture control of surface drainage, in addition to the subsurface flow system. In some cases, fault zones are also associated with bedrock depressions and thicker accumulations of unconsolidated sediment.     

Shear Wave Velocities in the Estimation of Earthquake Hazard Over Alluvium in Seismically Active Region

It is well known that the potential hazard during an earthquake is mainly in alluvium or alluvium filled basins; shear wave velocity plays a significant role in estimating the possible hazard during an earth quake in such an area. This paper presents shear wave velocity profile from Jabalpur, Central India mainly dominated by alluvial soil that was moderately affected by an earthquake of magnitude 6.5 in May, 1997. The acquired shear wave velocity by Multichannel Analysis of Surface Wave (MASW) in as many as 36 sites over alluvial soil ranges from 200 m/sec to 400 m/sec and in a few sites marginally less than 200 m/sec corresponding to a sub-surface depth of 30–35 m. Further, the computed N values vary as low as near zero to less than 25. The study is substantiated by the estimation of P-wave velocity by refraction seismic method at the same locations of MASW which ranges from 350 m/sec to 2200 m/sec. The results suggest that the damage during an earth quake appears to be highly unlikely in view of the marginally high Vs up to depth of 30 m. This study on seismic hazard is substantiated by the estimation of frequency of the ground as well as amplification which is found to be a maximum of about 2.5 in the frequency band of 2–6 Hz in west and north western portion of the study area.     

Subsurface characterization using seismic refraction and surface wave methods: a case of Lagos State University, Ojo, Lagos State

This study contains the finding of geophysical investigations conducted at the proposed science complex site at Lagos State University, Ojo, Lagos, Nigeria. Surface wave and seismic refraction tests are non-invasive seismic techniques and have been used to determine the shear wave velocity profile of soil deposits. The methods provide a simplified characterization of subsurface in two-dimensional (2D) (distance and depth) profiles. Seismic records obtained were processed/analyzed by Seis-Imager software to obtain one-dimensional shear wave velocity (Vs) distribution. Multiple Vs obtained were integrated and used to construct two-dimensional Vs map. The measured P- and S-wave velocities were also used to estimate Poisson’s ratio, rigidity modulus, and N-values. The study had shown that the area investigated composed mainly of loose sediments (clay formation) to the depth of 12 m with P-wave velocity ranging between 125 and 205 m/s and corresponding S-wave velocity between 60 and 100 m/s. The results presented in this study will be vital information for the engineers in construction of the proposed science complex.     

Evaluating the structural model from ambient vibration recordings and shallow seismic refraction techniques. A case study of Zagazig City,Egypt

This study uses the seismic refraction and noise measurements to investigate the velocity structure of the subsurface and emphasize the advantage of ambient vibration over the conventional seismic refraction technique. Field measurements were carried out at nine sites in and around Zagazig city. Shallow seismic refraction data were interpreted using the delay time method to obtain the two-dimension ground model at each site. Ambient vibration arrays are used to infer the one-dimensional compressional and shear wave velocity profiles through two main steps. The first step is to derive the dispersion curve from the recorded signals using the frequency-wavenumber method. The second is to invert the dispersion curve to obtain the site velocity profiles. The results of the compressional wave velocities obtained from seismic refraction technique showed that the subsurface consists of a number of layers ranging from two to four layers and give a good agreement with the results of the seismic wave velocities obtained from the ambient vibration arrays. The ambient vibration arrays gave a deeper depth of penetration than the other method, providing more information on the subsurface structure without any disturbance to the environment. This work provides reliable estimates of the seismic velocity structures of both shallow and deep sedimentary layers within the area of interest.