Small-scale layered structures at the inner core boundary.

The fine-scale seismic features near the inner core boundary (ICB) provide critical insights into the thermal, chemical, and geodynamical interactions between liquid and solid cores, and may shed light on the evolution mechanism of the Earth's core. Here, we utilize a dataset of pre-critical PKiKP waveforms to constrain the fine structure at the ICB, considering the influence of various factors such as source complexity, structural anomalies in the mantle, and properties at the ICB. Our modeling suggests a sharp ICB beneath Mongolia and most of Northeast Asia, but a locally laminated ICB structure beneath Central Asia, Siberia, and part of Northeast Asia. The complex ICB structure might be explained by either the existence of a kilometer-scale thickness of mushy zone, or the localized coexistence of bcc and hcp iron phase at the ICB. We infer that there may be considerable lateral variations in the dendrites growing process at ICB, probably due to the complicated thermochemical and geodynamical interaction between the outer and inner core.

Compositional and thermal state of the lower mantle from joint 3D inversion with seismic tomography and mineral elasticity.

The compositional and thermal state of Earth's mantle provides critical constraints on the origin, evolution, and dynamics of Earth. However, the chemical composition and thermal structure of the lower mantle are still poorly understood. Particularly, the nature and origin of the two large low-shear-velocity provinces (LLSVPs) in the lowermost mantle observed from seismological studies are still debated. In this study, we inverted for the 3D chemical composition and thermal state of the lower mantle based on seismic tomography and mineral elasticity data by employing a Markov chain Monte Carlo framework. The results show a silica-enriched lower mantle with a Mg/Si ratio less than ~1.16, lower than that of the pyrolitic upper mantle (Mg/Si = 1.3). The lateral temperature distributions can be described by a Gaussian distribution with a standard deviation (SD) of 120 to 140 K at 800 to 1,600 km and the SD increases to 250 K at 2,200 km depth. However, the lateral distribution in the lowermost mantle does not follow the Gaussian distribution. We found that the velocity heterogeneities in the upper lower mantle mainly result from thermal anomalies, while those in the lowermost mantle mainly result from compositional or phase variations. The LLSVPs have higher density at the base and lower density above the depth of ~2,700 km than the ambient mantle, respectively. The LLSVPs are found to have ~500 K higher temperature, higher Bridgmanite and iron content than the ambient mantle, supporting the hypothesis that the LLSVPs may originate from an ancient basal magma ocean formed in Earth's early history.

Improving seismic remote sensing of typhoon with a three-dimensional Earth model.

Typhoon-induced P-wave microseisms can be observed using seismological arrays and analyzed for the seismic monitoring of ocean storms. This paper presents a frequency-domain beamforming (FB) method that integrates a three-dimensional (3-D) Earth model to better capture the heterogeneities in the subsurface structure, and therefore yield more accurate ray-tracing and travel-time predictions. This method is applied to the Super Typhoon Lupit (2009) using seismological array observations from the Northeast China Extended Seismic Array (NECESSArray) and high-sensitivity seismograph network in Japan (Hi-net). The results show that the localized P-wave microseism source regions based on the 3-D model are in better agreement with the theoretical source regions and typhoon centers than those based on a conventional one-dimensional (1-D) model. The significance of using a 3-D model instead of a 1-D model in the FB method is further investigated by comparing the consistency of the localization results for the two different arrays, with the localized source regions being more mutually concordant when using the 3-D model. The results demonstrate that integrating the 3-D model into the FB method improves the accuracy of locating the typhoon-induced P-wave microseism source regions.

Seismic ocean thermometry.

More than 90% of the energy trapped on Earth by increasingly abundant greenhouse gases is absorbed by the ocean. Monitoring the resulting ocean warming remains a challenging sampling problem. To complement existing point measurements, we introduce a method that infers basin-scale deep-ocean temperature changes from the travel times of sound waves that are generated by repeating earthquakes. A first implementation of this seismic ocean thermometry constrains temperature anomalies averaged across a 3000-kilometer-long section in the equatorial East Indian Ocean with a standard error of 0.0060 kelvin. Between 2005 and 2016, we find temperature fluctuations on time scales of 12 months, 6 months, and ~10 days, and we infer a decadal warming trend that substantially exceeds previous estimates.

Author Correction: Velocity and density characteristics of subducted oceanic crust and the origin of lower-mantle heterogeneities.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Velocity and density characteristics of subducted oceanic crust and the origin of lower-mantle heterogeneities.

Seismic heterogeneities detected in the lower mantle were proposed to be related to subducted oceanic crust. However, the velocity and density of subducted oceanic crust at lower-mantle conditions remain unknown. Here, we report ab initio results for the elastic properties of calcium ferrite-type phases and determine the velocities and density of oceanic crust along different mantle geotherms. We find that the subducted oceanic crust shows a large negative shear velocity anomaly at the phase boundary between stishovite and CaCl2-type silica, which is highly consistent with the feature of mid-mantle scatterers. After this phase transition in silica, subducted oceanic crust will be visible as high-velocity heterogeneities as imaged by seismic tomography. This study suggests that the presence of subducted oceanic crust could provide good explanations for some lower-mantle seismic heterogeneities with different length scales except large low shear velocity provinces (LLSVPs).

Inferring Earth's discontinuous chemical layering from the 660-kilometer boundary topography.

