Linking spatial patterns of terrestrial herbivore community structure to trophic interactions

  1. Jakub Witold Bubnicki  Is a corresponding author
  2. Marcin Churski
  3. Krzysztof Schmidt
  4. Tom A Diserens
  5. Dries PJ Kuijper
  1. Mammal Research Institute, Polish Academy of Sciences, Poland
51 figures, 2 tables and 3 additional files

Figures

Maps of the study area (left) and camera trapping sampling design (right).

On the study area map we marked the major settlements connected by public roads open for cars.

https://doi.org/10.7554/eLife.44937.003
Maps of candidate landscape-scale covariates for the ungulate and carnivore models representing environmental and human-related landscape gradients.

All covariates were scaled and zero-centered prior to modelling. POIs – density of tourist infrastructure. When selecting covariates for the final models we ensured that the Pearson correlation …

https://doi.org/10.7554/eLife.44937.004
Spatial distribution of large carnivores is determined by human activity.

(a, d) Estimated parameters for the wolf and lynx models, (b, e) maps of the predicted variation in their landscape use (relative density surfaces) and (c, f) fitted spatial random effects (SRE). …

https://doi.org/10.7554/eLife.44937.005
Spatial distribution of each ungulate species was related to a unique combination of bottom-up and top-down factors.

(a) Estimated effects of covariates explaining the spatial variation in the parameter λ, which is the expected number of individuals using a given landscape grid cell (25 ha pixel) during the …

https://doi.org/10.7554/eLife.44937.006
Only red deer (the major wolf prey) showed a negative spatial association with parts of the landscape intensively used by wolves, whereas distributions of other species were shaped by bottom-up factors.

The maps show substantial variation in predicted landscape use for the five studied ungulate species (relative density surfaces). Specifically, this figure presents the spatial predictions for the …

https://doi.org/10.7554/eLife.44937.007
All ungulate species except red deer males showed some level of remaining spatial auto-correlation not explained by the raster covariates included in the models.

The maps show the fitted spatial random effects (SRE) for the models of five studied ungulate species. The SRE are deviations from 0 at log-scale. The SRE captured all spatial variation not …

https://doi.org/10.7554/eLife.44937.008
Factors influencing detection rates of ungulate species.

