Search

A long-term time series of ice sheet surface elevation change (SEC) is an essential parameter to assess the impact of climate change. In this study, we used an updated plane-fitting least-squares regression strategy to generate a 30-year surface elevation time series for the Greenland Ice Sheet (GrIS) at monthly temporal resolution and 5×5 km grid spatial resolution using ERS-1 (European Remote Sensing), ERS-2, Envisat, and CryoSat-2 satellite radar altimeter observations obtained between August 1991 and December 2020. The ingenious corrections for intermission bias were applied using an updated plane-fitting least-squares regression strategy. Empirical orthogonal function (EOF) reconstruction was used to supplement the sparse monthly gridded data attributable to poor observations in the early years. Validation using both airborne laser altimeter observations and the European Space Agency GrIS Climate Change Initiative (CCI) product indicated that our merged surface elevation time series is reliable. The accuracy and dispersion of errors of SECs of our results were 19.3 % and 8.9 % higher, respectively, than those of CCI SECs and even 30.9 % and 19.0 % higher, respectively, in periods from 2006–2010 to 2010–2014. Further analysis showed that our merged time series could provide detailed insight into GrIS SEC on multiple temporal (up to 30 years) and spatial scales, thereby providing an opportunity to explore potential associations between ice sheet change and climatic forcing. The merged surface elevation time series data are available at https://doi.org/10.11888/Glacio.tpdc.271658 (Zhang et al., 2021).

The Greenland Ice Sheet (GrIS) mass balance shows significant variabilities over a range of time scales. As geodetic records lengthen over time, it becomes insufficient to characterize the temporal evolution of the ice sheet by using a best-fit linear trend over a given observation period. This study investigates the joint analysis of laser and radar satellite altimeter measurements for estimating GrIS surface elevation changes (SECs) with a 30-day resolution. We first apply a crossover analysis to assess the precisions of the surface elevations measured by ICESat-1/2 laser altimeters and CryoSat-2 radar altimeter over the GrIS, which are needed for assigning weights for each data set in the joint analysis. Then, based on a modified repeat-track approach, we analyze the surface elevation measurements of ICESat-1/2 and CryoSat-2 to produce monthly SEC estimates for the past two decades, together with their associated uncertainties. The multimission SEC estimates are further assessed by using IceBridge airborne laser measurements, showing differences with a median value of -12 cm ± 60 cm. The monthly SEC time series reveal important variations over a range of time scales across different parts of the GrIS and would facilitate the investigation of complex spatiotemporal patterns of GrIS changes.

Ice sheet subsurface temperature is important for understanding glacier dynamics, yet existing methods to obtain the temperature of the ice sheet column are limited to in situ sources at present. The ultrawideband software-defined microwave radiometer (UWBRAD) has been developed to investigate the remote sensing of ice sheet internal temperatures. UWBRAD measures brightness temperature spectra from 0.5 to 2 GHz using 12 subchannels and employs a sophisticated algorithm for detection and mitigation of radio frequency interference (RFI). The instrument was deployed during a flight over northwestern Greenland in September 2017 and acquired the first wideband low-frequency brightness temperature spectra over the ice sheet and coastal regions. The results reveal strong spatial and spectral variations that correlate well with internal ice sheet temperature information. In this article, the section of the flight path ranging from the Camp Century to NEEM to NGRIP boreholes is used for subsurface temperature estimation. A ``partially coherent'' forward model is applied along with a Robin model for the temperature profile and a two-scale model of ice sheet density variations to describe measured brightness temperatures. Using this model, vertical temperature profiles are retrieved along the flight path using a sequential Bayesian estimator; borehole measurements at the three campsites are used to obtain Bayesian priors. The retrieved temperature profiles show reasonable behaviors and demonstrate the potential of ultrawideband microwave radiometry for remotely sensing internal ice sheet temperatures.

Surface melting of the Greenland Ice Sheet contributes a large amount to current and future sea-level rise. Increased surface melt, algae growth, debris, and dust deposition lower the reflectivity of the ice surface and thereby increase melt rates: the so-called melt-albedo feedback describes this potentially self-sustaining increase in surface melting. Here we present a simplified version of the diurnal Energy Balance Model (dEBM-simple) which is implemented as a surface melt module in the Parallel Ice Sheet Model (PISM). dEBM-simple is a modification of diurnal Energy Balance Model (dEBM), a surface melt scheme of intermediate complexity useful for simulations over centennial to multi-millennial timescales. dEBM-simple is computationally efficient, suitable for standalone ice-sheet modeling and includes a simple representation of the melt-albedo feedback. Using dEBM-simple and PISM, we find that this feedback increases ice loss until 2300 through surface warming by 60 % for the high-emission scenario RCP8.5. With an increase of 90 %, the effect is more pronounced for lower surface warming under RCP2.6. Furthermore, assuming an immediate darkening of the ice surface over all summer months, we estimate an upper bound for this effect to be +70 % in the RCP8.5 scenario and a more than fourfold increase under RCP2.6. With dEBM-simple implemented in PISM, we find that the melt-albedo feedback is an essential contributor to mass loss in dynamic simulations of the Greenland Ice Sheet under future warming.

Ice cores from polar ice sheets and glaciers are an important climate archive. Snow layers, consecutively deposited and buried, contain climatic information from the time of their formation. However, particularly low-accumulation areas are characterised by temporally intermittent precipitation, which can be further redistributed after initial deposition, depending on the local surface features at different spatial scales. Therefore, the accumulation conditions at an ice core site influence the quantity and quality of the recorded climate signal in proxy records. This study aims to characterise the local accumulation patterns and the evolution of the snow height to describe the contribution of the snow (re-)deposition to the overall noise level in climate records from ice cores. To this end, we applied a structure-from-motion photogrammetry approach to generate near-daily elevation models of the surface snow for a 195 m2 area in the vicinity of the deep drilling site of the East Greenland Ice-core Project in northeast Greenland. Based on the snow height information we derive snow height changes on a day-to-day basis throughout our observation period from May to August 2018 and find an average snow height increase of ∼ 11 cm. The spatial and temporal data set also allows an investigation of snow deposition versus depositional modifications. We observe irregular snow deposition and erosion causing uneven snow accumulation patterns, a removal of more than 60 % of the deposited snow, and a negative relationship between the initial snow height and the amount of accumulated snow. Furthermore, the surface roughness decreased by approximately a factor of 2 throughout the spring and summer season at our study site. Finally, our study shows that structure from motion is a relatively simple method to demonstrate the potential influences of depositional processes on proxy signals in snow and ice.

