Global multiple-frequency seismic tomography using teleseismic and core-diffracted body waves

Seismic tomography is the pre-eminent tool for imaging the Earth's interior. Since the advent of this method in the 1980's, the internal structure of Earth has been vastly sampled and imaged at a variety of scales, and the resulting models have served as the primary means to investigate the processes driving our planet. Significant recent advances in seismic data acquisition and computing power have drastically progressed tomographic methods. Broad-band seismic waveforms can now be simulated up to the highest naturally occurring frequencies and consequently, measurement techniques can exploit seismic waves in their entire usable spectrum and in multiple frequencies. This dissertation revolves around aspects of global multiple-frequency seismic tomography, from retrieving and processing of large seismological data sets to explore the multi-scale structure of the earth. The centrepiece of this work is an efficient processing strategy to assemble the largest possible data sets for waveform-based tomographic inversions. Motivated by the complex but loosely-constrained structure of the lowermost mantle, we aim to increase the spatial resolution and coverage of the mantle in all depths by extracting a maximum of information from observed seismograms. We first present a method that routinely measures finite-frequency traveltimes of Pdiff waves by cross-correlating observed waveforms with synthetic seismograms across the broad-band frequency range. Large volumes of waveform data of ~ 2000 earthquakes are retrieved and pre-processed using fully automatic software built for this purpose. Synthetic seismograms for these earthquakes are calculated by semi-analytical wave propagation through a spherically symmetric earth model, to 1 Hz dominant frequency. This way, we construct one of the largest core-diffracted P wave traveltime collections so far with a total of 479,559 traveltimes in frequency passbands ranging from 30.0 to 2.7 s dominant period. Projected onto their core-grazing ray segments, the Pdiff observations recover major structural, lower-mantle heterogeneities known from existing global mantle models. An inversion framework with adaptive parameterisation and locally-adjusted regularisation is developed to accurately map the information of this data set onto the desired model parameters. This broad-band waveform inversion seamlessly incorporates the Pdiff measurements alongside a very large data set of conventional teleseismic P and PP measurements. We obtain structural heterogeneities of considerable detail in all mantle depths. The mapped features confirm several previously imaged structures. At the same time, sharper outlines for several subduction systems (e.g., Tethyan, Aegean and Farallon slabs) and uprising mantle plumes (e.g., Iceland, Afar and Tristan da Cunha) appear in our model. We trace some of these features throughout the mantle to investigate their morphological characteristics in a large (whole-mantle) context. Moreover, we report the structural findings revealed by our model. This ranges from geometries of slab complexes and subdivisions of Large Low Shear Velocity Provinces at the root of the mantle to tomographic evidence to support the existence of deep-mantle plumes beneath Iceland and Tristan da Cunha.