Search

River networks represent the largest biogeochemical nexus between the continents, ocean and atmosphere. Our current understanding of the role of rivers in the global carbon cycle remains limited, which makes it difficult to predict how global change may alter the timing and spatial distribution of riverine carbon sequestration and greenhouse gas emissions. Here we review the state of river ecosystem metabolism research and synthesize the current best available estimates of river ecosystem metabolism. We quantify the organic and inorganic carbon flux from land to global rivers and show that their net ecosystem production and carbon dioxide emissions shift the organic to inorganic carbon balance en route from land to the coastal ocean. Furthermore, we discuss how global change may affect river ecosystem metabolism and related carbon fluxes and identify research directions that can help to develop better predictions of the effects of global change on riverine ecosystem processes. We argue that a global river observing system will play a key role in understanding river networks and their future evolution in the context of the global carbon budget.

Humans have drastically disrupted the global sediment cycle. Suspended sediment flux and concentration are key controls over both river morphology and river ecosystems. Our ability to understand sediment dynamics within river corridors is limited by observations. Here, we present RivSed, a database of satellite observations of suspended sediment concentration from 1984-2018 across 460 large (>60 m wide) US rivers that provides a new, spatially explicit, view of river sediment. We found that 32% of US rivers have a declining temporal trend in sediment concentration, with a mean reduction of 40% since 1984, whereas only 2% have an increasing trend. Most rivers (52%) show decreasing sediment concentration longitudinally moving downstream, typically due to a few large dams rather than the accumulated effect of many small dams. Comparing our observations with modeled “pre-dam” longitudinal suspended sediment concentrations, most rivers (53%) show different patterns. However, contemporary longitudinal patterns in concentration are remarkably stable from year to year since 1984, with more stability in large, highly managed rivers with less cropland. RivSed has broad applications for river geomorphology and ecology and highlights anthropogenic effects on river corridors across the US.

Rivers are among the most diverse, dynamic, and productive ecosystems on Earth. River flow regimes are constantly changing, but characterizing and understanding such changes have been challenging from a long-term and global perspective. By analyzing water extent variations observed from four-decade Landsat imagery, we here provide a global attribution of the recent changes in river regime to morphological dynamics (e.g., channel shifting and anabranching), expansion induced by new dams, and hydrological signals of widening and narrowing. Morphological dynamics prevailed in ~20% of the global river area. Booming reservoir constructions, mostly skewed in Asia and South America, contributed to ~32% of the river widening. The remaining hydrological signals were characterized by contrasting hotspots, including prominent river widening in alpine and pan-Arctic regions and narrowing in the arid/semi-arid continental interiors, driven by varying trends in climate forcing, cryospheric response to warming, and human water management. Our findings suggest that the recent river extent dynamics diverge based on hydroclimate and socio-economic conditions, and besides reflecting ongoing morphodynamical processes, river extent changes show close connections with external forcings, including climate change and anthropogenic interference.

Ice plays key roles in regulating hydrological, ecological, biogeochemical, and socioeconomic functions of lakes. Long-term in situ lake ice phenological records indicate that lake ice is trending toward later freeze-up, earlier breakup, and a shorter ice duration. Parallel to study of lake ice using in situ records and process-based models, satellite remote sensing can expand our understanding of lake ice change over large spatial scales. However, most remote sensing studies have focused on large lakes or short periods of time, which may not robustly represent changes over multidecadal time periods or in the much more numerous small lakes. Here, we present a random forest model, Sensitive Lake Ice Detection (SLIDE), to accurately extract ice conditions from Landsat TM, ETM+, and OLI images. We trained the model using a manually labeled lake ice dataset (1089 labeled areas over 995 lakes globally). Our results show that our model achieves accurate classification between ice/snow and water (accuracy: 97.8%, kappa coefficient: 95.5%). Comparing Landsat-derived ice cover with in situ ice conditions, we show that our model produces less bias, lower RMSE, and higher kappa than does the Landsat snow/ice flag from the quality assessment band. This is especially true during the transitional period surrounding the ice on and off dates reported from in situ (mean bias -7.3% from our model, -17.3% from the Landsat quality band). Our results demonstrate the feasibility of mining the rich Landsat archive to study lake ice dynamics and of better flagging ice-affected lake observations.

