Higher Antarctic ice sheet accumulation and surface melt rates revealed at 2 km resolution.

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.

Increasing extreme melt in northeast Greenland linked to foehn winds and atmospheric rivers.

The Greenland Ice Sheet has been losing mass at an increased rate in recent decades. In northeast Greenland, increasing surface melt has accompanied speed-ups in the outlet glaciers of the Northeast Greenland Ice Stream, which contain over one meter of sea level rise potential. Here we show that the most intense northeast Greenland melt events are driven by atmospheric rivers (ARs) affecting northwest Greenland that induce foehn winds in the northeast. Near low-elevation outlet glaciers, 80-100% of extreme (> 99th percentile) melt occurs during foehn conditions and 50-75% during ARs. These events have become more frequent during the twenty-first century, with 5-10% of total northeast Greenland melt in several recent summers occurring during the ~1% of times with strong AR and foehn conditions. We conclude that the combined AR-foehn influence on northeast Greenland extreme melt will likely continue to grow as regional atmospheric moisture content increases with climate warming.

Peak refreezing in the Greenland firn layer under future warming scenarios.

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.

Greenland Mass Trends From Airborne and Satellite Altimetry During 2011-2020.

We use satellite and airborne altimetry to estimate annual mass changes of the Greenland Ice Sheet. We estimate ice loss corresponding to a sea-level rise of 6.9 ± 0.4 mm from April 2011 to April 2020, with a highest annual ice loss rate of 1.4 mm/yr sea-level equivalent from April 2019 to April 2020. On a regional scale, our annual mass loss timeseries reveals 10-15 m/yr dynamic thickening at the terminus of Jakobshavn Isbræ from April 2016 to April 2018, followed by a return to dynamic thinning. We observe contrasting patterns of mass loss acceleration in different basins across the ice sheet and suggest that these spatiotemporal trends could be useful for calibrating and validating prognostic ice sheet models. In addition to resolving the spatial and temporal fingerprint of Greenland's recent ice loss, these mass loss grids are key for partitioning contemporary elastic vertical land motion from longer-term glacial isostatic adjustment (GIA) trends at GPS stations around the ice sheet. Our ice-loss product results in a significantly different GIA interpretation from a previous ice-loss product.

Accelerating Ice Loss From Peripheral Glaciers in North Greenland.

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.

Increased variability in Greenland Ice Sheet runoff from satellite observations.

Runoff from the Greenland Ice Sheet has increased over recent decades affecting global sea level, regional ocean circulation, and coastal marine ecosystems, and it now accounts for most of the contemporary mass imbalance. Estimates of runoff are typically derived from regional climate models because satellite records have been limited to assessments of melting extent. Here, we use CryoSat-2 satellite altimetry to produce direct measurements of Greenland's runoff variability, based on seasonal changes in the ice sheet's surface elevation. Between 2011 and 2020, Greenland's ablation zone thinned on average by 1.4 ± 0.4 m each summer and thickened by 0.9 ± 0.4 m each winter. By adjusting for the steady-state divergence of ice, we estimate that runoff was 357 ± 58 Gt/yr on average - in close agreement with regional climate model simulations (root mean square difference of 47 to 60 Gt/yr). As well as being 21 % higher between 2011 and 2020 than over the preceding three decades, runoff is now also 60 % more variable from year-to-year as a consequence of large-scale fluctuations in atmospheric circulation. Because this variability is not captured in global climate model simulations, our satellite record of runoff should help to refine them and improve confidence in their projections.

Steep Glacier Bed Knickpoints Mitigate Inland Thinning in Greenland.

