Human alteration of global surface water storage variability.

Knowing the extent of human influence on the global hydrological cycle is essential for the sustainability of freshwater resources on Earth1,2. However, a lack of water level observations for the world's ponds, lakes and reservoirs has limited the quantification of human-managed (reservoir) changes in surface water storage compared to its natural variability3. The global storage variability in surface water bodies and the extent to which it is altered by humans therefore remain unknown. Here we show that 57 per cent of the Earth's seasonal surface water storage variability occurs in human-managed reservoirs. Using measurements from NASA's ICESat-2 satellite laser altimeter, which was launched in late 2018, we assemble an extensive global water level dataset that quantifies water level variability for 227,386 water bodies from October 2018 to July 2020. We find that seasonal variability in human-managed reservoirs averages 0.86 metres, whereas natural water bodies vary by only 0.22 metres. Natural variability in surface water storage is greatest in tropical basins, whereas human-managed variability is greatest in the Middle East, southern Africa and the western USA. Strong regional patterns are also found, with human influence driving 67 per cent of surface water storage variability south of 45 degrees north and nearly 100 per cent in certain arid and semi-arid regions. As economic development, population growth and climate change continue to pressure global water resources4, our approach provides a useful baseline from which ICESat-2 and future satellite missions will be able to track human modifications to the global hydrologic cycle.

High Voltage: The Molecular Properties of Redox-Active Dissolved Organic Matter in Northern High-Latitude Lakes.

Redox-active functional groups in dissolved organic matter (DOM) are crucial for microbial electron transfer and methane emissions. However, the extent of aquatic DOM redox properties across northern high-latitude lakes and their relationships with DOM composition have not been thoroughly described. We quantified electron donating capacity (EDC) and electron accepting capacity (EAC) in lake DOM from Canada to Alaska and assessed their relationships with parameters from absorbance, fluorescence, and ultrahigh resolution mass spectrometry (FT-ICR MS) analyses. EDC and EAC are strongly tied to aromaticity and negatively related to aliphaticity and protein-like content. Redox-active formulae spanned a range of aromaticity, including highly unsaturated phenolic formulae, and correlated negatively with many aliphatic N and S-containing formulae. This distribution illustrates the compositional diversity of redox-sensitive functional groups and their sensitivity to ecosystem properties such as local hydrology and residence time. Finally, we developed a reducing index (RI) to predict EDC in aquatic DOM from FT-ICR MS spectra and assessed its robustness using riverine DOM. As the hydrology of the northern high-latitudes continues to change, we expect differences in the quantity and partitioning of EDC and EAC within these lakes, which have implications for local water quality and methane emissions.

Direct measurements of meltwater runoff on the Greenland ice sheet surface.

Meltwater runoff from the Greenland ice sheet surface influences surface mass balance (SMB), ice dynamics, and global sea level rise, but is estimated with climate models and thus difficult to validate. We present a way to measure ice surface runoff directly, from hourly in situ supraglacial river discharge measurements and simultaneous high-resolution satellite/drone remote sensing of upstream fluvial catchment area. A first 72-h trial for a 63.1-km2 moulin-terminating internally drained catchment (IDC) on Greenland's midelevation (1,207-1,381 m above sea level) ablation zone is compared with melt and runoff simulations from HIRHAM5, MAR3.6, RACMO2.3, MERRA-2, and SEB climate/SMB models. Current models cannot reproduce peak discharges or timing of runoff entering moulins but are improved using synthetic unit hydrograph (SUH) theory. Retroactive SUH applications to two older field studies reproduce their findings, signifying that remotely sensed IDC area, shape, and supraglacial river length are useful for predicting delays in peak runoff delivery to moulins. Applying SUH to HIRHAM5, MAR3.6, and RACMO2.3 gridded melt products for 799 surrounding IDCs suggests their terminal moulins receive lower peak discharges, less diurnal variability, and asynchronous runoff timing relative to climate/SMB model output alone. Conversely, large IDCs produce high moulin discharges, even at high elevations where melt rates are low. During this particular field experiment, models overestimated runoff by +21 to +58%, linked to overestimated surface ablation and possible meltwater retention in bare, porous, low-density ice. Direct measurements of ice surface runoff will improve climate/SMB models, and incorporating remotely sensed IDCs will aid coupling of SMB with ice dynamics and subglacial systems.

Efficient meltwater drainage through supraglacial streams and rivers on the southwest Greenland ice sheet.

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(-1)) 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(3)⋅d(-1) vs. ∼0.103 km(3)⋅d(-1)), 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.

Toward global mapping of river discharge using satellite images and at-many-stations hydraulic geometry.