Greenland Ice Sheet mass loss is impacting connected terrestrial and marine hydrologic systems with global consequences. Groundwater is a key component of water cycling in the Arctic, underlying the 1.7e6 km2 ice sheet and forming offshore freshwater reserves. However, despite its vast extent, the response of Greenland’s groundwater to ongoing ice sheet change is unknown. Here we present in-situ observations of deep groundwater conditions under the Greenland Ice Sheet, obtained in a 651-metre-long proglacial bedrock borehole angled under the ice sheet margin. We find that Greenland’s groundwater system responds rapidly and sensitively to relatively minor ice sheet forcing. Hydraulic head clearly varies over multi-annual, seasonal and diurnal timescales, which we interpret as a response to fluid pressure forcing at the ice/bed interface associated with changes in overlying ice loading and ice sheet hydrology. We find a systematic decline in hydraulic head over the eight-year observational period is linked primarily to ice sheet mass loss. Ongoing and future ice thinning will probably reduce groundwater discharge rates, with potential impacts to submarine freshwater discharge, freshwater delivery to fjords and biogeochemical fluxes in the Arctic. Greenland’s groundwater system responds rapidly to ice-sheet change, according to borehole observations from underneath the ice-sheet margin.

The physical interactions between ice sheets and their surroundings are major factors in determining the state of the climate system, yet many current Earth system models omit them entirely or approximate them in a heavily parameterised manner. In this work we have improved the snow and ice sheet surface physics in the FAMOUS climate model, with the aim of improving the representation of polar climate and implementing a bidirectional coupling to the Glimmer dynamic ice sheet model using the water and energy fluxes calculated by FAMOUS. FAMOUS and Glimmer are both low-resolution, computationally affordable models used for multi-millennial simulations. Glaciated surfaces in the new FAMOUS-ice are modelled using a multi-layer snow scheme capable of simulating compaction of firn and the percolation and refreezing of surface melt. The low horizontal resolution of FAMOUS compared to Glimmer is mitigated by implementing this snow model on sub-grid-scale tiles that represent different elevations on the ice sheet within each FAMOUS grid box. We show that with this approach FAMOUS-ice can simulate relevant physical processes on the surface of the modern Greenland ice sheet well compared to higher-resolution climate models and that the ice sheet state in the coupled FAMOUS-ice–Glimmer system does not drift unacceptably. FAMOUS-ice coupled to Glimmer is thus a useful tool for modelling the physics and co-evolution of climate and grounded ice sheets on centennial and millennial timescales, with applications to scientific questions relevant to both paleoclimate and future sea level rise.

Radar reflections from the interior of the Greenland ice sheet contain a comprehensive archive of past accumulation rates, ice dynamics, and basal melting. Combining these data with dynamic ice sheet models may greatly aid model calibration, improve past and future sea level estimates, and enable insights into past ice sheet dynamics that neither models nor data could achieve alone. Unfortunately, simulating the continental-scale ice sheet stratigraphy represents a major challenge for current ice sheet models. In this study, we present the first three-dimensional ice sheet model that explicitly simulates the Greenland englacial stratigraphy. Individual layers of accumulation are represented on a grid whose vertical axis is time so that they do not exchange mass with each other as the flow of ice deforms them. This isochronal advection scheme does not influence the ice dynamics and only requires modest input data from a host thermomechanical ice sheet model, making it easy to transfer to a range of models. Using an ensemble of simulations, we show that direct comparison with the dated radiostratigraphy data yields notably more accurate results than calibrating simulations based on total ice thickness. We show that the isochronal scheme produces a more reliable simulation of the englacial age profile than traditional age tracers. The interpretation of ice dynamics at different times is possible but limited by uncertainties in the upper and lower boundary conditions, namely temporal variations in surface mass balance and basal friction.

The Programme for Monitoring of the Greenland Ice Sheet (PROMICE) has been measuring climate and ice sheet properties since 2007. Currently, the PROMICE automatic weather station network includes 25 instrumented sites in Greenland. Accurate measurements of the surface and near-surface atmospheric conditions in a changing climate are important for reliable present and future assessment of changes in the Greenland Ice Sheet. Here, we present the PROMICE vision, methodology, and each link in the production chain for obtaining and sharing quality-checked data. In this paper, we mainly focus on the critical components for calculating the surface energy balance and surface mass balance. A user-contributable dynamic web-based database of known data quality issues is associated with the data products at https://github.com/GEUS-Glaciology-and-Climate/PROMICE-AWS-data-issues/ (last access: 7 April 2021). As part of the living data option, the datasets presented and described here are available at https://doi.org/10.22008/promice/data/aws (Fausto et al., 2019).

The Nioghalvfjerdsfjorden glacier (also known as the 79∘ North Glacier) drains approximately 8 % of the Greenland Ice Sheet. Supraglacial lakes (SGLs), or surface melt ponds, are a persistent summertime feature and are thought to drain rapidly to the base of the glacier and influence seasonal ice velocity. However, seasonal development and spatial distribution of SGLs in the north-east of Greenland are poorly understood, leaving a substantial error in the estimate of meltwater and its impacts on ice velocity. Using results from an automated detection of melt ponds, atmospheric and surface mass balance modelling, and reanalysis products, we investigate the role of specific climatic conditions in melt onset, extent, and duration from 2016 to 2019. The summers of 2016 and 2019 were characterised by above-average air temperatures, particularly in June, as well as a number of rainfall events, which led to extensive melt ponds to elevations up to 1600 m. Conversely, 2018 was particularly cold, with a large accumulated snowpack, which limited the development of lakes to altitudes less than 800 m. There is evidence of inland expansion and increases in the total area of lakes compared to the early 2000s, as projected by future global warming scenarios.

