Using the Species – Area Relationship to Set Baseline Targets for Conservation

This paper demonstrates how the power form of the Species–Area Relationship (SAR) can be used to set conservation targets for land classes using biodiversity survey data. The log-transformation of the power model is a straight line; therefore, if one knows the average number of species recorded per survey site and can estimate the true species number present in the land class, using EstimateS software, it is possible to calculate the slope of the curve, or z-value. The z-value is the exponent in the power model and it can then be used to estimate the proportion of area required to represent a given proportion of species present in any land class. This application of the SAR is explored using phytosociological relevé data from South Africa’s Succulent Karoo biome. We also provide suggestions for extrapolating the estimated z-values to other land classes within a bioregion that lack sufficient survey data, using the relationship between z-values and remotely determined landscape variables such as habitat diversity (topographic diversity) and geographic location (latitude and longitude). The SAR predicts that for most Succulent Karoo vegetation types a conservation target of 10% of the land area would not be sufficient to conserve the majority of species. We also demonstrate that not all land classes are equal from a plant biodiversity perspective, so applying one target to all land classes in a region will lead to significant gaps and inefficiencies in any reserve network based on this universal target.


INTRODUCTION
Targets are an integral part of contemporary conservation planning, implementation, and monitoring.Systematic conservation planning is dependent on explicitness, accountability, and defensibility in identifying priority conservation areas (Margules and Pressey 2000).Conservation targets underpin this process by providing a clear purpose for conservation decisions, lending them accountability and defensibility (Pressey et al. 2003).Targets are basically quantitative interpretations of broad conservation goals that are established in policy by experts, implementing agencies, or other stakeholders (Cowling et al. 1999b, Margules and Pressey 2000, Pressey et al. 2003).For example, an agency may specify that it wishes to conserve at least 10% of each vegetation type and three populations of endangered species within its jurisdiction.Consequently, targets also provide a benchmark against which to measure the success of conservation action.
Conservation targets can be divided into two broad categories based on the scale of biodiversity surrogate targeted (Noss 1996, Pressey et al. 2003).Coarse-filter approaches set targets for features such as vegetation types, ecosystems, or land classes.Fine-filter approaches use species or populations as the focal feature for conservation action (Noss 1996).Although both approaches are complementary, for most regions limitations in species distribution data sets obligate the use of coarse-filter surrogates (Lombard et al. 2003;P. Desmet, A. Lombard, and A. Cowling, unpublished manuscript).
Vegetation or land-class maps have the advantage of covering the entire landscape, thereby eliminating the inherent spatial and taxonomic bias of species data sets (Lombard et al. 2003, Pressey et al. 2003;P. Desmet, A. Lombard, and A. Cowling, unpublished manuscript).There are limitations with using such maps.Firstly, reserve selection using the coarse-filter approach is likely to protect many species yet to be discovered or for which records are deficient.However, unless complementary Ecology and Society 9(2): 11. http://www.ecologyandsociety.org/vol9/iss2/art11fine-filter information is incorporated in the process, other species, especially rarer ones, are likely to be missed (Kirkpatrick andBrown 1994, Lombard et al. 2003;P. Desmet, A. Lombard, and A. Cowling, unpublished manuscript).Secondly, the spatial, landclass compositional, or process requirements of certain species are unlikely to be satisfied unless they are specifically targeted (Pressey et al. 2003).Such taxa have been referred to as "focal species" (Lambeck 1997, Boshoff and Kerley 1999, Noss et al. 1999, Boshoff et al. 2001) or "landscape species" (Sanderson et al. 2002).These are essentially an alternative approach to interpreting the "umbrella" or "keystone" species concepts (Bond 1993, Mills et al. 1993, Launer et al. 1994, Fleisman et al. 2000).Thirdly, land classes do not explicitly target natural processes (Pressey et al. 2003).These need to be targeted if biodiversity is to persist (Cowling et al. 1999b, Cowling andPressey 2001).
These problems are compounded by problems related to scale.Targets framed as percentages of countries or regions can be achieved even though they fail to protect the natural features most urgently in need of protection (Pressey et al. 2003).Large regions are heterogeneous in terms of biodiversity and potential for anthropogenic transformation.Conservation areas have often been relegated to the least useable portions of regions, thereby avoiding areas where past impacts on biodiversity have been greatest and future threats are most serious (Noss et al. 1999, Pressey et al. 2000, Scott et al. 2001).This is true even of regions with formal protection overall equal to or greater than 10% (Armesto et al. 1998, Barnard et al. 1998, Soulé and Sanjayan 1998, Pressey et al. 2000, Rouget et al. 2003).Also, coarse-scale maps do not capture all possible land-class combinations, so even if vulnerability over the whole area is low, certain landscape biodiversity features (e.g., rare vegetation types or habitats) can fall through the cracks (P.Desmet, A. Lombard, and A. Cowling, unpublished manuscript).These issues can be addressed by mapping at finer scales and with improved mapping techniques (Ferrier 2002;P. Desmet, A. Lombard, and A. Cowling, unpublished manuscript).Targeting better mapped land types with classes that are more homogeneous, in terms of biodiversity and land use potential, limits the potential for conservation action to miss capturing all biodiversity (Bedward et al. 1992).
Considering these limitations, it must be accepted that, for the majority of areas on this planet, land-class maps of some sort will be the primary biodiversity feature used for conservation and land-use planning.Thus, how can biologically meaningful quantitative conservation targets for land classes be set?Although there are some studies dealing with a range of species (Margules et al. 1988, Saetersdal et al. 1993, Travaini et al. 1997), minimum viable population (Nunney and Campbell 1993, Boshoff et al. 2001, Burgman et al. 2001), meta population (Lindenmayer and Lacy 1995), genetic diversity (Lacy 1997, Ferguson et al. 1998), community (Prins et al. 1998), habitat (Turner et al. 1999, Calkin et al. 2002), or ecosystem (Turner et al. 1992, Noss 1996) targets, there is generally a paucity of work dealing specifically with targets for land classes.
The widely used 10% target, recommended by IUCN (International Union for the Conservation of Nature), when applied to land classes, implies that all are equal in terms of their species diversity, abundance, and distribution.This is certainly not the case.More questionable though, is the biological foundation for this target.Despite the potential arbitrariness of this target (Soulé andSanjayan 1998, Pressey et al. 2003), the origin is partly founded in the original work on the species-area relationship.A general "rule of thumb" first noted by Darlington (1957) and developed from early observations of species-area relationships was that a 10fold decrease in area resulted in a twofold decrease in species, or, alternatively, 10% of an area would conserve 50% of the species (Diamond and May 1976).Even if at some point this observation influenced the conception of the IUCN 10% target at the World Parks Congress in Caracas in 1993, the question arises as to whether saving 50% of the planet's terrestrial species is really adequate.General consensus in the literature is that 10% of area is not sufficient to represent the majority of biodiversity, assuming that the remainder of the landscape is cleared or not conserved (Soulé andSanjayan 1998, Rodrigues andGaston 2001).Is it possible to set ecologically meaningful targets for land classes?In an attempt to address this question, this paper returns to one of ecology's oldest observations, the species-area relationship (SAR), to find answers."You will find more species if you sample a larger area" (Rosenzweig 1995).This could also be stated as follows: "you will conserve more species if you conserve a larger area."This observation forms the basis of the SAR.Patterns in the SAR are well explored in the ecological literature across a range of spatial and temporal scales (Rosenzweig 1995, Lomolino 2000, Knowles 2001, Lomolino and Weiser 2001, Collins et al. 2002, Haila 2002).Attempts have http://www.ecologyandsociety.org/vol9/iss2/art11been made to develop a functional understanding of the SAR (e.g., Harte et al. 1999b).Although there is as yet no widely accepted ecological theory to explain the relationship, the basic pattern is real.Many mathematical models can and have been used to describe the SAR (Lomolino 2000), but the power relationship or Arrhenius equation is probably the most popular in the literature (Rosenzweig 1995): (1) In this relationship, k is a scaling factor that relates to the sample size used (Rosenzweig 1995, Lomolino 2000) and it can be ignored when comparing proportions or percentages of area and species (Vreugdenhil et al. 2003).The meaning of k is debatable (Gould 1979, Lomolino 2000, 2001), but we will not discuss this here.The rate at which species are encountered in a system is described by the parameter z (Rosenzweig 1995, Lomolino 2000).The SAR has been applied to a host of questions in ecology.Those with a direct conservation angle range from the "Single Large Or Several Small" debate (Diamond and May 1976, Bond et al. 1988, Deshaye and Morisse 1989, Cowling and Bond 1991) to predicting the loss of species from fragmented landscapes (Tilman et al. 1994, Kemper et al. 1999, Cornelius et al. 2000, Miller and Cale 2000, Magura et al. 2001, Acosta and Robertson 2002, Haila 2002, Tscharntke et al. 2002, Zurlini et al. 2002).
Fragmentation research has come very close to actually addressing the target issue by using Eq. 1 to estimate the proportion of species lost, given a reduction in area (Brooks et al. 1999, Kinzig andHarte 2000).However, the suitability of using the SAR to predict species loss due to habitat transformation has been questioned recently (Seabloom et al. 2002).From a conservation planning perspective, rather than asking how many species will be lost if a landscape is fragmented, the question can be turned around to ask how many species will be gained by a protected area network if more landscape is added to the network.Diamond and May (1976) used the power equation to predict the number of species remaining if a given percentage of a landscape was transformed.If the zvalue for a biota or land class is known, then by using Eq. 2 it is possible to predict the proportion of species remaining if the area is reduced by a given proportion. (2) Here S' and A' denote the proportion of species and area rather than absolute values.This equation can be reordered to address conservation targets to determine the proportion of area required to represent a given percentage of species: (3) For the SAR rule of thumb previously discussed (a 10fold decrease in area equates to a twofold decrease in species) the z-value here is approximately 0.3.Using this same z-value but changing the percentage of species targeted to 75%, it follows from Eq. 3 that the percentage of the area required increases from 10% to 38%.If the species target is increased to 95%, then the area target becomes 84%!These are quite significant changes in the area required to meet this basic biodiversity target of representing each species at least once.Published z-values for biotas range between approximately 0.1 and 0.4 (Rosenzweig 1995).Although this range in the exponent is small, the nature of the power equation means that for a species target of 75%, the area target can range from 5% to 48%, respectively, for z-values from 0.1 to 0.4.
It is possible to use the SAR to set conservation targets for land classes.The method for setting this target involves estimating the area of a land class that is required to represent a given proportion of the species occurring in the land class.Thus, it remains a question of being able to calculate the z-value of the SAR for a land class.To achieve this it is necessary to generate a SAR curve based on some form of inventory data.Because the curve of the power model is a straight line in log-log space and the slope is the z-value, it is possible to calculate the z-value without generating the actual curve by using species-area data that sample larger and larger proportions of a land class.There is no need to demonstrate the relationship by fitting the power curve, because the suitability of the model is assumed a priori.For the log transformation of the power model, the slope of the curve (hence, the zvalue) can be determined using the formula for http://www.ecologyandsociety.org/vol9/iss2/art11calculating the slope of a straight line: (4) Here z is the slope of the straight line, y 2 = log(total number of species in a land class); y 1 = log(average number of species per survey sample); x 2 = log(total area of land class); and, x 1 = log(average area of samples).When using inventory data, three of these variables are known; all that remains is to estimate the total number of species that occur in the vegetation type.Several nonparametric estimator functions are available for estimating the true number of species in an area based on a set of random samples (Colwell andCoddington 1995, Colwell 1997).
In this paper, phytosociological survey data are used to calculate the z-values for land classes (Succulent Karoo biome vegetation types) by estimating the true number of species per vegetation type using EstimateS software (Colwell 1997).Firstly, z-values for vegetation types with sufficient survey sites to allow for meaningful estimates to be made are calculated.Following from this, observed z-values are extrapolated to other vegetation types in the Succulent Karoo that do not have survey data.This is achieved by relating known z-values to landscape physical properties that act as surrogates for geographic species turnover and habitat diversity.Finally, some generalizations about land-class characteristics are made that can be applied to other conservation situations in which suitable survey data are lacking.

