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.

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.

We mapped daily streamflow from 1984 to 2018 in approximately 2.9 million rivers to assess recent changes to global river systems. We found that river outlets were dominated by significant decreases in flow, whereas headwaters were 1.7 times more likely to have significantly increased flow than decreased. These changes result in a significant upstream shift in streamflow experienced by about 29% of the global land surface. We found the most changes in the smallest steams in our study: increases in erosion potential (approximately 5% increase in stream power), flood frequency (approximately 42% increase in 100-year floods), and likely nutrient dynamics (altered seasonal flow regimes). We revealed these changes using “detail at scale” by mapping millions of individual rivers. Widely adopting this approach could reveal other changes to the hydrosphere.

Recent advances in remote sensing and the upcoming launch of the joint NASA/CNES/CSA/UKSA Surface Water and Ocean Topography (SWOT) satellite point toward improved river discharge estimates in ungauged basins. Existing discharge methods rely on “prior river knowledge” to infer parameters not directly measured from space. Here, we show that discharge estimation is improved by classifying and parameterizing rivers based on their unique geomorphology and hydraulics. Using over 370,000 in situ hydraulic observations as training data, we test unsupervised learning and an “expert” method to assign these hydraulics and geomorphology to rivers via remote sensing. This intervention, along with updates to model physics, constitutes a new method we term “geoBAM,” an update of the Bayesian At‐many‐stations hydraulic geometry‐Manning's (BAM) algorithm. We tested geoBAM on Landsat imagery over more than 7,500 rivers (108 are gauged) in Canada's Mackenzie River basin and on simulated hydraulic data for 19 rivers that mimic SWOT observations without measurement error. geoBAM yielded considerable improvement over BAM, improving the median Nash‐Sutcliffe efficiency (NSE) for the Mackenzie River from −0.05 to 0.26 and from 0.16 to 0.46 for the SWOT rivers. Further, NSE improved by at least 0.10 in 78/108 gauged Mackenzie rivers and 8/19 SWOT rivers. We attribute geoBAM improvement to parameterizing rivers by type rather than globally, but prediction accuracy worsens if parameters are misassigned. This method is easily mapped to rivers at the global scale and paves the way for improving future discharge estimates, especially when coupled with hydrologic models.

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.