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.     

In situ measurements of sound velocities and densities across the orthopyroxene → high-pressure clinopyroxene transition in MgSiO3 at high pressure

Using acoustic measurement interfaced with a large volume multi-anvil apparatus in conjunction with in situ X-radiation techniques, we are able to measure the density and elastic wave velocities (V_P and V_S) for both ortho- and high-pressure clino-MgSiO₃ polymorphs in the same experimental run. The elastic bulk and shear moduli of the unquenchable high-pressure clinoenstatite phase were measured within its stability field for the first time. The measured density contrast associated with the phase transition OEN → HP-CEN is 2.6-2.9% in the pressure of 7-9 GPa, and the corresponding velocity jumps are 3-4% for P waves and 5-6% for S waves. The elastic moduli of the HP-CEN phase are K_S=156.7(8) GPa, G = 98.5(4) GPa and their pressure derivatives are K_S′=5.5(3) and G′ = 1.5(1) at a pressure of 6.5 GPa, room temperature. In addition, we observed anomalous elastic behavior in orthoenstatite at pressure above 9 GPa at room temperature. Both elastic wave velocities exhibited softening between 9 and 13-14 GPa, which we suggest is associated with a transition to a metastable phase intermediate between OEN and HP-CEN.     

Low-velocity zone and topography as a source of site amplification effect on Tarzana hill,California

