Non-breeding waterbirds benefit from protected areas when adjusting their distribution to climate warming

Climate warming is driving changes in species distributions, although many species show a so-called climatic debt, where their range shifts lag behind the fast shift in temperature isoclines. Protected areas (PAs) may impact the rate of distribution changes both positively and negatively. At the cold edges of species distributions, PAs can facilitate species distribution changes by increasing the colonization required for distribution change. At the warm edges, PAs can mitigate the loss of species, by reducing the local extinction of vulnerable species. To assess the importance of PAs to affect species distribution change, we evaluated the changes in a non-breeding waterbird community as a response to temperature increase and PA status, using changes of species occurrence in the Western Palearctic over 25 years (97 species, 7,071 sites, 39 countries, 1993 2017). We used a community temperature index (CTI) framework based on species thermal affinities to investigate the species turn-over induced by temperature increase. In addition, we measured whether the thermal community adjustment was led by cold-dwelling species extinction and/or warm-dwelling species colonization, by modelling the change in standard deviation of the CTI (CTIsd). Using linear mixed-effects models, we investigated whether communities within PAs had lower climatic debt and different patterns of community change regarding the local PA surface. Thanks to the combined use of the CTI and CTIsd, we found that communities inside PAs had more species, higher colonization, lower extinction and the climatic debt was 16% lower than outside PAs. The results suggest the importance of PAs to facilitate warm-dwelling species colonization and attenuate cold-dwelling species extinction. The community adjustment was however not sufficiently fast to keep pace with the strong temperature increase in central and northeastern Western Palearctic regions. Our study underlines the potential of the combined CTI and CTIsd metrics to understand the colonization-extinction patterns driven by climate warming.