Ice-marginal lakes form at the edge of the Greenland Ice Sheet and its surrounding peripheral glaciers and ice caps (PGIC), where outflowing glacial meltwater is trapped by a moraine, or by the ice itself, and create a reservoir that is in contact with the adjacent ice. While glacial meltwater is typically assumed to flow directly into the ocean, ice-marginal lakes temporarily store a portion of this runoff, influencing glacier dynamics and ablation, ecosystems, and downstream hydrology. Their presence, and change in abundance and size, remain under-represented in projections of sea level change and glacier mass loss. Here, we present an eight-year (2016–2023) inventory of 2918 automatically classified ice-marginal lakes (≧0.05 km2) across Greenland, tracking changes in lake abundance, surface extent, and summer surface temperature over time. Fluctuations in lake abundance were most pronounced at the north (22 %) and northeast (14 %) PGIC margins and the southwest Ice Sheet margin (8 %). Over the study period, an increase in surface lake area was evident at 283 lakes, a decreasing trend was evident at 240 lakes, and 1373 remained stable (±0.05 km2). The northeast region contained the largest lakes, with a median size of 0.40 km2 at the ice sheet margin and 0.24 km2 at PGIC margins. Average summer surface temperatures fluctuated between 3.8 °C (2018) and 5.3 °C (2023), with spatial and temporal trends identified with possible links to lake setting and size. Validation against manually identified lakes showed 64 % agreement, yielding an error estimate of −809 lakes (36 %), while lake area uncertainty was ±5 %. Surface temperature estimates showed strong agreement with in situ measurements (r2=0.87, RMSE =1.68 °C, error ±1.2 °C). This dataset provides a crucial foundation for quantifying meltwater storage at ice margins and refining sea level contribution projections while supporting research on glacier-lake interactions, Arctic ecology, and environmental management. The inventory series is openly accessible on the GEUS Dataverse (https://doi.org/10.22008/FK2/MBKW9N, How et al., 2025) with full metadata and documentation, and a reproducible processing workflow.

Indoor radon and its decay products are the primary sources of the population's exposure to background ionizing radiation. Radon decay products are one of the leading causes of lung cancer, with a higher lung cancer risk for smokers due to the synergistic effects of radon decay products and cigarette smoking. A total of 459 year-long radon measurements in 257 detached and semi-detached residential homes in southwest and south Greenland were carried out, and a dose assessment for adults was performed. The annual arithmetic and geometric means of indoor radon concentrations was 10.5 +/- 0.2 Bq m- 3 and 8.0 +/- 2.3 Bq m- 3 in Nuuk, 139.0 +/- 1.0 Bq m- 3 and 97.3 +/- 2.1 Bq m- 3 in Narsaq, and 42.1 +/- 0.7 Bq m- 3 and 22.0 +/- 3.1 Bq m- 3 in Qaqortoq. Arithmetic and geometric mean radon concentration of 79.0 Bq m- 3 and 50.3 Bq m- 3 were estimated for adult, person-weighted living in south Greenland. The total number of detached and semi-detached residential homes where indoor radon is exceeding 100 Bq m- 3, 200 Bq m- 3, and 300 Bq m- 3 is 37 homes (15.0%), 13 homes (5.2%), and 8 homes (3.2%), respectively. A positive correlation between indoor air radon concentrations and underlying geology was observed.The indoor radon contribution to the annual inhalation effective dose to an average adult was 0.5 mSv in Nuuk, 6.5 mSv in Narsaq, 2.0 mSv in Qaqortoq, and 4.0 mSv for south Greenland adult person weighted. The estimated annual average effective dose to adults in Narsaq is higher than the world's average annual effective dose of 1.3 mSv due to inhalation of indoor radon. Cost-efficient mitigation methods exist to reduce radon in existing buildings, and to prevent radon entry into new buildings.

The Greenland Ice Sheet (GrIS) is losing mass as the climate warms through both increased meltwater runoff and ice discharge at marine-terminating sectors. At the ice sheet surface, meltwater runoff forms a dynamic supraglacial hydrological system which includes stream and river networks and large supraglacial lakes (SGLs). Streams and rivers can route water into crevasses or into supraglacial lakes with crevasses underneath, both of which can then hydrofracture to the ice sheet base, providing a mechanism for the surface meltwater to access the bed. Understanding where, when, and how much meltwater is transferred to the bed is important because variability in meltwater supply to the bed can increase ice flow speeds, potentially impacting the hypsometry of the ice sheet in grounded sectors, and iceberg discharge to the ocean. Here we present a new, physically based, supraglacial hydrology model for the GrIS that is able to simulate (a) surface meltwater routing and SGL filling; (b) rapid meltwater drainage to the ice sheet bed via the hydrofracture of surface crevasses both in and outside of SGLs; (c) slow SGL drainage via overflow in supraglacial meltwater channels; and, by offline coupling with a second model, (d) the freezing and unfreezing of SGLs from autumn to spring. We call the model the Supraglacial Hydrology Evolution and Drainage (or SHED) model. We apply the model to three study regions in southwest Greenland between 2015 and 2019 (inclusive) and evaluate its performance with respect to observed supraglacial lake extents and proglacial discharge measurements. We show that the model reproduces 80 % of observed lake locations and provides good agreement with observations in terms of the temporal evolution of lake extent. Modelled moulin density values are in keeping with those previously published, and seasonal and inter-annual variability in proglacial discharge agrees well with that which is observed, though the observations lag the model by a few days since they include transit time through the subglacial system, while the model does not. Our simulations suggest that lake drainage behaviours may be more complex than traditional models suggest, with lakes in our model draining through a combination of both overflow and hydrofracture and with some lakes draining only partially and then refreezing. This suggests that, in order to simulate the evolution of Greenland's surface hydrological system with fidelity, a model that includes all of these processes needs to be used. In future work, we will couple our model to a subglacial model and an ice flow model and thus use our estimates of where, when, and how much meltwater gets to the bed to understand the consequences for ice flow.

Abstract. We compile, analyse and map all available geothermal heat flow measurements collected in and around Greenland into a new database of 419 sites and generate an accompanying spatial map. This database includes 290 sites previously reported by the International Heat Flow Commission (IHFC), for which we now standardize measurement and metadata quality. This database also includes 129 new sites, which have not been previously reported by the IHFC. These new sites consist of 88 offshore measurements and 41 onshore measurements, of which 24 are subglacial. We employ machine learning to synthesize these in situ measurements into a gridded geothermal heat flow model that is consistent across both continental and marine areas in and around Greenland. This model has a native horizontal resolution of 55 km. In comparison to five existing Greenland geothermal heat flow models, our model has the lowest mean geothermal heat flow for Greenland onshore areas (44 mW m–2). Our model’s most distinctive spatial feature is pronounced low geothermal heat flow (< 40 mW m–2) across the North Atlantic Craton of southern Greenland. Crucially, our model does not show an area of elevated heat flow that might be interpreted as remnant from the Icelandic Plume track. Finally, we discuss the substantial influence of paleoclimatic and other corrections on geothermal heat flow measurements in Greenland. The in-situ measurement database and gridded heat flow model, as well as other supporting materials, are freely available from the GEUS DataVerse (https://doi.org/10.22008/FK2/F9P03L; Colgan and Wansing, 2021).

The majority of the freshwater input from Greenland stems from the Greenland Ice Sheet. Despite its importance in terms of freshwater totals, there is a much higher number of individual catchments disconnected from the ice sheet contributing on average about 26% of the total Greenland freshwater flux. Most of those catchments have local glacier cover, only very few of them are instrumented and little scientific literature exists. We present a dataset of 12 years of discharge of four catchments less than 15 km apart, that are different in size (between 7 and 32 km²), local glacier coverage (4-11%) and lake cover (0-5%). They all drain into Kobbefjord, a well-studied fjord in West Greenland, near Greenland’s capital Nuuk. We find that annual specific discharge totals vary greatly (between 1.2 and 1.9 m/yr on a 12-year average within 15 km) due to a general climatic gradient and different strengths of orographic shading. The seasonal cycle differs among the sites mainly due to different exposure to solar radiation as a driver for major snowmelt; small ice coverage in the catchments plays only a minor role in discharge variability. Dry years generally increase the magnitude of spatial gradients in specific discharge and no significant temporal trends have been found in the studied catchments. On the sub-daily scale, the presence and elevation of lakes determines the catchment’s response during sunny days, leading to a difference in the timing of maximum discharge of between 7 and 12 hours depending on the site and the time of the year. The response of discharge to major precipitation events is discussed, where uniform reaction is found for the catchments with no lakes near the gauge and a delay of between 5 and 7 hours in the catchment with low-lying lakes. A comparison with a recently published modelled discharge time series on individual catchment scale shows the model’s capability of reproducing both snowmelt and large-scale storm events; however, the strong spatial heterogeneity of discharge magnitude and timing as well as the presence and variability of base-flow is not captured. We discuss methods to combine observational data with existing model output in order to improve the potential of their combined usage on the Greenland-scale.

The majority of freshwater input from Greenland to the global ocean stems from the Greenland Ice Sheet. Currently, almost a quarter of the freshwater flowing from Greenland is derived from catchments that are disconnected from the Greenland Ice Sheet. Despite their importance to the total freshwater flux and influence on fjord geochemistry, there is relatively little monitoring data available for those catchments and therefore the drivers of regional differences in export are largely unknown. We present a dataset of 12 years of discharge of four catchments less than 15 km apart, that are different in size (between 7 and 32 km(2)), local glacier coverage (4%-11%) and lake cover (0%-5%). They all drain into Kobbefjord, a well-studied fjord in West Greenland, near Greenland's capital Nuuk. Between catchments, the magnitude of discharge varies at annual, seasonal and sub-daily timescales, due to differences in physical catchment properties as well as local climate variability. We find that annual specific discharges vary greatly (between 1.2 and 1.9 m/year on a 12-year average) due to a longitudinal precipitation gradient from West to East caused by different amount of orographic precipitation shading. The seasonal cycle of discharge (amplitude, timing and minimum flow) differs among the sites mainly due to different exposure to solar radiation as a driver for major snowmelt; the small ice coverage in the catchments plays only a minor role in discharge variability. Dry years generally increase the relative differences in annual specific discharge and no significant temporal trends have been found in the studied catchments. On a sub-daily timescale, the difference in timing of maximum discharge during fair-weather days (>80% maximum solar radiation and no precipitation) ranged between 7 and 12 h, which is attributed to differences between the presence and elevation of lakes among the catchments. The response of discharge to major precipitation events is discussed, where a delay of between 5 and 7 h is found for the catchments that do not contain lakes near the gauge.

Meltwater runoff from the Greenland Ice Sheet changes water levels in glacial lakes and can lead to glacial lake outburst flooding (GLOF) events that threaten lives and property. Icebergs produced at Greenland’s marine terminating glaciers drift into Baffin Bay and the North Atlantic, where they can threaten shipping and offshore installations. Thus, monitoring glacial lake water levels and the drift of icebergs can enhance safety and aid in the scientific studies of glacial hydrology and iceberg-ocean interactions. The Maker Buoy was originally designed as a low-cost and open source sensor to monitor surface ocean currents. The open source framework, low-cost components, rugged construction and affordable satellite data transmission capabilities make it easy to customize for environmental monitoring in remote areas and under harsh conditions. Here, we present two such Maker Buoy variants that were developed to monitor water level in an ice-infested glacial lake in southern Greenland and to track drifting icebergs and moorings in the Vaigat Strait (Northwest Greenland). We describe the construction of each design variant, methods to access data in the field without an internet connection, and deployments in Greenland in summer 2019. The successful deployments of each Maker Buoy variant suggest that they may also be useful in operational iceberg management strategies and in GLOF monitoring programs.

Proglacial streams deliver melt water and chemical weathering products, including nutrients and radiogenic isotopes, from continental ice sheets to the ocean. Weathering products are also delivered to the ocean in non-glacial streams that form following ice sheet retreat and are disconnected from ice sheet meltwater by hydrologic divides. If weathering reactions differ in non-glacial and proglacial stream catchments, the streams could deliver different types and magnitudes of solutes to the ocean, depending on relative discharge volumes. Unlike proglacial streams, however, little is known of non-glacial stream solute compositions or discharge. Here we show specific discharges are similar from a proglacial stream draining the Greenland Ice Sheet (GrIS) with several streams disconnected from the ice sheet. We also evaluate weathering reactions across a 170-km transect in western Greenland that contains one proglacial stream draining the GrIS, and two coastal (ice distal) and three inland (ice proximal) areas with non-glacial streams. Non-glacial streams exhibit solute compositions and offsets between dissolved and bedload Sr isotope ratios that indicate weathering increases toward the coast with exposure age and precipitation. Major element mass balance calculations show weathering reactions shift from predominately carbonic acid weathering of carbonate minerals inland near the ice sheet to predominately sulfuric acid weathering of carbonate minerals near the coast. Strontium concentrations and isotopic ratios of the proglacial stream reflect mixing of at least two subglacial sources and minor in-stream weathering that consumes CO2. About 5 times less CO2 is consumed per liter in the proglacial than inland non-glacial streams; however, arid conditions inland suggest limited discharge from the ungauged inland streams leads to less total CO2 weathering than proglacial stream. One coastal area consumes less CO2 per liter than the proglacial stream and another coastal area exhibits net CO2 production. These results indicate estimates for glacial foreland solute fluxes and CO2 weathering consumption and production should include estimates from both non-glacial and proglacial streams. Understanding weathering fluxes from these two types of streams will be important for evaluations of past ice sheet retreat and predictions of future solute and CO2 fluxes associated with continued ice sheet retreat.

Ice-sheet melting is the primary water source for the proglacial Watson River in southern west Greenland. Discharge from the large, approximately 12,000 km2 ice-sheet catchment draining through the Watson River has been monitored since 2006. While this record is of respectable length for a Greenland monitoring effort, it is too short to resolve climate signals. Therefore, we use observed Tasersiaq lake discharge and Kangerlussuaq air temperature to reconstruct annual Watson River discharge back to 1949. The resulting sixty-five-year record shows that average ice-sheet runoff since 2003 has roughly increased by 46 percent relative to the 1949–2002 period. The time series suggests that the five top-ranking discharge years occurred since 2003. The three top-ranking discharge years (2010, 2012, and 2016) are characterized by melt seasons that were both long and intense. Interannual variability more than doubled since 2003, which we speculate to be because of hypsometric runoff amplification enhanced by albedo decrease and decreased firn permeability. The reconstructed time series proves to be a valuable tool for long-term evaluation of Greenland Ice Sheet surface mass balance models. A comparison with freshwater fluxes calculated by a downscaled version of the regional climate model RACMO2 reveals high correlation (r = 0.89), and also shows that the model possibly underestimates runoff by up to 26 percent in above-average melt years.

The recent decades of accelerating mass loss of the Greenland ice sheet have arisen from an increase in both surface meltwater runoff and ice flow discharge from tidewater glaciers. Despite the role of the Greenland ice sheet as the dominant individual cryospheric contributor to sea level rise in recent decades, no observational record of its mass loss spans the 30-year period needed to assess its climatological state. We present for the first time a 40-year (19752014) time series of observed meltwater discharge from a >6500-km(2) catchment of the southwestern Greenland ice sheet. We find that an abrupt 80% increase in runoff occurring between the 1976-2002 and 2003-2014 periods is due to a shift in atmospheric circulation, with meridional exchange events occurring more frequently over Greenland, establishing the first observation-based connection between ice sheet runoff and climate change.

Freshwater runoff to fjords with marine-terminating glaciers along the Greenland Ice Sheet margin has an impact on fjord circulation and potentially ice sheet mass balance through increasing heat transport to the glacier front. Here, the authors use the high-resolution (5.5 km) HIRHAM5 regional climate model, allowing high detail in topography and surface types, to estimate freshwater input to Godthabsfjord in southwest Greenland. Model output is compared to hydrometeorological observations and, while simulated daily variability in temperature and downwelling radiation shows high correlation with observations (typically >0.9), there are biases that impact the results. In particular, overestimated albedo leads to underestimation of melt and runoff at low elevations.In the model simulation (1991-2012), the ice sheet experiences increasing energy input from the surface turbulent heat flux (up to elevations of 2000 m) and shortwave radiation (at all elevations). Southerly wind anomalies and declining cloudiness due to an increase in atmospheric pressure over north Greenland contribute to increased summer melt. This results in declining surface mass balance (SMB), increasing surface runoff, and upward shift of the equilibrium line altitude.SMB is reconstructed back to 1890 though regression between simulated SMB and observed temperature and precipitation, with added uncertainty in the period 1890-1952 because of possible inhomogeneity in the precipitation record. SMB as low as in recent years appears to have occurred before, most notably around 1930, 1950, and 1960. While previous low SMBs were mainly caused by low accumulation, those around 1930 and in the 2000s are mainly due to warming.

The Greenland ice sheet is experiencing dramatic melt that is likely to continue with rapid Arctic warming. However, the proportion of meltwater stored before reaching the global ocean remains difficult to quantify. We use NASA MODIS surface reflectance data to estimate river discharge from two West Greenland rivers - the Watson River near Kangerlussuaq and the Naujat Kuat River near Nuuk - over the summers of 2000-12. By comparison with in situ river discharge observations, 'inundation-discharge' relations were constructed for both rivers. MODIS-based total annual discharges agree well with total discharge estimated from in situ observations (86% of summer discharge in 2009 to 96% in 2011 at the Watson River, and 106% of total discharge in 2011 to 104% in 2012 at the Naujat Kuat River). We find, however, that a time-lapse camera, deployed at the Watson River in summer 2012, better captures the variations in observed discharge, benefiting from fewer data gaps due to clouds. The MODIS-derived estimates indicate that summer discharge has not significantly increased over the last decade, despite a strong warming trend. Also, meltwater runoff estimates derived from the regional climate model RACMO2/GR for the drainage basins are higher than our reconstructions of river discharge. These results provide indirect evidence for a considerable component of water storage within the glacio-hydrological system.

The drainage of ice‐dammed lakes in the form of outburst floods in Greenland is detected regularly by remote sensing, and these events are expected to occur more frequently in a warmer climate. However, their impact on ice sheet stability and neighboring water bodies is still unknown. In this interdisciplinary study, we investigate lake drainages from the Greenland Ice Sheet into a west Greenland fjord by analyzing simultaneous time series of satellite observations and direct hydrographic measurements of temperature and salinity in the fjord. Satellite images show that, in general, lake drainages have occurred quasiperiodically during the last decade. A particular sequence of drainage events was observed by satellite in 2009 and was analyzed together with the first direct hydrographic observations. Signs of ice‐dammed lake drainages were observed by a downstream mooring located just below the intertidal zone. The release of freshwater occurred at the fjord subsurface at a tidewater outlet glacier. The downstream in‐water sequence of property changes in relation to these drainage events was observed as an almost immediate decrease in surface layer temperature (~2°C) followed within a week by the arrival of a high‐saline pulse (~ +5 units) with elevated salinity lasting for several days during the passage. During lake drainages, large amounts of relatively warm and saline intermediate water are brought to the near‐surface layers by entrainment processes near the glacier front, and this influences the hydrography of the fjord but also impacts the ecosystem through upwelling of nutrient‐rich intermediate water.

AbstractWe assess the runoff and surface mass balance (SMB) of the Greenland ice sheet in the Nuuk region (southwest) using output of two regional climate models (RCMs) evaluated by observations. The region encompasses six glaciers that drain into Godthåbsfjord. RCM data (1960–2012) are resampled to a high spatial resolution to include the narrow (relative to the native grid spacing) glacier trunks in the ice mask. Comparing RCM gridded results with automatic weather station (AWS) point measurements reveals that locally models can underestimate ablation and overestimate accumulation by up to tens of per cent. However, comparison with lake discharge indicates that modelled regional runoff totals are more accurate. Model results show that melt and runoff in the Nuuk region have doubled over the past two decades. Regional SMB attained negative values in recent high-melt years. Taking into account frontal ablation of the marine-terminating glaciers, the region lost 10–20 km3 w.e. a–1 in 2010–12. If 2010 melting prevails during the remainder of this century, a low-end estimate of sea-level rise of 5 mm is expected by 2100 from this relatively small section (2.6%) of the ice sheet alone.