International audienceWe present forward and adjoint spectral-element simulations of coupled acoustic and (an)elastic seismic wave propagation on fully unstructured hexahedral meshes. Simulations benefit from recent advances in hexahedral meshing, load balancing and software optimization. Meshing may be accomplished using a mesh generation tool kit such as CUBIT, and load balancing is facilitated by graph partitioning based on the SCOTCH library. Coupling between fluid and solid regions is incorporated in a straightforward fashion using domain decomposition. Topography, bathymetry and Moho undulations may be readily included in the mesh, and physical dispersion and attenuation associated with anelasticity are accounted for using a series of standard linear solids. Finite-frequency Fr'echet derivatives are calculated using adjoint methods in both fluid and solid domains. The software is benchmarked for a layercake model. We present various examples of fully unstructured meshes, snapshots of wavefields and finite-frequency kernels generated by Version 2.0 'Sesame' of our widely used open source spectral-element package SPECFEM3D

No abstract

No abstract

We built a 3D seismic model of the Po Plain and neighboring regions of northern Italy, covering altogether an area about 600 km by 300 km with an approximately 1 km spaced grid. We started by collecting an extensive and diverse set of geological and geophysical data, including seismic reflection and refraction profiles, borehole logs, and available geological information. Major geological boundaries and discontinuities have thus been identified and mapped into the model. We used kriging to interpolate the geographically sparse information into continuous surfaces delimiting geological bodies with laterally varying thickness. Seismic-wave properties have been assigned to each unit using a rule-based system and, V P , V S , and ρ derived from other studies. Sedimentary strata, although with varying levels of compaction and hence material properties, may locally reach a thickness of 15 km and give rise to significant effects in seismic-wave propagation. We have used our new model to compute the seismic response for two recent earthquakes, to test its performance. Results show that the 3D model reproduces the large amplitude and the long duration of shaking seen in the observed waveforms recorded on sediments, whereas paths outside the basin may be well fit by more homogeneous (1D) hard rock structure. We conclude that the new model is suited for simulation of wave propagation, mostly for T > 3 s, and may serve well as a constraint for earthquake location and further improvements via body-or surface-wave inversion.

(2014), A 3-D spectral-element and frequency-wave number hybrid method for high-resolution seismic array imaging, Geophys. Res. Lett., 41, 7025-7034, doi:10.1002 Abstract We present a three-dimensional (3-D) hybrid method that interfaces the spectral-element method (SEM) with the frequency-wave number (FK) technique to model the propagation of teleseismic plane waves beneath seismic arrays. The accuracy of the resulting 3-D SEM-FK hybrid method is benchmarked against semianalytical FK solutions for 1-D models. The accuracy of 2.5-D modeling based on 2-D SEM-FK hybrid method is also investigated through comparisons to this 3-D hybrid method. Synthetic examples for structural models of the Alaska subduction zone and the central Tibet crust show that this method is capable of accurately capturing interactions between incident plane waves and local heterogeneities. This hybrid method presents an essential tool for the receiver function and scattering imaging community to verify and further improve their techniques. These numerical examples also show the promising future of the 3-D SEM-FK hybrid method in high-resolution regional seismic imaging based on waveform inversions of converted/scattered waves recorded by seismic array.