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

SignificanceMeltwater runoff is an important hydrological process operating on the Greenland ice sheet surface that is rarely studied directly. By combining satellite and drone remote sensing with continuous field measurements of discharge in a large supraglacial river, we obtained 72 h of runoff observations suitable for comparison with climate model predictions. The field observations quantify how a large, fluvial supraglacial catchment attenuates the magnitude and timing of runoff delivered to its terminal moulin and hence the bed. The data are used to calibrate classical fluvial hydrology equations to improve meltwater runoff models and to demonstrate that broad-scale surface water drainage patterns that form on the ice surface powerfully alter the timing, magnitude, and locations of meltwater penetrating into the ice sheet.

Millions of dollars are spent annually on revitalizing salmon spawning in riverbeds where redd building by female salmon is inhibited by sediment that is too big for fish to move. Yet the conditions necessary for productive spawning remain unclear. There is no gauge for quantifying how grain size influences the reproductive potential of coarse-bedded rivers. Hence, managers lack a quantitative basis for optimizing spawning habitat restoration for reproductive value. To overcome this limitation, we studied spawning by Chinook, sockeye, and pink salmon (Oncorhynchus tshawytscha, O. nerka, and O. gorbuscha) in creeks and rivers of California and the Pacific Northwest. Our analysis shows that coarse substrates have been substantially undervalued as spawning habitat in previous work. We present a field-calibrated approach for estimating the number of redds and eggs a substrate can accommodate from measurements of grain size and fish length. Bigger fish can move larger sediment and thus use more riverbed area for spawning. They also tend to have higher fecundity, and so can deposit more eggs per redd. However, because redd area increases with fish length, the number of eggs a substrate can accommodate is maximized for moderate-sized fish. This previously unrecognized tradeoff raises the possibility that differences in grain size help regulate river-to-river differences in salmon size. Thus, population diversity and species resilience may be linked to lithologic, geomorphic, and climatic factors that determine grain size in rivers. Our approach provides a tool for managing grain-size distributions in support of optimal reproductive potential and species resilience.

[1] Understanding river form and behavior requires an efficient means of measuring channel morphology. This study evaluated the potential to map the bathymetry of two clear-flowing, shallow (<3 m deep) gravel bed rivers <60 m wide from 2 m-pixel WorldView2 satellite images. Direct measurements of water column optical properties were used to quantify constraints on depth retrieval. The smallest detectable change in depth was 0.01-0.04 m and the maximum detectable depth was 5 m in green bands but <2 m in the near-infrared; lower sensor radiometric resolution yields less precise estimates over a smaller range. An algorithm for calibrating a band ratio X to field measurements of depth d proved effective when applied to spectra extracted from images (R 2 = 0.822 and 0.594 for the larger and smaller stream, respectively) or measured in the field (R 2 = 0.769 and 0.452). This procedure also identified optimal wavelength combinations, but different bands were selected for each site. Accuracy assessment of bathymetric maps produced using various calibration approaches and image types indicated that: 1) a linear d vs. X relation provided depth estimates nearly as accurate as a quadratic formulation; 2) panchromatic and pan-sharpened multispectral images with smaller 0.5 m pixels did not yield more reliable depth estimates than the original images; and 3) depth retrieval was less reliable in pools due to saturation of the radiance signal. This investigation thus demonstrated the feasibility, as well as the limitations, of measuring the bathymetry of clear, shallow gravel bed rivers from space.

The forthcoming Surface Water and Ocean Topography (SWOT) satellite mission will provide global measurements of the free surface of large rivers, providing new opportunities for remote sensing‐derived estimates of river discharge in gaged and ungaged basins. SWOT discharge algorithms have been developed and benchmarked using synthetic data but remain untested on real‐world swath altimetry observations. We present the first discharge estimates from AirSWOT, a SWOT‐like airborne Ka‐band radar, using 6 days of measurements over a 40‐km segment of the Willamette River in Oregon, USA. The three evaluated discharge algorithms estimated discharge with normalized root‐mean‐square errors of 10–31% when compared with in situ gage data but were sensitive to an initial estimate of mean annual discharge. Our results show that these discharge algorithms provide reliable discharge estimates on remotely sensed data at SWOT‐like spatial scales while highlighting the need for further algorithm sensitivity tests.