Optimization model to start harvesting in stochastic aquaculture system

Establishment of cost‐effective management strategy of aquaculture is one of the most important issues in fishery science, which can be addressed with bio‐economic mathematical modeling. This paper deals with the aforementioned issue using a stochastic process model for aquacultured non‐renewable fishery resources from the viewpoint of an optimal stopping (timing) problem. The goal of operating the model is to find the optimal criteria to start harvesting the resources under stochastic environment, which turns out to be determined from the Bellman equation (BE). The BE has a separation of variables type structure and can be simplified to a reduced BE with a fewer degrees of freedom. Dependence of solutions to the original and reduced BEs on parameters and independent variables is analyzed from both analytical and numerical standpoints. Implications of the analysis results to management of aquaculture systems are presented as well. Numerical simulation focusing on aquacultured Plecoglossus altivelis in Japan validates the mathematical analysis results. Copyright © 2017 John Wiley & Sons, Ltd.     

Optimal harvesting of a Gompertz population model with a marine protected area and interval‐value biological parameters

To protect fishery populations on the verge of extinction and sustain the biodiversity of the marine ecosystem, marine protected areas (MPA) are established to provide a refuge for fishery resource. However, the influence of current harvesting policies on the MPA is still unclear, and precise information of the biological parameters has yet to be conducted. In this paper, we consider a bioeconomic Gompertz population model with interval‐value biological parameters in a 2‐patch environment: a free fishing zone (open‐access) and a protected zone (MPA) where fishing is strictly prohibited. First, the existence of the equilibrium is proved, and by virtue of Bendixson‐dulac Theorem, the global stability of the nontrivial steady state is obtained. Then, the optimal harvesting policy is established by using Pontryagin's maximum principle. Finally, the results are illustrated with the help of some numerical examples. Our results show that the current harvesting policy is advantageous to the protection efficiency of an MPA on local fish populations.     

A JUVENILE‐ADULT DISCRETE‐TIME PRODUCTION MODEL OF EXPLOITED FISHERY SYSTEMS

Abstract. Previous mathematical modeling of the population dynamics of Georges Bank Atlantic cod fishery employed discrete‐time models without age‐structure. To make use of a much wider variety of data on fisheries and fish stocks than was possible with an unstructured model, we introduce a juvenile‐adult age‐structured production exploited fishery model with a very general recruitment function. We use the age‐structured model to study the interaction between fish exploitation levels and recruitment dynamics. As case studies, we use our model results and historical fish population data from Georges Bank to investigate the impact of recent harvesting levels on the sustainability of cod fishery. We show that a constant harvesting policy with the same harvesting rate of 2007 would lead to the recovery and sustainability of Georges Bank cod fishery.     

An optimal switching approach toward cost‐effective control of a stand‐alone photovoltaic panel system under stochastic environment

The operation of a stand‐alone photovoltaic (PV) system ultimately aims for the optimization of its energy storage. We present a mathematical model for cost‐effective control of a stand‐alone system based on a PV panel equipped with an angle adjustment device. The model is based on viscosity solutions to partial differential equations, which serve as a new and mathematically rigorous tool for modeling, analyzing, and controlling PV systems. We formulate a stochastic optimal switching problem of the panel angle, which is here a binary variable to be dynamically controlled under stochastic weather condition. The stochasticity comes from cloud cover dynamics, which is modeled with a nonlinear stochastic differential equation. In finding the optimal control policy of the panel angle, switching the angle is subject to impulsive cost and reduces to solving a system of Hamilton‐Jacobi‐Bellman quasi‐variational inequalities (HJBQVIs). We show that the stochastic differential equation is well posed and that the HJBQVIs admit a unique viscosity solution. In addition, a finite‐difference scheme is proposed for the numerical discretization of HJBQVIs. A demonstrative computational example of the HJBQVIs, with emphasis on a stand‐alone experimental system, is finally presented with practical implications for its cost‐effective operation.     

HARVESTING AN AGE‐STRUCTURED POPULATION AS BIOMASS: DOES IT WORK?

Abstract The economics of fisheries is based heavily on describing fish populations by the surplus production model. Both economists and ecologists have different opinions on whether this approach provides an adequate biological basis for economic analysis. This study takes an age‐structured population model and shows how, under equilibrium conditions, it determines the surplus production model. The surplus production model is then used to solve an optimal feedback policy for a generic optimal harvesting problem. Next, it is assumed that the fishery manager applies this feedback policy even though the fish population actually evolves according to the age‐structured model. This framework is applied to the widow rockfish, Atlantic menhaden, and Pacific halibut fisheries. Population age‐structure contains information on future harvest possibilities. The surplus production model neglects this information and may lead to major deviations between the expected and actual outcomes especially under multiple steady states and nonlinearities.     

