Comparison of spectral ratio techniques for estimation of site effects using microtremor data and earthquake motions recorded at the surface and in boreholes

The authors apply soil response estimation techniques employing accelerograms for fifteen earthquakes recorded at the Yokohama Strong Motion Network and its vertical array of nine sites, plus microtremor data recorded at all 150 sites. Assessment of the reliability of surface to reference site spectral ratios, and horizontal to vertical spectral ratios of S‐waves, coda, and microtremors, relied on cross‐validation with both surface to borehole spectral ratios, and theoretical soil response functions for vertically incident SH plane waves calculated from the vertical array recordings and logging data, respectively. The results indicate the following: (1) the reference site method provides more accurate soil response estimations; (2) the H/V of microtremors provides relatively better soil response estimations than its S‐wave and coda counterparts; and (3) the H/V of microtremors provides coincident soil response estimations at two‐thirds of the sites when considering no more than ±30% mismatch in the determination of the predominant periods. The matching of parameters obtained from the analysis of microtremors, combined with their relatively low economic cost and the possibility of recording without strict spatial or time restrictions turns microtremors into a particularly appealing approach. Copyright © 2003 John Wiley & Sons, Ltd.     

The Effect of Velocity Inversions on H/V

We analyzed the phenomenology of microtremor H/V curves under inversions in the shear-wave velocity (V_s) profile in the subsoil. Under no V_s inversion the spectral signature of the H/V peaks is found to be ‘eye-shaped’ with the horizontal components higher than the vertical. Conversely, under negative velocity gradients, numerous of differences emerge. I) A H/V ratio below 1 is observed for a wide range of frequencies, due to the decrease of the horizontal components below the vertical one. II) In the presence of persistent H/V < 1, small bumps in the H/V ratio given by local minima in the vertical spectral component may represent the relics of the peaks indicating resonances and stratigraphic discontinuities. As a consequence, in the presence of velocity inversions the H/V > 2 SESAME (2004) criterion fails but a stratigraphic interpretation may still be possible. III) The H/V curves should always be interpreted together with the single component spectra. IV) Microtremor H/V measurements for stratigraphic/microzonation purposes on stiff artificial soils, (asphalt, concrete, cement, pavements) should always be avoided since the latter often produce velocity inversions. This may have consequences in the intermediate to high frequency domain ( > 1 Hz) also in the application of reference site methods, like H_site/H_bedrock, to microtremor. Theoretical modeling confirms these experimental findings.     

Reliability of rayleigh wave dispersion curve obtained from f–k spectral analysis of microtremor array measurement

The paper utilizes previously developed microtremor simulation technique to evaluate the reliability of Rayleigh wave dispersion curve estimated by f–k spectral analysis of microtremor array measurement. The simulated microtremors are obtained for a fictitious (TEST) site. Attempt is also made to obtain the dispersion curve for two real sites (OHDATE and SKC) by inverse analysis of the microtremor array measurement using f–k spectra method. The estimated dispersion curve from simulated microtremors (TEST site) compares well with the theoretical dispersion curve, demonstrating the reliability of f–k spectra method and indicating that the estimated dispersion curve from microtremor measurements could be adequately used as the target for inverse analysis purposes. It is also demonstrated that the dispersion curve from microtremor measurements can be utilized to estimate the soil profile at OHDATE and SKC sites by inverse analysis. Results show that the theoretical dispersion curve of the fundamental mode of Rayleigh wave after the end of inverse analysis are in good agreement with the dispersion curve obtained by f–k spectral analysis of microtremor array measurement.     

Comparison between different techniques for evaluation of predominant periods using strong ground motion records and microtremors in Pereira Colombia

This paper presents a comparison between different techniques for evaluation of predominant periods in soft soil, for the urban area of Pereira city, Western Colombia. In this study we used microtremor and strong ground motion records obtained by a local array of seven accelerographs stations deployed in the city. Response spectra and spectral ratios have been calculated and compared with strong seismic events recorded in solid rock and soft soil stations. These observations allowed the determination of dominant response spectra for several sectors in the urban area. For the microtremor measurements and earthquake data, dominant periods were determined using interpretation of Fourier amplitude spectra and Nakamura's technique. A comparison between dominant periods obtained from strong ground motion records and those obtained from microtremor measurements show similarities, which is in the range 0.2–0.5 s. A preliminary version of a site response map for Pereira city was obtained from this analysis.     

