The impact of climate oscillations on the surface energy budget over the Greenland Ice Sheet in a changing climate

Climate change is particularly strong in Greenland primarily as a result of changes in advection of heat and moisture fluxes from lower latitudes. The atmospheric structures involved influence the surface mass balance and their pattern are largely explained by climate oscillations which describe the internal climate variability. Based on a clustering method, we combine the Greenland Blocking Index and the North Atlantic Oscillation index with the vertically integrated water vapor to analyze inter-seasonal and regional impacts of the North Atlantic influence on the surface energy components over the Greenland Ice Sheet. In comparison to the reference period (1959–1990), the atmosphere has become warmer and moister during recent decades (1991–2020) for contrasting atmospheric circulation patterns. Particularly in the northern regions, increases in tropospheric water vapor enhance incoming longwave radiation and thus contribute to surface warming. Surface warming is most evident in winter, although its magnitude and spatial extent depend on the prevailing atmospheric configuration. Relative to the reference period, increases in sensible heat flux in the summer ablation zone are found irrespective of the atmospheric circulation pattern. Especially in the northern ablation zone, these are explained by the stronger katabatic winds which are partly driven by the larger surface pressure gradients between the ice/snow-covered surface and adjacent seas, and by the larger temperature gradient between near-surface air and the air above. Increases in net shortwave radiation are mainly connected to high-pressure systems. Whereas in the southern part of Greenland the atmosphere has gotten optical thinner, thus allowing more incoming shortwave radiation to reach the surface, in the northern part the incoming shortwave radiation flux has changed little with respect to the reference period, but the surface albedo decreased due to the expansion of the bare ice area.

No evidence for large subglacial source of mercury from the southwestern margin of the Greenland Ice Sheet.

In the current Matters Arising we present results from verifying control measurements of dissolved mercury (Hg) in glacial meltwater from the Greenland Ice Sheet (GrIS), which significantly challenges the conclusions of the recent publication by Hawkings et al. (2021). By direct measurements of meltwater in the same glacial catchment area, we demonstrate that the input Hg concentration for the regional upscaling in Hawkings et al (2021) is likely vastly over-estimated with major implications for the validity of the asserted extreme yield of Hg from the GrIS. In addition, we present a plausible explanation for the high Hg concentration values in the study, namely hitherto unidentified cross-contamination of water samples by mercury chloride (HgCl2), which was present and used for other purposes during field work. Together, the result of our control study potentially invalidates the suggested implications of geologically sourced Hg under the southwestern margin of the GrIS on the Arctic ecosystem in both current and future climate conditions.

An empirical algorithm to map perennial firn aquifers and ice slabs within the Greenland Ice Sheet using satellite L-band microwave radiometry

