{"id":3856,"date":"2024-09-25T20:30:06","date_gmt":"2024-09-26T00:30:06","guid":{"rendered":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/?post_type=chapter&#038;p=3856"},"modified":"2025-04-07T14:39:35","modified_gmt":"2025-04-07T18:39:35","slug":"tutorial-policy-exploration-procedure","status":"publish","type":"chapter","link":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/tutorial-policy-exploration-procedure\/","title":{"raw":"Tutorial: Policy exploration procedure","rendered":"Tutorial: Policy exploration procedure"},"content":{"raw":"<div class=\"textbox shaded\">\r\n\r\n<strong>About this tutorial: <\/strong>It is intended for more complex applications of the fishing policy module rather than for introductory EwE courses.\u00a0There's a simpler tutorial to start with (<a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/tutorial-trade-offs-between-policy-objectives\/\">link<\/a>).\r\n\r\n<strong>Preparation<\/strong>: Read the Fishing policy optimization chapter (see <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/fishing-policy-exploration\/\">chapter<\/a>) and the User Guide interface description (see <a href=\"https:\/\/pressbooks.bccampus.ca\/eweguide\/chapter\/fishing-policy-exploration\/\">chapter<\/a>) before embarking on this tutorial.\r\n\r\n<\/div>\r\nThe EwE policy exploration module is a complex but capable beast, designed for policy <span style=\"text-decoration: underline\">exploration<\/span> of trade-offs, not for providing management advice for direct implementation. The policy advice it can produce is strategic rather than tactical (i.e. broad, directional policy advice rather than specific management advice). It can thus contribute at the table where policy discussions take place, in particular by providing options for and trade-offs in ecosystem-based management.\r\n\r\nIn this tutorial, we'll go through and explain details for how the module may be used for actual policy exploration. As part of this, we will outline, step by step, a procedure we find useful for conducting a more complete policy exploration that can be published and potentially can contribute to policy development.\r\n<h2>Model scope and behaviour<\/h2>\r\nWe assume that your model is indeed to be used for actual policy exploration, and advise that the model, while being predictive (rather than descriptive, see the <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/defining-the-ecosystem\/\">Defining the ecosystem<\/a> chapter) should include the fleets among which the policy module will seek to balance trade-offs \u2013 this may well be all fleets operating in the given ecosystem as trade-offs often are through food web interactions and bycatch. \u00a0Further, the target species for the fishery should be included in the model, along valued species that they take as bycatch, along with their prey and competitors, and where applicable top predators such as marine mammals and species of conservation concern.\r\n\r\nFor the model behaviour in response to proposed fisheries changes it is important that the model is fitted to time series data \u2013 this implies that density dependence related to <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/density-dependence-carrying-capacity-and-vulnerability-multipliers\/\">carrying capacity<\/a> has been considered (see <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/density-dependence-carrying-capacity-and-vulnerability-multipliers\/\">Density dependence and carrying capacity<\/a> chapter), and the vulnerability multipliers that affect population resilience have been modified accordingly, (see <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/vulnerability-and-vulnerability-multipliers\/\">Vulnerability and vulnerability multiplier<\/a> chapter).\r\n<div class=\"textbox shaded\">\r\n\r\nIt should hardly come as a surprise that we use Anchovy Bay for illustration in this tutorial. \u00a0If needed, you can download a version of the Anchovy Bay model that is fitted to time series data from <a href=\"https:\/\/ln5.sync.com\/dl\/552986540\/3k79yq3u-cmsf49gt-xkvbtqiw-5ikcwapm\">this link<\/a>. This is the version that is used for this tutorial.\r\n\r\nNote though that the time series fitting is rather incomplete in this version of the Anchovy Bay model, e.g., effort is lacking for several fleets.\r\n\r\n<\/div>\r\n<h3>Fleet considerations<\/h3>\r\nIf you have fleets in your system that it does not make sense to optimize for, e.g., optimizing for profit for a recreational fleet or a \"catch-all\" IUU[footnote]Illegal, unregulated and unreported[\/footnote] fleet, they can be considered in optimizations without being included in optimization searches. \u00a0For this, on the <em>Blocks<\/em> form in the interface select the first (black) block, then click the name of the fleet in the spreadsheet, and all years will be blackened out. When this is done, the Ecosim effort will remain as entered for that group, but the calculations of objectives will still include impacts of and on the blocked fleet(s).\r\n\r\nAs an example, the recreational fleet may be relying on a species that is also a target for a commercial fleet that is considered in the optimizations. Abundance changes caused by changing the commercial fleet effort will then impact the catch rates of the recreational fleet, which in turn will impact optimization measures that are affected by the recreational fleet.\r\n\r\nFor evaluating the <em>Social value (employment)<\/em> objective, it is important to set the relative <em>Jobs\/catch value<\/em> parameters in the <em>Ecosim &gt; Tools &gt; Fishing policy search<\/em> table. By default these are at unity (1) for all fleets, which is a reasonable assumption for Anchovy Bay (so we will leave all <em>Jobs\/catch value<\/em> at 1 for this tutorial), and indeed perhaps for many other places. Still, this parameters is important and it is best if you can use local or regional data.\r\n<h3>Time period: retrospective or forward?<\/h3>\r\nThere are two different approaches that have been used for policy optimizations, retrospective analysis vs. forward-looking scenarios, and we will describe both in the present extended tutorial.\r\n\r\nThe general advice is, where possible use forward-looking scenarios, but if for some reason this isn't feasible, the retrospective analysis is also a good option.\r\n<h4>Retrospective analysis<\/h4>\r\nFor this, the idea is to explore a what-if of the kind \"<em>how might the situation in this system have developed if the policy objectives had been implemented throughout the historic time period<\/em>\". This allows fitting of the model to time series data and to include environmental productivity patters. So, here we know what happened and we compare that to what might have happened if we had optimized the fisheries from early on.\r\n<h4>Forward-looking scenarios<\/h4>\r\nFor these, the idea is to use a model that has been fitted to time series data, running up to the present. We then explore what options we have looking forward to explore alternative policy optimization scenarios. \u00a0We may for instance add 20 years[footnote]Note that the policy optimization always runs for an additional 20 years without showing the results. This is to ensure that the optimizations do not result in an \"empty sea\" where fleets are encouraged to fish out the resources (which might happen if there was no tomorrow to consider).[\/footnote] to the runtime in Ecosim <em>(Ecosim &gt; Input &gt; Ecosim parameters &gt; Duration of simulation),\u00a0<\/em>then on the <em>Ecosim &gt; Tools &gt; Fishing policy, <\/em>click the black colour in the <em>Blocks<\/em> section, and block the historic time period for all fleets so that the optimizations will only be done for future years.\r\n\r\nForward-looking scenarios beg the question of how to handle a changing environment. The often-taken approaches for this are,\r\n<ol>\r\n \t<li><span style=\"text-align: initial;font-size: 1em\">Make the policy optimizations on the background of no change, i.e. keep the environmental productivity and patterns as they were at the end of the historic period,<\/span><\/li>\r\n \t<li><span style=\"text-align: initial;font-size: 1em\">If you have environmental productivity patterns, then repeat those pattens going forward,<\/span><\/li>\r\n \t<li><span style=\"text-align: initial;font-size: 1em\">Use output from climate models to drive the Ecosim environment.[footnote]This approach is used commonly for climate change scenarios, but may never have been used for policy optimizations.[\/footnote]\u00a0<\/span><\/li>\r\n<\/ol>\r\n<h2>Exploratory analysis<\/h2>\r\n<h3>Objective ranges<\/h3>\r\nPolicy explorations are often intended to explore less extreme, more balanced solutions for fleet tradeoffs. That calls for using weights on several policy objectives (see textbook <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/fishing-policy-exploration\/\">Policy exploration<\/a> chapter) \u2013 but what weights are needed to make the resulting fishing efforts \"balanced\"? Using the same weight (e.g., 1) for all objectives is not likely to results in a reasonable balance of performance measures. What then?\r\n\r\nWhen you are ready to explore the policy optimization, the first step is to evaluate the range of optimized fishing efforts that result for each objective. Open your model (or Anchovy Bay from <a href=\"https:\/\/ln5.sync.com\/dl\/0d48362b0\/ahwwesvy-2sz8peba-2x99usgn-ytcs4vrw\">this link<\/a>), load the Ecosim scenario you want to use, but do not load time series. Then run four policy searches with default settings varying only the objective weights. In the first search give the <em>Net economic value<\/em> a relative weight of 1, and leave the relative weight on all other objectives at 0. In the second run, give the <em>Social value (employment)<\/em> a value of 1, and all others 0.\u00a0 Then third and fourth runs are with only weights on <em>Ecosystem structure<\/em> and <em>Biodiversity<\/em>, respectively.\r\n\r\nThere is no need to include the <em>Mandated rebuilding<\/em> objective in the range of tested optimizations as this objective differs in behaviour, and serves a purpose different from other objectives. When invoked it is either, (1) a <u>forced rebuilding<\/u> or (2) providing a limit for the minimum acceptable biomass (<em>B<sub>lim<\/sub>\/B<sub>0<\/sub><\/em>)[footnote]<em>B<sub>lim<\/sub><\/em> is the lowest acceptable biomass and <em>B<sub>0<\/sub><\/em> the Ecopath biomass for the group. The objective is thus entered as a relative biomass term. [\/footnote]. This objective is thus intended to take precedence over all other objectives, and it should have a weight that trump all other objectives, (so very high, e.g., 100). The <em>Mandated rebuilding <\/em>is functional group specific, and only impacts the optimization when the biomass of a specified group falls below a user-defined reference level (<em>B<sub>lim<\/sub>\/B<sub>0<\/sub><\/em>). When the biomass is at or above <em>B<sub>lim<\/sub>\/B<sub>0<\/sub>, <\/em>the objective has no impact on the optimization.\r\n\r\nFor Anchovy Bay, we may get results as in Table 1 from the single objectives optimization runs.\r\n\r\n<strong>Table 1. Policy optimization objective ranges for four runs, each with weight on only one objective at the time. \u00a0<\/strong>\r\n\r\n[table id=14 \/]\r\n\r\nTable 1 shows the outcome from the four objective-by-objective runs. The range of objective values are indicated (ignoring negative <em>Net economic values)<\/em> and make it clear that the two economic objectives have the largest range, followed by <em>Ecosystem structure<\/em>. The\u00a0<em>Biodiversity\u00a0<\/em>objective has a much more narrow range. \u00a0There is no truth or absolute values coming out of this exploratory analysis, but it serves to illustrate that using equal weight on all objectives is unlikely to lead to a balanced solution. Instead, as a first estimate for weights to use across fleets one may be able to use the inverse of the ranges, see Table 1.\r\n<p style=\"font-weight: 400\">Using the inverse of the objective range <span style=\"text-decoration: underline\">may<\/span> lead to a more \u201cbalanced\u201d outcome during optimization exercises as this scaling better reflects the how \"easy\" it is for each objective to influence the objective function. This should, however, not be taken to mean that using the inverse range represents the truth, the whole truth and nothing but the truth: it is merely a more reasonable starting point than assuming the using the same relative weight for all objectives is a balanced solution.<\/p>\r\n&nbsp;\r\n<div>\r\n<div class=\"textbox shaded\">\r\n\r\n<strong>If a search is stuck with negative profit<\/strong>\r\n\r\nIt can happen when you start running optimizations with the \"limit cost &gt; earnings\" option, that both the net economic value and the social indicator (jobs) are negative already at the two initial runs. This indicates that a penalty has been applied in the search routine where the penalty increases rapidly as the ratio cost\/income increases toward and exceeds 1.0. Such a penalty is needed to make the optimization move away from fleet efforts that drive cost &gt; earnings.\r\n\r\nWhen it happens from the onset, it indicates that the baseline effort is unsustainable. In the case of Anchovy Bay, the culprit is the <em>sealers<\/em> fleet, which has a high, unsustainable effort in the base year. In Ecosim run, the fleet is shut down after a few years, but in the optimization, the high initial effort is maintained through the run. The optimization takes the cost and value at the baseline and sums the cost and value at the last year, and it multiplies that last year with a discounted value of what the last years catch would be worth if it were continued for an additional 20 years (i.e. there's a high weight on the end state relative to the baseline).\r\n\r\nIn some cases, the optimization routine can find its way out of the unsustainable fleet effort range, but not always. If the routine keeps producing negative indicators for the first two objectives, try making a run where you flatline the effort (1 throughout), and see which fleets end up having negative profit in the last year. Then reduce the effort for those fleets, and run the optimization again.\r\n\r\nSearches with \"limit cost &gt; earnings\" may have issues, and are not at all guaranteed to work. Consider if you can get by without using this option if there are problems with this in your model.\r\n\r\n<\/div>\r\nWhen running the policy optimization for Anchovy Bay with objective weights = the inverse of ranges from Table 1, the results in the table below are obtained.\r\n\r\n<strong>Table 2. Objective function values and fleet effort for a retrospective optimization for Anchovy Bay with weights set to inverse range of objectives. Optimization started with Ecopath base effort.<\/strong>\r\n\r\n[table id=16 \/]\r\n<h3>Introducing priors for weighting schedules<\/h3>\r\n<p style=\"font-weight: 400\">This type of weighting schedule perhaps best illustrates the outcomes of broad policy goals such as \u201c<em>securing the triple bottom line<\/em>\u201d or \u201c<em>promoting profitable, sustainable and just fisheries<\/em>\u201d.<\/p>\r\n<p style=\"font-weight: 400\">Granted that this discourse has become more mainstream, we are work under the assumption that fisheries stakeholders and the civil society are generally aiming for sustainable futures. However, different stakeholder groups will have different ways of weighting objectives, within and across fisheries systems. For example, fishers in some location may want to maximize fisheries rent twice as much as ecosystem structure, while recreational divers in another may want to maximize biodiversity 50% more than social benefits and 90% more than fisheries rent.<\/p>\r\n<p style=\"font-weight: 400\">Published literature using EwE\u2019s policy search routine has addressed this by building scenarios where different weighting schedules are applied to represent a \u201ccompromise\u201d among stakeholder objectives.[footnote]e.g., Natugonza et al. (2020) <a href=\"https:\/\/doi.org\/10.1016\/j.fishres.2020.105593\">https:\/\/doi.org\/10.1016\/j.fishres.2020.105593<\/a>, and Alms et al., (2022) <a href=\"https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349\">https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349<\/a>[\/footnote] However, grounding scenarios with field data is something that modellers should strive for. To adequately capture people\u2019s viewpoints, researchers should collect primary data (e.g., through semi-structured interviews or surveys) that allows them to construct alternative weighting schedules for different types of stakeholders operating in the system. This requires questions that rank optimization objectives and map the relative distance between them.<\/p>\r\n<p style=\"font-weight: 400\">For instance, let\u2019s say that a survey was implemented along the multiple communities that inhabit Anchovy Bay. The goal was to gather data on the policy preferences of coastal dwellers rather than those of the more \u201cusual suspects\u201d of fisheries systems (e.g., fishers, seafood processors, distributors, or sellers). People were asked to agree with certain statements (e.g., \u201c<em>Economic growth must be a priority for the my country, even if it affects the environment<\/em>\u201d or \u201c<em>Fishing effort should be restricted to prevent overfishing, even at the expense of losses in employment\"<\/em>), and the answers were arranged using a 5 point Likert scale with options ranging from \u201c<em>Strongly disagree<\/em>\u201d to \u201c<em>Strongly agree<\/em>\u201d. As questions addressed all permutations among objectives included in the policy search routine, we created a distance matrix for them and used it to rank them.<\/p>\r\n<p style=\"font-weight: 400\">Inhabitant of Anchovy Bay prioritized maximizing ecosystem structure as the central objective of local fisheries policies. This objective was closely followed by maximizing biodiversity, and more loosely followed by employment and revenue considerations. For the sake of simplicity, no questions addressed the mandated rebuilding objective. Inhabitants of Anchovy Bay are thus more pro-environment than the fishers, and certainly more in favour of conservation than local politicians believe their constituents to be. Using the responses we created an index where the objective with the highest value was given a value of 1, and all other objectives were scaled to it based on their relative importance, which leads to the weightings in Table 3. \u00a0These values can then be used as multipliers for the weights of the objectives used in the \u201cbalanced\u201d outcome run, to skew them towards their preferences and highlight the differences.<\/p>\r\n<strong>Table 3. Relative weighting as obtained from interviews with Anchovy Bay community members. The table also shows the inverse weighting from Table 1 and the resulting combined weight, which is the product of the community weighting and the inverse weighting.<\/strong>\r\n\r\n[table id=18 \/]\r\n\r\n<\/div>\r\n<strong>Table 4. Objective function values and fleet effort for a retrospective optimization for Anchovy Bay with weights set to the combined weigthings from Table 3. Optimization started with Ecopath base effort.<\/strong>\r\n\r\n[table id=19 \/]\r\n\r\nComparing Table 2 and Table 4, it's noticeable that adding the strong conversation-oriented view of the Anchovy Bay community has some socio-economic consequences. The economic rent and social benefits performance values were reduced with 20-25%, the ecosystem structure indicator increased ~30%, while the biodiversity measure (which is relative hard to impact) increased with a couple of percent.\r\n\r\nIf you want to explore the effect that the two optimizations have on fleet catches and values, and on group biomasses, you can extract that information (after each optimization) from the tables on the <em>Ecosim &gt; Output &gt; Ecosim results<\/em>\u00a0form.\r\n<div>\r\n<h3>Local maxima<\/h3>\r\nAs part of the exploratory analysis, it is important to check whether the maximization search is impacted by the start point, i.e. whether the optimization solutions are unique. By default the optimization routine will start with the fishing rates defined by the Ecopath baseline. It's possible, however, to instead using random fleet effort <em>(Random F's<\/em> in the policy interface) to check if the optimization routine is likely to get stuck at local maxima. \u00a0All optimization routines are impacted by this, the ones in EwE being no exceptions.\r\n\r\nWe illustrate this for Anchovy Bay by running 70 optimizations with the combined rel. weights from Table 3, and runs initialized with random F's. The outcome of that exploratory analysis is presented in Figure 1.\r\n\r\n<img class=\"alignnone size-full wp-image-4120\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay.png\" alt=\"\" width=\"1000\" height=\"400\" \/>\r\n\r\n<strong>Figure 1. Box plot[footnote]<\/strong>The R-code and CSV-file used for this figure can be downloaded from this <a href=\"https:\/\/ln5.sync.com\/dl\/3065219c0\/gb3t8j6n-pdmu9gh7-yqqgnusb-jh89zs6a\">link<\/a>.<strong>[\/footnote] showing minimum, first quartile, median, third quartile and max value of objective function value for indicators and relative effort for fleets for 70 policy optimizations run with random <em>Starting F<\/em> initialization. \u00a0It's clear that only the relative effort for the <em>bait boats<\/em> varies between runs (with one run way higher than the 69 others), and that the optimizations overall are robust to local maxima.<\/strong>\r\n\r\nFrom Figure 1, it is clear that for this version[footnote]Optimizations depend on fitting and model parameters. Different versions of the Anchovy Bay may well lead to different optimizations[\/footnote] of the Anchovy Bay model, there is very little tendency for the optimization to get stuck on local maxima. The variation in the objective estimates and effort patterns are very similar across all runs (apart perhaps from <em>bait boat<\/em> effort), with only a few runs indicating presence of local maxima.\r\n<div class=\"textbox shaded\">\r\n\r\nThat this Anchovy Bay model isn't prone to get stuck on local maxima should not invite complacency, but be seen as a warning that there <span style=\"text-decoration: underline\">may<\/span> be local maxima. We ran the optimizations 70 times for Figure 1 to illustrate this point \u2013 but also to show that it rarely happens. Bottom line is, run your optimization a number of times (maybe 10 to 20), check how consistent the output is. Then run the model once, check if the outcome corresponds to the majority of the random runs. If it does, you're good to go.\r\n\r\n<\/div>\r\nThe overall conclusion is that policy optimizations for Anchovy Bay are not very prone to get stuck on local maxima. This is also what we have found for many other ecosystem model optimizations, giving some comfort that the starting point isn't very critical. Still, this needs to be checked for all models, so including a search with random <em>Starting F's<\/em> should be included in all more serious policy explorations.\r\n<h2>Fleet trade-off analysis<\/h2>\r\n<img class=\"size-medium wp-image-4013 alignleft\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Anchovy-Bay-FleetTradeoff-e1733339668222-300x296.jpg\" alt=\"\" width=\"300\" height=\"296\" \/>\r\n\r\n<strong>Figure 2. Fleet trade-off analysis for Anchovy Bay showing impact a 10% reduction in effort for the fleet listed in rows has on the fleets listed above columns. Negative impacts are in red, and positive in blue.\u00a0<\/strong>\r\n\r\nA next step of exploratory analysis is the fleet trade-off analysis described in the Fishing policy chapter (<a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/fishing-policy-exploration#fleet tradeoff\">link to fleet trade-of<\/a>f). We refer to that section for description, including for code to produce plots.\r\n\r\nWe suggest that you perform the fleet trade-off analysis for your model and explore the outcome. For Anchovy Bay (Figure 2), the plot shows the impact that a 10% reduction for the fleets mentioned to the left, are predicted to have in the value of the catch for all fleets. <span style=\"color: #ff0000\">Red<\/span> circles indicate reduction in value, and <span style=\"color: #0000ff\">blue<\/span> increase. \u00a0The impacts are displayed so that the circle areas are proportional to the changes in value of the catch, and are thus comparable across fleets[footnote]The monetary value of the fleet-tradeoffs can be obtained from the CSV file used for producing the fleet trade-off plots[\/footnote].\r\n\r\nFor Anchovy Bay, the fleet trade-off analysis shows some both straightforward and more complex patterns. Notice for instance that a reduction in <em>bait boats<\/em>' effort will have a positive impact on <em>seiners<\/em>. That makes sense since the two fleets both catch anchovy. But conversely, a reduction in <em>seiners<\/em>' effort lead to a small reduction in landed value for the <em>bait boats<\/em>. Why? <em>Seiners<\/em> also catch mackerel, and the reduced effort will lead to more mackerel, which in turn will have a negative impact on anchovy, and hence on <em>bait boats<\/em>, (which only catch anchovy).\r\n\r\nThe strongest impact of effort reduction is for <em>trawlers<\/em> and <em>shrimpers<\/em>, the two fleets that have the highest value in the model base year. Reduction in <em>trawlers<\/em>' effort has a considerable negative impact on <em>trawlers<\/em>, but also an almost corresponding negative impact on <em>shrimpers<\/em>. That makes sense as the reduction should lead to more cod and whiting, both of which eat shrimp. \u00a0But conversely, reducing <em>shrimpers<\/em>' effort leads to an <span style=\"text-decoration: underline\">increase<\/span> in shrimp landings (so they must be overexploited in the baseline \u2013 given the baseline predator-prey conditions though). But why does this not lead to an increase in the value of <em>trawlers<\/em>' landings? \u00a0The reason is that shrimp have a positive impact on whiting and mackerel, but a negative impact on cod. Notably, the increase in whiting impacts cod negatively. \u00a0Those relationships becomes clearer if you check out the Mixed Trophic Impact analysis <em>(Ecopath &gt; Output &gt; Tools &gt; Network Analysis &gt; Mixed trophic impact &gt; Impact<\/em> data), which shows that <em>shrimpers<\/em> have opposite impact on cod <em>versus<\/em> whiting and mackerel.\r\n<h2>Your policy questions?<\/h2>\r\nHaving explored the behaviour of your model, e.g., as described above, the next issue is to clearly define what questions you are asking for your model. An example of this is Alms et al. 2022, [footnote]Alms V, Romagnoni G, Wolff M. Exploration of fisheries management policies in the Gulf of Nicoya (Costa Rica) using ecosystem modelling. Ocean and Coastal Management 230 (2022) 106349.\u00a0<a href=\"https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349\">https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349<\/a>[\/footnote] who compared output from three defined Ecosim scenarios, (1) ban on shrimp trawling, (2) gill net effort reduction of 25%, and (3) the combination of (1) and (2), and then compared these to the output of single- and multi-objective policy optimizations. Here, the multi-objective optimizations were designed to serve as more balanced solutions. \u00a0See the paper for details.\r\n<h3>Anchovy Bay<\/h3>\r\nWhile the complex patterns in the fleet tradeoff (Figure 2) can be explained as above, they raise some questions. There is a big negative impact of reducing <em>trawlers<\/em>' effort and a big positive impact of reducing <em>shrimpers<\/em>' effort. \u00a0Is this then what we should explore? Well, it's certainly interesting scenarios, but there are some complications, and it is in line with what actually happened in Anchovy Bay. Effort of <em>trawlers<\/em> indeed increased 2-3 times over time, but the shrimp effort increased by almost an order of magnitude. \u00a0Does this make sense when the fleet trade-off analysis indicate that shrimp catches would increase if <em>shrimpers<\/em>' effort was reduced? \u00a0The finding is not really wrong, but it doesn't not consider that the increase in the <em>trawlers<\/em>' effort reduced the abundance of cod and whiting, which in turn lead to many more shrimps being available for the <em>shrimpers<\/em>. Predator-release!\r\n\r\nWhat the above indicates is in essence that the Anchovy Bay ecosystem of today is very different from that of 1970. Therefore, we will next run the policy optimization for Anchovy Bay as a forward-looking scenario where we keep the run from 1970-2010 as it was in the fitting, and then conduct the policy optimization from 2011 (year 42) onwards. Figure 3 illustrates the setup.\r\n\r\nIf you change the \u00a0<em>Ecosim &gt; Tools &gt; Fishing policy search &gt; Base year<\/em> to the end of the time series, the routine will automatically black out the years of the time series. In that case, the economic data (cost and value) in <em>Ecopath &gt; Input &gt; Fishery &gt; Fleets<\/em> are assumed to represent your base year, not the Ecopath model year (Year 1 or 1970 for Anchovy Bay).\r\n\r\n<\/div>\r\n<img class=\"alignnone size-full wp-image-4123\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18.png\" alt=\"\" width=\"3416\" height=\"388\" \/>\r\n\r\n<strong>Figure 3. Fishing policy search interface set up to do a forward scenario for Anchovy Bay. The optimization will not change the effort for the first 41 years (in black), but search for one effort for each fleet for the last 20 years. If you cannot see all the years, try the zoom icons in the interface and widen the aspect ratio of the interface. \u00a0\u00a0<\/strong>\r\n\r\nWe run the forward-looking scenarios with <em>Random F<\/em> drawn to start the optimizations, which will ensure that we do not get stuck in a local maxima for all runs (even though we have found above that the risk of this is small). We here do not use the <em>\"limit cost &gt; earnings\" <\/em>as this tends to hang the optimizations for the Anchovy Bay model[footnote]This is likely because we for this tutorial are using very rough economic data. If you have this issue, do check your assumptions about cost vs value in your base year.[\/footnote].\r\n\r\n<img class=\"size-full wp-image-4139 aligncenter\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Objective-plot.png\" alt=\"\" width=\"600\" height=\"400\" \/>\r\n\r\n<strong>Figure 4. Indicator values from Anchovy Bay model policy optimizations with <em>inverse range<\/em> weighting <em>vs<\/em>. <em>community-based<\/em> weightings (from Table 3). The community places higher weights on <em>ecosystem structure<\/em> and <em>biodiversity<\/em>, and the optimizations indeed reflect this. It's clear that all of the objectives contributes to the total objective function score, which is the sum of the four objectives \u2013 actually with less variation for the <em>community-based <\/em>weighting. But the improved ecological score is offset by lower employment and profit[footnote]<\/strong>R-code and data files to produce this code are included in the Zip file that can be downloaded for Figure 5, see below.<strong>[\/footnote]. \u00a0\u00a0<\/strong>\r\n\r\n&nbsp;\r\n\r\n<img class=\"alignnone size-full wp-image-4130\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-results-X.png\" alt=\"\" width=\"1200\" height=\"800\" \/>\r\n\r\n<strong>Figure 5. Plot[footnote]<\/strong>Download R-code and CSV-files used for Figure 4 and 5 from this <a href=\"https:\/\/ln5.sync.com\/dl\/6be981c90\/djen2t7n-keic5kps-6u9scv7p-b48hnx2r\">link<\/a>.<strong>[\/footnote] comparing fishing policy optimizations for Anchovy Bay model with weight as defined based on inverse range and community input in Table 3. Value is proportional to jobs, so the <em>Jobs end<\/em> plot indicates around a 25% reduction in jobs with the community-based weighting. \u00a0<\/strong>\r\n\r\nThe fleet effort value, catch and biomass results for the Anchovy Bay model forward-looking policy scenario are presented in Figure 5, again comparing optimizations with the inverse range weightings vs the community opinion based weightings (both from Table 3). The results show the complex trade-off patterns between the two weighting schemes.\r\n<div>\r\n<h3>Your case<\/h3>\r\nWe suggest that you as a starting point go through this extended tutorial and conduct policy optimizations on a fitted version of your ecosystem model, working your way through the steps described in the tutorial.\r\n<div class=\"textbox shaded\">\r\n\r\n<strong>Pretty good yield<\/strong>\r\n\r\nWhen you run optimizations, you'll often see that the objective function may rather quickly get to more than 98% of the final objective function score. If you start with <em>Ecopath base F<\/em>, the effort changes may be rather limited when you pass the 98% mark, but the last few percent may cause big changes in effort. \u00a0It may be worth exploring those intermediate states (by stopping the optimizations before it finishes by itself).\r\n\r\nThis is to some degree related to the idea behind Ray Hilborn's <em>Pretty Good Yield[footnote]<a href=\"https:\/\/doi.org\/10.1016\/j.marpol.2009.04.013\">https:\/\/doi.org\/10.1016\/j.marpol.2009.04.013 <\/a>[\/footnote]<\/em> for MSY.\r\n\r\n<\/div>\r\n&nbsp;\r\n<div class=\"textbox textbox--examples\"><header class=\"textbox__header\">\r\n<p class=\"textbox__title\"><strong>Acknowledgement<\/strong><\/p>\r\n\r\n<\/header>\r\n<div class=\"textbox__content\">\r\n<div class=\"textbox__content\"><img class=\"alignright wp-image-3830 size-medium\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/EcoScope-logo-300x113.png\" alt=\"\" width=\"300\" height=\"113\" \/>This chapter was developed for the <a href=\"https:\/\/ecoscopium.eu\">EcoScope<\/a> project to guide implementation of the EwE Policy Search for the project case studies. EcoScope is funded from the <a href=\"https:\/\/ec.europa.eu\/programmes\/horizon2020\/en\/home\" target=\"_blank\" rel=\"noopener noreferrer\">European Commission\u2019s Horizon 2020 Research and Innovation programme<\/a> under grant agreement No 101000302. Project coordinator: Aristotle University of Thessaloniki, Greece. \u00a0Parts of the text are from the unpublished EwE User Guide: Christensen V, C Walters, D Pauly, R Forrest. Ecopath with Ecosim. User Guide. November 2008.<\/div>\r\n<\/div>\r\n<\/div>\r\n&nbsp;\r\n\r\n<\/div>","rendered":"<div class=\"textbox shaded\">\n<p><strong>About this tutorial: <\/strong>It is intended for more complex applications of the fishing policy module rather than for introductory EwE courses.\u00a0There&#8217;s a simpler tutorial to start with (<a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/tutorial-trade-offs-between-policy-objectives\/\">link<\/a>).<\/p>\n<p><strong>Preparation<\/strong>: Read the Fishing policy optimization chapter (see <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/fishing-policy-exploration\/\">chapter<\/a>) and the User Guide interface description (see <a href=\"https:\/\/pressbooks.bccampus.ca\/eweguide\/chapter\/fishing-policy-exploration\/\">chapter<\/a>) before embarking on this tutorial.<\/p>\n<\/div>\n<p>The EwE policy exploration module is a complex but capable beast, designed for policy <span style=\"text-decoration: underline\">exploration<\/span> of trade-offs, not for providing management advice for direct implementation. The policy advice it can produce is strategic rather than tactical (i.e. broad, directional policy advice rather than specific management advice). It can thus contribute at the table where policy discussions take place, in particular by providing options for and trade-offs in ecosystem-based management.<\/p>\n<p>In this tutorial, we&#8217;ll go through and explain details for how the module may be used for actual policy exploration. As part of this, we will outline, step by step, a procedure we find useful for conducting a more complete policy exploration that can be published and potentially can contribute to policy development.<\/p>\n<h2>Model scope and behaviour<\/h2>\n<p>We assume that your model is indeed to be used for actual policy exploration, and advise that the model, while being predictive (rather than descriptive, see the <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/defining-the-ecosystem\/\">Defining the ecosystem<\/a> chapter) should include the fleets among which the policy module will seek to balance trade-offs \u2013 this may well be all fleets operating in the given ecosystem as trade-offs often are through food web interactions and bycatch. \u00a0Further, the target species for the fishery should be included in the model, along valued species that they take as bycatch, along with their prey and competitors, and where applicable top predators such as marine mammals and species of conservation concern.<\/p>\n<p>For the model behaviour in response to proposed fisheries changes it is important that the model is fitted to time series data \u2013 this implies that density dependence related to <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/density-dependence-carrying-capacity-and-vulnerability-multipliers\/\">carrying capacity<\/a> has been considered (see <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/density-dependence-carrying-capacity-and-vulnerability-multipliers\/\">Density dependence and carrying capacity<\/a> chapter), and the vulnerability multipliers that affect population resilience have been modified accordingly, (see <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/vulnerability-and-vulnerability-multipliers\/\">Vulnerability and vulnerability multiplier<\/a> chapter).<\/p>\n<div class=\"textbox shaded\">\n<p>It should hardly come as a surprise that we use Anchovy Bay for illustration in this tutorial. \u00a0If needed, you can download a version of the Anchovy Bay model that is fitted to time series data from <a href=\"https:\/\/ln5.sync.com\/dl\/552986540\/3k79yq3u-cmsf49gt-xkvbtqiw-5ikcwapm\">this link<\/a>. This is the version that is used for this tutorial.<\/p>\n<p>Note though that the time series fitting is rather incomplete in this version of the Anchovy Bay model, e.g., effort is lacking for several fleets.<\/p>\n<\/div>\n<h3>Fleet considerations<\/h3>\n<p>If you have fleets in your system that it does not make sense to optimize for, e.g., optimizing for profit for a recreational fleet or a &#8220;catch-all&#8221; IUU<a class=\"footnote\" title=\"Illegal, unregulated and unreported\" id=\"return-footnote-3856-1\" href=\"#footnote-3856-1\" aria-label=\"Footnote 1\"><sup class=\"footnote\">[1]<\/sup><\/a> fleet, they can be considered in optimizations without being included in optimization searches. \u00a0For this, on the <em>Blocks<\/em> form in the interface select the first (black) block, then click the name of the fleet in the spreadsheet, and all years will be blackened out. When this is done, the Ecosim effort will remain as entered for that group, but the calculations of objectives will still include impacts of and on the blocked fleet(s).<\/p>\n<p>As an example, the recreational fleet may be relying on a species that is also a target for a commercial fleet that is considered in the optimizations. Abundance changes caused by changing the commercial fleet effort will then impact the catch rates of the recreational fleet, which in turn will impact optimization measures that are affected by the recreational fleet.<\/p>\n<p>For evaluating the <em>Social value (employment)<\/em> objective, it is important to set the relative <em>Jobs\/catch value<\/em> parameters in the <em>Ecosim &gt; Tools &gt; Fishing policy search<\/em> table. By default these are at unity (1) for all fleets, which is a reasonable assumption for Anchovy Bay (so we will leave all <em>Jobs\/catch value<\/em> at 1 for this tutorial), and indeed perhaps for many other places. Still, this parameters is important and it is best if you can use local or regional data.<\/p>\n<h3>Time period: retrospective or forward?<\/h3>\n<p>There are two different approaches that have been used for policy optimizations, retrospective analysis vs. forward-looking scenarios, and we will describe both in the present extended tutorial.<\/p>\n<p>The general advice is, where possible use forward-looking scenarios, but if for some reason this isn&#8217;t feasible, the retrospective analysis is also a good option.<\/p>\n<h4>Retrospective analysis<\/h4>\n<p>For this, the idea is to explore a what-if of the kind &#8220;<em>how might the situation in this system have developed if the policy objectives had been implemented throughout the historic time period<\/em>&#8220;. This allows fitting of the model to time series data and to include environmental productivity patters. So, here we know what happened and we compare that to what might have happened if we had optimized the fisheries from early on.<\/p>\n<h4>Forward-looking scenarios<\/h4>\n<p>For these, the idea is to use a model that has been fitted to time series data, running up to the present. We then explore what options we have looking forward to explore alternative policy optimization scenarios. \u00a0We may for instance add 20 years<a class=\"footnote\" title=\"Note that the policy optimization always runs for an additional 20 years without showing the results. This is to ensure that the optimizations do not result in an &quot;empty sea&quot; where fleets are encouraged to fish out the resources (which might happen if there was no tomorrow to consider).\" id=\"return-footnote-3856-2\" href=\"#footnote-3856-2\" aria-label=\"Footnote 2\"><sup class=\"footnote\">[2]<\/sup><\/a> to the runtime in Ecosim <em>(Ecosim &gt; Input &gt; Ecosim parameters &gt; Duration of simulation),\u00a0<\/em>then on the <em>Ecosim &gt; Tools &gt; Fishing policy, <\/em>click the black colour in the <em>Blocks<\/em> section, and block the historic time period for all fleets so that the optimizations will only be done for future years.<\/p>\n<p>Forward-looking scenarios beg the question of how to handle a changing environment. The often-taken approaches for this are,<\/p>\n<ol>\n<li><span style=\"text-align: initial;font-size: 1em\">Make the policy optimizations on the background of no change, i.e. keep the environmental productivity and patterns as they were at the end of the historic period,<\/span><\/li>\n<li><span style=\"text-align: initial;font-size: 1em\">If you have environmental productivity patterns, then repeat those pattens going forward,<\/span><\/li>\n<li><span style=\"text-align: initial;font-size: 1em\">Use output from climate models to drive the Ecosim environment.<a class=\"footnote\" title=\"This approach is used commonly for climate change scenarios, but may never have been used for policy optimizations.\" id=\"return-footnote-3856-3\" href=\"#footnote-3856-3\" aria-label=\"Footnote 3\"><sup class=\"footnote\">[3]<\/sup><\/a>\u00a0<\/span><\/li>\n<\/ol>\n<h2>Exploratory analysis<\/h2>\n<h3>Objective ranges<\/h3>\n<p>Policy explorations are often intended to explore less extreme, more balanced solutions for fleet tradeoffs. That calls for using weights on several policy objectives (see textbook <a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/fishing-policy-exploration\/\">Policy exploration<\/a> chapter) \u2013 but what weights are needed to make the resulting fishing efforts &#8220;balanced&#8221;? Using the same weight (e.g., 1) for all objectives is not likely to results in a reasonable balance of performance measures. What then?<\/p>\n<p>When you are ready to explore the policy optimization, the first step is to evaluate the range of optimized fishing efforts that result for each objective. Open your model (or Anchovy Bay from <a href=\"https:\/\/ln5.sync.com\/dl\/0d48362b0\/ahwwesvy-2sz8peba-2x99usgn-ytcs4vrw\">this link<\/a>), load the Ecosim scenario you want to use, but do not load time series. Then run four policy searches with default settings varying only the objective weights. In the first search give the <em>Net economic value<\/em> a relative weight of 1, and leave the relative weight on all other objectives at 0. In the second run, give the <em>Social value (employment)<\/em> a value of 1, and all others 0.\u00a0 Then third and fourth runs are with only weights on <em>Ecosystem structure<\/em> and <em>Biodiversity<\/em>, respectively.<\/p>\n<p>There is no need to include the <em>Mandated rebuilding<\/em> objective in the range of tested optimizations as this objective differs in behaviour, and serves a purpose different from other objectives. When invoked it is either, (1) a <u>forced rebuilding<\/u> or (2) providing a limit for the minimum acceptable biomass (<em>B<sub>lim<\/sub>\/B<sub>0<\/sub><\/em>)<a class=\"footnote\" title=\"Blim is the lowest acceptable biomass and B0 the Ecopath biomass for the group. The objective is thus entered as a relative biomass term.\" id=\"return-footnote-3856-4\" href=\"#footnote-3856-4\" aria-label=\"Footnote 4\"><sup class=\"footnote\">[4]<\/sup><\/a>. This objective is thus intended to take precedence over all other objectives, and it should have a weight that trump all other objectives, (so very high, e.g., 100). The <em>Mandated rebuilding <\/em>is functional group specific, and only impacts the optimization when the biomass of a specified group falls below a user-defined reference level (<em>B<sub>lim<\/sub>\/B<sub>0<\/sub><\/em>). When the biomass is at or above <em>B<sub>lim<\/sub>\/B<sub>0<\/sub>, <\/em>the objective has no impact on the optimization.<\/p>\n<p>For Anchovy Bay, we may get results as in Table 1 from the single objectives optimization runs.<\/p>\n<p><strong>Table 1. Policy optimization objective ranges for four runs, each with weight on only one objective at the time. \u00a0<\/strong><\/p>\n<table id=\"tablepress-14\" class=\"tablepress tablepress-id-14\">\n<thead>\n<tr class=\"row-1\">\n<th class=\"column-1\"><strong>Optimizing for \\ Objective<\/strong><\/th>\n<th class=\"column-2\"><strong>Econ.<\/strong><\/th>\n<th class=\"column-3\"><strong>Social<\/strong><\/th>\n<th class=\"column-4\"><strong>Ecosys.<\/strong><\/th>\n<th class=\"column-5\"><strong>Biodiversity<\/strong><\/th>\n<\/tr>\n<\/thead>\n<tbody class=\"row-striping row-hover\">\n<tr class=\"row-2\">\n<td class=\"column-1\">Econ.<\/td>\n<td class=\"column-2\"><strong>  1.84 <\/strong><\/td>\n<td class=\"column-3\">  1.15 <\/td>\n<td class=\"column-4\">  0.96 <\/td>\n<td class=\"column-5\">  1.02 <\/td>\n<\/tr>\n<tr class=\"row-3\">\n<td class=\"column-1\">Social<\/td>\n<td class=\"column-2\">  (2.26)<\/td>\n<td class=\"column-3\"><strong>  2.34 <\/strong><\/td>\n<td class=\"column-4\">  0.59 <\/td>\n<td class=\"column-5\">  0.88 <\/td>\n<\/tr>\n<tr class=\"row-4\">\n<td class=\"column-1\">Ecosys.<\/td>\n<td class=\"column-2\">  1.28 <\/td>\n<td class=\"column-3\">  0.68 <\/td>\n<td class=\"column-4\"><strong>  1.32 <\/strong><\/td>\n<td class=\"column-5\">  1.05 <\/td>\n<\/tr>\n<tr class=\"row-5\">\n<td class=\"column-1\">Biodiversity<\/td>\n<td class=\"column-2\">  0.71 <\/td>\n<td class=\"column-3\">  0.50 <\/td>\n<td class=\"column-4\">  1.08 <\/td>\n<td class=\"column-5\"><strong>  1.05 <\/strong><\/td>\n<\/tr>\n<tr class=\"row-6\">\n<td class=\"column-1\">Min. value<\/td>\n<td class=\"column-2\">  -   <\/td>\n<td class=\"column-3\">  0.50 <\/td>\n<td class=\"column-4\">  0.59 <\/td>\n<td class=\"column-5\">  0.88 <\/td>\n<\/tr>\n<tr class=\"row-7\">\n<td class=\"column-1\">Max. value<\/td>\n<td class=\"column-2\">  1.84 <\/td>\n<td class=\"column-3\">  2.34 <\/td>\n<td class=\"column-4\">  1.32 <\/td>\n<td class=\"column-5\">  1.05 <\/td>\n<\/tr>\n<tr class=\"row-8\">\n<td class=\"column-1\">Range<\/td>\n<td class=\"column-2\">  1.84 <\/td>\n<td class=\"column-3\">  1.83 <\/td>\n<td class=\"column-4\">  0.73 <\/td>\n<td class=\"column-5\">  0.17 <\/td>\n<\/tr>\n<tr class=\"row-9\">\n<td class=\"column-1\">1 \/ range<\/td>\n<td class=\"column-2\">  0.5 <\/td>\n<td class=\"column-3\">  0.5 <\/td>\n<td class=\"column-4\">  1.4 <\/td>\n<td class=\"column-5\">  5.9 <\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p><!-- #tablepress-14 from cache --><\/p>\n<p>Table 1 shows the outcome from the four objective-by-objective runs. The range of objective values are indicated (ignoring negative <em>Net economic values)<\/em> and make it clear that the two economic objectives have the largest range, followed by <em>Ecosystem structure<\/em>. The\u00a0<em>Biodiversity\u00a0<\/em>objective has a much more narrow range. \u00a0There is no truth or absolute values coming out of this exploratory analysis, but it serves to illustrate that using equal weight on all objectives is unlikely to lead to a balanced solution. Instead, as a first estimate for weights to use across fleets one may be able to use the inverse of the ranges, see Table 1.<\/p>\n<p style=\"font-weight: 400\">Using the inverse of the objective range <span style=\"text-decoration: underline\">may<\/span> lead to a more \u201cbalanced\u201d outcome during optimization exercises as this scaling better reflects the how &#8220;easy&#8221; it is for each objective to influence the objective function. This should, however, not be taken to mean that using the inverse range represents the truth, the whole truth and nothing but the truth: it is merely a more reasonable starting point than assuming the using the same relative weight for all objectives is a balanced solution.<\/p>\n<p>&nbsp;<\/p>\n<div>\n<div class=\"textbox shaded\">\n<p><strong>If a search is stuck with negative profit<\/strong><\/p>\n<p>It can happen when you start running optimizations with the &#8220;limit cost &gt; earnings&#8221; option, that both the net economic value and the social indicator (jobs) are negative already at the two initial runs. This indicates that a penalty has been applied in the search routine where the penalty increases rapidly as the ratio cost\/income increases toward and exceeds 1.0. Such a penalty is needed to make the optimization move away from fleet efforts that drive cost &gt; earnings.<\/p>\n<p>When it happens from the onset, it indicates that the baseline effort is unsustainable. In the case of Anchovy Bay, the culprit is the <em>sealers<\/em> fleet, which has a high, unsustainable effort in the base year. In Ecosim run, the fleet is shut down after a few years, but in the optimization, the high initial effort is maintained through the run. The optimization takes the cost and value at the baseline and sums the cost and value at the last year, and it multiplies that last year with a discounted value of what the last years catch would be worth if it were continued for an additional 20 years (i.e. there&#8217;s a high weight on the end state relative to the baseline).<\/p>\n<p>In some cases, the optimization routine can find its way out of the unsustainable fleet effort range, but not always. If the routine keeps producing negative indicators for the first two objectives, try making a run where you flatline the effort (1 throughout), and see which fleets end up having negative profit in the last year. Then reduce the effort for those fleets, and run the optimization again.<\/p>\n<p>Searches with &#8220;limit cost &gt; earnings&#8221; may have issues, and are not at all guaranteed to work. Consider if you can get by without using this option if there are problems with this in your model.<\/p>\n<\/div>\n<p>When running the policy optimization for Anchovy Bay with objective weights = the inverse of ranges from Table 1, the results in the table below are obtained.<\/p>\n<p><strong>Table 2. Objective function values and fleet effort for a retrospective optimization for Anchovy Bay with weights set to inverse range of objectives. Optimization started with Ecopath base effort.<\/strong><\/p>\n<table id=\"tablepress-16\" class=\"tablepress tablepress-id-16\">\n<thead>\n<tr class=\"row-1\">\n<th class=\"column-1\"><strong>Total<\/strong><\/th>\n<th class=\"column-2\"><strong>Econ.<\/strong><\/th>\n<th class=\"column-3\"><strong>Social<\/strong><\/th>\n<th class=\"column-4\"><strong>Ecosys.<\/strong><\/th>\n<th class=\"column-5\"><strong>Biodiversity<\/strong><\/th>\n<\/tr>\n<\/thead>\n<tbody class=\"row-striping row-hover\">\n<tr class=\"row-2\">\n<td class=\"column-1\">  3.70 <\/td>\n<td class=\"column-2\">  1.69 <\/td>\n<td class=\"column-3\">  1.08 <\/td>\n<td class=\"column-4\">  1.05 <\/td>\n<td class=\"column-5\">  1.02 <\/td>\n<\/tr>\n<tr class=\"row-3\">\n<td class=\"column-1\"><\/td>\n<td class=\"column-2\"><\/td>\n<td class=\"column-3\"><\/td>\n<td class=\"column-4\"><\/td>\n<td class=\"column-5\"><\/td>\n<\/tr>\n<tr class=\"row-4\">\n<td class=\"column-1\"><strong>Sealers<\/strong><\/td>\n<td class=\"column-2\"><strong>Trawlers<\/strong><\/td>\n<td class=\"column-3\"><strong>Seiners<\/strong><\/td>\n<td class=\"column-4\"><strong>Bait boats<\/strong><\/td>\n<td class=\"column-5\"><strong>Shrimpers<\/strong><\/td>\n<\/tr>\n<tr class=\"row-5\">\n<td class=\"column-1\">  0.89 <\/td>\n<td class=\"column-2\">  1.66 <\/td>\n<td class=\"column-3\">  0.72 <\/td>\n<td class=\"column-4\">  0.85 <\/td>\n<td class=\"column-5\">  0.59 <\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p><!-- #tablepress-16 from cache --><\/p>\n<h3>Introducing priors for weighting schedules<\/h3>\n<p style=\"font-weight: 400\">This type of weighting schedule perhaps best illustrates the outcomes of broad policy goals such as \u201c<em>securing the triple bottom line<\/em>\u201d or \u201c<em>promoting profitable, sustainable and just fisheries<\/em>\u201d.<\/p>\n<p style=\"font-weight: 400\">Granted that this discourse has become more mainstream, we are work under the assumption that fisheries stakeholders and the civil society are generally aiming for sustainable futures. However, different stakeholder groups will have different ways of weighting objectives, within and across fisheries systems. For example, fishers in some location may want to maximize fisheries rent twice as much as ecosystem structure, while recreational divers in another may want to maximize biodiversity 50% more than social benefits and 90% more than fisheries rent.<\/p>\n<p style=\"font-weight: 400\">Published literature using EwE\u2019s policy search routine has addressed this by building scenarios where different weighting schedules are applied to represent a \u201ccompromise\u201d among stakeholder objectives.<a class=\"footnote\" title=\"e.g., Natugonza et al. (2020) https:\/\/doi.org\/10.1016\/j.fishres.2020.105593, and Alms et al., (2022) https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349\" id=\"return-footnote-3856-5\" href=\"#footnote-3856-5\" aria-label=\"Footnote 5\"><sup class=\"footnote\">[5]<\/sup><\/a> However, grounding scenarios with field data is something that modellers should strive for. To adequately capture people\u2019s viewpoints, researchers should collect primary data (e.g., through semi-structured interviews or surveys) that allows them to construct alternative weighting schedules for different types of stakeholders operating in the system. This requires questions that rank optimization objectives and map the relative distance between them.<\/p>\n<p style=\"font-weight: 400\">For instance, let\u2019s say that a survey was implemented along the multiple communities that inhabit Anchovy Bay. The goal was to gather data on the policy preferences of coastal dwellers rather than those of the more \u201cusual suspects\u201d of fisheries systems (e.g., fishers, seafood processors, distributors, or sellers). People were asked to agree with certain statements (e.g., \u201c<em>Economic growth must be a priority for the my country, even if it affects the environment<\/em>\u201d or \u201c<em>Fishing effort should be restricted to prevent overfishing, even at the expense of losses in employment&#8221;<\/em>), and the answers were arranged using a 5 point Likert scale with options ranging from \u201c<em>Strongly disagree<\/em>\u201d to \u201c<em>Strongly agree<\/em>\u201d. As questions addressed all permutations among objectives included in the policy search routine, we created a distance matrix for them and used it to rank them.<\/p>\n<p style=\"font-weight: 400\">Inhabitant of Anchovy Bay prioritized maximizing ecosystem structure as the central objective of local fisheries policies. This objective was closely followed by maximizing biodiversity, and more loosely followed by employment and revenue considerations. For the sake of simplicity, no questions addressed the mandated rebuilding objective. Inhabitants of Anchovy Bay are thus more pro-environment than the fishers, and certainly more in favour of conservation than local politicians believe their constituents to be. Using the responses we created an index where the objective with the highest value was given a value of 1, and all other objectives were scaled to it based on their relative importance, which leads to the weightings in Table 3. \u00a0These values can then be used as multipliers for the weights of the objectives used in the \u201cbalanced\u201d outcome run, to skew them towards their preferences and highlight the differences.<\/p>\n<p><strong>Table 3. Relative weighting as obtained from interviews with Anchovy Bay community members. The table also shows the inverse weighting from Table 1 and the resulting combined weight, which is the product of the community weighting and the inverse weighting.<\/strong><\/p>\n<table id=\"tablepress-18\" class=\"tablepress tablepress-id-18\">\n<thead>\n<tr class=\"row-1\">\n<th class=\"column-1\">Objective<\/th>\n<th class=\"column-2\">Rel. weight<\/th>\n<th class=\"column-3\">Inverse weighting<\/th>\n<th class=\"column-4\">Combined rel. weight<\/th>\n<\/tr>\n<\/thead>\n<tbody class=\"row-striping row-hover\">\n<tr class=\"row-2\">\n<td class=\"column-1\">Economic rent<\/td>\n<td class=\"column-2\">0.56<\/td>\n<td class=\"column-3\">0.5<\/td>\n<td class=\"column-4\">0.3<\/td>\n<\/tr>\n<tr class=\"row-3\">\n<td class=\"column-1\">Social benefits<\/td>\n<td class=\"column-2\">0.75<\/td>\n<td class=\"column-3\">0.5<\/td>\n<td class=\"column-4\">0.4<\/td>\n<\/tr>\n<tr class=\"row-4\">\n<td class=\"column-1\">Ecosystem structure<\/td>\n<td class=\"column-2\">1<\/td>\n<td class=\"column-3\">1.4<\/td>\n<td class=\"column-4\">1.4<\/td>\n<\/tr>\n<tr class=\"row-5\">\n<td class=\"column-1\">Biodiversity<\/td>\n<td class=\"column-2\">0.94<\/td>\n<td class=\"column-3\">5.9<\/td>\n<td class=\"column-4\">5.5<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p><!-- #tablepress-18 from cache --><\/p>\n<\/div>\n<p><strong>Table 4. Objective function values and fleet effort for a retrospective optimization for Anchovy Bay with weights set to the combined weigthings from Table 3. Optimization started with Ecopath base effort.<\/strong><\/p>\n<table id=\"tablepress-19\" class=\"tablepress tablepress-id-19\">\n<thead>\n<tr class=\"row-1\">\n<th class=\"column-1\"><strong>Total<\/strong><\/th>\n<th class=\"column-2\"><strong>Econ.<\/strong><\/th>\n<th class=\"column-3\"><strong>Social<\/strong><\/th>\n<th class=\"column-4\"><strong>Ecosys.<\/strong><\/th>\n<th class=\"column-5\"><strong>Biodiversity<\/strong><\/th>\n<\/tr>\n<\/thead>\n<tbody class=\"row-striping row-hover\">\n<tr class=\"row-2\">\n<td class=\"column-1\">3.98<\/td>\n<td class=\"column-2\">1.33<\/td>\n<td class=\"column-3\">0.81<\/td>\n<td class=\"column-4\">1.38<\/td>\n<td class=\"column-5\">1.04<\/td>\n<\/tr>\n<tr class=\"row-3\">\n<td class=\"column-1\"><\/td>\n<td class=\"column-2\"><\/td>\n<td class=\"column-3\"><\/td>\n<td class=\"column-4\"><\/td>\n<td class=\"column-5\"><\/td>\n<\/tr>\n<tr class=\"row-4\">\n<td class=\"column-1\"><strong>Sealers<\/strong><\/td>\n<td class=\"column-2\"><strong>Trawlers<\/strong><\/td>\n<td class=\"column-3\"><strong>Seiners<\/strong><\/td>\n<td class=\"column-4\"><strong>Bait boats<\/strong><\/td>\n<td class=\"column-5\"><strong>Shrimpers<\/strong><\/td>\n<\/tr>\n<tr class=\"row-5\">\n<td class=\"column-1\">0.05<\/td>\n<td class=\"column-2\">1.21<\/td>\n<td class=\"column-3\">0.03<\/td>\n<td class=\"column-4\">0.46<\/td>\n<td class=\"column-5\">0.62<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p><!-- #tablepress-19 from cache --><\/p>\n<p>Comparing Table 2 and Table 4, it&#8217;s noticeable that adding the strong conversation-oriented view of the Anchovy Bay community has some socio-economic consequences. The economic rent and social benefits performance values were reduced with 20-25%, the ecosystem structure indicator increased ~30%, while the biodiversity measure (which is relative hard to impact) increased with a couple of percent.<\/p>\n<p>If you want to explore the effect that the two optimizations have on fleet catches and values, and on group biomasses, you can extract that information (after each optimization) from the tables on the <em>Ecosim &gt; Output &gt; Ecosim results<\/em>\u00a0form.<\/p>\n<div>\n<h3>Local maxima<\/h3>\n<p>As part of the exploratory analysis, it is important to check whether the maximization search is impacted by the start point, i.e. whether the optimization solutions are unique. By default the optimization routine will start with the fishing rates defined by the Ecopath baseline. It&#8217;s possible, however, to instead using random fleet effort <em>(Random F&#8217;s<\/em> in the policy interface) to check if the optimization routine is likely to get stuck at local maxima. \u00a0All optimization routines are impacted by this, the ones in EwE being no exceptions.<\/p>\n<p>We illustrate this for Anchovy Bay by running 70 optimizations with the combined rel. weights from Table 3, and runs initialized with random F&#8217;s. The outcome of that exploratory analysis is presented in Figure 1.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-4120\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay.png\" alt=\"\" width=\"1000\" height=\"400\" srcset=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay.png 1000w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-300x120.png 300w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-768x307.png 768w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-65x26.png 65w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-225x90.png 225w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-350x140.png 350w\" sizes=\"auto, (max-width: 1000px) 100vw, 1000px\" \/><\/p>\n<p><strong>Figure 1. Box plot<a class=\"footnote\" title=\"The R-code and CSV-file used for this figure can be downloaded from this link.\" id=\"return-footnote-3856-6\" href=\"#footnote-3856-6\" aria-label=\"Footnote 6\"><sup class=\"footnote\">[6]<\/sup><\/a> showing minimum, first quartile, median, third quartile and max value of objective function value for indicators and relative effort for fleets for 70 policy optimizations run with random <em>Starting F<\/em> initialization. \u00a0It&#8217;s clear that only the relative effort for the <em>bait boats<\/em> varies between runs (with one run way higher than the 69 others), and that the optimizations overall are robust to local maxima.<\/strong><\/p>\n<p>From Figure 1, it is clear that for this version<a class=\"footnote\" title=\"Optimizations depend on fitting and model parameters. Different versions of the Anchovy Bay may well lead to different optimizations\" id=\"return-footnote-3856-7\" href=\"#footnote-3856-7\" aria-label=\"Footnote 7\"><sup class=\"footnote\">[7]<\/sup><\/a> of the Anchovy Bay model, there is very little tendency for the optimization to get stuck on local maxima. The variation in the objective estimates and effort patterns are very similar across all runs (apart perhaps from <em>bait boat<\/em> effort), with only a few runs indicating presence of local maxima.<\/p>\n<div class=\"textbox shaded\">\n<p>That this Anchovy Bay model isn&#8217;t prone to get stuck on local maxima should not invite complacency, but be seen as a warning that there <span style=\"text-decoration: underline\">may<\/span> be local maxima. We ran the optimizations 70 times for Figure 1 to illustrate this point \u2013 but also to show that it rarely happens. Bottom line is, run your optimization a number of times (maybe 10 to 20), check how consistent the output is. Then run the model once, check if the outcome corresponds to the majority of the random runs. If it does, you&#8217;re good to go.<\/p>\n<\/div>\n<p>The overall conclusion is that policy optimizations for Anchovy Bay are not very prone to get stuck on local maxima. This is also what we have found for many other ecosystem model optimizations, giving some comfort that the starting point isn&#8217;t very critical. Still, this needs to be checked for all models, so including a search with random <em>Starting F&#8217;s<\/em> should be included in all more serious policy explorations.<\/p>\n<h2>Fleet trade-off analysis<\/h2>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-4013 alignleft\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Anchovy-Bay-FleetTradeoff-e1733339668222-300x296.jpg\" alt=\"\" width=\"300\" height=\"296\" srcset=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Anchovy-Bay-FleetTradeoff-e1733339668222-300x296.jpg 300w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Anchovy-Bay-FleetTradeoff-e1733339668222-1024x1009.jpg 1024w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Anchovy-Bay-FleetTradeoff-e1733339668222-768x757.jpg 768w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Anchovy-Bay-FleetTradeoff-e1733339668222-65x64.jpg 65w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Anchovy-Bay-FleetTradeoff-e1733339668222-225x222.jpg 225w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Anchovy-Bay-FleetTradeoff-e1733339668222-350x345.jpg 350w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Anchovy-Bay-FleetTradeoff-e1733339668222.jpg 1242w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/p>\n<p><strong>Figure 2. Fleet trade-off analysis for Anchovy Bay showing impact a 10% reduction in effort for the fleet listed in rows has on the fleets listed above columns. Negative impacts are in red, and positive in blue.\u00a0<\/strong><\/p>\n<p>A next step of exploratory analysis is the fleet trade-off analysis described in the Fishing policy chapter (<a href=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/chapter\/fishing-policy-exploration#fleet tradeoff\">link to fleet trade-of<\/a>f). We refer to that section for description, including for code to produce plots.<\/p>\n<p>We suggest that you perform the fleet trade-off analysis for your model and explore the outcome. For Anchovy Bay (Figure 2), the plot shows the impact that a 10% reduction for the fleets mentioned to the left, are predicted to have in the value of the catch for all fleets. <span style=\"color: #ff0000\">Red<\/span> circles indicate reduction in value, and <span style=\"color: #0000ff\">blue<\/span> increase. \u00a0The impacts are displayed so that the circle areas are proportional to the changes in value of the catch, and are thus comparable across fleets<a class=\"footnote\" title=\"The monetary value of the fleet-tradeoffs can be obtained from the CSV file used for producing the fleet trade-off plots\" id=\"return-footnote-3856-8\" href=\"#footnote-3856-8\" aria-label=\"Footnote 8\"><sup class=\"footnote\">[8]<\/sup><\/a>.<\/p>\n<p>For Anchovy Bay, the fleet trade-off analysis shows some both straightforward and more complex patterns. Notice for instance that a reduction in <em>bait boats<\/em>&#8216; effort will have a positive impact on <em>seiners<\/em>. That makes sense since the two fleets both catch anchovy. But conversely, a reduction in <em>seiners<\/em>&#8216; effort lead to a small reduction in landed value for the <em>bait boats<\/em>. Why? <em>Seiners<\/em> also catch mackerel, and the reduced effort will lead to more mackerel, which in turn will have a negative impact on anchovy, and hence on <em>bait boats<\/em>, (which only catch anchovy).<\/p>\n<p>The strongest impact of effort reduction is for <em>trawlers<\/em> and <em>shrimpers<\/em>, the two fleets that have the highest value in the model base year. Reduction in <em>trawlers<\/em>&#8216; effort has a considerable negative impact on <em>trawlers<\/em>, but also an almost corresponding negative impact on <em>shrimpers<\/em>. That makes sense as the reduction should lead to more cod and whiting, both of which eat shrimp. \u00a0But conversely, reducing <em>shrimpers<\/em>&#8216; effort leads to an <span style=\"text-decoration: underline\">increase<\/span> in shrimp landings (so they must be overexploited in the baseline \u2013 given the baseline predator-prey conditions though). But why does this not lead to an increase in the value of <em>trawlers<\/em>&#8216; landings? \u00a0The reason is that shrimp have a positive impact on whiting and mackerel, but a negative impact on cod. Notably, the increase in whiting impacts cod negatively. \u00a0Those relationships becomes clearer if you check out the Mixed Trophic Impact analysis <em>(Ecopath &gt; Output &gt; Tools &gt; Network Analysis &gt; Mixed trophic impact &gt; Impact<\/em> data), which shows that <em>shrimpers<\/em> have opposite impact on cod <em>versus<\/em> whiting and mackerel.<\/p>\n<h2>Your policy questions?<\/h2>\n<p>Having explored the behaviour of your model, e.g., as described above, the next issue is to clearly define what questions you are asking for your model. An example of this is Alms et al. 2022, <a class=\"footnote\" title=\"Alms V, Romagnoni G, Wolff M. Exploration of fisheries management policies in the Gulf of Nicoya (Costa Rica) using ecosystem modelling. Ocean and Coastal Management 230 (2022) 106349.\u00a0https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349\" id=\"return-footnote-3856-9\" href=\"#footnote-3856-9\" aria-label=\"Footnote 9\"><sup class=\"footnote\">[9]<\/sup><\/a> who compared output from three defined Ecosim scenarios, (1) ban on shrimp trawling, (2) gill net effort reduction of 25%, and (3) the combination of (1) and (2), and then compared these to the output of single- and multi-objective policy optimizations. Here, the multi-objective optimizations were designed to serve as more balanced solutions. \u00a0See the paper for details.<\/p>\n<h3>Anchovy Bay<\/h3>\n<p>While the complex patterns in the fleet tradeoff (Figure 2) can be explained as above, they raise some questions. There is a big negative impact of reducing <em>trawlers<\/em>&#8216; effort and a big positive impact of reducing <em>shrimpers<\/em>&#8216; effort. \u00a0Is this then what we should explore? Well, it&#8217;s certainly interesting scenarios, but there are some complications, and it is in line with what actually happened in Anchovy Bay. Effort of <em>trawlers<\/em> indeed increased 2-3 times over time, but the shrimp effort increased by almost an order of magnitude. \u00a0Does this make sense when the fleet trade-off analysis indicate that shrimp catches would increase if <em>shrimpers<\/em>&#8216; effort was reduced? \u00a0The finding is not really wrong, but it doesn&#8217;t not consider that the increase in the <em>trawlers<\/em>&#8216; effort reduced the abundance of cod and whiting, which in turn lead to many more shrimps being available for the <em>shrimpers<\/em>. Predator-release!<\/p>\n<p>What the above indicates is in essence that the Anchovy Bay ecosystem of today is very different from that of 1970. Therefore, we will next run the policy optimization for Anchovy Bay as a forward-looking scenario where we keep the run from 1970-2010 as it was in the fitting, and then conduct the policy optimization from 2011 (year 42) onwards. Figure 3 illustrates the setup.<\/p>\n<p>If you change the \u00a0<em>Ecosim &gt; Tools &gt; Fishing policy search &gt; Base year<\/em> to the end of the time series, the routine will automatically black out the years of the time series. In that case, the economic data (cost and value) in <em>Ecopath &gt; Input &gt; Fishery &gt; Fleets<\/em> are assumed to represent your base year, not the Ecopath model year (Year 1 or 1970 for Anchovy Bay).<\/p>\n<\/div>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-4123\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18.png\" alt=\"\" width=\"3416\" height=\"388\" srcset=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18.png 3416w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18-300x34.png 300w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18-1024x116.png 1024w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18-768x87.png 768w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18-1536x174.png 1536w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18-2048x233.png 2048w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18-65x7.png 65w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18-225x26.png 225w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Screenshot-2025-02-05-at-10.41.18-350x40.png 350w\" sizes=\"auto, (max-width: 3416px) 100vw, 3416px\" \/><\/p>\n<p><strong>Figure 3. Fishing policy search interface set up to do a forward scenario for Anchovy Bay. The optimization will not change the effort for the first 41 years (in black), but search for one effort for each fleet for the last 20 years. If you cannot see all the years, try the zoom icons in the interface and widen the aspect ratio of the interface. \u00a0\u00a0<\/strong><\/p>\n<p>We run the forward-looking scenarios with <em>Random F<\/em> drawn to start the optimizations, which will ensure that we do not get stuck in a local maxima for all runs (even though we have found above that the risk of this is small). We here do not use the <em>&#8220;limit cost &gt; earnings&#8221; <\/em>as this tends to hang the optimizations for the Anchovy Bay model<a class=\"footnote\" title=\"This is likely because we for this tutorial are using very rough economic data. If you have this issue, do check your assumptions about cost vs value in your base year.\" id=\"return-footnote-3856-10\" href=\"#footnote-3856-10\" aria-label=\"Footnote 10\"><sup class=\"footnote\">[10]<\/sup><\/a>.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"size-full wp-image-4139 aligncenter\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Objective-plot.png\" alt=\"\" width=\"600\" height=\"400\" srcset=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Objective-plot.png 600w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Objective-plot-300x200.png 300w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Objective-plot-65x43.png 65w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Objective-plot-225x150.png 225w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/Objective-plot-350x233.png 350w\" sizes=\"auto, (max-width: 600px) 100vw, 600px\" \/><\/p>\n<p><strong>Figure 4. Indicator values from Anchovy Bay model policy optimizations with <em>inverse range<\/em> weighting <em>vs<\/em>. <em>community-based<\/em> weightings (from Table 3). The community places higher weights on <em>ecosystem structure<\/em> and <em>biodiversity<\/em>, and the optimizations indeed reflect this. It&#8217;s clear that all of the objectives contributes to the total objective function score, which is the sum of the four objectives \u2013 actually with less variation for the <em>community-based <\/em>weighting. But the improved ecological score is offset by lower employment and profit<a class=\"footnote\" title=\"R-code and data files to produce this code are included in the Zip file that can be downloaded for Figure 5, see below.\" id=\"return-footnote-3856-11\" href=\"#footnote-3856-11\" aria-label=\"Footnote 11\"><sup class=\"footnote\">[11]<\/sup><\/a>. \u00a0\u00a0<\/strong><\/p>\n<p>&nbsp;<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-4130\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-results-X.png\" alt=\"\" width=\"1200\" height=\"800\" srcset=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-results-X.png 1200w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-results-X-300x200.png 300w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-results-X-1024x683.png 1024w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-results-X-768x512.png 768w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-results-X-65x43.png 65w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-results-X-225x150.png 225w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/policy-optim-anchovy-bay-results-X-350x233.png 350w\" sizes=\"auto, (max-width: 1200px) 100vw, 1200px\" \/><\/p>\n<p><strong>Figure 5. Plot<a class=\"footnote\" title=\"Download R-code and CSV-files used for Figure 4 and 5 from this link.\" id=\"return-footnote-3856-12\" href=\"#footnote-3856-12\" aria-label=\"Footnote 12\"><sup class=\"footnote\">[12]<\/sup><\/a> comparing fishing policy optimizations for Anchovy Bay model with weight as defined based on inverse range and community input in Table 3. Value is proportional to jobs, so the <em>Jobs end<\/em> plot indicates around a 25% reduction in jobs with the community-based weighting. \u00a0<\/strong><\/p>\n<p>The fleet effort value, catch and biomass results for the Anchovy Bay model forward-looking policy scenario are presented in Figure 5, again comparing optimizations with the inverse range weightings vs the community opinion based weightings (both from Table 3). The results show the complex trade-off patterns between the two weighting schemes.<\/p>\n<div>\n<h3>Your case<\/h3>\n<p>We suggest that you as a starting point go through this extended tutorial and conduct policy optimizations on a fitted version of your ecosystem model, working your way through the steps described in the tutorial.<\/p>\n<div class=\"textbox shaded\">\n<p><strong>Pretty good yield<\/strong><\/p>\n<p>When you run optimizations, you&#8217;ll often see that the objective function may rather quickly get to more than 98% of the final objective function score. If you start with <em>Ecopath base F<\/em>, the effort changes may be rather limited when you pass the 98% mark, but the last few percent may cause big changes in effort. \u00a0It may be worth exploring those intermediate states (by stopping the optimizations before it finishes by itself).<\/p>\n<p>This is to some degree related to the idea behind Ray Hilborn&#8217;s <em>Pretty Good Yield<a class=\"footnote\" title=\"https:\/\/doi.org\/10.1016\/j.marpol.2009.04.013\" id=\"return-footnote-3856-13\" href=\"#footnote-3856-13\" aria-label=\"Footnote 13\"><sup class=\"footnote\">[13]<\/sup><\/a><\/em> for MSY.<\/p>\n<\/div>\n<p>&nbsp;<\/p>\n<div class=\"textbox textbox--examples\">\n<header class=\"textbox__header\">\n<p class=\"textbox__title\"><strong>Acknowledgement<\/strong><\/p>\n<\/header>\n<div class=\"textbox__content\">\n<div class=\"textbox__content\"><img loading=\"lazy\" decoding=\"async\" class=\"alignright wp-image-3830 size-medium\" src=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/EcoScope-logo-300x113.png\" alt=\"\" width=\"300\" height=\"113\" srcset=\"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/EcoScope-logo-300x113.png 300w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/EcoScope-logo-65x25.png 65w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/EcoScope-logo-225x85.png 225w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/EcoScope-logo-350x132.png 350w, https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-content\/uploads\/sites\/1902\/2024\/09\/EcoScope-logo.png 701w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/>This chapter was developed for the <a href=\"https:\/\/ecoscopium.eu\">EcoScope<\/a> project to guide implementation of the EwE Policy Search for the project case studies. EcoScope is funded from the <a href=\"https:\/\/ec.europa.eu\/programmes\/horizon2020\/en\/home\" target=\"_blank\" rel=\"noopener noreferrer\">European Commission\u2019s Horizon 2020 Research and Innovation programme<\/a> under grant agreement No 101000302. Project coordinator: Aristotle University of Thessaloniki, Greece. \u00a0Parts of the text are from the unpublished EwE User Guide: Christensen V, C Walters, D Pauly, R Forrest. Ecopath with Ecosim. User Guide. November 2008.<\/div>\n<\/div>\n<\/div>\n<p>&nbsp;<\/p>\n<\/div>\n<div class=\"media-attributions clear\" prefix:cc=\"http:\/\/creativecommons.org\/ns#\" prefix:dc=\"http:\/\/purl.org\/dc\/terms\/\"><h2>Media Attributions<\/h2><ul><li >Original       <\/li><li >Original       <\/li><li >Ecosim &gt; Tools &gt; Policy search       <\/li><li >Original       <\/li><li >Original       <\/li><li >EcoScope logo       <\/li><\/ul><\/div><hr class=\"before-footnotes clear\" \/><div class=\"footnotes\"><ol><li id=\"footnote-3856-1\">Illegal, unregulated and unreported <a href=\"#return-footnote-3856-1\" class=\"return-footnote\" aria-label=\"Return to footnote 1\">&crarr;<\/a><\/li><li id=\"footnote-3856-2\">Note that the policy optimization always runs for an additional 20 years without showing the results. This is to ensure that the optimizations do not result in an \"empty sea\" where fleets are encouraged to fish out the resources (which might happen if there was no tomorrow to consider). <a href=\"#return-footnote-3856-2\" class=\"return-footnote\" aria-label=\"Return to footnote 2\">&crarr;<\/a><\/li><li id=\"footnote-3856-3\">This approach is used commonly for climate change scenarios, but may never have been used for policy optimizations. <a href=\"#return-footnote-3856-3\" class=\"return-footnote\" aria-label=\"Return to footnote 3\">&crarr;<\/a><\/li><li id=\"footnote-3856-4\"><em>B<sub>lim<\/sub><\/em> is the lowest acceptable biomass and <em>B<sub>0<\/sub><\/em> the Ecopath biomass for the group. The objective is thus entered as a relative biomass term.  <a href=\"#return-footnote-3856-4\" class=\"return-footnote\" aria-label=\"Return to footnote 4\">&crarr;<\/a><\/li><li id=\"footnote-3856-5\">e.g., Natugonza et al. (2020) <a href=\"https:\/\/doi.org\/10.1016\/j.fishres.2020.105593\">https:\/\/doi.org\/10.1016\/j.fishres.2020.105593<\/a>, and Alms et al., (2022) <a href=\"https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349\">https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349<\/a> <a href=\"#return-footnote-3856-5\" class=\"return-footnote\" aria-label=\"Return to footnote 5\">&crarr;<\/a><\/li><li id=\"footnote-3856-6\"><\/strong>The R-code and CSV-file used for this figure can be downloaded from this <a href=\"https:\/\/ln5.sync.com\/dl\/3065219c0\/gb3t8j6n-pdmu9gh7-yqqgnusb-jh89zs6a\">link<\/a>.<strong> <a href=\"#return-footnote-3856-6\" class=\"return-footnote\" aria-label=\"Return to footnote 6\">&crarr;<\/a><\/li><li id=\"footnote-3856-7\">Optimizations depend on fitting and model parameters. Different versions of the Anchovy Bay may well lead to different optimizations <a href=\"#return-footnote-3856-7\" class=\"return-footnote\" aria-label=\"Return to footnote 7\">&crarr;<\/a><\/li><li id=\"footnote-3856-8\">The monetary value of the fleet-tradeoffs can be obtained from the CSV file used for producing the fleet trade-off plots <a href=\"#return-footnote-3856-8\" class=\"return-footnote\" aria-label=\"Return to footnote 8\">&crarr;<\/a><\/li><li id=\"footnote-3856-9\">Alms V, Romagnoni G, Wolff M. Exploration of fisheries management policies in the Gulf of Nicoya (Costa Rica) using ecosystem modelling. Ocean and Coastal Management 230 (2022) 106349.\u00a0<a href=\"https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349\">https:\/\/doi.org\/10.1016\/j.ocecoaman.2022.106349<\/a> <a href=\"#return-footnote-3856-9\" class=\"return-footnote\" aria-label=\"Return to footnote 9\">&crarr;<\/a><\/li><li id=\"footnote-3856-10\">This is likely because we for this tutorial are using very rough economic data. If you have this issue, do check your assumptions about cost vs value in your base year. <a href=\"#return-footnote-3856-10\" class=\"return-footnote\" aria-label=\"Return to footnote 10\">&crarr;<\/a><\/li><li id=\"footnote-3856-11\"><\/strong>R-code and data files to produce this code are included in the Zip file that can be downloaded for Figure 5, see below.<strong> <a href=\"#return-footnote-3856-11\" class=\"return-footnote\" aria-label=\"Return to footnote 11\">&crarr;<\/a><\/li><li id=\"footnote-3856-12\"><\/strong>Download R-code and CSV-files used for Figure 4 and 5 from this <a href=\"https:\/\/ln5.sync.com\/dl\/6be981c90\/djen2t7n-keic5kps-6u9scv7p-b48hnx2r\">link<\/a>.<strong> <a href=\"#return-footnote-3856-12\" class=\"return-footnote\" aria-label=\"Return to footnote 12\">&crarr;<\/a><\/li><li id=\"footnote-3856-13\"><a href=\"https:\/\/doi.org\/10.1016\/j.marpol.2009.04.013\">https:\/\/doi.org\/10.1016\/j.marpol.2009.04.013 <\/a> <a href=\"#return-footnote-3856-13\" class=\"return-footnote\" aria-label=\"Return to footnote 13\">&crarr;<\/a><\/li><\/ol><\/div>","protected":false},"author":1909,"menu_order":10,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":["villy","santiago-de-la-puente"],"pb_section_license":""},"chapter-type":[49],"contributor":[66,60],"license":[],"class_list":["post-3856","chapter","type-chapter","status-publish","hentry","chapter-type-numberless","contributor-santiago-de-la-puente","contributor-villy"],"part":438,"_links":{"self":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapters\/3856","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":26,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapters\/3856\/revisions"}],"predecessor-version":[{"id":4194,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapters\/3856\/revisions\/4194"}],"part":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/parts\/438"}],"metadata":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapters\/3856\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/wp\/v2\/media?parent=3856"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/pressbooks\/v2\/chapter-type?post=3856"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/wp\/v2\/contributor?post=3856"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/ewemodel\/wp-json\/wp\/v2\/license?post=3856"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}