Does mutual interference stabilize prey-predator model with Bazykin-Crowley-Martin trophic function?

We investigated a system of ordinary differential equations that describes the dynamics of prey and predator populations, taking into account the Allee effect affecting the reproduction of the predator population, and mutual interference amongst predators, which is modeled with the Bazykin-Crowley-Martin (BCM) trophic function. Bifurcation analysis revealed a rich spectrum of bifurcations occurring in the system. In particular, analytical conditions for the saddle-node, Hopf, cusp, and Bogdanov-Takens bifurcations were derived for the model parameters, quantifying the strength of the predator interference, the Allee effect, and the predation efficiency. Numerical simulations verify and illustrate the analytical findings. The main purpose of the study was to test whether the mutual interference in the model with BCM trophic function provides a stabilizing or destabilizing effect on the system dynamics. The obtained results suggest that the model demonstrates qualitatively the same pattern concerning varying the interference strength as other predator-dependent models: both low and very high interference levels increase the risk of predator extinction, while moderate interference has a favorable effect on the stability and resilience of the prey-predator system.

Predator-Prey Interaction Between Xylocoris sordidus (Hemiptera: Anthocoridae) and Enneothrips enigmaticus (Thysanoptera: Thripidae)

The peanut thrips, Enneothrips enigmaticus (Thysanoptera: Thrypidae), is an important pest of the peanut (Arachis hypogaea) in South America. Due to concerns about the environment and human health induced by the extensive use of pesticides in the management control of pests, environmentally and friendlier tactics must be targeted. Thus, this study investigates, for the first time, the behavior of Xylocoris sordidus (Hemiptera: Anthocoridae) as a biological control agent for E. enigmaticus. The methodology included no-choice tests to assess whether the predation rate varies according to the developmental stage of the prey, as well as the predator's developmental stage with the highest predation capacity. Additionally, an analysis of the functional response of adult and 5th instar nymphs of X. sordidus exposed to different densities of E. enigmaticus nymphs (1, 2, 4, 8, 16, and 32) was conducted. The results confirm the predation of peanut thrips by X. sordidus, with a higher predation rate in the nymphal stages of the prey. There was no difference in predation capacity between predator nymphs and adults, and exhibiting a type II functional response. Therefore, the potential of X. sordidus as a biological control agent for E. enigmaticus is confirmed, showing the importance of adopting measures to preserve this predator in peanut crops.

Assessing the potential phytosanitary threat of the house cricket Acheta domesticus

Phytosanitary threats can pose substantial risks to global agriculture and ecological systems, affecting biodiversity, human well-being, and food security. Meanwhile, global warming is projected to exacerbate these threats in the future. One in Europe already widely distributed potential phytosanitary threat that may benefit from global warming is the house cricket Acheta domesticus. This study explored the potential of A. domesticus as a relevant non-native phytosanitary threat under changing climatic conditions by conducting a series of functional response experiments across a temperature gradient (20, 25, and 30 °C). Acheta domesticus exhibited comparable patterns of seed consumption and functional responses. Seed type (millet seeds, wheat grains) and temperature increase influenced the damage inflicted on seeds, with softer and smaller seeds being more susceptible to damage, further amplified by warmer temperatures. The study's outcomes underline the phytosanitary threat that A. domesticus may pose. Considering the species' established presence and adaptable nature in urban environments exacerbates the potential for A. domesticus to transition to rural and agricultural areas. Its increasing production as a food item, paired with the here-identified potential to damage seeds, emphasizes the need for proactive and science-based strategies to address emerging phytosanitary threats driven by non-native species under changing climatic conditions. As global temperatures continue to rise, the assessment and management of potential pest species like A. domesticus will be crucial for safeguarding agriculture productivity and ecological balance.

The inappropriate application of imidacloprid destroys the ability of predatory natural enemies to control pests in the food chain: A case study of the feeding behavior of Orius similis on Frankliniella occidentalis

Insecticides are an indispensable and important tool for agricultural production. However, the inappropriate application of insecticides can cause damage to the food chain and ecosystem. Orius similis is an important predatory and natural enemy of Frankliniella occidentalis. Imidacloprid is widely used to control pests, but will inevitably exert adverse effects on O. similis. In order to determine the effect of different imidacloprid treatments on the ability of O. similis to prey on the 2nd-instar nymphs of F. occidentalis, we determined the toxicity and predation of imidacloprid on different stages of O. similis under contact and ingestion treatments. In addition, we used the Holling disc equation to evaluate the ability of O. similis to search and exhibit predatory activity following contact and ingestion treatments. Analysis showed that the highest LC10 and LC20 values for imidacloprid contact and ingestion toxicity treatment were 17.06 mg/L and 23.74 mg/L, respectively. Both imidacloprid treatments led toa reduction in the predatory of O. similis on prey. The functional responses of the 3rd to 5th instar nymphs, along with female and male O. similis adults to the 2nd-instar nymphs of F. occidentalis were consistent with the Holling type II response following contact and ingestion with imidacloprid. However, following imidacloprid treatment, the handing time (Th) of O. similis with single F. occidentalis was prolonged and the instantaneous attack rate (a) was reduced after imidacloprid treatment. The predatory capacity (a/Th) of female O. similis adults when treated with the LC10 concentration of imidacloprid by ingestion was 52.85; this was lower than that of the LC10 concentration of imidacloprid in the contact treatment (57.67). The extent of predation of O. similis on the 2nd-instar nymphs of F. occidentalis was positively correlated with prey density, although the search effect was negatively correlated with prey density. The most extensive search effect was exhibited by adult O. similis females. Simulations with the Hessell-Varley interference model showed that an increase in the number of O. similis would reduce search efficiency regardless of whether they were treated with imidacloprid or not. Thus, O. similis, especially female adults, exhibited strong potential for controlling the 2nd-instar nymphs of F. occidentalis. The toxicity of ingestion following treatment with the same concentration of imidacloprid in O. similis was greater than that of contact treatment. When using O. similis to control F. occidentalis in the field, we should increase the number of female adults released, and prolong the interval between imidacloprid treatment and O. similis exposure. This strategy will improve the control ability of O. similis, coordinate both chemical and biological control, reduce the impact of pesticides on the environment, and improve the efficiency of agricultural production.

Functional response of Neoseiulus californicus (Acari: Phytoseiidae) to Tetranychus urticae (Acari: Tetranychidae) at different temperatures

Environmental factors like temperature have a great impact on the predation potential of biological control agents. In the present study, the functional response of the predatory mite Neoseiulus californicus (Acari: Phytoseiidae) to the pest mite Tetranychus urticae (Acari: Tetranychidae) at moderate to high temperatures under laboratory conditions was determined. The study aimed to understand the prey-predator interaction under different temperatures and prey densities. Five constant temperatures (24 °C, 27 °C, 30 °C, 33 °C, and 36 °C), and thirteen prey densities (4, 5, 8, 10, 12, 15, 16, 20, 24, 25, 30, 32, and 40) of each stage (adult, nymph, larvae, and egg stage) were employed in the experiment. Observations were made 24 h after the start of each experiment. Results revealed that the predatory mites showed type II functional response to adult females of T. urticae, whereas type I to other stages (nymphs, larvae, and eggs) of T. urticae. The predation capability of adult predatory mites on T. urticae was significant at 24-36 °C. The instantaneous attack rate (a) of N. californicus increased and the handling time (Th) decreased with an increase in temperature. The maximum attack rate was recorded at 36 °C (1.28) for the egg stage. The longest handling time was (0.78) for the larval stage of T. urticae at 30 °C. Daily consumption increased with increasing prey density. Maximum daily consumption was observed at 33 °C (30.00) at the prey density of 40. Searching efficiency decreased with the increase in prey density but was found to increase with the rise in temperature. N. californicus was found to be voracious on the larval and egg stages. Conclusively, the incorporation of N. californicus at earlier stages (larvae and eggs) of T. urticae would be beneficial under warm conditions because managing a pest at its initial stage will save the crop from major losses. The results presented in this study at various temperatures will be helpful in different areas with different temperature extremes. The results of the functional response can also be applied to mass rearing, quality testing, and integrated pest management programmes.

No sensitivity to functional forms in the Rosenzweig-MacArthur model with strong environmental stochasticity

