Search

The Arctic is warming rapidly, significantly reducing the Greenland ice sheet (GrIS) surface mass balance (SMB) and raising its contribution to global sea‐level rise. Since these trends are expected to continue, it is essential to explore the GrIS SMB response to projected climate warming. We compare projections from three polar regional climate models, RACMO, MAR, and HIRHAM, forced by the Community Earth System Model CESM2 under a high‐end warming scenario (SSP5‐8.5, 1970–2099). We reveal different modeled SMB by 2100, including a twofold larger annual surface mass loss in MAR (−1735 Gt/yr) and HIRHAM (−1698 Gt/yr) relative to RACMO (−964 Gt/yr). Discrepancies primarily stem from differences in projected runoff, triggering melt‐albedo positive feedback and subsequent modeled ablation zone expansion. In addition, we find different responses of modeled meltwater production to similar atmospheric warming. Our analysis suggests clear avenues for model developments to further improve SMB projections and contribution to sea‐level rise. Plain Language Summary: We explore how three different polar climate models predict the future surface mass balance (SMB) of the Greenland ice sheet (GrIS), a major contributor to global sea‐level rise. Our results show that SMB projections among these models differ significantly by the end of the century. Differences primarily stem from how models convert meltwater into surface runoff toward the ocean, a crucial SMB component that directly contributes to sea‐level rise. Another key factor is the response of these models to climate warming, substantially affecting future melting rates. Our research highlights the need for further improvements of these climate models. By identifying and understanding the drivers behind model differences, we can improve SMB predictions and better estimate future global sea‐level rise. Key Points: With identical forcing, Greenland Ice Sheet surface mass balance from 3 regional climate models shows a two‐fold difference by 2100Different runoff projections stem from substantial discrepancies in projected ablation zone expansion, and reciprocallyThe response of meltwater production to similar atmospheric warming differs significantly among regional climate models [ABSTRACT FROM AUTHOR]

At least three primary iceberg calving styles have been identified in Greenland: serac collapse, which produces falling icebergs tens of meters in length; slab capsize, which produces rotating icebergs hundreds of meters in length; and tabular rifting, which produces kilometer‐scale icebergs. However, calving styles are mostly undocumented across Greenland. Here, we develop a method to disentangle the sizes of individual calving events and map the dominant calving style at glaciers, using the characteristic properties of step retreats in satellite‐derived terminus positions. At glaciers known to frequently produce calving teleseisms, step retreats greater than 200 m account for >80% ${ >} 80\%$ of net calved length since 2018. In contrast, at glaciers known to calve by serac failure, 200 m step retreats account <20% ${< } 20\%$ of net calving. Thus, terminus change timeseries can offer promising insight into the dominant calving styles at marine‐terminating glaciers. Plain Language Summary: Iceberg calving is ubiquitous at the ocean‐bounded edges of ice sheets. However, not all iceberg calving is the same: at least three distinct mechanisms govern the discharge of solid ice from the Greenland Ice Sheet. We show these are associated with a characteristic range of iceberg sizes. Because calving constitutes a major component of Greenland's continued mass loss, it is important that these distinct mechanisms be understood. We develop a method to identify which calving mechanisms are important at various glaciers. We do so by using satellite observations to estimate the lengths of terminus retreats producing icebergs at individual glaciers. Key Points: Greenland outlet glaciers with different calving styles have different characteristic retreat magnitudes in consecutive satellite imagesAt glaciers which frequently produce calving teleseisms, step retreats in terminus position are dominated by retreats >200 mThe distinct processes underlying different calving styles will likely necessitate several distinct calving laws [ABSTRACT FROM AUTHOR]

Pleistocene Ice Ages display abrupt Dansgaard-Oeschger (DO) climate oscillations that provide prime examples of Earth System tipping points--abrupt transition that may result in irreversible change. Greenland ice cores provide key records of DO climate variability, but gas-calibrated estimates of the temperature change magnitudes have been limited to central and northwest Greenland. Here, we present ice-core δ15N-N2 records from south (Dye 3) and coastal east Greenland (Renland) to calibrate the local water isotope thermometer and provide a Greenland-wide spatial characterization of DO event magnitude. We combine these data with existing records of δ18O, deuterium excess, and accumulation rates to create a multiproxy "fingerprint" of the DO impact on Greenland. Isotope-enabled climate models have skill in simulating the observational multiproxy DO event impact, and we use a series of idealized simulations with such models to identify regions of the North Atlantic that are critical in explaining DO variability. Our experiments imply that wintertime sea ice variation in the subpolar gyre, rather than the commonly invoked Nordic Seas, is both a sufficient and a necessary condition to explain the observed DO impacts in Greenland, whatever the distal cause. Moisture-tagging experiments support the idea that Greenland DO isotope signals may be explained almost entirely via changes in the vapor source distribution and that site temperature is not a main control on δ18O during DO transitions, contrary to the traditional interpretation. Our results provide a comprehensive, multiproxy, data-model synthesis of abrupt DO climate variability in Greenland. [ABSTRACT FROM AUTHOR]

The coastal circulation around Southern Greenland transports fresh, buoyant water masses from the Arctic and Greenland Ice Sheet near regions of convection, sinking, and deep-water formation in the Irminger and Labrador Seas. Here, we track the pathways and fate of these fresh water masses by initializing synthetic particles in the East Greenland Coastal Current on the Southeast Greenland shelf and running them through altimetry-derived surface currents from 1993 to 2021. We report that the majority of waters (83%) remain on the shelf around the southern tip of Greenland. Variability in the shelf-basin exchange of the remaining particles closely follows the number of tip jet wind events on seasonal and interannual timescales. The probability of a particle exiting the shelf increases almost fivefold during a tip jet event. These results indicate that the number of tip jets is a close proxy of the shelf-basin exchange around Southern Greenland. [ABSTRACT FROM AUTHOR]

The surface elevation of the Greenland Ice Sheet is constantly changing due to the interplay between surface mass balance processes and ice dynamics, each exhibiting distinct spatiotemporal patterns. Here, we employ satellite and airborne altimetry data with fine spatial (1 km) and temporal (monthly) resolutions to document this spatiotemporal evolution from January 2003 to August 2023. To estimate elevation changes of the Greenland Ice Sheet (GIS), we utilize radar altimetry data from CryoSat-2 and EnviSat, laser altimetry data from the ICESat and ICESat-2, and laser altimetry data from NASA's Operation IceBridge Airborne Topographic Mapper. We produce continuous monthly ice surface elevation changes from January 2003 to August 2023 on a 1 km grid covering the entire GIS. We estimate cumulative ice loss of 4,352 Gt ± 315 Gt (12.1 ± 0.9 mm sea level equivalent) during this period, excluding peripheral glaciers. Between 2003 and 2023, the ice sheet land-terminating margin underwent a significant cumulative thinning of several meters. Ocean-terminating glaciers exhibited thinning between 20–40 m, with Jakobshavn Isbræ experiencing an exceptional thinning of nearly 70 m. This dataset of fine-resolution altimetry data in both space and time will support studies of ice mass loss and useful for GIS ice sheet modelling. To validate our monthly mass changes of the Greenland ice sheet, we use mass change from satellite gravimetry and mass change from the Input-Output method. On multiannual timescales, there is a strong correlation between the time series, with R values ranging from 0.88 to 0.92. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet is highly sensitive to climate change, leading to significant retreat along its edges. This rapid ice loss contributes to rising sea levels and impacts the Earth's climate stability. Understanding the extent of recent glacier retreat is crucial in order to determine if it is unprecedented or within ranges of natural variability. Palaeoenvironmental studies aim to identify past glacial phases and landscape changes using advanced dating methods such as cosmic ray exposure (CRE) dating. In NE Greenland, CRE dating has helped establish the timing of glacial oscillations, yet a comprehensive understanding of glacial fluctuations during specific periods still needs to be developed. This study aims to chronologically constrain the postglacial landscape evolution of two NE Greenland valleys from the Young Sund–Tyrolerfjord area (74°N, 20–25°E) from the onset of deglaciation and throughout the Holocene to better understand glacial and postglacial changes. The chronological framework relies on 27 10Be cosmic‐ray exposure ages that constrain our interpretation of the geomorphological features in both valleys. Inconsistencies were observed in the ages dataset, highlighting potential bias associated with nuclide inheritance and post‐glacial dynamics. Despite limitations, the CRE results confirm the general pattern observed in NE Greenland: (i) major deglaciation and disconnection of glaciers from the main glacial systems during the Lateglacial and Early Holocene with a rapid but not homogeneous deglaciation within the range from ~14.3 to 11.9 ka; (ii) no evidence of glacial activity during the Middle Holocene, probably associated with the withdrawn position of the ice masses’ fronts; and (iii) glacier expansion during the Late Holocene, with a Little Ice Age advance as the last significant period of glacial regrowth. [ABSTRACT FROM AUTHOR]