P- and S-wave velocities of consolidated sediments from a seafloor seismic survey in the North Celtic Sea Basin, offshore Ireland

A geophysical survey was conducted over a hydrocarbon prospect in the North Celtic Sea Basin using a small array of ocean‐bottom seismographs (OBSs). The purpose of this study was to determine the ratio of compressional (P)‐ to shear (S)‐wave velocity of consolidated sedimentary rocks in order to constrain possible subsurface variations in pore‐fluid content. The ratio of V_P and V_S is known to be particularly sensitive to lithology, porosity and pore‐fluid content, making it a useful parameter for evaluating hydrocarbon prospects. OBSs offer a relatively cheap and time‐effective means of acquiring multi‐component data compared with ocean‐bottom cables. In this contribution, we demonstrate the ability of an OBS survey comprising three pairs of two OBSs spaced at 1.6 km to recover lateral variations in the V_P/V_S ratio. A key requirement of this type of study is that S waves will be generated by mode conversions in the subsurface, since they cannot be generated in nor travel through fluids. In this survey, the contrast in physical properties of the hard seabed of the North Celtic Sea Basin provided a means of generating converted S waves. Two‐dimensional ray‐tracing and forward modelling was used to create both V_P and V_S models along a profile crossing the Blackrock prospect in the North Celtic Sea Basin. These models comprise four layers and extend to a maximum depth of 1.1 km. The observed northward decrease in the V_P/V_S ratio at depths of 500–1000 m below the seafloor in the study area is interpreted to represent lateral variation in the amount of gas present in the pore space of Upper Cretaceous chalks and shales overlying the prospective reservoir.     

Characterisation of ground motion recording stations in the Groningen gas field

The seismic hazard and risk analysis for the onshore Groningen gas field requires information about local soil properties, in particular shear-wave velocity (V_S). A fieldwork campaign was conducted at 18 surface accelerograph stations of the monitoring network. The subsurface in the region consists of unconsolidated sediments and is heterogeneous in composition and properties. A range of different methods was applied to acquire in situ V_S values to a target depth of at least 30 m. The techniques include seismic cone penetration tests (SCPT) with varying source offsets, multichannel analysis of surface waves (MASW) on Rayleigh waves with different processing approaches, microtremor array, cross-hole tomography and suspension P-S logging. The offset SCPT, cross-hole tomography and common midpoint cross-correlation (CMPcc) processing of MASW data all revealed lateral variations on length scales of several to tens of metres in this geological setting. SCPTs resulted in very detailed V_S profiles with depth, but represent point measurements in a heterogeneous environment. The MASW results represent V_S information on a larger spatial scale and smooth some of the heterogeneity encountered at the sites. The combination of MASW and SCPT proved to be a powerful and cost-effective approach in determining representative V_S profiles at the accelerograph station sites. The measured V_S profiles correspond well with the modelled profiles and they significantly enhance the ground motion model derivation. The similarity between the theoretical transfer function from the V_S profile and the observed amplification from vertical array stations is also excellent.     

Mapping of underground cavities by the passive seismic standing waves method: the case study of Barsukovskaya cave (Novosibirsk region,Russia)

This article presents the results of mapping a karst cave by the passive seismic standing waves method. Barsukovskaya cave is located about 100 km southeast of the city of Novosibirsk (Russia). The total length of the cave's passages and grottoes is estimated at about 200 m, the maximum depth from the earth's surface is about 19 m. The method for studying underground cavities used is based on the effect of the generation of standing waves by microtremor in the space between the earth's surface and the cave roof. The accumulation of amplitude spectra of a large number of microtremor records makes it possible to determine the frequencies of the first few modes of these waves. Areal passive seismic survey on the earth's surface above the cave made it possible to construct a map of the lowest mode frequency distribution over the cave roof. Since no standing waves were observed at other points, this map reflects the cave structure in plan, which confirms the comparison with the cave diagram drawn up earlier by one of the speleologists. A schematic map of the depth of the cave roof was constructed using the longitudinal wave velocity V_p = 3120 m/s determined by the rock samples selected near the entrance to the cave. This map at a qualitative level also agrees with the data of speleologists, which indicate that the cave, on average, gradually becomes deeper from the entrance to its dead‐end branches. The shallower depths in comparison with the data of speleologists are apparently explained by a very low estimate of the velocity determined from a rock sample taken near the entrance to the cave. The reliability of the obtained cave mapping results is confirmed by the numerical simulation results using the finite‐element method.     

