Search

Estimating the future evolution of the Antarctic Ice Sheet (AIS) is critical for improving future sea level rise (SLR) projections. Numerical ice sheet models are invaluable tools for bounding Antarctic vulnerability; yet, few continental-scale projections of century-scale AIS SLR contribution exist, and those that do vary by up to an order of magnitude. This is partly because model projections of future sea level are inherently uncertain and depend largely on the model's boundary conditions and climate forcing, which themselves are unknown due to the uncertainty in the projections of future anthropogenic emissions and subsequent climate response. Here, we aim to improve the understanding of how uncertainties in model forcing and boundary conditions affect ice sheet model simulations. With use of sampling techniques embedded within the Ice Sheet System Model (ISSM) framework, we assess how uncertainties in snow accumulation, ocean-induced melting, ice viscosity, basal friction, bedrock elevation, and the presence of ice shelves impact continental-scale 100-year model simulations of AIS future sea level contribution. Overall, we find that AIS sea level contribution is strongly affected by grounding line retreat, which is driven by the magnitude of ice shelf basal melt rates and by variations in bedrock topography. In addition, we find that over 1.2 m of AIS global mean sea level contribution over the next century is achievable, but not likely, as it is tenable only in response to unrealistically large melt rates and continental ice shelf collapse. Regionally, we find that under our most extreme 100-year warming experiment generalized for the entire ice sheet, the Amundsen Sea sector is the most significant source of model uncertainty (1032 mm 6σ spread) and the region with the largest potential for future sea level contribution (297 mm). In contrast, under a more plausible forcing informed regionally by literature and model sensitivity studies, the Ronne basin has a greater potential for local increases in ice shelf basal melt rates. As a result, under this more likely realization, where warm waters reach the continental shelf under the Ronne ice shelf, it is the Ronne basin, particularly the Evans and Rutford ice streams, that are the greatest contributors to potential SLR (161 mm) and to simulation uncertainty (420 mm 6spread).

Greenland's bed topography is a primary control on ice flow, grounding line migration, calving dynamics, and subglacial drainage. Moreover, fjord bathymetry regulates the penetration of warm Atlantic water (AW) that rapidly melts and undercuts Greenland's marine-terminating glaciers. Here we present a new compilation of Greenland bed topography that assimilates seafloor bathymetry and ice thickness data through a mass conservation approach. A new 150 m horizontal resolution bed topography/bathymetric map of Greenland is constructed with seamless transitions at the ice/ocean interface, yielding major improvements over previous data sets, particularly in the marine-terminating sectors of northwest and southeast Greenland. Our map reveals that the total sea level potential of the Greenland ice sheet is 7.42 ± 0.05 m, which is 7 cm greater than previous estimates. Furthermore, it explains recent calving front response of numerous outlet glaciers and reveals new pathways by which AW can access glaciers with marine-based basins, thereby highlighting sectors of Greenland that are most vulnerable to future oceanic forcing.

Basal sliding is a main control on glacier flow primarily driven by water pressure at the glacier base. The ongoing increase in surface melting of the Greenland Ice Sheet warrants an examination of its impact on basal water pressure and in turn on basal sliding. Here we examine the case of Russell Glacier, in West Greenland, where an extensive set of observations has been collected. These observations suggest that the recent increase in melt has had an equivocal impact on the annual velocity, with stable flow on the lower part of the drainage basin but accelerated flow above the Equilibrium Line Altitude (ELA). These distinct behaviors have been attributed to different evolutions of the subglacial draining system during and after the melt season. Here we use a high-resolution subglacial hydrological model forced by reconstructed surface runoff for the period 2008 to 2012 to investigate the cause of these distinct behaviors. We find that the increase in meltwater production at low elevation yields a more efficient drainage system compatible with the observed stagnation of the mean annual flow below the ELA. At higher elevation, the model indicates that the drainage system is mostly inefficient and is therefore strongly sensitive to an increase in meltwater availability, which is consistent with the observed increase in ice velocity.

