Planning for robust reserve networks using uncertainty analysis

Planning land-use for biodiversity conservation frequently involves computer-assisted reserve selection algorithms. Typically such algorithms operate on matrices of species presence–absence in sites, or on species-specific distributions of model predicted probabilities of occurrence in grid cells. There are practically always errors in input data—erroneous species presence–absence data, structural and parametric uncertainty in predictive habitat models, and lack of correspondence between temporal presence and long-run persistence. Despite these uncertainties, typical reserve selection methods proceed as if there is no uncertainty in the data or models. Having two conservation options of apparently equal biological value, one would prefer the option whose value is relatively insensitive to errors in planning inputs. In this work we show how uncertainty analysis for reserve planning can be implemented within a framework of information-gap decision theory, generating reserve designs that are robust to uncertainty. Consideration of uncertainty involves modifications to the typical objective functions used in reserve selection. Search for robust-optimal reserve structures can still be implemented via typical reserve selection optimization techniques, including stepwise heuristics, integer-programming and stochastic global search.     

Variance and Uncertainty in the Expected Number of Occurrences in Reserve Selection

Abstract:  Reserve selection often concerns the design of reserve networks for the long-term maintenance of biodiversity. We considered uncertainty in the context of three common reserve-selection formulations, the expected number of populations, proportional coverage of land-cover types, and the probability of having at least one population. By uncertainty, we mean variance in the outcome of any probability-based reserve selection formulation. A typical reserve-selection formulation might ask for the least expensive set of sites that contains n populations per species. It is implicit here that this requirement concerns the expected number of populations, which actually is obtained only with a 50% chance. If the requirement is changed to select the least expensive set of sites that gives n populations per species with a 95% probability, the number of sites required in the solution increases and the identity of the sites is changed toward sites that have high probabilities of persistence (or occurrence) and low associated binomial variance. Anthropogenic threat is one factor that may cause probabilistic uncertainty in the context of proportional area coverage.     

Identifying Conservation Areas on the Basis of Alternative Distribution Data Sets

Abstract:  Distribution data on biodiversity features is a major component of conservation planning that are often inaccurate; thus, the true distribution of each feature is commonly over- or underrepresented. The selection of distribution data sets may therefore lead to variability in the spatial configuration and size of proposed reserve networks and uncertainty regarding the extent to which these networks actually contain the biodiversity features they were identified to protect. Our goals were to investigate the impact on reserve selection of choosing different distribution data sets and to propose novel methods to minimize uncertainty about target attainment within reserves. To do so, we used common prioritization methods (richness mapping, systematic reserve design, and a novel approach that integrates multiple types of distribution data) and three types of data on the distribution of mammals (predicted distribution models, occurrence records, and a novel combination of the two) to simulate the establishment of regional biodiversity reserves for the state of Arizona (U.S.A.). Using the results of these simulations, we explored variability in reserve placement and size as a function of the distribution data set. Spatial overlap of reserve networks identified with only predicted distribution data or only occurrence distribution data never exceeded 16%. In pairwise comparisons between reserves created with all three types of distribution data, overlap never achieved 50%. The reserve size required to meet conservation targets also varied with the type of distribution data used and the conservation goal; the largest reserve system was 10 times the smallest. Our results highlight the impact of employing different types of distribution data and identify novel tools for application to existing distribution data sets that can minimize uncertainty about target attainment.     

Effectiveness of Vegetation Surrogates for Parasitoid Wasps in Reserve Selection

Abstract:  Selecting suitable nature reserves is a continuing challenge in conservation, particularly for target groups that are time-consuming to survey, species rich, and extinction prone. One such group is the parasitoid Hymenoptera, which have been excluded from conservation planning. If basic characteristics of habitats or vegetation could be used as reliable surrogates of specific target taxa, this would greatly facilitate appropriate reserve selection. We identified a range of potential habitat indicators of the species richness of pimpline parasitoid communities (Hymenoptera: Ichneumonidae: Pimplinae, Diacritinae, Poemeniinae) and tested their efficiency at capturing the observed diversity in a group of small woodlands in the agricultural landscape of the Vale of York (United Kingdom). Eight of the 18 vegetation-based reserve-selection strategies were significantly better at parasitoid species inclusion than random selection of areas. The best strategy maximized richness of tree species over the entire reserve network through complementarity. This strategy omitted only 2–3 species more (out of 38 captured in the landscape as a whole) than selections derived from the parasitoid survey data. In general, strategies worked equally well at capturing species richness and rarity. Our results suggest that vegetation data as a surrogate for species richness could prove an informative tool in parasitoid conservation, but further work is needed to test how broadly applicable these indicators may be.     

