Greenland runoff, from ice mass loss and increasing rainfall, is increasing. That runoff, as discharge, impacts the physical, chemical, and biological properties of the adjacent fjords. However, where and when the discharge occurs is not readily available in an open database. Here we provide data sets of high-resolution Greenland hydrologic outlets, basins, and streams, as well as a daily 1958 through 2019 time series of Greenland liquid water discharge for each outlet. The data include 24 507 ice marginal outlets and upstream basins and 29 635 land coast outlets and upstream basins, derived from the 100 m ArcticDEM and 150 m BedMachine. At each outlet there are daily discharge data for 22 645 d – ice sheet runoff routed subglacially to ice margin outlets and land runoff routed to coast outlets – from two regional climate models (RCMs; MAR and RACMO). Our sensitivity study of how outlet location changes for every inland cell based on subglacial routing assumptions shows that most inland cells where runoff occurs are not highly sensitive to those routing assumptions, and outflow location does not move far. We compare RCM results with 10 gauges from streams with discharge rates spanning 4 orders of magnitude. Results show that for daily discharge at the individual basin scale the 5 % to 95 % prediction interval between modeled discharge and observations generally falls within plus or minus a factor of 5 (half an order of magnitude, or +500 %$/-\mathrm{80}$ %). Results from this study are available at https://doi.org/10.22008/promice/freshwater (Mankoff, 2020 a) and code is available at http://github.com/mankoff/freshwater (last access: 6 November 2020) (Mankoff, 2020 b).

Calving glaciers are rapidly retreating in many regions under the influence of ice‐water interactions at the glacier front. In contrast to the numerous researches conducted on fjords in front of tidewater glaciers, very few studies have been reported on lakes in which freshwater calving glaciers terminate. To better understand ice‐water interactions at the front of freshwater calving glaciers, we measured lakewater temperature, turbidity, and bathymetry near Glaciar Perito Moreno, Upsala, and Viedma, large calving glaciers of the Southern Patagonia Icefield. The thermal structures of these lakes were significantly different from those reported in glacial fjords. There was no indication of upwelling subglacial meltwater; instead, turbid and cold glacial water discharge filled the region near the lake bottom. This was because water density was controlled by suspended sediment concentrations rather than by water temperature. Near‐surface wind‐driven circulation reaches a depth of ~180 m, forming a relatively warm isothermal layer (mean temperature of ~5–6°C at Perito Moreno, ~3–4°C at Upsala, and ~6–7°C at Viedma), which should convey heat energy to the ice‐water interface. However, the deeper part of the glacier front is in contact with stratified cold water, implying a limited amount of melting there. In the lake in front of Glaciar Viedma, the region deeper than 120 m was filled entirely with turbid and very cold water at pressure melting temperature. Our results revealed a previously unexplored thermal structure of proglacial lakes in Patagonia, suggesting its importance in the subaqueous melting of freshwater calving glaciers.

Despite the importance of glacial lake development in ice dynamics and glacier thinning, in situ and satellite-based measurements from lake-terminating glaciers are sparse in the Bhutanese Himalaya, where a number of proglacial lakes exist. We acquired in situ and satellite-based observations across lake- and land-terminating debris-covered glaciers in the Lunana region, Bhutanese Himalaya. A repeated differential global positioning system survey reveals that thickness change of the debris-covered ablation area of the lake-terminating Lugge Glacier ($-\mathrm{4.67}\pm \mathrm{0.07}$ m a⁻¹) is more than 3 times more negative than that of the land-terminating Thorthormi Glacier ($-\mathrm{1.40}\pm \mathrm{0.07}$ m a⁻¹) for the 2004–2011 period. The surface flow velocities decrease down-glacier along Thorthormi Glacier, whereas they increase from the upper part of the ablation area to the terminus of Lugge Glacier. Numerical experiments using a two-dimensional ice flow model demonstrate that the rapid thinning of Lugge Glacier is driven by both a negative surface mass balance and dynamically induced ice thinning. However, the thinning of Thorthormi Glacier is minimised by a longitudinally compressive flow regime. Multiple supraglacial ponds on Thorthormi Glacier have been expanding since 2000 and have merged into a single proglacial lake, with the glacier terminus detaching from its terminal moraine in 2011. Numerical experiments suggest that the thinning of Thorthormi Glacier will accelerate with continued proglacial lake development.

Underwater ice geometry at the front of calving glaciers provides crucial information for calving and underwater melting. In this study, we present ice geometry captured by operating a side‐scanning sonar near the front of Glaciar Grey, a freshwater calving glacier in Patagonia. The observations revealed ice projecting into the lake with a substantially different structure from that of known tidewater glaciers. Terrace‐like ice structures were found at several tens of meters below the water surface and extended up to 100 m from the aerial ice front. The structure depicted by the sonar was confirmed when the ice front was exposed by flotation during a major calving event. We infer that buoyant force acting on the submerged ice terrace acted as a driver of the calving event. Our study demonstrates the importance of the underwater ice geometry, which affects sizable calving at the front of freshwater calving glaciers.