During the last glacial period (c. 120–11 kyr BP), dramatic temperature swings, known as Dansgaard-Oeschger (D-O) events, are clearly manifest in high resolution oxygen isotope records from the Greenland ice sheet. Although variability in the Atlantic Meridional Overturning Circulation (AMOC) is often invoked, a unified explanation for what caused these "sawtooth shaped" climate patterns has yet to be accepted. Of particular interest is the most recent D-O shaped climate pattern that occurred from ∼ 14,600 to 11,500 years ago – the Bølling/Allerød (BA) warm interstadial and the subsequent Younger Dryas (YD) cold stadial. Unlike earlier D-O stadials, the YD is frequently considered a unique event, potentially resulting from a rerouting and/or flood of glacial meltwater into the North Atlantic, a meteorite impact, or a volcanic eruption. Yet, these mechanisms are seldom considered as the cause of the earlier stadials. Using a robust multivariate outlier detection scheme – a novel approach for traditional paleoclimate research – we show that the pattern of climate change during the BA/YD is not statistically different from the other D-O events in the Greenland record, and that it should not be considered unique when investigating the drivers of abrupt climate change. Our results thus raise important questions about the ability of glacial meltwater input and other "one off" events to trigger abrupt, centennial-to-millennial length, changes in climate.

The aerodynamic roughness of heat, moisture, and momentum of a natural surface are important parameters in atmospheric models, as they co-determine the intensity of turbulent transfer between the atmosphere and the surface. Unfortunately this parameter is often poorly known, especially in remote areas where neither high-resolution elevation models nor eddy-covariance measurements are available. In this study we adapt a bulk drag partitioning model to estimate the aerodynamic roughness length (z0m) such that it can be applied to 1D (i.e. unidirectional) elevation profiles, typically measured by laser altimeters. We apply the model to a rough ice surface on the K-transect (west Greenland Ice Sheet) using UAV photogrammetry, and we evaluate the modelled roughness against in situ eddy-covariance observations. We then present a method to estimate the topography at 1 m horizontal resolution using the ICESat-2 satellite laser altimeter, and we demonstrate the high precision of the satellite elevation profiles against UAV photogrammetry. The currently available satellite profiles are used to map the aerodynamic roughness during different time periods along the K-transect, that is compared to an extensive dataset of in situ observations. We find a considerable spatio-temporal variability in z0m, ranging between 10−4 m for a smooth snow surface and 10−1 m for rough crevassed areas, which confirms the need to incorporate a variable aerodynamic roughness in atmospheric models over ice sheets.

The Greenland Ice Sheet (GrIS) rapid mass loss is primarily driven by an increase in meltwater runoff, which highlights the importance of understanding the formation, evolution, and impact of meltwater features on the ice sheet. Buried lakes are meltwater features that contain liquid water and exist under layers of snow, firn, and/or ice. These lakes are invisible in optical imagery, challenging the analysis of their evolution and implication for larger GrIS dynamics and mass change. Here, we present a method that uses a convolutional neural network, a deep learning method, to automatically detect buried lakes across the GrIS. For the years 2018 and 2019 (which represent low- and high-melt years, respectively), we compare total areal extent of both buried and surface lakes across six regions, and we use a regional climate model to explain the spatial and temporal differences. We find that the total buried lake extent after the 2019 melt season is 56 % larger than after the 2018 melt season across the entire ice sheet. Northern Greenland has the largest increase in buried lake extent after the 2019 melt season, which we attribute to late-summer surface melt and high autumn temperatures. We also provide evidence that different processes are responsible for buried lake formation in different regions of the ice sheet. For example, in southwest Greenland, buried lakes often appear on the surface during the previous melt season, indicating that these meltwater features form when surface lakes partially freeze and become insulated as snowfall buries them. Conversely, in southeast Greenland, most buried lakes never appear on the surface, indicating that these features may form due to downward percolation of meltwater and/or subsurface penetration of shortwave radiation. We provide support for these processes via the use of a physics-based snow model. This study provides additional perspective on the potential role of meltwater on GrIS dynamics and mass loss.

Higher boreal summer insolation in the early to middle Holocene drove thousands of years of summer warming across the Arctic. Modern-day warming has distinctly different causes, but geologic data from this past warm period hold lessons for the future. We compile Holocene temperature reconstructions from ice, lake, and marine cores around Greenland, where summer temperatures are globally important due to their influence on ice sheet mass balance, ocean circulation, and sea ice. Highlighting and accounting for some key issues with proxy interpretation, we find that much of Greenland experienced summers 3 to 5°C warmer than the mid-twentieth century in the early Holocene—earlier and stronger warming than often presumed. Warmth had dramatic consequences: Many glaciers disappeared, perennial sea ice retreated, plants and animals migrated northward, the Greenland Ice Sheet shrank rapidly, and increased meltwater discharge led to strong marine water stratification and enhanced winter sea ice in some areas. ▪ Summer air temperatures and open ocean temperatures around much of Greenland peaked in the early Holocene in response to elevated summer insolation. ▪ Peak summer air temperatures ranged from 3 to 5°C warmer than the mid-twentieth century in northwest and central Greenland to perhaps 1 to 2°C warmer in south Greenland. ▪ Many differences between records can be explained by proxy seasonality, ice sheet elevation changes, vegetation analogs and lags, and the nearshore effects of ice sheet meltwater. ▪ Early Holocene warmth dramatically affected glaciers and the Greenland Ice Sheet; meltwater discharge, nearshore ocean salinity, and sea ice; and diverse flora and fauna.

