Wide-angle seismic velocities in heterogeneous crust

Seismic velocities measured by wide-angle surveys are commonly used to constrain material composition in the deep crust. Therefore, it is important to understand how these velocities are affected by the presence of multiscale heterogeneities. The effects may be characterised by the scale of the heterogeneity relative to the dominant seismic wavelength (λ); what is clear is that heterogeneities of all scales and strengths bias wide-angle velocities to some degree. Waveform modelling was used to investigate the apparent wide-angle P -wave velocities of different heterogeneous lower crusts. A constant composition (50 per cent felsic and 50 per cent ultramafic) was formed into a variety of 1- and 2-D heterogeneous arrangements and the resulting wide-angle seismic velocity was estimated. Elastic, 1-D models produced the largest velocity shift relative to the true average velocity of the medium (which is the velocity of an isotropic mixture of the two components). Thick (width > λ) horizontal layers, as a result of Fermat's Principle, provided the largest increase in velocity; thin (width ≪λ) vertical layers produced the largest decrease in velocity. Acoustic 2-D algorithms were shown to be inadequate for modelling the kinematics of waves in bodies with multiscale heterogeneities. Elastic, 2-D modelling found velocity shifts (both positive and negative) that were of a smaller magnitude than those produced by 1-D models. The key to the magnitude of the velocity shift appears to be the connectivity of the fast (and/or slow) components. Thus, the models with the highest apparent levels of connectivity between the fast phases, the 1-D layers, produced the highest-magnitude velocity shifts. To understand the relationship between measured seismic velocities and petrology in the deep crust it is clear that high-resolution structural information (which describes such connectivity) must be included in any modelling.     

Fully normalized spherical cap harmonics: application to the analysis of sea-level data from TOPEX/POSEIDON and ERS-1

We introduce two sets of fully normalized harmonics for the spectral analysis of functions defined on a spherical cap. The harmonics are the products of Fourier functions and the fully normalized associated Legendre functions of non-integer degree. Using Sturm-Liouville theory for boundary-value problems, we present two convenient and stable formulae for computing the zeros of the associated Legendre functions that form two sets of orthogonal functions. Formulae for the stable numerical evaluation of the fully normalized associated Legendre functions of non-integer degree that avoid the gamma function are also derived. The result from the expansions of sea-level anomaly from altimetry into Set 2 fully normalized cap harmonics shows fast convergence of the series, and the degree variances decay rapidly without aliasing effects. The zero-degree coefficients (Set 2) of sea-level anomaly from TOPEX/ POSEIDON (T/P) and ERS-1 indicate an El Niño event during 1993 January-1993 July, and a La Niño event during 1993 November-1994 July, although the ERS-1 result is less obvious. Ocean circulations over the South China Sea and the Kuroshio area are clearly identified with the low-degree expansions of sea-surface topography (SST) from T/P and ERS-1. A cold-core eddy of 4° in diameter centred at 17.5°N, 118 E was detected with the expansion of SST from T/P cycle 47, and a property of the cap harmonics is used to compute this eddy's kinetic energy. The kinetic energy is at a low in winter and high in summer, and its variation seems to be periodic with an amplitude of 0.4 m² S^-2.     

Phase-array decomposition of a seismic section