Topography, or depth variation, of certain interfaces in the solid Earth can provide important insights into the dynamics of our planet interior. Although the intermediate- and long-range topographic variation of the 660-kilometer boundary between Earth's upper and lower mantle is well studied, small-scale measurements are far more challenging. We found a surprising amount of topography at short length scale along the 660-kilometer boundary in certain regions using scattered P'P' seismic waves. Our observations required chemical layering in regions with high short-scale roughness. By contrast, we did not see such small-scale topography along the 410-kilometer boundary in the upper mantle. Our findings support the concept of partially blocked or imperfect circulation between the upper and lower mantle.

The Pawnee earthquake as a result of the interplay among injection, faults and foreshocks.

The Pawnee M5.8 earthquake is the largest event in Oklahoma instrument recorded history. It occurred near the edge of active seismic zones, similar to other M5+ earthquakes since 2011. It ruptured a previously unmapped fault and triggered aftershocks along a complex conjugate fault system. With a high-resolution earthquake catalog, we observe propagating foreshocks leading to the mainshock within 0.5 km distance, suggesting existence of precursory aseismic slip. At approximately 100 days before the mainshock, two M ≥ 3.5 earthquakes occurred along a mapped fault that is conjugate to the mainshock fault. At about 40 days before, two earthquakes clusters started, with one M3 earthquake occurred two days before the mainshock. The three M ≥ 3 foreshocks all produced positive Coulomb stress at the mainshock hypocenter. These foreshock activities within the conjugate fault system are near-instantaneously responding to variations in injection rates at 95% confidence. The short time delay between injection and seismicity differs from both the hypothetical expected time scale of diffusion process and the long time delay observed in this region prior to 2016, suggesting a possible role of elastic stress transfer and critical stress state of the fault. Our results suggest that the Pawnee earthquake is a result of interplay among injection, tectonic faults, and foreshocks.

Solid-liquid transitions of sodium chloride at high pressures.

We investigate solid-liquid transitions in NaCl at high pressures using molecular dynamics simulations, including the melting curve and superheating/supercooling as well as solid-solid and liquid-liquid transitions. The first-order B1-B2 (NaCl-CsCl type) transition in solid is observed at high pressures besides continuous liquid structure transitions, which are largely analogous to the B1-B2 transition in solid. The equilibrium melting temperatures (T(m)) up to megabar pressure are obtained from the solid-liquid coexistence technique and the superheating-supercooling hysteresis method. Lindemann's vibrational and Born's mechanical instabilities are found upon melting. The Lindemann frequency is calculated from the vibrational density of states. The Lindemann parameter (fractional root-mean-squared displacement) increases with pressure and approaches a constant asymptotically, similar to the Lennard-Jones system. However, the Lindemann melting relation holds for both B1 and B2 phases to high accuracy as for the Lennard-Jonesium. The B1 and B2 NaCl solids can be superheated by 0.18T(m) and 0.24T(m), and the NaCl liquid, supercooled by 0.22T(m) and 0.32T(m), respectively, at heating or cooling rates of 1 K/s and 1 K/ps. The amount of maximum superheating or supercooling and its weak pressure dependence observed for NaCl are in accord with experiments on alkali halides and with simulations on the Lennard-Jones system and Al.

Densification of silica glass at ambient pressure.

We show that densification of silica glass at ambient pressure as observed in irradiation experiments can be attributed to defect generation and subsequent structure relaxation. In our molecular dynamics simulations, defects are created by randomly removing atoms, by displacing atoms from their nominal positions in an otherwise intact glass, and by assigning certain atom excess kinetic energy (simulated ion implantation). The former forms vacancies; displacing atoms and ion implantation produce both vacancies and "interstitials." Appreciable densification is induced by these defects after equilibration of the defective glasses. The structural and vibrational properties of the densified glasses are characterized, displaying resembling features regardless of the means of densification. These results indicate that relaxation of high free-energy defects into metastable amorphous structures enriched in atomic coordination serves as a common mechanism for densification of silica glass at ambient pressure.

Rupture process of the 2004 Sumatra-Andaman earthquake.

The 26 December 2004 Sumatra-Andaman earthquake initiated slowly, with small slip and a slow rupture speed for the first 40 to 60 seconds. Then the rupture expanded at a speed of about 2.5 kilometers per second toward the north northwest, extending 1200 to 1300 kilometers along the Andaman trough. Peak displacements reached approximately 15 meters along a 600-kilometer segment of the plate boundary offshore of northwestern Sumatra and the southern Nicobar islands. Slip was less in the northern 400 to 500 kilometers of the aftershock zone, and at least some slip in that region may have occurred on a time scale beyond the seismic band.

Seismology: energy radiation from the Sumatra earthquake.

We determined the duration of high-frequency energy radiation from Indonesia's great Sumatra-Andaman earthquake (26 December 2004) to be about 500 seconds. This duration can be translated into a rupture length of about 1,200 km, which is more than twice as long as that inferred from body-wave analyses performed soon after the event. Our analysis was able rapidly to define the extent of rupture, thereby aiding the assessment of seismic hazard in the immediate future.

Sharp sides to the African superplume.

Beneath southern Africa is a large structure about 1200 kilometers across and extending obliquely 1500 kilometers upward from the core-mantle boundary with a shear velocity reduction of about 3%. Using a fortuitous set of SKS phases that travel along its eastern side, we show that the boundary of the anomaly appears to be sharp, with a width less than 50 kilometers, and is tilted outward from its center. Dynamic models that fit the seismic constraints have a dense chemical layer within an upwardly flowing thermal structure. The tilt suggests that the layer is dynamically unstable on geological time scales.