Deconvolution of teleseismic recordings for mantle structure

Delineation of detailed mantle structure frequently requires the separation of source signature and structural response from seismograms recorded at teleseismic distances. This deconvolution problem can be posed in a log-spectral domain where the operation of time-domain convolution is reduced to an additive form. The introduction of multiple events recorded at many stations leads to a system of consistency equations that must be honoured by both the source time functions and the impulse responses associated with propagation paths between sources and receivers. The system is inherently singular, and stabilization is accomplished through the supply of an initial estimate of the source time function. Although alternative choices exist, an effective estimate is derived from the eigenimage associated with the largest eigenvalue in a singular-value decomposition of the suite of aligned seismograms corresponding to a given event. The relation of the deconvolution scheme to simultaneous least-squares deconvolution is examined. Application of the methodology to broadband teleseismic P waveforms recorded on the Canadian National Seismograph Network demonstrates the retrieval of effective Green's functions including secondary phases associated with upper-mantle structure.     

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.     

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)     

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.