Estimation of sedimentary layer shear wave velocity using micro-tremor H/V ratio measurements for Bangalore city

For assessing earthquake hazard of metro cities, knowledge of soil amplification, thickness and properties of sedimentary layer are essential. In order to map the soil thickness using microtremor survey method, in Bangalore city, it is required to calibrate the relation between fundamental resonance frequency of the soil layer and its thickness for the region. For this purpose microtremor survey was carried out at 34 locations in the city where borehole log was available. The resonance frequency of the soil is evaluated from the microtremor recordings using the H/V ratio technique. A nonlinear regression relation between the thickness of sedimentary layer h (m), from the borehole logs, and the resonance frequency f_r (Hz), was derived as h=(58.3±8.8)f_r^{−(0.95)±0.1}. Using the model of shear wave velocity increasing with depth at these locations, the derived average shear wave velocity and the corresponding soil thickness were used, to get an empirical relation between V_S (m/s) and depth z(m), as V_s=(174±28)(1+z)^0.16±0.07. This relation also compares reasonably with the fit obtained between simulated V_S and depth from borehole logs for Bangalore city. The calibrated relations can be used at locations in Bangalore city where borehole logs are not available, for finding the thicknesses and shear wave velocities of the local soil layers at the survey locations.     

Characteristics of earthquake ground motions and the H/V of microtremors in the southwestern part of Taiwan

The relationships between the spectral characteristics of earthquake ground motions and those of micro‐tremors are investigated using the observed data from a dense strong‐motion network consisting of 108 stations in the Yun‐Li, Chia‐Yi and Tai‐Nan areas in southwestern Taiwan. Many high‐quality recordings, including those of the 921 Chi‐Chi earthquake (M_w=7.6), the 1022 Chia‐Yi mainshock (M_L=6.4), the 1022 major aftershock (M_L=6.0), as well as some weak motion events are selected to evaluate site responses. Microtremor measurements are also performed at most ground motion stations. With many stations in the area located on an alluvium structure, however, it is difficult to find good reference stations on rock sites, which therefore necessitates the calculation of single‐station H/V ratios. The predominant frequencies obtained from H/V ratios are consistent with those from spectral ratios. The site characteristics between the strong and weak events are different, however. This implies that a nonlinear effect probably occurred with the strong‐motion events. The main peak in the H/V spectra of the microtremors is in good agreement with the first peak obtained from the spectra of earthquake ground motions. It is reasonable to claim that the main peak reflects the deep underground structure. On the basis of the H/V ratios of the microtremors, it is concluded that the lower predominant frequencies appear in the plain area, while the higher values are near the mountainous region. Copyright © 2002 John Wiley & Sons, Ltd.     

Detecting underground cavities using microtremor data: physical modelling and field experiment

In this paper, we study the possibilities of the use of microtremor records in the detection and delineation of near‐surface underground cavities. Three‐dimensional physical modelling data showed that the averaging amplitude spectra of a large number of microtremor records makes it possible to evaluate the frequencies and amplitudes of compressional standing waves generated by microtremor in the space between the ground surface and underground cavities. We illustrate how these parameters can be used to estimate the shape of the underground cavity horizontal projection. If the compressional wave velocity in the enclosing rock is known, it is possible to evaluate the depth to the cavity top using the frequencies of the standing waves. The results of the field experiment confirmed the possibility of underground cavities detection using microtremor data.     

Determination of shale stiffness tensor from standard logs

