Hierarchical demographic approaches for assessing invasion dynamics of non-indigenous species: An example using northern snakehead (Channa argus)

Models of species’ demographic features are commonly used to understand population dynamics and inform management tactics. Hierarchical demographic models are ideal for the assessment of non-indigenous species because our knowledge of non-indigenous populations is usually limited, data on demographic traits often come from a species’ native range, these traits vary among populations, and traits are likely to vary considerably over time as species adapt to new environments. Hierarchical models readily incorporate this spatiotemporal variation in species’ demographic traits by representing demographic parameters as multi-level hierarchies. As is done for traditional non-hierarchical matrix models, sensitivity and elasticity analyses are used to evaluate the contributions of different life stages and parameters to estimates of population growth rate. We applied a hierarchical model to northern snakehead (Channa argus), a fish currently invading the eastern United States. We used a Monte Carlo approach to simulate uncertainties in the sensitivity and elasticity analyses and to project future population persistence under selected management tactics. We gathered key biological information on northern snakehead natural mortality, maturity and recruitment in its native Asian environment. We compared the model performance with and without hierarchy of parameters. Our results suggest that ignoring the hierarchy of parameters in demographic models may result in poor estimates of population size and growth and may lead to erroneous management advice. In our case, the hierarchy used multi-level distributions to simulate the heterogeneity of demographic parameters across different locations or situations. The probability that the northern snakehead population will increase and harm the native fauna is considerable. Our elasticity and prognostic analyses showed that intensive control efforts immediately prior to spawning and/or juvenile-dispersal periods would be more effective (and probably require less effort) than year-round control efforts. Our study demonstrates the importance of considering the hierarchy of parameters in estimating population growth rate and evaluating different management strategies for non-indigenous invasive species.     

Temporal limits to simulating the future spread pattern of invasive species: Buddleja davidii in Europe and New Zealand

The spread of invasive species is a major ecological and economic problem. Dynamic spread modelling is a potentially valuable tool to assist regional and central government authorities to monitor and control invasive species. To date a lack of suitable data has meant that most broad scale dispersal models have not been validated with independent datasets, and so their predictive ability and reliability has remained unscrutinised. A dynamic, stochastic dispersal model of the widely invasive plant Buddleja davidii was calibrated on European spread data and then used to project the temporal progression of B. davidii's distribution in New Zealand, starting from several different historical distributions. To assess the model's performance, we constructed an occupancy map based on the average number of simulation realisations that have a population present. The application of Receiver Operating Characteristic (ROC) curves to occupancy maps is introduced, but with specificity substituted by the proportion of available area used in a realisation. A derivative measure, the partial area under these curves when assessed through time (pAUC), is introduced and used to assess overall performance of the spread model. The model was able to attain a high level of model sensitivity, encompassing all of the known locations within the occupancy envelope. However, attempting to simulate the spread of this invasive species beyond a decade had very low model specificity. This is due to several factors, including the exponential process of spread (the further a population spreads the more sites exist from which it can spread stochastically), and the Markovian chain property of the stochastic system whereby differences between realisations compound through time. These features are seen in many reports of spread models, without being explicitly acknowledged. Our measure of pAUC through time allows a model's temporal performance and its specificity to be simultaneously assessed. While the rapid deterioration in model performance limits the utility of this type of modelling for forecasting long-term broad-scale strategic management of biological invasions, it does not necessarily limit its attractiveness for informing smaller scale and shorter term invasion management activities such as surveillance, containment and local eradication.     

Unveiling human-assisted dispersal mechanisms in invasive alien insects: Integration of spatial stochastic simulation and phenology models

