ALBERT

Description: The Yellowstone–East Snake River Plain hotspot track has been intensely studied since several decades and is widely considered to result from the interaction of a mantle plume with the North American plate. An integrated conclusive geodynamic interpretation of this extensive data set is however presently still lacking, and our knowledge of the dynamical processes beneath Yellowstone is patchy. It has been argued that the Yellowstone plume has delaminated the lower part of the thick Wyoming cratonic lithosphere. We derive an original dynamic model to quantify delamination processes related to mantle plume–lithosphere interactions. We show that fast (~300 ka) lithospheric delamination is consistent with the observed timing of formation of successive volcanic centres along the Yellowstone hotspot track and requires (i) a tensile stress regime within the whole lithosphere exceeding its failure threshold, (ii) a purely plastic rheology in the lithosphere when stresses reach this yield limit, (iii) a dense lower part of the 200 km thick Wyoming lithosphere and (iv) a decoupling melt horizon inside the median part of the lithosphere. We demonstrate that all these conditions are verified and that ~150 km large and ~100 km thick lithospheric blocks delaminate within 300 ka when the Yellowstone plume ponded below the 200 km thick Wyoming cratonic lithosphere. Furthermore, we take advantage of the available extensive regional geophysical and geological observation data sets to design a numerical 3-D upper-mantle convective model. We propose a map of the ascending convective sheets contouring the Yellowstone plume. The model further evidences the development of a counter-flow within the lower part of the lithosphere centred just above the Yellowstone mantle plume axis. This counter-flow controls the local lithospheric stress field, and as a result the trajectories of feeder dykes linking the partial melting source within the core of the mantle plume with the crust by crosscutting the lithospheric mantle. This counter-flow further explains the 50 km NE shift observed between the mantle plume axis and the present-day Yellowstone Caldera as well as the peculiar shaped crustal magma chambers.

Description: The representation of viscoelastic media in the time domain becomes more challenging with greater bandwidth of the propagating waves and number of travelled wavelengths. With the continuously increasing computational power, more extreme parameter regimes become accessible, which requires the reassessment and improvement of the standard ‘memory variable’ methods to implement attenuation in time-domain seismic wave-propagation methods. In this paper, we propose a method to minimize the error in the wavefield for a fixed complexity of the anelastic medium. This method consists of defining an appropriate misfit criterion based on a first-order analysis of how errors in the discretized medium propagate into errors in the wavefield and a simulated annealing optimization scheme to find the globally optimal parametrization. Furthermore, we derive an analytical time-stepping scheme for the memory variables that encode the strain history of the medium. Then we develop the coarse grained memory variable approach for the spectral element method (SEM) and benchmark it using the 2.5-D code AxiSEM for global body waves up to 1 Hz. Showing very good agreement with a reference solution, it also leads to a speedup of a factor of 5 in the anelastic part of the code (factor 2 in total) in this 2.5-D approach. A factor of 15 (3 in total) can be expected for the 3-D case compared to conventional implementations.

Description: We present a numerical method to compute 3-D elastic waves in fully anisotropic axisymmetric media. This method is based on a decomposition of the wave equation into a series of uncoupled 2-D equations for which the dependence of the wavefield on the azimuth can be solved analytically. Four independent equations up to quadrupole order appear as solutions for moment-tensor sources located on the symmetry axis while single forces can be accommodated by two separate solutions up to dipole order. This decomposition gives rise to an efficient solution of the 3-D wave equation in a 2-D axisymmetric medium. First, we prove the validity of the decomposition of the wavefield in the presence of general anisotropy. Then we use it to derive the reduced 2-D equations of motions and discretize them using the spectral element method. Finally, we benchmark the numerical implementation for global wave propagation at 1 Hz and consider inner core anisotropy as an application for high-frequency wave propagation in anisotropic media at frequencies up to 2 Hz.

Description: The topography of the core–mantle boundary (CMB) is directly linked to the dynamics of both the mantle and the outer core, although it is poorly constrained and understood. Recent studies have produced topography models with mutual agreement up to degree 2. A broad-band waveform inversion strategy is introduced and applied here, with relatively low computational cost and based on a first-order Born approximation. Its performance is validated using synthetic waveforms calculated in theoretical earth models that include different topography patterns with varying lateral wavelengths, from 600 to 2500 km, and magnitudes (~10 km peak-to-peak). The source–receiver geometry focuses mainly on the P diff , PKP , PcP and ScS phases. The results show that PKP branches, PcP and ScS generally perform well and in a similar fashion, while P diff yields unsatisfactory results. We investigate also how 3-D mantle correction influences the output models, and find that despite the disturbance introduced, the models recovered do not appear to be biased, provided that the 3-D model is correct. Using cross-correlated traveltimes, we derive new topography models from both P and S waves. The static corrections used to remove the mantle effect are likely to affect the inversion, compromising the agreement between models derived from P and S data. By modelling traveltime residuals starting from sensitivity kernels, we show how the simultaneous use of volumetric and boundary kernels can reduce the bias coming from mantle structures. The joint inversion approach should be the only reliable method to invert for CMB topography using absolute cross-correlation traveltimes.

Description: We present a new, computationally efficient numerical method to simulate global seismic wave propagation in realistic 3-D Earth models. We characterize the azimuthal dependence of 3-D wavefields in terms of Fourier series, such that the 3-D equations of motion reduce to an algebraic system of coupled 2-D meridian equations, which is then solved by a 2-D spectral element method (SEM). Computational efficiency of such a hybrid method stems from lateral smoothness of 3-D Earth models and axial singularity of seismic point sources, which jointly confine the Fourier modes of wavefields to a few lower orders. We show novel benchmarks for global wave solutions in 3-D structures between our method and an independent, fully discretized 3-D SEM with remarkable agreement. Performance comparisons are carried out on three state-of-the-art tomography models, with seismic period ranging from 34 s down to 11 s. It turns out that our method has run up to two orders of magnitude faster than the 3-D SEM, featured by a computational advantage expanding with seismic frequency.

Description: An important goal of computational seismology is to simulate dynamic earthquake rupture and strong ground motion in realistic models that include crustal heterogeneities and complex fault geometries. To accomplish this, we incorporate dynamic rupture modelling capabilities in a spectral element solver on unstructured meshes, the 3-D open source code SPECFEM3D, and employ state-of-the-art software for the generation of unstructured meshes of hexahedral elements. These tools provide high flexibility in representing fault systems with complex geometries, including faults with branches and non-planar faults. The domain size is extended with progressive mesh coarsening to maintain an accurate resolution of the static field. Our implementation of dynamic rupture does not affect the parallel scalability of the code. We verify our implementation by comparing our results to those of two finite element codes on benchmark problems including branched faults. Finally, we present a preliminary dynamic rupture model of the 2011 M w 9.0 Tohoku earthquake including a non-planar plate interface with heterogeneous frictional properties and initial stresses. Our simulation reproduces qualitatively the depth-dependent frequency content of the source and the large slip close to the trench observed for this earthquake.