The elevation of ice sheets changes due to climate change, and satellite altimetry is the preferred tool for measuring ice sheet-wide height changes. In situ validation is needed to ensure the quality of the observed elevation changes, but the cost often limits the amount of in situ data which can be collected. As more tourists are accessing the ice sheets, citizen science might provide in situ data in an environmentally friendly and cost-efficient way. Here, we investigate the opportunistic kinematic global positioning system (GPS) profiles across the Greenland ice sheet, collected by the American-Icelandic expedition on the Greenlandic icecap 2018. The collected GPS data are in good agreement with the widely used NASA’s Operation IceBridge Airborne LiDAR data measured within ± 10 days, with an average difference of 10.7 cm ± 11.7 cm. The main difference is attributed to changes in the compaction of the snow while driving and changes in the tires’ pressure. The kinematic GPS data are then compared with data from the European Space Agency’s CryoSat-2 mission. Here, an average bias of 92.3 cm ± 65.7 cm in the two records is observed between the spring CryoSat-2 and the truck GPS data obtained largely in the dry-snow zone. This suggests that the surface penetration of Ku-band radar on the Greenland ice sheet and the observed magnitude are consistent with the literature. Finally, we compared the 2018 GPS data to a profile obtained in 2005 near Kangerlussuaq, West Greenland. Here, the records show an average ice-elevation decrease of 9 m, with peaks at 26 m. These results show that the citizen science kinematic GPS data can provide high-resolution data necessary for the validation of satellite altimetry, with the added benefit of potential direct sampling properties of the surface and firn. Linking up with citizen-science expeditions is a beneficial way of providing cost-effective satellitevalidations and may also have a societal impact by involving more people in the climate monitoring of ice sheets.

Physical geographer Konrad “Koni” Steffen, lost 8 August 2020 in a crevasse on the Greenland ice sheet, was a pioneer in satellite remote sensing and field observations of the Greenland ice sheet. This Classics Revisited piece honors the memory of Koni Steffen and examines the impact of a work which laid the foundation for numerous studies that made the Greenland ice sheet and the man global icons of climate change.

Radiative transfer in snow and ice is often not modeled explicitly in regional climate models. In this study, we evaluate a new englacial radiative transfer scheme and assess the surface mass and energy budget for the Greenland ice sheet in the latest version of the regional climate model RACMO2, version 2.3p3. We also evaluate the modeled (sub)surface temperature and melt, as radiation penetration now enables internal heating. The results are compared to the previous model version and are evaluated against stake measurements and automatic weather station data of the K-transect and PROMICE projects. In addition, subsurface snow temperature profiles are compared at the K-transect, Summit, and southeast Greenland. The surface mass balance is in good agreement with observations, with a mean bias of −31 mm w.e. yr−1 (−2.67 %), and only changes considerably with respect to the previous RACMO2 version around the ice margins and near the percolation zone. Melt and refreezing, on the other hand, are changed more substantially in various regions due to the changed albedo representation, subsurface energy absorption, and meltwater percolation. Internal heating leads to higher snow temperatures in summer, in agreement with observations, and introduces a shallow layer of subsurface melt. Hence, this study shows the consequences and necessity of radiative transfer in snow and ice for regional climate modeling of the Greenland ice sheet.

Sea level contributions from the Greenland Ice Sheet are influenced by the rapid changes in glacial terminus positions. The documentation of these evolving calving front positions, for which satellite imagery forms the basis, is therefore important. However, the manual delineation of these calving fronts is time consuming, which limits the availability of these data across a wide spatial and temporal range. Automated methods face challenges that include the handling of clouds, illumination differences, sea ice mélange, and Landsat 7 scan line corrector errors. To address these needs, we develop the Calving Front Machine (CALFIN), an automated method for extracting calving fronts from satellite images of marine-terminating glaciers, using neural networks. The results are often indistinguishable from manually curated fronts, deviating by on average 86.76 ± 1.43 m from the measured front. Landsat imagery from 1972 to 2019 is used to generate 22 678 calving front lines across 66 Greenlandic glaciers. This improves on the state of the art in terms of the spatiotemporal coverage and accuracy of its outputs and is validated through a comprehensive intercomparison with existing studies. The current implementation offers a new opportunity to explore subseasonal and regional trends on the extent of Greenland's margins and supplies new constraints for simulations of the evolution of the mass balance of the Greenland Ice Sheet and its contributions to future sea level rise.

Sometime during the middle to late Holocene (8.2 ka to ∼ 1850–1900 CE), the Greenland Ice Sheet (GrIS) was smaller than its current configuration. Determining the exact dimensions of the Holocene ice-sheet minimum and the duration that the ice margin rested inboard of its current position remains challenging. Contemporary retreat of the GrIS from its historical maximum extent in southwestern Greenland is exposing a landscape that holds clues regarding the configuration and timing of past ice-sheet minima. To quantify the duration of the time the GrIS margin was near its modern extent we develop a new technique for Greenland that utilizes in situ cosmogenic 10Be – 14C – 26Al in bedrock samples that have become ice-free only in the last few decades due to the retreating ice-sheet margin at Kangiata Nunaata Sermia ( n=12 sites, 36 measurements; KNS), southwest Greenland. To maximize the utility of this approach, we refine the deglaciation history of the region with stand-alone 10Be measurements ( n=49 ) and traditional 14C ages from sedimentary deposits contained in proglacial–threshold lakes. We combine our reconstructed ice-margin history in the KNS region with additional geologic records from southwestern Greenland and recent model simulations of GrIS change to constrain the timing of the GrIS minimum in southwest Greenland and the magnitude of Holocene inland GrIS retreat, as well as to explore the regional climate history influencing Holocene ice-sheet behavior. Our 10Be – 14C – 26Al measurements reveal that (1) KNS retreated behind its modern margin just before 10 ka , but it likely stabilized near the present GrIS margin for several thousand years before retreating farther inland, and (2) pre-Holocene 10Be detected in several of our sample sites is most easily explained by several thousand years of surface exposure during the last interglaciation. Moreover, our new results indicate that the minimum extent of the GrIS likely occurred after ∼5 ka , and the GrIS margin may have approached its eventual historical maximum extent as early as ∼2 ka . Recent simulations of GrIS change are able to match the geologic record of ice-sheet change in regions dominated by surface mass balance, but they produce a poorer model–data fit in areas influenced by oceanic and dynamic processes. Simulations that achieve the best model–data fit suggest that inland retreat of the ice margin driven by early to middle Holocene warmth may have been mitigated by increased precipitation. Triple 10Be – 14C – 26Al measurements in recently deglaciated bedrock provide a new tool to help decipher the duration of smaller-than-present ice over multiple timescales. Modern retreat of the GrIS margin in southwest Greenland is revealing a bedrock landscape that was also exposed during the migration of the GrIS margin towards its Holocene minimum extent, but it has yet to tap into a landscape that remained ice-covered throughout the entire Holocene.

