Journal Information
Vol. 21. Issue 3.
Pages 216-223 (July - September 2023)
Share
Share
Download PDF
More article options
Visits
1012
Vol. 21. Issue 3.
Pages 216-223 (July - September 2023)
Research Letters
Full text access
Does agricultural intensification impact pest regulation service by frogs in a natural multi-trophic system?
Visits
1012
Deyatima Ghosha,b, Parthiba Basub,
Corresponding author
pbzoo@caluniv.ac.in

Corresponding author.
a College of Biology and the Environment, Nanjing Forestry University, No. 159, Longpan Road, Xuanwu District, Nanjing City, Jiangsu Province 210037, China
b Ecology Research Unit, Department of Zoology, University of Calcutta, 35, Ballygunge Circular Road, Kolkata, 700019, India
Highlights

  • Frog density below a certain threshold cannot regulate rice pest population irrespective of agricultural intensity.

  • Due to strong intraguild predation between amphibians and arthropod pest predators, pest regulation is suppressed.

  • At advanced crop growth stage increased foliage provides refuge to pests causing a switch in prey preference of the frogs towards more agile arthropod pest predators.

  • Improving habitat heterogeneity will increase the frog density in the landscape and improve pest regulation.

This item has received
Article information
Abstract
Full Text
Bibliography
Download PDF
Statistics
Figures (5)
Show moreShow less
Tables (3)
Table 1. List of crop pest predators and pests observed in sweep net along with their relative abundance.
Table 2. Summary of the results from ANOVA Wald Chi-square test for significance testing of fixed effects (showing Chi-square, df and p values) and GLMM results (showing estimates, standard error, and p-values) between levels of fixed effects and random effects. Significant results (p-value < 0.05) are in bold.
Table 3. Results of posthoc test for the fixed effects with estimates, standard error, and p-values. Significant results (p-value < 0.05) are in bold.
Show moreShow less
Additional material (2)
Abstract

Agricultural intensification is a major driver of biodiversity loss. However, the exact impacts of such loss of key ecosystem service (ES) provisions in agriculture require more scrutiny. We assessed the population loss impacts of a key ES-providing species in an agricultural landscape. We tested the hypothesis that intensive agriculture causes density reduction of frogs known as pest regulators and that negatively affects pest regulation. Different frog densities (high and low) observed in a previous study in low and high agricultural intensification areas were used as treatments in a semi-controlled field experiment. Functional response of the frog species was also studied. Neither high nor low frog density had any significant effect on respective pest populations. Limited feeding rate of the frog species might have contributed to reduced predation pressure on pests. Surprisingly, at low density, frogs significantly reduced the arthropod natural enemies, probably their preferred prey due to the latter’s agility. Unlike in low frog density treatment, increased intraspecific competition at high frog density made them seek out pests as alternative prey to the preferred arthropod pest enemies, whose populations were not affected by frog density. This study for the first time links the population loss of frogs, a potential bio-control agent to ES provision in a multi-trophic system.

Keywords:
Frogs
Biological pest regulation
Ecosystem service and disservice
Intraguild predation
Trophic cascade effect
Full Text
Introduction

Extensive use of agrochemicals and habitat loss due to expansion of the agricultural area for intensification of food production have been identified as major drivers of biodiversity loss in agricultural landscapes (Newbold et al., 2015, Dudley and Alexander, 2017). Such biodiversity losses have led to the loss of key ecosystem services (ES) like pollination or natural pest regulation in agriculture (Brook et al., 2008). Sustaining the ES can be crucial in bridging the yield gap between the attainable and the actual yield (Bommarco et al., 2013). Biological control of crop pests delivered by pests’ natural predators is a key ES (Gurr et al., 2003). However, despite years of research, our understanding of the underlying ecological processes that govern biological control, particularly the complex food web interactions that regulate the process, is still incomplete (Bengtsson, 2015). For example, notwithstanding the benefits natural enemies of pests provide, their impact on other non-target arthropod populations has been highlighted (Howarth, 1991). On the other hand, the behavior and efficacy of natural predators in a multi-trophic and multi-species system are still mostly unclear and warrant scrutiny.

Amphibians have long been recognized as potential natural predators of crop pests (Hamer et al., 2004; Loman and Lardner, 2006) being both generalist and opportunist predators (Mahan and Johnson, 2007). Khatiwada et al. (2016) reported a high proportion of rice pests in frogs’ diet and recommended the introduction of frogs for biological pest control in rice fields. A sizeable abundance of amphibians is regarded as being efficient in bringing down the rice pest population (Teng et al., 2016; Fang et al., 2019). However, intensive agriculture has been found to negatively impact frog populations (Davies et al., 2008). Ghosh and Basu (2020) reported a reduced density of amphibians in areas of high agricultural intensification compared to less intensive traditional agricultural areas (for details on the basis for classifying agricultural intensification see, Ghosh and Basu, 2020). A decline in frog population in an intensive agricultural landscape can therefore be expected to negatively impact ES delivery. Apart from this, there could be other factors working in synergy with diversity loss that can disrupt the overall pest regulation by a given frog species. For example, in a multi-species and multi-trophic system within a guild, predation of the arthropod pest predators by super-predator frogs, referred to as “intraguild predation”, may also influence pest regulation (Chailleux et al., 2014).

