Search

To describe the energy transport in the seismic coda, we introduce a system of radiative transfer equations for coupled surface and body waves in a scalar approximation. Our model is based on the Helmholtz equation in a half-space geometry with mixed boundary conditions. In this model, Green's function can be represented as a sum of body waves and surface waves, which mimics the situation on Earth. In a first step, we study the single-scattering problem for point-like objects in the Born approximation. Using the assumption that the phase of body waves is randomized by surface reflection or by interaction with the scatterers, we show that it becomes possible to define, in the usual manner, the cross-sections for surface-to-body and body-to-surface scattering. Adopting the independent scattering approximation, we then define the scattering mean free paths of body and surface waves including the coupling between the two types of waves. Using a phenomenological approach, we then derive a set of coupled transport equations satisfied by the specific energy density of surface and body waves in a medium containing a homogeneous distribution of point scatterers. In our model, the scattering mean free path of body waves is depth dependent as a consequence of the body-to-surface coupling. We demonstrate that an equipartition between surface and body waves is established at long lapse-time, with a ratio which is predicted by usual mode counting arguments. We derive a diffusion approximation from the set of transport equations and show that the diffusivity is both anisotropic and depth dependent. The physical origin of the two properties is discussed. Finally, we present Monte-Carlo solutions of the transport equations which illustrate the convergence towards equipartition at long lapse-time as well as the importance of the coupling between surface and body waves in the generation of coda waves.

A narrow rectilinear valley in the French Pyrenees, affected in the past by damaging earthquakes, has been chosen as a test site for soil response characterization. The main purpose of this initiative was to compare experimental and numerical approaches. A temporary network of 10 stations has been deployed along and across the valley during two years; parallel various experiments have been conducted, in particular ambient noise recording, and seismic profiles with active sources for structure determination at the 10 sites. Classical observables have been measured for site amplification evaluation, such as spectral ratios of horizontal or vertical motions between site and reference stations using direct S waves and S coda, and spectral ratios between horizontal and vertical (H/V) motions at single stations using noise and S-coda records. Vertical shear-velocity profiles at the stations have first been obtained from a joint inversion of Rayleigh wave dispersion curves and ellipticity. They have subsequently been used to model the H/V spectral ratios of noise data from synthetic seismograms, the H/V ratio of S-coda waves based on equipartition theory, and the 3D seismic response of the basin using the spectral element method. General good agreement is found between simulations and observations. The 3D simulation reveals that topography has a much lower contribution to site effects than sedimentary filling, except at the narrow ridge crests. We find clear evidence of a basin edge effect, with an increase of the amplitude of ground motion at some distance from the edge inside the basin and a decrease immediately at the slope foot.

This paper investigates the reconstruction of elastic Green's function from the cross-correlation of waves excited by random noise in the context of scattering theory. Using a general operator equation, -the resolvent formula-, Green's function reconstruction is established when the noise sources satisfy an equipartition condition. In an inhomogeneous medium, the operator formalism leads to generalized forms of optical theorem involving the off-shell $T$-matrix of elastic waves, which describes scattering in the near-field. The role of temporal absorption in the formulation of the theorem is discussed. Previously established symmetry and reciprocity relations involving the on-shell $T$-matrix are recovered in the usual far-field and infinitesimal absorption limits. The theory is applied to a point scattering model for elastic waves. The $T$-matrix of the point scatterer incorporating all recurrent scattering loops is obtained by a regularization procedure. The physical significance of the point scatterer is discussed. In particular this model satisfies the off-shell version of the generalized optical theorem. The link between equipartition and Green's function reconstruction in a scattering medium is discussed.

We describe a time-resolved monitoring technique for heterogeneous media. Our approach is based on the spatial variations of the cross-coherence of coda waveforms acquired at fixed positions but at different dates. To locate and characterize a weak change that occurred between successive acquisitions, we use a maximum likelihood approach combined with a diffusive propagation model. We illustrate this technique, called LOCADIFF, with numerical simulations. In several illustrative examples, we show that the change can be located with a precision of a few wavelengths and its effective scattering cross-section can be retrieved. The precision of the method depending on the number of source receiver pairs, time window in the coda, and errors in the propagation model is investigated. Limits of applications of the technique to real-world experiments are discussed., Comment: 11 pages, 14 figures, 1 table