Capturing the spread of biological invasions in heterogeneous landscapes is a complex modelling task where information on both dispersal and population dynamics needs to be integrated. Spatial stochastic simulation and phenology models have rarely been combined to assist in the study of human-assisted long-distance dispersal events.Here we develop a process-based spatially explicit landscape-extent simulation model that considers the spread and detection of invasive insects. Natural and human-assisted dispersal mechanisms are modelled with an individual-based approach using negative exponential and negative power law dispersal kernels and gravity models. The model incorporates a phenology sub-model that uses daily temperature grids for the prediction and timing of the population dynamics in each habitat patch. The model was applied to the study of the invasion by the important maize pest western corn rootworm (WCR) Diabrotica virgifera ssp. virgifera in Europe. We parameterized and validated the model using maximum likelihood and simulation methods from the historical invasion of WCR in Austria.WCR was found to follow stratified dispersal where international transport networks in the Danube basin played a key role in the occurrence of long-distance dispersal events. Detection measures were found to be effective and altitude had a significant effect on limiting the spread of WCR. Spatial stochastic simulation combined with phenology models, maximum likelihood methods and predicted versus observed regression showed a high degree of flexibility that captured the salient features of WCR spread in Austria. This modelling approach is useful because it allows to fully exploit and the often limited and heterogeneous information available regarding the population dynamics and dispersal of alien invasive insects.     

Longevity can buffer plant and animal populations against changing climatic variability

Both means and year-to-year variances of climate variables such as temperature and precipitation are predicted to change. However, the potential impact of changing climatic variability on the fate of populations has been largely unexamined. We analyzed multiyear demographic data for 36 plant and animal species with a broad range of life histories and types of environment to ask how sensitive their long-term stochastic population growth rates are likely to be to changes in the means and standard deviations of vital rates (survival, reproduction, growth) in response to changing climate. We quantified responsiveness using elasticities of the long-term population growth rate predicted by stochastic projection matrix models. Short-lived species (insects and annual plants and algae) are predicted to be more strongly (and negatively) affected by increasing vital rate variability relative to longer-lived species (perennial plants, birds, ungulates). Taxonomic affiliation has little power to explain sensitivity to increasing variability once longevity has been taken into account. Our results highlight the potential vulnerability of short-lived species to an increasingly variable climate, but also suggest that problems associated with short-lived undesirable species (agricultural pests, disease vectors, invasive weedy plants) may be exacerbated in regions where climate variability decreases.     

Climate change threatens polar bear populations: a stochastic demographic analysis

The polar bear (Ursus maritimus) depends on sea ice for feeding, breeding, and movement. Significant reductions in Arctic sea ice are forecast to continue because of climate warming. We evaluated the impacts of climate change on polar bears in the southern Beaufort Sea by means of a demographic analysis, combining deterministic, stochastic, environment-dependent matrix population models with forecasts of future sea ice conditions from IPCC general circulation models (GCMs). The matrix population models classified individuals by age and breeding status; mothers and dependent cubs were treated as units. Parameter estimates were obtained from a capture-recapture study conducted from 2001 to 2006. Candidate statistical models allowed vital rates to vary with time and as functions of a sea ice covariate. Model averaging was used to produce the vital rate estimates, and a parametric bootstrap procedure was used to quantify model selection and parameter estimation uncertainty. Deterministic models projected population growth in years with more extensive ice coverage (2001-2003) and population decline in years with less ice coverage (2004-2005). LTRE (life table response experiment) analysis showed that the reduction in lambda in years with low sea ice was due primarily to reduced adult female survival, and secondarily to reduced breeding. A stochastic model with two environmental states, good and poor sea ice conditions, projected a declining stochastic growth rate, log lambdas, as the frequency of poor ice years increased. The observed frequency of poor ice years since 1979 would imply log lambdas approximately - 0.01, which agrees with available (albeit crude) observations of population size. The stochastic model was linked to a set of 10 GCMs compiled by the IPCC; the models were chosen for their ability to reproduce historical observations of sea ice and were forced with "business as usual" (A1B) greenhouse gas emissions. The resulting stochastic population projections showed drastic declines in the polar bear population by the end of the 21st century. These projections were instrumental in the decision to list the polar bear as a threatened species under the U.S. Endangered Species Act.     