The classic Rosenzweig-MacArthur predator-prey model has been shown to exhibit, like other coupled nonlinear ordinary differential equations (ODEs) from ecology, worrying sensitivity to model structure. This sensitivity manifests as markedly different community dynamics arising from saturating functional responses with nearly identical shapes but different mathematical expressions. Using a stochastic differential equation (SDE) version of the Rosenzweig-MacArthur model with the three functional responses considered by Fussmann & Blasius (2005), I show that such sensitivity seems to be solely a property of ODEs or stochastic systems with weak noise. SDEs with strong environmental noise have by contrast very similar fluctuations patterns, irrespective of the mathematical formula used. Although eigenvalues of linearized predator-prey models have been used as an argument for structural sensitivity, they can also be an argument against structural sensitivity. While the sign of the eigenvalues' real part is sensitive to model structure, its magnitude and the presence of imaginary parts are not, which suggests noise-driven oscillations for a broad range of carrying capacities. I then discuss multiple other ways to evaluate structural sensitivity in a stochastic setting, for predator-prey or other ecological systems.

Bacillus thuringiensis-based bioinsecticides affect predation of Euborellia annulipes on diamondback moth larvae

Interactions between earwigs and entomopathogens, such as Bacillus thuringiensis (Bt), are still poorly understood. This study tested whether Bt-based bioinsecticides have any effect on the predation of Euborellia annulipes (Lucas) (Dermaptera: Anisolabididae) on Plutella xylostella (L.) (Lepidoptera: Plutellidae), one of the pests with the largest number of cases of use and resistance to Bt. Fourth instar larvae were Bt infected by feeding on collard green leaves treated with Dipel®WG and XenTari®WG at the manufacturer-recommended doses. We used one no-choice condition, in which the predator had access to uninfected or Bt-infected larvae separately, and four free-choice conditions: uninfected vs Dipel®-infected larvae, uninfected vs XenTari®-infected larvae, Dipel®-infected vs XenTari®-infected larvae, and uninfected vs Bt-infected larvae with both bioinsecticides. Uninfected larvae were less consumed than those infected by both Bt-bioinsecticides in the no-choice condition. There was a higher consumption of uninfected over Dipel®-infected larvae in the free-choice condition. Overall, uninfected larvae were preferred over both Bt-based bioinsecticides infected larvae. We also used six different prey densities. The ringlegged earwig's predation rate enhanced as the prey population density increased, but the functional response was not affected by Bt-infection, being type II. The predator invested a low amount of handling time on Bt-fed prey and increased the maximum predation rate. Bt-based bioinsecticides cause effects on E. annulipes predation by altering their feeding preference and some aspects of its predatory behavior. The results of our study provide an important background for understanding interactions between earwigs and Bt. In addition, they can be used for decision making during approaches to integrated P. xylostella management.

From the distributions of times of interactions to preys and predators dynamical systems

We consider a stochastic individual based model where each predator searches and then manipulates its prey or rests during random times. The time distributions may be non-exponential and density dependent. An age structure allows to describe these interactions and get a Markovian setting. The process is characterized by a measure-valued stochastic differential equation. We prove averaging results in this infinite dimensional setting and get the convergence of the slow-fast macroscopic prey predator process to a two dimensional dynamical system. We recover classical functional responses. We also get new forms arising in particular when births and deaths of predators are affected by the lack of food.

Resilient amphipods: Gammarid predatory behaviour is unaffected by microplastic exposure and deoxygenation

Microplastics are a ubiquitous and persistent form of pollution globally, with impacts cascading from the cellular to ecosystem level. However, there is a paucity in understanding interactions between microplastic pollution with other environmental stressors, and how these could affect ecological functions and services. Freshwater ecosystems are subject to microplastic input from anthropogenic activities (eg. wastewater), but are also simultaneously exposed to many other stressors, particularly reduced dissolved oxygen availability associated with climatic warming and pollutants, as well as biological invasions. Here, we employ the comparative functional response method (CFR; quantifying and comparing organism resource use as a function of resource density) to investigate the relative impact of different microplastic concentrations and oxygen regimes on predatory trophic interactions of a native and an invasive alien gammarid (Gammarus duebeni and Gammarus pulex). No significant effect on trophic interaction strengths was found from very high concentrations of microplastics (200 mp/L and 200,000 mp/L) or low oxygen (40 %) stressors on either species. Additionally, both gammarid species exhibited significant Type II functional responses, with attack rates and handling times not significantly affected by microplastics, oxygen or gammarid invasion status. Thus, both species showed resistance to the simultaneous effects of microplastics and deoxygenation in terms of feeding behaviour. Based on these findings, we suggest that the trophic function, in terms of predation rate, of Gammarus spp. may be sustained under acute bouts of microplastic pollution even in poorly‑oxygenated waters. This is the first study to investigate microplastic and deoxygenation interactions and to find no evidence for an interaction on a key invertebrate ecosystem service. We argue that our CFR methods can help understand and predict the future ecological ramifications of microplastics and other stressors across taxa and habitats.

Parasite infection but not chronic microplastic exposure reduces the feeding rate in a freshwater fish

Microplastics (plastics <5 mm) are an environmental contaminant that can negatively impact the behaviour and physiology of aquatic biota. Although parasite infection can also alter the behaviour and physiology of their hosts, few studies have investigated how microplastic and parasite exposure interact to affect hosts. Accordingly, an interaction experiment tested how exposure to environmentally relevant microplastic concentrations and the trophically transmitted parasite Pomphorhynchus tereticollis affected the parasite load, condition metrics and feeding rate of the freshwater fish final host chub Squalius cephalus. Microplastic exposure was predicted to increase infection susceptibility, resulting in increased parasite loads, whereas parasite and microplastic exposure were expected to synergistically and negatively impact condition indices and feeding rates. Following chronic (≈170 day) dietary microplastic exposure, fish were exposed to a given number of gammarids (4/8/12/16/20), with half of the fish presented with parasite infected individuals, before a comparative functional response experiment tested differences in feeding rates on different live prey densities. Contrary to predictions, dietary microplastic exposure did not affect parasite abundance at different levels of parasite exposure, specific growth rate was the only condition index that was lower for exposed but unexposed fish, with no single or interactive effects of microplastic exposure detected. However, parasite infected fish had significantly lower feeding rates than unexposed fish in the functional response experiment, with exposed but unexposed fish also showing an intermediate decrease in feeding rates. Thus, the effects of parasitism on individuals were considerably stronger than microplastic exposure, with no evidence of interactive effects. Impacts of environmentally relevant microplastic levels might thus be relatively minor versus other stressors, with their interactive effects difficult to predict based on their single effects.

When worlds collide: Invader-driven benthic habitat complexity alters predatory impacts of invasive and native predatory fishes

Interactions between multiple invasive alien species (IAS) might increase their ecological impacts, yet relatively few studies have attempted to quantify the effects of facilitative interactions on the success and impact of aquatic IAS. Further, the effect of abiotic factors, such as habitat structure, have lacked consideration in ecological impact prediction for many high-profile IAS, with most data acquired through simplified assessments that do not account for real environmental complexities. In the present study, we assessed a potential facilitative interaction between a predatory invasive fish, the Ponto-Caspian round goby (Neogobius melanostomus), and an invasive bivalve, the Asian clam (Corbicula fluminea). We compared N. melanostomus functional responses (feeding-rates under different prey densities) to a co-occurring endangered European native analogue fish, the bullhead (Cottus gobio), in the presence of increased levels of habitat complexity driven by the accumulation of dead C. fluminea biomass that persists within the environment (i.e. 0, 10, 20 empty bivalve shells). Habitat complexity significantly influenced predation, with consumption in the absence of shells being greater than where 10 or 20 shells were present. However, at the highest shell density, invasive N. melanostomus maximum feeding-rates and functional response ratios were substantially higher than those of native C. gobio. Further, the Relative Impact Potential metric, by combining per capita effects and population abundances, indicated that higher shell densities exacerbate the relative impact of the invader. It therefore appears that N. melanostomus can better tolerate higher IAS shell abundances when foraging at high prey densities, suggesting the occurrence of an important facilitative interaction. Our data are thus fully congruent with field data that link establishment success of N. melanostomus with the presence of C. fluminea. Overall, we show that invader-driven benthic habitat complexity can alter the feeding-rates and thus impacts of predatory fishes, and highlight the importance of inclusion of abiotic factors in impact prediction assessments for IAS.

Incorporating nonlinearity with generalized functional responses to simulate multiple predator effects