A technique allowing inversion of the shale stiffness tensor from standard logging data: sonic velocities, density, porosity and clay content is developed. The inversion is based on the effective medium theory. The testing of the technique on laboratory measurements of the elastic wave velocities in shale samples shows that the inversion makes it possible to predict the elastic wave velocities V_P, V_S1 and V_S2 in any direction within an error of a few per cent. The technique has been applied for the stiffness tensor inversion along a well penetrating a shale formation of the Mississippian age altered by thin layers of limestone. It is demonstrated that the symmetry of a stiffness tensor inverted at the sonic frequency (2 kHz) is slightly orthorhombic and taking into account the experimental errors, can be related to the vertical transverse isotropy symmetry. For the productive interval of the shale formation, the Thomsen parameters ?, γ, and δ average, respectively, 0.32, 0.25 and 0.21, which indicate anelliptic behaviour of the velocities in this shale. The coefficients of anisotropy of this shale interval are around 24% and 20% for the compressional and shear waves, respectively. The values of the inverted velocities in the bedding plane for this interval are in good agreement with the laboratory measurements. The technique also allows inversion of the water saturation of the formation (Sw) and the inverted values are in agreement with the Sw values available for this formation. A Backus‐like upscaling of the inverted stiffness tensors is carried out for the lower and upper bounds of the frequency band used in the crosswell tomography (100 Hz and 500 Hz). These results can serve as an initial velocity model for the microearthquake location during hydrofracking of the shale formation.     

Coherence methods in mapping AVO anomalies and predicting P-wave and S-wave impedances

Filters for migrated offset substacks are designed by partial coherence analysis to predict ‘normal’ amplitude variation with offset (AVO) in an anomaly free area. The same prediction filters generate localized prediction errors when applied in an AVO‐anomalous interval. These prediction errors are quantitatively related to the AVO gradient anomalies in a background that is related to the minimum AVO anomaly detectable from the data. The prediction‐error section is thus used to define a reliability threshold for the identification of AVO anomalies. Coherence analysis also enables quality control of AVO analysis and inversion. For example, predictions that are non‐localized and/or do not show structural conformity may indicate spatial variations in amplitude–offset scaling, seismic wavelet or signal‐to‐noise (S/N) ratio content. Scaling and waveform variations can be identified from inspection of the prediction filters and their frequency responses. S/N ratios can be estimated via multiple coherence analysis. AVO inversion of seismic data is unstable if not constrained. However, the use of a constraint on the estimated parameters has the undesirable effect of introducing biases into the inverted results: an additional bias‐correction step is then needed to retrieve unbiased results. An alternative form of AVO inversion that avoids additional corrections is proposed. This inversion is also fast as it inverts only AVO anomalies. A spectral coherence matching technique is employed to transform a zero‐offset extrapolation or near‐offset substack into P‐wave impedance. The same technique is applied to the prediction‐error section obtained by means of partial coherence, in order to estimate S‐wave velocity to P‐wave velocity (V_S/V_P) ratios. Both techniques assume that accurate well ties, reliable density measurements and P‐wave and S‐wave velocity logs are available, and that impedance contrasts are not too strong. A full Zoeppritz inversion is required when impedance contrasts that are too high are encountered. An added assumption is made for the inversion to the V_S/V_P ratio, i.e. the Gassmann fluid‐substitution theory is valid within the reservoir area. One synthetic example and one real North Sea in‐line survey illustrate the application of the two coherence methods.     

Seismic hazard evaluation in Anjar city area of western India: Microtremor array measurement

In western India during the Bhuj earthquake (Mw 7.6) on January 26, 2001, the Anjar City at ~30 km southwest of Bhuj experienced three types of damage scenario: severely damaged, less damaged and non-damaged. Similar damage patterns were also observed for the 1819 (Mw 7.8) and the 1956 (Mw 6.0) earthquakes. Microtremor array measurements were conducted in and around the Anjar city to examine the strength of soil structures and damage pattern. Significant differences are observed in frequencies and amplitudes in horizontal-to-vertical spectral ratio (HVSR) using microtremor measurements. The severely- damaged site shows two peak amplitudes: 2.8 at 1.2 Hz; and 4.0 at 8.0 Hz. The less-damaged site also shows two amplitudes: 2.5 and 2.1 at 1.4 Hz; and 2.0 Hz, respectively. The non-damaged site, on the other hand, shows that the HVSR curves become almost flatter. Similar results for three types of damage scenario based on analyses of earthquake records are also observed for the study area. The microtremor array measurements has revealed shear wave velocity V_s≥400 m/s at 18 m depth in the non-damaged, at 40 m in the less-damaged and at 60 m depth in the severely-damaged sites. The site amplitudes and the V_s values show a good correlation with the soil characteristics and damage pattern, suggesting that strength of soil layers at varying depths is a dictating factor for the estimate of the earthquake risk evaluation of the area under study.     

