# Prey and predator density-dependent interactions under different water volumes

Ross N Cuthbert, Tatenda Dalu, Ryan J Wasserman, Arnaud Sentis, Olaf L F Weyl, P. William

William Froneman, Amanda Callaghan, Jaimie T A Dick

► To cite this version:

Ross N Cuthbert, Tatenda Dalu, Ryan J Wasserman, Arnaud Sentis, Olaf L F Weyl, et al.. Prey and predator density-dependent interactions under different water volumes. Ecology and Evolution, 2021, 11 (11), pp.1-10. 10.1002/ece3.7503. hal-03321109

HAL Id: hal-03321109

https://hal.inrae.fr/hal-03321109

Submitted on 17 Aug 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Prey and predator density-dependent interactions under different water volumes

Ross N. Cuthbert<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> ID | Tatenda Dalu<fnref n="3" />,<fnref n="4" /> ID | Ryan J. Wasserman<fnref n="3" />,<fnref n="5" /> ID | Arnaud Sentis<fnref n="6" /> ID |
Olaf L. F. Weyl<fnref n="7" /> | P. William Froneman<fnref n="5" /> | Amanda Callaghan<fnref n="8" /> | Jaimie T. A. Dick<fnref n="2" />

**Correspondence**

Ross N. Cuthbert, GEOMAR, Helmholtz-Zentrum für Ozeanforschung Kiel, 24105 Kiel, Germany.

Email: rossnoalcuthbert@gmail.com

**Funding information**

Department for the Economy, Northern Ireland; Alexander von Humboldt Foundation; National Research Foundation, Grant/Award Number: 117700 and 110507

**Abstract**

Predation is a critical ecological process that directly and indirectly mediates population stabilities, as well as ecosystem structure and function. The strength of interactions between predators and prey may be mediated by multiple density dependences concerning numbers of predators and prey. In temporary wetland ecosystems in particular, fluctuating water volumes may alter predation rates through differing search space and prey encounter rates. Using a functional response approach, we examined the influence of predator and prey densities on interaction strengths of the temporary pond specialist copepod *Lovenula raynerae* preying on cladoceran prey, *Daphnia pulex*, under contrasting water volumes. Further, using a population dynamic modeling approach, we quantified multiple predator effects across differences in prey density and water volume. Predators exhibited type II functional responses under both water volumes, with significant antagonistic multiple predator effects (i.e., antagonisms) exhibited overall. The strengths of antagonistic interactions were, however, enhanced under reduced water volumes and at intermediate prey densities. These findings indicate important biotic and abiotic contexts that mediate predator-prey dynamics, whereby multiple predator effects are contingent on both prey density and search area characteristics. In particular, reduced search areas (i.e., water volumes) under intermediate prey densities could enhance antagonisms by heightening predator-predator interference effects.

**KEYWORDS**

antagonism, consumer-resource, functional response, multiple predator effects, temporary pond, zooplankton

## 1 | INTRODUCTION

Interspecific interactions mediate resource population stabilities, with implications for ecosystem structure and function (Brooks & Dodson, 1965; Paine, 1980; Wasserman & Froneman, 2013). Interaction strengths between consumers and resources, such as

