skip to main content


Title: A High‐Resolution, Daily Hindcast (1990–2021) of Alaskan River Discharge and Temperature From Coupled and Optimized Physical Models
Abstract

Water quality and freshwater ecosystems are affected by river discharge and temperature. Models are frequently used to estimate river temperature on large spatial and temporal scales due to limited observations of discharge and temperature. In this study, we use physically based river routing and temperature models to simulate daily discharge and river temperature for rivers in 138 basins in Alaska, including the entire Yukon River basin, from 1990–2021. The river temperature model was optimized for ice free months using a surrogate‐based model optimization method, improving model performance at uncalibrated river gages. A common statistical model relating local air and water temperature was used as a benchmark. The physically based river temperature model exhibited superior performance compared to the benchmark statistical model after optimization, suggesting river temperature model optimization could become more routine. The river temperature model demonstrated high sensitivity to air temperature and model parameterization, and lower sensitivity to discharge. Validation of the models showed a Kling‐Gupta Efficiency of 0.46 for daily river discharge and a root mean square error of 2.04°C for daily river temperature, improving on the non‐optimized physical model and the benchmark statistical model, which had root mean square errors of 3.24 and 2.97°C, respectively. The simulation shows that rivers in northern Alaska have higher maximum summer temperatures and more variability than rivers in the Central and Southern regions. Furthermore, this framework can be readily adapted for use across models and regions.

 
more » « less
NSF-PAR ID:
10504764
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
60
Issue:
4
ISSN:
0043-1397
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Many canyon‐bound rivers have been dammed and downstream flow and water temperatures modified. In some regions, climate change is expected to cause lower storage in reservoirs and warmer release temperatures, which may further alter downstream flow and thermal regimes. To anticipate potential future changes, we first need to understand the dominant heat transfer mechanisms in canyon‐bound river systems. Toward this end, we adapt a dynamic process‐based river routing and temperature model to account for complex shading and radiation characteristics found in canyon‐bound rivers. We apply the model to a 362 km segment of the Colorado River in Grand Canyon National Park, USA to simulate temperature over an 18‐year period. Extensive temperature and flow data sets from within the canyon were used to assess model performance. At the most downstream gaging location, root mean square errors of hourly flow routing and temperature predictions were 11.5 m3/s and 0.93°C, respectively. We found that heat fluxes controlling temperatures were highly variable over space and time, primarily due to shortwave radiation dynamics and hydropeaking flow conditions. Additionally, the large differences between air and water temperature during summer periods resulted in high sensible and latent heat fluxes. Sensitivity analyses indicate that reservoir release temperatures are most influential above the RM88 gage (141 km below Glen Canyon Dam), while a combination of discharge, shortwave radiation, and air temperature become more important farther downstream. This study illustrates the importance of understanding the spatial and temporal variability of topographic shading when predicting water temperatures in canyon‐bound rivers.

     
    more » « less
  2. Abstract

    Water temperature is a critical ecological indicator; however, few studies have statistically modeled century‐scale trends in riverine or estuarine water temperature, or their cause. Here, we recover, digitize, and analyze archival temperature measurements from the 1850s onward to investigate how and why water temperatures in the lower Columbia River are changing. To infill data gaps and explore changes, we develop regression models of daily historical Columbia River water temperature using time‐lagged river flow and air temperature as the independent variables. Models were developed for three time periods (mid‐19th, mid‐20th, and early 21st century), using archival and modern measurements (1854–1876; 1938–present). Daily and monthly averaged root‐mean‐square errors overall are 0.89°C and 0.77°C, respectively for the 1938–2018 period. Results suggest that annual averaged water temperature increased by 2.2°C ± 0.2°C since the 1850s, a rate of 1.3°C ± 0.1°C/century. Increased water temperatures are seasonally dependent. An increase of approximately 2.0°C ± 0.2°C/century occurs in the July–Dec time‐frame, while springtime trends are statistically insignificant. Rising temperatures change the probability of exceeding ecologically important thresholds; since the 1850s, the number of days with water temperatures over 20°C increased from ~5 to 60 per year, while the number below 2°C decreased from ~10 to 0 days/per year. Overall, the modern system is warmer, but exhibits less temperature variability. The reservoir system reduces sensitivity to short‐term atmospheric forcing. Statistical experiments within our modeling framework suggest that increased water temperature is driven by warming air temperatures (~29%), altered river flow (~14%), and water resources management (~57%).

     
    more » « less
  3. Abstract

    Bedrock rivers are the pacesetters of landscape evolution in uplifting fluvial landscapes. Water discharge variability and sediment transport are important factors influencing bedrock river processes. However, little work has focused on the sensitivity of hillslope sediment supply to precipitation events and its implications on river evolution in tectonically active landscapes. We model the temporal variability of water discharge and the sensitivity of sediment supply to precipitation events as rivers evolve to equilibrium over 106model years. We explore how coupling sediment supply sensitivity with discharge variability influences rates and timing of river incision across climate regimes. We find that sediment supply sensitivity strongly impacts which water discharge events are the most important in driving river incision and modulates channel morphology. High sediment supply sensitivity focuses sediment delivery into the largest river discharge events, decreasing rates of bedrock incision during floods by orders of magnitude as rivers are inundated with new sediment that buries bedrock. The results show that the use of river incision models in which incision rates increase monotonically with increasing river discharge may not accurately capture bedrock river dynamics in all landscapes, particularly in steep landslide prone landscapes. From our modeling results, we hypothesize the presence of an upper discharge threshold for river incision at which storms transition from being incisional to depositional. Our work illustrates that sediment supply sensitivity must be accounted for to predict river evolution in dynamic landscapes. Our results have important implications for interpreting and predicting climatic and tectonic controls on landscape morphology and evolution.

     
    more » « less
  4. Abstract

    The evasion of CO2from inland waters, a major carbon source to the atmosphere, depends on dissolved inorganic carbon (DIC) concentrations. Our understanding of DIC dynamics across gradients of climate, geology, and vegetation conditions however have remained elusive. To understand its large‐scale patterns and drivers, we collated instantaneous and mean (multiyear average) DIC concentrations from about 100 rivers draining minimally‐impacted watersheds in the contiguous United States. Within individual sites, instantaneous concentrations (C) measured at daily to seasonal scales exhibit a near‐universal response to changes in river discharge (Q) in a negative power law form. High concentrations occur at low discharge when DIC‐enriched groundwater dominates river discharge; low concentrations occur under high flow when relatively DIC‐poor shallow soil water predominates river discharge. Such patterns echo the widely observed increase of soil CO2and DIC with depth and the shallow‐and‐deep hypothesis that emphasizes the importance of flow paths and source water chemistry. Across sites, mean concentrations (Cm) decrease with increasing mean discharge (Qm), a long‐term climate measure, and reachs maxima at around 200 mm/yr. A parsimonious model reveals that high mean DIC arises from soil CO2accumulation when rates of DIC‐generating reactions are relatively high compared to its export fluxes in arid climates. Although instantaneous and mean DIC concentrations similarly decrease with increasing discharge, results here highlight their distinct drivers: daily to seasonal‐scale instantaneous concentration variations (C) are controlled by subsurface CO2distribution over depth (from below), whereas long‐term mean concentrations (Cm) are regulated by climate (from above). The results emphasize the significance of land‐river connectivity via subsurface flow paths. They also underscore the importance of characterizing subsurface CO2distribution to illuminate belowground processes in order to project the future of water and carbon cycles in a warming climate.

     
    more » « less
  5. Abstract

    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.

     
    more » « less