Regional thermal analysis approach: A management tool for predicting water temperature metrics relevant for thermal fish habitat.

• Published in: Ecological Informatics
• Main Authors: Abidi, Olfa, St-Hilaire, André, Ouarda, Taha B. M. J., Charron, Christian, Boyer, Claudine, Daigle, Anik
• Format: Article in Journal/Newspaper
• Language: English
• Published: 2022
• Subjects: river water temperature; regional thermal analysis; multiple linear regression model MLR; generalized additive model GAM; thermal homogeneous regions; HCA; Canada; Gam; Atlantic salmon; Salmo salar
• Online Access: https://espace.inrs.ca/id/eprint/12628/; https://espace.inrs.ca/id/eprint/12628/1/P4131.pdf; https://doi.org/10.1016/j.ecoinf.2022.101692

The field of water resource management, including fisheries, is facing new challenges associated with climate change. This study sheds light on the modeling of water temperature indices (metrics) that describe critical thermal maxima of the Atlantic salmon (salmo salar). These thermal metrics include MaxWaterTmax (interannual mean of maximum summer temperature), MaxNumDay (interannual mean of the number of consecutive days with maximum water temperature > 25 °C and minimum water temperature > 20 °C). The latter is an important indicator to evaluate thermal variability. Three other parameters of a Gaussian function fitted to the interannual daily mean temperatures characterizing the thermal regime of 146 stations located in Eastern Canada were estimated. These three parameters are Gaussian_a (maximum of interannual daily mean temperature), Gaussian_b (mean duration of the warm period), and Gaussian_c (date of occurrence of the interannual maximum temperature). The classical Multiple linear regression model (MLR) and the non-linear Generalized additive model (GAM) were tested and compared to estimate the five thermal metrics. The regression-based approaches involve the identification of thermally homogeneous regions based on three approaches: hierarchical clustering analysis (HCA), regions of influence (ROI) as well as canonical correlation analysis (CCA). Then, the regional MLR and GAM models were applied within the delineated homogenous regions. Also, the regional models were compared to models encompassing all stations (i.e., one region). For each regional estimation model and each thermal metric, a set of optimal explanatory variables were selected using a forward stepwise procedure. The database consisted of 22 environmental predictors related to physiography, topography, climate, land cover and surface deposits. To assess performance of the models, the following statistical metrics were used: coefficient of determination R², root mean square error (RMSE), bias, relative root mean square error (RRMSE) ...