Sliding and oscillations in fisheries with on–off harvesting and different switching times

In this paper, we propose a fishery model with a discontinuous on–off harvesting policy, based on a very simple and well known rule: stop fishing when the resource is too scarce, i.e. whenever fish biomass is lower than a given threshold. The dynamics of the one-dimensional continuous time model, represented by a discontinuous piecewise-smooth ordinary differential equation, converges to the Schaefer equilibrium or to the threshold through a sliding process. We also consider the model with discrete time impulsive on–off switching that shows oscillations around the threshold value. Finally, a discrete-time version of the model is considered, where on–off harvesting switchings are decided with the same discrete time scale of non overlapping reproduction seasons of the harvested fish species. In this case the border collision bifurcations leading to the creations and destruction of periodic oscillations of the fish biomass are studied.     

OPTIMAL FISH HARVESTING FOR A POPULATION MODELED BY A NONLINEAR PARABOLIC PARTIAL DIFFERENTIAL EQUATION

As the human population continues to grow, there is a need for better management of our natural resources in order for our planet to be able to produce enough to sustain us. One important resource we must consider is marine fish populations. We use the tool of optimal control to investigate harvesting strategies for maximizing yield of a fish population in a heterogeneous, finite domain. We determine whether these solutions include no‐take marine reserves as part of the optimal solution. The fishery stock is modeled using a nonlinear, parabolic partial differential equation with logistic growth, movement by diffusion and advection, and with Robin boundary conditions. The objective for the problem is to find the harvest rate that maximizes the discounted yield. Optimal harvesting strategies are found numerically.     

Partial differential integral equation model for pricing American option under multi state regime switching with jumps

In this paper, we consider a two dimensional partial differential integral equation (PDIE) model for pricing American option. A nonlinear rationality parameter function for two asset problems is introduced to deal with the free boundary. The rationality parameter function is added in the PDIEs used for pricing American option problems under multi-state regime switching with jumps. The resulting two dimensional nonlinear system of PDIE is then numerically solved. Based on real poles rational approximation, a strongly stable highly efficient and reliable method is developed to solve such complicated systems of PIDEs. The method is build in a predictor corrector style which makes it linearly implicit, therefore, avoids solving nonlinear systems of equations at each time step in all regimes. The method is seen to maintain the stability and convergence for large jump sizes and high volatility in each regime. The impact of regime switching on option prices corresponding to different values interest rate, volatility, and rationality parameter is computed, illustrated by graphs and given in the tables. Convergence results in each regime are presented and time evolution graphs are given to show the effectiveness and reliability of the method.     

Dynamics of a fishery system in a patchy environment with nonlinear harvesting

In this paper, a stock‐effort dynamical model with two fishing zones is discussed. The nonlinear harvesting function is assumed depending upon stock size as well as fishing effort. The migration of fish is considered between two zones. The harvesting vessels also move between zones to increase their revenue. The movements of fish and fishing vessels between zones are assumed to take place at a faster time scale as compared with processes involving growth and harvesting occurring at a slow time scale. The aggregated model is obtained for total fish stock and fishing effort. This aggregated (reduced) model is analyzed analytically as well as numerically. Biological and bionomic equilibria of the system are obtained, and criteria for local stability or instability of the system are derived. The impact of levels of taxation T on the fish population and on the revenue earned by the fishery is investigated. An optimal harvesting policy is also discussed using the Pontryagin's maximum principle. The aggregated model also exhibits Hopf and transcritical bifurcation with respect to the bifurcation parameter tax T. Numerical simulations are presented to illustrate the results.     

A BIOECONOMIC APPROACH TO THE FAUSTMANN–HARTMAN MODEL: ECOLOGICAL INTERACTIONS IN MANAGED FOREST

Abstract This paper develops a bioeconomic forestry model that makes it possible to take ecosystem services that are independent of the age structure of trees into account. We derive the Faustmann–Hartman optimal harvesting strategy as a special case. The bioeconomic model is then extended to account for the fact that forest harvesting decisions impact on other ecological resources, which provide benefits for the wider community. The paper focuses on impacts associated with disturbance caused by logging operations and habitat destruction due to tree removal. This enables us to explore the interactions between forest management and the dynamics of ecological resources. The optimal rotation rule is obtained as a variation on the traditional Faustmann–Hartman equation, where an additional term captures the potential benefits derived from the growth of the ecological resource valued at its shadow price. The steady‐state solutions to the problem and sensitivity to model parameter are identified using numerical analysis.     