predators and prey, may be influenced by numerous biotic and abiotic environmental contexts that can challenge predictions of interaction strengths due to emergent effects (Cuthbert, Wasserman, et al., 2020; Sentis et al., 2019; Wasserman, Alexander, Dalu, et al., 2016), such as multiple predator effects (Barrios-O'Neill et al., 2014; Sih et al., 1998). Predator effects can be additive (i.e., predator-predator independence), whereby consumption rates are

<footnote n="1">¹GEOMAR Helmholtz-Zentrum für Ozeanforschung Kiel, Kiel, Germany</footnote>

<footnote n="2">²Institute for Global Food Security, School of Biological Sciences, Queen's University Belfast, Belfast, UK</footnote>

<footnote n="3">³South African Institute for Aquatic Biodiversity, Makhanda, South Africa</footnote>

<footnote n="4">⁴School of Biology and Environmental Sciences, University of Mpumalanga, Nelspruit, South Africa</footnote>

<footnote n="5">⁵Department of Zoology and Entomology, Rhodes University, Makhanda, South Africa</footnote>

<footnote n="6">⁶INRAE, Aix Marseille University, UMR RECOVER, Aix-en-Provence, France</footnote>

<footnote n="7">⁷DSI/NRF Research Chair in Inland Fisheries and Freshwater Ecology, South African Institute for Aquatic Biodiversity, Makhanda, South Africa</footnote>

<footnote n="8">⁸Ecology and Evolutionary Biology, School of Biological Sciences, University of Reading, Reading, UK</footnote>

<footnote>This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.</footnote>

<footnote>© 2021 The Authors. *Ecology and Evolution* published by John Wiley & Sons Ltd.</footnote>
predictable based on feeding rates of individuals (Cuthbert et al., 2020). In contrast, multiple predator effects can manifest as antagonistic (e.g., via predator-predator interference) or synergistic (e.g., via predator-predator cooperation) outcomes, which lessen and heighten prey risk, respectively (Sih et al., 1998; Vance-Chalcraft & Soluk, 2005; Wasserman, Alexander, Dalu, et al., 2016). Evidence for the importance of trait-mediated interactions such as multiple predator effects in influencing population dynamics is increasing (Anton et al., 2020; Gabowski, 2004; Schmitz et al., 2004; Trussell et al., 2004). In particular, trait-mediated effects in aquatic ecosystems have been found to be pervasive due to the presence of water-borne cues that illicit responses over entire populations (Alexander et al., 2013; Peacor & Werner, 2001; Trussell et al., 2004, 2008).

Multiple predator effects have been classically quantified by measuring predatory interactions toward one or just a few prey densities, using experimental designs that considered interspecific but not intraspecific consumers (Soluk, 1993; but see Skalski & Gilliam, 2001; Rall et al., 2008). Accordingly, there has been a relative lack of examination of potential conspecific multiple predator effects, despite the importance of intraspecific interactions in certain habitat types (Cuthbert, Dalu, Wasserman, Weyl, et al., 2020). Indeed, studies have found intraspecific interference to be common across many study systems, and these tend to be intermediate in magnitude overall (-0.6 to -0.7 m), despite high variability among studies (DeLong & Vasseur, 2011).

In addition to predator density effects, predator-prey interactions are known to be inherently prey density-dependent (Dick et al., 2014; Holling, 1959). Differences in prey density dependence forms can be driven by prey refuge effects that can arise at low-density prey populations (Alexander et al., 2012). For example, predators can target prey even when relatively rare in environments, promoting their extirpation. Conversely, predators can avoid rare prey and instead predate abundant groups, helping to stabilize populations and minimize extirpations through frequency-dependent predation (Murdoch, 1969; Murdoch et al., 1975). Classically, the "functional response" has been used to quantify rates of resource consumption as a function of resource density (Holling, 1959). Functional responses are characterized into three forms: types I (i.e., linear), II (i.e., hyperbolic), and III (i.e., sigmoid) (Hassell, 1978). Type I functional responses are mechanistically restricted to filter feeders owing to an absence of resource search-related constraints (Jeschke et al., 2004). Type II functional responses can be stabilizing compared with type I due to consumption limitations at high-resource densities. Conversely, type II functional responses can be destabilizing compared with type III owing to high consumption rates at low prey densities, while consumption rates are lower in type III functional responses at low prey densities (Dick et al., 2014; Hassell, 1978). The type III response also prevents prey from reaching high densities where population fluctuations are more likely to occur. Furthermore, recent advances have applied functional responses for the discernment of multiple predator effects across a range of resource densities, with several methods applied (e.g., multiplicative risk model: Wasserman, Alexander, Dalu, et al., 2016; population dynamic

model: Sentis & Boukal, 2018). Indeed, functional responses are known to be both predator and prey density-dependent (Abrams & Ginzburg, 2000; Coblentz & DeLong, 2020). Fundamentally, examinations of multiple predator effects using functional responses enable quantification of effects of both predator and prey density dependences simultaneously, as well as potential interaction effects among these factors. Previous works concerning multiple predator effects have shown significant dependence on prey availability, as well as other factors such as climatic warming and body size (Dalal et al., 2020; Sentis et al., 2016).

Functional responses have been recently applied to quantify interaction strengths within temporary pond ecosystems in arid environments (Buxton et al., 2020; Cuthbert et al., 2018; Wasserman, Alexander, Barrios-O'Neill, et al., 2016). Trophic structuring in these systems is atypical and can be determined by phenologies of dormant eggs, which hatch when ponds fill with water and the hydroperiod begins (Wasserman, Alexander, Barrios-O'Neill, et al., 2016). Communities in these systems are highly biodiverse (Bird et al., 2019) and increasingly threatened by anthropogenic global change (Dalu, Wasserman, & Dalu, 2017). In austral environments, temporary ponds are poorly studied despite their suitability as model systems in ecology for testing theories (De Meester et al., 2005). Predation pressure in these systems is transient throughout the hydroperiod, with higher predators such as hexapods arriving late in the hydroperiod via aerial dispersal (Wasserman et al., 2018). Accordingly, for much of the early hydroperiod stages, or throughout during short wet phases, predator assemblages are dominated by abundant zooplankton species owing to rapid mass hatching events following inundation. Interspecific interactions in these systems among conspecific predatory zooplankton can thus be marked, as such consumers are under temporal pressures to utilize resources to enable rapid development and maturation to reproduce prior to the hydroperiod ending (i.e., pond drying) (Dalu, Wasserman, Vink, et al., 2017). However, despite high abundances of zooplankton in these systems and spatiotemporally changeable population dynamics, few works have examined the potential for conspecific multiple predator effects as predator and prey densities simultaneously shift (Cuthbert, Dalu, Wasserman, Weyl, et al., 2020).

Naturally, water volumes are extremely variable in temporary pond ecosystems owing to wetting and drying, which can mediate interaction strengths owing to differences in search areas and encounter rates with prey or other predators (Uiterwaal et al., 2018; Uiterwaal & DeLong, 2018). Moreover, anthropogenic processes such as water extraction could also alter the volume of ponds, particularly as water shortages become more severe (Dalu, Wasserman, & Dalu, 2017). Reductions in search area might heighten encounter rates with prey by improving search efficiencies, and yet could also drive greater interference or mutualistic effects between interacting predators (i.e., multiple predator effects). The present study thus examined the influence of water volume on emergent multiple predator effects between conspecifics of the calanoid copepod *Lovenula raynerae* Suárez-Morales et al. (2015), preying on representative cladoceran prey. This copepod species has been
identified as an important and abundant predator in temporary pond ecosystems, where cladocerans such as daphniids coexist with predatory copepods due to simultaneous resting egg hatching (Cuthbert et al., 2018; Wasserman, Alexander, Barrios-O'Neill, et al., 2016; Wasserman et al., 2018). We systematically altered predator abundances using an additive experimental design, while also varying water volumes and prey densities factorially. We then used a comparative functional response approach, combined with a population dynamic model, to predict intraspecific multiple predator effects as population densities shift.

## 2 | MATERIALS AND METHODS

Adult male *L. raynerae* (4.5–5.0 mm) (Suárez-Morales et al., 2015) and mature *Daphnia pulex* complex (2.2–2.6 mm) were collected from a temporary pond close to Makhanda in the Eastern Cape Province of South Africa (33°15'04.1"S 26°26'17.0"E) by hauling a zooplankton net through the water column. Zooplankton were transported in source water to the Department of Zoology and Entomology, Rhodes University, separated, and housed in a controlled environment room (21 ± 1°C; 12:12 light:dark phase) in 5-L tanks containing filtered (20 µm sieved) source water. Copepods were unfed for 24 hr before use to standardize levels of hunger.

Feeding rates of *L. raynerae* were quantified toward daphniid prey under two levels of water volume and four predator densities. Prey were introduced at one of six densities (2, 4, 8, 16, 32, and 64) into 100-ml experimental arenas (5.6 cm diameter) containing 40 or 80 ml of filtered source water. After settling for 2 hr, one of four predator densities (1, 2, 3, or 4 ind.) was introduced into experimental arenas, with predators allowed to feed for 18 hr. Five replicates were undertaken per experimental group (i.e., 240 units overall; 6 prey densities × 2 volumes × 4 predator densities × 5 replicates), and the design was fully randomized to avoid positional effects. Five replicates of predator-free controls were also run at each prey density and water volume to quantify background prey mortality rates. After the allocated feeding period, predators were removed and remaining live prey counted to quantify numbers killed. Predators were not used in more than one experimental trial to standardize prior predator experience and avoid the potentially confounding factor of predator identity.

Differences in the proportions of prey eaten were analyzed using quasi-binomial generalized linear models (owing to residual overdispersion) as a function of water volume, predator density and prey density, and their interactions. Nonsignificant terms were removed from the model such that the final model contained only significant predictors (Crawley, 2007). Post hoc pairwise tests were computed using Tukey's comparisons.

Binomial generalized linear models were additionally used to categorize functional response types for each depth treatment at the single-predator density (i.e., predator density = 1 predator; Juliano, 2001; Pritchard et al., 2017). A type II functional response was indicated through the presence of a significantly negative linear

coefficient in response to increasing prey density. Given that prey were not replaced following consumption over the course of the experiment, Rogers' random predator equation was used to model functional responses at the single-predator densities (Rogers, 1972):

$$N_e = N_0 (1 - \exp(a(N_e h - T))) \quad (1)$$

where $N_e$ is the number of prey eaten, $N_0$ is the initial density of prey, $a$ is the attack constant, $h$ is the handling time, and $T$ is the total experimental period. The Lambert W function was used to fit the model to the data (Bolker, 2008; Pritchard et al., 2017). The random predator equation is robust to prey depletion in parameter estimation (Cuthbert, Wasserman, et al., 2020). Total prey depletion occurred across all replicates in the following two treatments: at 1 predator, 40 ml volume, 2 prey; and at 3 predators, 80 ml volume, 2 prey. Attack rates and handling times were compared using the differences (delta) method outlined in Juliano (2001), pairwise between each water volume.

We quantified interaction strength (IS) as the proportion of prey killed at each predator density, water volume, and prey density by dividing the number of prey consumed by the initial prey density (Veselý et al., 2019):

$$IS(P, Z) = \frac{N_p - N_{P,Z}}{N_p} \quad (2)$$

where $N_p$ and $N_{P,Z}$ are the numbers of live prey at the beginning and end of the experiment, respectively. The proportion of prey killed (IS) includes both trophic interactions (i.e., feeding on prey) and nontrophic interactions that can enhance (e.g., facilitations among predators) or reduce (e.g., interference among predators) trophic interactions. To disentangle trophic (IS$_T$) and nontrophic (IS$_{NT}$) interactions, we next used a population dynamic approach that quantifies IS$_{NT}$ as the difference between the observed IS and the predicted IS from single-predator functional responses (Sentis et al., 2017). Accordingly, we used our attack rate and handling time estimates from single-predator functional responses to predict multiple predator feeding rates, which were then compared with observed multiple predator feeding rates. This was done separately using functional response parameters for each water volume using the corresponding single-predator functional response parameters. Estimations of IS$_T$ were calculated following McCoy et al. (2012) and Sentis and Boukal (2018):

$$\frac{dN}{dt} = - \frac{aN}{1 + ahN} P \quad (3)$$

where N is the prey population density, P is the predator population density, and a and h are the attack rate and handling time obtained from the single-predator functional response estimates. This model assumes no emergent multiple predator effects, and its predictions can be compared with multiple predator feeding trials to assess the sign and strength of multiple predator effects. To generate predictions of expected prey survival in the multipredator experiments, initial values of N and P are set at the experimental initial prey and predator densities corresponding to the experimental treatment. For each
predator density, water volume, and prey density, Equation 3 was then integrated over the full experimental time to obtain expected numbers of surviving prey. To estimate the variance around the predictions, we used a global sensitivity analysis that uses the 95% confidence intervals of each functional response parameter estimate and their variance-covariance matrix (covariance is assumed to be zero when unknown) to generate 100 random parameter sets using a Latin hypercube sampling algorithm (Soetaert & Petzoldt, 2010). For each parameter set (n = 100), Equation 3 is then integrated over time and expected prey survival calculated using the “sensRange” function in the R package “FME” (Soetaert & Petzoldt, 2010). In turn, these numbers of surviving prey predicted were used reciprocally to determine the numbers of prey killed at each multiple predator density, water volume and prey density, and thus the proportion of prey killed per treatment for comparison with IS (Equation 2). Multiple predator effects, including interference effects (i.e., $IS_{NT}$), were then calculated by subtracting the mean $IS_{T}$ (i.e., predictions) from IS (i.e., observations). We note that variance in our prediction of $IS_{NT}$ at predator densities of 1 was due to differences (i.e., error) between the modeled functional responses and underlying consumption data (see Figure 1). Furthermore, the $IS_{T}$ predictions were made on a continuous scale, whereas IS data (proportion killed) are based on discrete counts of prey killed. This may lead to some divergence between predictions and observations, especially at low prey densities where continuous and discrete scales are more likely to differ and thus be more sensitive to less accurate model fitting. Positive and negative values of $IS_{NT}$ correspond to prey risk enhancement and reduction, respectively. Multiple predator effects were analyzed using linear models. Homogeneity of variances and residual normality were examined. A backward step deletion process was employed to obtain the most parsimonious model through removal of nonsignificant terms and interactions, as before (Crawley, 2007). We

compared models with the prey density covariate included as linear and quadratic terms via AICc, to select the model that minimized information loss. Tukey's tests were used post hoc for comparison of levels within significant effects.

## 3 | RESULTS

Most control prey survived across predator-free treatments (96% survival, ±8% SD), and thus, no adjustments were made to experimental mortality rates. Predation rates differed significantly according to predator density and prey density (Table 1); there was no significant effect of water depth on total consumption rates, nor any significant two- or three-way interaction terms. With the exception of 2 versus 3 predators ($p > 0.05$), total feeding rates always significantly increased with increasing predator densities (all $p < 0.05$). Predation rates fell significantly with increasing prey density (Table 1).

The proportion of prey consumed was negatively related to initial prey density, and thus, functional responses were categorized as type II at the single-predator density for both water volume treatments (Table 2; Figure 1). Attack rates tended to be higher, and handling times shorter, in the low- compared with high-volume treatments. However, neither attack rates nor handling times differed significantly pairwise according to volume treatments at each predator density (attack rate, $z = 0.73$, $p > 0.05$; handling time, $z = 0.41$, $p > 0.05$).

Multiple predator effects were predominantly antagonistic (i.e., negative $IS_{NT}$) among conspecific *L. raynerae* overall (Figure 2), indicating prey risk reductions in multiple predator groupings. Multiple predator effects differed significantly among predator

<figure><img src="image_2.png" /><figcaption>**FIGURE 1** Functional responses of *Lovenula raynerae* under single-predator densities between water volumes. Points are means ± 1 standard error, and smaller points are raw data. Note that the curves were modeled using data from all predators within a given treatment. The solid lines show the functional responses modeled from the random predator equation</figcaption></figure>
densities ($F_{3,233} = 5.16$, *p* < 0.01), whereby predator densities of 3 and 4 were significantly more antagonistic than densities of 1 overall (both *p* < 0.05). The strength of negative multiple predator effects among *L. raynerae* was also significantly greater under reduced water volumes overall ($F_{1, 233} = 5.32$, *p* < 0.05), whereby greater single-predator attack rates at reduced volumes (Table 2) caused higher predictions, and thus observed/predicted divergence, at low prey densities. Inclusion of prey density as a quadratic term improved the model as compared to a linear term alone ($\Delta AICc = 3.88$; linear, *t* = 0.12, *p* > 0.05; quadratic, *t* = 2.43, *p* < 0.05), indicating that antagonistic multiple predator effects scaled unimodally at intermediate prey densities. The strength of antagonisms tended to be

greatest around prey densities of 8 (Figure 2). There were no significant interaction terms (all *p* > 0.05).

## 4 | DISCUSSION

The present study found density and water volume-dependent multiple predator effects to mediate predator-prey interactions between temporary pond specialists. While feeding rates were found to increase overall under increasing predator densities, negative multiple predator effects pointed to increasing interference effects among conspecific predators, and thus reduced per capita consumption. These negative effects were, however, most prevalent under reduced water volumes, indicating that smaller search areas could exacerbate antagonistic interactions that alleviate prey risk. In ephemeral ponds that are subject to changeable water volumes, this could result in changeable multiple predator effects as ponds dry. Furthermore, while feeding rates at low prey densities were not suppressed (i.e., type II functional responses), the magnitude of multiple predator effects peaked at intermediate prey densities, suggesting that the availability of prey in environments can further influence predator-predator interactions. This finding corroborates previous studies that have found multiple predator interactions to be prey density-dependent (Sentis et al., 2016). Centrally, as predators seldom occur in individually in ecosystems, understanding multiple predator effects is critical for predicting interaction strengths holistically (Cuthbert, Dalu, Wasserman, Weyl, et al., 2020; Griffin et al., 2013; Sih et al., 1998; Wasserman, Alexander, Dalu,

<table><caption>TABLE 1 Generalized linear model results considering feeding rates by multiple *Lovenula raynerae* as a function of predator density, water volume, and prey density. Significant predictors are emboldened</caption><thead><tr><th>Predictor</th><th>F-value (df)</th><th>p-value</th></tr></thead><tbody><tr><td>Predator density</td><td>18.79 (3, 236)</td><td>&lt;0.001</td></tr><tr><td>Volume</td><td>0.02 (1, 235)</td><td>0.90</td></tr><tr><td>Prey density</td><td>247.53 (1, 234)</td><td>&lt;0.001</td></tr><tr><td>Predator density × volume</td><td>0.50 (3, 231)</td><td>0.68</td></tr><tr><td>Predator density × prey density</td><td>1.40 (3, 228)</td><td>0.24</td></tr><tr><td>Volume × prey density</td><td>1.34 (1, 227)</td><td>0.24</td></tr><tr><td>Predator density × volume × prey density</td><td>0.03 (3, 224)</td><td>0.99</td></tr></tbody></table>

<table><caption>TABLE 2 Functional response linear coefficients (and *p*-values), attack rates, and handling times across predator density and water volume treatments alongside standard errors (SE)</caption><thead><tr><th>Predator density</th><th>Volume (ml)</th><th>Linear coefficient, p-value</th><th>Attack rate, SE</th><th>Handling time, SE</th></tr></thead><tbody><tr><td>1</td><td>40</td><td>-0.03, &lt;0.001</td><td>1.52, 0.93</td><td>0.25, 0.05</td></tr><tr><td>1</td><td>80</td><td>-0.03, &lt;0.001</td><td>0.79, 0.38</td><td>0.28, 0.06</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>FIGURE 2 Predicted nontrophic interaction strength (i.e., multiple predator effects) of (a) one, (b) two, (c) three, and (d) four *Lovenula raynerae* predators between water volumes and across prey densities. Negative values are indicative of prey risk reduction, while positive values indicate prey risk enhancement</figcaption></figure>
et al., 2016). For example, these quantifications could be used to better understand the implications of predator extirpation for trophic dynamics within ecosystems, or the effects of hydrological alterations following environmental change in aquatic systems. In tropical waters, increasing droughts may reduce water volumes in temporary ponds due to higher water demands and changing precipitation patterns, or could result in their increase via deepening through land modifications for water storage purposes (Dalu, Wasserman, & Dalu, 2017).

Studies concerning multiple predator effects have reported a lack of generality among consumer-resource systems, with examples of antagonistic, additive, and synergistic outcomes in recent decades (Barrios-O'Neill et al., 2014; Losey & Denno, 1998; Sih et al., 1998; Soluk, 1993; Vance-Chalcraft & Soluk, 2005; Wasserman, Alexander, Dalu, et al., 2016). The lack of consistency suggests that multiple predator effects are system-specific and may be mediated by context-dependent factors such as taxonomic grouping, population demographics, and habitat characteristics. However, the current sparsity of studies into multiple predator effects overall negates broader generalizations within or between such factors, particularly at the conspecific level. Recent works into another temporary pond specialist copepod, *Paradiaptomus lamellatus* Sars, 1985, has similarly found antagonistic interactions in calanoid copepods (Cuthbert, Dalu, Wasserman, Monaco, et al., 2020). This species is known to coexist with the focal species in the present study, *L. raynerae*, although it is less voracious (Wasserman, Alexander, Barrios-O'Neill, et al., 2016). Interspecific multiple predator effects among temporary pond copepods have, however, previously been found to combine additively (Cuthbert, Callaghan, et al., 2020). Moreover, Buxton et al. (2020) found multiple predator effects of a congeneric copepod to combine additively, if not synergistically, with higher-order predatory hexapods from temporary ponds. Accordingly, a variety of multiple predator effects have been reported in austral temporary pond ecosystems, with implications for prey population stability and "boom-bust" population dynamics (Wasserman et al., 2018). Despite these recent advances, the influence of predator densities beyond pairs, as well as effects in combination with water depth variations, had yet to be considered prior to this study. Nonetheless, the present study found the multiple predator effect to be strongest at predator densities beyond two, with relatively little further reduction in effects between densities of three and four. The increasing intensities of nontrophic interactions at reduced water volumes found in the present study likely emanate from greater predator-predator encounter rates, which interfere with predatory efficiencies, yet the mechanisms underlying this require further elucidation empirically at a broader range of volumes.

Prey density implications for multiple predator effects have been poorly studied until recently. Previous studies have identified intraguild predation (i.e., interactions among predators) to lessen under increasing prey densities (Lucas & Rosenheim, 2011; Sentis et al., 2013). In agreement with the finding of this study where multiple predator effects were highest at intermediate prey densities, Sentis et al. (2016) found unimodal scaling of multiple predator

effects. This patterning might arise from the interplays between predation rates and predator satiation. At low prey densities, most if not all prey are consumed, which masks any potential interference effects in nonreplacement experimental feeding designs. Conversely, at high prey densities, not all prey are consumed, and thus high predator satiation minimizes multiple predator interference with one another. The comparative functional response approach thus lends itself to studies of multiple predator effects, as they are able to capture effects across a range of resource densities, which were not considered in many previous studies (Griffen, 2006; Soluk, 1993; Tylianakis & Romo, 2010; Vance-Chalcraft & Soluk, 2005). This thus enhances the resolution of results across different population abundance scenarios.

In the present study, owing to the lack of interaction effects, this unimodal patterning where antagonisms peaked at intermediate prey levels was consistent irrespective of predator density or water volume. However, the latter also exacerbated multiple predator effects where water volumes were reduced. In the case of water volume, reductions in search area associated with reduced water levels likely increase encounter rates among predators, and particularly for actively pelagic consumers, such as *L. raynerae* in the present study. Copepods utilize hydromechanical cues, detected by mechanoreceptors on the antennules, for the detection of predators and prey (Hwang & Strickler, 2001). While considered the most important cue factor for copepod predator-prey interactions, and with evidence suggesting that copepods can differentiate between signals from predators and prey (Hwang & Strickler, 2001), conspecific hydromechanical cue recognition dynamics are unknown. Reduced proximities between conspecifics in association with increased individual numbers, or decreased water volumes, may result in increased hydromechanical signal reception. If these signals are indistinguishable from predator cues, crowding may result in increases in anti-predator responses, with implications for foraging time and efficacy at reduced water volumes. Alternatively, antagonistic effects might merely reflect an increase in predator-predator encounters and thus interference. These effects may be less pronounced for benthic ambush predators, which might exhibit lower predator-predator encounters. Indeed, Sentis et al. (2016) found multiple predator effects to differ markedly depending on the predator assemblage characteristics. In the context of functional responses, previous works have also found arena size to significantly alter the scaling of parameters (attack rate and handling time) (Uiterwaal et al., 2018; Uiterwaal & DeLong, 2018); thus, in a similar vein, we suggest that consideration for search area is critical in studies into multiple predator effects in aquatic systems.

Austral temporary ponds lack examination as they remain inconspicuous during dry periods. This is despite them being threatened by a number of anthropogenically mediated processes, including climate and land use changes (Dalu, Wasserman, & Dalu, 2017). In particular, increased droughts, pollutants and hydrological alternations *via* water extraction threaten these systems and the biodiversity they support (Mabidi et al., 2018). Many temporary ponds are not mapped let alone studied, and understandings of trophic dynamics
in these systems are lacking. Predatory copepods, such as *L. raynerae*, are often highly abundant in temporary ponds in the study region (Wasserman et al., 2018). These copepods hatch early in the hydroperiod from drought-resistant eggs, and thus have high potential to engage in predator-predator interactions with conspecifics. The focal prey in this study, *D. pulex*, also hatch from dormant eggs and cooccur with copepods, and thus, our predator-prey system is representative of multiple community composition scenarios considering predator and prey densities. This is particularly the case given that food webs in these systems can be highly simplified, particularly during the early stages of hydroperiod where recruitment is largely restricted to internal hatching events. Overall, the results from the present study advance our understandings of predator and prey density dependences of inter- and intraspecific interactions in these ecosystems, and also indicate how differences in water depth alter the nature of trophic dynamics. Reductions in water depth through the hydroperiod likely exacerbate predator-predator interactions, and thus may have implications for prey population stability. Future studies are required to examine other trait-mediated effects in these systems, as well as for considering more complex assemblages of trophic groups and additional environmental variables. Empirical work that considers a broader range of water volumes is also required to further our understandings of multiple predator effects in aquatic ecosystems.

## ACKNOWLEDGMENTS

This article is dedicated to Prof. Olaf L.F. Weyl, who passed away suddenly on the 14 November 2020. He provided immense support for work on ephemeral wetlands in Southern Africa, and other freshwater-waters. RNC is funded through a Postdoctoral Research Fellowship from the Alexander von Humboldt Foundation. TD acknowledges funding from the National Research Foundation (Grant No. 117700). This study received funding from the Department for the Economy, Northern Ireland. We thank Rhodes University for the provision of laboratory facilities. We acknowledge the use of infrastructure and equipment provided by the South African Institute for Aquatic Biodiversity (SAIAB) Research Platform and the funding channeled through the National Research Foundation—SAIAB Institutional Support system. This study was partly funded by the National Research Foundation—South African Research Chairs Initiative of the Department of Science and Innovation (Inland Fisheries and Freshwater Ecology, Grant No. 110507).

## CONFLICT OF INTEREST

There are no conflicts of interest relating to this study.

## AUTHOR CONTRIBUTION

**Ross N. Cuthbert:** Conceptualization (lead); Data curation (lead); Formal analysis (lead); Investigation (lead); Methodology (lead); Visualization (lead); Writing—original draft (lead). **Tatenda Dalu:** Conceptualization (equal); Writing—review & editing (equal). **Ryan J. Wasserman:** Conceptualization (equal); Writing—review & editing (equal). **Arnaud Sentis:** Formal analysis (equal); Software (equal);

Writing—review & editing (equal). **Olaf L. F. Weyl:** Project administration (equal); Resources (equal); Writing—review & editing (equal).

**P. William Froneman:** Resources (equal); Writing—review & editing (equal). **Amanda Callaghan:** Supervision (equal); Writing—review & editing (equal). **Jaimie T. A. Dick:** Supervision (equal); Writing—review & editing (equal).

## DATA AVAILABILITY STATEMENT

Underlying data are available on the Dryad digital repository (https://doi.org/10.5061/dryad.qz612jmf5).

## ORCID

Ross N. Cuthbert

<figure><img src="image_4.png" /></figure>

https://orcid.org/0000-0003-2770-254X

Tatenda Dalu

<figure><img src="image_5.png" /></figure>

https://orcid.org/0000-0002-9019-7702

Ryan J. Wasserman

<figure><img src="image_6.png" /></figure>

https://orcid.org/0000-0002-4162-1503

Arnaud Sentis

<figure><img src="image_7.png" /></figure>

https://orcid.org/0000-0003-4617-3620

## REFERENCES

Abrams, P. A., & Ginzburg, L. R. (2000). The nature of predation: Prey dependent, ratio dependent or neither? *Trends in Ecology and Evolution*, 15, 337–341. https://doi.org/10.1016/S0169-5347(00)01908-X

Alexander, M. E., Dick, J. T. A., & O'Connor, N. E. (2013). Trait-mediated indirect interactions in a marine intertidal system as quantified by functional responses. *Oikos*, 122, 1521–1531. https://doi.org/10.1111/j.1600-0706.2013.00472.x

Alexander, M. E., Dick, J. T. A., O'Connor, N. E., Haddaway, N. R., & Farnsworth, K. D. (2012). Functional responses of the intertidal amphipod *Echinogammarus marinus*: Effects of prey supply, model selection and habitat complexity. *Marine Ecology Progress Series*, 468, 191–202. https://doi.org/10.3354/meps09978

Anton, A., Geraldi, N. R., Ricciardi, A., & Dick, J. T. A. (2020). Global determinants of prey naiveté to exotic predators. *Proceedings of the Royal Society B: Biological Sciences*, 287, 20192978. https://doi.org/10.1098/rspb.2019.2978

Barrios-O'Neill, D., Dick, J. T. A., Emmerson, M. C., Ricciardi, A., MacIsaac, H. J., Alexander, M. E., & Bovy, H. C. (2014). Fortune favours the bold: A higher predator reduces the impact of a native but not an invasive intermediate higher predator reduces the impact of a native but not an invasive intermediate predator. *Journal of Animal Ecology*, 83, 693–701.

Bird, M. S., Mlambo, M. C., Wasserman, R. J., Dalu, T., Holland, A. J., Day, J. A., Villet, M. H., Bilton, D. T., Barber-James, H. M., & Brendonck, L. (2019). Deeper knowledge of shallow waters: Reviewing the invertebrate fauna of southern African temporary wetlands. *Hydrobiologia*, 827, 89–121. https://doi.org/10.1007/s10750-018-3772-z

Bolker, B. M. (2008). *emdbbook: Ecological models and data in R*. Princeton University Press.

Brooks, J. L., & Dodson, S. I. (1965). Predation, body size, and composition of plankton. *Science*, 150, 28–35. https://doi.org/10.1126/science.150.3692.28

Buxton, M., Cuthbert, R. N., Dalu, T., Nyamukondiwa, C., & Wasserman, R. J. (2020). Complementary impacts of heterospecific predators facilitate improved biological control of mosquito larvae. *Biological Control*, 144, 104216. https://doi.org/10.1016/j.biocontrol.2020.104216

Coblenzt, K. E., & DeLong, J. P. (2020). Predator-dependent functional responses alter the coexistence and indirect effects among prey that share a predator. *Oikos*, 129, 1404–1414.

Crawley, M. J. (2007). *The R book*. Wiley.

Cuthbert, R. N., Callaghan, A., Sentis, A., Dalal, A., & Dick, J. T. A. (2020). Additive multiple predator effects can reduce mosquito populations.
*Ecological Entomology*, 45, 243-450. https://doi.org/10.1111/
een.12791

Cuthbert, R. N., Dalu, T., Wasserman, R. J., Callaghan, A., Weyl, O. L. F., & Dick, J. T. A. (2018). Calanoid copepods: An overlooked tool in the biocontrol of disease vector mosquitoes. *Journal of Medical Entomology*, 55, 1656-1658.

Cuthbert, R. N., Dalu, T., Wasserman, R. J., Monaco, C. J., Callaghan, A., Weyl, O. L. F., & Dick, J. T. A. (2020). Assessing multiple predator, diurnal and search area effects on predatory impacts by ephemeral wetland specialist copepods. *Aquatic Ecology*, 54, 181-191. https://doi.org/10.1007/s10452-019-09735-y

Cuthbert, R. N., Dalu, T., Wasserman, R. J., Weyl, O. L. F., Froneman, P. W., Callaghan, A., & Dick, J. T. A. (2020). Examining intraspecific multiple predator effects across shifting predator sex ratios. *Basic and Applied Ecology*, 45, 12-21. https://doi.org/10.1016/j.baae.2020.03.002

Cuthbert, R. N., Wasserman, R. J., Dalu, T., Kaiser, H., Weyl, O. L. F., Dick, J. T. A., Sentis, A., McCoy, M. W., & Alexander, M. E. (2020). Influence of intra- and interspecific variations in predator-prey body size ratios on trophic interaction strengths. *Ecology and Evolution*, 10, 5946-5962.

Dalal, A., Cuthbert, R. N., Laverty, C., Dick, J. T. A., Sentis, A., Barrios-O'Neill, D., Ortiz-Perea, N., Callaghan, A., & Gupta, S. (2020). Prey size and predator density modify impacts by natural enemies towards mosquitoes. *Ecological Entomology*, 45, 423-433. https://doi.org/10.1111/een.12807

Dalu, T., Wasserman, R. J., & Dalu, M. T. B. (2017). Agricultural intensification and drought frequency increases may have landscape-level consequences for ephemeral ecosystems. *Global Change Biology*, 23, 983-985. https://doi.org/10.1111/gcb.13549

Dalu, T., Wasserman, R. J., Vink, T. J. F., & Weyl, O. L. F. (2017). Sex and species specific isotopic niche specialisation increases with trophic complexity: Evidence from an ephemeral pond ecosystem. *Scientific Reports*, 7, 43229. https://doi.org/10.1038/srep43229

De Meester, L., Declerck, S., Stoks, R., Louette, G., Van de Meutter, F., De Bie, T., Michels, E., & Brendonck, L. (2005). Ponds and pools as model systems in conservation biology, ecology and evolutionary biology. *Aquatic Conservation: Marine and Freshwater Ecosystems*, 15, 715-725. https://doi.org/10.1002/aqc.748

DeLong, J. P., & Vasseur, D. A. (2011). Mutual interference is common and mostly intermediate in magnitude. *BMC Ecology*, 11, 1. https://doi.org/10.1186/1472-6785-11-1

Dick, J. T. A., Alexander, M. E., Jeschke, J. M., Ricciardi, A., Maclsaac, H. J., Robinson, T. B., Kumschick, S., Weyl, O. L. F., Dunn, A. M., Hatcher, M. J., Paterson, R. A., Farnsworth, K. D., & Richardson, D. M. (2014). Advancing impact prediction and hypothesis testing in invasion ecology using a comparative functional response approach. *Biological Invasions*, 16, 735-753. https://doi.org/10.1007/s10530-013-0550-8

Gabowski, J. H. (2004). Habitat complexity disrupts predator - Prey interactions but not the trophic cascade on oyster reefs. *Ecology*, 85, 995-1004. https://doi.org/10.1890/03-0067

Griffen, B. (2006). Detecting emergent effects of multiple predator species. *Oecologia*, 148, 702-709. https://doi.org/10.1007/s00442-006-0414-3

Griffin, J. N., Byrnes, J. E., & Cardinale, B. J. (2013). Effects of predator richness on prey suppression: A meta-analysis. *Ecology*, 94, 2180-2187. https://doi.org/10.1890/13-0179.1

Hassell, M. P. (1978). *The dynamics of arthropod predator-prey systems*. Princeton University Press.

Holling, C. S. (1959). Some characteristics of simple types of predation and parasitism. *Canadian Entomologist*, 91, 385-398. https://doi.org/10.4039/Ent91385-7

Hwang, J. S., & Strickler, J. R. (2001). Can copepods differentiate prey from predator hydromechanically? *Zoological Studies*, 40, 1-6.

Jeschke, J. M., Kopp, M., & Tollrian, R. (2004). Consumer-food systems: Why type I functional responses are exclusive to filter feeders.

*Biological Reviews of the Cambridge Philosophical Society*, 79, 337-349. https://doi.org/10.1017/S1464793103006286

Juliano, S. A. (2001). Nonlinear curve fitting: Predation and functional response curves. In S. M. Scheiner, & J. Gurevitch (Eds.), *Design and analysis of ecological experiments* (pp. 178-196). Oxford University Press.

Losey, J. E., & Denno, R. F. (1998). Positive predator-predator interactions: Enhanced predation rates and synergistic suppression of aphid populations. *Ecology*, 79, 2143-2152.

Lucas, E., & Rosenheim, J. A. (2011). Influence of extraguild prey density on intraguild predation by heteropteran predators: A review of the evidence and a case study. *Biological Control*, 59, 61-67. https://doi.org/10.1016/j.biocontrol.2011.05.010

Mabidi, A., Bird, M. S., & Perissinotto, R. (2018). Increasing salinity drastically reduces hatching success of crustaceans from depression wetlands of the semi-arid Eastern Cape Karoo region, South Africa. *Scientific Reports*, 8, 5983. https://doi.org/10.1038/s41598-018-24137-0

McCoy, M. W., Stier, A. C., & Osenberg, C. W. (2012). Emergent effects of multiple predators on prey survival: The importance of depletion and the functional response. *Ecology Letters*, 15, 1449-1456. https://doi.org/10.1111/ele.12005

Murdoch, W. W. (1969). Switching in general predators: Experiments on predator specificity and stability of prey populations. *Ecological Monographs*, 39, 335-354. https://doi.org/10.2307/1942352

Murdoch, W. W., Avery, S., & Smyth, M. E. B. (1975). Switching in predatory fish. *Ecology*, 56, 1094-1105. https://doi.org/10.2307/1936149

Paine, R. T. (1980). Food webs: Linkage, interaction strength and community infrastructure. *Journal of Animal Ecology*, 49, 667-685. https://doi.org/10.2307/4220

Peacor, S. D., & Werner, E. E. (2001). The contribution of trait-mediated indirect effects to the net effects of a predator. *Proceedings of the National Academy of Sciences of the United States of America*, 98, 3904-3908. https://doi.org/10.1073/pnas.071061998

Pritchard, D. W., Paterson, R., Bovy, H. C., & Barrios-O'Neill, D. (2017). Frair: An R package for fitting and comparing consumer functional responses. *Methods in Ecology and Evolution*, 8, 1528-1534.

Rall, B. C., Guill, C., & Brose, U. (2008). Food web connectance and predator interference dampen the paradox of enrichment. *Oikos*, 117, 202-213. https://doi.org/10.1111/j.2007-0030-1299.15491.x

Rogers, D. (1972). Random search and insect population models. *Journal of Animal Ecology*, 41, 369-383. https://doi.org/10.2307/3747

Schmitz, O. J., Krivan, V., & Ovadia, O. (2004). Trophic cascades: The primacy of trait-mediated indirect interactions. *Ecology Letters*, 7, 153-163. https://doi.org/10.1111/j.1461-0248.2003.00560.x

Sentis, A., Bertram, R., Dardenne, N., Ramon-Portugal, F., Louit, I., Trionnaire, G. L., Simon, J.-C., Magro, A., Pujol, B., Hemptinne, J. L., & Danchin, E. (2019). Different phenotypic plastic responses to predators observed among aphid lineages specialized on different host plants. *Scientific Reports*, 9, 9017. https://doi.org/10.1038/s41598-019-45220-0

Sentis, A., Binzer, A., & Boukal, D. S. (2017). Temperature-size responses alter food chain persistence across environmental gradients. *Ecology Letters*, 20, 852-862. https://doi.org/10.1111/ele.12779

Sentis, A., & Boukal, D. S. (2018). On the use of functional responses to quantify emergent multiple predator effects. *Scientific Reports*, 8, 11787. https://doi.org/10.1038/s41598-018-30244-9

Sentis, A., Gémard, C., Jaugeon, B., & Boukal, D. S. (2016). Predator diversity and environmental change modify the strengths of trophic and nontrophic interactions. *Global Change Biology*, 23, 2629-2640. https://doi.org/10.1111/gcb.13560

Sentis, A., Hemptinne, J. L., & Brodeur, J. (2013). How functional response and productivity modulate intraguild predation. *Ecosphere*, 4, 46. https://doi.org/10.1890/ES12-00379.1
Sih, A., Englund, G., & Wooster, D. (1998). Emergent impacts of multiple predators on prey. *Trends in Ecology & Evolution, 13*, 350-355. https://doi.org/10.1016/S0169-5347(98)01437-2

Skalski, G. T., & Gilliam, J. F. (2001). Functional responses with predator interference: Viable alternatives to the Holling type II model. *Ecology*, 82, 3083-3092.

Soetaert, K., & Petzoldt, T. (2010). Inverse modelling, sensitivity and Monte Carlo analysis in R using package FME. *Journal of Statistical Software*, 33, 1-28.

Soluk, D. A. (1993). Multiple predator effects: Predicting combined functional response of stream fish and invertebrate predators. *Ecology*, 74, 219-255. https://doi.org/10.2307/1939516

Suárez-Morales, E., Wasserman, R. J., & Dalu, T. (2015). A new species of *Lovenula* Schmeil (Copepoda, Diaptomidae) from the Eastern Cape province of South Africa. *Crustaceana, 88*, 324-342.

Trussell, G. C., Ewanchuk, P. J., Bertness, M. D., & Silliman, B. R. (2004). Trophic cascades in rocky shore tide pools: Distinguishing lethal and nonlethal effects. *Oecologia*, 139, 427-432. https://doi.org/10.1007/s00442-004-1512-8

Trussell, G. C., Ewanchuk, P. J., & Matassa, C. M. (2008). Resource iden-
tity modifies the influence of predation risk on ecosystem function.
*Ecology*, 89, 2798-2807. https://doi.org/10.1890/08-0250.1

Tylianakis, J. M., & Romo, C. M. (2010). Natural enemy diversity and biological control: Making sense of the context-dependency. *Basic and Applied Ecology, 11*, 657-668. https://doi.org/10.1016/j.baae.2010.08.005

Uiterwaal, S., Dell, A. I., & DeLong, J. P. (2018). Arena size modulates functional responses via behavioral mechanisms. *Behavioral Ecology*, 30, 483-489. https://doi.org/10.1093/beheco/ary188

Uiterwaal, S. F., & DeLong, J. P. (2018). Multiple factors, including arena size, shape the functional responses of ladybird beetles. *Functional Ecology*, 55, 2429-2438.

Vance-Chalcraft, H. D., & Soluk, D. A. (2005). Multiple predator ef-
fects result in risk reduction for prey across multiple prey den-
sities. *Oecologia*, 144, 472-480. https://doi.org/10.1007/s0044
2-005-0077-5

Veselý, L., Boukal, D. S., Buřič, M., Kuklina, I., Fořt, M., Yazicioglu, B., Prchal, M., Kozák, P., Kouba, A., & Sentis, A. (2019). Temperature and prey density jointly influence trophic and non-trophic interactions in multiple predator communities. *Freshwater Biology*, 64, 1984-1993. https://doi.org/10.1111/fwb.13387

Wasserman, R. J., Alexander, M. E., Barrios-O'Neill, D., Weyl, O. L. F., & Dalu, T. (2016). Using functional responses to assess predator hatching phenology implications for pioneering prey in arid temporary pools. *Journal of Plankton Research*, 38, 154-158. https://doi.org/10.1093/plankt/fbv114

Wasserman, R. J., Alexander, M. E., Dalu, T., Ellender, B. R., Kaiser, H., & Weyl, O. L. F. (2016). Using functional responses to quantify interaction effects among predators. *Functional Ecology*, 30, 1988-1998. https://doi.org/10.1111/1365-2435.12682

Wasserman, R. J., & Froneman, P. W. (2013). Trophic level stability-
inducing effects of predaceous early juvenile fish in an estuarine me-
socosm study. *PLoS One*, 8, e61019.

Wasserman, R. J., Weston, M., Weyl, O. L. F., Froneman, P. W., Welch, R. J., Vink, T. J. F., & Dalu, T. (2018). Sacrificial males: The potential role of copulation and predation in contributing to copepod sex-skewed ratios. *Oikos*, 127, 970-980. https://doi.org/10.1111/oik.04832

SUPPORTING INFORMATION

Additional supporting information may be found online in the
Supporting Information section.

How to cite this article: Cuthbert RN, Dalu T, Wasserman RJ, et al. Prey and predator density-dependent interactions under different water volumes. *Ecol Evol.* 2021;11:6504–6512. https://doi.org/10.1002/ece3.7503