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.     

A quantitative study of the energy release in the aftershocks of the Bhuj earthquake, 2001, India, using Lg phase

The devastating earthquake on 26 January 2001 at Bhuj, India, resulted in large-scale death and destruction of properties of several million US dollars. The moment magnitude of the earthquake was 7.7 and its maximum focal intensity exceeded X in MM scale. The rate of aftershocks of this earthquake, recorded at Gauribidanur seismic array station (GBA), shows a monotonic decay with time superposed with oscillations. For the Indian continent the Lg phase is a prominent arrival at regional distances. The estimate of Lg amplitude is obtained by optimally fitting the Lg wave train to a exponential decay curve. The logarithm of these amplitudes and logarithm of root mean square (rms) value of actual amplitudes of the Lg are calibrated with USGS m_b to create a local m_bLg magnitude scale. The energy released from these aftershocks is calculated from the rms value of Lg phase. The plot of cumulative energy release with time follows the power law of the form t^p, superposed with oscillations. The exponent of the power law, p, is estimated both by a time-window scanning method and by an interpolation method. The value of p is 0.434 for time-window scanning method and 0.432 for the interpolation method. The predominant periods found in the oscillatory part of the cumulative energy, obtained by differencing the observed from the power law fit, are 10.6, 7.9, 5.4, 4.6 and 3.5 h for time-window scanning method. The corresponding periods for interpolation method are 13.4, 11.5, 7.4, 4.2, 3.5, 2.6 and 2.4 h.     

Seismic wave velocity and anisotropy of serpentinized peridotite in the Oman ophiolite

Shallow seismic measurements in harzburgite from the Oman ophiolite performed in a zone where the maximum horizontal anisotropy is expected (vertical foliation and horizontal lineation) point to a dominant dependence of seismic properties on fracturing.

Optical microscopy studies show that microcracks are guided by the serpentine (lizardite) penetrative network oriented subparallel to the harzburgite foliation and subperpendicular to the mineral lineation, and that serpentine (lizardite) vein filling has a maximum concentration of (001) planes parallel to the veins walls. The calculated elastic properties of the oriented alteration veins filled with serpentine in an anisotropic matrix formed by oriented crystals of olivine and orthopyroxene are compared with seismic velocities measured on hand specimens.

Laboratory ultrasonic data indicate that open microcracks are closed at 100 MPa pressure, e.g. (J. Geophys. Res. 65, (1960) 1083) and (Proc. ODP Sci. Results Leg 118, (1990) 227). Above this pressure, laboratory measurements and modeling show that P-compressional and S-shear wave velocities are mainly controlled by the mineral preferred orientation. Veins sealed with serpentine are effective in slightly lowering P and S velocities and increasing anisotropy. The penetrative lizardite network does not affect directly the geometry of seismic anisotropy, but contributes indirectly in the fact that this network controls the microcrack orientations.

Comparison between seismic measurements of peridotite and gabbro in the same conditions suggest that P- and S-waves anisotropies are a possible discriminating factor between the two lithologies in the suboceanic lithosphere.     

Architecture of fault zones determined from outcrop, cores, 3-D seismic tomography and geostatistical modeling: example from the Albalá Granitic Pluton, SW Iberian Variscan Massif

The 3-D seismic tomographic data are used together with field, core and well log structural information to determine the detailed 3-D architecture of fault zones in a granitic massif of volume 500×575×168 m at Mina Ratones area in the Albalá Granitic Pluton. To facilitate the integration of the different data, geostatistical simulation algorithms are applied to interpolate the relatively sparse structural (hard) control data conditioned to abundant but indirect 3-D (soft) seismic tomographic data. To effectively integrate geologic and tomographic data, 3-D migration of the velocity model from the time domain into the depth domain was essential. The resulting 3-D model constitutes an image of the fault zone architecture within the granitic massif that honours hard and soft data and provides an evaluation of the spatial variability of structural heterogeneities based on the computation of 3-D experimental variograms of Fracture Index (fault intensity) data. This probabilistic quantitative 3-D model of spatially heterogeneous fault zones is suitable for subsequent fluid flow simulations. The modeled image of the 3-D fault distribution is consistent with the fault architecture in the Mina Ratones area, which basically consists of two families of subvertical structures with NNE–SSW and ENE–WSW trends that displaces the surfaces of low-angle faults (North Fault) and follows their seismically detected staircase geometry. These brittle structures cut two subvertical dykes (27 and 27′ Dykes) with a NNE–SSW to N–S trend. The faults present high FI (FI>12) adjacent bands of irregular geometry in detail that intersect in space delimiting rhombohedral blocks of relatively less fractured granite (FI<6). Both structural domains likely correspond with the protolith and the damaged zone/fault core in the widely accepted model for fault zone architecture. Therefore, the construction of 3-D grids of the FI in granitic areas affected by brittle tectonics permits the quantitative structural characterization of the rock massif.     