In this study, we tested the hypothesis that population loss of ES-providing species led to ES reduction using an amphibian species as the model and pest regulation as the targeted ES. We build on the previously reported data on the density of amphibians by Ghosh and Basu (2020) as this reflects the real impact of agricultural intensification on amphibian density. We replicated the densities of a commonly found frog species in the rice fields (Fejervarya sp.) according to the density reported from areas of low and high agricultural intensification (Ghosh and Basu, 2020) in a semi-controlled field experiment to assess the impacts of two different frog densities on the rice pest population. Further, we assessed how trophic interactions between frogs (super-predator) and other natural insect pest enemies might impact the ES. Not only does this study bring out the impact of population losses of key ES providers on ES provision, it also highlights the importance of multi-trophic interaction on ES delivery which hitherto is unexplored especially in relation to amphibians as super-predators.

Materials and methodsExperiment site

The experiment was performed at the Agricultural Experimental Farm, University of Calcutta, Baruipur Campus in West Bengal, India in 2018 (22.3787°N, 88.4361°E). The farm is spread across an area of 85 ha with rice paddy being the major crop grown during the monsoon. Total experimental area was 7500 m2. The follow-up mesocosm experiment to study the nature of feeding response in the model amphibian was performed on the University campus in the year 2019.

Experimental design

Three blocks of 50 × 50 m area were selected for the study. Each of the blocks was divided into six experimental plots of 10 × 10 m where we installed our experimental units using a drift fence of length 10 m and height of 91.44 cm. We installed a total of 18 experimental plots across an area of 7500 m2 leaving a buffer of 5 m between each experimental plot. Based on the field realistic density of frogs in a previous study (Ghosh and Basu, 2020) we replicated the frog densities in this experiment. We randomly assorted the plots to three different frog treatments — one control without frogs, one with a treatment density of 10 frogs/100 m2 as observed in low agricultural intensification (we refer to this density as high frog density), and a third with a treatment of 5 frogs/100 m2 representing the density observed in high agricultural intensification zones (we refer to this density as low frog density).

Experimental model

We used adult (5–6 cm) paddy field frog species Fejervarya sp. as our model organism as they are the most abundant in the paddy fields of India and South Asia (Sumida et al., 2007). Fejervarya sp. is categorized as least concerned by IUCN and is also a generalist feeding mainly on insects (Chuang and Borzée, 2019) hence making it a potential bioresource to test their potential as pest regulators. This species belongs to a species complex, hence we have restrained to using only the genus name. We collected amphibians by active search in August 2018 from within the campus area for release in the experimental plots.

Experimental plot preparation

We selected paddy fields for their economic demand and for their potential in conserving amphibian diversity. For this study, we used a local rice cultivar “Patnai”, popularly grown in the study area. We prepared the seedbed with a density of 1916 seeds/m2 for a total seedbed area of 210 m2. Paddy saplings were transplanted around 30 days after preparing the seedbed at an average density of 762 plants/100 m2.

Each plot was secured using 3 ft high (from above the ground) drift fences on all four sides. Plots were secured with a mesh net of 4.8 cm gap size on the enclosure from the top as well as the sides of the fences (that allowed movements of arthropods) to exclude predation risk on both the pests and experimental frogs by birds as well as to prevent the escape of any frogs from the experimental enclosures by climbing on the fences (Fig. 1). Plots were left undisturbed for 15 days before the release of frogs. Frogs were released after sunset when the temperature was low to prevent desiccation. No pesticides, herbicides, or fertilizers were additionally applied. The study continued throughout the vegetative phase through the flowering stage from August to December.

Fig. 1.

Image of the experimental lay out, displaying the arrangement of the experimental units, showing the plot design with the drift fences, and mesh covering.

(0.68MB).
Arthropod sampling and identification

Sampling was divided into four phases- sample phase 1 (early vegetative phase), phase 2 (late vegetative phase), phase 3 (flowering phase) and phase 4 (fruiting phase). Insect sampling was conducted from 6.00 h to 13.00 h of the day. We used the sweep netting method with a standard-sized sweep net of 60.96 cm in length and an opening diameter of 30.48 cm (Fatahuddin et al., 2020). We used sweep netting owing to the accurate representation of the frog diet and the available prey as was also reported by Yousaf et al. (2010). Sweep netting was done 10 times along each of the four sides of the experimental plots leaving 3 m space on either side making a total of 40 sweeps per plot. All specimens collected were wet preserved in 70% alcohol for later identification and categorization. Table 1 provides a list of crop pests and arthropod pest predators along with their relative abundance. In most cases, we identified arthropods obtained in active sampling up to the genus level and categorized them into pests and non-pest natural predators (NICRA, 2011).

Table 1.

List of crop pest predators and pests observed in sweep net along with their relative abundance.

  Species  Family  Relative abundance 