Modelling S-wave radiation for SP-waves in isotropic and anisotropic media

Utilizing shear-wave (S-wave) data acquired with compressional waves (P-waves) is becoming more common as joint imaging and inversion techniques improve. Interest in S-waves radiated from vertical sources and buried explosives exploits conversion to P-waves as primary reflections (SP-waves) for reducing acquisition costs and for application to legacy data. However, recent investigations overstate the extent of SP-wave illumination and show isotropic processing results with narrow bandwidth frequency and wavenumber data. I demonstrate that illumination with SP-waves is limited in general to near vertical polar angles up to around 30° or 35° for V_P/V_S of 2 or 3, respectively. At greater angles, S-waves are typically in the P-wave evanescent range and cannot excite SP-wave reflections. Contrary to recent claims, these sources for P-wave do not radiate S_H-waves polarized in horizontal planes in all azimuths. I show these properties for isotropic media with radiation expressions for amplitude derived in vector slowness coordinates. Also, I extend these expressions to transversely isotropic media with a vertical symmetry axis to show agreement with synthetic seismic data that only quasi S_V-waves are radiated and become more narrowly focused towards 45°. Furthermore, in orthorhombic media, synthetic data show that fast S₁- and slow S₂-waves polarized parallel and perpendicular to fractures may appear as S_V- and S_H-waves. For the partially saturated fracture model studied here, S₁-wave radiation has broader azimuthal illumination than slow S₂-waves, which are more narrowly focused in azimuth. These produce SP-wave splitting signatures on vertical component reflection data that are nearly identical to PS-wave signatures on radial horizontal component data. Separating these fast and slow SP-waves is an additional processing challenge.     

Borehole study of compressional and shear attenuation of basalt flows penetrated by the Brugdan and William wells on the Faroes shelf

We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective Q_S/Q_P ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small.     

Estimation of Shallow S-Wave Velocity Structure in Damascus City,Syria, Using Microtremor Exploration

Array measurements of microtremors were carried out at thirty sites in Damascus city, Syria to estimate S-wave velocity structures of shallow soil formations for site effect analysis. The microtremor data were recorded by 6 vertical-component seismometers distributed along the circumferences of two circles as well as a 3-component seismometer deployed in the center. The phase velocities were estimated at each site from the vertical components of recorded microtremor data by using the Spatial Autocorrelation method. Then, Genetic Simulated Annealing Algorithm technique was applied for inversion of the phase velocities to estimate 1-D S-wave velocity structures beneath the sites. The inverted Vs profiles are not uniform in Damascus city and the results show that a shallow soft layer (∼200 m/s) appears in the eastern part of the city as well as the central part along Barada River. This layer controls the amplification distribution in the city with a high amplification mainly observed at the locations having this layer. The inversion results also show that the depth to the engineering bedrock (∼750 m/s) is very shallow along the foothills of Mt. Qasyoun in the north-west. Then the depth increases towards the east and the south. The maximum depth to the engineering bedrock (∼80 m) was observed in the southern part of Damascus. To validate the results of the inversions, the spectral ratios between the horizontal and vertical components (H/V) of the recorded microtremor data at the central seismometer were compared with the computed ellipticities of the fundamental-mode Rayleigh-waves based on the respective Vs structure. The results show a good agreement in a period range of 0.05 s to 0.5 s. In this period range, the dominant peaks of the H/V ratios are due to the overall effect of the velocity contrasts between the shallow layers representing the subsurface S-wave velocity structure. Moreover, the average S-wave velocity for the top 10 m of soils (V_S10) shows a better correlation with the averaged site amplification in a period range of 0.05 s to 0.5 s than V_S30 which indicates that V_S10 can be a better proxy for high-frequency site amplification in the case of Damascus city.     

Seismic microzonation of the greater Bangkok area using microtremor observations

