OPTIMAL HARVESTING AND SPATIAL PATTERNS IN A SEMIARID VEGETATION SYSTEM

We consider an infinite time horizon spatially distributed optimal harvesting problem for a vegetation and soil water reaction diffusion system, with rainfall as the main external parameter. By Pontryagin's maximum principle we derive the associated four component canonical system, and numerically analyze this and hence the optimal control problem in two steps. First we numerically compute a rather rich bifurcation structure of flat (spatially homogeneous) and patterned canonical steady states (FCSS and PCSS, respectively), in 1D and 2D. Then we compute time dependent solutions of the canonical system that connect to some FCSS or PCSS. The method is efficient in dealing with non-unique canonical steady states, and thus also with multiple local maxima of the objective function. It turns out that over wide parameter regimes the FCSS, i.e., spatially uniform harvesting, are not optimal. Instead, controlling the system to a PCSS yields a higher profit. Moreover, compared to (a simple model of) private optimization, the social control gives a higher yield, and vegetation survives for much lower rainfall. In addition, the computation of the optimal (social) control gives an optimal tax to incorporate into the private optimization.

[1]  Hannes Uecker,et al.  Numerical Results for Snaking of Patterns over Patterns in Some 2D Selkov-Schnakenberg Reaction-Diffusion Systems , 2013, SIAM J. Appl. Dyn. Syst..

[2]  Florian Wagener,et al.  Bifurcations of optimal vector fields in the shallow lake model , 2010 .

[3]  Jean-Pierre Raymond,et al.  Pontryagin's Principle for Time-Optimal Problems , 1999 .

[4]  Dieter Grass,et al.  Optimal management and spatial patterns in a distributed shallow lake model , 2015, 1503.05438.

[5]  From 0D to 1D spatial models using OCMat , 2015, 1505.03956.

[6]  Paul Nicholas,et al.  Pattern in(formation) , 2012 .

[7]  Jean-Pierre Raymond,et al.  Hamiltonian Pontryagin's Principles for Control Problems Governed by Semilinear Parabolic Equations , 1999 .

[8]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[9]  N. Rashevsky,et al.  Mathematical biology , 1961, Connecticut medicine.

[10]  Anastasios Xepapadeas,et al.  Diffusion-Induced Instability and Pattern Formation in Infinite Horizon Recursive Optimal Control , 2006 .

[11]  H. Prins,et al.  VEGETATION PATTERN FORMATION IN SEMI-ARID GRAZING SYSTEMS , 2001 .

[12]  Wolf-Jürgen Beyn,et al.  Dynamic optimization and Skiba sets in economic examples , 2001 .

[13]  J. A. Salvato.  John wiley & sons. , 1994, Environmental science & technology.

[14]  Carmen Camacho,et al.  Land use dynamics and the environment , 2014 .

[15]  Suzanne Lenhart,et al.  OPTIMAL HARVESTING OF A SPATIALLY EXPLICIT FISHERY MODEL , 2008 .

[16]  Thomas Slawig,et al.  A Smooth Regularization of the Projection Formula for Constrained Parabolic Optimal Control Problems , 2011 .

[17]  M. Neubert,et al.  Marine reserves and optimal harvesting , 2003 .

[18]  W. Brock,et al.  Energy balance climate models and general equilibrium optimal mitigation policies , 2013 .

[20]  A. Skiba,et al.  Optimal Growth with a Convex-Concave Production Function , 1978 .

[21]  Anastasios Xepapadeas,et al.  Pattern Formation, Spatial Externalities and Regulation in Coupled Economic-Ecological Systems , 2008 .

[22]  Yulong Xing,et al.  OPTIMAL FISH HARVESTING FOR A POPULATION MODELED BY A NONLINEAR PARABOLIC PARTIAL DIFFERENTIAL EQUATION , 2016 .

[23]  John T. Workman,et al.  Optimal Control Applied to Biological Models , 2007 .

[24]  Donato Trigiante,et al.  A Hybrid Mesh Selection Strategy Based on Conditioning for Boundary Value ODE Problems , 2004, Numerical Algorithms.

[25]  Nicolas Barbier,et al.  The global biogeography of semi‐arid periodic vegetation patterns , 2008 .

[26]  S. Carpenter,et al.  Early-warning signals for critical transitions , 2009, Nature.

[27]  Radu Strugariu,et al.  An optimal control problem for a two-prey and one-predator model with diffusion , 2014, Comput. Math. Appl..

[28]  Yuval R. Zelnik,et al.  Regime shifts in models of dryland vegetation. , 2013, Philosophical transactions. Series A, Mathematical, physical, and engineering sciences.

[29]  Herb Kunze,et al.  Optimal control and long-run dynamics for a spatial economic growth model with physical capital accumulation and pollution diffusion , 2013, Appl. Math. Lett..

[30]  M. Rietkerk,et al.  Self-Organized Patchiness and Catastrophic Shifts in Ecosystems , 2004, Science.

[31]  Sebastian Aniţa,et al.  An Introduction to Optimal Control Problems in Life Sciences and Economics: From Mathematical Models to Numerical Simulation with MATLAB® , 2010 .

[32]  M. Silber,et al.  Transitions between patterned states in vegetation models for semiarid ecosystems. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[33]  M. Pachter,et al.  Optimal control of partial differential equations , 1980 .

[34]  Suzanne Lenhart,et al.  OPTIMAL CONTROL OF HARVESTING IN A PARABOLIC SYSTEM MODELING TWO SUBPOPULATIONS , 2001 .

[35]  Daniel Wetzel,et al.  pde2path - A Matlab package for continuation and bifurcation in 2D elliptic systems , 2012, 1208.3112.

[36]  Anastasios Xepapadeas,et al.  The spatial dimension in environmental and resource economics , 2010, Environment and Development Economics.