Modeling of Rayleigh waves dispersion in the Sannio region (Southern Italy) from seismic active refraction data

Short period surface waves, recorded during a seismic refractionsurvey in the Sannio region (Southern Italy), have been modeled to infera shallow velocity model for the area. Based on the decrease of resolutionwith depth, due to the bias on group velocity estimates arising frominterference of the Rayleigh waves with higher modes, we carried out aprocedure of fitting, with synthetic seismograms, of selected filtered traceswith a gaussian filter, having a width at half height equal to 1 Hz and acentral frequency lying in the range [1,4] Hz. We estimated the likelihoodbetween synthetic and observed seismograms by measuring their semblance.In this way we were able to infer a more refined local velocity modelcharacterized by a high Vp and Vs vertical gradient in the sedimentarycover. Two ad hoc resolution studies, based on group velocity andamplitude data respectively, indicate that the local velocity model is a goodvelocity model also for the entire studied area. The increase in the numberof available data when using amplitude information allows us to make amore selective choice in the model parameter space (Vp and Vs of eachlayer) and to solve for the Vp/Vs ratio. The inferred Vp velocity in thehalf-space is equal to 2.8 km/s. This value is in excellent agreement withthat inferred by other authors (3 km/s) by modeling P-wave travel timevs. distance. The best-fit model furnish low Vp/Vs for the sedimentarycover so indicating a high degree of the sediment's compaction in thestudied area. The inferred shallow high-velocity gradient indicates thatthe shallow sedimentary layer in the area could trap and focus the energytraveling into it.     

Practical VTI approximations: a systematic anatomy

Transverse isotropy (TI) with a vertical symmetry axis (VTI) often provides an appropriate earth model for prestack imaging of steep-dip reflection seismic data. Exact P-wave and SV-wave phase velocities in VTI media are described by complicated equations requiring four independent parameters. Estimating appropriate multiparameter earth models can be difficult and time-consuming, so it is often useful to replace the exact VTI equations with simpler approximations requiring fewer parameters. The accuracy limits of different previously published VTI approximations are not always clear, nor is it always obvious how these different approximations relate to each other. Here I present a systematic framework for deriving a variety of useful VTI approximations. I develop first a sequence of well-defined approximations to the exact P-wave and SV-wave phase velocities. In doing so, I show how the useful but physically questionable heuristic of setting shear velocities identically to zero can be replaced with a more precise and quantifiable approximation. The key here to deriving accurate approximations is to replace the stiffness a₁₃ with an appropriate factorization in terms of velocity parameters. Two different specific parameter choices lead to the P-wave approximations of Alkhalifah (Geophysics 63 (1998) 623) and Schoenberg and de Hoop (Geophysics 65 (2000) 919), but there are actually an infinite number of reasonable parametrizations possible. Further approximations then lead to a variety of other useful phase velocity expressions, including those of Thomsen (Geophysics 51 (1986) 1954), Dellinger et al. (Journal of Seismic Exploration 2 (1993) 23), Harlan (Stanford Exploration Project Report 89 (1995) 145), and Stopin (Stopin, A., 2001. Comparison of v(θ) equations in TI medium. 9th International Workshop on Seismic Anisotropy). Each P-wave phase velocity approximation derived this way can be paired naturally with a corresponding SV-wave approximation. Each P-wave or SV-wave phase velocity approximation can then be converted into an equivalent dispersion relation in terms of horizontal and vertical slownesses. A simple heuristic substitution also allows each phase velocity approximation to be converted into an explicit group velocity approximation. From these, in turn, travel time or moveout approximations can also be derived. The group velocity and travel time approximations derived this way include ones previously used by Byun et al. (Geophysics 54 (1989) 1564), Dellinger et al. (Journal of Seismic Exploration 2 (1993) 23), Tsvankin and Thomsen (Geophysics 59 (1994) 1290), Harlan (89 (1995) 145), and Zhang and Uren (Zhang, F. and Uren, N., 2001. Approximate explicit ray velocity functions and travel times for P-waves in TI media. 71st Annual International Meeting, Society of Exploration Geophysicists, Expanded Abstracts, 106–109).     

Nonlinear Rate Dependent Model of High Damping Rubber Bearing

The accurate analysis of the response of isolated structures requires the application of appropriate models of isolation devices. The purpose of this paper is to analyse a nonlinear strain rate dependent model of a high damping rubber bearing which simulates the horizontal behaviour of the device under specified vertical load using a nonlinear elastic spring-dashpot element. The effectiveness of the model is checked by fitting the experimental data concerning three different rubber bearings. The results of the study show that the model can simulate the bearing behaviour over a wide shear strain range with small simulation errors. This revised version was published online in July 2006 with corrections to the Cover Date.     

Seismic evidence for deep low-velocity anomalies in the transition zone beneath West Antarctica

