Effects of acute and chronic temperature changes on the functional responses of the dogfish <Emphasis Type="Italic">Scyliorhinus canicula</Emphasis> (Linnaeus, 1758) towards amphipod prey <Emphasis Type="Italic">Echinogammarus marinus</Emphasis> (Leach, 1815)

Predation is a strong driver of population dynamics and community structure and it is essential to reliably quantify and predict predation impacts on prey populations in a changing thermal landscape. Here, we used comparative functional response analyses to assess how predator-prey interactions between dogfish and invertebrate prey change under different warming scenarios. The Functional Response Type, attack rate, handling time and maximum feeding rate estimates were calculated for Scyliorhinus canicula preying upon Echinogammarus marinus under temperatures of 11.3 °C and 16.3 °C, which represent both the potential daily variation and predicted higher summer temperatures within Strangford Lough, N. Ireland. A two x two design of “Predator Acclimated”, “Prey Acclimated”, “Both Acclimated”, and “Both Unacclimated” was implemented to test functional responses to temperature rise. Attack rate was higher at 11.3 °C than at 16.3 °C, but handling time was lower and maximum feeding rates were higher at 16.3 °C. Non-acclimated predators had similar maximum feeding rate towards non-acclimated and acclimated prey, whereas acclimated predators had significantly higher maximum feeding rates towards acclimated prey as compared to non-acclimated prey. Results suggests that the predator attack rate is decreased by increasing temperature but when both predator and prey are acclimated the shorter handling times considerably increase predator impact. The functional response of the fish changed from Type II to Type III with an increase in temperature, except when only the prey were acclimated. This change from population destabilizing Type II to more stabilizing Type III could confer protection to prey at low densities but increase the maximum feeding rate by Scyliorhinus canicula in the future. However, predator movement between different thermal regimes may maintain a Type II response, albeit with a lower maximum feeding rate. This has implications for the way the increasing population Scyliorhinus canicula in the Irish Sea may exploit valuable fisheries stocks in the future.     

Effects of prey-size and predator-instar on the predation of Daphnia by Notonecta

Abstract. 1. Attack rates and handling times are measured by a series of separate functional response experiments for each instar of Notonecta glauca attacking four size classes of Daphnia magna as prey. The resulting attack rate and handling time surfaces are complex, with maximum attack rates for small predators attacking small prey, and large predators attacking large prey. Adult Notonecta have lower attack rates than the two previous juvenile instars (4 and 5).
2. The literature on attack rates and handling times in other predator—prey interactions that involve a series of different predator and prey size or age classes is reviewed in the context of the Notonecta-Daphnia results. The data suggest that small predator instars will usually compete with large instars for food, unless there is spatial or temporal separation between them.
3. Complex attack rate and handling time surfaces are to be expected wherever a wide range of prey and predator sizes is involved.
4. Size related changes in attack rates and handling times can introduce very complex dynamics into predator-prey interactions.     

Not attackable or not crackable—How pre‐ and post‐attack defenses with different competition costs affect prey coexistence and population dynamics

It is well‐known that prey species often face trade‐offs between defense against predation and competitiveness, enabling predator‐mediated coexistence. However, we lack an understanding of how the large variety of different defense traits with different competition costs affects coexistence and population dynamics. Our study focusses on two general defense mechanisms, that is, pre‐attack (e.g., camouflage) and post‐attack defenses (e.g., weaponry) that act at different phases of the predator—prey interaction. We consider a food web model with one predator, two prey types and one resource. One prey type is undefended, while the other one is pre‐ or post‐attack defended paying costs either by a higher half‐saturation constant for resource uptake or a lower maximum growth rate. We show that post‐attack defenses promote prey coexistence and stabilize the population dynamics more strongly than pre‐attack defenses by interfering with the predator's functional response: Because the predator spends time handling “noncrackable” prey, the undefended prey is indirectly facilitated. A high half‐saturation constant as defense costs promotes coexistence more and stabilizes the dynamics less than a low maximum growth rate. The former imposes high costs at low resource concentrations but allows for temporally high growth rates at predator‐induced resource peaks preventing the extinction of the defended prey. We evaluate the effects of the different defense mechanisms and costs on coexistence under different enrichment levels in order to vary the importance of bottom‐up and top‐down control of the prey community.     

Microgeographic divergence of functional responses among salamanders under antagonistic selection from apex predators