Rivers provide critical water supply for many human societies and ecosystems, yet global knowledge of their flow rates is poor. We show that useful estimates of absolute river discharge (in cubic meters per second) may be derived solely from satellite images, with no ground-based or a priori information whatsoever. The approach works owing to discovery of a characteristic scaling law uniquely fundamental to natural rivers, here termed a river's at-many-stations hydraulic geometry. A first demonstration using Landsat Thematic Mapper images over three rivers in the United States, Canada, and China yields absolute discharges agreeing to within 20-30% of traditional in situ gauging station measurements and good tracking of flow changes over time. Within such accuracies, the door appears open for quantifying river resources globally with repeat imaging, both retroactively and henceforth into the future, with strong implications for water resource management, food security, ecosystem studies, flood forecasting, and geopolitics.

New Trans-Arctic shipping routes navigable by midcentury.

Recent historic observed lows in Arctic sea ice extent, together with climate model projections of additional ice reductions in the future, have fueled speculations of potential new trans-Arctic shipping routes linking the Atlantic and Pacific Oceans. However, numerical studies of how projected geophysical changes in sea ice will realistically impact ship navigation are lacking. To address this deficiency, we analyze seven climate model projections of sea ice properties, assuming two different climate change scenarios [representative concentration pathways (RCPs) 4.5 and 8.5] and two vessel classes, to assess future changes in peak season (September) Arctic shipping potential. By midcentury, changing sea ice conditions enable expanded September navigability for common open-water ships crossing the Arctic along the Northern Sea Route over the Russian Federation, robust new routes for moderately ice-strengthened (Polar Class 6) ships over the North Pole, and new routes through the Northwest Passage for both vessel classes. Although numerous other nonclimatic factors also limit Arctic shipping potential, these findings have important economic, strategic, environmental, and governance implications for the region.

Development of Ice-Shelf Estuaries Promotes Fractures and Calving.

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.

Dark zone of the Greenland Ice Sheet controlled by distributed biologically-active impurities.

Albedo-a primary control on surface melt-varies considerably across the Greenland Ice Sheet yet the specific surface types that comprise its dark zone remain unquantified. Here we use UAV imagery to attribute seven distinct surface types to observed albedo along a 25 km transect dissecting the western, ablating sector of the ice sheet. Our results demonstrate that distributed surface impurities-an admixture of dust, black carbon and pigmented algae-explain 73% of the observed spatial variability in albedo and are responsible for the dark zone itself. Crevassing and supraglacial water also drive albedo reduction but due to their limited extent, explain just 12 and 15% of the observed variability respectively. Cryoconite, concentrated in large holes or fluvial deposits, is the darkest surface type but accounts for <1% of the area and has minimal impact. We propose that the ongoing emergence and dispersal of distributed impurities, amplified by enhanced ablation and biological activity, will drive future expansion of Greenland's dark zone.

Methane bubbling from northern lakes: present and future contributions to the global methane budget.

Large uncertainties in the budget of atmospheric methane (CH4) limit the accuracy of climate change projections. Here we describe and quantify an important source of CH4 -- point-source ebullition (bubbling) from northern lakes -- that has not been incorporated in previous regional or global methane budgets. Employing a method recently introduced to measure ebullition more accurately by taking into account its spatial patchiness in lakes, we estimate point-source ebullition for 16 lakes in Alaska and Siberia that represent several common northern lake types: glacial, alluvial floodplain, peatland and thermokarst (thaw) lakes. Extrapolation of measured fluxes from these 16 sites to all lakes north of 45 degrees N using circumpolar databases of lake and permafrost distributions suggests that northern lakes are a globally significant source of atmospheric CH4, emitting approximately 24.2+/-10.5Tg CH4yr(-1). Thermokarst lakes have particularly high emissions because they release CH4 produced from organic matter previously sequestered in permafrost. A carbon mass balance calculation of CH4 release from thermokarst lakes on the Siberian yedoma ice complex suggests that these lakes alone would emit as much as approximately 49000Tg CH4 if this ice complex was to thaw completely. Using a space-for-time substitution based on the current lake distributions in permafrost-dominated and permafrost-free terrains, we estimate that lake emissions would be reduced by approximately 12% in a more probable transitional permafrost scenario and by approximately 53% in a 'permafrost-free' Northern Hemisphere. Long-term decline in CH4 ebullition from lakes due to lake area loss and permafrost thaw would occur only after the large release of CH4 associated thermokarst lake development in the zone of continuous permafrost.

Rapid early development of circumarctic peatlands and atmospheric CH4 and CO2 variations.

An analysis of 1516 radiocarbon dates demonstrates that the development of the current circumarctic peatlands began approximately 16.5 thousand years ago (ka) and expanded explosively between 12 and 8 ka in concert with high summer insolation and increasing temperatures. Their rapid development contributed to the sustained peak in CH4 and modest decline of CO2 during the early Holocene and likely contributed to CH4 and CO2 fluctuations during earlier interglacial and interstadial transitions. Given the decreased tempo of peatland initiation in the late Holocene and the transition of many from fens (which generated high levels of CH4) to ombrotrophic bogs, a neoglacial expansion of northern peatlands cannot explain the increase in atmospheric CH4 that occurred after 6 ka.