Seismic evidence for a 1000 km mantle discontinuity under the Pacific.

Seismic discontinuities in the mantle are indicators of its thermo-chemical state and offer clues to its dynamics. Ray-based seismic methods, though limited by the approximations made, have mapped mantle transition zone discontinuities in detail, but have yet to offer definitive conclusions on the presence and nature of mid-mantle discontinuities. Here, we show how to use a wave-equation-based imaging method, reverse-time migration of precursors to surface-reflected seismic body waves, to uncover both mantle transition zone and mid-mantle discontinuities, and interpret their physical nature. We observe a thinned mantle transition zone southeast of Hawaii, and a reduction in impedance contrast around 410 km depth in the same area, suggesting a hotter-than-average mantle in the region. Here, we furthermore reveal a 4000-5000 km-wide reflector in new images of the mid mantle below the central Pacific, at 950-1050 km depth. This deep discontinuity exhibits strong topography and generates reflections with polarity opposite to those originating at the 660 km discontinuity, implying an impedance reversal near 1000 km. We link this mid-mantle discontinuity to the upper reaches of deflected mantle plumes upwelling in the region. Reverse-time migration full-waveform imaging is a powerful approach to imaging Earth's interior, capable of broadening our understanding of its structure and dynamics and shrinking modeling uncertainties.

The Debiased Spatial Whittle likelihood.

We provide a computationally and statistically efficient method for estimating the parameters of a stochastic covariance model observed on a regular spatial grid in any number of dimensions. Our proposed method, which we call the Debiased Spatial Whittle likelihood, makes important corrections to the well-known Whittle likelihood to account for large sources of bias caused by boundary effects and aliasing. We generalize the approach to flexibly allow for significant volumes of missing data including those with lower-dimensional substructure, and for irregular sampling boundaries. We build a theoretical framework under relatively weak assumptions which ensures consistency and asymptotic normality in numerous practical settings including missing data and non-Gaussian processes. We also extend our consistency results to multivariate processes. We provide detailed implementation guidelines which ensure the estimation procedure can be conducted in O ( n log n ) operations, where n is the number of points of the encapsulating rectangular grid, thus keeping the computational scalability of Fourier and Whittle-based methods for large data sets. We validate our procedure over a range of simulated and realworld settings, and compare with state-of-the-art alternatives, demonstrating the enduring practical appeal of Fourier-based methods, provided they are corrected by the procedures developed in this paper.

The origin of secondary microseism Love waves.

The interaction of ocean surface waves produces pressure fluctuations at the seafloor capable of generating seismic waves in the solid Earth. The accepted mechanism satisfactorily explains secondary microseisms of the Rayleigh type, but it does not justify the presence of transversely polarized Love waves, nevertheless widely observed. An explanation for two-thirds of the worldwide ambient wave field has been wanting for over a century. Using numerical simulations of global-scale seismic wave propagation at unprecedented high frequency, here we explain the origin of secondary microseism Love waves. A small fraction of those is generated by boundary force-splitting at bathymetric inclines, but the majority is generated by the interaction of the seismic wave field with three-dimensional heterogeneity within the Earth. We present evidence for an ergodic model that explains observed seismic wave partitioning, a requirement for full-wave field ambient-noise tomography to account for realistic source distributions.

Imaging the Galápagos mantle plume with an unconventional application of floating seismometers.

We launched an array of nine freely floating submarine seismometers near the Galápagos islands, which remained operational for about two years. P and PKP waves from regional and teleseismic earthquakes were observed for a range of magnitudes. The signal-to-noise ratio is strongly influenced by the weather conditions and this determines the lowest magnitudes that can be observed. Waves from deep earthquakes are easier to pick, but the S/N ratio can be enhanced through filtering and the data cover earthquakes from all depths. We measured 580 arrival times for different raypaths. We show that even such a limited number of data gives a significant increase in resolution for the oceanic upper mantle. This is the first time an array of floating seismometers is used in seismic tomography to improve the resolution significantly where otherwise no seismic information is available. We show that the Galápagos Archipelago is underlain by a deep (about 1900 km) 200-300 km wide plume of high temperature, with a heat flux very much larger than predicted from its swell bathymetry. The decrease of the plume temperature anomaly towards the surface indicates that the Earth's mantle has a subadiabatic temperature gradient.