Consistent estimation of water surface area from remote sensing remains challenging in regions such as South Asia with vegetation, mountainous topography, and persistent monsoonal cloud cover. High-resolution optical imagery, which is often used for global inundation mapping, is highly impacted by clouds, while synthetic aperture radar (SAR) imagery is not impacted by clouds and is affected by both topographic layover and vegetation. Here, we compare and contrast inundation extent measurements from visible (Landsat-8 and Sentinel-2) and SAR (Sentinel-1) imagery. Each data type (wavelength) has complementary strengths and weaknesses which were gauged separately over selected water bodies in Bangladesh. High-resolution cloud-free PlanetScope imagery at 3-m resolution was used as a reference to check the accuracy of each technique and data type. Next, the optical and radar images were fused for a rule-based water area classification algorithm to derive the optimal decision for the water mask. Results indicate that the fusion approach can improve the overall accuracy by up to 3.8%, 18.2%, and 8.3% during the wet season over using the individual products of Landsat8, Sentinel-1, and Sentinel-2, respectively, at three sites, while providing increased observational frequency. The fusion-derived products resulted in overall accuracy ranging from 85.8% to 98.7% and Kappa coefficient varying from 0.61 to 0.83. The proposed SAR-visible fusion technique has potential for improving satellite-based surface water monitoring and storage changes, especially for smaller water bodies in humid tropical climate of South Asia.

Global demand for vegetable oil and high oil palm yield have driven the rapid expansion of oil palm plantations in tropical countries. The research literature widely recognizes the effects of forest cover conversion into oil palm on biodiversity, deforestation, and carbon dynamics. However, research on the ecohydrological impacts of oil palm plantations is sparse, despite studies indicating that oil palm development may reshape land-water interactions and the availability and movement of water at different spatial and temporal scales. We address this gap by conducting a systematic literature review on oil palm development and its relation to ecohydrological processes. We found 139 relevant papers up to the year 2021, addressing different ecohydrological processes related to oil palm. We reviewed their spatiotemporal scales, geographic distribution, oil palm species and age, and the effects of land conversion from forest, cropland, and pastures. We also incorporated societal aspects regarding community perceptions of water. Our review highlights the effects of oil palm plantations on three main components of the water cycle: (i) land-atmosphere, (ii) fluvial systems, and (iii) soils and groundwater. Most studies include analyses of the Indo-Malayan and Australasian biogeographical regions (113), followed by the Neotropics (49) and the Afrotropics (15). Compared to rainforests, oil palm monocultures are warmer and drier. They have higher evapotranspiration (ET) rates, lower runoff regulation and infiltration capacity, and lower soil organic carbon (SOC). Although less often implemented, alternative oil palm management practices, including oil palm agroforestry, can help to mitigate some of these effects. Forest to oil palm conversion is the most studied land transition, while conversions from croplands, pastures, and grasslands are less studied. Overall, we identify gaps in understanding the long-term effects of management on ecohydrological processes under different land conversions, especially in the Neotropics and the Afrotropics, precluding research-informed policy to manage impacts of this expanding crop.

Author(s): Durand, Michael; Barros, Ana; Dozier, Jeff; Adler, Robert; Cooley, Sarah; Entekhabi, Dara; Forman, Barton A; Konings, Alexandra G; Kustas, William P; Lundquist, Jessica D; Pavelsky, Tamlin M; Rodell, Matthew; Steele‐Dunne, Susan

