Transient ice loss in the Patagonia Icefields during the 2015-2016 El Niño event.

The Patagonia Icefields (PIF) are the largest non-polar ice mass in the southern hemisphere. The icefields cover an area of approximately 16,500 km2 and are divided into the northern and southern icefields, which are ~ 4000 km2 and ~ 12,500 km2, respectively. While both icefields have been losing mass rapidly, their responsiveness to various climate drivers, such as the El Niño-Southern Oscillation, is not well understood. Using the elastic response of the earth to loading changes and continuous GPS data we separated and estimated ice mass changes observed during the strong El Niño that started in 2015 from the complex hydrological interactions occurring around the PIF. During this single event, our mass balance estimates show that the northern icefield lost ~ 28 Gt of mass while the southern icefield lost ~ 12 Gt. This is the largest ice loss event in the PIF observed to date using geodetic data.

Centennial response of Greenland's three largest outlet glaciers.

The Greenland Ice Sheet is the largest land ice contributor to sea level rise. This will continue in the future but at an uncertain rate and observational estimates are limited to the last few decades. Understanding the long-term glacier response to external forcing is key to improving projections. Here we use historical photographs to calculate ice loss from 1880-2012 for Jakobshavn, Helheim, and Kangerlussuaq glacier. We estimate ice loss corresponding to a sea level rise of 8.1 ± 1.1 millimetres from these three glaciers. Projections of mass loss for these glaciers, using the worst-case scenario, Representative Concentration Pathways 8.5, suggest a sea level contribution of 9.1-14.9 mm by 2100. RCP8.5 implies an additional global temperature increase of 3.7 °C by 2100, approximately four times larger than that which has taken place since 1880. We infer that projections forced by RCP8.5 underestimate glacier mass loss which could exceed this worst-case scenario.

Months-long thousand-kilometre-scale wobbling before great subduction earthquakes.

Megathrust earthquakes are responsible for some of the most devastating natural disasters1. To better understand the physical mechanisms of earthquake generation, subduction zones worldwide are continuously monitored with geophysical instrumentation. One key strategy is to install stations that record signals from Global Navigation Satellite Systems2,3 (GNSS), enabling us to track the non-steady surface motion of the subducting and overriding plates before, during and after the largest events4-6. Here we use a recently developed trajectory modelling approach7 that is designed to isolate secular tectonic motions from the daily GNSS time series to show that the 2010 Maule, Chile (moment magnitude 8.8) and 2011 Tohoku-oki, Japan (moment magnitude 9.0) earthquakes were preceded by reversals of 4-8 millimetres in surface displacement that lasted several months and spanned thousands of kilometres. Modelling of the surface displacement reversal that occurred before the Tohoku-oki earthquake suggests an initial slow slip followed by a sudden pulldown of the Philippine Sea slab so rapid that it caused a viscoelastic rebound across the whole of Japan. Therefore, to understand better when large earthquakes are imminent, we must consider not only the evolution of plate interface frictional processes but also the dynamic boundary conditions from deeper subduction processes, such as sudden densification of metastable slab.

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.

Illuminating subduction zone rheological properties in the wake of a giant earthquake.

Deformation associated with plate convergence at subduction zones is accommodated by a complex system involving fault slip and viscoelastic flow. These processes have proven difficult to disentangle. The 2010 M w 8.8 Maule earthquake occurred close to the Chilean coast within a dense network of continuously recording Global Positioning System stations, which provide a comprehensive history of surface strain. We use these data to assemble a detailed picture of a structurally controlled megathrust fault frictional patchwork and the three-dimensional rheological and time-dependent viscosity structure of the lower crust and upper mantle, all of which control the relative importance of afterslip and viscoelastic relaxation during postseismic deformation. These results enhance our understanding of subduction dynamics including the interplay of localized and distributed deformation during the subduction zone earthquake cycle.

Observed rapid bedrock uplift in Amundsen Sea Embayment promotes ice-sheet stability.

The marine portion of the West Antarctic Ice Sheet (WAIS) in the Amundsen Sea Embayment (ASE) accounts for one-fourth of the cryospheric contribution to global sea-level rise and is vulnerable to catastrophic collapse. The bedrock response to ice mass loss, glacial isostatic adjustment (GIA), was thought to occur on a time scale of 10,000 years. We used new GPS measurements, which show a rapid (41 millimeters per year) uplift of the ASE, to estimate the viscosity of the mantle underneath. We found a much lower viscosity (4 × 1018 pascal-second) than global average, and this shortens the GIA response time scale from tens to hundreds of years. Our finding requires an upward revision of ice mass loss from gravity data of 10% and increases the potential stability of the WAIS against catastrophic collapse.