Predicting the combined effects of predators on shared prey has long been a focus of community ecology, yet quantitative predictions often fail. Failure to account for nonlinearity is one reason for this. Moreover, prey depletion in multiple predator effects (MPE) studies generates biased predictions in applications of common experimental and quantitative frameworks. Here, we explore additional sources of bias stemming from nonlinearities in prey predation risk. We show that in order to avoid bias, predictions about the combined effects of independent predators must account for nonlinear size-dependent risk for prey as well as changes in prey risk driven by nonlinear predator functional responses and depletion. Historical failure to account for biases introduced by well-known nonlinear processes that affect predation risk suggest that we may need to reevaluate the general conclusions that have been drawn about the ubiquity of emergent MPEs over the past three decades.

Effects of prey density and flow speed on plankton feeding by garden eels: A flume study

Feeding by zooplanktivorous fish depends on their foraging movements and the flux of prey to which they are exposed. While prey flux is a linear function of zooplankton density and flow speed, those two factors are expected to contribute differently to fish movements. Our objective was to determine the effects of these factors for garden eels, stationary fish that feed while anchored to the sandy bottom by keeping the posterior parts of their bodies inside a burrow. Using a custom-made flume with a sandy bottom, we quantified the effects of prey density and flow speed on feeding rates by spotted garden eels (Heteroconger hassi). Feeding rates increased linearly with prey density. However, feeding rates did not show a linear relationship with flow speed and decreased at 0.25 m s⁻¹. Using label-free tracking of body points and 3D movement analysis, we found that the reduction in feeding rates was related to modulation of the eel's movements, whereby the expected increase in energy expenditure was avoided by reducing exposure and drag. No effects of flow speed on strike speed, reactive distance or vectorial dynamic body acceleration (VeDBA) were found. A foraging model based on the body length extended from the burrow showed correspondence with observations. These findings suggest that as a result of their unique foraging mode, garden eels can occupy self-made burrows in exposed shelter-free sandy bottoms where they can effectively feed on drifting zooplankton.

Summary: Analysis of feeding by anchored garden eels in a flume shows their unique strategies, modulating dependence on burrows and foraging movements, to enable effective feeding at a wide range of flow speeds.

Plasticity of Foraging Strategies Adopted by the Painted Ghost Crab, Ocypode gaudichaudii, in Response to in situ Food Resource Manipulation Experiments

The feeding strategies of Ocypode gaudichaudii at two sandy beaches, Culebra Beach (CB) and Playa Venao (PV) in Panama, were studied via three experiments. Two separate manipulative in situ experiments were conducted to determine how the densities of food resources and the size of the supplemented food offered to the crabs can affect their diet and food handling behavior. The third experiment, a transplantation study, was also conducted to determine the plasticity of the feeding behavior of the displaced crabs. In the first experiment, freshly-emerged crabs showed different feeding modes when washed-sediment was seeded with different densities of diatoms and rove beetles, which suggests that they are optimal foragers. Crabs hoarded food in the second experiment when food augmentation was performed, in which small and large food pellets were placed around the burrows at the beginning and end of the crabs' feeding cycle. All freshly-emerged crabs from both sites foraged on the small pellets outside their burrows and did not cache food; when pellets were provided at the end of the feeding cycle, crabs from CB fed on some of the small pellets and returned to their burrows with the uneaten pellets left on the surface, whereas crabs at PV picked up all the small food pellets and transferred them into their burrows over several trips before plugging their burrow entrances. Only the crabs from PV carried the large food pellets supplemented at the start and end of the feeding cycle into their burrows. In contrast, the crabs at CB often left behind the partially-eaten pellets on the sand surface, probably due to the increased risk of predation associated with the prolonged handling time of the large food pellets. Excavation of the burrows of the crabs that hoarded food showed that all the pellets were deposited at the bend of the burrows, indicating that they were not consumed immediately. Crabs that fed in droves at PV stopped droving and foraged around their burrows after being transplanted to CB. This is the first documentation of food hoarding in a sandy beach macroinvertebrate at a resource-impoverished habitat. The plasticity of feeding strategies adopted by the painted ghost crab in response to different densities of food resources in the habitat could be an adaptation to the dynamic sandy beach environment.

Temperature, not salinity, drives impact of an emerging invasive species

Biological invasions are a growing ecological and socioeconomic problem worldwide. While robust predictions of impactful future invaders are urgently needed, understandings of invader impacts have been challenged by context-dependencies. In aquatic systems in particular, future climate change could alter the impacts of invasive non-native species. Widespread warming coupled with sea freshening may exacerbate ecological impacts of invaders in marine environments, compromising ecosystem structure, function and stability. We examined how multiple abiotic changes affect the potential ecological impact of an emerging invasive non-native species from the Ponto-Caspian region - a notorious origin hotspot for invaders, characterised by high salinity and temperature variation. Using a comparative functional response (feeding rates across prey densities) approach, the potential ecological impacts of the gammarid Pontogammarus maeoticus towards native chironomid prey were examined across a range of current and future temperature (18, 22 °C) and salinity (14, 10, 6, 2 ppt) regimes in a factorial design. Feeding rates of P. maeoticus on prey significantly increased with temperature (by 60%), but were not significantly affected by salinity regime. Gammarids displayed significant Type II functional responses, with attack rates not significantly affected by warming across all salinities. Handling times were, however, shortened by warming, and thus maximum feeding rates significantly increased, irrespective of salinity regime. Functional responses were significantly different following warming at high prey densities under all salinities, except under the ambient 10 ppt. Euryhalinity of invasive non-native species from the Ponto-Caspian region thus could allow sustained ecological impacts across a range of salinity regimes. These results corroborate high invasion success and field impacts of Ponto-Caspian gammarids in brackish through to freshwater ecosystems. Climate warming will likely worsen the potential ecological impact of P. maeoticus. With invasions growing worldwide, quantifications of how combined elements of climate change will alter the impacts of emerging invasive non-native species are needed.

Microplastics do not affect the feeding rates of a marine predator

Microplastics may affect the physiology, behaviour and populations of aquatic and terrestrial fauna through many mechanisms, such as direct consumption and sensory disruption. However, the majority of experimental studies have employed questionably high dosages of microplastics that have little environmental relevance. Predation, in particular, is a key trophic interaction that structures populations and communities and influences ecosystem functioning, but rarely features in microplastic research. Here, we quantify the effects of low (~65-114 MP/L) and high (~650-1140 MP/L) microplastic concentrations on the feeding behaviour of a ubiquitous and globally representative key marine predator, the shore crab, Carcinus maenas. We used a functional response approach (predator consumption across prey densities) to determine crab consumption rates towards a key marine community prey species, the blue mussel Mytilus edulis, under low and high microplastic concentrations with acute (8h) and chronic (120h) microplastic exposure times. For both the acute and chronic microplastic exposure experiments, proportional prey consumption by crabs did not differ with respect to microplastic concentration, but significantly decreased over increasing prey densities. The crabs thus displayed classical, hyperbolic Type II functional responses in all experimental groups, characterised by high consumption rates at low prey densities. Crab attack rates, handling times and maximum feeding rates (i.e. functional response curves) were not significantly altered under lower or higher microplastics concentrations, or by acute or chronic microplastic exposures. Here, we show that functional response analyses could be widely employed to ascertain microplastic impacts on consumer-resource interactions. Furthermore, we suggest that future studies should adopt both acute and chronic microplastic exposure regimes, using environmentally-relevant microplastic dosages and types as well as elevated future scenarios of microplastic concentrations.

Foraging behaviour of Scymnus coccivora Ayyar against cotton mealybug Phenacoccus solenopsis Tinsley