A predator''s functional response determines predator–prey interactions by describing the relationship between the number of prey available and the number eaten. Its shape and parameters fundamentally govern the dynamic equilibrium of predator–prey interactions and their joint abundances. Yet, estimates of these key parameters generally assume stasis in space and time and ignore the potential for local adaptation to alter feeding responses and the stability of trophic dynamics. Here, we evaluate if functional responses diverge among populations of spotted salamander (Ambystoma maculatum) larvae that face antagonistic selection on feeding strategies based on their own risk of predation. Common garden experiments revealed that spotted salamander from ponds with varying predation risks differed in their functional responses, suggesting an evolutionary response. Applying mechanistic equations, we discovered that the combined changes in attack rates, handling times and shape of the functional response enhanced feeding rate in environments with high densities of gape-limited predators. We suggest how these parameter changes could alter community equilibria and other emergent properties of food webs. Community ecologists might often need to consider how local evolution at fine scales alters key relationships in ways that alter local diversity patterns, food web dynamics, resource gradients and community responses to disturbance.     

Handling time promotes the coevolution of aggregation in predator-prey systems

Predators often have type II functional responses and live in environments where their life history traits as well as those of their prey vary from patch to patch. To understand how spatial heterogeneity and predator handling times influence the coevolution of patch preferences and ecological stability, we perform an ecological and evolutionary analysis of a Nicholson-Bailey type model. We prove that coevolutionarily stable prey and searching predators prefer patches that in isolation support higher prey and searching predator densities, respectively. Using this fact, we determine how environmental variation and predator handling times influence the spatial patterns of patch preferences, population abundances and per-capita predation rates. In particular, long predator handling times are shown to result in the coevolution of predator and prey aggregation. An analytic expression characterizing ecological stability of the coevolved populations is derived. This expression implies that contrary to traditional theoretical expectations, predator handling time can stabilize predator-prey interactions through its coevolutionary influence on patch preferences. These results are shown to have important implications for classical biological control.     

PREY ADAPTATION AS A CAUSE OF PREDATOR-PREY CYCLES

We analyze simple models of predator-prey systems in which there is adaptive change in a trait of the prey that determines the rate at which it is captured by searching predators. Two models of adaptive change are explored: (1) change within a single reproducing prey population that has genetic variation for vulnerability to capture by the predator; and (2) direct competition between two independently reproducing prey populations that differ in their vulnerability. When an individual predator's consumption increases at a decreasing rate with prey availability, prey adaptation via either of these mechanisms may produce sustained cycles in both species' population densities and in the prey's mean trait value. Sufficiently rapid adaptive change (e.g., behavioral adaptation or evolution of traits with a large additive genetic variance), or sufficiently low predator birth and death rates will produce sustained cycles or chaos, even when the predator-prey dynamics with fixed prey capture rates would have been stable. Adaptive dynamics can also stabilize a system that would exhibit limit cycles if traits were fixed at their equilibrium values. When evolution fails to stabilize inherently unstable population interactions, selection decreases the prey's escape ability, which further destabilizes population dynamics. When the predator has a linear functional response, evolution of prey vulnerability always promotes stability. The relevance of these results to observed predator-prey cycles is discussed.     

Consequences of size structure in the prey for predator-prey dynamics: the composite functional response

1. Current formulations of functional responses assume that the prey is homogeneous and independent of intraspecific processes. Most prey populations consist of different coexisting size classes that often engage in asymmetrical intraspecific interactions, including cannibalism, which can lead to nonlinear interaction effects. This may be important as the size structure with the prey could alter the overall density-dependent predation rates. 2. In a field experiment with damselfly and dragonfly larvae, 16 treatments manipulated the density of a small prey stage, the presence of large conspecific prey and the presence of heterospecific predators. 3. Size structure in the prey (i.e. when both prey stages were present) decreased the impact of the predator on overall prey mortality by 25-48% at mid and high prey densities, possibly due to density-dependent size-structured cannibalism in the prey. The predation rates on small prey stages were determined by the interaction of large prey and predators. Predation rates increased with prey density in the absence of large prey, but predation rates were constant across densities when large conspecifics were present. 4. The functional response for unstructured prey followed a Holling type III model, but the predation rate for size-structured prey was completely different and followed a complex pattern that could not be explained with any standard functional response. 5. Using additional laboratory experiments, a mortality model was developed and parameterized. It showed that the overall prey mortality of size-structured prey can be adequately predicted with a composite functional response model that modelled the individual functional responses of each prey stage separately and accounted for their cannibalistic interaction. 6. Thus, treating a prey population as a homogeneous entity will lead to erroneous predictions in most real-world food webs. However, if we account for the effects of size structure and the intraspecific interactions on functional responses by treating size classes as different functional groups, it is possible to reliably predict the dynamics of size-structured predator-prey systems.     

