Thermally incised meltwater channels that flow each summer across melt-prone surfaces of the Greenland ice sheet have received little direct study. We use high-resolution WorldView-1/2 satellite mapping and in situ measurements to characterize supraglacial water storage, drainage pattern, and discharge across 6,812 km 2 of southwest Greenland in July 2012, after a record melt event. Efficient surface drainage was routed through 523 high-order stream/river channel networks, all of which terminated in moulins before reaching the ice edge. Low surface water storage (3.6 ± 0.9 cm), negligible impoundment by supraglacial lakes or topographic depressions, and high discharge to moulins (2.54-2.81 cm·d) indicate that the surface drainage system conveyed its own storage volume every <2 d to the bed. Moulin discharges mapped inside ∼52% of the source ice watershed for Isortoq, a major proglacial river, totaled ∼41-98% of observed proglacial discharge, highlighting the importance of supraglacial river drainage to true outflow from the ice edge. However, Isortoq discharges tended lower than runoff simulations from the Modèle Atmosphérique Régional (MAR) regional climate model (0.056-0.112 km ), and when integrated over the melt season, totaled just 37-75% of MAR, suggesting nontrivial subglacial water storage even in this melt-prone region of the ice sheet. We conclude that (i) the interior surface of the ice sheet can be efficiently drained under optimal conditions, (ii) that digital elevation models alone cannot fully describe supraglacial drainage and its connection to subglacial systems, and (iii) that predicting outflow from climate models alone, without recognition of subglacial processes, may overestimate true meltwater export from the ice sheet to the ocean.Greenland ice sheet | supraglacial hydrology | meltwater runoff | mass balance | remote sensing M eltwater runoff from the Greenland ice sheet (GrIS) accounts for half or more of its total mass loss to the global ocean (1, 2) but remains one of the least-studied hydrologic processes on Earth. Each summer, a complex system of supraglacial meltwater ponds, lakes, streams, rivers, and moulins develops across large areas of the southwestern GrIS surface, especially below ∼1,300 m elevation (3-7), with supraglacial erosion driven by thermal and radiative processes (5). Digital elevation models (DEMs) suggest a poorly drained surface resulting from abundant topographic depressions, which computational flow routing models must artificially "fill" to allow hydrological flow paths extending from the ice sheet interior to its edge (8-11). The realism of such modeled flow paths remains largely untested by real-world observations. To date, most observational studies of GrIS supraglacial hydrology have focused on large lakes (∼1 km 2

No abstract

We document the existence of widespread firn aquifers in an elevation range of ~1200–2000 m, in the high snow‐accumulation regions of the Greenland ice sheet. We use NASA Operation IceBridge accumulation radar data from five campaigns (2010–2014) to estimate a firn‐aquifer total extent of 21,900 km2. We investigate two locations in Southeast Greenland, where repeated radar profiles allow mapping of aquifer‐extent and water table variations. In the upper part of Helheim Glacier the water table rises in spring following above‐average summer melt, showing the direct firn‐aquifer response to surface meltwater production changes. After spring 2012, a drainage of the firn‐aquifer lower margin (5 km) is inferred from both 750 MHz accumulation radar and 195 MHz multicoherent radar depth sounder data. For 2011–2014, we use a ground‐penetrating radar profile located at our Ridgeline field site and find a spatially stable aquifer with a water table fluctuating less than 2.5 m vertically. When combining radar data with surface topography, we find that the upper elevation edge of firn aquifers is located directly downstream of locally high surface slopes. Using a steady state 2‐D groundwater flow model, water is simulated to flow laterally in an unconfined aquifer, topographically driven by ice sheet surface undulations until the water encounters crevasses. Simulations suggest that local flow cells form within the Helheim aquifer, allowing water to discharge in the firn at the steep‐to‐flat transitions of surface topography. Supported by visible imagery, we infer that water drains into crevasses, but its volume and rate remain unconstrained.

In summer 2003 central Europe suffered an unusually severe heat wave, with air temperatures similar to those predicted for an average summer during the late 21st century. We use a unique set of over half a century of lake data from two lakes in Switzerland to determine the effect of the 2003 heat wave on water temperature and oxygen conditions in order to assess how temperate lakes will react when exposed to the increased ambient summer air temperatures that will be encountered in a generally warmer world and to test the predictions of relevant simulation models. In both lakes, surface temperature and thermal stability in summer 2003 were the highest ever recorded, exceeding the long-term mean by more than 2.5 standard deviations. The extremely high degree of thermal stability resulted in extraordinarily strong hypolimnetic oxygen depletion. These results are consistent with the predictions of the simulation models. Additionally, the results indicate that climatic warming will increase the risk of occurrence of deep-water anoxia, thus counteracting long-term efforts that have been undertaken to ameliorate the effects of anthropogenic eutrophication.

Past estimates of Greenland Ice Sheet accumulation rates have been multiyear climatologies based on ice/firn cores and coastal precipitation records. Existing annually resolved estimates have incompletely quantified uncertainty, owing primarily to incomplete spatial coverage. This study improves upon these shortcomings by calibrating annual (1958–2007) solid precipitation output from the Fifth Generation Mesoscale Model modified for polar climates (Polar MM5) using firn core and meteorological station data. The calibration employs spatial interpolation of regionally derived linear correction functions. Residual uncertainties exhibit coherent spatial patterns, which are modeled via spatial interpolation of root mean squared errors. Mean 1958–2007 Greenland Ice Sheet annual accumulation rate is 337 ± 48 mm/yr water equivalent (w.e.) or 591 ± 83 Gt/yr. Annual estimates contain one standard deviation uncertainties of 74 mm/yr w.e., 22%, or 129 Gt/yr. Accumulation rates in southeast Greenland are found to exceed 2000 mm/yr w.e. and to dominate interannual variability in Greenland Ice Sheet total accumulated mass, representing 31% of the whole. Accumulation rates in the southeast are of sufficient magnitude to affect the sign of Greenland mass balance during some years. The only statistically significant temporal change in total ice sheet accumulation in the 1958–2007 period occurred between 1960 and 1972, when a simultaneous accumulation increase and decrease occurred in west and east Greenland, respectively. No statistically significant uniform change in ice sheet‐wide accumulation is evident after 1972. However, regional changes do occur, including an accumulation increase on the west coast post‐1992. The high accumulation rates of 2002–2003 appear to be confined to the southeast.