Surface melt and mass loss of the Greenland Ice Sheet may play crucial roles in global climate change due to their positive feedbacks and large fresh-water storage. With few other regular meteorological observations available in this extreme environment, measurements from automatic weather stations (AWS) are the primary data source for studying surface energy budgets, and for validating satellite observations and model simulations. Station tilt, due to irregular surface melt, compaction and glacier dynamics, causes considerable biases in the AWS shortwave radiation measurements. In this study, we identify tilt-induced biases in the climatology of surface shortwave radiative flux and albedo, and retrospectively correct these by iterative application of solar geometric principles. We found, over all the AWS from the Greenland Climate Network (GC-Net), the Kangerlussuaq transect (K-transect) and the Programme for Monitoring of the Greenland Ice Sheet (PROMICE) networks, insolation on fewer than 40% of clear days peaks within ±0.5h of solar noon time, with the largest shift exceeding 3h due to tilt. Hourly absolute biases in the magnitude of surface insolation can reach up to 200W m-2, with respect to the well-understood clear-day insolation. We estimate the tilt angles and their directions based on the solar geometric relationship between the simulated insolation at a horizontal surface and the observed insolation by these tilted AWS under clear-sky conditions. Our adjustment reduces the root mean square error (RMSE) against references from both satellite observation and reanalysis by 16W m-2 (24%), and raises the correlation coefficients with them to above 0.95. Averaged over the whole Greenland Ice Sheet in the melt season, the adjustment in insolation to compensate station tilt is ∼ 11W m-2, enough to melt 0.24m of snow water equivalent. The adjusted diurnal cycles of albedo are smoother, with consistent semi-smiling patterns. The seasonal cycles and inter-annual variabilities of albedo agree better with previous studies. This tilt-corrected shortwave radiation data set derived using the Retrospective, Iterative, Geometry-Based (RIGB) method provide more accurate observations and validations for surface energy budgets studies on the Greenland Ice Sheet, including albedo variations, surface melt simulations and cloud radiative forcing estimates.

We compare four independent estimates of the mass balance of the Amundsen Sea Embayment of West Antarctica, an area experiencing rapid retreat and mass loss to the sea. We use ICESat and Operation IceBridge laser altimetry, Envisat radar altimetry, GRACE time-variable gravity, RACMO2.3 surface mass balance, ice velocity from imaging radars, and ice thickness from radar sounders. The four methods agree in terms of mass loss and acceleration in loss at the regional scale. Over 1992-2013, the mass loss is 83 ± 5 Gt/yr with an acceleration of 6.1 ± 0.7 Gt/yr2. During the common period 2003-2009, the mass loss is 84 ± 10 Gt/yr with an acceleration of 16.3 ± 5.6 Gt/yr2, nearly 3 times the acceleration over 1992-2013. Over 2003-2011, the mass loss is 102 ± 10 Gt/yr with an acceleration of 15.7 ± 4.0 Gt/yr2. The results reconcile independent mass balance estimates in a setting dominated by change in ice dynamics with significant variability in surface mass balance.

We use Gravity Recovery and Climate Experiment (GRACE) monthly gravity fields to determine the regional acceleration in ice mass loss in Greenland and Antarctica for 2003-2013. We find that the total mass loss is controlled by only a few regions. In Greenland, the southeast and northwest generate 70% of the loss (280±58 Gt/yr) mostly from ice dynamics, the southwest accounts for 54% of the total acceleration in loss (25.4±1.2 Gt/yr2) from a decrease in surface mass balance (SMB), followed by the northwest (34%), and we find no significant acceleration in the northeast. In Antarctica, the Amundsen Sea (AS) sector and the Antarctic Peninsula account for 64% and 17%, respectively, of the total loss (180±10 Gt/yr) mainly from ice dynamics. The AS sector contributes most of the acceleration in loss (11±4 Gt/yr2), and Queen Maud Land, East Antarctica, is the only sector with a significant mass gain due to a local increase in SMB (63±5 Gt/yr).

This study evaluates the impact of a recent upgrade in the physics package of the regional atmospheric climate model RACMO2 on the simulated surface mass balance (SMB) of the Antarctic ice sheet. The modelled SMB increases, in particular over the grounded ice sheet of East Antarctica (+44Gt a-1), with a small change in West Antarctica. This mainly results from an increase in precipitation, which is explained by changes in the cloud microphysics, including a new parameterization for ice cloud supersaturation, and changes in large-scale circulation patterns, which alter topographically forced precipitation. The spatial changes in SMB are evaluated using 3234 in situ SMB observations and ice-balance velocities, and the temporal variability using GRACE satellite retrievals. The in situ observations and balance velocities show a clear improvement of the spatial representation of the SMB in the interior of East Antarctica, which has become considerably wetter. No improvements are seen for West Antarctica and the coastal regions. A comparison of model SMB temporal variability with GRACE satellite retrievals shows no significant change in performance.

