Regional tomographic inversion of the amplitude and phase of Rayleigh waves with 2-D sensitivity kernels

In this study, we test the adequacy of 2-D sensitivity kernels for fundamental-mode Rayleigh waves based on the single-scattering (Born) approximation to account for the effects of heterogeneous structure on the wavefield in a regional surface wave study. The calculated phase and amplitude data using the 2-D sensitivity kernels are compared to phase and amplitude data obtained from seismic waveforms synthesized by the pseudo-spectral method for plane Rayleigh waves propagating through heterogeneous structure. We find that the kernels can accurately predict the perturbation of the wavefield even when the size of anomaly is larger than one wavelength. The only exception is a systematic bias in the amplitude within the anomaly itself due to a site response.
An inversion method of surface wave tomography based on the sensitivity kernels is developed and applied to synthesized data obtained from a numerical simulation modelling Rayleigh wave propagation over checkerboard structure. By comparing recovered images to input structure, we illustrate that the method can almost completely recover anomalies within an array of stations when the size of the anomalies is larger than or close to one wavelength of the surface waves. Surface wave amplitude contains important information about Earth structure and should be inverted together with phase data in surface wave tomography.     

Surface wave tomography for azimuthal anisotropy in a strongly reduced parameter space

Large scale seismic anisotropy in the Earth's mantle is likely dynamically supported by the mantle's deformation; therefore, tomographic imaging of 3-D anisotropic mantle seismic velocity structure is an important tool to understand the dynamics of the mantle. While many previous studies have focused on special cases of symmetry of the elastic properties, it would be desirable for evaluation of dynamic models to allow more general axis orientation. In this study, we derive 3-D finite-frequency surface wave sensitivity kernels based on the Born approximation using a general expression for a hexagonal medium with an arbitrarily oriented symmetry axis. This results in kernels for two isotropic elastic coefficients, three coefficients that define the strength of anisotropy, and two angles that define the symmetry axis. The particular parametrization is chosen to allow for a physically meaningful method for reducing the number of parameters considered in an inversion, while allowing for straightforward integration with existing approaches for modelling body wave splitting intensity measurements. Example kernels calculated with this method reveal physical interpretations of how surface waveforms are affected by 3-D velocity perturbations, while also demonstrating the non-linearity of the problem as a function of symmetry axis orientation. The expressions are numerically validated using the spectral element method. While challenges remain in determining the best inversion scheme to appropriately handle the non-linearity, the approach derived here has great promise in allowing large scale models with resolution of both the strength and orientation of anisotropy.     

Finite-frequency vectorial tomography: a new method for high-resolution imaging of upper mantle anisotropy

Amplitude measurements of the transverse component of SKS waves, the so-called splitting intensity, can be used to formulate a non-linear inverse problem to image the 3-D variations of upper mantle anisotropy. Assuming transverse isotropy (or hexagonal symmetry), one can parametrize anisotropy by two anisotropic parameters and two angles describing the orientation of the symmetry axis. These can also be written as two collinear pseudo-vectors. The tomographic process consists of retrieving the spatial distribution of these pseudo-vectors, and thus resembles surface wave vectorial tomography. Spatial resolution results from the sensitivity of low-frequency SKS waves to seismic anisotropy off the ray path. The expressions for the 3-D sensitivity kernels for splitting intensity are derived, including the near-field contributions, and validated by comparison with a full wave equation solution based upon the finite element method. These sensitivity kernels are valid for any orientation of the symmetry axis, and thus generalize previous results that were only valid for a horizontal symmetry axis. It is shown that both lateral and vertical subwavelength variations of anisotropy can be retrieved with a dense array of broad-band stations, even in the case of vertically propagating SKS waves.     

Finite-frequency sensitivity kernels for head waves