Given the long response time of ice sheets, simulations of the Greenland ice sheet typically exceed the availability of input climate data to reliably simulate the fast processes underlying surface mass balance. Strong feedback processes are known to make the mass balance sensitive to inter- and intra-annual variability. Even simulations with climate models do not always cover the full period of interest, motivating bridging these gaps using relatively coarsely resolved climate reconstructions or temporal interpolation methods. However, both of these approaches usually only provide information about the climatological average but not variability. We investigate how this simplification impacts the surface mass balance using the BErgen Snow SImulator. The model was run for up to 500 years using the same atmospheric climatology but different synthetic variabilities. While changing inter-annual variations has an impact of less than 5 % on the surface mass balance of the Greenland ice sheet, neglecting intra-annual variability by using a daily climatology causes a 40 % change in mass balance. Decomposing the total effect into contributions from different input variables, the biggest contributor is precipitation followed by temperature. Using a daily climatology, a small amount of snowfall every day overestimates the albedo and thus surface mass balance (SMB). We propose a correction that re-captures the effect of intermittent precipitation, reducing the SMB overestimation to 15 %–25 %. We conclude that simulations of the Greenland surface mass and energy balance should be forced with a transient climate, in particular for models that are calibrated with transient data. [ABSTRACT FROM AUTHOR]

As the only remaining ice sheet in the Northern Hemisphere, the Greenland ice sheet (GrIS) plays a crucial role in influencing atmospheric circulations, particularly with its rapid melting under global warming. In this paper, the influences of GrIS topography and surface thermal conditions are investigated by a series of aquaplanet experiments. The results show that the GrIS topography induces stationary waves and favors more blocking events through the generation of negative potential vorticity (PV) anomalies, while it tends to suppress local storm activities through the induced stationary waves. The surface cooling center of the GrIS is found to strengthen the jet streams by enhancing the meridional temperature gradient and thermal wind, while it causes the PV and static stability to increase during near-Greenland blocking days, thereby disfavoring blocking onset. Altogether, the topography and surface thermal effects of GrIS appear to compete with each other so that the net effect would determine the final response. Nevertheless, nonlinearity is found in both GrIS-topography alone and GrIS-surface temperature alone experiments, where nonlinear responses of atmospheric circulation are detected when the GrIS topography height or surface temperature exceeds their critical values, respectively. Hence, through this study, the response of the blocking in the vicinity of Greenland to the combined effects of topography and surface thermal conditions may shed light on comprehending the underlying mechanism of blocking alteration in a changing climate. Significance Statement: Although there have been numerous observation-based studies showing that the blocking around Greenland has increased over the past few decades, which is a predominant driver in accelerating the Greenland ice sheet (GrIS) melting, the reasons for this change are still unclear. In addition, the impact of GrIS on blocking still needs investigation. Through idealized aquaplanet simulations, this study examines the effect of GrIS's topography and surface thermal conditions upon blocking onset. It suggests the nonlinearity of the blocking response in the vicinity of Greenland to the combined effects in the variation of the height and surface temperature of GrIS. This study will enhance our understanding of possible changes in blocking due to global warming. [ABSTRACT FROM AUTHOR]

Objectives: As a sensitive indicator of global changes such as ice sheet mass loss, sea level rise and climate change, supraglacial lakes play an important role in the study of ice sheet stability. Methods: U-Net is a semantic segmentation model, which extracts supraglacial lakes from Landsat satellite images in different situations. Results: The extraction of supraglacial lakes by U-Net possesses the optimal performance and efficiency. And the supraglacial lakes near the Sermeq Avannarleq Glacier in Greenland ice sheet are extracted from 2000 to 2020, and the seasonal and interannual variations of the area are analyzed. Conclusions: The area of supraglacial lakes starts to increase in mid-May, reaches the peak in the midlate-July, and almost disappears in September in a year. For the long time series change, the area of supra‐ glacial lakes in this region shows an increasing trend from 2000 to 2020, indicating that the amount of ice sheet melting is gradually increasing. [ABSTRACT FROM AUTHOR]

The integration of data from multiple satellite altimetry missions, each offering unique observational characteristics, has enabled us to discern both short-term variability and long-term climate trends affecting Greenland's peripheral glaciers and the Greenland Ice Sheet (GrIS). Our methodology, adapted and informed by lessons learned from our similar efforts in Antarctica, ensures the consistency and reliability of the derived elevation change dataset. The data covers the years 1992–2023 and is made publicly available as part of the NASA's Making Earth System Data Records for Use in Research Environments (MEaSUREs) Inter-Mission Time Series of Land Ice Velocity and Elevation (ITS_LIVE) project. Our analysis of the dataset reveals significant patterns of mass loss across the GrIS. We find that the ice sheet and peripheral glaciers have experienced an average mass loss of -173 ± 19 Gt a-1 and -23 ± 5 Gt a-1 respectively over the 1992–2022 time period (given temporal availability of selected firn models), with notable variations over time. Specifically, the early years of the record exhibit a positive mass balance, likely due to anomalously positive surface mass balance. However, this trend shifts in subsequent years, with a pronounced increase in mass loss rates, highlighting the accelerating impact of a changing climate on ice sheet dynamics and surface mass balance. Moreover, our analysis underscores the importance of considering peripheral glaciers in addition to the continental ice sheet when assessing overall mass trends. By incorporating data from peripheral glaciers, we provide a more comprehensive understanding of Greenland's total contributions to global sea level rise. Our findings reveal not only the magnitude of mass loss but also its evolution over time, emphasizing the need for continued monitoring and research to better understand the impacts of climate change on Earth's cryosphere. [ABSTRACT FROM AUTHOR]

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]

Meltwater from the Greenland Ice Sheet can alter the continental shelf/slope circulation and cross-shelf freshwater fluxes and limit deep convection in adjacent basins through surface freshening. We explore the impacts on the West Greenland Current and eastern Labrador Sea with different vertical distributions of the meltwater forcing. In this study, we present the results from global coupled ocean/sea ice simulations, forced with atmospheric reanalysis, that are mesoscale eddy-active (∼2–3-km horizontal spacing) and eddy-permitting (∼6–7-km horizontal spacing) in the study region. We compare the West Greenland Current in mesoscale eddy-active and eddy-permitting without meltwater to highlight the role of small-scale features. The mesoscale eddy-active configuration is then used to assess the change in the eastern Labrador Sea when meltwater is added to the surface or vertically distributed to account for mixing within fjords. In both simulations with meltwater, the West Greenland and West Greenland Coastal Currents are faster than in the simulation with no meltwater; their mean surface speeds are the highest in the vertical distribution case. In the latter case, there is enhanced baroclinic conversion at the shelf break compared to the simulation with no meltwater. When meltwater is vertically distributed, there is an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy. In addition, in the eastern Labrador Sea, the salinity is lower and the meltwater volume is greater when meltwater is vertically distributed. Therefore, the West Greenland Current is sensitive to how meltwater is added to the ocean with implications for the freshening of the Labrador Sea. Significance Statement: Our goal is to understand how the flux of freshwater across the West Greenland continental slope into the Labrador Sea is modified by meltwater from the Greenland Ice Sheet. We compare the simulations of the ocean that capture key dynamics along the West Greenland continental slope that have no meltwater, meltwater added to the ocean surface, and meltwater distributed vertically to represent the mixing within fjords. When meltwater is added, the currents along the continental slope are faster, with the greatest increase when meltwater is vertically distributed. In that case, there is enhanced freshening of the Labrador Sea because modified density gradients generate more eddies. Proper representation of the vertical structure of meltwater is important for projecting the impact of freshwater on the subpolar North Atlantic. [ABSTRACT FROM AUTHOR]

Present understanding of Greenland's subglacial geology is derived mostly from interpolation of geologic mapping of its ice‐free margins and unconstrained by geophysical data. Here we refine the extent of its geologic provinces by synthesizing geophysical constraints on subglacial geology from seismic, gravity, magnetic and topographic data. North of 72°N, no province clearly extends across the whole island, leaving three distinct subglacial regions yet to be reconciled with margin geology. Geophysically coherent anomalies and apparent province boundaries are adjacent to the onset of faster ice flow at both Petermann Glacier and the Northeast Greenland Ice Stream. Separately, based on their subaerial expression, dozens of unusually long, straight and sub‐parallel subglacial valleys cross Greenland's interior and are not yet resolved by current syntheses of its subglacial topography. Plain Language Summary: The Greenland Ice Sheet obscures the rocks beneath 79% of Greenland. By necessity, scientists have relied mostly on studying the rocks exposed along Greenland's edge to understand the island's interior geology. We examine geophysical data from seismometers on the ground, satellites that measure Earth's gravity and magnetic fields and surface topography, and aircraft that measure those same properties and ice thickness. We draw a new map of Greenland's geology beneath the ice sheet by examining where those data show similar signals regarding the nature of the underlying rock, and where they could be related to mapped rock exposures. We also find evidence of some areas with geophysical expressions that are distinct from the rocks found at the island's edges. Some geologic structures, which are entirely covered by ice, may affect how ice flows from Greenland's vast interior toward its coast. Finally, we identify many valleys beneath the ice that are very long and often aligned with each other, but which are not yet fully captured in present maps of the topography beneath the ice. Key Points: We produced a new synthesis of subglacial boundaries for Greenland's geologic provinces from seismic, gravity, magnetic and topography dataThree subglacial regions in central and northern Greenland cannot yet be reconciled with surface‐exposed geologic provincesWe find evidence for a large subglacial valley network that is not fully resolved by subglacial topography syntheses for Greenland [ABSTRACT FROM AUTHOR]

