A Concept of Bayesian Regulation in Fisheries Management

Stochastic variability of biological processes and uncertainty of stock properties compel fisheries managers to look for tools to improve control over the stock. Inspired by animals exploiting hidden prey, we have taken a biomimetic approach combining catch and effort in a concept of Bayesian regulation (BR). The BR provides a real-time Bayesian stock estimate, and can operate without separate stock assessment. We compared the performance of BR with catch-only regulation (CR), alternatively operating with N-target (the stock size giving maximum sustainable yield, MSY) and F-target (the fishing mortality giving MSY) on a stock model of Baltic Sea herring. N-targeted BR gave 3% higher yields than F-targeted BR and CR, and 7% higher yields than N-targeted CR. The BRs reduced coefficient of variance (CV) in fishing mortality compared to CR by 99.6% (from 25.2 to 0.1) when operated with F-target, and by about 80% (from 158.4 to 68.4/70.1 depending on how the prior is set) in stock size when operated with N-target. Even though F-targeted fishery reduced CV in pre-harvest stock size by 19–22%, it increased the dominant period length of population fluctuations from 20 to 60–80 years. In contrast, N-targeted BR made the periodic variation more similar to white noise. We discuss the conditions when BRs can be suitable tools to achieve sustainable yields while minimizing undesirable fluctuations in stock size or fishing effort.