Pest predators  Micraspis sp.  Coccinellidae  0.265 
  Pholcus sp.  Pholcidae  0.006 
  Identified till order  Odonata (order)  0.288 
  Tetragnatha sp.  Tetragnathidae  0.143 
  Oxyopidae sp.  Oxyopidae  0.001 
  Ophionea sp.  Carabidae  0.038 
  Identified till order  Mantidae  0.003 
  Lycosidae sp.  Lycosidae  0.001 
  Identified till order  Arachnida (class)  0.232 
  Anisoptera sp.  Dipterocarpaceae  0.012 
  Paederus sp.  Staphylinidae  0.005 
Crop pests  Identified till order  Formicidae  0.095 
  Altica sp.  Chrysomelidae  0.196 
  Thoracites sp.  Calliphoridae  0.003 
  Mythimna sp.  Noctuidae  0.001 
  Nilaparvata sp.  Delphacidae  0.003 
  Identified till order  Caterpillar of Lepidoptera  0.007 
  Identified till order  Coccinallidae  0.036 
  Identified till order  Coleoptera (order)  0.007 
  Identified till order  Diptera (order)  0.012 
  Phymatostetha sp.  Cercopoidea  0.030 
  Hieroglyphus sp.  Acrididae  0.062 
  Nephotettix sp.  Cicadellidae  0.052 
  Identified till order  Hemiptera (order)  0.008 
  Cnaphalocrocis sp.  Crambidae  0.018 
  Oebalus sp.  Pentatomidae  0.118 
  Identified till order  Phoridae  0.001 
  Identified till order  Delphacidae  0.006 
  Aulacophora sp.  Chrysomelidae  0.004 
  Pyrrhochoris sp.  Pyrrhocoridae  0.006 
  Dicladispa sp.  Chrysomelidae  0.001 
  Echinocnemus sp.  Curculionidae  0.002 
  Identified till order  Sciomyzidae  0.022 
  Scirpophaga sp.  Pyraustidae  0.011 
  Naranga sp.  Noctuidae  0.006 
  Oxya sp.  Acrididae  0.178 
  Pelopidas sp.  Hesperiidae  0.017 
  Spodoptera sp.  Noctuidae  0.019 
  Aspidimorpha sp.  Chrysomelidae  0.001 
  Leptocorisa sp.  Alydidae  0.043 
  Sogatella sp.  Delphacidae  0.034 
Estimating frog survival and diet analysis

From 17.00 h to 21.00 h every day we counted the frogs and checked for any dead or desiccated frogs inside the experimental plot to ensure that the number of frogs remained constant. We performed a diet analysis study following the stomach flush technique which allowed us to select prey for the functional response study (Solé et al., 2005).

Experiment to study the feeding potential of the model species

We performed a mesocosm experiment with one size class of frogs (5–6 cm, adult) and three prey densities of 8, 33, and 43. We used medium-sized grasshoppers (Gesonula sp., family Acrididae) in this experiment, as representative of crop pests, as grasshoppers were found to frequently constitute frogs’ diet during stomach flush. We used this species as representative pest species (the pest species we observed in sweep nets belonged to the same family, refer to Table 1) as they can be raised and maintained in the laboratory for scientific purpose hence was outsourced and at the same time Gesonula sp. are essentially defoliators of crops as well. The densities of grasshoppers used were replicates of the observed total available pest densities in control treatments (no frog treatments) across sampling sessions (observed pest densities per 10 × 10-m experimental plots were 8, 33, 43, and 44 in four consecutive sampling sessions). A total of 12 frogs participated in this experiment. Four frogs were randomly allotted to each of the three different pest density treatments. The experiment was conducted from September to November 2019. On the day of the experiment, one predator was transferred to the experimental mesocosm of 1 × 1 m that was designed to mimic the natural environment of these predators. Health of the mesocosm was properly maintained. Total experimental time was 17 h. At the end of the experiment, we removed the predator and counted the remaining prey. All predators had the same starvation period before they were included in an experiment. Since they feed on live prey, we also maintained a separate terrarium for the prey at varying densities to check the mortality rate within 24 h and found no death during the period.

To collect data on the handling time and attack rate of frogs, we replicated the experiment with varying prey densities and recorded the feeding rates for 30 min using a video camera (Nikon 5200D) fitted outside the terrarium.

Statistical analysisEffect of frog density and paddy growth phase on crop pest and predator build-up

A generalized Linear Mixed Model (GLMM) was used to predict changes in pest and natural predator abundance with varying phases of paddy growth (represented by sampling sessions) and frog treatments (Ghosh and Basu, 2020). Sampling phases and frog treatments were the fixed effects, and treatment within blocks was included as random effects in the model. We tested the significance of the two fixed effects using Analysis of Variance (ANOVA) Wald Chi-Square test and performed a posthoc test on the significant fixed effects. Further, for testing the significance of the random effects, we compared the models with and without random effects using ANOVA.

Functional response of predator