We analyze the impact of a disintegrated Greenland Ice Sheet (GrIS) on the global climate through steady-state simulations with the MPI-ESM (Max Planck Institute for Meteorology Earth System Model). This advances our understanding of the intricate feedbacks between the GrIS and the full climate system. Sensitivity experiments enable the quantification of the individual contributions of altered Greenland surface elevation and properties (e.g., land cover) to the atmospheric and oceanic climate response. Removing the GrIS results in reduced mechanical atmospheric blocking, warmer air temperatures over Greenland and thereby changes in the atmospheric circulation. The latter alters the wind stress on the ocean, which controls the ocean-mass transport through the Arctic Gateways. Without the GrIS, the upper Nordic Seas are fresher, attenuating deep-water formation. In the Labrador Sea, deep-water formation is weaker despite a higher upper-ocean salinity, as the inflow of dense overflow from the Denmark Strait is reduced. Our sensitivity experiments show that the atmospheric response is primarily driven by the lower surface elevation, whereas altered Greenland surface properties mostly amplify but also counteract few of the changes. The lower Greenland elevation dominates the ocean response through wind-stress changes. Only in the Labrador Sea, altered Greenland surface properties dominate the ocean response, as this region stores excessive heat from the Greenland warming. The main drivers vary vertically: The elevation effect controls upper-ocean densities, while surface properties are important for the intermediate and deep ocean. Despite the confinement of most responses to the Arctic, a disintegrated GrIS also influences remote climates. The altered climate in response to a GrIS disintegration also constrains a potential ice-sheet regrowth to high-bedrock eastern Greenland. [ABSTRACT FROM AUTHOR]

We present an ice-penetrating radar data set acquired over the Greenland ice sheet by aircraft during the years 1971, 1972, 1974, 1978, and 1979. The data set comprises over 177 000 km of flight lines and contains a wealth of information on the state of the Greenland ice sheet, including information on ice thickness and englacial properties. During data collection in the 1970s, the data were recorded on optical film rolls, and in this paper, we document the digitization of these film rolls and their associated geographical information. Our data digitization enables interaction with and analysis of the data and facilitates comparison with modern-day radar observations. The complete data set in full resolution is available in the Stanford Digital Repository (10.25740/wm135gp2721;), with the associated technical reports. Part of the data set is available as low-resolution JPG files at the Technical University of Denmark's data repository (10.6084/m9.figshare.c.7235299.v1;), with associated technical reports and digitized geographical information. The Stanford Digital Repository serves as long-term storage, providing archival historic preservation in perpetuity, and is not intended as a primary data access point. The Technical University of Denmark (DTU) data repository serves as a primary entry point for data access, with files organized according to acquisition year and flight line in a simple folder structure. Here, we release the full data sets to enable the larger community to access and interact with the data. [ABSTRACT FROM AUTHOR]

Basal materials in ice cores hold information about paleoclimate conditions, glacial processes, and the timing of past ice-free intervals, all of which aid understanding of ice sheet stability and its contribution to sea level rise in a warming climate. Only a few cores have been drilled through ice sheets into the underlying sediment and bedrock, producing limited material for analysis. The last of three Camp Century ice cores, which the U.S. Army collected in northwestern Greenland from 1963–1966 CE, recovered about 3.5 m of subglacial material, including ice and sediment. Here, we document the scientific history of the Camp Century subglacial material. We present our recent core-cutting, sub-sampling, and processing methodology and results for this unique archive. In 1972 CE, curators at the Buffalo, New York, Ice Core Laboratory cut the original core sections into 32 segments that were each about 10 cm long. Since then, two segments were lost and are unaccounted for, two were thawed, and two were cut as pilot samples in 2019 CE. Except for the two thawed segments, the rest of the extant core has remained frozen since collection. In 2021 CE, we documented, described, and then cut each of the remaining frozen archived segments (n=26). We saved an archival half and cut the working half into eight oriented sub-samples under controlled temperature and light conditions for physical, geochemical, isotopic, sedimentological, magnetic, and biological analyses. Our approach was designed to maximize sample usage for multiproxy analysis, minimize contamination, and preserve archive material for future analyses of this legacy subglacial material. Grain size, bulk density, sedimentary features, magnetic susceptibility, and ice content, as well as pore ice pH and conductivity, suggest that the basal sediment contains five stratigraphic units. We interpret these stratigraphic units as representing different depositional environments in subglacial or ice-free conditions: from bottom to top, a diamicton with subhorizontal ice lenses (Unit 1), vertically fractured ice with dispersed fine-grained sediments (<20 % in mass) (Unit 2), a normally graded bed of pebbles to very fine sand in an icy matrix (Unit 3), bedded very fine to fine sand (Unit 4), and stratified medium to coarse sand (Unit 5). Plant macrofossils are present in all samples and are most abundant in Units 3 and 4; insect remains are present in some samples (Units 1, 3, and 5). Our approach provides a working template for future studies of ice core basal materials because it includes intentional planning of core sub-sampling, processing methodologies, and archiving strategies to optimize the collection of paleoclimate, glacial process, geochemical, geochronological, and sediment properties from archives of limited size. Our work benefited from a carefully curated and preserved archive, allowing the application of analytical techniques not available in 1966 CE. Preserving uncontaminated core material for future analyses that use currently unavailable tools and techniques is an important consideration for rare archive materials such as these from Camp Century. [ABSTRACT FROM AUTHOR]

Extreme melt and rainfall events can induce temporary acceleration of Greenland Ice Sheet motion, leading to increased advection of ice to lower elevations where melt rates are higher. In a warmer climate, these events are likely to become more frequent. In September 2022, seasonally unprecedented air temperatures caused multiple melt events over the Greenland Ice Sheet, generating the highest melt rates of the year. The scale and timing of the largest event overwhelmed the subglacial drainage system, enhancing basal sliding and increasing ice velocities by up to ∼240% relative to pre‐event velocities. However, ice motion returned rapidly to pre‐event levels, and the speed‐ups caused a regional increase in annual ice discharge of only ∼2% compared to when the effects of the speed‐ups were excluded. Therefore, although late melt‐season events are forecast to become more frequent and drive significant runoff, their impact on net mass loss via ice discharge is minimal. Plain Language Summary: Extreme melt and rainfall events can cause the flow of ice on the Greenland Ice Sheet to accelerate, potentially causing more ice to move to lower elevations, where temperatures are warmer and melt rates are higher. In September 2022, there were multiple unprecedented melt events. Their intensity caused some glaciers on the ice sheet to speed up by 240% relative to pre‐event speeds. Despite these accelerations, our analyses show that these events had only a minimal long‐term impact on how much ice was moved to lower elevations due to the short duration of the speed‐ups. As a result, while these melt‐induced speed‐ups are expected to become more common in a warmer climate, their effect on the amount of ice transported toward the ice margins is minimal. Key Points: September 2022 saw multiple melt events over the west Greenland Ice Sheet, with the largest daily runoff of any late melt‐season since 1950Large quantities of surface‐generated meltwater overwhelmed the subglacial drainage system causing brief increases in ice velocityLate‐season runoff‐induced speed‐ups have only minimal impact on the mass balance of the Greenland Ice Sheet via dynamical processes [ABSTRACT FROM AUTHOR]

The melting of the Greenland ice sheet (GrIS) in summer has rapidly and significantly increased in recent decades, especially for the northern GrIS. Circulation related to GrIS melting is important for understanding the contribution of the GrIS to the global sea level. In this paper, we used the SOM method to obtain three spatial patterns of GrIS melting based on model output data: overall melting, northern melting, and southern melting patterns. We also examined their linkages to the observed atmospheric circulation. GrIS melting is primarily related to Greenland blocking (GB), while different types of GB are responsible for different melting patterns. The overall and northern melting patterns are both associated with high-latitude GB, which is associated with the decay and westward movement of mid-latitude and high-latitude European blocking (EB), respectively. It is found that the generation of mid-latitude and high-latitude EBs are related to mid-latitude eastern North Atlantic warming and Greenland–Norskehavet–North Sea warming, respectively, while the movement of EB is related t upstream weakened zonal wind anomalies. Furthermore, the southern melting pattern is linked to mid-latitude GB, which is associated with the wave train from the East Pacific to Southern Greenland through North America and further related to the mid-latitude of East Pacific warming. [ABSTRACT FROM AUTHOR]