Changes in a river's width reflect natural and anthropogenic impacts on local and upstream/downstream hydraulic and hydrologic processes. Temporal variation of river width also impacts biogeochemical exchange and reflects geomorphologic evolution. However, while global maps of mean river width and dynamic water surface extent exist, there is currently no standardized global assessment of river widths that documents changes over time. Therefore, we made repeated width measurements from Landsat images for all rivers wider than 90 m collected from 1984 to 2020 (named Global LOng-term river Width, GLOW), which consists of ∼1.2 billion cross-sectional river width measurements, with an average of 3,000 width measurements per 10-km reach. With GLOW, we investigated the temporal variations of global river width, quantified by the interquartile range (IQR) and temporal trend. We found that 85% of global rivers have a width IQR 150 m. We also found that 37% of global river segments show significant temporal trends in width over the past 37 years, and this number is higher (46%) for human-regulated rivers. Further, we leveraged machine learning to identify the most important factors explaining river width variations and revealed that these driving factors are significantly different between free-flowing and non-free-flowing rivers. Specifically, the most important factor driving temporal variations in river width is the climate for free-flowing rivers, and is soil condition for human-impacted rivers. Finally, we anticipate that this study and the public release of GLOW will improve the understanding of river dynamics and catalyze additional interdisciplinary studies.

The wetted width of a river is one of the most important hydraulic parameters that can be readily measured using remote sensing. Remotely sensed river widths are used to estimate key attributes of river systems, including changes in their surface area, channel storage, and discharge. Although several published algorithms automate river network and width extraction from remote sensing images, they are limited by only being able to run on local computers and do not automatically manage cloudy images as input. Here we present RivWidthCloud, a river width software package developed on the Google Earth Engine cloud computing platform. RivWidthCloud automatically extracts river centerline and widths from optical satellite images with the ability to flag observations that are obstructed by features like clouds, cloud shadows, and snow based on existing quality band classification. Because RivWidthCloud is built on a popular cloud computing platform, it allows users to easily apply the algorithm to the platform’s vast archive of remote sensing images, thereby reducing the users’ overhead for computing hardware and data storage. By comparing RivWidthCloud-derived widths from Landsat images to in situ widths from the U.S. and Canada, we show that RivWidthCloud can estimate widths with high accuracy (root mean square error: 99 m; mean absolute error: 43 m; mean bias: −21 m). By making RivWidthCloud publicly available, we anticipate that it will be used to address both river science questions and operational applications of water resource management.

More than one-third of Earth’s landmass is drained by rivers that seasonally freeze over. Ice transforms the hydrologic1,2, ecologic3,4, climatic5 and socio-economic6–8 functions of river corridors. Although river ice extent has been shown to be declining in many regions of the world1, the seasonality, historical change and predicted future changes in river ice extent and duration have not yet been quantified globally. Previous studies of river ice, which suggested that declines in extent and duration could be attributed to warming temperatures9,10, were based on data from sparse locations. Furthermore, existing projections of future ice extent are based solely on the location of the 0-°C isotherm11. Here, using satellite observations, we show that the global extent of river ice is declining, and we project a mean decrease in seasonal ice duration of 6.10 ± 0.08 days per 1-°C increase in global mean surface air temperature. We tracked the extent of river ice using over 400,000 clear-sky Landsat images spanning 1984–2018 and observed a mean decline of 2.5 percentage points globally in the past three decades. To project future changes in river ice extent, we developed an observationally calibrated and validated model, based on temperature and season, which reduced the mean bias by 87 per cent compared with the 0-degree-Celsius isotherm approach. We applied this model to future climate projections for 2080–2100: compared with 2009–2029, the average river ice duration declines by 16.7 days under Representative Concentration Pathway (RCP) 8.5, whereas under RCP 4.5 it declines on average by 7.3 days. Our results show that, globally, river ice is measurably declining and will continue to decline linearly with projected increases in surface air temperature towards the end of this century. An analysis based on Landsat imagery shows that the extent of river ice has declined extensively over past decades and that this trend will continue under future global warming.

