“…While time-varying covariate effects can be represented using separate models for each season (Mohseni et al, 1998;Sohrabi et al, 2017), this may cause unnatural, abrupt changes at seasonal transitions. Time-varying coefficients, including those used in generalized additive models (GAMs) (Pedersen et al, 2019;Wood, 2017) use continuous functions that avoid these abrupt changes (Li et al, 2014;Jackson et al, 2018;Siegel & Volk, 2019). While GAMs have been used in daily stream temperature modeling for single-site prediction (Boudreault et al, 2019;Coleman et al, 2021;Glover et al, 2020;Laanaya et al, 2017;Siegel et al, 2022), spatiotemporal prediction (Jackson et al, 2018;Siegel & Volk, 2019), identifying extreme events (Georges et al, 2021), and trend assessment (Yang & Moyer, 2020), few studies have used GAMs to model seasonally varying flow effects or identify when stream temperatures are most affected by flow variation (Glover et al, 2020;Yang & Moyer, 2020).…”