Data from: The primacy of densitymediated indirect effects in a community of wolves, elk, and aspen
Data files
Aug 02, 2024 version files 3.02 MB

BriceEtAl2024_Data.xlsx

BriceEtAl2024_SupportingTables.xlsx

README.md
Sep 19, 2024 version files 3.09 MB

BriceEtAl2024_Data.xlsx

BriceEtAl2024_SupportingTables.xlsx

README.md
Abstract
The removal or addition of a predator in an ecosystem can trigger a trophic cascade, whereby the predator indirectly influences plants and/or abiotic processes via direct effects on its herbivore prey. A trophic cascade can operate through a densitymediated indirect effect (DMIE), where the predator reduces herbivore density via predation, and/or through a traitmediated indirect effect (TMIE), where the predator induces an herbivore trait response that modifies the herbivore’s effect on plants. Manipulative experiments suggest that TMIEs are an equivalent or more important driver of trophic cascades than are DMIEs. Whether this applies generally in nature is uncertain because few studies have directly compared the magnitudes of trait and densitymediated indirect effects on natural unmanipulated field patterns. A TMIE is often invoked to explain the textbook trophic cascade involving wolves (Canis lupus), elk (Cervus canadensis), and aspen (Populus tremuloides) in northern Yellowstone National Park. This hypothesis posits that wolves indirectly increase recruitment of young aspen into the overstory primarily through reduced elk browsing in response to spatial variation in wolf predation risk rather than through reduced elk population density. To test this hypothesis, we compared the effects of spatiotemporal variation in wolf predation risk and temporal variation in elk population density on unmanipulated patterns of browsing and recruitment of young aspen across 113 aspen stands over a 21year period (19992019) in northern Yellowstone National Park. Only two of ten indices of wolf predation risk had statistically meaningful effects on browsing and recruitment of young aspen, and these effects were 820 times weaker than the effect of elk density. To the extent that temporal variation in elk density was attributable to wolf predation, our results suggest that the wolfelkaspen trophic cascade was primarily densitymediated rather than traitmediated. This aligns with the alternative hypothesis that wolves and other actively hunting predators with broad habitat domains cause DMIEs to dominate whenever prey, such as elk, also have a broad habitat domain. For at least this type of predatorprey community, our study suggests that riskinduced trait responses can be abstracted or ignored while still achieving an accurate understanding of trophic cascades.
README: Data from: The primacy of densitymediated indirect effects in a community of wolves, elk, and aspen
https://doi.org/10.5061/dryad.2bvq83c0d
The file "BriceEtAl2024_Data.xlsx" includes all data needed to replicate models of aspen browse and recruitment probability as part of a test on density and traitmediated trophic cascades in Yellowstone National Park. The file "BriceEtAl2024_SupportingTables.xlsx" contains six supplemental tables from our analysis.
File Structure
BriceEtAl2024_Data.xlsx
This file has data on 26,012 individual young aspen collected between 19992019. Data were not collected in 2000 or 2015. The file has 3 tabs: one for data on aspen browse, and two for data on aspen recruitment. For the browsing data, each row represents an individual aspen at a specific site (i.e., "Plot", Column A), and the columns have data on whether that individual aspen was browsed by elk, and what it's height was. There are 10 columns for 10 separate measures of wolf predation risk of elk (described below and in methods), a column for annual elk abundance and annual elk density (abundance divided by 995 km2, i.e., study area), and columns for summer precipitation, winter snow water equivalent (SWE), and whether the site was within the 1988 Yellowstone fire boundary. Finally, there are three additional columns for height that were used in modeling: the first is a version of height that was scaled to have a mean of 0 and standard deviation of 1, and the second two are values of scaled height that are needed to model height as a spline.
In the tabs for aspen recruitment, each row represents a single site in a single year, with the total number of aspen individuals that were counted at that site that year, and how many of them have "recruited" into the overstory, meaning they have a height taller than 120cm. These tabs also have data for snow water equivalent, annual elk density, precipitation, and a single measure of wolf predation risk of elk (described below).
Below, we describe each column in the dataframes. Note that "winter" refers to Nov 1st – Apr 30th.
Browse
All data needed to reproduce models of browse probability in Brice et al. 2024. Columns are as follows:
 Plot (Column A): ID number for our 113 aspen plots in northern Yellowstone National Park.
 Year (Column B): the year that the data was collected.
 Browse (Column C): binary variable for whether an individual aspen was browsed by elk the previous winter. A value of 0 indicates that the aspen was not browsed, and a value of 1 indicates that the aspen was browsed.
 Height (Column D): height (in centimeters) of the tallest stem (i.e., leader stem) of individual aspen.
 WolfDensityAggregate (Column E): longterm average (i.e., 19992019) winter wolf density (wolves/30m2) estimated for each site. For example, the value of 0.0000903 for plot 1 in 1999 means that, within a 30m2 area, there was an estimated average of 0.0000903 wolves across all years of the study at plot 1.
 WolfDensityAnnual (Column F): winter wolf density (wolves/30m2) estimated at each site for the corresponding year (Col B). For example, the value of 0.0000429 for plot 1 in 1999 means that, within a 30m2 area, there was an estimated 0.0000429 wolves in 1999 at plot 1.
 AllElkKillsAggregate (Column G): longterm average (i.e., 19992019) density (kills/30m2) of all elk killed by wolves in winter at each site. For example, the value of 2.22E09 for plot 1 in 1999 means that, within a 30m2 area, there was an estimated average of 2.22E09 elk killed by wolves across all years of the study at plot 1.
 AllElkKillsAnnual (Column H): density (kills/30m2) of all elk killed by wolves in winter at each site for the corresponding year (Col B). For example, the value of 2.67E09 for plot 1 in 1999 means that, within a 30m2 area, there was an estimated of 2.67E09 elk killed by wolves in 1999 at plot 1.
 MaleElkKillsAggregate (Column I): longterm average (i.e., 19992019) density (kills/30m2) of adult and yearling male elk killed by wolves in winter at each site. This column can be interpreted the same way as Column G, but for just male elk instead of all elk.
 MaleElkKillsAnnual (Column J): density (kills/30m2) of adult and yearling male elk killed by wolves in winter at each site for the corresponding year (Col B).This column can be interpreted the same way as Column H, but for just male elk instead of all elk.
 FemaleElkKillsAggregate (Column K): longterm average (i.e., 19992019) density (kills/30m2) of adult female elk and elk calves killed by wolves in winter at each site. This column can be interpreted the same way as Column G, but for just female elk and elk calves instead of all elk.
 FemaleElkKillsAnnual (Column L): density (kills/30m2) of adult female elk and elk calves killed by wolves in winter at each site for the corresponding year (Col B). This column can be interpreted the same way as Column H, but for just female elk and elk calves instead of all elk.
 Openness (Column M): proportion of aspen plot that was not tree cover, estimated for 30m2 cells. A value of 0 indicates the landscape is completely covered by trees (i.e., all trees, 0% open), and a value of 1 indicates the landscape is completely open (i.e., no trees, 100% open)
 Smoothness (Column N): how smooth/flat the topography is. This value is a relative scale, with no units. A value of 0 indicates very rough terrain, and a value of 1 indicates very smooth/flat terrain.
 ElkAnnualCount (Column O): annual counts of the elk population for each year.
 Elk_Density (Column P): annual counts of the elk population divided by the study area to approximate elk density (elk/km2). For example, a value of 13.49 in 1999 at all stands indicates that there were an estimated 13.49 elk per km2 in 1999
 Precipitation_cm (Column Q): sum of daily precipitation (cm; sum of all forms converted to waterequivalent) from Apr 1st – July 31st at each stand each year. For example, a value of 17.04 at plot 1 in 1999 indicates that there was an estimated 17.04 cm of rain at plot 1 during the summer of 1999.
 SWE_tons (Column R): total winter snow water equivalent (tons/m2) at each stand each year. For example, a value of 14.49 at plot 1 in 1999 indicates that there was an estimated 14.49 tons of snow per m2 that accumulated during the winter before aspen was measured in 1999 at plot 1.
 Fire (Column S): a binary variable indicating whether a plot was within the perimeter of the 1988 Yellowstone fires. A value of 1 indicates that the plot was within the fire perimeter, and a value of 0 indicates the site was outside of the fire perimeter.
 HeightScaled (Column T): aspen heights (Col D) scaled to have a mean of 0 and standard deviation of 1. This column was used to create the values needed for the height spline (Cols U & V)
 HeightSpline1 (Column U): values for the first height spline, with a knot at 120cm height. Aspen stems that are greater than or equal to 120cm will have a value of 0.275. This column was used (along with Col V) to model aspen height as a nonlinear predictor of aspen browse.
 HeightSpline2 (Column V): values for the second height spline, with a knot at 120cm height. Aspen stems greater than 120cm will have a value of HeightScaled  0.275. Aspen stems that are less than 120cm will have a value of 0.