Predation is one of the significant biotic mortality factors reducing the insect pest population as functional response and the numerical response of the predator are the key factor regulating the population dynamics of predator prey species. This study is aimed to evaluate the functional response of all the developmental stages of Scymnus coccivora Ayyar (Coleoptera: Coccinellidae) against the different densities of cotton mealybug, Phenacoccus solenopsis (Tinsley) (Hemiptera: Pseudococcidae) and the numerical response of female predator. Experiments were carried out in controlled environment laboratory conditions at 25 ± 1 °C temperature, 60 ± 5% relative humidity and photoperiod of 16 h. Number of eggs consumed, number of eggs laid and the Efficiency of Conversion of Ingested food (ECI) were recorded daily. Results from the study revealed that all the developmental stages of S. coccivora exhibited a Type II response. Different parameters such as attack rate (a'), handling time (T_h) and the maximum rate of predation were estimated using Roger's random attack equation and Holling Disc equation in which Rogers random attack equation was found best fit. Female has shown the highest attack rate (a') followed by IV^th instar grub, male, III^rd, II^nd and I^st instar grub. With low handling time, IV^th instar grub has shown maximum predation rate of 76.40 per day followed by female (75.86), male (58.79), III^rd (22.84), II^nd (19.65) and I^st instar grub (15.39). The numerical response increase was curvilinearly related to different prey densities with the highest number of eggs (11.8 ± 3.44) produced at highest prey density (160). The Efficiency of Conversion of Ingested food (ECI) was highest (64.49 ± 8.03) at prey density of 10. Understanding the factors that lead to variation in functional response of predator in natural population will advance our understanding of the effects of predation on individual and the effectiveness of coccinellid predators as biocontrol agent against cotton mealybug.

Fitting stochastic predator-prey models using both population density and kill rate data

Most mechanistic predator-prey modelling has involved either parameterization from process rate data or inverse modelling. Here, we take a median road: we aim at identifying the potential benefits of combining datasets, when both population growth and predation processes are viewed as stochastic. We fit a discrete-time, stochastic predator-prey model of the Leslie type to simulated time series of densities and kill rate data. Our model has both environmental stochasticity in the growth rates and interaction stochasticity, i.e., a stochastic functional response. We examine what the kill rate data brings to the quality of the estimates, and whether estimation is possible (for various time series lengths) solely with time series of population counts or biomass data. Both Bayesian and frequentist estimation are performed, providing multiple ways to check model identifiability. The Fisher Information Matrix suggests that models with and without kill rate data are all identifiable, although correlations remain between parameters that belong to the same functional form. However, our results show that if the attractor is a fixed point in the absence of stochasticity, identifying parameters in practice requires kill rate data as a complement to the time series of population densities, due to the relatively flat likelihood. Only noisy limit cycle attractors can be identified directly from population count data (as in inverse modelling), although even in this case, adding kill rate data - including in small amounts - can make the estimates much more precise. Overall, we show that under process stochasticity in interaction rates, interaction data might be essential to obtain identifiable dynamical models for multiple species. These results may extend to other biotic interactions than predation, for which similar models combining interaction rates and population counts could be developed.

Can contamination by major systemic insecticides affect the voracity of the harlequin ladybird?

Systemic neurotoxic insecticides are widely used to control aphid pests worldwide and their potential non-target effects on aphid predators are often unknown. Behavioral responses linked to biological control services are crucial when assessing the compatibility of chemicals with biocontrol organisms. This is particularly relevant for insecticides at low and sublethal concentrations. We studied the acute toxicity and the sublethal effect on the voracity of the generalist predator Harmonia axyridis (Coleoptera: Coccinellidae) caused by the exposure to three systemic insecticides routinely used against aphids. The tested insecticide concentrations were the Lethal Concentration 50% (LC₅₀), 20% (LC₂₀) and 1% (LC₁) estimated for the target pest Aphis gossypii (Hemiptera: Aphididae) in a companion study. The survival and the voracity differed among the tested chemicals and concentrations, but only thiamethoxam at LC₅₀ caused a significant predator mortality, and individuals that survived showed a reduced predation rate. The predators showed a density independent functional response after the exposure to most of the insecticide-concentration combinations, while an inverse density dependence of the prey consumption rate was observed for coccinellids exposed to sulfoxaflor and thiamethoxam at their lowest tested concentration. The estimated parameters, i.e., the attack rate and the prey handling time, were affected at higher concentrations by both imidacloprid and sulfoxaflor. These findings stress the importance of carefully evaluating side effects of insecticides at very low concentrations on beneficial arthropods in the risk assessment schemes for sustainable pest control programmes.

Predatory functional response and fitness parameters of Orius strigicollis Poppius when fed Bemisia tabaci and Trialeurodes vaporariorum as determined by age-stage, two-sex life table

Background

The polyphagous predatory bug O. strigicollis is an active predator used to control thrips and aphids. The whitefly species Bemisia tabaci and Trialeurodes vaporariorum are voracious pests of different economic agricultural crops and vegetables.

Method

In this study, the Holling disc equation and the age-stage, two-sex life table technique were used to investigate the functional response and biological traits of third instar nymphs and adult female O. strigicollis when presented third instar nymphs of both whitefly species as prey.

Results

The results showed a type II functional response for each life stage of O. strigicollis when fed each whitefly species. The calculated prey handling time for different O. strigicollis life stages were shorter when fed T. vaporariorum than when fed B. tabaci nymphs. In contrast, the nymphal development of O. strigicollis was significantly shorter when fed B. tabaci than T. vaporariorum nymphs. Additionally, the total pre-oviposition period of adult females was statistically shorter when fed B. tabaci nymphs than T. vaporariorum nymphs. Furthermore, the survival rates and total fecundity of O. strigicollis were higher when fed B. tabaci than T. vaporariorum. There were no significant differences in any population parameters of O. strigicollis when fed either whitefly species. These results show that O. strigicollis could survive and maintain its populations on both species of whitefly and could therefore serve as a biological control agent in integrated pest management (IPM).

Global redistribution and local migration in semi-discrete host-parasitoid population dynamic models

Host-parasitoid population dynamics is often probed using a semi-discrete/hybrid modeling framework. Here, the update functions in the discrete-time model connecting year-to-year changes in the population densities are obtained by solving ordinary differential equations that mechanistically describe interactions when hosts become vulnerable to parasitoid attacks. We use this semi-discrete formalism to study two key spatial effects: local movement (migration) of parasitoids between patches during the vulnerable period; and yearly redistribution of populations across patches outside the vulnerable period. Our results show that in the absence of any redistribution, constant density-independent migration and parasitoid attack rates are unable to stabilize an otherwise unstable host-parasitoid population dynamics. Interestingly, inclusion of host redistribution (but not parasitoid redistribution) before the start of the vulnerable period can lead to stable coexistence of both species. Next, we consider a Type-III functional response (parasitoid attack rate increases with host density), where the absence of any spatial effects leads to a neutrally stable host-parasitoid equilibrium. As before, density-independent parasitoid migration by itself is again insufficient to stabilize the population dynamics and host redistribution provides a stabilizing influence. Finally, we show that a Type-III functional response combined with density-dependent parasitoid migration leads to stable coexistence, even in the absence of population redistributions. In summary, we have systematically characterized parameter regimes leading to stable/unstable population dynamics with different forms of spatial heterogeneity coupled to the parasitoid's functional response using mechanistically formulated semi-discrete models.

Interactions between different predator-prey states: a method for the derivation of the functional and numerical response

In this paper we introduce a formal method for the derivation of a predator's functional response from a system of fast state transitions of the prey or predator on a time scale during which the total prey and predator densities remain constant. Such derivation permits an explicit interpretation of the structure and parameters of the functional response in terms of individual behaviour. The same method is also used here to derive the corresponding numerical response of the predator as well as of the prey.

Functional response and size-selective clearance of suspended matter by an invasive mussel

Filter feeding activities link suspension feeders with their environment and underpin their impact on aquatic ecosystems. Despite their ecological and economic impacts, the functional response and size-selective capture of suspended particulates have not been well documented for the golden mussel Limnoperna fortunei. Here we demonstrated that golden mussels had a type I functional response, with an attack rate a = 0.085 and negligible handling time (h). Clearance rate ranged between 72.6 ± 27.0 and 305.5 ± 105.9 mL ind.^-1h^-1 (Mean ± S.E.), depending on food concentrations, which exhibited an inverse relationship with clearance rate. Presence of golden mussels suppressed chlorophyll a concentration in experimental mesocosms, the extent of which was dependent on mussel abundance. Concentration of suspended particles in experimental mesocosms experienced a sharp initial decline across all size categories (≤1->50 μm), though with increased final concentration of large particles (>25 μm), indicating packaging and egestion by golden mussels of fine particles (down to ≤1 μm). Capture efficiency of quantitatively-dominant suspended matter (≤1-50 μm) by golden mussels was inversely related to particle size. Animal abundance, particle size, and their interaction (abundance × particle size) determined the extent to which matter was removed from the water column. Presently L. fortunei occurs primarily in the southern end of the central route of South to North Water Diversion Project (China), but the species is spreading north; we anticipate that impacts associated with filtering of L. fortunei will correspond with local population abundance along this gradient.

Age and sex differences in numerical responses, dietary shifts, and total responses of a generalist predator to population dynamics of main prey

Fluctuations in the abundance of main prey species might shape animal communities, by inducing numerical responses and dietary shifts in predators. Whether numerical responses and dietary shifts differ among individuals of different age and sex has so far gained little attention. These differences could affect how much predators consume main and alternative prey, thus causing variation in predation pressure on main and alternative prey species. We studied the effect of fluctuating main prey abundance (voles) in autumn on the age and sex composition of a food-hoarding population of Eurasian pygmy owls Glaucidium passerinum (327 individuals), and on the species composition of their food stores in western Finland during 2003-2017 (629 food stores). Numbers of yearlings (< 1-year old) of both sexes and adult (+ 1-year old) females increased with increasing vole abundance. During low vole abundance, adult owls stored more small birds and less small mammals than yearlings. Females stored more small mammals than males and showed a tendency to store less birds. The amount of consumed birds (the most important alternative prey), and in particular of crested, willow, great, and blue tits, increased with low vole densities. Our results show that numerical, functional, and total responses of pygmy owls, and probably also other vertebrate predators, to the availability of the main prey in winter are shaped by the age and sex composition of the predator population, which both show large spatio-temporal variation in boreal forests.

Revisiting the "fallacy of averages" in ecology: Expected gain per unit time equals expected gain divided by expected time

Fitness is often defined as the average payoff an animal obtains when it is engaged in several activities, each taking some time. We point out that the average can be calculated with respect to either the time distribution, or to the event distribution of these activities. We show that these two averages lead to the same fitness function. We illustrate this result through two examples from foraging theory, Holling II functional response and the diet choice model, and one game-theoretic example of Hamilton's rule applied to the time-constrained Prisoner's dilemma (PD). In particular, we show that in these models, fitness defined as expected gain per unit time equals fitness defined as expected gain divided by expected time. We also show how these fitnesses predict the optimal outcome for diet choice and the prevalence of cooperation in the repeated PD game.

Are tintinnids picky grazers: Feeding experiments on a mixture of mixotrophic dinoflagellates and implications for red tide dynamics

To understand and predict the outbreak of red tides, which are often dominated by mixotrophic dinoflagellates (MTDs), the effects of "top-down" control by co-occurring predators on red-tide MTDs should be taken into consideration. We studied the numerical and functional responses of the tintinnid ciliate Favella ehrenbergii feeding on two red-tide MTDs, Scrippsiella trochoidea and Heterocapsa triquetra, under single and mixed prey conditions. Our results suggest that a mixed diet could support a better growth of predators compared to a monodiet. In addition, the predators preferred to graze S. trochoidea in the mixed diets, suggesting that predators may switch their feeding preference. The grazing by tintinnid predators could potentially inhibit the outbreaks of red tides dominated by MTDs. The findings in this study provide basic data and new insights for understanding the complex predator-prey relationships in marine microbial food webs, and the dynamics of red tides dominated by MTDs.

FMEM: Functional Mixed Effects Models for Longitudinal Functional Responses

The aim of this paper is to conduct a systematic and theoretical analysis of estimation and inference for a class of functional mixed effects models (FMEM). Such FMEMs consist of fixed effects that characterize the association between longitudinal functional responses and covariates of interest and random effects that capture the spatial-temporal correlations of longitudinal functional responses. We propose local linear estimates of refined fixed effect functions and establish their weak convergence along with a simultaneous confidence band for each fixed-effect function. We propose a global test for the linear hypotheses of varying coefficient functions and derive the associated asymptotic distribution under the null hypothesis and the asymptotic power under the alternative hypothesis are derived. We also establish the convergence rates of the estimated spatial-temporal covariance operators and their associated eigenvalues and eigenfunctions. We conduct extensive simulations and apply our method to a white-matter fiber data set from a national database for autism research to examine the finite-sample performance of the proposed estimation and inference procedures.

Eco-epidemiological interactions with predator interference and infection

Predator interference is a form of competition between predator individuals over access to their prey. There is broad empirical evidence for interference to exist in different strengths in various types of ecological communities. At the same time, parasites are increasingly recognized to alter food web structure and dynamics. In order to investigate the eco-epidemiological interplay between interference and infection, we develop and analyze mathematical models of a predator-prey system, where the predators are subject to both interference and infectious disease. In the absence of infection, equilibrium predator density is known to show a non-monotonic response to interference by first increasing and then decreasing with increasing interference levels. We show that predator infection can change this pattern into a monotonically decreasing predator response to interference, provided the transmissibility is large enough and the pathogenicity is moderate such that the impact of disease on host population density prevails over interference effects. This holds for both types of disease transmission studied here, density-dependent and frequency-dependent. For density-dependent transmission, we find that intermediate values of interference can facilitate disease persistence, whereas the disease would disappear for small or large interference levels. By contrast, for frequency-dependent transmission, disease emergence is independent of interference levels. These dynamic interactions may be important for the understanding of potential biocontrol measures and of spread patterns of zoonotic diseases.

The lacewing Ceraeochrysa caligata as a potential biological agent for controlling the red palm mite Raoiella indica

Background

Compared to chemical control, the use of naturally occurring biological agents to control invasive pests is less threatening to the environment and human health.

Objectives

Here, we assessed the ability of immature stages of the lacewing Ceraeochrysa caligata (Neuroptera: Chrysopidae) to prey upon different developmental stages of the red palm mite Raoiella indica (Acari: Tenuipalpidae), one of the most destructive invasive pests of palm trees in Neotropical regions.

Methods

Increasing densities of three stages of R. indica (eggs, immature stages, and adult females) were offered to C. caligata in coconut leaf arenas. The immature stages of C. caligata were less than 24 h old and were starved before being transferring to the arenas. The amount of prey consumed was recorded 6 h after releasing the C. caligata.

Results

Our results indicated that the ability of C. caligata to feed upon R. indica increased with the larval development of the predator. Higher feeding levels and shorter handling times were recorded for the first and second instars of C. caligata when preying upon the eggs and immature stages of R. indica. Furthermore, C. caligata individuals of different stages exhibited differential functional responses according to prey type (i.e., eggs, immatures, or adult females of R. indica). Ceraeochrysa caligata second instar individuals exhibited a sigmoid increase in consumption rate with increasing prey availability (i.e., a type III functional response) when preying upon immature stages of R. indica. However, when preying upon R. indica adult females, C. caligata second instar individuals exhibited a type II functional response (i.e., an increase in consumption rate with increasing prey availability, before reaching a plateau). Predator individuals of the first and third instar stages exhibited a type II functional response for all prey types.

Conclusions

Collectively, our findings demonstrate that C. caligata, especially at the second instar stage, has potential as a tool for ecological management of the red palm mite.

Foraging behaviour of top predators mediated by pollution of psychoactive pharmaceuticals and effects on ecosystem stability

Although pharmaceuticals are recognized as a major threat to aquatic ecosystems worldwide, little is known about their ecological effect on aquatic biota and ecosystems. Drug-induced behaviour changes could have a substantial impact on consumer-resource interactions influencing stability of the community and ecosystem. We combined laboratory experiments and functional response modelling to investigate effects of real wastewater treatment plant (WWTP) effluent, as well as environmentally relevant concentrations of the antidepressants citalopram and opioid pain medication tramadol, on trophic interactions. Our biological system consisted of dragonfly Aeshna cyanea larvae as predator of common carp Cyprinus carpio fry. Exposure to WWTP effluent significantly increased A. cyanea maximum feeding rate, while those parameters in tramadol and citalopram-exposed larvae were significantly lower from unexposed control group. This suggested the potential of all tested pollutants to have an effect on consumer-resource equilibrium in aquatic ecosystems. While WWTP effluent strengthened interaction strength (IS) of consumer-resource interaction dynamics making the food web more vulnerable to fluctuation and destabilization, tramadol and citalopram could inhibit the potential oscillations of the consumer-resource system by weakening the IS. Similar studies to reveal the potential of pervasive pharmaceuticals to change of consumer-resource interactions dynamics are needed, especially when real WWTP effluent consisting of mixture of various pharmaceuticals displayed very different effect from single compounds tested.

Functional consumers regulate the effect of availability of subsidy on trophic cascades in the Yellow River Delta, China

Understanding the environmental context where heterogeneous ecological processes affect biotic interactions is a key aim of ecological research. However, mechanisms underlying spatial variation in trophic interactions linked to resource availability across ecosystem gradients remains unclear. We experimentally manipulated the interactive effects of predator fish and quantitative gradient of leaf detritus on macroinvertebrates and benthic algae. We found that non-linear changes in the strength of trophic cascades were strongly linked to the retention rates of experimental leaf detritus and also determined by predatory consumers. Retention rate of leaf detritus influenced the recruitment of predatory invertebrates and foraging preference of predators, accounting for largely the variations in shift of strengthening and weakening trophic cascades. Our results highlight the importance to identify joint processes of recruitment and foraging responses of functional consumer in understanding the impacts of both anthropogenic and natural alterations in subsidy on trophic interaction of coastal food webs.

The influence of microplastics on trophic interaction strengths and oviposition preferences of dipterans

Microplastic (MP) pollution continues to proliferate in freshwater, marine and terrestrial environments, but with their biotic implications remaining poorly understood. Biotic interactions such as predation can profoundly influence ecosystem structuring, stability and functioning. However, we currently lack quantitative understandings of how trophic interaction strengths and associated behaviours are influenced by MP pollution, and how transference of MPs between trophic levels relates to consumptive traits. We also lack understanding of key life-history effects of MPs, for example, reproductive strategies such as oviposition. The present study examines the predatory ability of non-biting midge larvae, Chaoborus flavicans, towards larvae of Culex pipiens mosquitoes when the latter are exposed to MPs, using a functional response (FR) approach. Transfer of MPs occurred from larval mosquitoes to larval midges via predation. Microplastics transfer was significantly positively related to predation rates. Predation by C. flavicans followed a Type II FR, with average maximum feeding rates of 6.2 mosquito larvae per hour. These and other FR parameters (attack rates and handling times) were not significantly influenced by the presence of MPs. Further, C. pipiens adults did not avoid ovipositing in habitats with high concentrations of MPs. We thus demonstrate that MPs can move readily through freshwater food webs via biotic processes such as predation, and that uptake correlates strongly with consumption rates. Further, as MPs do not deter adult mosquitoes from ovipositing, our experiments reveal high potential for MP exposure and transference through ecosystems.

Sublethal and transgenerational effects of thiamethoxam on the demographic fitness and predation performance of the seven-spot ladybeetle Coccinella septempunctata L. (Coleoptera: Coccinellidae)

Seven-spot ladybird beetles, Coccinella septempunctata L., are critical aphidophagous predators in the agricultural environment. Thiamethoxam, a neonicotinoid insecticide, is commonly used for controlling pests but impairs their natural enemies at the same time. To improve effective IPM (integrated pest management) strategies, we evaluated the sublethal and transgenerational effects of thiamethoxam on C. septempunctata. Our results showed that thiamethoxam at doses of 0.1 × LC₁₀ (0.053 mg L^-1) and LC₁₀ (0.53 mg L^-1) significantly reduced adult emergence, fecundity and fertility of the parental generation. In unexposed progeny (F₁) of thiamethoxam-exposed parents, at the two doses 0.1 × LC₁₀ and LC₁₀, the larval stage was prolonged, and total longevity was decreased by 18.76 and 24.46%, respectively. The higher concentrations (0.1 × LC₁₀ and LC₁₀) also decreased the fecundity by 33.74 and 46.56%, respectively, and the oviposition period by 19.67 and 25.01%, respectively. In addition, demographic parameters including the intrinsic rate of increase (r), finite rate of increase (λ), net reproductive rate (R₀), and mean generation time (T) were significantly reduced when exposed to LC₁₀. Moreover, the predation activity of the F₁ generation was reduced by the transgenerational effects of LC₁₀. These results disclose negative influence of thiamethoxam at sublethal concentrations on this ladybird predator and its subsequent generation.

Hourly movement decisions indicate how a large carnivore inhabits developed landscapes

The ecology of wildlife living in proximity to humans often differs from that in more natural places. Animals may perceive anthropogenic features and people as threats, exhibiting avoidance behavior, or may acclimate to human activities. As development expands globally, changes in the ecology of species in response to human phenomena may determine whether animals persist in these changing environments. We hypothesize that American black bears (Ursus americanus) persist within developed areas by effectively avoiding risky landscape features. We test this by quantifying changes in the movements of adult females from a population living within exurban and suburban development. We collected hourly GPS data from 23 individuals from 2012 to 2014 and used step-selection functions to estimate selection for anthropogenic features. Females were more avoidant of roads and highways when with cubs than without and were more responsive to increased traffic volume. As bears occupied greater housing densities, selection for housing increased, while avoidance of roads and responsiveness to traffic increased. Behavioral flexibility allowed bears in highly developed areas to alter selection and avoidance for anthropogenic features seasonally. These findings support the hypothesis that black bears perceive human activity as risky, and effectively avoid these risks while inhabiting developed areas. We document a high amount of individual variation in selection of anthropogenic features within the study population. Our findings suggest that initially, wildlife can successfully inhabit developed landscapes by effectively avoiding human activity. However, variation among individuals provides the capacity for population-level shifts in behavior over time.

All you can eat: the functional response of the cold-water coral Desmophyllum dianthus feeding on krill and copepods

The feeding behavior of the cosmopolitan cold-water coral (CWC) Desmophyllum dianthus (Cnidaria: Scleractinia) is still poorly known. Its usual deep distribution restricts direct observations, and manipulative experiments are so far limited to prey that do not occur in CWC natural habitat. During a series of replicated incubations, we assessed the functional response of this coral feeding on a medium-sized copepod (Calanoides patagoniensis) and a large euphausiid (Euphausia vallentini). Corals showed a Type I functional response, where feeding rate increased linearly with prey abundance, as predicted for a tentaculate passive suspension feeder. No significant differences in feeding were found between prey items, and corals were able to attain a maximum feeding rate of 10.99 mg C h⁻¹, which represents an ingestion of the 11.4% of the coral carbon biomass per hour. These findings suggest that D. dianthus is a generalist zooplankton predator capable of exploiting dense aggregations of zooplankton over a wide prey size-range.

Size-dependent functional response of Xenopus laevis feeding on mosquito larvae

Predators can play an important role in regulating prey abundance and diversity, determining food web structure and function, and contributing to important ecosystem services, including the regulation of agricultural pests and disease vectors. Thus, the ability to predict predator impact on prey is an important goal in ecology. Often, predators of the same species are assumed to be functionally equivalent, despite considerable individual variation in predator traits known to be important for shaping predator–prey interactions, like body size. This assumption may greatly oversimplify our understanding of within-species functional diversity and undermine our ability to predict predator effects on prey. Here, we examine the degree to which predator–prey interactions are functionally homogenous across a natural range of predator body sizes. Specifically, we quantify the size-dependence of the functional response of African clawed frogs (Xenopus laevis) preying on mosquito larvae (Culex pipiens). Three size classes of predators, small (15–30 mm snout-vent length), medium (50–60 mm) and large (105–120 mm), were presented with five densities of prey to determine functional response type and to estimate search efficiency and handling time parameters generated from the models. The results of mesocosm experiments showed that type of functional response of X. laevis changed with size: small predators exhibited a Type II response, while medium and large predators exhibited Type III responses. Functional response data showed an inversely proportional relationship between predator attack rate and predator size. Small and medium predators had highest and lowest handling time, respectively. The change in functional response with the size of predator suggests that predators with overlapping cohorts may have a dynamic impact on prey populations. Therefore, predicting the functional response of a single size-matched predator in an experiment may misrepresent the predator’s potential impact on a prey population.

Comparative feeding rates of native and invasive ascidians

Ascidians have a recent history of species introductions globally, often with strong ecological impacts. Comparisons of per capita effects of invaders and comparable natives are useful to assess such impacts. Here, we explore ingestion rates (IR) and clearance rates (CR) of Ciona intestinalis and Ciona robusta, co-occurring native and non-native ascidians, respectively, from Brittany, France. IR was positively related to food concentration, with the invader responding more strongly to increasing food concentration. CR also differed by species, with the invader demonstrating higher values. C. robusta exhibited a higher functional response (Type I) than did C. intestinalis (Type II). Relative impact measured using seasonal abundance and IR revealed that C. robusta has a much greater impact than C. intestinalis at all food concentrations tested, though the former has a constrained distribution which limits its regional impact. Nevertheless, when abundant, we expect C. robusta to exert a greater impact on algal foods.

Shifting interaction strength between estuarine mysid species across a temperature gradient

In many coastal regions, mean coastal atmospheric and water temperatures are projected to shift as climate change ensues. Interaction strengths between organisms are likely to change along with environmental changes, given interspecific heterogeneity in responses to physico-chemical variables. Biological interaction outcomes have the potential to alter food web production and trophic level biomass distribution. This is particularly pertinent for key species that are either abundant or play disproportionately large roles in ecosystem processes. Using a functional response approach, we quantified the effects of shifting temperatures on interactions between key mysid species-sympatric in their distribution across a biogeographic transition zone along the east coast of South Africa. The Rhopalophthalmus terranatalis functional response type toward Mesopodopsis wooldridgei prey was independent of temperature, with all treatments producing Type II functional responses. Temperature effects on predator-prey dynamics were, however, evident as interaction strength was greatest at 21 °C, as measured by maximum feeding rates. Unlike maximum feeding rate, attack rates increased linearly with increasing temperature across the experimental treatments. Our findings suggest that interaction strength between the mysid shrimp species is likely to vary spatially along the current length of their sympatric distribution and temporally in certain regions where temperatures are projected to change. Such experimental interaction investigations are becoming increasingly important given our relatively poor understanding of the consequences of environmental change for effects on interactions among species and their wider ecosystem implications.

The Rosenzweig-MacArthur system via reduction of an individual based model

The Rosenzweig-MacArthur system is a particular case of the Gause model, which is widely used to describe predator-prey systems. In the classical derivation, the interaction terms in the differential equation are essentially derived from considering handling time vs. search time, and moreover there exist derivations in the literature which are based on quasi-steady state assumptions. In the present paper we introduce a derivation of this model from first principles and singular perturbation reductions. We first establish a simple stochastic mass action model which leads to a three-dimensional ordinary differential equation, and systematically determine all possible singular perturbation reductions (in the sense of Tikhonov and Fenichel) to two-dimensional systems. Among the reductions obtained we find the Rosenzweig-MacArthur system for a certain choice of small parameters as well as an alternative to the Rosenzweig-MacArthur model, with density dependent death rates for predators. The arguments to obtain the reductions are intrinsically mathematical; no heuristics are employed.

A study of the reproducibility and reliability of the musculo-articular stiffness of the ankle joint

The objective of this work was to evaluate the reproducibility, reliability and usefulness of the musculo-articular stiffness (MAS) of the ankle joint, measuring it by the free vibration technique. Seventeen (nine males and eight females) healthy university students were included in the study. Force (f), MAS (k) and unitary MAS (k_u) (defined as the ratio between the value of stiffness k obtained in the test (absolute terms) and the value of force (f)) were obtained. A test-retest protocol was designed and performed on the same day to determine the short-term reproducibility of f, k and k_u. Short-term reproducibility of k and k_u on 1 day in absolute terms (< 7% Coefficient of Variation (CV)) and relative reproducibility (Intraclass Correlation Coefficient (ICC) and Pearson ≥ 0.97) for both feet were obtained. The reliability of k and k_u in absol_ute terms (< 9% CV) and in relative terms (ICC and Pearson ≥0.93) based on repeating the protocol for 1 week was analysed for both feet. To analyse the usefulness, the Effect Size (ES) ratio = "Trivial" for all variables (for 1 day and 1 week) and the Smallest Worthwhile Change (SWC) ratio (Typical Error (TE)< SWC) = "GOOD" for k and k_u (1 day and 1 week) were considered. The Minimal Difference needed to be considered "real" (MD) for k_u ≅ 3.5% (1 day); k_u≅ 8.5% (1 week) (p < 0.05) was obtained. The statistical analysis carried out displayed the high reproducibility, reliability and usefulness of the MAS test, which was more consistent with k_u than k. Therefore, the unitary stiffness (k_u) proven to be representative of the mechanical response of the ankle joint obtained by free vibration techniques, which allows comparison between different subjects.

Immune system handling time may alter the outcome of competition between pathogens and the immune system

Predators may be limited in their ability to kill prey (i.e., have type II or III functional responses), an insight that has had far-reaching consequences in the ecological literature. With few exceptions, however, this possibility has not been extended to the behaviour of immune cells, which kill pathogens much as predators kill their prey. Rather, models of the within-host environment have tended to tacitly assume that immune cells have an unlimited ability to target and kill pathogens (i.e., a type I functional response). Here we explore the effects of changing this assumption on infection outcomes (i.e., pathogen loads). We incorporate immune cell handling time into an ecological model of the within-host environment that considers both the predatory nature of the pathogen-immune cell interaction as well as competition between immune cells and pathogens for host resources. Unless pathogens can preempt immune cells for host resources, adding an immune cell handling time increases equilibrium pathogen load. We find that the shape of the relationship between energy intake and pathogen load can change: with a type I functional response, pathogen load is maximised at intermediate inputs, while for a type II or III functional response, pathogen load is solely increasing. With a type II functional response, pathogen load can fluctuate rather than settling to an equilibrium, a phenomenon unobserved with type I or III functional responses. Our work adds to a growing literature highlighting the role of resource availability in host-parasite interactions. Implications of our results for adaptive anorexia are discussed.

Is structural sensitivity a problem of oversimplified biological models? Insights from nested Dynamic Energy Budget models

Many current issues in ecology require predictions made by mathematical models, which are built on somewhat arbitrary choices. Their consequences are quantified by sensitivity analysis to quantify how changes in model parameters propagate into an uncertainty in model predictions. An extension called structural sensitivity analysis deals with changes in the mathematical description of complex processes like predation. Such processes are described at the population scale by a specific mathematical function taken among similar ones, a choice that can strongly drive model predictions. However, it has only been studied in simple theoretical models. Here, we ask whether structural sensitivity is a problem of oversimplified models. We found in predator-prey models describing chemostat experiments that these models are less structurally sensitive to the choice of a specific functional response if they include mass balance resource dynamics and individual maintenance. Neglecting these processes in an ecological model (for instance by using the well-known logistic growth equation) is not only an inappropriate description of the ecological system, but also a source of more uncertain predictions.

Sensitivity of the dynamics of the general Rosenzweig-MacArthur model to the mathematical form of the functional response: a bifurcation theory approach

The equations in the Rosenzweig-MacArthur predator-prey model have been shown to be sensitive to the mathematical form used to model the predator response function even if the forms used have the same basic shape: zero at zero, monotone increasing, concave down, and saturating. Here, we revisit this model to help explain this sensitivity in the case of three response functions of Holling type II form: Monod, Ivlev, and Hyperbolic tangent. We consider both the local and global dynamics and determine the possible bifurcations with respect to variation of the carrying capacity of the prey, a measure of the enrichment of the environment. We give an analytic expression that determines the criticality of the Hopf bifurcation, and prove that although all three forms can give rise to supercritical Hopf bifurcations, only the Trigonometric form can also give rise to subcritical Hopf bifurcation and has a saddle node bifurcation of periodic orbits giving rise to two coexisting limit cycles, providing a counterexample to a conjecture of Kooji and Zegeling. We also revisit the ranking of the functional responses, according to their potential to destabilize the dynamics of the model and show that given data, not only the choice of the functional form, but the choice of the number and/or position of the data points can influence the dynamics predicted.

Transient effects following peak exposures towards pesticides - An explanation for the unresponsiveness of in situ measured functional variables

Invertebrate-mediated leaf litter decomposition is frequently used to assess stress-related implications in stream ecosystem integrity. In situ measures such as the mass loss from leaf bags or the feeding of caged invertebrates deployed for days or weeks may, however, fail to detect transient effects due to recovery or compensatory mechanisms. We assessed the relevance of transient effects using the peak exposure towards an insecticide (i.e., etofenprox) as a model scenario at three levels of complexity. These were 1) the assessment of the decomposition realised by invertebrate communities in stream mesocosms over 21 days via leaf bags, 2) 7-days lasting in situ bioassays quantifying the leaf consumption of Gammarus fossarum, and 3) a laboratory experiment determining the daily feeding rate of the same species over 7 days. Etofenprox did not trigger a significantly altered decomposition by invertebrate communities during the leaf bag assay, while in situ bioassays detected a significant reduction in gammarids' feeding rate at the highest tested concentration. The laboratory bioassay suggests that observed mismatches might be explained by recovery and post-exposure compensation. As leaf-shredding invertebrates are likely in a vulnerable state following transient effects, biomonitoring for implications of peak exposures and other pulsed stress events must happen at an adequate temporal resolution.