Bangkok, the capital city of Thailand, is located at a remote distance from seismic sources. However, it has a substantial risk from these distant earthquakes due to the ability of the underlying soft clay to amplify ground motions. It is therefore imperative to conduct a detailed seismic hazard assessment of the area. Seismic microzonation of big cities, like Bangkok, provides a basis for site‐specific hazard analysis, which can assist in systematic earthquake mitigation programs. In this study, a seismic microzonation map for the greater Bangkok area is constructed using microtremor observations. Microtremor observations were carried out at more than 150 sites in the greater Bangkok area. The predominant periods of the ground were determined from the horizontal‐to‐vertical (H/V) spectral ratio technique. A microzonation map was then developed for the greater Bangkok area based on the observations. Moreover, the transfer functions were calculated for the soil profile at eight sites, using the computer program SHAKE91, to validate the results from the microtremor analysis. The areas near the Gulf of Thailand, underlaid by a thick soft clay layer, were found to have long natural periods ranging from 0.8s to 1.2s. However, the areas outside the lower central plain have shorter predominant periods of less than 0.4s. The study shows that there is a great possibility of long‐period ground vibration in Bangkok, especially in the areas near the Gulf of Thailand. This may have severe effects on long‐period structures, such as high‐rise buildings and long‐span bridges. Copyright © 2003 John Wiley & Sons, Ltd.     

Application of microtremor measurements to the estimation of earthquake ground motions in Kushiro city during the Kushiro-Oki earthquake of 15 january 1993

Microtremor measurements were conducted in Kushiro City, Hokkaido Island of Japan. The thickness of alluvial deposits in the central area of Kushiro City abruptly changes from 0 to 80 m. Predominant ground motion at a microtremor measuring site is represented by the spectral ratio, the ratio of Fourier amplitude spectrum of microtremor in the horizontal direction to that in the vertical direction. The peak frequency of the spectral ratio corresponds well to the predominant frequency estimated from the thickness of alluvial layer at each site. Based on a hypothesis that the spectral ratio can be regarded identical with a half of the amplification factor from diluvial bed rock to the ground surface, the strong ground motion due to the Kushiro-oki earthquake at each site is estimated. The method of estimation is verified by the comparison of the synthesized accelerogram with the recorded one at West Port of Kushiro. The results satisfactorily explain major damages and vibrations felt in Kushiro City during the earthquake.     

Site Response Evaluation Using the H/V Ratio at the Yan-Liau Station in Hualien, Taiwan

—?Four days after the December 13, 1990 Hualien, Taiwan earthquake (M _L = 6.9), a temporary array of fifteen triaxial digital accelerographs was deployed in the epicenter area to monitor aftershocks. Approximately 600 earthquakes triggered this array during the three-month deployment period. The Yan-Liau station (S63) alone recorded 162 events. Most of the accelerograms at S63 exhibit resonance. We have estimated site responses at the Yan-Liau station using both the single-station spectral ratio (or H/V ratio) method and the traditional spectral ratio method that compares ratios at a soil site with those at a reference hard rock site. Based on site response analyses of S waves and coda waves of ground motion recordings, both types of waves show that the H/V ratio provide a good estimate at the resonant frequency although the site amplification factor is overestimated. In addition, the study of microtremor is also a good alternate for estimating the site predominant frequency. While the ground acceleration (or PGA) gradually increases, the resonant frequency shifts to lower frequencies.     

Analysis of Short-Period Rayleigh Waves Recorded in the Bohemian Massif Area During Celebration 2000 Experiment

The aim of this paper is to show the application of short-period surface waves recorded during deep seismic sounding experiment for constraining shallow velocity structure of the crust. Phase velocity of fundamental mode Rayleigh waves, observed along the CELEBRATION 2000 experiment profile CEL09, were obtained by a p-ω method and has been subsequently inverted for one-dimensional shear velocity models for the top 2 km. Multiple filter technique applied to one shot gather was used to carry out a joint inversion of phase and group velocity data and to provide γ_R data to be used for Q_β inversion. Validity of obtained V_S and Q_β models was confirmed by the reflectivity method. Noticeably, no clear dispersive wawes were observed in the Tepla-Barrandian Unit. Quasi-2D model based on the individual 1D V_S models is well correlated with the surface geology. Lower V_S are observed in the Saxothuringian Zone in comparison to the Moldanubian Zone. In the vicinity of the Central Bohemian and Moldanubian Plutons, the near-surface V_S values are relatively low, but below 1 km depth, they are higher than in surrounding areas. We interpret it as the result of the weathering and cracks within the granitoid rocks.