Satellite estimates of inland water quality have the potential to vastly expand our ability to observe and monitor the dynamics of large water bodies. For almost 50 years, we have been able to remotely sense key water quality constituents like total suspended sediment, dissolved organic carbon, chlorophyll a, and Secchi disk depth. Nonetheless, remote sensing of water quality is poorly integrated into inland water sciences, in part due to a lack of publicly available training data and a perception that remote estimates are unreliable. Remote sensing models of water quality can be improved by training and validation on larger data sets of coincident field and satellite observations, here called matchups. To facilitate model development and deeper integration of remote sensing into inland water science, we have built AquaSat, the largest such matchup data set ever assembled. AquaSat contains more than 600,000 matchups, covering 1984–2019, of ground-based total suspended sediment, dissolved organic carbon, chlorophyll a, and SDDSecchi disk depth measurements paired with spectral reflectance from Landsat 5, 7, and 8 collected within ±1 day of each other. To build AquaSat, we developed open source tools in R and Python and applied them to existing public data sets covering the contiguous United States, including the Water Quality Portal, LAGOS-NE, and the Landsat archive. In addition to publishing the data set, we are also publishing our full code architecture to facilitate expanding and improving AquaSat. We anticipate that this work will help make remote sensing of inland water accessible to more hydrologists, ecologists, and limnologists while facilitating novel data-driven approaches to monitoring and understanding critical water resources at large spatiotemporal scales.

Spatiotemporally continuous global river discharge estimates across the full spectrum of stream orders are vital to a range of hydrologic applications, yet they remain poorly constrained. Here we present a carefully designed modeling effort (Variable Infiltration Capacity land surface model and Routing Application for Parallel computatIon of Discharge river routing model) to estimate global river discharge at very high resolutions. The precipitation forcing is from a recently published 0.1° global product that optimally merged gauge‐, reanalysis‐, and satellite‐based data. To constrain runoff simulations, we use a set of machine learning‐derived, global runoff characteristics maps (i.e., runoff at various exceedance probability percentiles) for grid‐by‐grid model calibration and bias correction. To support spaceborne discharge studies, the river flowlines are defined at their true geometry and location as much as possible—approximately 2.94 million vector flowlines (median length 6.8 km) and unit catchments are derived from a high‐accuracy global digital elevation model at 3‐arcsec resolution (~90 m), which serves as the underlying hydrography for river routing. Our 35‐year daily and monthly model simulations are evaluated against over 14,000 gauges globally. Among them, 35% (64%) have a percentage bias within ±20% (±50%), and 29% (62%) have a monthly Kling‐Gupta Efficiency ≥0.6 (0.2), showing data robustness at the scale the model is assessed. This reconstructed discharge record can be used as a priori information for the Surface Water and Ocean Topography satellite mission's discharge product, thus named “Global Reach‐level A priori Discharge Estimates for Surface Water and Ocean Topography”. It can also be used in other hydrologic applications requiring spatially explicit estimates of global river flows., Key Points Thirty‐five‐year global reconstruction of river flows at unprecedented 2.94 million reaches extracted from 90‐m high‐accuracy DEMNovel bias correction (“sparse CDF matching”) is conducted grid‐by‐grid against machine learning‐derived global runoff characteristics mapsEvaluation against >14,000 gauges shows data robustness for hydrologic applications and utility in spaceborne discharge estimation efforts

