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.

Significance: Stream/river carbon dioxide (CO ₂ ) emission has significant spatial and seasonal variations critical for understanding its macroecosystem controls and plumbing of the terrestrial carbon budget. We relied on direct fluvial CO ₂ partial pressure measurements and seasonally varying gas transfer velocity and river network surface area estimates to resolve reach-level seasonal variations of the flux at the global scale. The percentage of terrestrial primary production (GPP) shunted into rivers that ultimately contributes to CO ₂ evasion increases with discharge across regions, due to a stronger response in fluvial CO ₂ evasion to discharge than GPP. This highlights the importance of hydrology, in particular water throughput, in terrestrial–fluvial carbon transfers and the need to account for this effect in plumbing the terrestrial carbon budget.

Arctic rivers drain ~15% of the global land surface and significantly influence local communities and economies, freshwater and marine ecosystems, and global climate. However, trusted and public knowledge of pan-Arctic rivers is inadequate, especially for small rivers and across Eurasia, inhibiting understanding of the Arctic response to climate change. Here, we calculate daily streamflow in 486,493 pan-Arctic river reaches from 1984-2018 by assimilating 9.18 million river discharge estimates made from 155,710 satellite images into hydrologic model simulations. We reveal larger and more heterogenous total water export (3-17% greater) and water export acceleration (factor of 1.2-3.3 larger) than previously reported, with substantial differences across basins, ecoregions, stream orders, human regulation, and permafrost regimes. We also find significant changes in the spring freshet and summer stream intermittency. Ultimately, our results represent an updated, publicly available, and more accurate daily understanding of Arctic rivers uniquely enabled by recent advances in hydrologic modeling and remote sensing.

All hydrological models need to be calibrated to obtain satisfactory streamflow simulations. Here we present a novel parameter regionalization approach that involves the optimization of transfer equations linking model parameters to climate and landscape characteristics. The optimization was performed in a fully spatially distributed fashion at high resolution (0.05°), instead of at lumped catchment scale, using an unprecedented database of daily observed streamflow from 4,229 headwater catchments (<5,000 km²) worldwide. The optimized equations were subsequently applied globally to produce parameter maps for the entire land surface including ungauged regions. The approach was evaluated using the Kling‐Gupta efficiency (KGE) and a gridded version of the hydrological model HBV. Tenfold cross validation was used to evaluate the generalizability of the approach and to obtain an ensemble of parameter maps. For the 4,229 independent validation catchments, the regionalized parameters yielded a median KGE of 0.46. The median KGE improvement (relative to uncalibrated parameters) was 0.29, and improvements were obtained for 88% of the independent validation catchments. These scores compare favorably to those from previous large catchment sample studies. The degree of performance improvement due to the regionalized parameters did not depend on climate or topography. Substantial improvements were obtained even for independent validation catchments located far from the catchments used for optimization, underscoring the value of the derived parameters for poorly gauged regions. The regionalized parameters—available via www.gloh2o.org/hbv—should be useful for hydrological applications requiring accurate streamflow simulations.

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.