P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation

P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.     

Estimation of in situ viscoelastic parameters of weak floor strata by plate-loading tests

Summary This paper explores an approach to estimation of the viscoelastic parameters of weak floor strata using plate-loading tests. Mathematical equations are derived to describe deformation-time history for the standard Burgers model under three-dimensional stress-strain conditions, which are verified with a finite element model. A number of time-dependent in-mine plate-loading tests were conducted and the values of the viscoelastic parameters were estimated based on the derived equation. It is shown that large discrepancies between the parameter values estimated from the plate-loading tests and those acquired from the in-mine convergence data are the results of neglecting adjacent pillar interaction. Finite element modelling was conducted to investigate the effects of adjacent pillar interactions upon the pillar settlement. Based on the finite element analyses, a correction curve was developed to adjust the viscoelastic parameter values to take into account the interaction between adjacent pillars. With the correction, the parameter values estimated from the plate-loading tests compared favourably with those acquired from the in-mine convergence results.     

Modelling of soil temperature profiles using data from the pilot experiment at Anand

We present some preliminary results from a multi-layer soil temperature finite-difference model using the data set obtained during the pilot experiment on land surface processes at Anand. We hope that the present results will prove useful during the final field experiments scheduled for the coming monsoon season.     

土的边界面本构模型研究进展

简述了土的边界面本构模型国内外研究状况，阐述了其基本思想。对边界面模型作了简要评述，提出了边界面模型研究的发展趋向。     

Recent Advances in Prediction and Processing of Strong Ground Motions

We present an overview of our recent results on utilizing small earthquakes in the earthquake engineering practice. Site-specific ground motion time-histories of large earthquakes can be successfully simulated using recordings of small earthquakes which are often referred to as 'empirical Green's functions' in seismology. Another important practical problem is whether and how these observations can be used in seismic risk studies which are based on empirical attenuation relations for ground motion parameters. We study a possibility of extrapolating attenuation relations for small earthquakes, to larger magnitudes using the data from the Garner Valley downhole array in Southern California. Finally we introduce efficient ground motion processing techniques in frequency- and time-domains and apply them to site response estimation.     

Seismic Ground Motion Expected for the Eastern District of Naples

The seismic ground motion of a test area in the eastern district of Naples is computed with a hybrid technique based on the mode summation and the finite difference methods. This technique allows us the realistic modelling of source and propagation effects, including local soil conditions. In the modelling, we consider the 1980 Irpinia earthquake, a good example of strong shaking for the area of Naples, which is located about 90 km from the epicenter.The detailed geological setting is reconstructed from a large number of drillings. The sub-soil is mainly formed by alluvial (ash, stratified sand and peat) and pyroclastic materials overlying a pyroclastic rock (yellow neapolitan tuff), representing the neapolitan bedrock. The detailed information available on mechanical properties of the sub-soil and its geometry warrants the application of the sophisticated hybrid technique.As expected, the sedimentary cover causes an increase of the signal's amplitudes and duration. If thin peat layers are present, the amplification effects are reduced, and the peak ground accelerations are similar to those observed for the bedrock model. This can be explained by the backscattering of wave energy at such layers, that tend to seismically decouple the upper from the lower part of the structure.For SH-waves, the influence of the variations of the S-wave velocities on the spectral amplification is studied, by considering locally measured velocities and values determined from near-by down-hole measurements. The comparison between the computed spectral amplifications confirms the key role of an accurate determination of the seismic velocities of the different layers.The comparison performed between a realistic 2-D seismic response and a standard 1-D response, based on the vertical propagation of waves in a plane layered structure, shows considerable difference, from which it is evident that serious caution must be taken in the modelling of expected ground motion at a specific site.     

龙首二级面板堆石坝三维真非线性地震反应分析和评价

基于土石料三维粘弹塑性动力本构模型,并采用新型三维各向异性有厚度薄单元来模拟面板和堆石的接触面特性,建立了高面板堆石坝地震反应分析的三维真非线性动力分析方法.利用坝料动力特性的大型三轴试验成果,分析计算了龙首二级(西流水)面板堆石坝的地震反应,主要包括加速度反应、堆石体应力反应及坝体单元抗震安全系数、面板应力反应和变形及接缝位移、高趾墙动力反应等.为大坝的抗震设计提供了有力的技术依据.