Predicting [corrected] extinction risk in spite of predator-prey oscillations.

Most population viability analyses (PVA) assume that the effects of species interactions are subsumed by population-level parameters. We examine how robust five commonly used PVA models are to violations of this assumption. We develop a stochastic, stage-structured predator-prey model and simulate prey population vital rates and abundance. We then use simulated data to parameterize and estimate risk for three demographic models (static projection matrix, stochastic projection matrix, stochastic vital rate matrix) and two time series models (diffusion approximation [DA], corrupted diffusion approximation [CDA]). Model bias is measured as the absolute deviation between estimated and observed quasi-extinction risk. Our results highlight three generalities about the application of single-species models to multi-species conservation problems. First, our collective model results suggest that most single-species PVA models overestimate extinction risk when species interactions cause periodic variation in abundance. Second, the DA model produces the most (conservatively) biased risk forecasts. Finally, the CDA model is the most robust PVA to population cycles caused by species interactions. CDA models produce virtually unbiased and relatively precise risk estimates even when populations cycle strongly. High performance of simple time series models like the CDA owes to their ability to effectively partition stochastic and deterministic sources of variation in population abundance.     

How sensitive are elasticities of long-run stochastic growth to how environmental variability is modelled?

A variable environment leaves a signature in a population's dynamics. Deriving statistical and mathematical models of how environmental variability affects population projections has - in the wake of reports of substantial climatic fluctuations - received much recent attention. If the model changes, then so too does the population projection. This is because a different model of environmental variability changes estimates of long-run stochastic growth, which is a function of demographic rates and their temporal sequence. Decomposing elasticities of long-run stochastic growth into constituent parts can assess the relative influence of different components. Here, we investigate the consequences of changing the environmental state definition, and therefore altering the shape of demographic rate distributions and their temporal sequence, by using age-structured matrix models to project vertebrate populations into the future under a range of environmental scenarios. The identity of the most influential demographic rate was consistent among all approaches that perturbed only the mean, but was not when only the variance was perturbed. Furthermore, the influence of each demographic rate fluctuated among projections by up to factors of six and two for changes to the variance and mean, respectively. These changes in influence depend in part upon how environmental variability - in particular, the color of environmental noise - is incorporated. In the light of predictions of increasing climatic variability in the future, these results suggest caution when drawing quantitative conclusions from stochastic population projections.     

Invasion, Disturbance, and Competition: Modeling the Fate of Coastal Plant Populations

Abstract:  Wetland habitats are besieged by biotic and abiotic disturbances such as invasive species, hurricanes, habitat fragmentation, and salinization. Predicting how these factors will alter local population dynamics and community structure is a monumental challenge. By examining ecologically similar congeners, such as Iris hexagona and I. pseudacorus (which reproduce clonally and sexually and tolerate a wide range of environmental conditions), one can identify life-history traits that are most influential to population growth and viability. We combined empirical data and stage-structured matrix models to investigate the demographic responses of native ( I. hexagona ) and invasive ( I. pseudacorus ) plant populations to hurricanes and salinity stress in freshwater and brackish wetlands. In our models I. hexagona and I. pseudacorus responded differently to salinity stress, and species coexistence was rare. In 82% of computer simulations of freshwater marsh, invasive iris populations excluded the native species within 50 years, whereas native populations excluded the invasive species in 99% of the simulations in brackish marsh. The occurrence of hurricanes allowed the species to coexist, and species persistence was determined by the length of time it took the ecosystem to recover. Rapid recovery (2 years) favored the invasive species, whereas gradual recovery (30 years) favored the native species. Little is known about the effects of hurricanes on competitive interactions between native and invasive plant species in marsh ecosystems. Our models contribute new insight into the relationship between environmental disturbance and invasion and demonstrate how influential abiotic factors such as climate change will be in determining interspecific interactions.     