Glacial isostatic adjustment (GIA) represents a source of uncertainty for ice sheet mass balance estimates from the Gravity Recovery and Climate Experiment (GRACE) time-variable gravity measurements. We evaluate Greenland GIA corrections from Simpson et al (2009 Quat. Sci. Rev. 28 1631-57), A et al (2013 Geophys. J. Int. 192 557-72) and Wu et al (2010 Nature Geosci. 3 642-6) by comparing the spatial patterns of GRACE-derived ice mass trends calculated using the three corrections with volume changes from ICESat (Ice, Cloud, and land Elevation Satellite) and OIB (Operation IceBridge) altimetry missions, and surface mass balance products from the Regional Atmospheric Climate Model (RACMO). During the period September 2003-August 2011, GRACE ice mass changes obtained using the Simpson et al (2009 Quat. Sci. Rev. 28 1631-57) and A et al (2013 Geophys. J. Int. 192 557-72) GIA corrections yield similar spatial patterns and amplitudes, and are consistent with altimetry observations and surface mass balance data. The two GRACE estimates agree within 2% on average over the entire ice sheet, and better than 15% in four subdivisions of Greenland. The third GRACE estimate corrected using the (Wu et al 2010 Nature Geosci. 3 642-6)) GIA shows similar spatial patterns, but produces an average ice mass loss for the entire ice sheet that is 64 - 67 Gt yr-1 smaller. In the Northeast the recovered ice mass change is 46-49 Gt yr -1 (245-270%) more positive than that deduced from the other two corrections. By comparing the spatial and temporal variability of the GRACE estimates with trends of volume changes from altimetry and surface mass balance from RACMO, we show that the Wu et al (2010 Nature Geosci. 3 642-6) correction leads to a large mass increase in the Northeast that is inconsistent with independent observations. © 2014 IOP Publishing Ltd.

The contribution to sea level to 2200 from the grounded, mainland Antarctic Peninsula ice sheet (APIS) was calculated using an ice-sheet model initialized with a new technique computing ice fluxes based on observed surface velocities, altimetry and surface mass balance, and computing volume response using a linearized method. Volume change estimates of the APIS resulting from surface massbalance anomalies calculated by the regional model RACMO2, forced by A1B and E1 scenarios of the global models ECHAM5 and HadCM3, predicted net negative sea-level contributions between -0.5 and -12mm sea-level equivalent (SLE) by 2200. Increased glacier flow due to ice thickening returned ~15% of the increased accumulation to the sea by 2100 and ~30% by 2200. The likely change in volume of the APIS by 2200 in response to imposed 10 and 20km retreats of the grounding line at individual large outlet glaciers in Palmer Land, southern Antarctic Peninsula, ranged between 0.5 and 3.5mm SLE per drainage basin. Ensemble calculations of APIS volume change resulting from imposed grounding-line retreat due to ice-shelf break-up scenarios applied to all 20 of the largest drainage basins in Palmer Land (covering ~40% of the total area of APIS) resulted in net sea-level contributions of 7-16mm SLE by 2100, and 10-25mm SLE by 2200. Inclusion of basins in the northern peninsula and realistic simulation of grounding-line movement for AP outlet glaciers will improve future projections.

Analysis of ICESat-1 data (2003-2008) shows significant surface lowering of Totten Glacier, the glacier discharging the largest volume of ice in East Antarctica, and less change on nearby Moscow University Glacier. After accounting for firn compaction anomalies, the thinning appears to coincide with fast-flowing ice indicating a dynamical origin. Here, to elucidate these observations, we apply high-resolution ice-ocean modelling. Totten Ice Shelf is simulated to have higher, more variable basal melting rates. We link this variability to the volume of cold water, originating in polynyas upon sea ice formation, reaching the sub-ice-shelf cavity. Hence, we propose that the observed increased thinning of Totten Glacier is due to enhanced basal melting caused by a decrease in cold polynya water reaching its cavity. We support this hypothesis with passive microwave data of polynya extent variability. Considering the widespread changes in sea ice conditions, this mechanism could be contributing extensively to ice-shelf instability.

Until quite recently, the mass balance (MB) of the great ice sheets of Greenland and Antarctica was poorly known and often treated as a residual in the budget of oceanic mass and sea level change. Recent developments in regional climate modelling and remote sensing, especially altimetry, gravimetry and InSAR feature tracking, have enabled us to specifically resolve the ice sheet mass balance components at a near-annual timescale. The results reveal significant mass losses for both ice sheets, caused by the acceleration of marine-terminating glaciers in southeast, west and northwest Greenland and coastal West Antarctica, and increased run-off in Greenland. At the same time, the data show that interannual variability is very significant, masking the underlying trends. © 2011 The Author(s).

High-resolution (∼11 km) regional climate modeling shows total annual precipitation on the Greenland ice sheet for 1958-2007 to be up to 24% and surface mass balance up to 63% higher than previously thought. The largest differences occur in coastal southeast Greenland, where the much higher resolution facilitates capturing snow accumulation peaks that past five-fold coarser resolution regional climate models missed. The surface mass balance trend over the full 1958-2007 period reveals the classic pattern expected in a warming climate, with increased snowfall in the interior and enhanced runoff from the marginal ablation zone. In the period 1990-2007, total runoff increased significantly, 3% per year. The absolute increase in runoff is especially pronounced in the southeast, where several outlet glaciers have recently accelerated. This detailed knowledge of Greenland's surface mass balance provides the foundation for estimating and predicting the overall mass balance and freshwater discharge of the ice sheet. Copyright 2009 by the American Geophysical Union.