The Surface Water and Ocean Topography (SWOT) satellite mission aims to improve the frequency and accuracy of global observations of river water surface elevations (WSEs) and slopes. As part of the SWOT mission, an airborne analog, AirSWOT, provides spatially-distributed measurements of WSEs for river reaches tens to hundreds of kilometers in length. For the first time, we demonstrate the ability of AirSWOT to consistently measure temporal dynamics in river WSE and slope. We evaluate data from six AirSWOT flights conducted between June 7–22, 2015 along a ~90 km reach of the Tanana River, AK. To validate AirSWOT measurements, we compare AirSWOT WSEs and slopes against an in situ network of 12 pressure transducers (PTs). Assuming error-free in situ data, AirSWOT measurements of river WSEs have an overall root mean square difference (RMSD) of 11.8 cm when averaged over 1 km2 areas while measurements of river surface slope have an RMSD of 1.6 cm/km for reach lengths >5 km. AirSWOT is also capable of recording accurate river WSE changes between flight dates, with an RMSD of 9.8 cm. Regrettably, observed in situ slope changes that transpired between the six flights are well below AirSWOT's accuracy, limiting the evaluation of AirSWOT's ability to capture temporal changes in slope. In addition to validating the direct AirSWOT measurements, we compare discharge values calculated via Manning's equation using AirSWOT WSEs and slopes to discharge values calculated using PT WSEs and slopes. We define or calibrate the remaining discharge parameters using a combination of in situ and remotely sensed observations, and we hold these remaining parameters constant between the two types of calculations to evaluate the impact of using AirSWOT versus the PT observations of WSE and slope. Results indicate that AirSWOT-derived discharge estimates are similar to the PT-derived discharge estimates, with an RMSD of 13.8%. Additionally, 42% of the AirSWOT-based discharge estimates fall within the PT discharge estimates' uncertainty bounds. We conclude that AirSWOT can measure multitemporal variations in river WSE and spatial variations in slope with both high accuracy and spatial sampling, providing a compelling alternative to in situ measurements of regional-scale, spatiotemporal fluvial dynamics.

Using river centerlines created with Landsat images and the Shuttle Radar Topography Mission digital elevation model, we created spatially continuous maps of mean annual flow river width, slope, meander wavelength, sinuosity, and catchment area for all rivers wider than 90 m located between 60°N and 56°S. We analyzed the distributions of these properties, identified their typical ranges, and explored relationships between river planform and slope. We found width to be directly associated with the magnitude of meander wavelength and catchment area. Moreover, we found that narrower rivers show a larger range of slope and sinuosity values than wider rivers. Finally, by comparing simulated discharge from a water balance model with measured widths, we show that power laws between mean annual discharge and width can predict width typically to −35% to +81%, even when a single relationship is applied across all rivers with discharge ranging from 100 to 50,000 m 3 /s.

The origins of the Surface Water and Ocean Topography (SWOT) mission date back to the mid-1970s with the launch of GOES-3 and SEASAT. These missions were then followed in 1992 by the Topex-Poseidon...

Groundwater modules are critically important to the simulation of low flows in physically based land surface models (LSMs) and conceptual rainfall-runoff models (HBV). Here, we develop a Groundwater for Ungauged Basins (GrUB) module that uses only physically based properties for which data are widely available, thus allowing its application without the need for calibration. GrUB is designed to be computationally simple and readily adaptable to a wide variety of LSMs and rainfall-runoff models. We assess the performance of GrUB in 84 United States watersheds by incorporating it into HBV, a popular rainfall-runoff model. We compare predictions of low flows by the native (calibrated) HBV groundwater module with those by the (uncalibrated) GrUB module and find that GrUB generates error metrics that are equivalent or superior to those generated by the (calibrated) HBV groundwater module. To assess whether predictions by GrUB are robust to changes in the structure and parameterization of the overlying hydrologic model, we run tests for two artificial scenarios: Slow Recharge with rates of percolation below 0.1 mm/day, and Fast Recharge with rates of percolation of up to 1,000 mm/day. GrUB proves to be robust to these extreme changes, with mean absolute error (MAE) of predictions of low flows only increasing by an average of up to 17%, while average MAE increases by up to 158% when the same tests are performed on HBV without the GrUB module. We suggest GrUB as a potential tool for improving predictions of low flows in LSMs, as well as rainfall-runoff models when calibration data are sparse.