Accelerating changes in ice mass within Greenland, and the ice sheet's sensitivity to atmospheric forcing.

From early 2003 to mid-2013, the total mass of ice in Greenland declined at a progressively increasing rate. In mid-2013, an abrupt reversal occurred, and very little net ice loss occurred in the next 12-18 months. Gravity Recovery and Climate Experiment (GRACE) and global positioning system (GPS) observations reveal that the spatial patterns of the sustained acceleration and the abrupt deceleration in mass loss are similar. The strongest accelerations tracked the phase of the North Atlantic Oscillation (NAO). The negative phase of the NAO enhances summertime warming and insolation while reducing snowfall, especially in west Greenland, driving surface mass balance (SMB) more negative, as illustrated using the regional climate model MAR. The spatial pattern of accelerating mass changes reflects the geography of NAO-driven shifts in atmospheric forcing and the ice sheet's sensitivity to that forcing. We infer that southwest Greenland will become a major future contributor to sea level rise.

Seismic monitoring in the oceans by autonomous floats.

Our understanding of the internal dynamics of the Earth is largely based on images of seismic velocity variations in the mantle obtained with global tomography. However, our ability to image the mantle is severely hampered by a lack of seismic data collected in marine areas. Here we report observations made under different noise conditions (in the Mediterranean Sea, the Indian and Pacific Oceans) by a submarine floating seismograph, and show that such floats are able to fill the oceanic data gap. Depending on the ambient noise level, the floats can record between 35 and 63% of distant earthquakes with a moment magnitude M≥6.5. Even magnitudes <6.0 can be successfully observed under favourable noise conditions. The serendipitous recording of an earthquake swarm near the Indian Ocean triple junction enabled us to establish a threshold magnitude between 2.7 and 3.4 for local earthquakes in the noisiest of the three environments.

Mapping Greenland's mass loss in space and time.

The melting of polar ice sheets is a major contributor to global sea-level rise. Early estimates of the mass lost from the Greenland ice cap, based on satellite gravity data collected by the Gravity Recovery and Climate Experiment, have widely varied. Although the continentally and decadally averaged estimated trends have now more or less converged, to this date, there has been little clarity on the detailed spatial distribution of Greenland's mass loss and how the geographical pattern has varied on relatively shorter time scales. Here, we present a spatially and temporally resolved estimation of the ice mass change over Greenland between April of 2002 and August of 2011. Although the total mass loss trend has remained linear, actively changing areas of mass loss were concentrated on the southeastern and northwestern coasts, with ice mass in the center of Greenland steadily increasing over the decade.

Probabilistic assessment of sea level during the last interglacial stage.

With polar temperatures approximately 3-5 degrees C warmer than today, the last interglacial stage (approximately 125 kyr ago) serves as a partial analogue for 1-2 degrees C global warming scenarios. Geological records from several sites indicate that local sea levels during the last interglacial were higher than today, but because local sea levels differ from global sea level, accurately reconstructing past global sea level requires an integrated analysis of globally distributed data sets. Here we present an extensive compilation of local sea level indicators and a statistical approach for estimating global sea level, local sea levels, ice sheet volumes and their associated uncertainties. We find a 95% probability that global sea level peaked at least 6.6 m higher than today during the last interglacial; it is likely (67% probability) to have exceeded 8.0 m but is unlikely (33% probability) to have exceeded 9.4 m. When global sea level was close to its current level (>or=-10 m), the millennial average rate of global sea level rise is very likely to have exceeded 5.6 m kyr(-1) but is unlikely to have exceeded 9.2 m kyr(-1). Our analysis extends previous last interglacial sea level studies by integrating literature observations within a probabilistic framework that accounts for the physics of sea level change. The results highlight the long-term vulnerability of ice sheets to even relatively low levels of sustained global warming.