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.

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.

Abstract. Assessing impacts of climate change on hydrologic systems is critical for developing adaptation and mitigation strategies for water resource management, risk control, and ecosystem conservation practices. Such assessments are commonly accomplished using outputs from a hydrologic model forced with future precipitation and temperature projections. The algorithms used for the hydrologic model components (e.g., runoff generation) can introduce significant uncertainties into the simulated hydrologic variables. Here, a modeling framework was developed that integrates multiple runoff generation algorithms with a routing model and associated parameter optimizations. This framework is able to identify uncertainties from both hydrologic model components and climate forcings as well as associated parameterization. Three fundamentally different runoff generation approaches, runoff coefficient method (RCM, conceptual), variable infiltration capacity (VIC, physically based, infiltration excess), and simple-TOPMODEL (STP, physically based, saturation excess), were coupled with the Hillslope River Routing model to simulate surface/subsurface runoff and streamflow. A case study conducted in Santa Barbara County, California, reveals increased surface runoff in February and March but decreased runoff in other months, a delayed (3 d, median) and shortened (6 d, median) wet season, and increased daily discharge especially for the extremes (e.g., 100-year flood discharge, Q100). The Bayesian model averaging analysis indicates that the probability of such an increase can be up to 85 %. For projected changes in runoff and discharge, general circulation models (GCMs) and emission scenarios are two major uncertainty sources, accounting for about half of the total uncertainty. For the changes in seasonality, GCMs and hydrologic models are two major uncertainty contributors (∼35 %). In contrast, the contribution of hydrologic model parameters to the total uncertainty of changes in these hydrologic variables is relatively small (<6 %), limiting the impacts of hydrologic model parameter equifinality in climate change impact analysis. This study provides useful information for practices associated with water resources, risk control, and ecosystem conservation and for studies related to hydrologic model evaluation and climate change impact analysis for the study region as well as other Mediterranean regions.

Incorporating coastal ecosystems in climate adaptation planning is needed to maintain the well-being of both natural and human systems. Our vulnerability study uses a multidisciplinary approach to evaluate climate change vulnerability of an urbanized coastal community that could serve as a model approach for communities worldwide, particularly in similar Mediterranean climates. We synthesize projected changes in climate, coastal erosion and flooding, watershed runoff and impacts to two important coastal ecosystems, sandy beaches and coastal salt marshes. Using downscaled climate models along with other regional models, we find that temperature, extreme heat events, and sea level are expected to increase in the future, along with more intense rainfall events, despite a negligible change in annual rainfall. Consequently, more droughts are expected but the magnitude of larger flood events will increase. Associated with the continuing rise of mean sea level, extreme coastal water levels will occur with increasingly greater magnitudes and frequency. Severe flooding will occur for both natural (wetlands, beaches) and built environments (airport, harbor, freeway, and residential areas). Adaptation actions can reduce the impact of rising sea level, which will cause losses of sandy beach zones and salt marsh habitats that support the highest biodiversity in these ecosystems, including regionally rare and endangered species, with substantial impacts occurring by 2050. Providing for inland transgression of coastal habitats, effective sediment management, reduced beach grooming and removal of shoreline armoring are adaptations that would help maintain coastal ecosystems and the beneficial services they provide.