Large uncertainties remain in the current and future contribution to sea level rise from Antarctica. Climate warming may increase snowfall in the continents interior, but enhance glacier discharge at the coast where warmer air and ocean temperatures erode the buttressing ice shelves. Here, we use satellite interferometric synthetic-aperture radar observations from 1992 to 2006 covering 85 of Antarcticas coastline to estimate the total mass flux into the ocean. We compare the mass fluxes from large drainage basin units with interior snow accumulation calculated from a regional atmospheric climate model for 1980 to 2004. In East Antarctica, small glacier losses in Wilkes Land and glacier gains at the mouths of the Filchner and Ross ice shelves combine to a near-zero loss of 4±61 Gt yr1. In West Antarctica, widespread losses along the Bellingshausen and Amundsen seas increased the ice sheet loss by 59 in 10 years to reach 132±60 Gt yr1 in 2006. In the Peninsula, losses increased by 140 to reach 60±46 Gt yr1 in 2006. Losses are concentrated along narrow channels occupied by outlet glaciers and are caused by ongoing and past glacier acceleration. Changes in glacier flow therefore have a significant, if not dominant impact on ice sheet mass balance. © 2008 Nature Publishing Group.

The Greenland ice sheet (GrIS) is at present the largest single contributor to global-mass-induced sea-level rise, primarily because of Arctic amplification on an increasingly warmer Earth 1-5 . However, the processes of englacial water accumulation, storage and ultimate release remain poorly constrained. Here we show that a noticeable amount of the summertime meltwater mass is temporally buffered along the entire GrIS periphery, peaking in July and gradually reducing thereafter. Our results arise from quantifying the spatiotemporal behaviour of the total mass of water leaving the GrIS by analysing bedrock elastic deformation measured by Global Navigation Satellite System (GNSS) stations. The buffered meltwater causes a subsidence of the bedrock close to GNSS stations of at most approximately 5 mm during the melt season. Regionally, the duration of meltwater storage ranges from 4.5 weeks in the southeast to 9 weeks elsewhere. We also show that the meltwater runoff modelled from regional climate models may contain systematic errors, requiring further scaling of up to about 20% for the warmest years. These results reveal a high potential for GNSS data to constrain poorly known hydrological processes in Greenland, forming the basis for improved projections of future GrIS melt behaviour and the associated sea-level rise 6 ., (© 2024. The Author(s).)

Antarctic ice sheet (AIS) mass loss is predominantly driven by increased solid ice discharge, but its variability is governed by surface processes. Snowfall fluctuations control the surface mass balance (SMB) of the grounded AIS, while meltwater ponding can trigger ice shelf collapse potentially accelerating discharge. Surface processes are essential to quantify AIS mass change, but remain poorly represented in climate models typically running at 25-100 km resolution. Here we present SMB and surface melt products statistically downscaled to 2 km resolution for the contemporary climate (1979-2021) and low, moderate and high-end warming scenarios until 2100. We show that statistical downscaling modestly enhances contemporary SMB (3%), which is sufficient to reconcile modelled and satellite mass change. Furthermore, melt strongly increases (46%), notably near the grounding line, in better agreement with in-situ and satellite records. The melt increase persists by 2100 in all warming scenarios, revealing higher surface melt rates than previously estimated., (© 2023. The Author(s).)

Mass loss from the West Antarctic Ice Sheet is dominated by glaciers draining into the Amundsen Sea Embayment (ASE), yet the impact of anomalous precipitation on the mass balance of the ASE is poorly known. Here we present a 25-year (1996-2021) record of ASE input-output mass balance and evaluate how two periods of anomalous precipitation affected its sea level contribution. Since 1996, the ASE has lost 3331 ± 424 Gt ice, contributing 9.2 ± 1.2 mm to global sea level. Overall, surface mass balance anomalies contributed little (7.7%) to total mass loss; however, two anomalous precipitation events had larger, albeit short-lived, impacts on rates of mass change. During 2009-2013, persistently low snowfall led to an additional 51 ± 4 Gt yr -1 mass loss in those years (contributing positively to the total loss of 195 ± 4 Gt yr -1 ). Contrastingly, extreme precipitation in the winters of 2019 and 2020 decreased mass loss by 60 ± 16 Gt yr -1 during those years (contributing negatively to the total loss of 107 ± 15 Gt yr -1 ). These results emphasise the important impact of extreme snowfall variability on the short-term sea level contribution from West Antarctica., (© 2023. The Author(s).)

Firn (compressed snow) covers approximately 90[Formula: see text] of the Greenland ice sheet (GrIS) and currently retains about half of rain and meltwater through refreezing, reducing runoff and subsequent mass loss. The loss of firn could mark a tipping point for sustained GrIS mass loss, since decades to centuries of cold summers would be required to rebuild the firn buffer. Here we estimate the warming required for GrIS firn to reach peak refreezing, using 51 climate simulations statistically downscaled to 1 km resolution, that project the long-term firn layer evolution under multiple emission scenarios (1850-2300). We predict that refreezing stabilises under low warming scenarios, whereas under extreme warming, refreezing could peak and permanently decline starting in southwest Greenland by 2100, and further expanding GrIS-wide in the early 22[Formula: see text] century. After passing this peak, the GrIS contribution to global sea level rise would increase over twenty-fold compared to the last three decades., (© 2022. The Author(s).)