Ice slabs are multi-meter-thick layers of refrozen ice that limit meltwater storage in firn, leading to enhanced surface runoff and ice sheet mass loss. To date, ice slabs have primarily been mapped using airborne ice-penetrating radar, which has limited spatial and temporal coverage. This makes it difficult to fully assess the current extent and continuity of ice slabs or to validate predictive models of ice slab evolution that are key to understanding their impact on Greenland's surface mass balance. Here, for the first time, we map the extent of ice slabs and superimposed ice facies across the entire Greenland Ice Sheet at 500 m resolution using dual-polarization Sentinel-1 (S-1) synthetic-aperture radar (SAR) data collected in winter 2016–2017. We do this by selecting empirical thresholds for the cross-polarized backscatter ratio and HV backscattered power that jointly optimize the agreement between airborne ice-penetrating radar data detections of ice slabs and the S-1 estimates of ice slab extent. Our results show that there is a strong correlation between C-band backscatter and the ice content of the upper ∼ 7 m of the firn column that enables ice slab mapping with S-1. Our new mapping shows that ice slabs are nearly continuous around the entire margin of the ice sheet. This includes regions in southwest Greenland where ice slabs have not been previously identified by ice-penetrating radar but where the S-1 -inferred ice slab extent shows strong agreement with the extent of visible runoff mapped from optical imagery. The algorithm developed here lays the groundwork for the long-term monitoring of ice slab expansion with current and future C-band satellite systems and highlights the potential added value of future L-band missions for near-surface studies in Greenland. [ABSTRACT FROM AUTHOR]

The increasing anthropogenic CO 2 forcing of the climate system calls for a better understanding of how polar ice sheets may respond to accelerating global warming. The sensitivity of the Greenland ice sheet to polar amplification, changes in ocean heat transport, and deteriorating perennial sea ice conditions makes the Northeast Greenland margin a pertinent location with respect to understanding the impact of climate change on ice sheet instability and associated sea level rise. Throughout the Cenozoic, ocean heat fluxes toward and along Northeast Greenland have been controlled by water mass exchanges between the Arctic and Atlantic oceans. A key element here is the current flow through oceanic gateways, notably the Fram Strait and the Greenland–Scotland Ridge. To gain a long-term (million-year) perspective of ice sheet variability in this region, it is essential to understand the broader context of ice–ocean–tectonic interactions. Coupling between the ice sheet, the subsurface, the ocean, and sea ice are readily observable today in Northeast Greenland, but geological records to illuminate long-term trends and their interplay with other parts of the global climate system are lacking. Consequently, the NorthGreen workshop was organized by the Geological Survey of Denmark and Greenland in collaboration with Aarhus (Denmark) and Stockholm (Sweden) universities in November 2022 to develop mission-specific platform (MSP) proposals for drilling the Northeast Greenland margin under the umbrella of the MagellanPlus Workshop Series Programme of the European Consortium for Ocean Research Drilling (ECORD). Seventy-one participants representing a broad scientific community discussed key scientific questions and primary targets that could be addressed through scientific drilling in Northeast Greenland. Three pre-proposals were initiated during the workshop targeting Morris Jesup Rise, the Northeast Greenland continental shelf, and Denmark Strait. [ABSTRACT FROM AUTHOR]

Synthetic aperture radar (SAR) data from the Sentinel-1 (S1) mission with its high temporal and spatial resolution allows for an automated detection of lake ice break-up timings from surface backscatter differences across south (S), southwest (SW), and northwest (NW) Greenland (<71 ° N latitude) during the period 2017 to 2021. Median break-up dates of the 563 studied lakes range between 8 June and 10 July, with the earliest being in 2019 and the latest in 2018. There is a strong correlation between the break-up date and elevation, while a weak relationship with latitude and lake area could be observed. Lake-specific median break-up timings for 2017–2021 increase (i.e., are later) by 3 d per 100 m elevation gain. When assuming an earlier break-up timing of 8 d which corresponds to the observed median variability of ±8 d, the introduced excess energy due to a changing surface albedo from snow-covered ice surface to water translates to melting 0.4 ± 0.1 m thick ice at the melting point or heating up a water depth down to 35 ± 3 m by 1 K across the entire surface area of each respective lake. Upscaling the results to 100 486 lakes across the S, SW, and NW regions, which correspond to 64.5 % of all lakes or 62.1 % of the overall lake area in Greenland, yields an estimate of 1.8 × 10 6 TJ additional energy input. This translates to melting 5.8 Gt ice at the melting point or warming 432.3 Gt water by 1 K. [ABSTRACT FROM AUTHOR]

Seasonal and interannual changes in the Greenland Sea ice age composition during winter months for seven standard age categories are analyzed. To estimate the amount of the ice of different age categories, the digital archive of the region's ice charts for the period 1997–2022 produced by the Arctic and Antarctic Research Institute is used. Statistical insignificance of linear trends has been revealed for the main age categories of the Greenland Sea ice. A comparison of the obtained estimates of the variability in the amount of the ice of different age categories and earlier studies shows that the quantitative changes in the Greenland Sea ice age composition started before 1997. The data on the ice age composition over a 25-year observation period alone are insufficient to draw an unambiguous conclusion about a decrease in the Greenland Sea ice thickness. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet (GIS) and Atlantic Meridional Overturning Circulation (AMOC) are considered tipping elements in the climate system, where global warming exceeding critical threshold levels in forcing can lead to large-scale and nonlinear reductions in ice volume and overturning strength, respectively. The positive–negative feedback loop governing their interaction with a destabilizing effect on the AMOC due to ice loss and subsequent freshwater flux into the North Atlantic as well as a stabilizing effect of a net cooling around Greenland with an AMOC weakening may determine the long-term stability of both tipping elements. Here we explore the potential dynamic regimes arising from this positive–negative tipping feedback loop in a physically motivated conceptual model. Under idealized forcing scenarios we identify conditions under which different kinds of tipping cascades can occur: herein, we distinguish between overshoot/bifurcation tipping cascades, leading to tipping of both GIS and AMOC, and rate-induced tipping cascades, where the AMOC, despite not having crossed its own intrinsic tipping point, tips nonetheless due to the fast rate of ice loss from Greenland. The occurrence of these different cascades is affected by the ice sheet disintegration time and thus eventually by the imposed forcing and its timescales. Our results suggest that it is necessary not only to avoid surpassing the respective critical levels of the environmental drivers for the Greenland Ice Sheet and Atlantic Meridional Overturning Circulation, but also to respect safe rates of environmental change to mitigate potential domino effects. [ABSTRACT FROM AUTHOR]

Measurements of oceanographic characteristics, currents, and tides in the Fram Strait and the Greenland Sea are analyzed using historical measurement data. Measurements on moorings across the Fram Strait made it possible to estimate the speeds of the East Greenland and West Spitsbergen currents. It is shown that the speeds of these currents are low (within 10 cm/s) and the amplitudes of the speed of tidal currents are two or more times greater. Tidal currents are intensifying on the shelf of Greenland and Svalbard. Direct measurements are compared with the University of Oregon barotropic tide model. [ABSTRACT FROM AUTHOR]

Seismic velocity models provide important constraints on Greenland's deep structure, which, in turn, has profound implications for our understanding of the tectonic history of this region. However, the resolution of seismic models has been limited by a sparse network, particularly in northern and central Greenland. We address these limitations by generating new high-resolution Rayleigh-wave phase velocity maps encompassing Greenland and northeastern Canada by processing over three decades of teleseismic earthquake records and incorporating recently added stations in Greenland and Arctic Canada. These phase velocity maps are sensitive to structure from the lower crust down to the sub-lithospheric mantle (25–185 s period). We find significant heterogeneity and a strong correlation between isotropic and anisotropic seismic velocities with inferred geological structure. High seismic velocities associated with cratonic lithosphere are broadly divided into two regions, with a belt of reduced velocity spanning central Greenland, which we interpret as lithospheric erosion resulting from interaction between the Greenland continental keel and the Iceland plume. Within each region, we identify tectonic subdivisions that suggest fundamental differences between the blocks that make up Precambrian Greenland. In the south, the North Atlantic craton (NAC) has a high-velocity keel exhibiting anisotropic stratification. Between the NAC and the cratonic lithosphere further north, the Proterozoic Nagssugtoqidian orogenic belt shows a distinct signature of reduced seismic velocity to |$\sim$| 75 s period, but then appears to pinch out at depth. The northern Greenland lithosphere exhibits significant isotropic heterogeneity, with a distinct core of high velocities in the northwest (⁠|$\sim$| 55–75 s period) giving way to a set of distinct east-west trending high-velocity belts at longer periods. At all periods sensitive to the lithospheric mantle in this region, anisotropic fast orientations are E–W, consistent with a north–south Precambrian assembly of the Greenland shield. In contrast to the NAC, there is no evidence of anisotropic stratification in the northern part of the cratonic keel. Based on both isotropic and anisotropic phase-velocity anomalies, we suggest that the Phanerozoic Caledonian and Ellesmerian-Franklinian fold belts are relatively thin-skinned features onshore Greenland, though the Caledonian belt may have a stronger signature off the east coast. At the longest periods, a prominent low-velocity anomaly initially centred on Iceland migrates northwards and spreads beneath central-eastern Greenland. Coupled with NW–SE trending anisotropy, this feature is interpreted as the effect of mantle flow radiating outward from the Iceland plume and interacting with the eroded Greenland lithosphere. [ABSTRACT FROM AUTHOR]