Recruitment
There are two tabs of data for the recruitment models: (1) Annual wolf density and (2) Annual male elk kills. The recruitment data was created by counting the number of stems at each site each year that were taller than or equal to the previously estimated peak browse threshold, which was 120cm in the annual wolf density model, and 123cm in the annual male elk kill model. Note that the two tabs have different numbers of rows; this difference is because outliers in the risk variable were removed prior to estimating the number of recruited stems (i.e., some plotyear combinations missing from each tab).
The columns for both tabs are all the same except for Column H, which is the risk variable. Descriptions of each column are as follows:
 Plot (Column A): ID number for our 113 aspen plots in northern Yellowstone National Park.
 Year (Column B): the year that the data was collected.
 N_stems (Column C): total number of aspen stems counted at each site each year (regardless of height). For example, at plot 1, there were 37 stems measured in 1999 (row 2, col C), and 42 stems measured in 2001 (row 3, col C).
 N_Recruits (Column D): total number of aspen stems taller than the peak browse threshold (i.e., 120cm for annual wolf density model, 123cm for annual male elk kill density model) at each site each year. For example, at plot 1, there were 0 stems taller than 120cm in 1999 (row 2, col C, "Annual wolf density" tab), and 15 stems taller than 120cm in 2010 (row 9, col C).
 SWE_tons (Column E): total winter (Nov 1st – Apr 30th) snow water equivalent (tons/m2) at each stand each year.
 ElkDensity (Column F): annual counts of the elk population divided by the study area to approximate elk density (elk/km2).
 Precip_cm (Column G): sum of daily precipitation (cm; sum of all forms converted to waterequivalent) from Apr 1st – July 31st at each stand each year
 Risk variable (Column H):
 WolfDensityAnnual: wolf density (wolves/30m2) estimated at each site for the corresponding year
 MaleElkKillsAnnual: density (kills/30m2) of adult and yearling male elk killed by wolves in winter at each site for the corresponding year
