Legacy effects of herbivory on treeline dynamics along an elevational gradient

Treelines are expected to expand into alpine ecosystems with global warming, but herbivory may delay this expansion. This study quantifies long-term effects of temporally varying sheep densities on birch recruitment and growth in the treeline ecotone. We examined treeline ecotone successional trajectories and legacy effects in a replicated experimental setup, where enclosures were present for 14 years with three different sheep densities (0, 25, 80 sheep km(−2)). Before and after the enclosures were present, the site had an ambient sheep density of 20–25 km(−2). We sampled field data 4 years after enclosure removal and compared these to data sampled 8 and 9 years after enclosure erection. We sampled data on birch browsing pressure, birch distribution across life-stages (recruits, saplings, and mature trees), and birch annual radial growth. Fourteen years of increased or decreased sheep density had observable legacy effects depending on birch life-stage. Birch recruit prevalence decreased in areas, where sheep were reintroduced after being absent for 14 years. For the same areas, sapling and mature tree prevalence increased, indicating that these areas have entered alternative successional trajectories compared to areas, where sheep were present the whole time. Birch annual radial growth showed a lag effect of 2 years after enclosure removal, with growth decreasing in areas where sheep had been absent for 14 years and increasing where sheep densities were high. Thus, decadal-scale absences of herbivores can leave legacy effects due to increased numbers of trees that have high resistance to later-introduced herbivore browsing. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00442-022-05125-8.