We present a three-dimensional (3D) SV-wave velocity model of the upper mantle beneath the Antarctic plate constrained by fundamental and higher mode Rayleigh waves recorded at regional distances. The good agreement between our results and previous surface wave studies in the uppermost 200 km of the mantle confirms that despite strong differences in data processing, modern surface wave tomographic techniques allow to produce consistent velocity models, even at regional scale. At greater depths the higher mode information present in our data set allows us to improve the resolution compared to previous regional surface wave studies in Antarctica that were all restricted to the analysis of the fundamental mode. This paper is therefore mostly devoted to the discussion of the deeper part of the model. Our seismic model displays broad domains of anomalously low seismic velocities in the asthenosphere. Moreover, we show that some of these broad, low-velocity regions can be more deeply rooted. The most remarkable new features of our model are vertical low-velocity structures extending from the asthenosphere down to the transition zone beneath the volcanic region of Marie Byrd Land, West Antarctica and a portion of the Pacific-Antarctic Ridge close to the Balleny Islands hotspot. A deep low-velocity anomaly may also exist beneath the Ross Sea hotspot. These vertical structures cannot be explained by vertical smearing of shallow seismic anomalies and synthetic tests show that they are compatible with a structure narrower than 200 km which would have been horizontally smoothed by the tomographic inversion. These deep low-velocity anomalies may favor the existence of several distinct mantle plumes, instead of a large single one, as the origin of volcanism in and around West Antarctica. These hypothetical deep plumes could feed large regions of low seismic velocities in the asthenosphere.     

A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

Attenuation of body waves in Southeastern Sicily (Italy)

The attenuation of P- and S-waves in Southeastern Sicily was estimated by applying two different methods in time and frequency domains. We analyzed waveforms from about 290 local events (0.6≤M_L≤4.6) recorded at a three-component digital network.By applying the pulse broadening method to the first P-wave pulse, we found an average Q_p value of ca. 140. The application of the frequency decay method provided a Q_p value of ca. 120, in the low-frequency band (3-9 Hz). Conversely, in the high frequency range (16-27 Hz) the average Q_p is significantly larger (ca. 640). The frequency decay method was also applied to S-waves spectra. In the low frequency range (2-5 Hz) the estimated average Q_s is ca. 190. As for Q_p, also Q_s, in the high frequency range (16-27 Hz), is larger (ca. 700). These results evidenced a frequency dependence of both the quality factors Q_p and Q_s, as commonly observed in tectonically active zones characterized by high degree of heterogeneity.     

平面散点集Delaunay三角剖分的一种高效方法

以平面散点集逐点插入的Delaunay三角化的方法为基础，在三角化过程中采用一定策略，将其改进成为一种简单易行而高效的方法，能够适应包括多岛、多连通域等复杂情况的各种边界，能够生成贴体的三角网，网格能够保证符合Delaunay法则。     

Estimating palaeowind strength from beach deposits

Abstract The geological record of past wind conditions is well expressed in the coarse gravel, cobble and boulder beach deposits of Quaternary palaeolakes in the Great Basin of the western USA and elsewhere. This paper describes a technique, using the particle‐size distribution of beach deposits, to reconstruct palaeowind conditions when the lakes were present. The beach particle technique (BPT) is first developed using coarse beach deposits from the 1986–87 highstand of the Great Salt Lake in Utah, combined with instrumental wind records from the same time period. Next, the BPT is used to test the hypothesis that wind conditions were more severe than at present during the last highstand of Lake Lahontan (≈ 13 ka), which only lasted a decade or two at most. The largest 50 beach clasts were measured at nine beach sites located along the north, west and south sides of Antelope Island in the Great Salt Lake, all of which formed in 1986–87. At these sites, the largest clast sizes range from 10 to 28 cm (b‐axis), and fetch lengths range from 25 to 55 km. Nearshore wave height was calculated by assuming that the critical threshold velocity required to move the largest clasts represents a minimum estimate of the breaking wave velocity, which is controlled by wave height. Shoaling transformations are undertaken to estimate deep‐water wave heights and, ultimately, wind velocity. Wind estimates for the nine sites, using the BPT, range from 6·5 to 17·4 m s^?1, which is in reasonable agreement with the instrumental record from Salt Lake City Airport. The same technique was applied to eight late Pleistocene beaches surrounding the Carson Sink sub‐basin of Lake Lahontan, Nevada. Using the BPT, estimated winds for the eight sites range from 9·7 to 27·1 m s^?1. The strongest winds were calculated for a cobble/boulder beach with a fetch of 25 km. Instrumental wind records for the 1992–99 period indicate that wind events of 9–12 m s^?1 are common and that the strongest significant wind event (≥ 9 m s^?1 for ≥ 3 h) reached an average velocity of 15·5 m s^?1. Based on this preliminary comparison, it appears that the late Pleistocene western Great Basin was a windier place than at present, at least for a brief time.     

Permeability Estimation with the Augmented Lagrangian Method for a Nonlinear Diffusion Equation

We consider numerical identification of the piecewise constant permeability function in a nonlinear parabolic equation, with the augmented Lagrangian method. By studying this problem, we aim at also gaining some insight into the potential ability of the augmented Lagrangian method to handle permeability estimation within the full two-phase porous-media flow setting. The identification is formulated as a constrained minimization problem. The parameter estimation problem is reduced to a coupled nonlinear algebraic system, which can be solved efficiently with the conjugate gradient method. The methodology is developed and numerical experiments with the proposed method are presented.     

偏心结构非线性地震反应分析的一种简化方法

本文建立了基于结构状态方程和偏心结构体系的层剪力-扭矩等效屈服EYST面概念的求解平扭耦联非线性地震反应的分离-予估-修正递推算法，编制了相应的计算程序。通过与三个单层偏心结构模型的破坏性地震模拟试验结果的对比，验证了小变形阶段平-扭耦联地震反应算法的实用性和可靠性。算例表明，该程序数据处理量小，运算快捷，便于应用。