Sensitivity of Systematic Reserve Selection to Decisions about Scale, Biological Data, and Targets: Case Study from Southern British Columbia

Abstract:  The identification of conservation areas based on systematic reserve-selection algorithms requires decisions related to both spatial and ecological scale. These decisions may affect the distribution and number of sites considered priorities for conservation within a region. We explored the sensitivity of systematic reserve selection by altering values of three essential variables. We used a 1:20,000–scale terrestrial ecosystem map and habitat suitability data for 29 threatened vertebrate species in the Okanagan region of British Columbia, Canada. To these data we applied a reserve-selection algorithm to select conservation sites while altering selection unit size and shape, features of biodiversity (i.e., vertebrate species), and area conservation targets for each biodiversity feature. The spatial similarity, or percentage overlap, of selected sets of conservation sites identified (1) with different selection units was ≤40%, (2) with different biodiversity features was 59%, and (3) with different conservation targets was ≥94%. Because any selected set of sites is only one of many possible sets, we also compared the conservation value (irreplaceability) of all sites in the region for each variation of the data. The correlations of irreplaceability were weak for different selection units (0.23 ≤ r ≤ 0.67), strong for different biodiversity features ( r = 0.84), and mixed for different conservation targets ( r = 0.16; 0.16; 1.00). Because of the low congruence of selected sites and weak correlations of irreplaceability for different selection units, recommendations from studies that have been applied at only one spatial scale must be considered cautiously.     

Incorporating consumer–resource spatial interactions in reserve design

The persistence of species in reserves depends in large part on the persistence of functional ecological interactions. Despite their importance, however, ecological interactions have not yet been explicitly incorporated into conservation prioritization methods. We develop here a general method for incorporating consumer–resource interactions into spatial reserve design. This method protects spatial consumer–resource interactions by protecting areas that maintain the connectivity between the distribution of consumers and resources. We illustrate our method with a conservation planning case study of a mammalian predator, American marten (Martes americana), and its two primary prey species, Red-backed vole (Clethrionomys rutilus) and Deer mouse (Peromyscus maniculatus). The conservation goal was to identify a reserve for marten that comprised 12% of a forest management unit in the boreal forest in Québec, Canada. We compared reserves developed using analysis variants that utilized different levels of information about predator and prey habitat distributions, species-specific connectivity requirements, and interaction connectivity requirements. The inclusion of consumer–resource interactions in reserve-selection resulted in spatially aggregated reserves that maintained local habitat quality for the species. This spatial aggregation was not induced by applying a qualitative penalty for the boundary length of the reserve, but rather was a direct consequence of modelling the spatial needs of the interacting consumer and resources. Our method for maintaining connectivity between consumers and their resources within reserves can be applied even under the most extreme cases of either complete spatial overlap or complete spatial segregation of consumer–resource distributions. The method has been made available via public software.     

The Value of Biodiversity in Reserve Selection: Representation, Species Weighting, and Benefit Functions

Abstract:  The limited availability of resources for conservation has led to the development of many quantitative methods for selecting reserves that aim to maximize the biodiversity value of reserve networks. In published analyses, species are often considered equal, although some are in much greater need of protection than others. Furthermore, representation is usually treated as a threshold: a species is either represented or not, but varying levels of representation over or under a given target level are not valued differently. We propose that a higher representation level should also have higher value. We introduce a framework for reserve selection that includes species weights and benefit functions for under- and overrepresentation (number of locations for each species). We applied the method to conservation planning for herb-rich forests in southern Finland. Our use of benefit functions and weighting changed the identity of about 50% of the selected sites at different funding levels and improved the representation of rare and threatened species. We also identified a small area of additional land that would substantially enhance the existing reserve network. We suggest that benefit functions and species weighting should be considered as standard options in reserve-selection applications.     

Avoiding Pitfalls of Using Species Distribution Models in Conservation Planning

