RESUMEN
Rapid changes in thermal environments are threatening many species worldwide. Thermal acclimatisation may partially buffer species from the impacts of these changes, but currently, the knowledge about the temporal dynamics of acclimatisation remains limited. Moreover, acclimatisation phenotypes are typically determined in laboratory conditions that lack the variability and stochasticity that characterise the natural environment. Through a distributed lag non-linear model (DLNM), we use field data to assess how the timing and magnitude of past thermal exposures influence thermal tolerance. We apply the model to two Scottish freshwater Ephemeroptera species living in natural thermal conditions. Model results provide evidence that rapid heat hardening effects are dramatic and reflect high rates of change in temperatures experienced over recent hours to days. In contrast, temperature change magnitude impacted acclimatisation over the course of weeks but had no impact on short-term responses. Our results also indicate that individuals may de-acclimatise their heat tolerance in response to cooler environments. Based on the novel insights provided by this powerful modelling approach, we recommend its wider uptake among thermal physiologists to facilitate more nuanced insights in natural contexts, with the additional benefit of providing evidence needed to improve the design of laboratory experiments.