Sensitivity of population viability to spatial and nonspatial parameters using GRIP.

Metapopulation dynamics are influenced by spatial parameters including the amount and arrangement of suitable habitat, yet these parameters may be uncertain when deciding how to manage species or their habitats. Sensitivity analyses of population viability analysis (PVA) models can help measure relative parameter influences on predictions, identify research priorities for reducing uncertainty, and evaluate management strategies. Few spatial PVAs, however, include sensitivity analyses of both spatial and nonspatial parameters, perhaps because computationally efficient tools for such analyses are lacking or inaccessible. We developed GRIP, a program to facilitate sensitivity analysis of spatial and nonspatial input parameters for PVAs created in RAMAS Metapop, a widely applied software program. GRIP creates random sets of input files by varying parameters specified in the PVA model including vital rates and their correlations among populations, the number and configuration of populations, dispersal rates, dispersal survival, initial population abundances, carrying capacities, and the probability, intensity, and spatial extent of catastrophes, while drawing on specified parameter distributions. We evaluated GRIP's performance as a tool for sensitivity analysis of spatial PVAs and explored the consequences of varying spatial input parameters for predictions of a published PVA model of the sand lizard (Lacerta agilis). We used GRIP output to generate standardized regression coefficients (SRCs) and nonparametric correlation coefficients as indices of the relative sensitivity of predicted conservation status to input parameters. GRIP performed well; with a single analysis we were able to rank the relative influence of input parameters identified as influential by the PVA's original author, S. A. Berglind, who used three separate forms of sensitivity analysis. Our analysis, however, also underscored the value of exploring the relative influence of spatial parameters on PVA predictions; both SRCs and correlation coefficients indicated that the most influential parameters in the sand lizard model were spatial in nature. We provide annotated code so that GRIP may be modified to reflect particular species biology, customized for more complex spatial PVA models, upgraded to incorporate features added in newer versions of RAMAS Metapop, used as a template to develop similar programs, or used as it is for computationally efficient sensitivity analyses in support of conservation planning.

[1]  H. Caswell.  Matrix population models : construction, analysis, and interpretation , 2001 .

[2]  Frithjof Lutscher,et al.  Average dispersal success: linking home range, dispersal, and metapopulation dynamics to reserve design. , 2006, Ecological applications : a publication of the Ecological Society of America.

[3]  R. Hilborn,et al.  The Ecological Detective: Confronting Models with Data , 1997 .

[4]  D. Garcelon,et al.  A population viability analysis for the Island Fox on Santa Catalina Island, California , 2005 .

[5]  D. Macdonald,et al.  Simulating the proposed reintroduction of the European beaver (Castor fiber) to Scotland , 2000 .

[6]  Steven R. Beissinger,et al.  Emerging Issues in Population Viability Analysis , 2002, Conservation biology : the journal of the Society for Conservation Biology.

[7]  Kathryn E Sieving,et al.  An Experimental Test of Matrix Permeability and Corridor Use by an Endemic Understory Bird , 2006, Conservation biology : the journal of the Society for Conservation Biology.

[8]  J. A. Jackson,et al.  IVORY-BILLED WOODPECKER (CAMPEPHILUS PRINCIPALIS): HOPE, AND THE INTERFACES OF SCIENCE, CONSERVATION, AND POLITICS , 2006 .

[9]  Rampal S. Etienne,et al.  On optimal choices in increase of patch area and reduction of interpatch distance for metapopulation persistence , 2004 .

[10]  Susan M. Haig,et al.  Modeling Approaches in Avian Conservation and the Role of Field Biologists , 2006 .

[11]  Ilkka Hanski,et al.  Metapopulation Dynamics: From Concepts and Observations to Predictive Models , 1997 .

[12]  M. McCarthy,et al.  Species conservation and management : case studies , 2004 .

[13]  Kyle J. Haynes,et al.  Interpatch movement and edge effects: the role of behavioral responses to the landscape matrix , 2006 .

[14]  T. Ricketts.  The Matrix Matters: Effective Isolation in Fragmented Landscapes , 2001, The American Naturalist.

[15]  Scott Ferson,et al.  Risk assessment in conservation biology , 1993 .

[16]  D. Hosmer,et al.  Applied Logistic Regression , 1991 .

[17]  G. Swartzman,et al.  Ecological Simulation Primer , 1988 .

[18]  Paul C. Cross,et al.  Using Logistic Regression to Analyze the Sensitivity of PVA Models: a Comparison of Methods Based on African Wild Dog Models , 2001 .

[19]  N. Schtickzelle,et al.  Metapopulation dynamics of the bog fritillary butterfly: modelling the effect of habitat fragmentation , 2002 .

[20]  David B. Lindenmayer,et al.  A landscape‐scale test of the predictive ability of a spatially explicit model for population viability analysis , 2001 .

[21]  D. Lindenmayer,et al.  Small mammals, habitat patches and PVA models: a field test of model predictive ability , 2002 .

[22]  N. Schumaker,et al.  Assessing the potential impacts of alternative landscape designs on amphibian population dynamics , 2004, Landscape Ecology.

[23]  Brendan A Wintle,et al.  Incorporating landscape stochasticity into population viability analysis. , 2007, Ecological applications : a publication of the Ecological Society of America.

[24]  H. Kroon,et al.  ELASTICITIES: A REVIEW OF METHODS AND MODEL LIMITATIONS , 2000 .

[25]  David R. Anderson,et al.  Modeling Survival and Testing Biological Hypotheses Using Marked Animals: A Unified Approach with Case Studies , 1992 .

[26]  K. E. Hodges,et al.  Resolving Critical Habitat Designation Failures: Reconciling Law, Policy, and Biology , 2006, Conservation biology : the journal of the Society for Conservation Biology.

[27]  Nicolas Schtickzelle,et al.  Metapopulation dynamics and conservation of the marsh fritillary butterfly: Population viability analysis and management options for a critically endangered species in Western Europe , 2005 .

[28]  D. Ludwig.  Is it meaningful to estimate a probability of extinction , 1999 .

[29]  Scott Ferson,et al.  Sensitivity analysis for models of population viability , 1995 .

[30]  John R Fieberg,et al.  Assessing uncertainty in ecological systems using global sensitivity analyses: a case example of simulated wolf reintroduction effects on elk , 2005 .