Habitat choice in predator-prey systems: spatial instability due to interacting adaptive movements

The role of habitat choice behavior in the dynamics of predator-prey systems is explored using simple mathematical models. The models assume a three-species food chain in which each population is distributed across two or more habitats. The predator and prey adjust their locations dynamically to maximize individual per capita growth, while the prey's resource has a low rate of random movement. The two consumer species have Type II functional responses. For many parameter sets, the populations cycle, with predator and prey "chasing" each other back and forth between habitats. The cycles are driven by the aggregation of prey, which is advantageous because the predator's saturating functional response induces a short-term positive density dependence in prey fitness. The advantage of aggregation in a patch is only temporary because resources are depleted and predators move to or reproduce faster in the habitat with the largest number of prey, perpetuating the cycle. Such spatial cycling can stabilize population densities and qualitatively change the responses of population densities to environmental perturbations. These models show that the coupled processes of moving to habitats with higher fitness in predator and prey may often fail to produce ideal free distributions across habitats.     

Experimental duration and predator satiation levels systematically affect functional response parameters

Empirical feeding studies where density‐dependent consumption rates are fitted to functional response models are often used to parameterize the interaction strengths in models of population or food‐web dynamics. However, the relationship between functional response parameter estimates from short‐term feeding studies and real‐world, long‐term, trophic interaction strengths remains largely unexamined. In a critical first step to address this void, we tested for systematic effects of experimental duration and predator satiation on the estimate of functional response parameters, namely attack rate and handling time. Analyzing a large data set covering a wide range of predator taxa and body masses, we show that attack rates decrease with increasing experimental duration, and that handling times of starved predators are consistently shorter than those of satiated predators. Therefore, both the experimental duration and the predator satiation level have a strong and systematic impact on the predictions of population dynamics and food‐web stability. Our study highlights potential pitfalls at the intersection of empirical and theoretical applications of functional responses. We conclude our study with some practical suggestions for how these implications should be addressed in the future to improve predictive abilities and realism in models of predator–prey interactions.     

Modeling changes in predator functional response to prey across spatial scales

Extrapolation of predator functional responses from laboratory observations to the field is often necessary to predict predation rates and predator-prey dynamics at spatial and temporal scales that are difficult to observe directly. We use a spatially explicit individual-based model to explore mechanisms behind changes in functional responses when the scale of observation is increased. Model parameters were estimated from a predator-prey system consisting of the predator Delphastus catalinae (Coleoptera: Coccinellidae) and Bemisia tabaci biotype B (Hemiptera: Aleyrodidae) on tomato plants. The model explicitly incorporates prey and predator distributions within single plants, the search behavior of predators within plants, and the functional response to prey at the smallest scale of interaction (within leaflets) observed in the laboratory. Validation revealed that the model is useful in scaling up from laboratory observations to predation in whole tomato plants of varying sizes. Comparing predicted predation at the leaflet scale, as observed in laboratory experiments, with predicted predation on whole plants revealed that the predator functional response switches from type II within leaflets to type III within whole plants. We found that the magnitude of predation rates and the type of functional response at the whole plant scale are modulated by (1) the degree of alignment between predator and prey distributions and (2) predator foraging behavior, particularly the effect of area-concentrated search within plants when prey population density is relatively low. The experimental and modeling techniques we present could be applied to other systems in which active predators prey upon sessile or slow-moving species.     

CONSEQUENCES OF ALTERNATIVE FUNCTIONAL RESPONSE FORMULATIONS IN MODELS EXPLORING WHALE-FISHERY INTERACTIONS

We evaluated the utility of Ecosim for exploring interactions between cetacean predators, their prey, and fisheries. We formulated six Ecosim parameterizations, representing alternative hypotheses of feeding interactions (functional response) between cetaceans and their main fish prey, and examined differences in the predicted responses to simulated harvesting regimes for minke whales and their prey. Regardless of the type of function response formulated, intense fishing on the main fish prey of minke whales had a longer-lasting negative impact on minke whales than when minke whale biomass was removed directly by harvesting. Consumption rate, biomass, feeding time and mortality of minke whales were all sensitive to the type of functional response specified. Inclusion of "handling time" limited minke whales consumption at high prey densities and predicted higher consumption at low prey densities; features characteristic of a type II functional response. Predicted decline and recovery rates of minke whales were slower than when consumption rates were not limited. Addition of "foraging time" adjustments resulted in more conservative estimates of decline and recovery. However, when "other mortality" was linked to time spent foraging, exposure to higher mortality at low prey densities, and reduced mortality at high prey densities resulted in dramatic differences in predicted biomass trajectory. Sensitivity to the "other mortality" assumption is important for cetaceans whose predation mortality is only a small proportion of total mortality. Differences in the feeding and biomass dynamics were also observed when prey availability to predators was represented by changes in prey vulnerability, confirming earlier reports that Ecosim predictions are sensitive to this parameter.     

