The contribution of Greenland Ice Sheet meltwater runoff to global sea-level rise is accelerating due to increased melting of its bare-ice ablation zone. There is growing evidence, however, that climate models overestimate runoff from this critical area of the ice sheet. Climate models traditionally assume that all bare-ice runoff enters the ocean, unlike porous firn, in which some meltwater is retained and/or refrozen. We used field measurements and numerical modeling to reveal that extensive retention and refreezing also occurs in bare glacier ice. We found that, from 2009 to 2018, meltwater refreezing in bare, porous glacier ice reduced runoff by an estimated 11-17 Gt a-1 in southwest Greenland alone, equivalent to 9-15% of this sector's annual meltwater runoff simulated by climate models. This mass retention explains evidence from prior studies of runoff overestimation on bare ice by current generation climate models and may represent an overlooked buffer on projected runoff increases. Inclusion of bare-ice retention and refreezing processes in climate models therefore has immediate potential to improve forecasts of ice sheet runoff and its contribution to sea-level rise.

The forthcoming Surface Water and Ocean Topography (SWOT) satellite and AirSWOT airborne instrument are the first imaging radar-altimeters designed with near-nadir, 35.75 GHz Ka-band InSAR for mapping terrestrial water storage variability. Remotely sensed surface water extents are crucial for assessing such variability, but are confounded by emergent and inundated vegetation along shorelines. However, because SWOT-like measurements are novel, there remains some uncertainty in the ability to detect certain land and water classes. We study the likelihood of misclassification between 15 land cover types and develop the Ka-band Phenomenology Scattering (KaPS) scattering model to simulate changes to radar backscatter as a result of changing surface water fraction and roughness. Using a separability metric, we find that water is five times more distinct compared with dry land classes, but has the potential to be confused with littoral zone and wet soil cover types. The KaPS scattering model simulates AirSWOT backscatter for incidence angles 1-27°, identifying the conditions under which open water is likely to be confused with littoral zone and wet soil cover types. A comparison of KaPS simulated backscatter with AirSWOT observed backscatter shows good overall agreement across the 15 classes (median r2=0.76). KaPS characterization of the sensitivity of near-nadir, Ka-band SAR to small changes in both wet area fraction and surface roughness enables more nuanced classification of inundation area. These results provide additional confidence in the ability of SWOT to classify water inundation extent, and open the door for novel hydrological and ecological applications of future Ka-band SAR missions.

Meltwater runoff from the Greenland ice sheet (GrIS) is an important contributor to global sea level rise, but substantial uncertainty exists in its measurement and prediction. Common approaches for estimating ice sheet runoff are in situ gauging of proglacial rivers draining the ice sheet and surface mass balance (SMB) modeling. To obtain hydrological and meteorological data sets suitable for both runoff stage characterization and, pending the establishment of stage–discharge curves, SMB model evaluation, we established an automated weather station (AWS) and a cluster of traditional and experimental river stage sensors on the Minturn River, the largest proglacial river draining Inglefield Land, NW Greenland. Secondary installations measuring river stage were installed in the Fox Canyon River and North River at Pituffik Space Base, NW Greenland. Proglacial runoff at these sites is dominated by supraglacial processes only, uniquely advantaging them for SMB studies. The three installations provide rare hydrological time series and an opportunity to evaluate experimental measurements of river stage from a harsh, little-studied polar region. The installed instruments include submerged vented and non-vented pressure transducers, a bubbler sensor, experimental bank-mounted laser rangefinders, and time-lapse cameras. The first 3 years of observations (2019 to 2021) from these stations indicate (a) a meltwater runoff season from late June to late August/early September that is roughly synchronous throughout the region; (b) the early onset (∼ 23 June to 8 July) of a strong diurnal runoff signal in 2019 and 2020, suggesting minimal meltwater storage in snow and/or firn; (c) 1 d lagged air temperature that displays the strongest correlation with river stage; (d) river stage that correlates more strongly with ablation zone albedo than with net radiation; and (e) the late-summer rain-on-ice events appear to trigger the region's sharpest and largest floods. The new gauging stations provide valuable in situ hydrological observations that are freely available through the PROMICE network (https://promice.org/weather-stations/, last access: 14 September 2023).