Perennial firn aquifers are subsurface meltwater reservoirs consisting of a meters-thick water-saturated firn layer that can form on spatial scales as large as tens of kilometers. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting and high snow accumulation. Widespread perennial firn aquifers have been identified within the Greenland Ice Sheet (GrIS) via field expeditions, airborne ice-penetrating radar surveys, and satellite microwave sensors. In contrast, ice slabs are nearly continuous ice layers that can also form on spatial scales as large as tens of kilometers as a result of surface and subsurface water-saturated snow and firn layers sequentially refreezing following multiple melting seasons. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting but in areas where snow accumulation is at least 25 % lower as compared to perennial firn aquifer areas. Widespread ice slabs have recently been identified within the GrIS via field expeditions and airborne ice-penetrating radar surveys, specifically in areas where perennial firn aquifers typically do not form. However, ice slabs have yet to be identified from space. Together, these two ice sheet features represent distinct, but related, sub-facies within the broader percolation facies of the GrIS that can be defined primarily by differences in snow accumulation, which influences the englacial hydrology and thermal characteristics of firn layers at depth. Here, for the first time, we use enhanced-resolution vertically polarized L-band brightness temperature (TVB) imagery (2015–2019) generated using observations collected over the GrIS by NASA's Soil Moisture Active Passive (SMAP) satellite to map perennial firn aquifer and ice slab areas together as a continuous englacial hydrological system. We use an empirical algorithm previously developed to map the extent of Greenland's perennial firn aquifers via fitting exponentially decreasing temporal L-band signatures to a set of sigmoidal curves. This algorithm is recalibrated to also map the extent of ice slab areas using airborne ice-penetrating radar surveys collected by NASA's Operation IceBridge (OIB) campaigns (2010–2017). Our SMAP-derived maps show that between 2015 and 2019, perennial firn aquifer areas extended over 64 000 km2, and ice slab areas extended over 76 000 km2. Combined together, these sub-facies are the equivalent of 24 % of the percolation facies of the GrIS. As Greenland's climate continues to warm, seasonal surface melting will increase in extent, intensity, and duration. Quantifying the possible rapid expansion of these sub-facies using satellite L-band microwave radiometry has significant implications for understanding ice-sheet-wide variability in englacial hydrology that may drive meltwater-induced hydrofracturing and accelerated ice flow as well as high-elevation meltwater runoff that can impact the mass balance and stability of the GrIS.

Impact of freshwater runoff from the southwest Greenland Ice Sheet on fjord productivity since the late 19th century

Climate warming and the resulting acceleration of freshwater discharge from the Greenland Ice Sheet are impacting Arctic marine coastal ecosystems, with implications for their biological productivity. To accurately project the future of coastal ecosystems, and place recent trends into perspective, paleo-records are essential. Here, we present late 19th century to present runoff estimates for a large sub-Arctic fjord system (Nuup Kangerlua, southwest Greenland) influenced by both marine- and land-terminating glaciers. We followed a multiproxy approach to reconstruct spatial and temporal trends in primary production from four sediment cores, including diatom fluxes and assemblage composition changes, biogeochemical and sedimentological proxies (total organic carbon, nitrogen, C / N-ratio, biogenic silica, δ13C, δ15N, grain size distribution). We show that an abrupt increase in freshwater runoff in the mid-1990’s is reflected by a 3-fold increase in biogenic silica fluxes in the glacier-proximal area of the fjord. In addition to increased productivity, freshwater runoff modulates the diatom assemblages and drives the dynamics and magnitude of the diatom spring bloom. Our records indicate that marine productivity is higher today than it has been at any point since the late 19th century and suggest that increased mass loss of the Greenland Ice Sheet is likely to continue promoting high productivity levels at sites proximal to marine-terminating glaciers. We highlight the importance of paleo-records in offering a unique temporal perspective on ice-ocean-ecosystem responses to climate forcing beyond existing remote sensing or monitoring time-series.

Filling and drainage of a subglacial lake beneath the Flade Isblink ice cap, northeast Greenland

The generation, transport, storage and drainage of meltwater beneath the ice sheet play important roles in the Greenland ice sheet (GrIS) system. Active subglacial lakes, common features in Antarctica, have recently been detected beneath GrIS and may impact ice sheet hydrology. Despite their potential importance, few repeat subglacial lake filling and drainage events have been identified under Greenland Ice Sheet. Here we examine the surface elevation change of a collapse basin at the Flade Isblink ice cap, northeast Greenland, which formed due to sudden subglacial lake drainage in 2011. We estimate the subglacial lake volume evolution using multi-temporal ArcticDEM data and ICESat-2 altimetry data acquired between 2012 and 2021. Our long-term observations show that the subglacial lake was continuously filled by surface meltwater, with basin surface rising by up to 55 m during 2012–2021 and we estimate 138.2 × 106 m3 of meltwater was transported into the subglacial lake between 2012 and 2017. A second rapid drainage event occurred in late August 2019, which induced an abrupt ice dynamic response. Comparison between the two drainage events shows that the 2019 drainage released much less water than the 2011 event. We conclude that multiple factors, e.g., the volume of water stored in the subglacial lake and bedrock relief, regulate the episodic filling and drainage of the lake. By comparing the surface meltwater production and the subglacial lake volume change, we find only ~64 % of the surface meltwater successfully descended to the bed, suggesting potential processes such as meltwater refreezing and firn aquifer storage, need to be further quantified.

