Uncovering the different scales in deer-forest interactions

Succession and climatic stochasticity induce long-term decline of a forest browser

Removal of predators and creation of early seral habitat have, in many systems, caused substantial population growth of herbivores. Hyperabundant herbivores, in turn, induce cascading ecosystem effects, but few studies have investigated long-term browser density trends in relation to succession and stochastic climate events. Here, we use annual, empirical population estimates of a forest browser to relate forest succession to long-term decline of an herbivore that prefers early seral habitat. From 2007-2021, concurrent with reduced timber harvest, we used line-transect distance sampling to document annual changes in Columbian black-tailed deer (Odocoileus hemionus columbianus) density on a mid-sized (17.3km2) predator-free island. We documented successional changes associated with forest aggradation and decreased forage quality for deer: early successional shrub/scrub habitat declined 3.8%/year; timber volume increased 4.5%/year; and canopy coverage increased 2.5%. In 2007-2008, deer densities were the greatest observed (~44/km2), but then an historic snowstorm reduced deer density by 39%. From 2010-2021, as forests continued to mature, deer density decreased 4.0%/year, declining to 20 deer/km2. Using a multivariate approach to combine habitat variables (i.e., early seral coverage, timber volume, and canopy closure) into a measure of forest maturation, we found a significant negative relationship between deer density and forest aggradation. Thus, consistent with predictions for bottom-up limited browsers, we observed significant annual declines in a deer population throughout an extended period of forest regrowth. Despite declines, deer density on the island exceeds mainland densities, and overbrowsing likely continues to disrupt ecosystem processes.

Comparative genomics and genome-wide SNPs of endangered Eld's deer provide breeder selection for inbreeding avoidance

Eld's deer, a conserved wildlife species of Thailand, is facing inbreeding depression, particularly in the captive Siamese Eld's deer (SED) subspecies. In this study, we constructed genomes of a male SED and a male Burmese Eld's deer (BED), and used genome-wide single nucleotide polymorphisms to evaluate the genetic purity and the inbreeding status of 35 SED and 49 BED with limited pedigree information. The results show that these subspecies diverged approximately 1.26 million years ago. All SED were found to be purebred. A low proportion of admixed SED genetic material was observed in some BED individuals. Six potential breeders from male SED with no genetic relation to any female SED and three purebred male BED with no relation to more than 10 purebred female BED were identified. This study provides valuable insights about Eld's deer populations and appropriate breeder selection in efforts to repopulate this endangered species while avoiding inbreeding.

Long-term exclusion of invasive ungulates alters tree recruitment and functional traits but not total forest carbon

Forests are major carbon (C) sinks, but their ability to sequester C and thus mitigate climate change, varies with the environment, disturbance regime, and biotic interactions. Herbivory by invasive, nonnative ungulates can have profound ecosystem effects, yet its consequences for forest C stocks remain poorly understood. We determined the impact of invasive ungulates on C pools, both above- and belowground (to 30 cm), and on forest structure and diversity using 26 paired long-term (>20 years) ungulate exclosures and adjacent unfenced control plots located in native temperate rainforests across New Zealand, spanning 36-41° S. Total ecosystem C was similar between ungulate exclosure (299.93 ± 25.94 Mg C ha^-1 ) and unfenced control (324.60 ± 38.39 Mg C ha^-1 ) plots. Most (60%) variation in total ecosystem C was explained by the biomass of the largest tree (mean diameter at breast height [dbh]: 88 cm) within each plot. Ungulate exclusion increased the abundance and diversity of saplings and small trees (dbh ≥2.5, <10 cm) compared with unfenced controls, but these accounted for ~5% of total ecosystem C, demonstrating that a few, large trees dominate the total forest ecosystem C but are unaffected by invasive ungulates at a timescale of 20-50 years. However, changes in understory C pools, species composition, and functional diversity did occur following long-term ungulate exclusion. Our findings suggest that, although the removal of invasive herbivores may not affect total forest C at the decadal scale, major shifts in the diversity and composition of regenerating species will have longer term consequences for ecosystem processes and forest C.

Top-down and bottom-up forces explain patch utilization by two deer species and forest recruitment

Ungulates play an important role in temperate systems. Through their feeding behaviour, they can respond to vegetation by selecting patches or modify vegetation composition by herbivory. The degree in which they interact with vegetation can either reinforce landscape heterogeneity by creating disturbance or reduce heterogeneity in case of overbrowsing. This study evaluates how bottom-up (patch quality, structure), top-down forces (hunting, distance to village, forest edge) and deer features (feeding type, abundance) mediate patch utilization in a temperate forest and assess the implications of patch utilization and light on forest recruitment. Theory predicts that animals seek to maximize their energetic gains by food intake while minimizing the costs associated to foraging, such as the energy required for avoiding predators and exploiting resources. We focused on two deer species with contrasting feeding type: a browser (C. capreolus) and a mixed feeder (C. elaphus). We paired camera traps to vegetation sub-plots in ten forest sites in the Netherlands that widely ranged in deer abundance and landscape heterogeneity. Results showed that patch utilization is simultaneously explained by bottom-up, top-down forces and by deer abundance, as predicted by the safety-in-numbers hypothesis. Yet, forces best explaining patch utilization differed between deer species. Overall, higher patch utilization came with higher browsing, lower tree diversity and a large difference in forest composition: from a mix of broadleaves and conifers towards only conifers. We conclude that these two deer species, although living in the same area and belonging to the same guild, differentially perceive, interact with and shape their surrounding landscape.

