Jeroen Tromp

Measurements and global models of surface wave propagation

A new technique for making single-station phase velocity measurements is developed and applied to a large number of globally recorded Rayleigh and Love waves in the period range 35–150 s. The method is based on phase-matched filter theory and iteratively suppresses the effect of interfering overtones by minimizing residual dispersion. The model surface wave signal is described by its amplitude and apparent phase velocity, both of which are parameterized in terms of smooth B-spline functions of frequency. A misfit function is constructed which represents the difference between the model and observed waveforms, and the optimal spline coefficients are estimated in an iterative misfit minimization algorithm. In order to eliminate cycle skips in the measurements of phase at short periods, the waveforms are first matched at long periods, and the frequency range is gradually extended to include higher frequencies. The application of the algorithm to records from the Global Seismographic Network, using earthquakes in the Harvard centroid-moment tensor catalog, results in the determination of more than 50,000 high-quality dispersion curves. The observed variations in measured dispersion for pairwise similar paths are used to estimate realistic uncertainties in the data. Phase delays at discrete periods are inverted for global maps of variations in phase velocity expanded in spherical harmonics up to degree 40. A realistic resolution test indicates that structures are well recovered up to at least degree 20. The new phase velocity maps explain 70–96% of the observed variance in phase residuals, reflecting the high internal consistency of the dispersion measurements.

Adjoint Tomography of the Southern California Crust

Crustal Details Revealed In seismic tomography, a large collection of data representing paths through Earth are inverted to provide an analysis of variation of density in which errors are minimized. Typically, the inversion starts with a simple layered model of the tomographic region. Tape et al. (p. 988) show how, starting with a three-dimensional model, based on synthetic seismograms, an improved iterative inversion approach can lead to a much more detailed view of a region. Using the rich data for Southern California, the model reveals details of the geologic history of the crust in this region. Analysis of seismic data using a more realistic crustal model reveals detailed variations in density beneath southern California. Using an inversion strategy based on adjoint methods, we developed a three-dimensional seismological model of the southern California crust. The resulting model involved 16 tomographic iterations, which required 6800 wavefield simulations and a total of 0.8 million central processing unit hours. The new crustal model reveals strong heterogeneity, including local changes of ±30% with respect to the initial three-dimensional model provided by the Southern California Earthquake Center. The model illuminates shallow features such as sedimentary basins and compositional contrasts across faults. It also reveals crustal features at depth that aid in the tectonic reconstruction of southern California, such as subduction-captured oceanic crustal fragments. The new model enables more realistic and accurate assessments of seismic hazard.

Simulations of Ground Motion in the Los Angeles Basin Based upon the Spectral-Element Method

We use the spectral-element method to simulate ground motion generated by two recent and well-recorded small earthquakes in the Los Angeles basin. Simulations are performed using a new sedimentary basin model that is constrained by hundreds of petroleum-industry well logs and more than 20,000 km of seismic reflection profiles. The numerical simulations account for 3D variations of seismic-wave speeds and density, topography and bathymetry, and attenuation. Simulations for the 9 September 2001 M_w 4.2 Hollywood earthquake and the 3 September 2002 M_w 4.2 Yorba Linda earthquake demonstrate that the combination of a detailed sedimentary basin model and an accurate numerical technique facilitates the simulation of ground motion at periods of 2 sec and longer inside the basin model and 6 sec and longer in the regional model. Peak ground displacement, velocity, and acceleration maps illustrate that significant amplification occurs in the basin.

Finite-Frequency Kernels Based on Adjoint Methods

We derive the adjoint equations associated with the calculation of Frechet derivatives for tomographic inversions based upon a Lagrange multiplier method. The Frechet derivative of an objective function χ(m), where m denotes the Earth model, may be written in the generic form δχ = ∫ K_m(x) δ ln m(x) d^3x, where δ ln m = δm/m denotes the relative model perturbation and K_m the associated 3D sensitivity or Frechet kernel. Complications due to artificial absorbing boundaries for regional simulations as well as finite sources are accommodated. We construct the 3D finite-frequency “banana-doughnut” kernel K_m by simultaneously computing the so-called “adjoint” wave field forward in time and reconstructing the regular wave field backward in time. The adjoint wave field is produced by using time- reversed signals at the receivers as fictitious, simultaneous sources, while the regular wave field is reconstructed on the fly by propagating the last frame of the wave field, saved by a previous forward simulation, backward in time. The approach is based on the spectral-element method, and only two simulations are needed to produce the 3D finite-frequency sensitivity kernels. The method is applied to 1D and 3D regional models. Various 3D shear- and compressional-wave sensitivity kernels are presented for different regional body- and surface-wave arrivals in the seismograms. These kernels illustrate the sensitivity of the observations to the structural parameters and form the basis of fully 3D tomographic inversions.

Simulation of anisotropic wave propagation based upon a spectral element method