In recent decades, Greenland's peripheral glaciers have experienced large-scale mass loss, resulting in a substantial contribution to sea level rise. While their total area of Greenland ice cover is relatively small (4%), their mass loss is disproportionally large compared to the Greenland ice sheet. Satellite altimetry from Ice, Cloud, and land Elevation Satellite (ICESat) and ICESat-2 shows that mass loss from Greenland's peripheral glaciers increased from 27.2 ± 6.2 Gt/yr (February 2003-October 2009) to 42.3 ± 6.2 Gt/yr (October 2018-December 2021). These relatively small glaciers now constitute 11 ± 2% of Greenland's ice loss and contribute to global sea level rise. In the period October 2018-December 2021, mass loss increased by a factor of four for peripheral glaciers in North Greenland. While peripheral glacier mass loss is widespread, we also observe a complex regional pattern where increases in precipitation at high altitudes have partially counteracted increases in melt at low altitude., (© 2022. The Authors.)

Compared to other Arctic ice masses, Svalbard glaciers are low-elevated with flat interior accumulation areas, resulting in a marked peak in their current hypsometry (area-elevation distribution) at  ~450 m above sea level. Since summer melt consistently exceeds winter snowfall, these low-lying glaciers can only survive by refreezing a considerable fraction of surface melt and rain in the porous firn layer covering their accumulation zones. We use a high-resolution climate model to show that modest atmospheric warming in the mid-1980s forced the firn zone to retreat upward by  ~100 m to coincide with the hypsometry peak. This led to a rapid areal reduction of firn cover available for refreezing, and strongly increased runoff from dark, bare ice areas, amplifying mass loss from all elevations. As the firn line fluctuates around the hypsometry peak in the current climate, Svalbard glaciers will continue to lose mass and show high sensitivity to temperature perturbations.

Since the early 1990s, the Greenland ice sheet (GrIS) has been losing mass at an accelerating rate, primarily due to enhanced meltwater runoff following atmospheric warming. Here, we show that a pronounced latitudinal contrast exists in the GrIS response to recent warming. The ablation area in north Greenland expanded by 46%, almost twice as much as in the south (+25%), significantly increasing the relative contribution of the north to total GrIS mass loss. This latitudinal contrast originates from a different response to the recent change in large-scale Arctic summertime atmospheric circulation, promoting southwesterly advection of warm air toward the GrIS. In the southwest, persistent high atmospheric pressure reduced cloudiness, increasing runoff through enhanced absorption of solar radiation; in contrast, increased early-summer cloudiness in north Greenland enhanced atmospheric warming through decreased longwave heat loss. This triggered a rapid snowline retreat, causing early bare ice exposure, amplifying northern runoff.

Surface mass balance (SMB) provides mass input to the surface of the Antarctic and Greenland Ice Sheets and therefore comprises an important control on ice sheet mass balance and resulting contribution to global sea level change. As ice sheet SMB varies highly across multiple scales of space (meters to hundreds of kilometers) and time (hourly to decadal), it is notoriously challenging to observe and represent in models. In addition, SMB consists of multiple components, all of which depend on complex interactions between the atmosphere and the snow/ice surface, large-scale atmospheric circulation and ocean conditions, and ice sheet topography. In this review, we present the state-of-the-art knowledge and recent advances in ice sheet SMB observations and models, highlight current shortcomings, and propose future directions. Novel observational methods allow mapping SMB across larger areas, longer time periods, and/or at very high (subdaily) temporal frequency. As a recent observational breakthrough, cosmic ray counters provide direct estimates of SMB, circumventing the need for accurate snow density observations upon which many other techniques rely. Regional atmospheric climate models have drastically improved their simulation of ice sheet SMB in the last decade, thanks to the inclusion or improved representation of essential processes (e.g., clouds, blowing snow, and snow albedo), and by enhancing horizontal resolution (5-30 km). Future modeling efforts are required in improving Earth system models to match regional atmospheric climate model performance in simulating ice sheet SMB, and in reinforcing the efforts in developing statistical and dynamic downscaling to represent smaller-scale SMB processes., (©2019. The Authors.)

The Canadian Arctic Archipelago contains >300 glaciers that terminate in the ocean, but little is known about changes in their frontal positions in response to recent changes in the ocean-climate system. Here, we examine changes in glacier frontal positions since the 1950s and investigate the relative influence of oceanic temperature versus atmospheric temperature. Over 94% of glaciers retreated between 1958 and 2015, with a region-wide trend of gradual retreat before ~2000, followed by a fivefold increase in retreat rates up to 2015. Retreat patterns show no correlation with changes in subsurface ocean temperatures, in clear contrast to the dominance of ocean forcing in western Greenland and elsewhere. Rather, significant correlations with surface melt indicate that increased atmospheric temperature has been the primary driver of the acceleration in marine-terminating glacier frontal retreat in this region.

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., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