Range expansion of a fouling species indirectly impacts local species interactions

We investigated how recent changes in the distribution and abundance of a fouling organism affected the strength of interactions between a commercially important foundation species and a common predator. Increases in the abundance of boring sponges that bioerode the calcified shells of oysters and other shelled organisms have been attributed to increased salinization of estuarine ecosystems. We tested the hypothesis that fouling by boring sponges will change the interaction strength between oysters and a common predator (stone crabs). We generated five oyster density treatments crossed with two sponge treatments (sponge and no sponge). We contrasted the interaction strength between stone crabs and fouled and non-fouled oysters by comparing the parameters of fitted functional response curves based on Rogers random predation model. We found that fouled oysters suffered higher predation from stone crabs, and that the increased predation risk stemmed from a reduction in the handling time needed to consume the fouled oysters. These findings highlight the importance of understanding the effects of abiotic changes on both the composition of ecological communities, and on the strengths of direct and indirect interactions among species. Global climate change is altering local ecosystems in complex ways, and the success of restoration, management, and mitigation strategies for important species requires a better appreciation for how these effects cascade through ecosystems.

How to be a good neighbour: Facilitation and competition between two co-flowering species

Empirical evidence suggests that co-flowering species can facilitate each other through shared pollinators. However, the extent to which one co-flowering species can relieve pollination limitation of another while simultaneously competing for abiotic resource has rarely been examined. Using a deterministic model we explored the demographic outcome for one ("focal") species of its co-occurrence with a species that shares pollinators and competes for both pollinator visitation and abiotic resources. In this paper we showed how the overall impact can be positive or negative, depending on the balance between enhanced fertilization versus increased competition. Our model could predict the density of co-flowering species that will maximize the pollination rate of the focal species by attracting pollinators. Because that density will also give rise to competitive effects, a lower density of co-flowering species is required for optimizing the trade-off between enhanced fertilization and competition so as to give the maximum possible facilitation of reproduction in the focal species. Results were qualitatively different when we considered attractiveness of the co-flowering species, as opposed to its density, because attractiveness, unlike density, had no effect on competition for abiotic resources. Whereas unattractive neighbours would not bring in pollinators, very attractive neighbours would captivate pollinators, not sharing them with the focal species. Thus optimal benefit to the focal species came at intermediate levels of attractiveness in the co-flowering species. This intermediate level of attractiveness in co-flowering species simultaneously maximized pollination and overall facilitation of reproduction for the focal species. The likelihood of facilitation was predicted to decline with the selfing rate of the focal species, revealing an indirect cost for an inbreeding mating system. Whether a co-flowering species can be facilitative depends on the way pollinators respond to the plant density: only a Type III functional response for visitation rate can result in facilitation. Our model provided both a conceptual framework and precise quantitative measures for determining the impacts of a neighbouring co-flowering species on reproduction.

