{"id":1912,"date":"2023-11-24T16:42:44","date_gmt":"2023-11-24T21:42:44","guid":{"rendered":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/?post_type=chapter&#038;p=1912"},"modified":"2025-10-30T08:56:52","modified_gmt":"2025-10-30T12:56:52","slug":"predator-satiation-and-foraging-time","status":"publish","type":"chapter","link":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/predator-satiation-and-foraging-time\/","title":{"raw":"Predator satiation and foraging time","rendered":"Predator satiation and foraging time"},"content":{"raw":"<h2>Predator satiation and handling time effects<\/h2>\r\nEcosim and Ecospace allow you to represent two factors that may limit prey consumption rates per predator (<em>Q\/B<\/em>):\r\n<ol>\r\n \t<li>foraging time adjustments related to predation risk and\/or satiation; and<\/li>\r\n \t<li>handling time effects.<\/li>\r\n \t<li>Parameters for both are specified via the Ecosim Group info form.<\/li>\r\n<\/ol>\r\nSatiation and\/or choices to forage for short times in order to avoid higher predation risk are represented by setting non-zero values for the <em>Ecosim &gt; Input &gt; Group info &gt; Feeding time adjustment rate<\/em> of a group: larger values of this rate represent more rapid adjustment of foraging time. Non-zero foraging time adjustment rates cause Ecosim\/Ecospace to update relative foraging time during each simulation so as to represent predators as trying to maintain <em>Q\/B<\/em> near the Ecopath input base rate. For some organisms (particularly marine mammals) this foraging time adjustment may represent animals always trying to feed to satiation (<em>Q\/B<\/em> from Ecopath the satiation feeding rate) and taking more or less time to reach satiation depending on prey densities (and possibly also facing higher predation risk when foraging times are longer). For other organisms, the Ecopath base <em>Q\/B<\/em> may represent a much lower feeding rate than the animal could achieve under \"safe\" laboratory conditions, and in this case we view the base <em>Q\/B<\/em> as an evolutionary \"target\" rate representing results of natural selection for balancing benefits from feeding with predation risk costs of spending more time feeding.\r\n\r\nHandling time effects represent the notion that predators have limited time available for foraging and this time can be used up by \"handling time\" (pursuit\/manipulation\/ingestion time per prey captured) rather than searching for prey, when prey densities are high. The <em>Ecosim &gt; Input &gt; Group info &gt; QB<sub>max<\/sub>\/QB<sub>0<\/sub> (for handling time) (&gt;1)<\/em> parameter allows you to set ratios of maximum to Ecopath base food consumption rates per individual (or per biomass). These ratios are set to large values (1000) by default, which allows predators to increase their feeding rates without limit as prey densities increase (i.e., not limited by time required to handle each prey). In most scenarios, limitation of prey vulnerability prevents this unreasonable assumption from having noticeable effect. But in scenarios where vulnerable prey densities of at least one type do increase greatly, setting a low value (e.g., 2 or 3) for the predator\u2019s maximum\/base feeding rate ratio allows you to represent limits on feeding rate associated\u00a0with time needed to handle each prey. Without such limits, your predictions of increase in predator <em>Q\/B<\/em>, and hence productivity, at low predator density (or high prey density) might be too optimistic and lead you to errors like overestimating sustainable harvest rate for the predator. Also, ignoring handling time effects when one prey type increases greatly can cause an underestimate of the \u2018buffering\u2019 effect that such increases can have on predation rates felt by other prey: if the predator consumes more of the abundant prey, and spends more time handling\/resting because of this, predation rates on other prey species should decrease.\r\n\r\nEcosim\/Ecospace calculates feeding rates of predators using the \"multispecies disc equation\", a generalization of Holling\u2019s type II functional response model for multiple prey types (see the <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/functional-response\/\">Holling functional response<\/a> chapter). Using the maximum\/base feeding rate ratio <em>R<\/em><sub>j<\/sub> from the <em>Ecosim &gt; Input &gt; Group info<\/em> form along with the Ecopath base food consumption rate per predator, the program calculates a maximum ration and effective handling time per prey biomass eaten (handling time = 1 \/ (maximum prey biomass eaten per time)). This handling time (Holling\u2019s <em>h<\/em>\u00a0parameter) is used to calculate the denominator in the disc equation formulation <em>Q<sub>ij<\/sub>\/B<sub>j<\/sub><\/em> , biomass of prey type <em>i<\/em> consumed per time per unit biomass of predator <em>j<\/em>, as <em>Q<sub>ij<\/sub>\/B<sub>j<\/sub><\/em>\u00a0 = <em>a<sub>ij<\/sub> V<sub>ij<\/sub> \/ (1 + h<sub>j<\/sub> \u2211<sub>k<\/sub> a<sub>kj<\/sub> V<sub>kj<\/sub>)<\/em> where <em>a<sub>ij<\/sub><\/em> is the rate of effective search by predator <em>j<\/em> for type <em>i<\/em> prey, <em>h<sub>j<\/sub><\/em> is the predator handling time parameter, <em>V<sub>ij<\/sub><\/em> is the instantaneous density of prey type <em>i<\/em> vulnerable to predator <em>j, <\/em>and the sum in the denominator is over all prey types k taken by the predator. A useful fact about the multispecies disc equation is that <em>D<\/em><sub>j<\/sub>, the proportion of time spent feeding (reactive to prey rather than handling), is given by <em>D<\/em>j=1\/<em>(1 + h<sub>j<\/sub> \u2211<sub>k<\/sub> a<sub>kj<\/sub> V<sub>kj<\/sub>).<\/em> For more information about how the disc equation <em>D<\/em><sub>j<\/sub> enters food consumption rate calculations along with other factors that influence feeding, see the <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/foraging-arena-theory\/\">Foraging arena<\/a> chapter. The solution for vulnerable prey densities <em>V<sub>ij<\/sub><\/em> needed in the disc equation calculation over time involves\u00a0 a numerical procedure that can occasionally cause annoying \"chatter\" in the Ecosim results when handling times are large (ratio of maximum\/base consumption rate small).\r\n\r\nA helpful fact about the <em>D<sub>j<\/sub><\/em> proportion of time spent feeding in the disc equation formulation is that it can be calculated simply from the user-provided ratio <em>R<\/em><sub>j<\/sub> of maximum to ecopath base ration, as just D<sub>j<\/sub>=R<sub>j<\/sub>\/(R<sub>j<\/sub>-1). This is used to initialize the <em>D<\/em><sub>j<\/sub> at the start of each Ecosim run.\r\n\r\nBioenergetics models for fish most often indicate that feeding rates are low compared to maximum ration; typical ratios of estimated to maximum ration (Hewett-Johnson[footnote]Hewett SW and BL Johnson. 1992. Fish bioenergetics model 2. Univ. of Wisconsin Sea Grant Institute. <a href=\"https:\/\/repository.library.noaa.gov\/view\/noaa\/35468\">https:\/\/repository.library.noaa.gov\/view\/noaa\/35468<\/a>[\/footnote] <em>P<\/em>\u00a0parameter) are around 0.3-0.4. These estimates imply <em>R<sub>j<\/sub><\/em> (maximum\/Ecopath base ration) values of at least 2-4. If you choose to use such realistic values instead of the default 1000, and if this causes Ecosim\/Ecospace to show oscillatory behaviour, you need to consider two possibilities:\r\n<ul>\r\n \t<li>The oscillatory behaviour may be a numerical artifact of the procedure used to update <i>D<\/i><em><sub>j<\/sub><\/em>; or<\/li>\r\n \t<li>The model\u2019s \"correct\" behaviour for the parameter combinations you have provided is indeed a predator-prey cycle.<\/li>\r\n<\/ul>\r\nIf the oscillation has a period of several time steps (months), it is very likely a predator-prey cycle. Persistent predator-prey cycles are commonly predicted by models that include handling time, along with strong top-down control (high vulnerability multipliers <em>v<sub>ij<\/sub><\/em> of prey to predators). If you think the cycle is unrealistic, you should adjust the prey vulnerabilities multipliers (<em>Ecosim &gt; Input &gt; Vulnerabilities<\/em>) to lower values (toward \"bottom up\", prey vulnerability control) rather than just setting high <em>R<sub>j<\/sub><\/em> values. If you see very short cycles indicating numerical instability in the <i>D<\/i><em><sub>j<\/sub><\/em> adjustment procedure (usually happens for fast turnover groups like micro-zooplankton), you should set higher <em>R<sub>j<\/sub><\/em> values for the offending groups. This amounts to recognizing that Ecosim may be limited in its ability to represent very fast dynamic changes in groups that turn over very rapidly.","rendered":"<h2>Predator satiation and handling time effects<\/h2>\n<p>Ecosim and Ecospace allow you to represent two factors that may limit prey consumption rates per predator (<em>Q\/B<\/em>):<\/p>\n<ol>\n<li>foraging time adjustments related to predation risk and\/or satiation; and<\/li>\n<li>handling time effects.<\/li>\n<li>Parameters for both are specified via the Ecosim Group info form.<\/li>\n<\/ol>\n<p>Satiation and\/or choices to forage for short times in order to avoid higher predation risk are represented by setting non-zero values for the <em>Ecosim &gt; Input &gt; Group info &gt; Feeding time adjustment rate<\/em> of a group: larger values of this rate represent more rapid adjustment of foraging time. Non-zero foraging time adjustment rates cause Ecosim\/Ecospace to update relative foraging time during each simulation so as to represent predators as trying to maintain <em>Q\/B<\/em> near the Ecopath input base rate. For some organisms (particularly marine mammals) this foraging time adjustment may represent animals always trying to feed to satiation (<em>Q\/B<\/em> from Ecopath the satiation feeding rate) and taking more or less time to reach satiation depending on prey densities (and possibly also facing higher predation risk when foraging times are longer). For other organisms, the Ecopath base <em>Q\/B<\/em> may represent a much lower feeding rate than the animal could achieve under &#8220;safe&#8221; laboratory conditions, and in this case we view the base <em>Q\/B<\/em> as an evolutionary &#8220;target&#8221; rate representing results of natural selection for balancing benefits from feeding with predation risk costs of spending more time feeding.<\/p>\n<p>Handling time effects represent the notion that predators have limited time available for foraging and this time can be used up by &#8220;handling time&#8221; (pursuit\/manipulation\/ingestion time per prey captured) rather than searching for prey, when prey densities are high. The <em>Ecosim &gt; Input &gt; Group info &gt; QB<sub>max<\/sub>\/QB<sub>0<\/sub> (for handling time) (&gt;1)<\/em> parameter allows you to set ratios of maximum to Ecopath base food consumption rates per individual (or per biomass). These ratios are set to large values (1000) by default, which allows predators to increase their feeding rates without limit as prey densities increase (i.e., not limited by time required to handle each prey). In most scenarios, limitation of prey vulnerability prevents this unreasonable assumption from having noticeable effect. But in scenarios where vulnerable prey densities of at least one type do increase greatly, setting a low value (e.g., 2 or 3) for the predator\u2019s maximum\/base feeding rate ratio allows you to represent limits on feeding rate associated\u00a0with time needed to handle each prey. Without such limits, your predictions of increase in predator <em>Q\/B<\/em>, and hence productivity, at low predator density (or high prey density) might be too optimistic and lead you to errors like overestimating sustainable harvest rate for the predator. Also, ignoring handling time effects when one prey type increases greatly can cause an underestimate of the \u2018buffering\u2019 effect that such increases can have on predation rates felt by other prey: if the predator consumes more of the abundant prey, and spends more time handling\/resting because of this, predation rates on other prey species should decrease.<\/p>\n<p>Ecosim\/Ecospace calculates feeding rates of predators using the &#8220;multispecies disc equation&#8221;, a generalization of Holling\u2019s type II functional response model for multiple prey types (see the <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/functional-response\/\">Holling functional response<\/a> chapter). Using the maximum\/base feeding rate ratio <em>R<\/em><sub>j<\/sub> from the <em>Ecosim &gt; Input &gt; Group info<\/em> form along with the Ecopath base food consumption rate per predator, the program calculates a maximum ration and effective handling time per prey biomass eaten (handling time = 1 \/ (maximum prey biomass eaten per time)). This handling time (Holling\u2019s <em>h<\/em>\u00a0parameter) is used to calculate the denominator in the disc equation formulation <em>Q<sub>ij<\/sub>\/B<sub>j<\/sub><\/em> , biomass of prey type <em>i<\/em> consumed per time per unit biomass of predator <em>j<\/em>, as <em>Q<sub>ij<\/sub>\/B<sub>j<\/sub><\/em>\u00a0 = <em>a<sub>ij<\/sub> V<sub>ij<\/sub> \/ (1 + h<sub>j<\/sub> \u2211<sub>k<\/sub> a<sub>kj<\/sub> V<sub>kj<\/sub>)<\/em> where <em>a<sub>ij<\/sub><\/em> is the rate of effective search by predator <em>j<\/em> for type <em>i<\/em> prey, <em>h<sub>j<\/sub><\/em> is the predator handling time parameter, <em>V<sub>ij<\/sub><\/em> is the instantaneous density of prey type <em>i<\/em> vulnerable to predator <em>j, <\/em>and the sum in the denominator is over all prey types k taken by the predator. A useful fact about the multispecies disc equation is that <em>D<\/em><sub>j<\/sub>, the proportion of time spent feeding (reactive to prey rather than handling), is given by <em>D<\/em>j=1\/<em>(1 + h<sub>j<\/sub> \u2211<sub>k<\/sub> a<sub>kj<\/sub> V<sub>kj<\/sub>).<\/em> For more information about how the disc equation <em>D<\/em><sub>j<\/sub> enters food consumption rate calculations along with other factors that influence feeding, see the <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/foraging-arena-theory\/\">Foraging arena<\/a> chapter. The solution for vulnerable prey densities <em>V<sub>ij<\/sub><\/em> needed in the disc equation calculation over time involves\u00a0 a numerical procedure that can occasionally cause annoying &#8220;chatter&#8221; in the Ecosim results when handling times are large (ratio of maximum\/base consumption rate small).<\/p>\n<p>A helpful fact about the <em>D<sub>j<\/sub><\/em> proportion of time spent feeding in the disc equation formulation is that it can be calculated simply from the user-provided ratio <em>R<\/em><sub>j<\/sub> of maximum to ecopath base ration, as just D<sub>j<\/sub>=R<sub>j<\/sub>\/(R<sub>j<\/sub>-1). This is used to initialize the <em>D<\/em><sub>j<\/sub> at the start of each Ecosim run.<\/p>\n<p>Bioenergetics models for fish most often indicate that feeding rates are low compared to maximum ration; typical ratios of estimated to maximum ration (Hewett-Johnson<a class=\"footnote\" title=\"Hewett SW and BL Johnson. 1992. Fish bioenergetics model 2. Univ. of Wisconsin Sea Grant Institute. https:\/\/repository.library.noaa.gov\/view\/noaa\/35468\" id=\"return-footnote-1912-1\" href=\"#footnote-1912-1\" aria-label=\"Footnote 1\"><sup class=\"footnote\">[1]<\/sup><\/a> <em>P<\/em>\u00a0parameter) are around 0.3-0.4. These estimates imply <em>R<sub>j<\/sub><\/em> (maximum\/Ecopath base ration) values of at least 2-4. If you choose to use such realistic values instead of the default 1000, and if this causes Ecosim\/Ecospace to show oscillatory behaviour, you need to consider two possibilities:<\/p>\n<ul>\n<li>The oscillatory behaviour may be a numerical artifact of the procedure used to update <i>D<\/i><em><sub>j<\/sub><\/em>; or<\/li>\n<li>The model\u2019s &#8220;correct&#8221; behaviour for the parameter combinations you have provided is indeed a predator-prey cycle.<\/li>\n<\/ul>\n<p>If the oscillation has a period of several time steps (months), it is very likely a predator-prey cycle. Persistent predator-prey cycles are commonly predicted by models that include handling time, along with strong top-down control (high vulnerability multipliers <em>v<sub>ij<\/sub><\/em> of prey to predators). If you think the cycle is unrealistic, you should adjust the prey vulnerabilities multipliers (<em>Ecosim &gt; Input &gt; Vulnerabilities<\/em>) to lower values (toward &#8220;bottom up&#8221;, prey vulnerability control) rather than just setting high <em>R<sub>j<\/sub><\/em> values. If you see very short cycles indicating numerical instability in the <i>D<\/i><em><sub>j<\/sub><\/em> adjustment procedure (usually happens for fast turnover groups like micro-zooplankton), you should set higher <em>R<sub>j<\/sub><\/em> values for the offending groups. This amounts to recognizing that Ecosim may be limited in its ability to represent very fast dynamic changes in groups that turn over very rapidly.<\/p>\n<hr class=\"before-footnotes clear\" \/><div class=\"footnotes\"><ol><li id=\"footnote-1912-1\">Hewett SW and BL Johnson. 1992. Fish bioenergetics model 2. Univ. of Wisconsin Sea Grant Institute. <a href=\"https:\/\/repository.library.noaa.gov\/view\/noaa\/35468\">https:\/\/repository.library.noaa.gov\/view\/noaa\/35468<\/a> <a href=\"#return-footnote-1912-1\" class=\"return-footnote\" aria-label=\"Return to footnote 1\">&crarr;<\/a><\/li><\/ol><\/div>","protected":false},"author":1909,"menu_order":9,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-1912","chapter","type-chapter","status-publish","hentry"],"part":411,"_links":{"self":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapters\/1912","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/wp\/v2\/users\/1909"}],"version-history":[{"count":21,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapters\/1912\/revisions"}],"predecessor-version":[{"id":3693,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapters\/1912\/revisions\/3693"}],"part":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/parts\/411"}],"metadata":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapters\/1912\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/wp\/v2\/media?parent=1912"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapter-type?post=1912"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/wp\/v2\/contributor?post=1912"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/wp\/v2\/license?post=1912"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}