The Greenland ice sheet (GrIS) is a growing contributor to global sea-level rise 1 , with recent ice mass loss dominated by surface meltwater runoff 2,3 . Satellite observations reveal positive trends in GrIS surface melt extent 4 , but melt variability, intensity and runoff remain uncertain before the satellite era. Here we present the first continuous, multi-century and observationally constrained record of GrIS surface melt intensity and runoff, revealing that the magnitude of recent GrIS melting is exceptional over at least the last 350 years. We develop this record through stratigraphic analysis of central west Greenland ice cores, and demonstrate that measurements of refrozen melt layers in percolation zone ice cores can be used to quantifiably, and reproducibly, reconstruct past melt rates. We show significant (P < 0.01) and spatially extensive correlations between these ice-core-derived melt records and modelled melt rates 5,6 and satellite-derived melt duration 4 across Greenland more broadly, enabling the reconstruction of past ice-sheet-scale surface melt intensity and runoff. We find that the initiation of increases in GrIS melting closely follow the onset of industrial-era Arctic warming in the mid-1800s, but that the magnitude of GrIS melting has only recently emerged beyond the range of natural variability. Owing to a nonlinear response of surface melting to increasing summer air temperatures, continued atmospheric warming will lead to rapid increases in GrIS runoff and sea-level contributions.

The freshwater budget of the Arctic and sub-polar North Atlantic Oceans has been changing due, primarily, to increased river runoff, declining sea ice and enhanced melting of Arctic land ice. Since the mid-1990s this latter component has experienced a pronounced increase. We use a combination of satellite observations of glacier flow speed and regional climate modeling to reconstruct the land ice freshwater flux from the Greenland ice sheet and Arctic glaciers and ice caps for the period 1958-2016. The cumulative freshwater flux anomaly exceeded 6,300 ± 316 km 3 by 2016. This is roughly twice the estimate of a previous analysis that did not include glaciers and ice caps outside of Greenland and which extended only to 2010. From 2010 onward, the total freshwater flux is about 1,300 km 3 /yr, equivalent to 0.04 Sv, which is roughly 40% of the estimated total runoff to the Arctic for the same time period. Not all of this flux will reach areas of deep convection or Arctic and Sub-Arctic seas. We note, however, that the largest freshwater flux anomalies, grouped by ocean basin, are located in Baffin Bay and Davis Strait. The land ice freshwater flux displays a strong seasonal cycle with summer time values typically around five times larger than the annual mean. This will be important for understanding the impact of these fluxes on fjord circulation, stratification, and the biogeochemistry of, and nutrient delivery to, coastal waters.

Rapid changes in thickness and velocity have been observed at many marine-terminating glaciers in Greenland, impacting the volume of ice they export, or discharge, from the ice sheet. While annual estimates of ice-sheet wide discharge have been previously derived, higher-resolution records are required to fully constrain the temporal response of these glaciers to various climatic and mechanical drivers that vary in sub-annual scales. Here we sample outlet glaciers wider than 1 km (N = 230) to derive the first continuous, ice-sheet wide record of total ice sheet discharge for the 2000-2016 period, resolving a seasonal variability of 6 %. The amplitude of seasonality varies spatially across the ice sheet from 5 % in the southeastern region to 9 % in the northwest region. We analyze seasonal to annual variability in the discharge time series with respect to both modelled meltwater runoff, obtained from RACMO2.3p2, and glacier front position changes over the same period. We find that year-to-year changes in total ice sheet discharge are related to annual front changes ( r 2 = 0.59, p = 10 -4 ) and that the annual magnitude of discharge is closely related to cumulative front position changes ( r 2 = 0.79), which show a net retreat of > 400 km, or an average retreat of > 2 km at each surveyed glacier. Neither maximum seasonal runoff or annual runoff totals are correlated to annual discharge, which suggests that larger annual quantities of runoff do not relate to increased annual discharge. Discharge and runoff, however, follow similar patterns of seasonal variability with near-coincident periods of acceleration and seasonal maxima. These results suggest that changes in glacier front position drive secular trends in discharge, whereas the impact of runoff is likely limited to the summer months when observed seasonal variations are substantially controlled by the timing of meltwater input., Competing Interests: Competing interests The authors declare that they have no conflict of interest.