Abstract:   Museum records have great potential to provide valuable insights into the vulnerability, historic distribution, and conservation of species, especially when coupled with species-distribution models used to predict species' ranges. Yet, the increasing dependence on species-distribution models in identifying conservation priorities calls for a more critical evaluation of model robustness. We used 11 bird species of conservation concern in Brazil's highly fragmented Atlantic Forest and data on environmental conditions in the region to predict species distributions. These predictions were repeated for five different model types for each of the 11 bird species. We then combined these species distributions for each model separately and applied a reserve-selection algorithm to identify priority sites. We compared the potential outcomes from the reserve selection among the models. Although similarity in identification of conservation reserve networks occurred among models, models differed markedly in geographic scope and flexibility of reserve networks. It is essential for planners to evaluate the conservation implications of false-positive and false-negative errors for their specific management scenario before beginning the modeling process. Reserve networks selected by models that minimized false-positive errors provided a better match with priority areas identified by specialists. Thus, we urge caution in the use of models that overestimate species' occurrences because they may misdirect conservation action. Our approach further demonstrates the great potential value of museum records to biodiversity studies and the utility of species-distribution models to conservation decision-making. Our results also demonstrate, however, that these models must be applied critically and cautiously.     

The Boundary-Quality Penalty: a Quantitative Method for Approximating Species Responses to Fragmentation in Reserve Selection

Abstract:  Aggregation of reserve networks is generally considered desirable for biological and economic reasons: aggregation reduces negative edge effects and facilitates metapopulation dynamics, which plausibly leads to improved persistence of species. Economically, aggregated networks are less expensive to manage than fragmented ones. Therefore, many reserve-design methods use qualitative heuristics, such as distance-based criteria or boundary-length penalties to induce reserve aggregation. We devised a quantitative method that introduces aggregation into reserve networks. We call the method the boundary-quality penalty (BQP) because the biological value of a land unit (grid cell) is penalized when the unit occurs close enough to the edge of a reserve such that a fragmentation or edge effect would reduce population densities in the reserved cell. The BQP can be estimated for any habitat model that includes neighborhood (connectivity) effects, and it can be introduced into reserve selection software in a standardized manner. We used the BQP in a reserve-design case study of the Hunter Valley of southeastern Australia. The BQP resulted in a more highly aggregated reserve network structure. The degree of aggregation required was specified by observed (albeit modeled) biological responses to fragmentation. Estimating the effects of fragmentation on individual species and incorporating estimated effects in the objective function of reserve-selection algorithms is a coherent and defensible way to select aggregated reserves. We implemented the BQP in the context of the Zonation method, but it could as well be implemented into any other spatially explicit reserve-planning framework .     

The Effect of Incremental Reserve Design and Changing Reservation Goals on the Long-Term Efficiency of Reserve Systems

Abstract:  Selecting reserve areas based on percentages, such as 10% or 12% of a bioregion, is common in conservation planning despite widespread admission that such percentages are arbitrary and likely to be inadequate for the conservation of all biodiversity. Reserve systems based on these relatively low percentage targets are likely to require expansion in the future, resulting in the assembly of reserve systems over many years (incremental reserve design). How then will incremental reserve design, such as increasing percentage targets over time, affect the long-term efficiency of marine reserve systems? We used South Australia as a case study to investigate how changing percentage targets affects the contribution of individual planning units to efficient reserve design. Selection frequency counts provided a measure of a planning unit's conservation value. For the majority of planning units, changing targets led to a change in their conservation value indicating, for example, that planning units identified as high-value sites at a low-percentage conservation target may be of lesser importance when targets are increased. Despite the variability in the value of individual planning units at different targets, there was no loss in efficiency from incremental design of reserve systems based on systematic methods compared with purpose-built reserve systems (i.e., the system is assembled in a single iteration). The exception was when incrementally designed systems were based on South Australia's existing marine reserve system—a system developed in an ad hoc method. The result was reserve systems that were less efficient, less compact, and larger in size. This suggests that systematic approaches have an important role for efficient reserve design when there is uncertainty about the target level of reservation .     

Where Should Nature Reserves Be Located in South Africa? A Snake's Perspective

