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.

Remote sensing of river discharge (RSQ) is a burgeoning field rife with innovation. This innovation has resulted in a highly non-cohesive subfield of hydrology advancing at a rapid pace, and as a result misconceptions, mis-citations, and confusion are apparent among authors, readers, editors, and reviewers. While the intellectually diverse subfield of RSQ practitioners can parse this confusion, the broader hydrology community views RSQ as a monolith and such confusion can be damaging. RSQ has not been comprehensively summarized over the past decade, and we believe that a summary of the recent literature has a potential to provide clarity to practitioners and general hydrologists alike. Therefore, we here summarize a broad swath of the literature, and find after our reading that the most appropriate way to summarize this literature is first by application area (into methods appropriate for gauged, semi-gauged, regionally gauged, politically ungauged, and totally ungauged basins) and next by methodology. We do not find categorizing by sensor useful, and everything from un-crewed aerial vehicles (UAVs) to satellites are considered here. Perhaps the most cogent theme to emerge from our reading is the need for context. All RSQ is employed in the service of furthering hydrologic understanding, and we argue that nearly all RSQ is useful in this pursuit provided it is properly contextualized. We argue that if authors place each new work into the correct application context, much confusion can be avoided, and we suggest a framework for such context here. Specifically, we define which RSQ techniques are and are not appropriate for ungauged basins, and further define what it means to be ‘ungauged’ in the context of RSQ. We also include political and economic realities of RSQ, as the objective of the field is sometimes to provide data purposefully cloistered by specific political decisions. This framing can enable RSQ to respond to hydrology at large with confidence and cohesion even in the face of methodological and application diversity evident within the literature. Finally, we embrace the intellectual diversity of RSQ and suggest the field is best served by a continuation of methodological proliferation rather than by a move toward orthodoxy and standardization.

The Surface Water and Ocean Topography (SWOT) mission will vastly expand measurements of global rivers, providing critical new data sets for both gaged and ungaged basins. SWOT discharge products (available approximately 1 year after launch) will provide discharge for all river that reaches wider than 100 m. In this paper, we describe how SWOT discharge produced and archived by the US and French space agencies will be computed from measurements of river water surface elevation, width, and slope and ancillary data, along with expected discharge accuracy. We present for the first time a complete estimate of the SWOT discharge uncertainty budget, with separate terms for random (standard error) and systematic (bias) uncertainty components in river discharge time series. We expect that discharge uncertainty will be less than 30% for two‐thirds of global reaches and will be dominated by bias. Separate river discharge estimates will combine both SWOT and in situ data; these “gage‐constrained” discharge estimates can be expected to have lower systematic uncertainty. Temporal variations in river discharge time series will be dominated by random error and are expected to be estimated within 15% for nearly all reaches, allowing accurate inference of event flow dynamics globally, including in ungaged basins. We believe this level of accuracy lays the groundwork for SWOT to enable breakthroughs in global hydrologic science.

The forthcoming Surface Water and Ocean Topography (SWOT) NASA satellite mission will measure water surface width, height, and slope of major rivers worldwide. The resulting data could provide an unprecedented account of river discharge at continental scales, but reliable methods need to be identified prior to launch. Here we present a novel algorithm for discharge estimation from only remotely sensed stream width, slope, and height at multiple locations along a mass‐conserved river segment. The algorithm, termed the Bayesian AMHG‐Manning (BAM) algorithm, implements a Bayesian formulation of streamflow uncertainty using a combination of Manning's equation and at‐many‐stations hydraulic geometry (AMHG). Bayesian methods provide a statistically defensible approach to generating discharge estimates in a physically underconstrained system but rely on prior distributions that quantify the a priori uncertainty of unknown quantities including discharge and hydraulic equation parameters. These were obtained from literature‐reported values and from a USGS data set of acoustic Doppler current profiler (ADCP) measurements at USGS stream gauges. A data set of simulated widths, slopes, and heights from 19 rivers was used to evaluate the algorithms using a set of performance metrics. Results across the 19 rivers indicate an improvement in performance of BAM over previously tested methods and highlight a path forward in solving discharge estimation using solely satellite remote sensing.

The Surface Water and Ocean Topography (SWOT) satellite has the potential to transform global hydrologic science by offering simultaneous and synoptic estimates of river discharge and other hydraulic variables. Discharge is estimated from SWOT observations of water surface elevation, width, and slope. A first assessment using just the highest quality SWOT measurements, over the first 15 months (March 2023–July 2024) of the mission evaluated at 65 gauged reaches shows results consistent with pre‐launch expectations. SWOT estimates track discharge dynamics without relying on any gauge information: median correlation is 0.73, with a correlation interquartile range of 0.51–0.89. SWOT estimates capture discharge magnitude correctly in some cases but are biased (median bias is 50%) in others. There are already a total of 11,274 ungauged global locations with highest quality SWOT measurements where SWOT discharge is expected to accurately track discharge variations: this value will increase as SWOT data record length grows, algorithms are refined and SWOT measurements are reprocessed. This first look indicates that SWOT discharge is performing as expected for SWOT data that achieve performance requirements, providing observed information on discharge variations in ungauged basins globally.

Conventional satellite platforms are limited in their ability to monitor rivers at fine spatial and temporal scales: suffering from unavoidable trade‐offs between spatial and temporal resolutions. CubeSat constellations, however, can provide global data at high spatial and temporal resolutions, albeit with reduced spectral information. This study provides a first assessment of using CubeSat data for river discharge estimation in both gauged and ungauged settings. Discharge was estimated for 11 Arctic rivers with sizes ranging from 16 to >1,000 m wide using the Bayesian at‐many‐stations hydraulic geometry‐Manning algorithm (BAM). BAM‐at‐many‐stations hydraulic geometry solves for hydraulic geometry parameters to estimate flow and requires only river widths as input. Widths were retrieved from Landsat 8 and Sentinel‐2 data sets and a CubeSat (the Planet company) data set, as well as their fusions. Results show satellite data fusion improves discharge estimation for both large (>100 m wide) and medium (40–100 m wide) rivers by increasing the number of days with a discharge estimation by a factor of 2–6 without reducing accuracy. Narrow rivers (<40 m wide) are too small for Landsat and Sentinel‐2 data sets, and their discharge is also not well estimated using CubeSat data alone, likely because the four‐band sensor cannot resolve water surfaces accurately enough. BAM technique outperforms space‐based rating curves when gauge data are available, and its accuracy is acceptable when no gauge data are present (instead relying on global reanalysis for discharge priors). Ultimately, we conclude that the data fusion presented here is a viable approach toward improving discharge estimates in the Arctic, even in ungauged basins.