Optimal harvesting in a predator-prey model with Allee effect and sigmoid functional response

This work deals with the determination of the optimal harvest policy in an open access fishery in which both prey and predator species are subjected to non-selective harvesting.The model is described by autonomous ordinary differential equation systems, the functional response of the predators is Holling type III and the prey growth is affected by the Allee effect. The catch-rate functions are based on the catch per unit effort (CPUE) or Schaefer’s hypothesis.The problem of determining the optimal harvest policy is solved by using Pontryagin’s maximal principle. The problem here studied is to maximize a cost function representing the present value of a continuous time-stream of revenue of the fishery.     

Optimal and feedback strategies for managing multicohort populations

A Beverton and Holt type linear cohort dynamics model is integrated and combined with a nonlinear stock-recruitment relationship to obtain a discrete-time multicohort harvesting model. Assuming that each age class is individually controllable, it is shown, subject to certain assumptions, that the optimal harvesting strategy is to drive the population to the maximum sustainable yield solution in one time step. In most fisheries, this controllability assumption is not met and harvesting is agewise nonselective. In this case, it may be preferable to implement a harvesting policy based on suboptimal constant effort or stock level feedback strategies, rather than implement a more complicated optimal policy. This question is addressed through numerical studies on the management of an anchovy fishery.Dedicated to G. LeitmannThe author would like to thank M. Mangel, W. Reed, P. Sullivan, and G. Swartzman for commenting on a draft of this paper.     

Geometrical methods for analyzing the optimal management of tipping point dynamics

Natural resources are not infinitely resilient and should not be modeled as being such. Finitely resilient resources feature tipping points and history dependence. This paper provides a didactical discussion of mathematical methods that are needed to understand the optimal management of such resources: viscosity solutions of Hamilton–Jacobi–Bellman equations, the costate equation and the associated canonical equations, exact root counting, and geometrical methods to analyze the geometry of the invariant manifolds of the canonical equations. Recommendations for Resource Managers

• Management of natural resources has to take into account the possible breakdown of resilience and induced regime shifts.
• Depending on the characteristics of the resource and on its present and future economic importance, either for all initial states the same kind of management policy is optimal, or the type of the optimal management policy depends on the initial state.
• Modeling should reflect the finiteness of the data.

     

Optimal fishery with coastal catch

In many spatial resource models, it is assumed that an agent is able to harvest the resource over the complete spatial domain. However, agents frequently only have access to a resource at particular locations at which a moving biomass, such as fish or game, may be caught or hunted. Here, we analyze an infinite time‐horizon optimal control problem with boundary harvesting and (systems of) parabolic partial differential equations as state dynamics. We formally derive the associated canonical system, consisting of a forward–backward diffusion system with boundary controls, and numerically compute the canonical steady states and the optimal time‐dependent paths, and their dependence on parameters. We start with some one‐species fishing models, and then extend the analysis to a predator–prey model of the Lotka–Volterra type. The models are rather generic, and our methods are quite general, and thus should be applicable to large classes of structurally similar bioeconomic problems with boundary controls. Recommedations for Resource Managers

• Just like ordinary differential equation‐constrained (optimal) control problems and distributed partial differential equation (PDE) constrained control problems, boundary control problems with PDE state dynamics may be formally treated by the Pontryagin's maximum principle or canonical system formalism (state and adjoint PDEs).
• These problems may have multiple (locally) optimal solutions; a first overview of suitable choices can be obtained by identifying canonical steady states.
• The computation of canonical paths toward some optimal steady state yields temporal information about the optimal harvesting, possibly including waiting time behavior for the stock to recover from a low‐stock initial state, and nonmonotonic (in time) harvesting efforts.
• Multispecies fishery models may lead to asymmetric effects; for instance, it may be optimal to capture a predator species to protect the prey, even for high costs and low market values of the predators.

     

MULTISPECIES EXPLOITATION WITH EVOLUTIONARY SWITCHING OF HARVESTING STRATEGIES