The present dispersion of nature reserves in South Africa is the historical result of a series of ad hoc decisions and may not be biologically optimal We have adopted a method to identify the optimal geography of nature reserves for the conservation of South Africa's snake fauna. The method of reserve selection operated on two tiers, and the spatial unit of analysis was a quarter-degree-square cell (∼625 km²). First, two scoring indices were used to rank cells with respect to species richness or to rarity. Second, two different iterative reserve-selection algorithms selected sets of cells (reserves), where each set represented all snake species at least once. Finally, the selected cells were examined for their present level of protection and their ranked scores. Depending on the algorithm chosen, only 23 or 29 cells were required to represent all species at least once; 72–78% of these cells already contained some level of protection; 59–70% of cells fell in areas of high species richness; and 72–91% of cells fell in areas with high rarity scores. Thus we conclude that most of the snake species in South Africa may be adequately protected with only modest acquisition of new reserves, and that the iterative algorithms appear to be efficient at siting cells in areas of high richness and rarity. We recommend that the reserve placement method outlined in this report be applied to as many other taxa as possible in the formulation of a national plan for an optimal reserve system for South Africa.     

Use of Inverse Spatial Conservation Prioritization to Avoid Biological Diversity Loss Outside Protected Areas

Globally expanding human land use sets constantly increasing pressure for maintenance of biological diversity and functioning ecosystems. To fight the decline of biological diversity, conservation science has broken ground with methods such as the operational model of systematic conservation planning (SCP), which focuses on design and on‐the‐ground implementation of conservation areas. The most commonly used method in SCP is reserve selection that focuses on the spatial design of reserve networks and their expansion. We expanded these methods by introducing another form of spatial allocation of conservation effort relevant for land‐use zoning at the landscape scale that avoids negative ecological effects of human land use outside protected areas. We call our method inverse spatial conservation prioritization. It can be used to identify areas suitable for economic development while simultaneously limiting total ecological and environmental effects of that development at the landscape level by identifying areas with highest economic but lowest ecological value. Our method is not based on a priori targets, and as such it is applicable to cases where the effects of land use on, for example, individual species or ecosystem types are relatively small and would not lead to violation of regional or national conservation targets. We applied our method to land‐use allocation to peat mining. Our method identified a combination of profitable production areas that provides the needed area for peat production while retaining most of the landscape‐level ecological value of the ecosystem. The results of this inverse spatial conservation prioritization are being used in land‐use zoning in the province of Central Finland.     

A guide to representing variability and uncertainty in biodiversity indicators

Biodiversity indicators are used to inform decisions and measure progress toward global targets, such as the United Nations Sustainable Development Goals. Indicators aggregate and simplify complex information, so underlying information influencing its reliability and interpretation (e.g., variability in data and uncertainty in indicator values) can be lost. Communicating uncertainty is necessary to ensure robust decisions and limit misinterpretations of trends, yet variability and uncertainty are rarely quantified in biodiversity indicators. We developed a guide to representing uncertainty and variability in biodiversity indicators. We considered the key purposes of biodiversity indicators and commonly used methods for representing uncertainty (standard error, bootstrap resampling, and jackknife resampling) and variability (quantiles, standard deviation, median absolute deviation, and mean absolute deviation) with intervals. Using 3 high-profile biodiversity indicators (Red List Index, Living Planet Index, and Ocean Health Index), we tested the use, suitability, and interpretation of each interval method based on the formulation and data types underpinning the indicators. The methods revealed vastly different information; indicator formula and data distribution affected the suitability of each interval method. Because the data underpinning each indicator were not normally distributed, methods relying on normality or symmetrical spread were unsuitable. Quantiles, bootstrapping, and jackknifing provided useful information about the underlying variability and uncertainty. We built a decision tree to inform selection of the appropriate interval method to represent uncertainty or variation in biodiversity indicators, depending on data type and objectives. Our guide supports transparent and effective communication of biodiversity indicator trends to facilitate accurate interpretation by decision makers.     

Habitat Suitability Models and the Shortfall in Conservation Planning for African Vertebrates