Recent observations revealed the existence of ice slabs and aquifers on the Greenland ice sheet (GrIS). Both affect the ice sheet's hydrology: ice slabs facilitate runoff and aquifers modulate drainage to the bed. However, their climatic drivers and history remain unclear, as most observations cover only two decades. Here, we present a model simulation of the evolution of GrIS ice slabs and aquifers (1980–2020), evaluated using radar measurements. The results show that accumulation, melt and rain rates are good predictors for the spatial distribution of ice slabs and aquifers. Both features were already present in the late 1980s, and their extent remained relatively constant until the beginning of this century, after which increased melt led to their expansion. We show that almost any transect from the coast to the ice‐sheet interior will cross either an ice slab region, or an aquifer, or both. Plain Language Summary: The Greenland ice sheet is covered by a layer of old, porous snow called firn, which has the capacity to prevent surface melt runoff. The firn can however, form thick ice slabs or store large amounts of water underground. These change the fate of surface melt by respectively enhancing runoff and draining water to the ice sheet's bed. While observations have shown where these phenomena occur, it is still unknown what the exact atmospheric conditions are that lead to their formation. Here, we tackle this problem using computer simulations. Our results show that the amount of snowfall, melt and rain determine whether ice slabs or aquifers can form. The results also suggest that ice slabs and aquifers are more abundant than previously assumed, and that their observed expansion started only in the early 2000s, and they will continue to expand if the climate gets warmer and melt and rain become more abundant. Key Points: Accumulation, melt and rain are good predictors of ice slab and firn aquifer locationsThe expansion of modeled ice slabs and aquifers on the Greenland ice sheet started in the early 2000sWe show that, in between the ablation zone and the accumulation zone, there are always ice slabs and/or aquifers present [ABSTRACT FROM AUTHOR]

