ALBERT

Description: The mantle transition zone is delineated by seismic discontinuities around 410 and 660 km, which are generally related to mineral phase transitions. Study of the topography of the discontinuities further constrains which phase transitions play a role and, combined with their Clapeyron slopes, what temperature variations occur. Here we use P-to-s converted seismic waves or receiver functions to study the topography of the mantle seismic discontinuities beneath Europe and the effect of subducting and ponding slabs beneath southern Europe on these features. We combine roughly 28,000 of the highest quality receiver functions into a common conversion point stack. In the topography of the discontinuity around 660 km, we find broadscale depressions of 30 km beneath central Europe and around the Mediterranean. These depressions do not correlate with any topography on the discontinuity around 410 km. Explaining these strong depressions by purely thermal effects on the dissociation of ringwoodite to bridgmanite and periclase requires unrealistically large temperature reductions. Presence of several wt % water in ringwoodite leads to a deeper phase transition, but complementary observations, such as elevated Vp/Vs ratio, attenuation and electrical conductivity, are not observed beneath central Europe. Our preferred hypothesis is the dissociation of ringwoodite into akimotoite and periclase in cold downwelling slabs at the bottom of the transition zone. The strongly negative Clapeyron slope predicted for the subsequent transition of akimotoite to bridgmanite explains the depression with a temperature reduction of 200-300 K and provides a mechanism to pond slabs in the first place.

Description: 〈span〉〈div〉Summary〈/div〉The mantle transition zone is the region between the globally observed major seismic velocity discontinuities around depths of 410 and 660 km and is important for determining the style of convection and mixing between the upper and the lower mantle. In this study, P-to-S converted waves, or receiver functions, are used to study these discontinuities beneath the Alaskan subduction zone, where the Pacific plate subducts underneath the North American plate. Previous tomographic models do not agree on the depth extent of the subducting slab, therefore improved imaging of the Earth structure underneath Alaska is required. We use 27,800 high quality radial receiver functions to make common conversion point stacks. Upper mantle velocity anomalies are accounted for by two recently published regional tomographic S-wave velocity models. Using these two tomographic models, we show that the discontinuity depths within our CCP stacks are highly dependent on the choice of velocity model, between which velocity anomaly magnitudes vary greatly. We design a quantitative test to show whether the anomalies in the velocity models are too strong or too weak, leading to over- or under-corrected discontinuity depths. We also show how this test can be used to rescale the 3D velocity corrections in order to improve the discontinuity topography maps. After applying the appropriate corrections, we find a localised thicker mantle transition zone and an uplifted 410 discontinuity, which show that the slab has clearly penetrated into the mantle transition zone. Little topography is seen on the 660 discontinuity, indicating that the slab has probably not reached the lower mantle. In the southwest, P410s arrivals have very small amplitudes or no significant arrival at all. This could be caused by water or basalt in the subducting slab, reducing the strength at the 410, or by topography on the 410 discontinuity, preventing coherent stacking. In the southeast of Alaska, a thinner mantle transition zone is observed. This area corresponds to the location of a slab window, and thinning of the mantle transition zone may be caused by hot mantle upwellings.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Seismic free oscillations, or normal modes, provide a convenient tool to calculate low-frequency seismograms in heterogeneous Earth models. A procedure called ‘full mode coupling’ allows the seismic response of the Earth to be computed. However, in order to be theoretically exact, such calculations must involve an infinite set of modes. In practice, only a finite subset of modes can be used, introducing an error into the seismograms. By systematically increasing the number of modes beyond the highest frequency of interest in the seismograms, we investigate the convergence of full-coupling calculations. As a rule-of-thumb, it is necessary to couple modes 1–2 mHz above the highest frequency of interest, although results depend upon the details of the Earth model. This is significantly higher than has previously been assumed. Observations of free oscillations also provide important constraints on the heterogeneous structure of the Earth. Historically, this inference problem has been addressed by the measurement and interpretation of splitting functions. These can be seen as secondary data extracted from low frequency seismograms. The measurement step necessitates the calculation of synthetic seismograms, but current implementations rely on approximations referred to as self- or group-coupling and do not use fully accurate seismograms. We therefore also investigate whether a systematic error might be present in currently published splitting functions. We find no evidence for any systematic bias, but published uncertainties must be doubled to properly account for the errors due to theoretical omissions and regularization in the measurement process. Correspondingly, uncertainties in results derived from splitting functions must also be increased. As is well known, density has only a weak signal in low-frequency seismograms. Our results suggest this signal is of similar scale to the true uncertainties associated with currently published splitting functions. Thus, it seems that great care must be taken in any attempt to robustly infer details of Earth's density structure using current splitting functions.〈/span〉

Description: We present the new model SP12RTS of isotropic shear-wave ( V S ) and compressional-wave ( V P ) velocity variations in the Earth's mantle. SP12RTS is derived using the same methods as employed in the construction of the shear-wave velocity models S20RTS and S40RTS, and the same data types. SP12RTS includes additional traveltime measurements of P -waves and new splitting measurements: 33 normal modes with sensitivity to the compressional-wave velocity and 9 Stoneley modes with sensitivity primarily to the lowermost mantle. Contrary to S20RTS and S40RTS, variations in V S and V P are determined without invoking scaling relationships. Lateral velocity variations in SP12RTS are parametrised using spherical harmonics up to degree 12, to focus on long-wavelength features of V S and V P and their ratio R . Large-low-velocity provinces (LLVPs) are observed for both V S and V P . SP12RTS also features an increase of R up to 2500 km depth, followed by a decrease towards the core–mantle boundary. A negative correlation between the shear-wave and bulk-sound velocity variations is observed for both the LLVPs and the surrounding mantle. These characteristics can be explained by the presence of post-perovskite or large-scale chemical heterogeneity in the lower mantle.