Geodetic measurements reveal similarities between post-Last Glacial Maximum and present-day mass loss from the Greenland ice sheet.

Accurate quantification of the millennial-scale mass balance of the Greenland ice sheet (GrIS) and its contribution to global sea-level rise remain challenging because of sparse in situ observations in key regions. Glacial isostatic adjustment (GIA) is the ongoing response of the solid Earth to ice and ocean load changes occurring since the Last Glacial Maximum (LGM; ~21 thousand years ago) and may be used to constrain the GrIS deglaciation history. We use data from the Greenland Global Positioning System network to directly measure GIA and estimate basin-wide mass changes since the LGM. Unpredicted, large GIA uplift rates of +12 mm/year are found in southeast Greenland. These rates are due to low upper mantle viscosity in the region, from when Greenland passed over the Iceland hot spot about 40 million years ago. This region of concentrated soft rheology has a profound influence on reconstructing the deglaciation history of Greenland. We reevaluate the evolution of the GrIS since LGM and obtain a loss of 1.5-m sea-level equivalent from the northwest and southeast. These same sectors are dominating modern mass loss. We suggest that the present destabilization of these marine-based sectors may increase sea level for centuries to come. Our new deglaciation history and GIA uplift estimates suggest that studies that use the Gravity Recovery and Climate Experiment satellite mission to infer present-day changes in the GrIS may have erroneously corrected for GIA and underestimated the mass loss by about 20 gigatons/year.

Recharge of a subglacial lake by surface meltwater in northeast Greenland.

In a warming climate, surface meltwater production on large ice sheets is expected to increase. If this water is delivered to the ice sheet base it may have important consequences for ice dynamics. For example, basal water distributed in a diffuse network can decrease basal friction and accelerate ice flow, whereas channelized basal water can move quickly to the ice margin, where it can alter fjord circulation and submarine melt rates. Less certain is whether surface meltwater can be trapped and stored in subglacial lakes beneath large ice sheets. Here we show that a subglacial lake in Greenland drained quickly, as seen in the collapse of the ice surface, and then refilled from surface meltwater input. We use digital elevation models from stereo satellite imagery and airborne measurements to resolve elevation changes during the evolution of the surface and basal hydrologic systems at the Flade Isblink ice cap in northeast Greenland. During the autumn of 2011, a collapse basin about 70 metres deep and about 0.4 cubic kilometres in volume formed near the southern summit of the ice cap as a subglacial lake drained into a nearby fjord. Over the next two years, rapid uplift of the floor of the basin (which is approximately 8.4 square kilometres in area) occurred as surface meltwater flowed into crevasses around the basin margin and refilled the subglacial lake. Our observations show that surface meltwater can be trapped and stored at the bed of an ice sheet. Sensible and latent heat released by this trapped meltwater could soften nearby colder basal ice and alter downstream ice dynamics. Heat transport associated with meltwater trapped in subglacial lakes should be considered when predicting how ice sheet behaviour will change in a warming climate.

Bedrock displacements in Greenland manifest ice mass variations, climate cycles and climate change.

The Greenland GPS Network (GNET) uses the Global Positioning System (GPS) to measure the displacement of bedrock exposed near the margins of the Greenland ice sheet. The entire network is uplifting in response to past and present-day changes in ice mass. Crustal displacement is largely accounted for by an annual oscillation superimposed on a sustained trend. The oscillation is driven by earth's elastic response to seasonal variations in ice mass and air mass (i.e., atmospheric pressure). Observed vertical velocities are higher and often much higher than predicted rates of postglacial rebound (PGR), implying that uplift is usually dominated by the solid earth's instantaneous elastic response to contemporary losses in ice mass rather than PGR. Superimposed on longer-term trends, an anomalous 'pulse' of uplift accumulated at many GNET stations during an approximate six-month period in 2010. This anomalous uplift is spatially correlated with the 2010 melting day anomaly.

Crustal dilatation observed by GRACE after the 2004 Sumatra-Andaman earthquake.

We report the detection of an earthquake by a space-based measurement. The Gravity Recovery and Climate Experiment (GRACE) satellites observed a +/-15-microgalileo gravity change induced by the great December 2004 Sumatra-Andaman earthquake. Coseismic deformation produces sudden changes in the gravity field by vertical displacement of Earth's layered density structure and by changing the densities of the crust and mantle. GRACE's sensitivity to the long spatial wavelength of gravity changes resulted in roughly equal contributions of vertical displacement and dilatation effects in the gravity measurements. The GRACE observations provide evidence of crustal dilatation resulting from an undersea earthquake.