Introduction
A fundamental goal of ecology is to understand mechanisms that regulate the dynamics of populations and communities. In wet Neotropical forest ecosystems, a diverse assemblage of small terrestrial frogs occupies the leaf-litter layer on the forest floor, and a large body of literature has suggested that the ‘leaf-litter frog assemblage’ is regulated bottom-up (sensu Hunter & Price, Reference Hunter and Price1992) by the abundance of leaf litter. Species in the leaf-litter assemblage use litter to obtain abundant arthropod food resources (Toft Reference Toft1980, Lieberman Reference Lieberman1986, Whitfield & Donnelly Reference Whitfield and Donnelly2006), to select habitat to limit desiccation across life stages (including during oviposition; Seebacher and Alford Reference Seebacher and Alford2002; Schlaepfer Reference Schlaepfer2003; Socci et al. Reference Socci, Schlaepfer and Gavin2005) and to take refuge from a diverse array of predators (Talbot Reference Talbot1979, Greene Reference Greene, Almeda and Pringle1988, Cooper et al. Reference Cooper, Caldwell and Vitt2008a, Reference Cooper, Caldwell and Vitt2008b). Several studies have documented a positive relationship between standing leaf litter and abundance of terrestrial frogs and anoles in mainland ecosystems (Scott Reference Scott1976, Lieberman Reference Lieberman1986, Guyer Reference Guyer1988, Fauth et al. Reference Fauth, Crother and Slowinski1989, Heinen Reference Heinen1992, Whitfield et al. Reference Whitfield, Bell, Philippi, Sasa, Bolaños, Chaves, Savage and Donnelly2007), and Whitfield et al. (Reference Whitfield, Reider, Greenspan and Donnelly2014) provided the first experimental demonstration of standing leaf litter as a density-limiting factor. The leaf-litter frog assemblage declined 75% in abundance over a 35-year period at La Selva Biological Station in Costa Rica, and a climate-driven decline in standing leaf-litter mass was identified as the most likely source of this enigmatic amphibian decline (Whitfield et al., Reference Whitfield, Bell, Philippi, Sasa, Bolaños, Chaves, Savage and Donnelly2007, Reference Whitfield, Reider, Greenspan and Donnelly2014). Thus, a common thread among studies is the identification of leaf litter as a fundamental resource regulating the abundance of the leaf-litter frog assemblage, in ways consistent with bottom-up regulation (sensu Hunter and Price Reference Hunter and Price1992).
While a large literature has studied resource limitation and bottom-up regulation of terrestrial frogs in the Neotropics (e.g., Whitfield et al. Reference Whitfield, Reider, Greenspan and Donnelly2014, Reference Whitfield, Lips and Donnelly2016), much less consideration has been given how predators can exert top-down effects on frogs in the leaf-litter assemblage. The lack of information on top-down effects is surprising because the leaf-litter frog assemblage is prey to a diverse assemblage of both vertebrate and invertebrate predators (Greene Reference Greene, Almeda and Pringle1988). For example, La Selva Biological Station in Costa Rica (hereafter, La Selva) is home to an extremely species-rich community of vertebrate predators, but one of the most frequent predators of small frogs may actually be invertebrates. In particular, large Wandering Spiders (Ctenidae) are a common generalist predator occupying terrestrial environments in lowland Neotropical wet forests of Costa Rica. Ctenids attack and consume terrestrial and arboreal frogs in the families Centrolenidae, Craugastoridae and Hylidae (Hayes Reference Hayes1983, Szelistowski Reference Szelistowski1985, Donnelly & Guyer Reference Donnelly and Guyer1994, Lapinski & Tschapka Reference Lapinski and Tschapka2013, Folt & Lapinski Reference Folt and Lapinski2017) but frequently reject the poisonous frog Oophaga pumilio (Szelistowski Reference Szelistowski1985, Murray et al. Reference Murray, Bolton, Berg and Saporito2016, Folt & Lapinski Reference Folt and Lapinski2017). Ctenids are also capable of consuming small reptiles, being the most commonly observed source of predation during detailed studies of the anole Norops humilis (Guyer Reference Guyer1988), and ctenids are dominant predators of metamorphic hylid frogs emerging from ephemeral swamps (Donnelly & Guyer Reference Donnelly and Guyer1994). The latter study observed a strong pulse in metamorphosis of frogs and hypothesised that this synchronous emergence may have evolved as a mechanism to satiate terrestrial spider predators (Donnelly & Guyer Reference Donnelly and Guyer1994, Guyer & Donnelly Reference Guyer, Donnelly, Donnelly, Crother, Guyer, Wake and White2005). Recent studies from Central and South America have also found that ctenids are opportunistic predators of small vertebrates (Menin et al. Reference Menin, de Jesus Rodrigues and Salette de Azevedo2005, Maffei et al. Reference Maffei, Ubaid and Jim2010, Nyffeler & Altig Reference Nyffeler and Altig2020, Reyes-Olivares et al. Reference Reyes-Olivares, Guajardo-Santibáñez, Segura, Zañartu, Penna and Labra2020, Meneses et al. Reference Meneses, Corrêa, Fernades, Lopes, Citeli and Brandão2021), and the accumulating literature suggests that spiders may play an important role as vertebrate predators across Neotropical wet forests (Hayes Reference Hayes1983, Guyer Reference Guyer1988, Donnelly & Guyer Reference Donnelly and Guyer1994, Menin et al. Reference Menin, de Jesus Rodrigues and Salette de Azevedo2005). However, no study has examined the extent to which ctenid spiders influence terrestrial frogs or other small vertebrates in a top-down, predator–prey context.
In recent decades, amphibian populations have declined significantly across the world, and declines have been linked to habitat change, invasive species, climate change, pollution, contaminants and emerging infectious diseases (Stuart et al. Reference Stuart, Chanson, Cox, Young, Rodrigues, Fischman and Waller2004, Mendelson et al. Reference Mendelson, Lips, Gagliardlo, Rabb, Collins, Diffendorfer, Daszak, Roberto Ibáñz, Zippel, Lawson, Wright, Stuart, Gascon, da Silva, Burrowes, Joglar, La Marca, Lötters, du Preez, Weldon, Hyatt, Rodriguez-Mahecha, Hunt, Robertson, Lock, Raxworthy, Frost, Lacy, Alford, Campbell, Parra-Olea, Bolaños, Domingo, Halliday, Murphy, Wake, Coloma, Kuzmin, Price, Howell, Lau, Pethiyagoda, Boone, Lannoo, Blaustein, Dobson, Griffiths, Crump, Wake and Brodie2006, Whitfield et al. Reference Whitfield, Lips and Donnelly2016, Scheele et al. Reference Scheele, Pasmans, Skerrat, Berger, Martel, Beukema, Acevedo, Burrowes, Carvalho, Catenazzi, de la Riva, Fisher, Flechas, Foster, Frías-Álvarez, Garner, Gratwicke, Guayasamin, Hirschfield, Kolby, Kosch, La Marca, Lindenmayer, Lips, Longo, Maneyra, Mcdonald, Mendelson, Palacios-Rodriguez, Parra-Olea, Richards-Zawacki, Rödel, Rovito, Soto-Azat, Felipe Toledo, Voyles, Weldon, whitfield, Wilkinson, Zamudio and Canessa2020). While the current paradigm for the leaf-litter frog assemblage suggests that it is bottom-up regulated by leaf litter, understanding top-down effects may be particularly important because predation pressure may interact with resources in ways important for understanding the amphibian decline crisis in Central America and elsewhere. Because leaf litter has declined strongly at sites like La Selva (Whitfield et al., Reference Whitfield, Bell, Philippi, Sasa, Bolaños, Chaves, Savage and Donnelly2007) and the terrestrial frog assemblage uses leaf litter to hide from predators (Talbot Reference Talbot1979, Greene Reference Greene, Almeda and Pringle1988, Cooper et al. Reference Cooper, Caldwell and Vitt2008a, Reference Cooper, Caldwell and Vitt2008b), decreased litter availability may have caused predator–prey dynamics to shift such that frogs are now experiencing increased interactions with predators. Thus, understanding interactions between top-down and bottom-up regulatory factors might also provide insights into novel mechanisms for population regulation of declining terrestrial frogs in the Neotropics.
To this end, we sought to better understand top-down effects on the ecology of leaf-litter frogs in lowland wet Neotropical forests where ctenid spiders are frequent predators on terrestrial anurans. We used multi-species occupancy models to quantify patterns of site occupancy and detection for two small terrestrial frog species (Craugastor bransfordii and Oophaga pumilio) at La Selva Biological Station, and we modelled how frog occupancy and detection probability were influenced top-down by predatory spiders and bottom-up by resources (arthropod prey and leaf-litter microhabitat). Given the hypothesis that ctenid spiders are important predators of terrestrial frogs in the Neotropics (Menin et al. Reference Menin, de Jesus Rodrigues and Salette de Azevedo2005, Maffei et al. Reference Maffei, Ubaid and Jim2010, Folt & Lapinski Reference Folt and Lapinski2017, von May et al. Reference von May, Biggi, Cárdenas, Isabel Diaz, Alarcón, Herrera, Santa-Cruz, Tomasinelli, Westeen, Sánchez-Paredes, Larson, Title, Grundler, Grundler, Davis Rabosky and Rabosky2019, Nyffeler & Altig Reference Nyffeler and Altig2020, Reyes-Olivares et al. Reference Reyes-Olivares, Guajardo-Santibáñez, Segura, Zañartu, Penna and Labra2020, Meneses et al. Reference Meneses, Corrêa, Fernades, Lopes, Citeli and Brandão2021), we predicted that small terrestrial frogs would experience elevated occupancy in areas where predatory spiders are absent relative to where predators are present. To understand potential interactions between top-down and bottom-up factors, we also evaluated how top-down effects interacted with resources, by modelling how predator–prey interactions varied along a leaf-litter resource gradient.
Materials and methods
Study site and taxa
La Selva Biological Station (hereafter, La Selva) is a private reserve owned by the Organization for Tropical Studies (OTS) in the Caribbean lowlands of northeastern Costa Rica (10.42˚N, 84.02˚W; WGS 84). The site is characterised by an average temperature of 25.8˚C, receives ca. 4 m of precipitation/year and is classified as a Tropical Wet Forest in the Holdridge life zone system (McDade et al. Reference Mcdade, Bawa, Hespenheide and Hartshorn1994).
We studied two common terrestrial frogs and an assemblage of ctenid spiders from La Selva as the focal taxa (Figure 1). Craugastor bransfordii is an abundant frog species (Craugastoridae) that is highly variable in colouration; variation has been categorised into at least four distinct morphs defined by colour and dorsal ridging patterns, which are thought to decrease predation in the terrestrial leaf-litter environment (Savage & Emerson Reference Savage and Emerson1970, Cooper et al. Reference Cooper, Caldwell and Vitt2008b). Oophaga pumilio is an aposematic poison frog (Dendrobatidae) that also occupies terrestrial habitats; the species primarily eats ants and mites (Donnelly Reference Donnelly1991), a diet which provides alkaloid compounds that are sequestered into poison glands in the frog’s skin and confers an antipredatory defence (Saporito et al. Reference Saporito, Garraffo, Donnelly, Edwards, Longino and Daly2004, Reference Saporito, Donnelly, Norton, Garraffo, Spande and Daly2007a, Stynoski et al. Reference Stynoski, Torres-Mendoza, Sasa-Marin and Saporito2014b). Across its geographic distribution, O. pumilio is brightly coloured, which serves as an aposematic signal to reduce predation (Saporito et al. Reference Saporito, Zuercher, Roberts, Kenneth and Donnelly2007b). Ctenid spiders are common generalist predators occupying terrestrial environments in lowland Caribbean forests of Costa Rica. Three species in the genus Ctenus are conspicuous residents at La Selva (Ctenus curvipes, Ctenus sinuatipes and an undescribed Ctenus species; W. Lapinski, pers. comm.) which consume small vertebrates (Szelistowski Reference Szelistowski1985, Guyer Reference Guyer1988, Lapinski & Tschapka Reference Lapinski and Tschapka2013). In primary forest habitats like La Selva, these ctenid species are visual ambush predators that emerge from hidden sites in soil and leaf litter (Lapinski & Tschapka Reference Lapinski and Tschapka2013) to attack prey (e.g., Folt and Lapinski Reference Folt and Lapinski2017). Identification of Ctenus species is difficult in the field because the species are morphologically similar and also exhibit high intraspecific variation. Thus, we considered all ctenids as a single group in our analyses.
Two species of abundant terrestrial frogs and a predatory spider (Ctenidae) from lowland Caribbean forests of Costa Rica: (A) Craugastor bransfordii, (B) Oophaga pumilio and (C) Cupiennius sp. consuming a rain frog (Pristimantis ridens). See Folt and Lapinski (Reference Folt and Lapinski2017) for other examples of ctenid spiders consuming various frog species from the Caribbean lowlands of Costa Rica. Photographs by J. Folt (A, B) and W. Lapinski (C).
Because the focal frogs are relatively mobile compared to sit-and-wait predatory ctenids (B. Folt, pers. obs.), prey species should respond to predation threat through shifts in space, habitat use or behaviour (i.e., the prey response dominates; Sih Reference Sih1984). Therefore, if spider predation exerts a strong ecological pressure on frogs, we predicted that, relative to sites where spiders were absent, the presence of spiders would decrease patterns of occupancy and detection of frogs (Andrews Reference Andrews1979, Losos Reference Losos2009). However, because O. pumilio is chemically defended and ctenids have rejected O. pumilio during feeding trials (Szelistowski Reference Szelistowski1985, Murray et al. Reference Murray, Bolton, Berg and Saporito2016), we also predicted that O. pumilio would occupy sites independent of spiders as a result of a lack of perceived predation threat.
Data collection
We established replicate gridded plots around individual canopy trees (N = 14) in primary forest habitat; plot grids consisted of 21 cells that were each 3 m × 3 m in area (189 m2 total) and were marked using PVC tubing (Figure S1; Folt Reference Folt2017). We used plots of this area because comparable-sized plots were successfully used in other studies of amphibian and reptile population ecology at La Selva (Guyer Reference Guyer1988, Donnelly Reference Donnelly1989, Whitfield et al. Reference Whitfield, Reider, Greenspan and Donnelly2014).
All 21 plot cells were sampled for frogs and spiders in February, May, June, August, September, October and November 2014. We surveyed animals using diurnal visual encounter searches of the forest floor and other substrates up to 2 m above ground. We walked methodically through plots to search all cells and carefully prodded through the litter and other structures with a 1.2-m pole. Each plot cell was surveyed repeatedly (N = 3) during each month. We conducted the first survey on a randomly selected day and generally repeated surveys on the two consecutive days. Our sampling method focused on counting frogs, lizards and spiders present on the forest floor, or immediately adjacent vertical structures up to 1 m above the ground. However, the vast majority of samples were taken from the forest floor.
We quantified abundance of leaf litter and terrestrial arthropods in four randomly selected plot cells in each month. For leaf-litter and arthropod samples, we randomly selected a corner of each plot cell and used a 0.30 m × 0.30 m PVC frame to delimit samples of litter, cut the perimeter of the litter with a utility blade and collected the sample into a plastic bag. Before cutting out the litter, we measured leaf-litter depth (cm) and the number of leaves in the litter column in each corner of the PVC frame. We brought the samples into the laboratory and dried them in Berlese funnels beneath heated incandescent bulbs for 24–48 hours or until no moisture appeared to remain in the samples; during the drying process, arthropods were driven downward from within the drying litter and collected in vials containing 95% ethanol. We measured dried leaf-litter mass (g) and used a microscope to sort arthropods into taxonomic groups that are important components of the diets of the study species (Whitfield & Donnelly Reference Whitfield and Donnelly2006, Folt Reference Folt2017). Because leaf-litter and arthropod samples were collected from the corner of plot cells (Figure S1), we used the frog and spider survey data from each plot cell with resources covariates and any plot cells immediately adjacent to a resource covariate sample to model predator–prey occupancy. Our additional estimates of leaf-litter abundance (depth and number of leaves) were characterised by strong positive relationships with dried leaf-litter mass (B. Folt, unpubl. results); however, we chose to use estimates of dried leaf-litter mass as our primary assay of leaf-litter abundance because this estimate controls for moisture, which is highly variable in space and time due to rainfall in tropical wet forests like La Selva, and was paired with our estimate of arthropod food resources in the samples.
Statistical analysis
When species are detected imperfectly (detection probability < 1.0), analyses of co-occurrence patterns can be misleading and result in erroneous inferences (MacKenzie et al. Reference Mackenzie, Nichols, Royle, Pollock, Bailey and Hines2006). We took advantage of recent advances in quantitative ecology that have developed statistical methods to account for imperfect detection of individuals while estimating species occupancy or abundance at sites (occupancy modelling and derivations; MacKenzie et al. Reference Mackenzie, Nichols, Royle, Pollock, Bailey and Hines2006). In particular, such occupancy models have been expanded to include parameterisations that investigate relationships among co-occurring and potentially interacting species, such as predator–prey relationships (MacKenzie et al. Reference Mackenzie, Bailey and Nichols2004); these multi-species occupancy models are flexible and can include covariates, to evaluate support for competing hypotheses, such as resources, in addition to interspecific interactions (Richmond et al. Reference Richmond, Hines and Beissinger2010). Thus, we modelled how the spatial and behavioural ecology of frogs was influenced by predatory spiders and resources (leaf litter and arthropod food resources) using the ΨBa parameterisation of the multi-species, single-season occupancy model (MacKenzie et al. Reference Mackenzie, Bailey and Nichols2004, Richmond et al. Reference Richmond, Hines and Beissinger2010). The ΨBa parameterisation estimates the probability of occupancy of a dominant species (A) and a subordinate species (B). We assumed that spiders were dominant over frog species because we hypothesized that prey frogs would avoid predation risk and respond predictably to the presence of predatory spiders. Using these models, we estimated the probability of occupancy of the dominant predatory organisms (ΨA) and the occupancy of each subordinate taxon when dominant predators were present (ΨBA) and absent (ΨBa).
We evaluated co-occurrence patterns between ctenids and each focal frog species separately. Preliminary analyses indicated that spiders and frogs consistently differed in patterns of occupancy and detection for the two predator–prey pairs (ctenids + C. bransfordii and ctenids + O. pumilio). Thus, we used a model-building process that began with creating an initial model describing variance in occupancy and detection differing between dominant and subordinant species and then built six models that included all combinations of Ψ parameters (Table 1) describing hypotheses of how predatory spiders influence occupancy and/or detection of prey species (Table 2). Using this model-building process, we evaluated whether occupancy of the subordinate species was affected by the presence of the dominant species (ΨBA ≠ ΨBa) and whether detection probability of the subordinate species was influenced by detection (r BA) and/or occupancy (r Ba) of the dominant species (Table 2). Our model-building process generated a balanced set of six models with even representation of predator effects (ΨBA ≠ ΨBa, r BA and r Ba) either present or absent in models.
Description of model parameters used to evaluate hypotheses describing two-species co-occurrence patterns, following Richmond et al. (Reference Richmond, Hines and Beissinger2010). Species A are predators and are assumed to be dominant over the prey species B.
Six models describing hypotheses for co-occurrence patterns between predatory spiders and prey frogs in forest plots at La Selva, Costa Rica. Species A is assumed to be dominant over species B. All models include terms describing species A and B differing in probability of occupancy and detection.
To evaluate the relative roles of resources and predators contributing to frog ecology, we built 54 additional models to evaluate how frog occupancy and detection are explained by predators and resources. We expanded the initial six models (Table 2) to include covariate effects (litter mass and arthropod abundance) on occupancy, detection, and both occupancy and detection, in the presence and absence of predator effects. Because rainfall is strongly seasonal at La Selva and seasonal rainfall is thought to influence patterns of leaf-litter deposition, arthropod abundance, and abundance of terrestrial frogs and lizards (Frankie et al. Reference Frankie, Baker and Opler1974, Lieberman & Dock Reference Lieberman and Dock1982, Guyer Reference Guyer1988), we also built models evaluating how seasonality influenced patterns of occupancy and detection. The second model-building process generated 54 additional models with different combinations of predator, food resources, seasonal effects and null models in a completely balanced design (N = 60 models; Supplementary Table 1). We modelled the arthropod covariate differently for each species depending on its diet (Lieberman Reference Lieberman1986, Whitfield & Donnelly Reference Whitfield and Donnelly2006): the sum of Acari, Araneae, Coleoptera and Isopoda abundance for C. bransfordii, and Acari and Formicidae for O. pumilio. We transformed all covariates by log(x+1) to improve model convergence. We classified seasonality by categorising months as ‘dry’ (February–May) and ‘wet’ (June–November) based on historical mean monthly precipitation values (Guyer Reference Guyer1988) and by heuristically assessing rainfall data from La Selva during 2014 (B. Folt, pers. obs.).
We evaluated how well each model fit the data using Akaike’s Information Criterion corrected for small sample size (AICc; Hurvich and Tsai Reference Hurvich and Tsai1989) and calculated the probability a given model within a set is the best approximating model (model weight; w m). We then used model weights to calculate model-averaged coefficients (± unconditional SE) and model-averaged parameter weights (w p) across all models with w m > 0.025. We considered parameters as explaining significant variation when they occurred in a model in the top model set (ΔAICc < 2.00), unconditional SE estimates did not overlap, and w p > 0.70 (Burnham & Anderson Reference Burnham and Anderson2010).
Our analyses indicated significant predator–prey interactions in the top-model set by supporting the ΨBa parameter in the top model; we therefore calculated a species interaction factor (SIF). The SIF describes if the predator–prey pairs co-occur more frequently (SIF > 1) or less frequently (SIF < 1) than expected if species did not interact and occupancy probabilities were independent of each other (SIF ˜ 1). We calculated the SIF using model-averaged parameter estimates and by following the SIF equation described by Richmond et al. (Reference Richmond, Hines and Beissinger2010):
$${\rm{SIF}} = \;{{{\Psi ^{\rm{A}}}{\rm{*}}{\Psi ^{{\rm{BA}}}}} \over {{\Psi ^{\rm{A}}}\left( {{\Psi ^{\rm{A}}}{\rm{*}}{\Psi ^{{\rm{BA}}}} + \left( {1 - {\Psi ^{\rm{A}}}} \right){\rm{*}}{\Psi ^{{\rm{Ba}}}}} \right)}}$$
We performed all analyses in the statistical Program R (R; R Core Team 2018). We built single-season, two-species occupancy models using the function occmod() in the package ‘RPresence’ (MacKenzie & Hines Reference Mackenzie and Hines2016). In our analysis, we treated plot cells with covariates within each month as ‘sites’ (sensu MacKenzie et al. Reference Mackenzie, Nichols, Royle, Pollock, Bailey and Hines2006), such that individual plot cells sampled repeatedly among months were considered independent replicate samples. Our data, model code, and software packages are available on GitHub (https://github.com/brianfolt/predator-prey-models.git).
Results
During February–November 2014, we made 946 detections of C. bransfordii, 343 detections of O. pumilio and 271 detections of ctenid spiders during 3144 surveys of 1048 plot cells. Top-model sets of co-occurrence models for each predator–prey pair contained parameters describing interactions between predators and prey (ΨBa, r BA and r Ba) and effects of leaf litter on occupancy and detection (ΨLL and p LL; Table 3).
Number of parameters (k), AICc, ΔAICc and model weight (w m) for the top model set (w m > 0.10) among 60 models built to describe how occupancy and detection patterns of the frogs Craugastor bransfordii and Oophaga pumilio are influenced by predatory ctenid spiders and resources (leaf litter, arthropods) at La Selva, Costa Rica. See Table 1 for explanations of model parameters.
The top model for both frog species included the terms ΨBa, ΨLL, r BA and r Ba (Table 3), a model which indicated effects of predators (ΨBa) and leaf litter (ΨLL) on frog occupancy and a predator effect on frog detection (r BA). For C. bransfordii, the top model (w m = 0.41) was over one and a half times more likely than other models. For O. pumilio, the top model (w m = 0.35) was rivaled by an equivocally supported model, which included r Ba (w m = 0.33). Model averaging of parameters strongly supported three parameters for both frogs: ΨBa (C. bransfordii, wp = 0.99; O. pumilio, wp = 0.99), ΨLL (C. bransfordii, w p = 0.98; O. pumilio, wp = 0.97) and r BA (C. bransfordii, wp = 0.99; O. pumilio, wp = 0.99; Table 4 and 5). While the parameter r Ba was included in top-model sets for both species, it received low support by model averaging (w p < 0.70) for both frog species.
Model-averaged coefficient values (β), unconditional standard error (SE) and parameter weights (w p) generated by averaging 60 co-occurrence models (Supplementary Table 1) for Craugastor bransfordii and Oophaga pumilio and ctenid spider predators at La Selva, Costa Rica. See Table 1 for explanations of model parameters. Parameter weights were not provided for ΨA and ΨBA because these parameters were included in all models (w p = 1.00); parameters are not reported for ΨArthropods and ΨSeason because they received no support.
Model-averaged coefficient values (β), unconditional standard error (SE) and Akaike parameter weights (w p) generated by averaging 60 co-occurrence models (Supplementary Table 1) describing patterns of detection for the frogs Craugastor bransfordii and Oophaga pumilio and predatory ctenid spiders at La Selva, Costa Rica. See Table 1 for explanations of model parameters. Parameter weights were not provided for p A and p B, because they were included in all models (w p = 1.00); parameters are not reported for ΨArthropods and ΨSeason because they received no support.
Model-averaged coefficients described frog occupancy as increasing with leaf-litter mass and, relative to when spiders were present, frog occupancy was elevated when spiders were absent (Figure 2). Occupancy was more strongly predicted by the absence of predators (ΨBa) than the abundance of leaf litter (ΨLL) or arthropods (ΨArthropods) for both C. bransfordii and O. pumilio (Table 4). These patterns generated SIFs for frogs that described predator avoidance at sites with low leaf litter; however, at sites with increased leaf litter, co-occurrence patterns were independent between frogs and spiders (Figure 3). Detection probability of spiders (p A = 0.10–0.12) was lower than the detection of C. bransfordii (p B = 0.29) and O. pumilio (p B = 0.25), but detection probability of both species increased at sites where spiders were occupants and were detected (r BA = 0.45, 0.40, respectively; Figure 4). Arthropod abundance and seasonality received zero support as effects on frog occupancy and detection (Table 4 and 5).
Site occupancy of Craugastor bransfordii (A) and Oophaga pumilio (B) as a function of leaf-litter mass at sites conditionally occupied (ΨBA) and unoccupied (ΨBa) by predatory spiders (Ctenidae). ΨA is the occupancy of spiders. The black arrow indicates the median value of leaf-litter mass for both panels. Results were model-averaged across a balanced design of 60 models (Supplementary Table 1).
Species interaction factors (SIF) between predatory spiders (Ctenidae) and Craugastor bransfordii (A) and Oophaga pumilio (B) in relation to leaf-litter mass in forest plots at La Selva, Costa Rica. Ctenid spiders are assumed to be dominant over prey frogs. SIF values < 1.0 indicate avoidance of the subordinate species, values > 1.0 indicate aggregated with dominate species and values ˜ 1.0 (grey dotted line) indicates independent patterns of co-occupancy between spiders and frogs. The black arrows indicates the median value of leaf-litter mass.
Conditional detection probability of predatory spiders (dominant species; Ctenus sps.) and prey frogs (subordinate species; A – Craugastor bransfordii, B – Oophaga pumilio) in terrestrial habitats at La Selva, Costa Rica. p A – detection probability of dominant species, given the subordinate species is absent; p B – detection of the subordinate species, given the dominant species is absent; r BA – detection of subordinate species, given the dominant species is present and detected.
Discussion
In our study, site occupancy of both C. bransfordii and O. pumilio in plot cells was more strongly predicted by the presence of spiders than other resources. Specifically, both frogs exhibited decreased occupancy when spiders were present, suggesting that frogs selected habitat to reduce co-occurrence with spider predators. Occupancy of both frogs and spiders was also strongly influenced by leaf-litter abundance at sites, and we observed an effect in which predator avoidance by frogs decreased as litter abundance increased. Variance in frog detection probability was also consistent with predator avoidance: relative to when spiders were absent, detection of C. bransfordii and O. pumilio was elevated when spiders were present and detected within plots, suggesting that frogs exhibit more vigilant behaviour when spider predators were present and conspicuous in plots. These results were consistent with the hypothesis that spider predators exert strong ecological effects on prey vertebrates in Central American forests.
Guyer (Reference Guyer1988) observed seasonal abundance cycles for populations of a terrestrial anole (Norops humilis) in cacao plantations that oscillated in tandem with the leaf-drop phenology of cacao trees. Because seasonality of leaf litter drives arthropod abundance (Lieberman & Dock Reference Lieberman and Dock1982) and N. humilis was experimentally determined to be food limited, Guyer (Reference Guyer1988) hypothesised that leaf-fall events at the scale of single canopy trees regulate abundances of terrestrial anoles and similar vertebrates, such as terrestrial frogs. The Guyer (Reference Guyer1988) hypothesis assumes that populations are regulated bottom-up by leaf-litter and/or food resources and predicts that (1) anole and frog abundance through primary forest is a mosaic of sites at different stages depending on the leaf-drop phenology of the most proximate tree, and (2) low-resource patches should be vacated seasonally by individuals in search of areas of greater resource quality (hereafter, the litter-mosaic model; Guyer Reference Guyer1988). Our results have implications for reinterpreting and revising the litter-mosaic model in two ways. First, site occupancy of both frogs and N. humilis (Folt Reference Folt2017) was positively influenced by abundance of leaf litter, but not by the abundance of arthropods used as food. This result suggests that leaf litter alone may drive variance in the litter-mosaic model (i.e., Whitfield et al. Reference Whitfield, Reider, Greenspan and Donnelly2014). Second, because frog occupancy was more greatly influenced by predators than by litter abundance, predators appear to be an underemphasised control on anuran ecology in a patch-mosaic context. Specifically, our results suggested that frog demography is most sensitive to the presence of predators and only secondarily influenced by the availability of resources, such as leaf litter. In the context of the litter-mosaic model, frogs may still be incentivised to occupy sites of diminished leaf litter or vacate sites with abundant litter, depending on the presence of predators. Thus, our results revise the litter-mosaic model by placing greater emphasis on leaf litter and predators as mechanisms influencing site occupancy of small terrestrial vertebrates in a patch-mosaic model.
We observed an interaction between the predators and resources, where predator avoidance by frogs decreased with increasing abundance of leaf litter. The observed interaction between predators and litter suggests that top-down and bottom-up effects on terrestrial frogs do not occur independently in this ecosystem, but rather that habitat characteristics influence species’ interactions. Further, the relationship between leaf litter and predator–prey interactions also has implications for understanding seasonal abundance cycles and amphibian declines at La Selva and other comparable wet tropical forests. Because terrestrial leaf-litter abundance is strongly seasonal in lowland Caribbean wet forests like La Selva, our results suggest that food web interactions between spiders and frogs may vary temporally at the seasonal scale (Tunney et al. Reference Tunney, Mccann, Lester and Shuter2012). Leaf-drop phenology of most trees causes standing leaf litter to reach its greatest depth at the end of the dry season (Frankie et al. Reference Frankie, Baker and Opler1974, Levings & Windsor Reference Levings and Windsor1984, Folt Reference Folt2017); as the wet season begins and advances, decreased leaf drop, increased decomposition, and mechanical action of rainfall and shredding action of insects cause standing litter to decrease. While models suggested that predator–prey interactions vary across a gradient of litter abundance, ctenid spiders might exert greater effects on terrestrial frogs during periods of low leaf litter in the wet season. In these situations, decreased occupancy of frogs in the presence of spiders may be driven by reduced refugia in litter, increased risk of predation and/or greater incentive to vacate sites for areas with higher leaf litter. This effect can be extended with implications for models of amphibian population decline at La Selva, a site that experienced a 75% decline over a 35-year period that was hypothesised as a result of climate-driven reductions in the abundance of standing leaf litter (Whitfield et al. Reference Whitfield, Bell, Philippi, Sasa, Bolaños, Chaves, Savage and Donnelly2007). Given observed increase in predator–prey interactions at lower levels of litter abundance at La Selva, climate-driven declines in litter could have driven increased predator–prey interactions and a shift in community composition towards decreased abundance of frogs and increased abundance of spider predators. Thus, models of resource-driven amphibian declines (Whitfield et al. Reference Whitfield, Lips and Donnelly2016) should consider how altered predator–prey relationships may influence community dynamics in these ecosystems.
Because ctenid spiders are thought to avoid depredating O. pumilio, we predicted that O. pumilio would exhibit patterns of occupancy and detection independent of spiders. Contrary to this prediction, occupancy and detection patterns were consistent with spider avoidance and similar to those observed for a species thought to be frequently consumed by ctenids, C. bransfordii. The spider avoidance effect may be driven by (1) predation attempts by ctenids on O. pumilio that occasionally result in frog consumption (e.g., Summers Reference Summers1999), including at La Selva (M. Chaves, pers. comm.), (2) ctenid predation of O. pumilio eggs or tadpoles that have reduced chemical protection than adults (Stynoski et al. Reference Stynoski, Shelton and Stynoski2014a) and/or (3) an innate, general avoidance of spiders by frogs, because anurans in general have shared a long evolutionary history with spider predators. Because O. pumilio responded to spiders in ways consistent with predictions of a predator–prey relationship, our modelling results are at odds with behavioural studies of spider–O. pumilio interactions describing spiders as unimportant predators (Szelistowski Reference Szelistowski1985, Murray et al. Reference Murray, Bolton, Berg and Saporito2016).
In conclusion, our results provide evidence for strong effects of predatory spiders on prey frogs in Central America. Our results supported leaf litter as a bottom-up factor regulating the terrestrial vertebrate assemblage in lowland wet tropical forests but described a stronger and novel top-down effect of predators on anurans that has significant implications for improving models of patch dynamics, seasonal population cycles and amphibian declines (Guyer Reference Guyer1988, Whitfield et al. Reference Whitfield, Reider, Greenspan and Donnelly2014). While a large literature has suggested or assumed dominant bottom-up regulation of the terrestrial frogs and lizards in this system (Scott Reference Scott1976, Toft Reference Toft1980, Guyer Reference Guyer1988, Reference Guyer, McDade, Bawa, Hespenheide and Hartshorn1994, Fauth et al. Reference Fauth, Crother and Slowinski1989, Heinen Reference Heinen1992, Whitfield et al. Reference Whitfield, Bell, Philippi, Sasa, Bolaños, Chaves, Savage and Donnelly2007, Reference Whitfield, Reider, Greenspan and Donnelly2014, Folt & Reider Reference Folt and Reider2013), we provided the first evidence supporting a dominant effect of predators regulating this important anuran assemblage, a factor which appears habitat-dependent.
Acknowledgements
Funding support for this project was provided by the Organization for Tropical Studies (Christiane and Christopher Tyson Fellowship, the Dole Food Fellowship), the National Geographic Society (Young Explorer’s Grant), the Society for the Study of Amphibians and Reptiles (Andrew H. Price Field Research Grant), The Herpetologists’ League (E.E. Williams Research Grant), Auburn University (Thesis/Dissertation Research Award) and the Auburn University Museum of Natural History (Cunningham Foundation Graduate Research Assistantship). The manuscript benefited from helpful comments from M. Donnelly, S. Hermann, J. Feminella, C. McGowan, J. Hall, A. Tucker and C. McGowan. The study complied with permission from Auburn University IACUC #2013–1827 and the Costa Rican government (MINAET Resolución #009–2014). This paper is contribution no. 948 of the Auburn University Museum of Natural History.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S0266467421000274
Open access
