Integral Projection Models for trees: a new parameterization method and a validation of model output

Summary 1. Matrix models are popular tools for plant demographic studies, but their application to long-lived, slow-growing species is hampered by the fact that (i) model output is highly sensitive to category width and (ii) growth variation between individuals can only be partially accounted for. Integral Projection Models (IPMs) – an extension of matrix models – offer a solution to these problems. 2. Here, we introduce a new method to parameterize IPMs for trees – the ‘integration method’ – which allows constructing IPMs for long-lived, slow-growing species. This approach is more suitable than the ‘midpoint rule’, which is customarily used. 3. We built IPMs for six tree species from Vietnamese (sub)tropical forests. For four of these species, population growth rate (k) was highly sensitive to the number of categories in the transition matrix. Population growth stabilized for IPMs with 100–1000 categories, corresponding to categories of 0.1–1 cm in trunk diameter. This preferred width is much narrower than the 10-cm-wide categories customarily used in tree models. 4. The distribution of elasticity values over transition types (stasis, progression to next and further categories) is also highly sensitive to matrix dimension in IPMs. In addition, elasticity distribution is influenced by including or excluding growth variation. 5. Age estimates obtained from IPMs were also highly sensitive to matrix dimension: an IPM with 1000 size categories yielded 2–4 times higher age estimates for large trees than one with 10 size categories. Observed ages obtained from tree ring analyses for four of the study species allowed validating these estimates. IPMs with 10 categories strongly underestimated age, while those with 1000 categories yielded slight age overestimates. Underestimating age in small matrices is caused by the occurrence of unrealistically fast pathways through the life cycle and is probably widespread among tree models with broad categories. Overestimating ages in IPMs with narrow categories may be due to temporally autocorrelated growth or errors in fitting growth curves. 6. Synthesis. IPMs are highly suitable tools to analyse tree demography. We recommend that tree IPMs (and classical matrix models) apply narrow diameter categories (0.1–1 cm width) to obtain reliable model output.

[1]  E. Álvarez-Buylla,et al.  How old are tropical rain forest trees , 1998 .

[2]  P. Zuidema,et al.  Autocorrelated growth of tropical forest trees: Unraveling patterns and quantifying consequences , 2006 .

[3]  S. Ellner,et al.  SIZE‐SPECIFIC SENSITIVITY: APPLYING A NEW STRUCTURED POPULATION MODEL , 2000 .

[4]  Stephen P. Ellner,et al.  Stochastic stable population growth in integral projection models: theory and application , 2007, Journal of mathematical biology.

[5]  Mark Rees,et al.  Integral projection models perform better for small demographic data sets than matrix population models: a case study of two perennial herbs. , 2009 .

[6]  P. Zuidema,et al.  Lifetime growth patterns and ages of Bolivian rain forest trees obtained by tree ring analysis , 2006 .

[7]  R. Salguero‐Gómez,et al.  Keeping plant shrinkage in the demographic loop , 2010 .

[8]  John Sabo,et al.  Morris, W. F., and D. F. Doak. 2003. Quantitative Conservation Biology: Theory and Practice of Population Viability Analysis. Sinauer Associates, Sunderland, Massachusetts, USA , 2003 .

[9]  M. J. Adams.  Graph decompositions for demographic loop analysis , 2008, Journal of mathematical biology.

[10]  P. Zuidema,et al.  Do Persistently Fast‐Growing Juveniles Contribute Disproportionately to Population Growth? A New Analysis Tool for Matrix Models and Its Application to Rainforest Trees , 2009, The American Naturalist.

[11]  Stephen P. Ellner,et al.  Simple Methods for Calculating Age‐Based Life History Parameters for Stage‐Structured Populations , 1992 .

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

[13]  S. Ellner,et al.  Integral Projection Models for Species with Complex Demography , 2006, The American Naturalist.

[14]  W. Platt,et al.  DEMOGRAPHY OF A SHADE-TOLERANT TREE (FAGUS GRANDIFOLIA) IN A HURRICANE-DISTURBED FOREST , 1998 .

[15]  Kiwako S. Araki,et al.  Matrix models using fine size classes and their application to the population dynamics of tree species: Bayesian non-parametric estimation , 2007 .

[16]  R. Boot,et al.  Demography of the Brazil nut tree (Bertholletia excelsa) in the Bolivian Amazon: impact of seed extraction on recruitment and population dynamics , 2002, Journal of Tropical Ecology.

[17]  D. Peart,et al.  EARLY PERFORMANCE PREDICTS CANOPY ATTAINMENT ACROSS LIFE HISTORIES IN SUBALPINE FOREST TREES , 2005 .

[18]  P. Holgate,et al.  Matrix Population Models. , 1990 .

[19]  R. Bernal.  Demography of the vegetable ivory palm Phytelephas seemannii in Colombia, and the impact of seed harvesting , 1998 .

[20]  Paulette Bierzychudek,et al.  LOOKING BACKWARDS: ASSESSING THE PROJECTIONS OF A TRANSITION MATRIX MODEL , 1999 .

[21]  C. Pfister,et al.  INDIVIDUAL VARIATION AND ENVIRONMENTAL STOCHASTICITY: IMPLICATIONS FOR MATRIX MODEL PREDICTIONS , 2003 .

[22]  Jennifer L. Williams,et al.  The impact of invasive grasses on the population growth of Anemone patens, a long-lived native forb. , 2006, Ecology.

[23]  M. Franco,et al.  Integrating vital rate variability into perturbation analysis: an evaluation for matrix population models of six plant species , 2001 .

[24]  P. Zuidema,et al.  Conservation prospects for threatened Vietnamese tree species: results from a demographic study , 2008, Population Ecology.

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

[26]  P. Zuidema.  Demography of exploited tree species in the Bolivian Amazon , 2000 .

[27]  Carol C. Horvitz,et al.  Spatiotemporal Variation in Demographic Transitions of a Tropical Understory Herb: Projection Matrix Analysis , 1995 .

[28]  Stephen P. Ellner,et al.  Evolutionary demography of long‐lived monocarpic perennials: a time‐lagged integral projection model , 2008 .

[29]  Mark A. Burgman,et al.  Spatial Ecology of Desert Rodent Communities , 1999, Adaptations of Desert Organisms.

[30]  Shripad Tuljapurkar,et al.  A time to grow and a time to die: a new way to analyze the dynamics of size, light, age, and death of tropical trees. , 2009, Ecology.

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

[32]  P. Chien.  Demography of Threatened Tree Species in Vietnam , 2006 .

[33]  Burak Güneralp,et al.  An improved formal approach to demographic loop analysis. , 2007 .

[34]  J. Silvertown,et al.  A COMPARATIVE DEMOGRAPHY OF PLANTS BASED UPON ELASTICITIES OF VITAL RATES , 2004 .

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

[36]  Takuya Kubo,et al.  Effect of temporal autocorrelation on apparent growth rate variation in forest tree census data and an alternative distribution function of tree growth rate , 2004, Ecological Research.

[37]  P. Zuidema,et al.  Tropical tree rings reveal preferential survival of fast-growing juveniles and increased juvenile growth rates over time. , 2010, The New phytologist.

[38]  S. Ramula,et al.  Matrix dimensionality in demographic analyses of plants : when to use smaller matrices? , 2005 .

[39]  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.

[40]  D. Doak,et al.  Book Review: Quantitative Conservation biology: Theory and Practice of Population Viability analysis , 2004, Landscape Ecology.