How patch size and refuge availability change interaction strength and population dynamics: a combined individual- and population-based modeling experiment

Knowledge on how functional responses (a measurement of feeding interaction strength) are affected by patch size and habitat complexity (represented by refuge availability) is crucial for understanding food-web stability and subsequently biodiversity. Due to their laborious character, it is almost impossible to carry out systematic empirical experiments on functional responses across wide gradients of patch sizes and refuge availabilities. Here we overcame this issue by using an individual-based model (IBM) to simulate feeding experiments. The model is based on empirically measured traits such as body-mass dependent speed and capture success. We simulated these experiments in patches ranging from sizes of petri dishes to natural patches in the field. Moreover, we varied the refuge availability within the patch independently of patch size, allowing for independent analyses of both variables. The maximum feeding rate (the maximum number of prey a predator can consume in a given time frame) is independent of patch size and refuge availability, as it is the physiological upper limit of feeding rates. Moreover, the results of these simulations revealed that a type III functional response, which is known to have a stabilizing effect on population dynamics, fitted the data best. The half saturation density (the prey density where a predator consumes half of its maximum feeding rate) increased with refuge availability but was only marginally influenced by patch size. Subsequently, we investigated how patch size and refuge availability influenced stability and coexistence of predator-prey systems. Following common practice, we used an allometric scaled Rosenzweig–MacArthur predator-prey model based on results from our in silico IBM experiments. The results suggested that densities of both populations are nearly constant across the range of patch sizes simulated, resulting from the constant interaction strength across the patch sizes. However, constant densities with decreasing patch sizes mean a decrease of absolute number of individuals, consequently leading to extinction of predators in the smallest patches. Moreover, increasing refuge availabilities also allowed predator and prey to coexist by decreased interaction strengths. Our results underline the need for protecting large patches with high habitat complexity to sustain biodiversity.

