Forecasting extinction risk with nonstationary matrix models.

Matrix population growth models are standard tools for forecasting population change and for managing rare species, but they are less useful for predicting extinction risk in the face of changing environmental conditions. Deterministic models provide point estimates of lambda, the finite rate of increase, as well as measures of matrix sensitivity and elasticity. Stationary matrix models can be used to estimate extinction risk in a variable environment, but they assume that the matrix elements are randomly sampled from a stationary (i.e., non-changing) distribution. Here we outline a method for using nonstationary matrix models to construct realistic forecasts of population fluctuation in changing environments. Our method requires three pieces of data: (1) field estimates of transition matrix elements, (2) experimental data on the demographic responses of populations to altered environmental conditions, and (3) forecasting data on environmental drivers. These three pieces of data are combined to generate a series of sequential transition matrices that emulate a pattern of long-term change in environmental drivers. Realistic estimates of population persistence and extinction risk can be derived from stochastic permutations of such a model. We illustrate the steps of this analysis with data from two populations of Sarracenia purpurea growing in northern New England. Sarracenia purpurea is a perennial carnivorous plant that is potentially at risk of local extinction because of increased nitrogen deposition. Long-term monitoring records or models of environmental change can be used to generate time series of driver variables under different scenarios of changing environments. Both manipulative and natural experiments can be used to construct a linking function that describes how matrix parameters change as a function of the environmental driver. This synthetic modeling approach provides quantitative estimates of extinction probability that have an explicit mechanistic basis.

[1]  N. Gotelli.  Demographic Models for Leptogorgia Virgulata, A Shallow‐Water Gorgonian , 1991 .

[2]  A. Ellison,et al.  Nitrogen availability alters the expression of carnivory in the northern pitcher plant, Sarracenia purpurea , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[3]  D. Schnell.  Carnivorous Plants of the United States and Canada , 1976 .

[4]  Brian Dennis,et al.  ESTIMATING POPULATION TREND AND PROCESS VARIATION FOR PVA IN THE PRESENCE OF SAMPLING ERROR , 2004 .

[5]  Elena R. Alvarez-Buylla,et al.  Sustainable Harvesting of Tropical Trees: Demography and Matrix Models of Two Palm Species in Mexico , 1995 .

[6]  J. Halley.  Parameter drift stabilizes long‐range extinction forecasts , 2003 .

[7]  David R. Anderson,et al.  Some specialized risk assessment methodologies for vertebrate populations , 1995, Environmental and Ecological Statistics.

[8]  Eric S. Menges,et al.  Interpretation of elasticity matrices as an aid to the management of plant populations for conservation , 1996 .

[9]  A. Ellison,et al.  NITROGEN DEPOSITION AND EXTINCTION RISK IN THE NORTHERN PITCHER PLANT, SARRACENIA PURPUREA , 2002 .

[10]  Aaron M. Ellison,et al.  A Primer of Ecological Statistics , 2004 .

[11]  J. Silvertown,et al.  Comparing plant life histories using elasticity analysis: the importance of life span and the number of life-cycle stages , 1995, Oecologia.

[12]  M. Wisdom,et al.  Sensitivity analysis to guide population recovery : Prairie-chickens as an example , 1997 .

[13]  L. Crowder,et al.  population viability analysis for red-cockaded woodpeckers using an individual-based model , 2002 .

[14]  H. Caswell.  Analysis of life table response experiments I. Decomposition of effects on population growth rate , 1989 .

[15]  A. Grant,et al.  How to Keep Fit in the Real World: Elasticity Analyses and Selection Pressures on Life Histories in a Variable Environment , 1996, The American Naturalist.

[16]  S. Ellner.  When does parameter drift decrease the uncertainty in extinction risk estimates , 2003 .

[17]  A. Grant,et al.  ELASTICITY ANALYSIS FOR DENSITY‐DEPENDENT POPULATIONS IN STOCHASTIC ENVIRONMENTS , 2000 .

[18]  Janusz Uchmański,et al.  What promotes persistence of a single population: an individual-based model , 1999 .

[19]  J. Brewer.  A demographic analysis of fire-stimulated seedling establishment of Sarracenia alata (Sarraceniaceae). , 2001, American journal of botany.

[20]  A. Łomnicki.  Individual-based models and the individual-based approach to population ecology , 1999 .

[21]  C. Horvitz,et al.  HURRICANE DISTURBANCE AND THE POPULATION DYNAMICS OF A TROPICAL UNDERSTORY SHRUB: MEGAMATRIX ELASTICITY ANALYSIS , 1998 .

[22]  Brian Dennis,et al.  Joint effects of density dependence and rainfall on abundance of San Joaquin kit fox , 2000 .