Temperature sensors installed in a grid of nine full-depth boreholes drilled in the southwestern ablation zone of the Greenland Ice Sheet recorded cooling in discrete sections of ice over time within the lowest third of the ice column in most boreholes. Rates of temperature change outpace cooling expected from vertical conduction alone. Additionally, observed temperature profiles deviate significantly from the site-average thermal profile that is shaped by all thermomechanical processes upstream. These deviations imply recent, localized changes to the basal thermal state in the boreholes. Although numerous heat sources exist to add energy and warm ice as it moves from the central divide towards the margin such as strain heat from internal deformation, latent heat from refreezing meltwater, and the conduction of geothermal heat across the ice–bedrock interface, identifying heat sinks proves more difficult. After eliminating possible mechanisms that could cause cooling, we find that the observed cooling is a manifestation of previous warming in near-basal ice. Thermal decay after latent heat is released from freezing water in basal crevasses is the most likely mechanism resulting in the transient evolution of temperature and the vertical thermal structure observed at our site. We argue basal crevasses are a viable englacial heat source in the basal ice of Greenland's ablation zone and may have a controlling influence on the temperature structure of the near-basal ice.

High-resolution shallow seismic information supported by sedimentological and geochemical data from sediment cores has been used for reconstruction of Holocene sea-level changes along the northern coast of the Persian Gulf. These investigations in the area of the North Qeshm Island Waterway (NQIW) near the Strait of Hormuz reveal a series of five prominent sea-level fluctuations starting with the second phase of regional postglacial sea-level rise recorded between 15.0 ka and 9.5 ka (ka = thousands of years before present). Based on our study, sea-level had reached −22 m before 9.1 ka. Postglacial sea level high-stand in the study area has been dated to have occurred around 6.0 ka, which was followed by a relative sea-level fall coeval with a regional climate change from humid to arid conditions at about 5.4 ka. After a renewed, but relatively short period of sea level rise associated with climate changes around 4.35 ka, sea-level fell again at about 3.3 ka and remained relatively stable from 3.2 ka until the present. We conclude that major sea-level changes in the Persian Gulf coincided with sea-level fluctuations in the Indian Ocean and observed changes in North Atlantic Ocean and atmospheric circulation and Greenland Ice Sheet melting. Thus, despite of regional tectonic instability, regional sea-level stand in the Persian Gulf has been controlled mainly by eustatic sea-level changes of the global ocean.

We combine satellite observations and numerical models to show that Earth lost 28 trillion tonnes of ice between 1994 and 2017. Arctic sea ice (7.6 trillion tonnes), Antarctic ice shelves (6.5 trillion tonnes), mountain glaciers (6.1 trillion tonnes), the Greenland ice sheet (3.8 trillion tonnes), the Antarctic ice sheet (2.5 trillion tonnes), and Southern Ocean sea ice (0.9 trillion tonnes) have all decreased in mass. Just over half (58 %) of the ice loss was from the Northern Hemisphere, and the remainder (42 %) was from the Southern Hemisphere. The rate of ice loss has risen by 57 % since the 1990s – from 0.8 to 1.2 trillion tonnes per year – owing to increased losses from mountain glaciers, Antarctica, Greenland and from Antarctic ice shelves. During the same period, the loss of grounded ice from the Antarctic and Greenland ice sheets and mountain glaciers raised the global sea level by 34.6 ± 3.1 mm. The majority of all ice losses were driven by atmospheric melting (68 % from Arctic sea ice, mountain glaciers ice shelf calving and ice sheet surface mass balance), with the remaining losses (32 % from ice sheet discharge and ice shelf thinning) being driven by oceanic melting. Altogether, these elements of the cryosphere have taken up 3.2 % of the global energy imbalance.

Majority areas of Antarctica and Greenland are under the thick ice sheet and characterized by evolving cryosphere surroundings. In the polar region, associated with the recent trend on climate change such as global warming, glacier relating earthquakes are increasing during this 21st century. In this paper, a decade of progress in “Cryoseismology” at bi-polar regions is re-viewed by focusing on the contribution from Japanese researchers. In particular, the specific cryoseismic events are treated, which occurred in the coastal area of East Antarctica, around the Lutzow-Holm Bay, together with the coast and whole inland area of Greenland. As the major scientific results, frequency-overtone signals in the harmonic cryoseismic tremors were analyzed by assuming constant sources, suggesting inter-glacial asperities that generate characteristic tremors. Infrasound source locations were also determined by using the array deployment at the coastal regions in the Antarctic. In contrast, characteristics of glacial earthquakes and seismic interferometry approach have been conducted so as to achieve the fine structure of the Greenland Ice Sheet (GrIS) in particular the basal condition beneath the ice sheet.

Changes in the Greenland ice sheet (GIS) affect global sea level. Greenland stable water isotope (δ18O) records from ice cores offer information on past changes in the surface of the GIS. Here, we use the isotope-enabled Hadley Centre Coupled Model version 3 (HadCM3) climate model to simulate a set of last interglacial (LIG) idealised GIS surface elevation change scenarios focusing on GIS ice core sites. We investigate how δ18O depends on the magnitude and sign of GIS elevation change and evaluate how the response is altered by sea ice changes. We find that modifying GIS elevation induces changes in Northern Hemisphere atmospheric circulation, sea ice and precipitation patterns. These climate feedbacks lead to ice-core-averaged isotopic lapse rates of 0.49 ‰ (100 m)−1 for the lowered GIS states and 0.29 ‰ (100 m)−1 for the enlarged GIS states. This is lower than the spatially derived Greenland lapse rates of 0.62–0.72 ‰ (100 m)−1. These results thus suggest non-linearities in the isotope–elevation relationship and have consequences for the interpretation of past elevation and climate changes across Greenland. In particular, our results suggest that winter sea ice changes may significantly influence isotope–elevation gradients: winter sea ice effect can decrease (increase) modelled core-averaged isotopic lapse rate values by about −19 % (and +28 %) for the lowered (enlarged) GIS states, respectively. The largest influence of sea ice on δ18O changes is found in coastal regions like the Camp Century site.