Head waves are extremely important in determining the structure of the predominantly layered Earth. While several recent studies have shown the diffractive nature and the 3-D Fréchet kernels of finite-frequency turning waves, analogues of head waves in a continuous velocity structure, the finite-frequency effects and sensitivity kernels of head waves are yet to be carefully examined. We present the results of a numerical study focusing on the finite-frequency effects of head waves. Our model has a low-velocity layer over a high-velocity half-space and a cylindrical-shaped velocity perturbation placed beneath the interface at different locations. A 3-D finite-difference method is used to calculate synthetic waveforms. Traveltime and amplitude anomalies are measured by the cross-correlation of synthetic seismograms from models with and without the velocity perturbation and are compared to the 3-D sensitivity kernels constructed from full waveform simulations. The results show that the head wave arrival-time and amplitude are influenced by the velocity structure surrounding the ray path in a pattern that is consistent with the Fresnel zones. Unlike the 'banana–doughnut' traveltime sensitivity kernels of turning waves, the traveltime sensitivity of the head wave along the ray path below the interface is weak, but non-zero. Below the ray path, the traveltime sensitivity reaches the maximum (absolute value) at a depth that depends on the wavelength and propagation distance. The sensitivity kernels vary with the vertical velocity gradient in the lower layer, but the variation is relatively small at short propagation distances when the vertical velocity gradient is within the range of the commonly accepted values. Finally, the depression or shoaling of the interface results in increased or decreased sensitivities, respectively, beneath the interface topography.     

Surface wave tomography: global membrane waves and adjoint methods

We implement the wave equation on a spherical membrane, with a finite-difference algorithm that accounts for finite-frequency effects in the smooth-Earth approximation, and use the resulting 'membrane waves' as an analogue for surface wave propagation in the Earth. In this formulation, we derive fully numerical 2-D sensitivity kernels for phase anomaly measurements, and employ them in a preliminary tomographic application. To speed up the computation of kernels, so that it is practical to formulate the inverse problem also with respect to a laterally heterogeneous starting model, we calculate them via the adjoint method, based on backpropagation, and parallelize our software on a Linux cluster. Our method is a step forward from ray theory, as it surpasses the inherent infinite-frequency approximation. It differs from analytical Born theory in that it does not involve a far-field approximation, and accounts, in principle, for non-linear effects like multiple scattering and wave front healing. It is much cheaper than the more accurate, fully 3-D numerical solution of the Earth's equations of motion, which has not yet been applied to large-scale tomography. Our tomographic results and trade-off analysis are compatible with those found in the ray- and analytical-Born-theory approaches.     

Formalism for traveltime inversion with finite frequency effects

A simple modification of the waveform inversion formula, based on the normal mode perturbation theory, is shown to lead to a formula for traveltime anomalies. The kernel which is derived can be used for traveltime inversion with automatic inclusion of finite frequency effects. Inversion for Earth structure with such kernels will lead to better resolution estimates than ray-theoretical traveltime inversion. Examples of kernels for transverse component seismograms are shown for direct S waves, ScS , Love waves and diffracted S waves. A measure of finite frequency effects is also proposed by comparing our formula with the one from ray theory. A quantity which should be 1 in the case of ray theory is computed for the finite frequency kernels and is shown to have deviations up to about 30 per cent from 1. Therefore, the use of ray theory for long-period body waves applies incorrect weight along a ray path and may introduce a small bias to an earth model.     

Full three-dimensional tomography: a comparison between the scattering-integral and adjoint-wavefield methods

This paper analyses the computational issues of full 3-D tomography, in which the starting model as well as the model perturbation is 3-D and the sensitivity (Fréchet) kernels are calculated using the full physics of 3-D wave propagation. We compare two formulations of the structural inverse problem: the adjoint-wavefield (AW) method, which back-propagates the data from the receivers to image structure, and the scattering-integral (SI) method, which sets up the inverse problem by calculating and storing the Fréchet kernels for each data functional. The two inverse methods are closely related, but which one is more efficient depends on the overall problem geometry, particularly on the ratio of sources to receivers, as well as trade-offs in computational resources, such as the relative costs of compute cycles to data storage. We find that the SI method is computationally more efficient than the AW method in regional waveform tomography using large sets of natural sources, although it requires more storage.     

Finite-frequency sensitivity of body waves to anisotropy based upon adjoint methods