Study area
The Succulent Karoo biome covers an area of approximately 116 000 km 2 located in the southwestern corner of southern Africa (Fig. 1).It forms the arid ecotone (rainfall <350 mm/yr) between the winter-rainfall Fynbos biome and the summerrainfall biomes of the region (Cowling et al. 1999a, Dean andMilton 1999).Floristically, it forms part of the Cape Floral Kingdom (J.Born, H. P. Linder, and P. Desmet, unpublished data).The hallmark of the Succulent Karoo is its exceptionally diverse and endemic-rich flora, especially succulents and bulbs.The biome is home to 6356 plant species, 40% of which are endemic and 936 (17%) of which are Red Data Listed (Driver et al. 2003).In addition to its floral diversity, the biome is a center of diversity for reptiles and many groups of invertebrates.The Succulent Karoo biome is one of 25 internationally recognized biodiversity hotspots and is the world's only arid hotspot (Mittermeier et al. 1998).However, only 3.5% is formally conserved, and the biome's biodiversity is under pressure from a range of human impacts, especially mining, crop agriculture, ostrich farming, overgrazing, illegal collection of fauna and flora, and anthropogenic climate change (Driver et al. 2003).Hence the urgency to address the question of conservation targets.

Estimating z-values for selected Succulent Karoo vegetation types
Phytosociological relevé data from the Succulent Karoo biome were used to estimate true species numbers for Succulent Karoo vegetation types.The relevé database contained 5491 geo-referenced survey sites from 25 different studies conducted in the biome over the last 20 years.The Succulent Karoo was expertly classified and mapped into 132 vegetation types (Driver et al. 2003).Each relevé was assigned to a vegetation type based on its geographic location.In total, 42 vegetation types, each with more than 30 relevés, were used for the analyses.The 30-sample limit was determined in exploratory analyses as being the smallest average number of samples for which the standard deviation of the final estimate was less than 5% of the estimate.
The software package EstimateS (Colwell 1997) was used to generate estimates of the total number of species present in each vegetation type based on the samples of survey plots.EstimateS provides the user with 10 nonparametric statistical functions to estimate the true number of species based on species incidence in a sample of survey sites.These estimators are discussed in Colwell and Coddington (1995).Only seven estimators (ACE, ICE, Chao1, Chao2, Jack-knife1, Jack-knife2, and Bootstrap) were used in the comparative analyses, as the others (Michaelis-Menten runs, Michaelis-Menten means, and Cole) have not been widely tested for this purpose in the literature (Brose 2002, Herzog et al. 2002, Chiarucci et al. 2003, Petersen et al. 2003, Turnipseed et al. 2003).The estimated true species number was used in Eq. 4 to calculate the z-value of the SAR for each vegetation type.http://www.ecologyandsociety.org/vol9/iss2/art11

