Daily stream temperature predictions for free-flowing streams in the Pacific Northwest, USA

Supporting sustainable lotic ecosystems and thermal habitats for cold-water species like salmonids requires estimates of stream temperature that are high in scope and resolution across space and time. We combined and enhanced elements of existing stream temperature models to produce a new statistical model to address this need. Contrasting with previous models that estimated coarser metrics such as monthly or seasonal stream temperature or focused on individual watersheds, here we model and predict daily stream temperature across the entire calendar year for a broad geographic region. This model reflects mechanistic processes using publicly available climate and landscape covariates in a Generalized Additive Model (GAM) framework. We allowed covariates to interact while accounting for nonlinear relationships between temporal and spatial covariates to better capture seasonal patterns. Additionally, to represent variation in sensitivity to climate, we used a moving average of antecedent air temperatures over a variable duration linked to area-standardized streamflow. The moving average window size was longer for reaches classified as having a snow-dominated hydrology, especially at higher flows, whereas window size was relatively constant and low for reaches having rain-dominated hydrology. Our model’s ability to capture the temporally variable impact of snowmelt on stream sensitivity helped improve its capacity to predict stream temperature across diverse geography for multiple years. We fit the model to stream temperature data from 1993-2013 and used the model to predict daily stream temperatures for ~261,200 free-flowing stream reaches across the Pacific Northwest from 1990-2021. Our daily model fit well (RMSE = 1.78; MAE = 1.33 °C). Spatial and temporal cross-validation suggested that the model produced useful predictions at un-sampled locations and across diverse landscapes and climate conditions. We produced stream temperature predictions that will be immediately useful to natural resource practitioners in the Pacific Northwest, USA, especially for effective conservation planning in lotic ecosystems and for managing species such as Pacific salmon. Our approach is straightforward and can be easily adapted to new spatial regions, time periods, or scenarios such as the anticipated decline in snowpack with climate change.Model output - daily stream temperature predictions for each NHD 1:100k reach - are archived with Zenodo in the Riverscapes Consortium community. Prediction data is available in zip (7z) archives by 6-digit HUC for CONUS hydroregion 17 for the time period 1990-2021. Archives are 10-digit HUC based csv format files (COMID, year, day, stream temperature prediction), bundled by 6-digit HUC. All code is available at https://github.com/aimeefullerton/daily-st-pnw under the CC-BY-4.0 license. All model input data will be available as a riverscape compliant project on the Riverscape Exchange data portal and web-based explorer