We present continuous estimates of snow and firn density, layer depth and accumulation from a multi-channel, multi-offset, ground-penetrating radar traverse. Our method uses the electromagnetic velocity, estimated from waveform travel-times measured at common-midpoints between sources and receivers. Previously, common-midpoint radar experiments on ice sheets have been limited to point observations. We completed radar velocity analysis in the upper ~2 m to estimate the surface and average snow density of the Greenland Ice Sheet. We parameterized the Herron and Langway (1980) firn density and age model using the radar-derived snow density, radar-derived surface mass balance (2015–2017) and reanalysis-derived temperature data. We applied structure-oriented filtering to the radar image along constant age horizons and increased the depth at which horizons could be reliably interpreted. We reconstructed the historical instantaneous surface mass balance, which we averaged into annual and multidecadal products along a 78 km traverse for the period 1984–2017. We found good agreement between our physically constrained parameterization and a firn core collected from the dry snow accumulation zone, and gained insights into the spatial correlation of surface snow density.

The presence of a grounded Greenland Ice Sheet on the northeastern part of the Greenland continental shelf during the Last Glacial Maximum is supported by new swath bathymetry and high-resolution seismic data, supplemented with multi-proxy analyses of sediment gravity cores from Store Koldewey Trough. Subglacial till fills the trough, with an overlying drape of maximum 2.5 m thick glacier-proximal and glacier-distal sediment. The presence of mega-scale glacial lineations and a grounding zone wedge in the outer part of the trough, comprising subglacial till, provides evidence of the expansion of fast-flowing, grounded ice, probably originating from the area presently covered with the Storstrømmen ice stream and thereby previously flowing across Store Koldewey Island and Germania Land. Grounding zone wedges and recessional moraines provide evidence that multiple halts and/or readvances interrupted the deglaciation. The formation of the grounding zone wedges is estimated to be at least 130 years, while distances between the recessional moraines indicate that the grounding line locally retreated between 80 and 400 m yr−1 during the deglaciation, assuming that the moraines formed annually. The complex geomorphology in Store Koldewey Trough is attributed to the trough shallowing and narrowing towards the coast. At a late stage of the deglaciation, the ice stream flowed around the topography on Store Koldewey Island and Germania Land, terminating the sediment input from this sector of the Greenland Ice Sheet to Store Koldewey Trough.

Determining the sensitivity of the Greenland Ice Sheet (GrIS) to Holocene climate changes is a key prerequisite for understanding the future response of the ice sheet to global warming. In this study, we present new information on the Holocene glacial history of the GrIS in Inglefield Land, north Greenland. We use 10 Be and in situ 14 C exposure dating to constrain the timing of deglaciation in the area and radiocarbon dating of reworked molluscs and wood fragments to constrain when the ice sheet retreated behind its present-day extent. The 10 Be ages are scattered ranging from ca. 92.7 to 6.8 ka, whereas the in situ 14 C ages range from ca. 14.2 to 6.7 ka. Almost half of the apparent 10 Be ages predate the Last Glacial Maximum and up to 89 % are to some degree affected by nuclide inheritance. Based on the few reliable 10 Be ages, the in situ 14 C ages and existing radiocarbon ages from Inglefield Land, we find that the deglaciation along the coast commenced at ca. 8.6–8.3 ka cal BP in the western part and ca. 7.9 ka in the central part, following the opening of Nares Strait and arrival of warm waters. The ice margin reached its present-day position at ca. 8.2 ka at the Humboldt Glacier and ca. 6.7 ka in the central part of Inglefield Land. Radiocarbon ages of reworked molluscs and wood fragments show that the ice margin was behind its present-day extent from ca. 5.8 to 0.5 ka cal BP. After 0.5 ka cal BP, the ice advanced towards its Little Ice Age position. Our results emphasize that the slowly eroding and possibly cold-based ice in north Greenland makes it difficult to constrain the deglaciation history based on 10 Be ages alone unless they are paired with in situ 14 C ages. Further, combining our findings with those of recently published studies reveals distinct differences between deglaciation patterns of northwest and north Greenland. Deglaciation of the land areas in northwest Greenland occurred earlier than in north Greenland, and periods of restricted ice extent were longer, spanning the Middle and Late Holocene. Overall, this highlights past ice sheet sensitivity to Holocene climate changes in an area where little information was available just a few years ago.

Close-range (CR) L-band radiometry and quasi-simultaneous in situ snow characterizations were conducted in May 2019 at the Swiss Camp research site in the ablation zone of the western Greenland ice sheet. Snow liquid-water and its melt/refreeze are retrieved from microwave antenna temperatures measured with the ground-based L-band radiometer ELBARA-III. The emission model (EM) used in the retrieval is a two-layer configuration of the ‘L-Band Specific Microwave Emission Model of Layered Snowpack’. Consistent snow wetness retrievals were achieved from both single- and multi-angle CR observations of L-band antenna temperatures. This suggests that multi-angle observation is not a pre-requisite for snow wetness retrieval. Therefore, in addition to soil moisture and ocean salinity (SMOS) multi-angle measurements, snow wetness can be estimated from spaceborne L-band brightness temperatures measured at a single observation angle, such as from NASA's SMAP satellite. Our results provide partial validation of a recently presented snow wetness retrieval approach based on the same EM and applied over Greenland using multi-angle SMOS brightness temperatures. Agreement between measured CR antenna temperatures and SMOS brightness temperatures is found to be within the 95% confidence intervals of ELBARA-III and SMOS measurement uncertainties. Our measurements confirm the modeled response of antenna temperatures to diurnal variations of snow wetness.

