Simple Bayesian reconstruction and forecasting of stream water temperature for ecologists—A tool using air temperature, optionally flow, in a time series decomposition approach

Mitigating the impacts of global warming on wildlife entails four practical steps. First, we need to study how processes of interest vary with temperature. Second, we need to build good temperature scenarios. Third, processes can be forecast accordingly. Only then can we perform the fourth step, testing mitigating measures. While having good temperature data is essential, this is not straightforward for stream ecologists and managers. Water temperature (WT) data are often short and incomplete and future projections are currently not routinely available. There is a need for generic models which address this data gap with good resolution and current models are partly lacking. Here, we expand a previously published hierarchical Bayesian model that was driven by air temperature (AT) and flow (Q) as a second covariate. The new model can hindcast and forecast WT time series at a daily time step. It also allows a better appraisal of real uncertainties in the warming of water temperatures in rivers compared to the previous version, stemming from its hybrid structure between time series decomposition and regression. This model decomposes all-time series using seasonal sinusoidal periodic signals and time varying means and amplitudes. It then links the contrasted frequency signals of WT (daily and six month) through regressions to that of AT and optionally Q for better resolution. We apply this model to two contrasting case study rivers. For one case study, AT only is available as a covariate. This expanded model further improves the already good fitting and predictive capabilities of its earlier version while additionally highlighting warming uncertainties. The code is available online and can easily be run for other temperate rivers.