Abstract. Meltwater runoff from the Greenland Ice Sheet (GrIS) is an important contributor to global sea level rise, but substantial uncertainty exists in its measurement and prediction. Common approaches for estimating ice sheet runoff are in situ gauging of proglacial rivers draining the ice sheet, and surface mass balance (SMB) modeling. To obtain hydrological and meteorological datasets suitable for both runoff characterization and SMB model validation, we established an automated weather station (AWS) and cluster of traditional and experimental river stage sensors on the Minturn River, the largest proglacial river draining Inglefield Land, NW Greenland. Secondary installations measuring river stage were installed in the Fox Canyon River and North River at Thule Air Base, NW Greenland. Proglacial runoff at these sites is dominated by supraglacial processes only, uniquely advantaging them for SMB studies. The three installations provide rare hydrological time-series and an opportunity to evaluate experimental measurements of river stage from a harsh, little-studied polar region. The installed instruments include submerged vented and non-vented pressure transducers, a bubbler sensor, experimental bank-mounted laser rangefinders, and time-lapse cameras. The first three years of observations (2019 to 2021) from these stations indicate a) a meltwater runoff season from late June to late August/early September, roughly synchronous throughout the region; b) early onset (~ June 23 to July 8) of a strong diurnal runoff signal in 2019 and 2020, suggesting minimal meltwater storage in snow/firn; c) one-day lagged air temperature displays the strongest correlation with river stage; d) river stage correlates more strongly with ablation zone albedo than with net radiation; and e) late-summer rain-on-ice events appear to trigger the region’s sharpest and largest floods. The new gauging stations provide valuable in situ hydrological observations from a little-studied, rapidly changing area and are freely available through the PROMICE network (https://promice.org/weather-stations/).

The Greenland Ice Sheet is a leading source of global sea level rise, due to surface meltwater runoff and glacier calving. However, given a scarcity of proglacial river gauge measurements, ice sheet runoff remains poorly quantified. This lack of in situ observations is particularly acute in Northwest Greenland, a remote area releasing significant runoff and where traditional river gauging is exceptionally challenging. Here, we demonstrate that georectified time-lapse camera images accurately retrieve stage fluctuations of the proglacial Minturn River, Inglefield Land, over a 3 year study period. Camera images discern the river’s wetted shoreline position, and a terrestrial LiDAR scanner (TLS) scan of riverbank microtopography enables georectification of these positions to vertical estimates of river stage. This non-contact approach captures seasonal, diurnal, and episodic runoff draining a large (∼2,800 km 2 ) lobe of grounded ice at Inglefield Land with good accuracy relative to traditional in situ bubble-gauge measurements ( r 2 = 0.81, Root Mean Square Error (RMSE) ±0.185 m for image collection at 3-h frequency; r 2 = 0.92, RMSE ±0.109 m for resampled average daily frequency). Furthermore, camera images effectively supplement other instrument data gaps during icy and/or low flow conditions, which challenge bubble-gauges and other contact-based instruments. This benefit alone extends the effective seasonal hydrological monitoring period by ∼2–4 weeks each year for the Minturn River. We conclude that low-cost, non-contact time-lapse camera methods offer good promise for monitoring proglacial meltwater runoff from the Greenland Ice Sheet and other harsh polar environments.

Mass loss from the Greenland Ice Sheet (GrIS) is a primary contributor to sea level rise, but substantial uncertainty exists in estimates of future ice sheet losses. Surface mass balance (SMB) models, the current leading approach to sea level rise projection, anticipate continued dominance of runoff as a mass loss pathway. Despite their preeminence, SMB models in vulnerable northern environments lack adequate field validation, particularly for error-sensitive runoff estimates. We have installed a cluster of high quality field instruments at the Minturn Elv, a proglacial river site in Inglefield Land, NW Greenland to provide discharge and weather datasets for the validation and refinement of climate/SMB runoff models. The instrument cluster has meteorological, hydrological, and time lapse camera instrumentation, including a vented water level stage recorder, single shot and scanning lidars, time lapse cameras, and in situ ADCP discharge and terrestrial scanning lidar measurements. The instrument suite provides novel flow and weather datasets with the opportunity to evaluate experimental approaches to stage measurement in adverse, high-latitude areas. Inglefield is a uniquely advantaged location because proglacial runoff is dominated by SMB processes operating on the ice surface without interference from subglacial hydrology. Overall, our hydrometeorological instrument cluster at Inglefield Land will provide one of the few validation datasets for regional climate models outside of Southwest Greenland.

Areas of lakes that support emergent aquatic vegetation emit disproportionately more methane than open water but are under‐represented in upscaled estimates of lake greenhouse gas emissions. These shallow areas are typically less than ∼1.5 m deep and can be detected with synthetic aperture radar (SAR). To assess the importance of lake emergent vegetation (LEV) zones to landscape‐scale methane emissions, we combine airborne SAR mapping with field measurements of vegetated and open‐water methane flux. First, we use Uninhabited Aerial Vehicle SAR data from the NASA Arctic‐Boreal Vulnerability Experiment to map LEV in 4,572 lakes across four Arctic‐boreal study areas and find it comprises ∼16% of lake area, exceeding previous estimates, and exhibiting strong regional differences (averaging 59 [50–68]%, 22 [20–25]%, 1.0 [0.8–1.2]%, and 7.0 [5.0–12]% of lake areas in the Peace‐Athabasca Delta, Yukon Flats, and northern and southern Canadian Shield, respectively). Next, we account for these vegetated areas through a simple upscaling exercise using paired methane fluxes from regions of open water and LEV. After excluding vegetated areas that could be accounted for as wetlands, we find that inclusion of LEV increases overall lake emissions by 21 [18–25]% relative to estimates that do not differentiate lake zones. While LEV zones are proportionately greater in small lakes, this relationship is weak and varies regionally, underscoring the need for methane‐relevant remote sensing measurements of lake zones and a consistent criterion for distinguishing wetlands. Finally, Arctic‐boreal lake methane upscaling estimates can be improved with more measurements from all lake zones.

Surface melting impacts ice sheet sliding by supplying water to the bed, but subglacial processes driving ice accelerations are complex. We examine linkages between surface runoff, transient subglacial water storage, and short-term ice motion from 168 consecutive hourly measurements of meltwater discharge (moulin input) and GPS-derived ice surface motion for Rio Behar, a moulin-terminating supraglacial river catchment on the southwest Greenland Ice Sheet. Short-term accelerations in ice speed correlate strongly with lag-corrected measures of supraglacial river discharge ( r = 0.9, τ = 0.7, p < 0.01). Though our 7 days record cannot address seasonal-scale forcing, diurnal ice accelerations align with normalized differenced supraglacial and proglacial discharge, a proxy for subglacial storage change, better than GPS-derived ice surface uplift. These basal in absolute drive accelerations. present of in situ in a large supraglacial river draining the ice sheet surface, just of it plummets into a GPS of ice surface motion record brief accelerations in ice sliding speed that follow daily cycles in meltwater entering the moulin. By comparing these measurements with proglacial river discharges leaving the ice sheet, we identify daily fluctuations in subglacial water storage that track short-term accelerations in ice motion. These findings affirm the importance of supraglacial rivers to subglacial water pressure and ice dynamics, even in relatively thick ice >40 km inland from the ice terminus.

Glacial environments offer the opportunity to study the incipient stages of chemical weathering due to the high availability of finely ground sediments, low water temperatures, and typically short rock-water interaction times. In this study we focused on the geochemical behavior of germanium (Ge) in west Greenland, both during subglacial weathering by investigating glacier-fed streams, as well as during a batch reactor experiment by allowing water-sediment interaction for up to 2 years in the laboratory. Sampled in late August 2014, glacial stream Ge and Si concentrations were low, ranging between 12–55 pmol/L and 7–33 µmol/L, respectively (Ge/Si = 0.9–2.2 µmol/mol, similar to parent rock). As reported previously, the dissolved stable Ge isotope ratio (δ 74 Ge) of the Watson River was 0.86 ± 0.24‰, the lowest among global rivers and streams measured to date. This value was only slightly heavier than the suspended load (0.48 ± 0.23‰), which is likely representative of the bulk parent rock composition. Despite limited Ge/Si and δ74Ge Ge fractionation, both Ge and Si appear depleted relative to Na during subglacial weathering, which we interpret as the relatively congruent uptake of both phases by amorphous silica (aSi). Continued sediment-water interaction over 470–785 days in the lab produced a large increase in dissolved Si concentrations (up to 130–230 µmol/L), a much smaller increase in dissolved Ge (up to ∼70 pmol/L), resulting in a Ge/Si decrease (to 0.4–0.5 µmol/mol) and a significant increase in δ 74 Ge (to 1.9–2.2‰). We argue that during the experiment, both Si and Ge are released by the dissolution of previously subglacially formed aSi, and Ge is then incorporated into secondary phases (likely adsorbed to Fe oxyhydroxides), with an associated Δ 74 Ge secondary−dissolved fractionation factor of −2.15 ± 0.46‰. In summary, we directly demonstrate Ge isotope fractionation during the dissolution-precipitation weathering reactions of natural sediments in the absence of biological Ge and Si uptake, and highlight the significant differences in Ge behavior during subglacial and non-glacial weathering.

Recent work has identified complex perennial supraglacial stream and river networks in areas of the Greenland Ice Sheet (GrIS) ablation zone. Current surface mass balance (SMB) models appear to overestimate meltwater runoff in these networks compared to in-channel measurements of supraglacial discharge. Here, we constrain SMB models using the hillslope river routing model (HRR), a spatially explicit flow routing model used in terrestrial hydrology, in a 63 km2 supraglacial river catchment in southwest Greenland. HRR conserves water mass and momentum and explicitly accounts for hillslope routing (i.e., flow over ice and/or firn on the GrIS), and we produce hourly flows for nearly 10 000 channels given inputs of an ice surface digital elevation model (DEM), a remotely sensed supraglacial channel network, SMB-modeled runoff, and an in situ discharge dataset used for calibration. Model calibration yields a Nash–Sutcliffe efficiency as high as 0.92 and physically realistic parameters. We confirm earlier assertions that SMB runoff exceeds the conserved mass of water measured in this catchment (by 12 %–59 %) and that large channels do not dewater overnight despite a diurnal shutdown of SMB runoff production. We further test hillslope routing and network density controls on channel discharge and conclude that explicitly including hillslope flow and routing runoff through a realistic fine-channel network (as opposed to excluding hillslope flow and using a coarse-channel network) produces the most accurate results. Modeling complex surface water processes is thus both possible and necessary to accurately simulate the timing and magnitude of supraglacial channel flows, and we highlight a need for additional in situ discharge datasets to better calibrate and apply this method elsewhere on the ice sheet.

As the global climate warms, increased surface meltwater production on ice shelves may trigger ice-shelf collapse and enhance global sea-level rise. The formation of surface rivers could help prevent ice-shelf collapse if they can efficiently evacuate meltwater. Here, we present observations of the evolution of a surface river into an ice-shelf estuary atop the Petermann Ice Shelf in northwest Greenland, and identify a second estuary at the nearby Ryder Ice Shelf. This surface hydrology process can foster fracturing and enhance calving. At the Petermann estuary, sea ice was observed converging at the river mouth upstream, indicating a flow reversal. Seawater persists in the estuary, after the surrounding icescape is frozen. Along the base of Petermann estuary, linear fractures were initiated at the calving front and propagated upstream along the channel. Similar fractures along estuary channels shaped past large rectilinear calving events at the Petermann and Ryder Ice Shelves. Increased surface melting in a warming world will enhance fluvial incision, promoting estuary development, longitudinal fracturing orthogonal to ice-shelf fronts, and increase rectilinear calving. Estuaries could develop in Antarctica within the next half-century, resulting in increased calving and accelerating both ice loss and global sea-level rise.

Abstract The Greenland Ice Sheet’s contribution to global sea-level rise is accelerating1 due to increased melting of its bare-ice ablation zone2–6, but there is growing evidence that climate models overestimate runoff from this critical area of the ice sheet7–12. Current climate models assume all bare ice runoff escapes to the ocean, unlike snow covered areas where some fraction of runoff is retained and/or refrozen in porous firn13–15. Here we use in situ measurements and numerical modeling to reveal extensive retention and refreezing of liquid meltwater in bare glacial ice, explaining chronic runoff overestimation by climate models. From 2009–2018, refreezing of liquid meltwater in bare, porous glacial ice reduced meltwater runoff by 11–23 Gt a-1 in southwest Greenland alone, equivalent to 10–20% of annual meltwater production. This mass retention is commensurate with current estimates of climate model ice sheet meltwater runoff uncertainty, and may represent an overlooked buffer on projected runoff increases for the coming century16. Inclusion of bare-ice retention and refreezing processes in climate models therefore has immediate potential to improve forecasts of ice sheet runoff and its contribution to global sea-level rise.

Areas of lakes that support emergent aquatic vegetation emit disproportionately more methane than open water but are underrepresented in upscaled estimates of lake greenhouse gas emissions. These shallow areas are typically less than ̃1.5 m deep and can be estimated through synthetic aperture radar (SAR) mapping. To assess the importance of lake emergent vegetation (LEV) zones to landscape-scale methane emissions, we combine airborne SAR mapping with field measurements of vegetated and open-water methane flux. First, we use Uninhabited Aerial Vehicle SAR (UAVSAR) data from the NASA Arctic-Boreal Vulnerability Experiment (ABoVE) to map LEV in 4,572 lakes across four Arctic-boreal study areas and find it comprises ̃16% of lake area, exceeding previous estimates, and exhibiting strong regional differences (averaging 59 [50–68]%, 22 [20-25]%, 1.0 [0.8-1.2]%, and 7.0 [5.0-12]% of lake areas in the Peace-Athabasca Delta, Yukon Flats, and northern and southern Canadian Shield, respectively). Next, we account for these vegetated areas through a simple upscaling exercise using paired methane fluxes from regions of open water and LEV. After excluding vegetated areas that could be accounted for as wetlands, we find that inclusion of LEV increases overall lake emissions by 21 [18-25]% relative to estimates that do not differentiate lake zones. While LEV zones are proportionately greater in small lakes, this relationship is weak and varies regionally, underscoring the need for methane-relevant remote sensing measurements of lake zones and a consistent criterion for distinguishing wetlands. Finally, Arctic-boreal lake methane upscaling estimates can be improved with more measurements from all lake zones.

Supraglacial lakes and rivers dominate the storage and transport of meltwater on the southwest Greenland Ice Sheet (GrIS) surface. Despite functioning as interconnected hydrologic networks, supraglacial lakes and rivers are commonly studied as independent features, resulting in an incomplete understanding of their collective impact on meltwater storage and routing. We use Landsat 8 satellite imagery to assess the seasonal evolution of supraglacial lakes and rivers on the southwest GrIS during the 2015 melt season. Remotely sensed meltwater areas and volumes are compared with surface runoff simulations from three climate models (MERRA-2, MAR 3.6 and RACMO 2.3), and with in situ observations of proglacial discharge in the Watson River. We find: (1) at elevations >1600 m, 21% of supraglacial lakes and 28% of supraglacial rivers drain into moulins, signifying the presence of high-elevation surface-to-bed meltwater connections even during a colder-than-average melt season; (2) while supraglacial lakes dominate instantaneous surface meltwater storage, supraglacial rivers dominate total surface meltwater area and discharge; (3) the combined surface area of supraglacial lakes and rivers is strongly correlated with modeled surface runoff; and (4) of the three models examined here, MERRA-2 runoff yields the highest overall correlation with observed proglacial discharge in the Watson River.

In situ river discharge estimation is a critical component of studying rivers. A dominant method for establishing discharge monitoring in situ is a temporary gauge, which uses a rating curve to relate stage to discharge. However, this approach is constrained by cost and the time to develop the stage‐discharge rating curve, as rating curves rely on numerous flow measurements at high and low stages. Here, we offer a novel alternative approach to traditional temporary gauges: estimating Discharge via Arrays of Pressure Transducers (DAPT). DAPT uses a Bayesian discharge algorithm developed for the upcoming Surface Water Ocean Topography satellite (SWOT) to estimate in situ discharge from automated water surface elevation measurements. We conducted sensitivity tests over 4,954 model runs on five gauged rivers and conclude that the DAPT method can robustly reproduce discharge with an average Nash‐Sutcliffe Efficiency (NSE) of 0.79 and Kling‐Gupta Efficiency of 0.78. Further, we find that the DAPT method estimates discharge similarly to an idealized temporary gauge created from the same input data (NSE differences of less than 0.1), and that results improve significantly with accurate priors. Finally, we test the DAPT method in nine poorly gauged rivers in a realistic and complex field setting in the Peace‐Athabasca Delta, and show that the DAPT method largely outperforms a temporary gauge in this time and budget constrained setting. We therefore recommend DAPT as an effective tool for in situ discharge estimation in cases where there is not enough time or resources to develop a temporary gauge.

This dataset contains georeferenced three-band orthomosaics of green, red, and near-infrared (NIR) digital imagery at 1m resolution collected over selected surface waters across Alaska and Canada between July 9 and August 17, 2017. The orthomosaics were generated from individual images collected by a Cirrus Designs Digital Camera System (DCS) mounted on a Beechcraft Super King Air B200 aircraft from approximately 8-11 km altitude. Flights were over the following areas: Saskatchewan River, Saskatoon, Inuvik, Yukon River including Yukon Flats, Sagavanirktok River, Arctic Coastal Plain, Old Crow Flats, Peace-Athabasca Delta, Slave River, Athabasca River, Yellowknife, Great Slave Lake, Mackenzie River and Delta, Daring Lake, and other selected locations. Most locations were imaged twice during two flight campaigns in Canada and Alaska extending roughly SE-NW then NW-SE up to a month apart. The data were georeferenced using 303 ground control points (GCPs) across the study region.

Meltwater runoff from the Greenland Ice Sheet (GrIS) significantly contributes to sea level rise and is the dominant driver of enhanced mass loss. While most melt occurs during summer, little is known about its seasonal and/or interannual retention within the GrIS. Here, we document evidence of runoff during winter, ~4 months after summer melt. Ground‐penetrating radar and borehole surveys in the proglacial Isortoq River reveal slowly flowing water beneath >0.5 m of river ice. Geochemical analysis of this water indicates previous contact with the ice sheet bed. Comparable surveys in proglacial rivers draining four neighboring catchments found no winter drainage, despite a brief surface melt event ~10 days prior. We attribute the observed runoff to residual meltwater storage and release enabled by a >600 m deep trough beneath Isunguata Sermia, but not neighboring glaciers. We conclude that the GrIS bed can stay wet and drain small amounts of meltwater year‐round.

AirSWOT is an experimental airborne Ka-band radar interferometer developed by NASA-JPL as a validation instrument for the forthcoming NASA Surface Water and Ocean Topography (SWOT) satellite mission. In 2017, AirSWOT was deployed as part of the NASA Arctic Boreal Vulnerability Experiment (ABoVE) to map surface water elevations across Alaska and western Canada. The result is the most extensive known collection of near-nadir airborne Ka-band interferometric synthetic aperture radar (InSAR) data and derivative high-resolution (3.6 m pixel) digital elevation models to produce water surface elevation (WSE) maps. This research provides a synoptic assessment of the 2017 AirSWOT ABoVE dataset to quantify regional WSE errors relative to coincidentin situfield surveys and LiDAR data acquired from the NASA Land, Vegetation, and Ice Sensor (LVIS) airborne platform. Results show that AirSWOT WSE data can penetrate cloud cover and have nearly twice the swath-width of LVIS as flown for ABoVE (3.2 km vs. 1.8 km nominal swath-width). Despite noise and biases, spatially averaged AirSWOT WSEs can be used to estimate sub-seasonal hydrologic variability, as confirmed with field GPS surveys andin situpressure transducers. This analysis informs AirSWOT ABoVE data users of known sources of measurement error in the WSEs as influenced by radar parameters including incidence angle, magnitude, coherence, and elevation uncertainty. The analysis also provides recommended best practices for extracting information from the dataset by using filters for these four parameters. Improvements to data handing would significantly increase the accuracy and spatial coverage of future AirSWOT WSE data collections, aiding scientific surface water studies, and improving the platform's capability as an airborne validation instrument for SWOT.

Accurate, transparent knowledge of global reservoir levels is a prerequisite for effective management of water resources. However, no complete database exists because gauge data are not globally available and the current generation of satellite radar altimeters resolves only the world's largest reservoirs. Here, we investigate water level changes in global reservoirs using ICESat‐2, National Aeronautics and Space Administration (NASA)'s new satellite laser altimetry mission. In just the first 12 months of the mission, we find that ICESat‐2 accurately (±14.1 cm) retrieved water level changes for 3,712 global reservoirs having surface areas ranging from <1 to >10,000 km2. From this new global data set, we identify distinct regional patterns in reservoir level change that can be attributed to both water availability and management strategy. Our findings demonstrate that ICESat‐2 will form a crucial component of any global reservoir level inventory and enable new insight into how reservoir management responds to climatic variability and increasing human demand.