Projections of the contribution of the Greenland ice sheet to future sea-level rise include uncertainties primarily due to the imposed climate forcing and the initial state of the ice sheet model. Several state-of-the-art ice flow models are currently being employed on various grid resolutions to estimate future mass changes in the framework of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). Here we investigate the sensitivity to grid resolution of centennial sea-level contributions from the Greenland ice sheet and study the mechanism at play.We employ the finiteelement higher-order Ice-sheet and Sea-level System Model (ISSM) and conduct experiments with four different horizontal resolutions, namely 4, 2, 1 and 0.75 km. We run the simulation based on the ISMIP6 core climate forcing from the MIROC5 global circulation model (GCM) under the highemission Representative Concentration Pathway (RCP) 8.5 scenario and consider both atmospheric and oceanic forcing in full and separate scenarios. Under the full scenarios, finer simulations unveil up to approximately 5% more sea-level rise compared to the coarser resolution. The sensitivity depends on the magnitude of outlet glacier retreat, which is implemented as a series of retreat masks following the ISMIP6 protocol. Without imposed retreat under atmosphereonly forcing, the resolution dependency exhibits an opposite behaviour with approximately 5% more sea-level contribution in the coarser resolution. The sea-level contribution indicates a converging behaviour below a 1 km horizontal resolution. A driving mechanism for differences is the ability to resolve the bedrock topography, which highly controls ice discharge to the ocean. Additionally, thinning and acceleration emerge to propagate further inland in high resolution for many glaciers. A major response mechanism is sliding, with an enhanced feedback on the effective normal pressure at higher resolution leading to a larger increase in sliding speeds under scenarios with outlet glacier retreat., SCOPUS: ar.j, info:eu-repo/semantics/published

We used a deep-ultraviolet fluorescence mapping spectrometer, coupled to a drill system, to scan from the surface to 105 m depth into the Greenland ice sheet. The scan included firn and glacial ice...

Enhanced-resolution L-band brightness temperature (TB) image time series generated from observations collected over the Greenland Ice Sheet by NASA's Soil Moisture Active Passive (SMAP) satellite are used to map Greenland's perennial firn aquifers from space. Exponentially decreasing L-band TB signatures are correlated with perennial firn aquifer areas identified via the Center for Remote Sensing of Ice Sheets (CReSIS) Multi-Channel Coherent Radar Depth Sounder (MCoRDS) that was flown by NASA's Operation IceBridge (OIB) campaign. An empirical algorithm to map extent is developed by fitting these signatures to a set of sigmoidal curves. During the spring of 2016, perennial firn aquifer areas are found to extend over ∼66 000 km2.

The Greenland ice sheet is one of the largest contributors to global meansea-level rise today and is expected to continue to lose mass as the Arcticcontinues to warm. The two predominant mass loss mechanisms are increasedsurface meltwater run-off and mass loss associated with the retreat ofmarine-terminating outlet glaciers. In this paper we use a large ensemble ofGreenland ice sheet models forced by output from a representative subset ofthe Coupled Model Intercomparison Project (CMIP5) global climate models to project ice sheet changes and sea-level risecontributions over the 21st century. The simulations are part of theIce Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We estimate thesea-level contribution together with uncertainties due to future climateforcing, ice sheet model formulations and ocean forcing for the twogreenhouse gas concentration scenarios RCP8.5 and RCP2.6. The resultsindicate that the Greenland ice sheet will continue to lose mass in bothscenarios until 2100, with contributions of 90±50 and 32±17 mm to sea-level rise for RCP8.5 and RCP2.6, respectively. The largestmass loss is expected from the south-west of Greenland, which is governed bysurface mass balance changes, continuing what is already observed today.Because the contributions are calculated against an unforced controlexperiment, these numbers do not include any committed mass loss, i.e. massloss that would occur over the coming century if the climate forcingremained constant. Under RCP8.5 forcing, ice sheet model uncertaintyexplains an ensemble spread of 40 mm, while climate model uncertainty andocean forcing uncertainty account for a spread of 36 and 19 mm,respectively. Apart from those formally derived uncertainty ranges, thelargest gap in our knowledge is about the physical understanding andimplementation of the calving process, i.e. the interaction of the ice sheetwith the ocean., info:eu-repo/semantics/published

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_discharge (Mankoff, 2020a) and at https://github.com/mankoff/ice_discharge (last access: 6 June 2020, Mankoff, 2020e).

One-dimensional simulations of firn evolution neglect horizontal advection from ice flow, which transports the firn column across climate gradients as it is buried by accumulation. Using a suite of model runs, we demonstrate the impacts of horizontal advection on the development of firn density, temperature, and the stratigraphy of melt features through the Greenland ice sheet percolation zone. The simulations isolate processes in synthetic runs and investigate four specific transects and an ice core site. Relative to one-dimensional simulations, the horizontal advection process tends to increase the pore close-off depth, reduce the heat content, and decrease the frequency of melt features with depth by emplacing firn sourced from higher locations under increasingly warm and melt-affected surface conditions. Preservation of the advected pore space and cold content is strongly dependent upon the depth of meltwater infiltration. Horizontal ice flow interacts with topography, climate gradients, and meltwater infiltration to influence the evolution of the firn column structure; the interaction between these variables modulates the impact of horizontal advection on firn at locations around Greenland. Pore close-off and firn temperature are mainly impacted in the lowermost 20–30 km of the percolation zone, which may be relevant to migration of the lower percolation zone. Relatively high in the percolation zone, however, the stratigraphy of melt features can have an advection-derived component that should not be conflated with changing climate.