We introduce a numerical approach for modeling elastic wave propagation in 2-D and 3-D fully anisotropic media based upon a spectral element method. The technique solves a weak formulation of the wave equation, which is discretized using a high-order polynomial representation on a finite element mesh. For isotropic media, the spectral element method is known for its high degree of accuracy, its ability to handle complex model geometries, and its low computational cost. We show that the method can be extended to fully anisotropic media. The mass matrix obtained is diagonal by construction, which leads to a very efficient fully explicit solver. We demonstrate the accuracy of the method by comparing it against a known analytical solution for a 2-D transversely isotropic test case, and by comparing its predictions against those based upon a finite difference method for a 2-D heterogeneous, anisotropic medium. We show its generality and its flexibility by modeling wave propagation in a 3-D transversely isotropic medium with a symmetry axis tilted relative to the axes of the grid.

The global seismographic network surpasses its design goal

This year, the Global Seismographic Network (GSN) surpassed its 128-station design goal for uniform worldwide coverage of the Earth. A total of 136 GSN stations are now sited from the South Pole to Siberia, and from the Amazon Basin to the sea floor of the northeast Pacific Ocean—in cooperation with over 100 host organizations and seismic networks in 59 countries worldwide (Figure 1). Established in 1986 by the Incorporated Research Institutions for Seismology (IRIS) to replace the obsolete, analog Worldwide Standardized Seismograph Network (WWSSN),the GSN continues a tradition in global seismology that dates back more than a century to the network of Milne seismographs that initially spanned the globe. The GSN is a permanent network of state-of-the-art seismological and geophysical sensors connected by available telecommunications to serve as a multi-use scientific facility and societal resource for scientific research, environmental monitoring, and education for our national and international community.

Wave propagation near a fluid‐solid interface: A spectral‐element approach

We introduce a spectral-element method for modeling wave propagation in media with both fluid (acoustic) and solid (elastic) regions, as for instance in offshore seismic experiments. The problem is formulated in terms of displacement in elastic regions and a velocity potential in acoustic regions. Matching between domains is implemented based upon an interface integral in the framework of an explicit prediction-multicorrection staggered time scheme. The formulation results in a mass matrix that is diagonal by construction. The scheme exhibits high accuracy for a 2-D test case with known analytical solution. The method is robust in the case of strong topography at the fluid-solid interface and is a good alternative to classical techniques, such as finite differencing.

Spectral-Element Moment Tensor Inversions for Earthquakes in Southern California

We have developed and implemented an automated moment tensor inversion procedure to determine source parameters for southern California earthquakes. The method is based upon spectral-element simulations of regional seismic wave propagation in an integrated 3D southern California velocity model. Sensitivity to source parameters is determined by numerically calculating the Frechet derivatives required for the moment tensor inversion. We minimize a waveform misfit function, and allow limited time shifts between data and corresponding synthetics to accommodate additional 3D heterogeneity not included in our model. The technique is applied to three recent southern California earthquakes: the 9 September 2001, M_L 4.2 Hollywood event, the 22 February 2003, M_L 5.4 Big Bear event, and the 14 December 2001, M_L 4.0 Diamond Bar event. Using about half of the available three-component data at periods of 6 sec and longer, we obtain focal mechanisms, depths, and moment magnitudes that are generally in good agreement with estimates based upon traditional body-wave and surface-wave inversions.

The Spectral‐Element Method in Seismology

We present the main properties of the spectral-element method, which is well suited for numerical calculations of synthetic seismograms for three-dimensional Earth models. The technique is based upon a weak formulation of the equations of motion and combines the flexibility of a finite-element method with the accuracy of a pseudospectral method. The mesh is composed of hexahedral elements and honors the main discontinuities in the Earth model. The displacement vector is expressed in each element in terms of high-degree Lagrange interpolants, and integrals are computed based upon Gauss-Lobatto-Legendre quadrature, which leads to an exactly diagonal mass matrix and therefore drastically simplifies the algorithm. We use a fluid-solid coupling formulation that does not require iterations at the core-mantle or inner-core boundaries. The method is efficiently implemented on parallel computers with distributed memory based upon a message-passing methodology. We present two large-scale simulations for a realistic three-dimensional Earth model computed on the Japanese Earth Simulator at periods of 5 s and longer.

A 14.6 billion degrees of freedom, 5 teraflops, 2.5 terabyte earthquake simulation on the Earth Simulator

We use 1944 processors of the Earth Simulator to model seismic wave propagation resulting from large earthquakes. Simulations are conducted based upon the spectral-element method, a high-degree finite-element technique with an exactly diagonal mass matrix. We use a very large mesh with 5.5 billion grid points (14.6 billion degrees of freedom). We include the full complexity of the Earth, i.e., a three-dimensional wave-speed and density structure, a 3-D crustal model, ellipticity as well as topography and bathymetry. A total of 2.5 terabytes of memory is needed. Our implementation is purely based upon MPI, with loop vectorization on each processor. We obtain an excellent vectorization ratio of 99.3%, and we reach a performance of 5 teraflops (30% of the peak performance) on 38% of the machine. The very high resolution of the mesh allows us to perform fully three-dimensional calculations at seismic periods as low as 5 seconds.