In this paper, we present a nonparametric method to estimate the heterogeneity of a random medium from the angular distribution of intensity transmitted through a slab of random material. Our approach is based on the modeling of forward multiple scattering using Compound Poisson Processes on compact Lie groups. The estimation technique is validated through numerical simulations based on radiative transfer theory., Comment: 23 pages, 8 figures, 21 references

We have studied the partition of shear, compressional and kinetic energies in the coda of ten earthquakes recorded on a dense array, located at Pinyon Flats Observatory (PFO), California. We observe a clear stabilization of the shear to compressional ($W^s/W^p$) energy ratio in the coda, with an average value of about 2.8. The ratio between the vertical and horizontal kinetic energies ($V^2/H^2$) can be measured from 5 to 25Hz and shows an abrupt transitionfrom 0.1 in the 5-10Hz band, to about 0.8 in the 15-25Hz band. These measured values are in sharp contrast with the theoretical prediction for equipartitioned elastic waves in a homogeneous half-space. To explain these observations, we have developed a theory of equipartition in a layered elastic half-space. Using a rigorous spectral decomposition of the elastic wave equation, we define equipartition as a white noise distributed over the complete set of eigenfunctions. The theory predicts that close to the resonance frequency of a low-velocity layer, the ratio between shear and compressional energies strongly decreases. Using a detailed model of the subsurface at PFO, this conterintuitive result is found to be in good qualitative and quantitative agreement with the observations.

We propose to model the uppermost inner core as an aggregate of randomly oriented anisotropic ``patches''. A patch is defined as an assemblage of a possibly large number of crystals with identically oriented crystallographic axes. This simple model accounts for the observed velocity isotropy of short period body waves, and offers a reasonable physical interpretation for the scatterers detected at the top of the inner core. From rigorous multiple scattering modeling of seismic wave propagation through the aggregate, we obtain strong constraints on both the size and the elastic constants of iron patches. We perform a systematic search for iron models compatible with measured seismic velocities and attenuations. An iron model is characterized by its symmetry (cubic or hexagonal), elastic constants, and patch size. Independent of the crystal symmetry, we infer a most likely size of patch of the order of 400 m. Recent {\it bcc} iron models from the literature are in very good agreement with the most probable elastic constants of cubic crystals found in our inversion. Our study (1) suggests that the presence of melt may not be required to explain the low shear wavespeeds in the inner core and (2) supports the recent experimental results on the stability of cubic iron in the inner core, at least in its upper part.

The spectral decomposition of the elastic wave operator in a layered isotropic half-space is derived by means of standard functional analytic methods. Particular attention is paid to the coupled $P$-$SV$ waves. The problem is formulated directly in terms of displacements which leads to a $2 \times 2$ Sturm-Liouville system. The resolvent kernel (Green function) is expressed in terms of simple plane-wave solutions. Application of Stone's formula leads naturally to eigenfunction expansions in terms of generalized eigenvectors with oscillatory behavior at infinity. The generalized normal mode expansion is employed to define a diffuse field as a white noise process in modal space. By means of a Wigner transform, we calculate vertical to horizontal kinetic energy ratios in layered media, as a function of depth and frequency. Several illustrative examples are considered including energy ratios near a free surface, in the presence of a soft layer. Numerical comparisons between the generalized eigenfunction summation and a classical locked-mode approximation demonstrate the validity of the approach. The impact of the local velocity structure on the energy partitioning of a diffuse field is illustrated.

When averaged over sources or disorder, cross-correlation of diffuse fields yield the Green's function between two passive sensors. This technique is applied to elastic ultrasonic waves in an open scattering slab mimicking seismic waves in the Earth's crust. It appears that the Rayleigh wave reconstruction depends on the scattering properties of the elastic slab. Special attention is paid to the specific role of bulk to Rayleigh wave coupling, which may result in unexpected phenomena like a persistent time-asymmetry in the diffuse regime.