More viability analysis. Let's review what we've done so far (this was a bit confusing, so it's good to review): We've looked at single population of animals, and have come up with an estimate of how big a reserve should be in order to make sure that this population survives. How? Simplified, what we did is as follows: 1) Get information on rain fall 2) From this, estimate plant growth 3) Determine how this affects animal growth 4) Since #1 fluctuates considerably, it is logical that #3 will fluctuate 4a) plug random (but realistic) numbers for rainfall into the model 4b) see impact on #3 4c) Repeat 4a until animal numbers drop below a viable population level 4d) Record this number as "time to extinction" (generally in years) 5) Repeat # 4 100 (or more) times to get "average time to extinction". 6) Repeat 4 and 5 for different reserve sizes. 7) Since distribution is skewed toward shorter times, use "median time to extinction" instead. This can be done analytically (text), or simply by getting the median from the simulation 8) Once we know the "median time to extinction" for different reserve sizes, we can calculate the probability that the population will go extinct in a given number of years for different reserve sizes. 9) Fudge # 8 to take unexpected events into account. Can we generalize this procedure? In other words, for red kangaroos, we followed the above outline. But what if we're dealing with fruit flies? Bald eagles? Obviously, the above will need to be modified. Unfortunately, many generalized approaches use a single species model (only one species is modeled). As will be seen, there are some serious shortcomings with this approach. There are two methods available. We won't delve too deeply into this, since the math is very messy, and, in any case, there are no examples in the text. 1) Analytical. If one is willing to make certain assumptions (e.g., about mortality, fecundity, carrying capacity, etc.), an analytical approach can be used. Depending on the exact assumptions and how good your estimates are, the equations vary somewhat. They're presented on page 207 if you're interested, but we won't go into these. 2) Computer based. This is similar to what was done with the red kangaroos. One can either build one's own model, or use a "pre-packaged" model to do Population viability analysis (PVA). - the text mentions VORTEX which does this, and is available for free (!) on the internet (http://pw1.netcom.com/~rlacy/vortex.html). - careful - it's about 12 megs in size (I haven't had time to explore it, so I can't say too much about it). - this is still maintained, so I suspect it has changed considerably from the description in the book. - sample data sets are also available. - the basic idea is to plug in a number of parameters, and then ask for an estimate of the probability of extinction over different time intervals. - even though the program has undoubtedly changed, what the authors mention in the text is still true: - the program will only do what the operator asks. The operator is responsible for: - entering everything correctly (often we only know "guesses" as to some of these parameters - obviously this affects the output!!!) - making sure that he/she knows enough biology to understand what the program is doing - making sure the answers make sense, and even if they do, treating these very carefully (perhaps as an "initial" guess as to what might happen, hopefully to be bolstered by more research). This leads into the next topic. Assessing generalized PVA's. In particular, single species PVA's. In a single species PVA, generally the only factor considered is density. We ignore the effects of food availability, etc. - the assumption is that at higher densities, food (or whatever) becomes limiting, so "K", the carrying capacity should reflect this. Two models are in common use: - truncated exponential - logistic 1) Truncated exponential - main problem here is that both of the parameters needed have little to do with the actual ecology of the species and are hard to estimate. The text provides a few more details, but this will do. 2) logistic - this is used a lot in ecology, though often in a slightly different context (we're trying to estimate how long a population will survive). - since I think we've all had just about enough math, let's keep it simple. - the parameters are a little easier to estimate - note: just because someone puts "r" into an equation, doesn't mean it's the same "r" we've been talking about, or that we can just plug in an "r" that we've estimated. Be a little careful here. - cutting through all the "stuff", what we want to do is run a logistic growth model on our population of kangaroos, and see how it compares with the more sophisticated model described above. - We need to estimate our maximum rate of increase (this as .4 from our equation before) - We also need to estimate carrying capacity. Since this is hard to estimate (the authors mention that it appears from field data that the population has never hit it's "maximum", so we don't know what it is), it was chosen in such a way to "make the parameters align" (this isn't explained terribly well). - Essentially a density of 55 kangaroos/square kilometer was picked as an "average" carrying capacity. The actual carrying capacity varied and was adjusted once a year - (it is not quite clear from the description in the book, but it appears this was picked randomly once a year from a normal distribution with mean 55 and standard deviation 14.5) - so using the logistic equation, and proceeding with the same types of computer games as above: - tweaking r every three months, then tweaking carrying capacity once a year, getting growth information, and keeping the model going until it crashes. - then repeating this numerous times to get the same types of "average or median time of survival". - we then get the same curve as before "probability of surviving 100 years" (figure 7.10). - end result is that the logistic model implies that an area of 17 ha. is sufficient to ensure the survival of red kangaroos for 100 years. - remember, the previous estimate was 69 ha. - so which is better? - the authors argue strongly that the multi-species model is better (it models the plants as well as the kangaroos). - this also would appear to be obvious. - note also figure 7.8 - the first graph shows kangaroo growth rates as responding to density (using the 100 year simulation). - the second graph shows kangaroo growth rates as responding to biomass. - which shows a clearer trend? The bottom graph (in the top, the lines go all over the place). (the graphs might have been a little clearer if they'd just plotted points instead of lines). - another problem is that there are several different ways of applying the logistic, and every method has it's proponents. - it's a bit confusing, and I don't think conservation biologists are doing themselves a favor by these arguments. - if one person comes out with "it'll be dead in 20 years" and another "it'll be dead in 50 years", what do policy makers do? - bottom line for logistic: - if there's nothing else you can do (you don't know anything about the plants), this might be a good initial guess. Do more biology. - we said the same thing above, but now it's even more urgent. So what about PVA? - single species PVA is dangerous primarily because, in general, environmental factors are left out of the model. Both the logistic and truncated exponential consider only density and growth rates (but see below). - what about environmental fluctuations? This is only indirectly captured by carrying capacity? - but to be positive, it has allowed several first estimates for reserve sizes for such species as the Grizzly or African elephant. - presumably (this isn't explained very well) by tweaking carrying capacity, the effects of random events can be simulated on populations, and one can assess how a population might respond to certain events. - A VERY important point. - PVA assumes that the problem is low numbers. - It never explains why the species declined, or what can be done to get the numbers back up (other than "making the reserve bigger") - we should know from our discussions in here that bigger reserve size isn't always the answer (e.g., pollution, exotics, hunting, etc.). - The Mediterranean Monk seal is a great example - a study indicated that if certain population parameters are crossed, then every effort should be made to ensure juvenile survival. - but the population is already in trouble!! What happened?? - this is the biggest shortcoming of PVA - it generally does nothing about diagnosis or management considerations. - it is most useful in: - showing the vulnerability of small populations to random events - providing an initial estimate of reserve size.