(a) Estimated effects of covariates explaining the variation in the parameter γ, which is the expected (daily) detection (trapping) rate at a single camera trap site (see the general and formal …

https://doi.org/10.7554/eLife.44937.009
There was substantial spatial variation in the structure of the ungulate community, as indicated by the weak spatial associations between species.

The pairwise spatial overlap between the studied ungulate species is presented as the Pearson correlation heatmap of their relative density surfaces predicted by the models (a) and fitted spatial …

https://doi.org/10.7554/eLife.44937.010
Maps of estimated total biomass [kg/pixel] and functional diversity (FDis).

Both variables, together with species-specific biomass raster layers, were used in the hierarchical clustering analysis. Interestingly, the lowest values of both parameters were largely associated …

https://doi.org/10.7554/eLife.44937.011
Spatial variation and functional diversity of a herbivore community within a temperate forest herbivome, decomposed into landscape-scale herbivory regimes – herbiscapes.

(a) The spatial distribution of the five identified clusters (‘herbiscapes’) based on hierarchical clustering on the principal components analysis. (b) The distribution of identified clusters in the …

https://doi.org/10.7554/eLife.44937.012
The patterns of browsing intensity and the composition of regenerating tree species differ between herbiscapes and between the protected and managed parts of the BF landscape.

(a, b) The mean differences between the two tree height classes (saplings < 30 cm and ≥30 cm) in the proportional shares of Acer platanoides (Acer) and Carpinus betulus (Carpinus). The numbers on …

https://doi.org/10.7554/eLife.44937.013
Graphical abstract presenting the steps taken to reveal the relationships uncovered in this study.

Image credit: Lisa Sánchez Aguilar

https://doi.org/10.7554/eLife.44937.014
Ranked relative densities predicted for all ungulate species and for each 25 ha landscape pixel.

All plots follow a characteristic shape of a hollow or sigmoidal curve (Brown et al., 1995; Rocchini and Neteler, 2012) indicating a non-uniform, heterogenous distribution of all studied species.

https://doi.org/10.7554/eLife.44937.015
Scatter plot showing how wolves, by affecting the spatial distribution of red deer, re-structured the composition of the entire community of large herbivores.

The predicted biomass of red deer (both sexes, X axis) is shown as a share of the total biomass of the entire community of ungulates. Each dot on the plot is a 25 ha landscape pixel. In color, we …

https://doi.org/10.7554/eLife.44937.016
Appendix 1—figure 1
Heatmap of pairwise Pearson correlation values of all candidate landscape-scale covariates for the ungulates models; wolf: predicted use of the landscape by the wolf, lynx: predicted use of the landscape by the lynx, open: landscape openness, conif: % share of coniferous tree species, decid: % share of deciduous tree species, fheight: forest height, reserv: density of protected areas, elev: elevation, dsetta: distance to all settlements, dsettm: distance to major settlements, droads: distance to main roads, dedge: distance to forest edge, pois: density of tourist infrastructure.

When selecting covariates for the final models we ensured that the Pearson correlation value for all pairs of included covariates was lower than 0.7.

https://doi.org/10.7554/eLife.44937.023
Appendix 1—figure 2
The results of hierarchical clustering on the principal component analysis.

The five clusters were identified by assessing the inertia (i.e. change in within cluster homogeneity) gained by cutting the tree at different levels and the ecological interpretability of the …

https://doi.org/10.7554/eLife.44937.024
Appendix 1—figure 3
Spatial patterns of human activity in Białowieża Forest derived from the publicly available world-wide data provided by STRAVA (https://www.strava.com/heatmap) combined with (A) predicted use of the landscape by the wolf, (B) predicted use of the landscape by the lynx, (C) predicted use of landscape by the red deer (both sexes) and (D) distance to major settlements used in the models for both large carnivores.

The dark-red thick lines indicate high human activity and the light-blue thin lines indicate low human activity. For (A), (B) and (C) the color gradient from blue to red corresponds to the gradient …

https://doi.org/10.7554/eLife.44937.025
Appendix 1—figure 4
Comparison of predictions of wolf space use with existing radio-tracking data from collared wolves in Białowieża Forest (period 1994–1999, Jędrzejewski et al., 2007) and observations of known den locations (period 1993–2007, MRI PAS unpublished data).

In the background we plotted the raster map of wolf space use predicted by our model together with the raw camera trapping data (‘bubbles’ with trapping rates at each location). On top of it we …

https://doi.org/10.7554/eLife.44937.026
Appendix 1—figure 5
Bubble plot presenting raw camera trapping data for lynx collected during the carnivore survey.

Each bubble is a camera trap location and its size is proportional to the daily trapping rate (i.e. no. of individuals observed/no. of days).

https://doi.org/10.7554/eLife.44937.027
Appendix 1—figure 6
Bubble plot presenting raw camera trapping data for red deer female collected during the ungulate survey.

Each bubble is a camera trap location and its size is proportional to the daily trapping rate (i.e. no. of individuals observed/no. of days).

https://doi.org/10.7554/eLife.44937.028
Appendix 1—figure 7
Bubble plot presenting raw camera trapping data for red deer male collected during the ungulate survey.

Each bubble is a camera trap location and its size is proportional to the daily trapping rate (i.e. no. of individuals observed/no. of days).

https://doi.org/10.7554/eLife.44937.029
Appendix 1—figure 8
Bubble plot presenting raw camera trapping data for wild boar collected during the ungulate survey.

Each bubble is a camera trap location and its size is proportional to the daily trapping rate (i.e. no. of individuals observed/no. of days).

https://doi.org/10.7554/eLife.44937.030
Appendix 1—figure 9
Bubble plot presenting raw camera trapping data for roe deer collected during the ungulate survey.

Each bubble is a camera trap location and its size is proportional to the daily trapping rate (i.e. no. of individuals observed/no. of days).

https://doi.org/10.7554/eLife.44937.031
Appendix 1—figure 10
Bubble plot presenting raw camera trapping data for bison collected during the ungulate survey.

Each bubble is a camera trap location and its size is proportional to the daily trapping rate (i.e. no. of individuals observed/no. of days).

https://doi.org/10.7554/eLife.44937.032
Appendix 1—figure 11
Bubble plot presenting raw camera trapping data for moose collected during the ungulate survey.

Each bubble is a camera trap location and its size is proportional to the daily trapping rate (i.e. no. of individuals observed/no. of days).

https://doi.org/10.7554/eLife.44937.033
Appendix 1—figure 12
Basic meteorological data for the study period (05/2012 – 09/2014): 1) averaged daily temperature [C], 2) precipitation [mm] and 3) snow cover [cm].
https://doi.org/10.7554/eLife.44937.034
Appendix 1—figure 13
Distribution of camera trap locations grouped by quality of view in front of a camera and different UNESCO zones in Białowieża Forest (1: Strict Reserve of BNP, 2: BNP and nature reserves, 3: valuable unprotected forest stands, 4: managed forest stands).