Predicting predatory impact of juvenile invasive lionfish (<Emphasis Type="Italic">Pterois volitans</Emphasis>) on a crustacean prey using functional response analysis: effects of temperature,habitat complexity and light regimes

The ecological implications of biotic interactions, such as predator-prey relationships, are often context-dependent. Comparative functional responses analysis can be used under different abiotic contexts to improve understanding and prediction of the ecological impact of invasive species. Pterois volitans (Lionfish) [Linnaeus 1758] is an established invasive species in the Caribbean and Gulf of Mexico, with a more recent invasion into the Mediterranean. Lionfish are generalist predators that impact a wide range of commercial and non-commercial species. Functional response analysis was employed to quantify interaction strength between lionfish and a generic prey species, the shrimp (Paleomonetes varians) [Leach 1814], under the contexts of differing temperature, habitat complexity and light wavelength. Lionfish have prey population destabilising Type II functional responses under all contexts examined. Significantly more prey were consumed at 26 °C than at 22 °C. Habitat complexity did not significantly alter the functional response parameters. Significantly more prey were consumed under white light and blue light than under red light. Attack rate was significantly higher under white light than under blue or red light. Light wavelength did not significantly change handling times. The impacts on prey populations through feeding rates may increase with concomitant temperature increase. As attack rates are very high at low habitat complexity this may elucidate the cause of high impact upon degraded reef ecosystems with low-density prey populations, although there was little protection conferred through habitat complexity. Only red light (i.e. dark) afforded any reduction in predation pressure. Management initiatives should account for these environmental factors when planning mitigation and prevention strategies.     

How predator functional responses and Allee effects in prey affect the paradox of enrichment and population collapses

In Rosenzweig-MacArthur models of predator-prey dynamics, Allee effects in prey usually destabilize interior equilibria and can suppress or enhance limit cycles typical of the paradox of enrichment. We re-evaluate these conclusions through a complete classification of a wide range of Allee effects in prey and predator's functional response shapes. We show that abrupt and deterministic system collapses not preceded by fluctuating predator-prey dynamics occur for sufficiently steep type III functional responses and strong Allee effects (with unstable lower equilibrium in prey dynamics). This phenomenon arises as type III functional responses greatly reduce cyclic dynamics and strong Allee effects promote deterministic collapses. These collapses occur with decreasing predator mortality and/or increasing susceptibility of the prey to fall below the threshold Allee density (e.g. due to increased carrying capacity or the Allee threshold itself). On the other hand, weak Allee effects (without unstable equilibrium in prey dynamics) enlarge the range of carrying capacities for which the cycles occur if predators exhibit decelerating functional responses. We discuss the results in the light of conservation strategies, eradication of alien species, and successful introduction of biocontrol agents.     

Body masses,functional responses and predator–prey stability

The stability of ecological communities depends strongly on quantitative characteristics of population interactions (type‐II vs. type‐III functional responses) and the distribution of body masses across species. Until now, these two aspects have almost exclusively been treated separately leaving a substantial gap in our general understanding of food webs. We analysed a large data set of arthropod feeding rates and found that all functional‐response parameters depend on the body masses of predator and prey. Thus, we propose generalised functional responses which predict gradual shifts from type‐II predation of small predators on equally sized prey to type‐III functional‐responses of large predators on small prey. Models including these generalised functional responses predict population dynamics and persistence only depending on predator and prey body masses, and we show that these predictions are strongly supported by empirical data on forest soil food webs. These results help unravelling systematic relationships between quantitative population interactions and large‐scale community patterns.     

Trait‐mediated indirect interactions in a marine intertidal system as quantified by functional responses

