27 Bout feeding

Many predators do not feed continuously over time as assumed in derivation of the vulnerable abundance Va in the standard foraging arena equation (Eq. 1 below).  Rather, they obtain most of their food intake in short, intensive feeding “bouts”, typically at dawn and dusk when light levels are changing rapidly[1] [2].

Basic foraging arena equation

[latex]V_a=\frac{v_a \cdot B_{i(a)}}{v_a+v_a' +\alpha_a \cdot P_{j(a)}}\tag{1}[/latex]

Shared foraging arena extension (see preceding chapter)

[latex]V_a=\frac{v_a \cdot B_{i(a)}}{v_a+v_a'+ \sum \limits_k \alpha_{ak} \cdot P_k}\tag{2}[/latex]

Particularly when predators such as juvenile fish have severely restricted habitat use as a tactic for managing predation risk (hiding, schooling), only a small fraction of the system-scale prey biomass is available to them in the foraging arenas that they use during each feeding bout. As an example, juvenile Atlantic salmon have been shown to restrict the time they spend feeding rather than maximizing their growth when food is abundant[3].

Here, we show that overall trophic flow rates Qak (for predators k feeding in foraging arena a) over longer time scales can still be closely approximated by a continuous rate equation of the mass-action form Qak = αakVi(a)Pk(a), where the bout search rates α* and mean vulnerable prey densities per feeding bout are comparable to (but differ numerically from) the α,V predictions for continuous feeding.

Consider a single feeding bout in arena a of duration d (d << one day), during which an initial prey density Va(0) is depleted by predators k(a).  Assume that d is short enough that prey renewal and loss during the bout, (e.g., due to prey spatial movement and other mortality sources) can be safely ignored.  Assume that Va(0) is a proportion fa of total prey biomass Bi(a) and that renewal mechanisms between bouts make Va(0) = fa Bi(a).  Note that when used over multiple bouts, this prediction of Va(0) for each bout requires that arena prey abundance be independent of predation effects in previous bouts except through effects on Bi(a), i.e. that there are no carryover effects from previous bouts (extreme opposite of continuous feeding assumption).  Then if predators k search randomly within the arena, vulnerable prey density Va(t) will change during the bout according to the simple rate equation,

[latex]\frac{dV_a(t)}{dt}=-V_a(t)\sum \limits_k\alpha_{ak}P_k\tag{3}[/latex]

where the αak are predator rates of effective search with the same interpretation as for continuous feeding.

Integrating Eq. 3 over the bout duration d leads to the familiar exponential exploitation equation Va(d) = Va(0) exp(-d ∑k αak Pk) and to predicted total prey consumption per bout Qakbout by each predator k,

[latex]Q_{ak}^{bout}=\frac{\alpha_{ak}P_k}{\sum\limits_k\alpha_{ak}P_k}\cdot f_aB_{i(a)}\cdot [1-\exp(-d\sum\limits_k\alpha_{ak}P_k)]\tag{4}[/latex]

The first term of Eq. 4 simply apportions total prey consumption Va(0) – Va(d) over the bout among competing predators.  Further, the mean prey density Va* during the bout is given by the integral of V over the bout divided by bout duration d. This mean is just,

[latex]V_a^*=f_aB_{i(a)}\cdot \frac{1-\exp(-d\sum\limits_k\alpha_{ak}P_k)}{d\sum\limits_k\alpha_{ak}P_k}\tag{5}[/latex]

Expressed in terms of this mean arena prey density, consumption per bout Eq. 4 can be expressed more simply as

[latex]Q_{ak}^{bout}=\alpha_{ak}P_kV_a^*\tag{6}[/latex]

We could use this formula directly in a complex simulation model that steps forward in time by the interval Δt between feeding bouts, adding in other components of prey and predator abundance change over each such short interval.  Fortunately, such a tedious calculation is generally unnecessary.

Consider the component of overall prey biomass change caused by each feeding bout, where there are nb = 1/Δt bouts per year.  That (typically small) change in Bi(a) per bout is given by the sum of Eq. 4 terms over predators k, i.e.,

[latex]\Delta B_{i(a)} = f_aB_{i(a)}\cdot[1-\exp(-d\sum\limits_k\alpha_{ak}P_k)]\tag{7}[/latex]

Dividing this by the bout duration Δt gives a discrete-time component of the prey rate of change,

[latex]\frac {\Delta B_{i(a)}}{\Delta t}=\frac{1}{\Delta t}f_aB_{i(a)}\cdot [1-\exp(-d\sum\limits_k\alpha_{ak}P_k)] =n_bf_aB_{i(a)}\cdot [1-\exp(-d\sum\limits_k\alpha_{ak}P_k)]\tag{8}[/latex]

Since the time Δt between bouts is typically very short (nb is typically of the order of several hundred bouts per year), we can approximate Eq. 8 very accurately by treating it as a continuous rate component dBi(a)/dt.  This approximation leads immediately to a continuous rate equation for Qak comparable to the continuous feeding case where Qak = αakPkVa, namely,

[latex]Q_{ak}=\alpha_{ak}^* \ P_k \ v_a^* \ B_{i(a)}\cdot \frac{1-\exp(-\sum\limits_k\alpha_{ak}^*P_k)}{\sum\limits_k\alpha_{ak}^*P_k}=\alpha_{ak}^* \ P_k \ V_a^*\tag{9}[/latex]

Where αak* are the duration-weighted search rates αak* = αakd. va* = nb fa represents a total prey “fraction” that would become vulnerable over a one-year time scale, and  Va* (comparable to Eq 2) is given by

[latex]V_a^*=v_a^* \ B_{i(a)} \cdot \frac{1-\exp(-\sum\limits_k\alpha_{ak}^* \ P_k)}{\sum\limits_k\alpha_{ak}^* \ P_k}\tag{10}[/latex]

This model for vulnerable prey density obviously exhibits the same “ratio dependence” of available prey density on predator abundance as does Eq 2, but with the ratio effect 1/(v+v’+∑k αak Pk) replaced by a negative exponential effect.   At high predator abundances it also implies an upper bound Bi(a) on total removal rate Qa and hence on total instantaneous predation mortality rate Qa/Bi(a).

We can parameterize the continuous approximation to bout feeding from Ecopath inputs and assumed maximum predation rates in the same way as described in the previous section for continuous arena feeding.  That is, we set va* = Ka Ma(0) where Ka as above is a defined input ratio of maximum to Ecopath baseline predation rate.  We calculate base mean prey density per bout Va*(0) by substituting Ecopath base prey and predator abundances Bi(a)(0) and Pk(0) into Eq. 10 along with ∑k α Pk(0) = Qa(0) / Va(0), (where Qa(0) is the base total consumption rate summed over predators k), and solving for Va*, to give

[latex]V_{a}^*(0)=-\frac{Q_a(0)}{\ln(1-\frac{1}{K_a})}\tag{11}[/latex]

Then we simply calculate the αak* as

[latex]\alpha_{ak}^*=-\frac{Q_{ak}(0)}{P_k(0)\cdot V_a^*(0)}\tag{12}[/latex]

whereas above the arena-specific base consumption rate is calculated using assumed arena feeding proportions pak as Qak(0) = pakQi(a),k(0), and Qi(a),k(0) is the Ecopath base total consumption rate of prey i(a) by predator k.

 

Figure 1. Comparison of instantaneous mortality rates expressed relative to Ecopath baseline predation mortality rates for the original model formulation (“Continuous”) and compared to bout feeding with va* = va in case Bout A, and with va* < va set to give same limiting maximum consumption per predator, Q/P in case Bout B.

It is instructive to compare the predictions of instantaneous prey mortality rate M = Qa/Bi(a) from Eq. 8 to those of the continuous model defined by Eq. 2  (and Qa = ∑k αak Pk Va), for varying predator abundances Pk while holding prey biomasses Bi(a) constant (Figure 1). If we set va* = va, i.e. use the same Ka to calculate va* as we would for va in the continuous case (Bout A in Figure 1), the exponential term in the bout feeding model generally predicts steeper variation in M than the continuous model, i.e., it predicts that M will drop off more rapidly if P decreases from P(0) than does the continuous model.  This leads to weaker “compensation” measured in terms of increase in potential Q/P as P declines.  But if we set va* smaller than va, so as to predict the same limiting maximum consumption per predator (Q/P) at very low predator densities (Bout B in Figure 1), the two arena models give predicted patterns of variation in M that are the opposite, i.e. bout feeding predicts saturation of M at lower P than the continuous case.  This means that in Ecosim cases where Ka has been estimated by fitting the continuous arena model to time series data (the only option before the inclusion of bout feeding in EwE), and where feeding in reality has been of the bout type, the fitted Ka estimates have probably been somewhat too large, i.e Ka is in reality closer to 1.0 and predators are already causing (in the Ecopath base situation) what may be close to their maximum possible predation rates from bout feeding.

For Ecosim models that include multi-stanza population dynamics, a critically important capability is to represent adjustments in foraging time, particularly for juvenile stanzas.  Such adjustments allow juvenile fish to translate increases in potential feeding rate Q/P into reduced foraging time and predation risk when competitor abundance P decreases[4], leading to compensatory changes in juvenile mortality rates and emergent stock-recruitment relationships of the Beverton-Holt form[5].

Foraging time adjustments are modeled in Ecosim by including a dynamic variable Ti for each biomass type, with Ti at time zero set to 1.0.  Then Ti is varied over time so as to try and maintain Ecopath base feeding rate per predator (Q/P), by multiplying all search rates of type i for its prey by Ti, and all vulnerability exchange rates of type i into arenas where predators take it by Ti.  In the bout foraging representation, this means simply that (1) search parameters for type i as a predator are adjusted by varying bout durations d in proportion to Ti (i.e. setting α*(t) = α*(0) Ti, with Ti defined as the relative bout duration d(t)/d(0)) and (2) the vulnerable fraction f that define v* of i to its predators are also treated as being proportional to d by setting v*(t) = v*(0) Ti.

As a simple test of whether bout feeding is likely to make much difference to the ability of Ecosim models to fit historical time series data, we examined changes in a simple fitting criterion (sum of squared deviations, SS, from historical data, summed over all time series used in model fitting) for a collection of models that had been fitted to data using the continuous arena equations, when all trophic linkages were simply reset to assume bout feeding without correction or refitting of the Ka parameters (Figure 2,[6]). Ability of several of these models to fit historical data are reviewed in Walters and Martell (2004, Figure 12.6[7]), and these were mostly the same models used in single-species versus multispecies MSY comparisons by Walters et al.[8].

Figure 2. Changes in sum of squares goodness-of-fit criterion (SS) for an assortment of models that have been fit to historical data using the continuous feeding arena equations of Ecosim, when bout feeding is assumed instead for all trophic linkages initially without re-estimation of vulnerability exchange multipliers Ka.  (A): represents SS from original model.  (B): SS from original model but with bout feeding.  (C): SS for bout feeding, after fitting by varying 20 most important Ka.  (D): SS for continuous feeding after fitting.

Surprisingly, there was little change in the fitting criterion for many of the models, and one (Central North Pacific) even gave a better fit immediately.  For those models where there was a substantial increase in SS, it was easy to remedy the poor fits by re-estimating Ka under the global bout arena assumption.  When we refitted the models under both feeding assumptions (by nonlinear estimation search over the 20 Ka values with largest contributions to the sum of squares), we were easily able to find fits at least as good under the bout feeding assumption for most cases, and qualitatively as good for all cases.

Discussion

The equations introduced here obviously give considerable flexibility to represent trophic interactions in the Ecosim model more realistically than previously possible[9].  It is particularly comforting to see that the much more realistic assumption of bout rather than continuous feeding leads to very similar predictions of how prey mortality rates should vary with predator abundances as have been assumed in previous Ecosim models based on the unrealistic but mathematically convenient assumption of continuous feeding with rapid equilibration of vulnerable prey densities.

We recommend extreme care in using either the continuous or bout feeding equations to represent feeding by multiple predators in a relatively small number of arenas.  As noted above, the intense inter-specific competition implied by such concentration of feeding has very likely driven natural selection for differentiation in feeding behavior (use of different fine-scale arenas) as well as in diet composition.  See for instance Berec et al. (2006) for an experiment illustrating this. If such differentiation is excluded from the model parameterization, the Ecosim user risks building a model that will not retain observed biodiversity over time.

A few authors have referred to the basic Ecosim equation for predicting total flow rates Q =(α v B P) / (v + v’ + α P) as though it were a functional response equation comparable to assuming mass-action encounters and type II predation, e.g., Q =(α B P) / (1 + h B);[10].  Such comparisons reflect a misunderstanding about a basic proposition of foraging arena theory, namely that predators very generally encounter their prey in space-time restricted circumstances (foraging arenas), such that it is almost never appropriate to predict Q from the ecosystem-scale mean prey density B when trying to account for effects such as handling time and switching (changes in α).  We would argue that it is sometimes appropriate to account for handling time effects, but only if these are predicted using arena-scale vulnerable prey densities V, i.e. Q = α V P / (1 + hV), where V is adjusted away from the system-scale average B using assumptions about localization of foraging (effects of vulnerability exchanges v’s and/or available prey fractions per bout f’s in the arena equations). We explicitly allow switching in Ecosim, but again caution that it should be used in conjunction with predictions of vulnerable, rather than overall, prey densities.

It would be ignorant to assert that the equations presented in this chapter are the only or best way to represent differentiation of vulnerable prey biomasses V from system-scale average prey biomasses B in prediction of trophic interaction rates.  They do not for example account explicitly for some very gross system-scale effects that occur in highly disturbed systems, such as changes in overall prey and predator distributions and overlap patterns, (e.g., due to range contractions), and changes in spatial arena structure due to obvious habitat changes like growth and destruction of biogenic spatial refuges[11].  Some such changes can be accounted for in Ecosim through trophic ‘mediation functions’ that link v’s and α’s to abundances of species besides those engaged directly as predators and prey, (e.g., one can make α’s and v’s for juvenile fish that hide in macrophyte beds dependent on macrophyte biomass).  But there is still a long way to go in development of fully-defensible predictions of V for systems that are massively disturbed.

One option for dealing with the prediction of V would be to construct very detailed spatial models (with habitat and its use modelled at fine scale, maybe of a few m2) running on very short (bout) time scales (time steps of one hour or less).  But such models may be plagued by lack of detailed spatial data, lack of understanding of how organisms move and concentrate their activities at such fine scales, and risk of cumulative divergence of predictions from reality simply due to explosions over simulated time and space of small errors in behavioural movement predictions.

A key advantage of the relatively simple foraging arena equations for Q prediction is that we can easily force them to agree with baseline “observations” or estimates of system-scale abundances (B’s, P’s), feeding rates, and diet compositions (Qji) as summarized in static (point-in-time) mass-balance assessments like Ecopath.  But this is also a disadvantage, in the sense that the rate parameter estimates then become dependent on the often incomplete and possibly biased estimates entered as Ecopath inputs.  It is clear that Ecosim-type dynamic predictions are sensitive to those baseline inputs, and that this represents an especially severe issue for interactions involving small fish as prey (where the small fish typically represent only trivial and often overlooked proportions of their predators’ diets).

Even absent difficulties with Ecopath inputs, i.e. empirical knowledge of baseline ecosystem biomass flow rates and states, the most troublesome parameters for Ecosim users to specify have been the “vulnerability multipliers” Ka representing ratios of maximum to Ecopath base predation mortality rates.  One source of trouble is obviously that Ecopath inputs provide no information about the Ka, and such information can only come from either fine-scale analysis of spatial arena structures, from data collected at different times and/or places about how Q’s have varied with predator and prey abundances, or from assumptions about or estimates of where populations are relative to their carrying capacity.  Indeed, this is why we emphasize the importance of fitting Ecosim models to time series data by varying the Ka parameters.

Another, and important aspect is that the Ka are not purely “behavioural” or ecological parameters; rather, they depend as well on how large the Ecopath initial predator abundances Pk are compared to what the ecosystem might naturally support (see Density dependence chapter).  So for example a model that includes Atlantic cod stocks off Newfoundland, and uses the current low stock size as the Ecopath base, must have very high Ka values (1000+) for interactions between cod and its prey, else the model will not make enough prey available to the simulated cod stock for it to recover to anywhere near its historical abundance when simulated fishing is removed.

We can provide some guidance about reasonable ecological Ka values (corrected for effects of historical depletion on biomasses) from meta-analysis of Ka estimates for many fitted models.  One pattern that is becoming broadly evident from cases like those in Figure 2 is that fitted Ka values tend to be small (<2.0) for most trophic linkages in temperate and tropical systems, and for feeding by juvenile stanzas in all systems.  In contrast, fitted Ka values tend to be much larger for most interactions (except juvenile stanzas of demersal fish species) in high-latitude ecosystems like the Bering Sea.  The low Ka values are easily explained for juvenile fish and reef-associated older fish, as a consequence of severe spatial restriction in habitat use leading to low proportions of prey populations being available the fish at any time[12].  High Ka values in high-latitude systems likely reflect the wider spatial movement characteristic of northern fish, and tactics such as diel vertical migration that bring high proportions of widely distributed predators and prey into daily contact with one another, (i.e., high f’s for bout feeding during periods of diurnal contact[13]).

The Ka vulnerability multipliers are discussed in details in the vulnerability multiplier chapter.

Attribution This chapter is based on Walters and Christensen (2007)[14], used with permission from Elsevier, Licence Numbers 5663310244809 and 5663310474242.


  1. Helfman, G. S., 1993. Fish behaviour by day, night and twilight. In: T. J. Pitcher (Editor) Behaviour of Teleost Fishes. Chapman & Hall, London, Vol. 2. pp. 479-512.
  2. Rickel, S. and Genin, A., 2005. Twilight transitions in coral reef fish: the input of light-induced changes in foraging behaviour. Animal Behaviour, 70:133-144. https://doi.org/10.1016/j.anbehav.2004.10.014
  3. Orpwood, J. E., Griffiths, S. W. and Armstrong, J. D., 2006. Effects of food availability on temporal activity patterns and growth of Atlantic salmon. Journal of Animal Ecology, 75:677-685. https://doi.org/10.1111/j.1365-2656.2006.01088.x
  4. see, e.g., Orpwood, J. E., Griffiths, S. W. and Armstrong, J. D., 2006. Effects of food availability on temporal activity patterns and growth of Atlantic salmon. Journal of Animal Ecology, 75:677-685. https://doi.org/10.1111/j.1365-2656.2006.01088.x
  5. Walters, C. and Korman, J., 1999. Linking recruitment to trophic factors: revisiting the Beverton-Holt recruitment model from a life history and multispecies perspective. Reviews in Fish Biology and Fisheries, 9:187-202. https://doi.org/10.1023/A:1008991021305
  6. for details about the models, see Walters, C and V. Christensen. 2007. Adding realism to foraging arena predictions of trophic flow rates in Ecosim ecosystem models: shared foraging arenas and bout feeding. Ecological Modelling 209:342-350. https://doi.org/10.1016/j.ecolmodel.2007.06.025
  7. Walters, C. J. and Martell, S. J. D., 2004. Fisheries ecology and management. Princeton University Press, Princeton. 399 pp.
  8. Walters, C. J., Christensen, V., Martell, S. J. and Kitchell, J. F., 2005. Possible ecosystem impacts of applying MSY policies from single-species assessment. ICES Journal of Marine Science, 62:558-568. https://doi.org/10.1016/j.icesjms.2004.12.005
  9. Walters, C., Pauly, D., Christensen, V. and Kitchell, J. F., 2000. Representing density dependent consequences of life history strategies in aquatic ecosystems: EcoSim II. Ecosystems, 3:70-83. https://doi.org/10.1007/s100210000011
  10. see, e.g., Koen-Alonso, M. and Yodzis, P., 2005. Multispecies modelling of some components of the marine community of northern and central Patagonia, Argentina. Canadian Journal of Fisheries and Aquatic Sciences, 62:1490-1512. https://doi.org/10.1139/f05-087
  11. e.g., Rodriguez, C. F., Becares, E., Fernandez-Alaez, M. and Fernandez-Alaez, C., 2005. Loss of diversity and degradation of wetlands as a result of introducing exotic crayfish. Biological Invasions, 7:75-85. https://doi.org/10.1007/s10530-004-9636-7
  12. e.g., Gonzalez, M. and Tessier, A., 1997. Habitat segregation and interactive effects of multiple predators on a prey assemblage. Freshwater Biology, 38:179-191.
  13. see, e.g., Hrabik, T. R., Jensen, O. P., Martell, S. J. D., Walters, C. J. and Kitchell, J. F., 2006. Diel vertical migration in the Lake Superior pelagic community. I. Changes in vertical migration of coregonids in response to varying predation risk. Canadian Journal of Fisheries and Aquatic Sciences, 63:2286-2295. https://doi.org/10.1139/f06-12
  14. Walters, C and V. Christensen. 2007. Adding realism to foraging arena predictions of trophic flow rates in Ecosim ecosystem models: shared foraging arenas and bout feeding. Ecological Modelling 209:342-350. https://doi.org/10.1016/j.ecolmodel.2007.06.025

License

Icon for the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License

Ecosystem Modelling with EwE Copyright © 2024 by Ecopath International Initiative is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License, except where otherwise noted.

Share This Book