Search

The Northeast Greenland Ice Stream (NEGIS) has experienced substantial dynamic thinning in recent years. Here, we examine the evolving behavior of NEGIS, with focus on summer speedup at Zachariae Isstrøm, one of the NEGIS outlet glaciers, which has exhibited rapid retreat and acceleration, indicative of its vulnerability to changing climate conditions. Through a combination of Sentinel‐1 data, in‐situ GPS observations, and numerical ice flow modeling from 2007, we investigate the mechanisms driving short‐term changes. Our analysis reveals a summer speedup in ice flow both near the terminus and inland, with satellite data detecting changes up to 60 km inland, while GPS data capture changes up to 190 km inland along the glacier center line. We attribute this summer speedup to variations in subglacial hydrology, where surface meltwater runoff influences basal friction over the melt season. Incorporating subglacial hydrology into numerical models makes it possible to replicate observed ice velocity patterns. Plain Language Summary: The Northeast Greenland Ice Stream (NEGIS), a crucial part of the Greenland Ice Sheet, has been experiencing significant dynamic thinning recently. This study focuses on the summer speedup of Zachariae Isstrøm (ZI), one of NEGIS's outlet glaciers, which is rapidly retreating and accelerating, highlighting its sensitivity to climate change. Utilizing Sentinel‐1 satellite data, in‐situ GPS observations, and numerical ice flow modeling, we explore the mechanisms behind short‐term dynamic changes. We find that satellite data reveals short‐term summer (June to August) fluctuations in ice flow speed near the glacier terminus and up to 50–70 km inland. However, GPS data shows that this speedup extends further inland, up to at least 190 km along the main flow line. Only GPS data can detect the smaller‐scale summer speedups in these inland regions, providing critical observations for validating ice flow models. We determine that the seasonal acceleration of ice velocity at Zachariae Isstrøm is due to variations in subglacial hydrology, where surface meltwater runoff reduces basal friction by altering the subglacial hydrologic system during the melt season. Additionally, our study highlights that these findings are applicable beyond NEGIS, with similar speedup patterns observed in other Greenland glaciers. Key Points: GPS data reveal summer speed up at least 190 km inland along the main flowline of Zachariae IsstrømSubglacial hydrology is the main driver of the summer speedup near the terminus and deep inlandRecord high warming in 2019 led to a more intense and longer duration of the summer speedup [ABSTRACT FROM AUTHOR]

Monitoring the Greenland mass balance (GMB) is crucial in the context of global sea level rise. Currently, three main methods are used to measure GMB, with the primary source of uncertainty arising from the glacial isostatic adjustment (GIA) contribution. Here, we propose a novel approach based on a simple methodology that uses the entire Greenland GNSS network (GNET) as an instrument to monitor the present‐day mass changes. Our method is validated against GRACE‐derived GMB, and we find a very good agreement. This leads to an independent methodology for estimating present‐day mass changes from GNSS, bridging the gap between GRACE and GRACE‐FO in GMB estimates. Through a combined analysis of GMB from GRACE and GNET, we identify a consistency relation between the gravity and uplift signature of GIA, providing a new robust constraint for GIA models. Plain Language Summary: The rate of Greenland ice mass change is constantly monitored by satellites measuring gravity and surface elevation variations, airborne surveys, and a growing network of GNSS stations on the ground. Each technique has its strengths and weaknesses and spatial/temporal resolution. We developed an approach that uses the whole GNSS network as a "virtual instrument", capable of measuring the total ice loss from Greenland. Comparing this with more standard mass balance estimate techniques, we validate the method, and, by combining these techniques, we find a novel constraint on both the mass variations and glacial isostatic adjustment models. Key Points: New methodology for monitoring daily estimates of the Greenland Mass Balance (GMB) derived from the GNSS networkNew glacial isostatic adjustment constraint from combined GRACE and GNSS observationsFilling the gap in the GMB between GRACE and GRACE‐FO using GNSS data [ABSTRACT FROM AUTHOR]

In Greenland, 87% of the glacierized area terminates in the ocean, but mass lost at the ice‐ocean interface, or frontal ablation, has not yet been fully quantified. Using measurements and models we calculate frontal ablation of Greenland's 213 outlet and 537 peripheral glaciers and find a total frontal ablation of 481.8 ± 24.0 for 2000–2010 and 510.2 ± 18.6 Gt a−1 for 2010–2020. Ice discharge accounted for ∼90% of frontal ablation during both periods, while mass loss due to terminus retreat comprised the remainder. Only 16 glaciers were responsible for the majority (>50%) of frontal ablation from 2010 to 2020. These estimates, along with the climatic‐basal balance, allow for a more complete accounting of Greenland Ice Sheet and peripheral glacier mass balance. In total, Greenland accounted for ∼90% of Northern Hemisphere frontal ablation for 2000–2010 and 2010–2020. Plain Language Summary: We estimate the mass of ice lost from all Greenland glaciers that entered the ocean during each of the last two decades. This ice loss at the front of these marine‐terminating glaciers is called frontal ablation and is approximately equal to the mass of icebergs entering the ocean. Frontal ablation is important because 87% of glacier area in Greenland ends in the ocean, through 750 outlets, and previous work has only approximated frontal ablation. This study quantifies it for the first time, helping to close the mass budget for the Greenland Ice Sheet and better partition its mass balance into components. We find that Greenland accounts for ∼90% of all Northern Hemisphere frontal ablation and, of that contribution, just 17 glaciers for 2000–2010 and 16 glaciers for 2010–2020 account for more than half of total Greenland frontal ablation. Key Points: Frontal ablation of the Greenland Ice Sheet averaged 510.2 ± 18.6 Gt a−1 for 2010–2020, ∼90% of which came from ice dischargeThe frontal ablation we measured is larger than the total mass loss from the ice sheet, indicating a positive climatic‐basal balanceOnly 16 glaciers account for 50% of the total frontal ablation from the Greenland Ice Sheet [ABSTRACT FROM AUTHOR]