Abstract:  Ongoing loss of biodiversity requires identifying large-scale conservation priorities, but the detailed information on the distribution of species required for this purpose is often missing. We present a systematic reserve selection for 1223 African mammals and amphibians in which habitat suitability models are used as estimates of the area occupied by species. In the framework of the World Conservation Union (IUCN) Global Amphibian Assessment and IUCN Global Mammal Assessment, we collected the geographic range (extent of occurrence) and habitat preferences for each species. We used the latter to build species-specific habitat suitability models inside geographic ranges, and for 181 species we verified the models by comparing suitability levels to presence-absence data collected in the field. We then used the suitable areas as estimators of the area of occupancy and compared the results of systematic reserve selection based on geographic ranges to those based on estimated areas of occupancy. Our results showed that the reserve system would need a 30-100% expansion to achieve minimal conservation targets, concentrated in the tropics, where species richness reaches a maximum. Comparative analyses revealed that using geographic ranges, which overestimate the area occupied by species, underestimates the total amount of area that needs to be conserved. The area selected for conservation doubled when we used the estimated area of occupancy in place of the geographic ranges. This happened because the suitable areas potentially occupied by each species overlapped less than their geographic ranges. As a result, any given protected area contained fewer species than predicted by the analysis of ranges. Because species are more specialized than our estimates of distribution based on extent of occurrence suggest, we propose that this is a general effect in systematic conservation planning.     

Confronting Uncertainty and Missing Values in Environmental Value Transfer as Applied to Species Conservation

Abstract: The nonuse (or passive) value of nature is important but time‐consuming and costly to quantify with direct surveys. In the absence of estimates of these values, there will likely be less investment in conservation actions that generate substantial nonuse benefits, such as conservation of native species. To help overcome decisions about the allocation of conservation dollars that reflect the lack of estimates of nonuse values, these values can be estimated indirectly by environmental value transfer (EVT). EVT uses existing data or information from a study site such that the estimated monetary value of an environmental good is transferred to another location or policy site. A major challenge in the use of EVT is the uncertainty about the sign and size of the error (i.e., the percentage by which transferred value exceeds the actual value) that results from transferring direct estimates of nonuse values from a study to a policy site, the site where the value is transferred. An EVT is most useful if the decision‐making framework does not require highly accurate information and when the conservation decision is constrained by time and financial resources. To account for uncertainty in the decision‐making process, a decision heuristic that guides the decision process and illustrates the possible decision branches, can be followed. To account for the uncertainty associated with the transfer of values from one site to another, we developed a risk and simulation approach that uses Monte Carlo simulations to evaluate the net benefits of conservation investments and takes into account different possible distributions of transfer error. This method does not reduce transfer error, but it provides a way to account for the effect of transfer error in conservation decision making. Our risk and simulation approach and decision‐based framework on when to use EVT offer better‐informed decision making in conservation.     

Finding all optimal solutions to the reserve site selection problem: formulation and computational analysis

The problem of selecting nature reserves has received increased attention in the literature during the past decade, and a variety of approaches have been promoted for selecting those sites to include in a reserve network. One set of techniques employs heuristic algorithms and thus provides possibly sub-optimal solutions. Another set of models and accompanying algorithms uses an integer programming formulation of the problem, resulting in an optimization problem known as the Maximal Covering Problem, or MCP. Solution of the MCP provides an optimal solution to the reserve site selection problem, and while various algorithms can be employed for solving the MCP they all suffer from the disadvantage of providing a single optimal solution dictating the selection of areas for conservation. In order to provide complete information to decision makers, the determination of all alternate optimal solutions is necessary. This paper explores two procedures for finding all such solutions. We describe the formulation and motivation of each method. A computational analysis on a data set describing native terrestrial vertebrates in the state of Oregon illustrates the effectiveness of each approach.     

Combining Strategies to Select Reserves in Fragmented Landscapes

Abstract:  In the identification of reserve networks in fragmented landscapes with limited species-specific data at hand, one approach is to use selection criteria, such as patch size, to rank the habitat patches' conservation value and evaluate reserve-network alternatives. These criteria are assumed to be reasonable surrogates for the true network objectives. Caution is warranted, however, because the relationships between the selection criteria and the reserve-network objectives may be inconsistent. Conflicts are also likely to arise because no single reserve network will be optimal with respect to multiple objectives (or selection criteria) simultaneously. Instead, reserve planners must compromise between conflicting demands. We field tested the relationships between a variety of selection criteria and the objectives of a reserve network for the sandplain natural communities on Martha's Vineyard Island, Massachusetts (U.S.A.). Selection criteria that correlated with the reserve-network objectives were used in a multi-objective integer program to identify the 10-patch reserve networks that were optimal with each objective independently and those that offered optimal tradeoffs between the reserve-network objectives. From these 10-patch networks, one can select a final reserve network that provides the preferred compromise between the objectives.     

Effect of risk aversion on prioritizing conservation projects