BriceEtAl2024_SupportingTables.xlsx
This excel file contains six tabs of supplementary results, each with it's own table. Each tab contains a text box with a description of the data. The tabs are as follows:
 Table S1_Random Effects: when modeling aspen browse probability, we used generalized linear mixed models that contained random effects for plot and year. We also tested various combinations of random slopes to determine what random effect structure was the best fit for the model. This table shows the results from the model selection process.
 Model (Column A): explains the random effect structure being tested. We tested various structures for each predation risk metric
 N (Column B): Sample size for each model
 Constant (Column C): the number of constants (i.e., intercepts) in the model. There was one constant in each model.
 Rand. Intercept (Column D): the number of random intercepts included in the random effect structure. Two random intercepts (one for plot, one for year) were included in each model.
 Rand. Slopes (Column E): the number of random slopes (and associated covariances) being tested in each model. A value of 0 indicates that there were no random slopes included in the model. In a model where we test the effect of including a random slope for elk density that varies by plot, we would have a random slope value of 2. This number includes the variance of the random slope for elk, as well as the covariance between elk and plot. In a model where we we test the effect of including a random slope for elk density that varies by plot AND a random slope for predation risk that varies by plot, we would have a random slope value of 5. This number includes the variance of the random slope for elk density, the variance of the random slope of risk, the covariance between elk density and plot, the covariance between risk and plot, and the covariance between elk density and risk.
 Covariate (Column F): the number of covariates, or fixed effects, in the model. Each model included 5 fixed effects.
 K (Column G): the total number of parameters in the model, which is the sum of columns CF.
 LL (Column H): the log likelihood of the model.
 AICc (Column I): the AICc value for the model. The model with the lowest value is the best fit for the data.
 ∆AICc (Column J): the difference in AICc between each model and the best model in the set. A value of 0 indicates that model is the best fit for the data.
 Wi (Column K): the Akaike weight of the model. The Akaike weight of model i can be interpreted as the probability that model i is the best model given the collection of models considered.
 TableS2_Height Spline: when modeling aspen browse probability, we included aspen height as a predictor. However, the effect of height on browse is nonlinear, with browse increasing with height up to a certain point, and then decreasing with height. To determine at which height browse switches from increasing to decreasing, we modeled height as a piecewise linear spline. This table shows the model selection results from this process. Each row is for a model with a different height threshold being tested. Cells highlighted in green are all within 2 AICc of the best model, which is indicated with bold font. The columns are similar to Table S1:
 Knot (Column A): indicates what height was used as the point at which browsing switches from increasing with height to decreasing.
 N (Column B): Sample size for each model.
 Constant (Column C): the number of constants (i.e., intercepts) in the model. There was one constant in each model.
 Rand. Intercept (Column D): the number of random intercepts included in the random effect structure. Two random intercepts (one for plot, one for year) were included in each model.
 Rand. Slope (Column E): the number of random slopes (and associated covariances) being tested in each model. See item 1.5 above for further detail.
 Covariate (Column F): the number of covariates, or fixed effects, in the model. Each model included 5 fixed effects.
 K (Column G): the total number of parameters in the model, which is the sum of columns CF.
 LL (Column H): the log likelihood of the model
 AICc (Column I): the AICc value for the model. The model with the lowest value is the best fit for the data
 ∆AICc (Column J): the difference in AICc between each model and the best model in the set. A value of 0 indicates that model is the best fit for the data.
 Wi (Column K): the Akaike weight of the model.
 RE Structure (Column L): the random effect structure included in the model.
 TableS3_Risk Shape: when modeling aspen browse probability, we wanted to know whether the effect of predation risk on browsing was linear, or if there was a threshold of risk, wherein there was little effect of risk up to a point, after which browsing decreased with risk. To test for a threshold, we compared models with a linear effect of predation risk to models with thresholds at low, moderate, or high risk, which were represented with piecewise linear splines. This table shows the model selection results from this process. Each row is for a model with a different risk threshold being tested. The columns are similar to Table S1:
 Model (Column A): indicates the functional form of risk (i.e., linear, spline) being tested.
 N (Column B): Sample size for each model
 Covariate (Column C): the number of covariates, or fixed effects, in the model. Models with a linear form of risk have 5 fixed effects, while models testing for a threshold have two covariates for risk, and therefore have 6 fixed effects.
 K (Column D): the total number of parameters in the model.
 LL (Column E): the log likelihood of the model.
 AICc (Column F): the AICc value for the model. The model with the lowest value is the best fit for the data.
 ∆AICc (Column G): the difference in AICc between each model and the best model in the set.
 Wi (Column H): the Akaike weight of the model.
 TableS4_Elk Shape: when modeling aspen browse probability, we wanted to know whether the effect of elk density on browsing was linear, or if there was a threshold of elk density, below which there was a constant effect of elk density. To test for a threshold, we compared models with a linear effect of elk density to models with a log effect of density, and thresholds at 26 elk/km2, modeled with piecewise linear splines. This table shows the model selection results from this process. Each row is for a model with a different elk density threshold being tested. The columns are similar to Table S1:
 Model (Column A): indicates the functional form of elk density (i.e., linear, log, spline) being tested.
 N (Column B): Sample size for each model.
 Rand. Intercept (Column C): the number of random intercepts included in the random effect structure. Two random intercepts (one for plot, one for year) were included in each model.
 Rand. Slope (Column D): the number of random slopes (and associated covariances) being tested in each model.
 Covariate (Column E): the number of covariates, or fixed effects, in the model. Models with a linear form of elk density have 5 fixed effects, while models testing for a threshold have two covariates for elk density, and therefore have 6 fixed effects.
 K (Column F): the total number of parameters in the model, which is the sum of columns CF.
 LL (Column G): the log likelihood of the model.
 AICc (Column H): the AICc value for the model. The model with the lowest value is the best fit for the data.
 ∆AICc (Column I): the difference in AICc between each model and the best model in the set.
 Wi (Column J): the Akaike weight of the model.
 RE Structure (Column K): the random effect structure included in the model.
 TableS5_Browse Models: coefficient estimates from the final models of browse probability as a function of height, SWE, elk density, and predation risk. Height was represented as a piecewise linear spline, with separate effects for stems shorter than ~120cm than those taller than ~120cm. The heights used for the spline are specific to each risk variable (see Table S4).
 Risk metric (Column A): the predation risk metric included in the model.
 Parameter (Column B): the parameters (i.e., predictors) included in the model.
 β (Column C): the beta coefficient estimated for each parameter in the model.
 95% Confidence Interval (Columns D & E): the lower (col D) and upper (col E) limits of the 95% confidence interval for the parameter estimate.
 SE (Column F): the standard error estimated for the parameter.
 P (Column G): the pvalue associated with the parameter.
 RE Structure (Column H): the random effects structure included in the model.
 TableS6_DensityRisk Interaction: coefficient estimates from models of browse probability as a function of height, SWE, elk density, predation risk, and an interaction between elk density and risk. Height was represented as a piecewise linear spline, with separate effects for stems shorter than ~120cm than those taller than ~120cm. The heights used for the spline are specific to each risk variable (see Table S4).
 Risk metric (Column A): the predation risk metric included in the model.
 Parameter (Column B): the parameters (i.e., predictors) included in the model.
 β (Column C): the beta coefficient estimated for each parameter in the model.
 95% Confidence Interval (Columns D & E): the lower (col D) and upper (col E) limits of the 95% confidence interval for the parameter estimate.
 SE (Column F): the standard error estimated for the parameter.
 P (Column G): the pvalue associated with the parameter.
 RE Structure (Column H): the random effects structure included in the model.