Exploring metapopulation-scale suppression alternatives for a global invader in a river network experiencing climate change

Invasive species can dramatically alter ecosystems, but eradication is difficult, and suppression is expensive once they are established. Uncertainties in the potential for expansion and impacts by an invader can lead to delayed and inadequate suppression, allowing for establishment. Metapopulation viability models can aid in planning strategies to improve responses to invaders and lessen invasive species’ impacts, which may be particularly important under climate change. We used a spatially explicit metapopulation viability model to explore suppression strategies for ecologically damaging invasive brown trout (Salmo trutta), established in the Colorado River and a tributary in Grand Canyon National Park. Our goals were to estimate the effectiveness of strategies targeting different life stages and subpopulations within a metapopulation; quantify the effectiveness of a rapid response to a new invasion relative to delaying action until establishment; and estimate whether future hydrology and temperature regimes related to climate change and reservoir management affect metapopulation viability and alter the optimal management response. Our models included scenarios targeting different life stages with spatially varying intensities of electrofishing, redd destruction, incentivized angler harvest, piscicides, and a weir. Quasi-extinction (QE) was obtainable only with metapopulation-wide suppression targeting multiple life stages. Brown trout population growth rates were most sensitive to changes in age 0 and large adult mortality. The duration of suppression needed to reach QE for a large established subpopulation was 12 years compared with 4 with a rapid response to a new invasion. Isolated subpopulations were vulnerable to suppression; however, connected tributary subpopulations enhanced metapopulation persistence by serving as climate refuges. Water shortages driving changes in reservoir storage and subsequent warming would cause brown trout declines, but metapopulation QE was achieved only through refocusing and increasing suppression. Our modeling approach improves understanding of invasive brown trout metapopulation dynamics, which could lead to more focused and effective invasive species suppression strategies and, ultimately, maintenance of populations of endemic fishes.     

A hierarchical Bayesian non-linear spatio-temporal model for the spread of invasive species with application to the Eurasian Collared-Dove

The spread of invasive species is a long studied subject that garners much interest in the ecological research community. Historically the phenomenon has been approached using a purely deterministic mathematical framework (usually involving differential equations of some form). These methods, while scientifically meaningful, are generally highly simplified and fail to account for uncertainty in the data and process, of which our knowledge could not possibly exist without error. We propose a hierarchical Bayesian model for population spread that accommodates data sources with errors, dependence structures between population dynamics parameters, and takes into account prior scientific understanding via non-linear relationships between model parameters and space-time response variables. We model the process (i.e., the bird population in this case) as a Poisson response with spatially varying diffusion coefficients as well as a logistic population growth term using a common reaction-diffusion equation that realistically mimics the ecological process. We focus the application on the ongoing invasion of the Eurasian Collared-Dove.     

Relating chronic toxicity responses to population-level effects: A comparison of population-level parameters for three salmon species as a function of low-level toxicity