Studies of trait‐mediated indirect interactions (TMIIs) typically focus on effects higher predators have on per capita consumption by intermediate consumers of a third, basal prey resource. TMIIs are usually evidenced by changes in feeding rates of intermediate consumers and/or differences in densities of this third species. However, understanding and predicting effects of TMIIs on population stability of such basal species requires examination of the type and magnitude of the functional responses exhibited towards them. Here, in a marine intertidal system consisting of a higher‐order fish predator, the shanny Lipophrys pholis, an intermediate predator, the amphipod Echinogammarus marinus, and a basal prey resource, the isopod Jaera nordmanni, we detected TMIIs, demonstrating the importance of habitat complexity in such interactions, by deriving functional responses and exploring consequences for prey population stability. Echinogammarus marinus reacted to fish predator diet cues by reducing activity, a typical anti‐predator response, but did not alter habitat use. Basal prey, Jaera nordmanni, did not respond to fish diet cues with respect to activity, distribution or aggregation behaviour. Echinogammarus marinus exhibited type II functional responses towards J. nordmanni in simple habitat, but type III functional responses in complex habitat. However, while predator cue decreased the magnitude of the type II functional response in simple habitat, it increased the magnitude of the type III functional response in complex habitat. These findings indicate that, in simple habitats, TMIIs may drive down consumption rates within type II responses, however, this interaction may remain de‐stabilising for prey populations. Conversely, in complex habitats, TMIIs may strengthen regulatory influences of intermediate consumers on prey populations, whilst potentially maintaining prey population stability. We thus highlight that TMIIs can have unexpected and complex ramifications throughout communities, but can be unravelled by considering effects on intermediate predator functional response types and magnitudes. Synthesis Higher‐order predators and habitat complexity can influence behaviour of intermediate species, affecting their consumption of prey through trait‐mediated indirect interactions (TMIIs). However, it is not clear how these factors interact to determine prey population stability. Using functional responses (FRs), relating predator consumption to prey density, we detected TMIIs in a marine system. In simple habitats, TMIIs reduced consumption rates, but FRs remained de‐stabilising for prey populations. In complex habitats, TMIIs strengthened prey regulation with population stabilizing FRs. We thus demonstrate that FRs can assess interactions of environmental and biological cues that result in complex and unexpected outcomes for prey populations.     

Parasite modification of predator functional response

Parasite alteration of the host (predator) functional response provides a mechanism by which parasites can alter predator–prey population dynamics and stability. We tested the hypothesis that parasitic infection of a crab (Eurypanopeus depressus) by a rhizocephalan barnacle (Loxothylacus panopei) can modify the crab’s functional response to mussel (Brachidontes exustus) prey and investigated behavioral mechanisms behind a potential change in the response. Infection dramatically reduced mussel consumption by crabs across mussel densities, resulting in a decreased attack rate parameter and a nearly eightfold reduction in maximum consumption (i.e. the asymptote, or inverse of the handling time parameter) in a type II functional response model. To test whether increased handling time of infected crabs drove the decrease in maximum consumption rate, we independently measured handling time through observation. Infection had no effect on handling time and thus could not explain the reduction in consumption. Infection did, however, increase the time that it took crabs to begin handling prey after the start of the handling time experiment. Furthermore, crabs harboring relatively larger parasites remained inactive longer before making contact with prey. This behavioral modification likely contributed to the reduced mussel consumption of infected crabs. A field survey revealed that 20 % of crabs inhabiting oyster reefs at the study site (North Inlet estuary, Georgetown, South Carolina, USA) are infected by the barnacle parasite, indicating that parasite infection could have a substantial effect on the population level crab-mussel interaction.     

Examining intraspecific multiple predator effects across shifting predator sex ratios

Predator-predator interactions, or “multiple predator effects” (MPEs), are pervasive in the structuring of communities and complicate predictive quantifications of ecosystem dynamics. The nature of MPEs is also context-dependent, manifesting differently among species, prey densities and habitat structures. However, there has hitherto been a lack of consideration for the implications of intraspecific demographic variation within populations for the strength of MPEs. The present study extends MPE concepts to examine intraspecific interactions among male and female predators across differences in prey densities using a functional response approach. Focusing on a copepod-mosquito model predator-prey system, interaction strengths of different sex ratio pairs of Lovenula raynerae were quantified towards larval Culex pipiens complex prey, with observations compared to both additive and substitutive model predictions. Copepods exhibited destabilising Type II functional responses in all treatments, with female copepods significantly more voracious than males under multiple predator groups. Lovenula raynerae exhibited significantly negative MPEs overall which resulted in prey risk reductions. However, whilst not statistically clear, the magnitude of antagonistic interactions subtly differed among predator-predator compositions and prey densities, with female-female antagonisms generally prevalent at low prey densities, and male-male negative interactions greater under high prey densities. Mixed-sex copepod group predation was predicted by both additive and substitutive models, and additive models generated significantly higher consumption estimates than substitutive equivalents given that their predictions were based on the absence of antagonistic non-trophic interactions. We propose the importance of internal population sex demographics as a further context-dependency which influences the nature of MPEs, with demographic implications requiring investigation across other taxonomic and trophic groups.     