A seismic re fraction/wide-angle reflection profile is analysed for the presence of correlated events ('phases'). The correlation problem is formulated in terms of temporally, spatially and frequency-local complex covariances. For robustness, the method concentrates on phase rather than amplitude information. This allows a computationally efficient algorithm that can make allowance for signal correlation length and can model curved wavefronts. A statistical test based on residual phase misfit across the analysed subarray is used to assess the probability that a detected event represents a real correlated signal.
With our chosen analysis parameters and confidence level (over 99.9 per cent). 1222 events were detected in the data. Using simple techniques based on 1-D earth models, detected events are associated with a small number of particular wave types. In this way, we have succeeded in classifying almost 95 per cent of the detected events. Those that remain describe those components of the data that are inconsistent with our simple ray paths in the 1-D assumption and with our prescribed tolerance. These include reverberations, near-surface guided waves and reflected waves from strongly laterally inhomogeneous structures. According to our modelling, about 25 per cent of the detected events are consistent with simple P -wave reflected energy, and these are to a very large extent (over 85 per cent) distinct from all the other wave-type models we have used. A direct mapping of the detected events into the offset-depth domain reveals dear internal and external consistencies among the detections for the various wave types. Estimated earth structure is consistent with models from previous analyses based on much larger data sets.
We have thus succeeded in extracting correlated events from the data and decomposing these, approximately but meaningfully, into distinct classes (ray paths)     

利用不纯净十水芒硝与纯净氯化钾制取硫酸钾的研究

本文通过相图分析，对利用不纯净十水芒硝与纯净氯化钾通过二步转化制取硫酸钾进行了物料配比计算．给出了最佳物料配比及十水芒硝中氯离子的最大含量，对实际生产有指导意义。     

三峡工程永久船闸边坡岩体卸荷特性

从地貌学角度,分析了三峡工程永久船闸边坡岩体的自然及人为卸荷特性,指出了地表侵蚀、河流下切、开挖船闸等自然及人为地貌过程诱发了岩体的基本卸荷过程,构筑了一个三峡工程永久船闸边坡岩体卸荷特性的宏观模式。实地调砥及数值计算成果均表明,地貌分析法在分析实际岩体工程中,岩体卸荷特性方面具有较强的应用价值。     

基于GIS的山区年均温分布模拟与制图

针对某些下垫面复杂而气象站(哨)有限的山区,用地理信息系统(GIS)和气温推算数学模型相结合,来获取全域的年均温分布栅格图和等值线矢量图。以龙游县为例,说明了该方法从数学模型、因素图输入、图层格式转换、图幅叠置分析等一系列技术步骤。     

Noise reduction and detection of weak, coherent signals through phase-weighted stacks

We present a new tool for efficient incoherent noise reduction for array data employing complex trace analysis. An amplitude-unbiased coherency measure is designed based on the instantaneous phase, which is used to weight the samples of an ordinary, linear stack. The result is called the phase-weighted stack (PWS) and is cleaned from incoherent noise. PWS thus permits detection of weak but coherent arrivals. The method presented can easily be extended to phase-weighted cross-correlations or be applied in the τ p domain. We illustrate and discuss the advantages and disadvantages of PWS in comparison with other coherency measures and present examples. We further show that our non-linear stacking technique enables us to detect a weak lower-mantle P -to- S conversion from a depth of approximately 840 km on array data. Hints of an 840 km discontinuity have been reported; however, such a discontinuity is not yet established due to the lack of further evidence.     

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.     

Candidate precursors: pulse-like geoelectric signals possibly related to recent seismic activity in Japan

Telemetric network observations of pulse-like geoelectric charge signals using a vertical dipole buried under the ground were undertaken at various observation sites over a wide area in Japan from April 1996. From continuous records of the signals during the six months following that, we attempted to select anomalous signals, possibly related to seismic electric activity. Special attention was paid to the elimination of noise resulting from industrial and meteorological electric activity, comparison with other electromagnetic signals in the VLF band and the relation between the precursor period and the distance from the eventual main shock that occurred in Japan. Four candidate precursor electric signals, which were not contaminated by industrial and meteorological electric activity, were then selected, of which the second appeared before the Akita-ken Nairiku-nanbu earthquake swarm, for which the maximum M = 5.9 occurred on 1996 August 11, and the third and fourth before the Chiba-ken Toho-oki earthquake, M = 6.6, on 1996 September 11. A tentative qualitative model explaining why the candidate precursory signal is related to stress building up before an earthquake is presented in terms of the electrification of gases released from fracturing rocks immediately prior to the main shock.