Dark (low albedo) surface ice on the Greenland Ice Sheet enhances melting and subsequent runoff, a major mass loss contributor during the ablation season. The accumulation of both biological (e.g., glacier ice algae) and abiotic (e.g., mineral dust) light‐absorbing particulates are important darkening factors, that are potentially influenced by the duration of snow‐free, bare ice (a phenological factor), and other geo‐topographical factors such as elevation, slope, aspect and the distance from the ice margin. Here, we present the first medium‐resolution (30 m) analysis of the phenological and geo‐topographical controls on the distribution of dark ice in SE and SW Greenland from statistical analysis of data derived from a harmonized satellite albedo product and ArcticDEM. The duration of bare ice primarily controls the distribution of dark surface ice, allowing for algae growth on inland ice surfaces in particular, whereas geo‐topographical factors are only secondary controls. Plain Language Summary: Meltwater that runs off the Greenland Ice Sheet is currently the largest single glaciological contributor to global sea level rise. Recent darkening of the surface ice has contributed significantly to increased melt rates. We used medium resolution (30m) satellite data to link ice surface darkness (albedo) to the topography (elevation, slope, and aspect) of the ice sheet surface, the distance from the ice margin, and the duration of bare or snow‐free surface ice. Our analysis shows that the duration of bare ice or snow‐free conditions was the most significant factor in driving the darkening of ice. We also found that the impact of topographic factors on the ice darkening varied depending on whether abiotic or biological processes were driving the darkening. This information helps us to better explain where the ice sheet is melting fastest and to predict the loss of ice mass, particularly considering the expected increase in the area of bare, snow‐free ice in our warming climate. Key Points: The duration of snow‐free, bare ice is the primary control on the distribution of dark iceLonger durations of bare ice (median = 40 days) are prerequisites for ice to become darkSurface slope and aspect influence darkening by modulating the duration of bare ice and glacier ice algae growth [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet is a key contributor to contemporary global sea level rise, but its long-term history and response to episodes of warming in Earth's geological past remain uncertain. The terrain covered by the ice sheet comprises ∼ 79 % of Greenland and ∼ 1.1 % of the Earth's land surface and contains geomorphological records that may provide valuable insights into past ice-sheet behaviour. Here we use ice surface morphology and radio-echo sounding data to identify ice-covered valleys within the highlands of southern and eastern Greenland and use numerical ice-sheet modelling to constrain the climatological and glaciological conditions responsible for valley incision. Our mapping reveals intricate subglacial valley networks with morphologies that are indicative of substantial glacial modification of an inherited fluvial landscape, yet many of these valleys are presently situated beneath cold-based, slow-moving (i.e. non-erosive) ice. We use the morphology of the valleys and our simple ice-sheet model experiments to infer that incision likely occurred beneath erosive mountain valley glaciers during one or more phases of Greenland's glacial history when ice was restricted to the southern and eastern highlands and when Greenland's contribution to barystatic sea level was up to + 7 m relative to today. We infer that this valley incision primarily occurred prior to the growth of a continental-scale ice sheet, most likely during the late Miocene (ca. 7–5 Ma) and/or late Pliocene (ca. 3.6–2.6 Ma). Our findings therefore provide new data-based constraints on early Greenland Ice Sheet extent and dynamics that can serve as valuable boundary conditions in models of regional and global palaeoclimate during past warm periods that are important analogues for climate change in the 21st century and beyond. [ABSTRACT FROM AUTHOR]

We present an ice-penetrating radar data set acquired over the Greenland ice sheet by aircraft during the years 1971, 1972, 1974, 1978, and 1979. The data set comprises over 177,000 km of flight lines and contains a wealth of information on the state of the Greenland ice sheet including information on ice thickness and englacial properties. During data collection in the 1970s, the data were recorded on optical film rolls and in this manuscript, we document the digitization of these film rolls and their associated geographical information. Our digitization of the data enables interaction with and analysis of the data and facilitates comparison with modern-day radar observations. The complete data set in full resolution is available at the Stanford Digital Repository (https://doi.org/10.25740/wm135gp2721) with the associated technical reports (Karlsson et al., 2023). Part of the data set has been quality checked and is available as low-resolution JPG files at The Technical University of Denmark's data repository, DTU Data, with associated geographical information (DOI pending). Here, we release the full data sets to enable the larger community to access and interact with the data. [ABSTRACT FROM AUTHOR]

Microbial blooms colonize the Greenland Ice Sheet bare ice surface during the ablation season and significantly reduce its albedo. On the ice surface, microbes are exposed to high levels of irradiance, freeze–thaw cycles, and low nutrient concentrations. It is well known that microorganisms secrete metabolites to maintain homeostasis, communicate with other microorganisms, and defend themselves. Yet, the exometabolome of supraglacial microbial blooms, dominated by the pigmented glacier ice algae Ancylonema alaskanum and Ancylonema nordenskiöldii, remains thus far unstudied. Here, we use a high‐resolution mass spectrometry‐based untargeted metabolomics workflow to identify metabolites in the exometabolome of microbial blooms on the surface of the southern tip of the Greenland Ice Sheet. Samples were collected every 6 h across two diurnal cycles at 5 replicate sampling sites with high similarity in community composition, in terms of orders and phyla present. Time of sampling explained 46% (permutational multivariate analysis of variance [PERMANOVA], pseudo‐F = 3.7771, p = 0.001) and 27% (PERMANOVA, pseudo‐F = 1.8705, p = 0.001) of variance in the exometabolome across the two diurnal cycles. Annotated metabolites included riboflavin, lumichrome, tryptophan, and azelaic acid, all of which have demonstrated roles in microbe–microbe interactions in other ecosystems and should be tested for potential roles in the development of microbial blooms on bare ice surfaces. [ABSTRACT FROM AUTHOR]

Surface melt on the Greenland ice sheet has been increasing in intensity and extent over the last decades due to Arctic atmospheric warming. Surface melt depends on the surface energy balance, which includes the atmospheric forcing but also the thermal budget of the snow, firn and ice near the ice sheet surface. The temperature of the ice sheet subsurface has been used as an indicator of the thermal state of the ice sheet's surface. Here, we present a compilation of 4612 measurements of firn and ice temperature at 10 m below the surface (T10m) across the ice sheet, spanning from 1912 to 2022. The measurements are either instantaneous or monthly averages. We train an artificial neural network model (ANN) on 4597 of these point observations, weighted by their relative representativity, and use it to reconstruct T10m over the entire Greenland ice sheet for the period 1950–2022 at a monthly timescale. We use 10-year averages and mean annual values of air temperature and snowfall from the ERA5 reanalysis dataset as model input. The ANN indicates a Greenland-wide positive trend of T10m at 0.2 ∘ C per decade during the 1950–2022 period, with a cooling during 1950–1985 (-0.4 ∘ C per decade) followed by a warming during 1985–2022 (+0.7 ∘ per decade). Regional climate models HIRHAM5, RACMO2.3p2 and MARv3.12 show mixed results compared to the observational T10m dataset, with mean differences ranging from -0.4 ∘ C (HIRHAM) to 1.2 ∘ C (MAR) and root mean squared differences ranging from 2.8 ∘ C (HIRHAM) to 4.7 ∘ C (MAR). The observation-based ANN also reveals an underestimation of the subsurface warming trends in climate models for the bare-ice and dry-snow areas. The subsurface warming brings the Greenland ice sheet surface closer to the melting point, reducing the amount of energy input required for melting. Our compilation documents the response of the ice sheet subsurface to atmospheric warming and will enable further improvements of models used for ice sheet mass loss assessment and reduce the uncertainty in projections. [ABSTRACT FROM AUTHOR]

Numerical modeling of ice sheet motion and hence projections of global sea level rise require information about the evolving subglacial environment, which unfortunately remains largely unknown due to its difficulty of access. Here we advance such subglacial observations by reporting multi‐year observations of seismic tremor likely associated with glacier sliding at Helheim Glacier. This association is confirmed by correlation analysis between tremor power and multiple environmental forcings on different timescales. Variations of the observed tremor power indicate that different factors affect glacial sliding on different timescales. Effective pressure may control glacial sliding on long (seasonal/annual) timescales, while tidal forcing modulates the sliding rate and tremor power on short (hourly/daily) timescales. Polarization results suggest that the tremor source comes from an upstream subglacial ridge. This observation provides insights on how different factors should be included in ice sheet modeling and how their timescales of variability play an essential role. Plain Language Summary: The Greenland Ice Sheet has lost ice increasingly in the past few decades, contributing significantly to global sea level rise. However, large uncertainties remain in computer simulations of such ice mass loss and hence sea level rise. This uncertainly is mainly attributed to the lack of information on what is happening at the bottom of the ice sheet and how that affects ice sheet movement. For example, how much water is there and whether it is in isolated pockets or is well distributed as a thin sheet can have an important effect on ice movement. In this paper, we observe small ground motions that are generated during glacier movement at a marine‐terminating glacier in Greenland. By analyzing the energy variation of the small ground motions over multiple years as well as observing the reasons that cause the variations, we learn that subglacial water pressure may control ice flow speeds on a seasonal/annual scale, and that ocean tides can change ice flow on an hourly/daily scale. This observation provides important constraints for computer simulations of future sea level rise in terms of the impacting factors and their respective timescales. Key Points: Multi‐year seismic records reveal glacial slip and tremor variability at Helheim GlacierTremor correlates with high effective pressure at tidal timescales, opposite to the expectation at longer timescalesThe tremor source points to an upstream subglacial ridge [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet (GrIS) has been losing mass since the 1990s as a direct consequence of rising temperatures and has been projected to continue to lose mass at an accelerating pace throughout the 21st century, making it one of the largest contributors to future sea-level rise. The latest Coupled Model Intercomparison Project Phase 6 (CMIP6) models produce a greater Arctic amplification signal and therefore also a notably larger mass loss from the GrIS when compared to the older CMIP5 projections, despite similar forcing levels from greenhouse gas emissions. However, it is also argued that the strength of regional factors, such as melt–albedo feedbacks and cloud-related feedbacks, will partly impact future melt and sea-level rise contribution, yet little is known about the role of these regional factors in producing differences in GrIS surface melt projections between CMIP6 and CMIP5. In this study, we use high-resolution (15 km) regional climate model simulations over the GrIS performed using the Modèle Atmosphérique Régional (MAR) to physically downscale six CMIP5 Representative Concentration Pathway (RCP) 8.5 and five CMIP6 Shared Socioeconomic Pathway (SSP) 5-8.5 extreme high-emission-scenario simulations. Here, we show a greater annual mass loss from the GrIS at the end of the 21st century but also for a given temperature increase over the GrIS, when comparing CMIP6 to CMIP5. We find a greater sensitivity of Greenland surface mass loss in CMIP6 centred around summer and autumn, yet the difference in mass loss is the largest during autumn with a reduction of 27.7 ± 9.5 Gt per season for a regional warming of +6.7 ∘ C and 24.6 Gt per season more mass loss than in CMIP5 RCP8.5 simulations for the same warming. Assessment of the surface energy budget and cloud-related feedbacks suggests a reduction in high clouds during summer and autumn – despite enhanced cloud optical depth during autumn – to be the main driver of the additional energy reaching the surface, subsequently leading to enhanced surface melt and mass loss in CMIP6 compared to CMIP5. Our analysis highlights that Greenland is losing more mass in CMIP6 due to two factors: (1) a (known) greater sensitivity to greenhouse gas emissions and therefore warmer temperatures and (2) previously unnotified cloud-related surface energy budget changes that enhance the GrIS sensitivity to warming. [ABSTRACT FROM AUTHOR]

There is currently poor scientific agreement on whether the ice–bed interface is frozen or thawed beneath approximately one third of the Greenland ice sheet. This disagreement in basal thermal state results, at least partly, from differences in the subglacial geothermal heat-flow basal boundary condition used in different ice-flow models. Here, we employ seven widely used Greenland geothermal heat-flow maps in 10 000-year spin-ups of the Community Ice Sheet Model (CISM). We perform two spin-ups: one nudged toward thickness observations and the other unconstrained. Across the seven heat-flow maps, and regardless of unconstrained or nudged spin-up, the spread in basal ice temperatures exceeds 10 ∘ C over large areas of the ice–bed interface. For a given heat-flow map, the thawed-bed ice-sheet area is consistently larger under unconstrained spin-ups than nudged spin-ups. Under the unconstrained spin-up, thawed-bed area ranges from 33.5 % to 60.0 % across the seven heat-flow maps. Perhaps counterintuitively, the highest iceberg calving fluxes are associated with the lowest heat flows (and vice versa) for both unconstrained and nudged spin-ups. These results highlight the direct, and non-trivial, influence of the heat-flow boundary condition on the simulated equilibrium thermal state of the ice sheet. We suggest that future ice-flow model intercomparisons should employ a range of basal heat-flow maps, and limit direct intercomparisons with simulations using a common heat-flow map. [ABSTRACT FROM AUTHOR]

Subglacial environments provide conditions suitable for the microbial production of methane, an important greenhouse gas, which can be released from beneath the ice as a result of glacial melting. High gaseous methane emissions have recently been discovered at Russell Glacier, an outlet of the southwestern margin of the Greenland Ice Sheet, acting not only as a potential climate amplifier but also as a substrate for methane consuming microorganisms. Here, we describe the composition of the microbial assemblage exported in meltwater from the methane release hotspot at Russell Glacier and its changes over the melt season and as it travels downstream. We found that a substantial part (relative abundance 27.2% across the whole dataset) of the exported assemblage was made up of methylotrophs and that the relative abundance of methylotrophs increased as the melt season progressed, likely due to the seasonal development of the glacial drainage system. The methylotrophs were dominated by representatives of type I methanotrophs from the Gammaproteobacteria; however, their relative abundance decreased with increasing distance from the ice margin at the expense of type II methanotrophs and/or methylotrophs from the Alphaproteobacteria and Betaproteobacteria. Our results show that subglacial methane release hotspot sites can be colonized by microorganisms that can potentially reduce methane emissions. [ABSTRACT FROM AUTHOR]

The microbiome of Greenland Ice Sheet supraglacial habitats is still underinvestigated, and as a result there is a lack of representative genomes from these environments. In this study, we investigated the supraglacial microbiome through a combination of culturing-dependent and -independent approaches. We explored ice, cryoconite, biofilm, and snow biodiversity to answer: (1) how microbial diversity differs between supraglacial habitats, (2) if obtained bacterial genomes reflect dominant community members, and (3) how culturing versus high throughput sequencing changes our observations of microbial diversity in supraglacial habitats. Genomes acquired through metagenomic sequencing (133 high-quality MAGs) and whole genome sequencing (73 bacterial isolates) were compared to the metagenome assemblies to investigate abundance within the total environmental DNA. Isolates obtained in this study were not dominant taxa in the habitat they were sampled from, in contrast to the obtained MAGs. We demonstrate here the advantages of using metagenome SSU rRNA genes to reflect whole-community diversity. Additionally, we demonstrate a proof-of-concept of the application of in situ culturing in a supraglacial setting. [ABSTRACT FROM AUTHOR]

It is highly uncertain how the humidity flux between the snow surface and the atmosphere contributes to the surface mass balance (SMB) of the interior Greenland Ice Sheet (GrIS). Due to sparse observations, evaluations of the simulated humidity flux are limited. Model-based estimates of the humidity flux contribution to the SMB are, therefore, unconstrained and even disagree in magnitude and sign. In this study, we evaluate the regional climate model MAR at the EGRIP (East Greenland Ice-Core Project) site in the accumulation zone of the GrIS. We use a combined dataset of continuous one-level bulk estimates of the humidity flux covering the period of May 2016–August 2019 and eddy-resolving eddy-covariance humidity flux measurements from all four summer seasons. In summer, we document a bias of too little sublimation (- 1.3 W m −2 , - 1.65 mm w.e.) caused by a cold bias in both air and surface temperature, leading to a reduced humidity gradient. In winter, MAR overestimates vapor deposition by about 1 order of magnitude. This is a consequence of an overestimated temperature gradient in too stable atmospheric conditions compared to observations. Both systematic errors cause a large discrepancy in the annual net humidity flux between the model and observations of - 9 mm w.e. yr −1. Remarkably, the simulated net annual humidity flux contributes positively to the SMB, contrary to observations documenting a net sublimation flux. We correct the systematic errors by applying a simple but effective correction function to the simulated latent heat flux. Using this correction, we find that 5.1 % of the annual mass gain at the EGRIP site sublimates again, and 4.3 % of the total mass gain is deposited vapor from the near-surface air. The estimated net humidity flux contribution to the annual SMB is about - 1 % (net sublimation) compared to + 5.6 % for the uncorrected simulation. In summer, the corrected MAR simulation shows that vapor deposition accounts for 9.6 % of the total mass gain and that 31 % of the total mass gain at the EGRIP site sublimates again. The net fluxes contribute to - 32 % of the summer SMB. These results demonstrate that the humidity flux is a major driver of the summer SMB in the accumulation zone of the GrIS and highlight that even small changes could increase its importance for the annual SMB in a warming climate. [ABSTRACT FROM AUTHOR]

The negative surface mass balance of glaciers and ice caps under a warming climate impacts local ecosystems, influencing the volume and timing of water flow in local catchments while also contributing to global sea level rise. Peripheral glaciers distinct to the Greenland ice sheet respond faster to climate change than the main ice sheet. Accurate assessment of surface mass balance depends on in situ observations of near-surface climate and ice ablation, but very few in situ observations of near-surface climate and ice ablation are freely available for Greenland's peripheral glaciers. The transect of three automated weather and ablation stations on the peripheral A. P. Olsen Ice Cap in northeast Greenland is an example of these much needed data. The transect has been monitored since 2008, and in 2022, the old weather and ablation stations were replaced by a new standardized setup. In order to ensure comparable data quality of the old and new monitoring station setups, it is necessary to re-evaluate the data collected between 2008 and 2022. This paper presents the fully reprocessed near-surface climate and ablation data from the A. P. Olsen Ice Cap transect from 2008 to 2022, with a focus on data quality and the usability in ice ablation process studies. The usability of the data is exemplified by the data in an energy balance melt model for two different years. We show that the inherent uncertainties in the data result in an accurate reproduction of ice ablation for just one of the two years. A transect of three automatic ablation and weather stations of this length is unique to Greenland's peripheral glaciers, and it has a broad scale of usage from input to climate reanalysis and detailed surface ablation studies. The dataset can be downloaded at 10.22008/FK2/X9X9GN. [ABSTRACT FROM AUTHOR]

We document the isotopic evolution of near-surface snow at the East Greenland Ice Core Project (EastGRIP) ice core site in northeast Greenland using a time-resolved array of 1 m deep isotope (δ18O , δD) profiles. The snow profiles were taken from May–August during the 2017–2019 summer seasons. An age–depth model was developed and applied to each profile, mitigating the impacts of stratigraphic noise on isotope signals. Significant changes in deuterium excess (d) are observed in surface snow and near-surface snow as the snow ages. Decreases in d of up to 5 ‰ occur during summer seasons after deposition during two of the three summer seasons observed. The d always experiences a 3 ‰–5 ‰ increase after aging 1 year in the snow due to a broadening of the autumn d maximum. Models of idealized scenarios coupled with prior work indicate that the summertime post-depositional changes in d (Δd) can be explained by a combination of surface sublimation, forced ventilation of the near-surface snow down to 20–30 cm , and isotope-gradient-driven diffusion throughout the column. The interannual Δd is also partly explained with isotope-gradient-driven diffusion, but other mechanisms are at work that leave a bias in the d record. Thus, d does not just carry information about source-region conditions and transport history as is commonly assumed, but also integrates local conditions into summer snow layers as the snow ages through metamorphic processes. Finally, we observe a dramatic increase in the seasonal isotope-to-temperature sensitivity, which can be explained solely by isotope-gradient-driven diffusion. Our results are dependent on the site characteristics (e.g., wind, temperature, accumulation rate, snow properties) but indicate that more process-based research is necessary to understand water isotopes as climate proxies. Recommendations for monitoring and physical modeling are given, with special attention to the d parameter. [ABSTRACT FROM AUTHOR]

We analyze the temperature time series of the EPICA Dome C ice cores in Antarctica and of the Greenland project, Summit, with durations of 800 000 and 248 000 years, respectively, with a recent mathematical tool defined through the Fourier phases of the series, known as the J -index. This data driven index can differentiate between purely random dynamics and dynamics with a deterministic component. It is sensitive to nonlinear components and robust to the presence of noise. Our J -index data analysis shows that both Greenland and Antarctica climatic fluctuations possess deterministic traits and suggests the presence of an underlying nonlinear dynamics. Furthermore, in both regions, it reveals the simultaneous occurrence of an important global event known as the "Pelukian transgression." For Antarctica, it also detects the marine isotopic stage 11. Additionally, our calculation of the time series Hurst exponents and our detrended fluctuation analysis show the presence of long-range persistent correlations for Antarctica and anti-persistent correlations for Greenland. For the latter case, our fractal dimension determinations are indicative of a more complex climatic dynamics in Greenland with respect to Antarctica. Our results are encouraging for further development of climate variability deterministic models for these regions. [ABSTRACT FROM AUTHOR]

The mass balance of the Greenland Ice Sheet is strongly influenced by the dynamics of its outlet glaciers. Therefore, it is of paramount importance to accurately and continuously monitor these glaciers, especially the variation in their frontal positions. A temporally comprehensive parameterization of glacier calving is essential for understanding dynamic changes and constraining ice sheet modeling. However, many current calving front records are limited in terms of temporal resolution as they rely on manual delineation, which is laborious and not appropriate considering the increasing amount of satellite imagery available. In this contribution, we address this problem by applying an automated method to extract calving fronts from optical satellite imagery. The core of this workflow builds on recent advances in the field of deep learning while taking full advantage of multispectral input information. The performance of the method is evaluated using three independent test datasets. For the three datasets, we calculate mean delineation errors of 61.2, 73.7, and 73.5 m, respectively. Eventually, we apply the technique to Landsat-8 imagery. We generate 9243 calving front positions across 23 outlet glaciers in Greenland for the period 2013–2021. Resulting time series not only resolve long-term and seasonal signals but also resolve subseasonal patterns. We discuss the implications for glaciological studies and present a first application for analyzing the effect of bedrock topography on calving front variations. Our method and derived results represent an important step towards the development of intelligent processing strategies for glacier monitoring, opening up new possibilities for studying and modeling the dynamics of Greenland's outlet glaciers. [ABSTRACT FROM AUTHOR]

We present the first continuous mooring records of the West Greenland Coastal Current (WGCC), a conduit of fresh, buoyant outflow from the Arctic Ocean and the Greenland Ice Sheet. Nearly two years of temperature, salinity, and velocity data from 2018 to 2020 demonstrate that the WGCC on the southwest Greenland shelf is a well-formed current distinct from the shelfbreak jet but exhibits strong chaotic variability in its lateral position on the shelf, ranging from the coastline to the shelf break (50 km offshore). We calculate the WGCC volume and freshwater transports during the 35% of the time when the mooring array fully bracketed the current. During these periods, the WGCC remains as strong (0.83 ± 0.02 Sverdrups; 1 Sv ≡ 106 m3 s−1) as the East Greenland Coastal Current (EGCC) on the southeast Greenland shelf (0.86 ± 0.05 Sv) but is saltier than the EGCC and thus transports less liquid freshwater (30 × 10−3 Sv in the WGCC vs 42 × 10−3 Sv in the EGCC). These results indicate that a significant portion of the liquid freshwater in the EGCC is diverted from the coastal current as it rounds Cape Farewell. We interpret the dominant spatial variability of the WGCC as an adjustment to upwelling-favorable wind forcing on the West Greenland shelf and a separation from the coastal bathymetric gradient. An analysis of the winds near southern Greenland supports this interpretation, with nonlocal winds on the southeast Greenland shelf impacting the WGCC volume transport more strongly than local winds over the southwest Greenland shelf. [ABSTRACT FROM AUTHOR]

Satellite data have revealed that the Greenland and Antarctic Ice Sheets are changing rapidly due to warming air and ocean temperatures. Crucially, Earth Observations can now be used to measure ice sheet mass balance at the continental scale, which can help reduce uncertainties in the ice sheets' past, present, and future contributions to global mean sea level. The launch of satellite missions dedicated to the polar regions led to great progress towards a better assessment of the state of the ice sheets, which, in combination with ice sheet models, have furthered our understanding of the physical processes leading to changes in the ice sheets' properties. There is now a three-decade-long satellite record of Antarctica and Greenland mass changes, and new satellite missions are planned to both continue this record and further develop our observational capabilities, which is critical as the ice sheets remain the most uncertain component of future sea-level rise. In this paper, we review the mechanisms leading to ice sheets' mass changes and describe the state of the art of the satellite techniques used to monitor Greenland's and Antarctica's mass balance, providing an overview of the contributions of Earth Observations to our knowledge of these vast and remote regions. [ABSTRACT FROM AUTHOR]

As the main contributor to global sea-level rise, the Greenland Ice Sheet (GrIS) has undergone significant mass change over the last two decades. The satellite mission of GRACE (Gravity Recovery And Climate Experiment) and its follow-on mission (GRACE-FO) provide accurate observations but low-spatial resolution. In contrast, satellite altimetry provides observations at a high-spatial resolution but with large uncertainties, limiting the understanding of glacier-scale mass change. To derive accurate and high-spatial resolution mass change estimates from GRACE/GRACE-FO observations, we present a novel constraint mascon method in which the regularization matrix is constructed with the signal variances from satellite altimetry. Based on the proposed method, we derive a series of high-resolution (25 km × 25 km) monthly mascon solutions from 2002 April to September. The glacier-scale estimates from the input–output method agree better with those from our mascon solutions than those from the global mascons of CSR (Center for Space Research, the University of Texas), JPL (Jet Propulsion Laboratory), and GSFC (Goddard Space Flight Center), with a higher linear regression coefficient of 0.71. Benefitting from the greatly improved spatial resolution, our estimates provide the first accurate monthly glacier-scale mass change estimates from GRACE/GRACE-FO observations over the GrIS, to our knowledge. The results show that 20 of the 260 glaciers contributed to more than 42 per cent of the ice loss in the GrIS from 2002 to 2022. Most strikingly, the mass loss of Jakobshavn Isbrae was the most significant at –18.7 ± 0.05 Gt yr−1, accounting for 7.4 per cent of the total in the GrIS during the study period. Furthermore, we find that the SMB (surface mass balance) and ice-dynamics-related mass changes contribute nearly equally to the observed mass changes, but the corresponding spatiotemporal characteristics differ. SMB contributed the most to the mass change in 2019, while ice dynamics played the most dominant role in 2018. Moreover, the SMB effect is widespread, and ice-dynamics-related mass loss is mainly concentrated in five small areas. [ABSTRACT FROM AUTHOR]

The occurrence of Supraglacial Lakes (SGLs) may influence the signals acquired with microwave radiometers, which may result in a degree of uncertainty when employing microwave radiometer data for the detection of surface melt. Accurate monitoring of surface melting requires a reasonable assessment of this uncertainty. However, there is a scarcity of research in this field. Therefore, in this study, we computed surface melt in the vicinity of Automatic Weather Stations (AWSs) by employing Defense Meteorological Satellite Program (DMSP) Ka-band data and Soil Moisture and Ocean Salinity (SMOS) satellite L-band data and extracted SGL pixels by utilizing Sentinel-2 data. A comparison between surface melt results derived from AWS air temperature estimates and those obtained with remote sensing inversion in the two different bands was conducted for sites below the mean snowline elevation during the summers of 2016 to 2020. Compared with sites with no SGLs, the commission error (CO) of DMSP morning and evening data at sites where these water bodies were present increased by 36% and 30%, respectively, and the number of days with CO increased by 12 and 3 days, respectively. The omission error (OM) of SMOS morning and evening data increased by 33% and 32%, respectively, and the number of days with OM increased by 17 and 21 days, respectively. Identifying the source of error is a prerequisite for the improvement of surface melt algorithms, for which this study provides a basis. [ABSTRACT FROM AUTHOR]

Introduction: Noncommunicable diseases (NCDs) constitute a massive global burden and are the leading cause of death and disability worldwide. In Greenland, the prevalence of NCDs has historically been low. However, during the past approximately 70 years, life circumstances have changed dramatically resulting in increased life expectancy. Today, the proportion of inhabitants in Greenland ≥65 years has nearly tripled since the 1980s, and the prevalence of obesity and diabetes has increased rapidly within the past decades. The aim of this study was to describe the burden of selected NCDs in a primary care setting in Nuuk and compare it to a modern westernized suburban general practice in Denmark. Methods: The study was performed as a cross sectional register-based study using data extracted from the electronic medical records (EMR) based on diagnosis codes from inhabitants living in Nuuk, Greenland, and a suburb in Denmark. Estimates of prevalence were age-standardized to the WHO world standard population. Results: In both Nuuk and the Danish suburb, the highest prevalence was observed for hypertension (13.2% for both populations), followed by asthma (4.4 and 9.5%, respectively) and diabetes (4.3 and 2.9%, respectively). The age-standardized prevalences of diabetes, COPD, atrial fibrillation, and heart failure, were significantly higher in Nuuk, while seven NCDs including asthma, ischemic heart disease, arthritis urica, psoriasis, hyperthyreosis, hypothyreosis and osteoporosis were significantly higher in the Danish suburb. Conclusion: In contrast to the disease pattern observed in Greenland in the last century, the prevalence of diagnosed NCDs in Nuuk is no longer rare. Thus, the overall prevalence of NCDs in the population of Nuuk is now comparable to or even higher than in the suburb in Denmark. This calls for increased focus on all NCDs in the primary healthcare system in Greenland and adaption of the primary healthcare services to a changed disease spectrum. [ABSTRACT FROM AUTHOR]

Heavily pigmented glacier ice algae Ancylonema nordenskiöldii and Ancylonema alaskanum (Zygnematophyceae, Streptophyta) reduce the bare ice albedo of the Greenland Ice Sheet, amplifying melt from the largest cryospheric contributor to eustatic sea-level rise. Little information is available about glacier ice algae interactions with other microbial communities within the surface ice environment, including fungi, which may be important for sustaining algal bloom development. To address this substantial knowledge gap and investigate the nature of algal-fungal interactions, an ex situ co-cultivation experiment with two species of fungi, recently isolated from the surface of the Greenland Ice Sheet (here proposed new species Penicillium anthracinoglaciei Perini, Frisvad and Zalar, Mycobank (MB 835602), and Articulospora sp.), and the mixed microbial community dominated by glacier ice algae was performed. The utilization of the dark pigment purpurogallin carboxylic acid-6-O-β-D-glucopyranoside (C18H18O12) by the two fungi was also evaluated in a separate experiment. P. anthracinoglaciei was capable of utilizing and converting the pigment to purpurogallin carboxylic acid, possibly using the sugar moiety as a nutrient source. Furthermore, after 3 weeks of incubation in the presence of P. anthracinoglaciei, a significantly slower decline in the maximum quantum efficiency (Fv/Fm, inverse proxy of algal stress) in glacier ice algae, compared to other treatments, was evident, suggesting a positive relationship between these species. Articulospora sp. did uptake the glycosylated purpurogallin, but did not seem to be involved in its conversion to aglycone derivative. At the end of the incubation experiments and, in conjunction with increased algal mortality, we detected a substantially increasing presence of the zoosporic fungi Chytridiomycota suggesting an important role for them as decomposers or parasites of glacier ice algae. [ABSTRACT FROM AUTHOR]

Atmospheric blocking is a phenomenon that can lead to extreme weather events over a large region, yet its causes are not fully understood. Global climate models show limitations in representing Northern Hemisphere blocking, especially its frequency, and decadal variability in Greenland blocking in summer in the recent decades. In this study we evaluate the ability of high-resolution (HighResMIP) Earth System Models (ESMs) to simulate summer blocking over the Greenland area, using different but complementary methods to describe the characteristics of blocking. We find that the HighResMIP ensemble can reproduce the spatial pattern of Greenland blocking events, albeit with systematic biases, and capture the relative frequencies of the main blocking patterns: namely the wave breaking structure, North Atlantic ridge, and omega-type blocking. However, the HighResMIP ensemble fails to simulate the observed temporal variations of Greenland blocking index (GB2) and the extremely high values of daily GB2 observed in recent decades. In addition, we do not find clearly superior representation of blocking features from higher-resolution in HighResMIP models compared with lower-resolution models. We also find large sea surface temperature (SST) biases over the North Atlantic and seas surrounding Greenland, and biases in moisture transport over the North Atlantic toward Greenland, especially over the western flank of blocking areas, which might together contribute to model biases in the representation of blocking magnitude. [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 recent decades, the Greenland ice sheet has been losing mass through glacier retreat and ice flow acceleration. This mass loss is linked with variations in submarine melt, yet existing ocean models are either coarse global simulations focused on decadal‐scale variability or fine‐scale simulations for process‐based investigations. Here, we unite these scales with a framework to downscale from a global state estimate (15 km) into a regional model (3 km) that resolves circulation on the continental shelf. We further downscale into a fjord‐scale model (500 m) that resolves circulation inside fjords and quantifies melt. We demonstrate this approach in Scoresby Sund, East Greenland, and find that interannual variations in submarine melt at Daugaard‐Jensen glacier induced by ocean temperature variability are consistent with the decadal changes in glacier ice dynamics. This study provides a framework by which coarse‐resolution models can be refined to quantify glacier submarine melt for future ice sheet projections. Plain Language Summary: Over the past several decades, the Greenland ice sheet has been losing ice and contributing to sea‐level rise. About half of this ice loss is induced by melt that occurs where glaciers meet the ocean. Using coarse‐scale ocean models that simulate circulation around the globe, previous studies have noted a strong link between ocean temperature and enhanced glacier ice loss. However, due to the small scale of Greenland's fjords, coarse models are unable to directly quantify circulation in these fjords and melt on submerged glaciers. In this study, we develop a new framework to "zoom in" on a fjord, using high‐resolution models driven by larger coarse‐resolution models. In this approach, we simulate melt on one of Greenland's biggest glaciers and find that periods of higher melt coincide with more ice loss as observed from satellites. Since this framework is adaptable to other regions, it could also be used to simulate melt on other glaciers and support estimates of future sea‐level rise. Key Points: Subsurface temperature variability is simulated in a narrow fjord network using regional models downscaled from a global state estimateModeled increases in ocean melt at Daugaard‐Jensen glacier coincide with the onset of acceleration in 2005 and retreat and thinning in 2011Model variations in shelf‐to‐fjord ocean properties match with observations, providing a basis to estimate ocean forcing in ice projections [ABSTRACT FROM AUTHOR]