We investigate the sensitivity of finite-frequency body-wave observables to mantle anisotropy based upon kernels calculated by combining adjoint methods and spectral-element modelling of seismic wave propagation. Anisotropy is described by 21 density-normalized elastic parameters naturally involved in asymptotic wave propagation in weakly anisotropic media. In a 1-D reference model, body-wave sensitivity to anisotropy is characterized by 'banana–doughnut' kernels which exhibit large, path-dependent variations and even sign changes. P -wave traveltimes appear much more sensitive to certain azimuthally anisotropic parameters than to the usual isotropic parameters, suggesting that isotropic P -wave tomography could be significantly biased by coherent anisotropic structures, such as slabs. Because of shear-wave splitting, the common cross-correlation traveltime anomaly is not an appropriate observable for S waves propagating in anisotropic media. We propose two new observables for shear waves. The first observable is a generalized cross-correlation traveltime anomaly, and the second a generalized 'splitting intensity'. Like P waves, S waves analysed based upon these observables are generally sensitive to a large number of the 21 anisotropic parameters and show significant path-dependent variations. The specific path-geometry of SKS waves results in favourable properties for imaging based upon the splitting intensity, because it is sensitive to a smaller number of anisotropic parameters, and the region which is sampled is mainly limited to the upper mantle beneath the receiver.     

Near-field waveforms from an arbitrarily expanding, transparent spherical cavity in a prestressed medium

Summary. A simple, approximate ('transparent') solution is derived for the near-field radiation emitted by a spherical cavity expanding in an initial pure shear prestress field. Near-field terms, their propagation and decay are discussed for a variety of growth histories, and are shown to be rather insensitive to the detailed variations of rupture velocity. The transparency approximation is shown to be adequate in the near field as well as in the far field; the main effect is a slight narrowing of far-field pulses. Time domain moment estimators at close range are more reliable for the S wave than for the P wave since transverse pulses are not as strongly contaminated by near-field effects.     

The Backus–Gilbert approach to the 3-D structure in the upper mantle – II. SH and SV velocity

Summary. Phase velocity variations obtained in the previous paper are inverted by the Backus–Gilbert method for the velocity structure of the upper mantle. Spheroidal modes and toroidal modes in the period range of 125–260 s are used in the inversion. The data cannot constrain all six parameters in a transversely isotropic medium and we chose to perturb only two parameters, SH and SV velocities. SV velocities are resolved between the depths of about 200 and 400 km and SH velocities between 0 and 200 km. Resolution kernels have half-peak widths of about 200–300 km in depth, becoming broader for deeper target depths. SV velocity kernels show secondary peaks near the surface of the Earth, with widths varying from 50 to 100 km. The deeper the target depths, the wider the secondary peaks near the surface. SH velocity kernels do not possess such secondary peaks. The trade-off between SV and SH velocities is small. SV velocity is essentially determined by spheroidal modes and SH velocity by toroidal modes. Because of the broad width of the resolution kernels, the structure in the resolved region is difficult to detect from our data set; for example the differences in SV velocity structure between 250 and 350 km or the differences in SH velocity between 100 and 200 km are difficult to distinguish. Considering the horizontal resolution of about 2000 km, obtained in the previous paper, averaging kernels for 3-D structure are quite elongated in the horizontal dimension.     

Spectral aspects of the Gaussian beam method: reflection from a homogeneous half-space

Summary. Use of paraxially approximated Gaussian beams continues to be actively pursued for construction of synthetic seismograms in complicated environments. How to select the beams in the stack remains a source of difficulty which has primarily been addressed by semi-heuristic considerations. In this paper, the classical example of line-source field reflection from a homogeneous half-space that can sustain a head wave is examined from a plane-wave spectral point of view. The individual beam fields are modelled exactly by the complex source point technique, which emphasizes the complex spectral content of these wave objects. The quality of the paraxial approximation of a typical reflected (Gaussian) beam characterized by different parameters is examined from this perspective, and is compared with uniform and non-uniform asymptotics generated from the exact beam field spectral integral. With this information as background, the reflected field for a real line-source is synthesized by beam superposition. Except for the immediate vicinity of the critical reflection angle, the well-known failure of narrow paraxial beams, no matter how densely stacked, to reproduce the head wave effects is shown to be due to the inadequate spectral content of these beams and not to the failure of beam stacking per se. When the rigorous solutions are used for the narrow-waist beams, even relatively few suffice to yield agreement with the exact solution. This circumstance emphasizes the importance of fully understanding the spectral implications of various beam stacking schemes.     

Profile response of a lacustrine multiple barred nearshore to a sequence of storm events