Should "handled" prey be considered? Some consequences for functional response, predator-prey dynamics and optimal foraging theory

Predator-prey models consider those prey that are free. They assume that once a prey is captured by a predator it leaves the system. A question arises whether in predator-prey population models the variable describing prey population shall consider only those prey which are free, or both free and handled prey together. In the latter case prey leave the system after they have been handled. The classical Holling type II functional response was derived with respect to free prey. In this article we derive a functional response with respect to prey density which considers also handled prey. This functional response depends on predator density, i.e., it accounts naturally for interference. We study consequences of this functional response for stability of a simple predator-prey model and for optimal foraging theory. We show that, qualitatively, the population dynamics are similar regardless of whether we consider only free or free and handled prey. However, the latter case may change predictions in some other cases. We document this for optimal foraging theory where the functional response which considers both free and handled prey leads to partial preferences which are not observed when only free prey are considered.     

Trait-mediated interactions: influence of prey size, density and experience

1. The role of non-consumptive predator effects in structuring ecological communities has become an important area of study for ecologists. Numerous studies have shown that adaptive changes in prey in response to a predator can improve survival in subsequent encounters with that predator. 2. Prey-mediated changes in the shapes of predators' functional response surfaces determine the qualitative predictions of theoretical models. However, few studies have quantified the effects of adaptive prey responses on the shape of predator functional responses. 3. This study explores how prey density, size and previous predator experience interact to change the functional response curves of different-sized predators. 4. We use a response surface design to determine how previous exposure to small or large odonate predators affected the short-term survival of squirrel tree frog (Hyla squirella) tadpoles across a range of sizes and densities (i.e. the shape of odonate functional response curves). 5. Predator-induced tadpoles in a given size class did not differ in shape, although induction changed tadpole behaviour significantly. Induced tadpoles survived better in lethal encounters with either predator than did similar-sized predator-naive tadpoles. 6. Induction by either predator resulted in increased survival with both predators at a given size. However, different mechanisms led to increased survival for induced tadpoles. Attack rate for the small predators, whereas handling time increased for the large predators.     

Foraging behavior of the coccinellid Nephus includens (Coleoptera: Coccinellidae) in response to Aphis gossypii (Hemiptera: Aphididae) with particular emphasis on larval parasitism

This study assessed the effect of parasitism of Nephus includens (Col.: Coccinellidae) larvae by Homalotylus flaminius (Hym.: Encyrtidae) on the predation rates of the predator on the cotton aphid, Aphis gossypii Glover (Hem.: Aphididae) by deriving functional responses for second- and fourth-instar predators at prey densities ranging from 10 to 80 aphids per arena. The relationship between the functional and numerical responses of adult females of N. includens also was determined for prey densities ranging from 10 to 140 aphids per arena. Predation rates of unparasitized and parasitized second-instar N. includens were both fit by a type II functional response model with parameters as follows: unparasitized (a = 0.0768 hours(-1) and T(h) = 0.975 h) and parasitized (a = 0.0787 hours(-1) and T(h) = 0.8823 hours). Predation rates of unparasitized and parasitized fourth-instar N. includens were fit by type III and II models, respectively, with the following parameters: unparasitized (b = 0.1702 hours(-1) and T(h) = 0.2369 hours) and parasitized (a = 0.038 h(-1) and T(h) = 0.539 h). The unparasitized fourth-instar was the most voracious stage, having the highest attack rate and lowest handling time. Considering these attributes, it would seem to be the most effective stage of this predator against A. gossypii. Adult female lady beetles (N. includens) showed a type III functional response and their numerical response increase to prey density was curvilinearly related to prey density, with the highest number of eggs being produced at highest prey densities. The maximum saturation level for both predation and egg production for adult females occurred at a prey density of 120 aphids. Thus, a ratio 1:120 (predator:prey) should be used when releasing this species for augmentative biological control. Release of either fourth-instar or adult stage N. includens should be minimized the potentially negative effect of parasitism by H. flaminius on early developmental stages, and hence increase its efficiency in biocontrol programs.