Standard laboratory toxicity tests assess the physiological responses of individual organisms to exposure to toxic substances under controlled conditions. Time and space restrictions often prevent the assessment of population-level responses to a toxic substance. Contaminants can affect various biological functions (e.g. growth, fecundity or behavior), which may alter different demographic traits, leading to population-level impacts. In this study, immune suppression, reproductive dysfunction and somatic growth impairment were examined using life history matrix models for coho salmon (Oncorhynchus kisutch), sockeye salmon (Oncorhynchus nerka) and chinook salmon (Oncorhynchus tshawytscha). Our intent was to gauge the relative magnitude of response to toxic effects among species and between life history stages, not provide a specific estimate of population growth rate or abundance. Effects due to immune suppression were modeled as reductions in age-specific survival. Toxic impacts on reproductive function were modeled as a 10% reduction in reproductive contribution for all reproductively mature age groups. Model runs that examined the effect of somatic growth reduction on population parameters incorporated both survival and reproductive impacts. All impacts were modeled as 10% reductions in the affected population demographic parameters. First-year survival and reproductive impacts produced similar population growth rates (λ), but resulted in different sensitivity and stable age distributions. Modeled somatic growth reduction produced additive effects on survival and reproduction. Toxic stressors producing similar changes in λ did not necessarily produce similar changes in the age distributions. Sensitivity and elasticity analyses demonstrated that changes to the first-year survival rate produced the greatest per-unit effect on λ for each species. Alteration in abundance of mature females also strongly influenced λ. Differences observed between species showed that the number of reproductive ages and time to reproductive maturity were important components for population-level responses. These results emphasize the importance of linking toxicity responses at low concentrations to the demographic traits they affect, and help to highlight the toxicity tests that are more suitable for assessing impacts on the focal species. Additionally, life history modeling is a useful tool for developing testable hypotheses regarding impacts on specific populations as well as for conducting comparisons between populations.     

Uncertainty analysis of transient population dynamics

Two types of demographic analyses, perturbation analysis and uncertainty analysis, can be conducted to gain insights about matrix population models and guide population management. Perturbation analysis studies how the perturbation of demographic parameters (survival, growth, and reproduction parameters) may affect the population projection, while uncertainty analysis evaluates how much uncertainty there is in population dynamic predictions and where the uncertainty comes from. Previously, both perturbation analysis and uncertainty analysis were conducted on the long-term population growth rate. However, the population may not reach its equilibrium state, especially when there is management by harvesting or hunting. Recently, there has been an increased interest in short-term transient dynamics, which can differ from asymptotic long-term dynamics. There are currently techniques to conduct perturbation analyses of short-term transient dynamics, but no techniques have been proposed for uncertainty analysis of such dynamics. In this study, we introduced an uncertainty analysis technique, the general Fourier Amplitude Sensitivity Test (FAST), to study uncertainties in transient population dynamics. The general FAST is able to identify the amount of uncertainty in transient dynamics and contributions by different demographic parameters. We applied the general FAST to a mountain goat (Oreamnos americanus) matrix population model to give a clear illustration of how uncertainty analysis can be conducted for transient dynamics arising from matrix population models.     

Effects of temperature variation on TSD in turtle (C. picta) populations

We formulate a two-sex model of temperature-dependent sex determination (TSD) for a freshwater turtle (C. picta) population. The aim is to understand how environmental temperature variations and nest heat conduction properties affect the long term dynamics of the population. This is a key to understanding how global temperature changes may affect their survival. With stochastic inputs of ambient temperature and solar radiation, the model uses the heat equation to determine the temperature in the egg layer in the nest; in turn, this determines the sex ratio in the egg clutch using a variable degree-day model. Finally, a nonlinear Leslie type, stage-based, two-sex model, is used to determine the long term male and female populations. A two-sex model is required because of different development rates for males and females. The model is flexible enough to enable other researchers to examine the effects of temperature variation variations on other species with TSD, e.g., crocodilians, reptilians, as well as other turtle species. It can be adapted to study effects of nest location, soil type, rain events, different incubation periods, and density effects, for example, the dependence of the mating function on the ratio of males to females and each’s contribution to the sex of hatchlings. Modifications can be easily made to fit a specific life history traits. The model is a beginning step in understanding the long term, high fitness shown by many reptile species with TSD, and it may suggest to experimentalists what data may be relevant to these issues; it can also be useful to wildlife managers in developing strategies for intervention if needed. Among the principal findings are that temperature variability and detailed nest heat conduction properties may buffer projected negative effects on a population.     

Joint modelling of breeding and survival in the kittiwake using frailty models