Profile change in a lacustrine multiple-barred nearshore was investigated over the ice-free season of 2001/2002 at Burley Beach on the southeastern shore of Lake Huron in order to identify the feedback mechanisms between the pre-existing morphology and the wave forcing and the consequence of those feedbacks to the behaviour of the nearshore environment. The characteristics of the offshore wave field were monitored using a Falmouth Scientific combined 3D-ACM wave recorder and pressure transducer. Supplemental wave data were downloaded from a 3-m discus buoy operated by the National Data Buoy Center, 75 km to the NW of the study site. The three nearshore bars were in a quasi-equilibrium state through a large part of the ice-free season, with dramatic changes occurring during relatively moderate storm events that followed much larger storms in late October. A comparison of the incident wave field with changes in the nearshore profile through canonical correlation analysis indicates that the morphology responds to the distribution of the significant, root-mean-square (rms) and average wave heights between surveys. The threshold between bar decay and onshore bar migration and growth is associated with the onset of breaking of the rms wave at the bar crest (H_rms h_cr⁻¹ ≈0.3–0.4). The threshold between onshore and offshore migration is associated with the onset of breaking of the average wave at the bar crest (H_avg h_cr⁻¹ ≈0.3–0.4), coincident with complete dissipation of the significant wave over the lakeward slope of the bar (H_s h_cr⁻¹ >0.6). Inshore wave data collected during an instrumented study at the same site revealed that the middle and inner bars remained at the threshold of onshore and offshore migration over a wide range of offshore significant wave heights (0.8 to 2.4 m) prior to the October storms. This self-organised equilibrium is a result of changes to the incident wave distribution through breaking on the outer bar. It is concluded that the prediction of bar response requires an understanding of the feedback between the bar and the local wave distribution in addition to an understanding of the feedback associated with the profile as a whole.     

Exact synthetic seismograms for an inhomogeneity in a layered elastic half-space

Summary. The propagation of a pulsed elastic wave in the following geometry is considered. An elastic half-space has a surface layer of a different material and the layer furthermore contains a bounded 3-D inhomogeneity. The exciting source is an explosion, modelled as an isotropic pressure point source with Gaussian behaviour in time.
The time-harmonic problem is solved using the null field approach (the T matrix method), and a frequency integral then gives the time-domain response. The main tools of the null field approach are integral representations containing the free space Green's dyadic, expansions in plane and spherical vector wave functions, and transformations between plane and spherical vector wave functions. It should be noted that the null field approach gives the solution to the full elastodynamic equations with, in principle, an arbitrarily high accuracy. Thus no ray approximations or the like are used. The main numerical limitation is that only low and intermediate frequencies, in the sense that the diameter of the inhomogeneity can only be a few wavelengths, can be considered.
The numerical examples show synthetic seismograms consisting of data from 15 observation points at increasing distances from the source. The normal component of the velocity field is computed and the anomalous field due to the inhomogeneity is sometimes shown separately. The shape of the inhomogeneity, the location and depth of the source, and the material parameters are all varied to illustrate the relative importance of the various parameters. Several specific wave types can be identified in the seismograms: Rayleigh waves, direct and reflected P -waves, and head waves.     

On elastic wave calculations in a sphere using moment tensor sources

Summary. Propagator matrix solutions to the elastic equations of motion in spherically symmetric, inhomogeneous media with moment tensor sources are recast into a simple and intuitively satisfying form which is applicable to both exact and approximate calculations. The transformed expression benefits from the analogous equations of normal mode excitation, while clearly distinguishing the finer partitions of the displacement field and the more flexible boundary conditions that body wave formulations provide. I believe that this new representation, because of its many advantages, should be favoured as the foundation for elastic wave calculations in a sphere.     

Inversion of Love-wave data for velocity and anelasticity using exact kernels

Summary. The general problem of inverting Love-wave dispersion and amplitude data to obtain a velocity and Q_s structure is considered. A formulation is used which incorporates attenuation into the Haskell-Thompson matrix method in an exact manner and thus retains the inherent non-linearity in the anelasticity. The resulting exact inversion kernels allow simultaneous inversion for velocity and intrinsic attenuation parameters. The method is applied to synthetic data which allows a comparison to be made with inexact kernels. The results indicate that the use of inexact kernels may introduce spurious oscillations into the Q_s structure and that a simultaneous inversion can be more stable than inverting for velocity alone.     