Effect of temperature and search area on the functional response of Anisops sardea (Hemiptera: Notonectidae) against Anopheles stephensi in laboratory bioassay

Present study was carried out to establish the influence of abiotic factors on foraging activities of a predatory hemipteran insect Anisops sardea against Anopheles stephensi larvae. The functional response of A. sardea was evaluated in variable density of prey items with variation in search area (100, 250, 500 and 1000ml water volume) and temperatures (20, 25 and 30°C). The results of laboratory bioassay revealed that prey consumption rate of predator species was positively related with increasing temperature and inversely related with increasing search area. Polynomial logistic regression equations and associated parameters showed that A. sardea exhibited a type II functional response in variable search area and type-III response at variable temperatures. Related response specific attack rates and handling times were also evaluated in presence of specific abiotic factors.

Diverse foraging opportunities drive the functional response of local and landscape-scale bear predation on Pacific salmon

The relationship between prey abundance and predation is often examined in single habitat units or populations, but predators may occupy landscapes with diverse habitats and foraging opportunities. The vulnerability of prey within populations may depend on habitat features that hinder predation, and increased density of conspecifics in both the immediate vicinity and the broader landscape. We evaluated the relative effects of physical habitat, local, and neighborhood prey density on predation by brown bears on sockeye salmon in a suite of 27 streams using hierarchical Bayesian functional response models. Stream depth and width were inversely related to the maximum proportion of salmon killed, but not the asymptotic limit on total number killed. Interannual variation in predation was density dependent; the number of salmon killed increased with fish density in each stream towards an asymptote. Seven streams in two geographical groups with ≥23 years of data in common were then analyzed for neighborhood density effects. In most (12 of 18) cases predation in a stream was reduced by increasing salmon abundance in neighboring streams. The uncertainty in the estimates for these neighborhood effects may have resulted from interactions between salmon abundance and habitat that influenced foraging by bears, and from bear behavior (e.g., competitive exclusion) and abundance. Taken together, the results indicated that predator-prey interactions depend on density at multiple spatial scales, and on habitat features of the surrounding landscape. Explicit consideration of this context dependency should lead to improved understanding of the ecological impacts of predation across ecosystems and taxa.