Conservation outcomes are uncertain. Agencies making decisions about what threat mitigation actions to take to save which species frequently face the dilemma of whether to invest in actions with high probability of success and guaranteed benefits or to choose projects with a greater risk of failure that might provide higher benefits if they succeed. The answer to this dilemma lies in the decision maker's aversion to risk—their unwillingness to accept uncertain outcomes. Little guidance exists on how risk preferences affect conservation investment priorities. Using a prioritization approach based on cost effectiveness, we compared 2 approaches: a conservative probability threshold approach that excludes investment in projects with a risk of management failure greater than a fixed level, and a variance‐discounting heuristic used in economics that explicitly accounts for risk tolerance and the probabilities of management success and failure. We applied both approaches to prioritizing projects for 700 of New Zealand's threatened species across 8303 management actions. Both decision makers’ risk tolerance and our choice of approach to dealing with risk preferences drove the prioritization solution (i.e., the species selected for management). Use of a probability threshold minimized uncertainty, but more expensive projects were selected than with variance discounting, which maximized expected benefits by selecting the management of species with higher extinction risk and higher conservation value. Explicitly incorporating risk preferences within the decision making process reduced the number of species expected to be safe from extinction because lower risk tolerance resulted in more species being excluded from management, but the approach allowed decision makers to choose a level of acceptable risk that fit with their ability to accommodate failure. We argue for transparency in risk tolerance and recommend that decision makers accept risk in an adaptive management framework to maximize benefits and avoid potential extinctions due to inefficient allocation of limited resources. El Efecto de la Aversión de Riesgo sobre la Priorización de Proyectos de Conservación     

Sensitivity of species-distribution models to error,bias, and model design: An application to resource selection functions for woodland caribou

Models that predict distribution are now widely used to understand the patterns and processes of plant and animal occurrence as well as to guide conservation and management of rare or threatened species. Application of these methods has led to corresponding studies evaluating the sensitivity of model performance to requisite data and other factors that may lead to imprecise or false inferences. We expand upon these works by providing a relative measure of the sensitivity of model parameters and prediction to common sources of error, bias, and variability. We used a one-at-a-time sample design and GPS location data for woodland caribou (Rangifer tarandus caribou) to assess one common species-distribution model: a resource selection function. Our measures of sensitivity included change in coefficient values, prediction success, and the area of mapped habitats following the systematic introduction of geographic error and bias in occurrence data, thematic misclassification of resource maps, and variation in model design. Results suggested that error, bias and model variation have a large impact on the direct interpretation of coefficients. Prediction success and definition of important habitats were less responsive to the perturbations we introduced to the baseline model. Model coefficients, prediction success, and area of ranked habitats were most sensitive to positional error in species locations followed by sampling bias, misclassification of resources, and variation in model design. We recommend that researchers report, and practitioners consider, levels of error and bias introduced to predictive species-distribution models. Formal sensitivity and uncertainty analyses are the most effective means for evaluating and focusing improvements on input data and considering the range of values possible from imperfect models.     

Making conservation decisions under uncertainty for the persistence of multiple species.

Population models for multiple species provide one of the few means of assessing the impact of alternative management options on the persistence of biodiversity, but they are inevitably uncertain. Is it possible to use population models in multiple-species conservation planning given the associated uncertainties? We use information-gap decision theory to explore the impact of parameter uncertainty on the conservation decision when planning for the persistence of multiple species. An information-gap approach seeks robust outcomes that are most immune from error. We assess the impact of uncertainty in key model parameters for three species, whose extinction risks under four alternative management scenarios are estimated using a metapopulation model. Three methods are described for making conservation decisions across the species, taking into account uncertainty. We find that decisions based on single species are relatively robust to uncertainty in parameters, although the estimates of extinction risk increase rapidly with uncertainty. When identifying the best conservation decision for the persistence of all species, the methods that rely on the rankings of the management options by each species result in decisions that are similarly robust to uncertainty. Methods that depend on absolute values of extinction risk are sensitive to uncertainty, as small changes in extinction risk can alter the ranking of the alternative scenarios. We discover that it is possible to make robust conservation decisions even when the uncertainties of the multiple-species problem appear overwhelming. However, the decision most robust to uncertainty is likely to differ from the best decision when uncertainty is ignored, illustrating the importance of incorporating uncertainty into the decision-making process.