基于SWAN模式的近10年南海北部波浪能资源研究

利用QN(QuikSCAT/NCEP)混合风场驱动第三代海浪数值模式SWAN(SimulatingWaves Nearshore),得到中国南海北部1999年8月—2009年7月的海浪场,对该海域的波浪能资源进行综合分析,研究发现中国南海北部海域蕴藏着丰富的、适宜开发的波浪能资源.     

Electromagnetic fields in a halfspace scattered by a buried sphere by the method of transformation of local elements

Summary. A problem in modelling electromagnetic fields used in exploration geophysics is treated mathematically. Analytical expressions are obtained for the electric field due to a harmonic current in a horizontal loop on or above a conducting ground in which is buried a conductive and permeable sphere (ore body). The loop is coaxial with the sphere. For a general time-varying current in the loop, the analysis is carried to the stage where a Fourier inversion can be used to obtain readily the electric field in the time-domain. A new relationship between spherical and cylindrical wave functions is obtained as a transformation of local elements.
Solution of this problem has not been presented before in this form. Lee's solution of 1975 which uses an integral-equation formulation treats a similar problem without taking account of differences in magnetic permeability. The effects of magnetic permeability may have important and useful implications for geophysical explorations.     

Global multiresolution models of surface wave propagation: comparing equivalently regularized Born and ray theoretical solutions

I invert a large set of teleseismic phase-anomaly observations, to derive tomographic maps of fundamental-mode surface wave phase velocity, first via ray theory, then accounting for finite-frequency effects through scattering theory, in the far-field approximation and neglecting mode coupling. I make use of a multiple-resolution pixel parametrization which, in the assumption of sufficient data coverage, should be adequate to represent strongly oscillatory Fréchet kernels. The parametrization is finer over North America, a region particularly well covered by the data. For each surface-wave mode where phase-anomaly observations are available, I derive a wide spectrum of plausible, differently damped solutions; I then conduct a trade-off analysis, and select as optimal solution model the one associated with the point of maximum curvature on the trade-off curve. I repeat this exercise in both theoretical frameworks, to find that selected scattering and ray theoretical phase-velocity maps are coincident in pattern, and differ only slightly in amplitude.     

Morphodynamic evolution of a coastal lagoon entrance during swash overwash

New field observations of beach berm growth resulting in the sub-aerial closure of an intermittently closed and open lake or lagoon are presented. These perched estuarine systems frequently exhibit ephemeral entrances that respond to process forcing almost instantaneously on a geomorphic timescale, with closure by a supra-tidal berm of depositional origin. The observations were made following a mechanical opening of the entrance and show very rapid vertical growth of an initially low beach berm as a result of swash overtopping and sediment overwash. Very little simultaneous seaward progradation was observed. Sediment overwash volumes were obtained from survey data over a period of 12 days either side of spring tide and the observations also provide an accurate measure of swash sediment transport at the berm crest. A process-based parametric model is developed through robust parameterisations of wave run-up, wave run-up distributions and sediment transport and is applied to predict the total overwash transport into and infilling the estuary entrance. The model is tested against the field data and compares well with the observations. While the field data are site-specific, the proposed modelling framework represents a first step in modelling the complex growth and recovery of natural beach berms in broad-scale morphodynamic models.     

The Backus–Gilbert theory for piecewise continuous structures with variable discontinuity levels and its application to the magnetotelluric inverse problem

Summary. The Backus–Gilbert theory is extended to the case when the models are piecewise continuous vector functions of depth with variable discontinuity locations. In addition, some of the layers may be represented by linear combinations of known functions. For such layers only a finite number of discrete parameters is to be determined. The iterative process for obtaining a model satisfying the data is convergent, the numerical procedure by which the iterations are performed being equivalent to the method of spectral decomposition for continuous structures. The method of obtaining the Fréchet kernels by using the first perturbation of the differential system satisfied by the corresponding functionals is shown to be valid.
The theory is applied to the magnetotelluric problem, Fréchet kernels being calculated for isotropic and non-isotropic structures.
A few numerical examples are described.