Millennial-scale migration of the frozen/melted basal boundary, western Greenland ice sheet

The geometry and thermal structure of western Greenland ice sheet are known to have undergone relatively substantial change over the Holocene. Evolution of the frozen and melted fractions of the bed associated with the ice-sheet retreat over this time frame remains unclear. We address this question using a thermo-mechanically coupled flowline model to simulate a 11 ka period of ice-sheet retreat in west central Greenland. Results indicate an episode of ~100 km of terminus retreat corresponded to ~16 km of upstream frozen/melted basal boundary migration. The majority of migration of the frozen area is associated with the enhancement of the frictional and strain heating fields, which are accentuated toward the retreating ice margin. The thermally active bedrock layer acts as a heat sink, tending to slow contraction of frozen-bed conditions. Since the bedrock heat flux in our region is relatively low compared to other regions of the ice sheet, the frozen region is relatively greater and therefore more susceptible to marginward changes in the frictional and strain heating fields. Migration of melted regions thus depends on both geometric changes and the antecedent thermal state of the bedrock and ice, both of which vary considerably around the ice sheet.

Concurrent Changepoints in Greenland Ice Core δ18O records and the North Atlantic Oscillation over the Past Millennium

Structural changes, or changepoints, coinciding in multiple ice core records over the Greenland Ice Sheet (GrIS) may reflect a widespread response of the GrIS to atmospheric forcing. Thus, to better understand how atmospheric circulation may regulate sudden changes in δ18O of Greenland precipitation, we seek synchronous changepoints occurring in ice core-derived δ18O time series across the GrIS and in the North Atlantic Oscillation (NAO) over the past millennium. By utilizing a Bayesian changepoint detection method, four changepoint horizons were revealed: at the beginning of the 20th century, in the late-15th century, and around the turn of the 11th and 10th centuries. Although the changepoints in ice core δ18O records exhibited distinctive spatial arrangements in each horizon, all corresponded to changepoints in the NAO, indicative of a consistent atmospheric influence on GrIS surface changes over the past millennium.

Microporphyritic and microspherulitic melt grains, Hiawatha crater, Northwest Greenland: Implications for post-impact cooling rates, hydration, and the cratering environment

Sand-sized impactite melt grains hand-picked from a glaciofluvial sample proximal to the Hiawatha impact crater in Northwest Greenland contain new information about the crystallization and cooling history of this impact structure, which is concealed by the Greenland Ice Sheet. Of course, the original locations of the individual sand grains are unknown, but this is offset by the substantial number and wide variety of impactite grains available for study. A detailed investigation of 16 melt grains shows that post-cratering crystallization took place under very variable conditions of strong undercooling with temperatures that dropped rapidly from high above their solidus to far below. A distinct event of near-isochemical hydration at above or ∼250 °C is recorded by intense perlitic fracturing and the growth of closely packed mordenite spherulites only 1−3 μm across in felsic melt grains, which was followed by lower temperature hydrothermal alteration along the pre-existing perlitic fractures. The formation of abundant mordenite microspherulites appears to be very rare or not previously recorded in impactite melts and suggests the rapid infilling of the Hiawatha crater by a hydrous source. The infilling did not occur immediately after the impact as in submarine impacts, but soon thereafter, and before the establishment of a low-temperature hydrothermal alteration system common to the waning stage of cooling in many impact structures. These observations and previous documentation of terrestrial organic matter in the impactites are consistent with an impact into a water-rich terrestrial environment, such as through the Greenland Ice Sheet or into a forested, lacustrine−fluvial region.