Each camera location was manually tagged with a 4-level categorical label (‘Exclude’, ‘Acceptable’, ‘Good’, ‘Perfect’) describing the quality of view in front of the camera. Locations labeled with …

https://doi.org/10.7554/eLife.44937.035
Appendix 1—figure 14
Red deer female – landscape use; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.036
Appendix 1—figure 15
Red deer female – detection rate; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.037
Appendix 1—figure 16
Red deer male – landscape use; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.038
Appendix 1—figure 17
Red deer male – detection rate; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.039
Appendix 1—figure 18
Roe deer – landscape use; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.040
Appendix 1—figure 19
Roe deer – detection rate; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.041
Appendix 1—figure 20
Wild boar – landscape use; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.042
Appendix 1—figure 21
Wild boar – detection rate; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.043
Appendix 1—figure 22
European bison – landscape use; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.044
Appendix 1—figure 23
European bison – detection rate; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.045
Appendix 1—figure 24
Moose – landscape use; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.046
Appendix 1—figure 25
Moose – detection rate; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.047
Appendix 1—figure 26
Wolf – landscape use; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.048
Appendix 1—figure 27
Eurasian lynx – landscape use; the posterior distributions of all parameters.
https://doi.org/10.7554/eLife.44937.049
Appendix 1—figure 28
Red deer female – model diagnostic plots.
https://doi.org/10.7554/eLife.44937.050
Appendix 1—figure 29
Red deer male – model diagnostic plots.
https://doi.org/10.7554/eLife.44937.051
Appendix 1—figure 30
Roe deer – model diagnostic plots.
https://doi.org/10.7554/eLife.44937.052
Appendix 1—figure 31
Wild boar – model diagnostic plots.
https://doi.org/10.7554/eLife.44937.053
Appendix 1—figure 32
European bison – model diagnostic plots.
https://doi.org/10.7554/eLife.44937.054
Appendix 1—figure 33
Moose – model diagnostic plots.
https://doi.org/10.7554/eLife.44937.055
Appendix 1—figure 34
Wolf – model diagnostic plots.
https://doi.org/10.7554/eLife.44937.056
Appendix 1—figure 35
Eurasian lynx – model diagnostic plots.
https://doi.org/10.7554/eLife.44937.057
Appendix 1—figure 36
Estimated effects of the covariates from the model for all ungulate species using only a subset of camera trap data covering a three month period (August - October) closely matching the period of the carnivore survey.

The covariates explain the spatial variation in the parameter λ, which is the expected number of individuals using a given landscape grid cell (25 ha pixel) during the sampling period (i.e. the …

https://doi.org/10.7554/eLife.44937.058
Appendix 1—figure 37
The spatial predictions for the parameter λ, which is the expected number of individuals using a given landscape grid cell (25 ha pixel) during the sampling period (see the general and formal description of the model in the Materials and methods section).

The predictions are based on the model for all ungulate species using only a subset of the camera trap data covering a three month period (August - October) closely matching the period of the …

https://doi.org/10.7554/eLife.44937.059

Tables

Appendix 1—table 1
Summary of the two camera trapping surveys conducted during the study period.

Days – the total number of trapping days (effort), Counts – the total number of individuals recorded during independent visits (events), Trate – trapping rate. Notice the dramatic differences in the …

https://doi.org/10.7554/eLife.44937.021
SpeciesDaysCountsTrate (daily)Trate sdTrate max
Ungulates survey (05.2012–05.2014)
Wild Boar981348180.495070.8788110.16667
Red Deer981320250.202550.363403.90909
Red Deer Female98139730.096480.203912.00000
Red Deer Male98135470.055440.128671.25000
European Bison98134400.045220.232124.50000
Roe Deer98131940.020430.066110.54545
Eurasian Elk9813820.008660.047140.66667
Wolf9813510.005410.029600.36364
Eurasian Lynx981330.000270.004580.08333
Carnivores survey (09–10.2015)
Wolf30934710.191290.495346.00000
Eurasian Lynx3093350.012680.049970.33333
Appendix 1—table 2
Matrix of the foraging-related traits used for the estimation of the functional diversity index (FDis) of the entire ungulate community: body mass (based on Jedrzejewska et al., 1994) and Borowik et al., 2016), diet type (based on Hofmann, 1989); B: browsers-concentrate selectors, BG: browsers/grazers-intermediate types, G: grazers-grass roughage eaters, O: omnivores) and gut type (R: ruminants, NR: non-ruminants).
https://doi.org/10.7554/eLife.44937.022
Roe deerWild boarRed deer femaleRed deer maleMooseEuropean bison
Body mass208090150200400
Diet typeBOBGBGBG
Gut typeRNRRRRR

Additional files

Download links