We develop and apply a full waveform inversion method that incorporates seismic data on a wide range of spatio-temporal scales, thereby constraining the details of both crustal and uppermantle structure. This is intended to further our understanding of crust-mantle interactions that shape the nature of plate tectonics, and to be a step towards improved tomographic models of strongly scale-dependent earth properties, such as attenuation and anisotropy.The inversion for detailed regional earth structure consistently embedded within a largescale model requires locally refined numerical meshes that allow us to (1) model regional wave propagation at high frequencies, and (2) capture the inferred fine-scale heterogeneities. The smallest local grid spacing sets the upper bound of the largest possible time step used to iteratively advance the seismic wave field. This limitation leads to extreme computational costs in the presence of fine-scale structure, and it inhibits the construction of full waveform tomographic models that describe earth structure on multiple scales. To reduce computational requirements to a feasible level, we design a multigrid approach based on the decomposition of a multiscale earth model with widely varying grid spacings into a family of single-scale models where the grid spacing is approximately uniform. Each of the single-scale models contains a tractable number of grid points, which ensures computational efficiency. The multi-to-singlescale decomposition is the foundation of iterative, gradient-based optimization schemes that simultaneously and consistently invert data on all scales for one multi-scale model.We demonstrate the applicability of our method in a full waveform inversion for Eurasia, with a special focus on Anatolia where coverage is particularly dense. Continental-scale structure is constrained by complete seismic waveforms in the 30-200 s period range. In addition to the well-known structural elements of the Eurasian mantle, our model reveals a variety of subtle features, such as the Armorican Massif, the Rhine Graben and the Massif Central. Anatolia is covered by waveforms with 8-200 s period, meaning that the details of both crustal and mantle structure are resolved consistently. The final model contains numerous previously undiscovered structures, including the extension-related updoming of lower-crustal material beneath the Menderes Massif in western Anatolia.Furthermore, the final model for the Anatolian region confirms estimates of crustal depth from receiver function analysis, and it accurately explains cross-correlations of ambient seismic noise at 10 s period that have not been used in the tomographic inversion. This provides strong independent evidence that detailed 3-D structure is well resolved.

No abstract

Here we present first-order results detailing the Anatolian crustal from receiver function analysis of data from approximately 300 stations within Turkey. Seismic data from the Kandilli Observatory array (KOERI; KO), the National Seismic Network of Turkey (AFAD-DAD; TU) and available IRIS data from the Northern Anatolian Fault experiment (YL) for the period between 2005 and 2010 is analysed. We calculate receiver functions in the frequency domain using water-level deconvolution. The results are analysed using a combination of H-K stacking and depth stacking to determine robust Moho conversion depths and V P /V s ratios across Anatolia. We detect a deep Moho in eastern Anatolia of up to ∼55 km, a generally normal Moho in Central Anatolia of ∼37-47 km and a thinned Moho in western Anatolia and Cyprus of ∼30 km. The V P /V s ratio across the Anatolian Plate is generally slightly elevated; regions of extremely high V P /V s ratio (>1.85) can be associated with recent volcanism in eastern and central Anatolia. High V P /V s ratio measurements (>1.85) in western Anatolia may be indicative of partial melt in the lower crust associated with regional extension.

SUMMARY Since 2004 more than 7000 km of full‐crustal reflection profiles have been collected across Australia to give a total of more than 11 000 km, providing valuable new constraints on crustal structure. A further set of hitherto unexploited results comes from 150 receiver functions distributed across the continent, mostly from portable receiver sites. These new data sets provide a dramatic increase in data coverage compared with previous studies, and reveal the complex structure of the Australian continent in considerable detail. A new comprehensive model for Moho depth across Australia and its immediate environment is developed by utilizing multiple sources of information. On‐shore and off‐shore refraction experiments are supplemented by receiver functions from a large number of portable stations and the recently augmented set of permanent stations, and Moho picks from the full suite of reflection transects. The composite data set provides a much denser sampler of most of the continent than before, though coverage remains low in the remote areas of the Simpson and Great Sandy deserts. The various data sets provide multiple estimates of the depth to Moho in many regions and the consistency between the different techniques is high. In a number of instances, differences in estimates of Moho depth can be associated with the aspects of the structure highlighted by the particular methods. The new results allow considerable refinement of the patterns of Moho depth across the continent. Some of the thinnest crust lies beneath the Archean cratons in the Pilbara and the southern part of the Yilgarn. Thick crust is encountered beneath parts of the Proterozoic in Central Australia, and beneath the Palaeozoic Lachlan fold belt in southeastern Australia. The refined data indicate a number of zones of sharp contrast in depth to Moho, notably in the southern part of Central Australia.

Stations on the Australian continent receive a rich mixture of continuous ground motion with ambient seismic noise from the surrounding oceans, and numerous small earthquakes in the earthquake belts to the north in Indonesia, and east in Tonga-Kermadec, as well as more distant source zones. The ground motion at a seismic station contains information about the structure in the vicinity of the site, and this can be exploited by applying an autocorrelation procedure to the continuous records. By creating stacked autocorrelograms of the ground motion at a single station, information on crust properties can be extracted in the form of a signal that includes the crustal reflection response convolved with the autocorrelation of the combined effect of source excitation and the instrument response. After applying suitable high-pass filtering, the reflection component can be extracted to reveal the most prominent reflectors in the lower crust, which often correspond to the reflection at the Moho. Because the reflection signal is stacked from arrivals from a wide range of slownesses, the reflection response is somewhat diffuse, but still sufficient to provide useful constraints on the local crust beneath a seismic station.Continuous vertical component records from 223 stations (permanent and temporary) across the continent have been processed using autocorrelograms of running windows 6 hr long with subsequent stacking. A distinctive pulse with a time offset between 8 and 30 s from zero is found in the autocorrelation results, with frequency content between 1.5 and 4 Hz, suggesting P-wave multiples trapped in the crust. Synthetic modelling, with control of multiple phases, shows that a local p m p phase can be recovered with the autocorrelation approach. This identification enables us to make out the depth to the most prominent crustal reflector across the continent. We obtain results that largely conform to those from previous studies using a combination of data from refraction, reflection profiles and receiver functions. This approach can be used for crustal property extraction using just vertical component records, and effective results can be obtained with temporary deployments of just a few months.