Assessment of population dynamics is central to population dynamics and conservation. In structured populations, matrix population models based on demographic data have been widely used to assess such dynamics. Although highlighted in several studies, the influence of heterogeneity among individuals in demographic parameters and of the possible correlation among these parameters has usually been ignored, mostly because of difficulties in estimating such individual-specific parameters. In the kittiwake (Rissa tridactyla), a long-lived seabird species, differences in survival and breeding probabilities among individual birds are well documented. Several approaches have been used in the animal ecology literature to establish the association between survival and breeding rates. However, most are based on observed heterogeneity between groups of individuals, an approach that seldom accounts for individual heterogeneity. Few attempts have been made to build models permitting estimation of the correlation between vital rates. For example, survival and breeding probability of individual birds were jointly modelled using logistic random effects models by [Cam, E., Link, W.A., Cooch, E.G., Monnat, J., Danchin, E., 2002. Individual covariation in life-history traits: seeing the trees despite the forest. Am. Naturalist, 159, in press]. This is the only example in wildlife animal populations we are aware of. Here we adopt the survival analysis approaches from epidemiology. We model the survival and the breeding probability jointly using a normally distributed random effect (frailty). Conditionally on this random effect, the survival time is modelled assuming a lognormal distribution, and breeding is modelled with a logistic model. Since the deaths are observed in year-intervals, we also take into account that the data are interval censored. The joint model is estimated using classic frequentist methods and also MCMC techniques in Winbugs. The association between survival and breeding attempt is quantified using the standard deviation of the random frailty parameters. We apply our joint model on a large data set of 862 birds, that was followed from 1984 to 1995 in Brittany (France). Survival is positively correlated with breeding indicating that birds with greater inclination to breed also had higher survival.     

Extinction-effective population index: incorporating life-history variations in population viability analysis

Viability status of populations is a commonly used measure for decision-making in the management of populations. One of the challenges faced by managers is the need to consistently allocate management effort among populations. This allocation should in part be based on comparison of extinction risks among populations. Unfortunately, common criteria that use minimum viable population size or count-based population viability analysis (PVA) often do not provide results that are comparable among populations, primarily because they lack consistency in determining population size measures and threshold levels of population size (e.g., minimum viable population size and quasi-extinction threshold). Here I introduce a new index called the "extinction-effective population index," which accounts for differential effects of demographic stochasticity among organisms with different life-history strategies and among individuals in different life stages. This index is expected to become a new way of determining minimum viable population size criteria and also complement the count-based PVA. The index accounts for the difference in life-history strategies of organisms, which are modeled using matrix population models. The extinction-effective population index, sensitivity, and elasticity are demonstrated in three species of Pacific salmonids. The interpretation of the index is also provided by comparing them with existing demographic indices. Finally, a measure of life-history-specific effect of demographic stochasticity is derived.     

Effects of Harvest of Nontimber Forest Products and Ecological Differences between Sites on the Demography of African Mahogany

Abstract: The demographic impacts of harvesting nontimber forest products (NTFP) have been increasingly studied because of reports of potentially unsustainable harvest. Nevertheless, our understanding of how plant demographic response to harvest is altered by variation in ecological conditions, which is critical for developing realistic sustainable‐use plans, is limited. We built matrix population models to test whether and how variation in ecological conditions affects population responses to harvest. In particular, we examined the effect of bark and foliage harvest on the demography of populations of African mahogany (Khaya senegalensis) in two contrasting ecological regions of Benin, West Africa. K. senegalensis bark and foliage harvest significantly reduced its stochastic population growth rates, but ecological differences between regions had a greater effect on population growth rates than did harvest. The effect of harvest on population growth rates (Δλ) was slightly stronger in the moist than in the drier region. Life‐table response experiments revealed that the mechanism by which harvesting reduced λ differed between ecological regions. Lowered stasis (persistence) of larger life stages lead to a reduction in λ in the drier region, whereas lowered growth of all life stages lowered λ in moist region. Potential strategies to increase population growth rates should include decreasing the proportion of individuals harvested, promoting harvester‐owned plantations of African mahogany, and increasing survival and growth by promoting no‐fire zones in gallery forests. Our results show how population responses to harvest of NTFP may be altered by ecological differences across sites and emphasize the importance of monitoring populations over the climatic range in which they occur to develop more realistic recommendations for conservation.     