In land surface models (LSMs), the hydraulic properties of the subsurface are commonly estimated according to the texture of soils at the Earth's surface. This approach ignores macropores, fracture flow, heterogeneity, and the effects of variable distribution of water in the subsurface on effective watershed‐scale hydraulic variables. Using hydrograph recession analysis, we empirically constrain estimates of watershed‐scale effective hydraulic conductivities (K) and effective drainable aquifer storages (S) of all reference watersheds in the conterminous United States for which sufficient streamflow data are available (n = 1,561). Then, we use machine learning methods to model these properties across the entire conterminous United States. Model validation results in high confidence for estimates of log(K) (r2 > 0.89; 1% 0.83; −70%

The volume of water stored in seasonal wetlands is a fundamental but difficult to measure variable for developing a physical understanding of wetland behavior. For seasonal wetlands that are a major source of water for rice and fish production, this physical understanding is key to planning for water-food security and ecosystem services. This study quantified variations in volumetric storage for the numerous seasonal wetlands of northeastern Bangladesh, locally known as “haors.” These haors receive transboundary runoff from densely vegetated and mountainous terrain in India and face persistent monsoonal cloud cover as they become full. We estimated volumetric storage for 13 haors by using extensive remote sensing data on water surface extent and elevation that was complemented with citizen-contributed gauge data. Assuming a trapezoidal bathymetry, an area-volume relationship was developed for selected haors. This relationship was assumed to be valid for extrapolating volumetric estimations over all the haors in the region. Results suggested that as haors get filled with the onset of monsoon rains, total estimated storage relative to the lowest observable level varied from 6.5 (±0.4) km3 in May to 30.9 (±2.0) km3 in July (peak of monsoon). Choosing a rectangular bathymetry can lead to 47% higher estimates compared to trapezoidal cross section. Estimating this intra-annual/interannual increase in storage is important for the region to plan water management policies that balance the human and ecosystem needs. Our analytical approach has potential for application to wetlands worldwide in light of the upcoming Surface Water and Ocean Topography (SWOT) mission.

Surface water resources are vital for biological, physical, and socioeconomic processes on the Earth’s surface. To advance monitoring capabilities of these resources, new remote sensing technologies including the forthcoming Surface Water and Ocean Topography (SWOT) satellite (expected launch 2022) and its experimental airborne prototype AirSWOT are being developed to repeatedly map water surface elevation (WSE) and slope (WSS) of the world’s rivers, lakes, and reservoirs. However, the vertical accuracies of these novel technologies are largely unverified thus standard and repeatable field procedures to validate remotely sensed WSE and WSS are needed. To that end, we designed, engineered and operationalized a Water Surface Profiler (WaSP) system that efficiently and accurately surveys WSE and WSS in a variety of surface water environments using Global Navigation Satellite Systems (GNSS) time-averaged measurements with Precise Point Positioning (PPP) corrections. Here, we present WaSP construction, deployment, and a data processing workflow. We demonstrate WaSP data collections from repeat field deployments in the North Saskatchewan River and three prairie pothole lakes near Saskatoon, Saskatchewan, Canada. We find that WaSP reproducibly measures WSE and WSS with vertical accuracies similar to standard field survey methods (WSE RMSD ~8 cm, WSS RMSD ~1.3 cm/km), and that repeat WaSP deployments accurately quantify water level changes (RMSD ~3 cm). Collectively, these results suggest that WaSP is an easily-deployed, self-contained system with sufficient accuracy for validating the decimeter level expected accuracies of SWOT and AirSWOT. We conclude by discussing the utility of WaSP for validating airborne and spaceborne WSE mappings, present 63 WaSP in-situ lake WSE measurements collected in support of NASA’s Arctic Boreal and Vulnerability Experiment (ABoVE), highlight routine deployment in support of the Lake Observation by Citizen Scientists and Satellites (LOCSS) project, and explore WaSP utility for validating a novel GNSS interferometric reflectometry (GNSS-IR) LArge Wave Warning System (LAWWS).