To understand how feeding rate might be influencing biological pest regulation we performed a diagnostic test for the functional response (Pritchard et al., 2017) with prey eaten as a function of prey provided. Based on the result, we developed a final model with attack rate and handling time.

All the analyses were performed with the software R (version 3.5.2) (R Core Team, 2018) using packages nlme (Pinheiro and Bates, 2023), Mass (Venables and Ripley, 2002), lmtest (Zeileis and Hothorn, 2002), multcomp (Hothorn et al., 2008), lme4 (Bates et al., 2015), car (Fox and Weisberg, 2019), and frair (Paterson et al., 2014).

ResultsEffects of frog abundance on crop pests and beneficial arthropod predators

Results from the GLMM showed that pest abundance was significantly affected by the paddy growth phases (Chi square = 261.590; df = 3; p < 0.001), especially by the 2nd (p < 0.001), 3rd (p < 0.001), and 4th (p < 0.001) phases. There was no significant effect of frog predation on pests at either high or low frog densities (Chi-square = 3.216; df = 2; p = 0.200; 10 frog treatment p = 0.369; 5 frog treatment p = 0.077). Arthropod pest predator abundance also was significantly affected by different paddy growth phases (Chi-square = 100.115; df = 3, p < 0.001), especially the 2nd (p < 0.001), 3rd (p < 0.001), and 4th (p = 0.006) phases. Interestingly, the abundance of beneficial arthropods was significantly affected by frog density (Chi-square = 6.442; df = 2; p = 0.039), especially by the lower density of frogs (5 frog treatment p = 0.016) and the effect was negative. However, abundance of beneficial arthropods was not affected by higher frog density (10 frog treatment p = 0.16). Table 2 provides a summary of the results we obtained from the ANOVA and GLMM. Tukey’s posthoc test suggests that 1st phase had lower values than the remaining phases, 2nd had lower values than the 3rd, and the 3rd and 4th did not differ for pests. For arthropod pest predators, 1st phase had lower values than the remaining phases, 2nd had lower values than the 4th, whereas 2nd and 3rd, and 3rd and 4th phases did not differ (Table 3). There was no significant difference between the models with random effects and without random effects (for pest abundance model: AIC value for model with random effects = 296, AIC value for model without random effects = 297; p = 0.08; for arthropod predator abundance: AIC value for model with random effects = 268, AIC value for model without random effects = 266, p = 0.35).

Table 2.

Summary of the results from ANOVA Wald Chi-square test for significance testing of fixed effects (showing Chi-square, df and p values) and GLMM results (showing estimates, standard error, and p-values) between levels of fixed effects and random effects. Significant results (p-value < 0.05) are in bold.

  Pest abundancePredator abundance
Fixed effects (Wald Chisq | Df |p-value)         
Sampling phase  261.59 3  <0.001  100.115 3  < 0.001 
Treatment  3.216 2  0.200  6.442 2  0.039 
Fixed effects (estimate ± SE | p-value)         
Intercept  2.225 ± 0.185  <0.001  2.733 ± 0.149  <0.001 
Second sampling phase  1.172 ± 0.114  <0.001  1.024 ± 0.142  <0.001 
Third sampling phase  1.706 ± 0.108  <0.001  0.881 ± 0.143  <0.001 
Fourth sampling phase  1.355 ± 0.112  <0.001  0.414 ± 0.151  <0.001 
Low density of frogs  0.276 ± 0.156  0.077  -0.377 ± 0.157  0.016 
High density of frogs  0.141 ± 0.157  0.369  −0.2184 ± 0.155  0.159 
Random effects (estimate ± SE | p-value)         
Treatment: Block  0.029 ± 0.172    0.01 ± 0.1   
Block  0.037 ± 0.194    0.00 ± 0.00   
Table 3.

Results of posthoc test for the fixed effects with estimates, standard error, and p-values. Significant results (p-value < 0.05) are in bold.

  Pest abundancePredator abundance
Fixed effects (estimate ± SE | p-value)         
Treatment         
High-control  0.289 ± 0.173  0.216  −0.366 ± 0.158  0.05 
Low-control  0.142 ± 0.173  0.690  −0.175 ± 0.155  0.494 
Low-high  -0.146 ± 0.172  0.670  0.190 ± 0.159  0.456 
Sampling phase         
2nd–1st  1.202 ± 0.181  <0.001  1.018 ± 0.147  <0.001 
3rd–1st  1.662 ± 0.178  <0.001  0.876 ± 0.148  <0.001 
4th–1st  1.358 ± 0.182  <0.001  0.553 ± 0.153  0.0015 
3rd–2nd  0.459 ± 0.157  0.01  −0.142 ± 0.129  0.688 
4th–2nd  0.155 ± 0.160  0.765  −0.464 ± 0.134  0.003 
4th–3rd  -0.304 ± 0.159  0.221  −0.322 ± 0.136  0.082 