Here, we present a compilation of 95 ice temperature profiles from 85 boreholes from the Greenland ice sheet and peripheral ice caps, as well as local ice caps in the Canadian Arctic. Profiles from only 31 boreholes (36 %) were previously available in open-access data repositories. The remaining 54 borehole profiles (64 %) are being made digitally available here for the first time. These newly available profiles, which are associated with pre-2010 boreholes, have been submitted by community members or digitized from published graphics and/or data tables. All 95 profiles are now made available in both absolute (meters) and normalized (0 to 1 ice thickness) depth scales and are accompanied by extensive metadata. These metadata include a transparent description of data provenance. The ice temperature profiles span 70 years, with the earliest profile being from 1950 at Camp VI, West Greenland. To highlight the value of this database in evaluating ice flow simulations, we compare the ice temperature profiles from the Greenland ice sheet with an ice flow simulation by the Parallel Ice Sheet Model (PISM). We find a cold bias in modeled near-surface ice temperatures within the ablation area, a warm bias in modeled basal ice temperatures at inland cold-bedded sites, and an apparent underestimation of deformational heating in high-strain settings. These biases provide process level insight on simulated ice temperatures. [ABSTRACT FROM AUTHOR]

The magnitude and azimuth of horizontal ice flow at Camp Century, Greenland have been measured several times since 1963. Here, we provide a further two independent measurements over the 2017–21 period. Our consensus estimate of horizontal ice flow from four independent satellite-positioning solutions is 3.65 ± 0.13 m a−1 at an azimuth of 236 ± 2°. A portion of the small, but significant, differences in ice velocity and azimuth reported between studies likely results from spatial gradients in ice flow. This highlights the importance of restricting inter-study comparisons of ice flow estimates to measurements surveyed within a horizontal distance of one ice thickness from each other. We suggest that ice flow at Camp Century is stable on seasonal to multi-decadal timescales. The airborne and satellite laser altimetry record indicates an ice thickening trend of 1.1 ± 0.3 cm a−1 since 1994. This thickening trend is qualitatively consistent with previously inferred ongoing millennial-scale ice thickening at Camp Century. The ice flow divide immediately north of Camp Century may now be migrating southward, although the reasons for this divide migration are poorly understood. The Camp Century flowlines presently terminate in the vicinity of Innaqqissorsuup Oqquani Sermeq (Gade Gletsjer) on the Melville Bay coast. [ABSTRACT FROM AUTHOR]

The largest floating tongue of Greenland's ice sheet, Nioghalvfjerdsbræ, has been relatively stable with respect to areal retreat until 2022. Draining more than 6 % of the ice sheet, a disintegration of Nioghalvfjerdsbræ's floating tongue and subsequent acceleration due to loss in buttressing are likely to lead to sea level rise. Therefore, the stability of the floating tongue is a focus of this study. We employed a suite of observational methods to detect recent changes at the calving front. We found that the calving style has changed since 2016 at the southern part of the eastern calving front, from tongue-type calving to a crack evolution initiated at frontal ice rises reaching 5–7 km and progressing further upstream compared to 2010. The calving front area is further weakened by an area upstream of the main calving front that consists of open water and an ice mélange that has substantially expanded, leading to the formation of a narrow ice bridge. These geometric and mechanical changes may be a precursor of instability of the floating tongue. We complement our study by numerical ice flow simulations to estimate the impact of future ice-front retreat and complete ice shelf disintegration on the discharge of grounded ice. These idealized scenarios reveal that a loss of the south-eastern area of the ice shelf would lead to a 0.2 % increase in ice discharge at the grounding line, while a sudden collapse of the frontal area (46 % of the floating tongue area) will enhance the ice discharge by 5.1 % due to loss in buttressing. Eventually, a full collapse of the floating tongue increases the grounding line flux by 166 %. [ABSTRACT FROM AUTHOR]

Geothermal heat flow is key to unraveling several large-scale geophysical systems, including the inner workings of the Greenlandic ice sheet, and by extension, the possibility of understanding the past and prior global climate. Similarly, it could provide insight into the paleo-trace of the Icelandic mantle plume, which in turn is integral in answering long-standing questions on the origin of mountains in western and eastern Greenland and in Norway. This study documents the results from an intra-scientific field approach, which combines geological, petrophysical, and satellite magnetic field data in a nonlinear probabilistic inversion. These results include Curie depths with associated uncertainties and Geothermal Heat Flux estimates. While baselines remain challenging to evaluate due to the strong nonlinearity of the problem posed, stress testing reveals a high robustness of the predicted spatial variations, which largely disagree with the classic straightforward northwest–southeast or east–west plume trace across Greenland. Instead, our results indicate a complex heat flux pattern, including a localized region with anomalously heightened heat flux near the origin of the Northeast Greenland Ice Stream. [ABSTRACT FROM AUTHOR]