Lakes are often defined by seasonal cycles. The seasonal timing, or phenology, of many lake processes are changing in response to human activities. However, long‐term records exist for few lakes, and extrapolating patterns observed in these lakes to entire landscapes is exceedingly difficult using the limited number of available in situ observations. Limited landscape‐level observations mean we do not know how common shifts in lake phenology are at macroscales. Here, we use a new remote sensing data set, LimnoSat‐US, to analyze U.S. summer lake color phenology between 1984 and 2020 across more than 26,000 lakes. Our results show that summer lake color seasonality can be generalized into five distinct phenology groups that follow well‐known patterns of phytoplankton succession. The frequency with which lakes transition from one phenology group to another is tied to lake and landscape level characteristics. Lakes with high inflows and low variation in their seasonal surface area are generally more stable, while lakes in areas with high interannual variations in climate and catchment population density show less stability. Our results reveal previously unexamined spatiotemporal patterns in lake seasonality and demonstrate the utility of LimnoSat‐US, which, with over 22 million remote sensing observations of lakes, creates novel opportunities to examine changing lake ecosystems at a national scale., Key Points Summer lake color phenology can be generalized into five distinct seasonal patterns of greening and blueing eventsSince the mid‐1990s, the number of lakes with color patterns corresponding to eutrophic waterbodies has been increasingWe observe these patterns using a new U.S. lake remote sensing data set that contains over 22 million lake observations

River channels store large volumes of water globally, critically impacting ecological and biogeochemical processes. Despite the importance of river channel storage, there is not yet an observationa...

Remote sensing of hydrology provides the key to one of the longest-standing problems in hydrology: ungaged basins. The global network of river discharge measurements available publicly accurately captures overall streamflow into global ocean, but coverage for smaller tributaries upstream, is vanishingly small. The forthcoming Surface Water and Ocean Topography (SWOT) mission will vastly expand measurements of global rivers. However, converting the SWOT measurements of height, width and slope into river discharge is non-trivial. SWOT discharge algorithms will provide discharge for all SWOT-observed rivers, but at lower accuracy than what is possible at a gage. In this paper, we answer the question of where and how SWOT discharge will be valuable. We show that median river discharge accuracy will be on the order of 50%, in terms of normalized root mean squared error (nRMSE), whereas unbiased nRMSE will be on the order of 10-20%. We discuss how SWOT discharge will contribute to global hydrologic science in the coming decade.

One important goal of the Surface Water and Ocean Topography (SWOT) satellite mission is to estimate global river discharges from observations of river width, height, and slope. While a range of algorithms have been developed and intercompared for SWOT (e.g., Durand et al. 2016), our understanding of the algorithm accuracy has been confined to tens of rivers globally, due to the limited SWOT-like observations currently collected from hydraulic model outputs.To scale up the assessment of SWOT discharge algorithms, this study will first collect discharge observations at thousands of global gauges whose river widths are wider than 50 m (i.e., observable by SWOT), to provide the most comprehensive observations to evaluate discharge estimations. Then at those gauges, all available Landsat images from 2010 to 2017 (8 years) will be collected to extract river widths with an automatic Google Earth Engine tool called RivWidthCloud (Yang et al. 2019). The extracted river width time series (temporally intermittent) will provide SWOT-like observations, which can be used to derive discharge using the Bayesian AMHG-Manning (BAM) algorithm (Hagemann et al. 2017; Feng et al. 2019). The prior discharge information needed by the BAM algorithm will come from an updated global discharge modeling database (Lin et al. 2019). These datasets, collectively, will provide a critical assessment of the SWOT discharge algorithms.This study is expected to provide the first geographically explicit assessment of the BAM algorithm at thousands of global locations, and the insights gained may also help the global hydrologic modeling community with their data assimilation efforts.