Description: Journal of the American Chemical Society DOI: 10.1021/jacs.5b03693

Description: As we strive to understand the most remote region of our planet, one critical area of investigation is the uppermost inner core since its structure is related to solidification of outer core material at the inner core boundary (ICB). Previous seismic studies have used body waves to show that the top ~ 100km of the inner core is isotropic. However, radial anisotropy cannot be uniquely determined by body wave observations. Alternatively, normal mode centre frequencies are sensitive to spherically symmetric Earth structure, therefore may provide a constraint on the existence of radial anisotropy in the inner core. Here we show that normal mode centre frequency measurements are compatible with 2-5% radial anisotropy in the top ~100km of the inner core with a fast direction radially outwards and a slow direction along the ICB. Given the uncertainties in the mineral physics and processes that produce anisotropy, the observed radial anisotropy may be reconciled with predictions based on either solidification processes or from texturing due to anisotropic growth.

Description: The inner core boundary marks the phase transition between the solid inner core and the fluid outer core. As the site of inner core solidification, the boundary provides insight into the processes generating the seismic structures of the inner core. In particular, it may hold the key to understanding the previously observed hemispherical asymmetry in inner core seismic velocity, anisotropy and attenuation. Here, we use a large PKiKP-PcP amplitude ratio and travel time residual dataset to investigate velocity and density contrast properties near the inner core boundary. Although hemispherical structure at the boundary has been proposed by previous inner core studies, we find no evidence for hemispheres in the amplitude ratios or travel time residuals. In addition, we find that the amplitude ratios are much larger than can be explained by variations in density contrast at the inner core boundary or core-mantle boundary. This indicates that PKiKP is primarily observed when it is anomalously large, due to focussing along its ray path. Using data in which PKiKP is not detected above the noise level, we calculate an upper estimate for the ICB density contrast of 1.2 kg m −3 . The travel time residuals display large regional variations, which differ on long and short length scales. These regions may be explained by large scale velocity variations in the F-layer just above the inner core boundary, and/or small scale topography of varying magnitude on the ICB, which also causes the large amplitudes. Such differences could arise from localised freezing and melting of the inner core.

Description: Reconciling the hemispherical structure of Earth’s inner core with its super-rotation Nature Geoscience 4, 264 (2011). doi:10.1038/ngeo1083 Authors: Lauren Waszek, Jessica Irving & Arwen Deuss Earth’s solid inner core grows through solidification of material from the fluid outer core onto its surface at rates of about 1 mm per year, freezing in core properties over time and generating an age–depth relation for the inner core. A hemispherical structure of the inner core is well-documented: an isotropic eastern hemisphere with fast seismic velocities contrasts with a slower, anisotropic western hemisphere. Independently, the inner core is reported to super-rotate at rates of up to 1° per year. Considering the slow growth, steady rotation rates of this magnitude would erase ’frozen-in’ regional variation and cannot coexist with hemispherical structure. Here, we exploit the age–depth relation, using the largest available PKIKP–PKiKP seismic travel time data set, to confirm hemispherical structure in the uppermost inner core, and to constrain the locations of the hemisphere boundaries. We find consistent eastward displacement of these boundaries with depth, from which we infer extremely slow steady inner core super-rotation of 0.1°–1° per million years. Our estimate of long-term super-rotation reconciles inner core rotation with hemispherical structure, two properties previously thought incompatible. It is in excellent agreement with geodynamo simulations, while not excluding the possibility that the much larger rotation rates inferred earlier correspond to fluctuations in inner core rotation on shorter timescales.

Description: SUMMARY Normal mode observations play an important role in studying broad-scale lateral variations in the Earth. Such studies require the calculation of accurate synthetic spectra in realistic earth models, and this remains a computationally challenging problem. Here, we describe a new implementation of the direct solution method for calculating normal mode spectra in laterally heterogeneous earth models. In this iterative direct solution method , the mode-coupling equations are solved in the frequency-domain using the preconditioned biconjugate gradient algorithm, and the time-domain solution is recovered using a numerical inverse Fourier transform. A number of example calculations are presented to demonstrate the accuracy and efficiency of the method for performing large ‘full coupling’ calculations as compared to methods based on matrix diagonalization and the traditional direct solution method.

Description: SUMMARY Normal mode observations play an important role in studying broad-scale lateral variations in the Earth. Such studies require the calculation of accurate synthetic spectra in realistic earth models, and this remains a computationally challenging problem. Here, we describe a new implementation of the direct solution method for calculating normal mode spectra in laterally heterogeneous earth models. In this iterative direct solution method , the mode-coupling equations are solved in the frequency-domain using the preconditioned biconjugate gradient algorithm, and the time-domain solution is recovered using a numerical inverse Fourier transform. A number of example calculations are presented to demonstrate the accuracy and efficiency of the method for performing large ‘full coupling’ calculations as compared to methods based on matrix diagonalization and the traditional direct solution method.