Demographic PVAs

1. Demographic PVAs Simulating Demographic Stochasticity and Density Dependence

2. Demographic stochasticity • Simulated by performing so-called Monte Carlo simulations: the fate of each individual in a certain class and a certain year is decided by a set of independent random choices, all of which are based on the same set of mean vital rates • However this can greatly slow a program

3. Variability caused by demographic stocahsticity in binomial vital rates

4. Techniques to how to perform Monte Carlo simulations • For vital rates that are probabilities: Pick a uniform random number and compare its value to the probabilities of different fates an individual might experience

5. How to use a uniform random number to decide between fates in a Monte Carlo simulation 0 Survive and shrink to class 3 a34 Survive and stay to class 4 a34+a44 Survive and grow to class 5 a34+a44+a54 Survive and grow to class 6 a34+a44+a54+a64 Die 1

6. Adding demographic stochasticity to reproduction • Determine if the individual lives • Use a Poisson or another discrete distributions to obtain the individual fertility

7. Number of flowers Population 1; 2000 Sampling from Poisson distributions to estimate flower production Population 2; 2000 P (y=r) = (e –μμr) / r!

8. Number of seeds per flower Population 1; 2000 Sampling from normal distributions to estimate seed production

9. Including density dependence • Two factors make it more difficult to account for density dependence in demographic PVAs • We will rarely have as many years of vital rate estimates from a demographic study • There are many more variable that are potentially density dependent

10. Three questions we must consider: • Which vital rates are density-dependent? • How do those rates change with density? • Which classes contribute to the density that each vital rate “feels”?

11. Two more limited approaches to building density-dependent projection matrix models 1. Assume that there is a maximum number or density of individuals in one or more classes, or of the population as a whole , that the available resources can support, and construct a simulation model that prevents the population vector from exceeding this limit.

12. Two more limited approaches to building density-dependent projection matrix models 2. Choose one or at most a few vital rates that, on the basis some evidence are suspected to be strongly density dependent

13. Placing a limit on the size of one or more classes Caps or ceilings on population density are most often used to introduce density dependence when the focal species is territorial

14. The Iberian Lynx

15. Gaona et al. (1998)lynx population in Doñana National Park post breeding census birth-pulse population b = probability that a territory-holding female will breed in a given year c = number of cubs produced by females that do breed p = proportion of cubs that are female

16. Gaona et al. (1998)lynx population in Doñana National Park post breeding census birth-pulse population • Density dependence acts on g, which represent the probability that a surviving juvenile or floater will acquire a territory next year • Just before the birth pulse that precedes census t+1, there will be s4n4(t) breeding females still in possession of a territory and K- s4n4(t) vacant territories available.

17. Gaona et al. (1998)lynx population in Doñana National Park post breeding census birth-pulse population • g= [K-s4n4(t)]/[s2n2(t)+ s3n3(t)]

19. Density dependent functions The Ricker function

20. Density dependent functions Beverton-Holt model

21. Density dependent functions