Extrapolating z-values to all Succulent Karoo vegetation types
The effect of area has already been accounted for in the z-value.Species diversity also relates strongly to habitat diversity (Rosenzweig 1995).Therefore, it makes sense to introduce a variable that can act as a surrogate for habitat diversity.Z-values are related to two independent measures of land-class habitat diversity.The first is a simplistic measure of topographic diversity summarized in two ways: (1) the standard deviation of the mean altitude for a vegetation type as determined from a 100-m grid cell resolution digital elevation model; and (2) the ratio of volumetric to planimetric surface area per vegetation type, determined from the same elevation model.The hypothesis is that more topographically heterogeneous land classes would have more habitats and, consequently, would support more species per unit area and have higher z-values.
The ultimate goal of this exercise is to generate conservation targets for all land classes within a biome or planning domain.The problem faced in the Succulent Karoo, as well as most of the rest of the world, is that there are inadequate survey data for all land classes.For this Succulent Karoo example, 42 out of 132 vegetation types had 30 or more survey sites, with only nine having more than 100 relevés.Therefore it is necessary to extrapolate z-values (or targets) to all vegetation types within a biogeographic province based on an observed relationship between the calculated z-values for vegetation types or land classes and some remotely measurable land-class properties that could explain macro patterns of diversity.
The second measure of habitat diversity is a count of the number of ecological land classes present within each http://www.ecologyandsociety.org/vol9/iss2/art11vegetation type.Eecological zones were developed using a generalized dissimilarity modeling technique discussed in Ferrier ( 2002) and applied by S. Ferrier, G. Manion, and P. Desmet (unpublished manuscript) for the Namaqualand region of the Succulent Karoo.Each zone is determined as a function of remotely determined topographic, edaphic, and climatic variables scaled according to observed patterns of plant diversity.Each zone can be considered environmentally and biologically homogeneous relative to other zones; thus a count of the number of zones represented in a vegetation type can be regarded as an alternative, albeit crude, measure of habitat diversity.
In addition to the three habitat or beta diversity variables, latitude and longitude were introduced as geographic or gamma diversity variables.For the analyses, the geographic centroid of each vegetation type was used.These variables are useful at the landscape scale, because patterns of gamma diversity relate strongly to distance between areas (Ferrier et al. 2002).
To develop the relationships between z-values and the landscape variables, the 42 vegetation types were grouped into eight higher order vegetation-type categories, i.e., categories between the lower level of vegetation type and upper level of biome (Table 1).The groups were based on known biogeographic and physiognomic similarities in vegetation types.Regression models relating z-values to the environmental variables were built using linear, nonlinear, and generalized additive modeling methods.The S-Plus statistical package (http://www.mathsoft.com/splus/)was used to perform analyses.

Estimating z-values for vegetation types with sufficient survey data
Estimates of true species number for Succulent Karoo vegetation types are presented in Table 2.As expected, there is considerable variation in the predictions of the different estimators.For this analysis, we decided to use only the Bootstrap estimator to calculate the zvalue rather than, for example, using a mean of all estimators.Exploratory analyses using model data indicated that this estimator provided the most consistent estimates under a range of modeled data conditions (P.Desmet, unpublished data).The Bootstrap estimator is consistently the lowest of the seven estimates of true species number (Table 2).This estimate is, for the most part, within or close to the lower standard deviation limit of the other estimates.The z-value was calculated for each vegetation type using Eq. 4. The average relevé size was taken to be 10 x 10 m or 100 m 2 .Varying the relevé area by an order of magnitude (i.e., 10 to 1000 m 2 ) changed the resultant z-value by less than 0.01%.

Extrapolating z-values to vegetation types with insufficient data
In order to extrapolate the observed patterns in zvalues to all vegetation types in the Succulent Karoo, the z-values were related to the five landscape environmental variables (Fig. 2).In Fig. 2, the z-values are classified into the higher order vegetation categories presented in Table 1.There appear to be relationships between observed z-values per vegetation type category and the landscape variables examined.
Regression models were built for the Namaqualand Succulent Karoo (Fig. 3); Namaqualand quartz patches (Fig. 4 and Fig. 5); and Sandveld (Fig. 6) vegetation categories (see Appendix 1).These were the only models that were significant, given the number of data points (i.e., vegetation types with calculated z-values) available fore each model.For all three cases, the basic model used to extrapolate values expresses the z-value, firstly, as a function of geographic location and, secondly, as a function of topographic diversity.
Comparing the response of z-values in the three models, it is clear that patterns of plant biodiversity within different higher order vegetation-type categories respond differently in relation to the landscape variables examined.Although geographic distance and habitat diversity are significant variables explaining the pattern in z-values, the observed pattern is different in each of the three vegetation groups.Thus, it was not possible to develop a single general model to extrapolate z-values for all Succulent Karoo vegetation types.For the remainder of vegetation types that did not fall into one of the three vegetation categories with models, we awarded the observed zvalue for the geographically nearest vegetation type within their respective higher order vegetation category.

Calculating conservation targets from z-values
Conservation targets were calculated for vegetation types using Eq. 3. Examples of the range in conservation target values, expressed as the percentage of vegetation type required to represent a given percentage of plant species occurring in that vegetation type, are calculated for the observed range in z-values (,Table 3).
If the conservation objective is to represent the majority of biodiversity (e.g., between 70% and 80% of species) within the formal reserve network, then the SAR would predict that well in excess of 10% would be required for most land classes.For the Succulent Karoo vegetation types, where a conservative average estimate of the z-value is 0.18, this would translate into a target of between 14% and 30% of the land area required to represent between 70% and 80% of the species, respectively.For the Succulent Karoo, a representation target of 75% of species was set based on regional-level discussion with stakeholders.The complete set of targets for Succulent Karoo vegetation types can be found in Driver et al. (2003).

Is 10% enough?
The approach to setting targets discussed in this paper does provide an ecological framework for testing the validity of the widely used 10% conservation target.
The SAR would predict that, for most Succulent Karoo vegetation types, a conservation target of 10% of the land area would not be sufficient to conserve the majority of species.It is likely that this is the same for most land classes elsewhere on Earth.A 10% target may be valid for only the most species-poor land classes.
Another important finding of this study is that not all land classes are equal from a biodiversity perspective.Just as 10% is not enough to represent most species, applying one target to all land classes will lead to significant gaps and inefficiencies in any resultant reserve network (Soulé andSanjayan 1998, Pressey et al. 2003).

Assumptions of the power model
The approach here assumes that the power model best describes the species-area relationship.This model rests on the assumption that species-abundance within in a land class follows a log-normal distribution (Preston 1948, May 1975, Rosenzweig 1995).It is questionable whether this assumption holds for most species.It is to be expected that as species-abundance distributions deviate from this distribution, so the difference between predicted and actual targets widens.Harte et. al. (1999a) have proposed an alternative derivation of the power model based on self-similarity in the distribution and abundance of species.This derivation does not assume a log-normal species-abundance distribution.Thus, the power model may still be valid for setting conservation targets if it is not dependent on the species-abundance distribution.
In addition, the nature of the power function means that the curve only reaches 100% of species at 100% of area.How valid this is in reality is debateable (Lomolino 2002).It is probably possible to represent 100% of species within less than 100% of area.Does this represent a breakdown in the validity of the power model as the curve approaches the asymptote, or a defect in the model as a whole?This artefact should not stop conservationists from targeting 100% of species.However, the power form of the SAR cannot be used to predict what the actual area will be for achieving this target.Ecology and Society 9(2): 11. http://www.ecologyandsociety.org/vol9/iss2/art11

Limitations of using z-values to set targets
The most important limitation of using z-values to set conservation targets is that it says nothing about where species are located in the landscape.It only provides an indication of the rate at which species are likely to be accumulated.Consequently, it says nothing about which 20% of the land class is required to represent the 75% of species being targeted.If species are distributed randomly in a land class, then reserving any 20% will capture roughly all of the predicted percentage of species targeted.Unfortunately, species are not distributed randomly and it is unlikely that the location of every species in the landscape will be known.Real-world applications of these targets will capture larger areas than those predicted by the SAR target due to inefficiencies in adding areas to a conservation network.This stresses the need for at least some species point-locality data to help guide conservation decisions.
Secondly, the z-value describes the accumulation of species based on a single occurrence of a species.A target set using z is analogous to saying, "select one occurrence of each species" in a minimum set.However, it is more than likely that numerous occurrences of common species may actually be incorporated into the notional or real system by the target.This is good, because common species often require larger populations for persistence (Lawton 1988).Also, no vegetation types in the Succulent Karoo have exclusive species complements.Consequently, the same species will be targeted in many vegetation types.Unfortunately, rare or very patchy (habitat-specific) species, i.e., the other 25% of species not targeted, are likely to be missed.Pointlocality data for species that are good surrogates for this group (e.g., rare habitats) are necessary in the conservation planning process.
If the conservation goal is to select at least three occurrences of each species, then the target will have to be increased to accommodate especially rare species.How the SAR method can be used to achieve this needs to be explored.
Lastly, z-values will increase as a result of land-class fragmentation (Rosenzweig 1995).Archipelagos typically have higher z-values than mainland areas.This is generally ascribed to different rates of immigration, extinction, and in situ speciation that occur in island or naturally fragmented biotas (Diamond and May 1976, Bond et al. 1988, Brown and Dinsmore 1988, Rosenzweig 1995, Lomolino and Perault 2001, Haila 2002).Targets derived here assume that a land class is untransformed.Under anthropogenic transformation, however, a larger area than predicted by the model will be required to achieve the same species target.This is a crucial point that needs to be borne in mind when applying this approach.Species relaxation in fragmented landscapes results in a net loss in the original number of species present, as species go extinct from habitat patches over time (Brooks et al. 1999, Robinson 1999, Debinski and Holt 2000, Gonzalez 2000, Kelt 2001).This effectively increases the z-value as a larger area is required to represent the same given percentage of the original species complement.How to adapt SAR targets to landscapes under contemporary transformation needs further investigation.

Which estimator?
Probably the largest source of error in this approach lies in the estimation of the true species number for a vegetation type.There is no consensus in the literature as to which is the best estimator to use (Colwell and Coddington 1995, Chiarucci et al. 2003, Petersen et al. 2003).There is agreement, however, that the Bootstrap estimator is the most conservative (Colwell and Coddington 1995), but no consensus as to which estimator provides the best predictions (Brose 2002, Herzog et al. 2002, Chiarucci et al. 2003, Petersen et al. 2003, Turnipseed et al. 2003).In this analysis it was observed that the Bootstrap estimator was consistently the lowest estimate of the seven estimators calculated.Therefore, the estimates of zvalues and targets calculated in this paper should to be regarded as conservative and probably underestimates of true targets.The rationale for using the Bootstrap technique here is based solely on the observed patterns in estimation error for modeled species data sets (P.Desmet, unpublished data).A better approach may be to use the average of several or all of the estimators.Another approach may be to calculate z-values using all estimators and then deriving a target range.This error, however, does not detract from the utility of the SAR for setting targets, but is rather a source of error in predictions of the model.
The best means of eliminating this error would be not to use the estimators to calculate the true species number, but instead to use an alternative technique for estimating z-values.Two techniques in the literature http://www.ecologyandsociety.org/vol9/iss2/art11hold promise in this regard.Firstly, Harte et al. (1999b) have developed a method of calculating the SAR using species spatial-turnover data.D. P. Faith (personal communication) has proposed a method for calculating z-values using the environmental diversity index (Faith 2003).Both of these techniques can use the same inventory data, but eliminate the need to estimate true species number.
In regions where there are no survey data, but there are inventory data of some form (such as museum collections), it would be useful to explore determining z-values directly from the species-abundance distribution in the pooled inventory data for a land class.Wright (1988) showed that the z-value can be determined directly from this distribution.The assumption here would be that the number of times that a species is recorded in an inventory would be indicative of the species' relative abundance in the land class.This would obviate the need for area-based survey sites, and would also create a novel and very important use for museum data.
As a point of clarification, the z-value cannot be calculated directly by generating a speciesaccumulation curve for a sample of survey sites.This curve is not a species-area curve; rather, it is a collector's curve.A species-area curve is constructed by adding successively larger sampled areas to the data pool until one has sampled the entire land class.The accumulation of species with this progression is then plotted.Survey data are generally sampled at the same spatial scale; hence, simply generating the collector's curve as one randomly adds sites to the data pool would result in significant errors in the z-value.In addition, the collector's curve rarely fits the power model.

Further sources of error
Further bias arises as a result of errors in the survey data.Errors in this type of data are a perennial problem.These include sampling biases leading to uneven sampling of land classes; geo-referencing errors; omission of cryptic species; species identification errors; and data capture and archival errors.No effort is made here to control for these errors beyond the normal checks and balances, such as checking spelling, involved with collating and curating a large biological survey database.Also, the survey data cover a range of projects that span 30 years of research in the Succulent Karoo and involve tens if not hundreds of workers.None of these projects was aimed at landscape-level biodiversity inventory, although phytosociological studies do tend to target all observed plant communities within their respective study areas.The potential for taxonomic errors is high, especially because identification of Mesembryanthema (Aizoaceae), the second largest family in the biome, is notoriously difficult (Smith et al. 1998).No attempts have been made, as yet, to estimate the degree of error in this data set.
Using survey data from a variety of projects that used different relevé sizes is not a significant source of error.Varying the relevé area by an order of magnitude either way (i.e., 10 to 1000 m 2 ) changed the z-value by less than 0.01%.Therefore, knowing the size of the sample relevé does not underpin the application of this technique.Consequently, variable survey area size is not a constraint to using this method.As long as the sample areas are within an order of magnitude of each other, they can be combined for the purposes of estimating targets.
These problems highlight the great need for systematic data collection over a variety of scales, to allow for proper comparisons of z-values at local, regional, and global scales.It is imperative that all of this is done in a highly standardized manner, and having in mind comparisons at the scale considered.Once such data would be available, they would constitute a great starting point for systematic and biologically meaningful target setting.
Another source of error lies in the delimitation of landclass boundaries.The vegetation types used here were mapped using expert knowledge.Errors in where the "true" boundaries of the vegetation types lie can lead to over-or underestimation of z-values.For example, the boundaries of the Eastern Bushmanland Quartz and Gravel Patch vegetation type, as used in this study, are incorrect.In revised versions of the South African vegetation map, this vegetation type has been divided between three vegetation types (one new and two existing; L. Mucina, personal communication).The resultant vegetation map not only agrees better with expert assessment, but also the vegetation types are more homogeneous and better reflect landscape-level vegetation patterns.The consequences for the targets are that they will have to be revisited for this area, and relevés will have to be reassigned to vegetation types according to the new boundaries.The calculated target http://www.ecologyandsociety.org/vol9/iss2/art11General cannot be extrapolated from the old vegetation type to the new ones.Another problem that arises as a result of incorrect vegetation-type boundaries is that relevés are incorrectly assigned to a vegetation type.Realistically, such errors can be controlled only through wide expert involvement in the delimitation of land-class boundaries, whether using expert mapping or modeling techniques.
The idea of using the SAR for setting conservation targets has been around for some time (May 1975, Diamond and May 1976, Rosenzweig 1995, Vreugdenhil et al. 2003).This work, however, represents the first attempt to quantitatively determine conservation targets for land classes using the SAR based on actual biodiversity data.Exciting as this advance is, the limitations of this approach, in terms of input data requirements, data error, and model assumptions, must not be forgotten.

Extrapolating z-values
Within a biogeographic province (i.e., vegetation-type group) there is considerable variation in z-values.There is, however, a generally accepted strong relationship between species diversity and habitat diversity (Holland 1978, Miller et al. 1995, Rosenzweig 1995), which was confirmed in this study.Using topographic diversity as a surrogate for habitat diversity, we constructed a model that relates z-values to an independent land-class metric that can be generated from remotely derived GIS data.The advantages of taking this step are substantial.It is now possible to approximate a z-value for all land classes in the Succulent Karoo, based on a measure of the diversity of habitats.Geographic location is important in explaining the pattern in z-values.This may reflect the historical influence of dynamic environments on the evolution of regional floras in and away from local centers of diversity or refugia.
Conservation practitioners also need to bear in mind that applying SAR targets is only one of many types of conservation targets.This target is based on the hypothesized accumulation of species in a sample of conservation areas.It does not explicitly take into account multiple occurrences of species, nor does it tell us anything about where, within a land class, the target should be achieved to conserve the target percentage of species.Also, these targets do not tell us anything about requirements for ecological processes.It is important to remember that the SAR target does not replace other approaches to setting targets that focus, for example, on minimum viable populations, meta-population dynamics, or ecological processes.
From a practical perspective, the two major limitations of binary or fixed conservation targets need to be stressed.Firstly, fixed targets distort the effectiveness of conservation implementation.Using fixed targets, it is possible to achieve targets for land classes that are just below target, and that probably require minimal effort to achieve, while ignoring those land classes that are far below target and that require significantly more effort to achieve their targets.This approach of picking the low-hanging fruits, although not always the case, exaggerates the success of conservation implementation while exposing the most vulnerable components, i.e., those least conserved, to potentially greater risk through being sidelined by the implementation process.
For the Succulent Karoo vegetation types, where there are not enough vegetation types in a category to build a significant model incorporating topographic diversity, the equally strong observed relationships between z-values and geographic location support the approach of assigning z-values to vegetation types in the same group, using the nearest neighbor principle.This approach would hold at least within vegetation groups within a biome.
It is difficult to make predictions for other biomes.It is almost certain that other biomes will have a range in zvalues and, hence, targets.This study makes it very clear that no single target will be suitable for all land classes within a region.As a very general rule, land classes with large numbers of range-restricted species will have higher targets.Whether there is a relationship between endemicity or rarity and z-values will need to be determined before more empirical statements can be made.Also, more topographically diverse land classes will have relatively higher zvalues.
Secondly, fixed targets promote the land-use philosophy of "clearing down to target."This is a dangerous philosophy, as it is generally accepted that one requires more than just each species represented in a reserve to conserve biodiversity.The SAR target approach applied here does not take into account ecological processes.http://www.ecologyandsociety.org/vol9/iss2/art11Survey plots are little more than slightly-less-thanrandom-samples of the complete biodiversity present at any point in space.These data, however, have formed the basis of much of the ecological research into how terrestrial systems are structured and work.It is important to make the best use of available information rather than wait for better data.Biologically informed decisions need to be made regarding conservation action and landscape management.
The methodology for setting conservation targets presented here is by no means a "save-all" solution to the problem of setting targets.It should be viewed as a tool that complements, rather than replaces, existing empirical or expert-based species, population, habitat, or ecosystem targets.This method is fraught with methodological and data assumptions that need to be addressed, but in the meantime, it would be wise to apply the method, being mindful of its limitations, rather than wait until these problems have been resolved.
It must reiterated that this work makes it clear that the IUCN 10% target is inadequate for capturing the majority of plant diversity within the Succulent Karoo Biome.This trend is probably true for many other terrestrial ecosystems.Further, land classes are not all equal from a biodiversity perspective, and setting a single target for all land classes does not make good conservation sense.

Fig. 1 .
Fig. 1.The location of the study area on the southwestern corner of the African continent.

Fig. 2 .
Fig. 2. The relationship between estimated z-values and: (A) the standard deviation of mean vegetation type altitude (ALT); (B) the ratio between planimetric and surface area of each vegetation type (RATIO); (C) longitude of vegetation type centroid; (D) latitude of vegetation type centroid; and (E) the number of modelled land classes per vegetation type.The Namaqualand and Little Karoo regions of the Succulent Karoo are separated by a dashed line in the longitude (C, Namaqualand to the left) and latitude (D, Namaqualand to the right) graphs.For this figure only the Namaqualand Succulent Karoo vegetation group has been divided into three subgroups: Kamiesberg Klipkop, Lowland Succulent Karoo, and Mountain Succulent Karoo.

Fig. 3 .
Fig. 3. Regression model outputs for the Namaqualand Succulent Karoo vegetation types group using a generalized additive model (Appendix 1).Click on the figure for a larger view.

Table 1 .
The classification of vegetation types into higher order vegetation categories or groups.Vegetation types in each category were used to develop regression models relating observed z-values for vegetation types within each category to remotely determined landscape variables.

Table 2 .
Estimates of true species number for 42 Succulent Karoo vegetation types.Values are rounded to the nearest whole number.Only mean values for ACE, ICE, Jack-knife 2 and Bootstrap are calculated by EstimateS.Values for number of species per plot and Chao1, Chao2, and Jack1 are means (with 1 SD in parentheses).The z-value was calculated using the Bootstrap estimator.

Table 3 .
Targets for Succulent Karoo vegetation types expressed as the percentage area of a vegetation type required to represent a given proportion of the total estimated number of plant species that occur in that vegetation type.Targets are calculated for the range of observed z-values and are rounded off to the nearest whole number.
to this article can be read online at: http://www.ecologyandsociety.org/vol9/iss2/art11/responses/index.htmlhttp://www.ecologyandsociety.org/vol9/iss2/art11df for terms and F values for nonparametric effects Residual standard error: 0.01238 on 3 degrees of freedom Multiple R 2 : 0.7792 F statistic: 10.59 on 1 and 3 degrees of freedom, the P value is 0.04735 2. Model with only the four strict Namaqualand Quartz Patch vegetation types Residual standard error: 0.007422 on 2 degrees of freedom Multiple R 2 : 0.9197 F statistic: 22.92 on 1 and 2 degrees of freedom, P = 0.04097 Excluding the Eastern Bushmanland Quartz and Gravel Patches improves the fit of the model despite the loss of one degree of freedom.The Bushmanland quartz patches should probably be considered a separate biogeographic region, as are the Little Karoo quartz patches {Schmiedel and Jürgens 1999, Schmiedel 2002}.Neither of the topographic diversity variables demonstrated any significant relationship with z-values for quartz patches.This is not surprising, because quartz patch landscapes are characteristically flat areas.Residual standard error: 0.00315 on 3 degrees of freedom Multiple R 2 : 0.9866 F statistic: 73.4 on 3 and 3 degrees of freedom, P = 0.002635