Current sea-level rise partly stems from increased surface melting and meltwater runoff from the Greenland ice sheet. Multi-year snow, also known as firn, covers about 80% of the ice sheet and retains part of the surface meltwater. Since the firn cold content integrates its physical and thermal characteristics, it is a valuable tool for determining the meltwater-retention potential of firn. We use gap-filled climatological data from nine automatic weather stations in the ice-sheet accumulation area to drive a surface-energy-budget and firn model, validated against firn density and temperature observations, over the 1998-2017 period. Our results show a stable top 20 m firn cold content (CC20) at most sites. Only at the lower-elevation Dye-2 site did CC20 decrease, by 24% in 2012, before recovering to its original value by 2017. Heat conduction towards the surface is the main process feeding CC20 at all nine sites, while CC20 reduction occurs through low-cold-content fresh-snow addition at the surface during snowfall and latent-heat release when meltwater refreezes. Our simulations suggest that firn densification, while reducing pore space for meltwater retention, increases the firn cold content, enhances near-surface meltwater refreezing and potentially sets favourable conditions for ice-slab formation., SCOPUS: ar.j, info:eu-repo/semantics/published

We present changes of the ice margin in Northwest Greenland at the Eastern part of the Nunatarssuaq Ice Cap (NIC) over six decades. The ice margin in this area terminates as a near-vertical ice cliff of between 9 and 33 m thickness. During the years 1954–1957 and in 1965 multi-disciplinary studies were performed. We digitise and orthorectify material, that is often difficult to access, in order to use the historical data as an absolute starting point of our change assessment. We compare the cliff morphology of the mid-1950s and the mid-1960s with various time-steps between 1985 and 2017. The studied ice margin remained remarkably constant with very subtle changes of changing sign: rather slow advance rates are reported from the 1950s and 1960s that accelerated until 1985 and were followed by a general retreat until 2012 and a subsequent advance until 2017. Thickness changes are negative throughout the entire time-period, however, different rates of thinning are shown and there is a positive relationship with air temperature anomalies. Compared to similar elevations on the adjacent Greenland ice sheet, we find significantly weaker thinning rates at the NIC.

Seismic studies of glaciers yield insights into spatio-temporal processes within and beneath glaciers on scales relevant to flow and deformation of the ice. These methods enable direct monitoring of the bed in ways that complement other geophysical techniques, such as geodetic or ground penetrating radar observations. In this work, we report on the analysis of passive seismic data collected from the interior of the North East Greenland Ice Stream, the Greenland ice sheet's largest outlet glacier. We record thousands of basal earthquakes, many of which repeat with nearly identical waveforms. We also record many long-duration glacial tremor episodes that migrate across the seismic network with slow velocities (e.g. ~4–12 m s−1). Analysis of the basal earthquakes indicates a transition between times of individual event activity and times of tremor activity. We suggest that both processes are produced by shear slip at localized asperities along the bed. The transition between discrete and quasi-continuous slipping modes may be driven by pore-water pressure transients or heterogeneous strain accumulation in the ice due to strength contrasts of the underlying till.

Airborne radio-echo sounding (RES) surveys are widely used to measure ice-sheet bed topography. Measuring bed topography as accurately and widely as possible is of critical importance to modelling ice dynamics and hence to constraining better future ice response to climate change. Measurement accuracy of RES surveys is influenced both by the geometry of bed topography and the survey design. Here we develop a novel approach for simulating RES surveys over glaciated terrain, to quantify the sensitivity of derived bed elevation to topographic geometry. Furthermore, we investigate how measurement errors influence the quantification of glacial valley geometry. We find a negative bias across RES measurements, where off-nadir return measurement error is typically −1.8 ± 11.6 m. Topographic highlands are under-measured an order of magnitude more than lowlands. Consequently, valley depth and cross-sectional area are largely under-estimated. While overall estimates of ice thickness are likely too high, we find large glacier valley cross-sectional area to be under-estimated by −2.8 ± 18.1%. Therefore, estimates of ice flux through large outlet glaciers are likely too low when this effect is not taken into account. Additionally, bed mismeasurements potentially impact our appreciation of outlet-glacier stability.

We investigate changing precipitation patterns in the Kangerlussuaq region of western central Greenland during the Holocene thermal maximum (HTM), using a new chronology of ice sheet terminus position through the Holocene and a novel inverse modeling approach based on the unscented transform (UT). The UT is applied to estimate changes in annual precipitation in order to reduce the misfit between modeled and observed terminus positions. We demonstrate the effectiveness of the UT for time-dependent data assimilation, highlighting its low computational cost and trivial parallel implementation. Our results indicate that Holocene warming coincided with elevated precipitation, without which modeled retreat in the Kangerlussuaq region is more rapid than suggested by observations. Less conclusive is whether high temperatures during the HTM were specifically associated with a transient increase in precipitation, as the results depend on the assumed temperature history. Our results highlight the important role that changing precipitation patterns had in controlling ice sheet extent during the Holocene.

The melting of glaciers and ice sheets is nowadays considered a symbol of climate change. Many complex mechanisms are involved in the melting of ice, and, among these processes, surface darkening due to organic material on bare ice has recently received attention from the scientific community. The presence of microbes on glaciers has been shown to decrease the albedo of ice and promote melting. Despite several studies from the Himalaya, Greenland, Andes, and Alaska, no quantitative studies have yet been conducted in the European Alps. In this paper, we made use of DNA sequencing, microscopy and field spectroscopy to describe the nature of glacier algae found at a glacier (Vadret da Morteratsch) of the European Alps and to evaluate their effect on the ice-albedo feedback. Among different algal species identified in the samples, we found a remarkable abundance of Ancylonema nordenskioeldii, a species that has never previously been quantitatively documented in the Alps and that dominates algal blooms on the Greenland Ice Sheet. Our results show that, at the end of the ablation season, the concentration of Ancylonema nordenskioeldii on the glacier surface is higher than that of other algal species (i.e. Mesotaenium berggrenii). Using field spectroscopy data, we identified a significant correlation between a reflectance ratio (750 nm/650 nm) and the algae concentration. This reflectance ratio could be useful for future mapping of glacier algae from remote sensing data exploiting band 6 (740 nm) and band 4 (665 nm) of the MultiSpectral Instrument (MSI) on board Sentinel-2 satellite. Here we show that the biological darkening of glaciers (i.e. the bioalbedo feedback) is also occurring in the European Alps, and thus it is a global process that must be taken into account when considering the positive feedback mechanisms related to glacier melting.