In this paper, we propose a bioeconomic model which describes a fishery in which each of two noninteracting species is harvested by a given group of fishers during a defined time period. Then the Fishing Regulatory Authority allows each fisher to reconsider the harvesting decision at fixed (discrete) periods of time. The model derives from an Italian fisheries management experience in the Northern Adriatic Sea, where this kind of “self‐adjusting” fishing policy has been proposed to regulate harvesting of two shellfish species. The proposed dynamic model assumes the form of a hybrid system, as the natural growth functions of the two species (in continuous time) are coupled with a discrete time adaptive system that regulates how agents switch from one harvesting strategy to the other period by period according to an evolutionary mechanism based on profit comparison. In order to obtain some insights into the basic mechanisms of the system, some relevant benchmark cases are analyzed before tackling (mainly numerically) the complete hybrid model. Our results suggest that, for proper sets of parameters, this kind of myopic and adaptive self‐regulation may ensure a virtuous trade‐off between profit maximization and resource conservation, driven by cost externalities and market pressure.     

THE INCOMPLETE‐INFORMATION SPLIT‐STREAM FISH WAR: EXAMINING THE IMPLICATIONS OF COMPETING RISKS

ABSTRACT. . It is now widely recognized that climactic regime shifts, which aperiodically alter a harvested fish stock's biomass and spatial distribution, may lead to distorted fisheries management decisions which negatively impact the fishery, both biologically and economically. This is particularly true for trans‐boundary migratory stocks, where optimal management relies on coordination among independent nation‐states. Unanticipated changes in stock distribution and abundance can upset expectations of national authorities, leading them to sanction inappropriate harvesting levels by their separately managed fleets targeting the same breeding fish stock. Our theoretical studies are based on a spatially‐distributed stochastic model, which we have called the “split‐stream model,‘ where two separately managed fleets harvest simultaneously at two separate sites. Our key assumption is that competing fleet managers, when harvesting noncooperatively, hold incomplete and asymmetric private information of current stock recruitment and spatial distribution. When subsequently negotiating to coordinate their harvests, they agree that they will share their information and then bargain over partition of the gains from their cooperation. This bargaining process takes into account the fleet's relative competitive strengths, particularly due to private information asymmetries. In this present article we introduce a more complex information structure than had been assumed in our earlier work (McKelvey and Golubtsov [2002], McKelvey, Miller and Golubtsov [2003], Mckelvey et al. [2004]). Specifically, both stock‐growth and stock‐split parameters vary stochastically and asynchronously. Thus, when harvesting noncooperatively, each fleet may possess private knowledge which is unavailable to the other. We examine the interplay of the harvesting game's information structure with other fishery characteristics, such as the fleets' economics and operating characteristics and their attitudes toward risk, to determine the implications of such structure for the outcome of the harvesting game. All of these changes are made to capture new conceptual phenomena and expand the range of applicability of the model.     

Application of Exp‐function method to EW‐Burgers equation

In this article, Exp‐function method is used to obtain an exact solution of the equal‐width wave‐Burgers equation (EW‐Burgers). The method is straightforward and concise, and its applications are promising. It is shown that Exp‐function method, with the help of symbolic computation, provides a very effective and powerful mathematical tool for solving EW‐Burgers equation. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Optimal impulsive harvesting policy for single population

In this paper, we established the exploitation of impulsive harvesting single autonomous population model by Logistic equation. By some special methods, we analysis the impulsive harvesting population equation and obtain existence, the explicit expression and global attractiveness of impulsive periodic solutions for constant yield harvest and proportional harvest. Then, we choose the maximum sustainable yield as management objective, and investigate the optimal impulsive harvesting policies respectively. The optimal harvest effort that maximizes the sustainable yield per unit time, the corresponding optimal population levels are determined. At last, we point out that the continuous harvesting policy is superior to the impulsive harvesting policy, however, the latter is more beneficial in realistic operation.     

Comparison between homotopy analysis method and optimal homotopy asymptotic method for nth‐order integro‐differential equation

This paper presents general framework for solving the nth‐order integro‐differential equation using homotopy analysis method (HAM) and optimal homotopy asymptotic method (OHAM). OHAM is parameter free and can provide better accuracy over the HAM at the same order of approximation. Furthermore, in OHAM the convergence region can be easily adjusted and controlled. Comparison, via two examples, between our solution using HAM and OHAM and the exact solution shows that the HAM and the OHAM are effective and accurate in solving the nth‐order integro‐differential equation. Copyright © 2011 John Wiley & Sons, Ltd.     

基于Leslie-Gower捕食-被食模型的渔业收获问题

考虑了基于Leslie-Gower模型的生物经济学捕获问题,通过对税收政策的控制影响渔业生态系统,研究了系统的动力学行为,并通过Pontryagins最大值原理考虑了最优的税收政策,最后给出了系统仿真.