Meltwater runoff from the Greenland ice sheet surface influences surface mass balance (SMB), ice dynamics, and global sea level rise, but is estimated with climate models and thus difficult to validate. We present a way to measure ice surface runoff directly, from hourly in situ supraglacial river discharge measurements and simultaneous high-resolution satellite/drone remote sensing of upstream fluvial catchment area. A first 72-h trial for a 63.1-km 2 moulin-terminating internally drained catchment (IDC) on Greenland's midelevation (1,207-1,381 m above sea level) ablation zone is compared with melt and runoff simulations from HIRHAM5, MAR3.6, RACMO2.3, MERRA-2, and SEB climate/SMB models. Current models cannot reproduce peak discharges or timing of runoff entering moulins but are improved using synthetic unit hydrograph (SUH) theory. Retroactive SUH applications to two older field studies reproduce their findings, signifying that remotely sensed IDC area, shape, and supraglacial river length are useful for predicting delays in peak runoff delivery to moulins. Applying SUH to HIRHAM5, MAR3.6, and RACMO2.3 gridded melt products for 799 surrounding IDCs suggests their terminal moulins receive lower peak discharges, less diurnal variability, and asynchronous runoff timing relative to climate/SMB model output alone. Conversely, large IDCs produce high moulin discharges, even at high elevations where melt rates are low. During this particular field experiment, models overestimated runoff by +21 to +58%, linked to overestimated surface ablation and possible meltwater retention in bare, porous, low-density ice. Direct measurements of ice surface runoff will improve climate/SMB models, and incorporating remotely sensed IDCs will aid coupling of SMB with ice dynamics and subglacial systems., Competing Interests: The authors declare no conflict of interest., (Copyright © 2017 the Author(s). Published by PNAS.)

We provide the first comprehensive analysis of the relationships between large-scale patterns of Southern Hemisphere climate variability and the detailed structure of Antarctic precipitation. We examine linkages between the high spatial resolution precipitation from a regional atmospheric model and four patterns of large-scale Southern Hemisphere climate variability: the southern baroclinic annular mode, the southern annular mode, and the two Pacific-South American teleconnection patterns. Variations in all four patterns influence the spatial configuration of precipitation over Antarctica, consistent with their signatures in high-latitude meridional moisture fluxes. They impact not only the mean but also the incidence of extreme precipitation events. Current coupled-climate models are able to reproduce all four patterns of atmospheric variability but struggle to correctly replicate their regional impacts on Antarctic climate. Thus, linking these patterns directly to Antarctic precipitation variability may allow a better estimate of future changes in precipitation than using model output alone.

Melting of the Greenland ice sheet (GrIS) and its peripheral glaciers and ice caps (GICs) contributes about 43% to contemporary sea level rise. While patterns of GrIS mass loss are well studied, the spatial and temporal evolution of GICs mass loss and the acting processes have remained unclear. Here we use a novel, 1 km surface mass balance product, evaluated against in situ and remote sensing data, to identify 1997 (±5 years) as a tipping point for GICs mass balance. That year marks the onset of a rapid deterioration in the capacity of the GICs firn to refreeze meltwater. Consequently, GICs runoff increases 65% faster than meltwater production, tripling the post-1997 mass loss to 36±16 Gt -1 , or ∼14% of the Greenland total. In sharp contrast, the extensive inland firn of the GrIS retains most of its refreezing capacity for now, buffering 22% of the increased meltwater production. This underlines the very different response of the GICs and GrIS to atmospheric warming.

We propose a new ice sheet model validation framework - the Cryospheric Model Comparison Tool (CmCt) - that takes advantage of ice sheet altimetry and gravimetry observations collected over the past several decades and is applied here to modeling of the Greenland ice sheet. We use realistic simulations performed with the Community Ice Sheet Model (CISM) along with two idealized, non-dynamic models to demonstrate the framework and its use. Dynamic simulations with CISM are forced from 1991 to 2013 using combinations of reanalysis-based surface mass balance and observations of outlet glacier flux change. We propose and demonstrate qualitative and quantitative metrics for use in evaluating the different model simulations against the observations. We find that the altimetry observations used here are largely ambiguous in terms of their ability to distinguish one simulation from another. Based on basin- and whole-ice-sheet scale metrics, we find that simulations using both idealized conceptual models and dynamic, numerical models provide an equally reasonable representation of the ice sheet surface (mean elevation differences of <1 m). This is likely due to their short period of record, biases inherent to digital elevation models used for model initial conditions, and biases resulting from firn dynamics, which are not explicitly accounted for in the models or observations. On the other hand, we find that the gravimetry observations used here are able to unambiguously distinguish between simulations of varying complexity, and along with the CmCt, can provide a quantitative score for assessing a particular model and/or simulation. The new framework demonstrates that our proposed metrics can distinguish relatively better from relatively worse simulations and that dynamic ice sheet models, when appropriately initialized and forced with the right boundary conditions, demonstrate predictive skill with respect to observed dynamic changes occurring on Greenland over the past few decades. An extensible design will allow for continued use of the CmCt as future altimetry, gravimetry, and other remotely sensed data become available for use in ice sheet model validation.