You can hide but you can't run: apparent competition, predator responses and the decline of Arctic ground squirrels in boreal forests of the southwest Yukon

Throughout much of North America's boreal forest, the cyclical fluctuations of snowshoe hare populations (Lepus americanus) may cause other herbivores to become entrained in similar cycles. Alternating apparent competition via prey switching followed by positive indirect effects are the mechanisms behind this interaction. Our purpose is to document a change in the role of indirect interactions between sympatric populations of hares and arctic ground squirrels (Urocitellus parryii plesius), and to emphasize the influence of predation for controlling ground squirrel numbers. We used mark-recapture to estimate the population densities of both species over a 25-year period that covered two snowshoe hare cycles. We analysed the strength of association between snowshoe hare and ground squirrel numbers, and the changes to the seasonal and annual population growth rates of ground squirrels over time. A hyperbolic curve best describes the per capita rate of increase of ground squirrels relative to their population size, with a single stable equilibrium and a lower critical threshold below which populations drift to extinction. The crossing of this unstable boundary resulted in the subsequent uncoupling of ground squirrel and hare populations following the decline phase of their cycles in 1998. The implications are that this sustained Type II predator response led to the local extinction of ground squirrels. When few individuals are left in a colony, arctic ground squirrels may also have exhibited an Allee effect caused by the disruption of social signalling of approaching predators.