USP. Surplus-Production Models

Similar documents
) were both constant and we brought them from under the integral.

N t = N o e rt. Plot N t+1 vs. N t. Recruits Four Desirable Properties of S-R Relationships

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Vehicle Arrival Models : Headway

5.2. The Natural Logarithm. Solution

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Section 7.4 Modeling Changing Amplitude and Midline

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Lecture 3: Exponential Smoothing

ADVANCED MATHEMATICS FOR ECONOMICS /2013 Sheet 3: Di erential equations

An introduction to the theory of SDDP algorithm

Advanced Organic Chemistry

Linear Response Theory: The connection between QFT and experiments

Lab #2: Kinematics in 1-Dimension

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

KINEMATICS IN ONE DIMENSION

Some Basic Information about M-S-D Systems

5.1 - Logarithms and Their Properties

PET467E-Analysis of Well Pressure Tests/2008 Spring Semester/İTÜ Midterm Examination (Duration 3:00 hours) Solutions

20. Applications of the Genetic-Drift Model

Sterilization D Values

Institute, Guangzhou , China. Technology, China. Chinese Academy of Fishery Sciences, Qingdao , China. *[

1.6. Slopes of Tangents and Instantaneous Rate of Change

Kalman filtering for maximum likelihood estimation given corrupted observations.

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin

The average rate of change between two points on a function is d t

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

N H. be the number of living fish outside area H, and let C be the cumulative catch of fish. The behavior of N H

Overview. COMP14112: Artificial Intelligence Fundamentals. Lecture 0 Very Brief Overview. Structure of this course

1 Differential Equation Investigations using Customizable

INDEX. Transient analysis 1 Initial Conditions 1

Age (x) nx lx. Age (x) nx lx dx qx

Perfect selectivity vs. biomass harvesting in an age- structured fishery population model. What is the gain?

Outline. Math Elementary Differential Equations. Solution of Linear Growth and Decay Models. Example: Linear Decay Model

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Understanding the asymptotic behaviour of empirical Bayes methods

Final Spring 2007

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

4.1 - Logarithms and Their Properties

STATE-SPACE MODELLING. A mass balance across the tank gives:

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

FITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

The Simple Linear Regression Model: Reporting the Results and Choosing the Functional Form

Air Traffic Forecast Empirical Research Based on the MCMC Method

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Final Exam. Tuesday, December hours

Solutions to Odd Number Exercises in Chapter 6

not to be republished NCERT MATHEMATICAL MODELLING Appendix 2 A.2.1 Introduction A.2.2 Why Mathematical Modelling?

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

Robust estimation based on the first- and third-moment restrictions of the power transformation model

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18

Math 10B: Mock Mid II. April 13, 2016

Echocardiography Project and Finite Fourier Series

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Class Meeting # 10: Introduction to the Wave Equation

Two Coupled Oscillators / Normal Modes

Affine term structure models

Testing for a Single Factor Model in the Multivariate State Space Framework

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED

Reliability of Technical Systems

SOFTWARE RELIABILITY GROWTH MODEL WITH LOGISTIC TESTING-EFFORT FUNCTION CONSIDERING LOG-LOGISTIC TESTING-EFFORT AND IMPERFECT DEBUGGING

Failure of the work-hamiltonian connection for free energy calculations. Abstract

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model:

CHAPTER 17: DYNAMIC ECONOMETRIC MODELS: AUTOREGRESSIVE AND DISTRIBUTED-LAG MODELS

THE CATCH PROCESS (continued)

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

You must fully interpret your results. There is a relationship doesn t cut it. Use the text and, especially, the SPSS Manual for guidance.

Lecture 33: November 29

Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Method

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

Dynamic models for largedimensional. Yields on U.S. Treasury securities (3 months to 10 years) y t

18 Biological models with discrete time

6. Solve by applying the quadratic formula.

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

Properties of Autocorrelated Processes Economics 30331

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

CSE 3802 / ECE Numerical Methods in Scientific Computation. Jinbo Bi. Department of Computer Science & Engineering

Comparing Means: t-tests for One Sample & Two Related Samples

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Fishing limits and the Logistic Equation. 1

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product

1. Consider a pure-exchange economy with stochastic endowments. The state of the economy

, where P is the number of bears at time t in years. dt (a) If 0 100, lim Pt. Is the solution curve increasing or decreasing?

1. VELOCITY AND ACCELERATION

Seminar 4: Hotelling 2

FITTING EQUATIONS TO DATA

22. Inbreeding. related measures: = coefficient of kinship, a measure of relatedness of individuals of a population; panmictic index, P = 1 F;

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

5. Stochastic processes (1)

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints

Transcription:

USP Surplus-Producion Models

2 Overview Purpose of slides: Inroducion o he producion model Overview of differen mehods of fiing Go over some criique of he mehod Source: Haddon 2001, Chaper 10 Hilborn and Walers 1992, Chaper 8

3 Inroducion Nomenclaure: surplus producion models = biomass dynamic models Simples analyical models available Pool recruimen, moraliy and growh ino a single producion funcion The sock is hus an undiffereniaed mass Minimum daa requiremens: Time series of index of relaive abundance Survey index Cach per uni of effor from he commercial fisheries Time series of cach

4 General form of surplus producion Are relaed direcly o Russels mass-balance formulaion: B+1 = B + R + G - M C B+1 = B + P C B+1 = B + f(b) C B+1: Biomass in he beginning of year +1 (or end of ) B: Biomass in he beginning of year P: Surplus producion he difference beween producion (recruimen + growh) and naural moraliy f(b): Surplus producion as a funcion of biomass in he sar of he year C: Biomass (yield) caugh during year

5 Funcional forms of surplus producions Classic Schaefer (logisic) form: f B rb B K 1 The more general Pella & Tomlinson form: f B r p B 1 B K p Noe: when p=1 he wo funcional forms are he same

Sock biomass B 6 Populaion rajecory according o Schaefer model 100 90 80 70 60 50 40 30 20 10 0 K: Asympoic carrying capaciy B 1 B rb 1 B K 0 5 10 15 20 25 Time Unharvesed curve

Surplus producion 7 Surplus producion as a funcion of sock size rb B 14 K 12 1 Maximum producion 10 8 6 4 2 K /2 K 0 0 10 20 30 40 50 60 70 80 90 100 Sock biomass B

8 Reference poins Once he parameers have been esimaed, fishery performance indicaors useful o fisheries managemen can be calculaed. Biomass giving maximum susainable yield: BMSY K 2 Maximum susainable yield: MSY rk 4 Effor ha should lead o MSY: E r 2q MSY Fishing moraliy rae a MSY: FMSY r 2

USP Mehods of fiing

10 Esimaion procedures All mehods rely on he assumpion ha: Normally assume: Equilibrium mehods CPUE Never use equilibrium mehods! C qb CPUE Ofen use: U f B CPUE E Regression (linear ransformaion) mehods Compuaionally quick Ofen make odd assumpion abou error srucure Time series fiing Currenly considered o be bes mehod available

11 The equilibrium model 1 The model: B B rb B K C 1 1 The equilibrium assumpion Cach = Producion CPUE B 1 B and hus C rb 1 B K Each years cach and effor daa represen a an equilibrium siuaion where he cach is equal o he surplus producion a ha level of effor If he fishing regime is changed (effor or harves rae) he sock is assumed o move insananeously o a differen sable equilibrium biomass wih associaed surplus producion. The ime series naure of he daa is hus ignored. THIS IS PATENTLY WRONG, SO NEVER FIT THE DATA USING THE EQUILIBRIUM ASSUMPTION C E qb

12 The equilibrium model 2 Given his assumpion: I can be shown ha: B 1 B rb 1 B K Y Can be simplified o C rb 1 B K CPUE C E qb Where C a be E E MSY = a/(2b) MSY = (a/2) 2 /b 2 C ae be

13 The equilibrium model 3

Yield 14 Equilibrium model 4 20 18 16 14 12 10 8 6 4 2 0 If effor increases as he fisheries develops we may expec o ge daa like he blue rajecory Observaions and he fi The rue equilibrium 0 10 20 30 40 50 60 70 80 90 100 Effor

15 The regression model 1 Nex biomass = his biomass + surplus producion - cach 1 2 B B B 1 B rb (1 ) qe B K CPUE q U Subsiuing 2 in 1 gives: U 1 U U U r (1 ) q q q qb q EU U ru U r 1 r E q 1 r U qe U qk U qk 1 1.and dividing by Y F B qe B U q

16 The regression model 2 rae of change in B U U 1 1 inrinsic growh rae r r U qe qk densiy dependen reducion fishing his means: moraliy r M F 0 Regress: U U 1 1 on U and E which is a muliple regression of he form: Y b b X b X where X U and X E 0 1 1 2 2 1 2 r b0 r, b1, b2 q K qk b0 b b 1 2

17 Time series fiing 1 The populaion model: B B B rb Y K 1 1 The observaion model: Uˆ qb or Uˆ qb e i i The saisical model min 2 2 SS U Uˆ or SS lnu lnuˆ min 0 0

18 Time series fiing 2 The populaion model: Given some iniial guesses of he parameers r and K and an iniial saring value B0, and given he observed yield (Y) a series of expeced B s, can be produced. The observaion model: The expeced value of B s is used o produce a prediced series of Û (CPUE) by muliplying B wih a cachabiliy coefficien (q) The saisical model: The prediced series of Û are compared wih hose observed (U ) and he bes se of parameer r,k,b0 and q obained by minimizing he sums of squares

19 Example: N-Ausralian iger prawn fishery 1 40000 35000 30000 25000 20000 15000 10000 5000 0 Effor 1970 1975 1980 1985 1990 1995 0.30 0.25 0.20 0.15 0.10 0.05 0.00 CPUE 1970 1975 1980 1985 1990 1995 7000 6000 5000 4000 3000 2000 1000 0 Cach 1970 1975 1980 1985 1990 1995 Hisory of fishery: 1970-76: Effor and cach low, cpue high 1976-83: Effor increased 5-7fold and cach 5fold, cpue declined by 3/4 1985-: Effor declining, cach inermediae, cpue gradually increasing. Objecive: Fi a sock producion model o he daa o deermine sae of sock and fisheries in relaion o reference values.

20 Example: N-Ausralian iger prawn fishery 2 0.30 0.25 0.20 0.15 0.10 0.05 0.00 CPUE 1970 1975 1980 1985 1990 1995 2.0 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 B/Bmsy 1970 1975 1980 1985 1990 1995 1.2 1.0 0.8 0.6 0.4 0.2 0.0 E/Emsy 1970 1975 1980 1985 1990 1995 The observe cach raes (red) and model fi (blue). Parameers: r= 0.32, K=49275, Bo=37050 Reference poins: MSY: 3900 Emsy: 32700 Bcurren/Bmsy: 1.4 The analysis indicaes ha he curren saus of he sock is above Bmsy and ha effor is below Emsy. Quesion is how informaive are he daa and how sensiive is he fi.

21 Adding auxiliary informaion In addiion o yield and biomass indices daa may have informaion on: Sock biomass may have been un-fished a he beginning of he ime series. I.e. B0 = K. Thus possible o assume ha U0 = qb0=qk However daa mos ofen no available a he beginning of he fisheries Absolue biomass esimaes from direc couns, acousic or rawl surveys for one or more years. Those years may be used o obain esimaes of q. Prior esimaes of r, K or q q: agging sudies r: basic biology or oher similar populaions In he laer case make sure ha he esimaes are sound K: area or habia available basis

USP Some word of cauion

23 One way rip One way rip Increase in effor and decline in CPUE wih ime A lo of cach and effor series fall under his caegory. Effor CPUE Time Time Cach Time

24 One way rip 2 Fi # r K q MSY Emsy SS R 1.24 3 10 5 10 10-5 100,000 0.6 10 6 ---- T1 0.18 2 10 6 10 10-7 100,000 1.9 10 6 3.82 T2 0.15 4 10 6 5 10-7 150,000 4.5 10 6 3.83 T3 0.13 8 10 6 2 10-7 250,000 9.1 10 6 3.83 Idenical fi o he CPUE series, however: K doubles from fi 1 o 2 o 3 Thus no informaion abou K from his daa se Thus esimaes of MSY and Emsy is very unreliable

25 The imporance of perurbaion hisories 1 Principle: You can no undersand how a sock will respond o exploiaion unless he sock has been exploied. (Walers and Hilborn 1992). Ideally, o ge a good fi we need hree ypes of siuaions: Sock size Effor ge parameer low low r high ( K) low K (given we know q) high/low high q (given we know r) Due o ime series naure of sock and fishery developmen i is virually impossible o ge hree such divergen & informaive siuaions

26 The imporance of perurbaion hisories 2 db d U U 1 1 r r U qe qk r = max when B is low and E is low This means overfishing When C/E = Kq, hen r = 0 E = rq, r = 0 Biomass Effor Usual siuaion, CPUE (C/E) declines as effor (E) increases. Very lile info on r

27 One more word of cauion Mos abused sock-assessmen echnique! Published applicaions based on he assumpion of equilibrium should be ignored Problem is ha equilibrium based models almos always produce workable managemen advice. Non-equilibrium model fiing may however reveal ha here is no informaion in he daa. Laer no useful for managemen, bu i is scienific Efficiency is likely o increase wih ime, hus may have: q q q Commercial CPUE indices may no be proporional o abundance. I.e. Is he relaionship is ruly proporional? CPUE 0 incr qb

CPUE 28 CPUE = f(b): Wha is he rue relaionship?? 10 9 8 7 6 5 4 3 2 1 0 0 5 10 Biomass Hypersabiliy Proporional Hyperdepleion

29 Accuracy = Precision + Bias No accurae and no precise Accurae bu no precise (Vaguely righ) Precise bu no accurae (Precisely wrong) Accurae and precise I is beer o be vaguely righ han precisely wrong!