Greenland's outlet glaciers have been a leading source of mass loss and accompanying sea-level rise from the Greenland Ice Sheet (GrIS) over the last 25 years. The dynamic component of outlet glacier mass loss depends on both the ice flux through the terminus and the inland extent of glacier thinning, initiated at the ice-ocean interface. Here, we find limits to the inland spread of thinning that initiates at glacier termini for 141 ocean-terminating outlet glaciers around the GrIS. Inland diffusion of thinning is limited by steep reaches of bed topography that we call "knickpoints." We show that knickpoints exist beneath the majority of outlet glaciers but they are less steep in regions of gentle bed topography, giving glaciers in gentle bed topography the potential to contribute to ongoing and future mass loss from the GrIS by allowing the diffusion of thinning far into the ice sheet interior.

Ocean forcing drives glacier retreat in Greenland.

The retreat and acceleration of Greenland glaciers since the mid-1990s have been attributed to the enhanced intrusion of warm Atlantic Waters (AW) into fjords, but this assertion has not been quantitatively tested on a Greenland-wide basis or included in models. Here, we investigate how AW influenced retreat at 226 marine-terminating glaciers using ocean modeling, remote sensing, and in situ observations. We identify 74 glaciers in deep fjords with AW controlling 49% of the mass loss that retreated when warming increased undercutting by 48%. Conversely, 27 glaciers calving on shallow ridges and 24 in cold, shallow waters retreated little, contributing 15% of the loss, while 10 glaciers retreated substantially following the collapse of several ice shelves. The retreat mechanisms remain undiagnosed at 87 glaciers without ocean and bathymetry data, which controlled 19% of the loss. Ice sheet projections that exclude ocean-induced undercutting may underestimate mass loss by at least a factor of 2.

Low elevation of Svalbard glaciers drives high mass loss variability.

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.

Rapid ablation zone expansion amplifies north Greenland mass loss.

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.

Forty-six years of Greenland Ice Sheet mass balance from 1972 to 2018.

We reconstruct the mass balance of the Greenland Ice Sheet using a comprehensive survey of thickness, surface elevation, velocity, and surface mass balance (SMB) of 260 glaciers from 1972 to 2018. We calculate mass discharge, D, into the ocean directly for 107 glaciers (85% of D) and indirectly for 110 glaciers (15%) using velocity-scaled reference fluxes. The decadal mass balance switched from a mass gain of +47 ± 21 Gt/y in 1972-1980 to a loss of 51 ± 17 Gt/y in 1980-1990. The mass loss increased from 41 ± 17 Gt/y in 1990-2000, to 187 ± 17 Gt/y in 2000-2010, to 286 ± 20 Gt/y in 2010-2018, or sixfold since the 1980s, or 80 ± 6 Gt/y per decade, on average. The acceleration in mass loss switched from positive in 2000-2010 to negative in 2010-2018 due to a series of cold summers, which illustrates the difficulty of extrapolating short records into longer-term trends. Cumulated since 1972, the largest contributions to global sea level rise are from northwest (4.4 ± 0.2 mm), southeast (3.0 ± 0.3 mm), and central west (2.0 ± 0.2 mm) Greenland, with a total 13.7 ± 1.1 mm for the ice sheet. The mass loss is controlled at 66 ± 8% by glacier dynamics (9.1 mm) and 34 ± 8% by SMB (4.6 mm). Even in years of high SMB, enhanced glacier discharge has remained sufficiently high above equilibrium to maintain an annual mass loss every year since 1998.

Atmospheric forcing of rapid marine-terminating glacier retreat in the Canadian Arctic Archipelago.

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.

Nonlinear rise in Greenland runoff in response to post-industrial Arctic warming.

The Greenland ice sheet (GrIS) is a growing contributor to global sea-level rise1, with recent ice mass loss dominated by surface meltwater runoff2,3. Satellite observations reveal positive trends in GrIS surface melt extent4, 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 rates5,6 and satellite-derived melt duration4 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.

Seasonal to decadal variability in ice discharge from the Greenland Ice Sheet.

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.

Direct measurements of meltwater runoff on the Greenland ice sheet surface.

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-km2 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.

Distinct patterns of seasonal Greenland glacier velocity.

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.