Geoid determination using adapted reference field, seismic Moho depths and variable density contrast

 The traditional remove-restore technique for geoid computation suffers from two main drawbacks. The first is the assumption of an isostatic hypothesis to compute the compensation masses. The second is the double consideration of the effect of the topographic–isostatic masses within the data window through removing the reference field and the terrain reduction process. To overcome the first disadvantage, the seismic Moho depths, representing, more or less, the actual compensating masses, have been used with variable density anomalies computed by employing the topographic–isostatic mass balance principle. In order to avoid the double consideration of the effect of the topographic–isostatic masses within the data window, the effect of these masses for the used fixed data window, in terms of potential coefficients, has been subtracted from the reference field, yielding an adapted reference field. This adapted reference field has been used for the remove–restore technique. The necessary harmonic analysis of the topographic–isostatic potential using seismic Moho depths with variable density anomalies is given. A wide comparison among geoids computed by the adapted reference field with both the Airy–Heiskanen isostatic model and seismic Moho depths with variable density anomaly and a geoid computed by the traditional remove–restore technique is made. The results show that using seismic Moho depths with variable density anomaly along with the adapted reference field gives the best relative geoid accuracy compared to the GPS/levelling geoid. Received: 3 October 2001 / Accepted: 20 September 2002 Correspondence to: H.A. Abd-Elmotaal     

地下水"恒流式"单井综合观测方法

以探索水位与水化学单井综合观测新方法为目的，建立了改性动水位概念，即恒定小流量状态下的井孔动水位。理论分析和实验表明，改性动水位与同一井孔对应的静水位之间相差一恒定水位落差。根据这一结果，可将观测得到的改性动水位值订正成相对应的静水位值。由此建立了以改性动水位和传统水化学为观测对象的“恒流式”单井综合观测方法，从而实现了静水位与传统水化学的单井综合观测。     

慢地震慢前兆的机制研究

通过岩石高温高压破裂实验与有机玻璃破裂试验的应力—应变曲线对比、地震波形记录(观测资料)与试样破裂波形(实验资料)对比、地震前后定点应变与主破裂前后应变形态变化的对比，认为慢地震是临界或亚临界破裂或预先滑动所致，是低频应变波动。它是材料屈服、弱化或塑性变化的必然结果。而慢前兆则是临近大地震(大破裂)前出现的诸如形变、低频地震波(破裂弹性波)等的短期及临震前兆现象，它比较可靠。但由于许多岩石主破裂前并无明显的临界破裂或预滑动现象，也并非所有岩石都一定出现明显屈服，所以也并非每个地震前都有慢地震，因此，也就不一定都出现低频波动。从而，也并不一定都出现相同的短期、临震前兆现象。     

震级加权四指标Blade算法及在地震带识别中的应用

Blade算法的产生是根据地震点的分布寻找地震断层，其优势在于能从较稀疏的点分布中找出线型条带模式，但不足的是完全不考虑地震的震级屑性。针对Blade算法的几个不足之处逐步改进，提出了震级加权四指标Blade算法。算法主要的改进体现在：①引入以震级为变量的权重函数，体现了高震级地震的控制作用，通过模拟与实例比较了三种权重函数的效果；②改进原方法中所挑出的Blade中心必须有地震点，引入四指标评价体系代替原方法的二指标评价体系，同时还根据指标值将挑选出的线型条带修正到更淮确的位置。将原算法与改进算法应用于华北地区4级以上的地震数据集以按寻地震带，结果表明，震级加权四指标Blade算法较Blade算法有显改进，取得了较好的应用效果。     

丽江7.0级地震重力前兆模式研究<

为了研究1996年丽江MＳ7.0地震前地震孕育过程或前兆表现，利用滇西地震实验场内高精度重复重力观测数据，结合地质调查和地球物理推断结果，并考虑到观测数据和模型粗差， 采用稳健或抗差-贝叶斯最小二乘算法和多断层位错模型，首次初步反演获得了研究区主要活动断裂滑动的时间变化分布．结果表明，1990～1997年断层运动的时间变化，较好地反映了1996年丽江MＳ7.0地震孕育过程．其主要前兆模式图象具有主震余震型特征，遵循地壳内部密度和地壳形变耦合运动模式（简称DD耦合运动模式）．      

利用地震层析成像数据计算地幔对流新模型的探讨

假设地幔地震层析成像数据对应的地幔横向不均匀结构是地幔热对流的结果. 将地震层析成像数据转化为地幔温度(或密度)不均匀分布，考虑热流体动力学的三个基本方程，顾及热输运方程中的非线性项，直接将地震层析成像转化的地幔温度不均匀分布作为内部荷载引入基本方程, 反演计算地幔对流. 本文在利用地震层析成像数据计算地幔对流模型的新理论和方法的基础上，用SH12WM13地震层析成像模型数据，计算了全球地幔对流格局. 结果表明，对流格局不仅依赖地震层析成像数据，而且在很大程度上受地幔动力学框架、热动力参数和边界条件的所确定的系统响应函数的影响. 显示了地幔中复杂的对流格局，特别是区域性层状对流以及多层对流环可能在地幔中存在的现象.