To further explain this, we explored the pest-predator build-up across the three treatments (Fig. 2). From the graph it is clear that pest abundance decreases substantially in high frog density (i.e., low intensification) compared to the control plot without frogs, however, this difference is not evident between control plots and low frog density (i.e., in high intensification). Especially, in the third phase of sampling, graphical representation revealed a prominent variation in pest-predator build-up pattern between the treatments. During the third phase, the availability of pests per arthropod predator was 3.9 at lower frog density. In contrast, at higher frog density, the availability of pests per arthropod predator was 1.7 while for the control treatment, it was 1.1.

Fig. 2.

Crop pest and natural arthropod pest predator build-up in the three treatments (grey-control; green-high frog density/low AI; yellow-low frog density/high AI) across four sampling phases.

(0.26MB).
Functional response

Diagnostic tests indicated the feeding response pattern in our focal frog species to be Holling’s Type II. The final model with attack rate and handling time, and with the exponential co-efficient “q” fixed at 0, showed a significant effect on handling time (p < 0.0001). This validates that the functional response is of Hollings’s Type II (Fig. 3).

Fig. 3.

Functional response output showing Holling’s Type II feeding nature in focal frog species.

(0.06MB).
Discussion

This study clearly brings out the importance of maintaining a threshold density of beneficial organisms in a natural multi-trophic system, below which not only does the ES provider fail to provide the desirable service it may be also counter acting. This was suggested by the fact that amphibians at lower density (in high-intensive farming) would feed on beneficial arthropods reducing the overall pest regulation. Further, it shows how population loss of key ES providing biodiversity elements can translate into the loss of ecosystem service provisioning using amphibians as model bioregulators of rice pests.

In our study, although frogs consumed pests, their existing field densities (Ghosh and Basu, 2020) were insufficient for significant pest regulation in areas of both low and high agricultural intensification. Even the highest observed frog density (10 frogs/100 m2) in the low-intensive farmlands, appears to fall short of the threshold at which they can significantly regulate pests.

There can be three possible reasons why frogs failed to regulate pests. Firstly, the existing density of frogs even in the low intensification areas is not adequate to control the pests. Teng et al. (2016) reported that an experimental density of 150 frogs/100 m2 could effectively control pests. There is a wealth of literature that shows that traditional or low-intensity farmlands harbor more biodiversity and density than conventional intensive farming areas (Gonthier et al., 2014; Stein-Bachinger et al., 2022). However, whether this higher diversity and density of an ecosystem service provider is sufficient for effective service delivery has hardly been questioned. The current density of amphibians in low-intensive agriculture is evidently below the threshold that can be expected to provide a considerable service to the ecosystem. However, reduction of intensification pressure has been reported to increase frog density (Ghosh and Basu, 2020) and as the present study shows increased frog density does impact pest regulation due to the cascading effect (of increased insect predators) in the multi-trophic system.

The second possible reason for the unaltered pest population at both high and low frog densities (in low and high-intensive farming respectively) could be due to trophic interactions interfering with pest regulation. As the results show arthropod predators were regulated by the frogs at a low density due to intraguild predation, which could have reduced the predation pressure by arthropod predators on pests. In a multi-species and multi-trophic environment with a super-predator like frog, natural pest regulation may remain unrealized. It is well established that generalist predators may have an antagonistic effect on biological control (Perez-Alvarez et al., 2019) especially because of diffused predation pressure (Wells, 2007). However, if frogs feed on crop pests only, the growth rate of beneficial arthropods would decrease as they are robbed off their resource, and eventually, this could further lead to an increase in crop pest abundance. Again, if the harvesting is shifted completely towards arthropod predators there would be an increase in crop pest abundance due to the release of interspecific interaction from the pest predator population. Therefore, the presence of frogs would cause negative feedback on the crop pest predator population, consequently increasing crop pest abundance. Our study shows frogs exhibit a preference for arthropod predators over pests, especially at a lesser density (Table 2, Fig. 4). The diet of amphibians in croplands has been found to be dominated by non-pests (94% of diet) (Brown, 1974; Khatiwada et al. 2016). Lisiecki (2019) also highlights the absence of any significant effect of lizards on crop damage due to intraguild predation on the beneficial arthropods. Though generalist species have a wide range of prey choices, they often weaken the biocontrol efficiency as they kill and eat natural enemies as well (Snyder et al., 2022). Therefore, intraguild predation can be a cause for disrupting pest regulation in a multitrophic system.

Fig. 4.

Impact of prey preference shift in super-predator frog resulting in disrupted pest regulation. The preference of frogs for crop pests (black arrow) reduces the arthropod pest predator density due to resource depletion (broken arrow) lowering the predation pressure on crop pests by the pest predators (grey arrow) (A). The preference of frogs for the arthropod pest predators (black arrow) causes an increase in crop pests due to lack of predation pressure from the predators (arrow in grey) (B). Width of arrows signify feeding intensity.

(0.3MB).

The third possible factor that may affect pest regulation by frogs could be their feeding behavior. The model frog species seems to follow a Type II functional response pattern, however, this conclusion is based on a small number of replicates. This means that the feeding rate remains constant beyond a certain prey density and is dependent upon the prey handling time. Handling time (Th) relates to every prey item consumed and is indifferent to the prey density. Even with a higher density of grasshoppers, the maximum number of prey consumed was determined by the time spent in prey handling (T/Th). This means that the rate of prey consumption would decline at higher prey densities (that includes both crop pests as well as arthropod pest predators) due to handling constraints (Thorp et al., 2018).