Estimating the functional form for the density dependence from life history data

Two contrasting approaches to the analysis of population dynamics are currently popular: demographic approaches where the associations between demographic rates and statistics summarizing the population dynamics are identified; and time series approaches where the associations between population dynamics, population density, and environmental covariates are investigated. In this paper, we develop an approach to combine these methods and apply it to detailed data from Soay sheep (Ovis aries). We examine how density dependence and climate contribute to fluctuations in population size via age- and sex-specific demographic rates, and how fluctuations in demographic structure influence population dynamics. Density dependence contributes most, followed by climatic variation, age structure fluctuations and interactions between density and climate. We then simplify the density-dependent, stochastic, age-structured demographic model and derive a new phenomenological time series which captures the dynamics better than previously selected functions. The simple method we develop has potential to provide substantial insight into the relative contributions of population and individual-level processes to the dynamics of populations in stochastic environments.     

A stage-structured, Aedes albopictus population model

Aedes albopictus has been the fastest spreading invasive animal species in the world from the mid-1980s until the mid-2000s. In areas it infests, it disrupts native mosquito ecology and can potentially vector up to 21 viruses. To better understand the population dynamics of this species, we created a temperature dependent population model. A stage-structured model was chosen to allow each life-stage to have different temperature dependent mortality and development rates, and each stage was modeled with an ordinary differential equation. Model parameters and distributions were based upon literature values. Initially, a basic model was constructed. This model then had parameters that were forced based upon daily average temperatures. Several criteria were used to evaluate the model, including a comparison to field data from Lubbock, TX. In a stochastic version of the model, a 95% confidence limit contained 70.7% of the field data points. Based upon these results, we feel reasonably confident that we have captured the role of temperature in driving the population dynamics of Ae. albopictus.     

A mathematical model of field population dynamics of the southern pine beetle, dendroctonus frontalis

The southern pine beetle, Dendroctonus frontalis Zimmermann, is a major insect pest of pines in the southern U.S.A. A population dynamic model of D. frontalis within an infestation was developed. The model is based on the integration of several component submodels that describe individual life stage processes of the beetle. The inclusion of variable forest stand density and microclimate conditions results in a model response which closely follows actual field responses. A sensitivity analysis was performed to investigate effects of weather conditions and the importance of the spatial patterns within the forest.     

Modeling the population dynamics of capybara Hydrochaeris hydrochaeris: a first step towards a management plan

Capybara (Hydrochaeris hydrochaeris) is the largest living rodent, widely distributed in South America, with a considerable potential as an economic resource. There are large populations in the region of Esteros del Ibera that could be exploited in a sustainable manner and contribute positively to the development of the region. Such exploitation requires to be done responsibly and following a specially designed management plan. In the particular case of wildlife management, a sustainable exploitation of a given species requires knowledge of its dynamics, its density in the area, and its vital rates, as well as understanding how the natural environment will respond to the proposed manipulation. It is also important to determine which environmental variables lead the population dynamics in each region. A mathematical model is a very convenient tool for analyzing possible strategies.This work presents a matrix population model structured in five stages. Sensitivity and elasticity analyses of the model provide an understanding of the effects of differentiated class-dependent survival and growth probabilities on the life cycle. A harvest term has been included in the model with the purpose of evaluating the sustainability of a given exploitation scenario or comparing different harvest strategies. The different strategies depend on the exploitation pattern, the exploitation intensity, and the season in which the harvest takes place.The purpose is to formulate a simple tool that can be understood by a manager without specific training. Hence, the software interface was specially designed for a user with elementary knowledge of vital parameters of the species.