The Atlantic Meridional Overturning Circulation (AMOC) is an important component of ocean thermohaline circulation. Melting of Greenland's ice sheet is freshening the North Atlantic; however, whether the augmented freshwater flux is disrupting the AMOC is unclear. Dense Labrador Sea Water (LSW), formed by winter cooling of saline North Atlantic water and subsequent convection, is a key component of the deep southward return flow of the AMOC. Although LSW formation recently decreased, it also reached historically high values in the mid-1990s, making the connection to the freshwater flux unclear. Here we derive a new estimate of the recent freshwater flux from Greenland using updated GRACE satellite data, present new flux estimates for heat and salt from the North Atlantic into the Labrador Sea and explain recent variations in LSW formation. We suggest that changes in LSW can be directly linked to recent freshening, and suggest a possible link to AMOC weakening.

The Greenland ice sheet has become one of the main contributors to global sea level rise, predominantly through increased meltwater runoff. The main drivers of Greenland ice sheet runoff, however, remain poorly understood. Here we show that clouds enhance meltwater runoff by about one-third relative to clear skies, using a unique combination of active satellite observations, climate model data and snow model simulations. This impact results from a cloud radiative effect of 29.5 (±5.2) W m(-2). Contrary to conventional wisdom, however, the Greenland ice sheet responds to this energy through a new pathway by which clouds reduce meltwater refreezing as opposed to increasing surface melt directly, thereby accelerating bare-ice exposure and enhancing meltwater runoff. The high sensitivity of the Greenland ice sheet to both ice-only and liquid-bearing clouds highlights the need for accurate cloud representations in climate models, to better predict future contributions of the Greenland ice sheet to global sea level rise.

Growing evidence has demonstrated the importance of ice shelf buttressing on the inland grounded ice, especially if it is resting on bedrock below sea level. Much of the Southern Antarctic Peninsula satisfies this condition and also possesses a bed slope that deepens inland. Such ice sheet geometry is potentially unstable. We use satellite altimetry and gravity observations to show that a major portion of the region has, since 2009, destabilized. Ice mass loss of the marine-terminating glaciers has rapidly accelerated from close to balance in the 2000s to a sustained rate of -56 ± 8 gigatons per year, constituting a major fraction of Antarctica's contribution to rising sea level. The widespread, simultaneous nature of the acceleration, in the absence of a persistent atmospheric forcing, points to an oceanic driving mechanism., (Copyright © 2015, American Association for the Advancement of Science.)

We present a new record of ice thickness change, reconstructed at nearly 100,000 sites on the Greenland Ice Sheet (GrIS) from laser altimetry measurements spanning the period 1993-2012, partitioned into changes due to surface mass balance (SMB) and ice dynamics. We estimate a mean annual GrIS mass loss of 243 ± 18 Gt ⋅ y(-1), equivalent to 0.68 mm ⋅ y(-1) sea level rise (SLR) for 2003-2009. Dynamic thinning contributed 48%, with the largest rates occurring in 2004-2006, followed by a gradual decrease balanced by accelerating SMB loss. The spatial pattern of dynamic mass loss changed over this time as dynamic thinning rapidly decreased in southeast Greenland but slowly increased in the southwest, north, and northeast regions. Most outlet glaciers have been thinning during the last two decades, interrupted by episodes of decreasing thinning or even thickening. Dynamics of the major outlet glaciers dominated the mass loss from larger drainage basins, and simultaneous changes over distances up to 500 km are detected, indicating climate control. However, the intricate spatiotemporal pattern of dynamic thickness change suggests that, regardless of the forcing responsible for initial glacier acceleration and thinning, the response of individual glaciers is modulated by local conditions. Recent projections of dynamic contributions from the entire GrIS to SLR have been based on the extrapolation of four major outlet glaciers. Considering the observed complexity, we question how well these four glaciers represent all of Greenland's outlet glaciers.

Predicting Greenland Ice Sheet mass loss due to ice dynamics requires a complete understanding of spatiotemporal velocity fluctuations and related control mechanisms. We present a 5 year record of seasonal velocity measurements for 55 marine-terminating glaciers distributed around the ice sheet margin, along with ice-front position and runoff data sets for each glacier. Among glaciers with substantial speed variations, we find three distinct seasonal velocity patterns. One pattern indicates relatively high glacier sensitivity to ice-front position. The other two patterns are more prevalent and appear to be meltwater controlled. These patterns reveal differences in which some subglacial systems likely transition seasonally from inefficient, distributed hydrologic networks to efficient, channelized drainage, while others do not. The difference may be determined by meltwater availability, which in some regions may be influenced by perennial firn aquifers. Our results highlight the need to understand subglacial meltwater availability on an ice sheet-wide scale to predict future dynamic changes., Key Points: First multi-region seasonal velocity measurements show regional differencesSeasonal velocity fluctuations on most glaciers appear meltwater controlledSeasonal development of efficient subglacial drainage geographically divided.