Inferring the elastic structure of the Earth's mantle using the spectral element method.

摘要

Mapping the elastic and anelastic structure of the Earth's mantle is crucial for understanding the temperature, composition and dynamics of our planet. Extracting the information contained in seismic waveforms is the key to constraining the elastic and anelastic structure within the Earth, and is the goal of our work. In the past quarter century, global tomography based on ray theory and first-order perturbation methods has imaged long-wavelength velocity heterogeneities of the Earth's mantle. However, the approximate techniques upon which global tomographers have traditionally relied become inadequate when dealing with crustal structure, as well as short-wavelength or large amplitude mantle heterogeneity. The spectral element method, on the other hand, permits accurate calculation of wave propagation through highly heterogeneous structures, and is computationally economical when coupled with a normal mode solution and applied to a restricted region of the earth such as the upper mantle (SEM: Capdeville et al., 2003). Importantly, SEM allows a dramatic improvement in accounting for the effects of crustal structure.;Here, we develop and apply a new hybrid method of tomography, which allows us to leverage the accuracy of SEM to model fundamental and high-mode long period (>60s) waveforms. We then present the first global model of upper mantle velocity and radial anisotropy developed using SEM. Our model, SEMum, confirms that the long-wavelength mantle structure imaged using approximate semi-analytic techniques is robust and representative of the Earth's true structure. Furthermore, it reveals structures in the upper mantle that were not clearly seen in previous global tomographic models, providing new constraints on the temperature, composition as well as flow in the mantle. We show that applying a clustering analysis to the absolute shear wave-speed profiles offers a powerful new way of exploring the relationship between surface expressions of tectonics and their elastic signature in the upper mantle. We note that this new hybrid approach to tomography can be applied to a bigger and higher-frequency dataset in order to gain new insights into the structure of the lower mantle and more robustly map seismic structure at the regional and smaller scales.