# (p.240) Appendix C

# (p.240) Appendix C

# Predicting the SAR: An alternative approach

Here we start with the METE results for the spatial abundance distribution at scale *A* _{0}, and the spatial abundance distribution *Π*(*n*|*A*,*n*,*A* _{0}), and simply substitute into Eq. 3.13 to get the desired answer at any scale:

The subscript 0 on *β* is to remind you that it is evaluated at scale *A* _{0}, and *β* _{0} is given by the solution to Eq. 7.27 as a function of the values of *S*̄ (*A* _{0}) and *N*(*A* _{0}). *Π*(0|*A*,*n*,*A* _{0}) in Eq. C.1 is determined from Eqs 7.48–7.50.

The predicted SAR can then be determined by numerical evaluation of Eq. C.1. For the case in which down-scaling of species' richness to determine *S*̄(*A*) is the goal, and the starting point is knowledge of *S*(*A* _{0}) and *N*(*A* _{0}), the procedure is straightforward. First solve for *β* _{0} using Eq. 7.27 (or the simpler Eq. 7.30, if *β* _{0} *N*(*A* _{0}) ≫ 1 and *β* _{0} ≪ 1), solve for *Π*(0) using Eqs 7.48–7.50, and then substitute that information in to Eq. 7.48 to determine *S*̄(*A*). The summation in Eq. C.1 can be carried out with many software packages, including Excel, MathCad, Matlab, and Mathematica.

To use Eq. C.1 to up-scale species' richness to scale *A* _{0}, starting with an empirical estimate of *S*̄(*A*), the procedure is more complicated because Eq. C.1 contains two unknowns, *β* _{0} and *S*̄(*A* _{0}), and so does Eq. 7.27. Now, just as with method 1, the two equations have to be solved simultaneously, rather than sequentially.

It might be objected that Eq. C.1 really contains three unknowns, *β* _{0}, *S*̄(*A* _{0}), and *N*(*A* _{0}), with only two equations (C.1 and 7.27) to determine them. As before, however, up-scaling abundance *N* from scale *A* to *A* _{0} is not a problem. Because we are examining the complete nested SAR, an average of the measured *N*-values in each of the plots of area *A* is an estimate of the average density of individuals in *A* _{0}. Hence, we can assume:

# (p.241) C.1 A useful approximation for method 2

While Eqs 7.27, 7.48–7.50, and C.1 can be solved numerically to either up-scale or down-scale species' richness, it is useful to examine some simplifying approximations that will often be applicable to real datasets. With these simplifications we can obtain analytically tractable solutions that do not require numerical solutions and that provide considerable insight.

We assume that *A* ≪ *A* _{0}, so that 1 − *Π*(0) can be approximated by Eq. 7.54. In practice, this means that we are considering ratios *A* _{0}/*A* > 32. Particularly for up-scaling applications, this assumption is not very constraining because we are often interested in up-scaling to areas much larger than the area of censused plots. As always, we assume *S* _{0}≫ 1 so that exp(–*β* _{0} *N*(*A* _{0})) ≪ 1.

With *A* « *A* _{0}, we can replace Eq. C.1 with:

Now wo cases separately: *β* _{0} *A* _{0}/*A* ≪ 1 and *β* _{0} *A* _{0}/*A* ~ 1 (we will show that the case *β* _{0} *A* _{0}/*A* ≫ 1 is unlikely to arise).

Box C.1 shows that if *β* _{0} *A* _{0}/*A* ≪ 1, the summation in Eq. 7.63 results in:

where *γ* is Euler's constant, 0.5772.…

If *β* _{0} *A* _{0}/*A* is ~ 1, then a more complicated result, also shown in Box C.1, is obtained.

The SAR given in Eq. C.4 can be applied either to upscaling or down-scaling species' richness. As before, for down-scaling, Eq. 7.27 is used to determine *β* _{0}, which is then substituted into Eq. C.4 to determine *S*̄(*A*). For up-scaling, Eqs C.4 and 7.27 must be solved simultaneously. We see from the form of Eq. C.4 that for *A* _{0} ≫ *A* and *β* _{0} *A* _{0}/*A* ≪ 1, the SAR is of the form:

where:

From Eqs C.11–C.13, we can extract further information about the shape of the species–area relationship. For the case of down-scaling, *β* _{0} is a fixed parameter determined using Eq.7.27 from the values of *S*̄(*A* _{0}) and *N* _{0} at scale *A* _{0}. Hence, *S*̄(*A*) depends logarithmically on area.

For up-scaling species' richness, we assume *S*̄(*A*) is known and express *S*̄(*A* _{0}), with *A* _{0} ≫ *A*, as:

Using Eq. 7.27, this can be rewritten as:

Rearranging terms, we derive:

At first glance, Eq. C.16 does not seem very helpful because the coefficient *β* _{0} in front of the log(*A* _{0}) term itself depends on *S*̄(*A* _{0}). However, numerical evaluation of the up-scaling SAR from simultaneous solution of Eqs 7.27 and C.14 indicates that *β* _{0} *N*(*A* _{0}) is nearly exactly constant over spatial scales for which our assumption *A* _{0}/*A* ≫ 1 holds. In particular, with each doubling of area, *A* _{0} → 2*A* _{0}, *N* doubles and *β* is approximately halved. At anchor scale, *β*(*A*)*N*(*A*) is only slightly smaller than *β* _{0} *N*(*A* _{0}) because while *N* exactly doubles with each area doubling, *β* decreases by slightly less than a factor of two for areas *A* _{0} only a little bigger than *A*. The bottom line is that the coefficient *β* _{0} *N*(*A* _{0}) in front of the log(*A* _{0}/*A*) term in Eq. C.16 is approximated by the value of *βN* at anchor scale. Numerical solutions to the method 2 up-scaling equations (7.27 and C.4) for a variety of anchor scale boundary conditions confirm the accuracy ot the approximations that lead to Eq. C.16.