[23]  M. Boyce.  Population Viability Analysis , 1992 .

[24]  Thomas J. Givnish,et al.  Carnivory in the Bromeliad Brocchinia reducta, with a Cost/Benefit Model for the General Restriction of Carnivorous Plants to Sunny, Moist, Nutrient-Poor Habitats , 1984, The American Naturalist.

[25]  B. Gilland.  World population and food supply: can food production keep pace with population growth in the next half-century? , 2002 .

[26]  Shripad Tuljapurkar,et al.  The Many Growth Rates and Elasticities of Populations in Random Environments , 2003, The American Naturalist.

[27]  N. Picard,et al.  Modeling forest dynamics with a combined matrix/individual-based model , 2002 .

[28]  William F. Morris,et al.  DETECTING POPULATION‐LEVEL CONSEQUENCES OF ONGOING ENVIRONMENTAL CHANGE WITHOUT LONG‐TERM MONITORING , 1999 .

[29]  C. Pfister,et al.  Patterns of variance in stage-structured populations: evolutionary predictions and ecological implications. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Christopher B. Field,et al.  Nitrogen and Climate Change , 2003, Science.

[31]  T. O'Shea,et al.  Population Viability Analysis of the Florida Manatee (Trichechus manatus latirostris), 1976–1991 , 1997 .

[32]  Paulette Bierzychudek,et al.  The Demography of Jack‐in‐the‐Pulpit, a Forest Perennial that Changes Sex , 1982 .

[33]  Paul Beier,et al.  Population viability analysis. , 2016 .

[34]  A. Ellison,et al.  PREY ADDITION ALTERS NUTRIENT STOICHIOMETRY OF THE CARNIVOROUS PLANT SARRACENIA PURPUREA , 2005 .

[35]  B. Kendall,et al.  Unstructured Individual Variation and Demographic Stochasticity , 2003 .

[36]  Hal Caswell,et al.  Elasticity: The Relative Contribution of Demographic Parameters to Population Growth Rate , 1986 .

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

[38]  Michael C. Dietze,et al.  COEXISTENCE: HOW TO IDENTIFY TROPHIC TRADE-OFFS , 2003 .

[39]  Brian D. Inouye,et al.  RESPONSE SURFACE EXPERIMENTAL DESIGNS FOR INVESTIGATING INTERSPECIFIC COMPETITION , 2001 .

[40]  E. Cowling,et al.  Reactive Nitrogen and The World: 200 Years of Change , 2002, Ambio.

[41]  P. V. Tienderen,et al.  Life Cycle Trade-Offs in Matrix Population Models , 1995 .

[42]  B. Craig,et al.  BAYESIAN ESTIMATION OF A DEMOGRAPHIC MATRIX MODEL FROM STAGE-FREQUENCY DATA , 2002 .

[43]  K. Gross.  Efficient data collection for estimating growth rates of structured populations , 2002 .

[44]  Shripad Tuljapurkar,et al.  Loop analysis: evaluating life history pathways in population projection matrices. , 1994 .

[45]  W. Hoffmann.  FIRE AND POPULATION DYNAMICS OF WOODY PLANTS IN A NEOTROPICAL SAVANNA: MATRIX MODEL PROJECTIONS , 1999 .

[46]  B. Efron.  Nonparametric estimates of standard error: The jackknife, the bootstrap and other methods , 1981 .

[47]  Larry B. Crowder,et al.  Models to Evaluate Headstarting as a Management Tool for Long‐Lived Turtles , 1996 .

[48]  E. Cortés.  Incorporating Uncertainty into Demographic Modeling: Application to Shark Populations and Their Conservation , 2002 .

[49]  D. Pyke,et al.  THE EFFECT OF STOCHASTIC TECHNIQUE ON ESTIMATES OF POPULATION VIABILITY FROM TRANSITION MATRIX MODELS , 2003 .

[50]  Michael J. Wisdom,et al.  Life Stage Simulation Analysis: Estimating Vital-Rate Effects on Population Growth for Conservation , 2000 .

[51]  J. Brewer.  Short-term effects of fire and competition on growth and plasticity of the yellow pitcher plant, Sarracenia alata (Sarraceniaceae). , 1999, American journal of botany.

[52]  Kirk A. Moloney,et al.  A generalized algorithm for determining category size , 1986, Oecologia.

[53]  John Vandermeer,et al.  Choosing category size in a stage projection matrix , 2004, Oecologia.

[54]  Dustin J Penn,et al.  The Evolutionary Roots of Our Environmental Problems: Toward a Darwinian Ecology , 2003, The Quarterly Review of Biology.