Introduction
Fisheries managers are challenged with two widely permeated properties of their study system, uncertainty [1][2][3] and variability [4][5][6].These can confound each other, e.g.imprecise spawning stock size estimates can generate apparent variability in the stock 2 recruitment relationship [7].Reduced inter-annual variability in effort has socio-economic benefits with yields matching the capacity of the processing industry [8], a more stable job market, and fewer years with over-dimensioned fleets [9].Often there is a trade-off between maximizing yield and stabilizing yield and fishing effort that makes management objectives ambiguous [10][11][12].There are also concerns that fishing increases the temporal variability of harvested stocks, [4,[13][14][15].
The problem with temporal variability in fisheries has led to an increased interest in the performance of alternative harvest control rules [16][17][18].Stephenson et al. [19] argue that the total allowable catch (TAC) could be used to prevent unsustainable use, but is insufficient to control spatial and temporal variability, and cannot be used to achieve socioeconomic objectives.In a reflection over proposed alternative harvest controls, May et al. [20] conclude that further mathematical refinement is probably not as important as developing ''robustly self-correcting strategies that can operate with only fuzzy knowledge about stock levels and recruitment curves''.If our belief of the stock size is a probability function (in contrast to a point estimate), Bayes' theorem postulates that harvesting information can be used to calculate a conditional probability function [21].Here we propose Bayesian regulation (BR) that uses catch-and effort data from the ongoing fisheries to make real-time estimates of stock size and fishing mortality.The years can be linked by using the posterior distribution as a prior in the sequential year, and hence exploitation is combined with stock size assessment.Real-time assessment of BR would be advantageous to management routines relying on forecasts of stock abundances.Given that most commercial fisheries with CR-based control of biomass and fishing mortality (e.g.TAC) routinely monitor effort and catch, this additional information is surprisingly poorly utilized during the CR fishing season.We suggest a methodology where this information can be used for updating population size estimates to facilitate in-season management decisions.
The BR is derived from Bayesian foraging theory in behavioral ecology, which describes a giving-up rule for patch foraging animals [22][23][24].The rule is a relaxation of the full information assumption of the marginal value theorem [25], with search time and number of prey caught as information variables.It applies to patches in which the prey are hidden, such as a woodpecker feeding on pupae under bark [26,27].Given some prior information available to the forager, e.g.experience from foraging in patches of its territory, a Bayesian posterior distribution can represent the forager's continuously changing belief about the prey density in the current patch [28].The animal maximizing its intake rate would thus leave the patch for a new one when the anticipated intake rate in the current patch drops below the one expected from patches on average in the territory.The expected number of remaining prey in a patch can be expressed with a fairly simple equation derived by Iwasa et al. [29], but the decision to leave should include a discounting of the value of further bits of information [30].
For fisheries applications, we have modified Iwasa's equation for Bayesian prey density estimation for recurrent exploitation of one population.Both are cases of Bayesian update of mean and variance of a population size estimate.In the original scenario, the exploiter is informed by the mean and variance of several exploited patches within its foraging area.In the fishery scenario, the exploiter (manager) is informed by previous exploitations and surveys in the same area, followed by quantitative assessment and stochastic forecast simulations.We compare CR with BR, and show how they perform in relation to management objectives and targets, by simulating fishing on a model of the main basin Baltic Sea herring (Appendix S2) [31].We present levels and temporal variability in yield, fishing mortality, stock abundance, spawning stock biomass (SSB), and finally how frequently it surpasses the MSY related reference point B trigger [32].

Methods
We look at three levels of decision-making in fisheries management as they are described in the common fisheries policy of the Council of the European Union [33]: (i) The management objective, which is the ultimate goal of fisheries management.The objective can be simple or more complex, for example weighing incompatible goals [34,35].We have chosen the maximum sustainable yield (MSY) because it is the current objective of fisheries management in the EU [36].(ii) The target of The subscript denotes type of target and the superscript the existence of supervision.The target column denotes the targets used: F MSY is the fishing mortality giving highest yield in the MSY-analysis of the SOM, and N' MSY is the corresponding post-harvest stock size.The management is either supervised with a separate assessment or unsupervised (see methods).For unsupervised management we use the mean, m, and the variance, v, of the pre-harvest stock size from MSY-analyses: N MSY .In supervised management m is instead the actual pre-harvest stock size, N, with an added randomly generated error term e drawn from the normal distribution, (mean = 0, var = var(N MSY ) where var(N MSY ) is taken from an MSY-analysis).doi:10.1371/journal.pone.0111614.t001 Figure 1.(A) Bayesian regulation curves indicating the combined total catch (y-axis) and total effort (x-axis) when fishing should be terminated for the season.The curves are given by equation 1 when r is set to the target number of fish, and solved for total catch.Hence, the curves show when the fishery target is achieved in catch-effort space, given by the ratio of variance (v) to mean stock size (m).There are four qualitative cases of the v-m ratio: when the ratio is zero, the fishery is regulated by catch only because there is no uncertainty around the mean.When the ratio equals 1, regulation is by effort only.Apart from these special cases, the fishery should be regulated by both catch and effort.If the variance is lower than the mean, the posterior mean will decrease with accumulating catches and effort.If the variance is higher than the mean, the posterior mean will increase with increasing catches but decrease with effort.These effects can be deduced from analyzing equation 1. (B) Here the solid curve denotes the Bayesian regulation when v.m, which is determined prior to the fishing season (see Methods).The grey curves are idealized trajectories of how cumulative catch and effort develops from the origin during the fishing season.When the trajectories cross the regulation curve, the Bayesian posterior is on target, and the fishery should be closed for the season.The hatched grey line indicates when the initial stock size is at the prior mean (m), the solid grey line when the stock size is at the mean plus one standard deviation (see F MSY -analysis).The BR thus allows larger catches when the stock size has been underestimated, and vice versa.doi:10.1371/journal.pone.0111614.g001 exploitation, which should be given in a quantitative unit that relates to the management objective.We compare the efficiency of two targets: fishing mortality (F MSY ) and post-harvest stock size (N MSY ).(iii) The harvest rules, which ''[lay] down the manner in which annual catch and/or fishing effort limits are to be calculated and provide for other specific management measures, taking account also of the effect on other species'' [33].On this level we choose to use other terms.We use the term regulation type to denote the use of catch and/or effort regulation, whereas the term harvest control rule (HCR) is used for the combined objective, target and regulation type.
The harvest control rules under development by the International Council for the Exploration of the Sea (ICES) involve a fishing mortality target (F*) and the spawning stock biomass reference point B trigger [32,37].If the SSB goes below B trigger the fishing mortality is reduced below target, a situation that should warrant a reinvestigation of the stock condition and the harvest control [32].Based on a stock assessment and a short-term forecast, the fishing mortality target is recalculated as a TAC, which is how the recommendation is presented to the European Commission [38].ICES harvest control rules use precautionary reference points below which fishing mortality is reduced [32], something we do not use in our simulations.We have focused on the principle differences between CR and BR, e.g. the elaborate data collection and assessment procedure is simplified by using the actual size of the simulated stock with a random error.

Harvest control
Altogether we compared eight different harvest controls achieved by the combinations of two exploitation targets, two regulation types and two stock assessment modes (Table 1).The Bayesian regulation BR(c, f|m, n) is a function of catch (c) and effort (f) given the prior information of the estimated stock size before harvest, m, and the uncertainty of that estimate given as variance, n.The Bayesian estimate of the number of remaining individuals (r) in a harvested stock is: where q, the catchability, is defined as the fraction of the population captured by one unit of fishing effort and is the scalar  between catch per unit effort indices and average population abundance, hence q~c f N .When fishing is closed by fulfillment of the target condition, f = f *, then the fishing mortality F = q f *.This is a generalization of Iwasa's [29] Bayesian estimates of remaining prey population, see Appendix S1 for details on how Equation ( 1) is derived from some specific probability distributions.With a known prior probability distribution and a given catch and effort, a posterior distribution can be calculated with r being the mean.It is very unlikely that the prior mean is equal to the actual population size.The posterior mean will approach the actual value with accumulating catch and effort, but there will always be a bias towards the prior [39].For large deviations of the prior mean in relation to the real population size, it may take more than one year to track the population size more closely.Hence, there can be temporal correlations in estimate biases.This depends on the random perturbations displacing the population size from the prior mean and the harvesting information making the posterior mean approach the real value.When Equation ( 1) is solved for catch as a function of f, and r is given a target value, it can be visualized as harvest control curves (Fig. 1A).These curves can be calculated prior to the fishing season and define the combination of total catch and total effort at which the Bayesian information indicates that the fishery has reached its target.The real-time accumulation of catch and effort during the fishing season is responsive to over-and underestimation of stock size and will develop different trajectories.The intersection of the real-time accumulation of catch and effort with the BR-curve gives the total catches when an initial over-or under-estimation is compensated for (Fig. 1B).Note that if we have no uncertainty in our prior estimate, i.e. v = 0, Equation ( 1) simplifies to r = m2c, in other words the remaining individuals in the population is our prior estimate minus the catch size.This means that catch-only regulation is a special case of BR (Fig. 1A).To be used in fisheries, Equation (1) needs to account for the simultaneous removal by predators (M), which is done by introducing two scaling parameters a and b: where parameter a~1zM=q q simply scales up the effort with the total mortality in proportion to q q.The notation q q denotes the assumed value of the catchability.As default q q is constant and unbiased, but we make separate runs to explore the impact of biases and stochastic errors in q q.Catchability is difficult to determine accurately and can be affected by aggregation behavior of fish and technical enhancement of fishing gear [40].Parameter b scales the catch c to the total number of casualties from fishing and natural mortality: Note that b depends on effort such that Equation ( 2) can work as a real-time estimator of the population size.The Bayesian estimator is used with the different targets (Table 1).When the target is formulated as a post-harvest population size (N' MSY ), harvesting should stop when the Bayesian estimate of remaining stock size (r) equals target N' MSY : The target N' MSY denotes the stock size at the end of the year, we therefore need to take into account the removal of individuals due to natural mortality for the remaining part of the year after fishing has ended, e { 1{f ð Þ M .When our target is defined as fishing mortality, F MSY , harvesting is aborted when.
Note that the prior m cannot be used in the denominator, because the estimated size of the pre-harvest population changes as we receive information from catch and effort.When the conditions in Equation (4), or (5) are fulfilled, we extract the fishing effort (f*) from Equation (2).Equations ( 2), ( 4) and ( 5) together define the control rules both for CR and BR.

Assessment and priors
Estimates of the pre-harvest stock size, the prior m, are produced in two ways: by supervised assessment and unsupervised assessment (Table 1).Supervised assessment represents the current annually revised assessment practiced by ICES and other stock assessors.In this case, m is calculated by adding an assessment error to the actual pre-harvest stock size, N. The error is a random value drawn from a normal distribution with mean = 0, and variance v which is equal to the variance of the pre-harvest stock size in the F MSY -analysis (see below).The coefficient of variance in the F MSY -analysis is 21%.This represents an ideal assessment where the error solely stems from the yearly variation of the population when fished at F MSY .The unsupervised assessment uses m equal to the average, and v the variance of the pre-harvest stock size from the F MSY -analysis (i.e. when fishing at F MSY ; see below).The rationale of using a constant m is that the Bayesian information of stock size when a fishery has been closed in the foregoing year indicates N' MSY at the end of the year when using N-target (Equation 4).Similarly for an F-targeted fishery, fishing is closed every year when at estimated F MSY , which is associated with the mean N' = N' MSY .The variance is used for the BR, whereas v = 0 for the CR.

The stochastic operating model
We have used a stock model of the herring population in the main basin of the Baltic Sea (ICES catch area subdivisions 25-27, 28.2, 29 and 32) for harvest control evaluation.The model is a stochastic operating model (SOM) parameterized from statistical analyses of ICES catch data and outputs from ICES XSA runs (Appendix S2).

MSY analysis
We performed an MSY-analysis on the SOM by stepping the fishing mortality in steps of 0.01, and for each F-value simulating 40,000 years and rejecting the first 500 to minimize the effects of initiation values.In contrast to simulated management, this algorithm executes perfectly-controlled constant fishing mortality.The relationship between yield and fishing mortality was used to identify the F MSY and its associated average post-harvest stock size.These were used as targets for F-targeted and N-targeted management, respectively.The average pre-harvest stock size, given fishing mortality F MSY , was used as the prior (m) in the unsupervised simulations, and the variance in the pre-harvest stock size was used as a measure of the uncertainty (v) of m in the BR.The lower 2.5% percentile of the associated SSB was used as the B trigger .Targets, priors and assessment modes were used in  defining the eight different HCRs, all aiming at achieving MSY (Table 1).

Evaluation of harvest control rules
The eight HCRs were applied to Monte Carlo simulations of the SOM.A simulation started with a number of years of controlled fishing with F MSY to get away from the initial population size and age-structure.After the initial period, data was collected for a number of years of applied HCR.We ran shorter time series with an initiation period of 1,000 years and 100 years of applied HCR.The time period of 100 years reflects a reasonably long time-horizon for management.We also ran longer series with an initiation period of 100 years and 20,000 years of applied HCR, in order to establish more accurate estimates of mean yields.We collected annual data on yield, fishing mortality, pre-harvest and post-harvest population size, and SSB.Cod SSB and year-specific growth were kept at a mean level (C y = 100,000 tons, k y = 0) during the simulations with the addition of random noise, the size of which was extracted from historical data after removing long-term changes [31].Technically, fishing on the SOM is performed by applying the effort, f*, from the solution of the quitting rules in Equation ( 4) and ( 5).f* is calculated by halved-distance iterations until the estimate (the left side of Equation 4and 5) differed from the target with less than 0.1%.The population is harvested by stepping the SOM one year with the catch being determined by Baranov's catch equation: in which we use the actual q of the SOM.

Sensitivity to error in catchability
We keep the estimated catchability q q in Equation ( 2) and (3) constant and unbiased in our base runs, but fishery dependent catchability may in reality be biased and vary over time.Such uncorrelated errors and biases in catchability will contribute to error or bias in the assessment.Our base model already has assessment error (CV = 21%, as described in Assessment and priors) for supervised HCRs, but unsupervised BRs do not (Table 1).Unsupervised BRs do not use separate assessments, but on the other hand they will be affected by error or bias in q q.We chose to compare supervised F-targeted CR and unsupervised F-targeted BR because the former is only affected by the error in  the assessed stock size, m, whereas the latter is affected only by the error in q q.We also explored the sensitivity to correlated biases in q qand m by changing their values 620%.Simulations ran for 41,000 generations and data from the last 40,000 were used for the analyses.

Statistical analyses
Unsupervised CR expectedly led to population crashes within a few decades or less.We therefore excluded them in the statistical analyses, which were performed using STATISTICA software (Statsoft, Tulsa, Oklahoma).General linear models (GLM) were used for three-way ANOVAs, testing the effects of supervision mode, regulation type and target type on yield, fishing mortality, spawning stock biomass, and post-harvest population size.Fourier analyses were performed on 1,000 year data series, tapered by 15% and padded to the length power of 2. The period length in years with the highest spectral density was identified after applying Hamming weights with a data window of size 7. Test of significant deviation from white noise was performed using Kolmogorov-Smirnov deviation (d) statistics, Table Y, r = 0.5, n.100 in Rohlf & Sokal [41].

Yield
Mean yield in 100-year simulations is significantly higher with BR than with CR although the differences are small (Fig. 2, Table 2, 3).With an increasing length of the time series the confidence intervals narrow (Fig. 3).Analysis of 20,000 years reveals more clearly that BRs operated with N-target give the highest average yields, even though the largest differences in means are less than 10% (Table 2, 4).Supervised and unsupervised BRs give higher yields (262 thousand tons) than the reference MSY from the F MSY -analysis (254 thousand tons; Table 2, 3, 4, 5).Contrary to the small differences in mean yields, the CV in yields is clearly affected by the choice of HCR.In general, using Ftargets gave less temporal variation than using N-targets (Fig. 2, Table 2, 3, 4, 5).There is also an effect of regulation type with BR, giving less variation than CR (Fig. 2, Table 3, 4).However, supervised BR and unsupervised BR have the same CV regardless of target (Table 5).

Fishing mortality
There is no significant difference in mean fishing mortality (F) due to any effect over 100 years (Table 6, Table 7).All HCRs reach the F-target of 0.17 as their mean, except supervised CR with N-target, which has a mean F of 0.20.Although the mean Fs are very similar across HCRs, the differences in CV are more pronounced.Supervised CR with N-target exhibits the highest CV of 158% (Table 6).F-targeted supervised CR has a CV of 25% (Table 6, Fig. 4).Since catchability is constant here, F is proportional to fishing effort.It is therefore not surprising that F-targeted HCRs result in less variation than N-targeted HCRs.F-targeted BRs, supervised and unsupervised, show very high precision in reaching the target (Fig. 4), the CV being as low as 0.1% (Table 6).As expected, N-targeted BRs exhibit higher CVs (68% unsupervised and 70% supervised; Table 6).In addition, it is only the CV between supervised and unsupervised BRs that does not differ significantly from each other (Table 8).

Post-harvest population size and SSB
Looking at post-harvest population size (N'), the harvest controls operating on N-targets are aiming for 41.5 billion individuals at the end of the year.The two N-targeted BRs reach this target with little variation between years (4% CV; Fig. 5, Table 9, 10, 11), whereas N-targeted CR leads to over-fishing in the sense that N' is on average 17% below target.In contrast to the N-targeted regulation types, there is a small difference in N' between Ftargeted BR and CR.This explains the significant interaction between regulation type and target type in the ANOVA (Table 10).Supervision has no significant effect on the mean N' (Table 10), nor the CV (Table 11).
The mean SSBs differ significantly by regulation type, harvest type and supervision (Table 12).BR SSBs are on average larger than CR SSBs (Fig. 6, Table 13, 13).Supervised CR with Ntarget has the lowest mean and is below B trigger in 50% of the years (Table 13, Fig. 6).Supervised BRs has on average larger SSBs than unsupervised, and the BRs operating with N-targets surpass the B trigger less frequently than when operating with F-targets (Table 13, Fig. 6).The lower CV for BRs operating with N-target (Table 13, 14) seems to be the main reason for this result.

Bias and error in harvest rates
So far we can summarize the difference between CR and BR as being BR's ability to reduce the variance in the target variable (F or N'; Tables 3, 4).Yield is partially dependent on F, and SSB on N', and thus the variance of these variables is also reduced, although to a lesser degree (Tables 2, 5).These results are based on the assumption of a constant and un-biased estimate of catchability (q q).When random error is added to the catchability (q) of the same magnitude as the error in m, the CV in c and F for the unsupervised F-targeted BR become similar to supervised Ftargeted CR (Fig. 7).The two F-targeted harvest rules also perform very similar in response to the off-set bias in m and q q.Hence, supervised CR and unsupervised BR behave similarly when error and bias in m and q are of the same relative magnitude.

Temporal variation in population size
We found that all HCRs decrease the amplitude of temporal variation in pre-and post-harvest population numbers, except for post-harvest numbers under supervised N-targeted CR, at which the variation is indifferent from the unexploited population (Table 15).The BR is effective in reaching the N-target, and hence reduces the CV markedly in the post-harvest population size (Fig. 8, Table 15).The effect is still evident in the pre-harvest population size, in which the inclusion of recruits from the foregoing year increases the variation (Table 15).
The unexploited population exhibited a very pronounced periodicity of about 20 years (Fig. 8, Table 16).There is only uncorrelated variation in our operating model, which means that any periodicity is due to the demographic parameters and structure of the population.The BRs with N-target remove much of this periodicity and puts it closer to white noise (Lower K-S dvalues, Table 16).HCRs with F-target increase the dominating period length to about 70-80 years (Table 16).Supervised CR with F-target, a simplified version of the HCR currently used,  reduces the amplitude of variation in both pre-and post-harvested population size (Table 15).At the same time, it increases the dominating period length from 20 years in the unexploited population to about 80 years in the exploited (pre-harvest population size, Table 16).

Discussion
We have suggested a new approach to fisheries regulation, which we name Bayesian regulation (BR), and made a basic evaluation of its properties in relation to catch-only regulation.The idea of BR is adopted from Bayesian foraging theory, in which patch-foragers are assumed to rely on a giving-up rule based on prior information and sampling information [29,30,42].The Bayesian foraging theory has been studied extensively, including the derivation of the probability function used and that cumulative catch and effort are sufficient information variables [29].The BR uses the uncertainty of the stock size estimate as a parameter, which is rarely the case in current management [43][44][45], and which has been requested for some time [46].The general applicability of the BR is emphasized by the fact that catch-only regulation, (CR), is a special case of BR _ i.e. when there is no uncertainty in the prior estimate of the stock size from an assessment or survey.We have analyzed three aspects of harvest controls: (i) the relevance of supervision (with or without stock size assessment), (ii) operational target type (constant fishing mortality, F, or constant stock size, N), (iii) type of regulation (catch only = CR, or Bayesian evaluation of combined catch and effort = BR), and their effect on mean and variance in yield and stock size.

The effect of fishing on stock dynamics
There is an ongoing discussion about whether fishing leads to increased temporal variability in stock size compared to unexploited populations [4,14].Observations from California, USA, corroborate the view that fishing magnifies fluctuations in fish abundance, driven by increased intrinsic growth rates as a response to higher mortality [4].Models without age-structure show an increased variability in population abundance as a consequence of exploitation [14].In a structured model, fishing can either increase or decrease the stock size variability [15].We find that the CV of population size has been reduced in all investigated HCRs, and especially with N-targeted BR.The different effects of fishing between our and the previous studies   Figure 7.A comparison of yield, fishing mortality (F), post-harvest stock size (N') and spawning stock biomass (SSB) when random errors of the same CV are applied to both the stock size assessment (m) and the estimated catchability parameter (q q).Relative positive and negative biases applied to q q and m are denoted on the x-axis.Circles denote results from F-targeted unsupervised BR and squares F-targeted supervised CR.Black symbols denote the value (left axis) and grey symbols the CV of the effect-variable (right axis).doi:10.1371/journal.pone.0111614.g007may depend on our implementation of errors in assessment, whereas earlier studies implemented perfect fishing mortality.

The comparison of Bayesian-and catch-regulation
We expected N-targeted BR to give the highest mean yields, because of its ability to maintain the population in the most productive state.It turned out to return the highest mean yield, but the differences between HCRs are small.The CV of fishing mortality and the population size is, however, much more affected using BR compared to CR.The advantage of F-targeted BR is seen in a reduced variation in fishing mortality and yield.Similarly, the N-targeted BR is associated with less variation in post-harvest population size and SSB.Under some conditions CR can be a blunt harvest control instrument; populations inevitably go extinct if not supervised, and it exhibits the largest CV in yield, F and SSB of all harvest controls when operated with N-target.For the better, F-target supervised CR is a commonly applied harvest control, e.g.within the European Union, and exhibits a smaller CV than N-targeted CR in all fisheries and stock variables investigated.Choosing BR, a manager has to consider the tradeoff between the pros and cons of an Fversus an N-targeted fishery, since targeting one variable pushes the variability to the other [47].From a fisheries perspective, F-targeted BR can appear as the most preferable HCR: low variability in effort (costs) and yield (income) enables a more stable economy, and an even process load is more manageable for the food industry [47].At the same time, single-stock F-targeted BR can potentially yield a conflict of interest between fishers, and may be less desirable when taking into account the ecological interaction between different species of fish.N-targeted BR, on the other hand, gives low variation in population size, which means lower risks of extinction and smaller indirect effects of variability on prey and predator species, which may also be subject to commercial fishing.The stable post-harvest population size of N-targeted BR, in combination with the inherent variation in the productivity of the population causes a highly variable yield.Our results show that F-targeted BRs come close to resolving the trade-off between a high long-term yield and the stable fishing effort and yield acknowledged by some authors [20,35].Although giving a slightly lower average yield than BRs with N-target, the yield is higher than with CR and the fluctuations in the yield are reduced considerably.
The ability of the BR to control the target variable depends on the reliability of the estimated fishery-dependent catchability.We have shown, as a rule of thumb, that if the CV of the error in catchability is as large as the error CV in the stock size assessment, the unsupervised BR performs similar to the supervised CR.Lack of knowledge makes it difficult to say in general which of the errors is the largest, in fishery-dependent catchability or in assessments.In the case of the Newfoundland cod, CPUE increased when stock abundance declined due to hyper-aggregation of the fish [40].Such deceitful behavior of CPUE warrants close attention, and neither BRs nor CRs are immune to it.It is possible to estimate catchability accurately when we know the fishing gear, the fishers' behavior, and the fish dispersal pattern.One could say that a transfer from CR to BR would require a shift of focus by managers from surveys and assessments to gear efficiency, behavior of fish and fishing fleets.

HCRs' effects on period lengths of stock size variation
Not only the amplitude of population variation, but also the frequency of the variation has consequences for the population [48].The induced periodic population variability in predator-prey interactions is well known in ecology [49], but the result of the dynamic interaction between the fish stock and the harvest control is much less studied [15].It has been hypothesized that populations truncated by fishing have a shorter inherent periodicity, and track environmental fluctuations more closely.Given the assumption that unharvested populations fluctuate at longer frequencies than the environment, this could be the reason for increased amplitude of population variation [4].Since populations have inherent dynamics controlled by age-specific demographic parameters, one can conclude that populations truncated by harvest should exhibit fluctuations of shorter period length than unexploited ones.In our study, this is true for the N-targeted BRs, but for the other harvest controls the period length of the fluctuations increase.For F-targeted fisheries, stock size deviations from equilibrium delays return times because of a positive feedback: enlarged populations face reduced mortality and reduced populations are subject to increased mortality.If stocks exhibit long-wave fluctuations with period lengths of 60-80 years due to fishing as indicated with F-targeted CR fisheries, the stress it imposes on ecosystems warrants more attention in future studies and consideration in the choice of method for harvest control.

Potential application of Bayesian regulation
The Bayesian regulation of fisheries modeled here is highly idealized.Real, well-informed management processes are far more complex, and may have sources of error not accounted for.We want to emphasize the potential of BR and stimulate ideas in how it might be implemented.The concept of BR can be developed further, for example to include a Bayesian estimate of the catchability parameter in the probability function (Equation 1).BR has the advantage that, within a year, it can give much higher precision in reaching the management target compared to CR.For instance, stock control with quantitative assessment and forecast can utilize a Bayesian HCR for the harvest year.New assessment methods that calculate a confidence interval for stock size estimate, such as SAM (Nielsen and Berg unpublished), is ideal as they provide priors for a Bayesian regulation.
Another useful property of BR is the ability to work without other assessments, since it is an assessment in itself.It can be an option for the more than 99% of fished species which lack assessments because of the costs and other requirements for collecting the necessary data [50].An application of BR would still require daily reporting of catches, effort and gear used.But with such a system in place, Bayesian estimates of stock size and fishing mortality could be provided daily, with fishers and managers more up-to-date with current stock status.
predator spawning stock biomass, and other parameters are set in this file.(R)

Figure 2 .
Figure 2. Yield over 100 years for different HCRs.The dashed reference line is the MSY = 254.4thousand tons obtained from the F MSY -analysis during controlled F. doi:10.1371/journal.pone.0111614.g002

Figure 3 .
Figure 3. Average yields and confidence intervals as functions of increasing length of simulation in years.Filled symbols denote supervised and open symbols unsupervised harvest control.Squares denote CR and circles BR.Solid lines denote N-target management whereas hatched lines denote F-targeted management.doi:10.1371/journal.pone.0111614.g003

Figure 4 .
Figure 4. Fishing mortality (F) for 100 years for different HCRs.The dashed reference line is the F MSY = 0.17 obtained from the F MSY -analysis during controlled F. doi:10.1371/journal.pone.0111614.g004

Figure 5 .
Figure 5. Post-harvest N-values for different HCRs in simulations over 100 years.The dashed reference line is the N-target = 41,491 million obtained from the F MSY analysis during controlled F. doi:10.1371/journal.pone.0111614.g005

Figure 6 .
Figure 6.Spawning stock biomass (SSB) for different HCRs in simulations over 100 years.The dashed reference line is the B trigger = 988 thousand tons.It is the 2.5% lower percentile of SSB obtained from the F MSY -analysis during controlled F. doi:10.1371/journal.pone.0111614.g006

Figure 8 .
Figure 8. Spectral densities of Fourier analysis of 1,000 years of pre-harvest population size under various HCRs.Period lengths from 0 to 100 years are shown (depicted on the abscissa).doi:10.1371/journal.pone.0111614.g008

Table 1 .
Characteristics and settings for the evaluation of eight different harvest controls: CR being a simple catch-only regulation; BR is the proposed regulation combining catch and effort.

Table 2 .
Mean, standard deviation and CV (%) of annual yield.

Table 3 .
ANOVA-table showing effects of supervision, regulation type, target type, and target type interacting with supervision/regulation type on annual yield.

Table 4 .
ANOVA-table showing effects of supervision, regulation type, target type, and target type interacting with supervision/regulation type on annual yield.

Table 5 .
Pair -wise test of CV in yield for the 100 year series, applying the F-statistics of testing differences in the variance of annual yield.F-values are presented in the lower left triangle and p-values in the upper right.doi:10.1371/journal.pone.0111614.t005

Table 7 .
ANOVA-table showing effects of supervision, regulation type, target type, and target type interacting with supervision/ regulation type on annual fishing mortality. doi:10.1371/journal.pone.0111614.t007

Table 6 .
Mean, standard deviation and CV (%) of annual fishing mortality (F) of 100 years of simulated fishing.

Table 8 .
Pair-wise test of CV-values applying the F-statistics of testing differences in the variance of annual fishing mortality.
F-values are presented in the lower left triangle and p-values in the upper right.doi:10.1371/journal.pone.0111614.t008

Table 9 .
Mean (Million), standard deviation and CV (%) of annual post-harvest population size, N', of 100 years of simulated fishing.

Table 10 .
ANOVA-table showing effects of supervision, regulation type, target type, and target type interacting with supervision/regulation type on annual post-harvest population size.

Table 11 .
Pair-wise test of CV-values applying the F-statistics of testing differences in the variance of annual post-harvest population size.
F-values are presented in the lower left triangle and p-values in the upper right.doi:10.1371/journal.pone.0111614.t011

Table 12 .
ANOVA-table showing effects of supervision, regulation type, target type, and target type interacting with supervision/regulation type on annual SSB.

Table 13 .
Mean in thousand tons, standard deviation and CV (%) of annual SSB of 100 years of simulated fishing.The percentage of years the SSB surpassed the B trigger (B trigger = 998 thousand tons), is given in the rightmost column. doi:10.1371/journal.pone.0111614.t013

Table 14 .
Pair-wise test of CV in yield for the 100 year series, applying the F-statistics of testing differences in the variance of annual SSB. are presented in the lower left triangle and p-values in the upper right.doi:10.1371/journal.pone.0111614.t014

Table 15 .
Coefficient of variation (CV) of pre-and post-harvest population numbers compared with the unexploited population.