Glacial lake outburst floods (GLOFs) or jökulhlaups from ice-dammed lakes are frequent in Greenland and can influence local ice dynamics and bedrock motion, cause geomorphological changes, and pose flooding hazards. Multidecadal time series of lake drainage dates, volumes, and flood outlets are extremely rare. However, they are essential for determining the scale and frequency of future GLOFs, for identifying drainage mechanisms, and for mitigating downstream flood effects. In this study, we use high-resolution digital elevation models (DEMs) and orthophotos (0.1 × 0.1 m) generated from uncrewed-aerial-vehicle (UAV) field surveys, in combination with optical satellite imagery. This allows us to reconstruct robust lake volume changes associated with 14 GLOFs between 2007 and 2021 at Russell Glacier, West Greenland. As a result, this is one of the most comprehensive and longest records of ice-dammed lake drainages in Greenland to date. Importantly, we find a mean difference of ∼ 10 % between our lake drainage volumes when compared with estimates derived from a gauged hydrograph 27 km downstream. Due to thinning of the local ice dam, the potential maximum drainage volume in 2021 is ∼ 60 % smaller than that estimated to have drained in 2007. Our time series also reveals variations in the drainage dates ranging from late May to mid-September and drainage volumes ranging between 0.9 and 37.7 Mm 3. We attribute these fluctuations between short periods of relatively high and low drainage volumes to a weakening of the ice dam and an incomplete sealing of the englacial tunnel following the large GLOFs. This syphoning drainage mechanism is triggered by a reduction in englacial meltwater, likely driven by late-season drainage and sudden air temperature reductions, as well as annual variations in the glacial drainage system. Furthermore, we provide geomorphological evidence of an additional drainage route first observed following the 2021 GLOF, with a subglacial or englacial flow pathway, as well as supraglacial water flow across the ice margin. It seems probable that the new drainage route will become dominant in the future. This will drive changes in the downstream geomorphology and raise the risk of flooding-related hazards as the existing buffering outlet lakes will be bypassed. [ABSTRACT FROM AUTHOR]

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. Plain Language Summary: The Arctic is warming more rapidly than the rest of the world. This warming has had an especially profound impact on Greenland's ice cover. Only 4% of Greenland's ice cover are small peripheral glaciers that are distinct from the ice sheet proper. Despite comprising this relatively small area, these small peripheral glaciers are responsible for 11% of the ice loss associated with Greenland's recent sea level rise contribution. Using the satellite laser platforms Ice, Cloud, and land Elevation Satellite (ICESat) and ICESat‐2, we estimate that ice loss from these Greenland glaciers increased from 27 ± 6 Gt/yr (2003–2009) to 42 ± 6 Gt/yr (2018–2021). We find that the largest acceleration in ice loss is in North Greenland, where we observe ice loss to increase by a factor of four between 2003 and 2021. In some areas, it appears that recent increases in snowfall at high altitudes have partially counteracted recent increases in melt at low altitudes. While many recent Greenland ice loss assessments have focused on only the ice sheet, the recent sharp increase in ice loss from small peripheral glaciers highlights the importance of accurately monitoring Greenland's small peripheral glaciers. These small peripheral glaciers appear poised to play an outsized role in Greenland ice loss for decades to come. Key Points: Greenland peripheral glacier mass loss increased from 27.2 Gt/yr (February 2003–October 2009) to 42.3 Gt/yr (October 2018–December 2021Mass loss increased by a factor of four for peripheral glaciers in North Greenland during February 2003–December 2021Enhanced precipitation in northeast Greenland during October 2018–December 2021 resulted in ice thickening at higher elevations [ABSTRACT FROM AUTHOR]

Studies have shown that geodetic Global Navigation Satellite System (GNSS) stations can be used to measure snow depths using GNSS interferometric reflectometry (GNSS-IR). Here, we study the results from a customized GNSS setup installed in March through August 2020 at the Programme for Monitoring of the Greenland Ice Sheet (PROMICE) automatic weather station NUK-K located on a small glacier outside Nuuk, Greenland. The setup is not optimized for reflectometry purposes. The site is obstructed between 85 and 215 degrees, and as the power supply is limited due to the remote location, the logging time is limited to 3 h per day. We estimate reflector heights using GNSS-IR and compare the results to a sonic ranger also placed on the weather station. We find that the snow melt measured by GNSS-IR is comparable to the melt measured by the sonic ranger. We expect that a period of up to 45 cm difference between the two is likely related to the much larger footprint GNSS-IR and the topography of the area. The uncertainty on the GNSS-IR reflector heights increase from approximately 2 cm for a snow surface to approximately 5 cm for an ice surface. If reflector height during snow free periods are part of the objective of a similar setup, we suggest increasing the logging time to reduce the uncertainty on the daily estimates. [ABSTRACT FROM AUTHOR]

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. Plain Language Summary: Greenland ice loss has accelerated over the last three decades. The Greenland Ice Sheet is currently one of the largest contributors to global sea‐level rise. We combine airborne and satellite altimetry measurements to make annual digital elevation models of the ice sheet during 2011–2020. Over this period, we find that the ice sheet lost an ice volume corresponding to 6.9 ± 0.4 mm of global sea‐level equivalent. The peak loss year was April 2019 to April 2020, when the ice sheet lost 1.4 mm of global sea‐level equivalent. This peak loss rate is equivalent to losing 15,850 tonnes of ice per second for 12 months. We also find that the acceleration of ice loss differs across different ice‐sheet sectors. We suggest that these regional ice loss trends may be a good target for the ice‐sheet models used to project future ice loss. Key Points: The Greenland Ice Sheet contributed 6.9 ± 0.4 mm to sea‐level from April 2011 to April 2020Satellite altimetry suggests a peak annual ice loss of 498 ± 45 Gt from April 2019 to April 2020The terminus of Jakobshavn Isbræ is once again dynamically thinning, following a period of dynamic thickening during 2016–2018 [ABSTRACT FROM AUTHOR]

Northeast Greenland is the place where the Greenland Ice Sheet (GrIS) experienced the largest areal changes since the Last Glacial Maximum. However, the age constraints of the last deglaciation are in some areas sparse. In this study, we use forty-seven new 10Be cosmogenic exposure ages to constrain the deglaciation of the present-day ice-free areas in Northeast Greenland. Our results show that the outer coast region was deglaciated between 12.8 ± 0.6 and 11.5 ± 0.2 ka and the region close to the present ice margin was deglaciated 2 to 4 ka later between 9.2 ± 0.3 to 8.6 ± 0.3 ka. By combining our new results with previously published data from North and Northeast Greenland, we find that the ice sheet advanced to the shelf edge between 26 and 20 cal. ka BP. The outer coast was deglaciated between 12.8 and 9.7 ka and the present ice extent was reached between 10.8 to 5.8 ka. The ice margin continued to retreat farther inland during the Middle Holocene before it readvanced toward its Little Ice Age position. The deglaciation was probably forced by a combination of increased atmospheric and ocean temperatures, but local topography also played an important role. These results add to the growing knowledge about the glaciation history of the GrIS and add useful constraints for future ice sheet models. [ABSTRACT FROM AUTHOR]

Greenland glaciers exhibit variable seasonal velocity signals that may reflect differences in subglacial hydrology. Here, we conduct a first GrIS-wide glacier classification based on seasonal velocity patterns derived from 2017 Sentinel-1 radar data. Our classification focuses on two distinct seasonal ice velocity patterns, with the first (type-2 from Moon and others, 2014) showing periods of both speedup and slowdown during the melt season, and the second (type-3) instead showing a longer period of slowdown from elevated velocities in the winter and spring. We analyze 221 glaciers in 2017 and show that 48 exhibit type-2 behavior, and 72 exhibit type-3 behavior. We extend the classification to 2018 and 2019 and find that while the glaciers meeting each criterion vary year to year, type-2 is consistently more common in the northern regions and type-3 is more common in the south. Our results highlight the varied impact of meltwater on subglacial drainage systems and glacier flow in Greenland. [ABSTRACT FROM AUTHOR]

We present a novel method to estimate dynamic ice loss of Greenland's three largest outlet glaciers: Jakobshavn Isbræ, Kangerlussuaq Glacier, and Helheim Glacier. We use Global Navigation Satellite System (GNSS) stations attached to bedrock to measure elastic displacements of the solid Earth caused by dynamic thinning near the glacier terminus. When we compare our results with discharge, we find a time lag between glacier speedup/slowdown and onset of dynamic thinning/thickening. Our results show that dynamic thinning/thickening on Jakobshavn Isbræ occurs 0.87 ± 0.07 years before speedup/slowdown. This implies that using GNSS time series we are able to predict speedup/slowdown of Jakobshavn Isbræ by up to 10.4 months. For Kangerlussuaq Glacier the lag between thinning/thickening and speedup/slowdown is 0.37 ± 0.17 years (4.4 months). Our methodology and results could be important for studies that attempt to model and understand mechanisms controlling short‐term dynamic fluctuations of outlet glaciers in Greenland. Plain Language Summary: A wide range of sensors and methods have been used to study the changes of the Greenland Ice Sheet, including satellite gravimetry, altimetry, and the input‐output method. Here, we present a novel fourth method to estimate dynamic ice loss of Greenland's three largest outlet glaciers: Jakobshavn Isbræ, Kangerlussuaq Glacier, and Helheim Glacier. We use Global Navigation Satellite System (GNSS) stations attached to bedrock to measure rise of land masses caused by ongoing ice mass loss near the glacier terminus. When we compare our results with ice discharge, we find a time lag between glacier speedup/slowdown and onset of dynamic induced thinning/thickening. Our results show that dynamic thinning/thickening on Jakobshavn Isbræ occurs 0.87 ± 0.07 years before speedup/slowdown. This implies that using GNSS uplift time series we are able to predict ice flow speedup/slowdown of Jakobshavn Isbræ by up to 10 months. For Kangerlussuaq Glacier and Helheim Glacier the lag between thinning/thickening and speedup/slowdown is 0.37 ± 0.17 years (4.4 months) and 0.03 ± 0.16 years, respectively. Our methodology and results could be important for studies that attempt to model and understand mechanisms controlling short‐term dynamic fluctuations of outlet glaciers in Greenland. Key Points: A novel method to estimate dynamic ice loss of Greenland's three largest outlet glaciers, Jakobshavn, Kangerlussuaq, and Helheim glacierDynamic thinning/thickening occurs 0.87 ± 0.07 years before speedup/slowdown at Jakobshavn IsbræA similar time lag between change in uplift rate and flow speed change allows us to predict future ice discharge from past uplift [ABSTRACT FROM AUTHOR]

Vertical land motion (VLM) from past and ongoing glacial changes can amplify or mitigate ongoing relative sea level change. We present a high‐resolution VLM model for the wider Arctic, that includes both present‐day ice loading (PDIL) and glacial isostatic adjustment (GIA). The study shows that the nonlinear elastic uplift from PDIL is significant (0.5–1 mm yr−1) in most of the wider Arctic and exceeds GIA at 15 of 54 Arctic GNSS sites, including sites in nonglaciated areas of the North Sea region and the east coast of North America. Thereby the sea level change from PDIL (1.85 mm yr−1) is significantly mitigated from VLM caused by PDIL. The combined VLM model was consistent with measured VLM at 85% of the GNSS sites (R=0.77) and outperformed a GIA‐only model (R=0.64). Deviations from GNSS‐measured VLM can be attributed to local circumstances causing VLM. Plain Language Summary: From 2003 to 2015, the Northern Hemisphere lost more than 6,000 gigatons of land ice, which added nearly 18 mm to the global mean sea level rise. Loss of land‐based ice results in the vertical deformation of the Earth's surface. Ongoing rebounding or subsidence caused by the end of the last ice age is often assumed to govern vertical deformation. However, present‐day ice loss from Greenland and Arctic glaciers also cause an immediate vertical deformation. By using an vertical deformation model, that includes both components, we can explain GPS‐measured deformation occurring in the Arctic. Our results show that the present‐day Arctic ice loss contribution to vertical deformation is ∼0.5 to 1 mm yr−1 in the wider northern region. This exceeds deformation caused by the disappearance of the last ice ages at many coastal regions, including the North Sea region and the North American Atlantic coast. The Arctic present‐day ice loss included in the VLM model equals a global sea level rise of 1.5 mm yr−1, which means that 30–80% of the sea level rise caused by Arctic ice loss is mitigated by surface uplift caused by the same ice loss. Key Points: Elastic VLM from present‐day ice loss in the Arctic causes significant uplift of coastlines in North America and Northern EuropeA combined VLM model that includes GIA and elastic VLM yields good agreement with GNSS stations in the wider ArcticResiduals between GNSS and modeled VLM provide an approximation of extraordinary VLM caused by local circumstances [ABSTRACT FROM AUTHOR]

The Northeast Greenland Ice Stream (NEGIS) extends around 600 km upstream from the coast to its onset near the ice divide in interior Greenland. Several maps of surface velocity and topography of interior Greenland exist, but their accuracy is not well constrained by in situ observations. Here we present the results from a GPS mapping of surface velocity in an area located approximately 150 km from the ice divide near the East Greenland Ice-core Project (EastGRIP) deep-drilling site. A GPS strain net consisting of 63 poles was established and observed over the years 2015–2019. The strain net covers an area of 35 km by 40 km, including both shear margins. The ice flows with a uniform surface speed of approximately 55 m a -1 within a central flow band with longitudinal and transverse strain rates on the order of 10 -4 a -1 and increasing by an order of magnitude in the shear margins. We compare the GPS results to the Arctic Digital Elevation Model and a list of satellite-derived surface velocity products in order to evaluate these products. For each velocity product, we determine the bias in and precision of the velocity compared to the GPS observations, as well as the smoothing of the velocity products needed to obtain optimal precision. The best products have a bias and a precision of ∼0.5 m a -1. We combine the GPS results with satellite-derived products and show that organized patterns in flow and topography emerge in NEGIS when the surface velocity exceeds approximately 55 m a -1 and are related to bedrock topography. [ABSTRACT FROM AUTHOR]

We present a 1986 through March 2020 estimate of Greenland Ice Sheet ice discharge. Our data include all discharging ice that flows faster than 100 m yr -1 and are generated through an automatic and adaptable method, as opposed to conventional handpicked gates. We position gates near the present-year termini and estimate problematic bed topography (ice thickness) values where necessary. In addition to using annual time-varying ice thickness, our time series uses velocity maps that begin with sparse spatial and temporal coverage and end with near-complete spatial coverage and 12 d updates to velocity. The 2010 through 2019 average ice discharge through the flux gates is ∼487±49 Gt yr -1. The 10 % uncertainty stems primarily from uncertain ice bed location (ice thickness). We attribute the ∼50 Gt yr -1 differences among our results and previous studies to our use of updated bed topography from BedMachine v3. Discharge is approximately steady from 1986 to 2000, increases sharply from 2000 to 2005, and then is approximately steady again. However, regional and glacier variability is more pronounced, with recent decreases at most major glaciers and in all but one region offset by increases in the northwest region through 2017 and in the southeast from 2017 through March 2020. As part of the journal's living archive option and our goal to make an operational product, all input data, code, and results from this study will be updated as needed (when new input data are available, as new features are added, or to fix bugs) and made available at https://doi.org/10.22008/promice/data/ice%5fdischarge and at https://github.com/mankoff/ice%5fdischarge (last access: 6 June 2020,). [ABSTRACT FROM AUTHOR]

The GRACE mission has had a revolutionary impact on the study of Earth system processes, but it provides a band‐limited representation of mass changes. This is particularly problematic when studying mass changes that tend to be concentrated in fairly narrow zones near the edges of the Greenland Ice Sheet. In this study, coarse‐resolution estimates of the mass change derived from GRACE have been enhanced by the introduction of heuristic scaling factors applied to model surface mass balance and observed surface elevation change. Corresponding results indicate large spatial heterogeneity in the gridded scaling factors at the 0.5° × 0.5° scale, reflecting significant mass losses concentrated along the ice sheet margin and relatively small internal ice sheet changes at higher elevations. The scaled GRACE‐derived vertical displacements are in the range from −2 to 14 mm/year from 2003 to 2015. The Greenland GPS network was used to examine the downscaling GRACE predictions of the crustal displacements. The results show consistency of scaled GRACE‐predicted and GPS‐observed seasonal and long‐term uplift in major drainage basins of Greenland. Our results indicate that GRACE predictions underestimate vertical displacements at sites located in regions characterized by concentrated loads, but perform well in other regions. Differences between predicted and observed uplift rates are mainly caused by the sensitivity kernels, because GPS and GRACE estimates are based on weighted averages of mass loss in different sensitivity ranges. Moreover, a large uncertainty in the glacial isostatic adjustment correction may also cause errors in the GPS‐to‐GRACE ratio. Key Points: GRACE results can be improved by using scaling factors estimated from two models of surface mass fields in GreenlandScaled GRACE‐derived and GPS‐observed vertical displacements were compared at 53 Greenland GPS Network (GNET) stationsGRACE‐based estimates of accurate uplift rates in Greenland needs to consider the spatial pattern of ice loss and reliable GIA correction [ABSTRACT FROM AUTHOR]

We present a 1986 through 2017 estimate of Greenland Ice Sheet ice discharge. Our data include all discharging ice that flows faster than 100 m yr−1 and are generated through an automatic and adaptable method, as opposed to conventional hand-picked gates. We position gates near the present-year termini and estimate problematic bed topography (ice thickness) values where necessary. In addition to using annual time-varying ice thickness, our time series uses velocity maps that begin with sparse spatial and temporal coverage and end with near-complete spatial coverage and 6 d updates to velocity. The 2010 through 2017 average ice discharge through the flux gates is ∼ 488 ± 48Gt yr−1. The 10 % uncertainty stems primarily from uncertain ice bed location (ice thickness). We attribute the ∼50 Gt yr−1 differences among our results and previous studies to our use of updated bed topography from BedMachine v3. Discharge is approximately steady from 1986 to 2000, increases sharply from 2000 to 2005, and then is approximately steady again. However, regional and glacier variability is more pronounced, with recent decreases at most major glaciers and in all but one region offset by increases in the NW (northwestern) region. As part of the journal's living archive option, all input data, code, and results from this study will be updated when new input data are accessible and made freely available at https://doi.org/10.22008/promice/data/ice_discharge. [ABSTRACT FROM AUTHOR]

The Greenland ice sheet has been one of the largest contributors to global sea-level rise over the past 20 years, accounting for 0.5 mm yr−1 of a total of 3.2 mm yr−1. A significant portion of this contribution is associated with the speed-up of an increased number of glaciers in southeast and northwest Greenland. Here, we show that the northeast Greenland ice stream, which extends more than 600 km into the interior of the ice sheet, is now undergoing sustained dynamic thinning, linked to regional warming, after more than a quarter of a century of stability. This sector of the Greenland ice sheet is of particular interest, because the drainage basin area covers 16% of the ice sheet (twice that of Jakobshavn Isbræ) and numerical model predictions suggest no significant mass loss for this sector, leading to an under-estimation of future global sea-level rise. The geometry of the bedrock and monotonic trend in glacier speed-up and mass loss suggests that dynamic drawdown of ice in this region will continue in the near future. [ABSTRACT FROM AUTHOR]

Abstract: We estimate the mass loss rate of Upernavik Isstrøm (UI) using surface elevation changes between a SPOT 5 Digital Elevation Model (DEM) from 2008 and NASA''s Airborne Topographic Mapper (ATM) data from 2010. To assess the validity of our mass loss estimate, we analyze GPS data between 2007 and 2011 from two continuous receivers, UPVK and SRMP which are established on bedrock and located and from the front of UI, respectively. We construct along-track elevation changes on UI for several time intervals during 2005–2011, based on ATM, SPOT 5 and Ice, Cloud, and land Elevation Satellite (ICESat) data to assess temporal changes of UI. We estimate a mass loss rate of −6.7±4.2Gt/yr, over an area of . The ice mass loss occurs primarily over the northern glacier of UI. This pattern is also observed upstream, where we observe glacier thinning at a rate of −1.6±0.3m/yr across the northern portion of UI and −0.5±0.1m/yr across the southern portion. GPS measurements suggest bedrock uplift rates of 7.6±0.6mm/yr (UPVK) and 16.2±0.6mm/yr (SRMP). The modeled ice mass loss of UI causes bedrock uplift rates of 1.3±0.6mm/yr (UPVK) and 8.3±4.2mm/yr (SRMP). Including additional contributions from ice mass changes outside UI and from Glacial Isostatic Adjustment (GIA), we obtain total modeled uplift rates of 4.7±0.6mm/yr (UPVK) and 13.8±4.2mm/yr (SRMP). The modeled uplift rates from our UI ice mass loss are substantially lower, indicating that additional mass loss is taking place outside of UI. We obtain a difference of 0.6mm/yr between the modeled and observed relative uplift rates (SRMP relative to UPVK), suggesting that the mass loss of UI is well captured in the model. We observe elevation changes from −15 to −40m/yr near the front during the period 2005–2011, indicating that UI undergoes large variations in thinning pattern over short time spans. [Copyright &y& Elsevier]

Global warming is predicted to have a profound impact on the Greenland Ice Sheet and its contribution to global sea-level rise. Recent mass loss in the northwest of Greenland has been substantial. Using aerial photographs, we produced digital elevation models and extended the time record of recent observed marginal dynamic thinning back to the mid-1980s. We reveal two independent dynamic ice loss events on the northwestern Greenland Ice Sheet margin: from 1985 to 1993 and 2005 to 2010, which were separated by limited mass changes. Our results suggest that the ice mass changes in this sector were primarily caused by short-lived dynamic ice loss events rather than changes in the surface mass balance. This finding challenges predictions about the future response of the Greenland Ice Sheet to increasing global temperatures. [ABSTRACT FROM AUTHOR]

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 GNETstations during an approximate six-month period in 2010. This anomalous uplift is spatially correlated with the 2010 melting day anomaly. [ABSTRACT FROM AUTHOR]

SUMMARY Jakobshavn Isbræ in west Greenland has been undergoing dramatic thinning since 1997. Applying the interferometric synthetic aperture radar (InSAR) technique to Radarsat-1 SAR data, we measure crustal uplift near Jakobshavn Isbræ caused by recent ice mass loss. The crustal uplift is predominantly at long spatial wavelengths (larger than 10 km), and thus is difficult to separate from InSAR orbit errors. We reduce the effects of orbit errors by removing long-wavelength deformation signals using conventional InSAR baseline fitting methods. We find good agreement between the remaining short-scale InSAR-estimated deformation rates during 2004-2008 and the corresponding short-scale components of a deformation model that is based on changes in ice elevation measured by NASA's Airborne Topographic Mapper (ATM). We are also able to use the InSAR-measured deformation to invert for the spatial pattern of ice thinning. Overall, our results suggest that despite the inherent difficulties of working with a signal that has significant large-scale components, InSAR-measured crustal deformation can be used to study the ice mass loss of a rapidly thinning glacier and its surrounding catchment, providing both a constraint on any existing model of ice mass loss and a data source that can be used to invert for ice mass loss. These new applications of InSAR can help to better understand a glacier's rapid response to a warming climate. [ABSTRACT FROM AUTHOR]

Studies of global sea level often exclude Tide Gauges (TGs) in glaciated regions due to vertical land movement. Recent studies show that geodetic GNSS stations can be used to estimate sea level by taking advantage of the reflections from the ocean surface using GNSS Interferometric Reflectometry (GNSS-IR). This method has the immediate benefit that one can directly correct for bedrock movements as measured by the GNSS station. Here we test whether GNSS-IR can be used for measurements of inter annual sea level variations in Thule, Greenland, which is affected by sea ice and icebergs during much of the year. We do this by comparing annual average sea level variations using the two methods from 2008–2019. Comparing the individual sea level measurements over short timescales we find a root mean square deviation (RMSD) of 13 cm, which is similar to other studies using spectral methods. The RMSD for the annual average sea level variations between TG and GNSS-IR is large (18 mm) compared to the estimated uncertainties concerning the measurements. We expect that this is in part due to the TG not being datum controlled. We find sea level trends from GNSS-IR and TG of −4 and −7 mm/year, respectively. The negative trend can be partly explained by a gravimetric decrease in sea level as a result of ice mass changes. We model the gravimetric sea level from 2008–2017 and find a trend of −3 mm/year. [ABSTRACT FROM AUTHOR]

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. The Greenland Ice Sheet is the largest land ice contributor to sea level rise and understanding the long-term glacier response to external forcing is key to improved projections. Here the authors show Greenland's three largest outlet glaciers will likely exceed current worst-case scenario [ABSTRACT FROM AUTHOR]

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. We provide 3D displacement time series from GPS available at (https://catalog-intaros.nersc.no/dataset/ice-mass-change-of-the-greenland-ice-sheet) to study uplift in response to present-day changes in ice mass. We compare present-day 3D deformations with modeled results. To retrieve 3D elastic displacements from GPS time series, we correct our observations for glacial-isostatic adjustment and tectonic plate motion. To model 3D elastic displacements, we first estimate mass loss using 1995–2014 NASA's Airborne Topographic Mapper (ATM) flights derived altimetry, supplemented with laser altimetry observations from the Ice, Cloud, and Land Elevation Satellite (ICESat) during 2003–2009; the airborne Land, Vegetation, and Ice Sensor (LVIS) instrument during 2007–2013; radar altimetry from the CryoSat-2 satellite during 2010–2018; and European Remote-Sensing Satellite–1 (ERS-1) and ERS-2 data during 1995–2003. We convert the volume loss rate into a mass loss rate accounting for firncompaction. We predict the elastic displacements by convolving mass loss estimates with Green's functions for vertical and horizontal displacements. [ABSTRACT FROM AUTHOR]

• Geodetic data and model reveal transient mass changes in Greenland. • NW Greenland exhibited different change patterns than the rest subregions. • Precipitation anomalies were opposite in signs in east & west Greenland. • Atmospheric circulation anomalies may cause most of the transient mass changes. Much of the research to understand the ice mass changes of Greenland ice sheet (GrIS) has focused on detecting linear rates and accelerations at decadal or longer periods. The transient (short-term, non-secular) mass changes show large variability, and if not properly accounted for, can introduce significant biases into estimates of long-term ice mass loss rates and accelerations. Despite the growing number of geodetic observations, in terms of spatial coverage, types of observables, and the extent of the time series, studies of the transient mass changes over GrIS are lacking. To address this limitation, we apply multi-channel singular spectral analysis to the Gravity Recovery and Climate Experiment (GRACE) mass concentrations (mascon), surface mass balance (SMB) model output, and ice discharge data, to determine the transient mass changes over Greenland over the decade (2007 to 2017). The goal of this analysis is to elucidate the spatio-temporal variability of the ice mass change. For the entire GrIS, both the mascon and SMB transient mass changes are characterized by a sustained mass gain from late 2007 to early 2010, a sustained mass loss from early 2010 to early 2013, and a mass gain from early 2013 to mid-2015. Global Positioning System sites deployed along the coast of Greenland showed uplift from early 2010 to early 2013 and subsidence from early 2013 to 2015, consistent with the corresponding ice mass loss and gain of the entire GrIS. The peak-to-peak amplitude of the transient mass change was estimated to be −294 ± 27 Gt from GRACE mascons and -252 ± 16 Gt from the SMB where the latter value includes the effect of ice discharge. The transient mass change due to ice discharge accounted for less than 10% of the total transient mass change. Our regional assessment reveals that the central-west, southwest, northeast, and southeast regions display similar time-varying patterns as we found for the entire GrIS, but the north and northwest regions show different patterns. Atmospheric circulation anomalies as measured by the Greenland Blocking Index (GBI) are able to explain most of these transient anomalies. More specifically, high-GBI-associated high temperature was one of the main reasons for the transient mass loss of the entire GrIS during 2010-2012 while low GBI can explain the transient mass gain during 2013-2015. Contrasting behaviors of precipitation anomalies in east and west Greenland under abnormally high or low GBI conditions may explain the different patterns of the transient mass change in the northwest and the rest of Greenland. [ABSTRACT FROM AUTHOR]

We investigate three annual mass-balance cycles on Helheim Glacier in south-east Greenland using TanDEM-X interferometric digital elevation models (DEMs), bedrock GPS measurements, and ice velocity from feature-tracking. The DEMs exhibit seasonal surface elevation cycles at elevations up to 800 m.a.s.l. with amplitudes of up to 19 m, from a maximum in July to a minimum in October or November, concentrated on the fast-flowing areas of the glacier indicating that the elevation changes have a mostly dynamic origin. By modelling the detrended bedrock loading/unloading signal we estimate a mean density for the loss of 671 ± 70 kg m − 3 and calculate that total water equivalent volume loss from the active part of the glacier (surface flow speeds >1 m day −1 ) ranges from 0.5 km 3 in 2011 to 1.6 km 3 in 2013. A rough ice-flux divergence analysis shows that at lower elevations (<200 m) mass loss by dynamic thinning fully explains seasonal elevation changes. In addition, surface elevations decrease by a greater amount than field observations of surface ablation or surface-energy-balance modelling predict, emphasising the dynamic nature of the mass loss. We conclude, on the basis of ice-front position observations through the time series, that melt-induced acceleration is most likely the main driver of the seasonal dynamic thinning, as opposed to changes triggered by retreat. [ABSTRACT FROM AUTHOR]

The article offers information about the Late Paleocene age for Greenland's Hiawatha impact structure. It mentions that Hiawatha structure, located beneath Hiawatha Glacier in northwestern Greenland, has been proposed as an impact structure that may have formed after the Pleistocene inception of the Greenland Ice Sheet.

Abstract: Remains from Paleo-Eskimo cultures are well-documented, but complete preservation is rare. Two kitchen middens in Greenland are known to hold extremely well-preserved organic artefacts. Here, we assess the fate of the Qajaa site in Western Greenland under future climate conditions based on site characteristics measured in situ and from permafrost cores. Measurements of thermal properties, heat generation, oxygen consumption and CO2 production show that the kitchen midden can be characterized as peat but produces 4–7 times more heat than natural sediment. An analytical model from permafrost research has been applied to assess future thawing of the midden. Results show that the preservation conditions are controlled by freezing temperatures and a high water/ice content limiting the subsurface oxygen availability. Threats to the future preservation are related to thawing followed by drainage and increasing subsurface oxygen availability and heat generation. The model predicts that the unique 4000-year-old Saqqaq layer below more than 1 m of peat is adequately protected against thawing for the next 70 years. [Copyright &y& Elsevier]

In the Northern Hemisphere, an insolation driven Early to Middle Holocene Thermal Maximum was followed by a Neoglacial cooling that culminated during the Little Ice Age (LIA). Here, we review the glacier response to this Neoglacial cooling in Greenland. Changes in the ice margins of outlet glaciers from the Greenland Ice Sheet as well as local glaciers and ice caps are synthesized Greenland-wide. In addition, we compare temperature reconstructions from ice cores, elevation changes of the ice sheet across Greenland and oceanographic reconstructions from marine sediment cores over the past 5,000 years. The data are derived from a comprehensive review of the literature supplemented with unpublished reports. Our review provides a synthesis of the sensitivity of the Greenland ice margins and their variability, which is critical to understanding how Neoglacial glacier activity was interrupted by the current anthropogenic warming. We have reconstructed three distinct periods of glacier expansion from our compilation: two older Neoglacial advances at 2,500 – 1,700 yrs. BP (Before Present = 1950 CE, Common Era) and 1,250 – 950 yrs. BP; followed by a general advance during the younger Neoglacial between 700-50 yrs. BP, which represents the LIA. There is still insufficient data to outline the detailed spatio-temporal relationships between these periods of glacier expansion. Many glaciers advanced early in the Neoglacial and persisted in close proximity to their present-day position until the end of the LIA. Thus, the LIA response to Northern Hemisphere cooling must be seen within the wider context of the entire Neoglacial period of the past 5,000 years. Ice expansion appears to be closely linked to changes in ice sheet elevation, accumulation, and temperature as well as surface-water cooling in the surrounding oceans. At least for the two youngest Neoglacial advances, volcanic forcing triggering a sea-ice /ocean feedback, could explain their initiation. There are probably several LIA glacier fluctuations since the first culmination close to 1250 CE (Common Era) and available data suggests ice culminations in the 1400s, early to mid-1700s and early to mid-1800s CE. The last LIA maxima lasted until the present deglaciation commenced around 50 yrs. BP (1900 CE). The constraints provided here on the timing and magnitude of LIA glacier fluctuations delivers a more realistic background validation for modelling future ice sheet stability. [ABSTRACT FROM AUTHOR]

• A new model with a stochastic term was proposed to describe Greenland mass changes. • Interannual ice mass variations have strong impacts on secular trend estimations. • Conventional method estimates lower bound of acceleration and upper bound of rate. • New method improves the model parameter estimates, esp. the rate and acceleration. • Rate and acceleration estimates are improved with longer observation records. The irregular interannual variations observed in the Greenland ice sheet (GrIS) mass balance can be interpreted as stochastic. These variations often have large amplitudes, and, if not accounted for correctly in the mass change model parameterization, could have profound impacts on the estimate of the secular trend and acceleration. Here we propose a new mass trajectory model that includes both the conventional deterministic components and a stochastic component. This new model simultaneously estimates the secular rate and acceleration, seasonal components, and the stochastic component of mass changes. Simulations show that this new model improves estimates of model parameters, especially accelerations, over the conventional model without stochastic component. Using this new model, we estimate an acceleration of −1.6 ± 1.3 Gt/yr2 in mass change (minus means mass loss) for 2003-2017 using the Gravity Recovery and Climate Experiment (GRACE) data and an acceleration of −1.1 ± 1.3 Gt/yr2 using the modeled surface mass balance plus observed ice discharge. The corresponding rates are estimated to be −288.2 ± 12.7 Gt/yr and −274.9 ± 13.0 Gt/yr. The greatest discrepancies between the new and the conventional model parameter determinations are found in the acceleration estimates, −1.6 Gt/yr2 vs. −7.5 Gt/yr2 from the GRACE data. The estimated accelerations using the new method are apparently smaller than those estimated by other studies in terms of mass loss. Our quantitative analysis elucidates that the acceleration estimate using the conventional method is the lower bound (i.e., −7.5 Gt/yr2 for 2003–2017) while the acceleration estimated by the new method lies in the middle of the possible ranges. It is also found that these discrepancies between the new and the conventional methods diminish with sufficiently long (>20 yr) observation records. [ABSTRACT FROM AUTHOR]

Knowledge about the future response of the Greenland Ice Sheet to global climate change, including ice sheet contributions to sea level rise, is important for understanding the impact of climate change on society. Such studies rely in ice sheet model predictions and improved chronological constraints of past ice sheet extents and paleoclimatic trends. Many regions in Greenland are well studied, but northwest Greenland and especially Melville Bay, being one of the most important regions in terms of dynamical ice mass loss, lack a firm chronology of Holocene ice marginal fluctuations. In this study, we present the first comprehensive chronology for Melville Bay spanning 73.1–75.7°N based on 36 new 10Be exposure ages of boulders and 39 new radiocarbon ages of marine molluscs in Little Ice Age moraines. From weighted mean 10Be exposure ages, excluding 6 outliers, we find that the outer coast in Melville Bay was deglaciated ∼11.6 ± 0.3 ka (n = 15) and the ice margin reached its present-day position 40 km farther inland ∼11.5 ± 0.3 ka (n = 15). Our results suggest an interval of rapid ice-marginal retreat (i.e. collapse) of the northwest GrIS in Melville Bay, most likely triggered by rapidly rising atmospheric temperatures in early Holocene. Additionally, combining the comprehensive dataset of new radiocarbon ages with 26 radiocarbon ages from previous studies shows a restricted ice sheet extent from 9.1 ± 0.2 to 0.4 ± 0.1 cal ka BP, which coincides with increased sea surface temperatures. Our results highlight past ice sheet sensitivity towards climate changes in one of the least explored and most vulnerable regions of Greenland. Furthermore, comparing our new results to already existing ice sheet models (Huy3 and Huy3b) emphasize the proximal relevance of the Agassiz ice core temperature reconstruction for Melville Bay, which indicates the possible sensitivity of the ice sheet to a warming climate and place improved constraints on ice sheet simulations. • The northwest Greenland Ice Sheet was very sensitive to past climate changes. • The marine-based part of the ice sheet collapsed in Early Holocene. • The ice sheet was smaller than present throughout middle and late Holocene. • Rising atmospheric and sea surface temperatures drove deglaciation. • Ice sheet models do not fully capture this behaviour and improvements are needed. [ABSTRACT FROM AUTHOR]

The response of the Greenland Ice Sheet (GIS) to changes in temperature during the twentieth century remains contentious, largely owing to difficulties in estimating the spatial and temporal distribution of ice mass changes before 1992, when Greenland-wide observations first became available. The only previous estimates of change during the twentieth century are based on empirical modelling and energy balance modelling. Consequently, no observation-based estimates of the contribution from the GIS to the global-mean sea level budget before 1990 are included in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Here we calculate spatial ice mass loss around the entire GIS from 1900 to the present using aerial imagery from the 1980s. This allows accurate high-resolution mapping of geomorphic features related to the maximum extent of the GIS during the Little Ice Age at the end of the nineteenth century. We estimate the total ice mass loss and its spatial distribution for three periods: 1900-1983 (75.1 ± 29.4 gigatonnes per year), 1983-2003 (73.8 ± 40.5 gigatonnes per year), and 2003-2010 (186.4 ± 18.9 gigatonnes per year). Furthermore, using two surface mass balance models we partition the mass balance into a term for surface mass balance (that is, total precipitation minus total sublimation minus runoff) and a dynamic term. We find that many areas currently undergoing change are identical to those that experienced considerable thinning throughout the twentieth century. We also reveal that the surface mass balance term shows a considerable decrease since 2003, whereas the dynamic term is constant over the past 110 years. Overall, our observation-based findings show that during the twentieth century the GIS contributed at least 25.0 ± 9.4 millimetres of global-mean sea level rise. Our result will help to close the twentieth-century sea level budget, which remains crucial for evaluating the reliability of models used to predict global sea level rise.