The high prey availability per arthropod pest predator especially for high frog density compared to control can be explained by the following:

  • i

    Availability of pest refuge — at sampling phase 3 (30–40 days after transplanting) the crops have maximum vegetative mass (leaves and tillers) (http://www.agritech.tnau.ac.in/expert_system/paddy/Botany.html). The pest build-up was also the maximum at this phase. The pests found refuge in increased vegetative mass helping them to evade the frogs. Prey activity, crypsis, and refuge are important components that determine the outcome of prey-predator interactions (Sih and Moore, 1990; Chakraborty et al., 2017) and can destabilize the prey–predator cycle. Increase in search time also influences the profitability of the prey causing the predators to attack the more conspicuous prey, in our case the arthropod predators.

  • ii

    High predation pressure of frogs on the arthropod predators - our model species is an ambush forager (field observation) that forages using visual cues of prey movement (Cooper, 1999). Behaviourally predatory arthropods are more agile than the herbivore pests (Ahmed et al., 2016) who remain increasingly well hidden in the higher available vegetative mass and this makes the natural predator easier prey for the super-predator frog.

What still remains to be explained is why higher frog density did not have any significant impact on the arthropod pest predators compared to lower frog density. A possible reason could be increased resource competition among frogs at higher frog density that compels them to seek out other food options i.e., crop pests. However, more conclusive results would require further experimental scrutiny.

Conserving frogs in intensive agricultural landscapes for effective pest regulation

The vulnerability of amphibians to agricultural intensification is globally acknowledged (Ghosh and Basu, 2020). Agricultural intensification is only going to increase in the coming years more so in developing countries to feed the growing population. Our study brings out the need to conserve frogs in the intensive agricultural landscape for effective natural pest regulation service. Scientific management of appropriate landscape structures at the patch, and landscape levels, like semi-natural habitats, hedgerows, ephemeral water bodies, degree of connectivity with adjoining remnant natural vegetation, and habitat heterogeneity (Ghosh and Basu, 2020), will help in maintaining a healthy frog community. Vegetated buffer strips that include vegetation along ditches, and levees connecting various crop areas are important for maintaining ecosystem health (Marshall and Moonen, 2002), at the same time, they also provide habitat refugia for the amphibians. Further, the presence of vegetated strips provides continuous prey availability for amphibians (Martori et al., 2005). Such linear features and their degree of connectivity to adjoining natural or semi-natural vegetation would eventually increase the abundance of herpetofauna in general. The consequent environmental heterogeneity between natural and croplands would increase species co-existence among natural arthropods as well as the amphibian assemblage (Da Silva et al., 2012). Alongside bolstering the frog density and diversity, complex landscapes can weaken the strength of intraguild predatory effects by reducing niche overlap through spatial separation or even alternate resource availability (Perez-Alvarez et al., 2019).

As our study shows, despite the limitations, we discuss in this study, amphibians can still make good pest regulators provided they are maintained in a threshold abundance. Conservation of species or groups of species and maintaining their population density above specific threshold levels in an agricultural landscape at which they can be effective is the only desirable pathway towards developing an ecologically healthy sustainable agricultural landscape.

Authors’ contribution

DG and PB conceived the idea, DG worked on the methodology, and collected and analysed data; DG prepared the first draft and DG and PB wrote the final manuscript. Both authors approve the publication of the manuscript.

Declaration of competing interest

The authors report no declarations of interest.

Acknowledgment

We thank Prof. Santanu Jha and his team from Bidhan Chandra Krishi Viswavidyalaya for their help in insect and pest infestation identification. We acknowledge Arnab Chakraborty for his support during the fieldwork.

Appendix A
Supplementary data

The following are the supplementary data to this article:

References
[Ahmed et al., 2016]
K.S. Ahmed, M.Z. Majeed, A.A. Haidary, N. Haider.
Integrated pest management tactics and predatory coccinellids: a review.
J. Entomol. Zool. Stud., 4 (2016), pp. 591-600
[Bates et al., 2015]
D. Bates, M. Mächler, B. Bolker, S. Walker.
Fitting linear mixed-effects models using lme4.
J. Stat. Softw., 67 (2015), pp. 1-48
[Bengtsson, 2015]
J. Bengtsson.
Biological control as an ecosystem service: partitioning contributions of nature and human inputs to yield.
Ecol. Entomol., 40 (2015), pp. 45-55
[Bommarco et al., 2013]
R. Bommarco, D. Kleijn, S.G. Potts.
Ecological intensification: harnessing ecosystem services for food security.
Trends Ecol. Evol., 28 (2013), pp. 230-238
[Brook et al., 2008]
B. Brook, N. Sodhi, C. Bradshaw.
Synergies among extinction drivers under global change.
Trends Ecol. Evol., 23 (2008), pp. 453-460
[Brown, 1974]
R.L. Brown.
Diets and habitat preferences of selected anurans in Southeastern Arkansas.
Am. Midl. Nat., (1974), pp. 468-473
[Chakraborty et al., 2017]
S. Chakraborty, P.K. Tiwari, S.K. Sasmal, S. Biswas, S. Bhattacharya, J. Chattopadhyay.
Interactive effects of prey refuge and additional food for predator in a diffusive predator–prey system.
Appl. Math. Modell., 47 (2017), pp. 128-140
[Chailleux et al., 2014]
A. Chailleux, E. Mohl, A. Teixeira, G. Messelink, N. Desneux.
Natural enemy-mediated indirect interactions among prey species: current theories and potential for enhancing biocontrol services in agro-ecosystems.
Pest Manag. Sci., 70 (2014), pp. 1769-1779
[Chuang and Borzée, 2019]
M-F. Chuang, A. Borzée.
Fejervarya limnocharis (Asian rice frog).
Diet. Herpetol. Rev., 50 (2019), pp. 336
[Cooper, 1999]
W.E. Cooper.
Supplementation of phylogenetically correct data by two species comparison: support for correlated evolution of foraging mode and prey chemical discrimination in lizards extended by first intrageneric evidence.
Oikos, 86 (1999), pp. 97-104
[Davies et al., 2008]
B. Davies, J. Biggs, P. Williams, M. Whitfield, P. Nicolet, D. Sear, S. Bray, S. Maund.
Comparative biodiversity of aquatic habitats in the European agricultural landscape.
Agric. Ecosyst. Environ., 125 (2008), pp. 1-8
[Dudley and Alexander, 2017]
N. Dudley, S. Alexander.
Agriculture and biodiversity: a review.
Biodiversity, 18 (2017), pp. 45-49
[Fang et al., 2019]
K. Fang, X. Yi, W. Dai, H. Gao, L. Cao.
Effects of integrated rice-frog farming on paddy field greenhouse gas emissions.
Int. J. Environ. Res. Public Health., 16 (2019), pp. 1930
[Fatahuddin et al., 2020]
F. Fatahuddin, M. Melina, F. Fatmawaty, S. Sulaeha.
The population density of arthropods in the rice field ecosystem with insecticide application.
IOP Conf. Ser.: Earth Environ. Sci., 486 (2020),
[Fox and Weisberg, 2019]
John Fox, Sanford Weisberg.
An {R} Companion to Applied Regression.
[Ghosh and Basu, 2020]
D. Ghosh, P. Basu.
Factors influencing herpetofauna abundance and diversity in a tropical agricultural landscape mosaic.
Biotropica, 52 (2020), pp. 927-937
[Gonthier et al., 2014]
D.J. Gonthier, K.K. Ennis, S. Farinas, H.Y. Hsieh, A.L. Iverson, P. Batáry, J. Rudolphi, T. Tscharntke, B.J. Cardinale, I. Perfecto.
Biodiversity conservation in agriculture requires a multi-scale approach.
Proc. R. Soc. B: Biol. Sci., 281 (2014),
[Gurr et al., 2003]
G. Gurr, S. Wratten, J. Luna.
Multi-function agricultural biodiversity: pest management and other benefits.
Basic Appl. Ecol., 4 (2003),
[Hamer et al., 2004]
A.J. Hamer, J.A. Makings, S.J. Lane, M.J. Mahony.
Amphibian decline and fertilizers used on agricultural land in south-eastern Australia.
Agric. Ecosyst. Environ., 102 (2004), pp. 299-305
[Hothorn et al., 2008]
T. Hothorn, F. Bretz, P. Westfall.
Simultaneous inference in general parametric models.
Biometr. J., 50 (2008), pp. 346-363
[Howarth, 1991]
F.G. Howarth.
Environmental impacts of classical biological control.
Annu. Rev. Entomol., 36 (1991), pp. 489-509
[Khatiwada et al., 2016]
J.R. Khatiwada, S. Ghimirec, S.P. Khatiwad, B. Paudel, R. Bischo, J. Jianga, T. Haugaasen.
Frogs as potential biological control agents in the rice fields of Chitwan, Nepal.
Agric. Ecosyst. Environ., 230 (2016), pp. 307-314
[Lisiecki, 2019]
C. Lisiecki.
Efficacy of the Aegean Wall Lizard (Podarcis erhardii) as a Potential Biological Control Agent in Mediterranean Agroecosystems (Doctoral dissertation).
(2019),
[Loman and Lardner, 2006]
J. Loman, B. Lardner.
Does pond quality limit frogs Rana arvalis and Rana temporaria in agricultural landscapes? A field experiment.
J. Appl. Ecol., 43 (2006), pp. 690-700
[Mahan and Johnson, 2007]
R.D. Mahan, J.R. Johnson.
Diet of the gray treefrog (Hyla versicolor) in relation to foraging site location.
[Marshall and Moonen, 2002]
E.J.P. Marshall, A.C. Moonen.
Field margins in northern Europe: their functions and interactions with agriculture.
Agri. Ecosyst. Environ., 89 (2002), pp. 5-21
[Martori et al., 2005]
R. Martori, L. Aun, F.L. Gallego, C. Rozzi Giménez.
Temporal variation and size class distribution in a herpetological assemblage from Córdoba, Argentina.
Cuadernos de Herpetología, 19 (2005),
[Newbold et al., 2015]
T. Newbold, et al.
Global effects of land use on local terrestrial biodiversity.
Nature, 520 (2015), pp. 45-50
[NICRA team of rice pest surveillance, 2011]
NICRA team of rice pest surveillance.
Manual for Rice Pest Surveillance.
Jointly published by National Centre for Integrated Pest Management, New Delhi, Central Research Institute for Dryland Agriculture and Directorate of Rice Research, Hyderabad, (2011),
[Paterson et al., 2014]
R.A. Paterson, J.T.A. Dick, D.W. Pritchard, M. Ennis, M.J. Hatcher, A.M. Dunn.
Predicting invasive species impacts: community module functional response experiments reveal context dependencies.
J. Anim. Ecol., 84 (2014), pp. 453-463
[Perez-Alvarez et al., 2019]
R. Perez-Alvarez, B.A. Nault, K. Poveda.
Effectiveness of augmentative biological control depends on landscape context.
[Pinheiro and Bates, 2023]
J. Pinheiro, D. Bates.
R Core Team. nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1-162.
[Pritchard et al., 2017]
D. Pritchard, R. Paterson, H. Bovy, D. Barrios-O’Neill.
Frair: an R package for fitting and comparing consumer functional responses.
Methods Ecol. Evol., 8 (2017),
[Sumida et al., 2007]
M. Sumida, M. Kotaki, M.M. Islam, T.D. Hon, T. Igawa, Y. Kondo, M. Matsui, A. Silva, W. Khonsue, M. Nishioka.
Evolutionary relationships and reproductive isolating mechanisms in the rice frog (Fejervarya limnocharis) species complex from Sri Lnaka, Thailand, Taiwan and Japan, inferred from mtDNA gene sequences, allozymes, and crossing experiments.
Zool. Sci., 24 (2007), pp. 547-562
[Thorp et al., 2018]
C.J. Thorp, M.E. Alexander, J.R. Vonesh, J. Measey.
Size-dependent functional response of Xenopus laevis feeding on mosquito larvae.
PeerJ, 6 (2018), pp. e5813
[R Core Team, 2018]
R Core Team.
R: A Language and Environment for Statistical Computing.
[Sih and Moore, 1990]
A. Sih, R.D. Moore.
Interacting effects of predator and prey behavior in determining diets.
[Da Silva et al., 2012]
F.R. Da Silva, C.P Candeira, D. de Cerqueira Rossa-Feres.
Dependence of anuran diversity on environmental descriptors in farmland ponds.
Biodivers. Conserv., 21 (2012), pp. 1411-1424
[Solé et al., 2005]
M. Solé, O. Beckmann, B. Pelz, A. Kwet, W. Engels.
Stomach-flushing for diet analysis in anurans: an improved protocol evaluated in a case study in Araucaria forests, southern Brazil.
Stud. Neotrop. Fauna Environ., 40 (2005), pp. 23-28
[Snyder et al., 2022]
G.B. Snyder, O.M. Smith, E.G. Chapman, M.S. Crossley, D.W. Crowder, Z. Fu, J.D. Harwood, A.S. Jensen, K.L. Krey, C.A. Lynch, W.E. Snyder.
Alternative prey mediate intraguild predation in the open field.
Pest Manag. Sci., 78 (2022), pp. 3939-3946
[Stein-Bachinger et al., 2022]
K. Stein-Bachinger, S. Preißel, S. Kühne, M. Reckling.
More diverse but less intensive farming enhances biodiversity.
Trends Ecol. Evol., 37 (2022), pp. 395-396
[Teng et al., 2016]
Q. Teng, X. Hu, F. Luo, C. Cheng, X. Ge, M. Yang, L. Liu.
Influences of introducing frogs in the paddy fields on soil properties and rice growth.
J. Soils Sediments, 16 (2016), pp. 51-61
[Venables and Ripley, 2002]
W.N. Venables, B.D. Ripley.
Modern Applied Statistics with S.
fourth edition, Springer, (2002),
[Wells, 2007]
K. Wells.
The Ecology and Behavior of Amphibians.
Bibliovault OAI Repository, The University of Chicago Press, (2007), http://dx.doi.org/10.7208/chicago/9780226893334.001.0001
[Yousaf et al., 2010]
S. Yousaf, T.M. Qureshi, M. Rias, I. Zia.
Population variation and food habits of ranid frogs in rice the rice-based cropping system n Gujranwala region in Pakistan.
Asian Herpetol. Res., 1 (2010), pp. 122-130
[Zeileis and Hothorn, 2002]
A. Zeileis, T. Hothorn.
Diagnostic checking in regression relationships.
R News, 2 (2002), pp. 7-10
Perspectives in Ecology and Conservation
Article options
Tools
Supplemental materials