Tarzana station is located in the foothills of the Santa Monica Mountains in California near the crest of a low (<20 m) natural hill with gentle slopes. The hill is about 500 m in length by 130 m in width and is formed of extremely weathered shale at the surface to fresh at depth. Average S-wave is about 250 m/s in the top 17–18 m, and S- and P-wave velocities significantly increase below this depth. According to the NEHRP classification based on V_S30∼300 m/s it is a site class D. Strong-motion instrumentation at Tarzana consisted of an accelerograph at the top of the hill, a downhole instrument at 60 m depth, and an accelerograph at the base of the hill. More than 20 earthquakes were recorded by at least three instruments at Tarzana from 1998 till 2003. Comparisons of recordings and Fourier spectra indicate strong directional resonance in a direction perpendicular to the strike of the hill. The dominant peaks in ground motion amplification on the top of the hill relative to the base are at frequencies ∼3.6 and 8–9 Hz for the horizontal components. Our hypothesis is that the hill acts like a wave trap. This results in an amplification at predominant frequencies f=V/4 h (h is layer's thickness) at f∼3.6 Hz for S-waves (using average V_S17=246 m/s and h=17 m) and f∼7.9 Hz for P-waves (using average V_P17=535 m/s and h=17 m). As was shown by Bouchon and Barker [Seismic response of a hill: the example of Tarzana, California. Bull Seism Soc Am 1996;86(1A):66–72], topography of this hill amplifies and polarizes ground motion in the frequency range of 3–5 Hz. Hill acts as a magnifying polarizing glass: It polarizes ground motion in the direction perpendicular to the strike of the hill and also amplifies ground motions that had been also amplified by a low-velocity layer.     

Near surface velocity and QS structure of the Quaternary sediment in Bohai basin,China

Heavily populated by Beijing and Tianjin cities, Bohai basin is a seismically active Cenozoic basin suffering from huge lost by devastating earthquakes, such as Tangshan earthquake. The attenuation (Q_P and Q_S) of the surficial Quaternary sediment has not been studied at natural seismic frequency (1?10 Hz), which is crucial to earthquake hazards study. Borehole seismic records of micro earthquake provide us a good way to study the velocity and attenuation of the surficial structure (0?500 m). We found that there are two pulses well separated with simple waveforms on borehole seismic records from the 2006 M_W4.9 Wen'an earthquake sequence. Then we performed waveform modeling with generalized ray theory (GRT) to confirm that the two pulses are direct wave and surface reflected wave, and found that the average v_P and v_S of the top 300 m in this region are about 1.8 km/s and 0.42 km/s, leading to high v_P/v_S ratio of 4.3. We also modeled surface reflected wave with propagating matrix method to constrain QS and the near surface velocity structure. Our modeling indicates that Q_S is at least 30, or probably up to 100, much larger than the typically assumed extremely low Q (~10), but consistent with Q_S modeling in Mississippi embayment. Also, the velocity gradient just beneath the free surface (0?50 m) is very large and velocity increases gradually at larger depth. Our modeling demonstrates the importance of borehole seismic records in resolving shallow velocity and attenuation structure, and hence may help in earthquake hazard simulation.     

Seismic Characterization of Neapolitan Soils

— Detailed shear-wave velocity profiles versus depth have been obtained in typical lithostratigraphies of Napoli. FTAN and hedgehog methods have been applied to Rayleigh surface waves recorded in refraction seismic surveys. The comparison with literature measurements shows good agreement with nearby down- and cross-hole tests. The pumiceous and lapilli content, and the different welding and alteration degree of the Neapolitan pyroclastic soils cause a strong scattering of the shear wave velocities (V_S) from bore-hole measurements, even for the same formation. Surface measurements, based on FTAN-hedgehog methods, determine average V_S along travel paths of about 100 m, give results that are comparable with down- and cross-hole velocity profiles, and have the additional advantage of being less scattered, and thus more representative of average properties than bore-hole measurements. The results of surface measurements should be preferred in the computation of realistic seismograms and are particularly suitable in urban areas, as they are not destructive and need just one receiver.     

Crustal thickness,VP/VS ratio,and shear wave velocity structures beneath Myanmar and their tectonic implications

The crustal structure in Myanmar can provide valuable information for the eastern margin of the ongoing Indo-Eurasian collision system. We successively performed H–k stacking of the receiver function and joint inversion of the receiver function and surface wave dispersion to invert the crustal thickness (H), shear wave velocity (V_S), and the V_P/V_S ratio (k) beneath nine permanent seismic stations in Myanmar. H was found to increase from 26 ?km in the south and east of the study area to 51 ?km in the north and west, and the V_P/V_S ratio was complex and high. Striking differences in the crust were observed for different tectonic areas. In the Indo-Burma Range, the thick crust (H ?～ ?51 ?km) and lower velocities may be related to the accretionary wedge from the Indian Plate. In the Central Myanmar Basin, the thin crust (H ?= ?26.9–35.5 ?km) and complex V_P/V_S ratio and V_S suggest extensional tectonics. In the Eastern Shan Plateau, the relatively thick crust and normal V_P/V_S ratio are consistent with its location along the western edge of the rigid Sunda Block.     

Shear Wave Velocity Estimates through Combined Use of Passive Techniques in a Tectonically Active Area

We made an attempt to assess the shear wave velocity values V _S and, to a lesser extent, the V _P values from ambient noise recordings in an array configuration. Five array sites were situated in the close proximity to borehole sites. Shear wave velocity profiles were modeled at these five array sites with the aid of two computational techniques, viz. spatial autocorrelation (SPAC) and H/V ellipticity. Out of these five array sites, velocity estimates could be reliably inferred at three locations. The shear wave velocities estimated by these methods are found to be quite consistent with each other. The computed V _S values up to 30 m depth are in the range from 275 to 375 m/s in most of the sites, which implies prevalence of a low velocity zone at some pocket areas. The results were corroborated by evidence of site geology as well as geotechnical information.     

Site Specific Ground Motion Modeling and Seismic Response Analysis for Microzonation of Baku,Azerbaijan

We investigated ground response for Baku (Azerbaijan) from two earthquakes of magnitude M6.3 occurred in Caspian Sea (characterized as a near event) and M7.5 in Shamakhi (characterized as a remote extreme event). S-wave velocity with the average shear wave velocity over the topmost 30 m of soil is obtained by experimental method from the V _P values measured for the soils. The downtown part of Baku city is characterized by low V_S30 values (< 250 m/s), related to sand, water-saturated sand, gravel-pebble, and limestone with clay. High surface PGA of 240 gal for the M7.5 event and of about 190 gal for the M6.3 event, and hence a high ground motion amplification, is observed in the shoreline area, through downtown, in the north-west, and in the east parts of Baku city with soft clays, loamy sands, gravel, sediments.     

Gravelly soils that liquefied during 2008 Wenchuan,China earthquake, Ms=8.0

Field investigations following the 2008 Wenchuan earthquake (M_s=8.0) identified 118 liquefaction sites nearly all of which are underlain by gravelly sediment in the Chengdu Plain and adjacent Mianyang area. Field studies, including core drilling, dynamic penetration tests (DPT), and multiple channel analysis of surface wave velocity tests (MASW) for measurement of shear wave velocities, reveal the following: (1) Sand boils and ground fissures, indicative of liquefaction, occurred across hundreds of square kilometers affecting 120 villages, 8 schools and 5 factories. (2) The Chengdu plain is underlain by sandy gravels ranging in thickness up to 540 m; loose upper layers within the gravels beds liquefied. (3) Mean grain sizes for gravelly layers that liquefied range from 1 mm to more than 30 mm. (4) Shear wave velocities in gravels that liquefied range up to 250 m/s. (5) A 50% probability curve, developed from logistic procedures, correctly bounds all but four data points for the 47 compiled V_s data.     

The effect of wave-induced turbulence on intertidal mudflats: Impact of boat traffic and wind

Semi-diurnal and fortnightly surveys were carried out to quantify the effects of wind- and navigation-induced high-energy events on bed sediments above intertidal mudflats. The mudflats are located in the upper fluvial part (Oissel mudflat) and at the mouth (Vasière Nord mudflat) of the macrotidal Seine estuary. Instantaneous flow velocities and mudflat bed elevation were measured at a high frequency and high resolution with an acoustic doppler velocimeter (ADV) and an ALTUS altimeter, respectively. Suspended particulate matter concentrations were estimated by calibrating the ADV acoustic backscattered intensity with bed sediments collected at the study sites. Turbulent bed shear stress values were estimated by the turbulent kinetic energy method, using velocity variances filtered from the wave contribution. Wave shear stress and maximum wave–current shear stress values were calculated with the wave–current interaction (WCI) model, which is based on the bed roughness length, wave orbital velocities and the wave period (T_S). In the fluvial part of the estuary, boat passages occurred unevenly during the surveys and were characterized by long waves (T_S>50 s) induced by the drawdown effect and by short boat-waves (T_S<10 s). Boat waves generated large bottom shear stress values of 0.5 N m⁻² for 2–5 min periods and, in burst of several seconds, larger bottom shear stress values up to 1 N m⁻². At the mouth of the estuary, west south-west wind events generated short waves (T_S<10 s) of H_S values ranging from 0.1 to 0.3 m. In shallow-water environment (water depth <1.5 m), these waves produced bottom shear stress values between 1 and 2 N m⁻². Wave–current shear stress values are one order of magnitude larger than the current-induced shear stress and indicate that navigation and wind are the dominant hydrodynamic forcing parameters above the two mudflats. Bed elevation and SPM concentration time series showed that these high energy events induced erosion processes of up to several centimetres. Critical erosion shear stress (τ_ce) values were determined from the SPM concentration and bed elevation measurements. Rough τ_ce values were found above 0.2 N m⁻² for the Oissel mudflat and about 1 N m⁻² for the Vasière Nord mudflat.     

Soil characterization of Tınaztepe region (İzmir/Turkey) using surface wave methods and nakamura (HVSR) technique

To determine the shear wave velocity structure and predominant period features of T?naztepe in ?zmir, Turkey, where new building sites have been planned, active–passive surface wave methods and single-station microtremor measurements are used, as well as surface acquisition techniques, including the multichannel analysis of surface waves (MASW), refraction microtremor (ReMi), and the spatial autocorrelation method (SPAC), to pinpoint shallow and deep shear wave velocity. For engineering bedrock (V _s > 760 m/s) conditions at a depth of 30 m, an average seismic shear wave velocity in the upper 30 m of soil (AV_s30) is not only accepted as an important parameter for defining ground behavior during earthquakes, but a primary parameter in the geotechnical analysis for areas to be classified by V _s30 according to the National Earthquake Hazards Reduction Program (NEHRP). It is also determined that Z1.0, which represents a depth to V _s = 1000 m/s, is used for ground motion prediction and changed from 0 to 54 m. The sediment–engineering bedrock structure for T?naztepe that was obtained shows engineering bedrock no deeper than 30 m. When compared, the depth of engineering bedrock and dominant period map and geology are generally compatible.     

四川、甘肃地区VS30经验估计研究

目前我国建筑工程抗震设计规范中对于工程场地条件的判断依据主要是地表以下20m深度范围内土层的等效剪切波速,简称VS20。相比之下,国外应用较广的是地表以下30m深度范围内的等效剪切波速,简称VS30。这种差别导致国内科研工作者在应用国外的地震工程、工程抗震模型时经常遇到对场地条件描述不准确的困难。为了解决这个问题,本文根据147个四川、甘肃地区国家强震动台站20m左右深度的钻孔剪切波速数据,利用延拓方法、场地分类统计方法以及基于地形特征的VS30估计方法研究各台站VS30与VS20的经验关系,对比发现基于速度梯度延拓的结果最为可取。参考国际上通用的Geomatrix Classification场地分类标准,最终得到四川、甘肃地区各类场地的平均VS30,此结果可以为缺乏钻孔数据的工程场地的VS30估计提供参考。     

Fine S wave velocity structure beneath Iwate volcano, northeastern Japan, as derived from receiver functions and travel times

A genetic algorithm inversion of receiver functions derived from a dense seismic network around Iwate volcano, northeastern Japan, provides the fine S wave velocity structure of the crust and uppermost mantle. Since receiver functions are insensitive to an absolute velocity, travel times of P and S waves propagating vertically from earthquakes in the subducting slab beneath the volcano are involved in the inversion. The distribution of velocity perturbations in relation to the hypocenters of the low-frequency (LF) earthquakes helps our understanding of deep magmatism beneath Iwate volcano. A high-velocity region (dV_S/V_S=10%) exists around the volcano at depths of 2–15 km, with the bottom depth decreasing to 11 km beneath the volcano’s summit. Just beneath the thinning high-velocity region, a low-velocity region (dV_S/V_S=−10%) exists at depths of 11–20 km. Intermediate-depth LF (ILF) events are distributed vertically in the high-velocity region down to the top of the low-velocity region. This distribution suggests that a magma reservoir situated in the low-velocity region supplies magma to a narrow conduit that is detectable by the hypocenters of LF earthquakes. Another broad low-velocity region (dV_S/V_S=−5 to −10%) occurs at depths of 17–35 km. Additional clusters of deep LF (DLF) events exist at depths of 32–37 km in the broad low-velocity zone. The DLF and ILF events are the manifestations of magma movement near the Moho discontinuity and in the conduit just beneath the volcano, respectively.     

Stepwise joint inversion of surface wave dispersion,Rayleigh wave ZH ratio,and receiver function data for 1D crustal shear wave velocity structure

Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio (i.e., ellipticity), and receiver function data to better resolve 1D crustal shear wave velocity (v _S) structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensitivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast (e.g., due to the existence of crustal interfaces) and v _P/v _S ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the v _S model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute v _S model and then incorporate receiver function data in the joint inversion to obtain a finer v _S model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal v _S structures and with little initial model dependency.     

Shear and dilatational wave velocities for unsaturated soils

The primary objective of this study is for presenting some simple-to-use expressions relating the shear and dilatational wave velocities (V_S and V_P) to some physical and constitutive parameters of unsaturated soils. To this purpose, a simplified formulation is developed using the theory of linear poroelasticity in conjunction with some constitutive parameters widely used in geotechnical engineering. The derived expressions are of practical interest in view of the fact that they could be employed for evaluating the involved soil parameters from V_S and V_P measurements by in-situ or laboratory geophysical tests.     

Integrative assessment of coastal pollution: Development and evaluation of sediment quality criteria from chemical contamination and ecotoxicological data

Elutriate embryo-larval bioassays with sea-urchins (Paracentrotus lividus) were conducted concurrently with chemical analyses of sediments and biota as part of an integrative assessment of pollution in highly productive coastal regions. High metal contents and organic compounds in sediments and mussels were found in localised areas from the inner part of the estuaries indicating a clear anthropogenic influence. In particular, average maximum concentrations of 2803 mg Cu/kg dw, 776 mg Pb/kg dw, 2.5 mg Hg/kg dw and 5803 μg ∑₇PAHs/kg dw were measured in sediments from the most polluted sites. Significant correlations were observed between sediment chemistry and toxicity bioassays. Moreover, the Mantel test revealed a significant correlation (r_M=0.80; p<0.01) between sediment pollutant concentrations and toxicity data profiles. In addition, sediment quality criteria were used to help in the ecological interpretation of sediment chemistry data and to identify pollutants of concern. The toxicity bioassays identified polluted sites and quantified the level of toxicity, providing a cost-effective tool to complement the routine chemical monitoring currently conducted in European coastal waters with ecologically relevant information. This is in line with the recent European legislation that advocates the use of biological tools with the ultimate aim of protecting marine resources from anthropogenic substances that will affect their sensitive early life stages.     

Evaluation of shear wave velocity profile using SPT based uphole method

The uphole method is a field seismic test which uses receivers on the ground surface and an underground source. A modified form of the uphole method is introduced in order to obtain efficiently the shear wave velocity (V_S) profile of a site. This method is called the standard penetration test (SPT)-uphole method because it uses the impact energy of the split spoon sampler in the SPT test as a source. Since the SPT-uphole method can be performed simultaneously with the SPT test it is economical and not labor intensive compared to the original uphole methods which use small explosives or a mechanical source. Field testing and interpretation procedures for the proposed method are described. To obtain reliable travel time information of the shear wave, the first peak point of the shear wave using two component geophones is recommended. Through a numerical study using the finite element method (FEM), the procedure of the proposed method was verified. Finally, the SPT-uphole method was performed at several sites, and the field applicability of the proposed method was verified by comparing the V_S profiles determined by the SPT-uphole method with the profiles determined by the downhole, the spectral analysis of surfaces waves (SASW) method and from the SPT-N values.     

Evolution of droplets in subsea oil and gas blowouts: Development and validation of the numerical model VDROP-J

The droplet size distribution of dispersed phase (oil and/or gas) in submerged buoyant jets was addressed in this work using a numerical model, VDROP-J. A brief literature review on jets and plumes allows the development of average equations for the change of jet velocity, dilution, and mixing energy as function of distance from the orifice. The model VDROP-J was then calibrated to jets emanating from orifices ranging in diameter, D, from 0.5 mm to 0.12 m, and in cross-section average jet velocity at the orifice ranging from 1.5 m/s to 27 m/s. The d₅₀/D obtained from the model (where d₅₀ is the volume median diameter of droplets) correlated very well with data, with an R² = 0.99. Finally, the VDROP-J model was used to predict the droplet size distribution from Deepwater Horizon blowouts. The droplet size distribution from the blowout is of great importance to the fate and transport of the spilled oil in marine environment.     

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.     

On the statistical distribution of seismic velocities in Earth's deep mantle

The existence of uncoupled shear (S) and compression (P) wave velocity variations in Earth's mantle is a characteristic that might only be explained by the presence of significant chemical and/or phase heterogeneity, with important implications for the dynamics and evolution of Earth's interior. While making a one-to-one comparison between tomographic models for P and S velocity (V_P and V_S) variations for a particular geographic region is ill-posed, their global statistical distributions reveal several robust characteristics indicative of the nature of uncoupled V_P and V_S in the deep mantle. We find that all of the V_P and V_S model distributions at a given depth are Gaussian-like throughout the lowermost mantle. However, a distinct low velocity feature is present in V_S distributions below ≈ 2200 km depth that is not present or is relatively weak in V_P models. The presence of anomalously low V_S material cannot be explained as an artifact, nor can the absence of a similarly strong feature in P models be ascribed to under-resolution. We propose that this feature can be partly explained by laterally variable occurrences of post-perovskite (pPv) lenses in the D″ layer, however, the persistence of significantly slow V_S regions at heights up to ≈ 700 km or more above the core–mantle boundary is likely to be incompatible with a pPv origin and might only be explained by the presence of a laterally discontinuous layer of chemically distinct material and/or some other kind of phase heterogeneity. There also exist significant discrepancies between tomographic models with respect to the width of the distributions as well as differences between the modeled peak values. We propose a scheme for comparison between different seismic models in which the widths of the dominant features in their statistical distributions is exploited.