Divergent roles of herbivory in eutrophying forests

Ungulate populations are increasing across Europe with important implications for forest plant communities. Concurrently, atmospheric nitrogen (N) deposition continues to eutrophicate forests, threatening many rare, often more nutrient-efficient, plant species. These pressures may critically interact to shape biodiversity as in grassland and tundra systems, yet any potential interactions in forests remain poorly understood. Here, we combined vegetation resurveys from 52 sites across 13 European countries to test how changes in ungulate herbivory and eutrophication drive long-term changes in forest understorey communities. Increases in herbivory were associated with elevated temporal species turnover, however, identities of winner and loser species depended on N levels. Under low levels of N-deposition, herbivory favored threatened and small-ranged species while reducing the proportion of non-native and nutrient-demanding species. Yet all these trends were reversed under high levels of N-deposition. Herbivores also reduced shrub cover, likely exacerbating N effects by increasing light levels in the understorey. Eutrophication levels may therefore determine whether herbivory acts as a catalyst for the "N time bomb" or as a conservation tool in temperate forests.

Ecological resilience of physical plant-soil feedback to chronic deer herbivory: Slow, partial, but functional recovery

Ecological resilience to ungulate overbrowsing is an important issue in forest ecosystems. After chronic herbivory, the recovery rate of understory vegetation and its related functions can be slow even with decreasing grazing intensity; thus, detecting elasticity during alternative successional trajectories is fundamental to understanding state perturbations. In this context, we focused on physical plant-soil feedback (functional interactions between plant growth and soil physical conditions) and evaluated elasticity and recovery processes according to deer density. The effects of 40 years of chronic herbivory by sika deer (average density 14.7 individuals km^-2 ) on the recovery of understory plant communities and associated changes in soil physical properties in headwater catchments were assessed. Using 8 years of catchment-wide exclusion (fenced) and reduction (only culled; average 4.3 individuals km^-2 ) treatments, plot sampling was conducted in 2010 (before treatment) and 2018 (after treatment). The recovery of vegetation and soil physical properties were evaluated, and functional plant-soil relationships and spatial variability were assessed to detect recovery processes during alternative successional trajectory. Woody species increased only under the exclusion treatment and the average soil bulk density was lower than that under reduction treatments. Soil bulk density was negatively correlated with root biomass in the fenced catchment, and root biomass was positively associated with woody species richness. Reduced soil bulk density (~0.5 g cm^-3 ) was observed with greater root biomass and woody species richness on upper hillslopes in the deer-excluded catchment where plant coverage was minimal. Successional failure under the reduction treatment suggested slow recovery with a depressed threshold according to deer density, indicating a clockwise hysteretic response to deer density. Unlike plant coverage during the earlier period of overbrowsing, woody species root development led the recovery of functional physical plant-soil feedback; however, this was probably limited by the higher soil erosion rate in riparian areas and an under-developed herb layer. Our results highlight an alternative recovery trajectory of physical plant-soil feedback driven by an alternative plant element (woody roots) to recovery trajectory with increasing plant cover. However, riparian erosion and herb layer would still suppress recovery. Therefore, recovery might be slower at the landscape scale.

Density dependence of daily activity in three ungulate species

Daily activity in herbivores reflects a balance between finding food and safety. The safety-in-numbers theory predicts that living in higher population densities increases safety, which should affect this balance. High-density populations are thus expected to show a more even distribution of activity-that is, spread-and higher activity levels across the day. We tested these predictions for three ungulate species; red deer (Cervus elaphus), roe deer (Capreolus capreolus), and wild boar (Sus scrofa). We used camera traps to measure the level and spread of activity across ten forest sites at the Veluwe, the Netherlands, that widely range in ungulate density. Food availability and hunting levels were included as covariates. Daily activity was more evenly distributed when population density was higher for all three species. Both deer species showed relatively more feeding activity in broad daylight and wild boar during dusk. Activity level increased with population density only for wild boar. Food availability and hunting showed no correlation with activity patterns. These findings indicate that ungulate activity is to some degree density dependent. However, while these patterns might result from larger populations feeling safer as the safety-in-numbers theory states, we cannot rule out that they are the outcome of greater intraspecific competition for food, forcing animals to forage during suboptimal times of the day. Overall, this study demonstrates that wild ungulates adjust their activity spread and level based on their population size.