Disentangling the effects of heterogeneity, stochastic dynamics and sampling in a community of aquatic insects

A dynamic and heterogeneous species abundance model generating the lognormal species abundance distribution is fitted to time series of species data from an assemblage of stoneflies and mayflies (Plecoptera and Ephemeroptera) of an aquatic insect community collected over a period of 15 years. In each year except one, we analyze 5 parallel samples taken at the same time of the season giving information about the over-dispersion in the sampling relative to the Poisson distribution. Results are derived from a correlation analysis, where the correlation in the bivariate normal distribution of log abundance is used as measurement of similarity between communities. The analysis enables decomposition of the variance of the lognormal species abundance distribution into three components due to heterogeneity among species, stochastic dynamics driven by environmental noise, and over-dispersion in sampling, accounting for 62.9, 30.6 and 6.5% of the total variance, respectively. Corrected for sampling the heterogeneity and stochastic components accordingly account for 67.3 and 32.7% of the among species variance in log abundance. By using this method, it is possible to disentangle the effect of heterogeneity and stochastic dynamics by quantifying these components and correctly remove sampling effects on the observed species abundance distribution.     

Incorporating low-resolution historic species location data decreases performance of distribution models

Developing robust species distribution models is important as model outputs are increasingly being incorporated into conservation policy and management decisions. A largely overlooked component of model assessment and refinement is whether to include historic species occurrence data in distribution models to increase the data sample size. Data of different temporal provenance often differ in spatial accuracy and precision. We test the effect of inclusion of historic coarse-resolution occurrence data on distribution model outputs for 187 species of birds in Australian tropical savannas. Models using only recent (after 1990), fine-resolution data had significantly higher model performance scores measured with area under the receiver operating characteristic curve (AUC) than models incorporating both fine- and coarse-resolution data. The drop in AUC score is positively correlated with the total area predicted to be suitable for the species (R² = 0.163-0.187, depending on the environmental predictors in the model), as coarser data generally leads to greater predicted areas. The remaining unexplained variation is likely to be due to the covariate errors resulting from resolution mismatch between species records and environmental predictors. We conclude that decisions regarding data use in species distribution models must be conscious of the variation in predictions that mixed-scale datasets might cause.     

Spatial statistics for modeling of abundance and distribution of wildlife species in the Masai Mara ecosystem, Kenya

This study illustrates the use of modern statistical procedures for better wildlife management by addressing three key issues: determination of abundance, modeling of animal distributions and variability of diversity in space and time. Prior information in Markov Chain Monte Carlo (MCMC) methods is used to improve estimates of abundance. Measures of autocorrelation are included when modeling distributions of animal counts, and a diversity index to indicate species abundance and richness for large herbivores is developed. Data from the Masai Mara ecosystem in Kenya are used to develop and demonstrate these procedures. The new abundance estimates are up to 35% more accurate than those obtained by existing methods. Significant temporal changes in spatial patterns are found from a space-time analysis of elephant counts over a 20-year period, with strong interactions over 5 km and 6 months space and time separations, respectively. The new diversity index is sensitive to both high abundance and species richness and is also able to capture year to year variation. It indicates an overall marginal decrease in diversity for large herbivores in the Mara ecosystem. The space-time analyses and diversity index can easily be computed thereby providing tools for rapid decision making.     

Modelling the abundance of three key plant species in New Zealand hill-pasture using a decision tree approach

Decision tree models were developed to investigate and predict the relative abundance of three key pasture plants [ryegrass (Lolium perenne), browntop (Agrostis capillaris), and white clover (Trifolium repens)] with integration of a geographical information system (GIS) in a naturalised hill-pasture in the North Island, New Zealand, and were compared with regression models with respect to model fit and predictive accuracy. The results indicated that the decision tree models had a better model fit in terms of average squared error (ASE) and a higher percentage of adequately predicted cases in model validation than the corresponding regression models. These decision tree models clearly revealed the relative importance of environmental and management variables in influencing the abundance of these three species. Hill slope was the most significant environmental factor influencing the abundance of ryegrass while soil Olsen P and annual P fertilizer input were the most significant factors influencing the abundance of browntop, and white clover, respectively. Soil Olsen P of approximately 10 μg/g, or a slope of about 10.5° was critical points where the competition between ryegrass and browntop tended to come to an equilibrium. Integrating the decision tree models with a GIS in this study not only facilitated the model development and analyses, but also provided a useful decision support tool in pasture management such as in assisting precision fertilizer placement. The insights obtained from the decision tree models also have important implications for pasture management, for example, it is important to maintain a soil Olsen P higher than 10 μg/g in order to keep the dominance of ryegrass in the hill-pasture.     

南亚热带森林3层次物种多度之广义泊松分布模型

生物群落中物种多度分布(species abundance distribution)呈典型的倒J形,即其中存在许多稀有种、少量常见种.物种多度分布模型研究有助于解决森林生态恢复中的物种配置等实际问题.本研究考察了一种过分散(over-dispersion,或称超分布,即方差大于均值)的离散型分布,即具有λ和α两个参数...     

Estimating marine biogeochemical rates of the carbonate pH system—A Kalman filter tested

Oxygen (O₂), nitrate (NO₃), dissolved inorganic carbon (DIC) or pCO₂, and pH or total alkalinity (TA), are useful indices of marine chemical, physical and biological processes operating on varying time-scales. Although these properties are increasingly being monitored at high frequency, they have not been extensively used for studying ecosystem dynamics. We test whether we can estimate time-evolving biogeochemical rates (e.g. primary production, respiration, calcification and carbonate dissolution, and nitrification) from synthetic high frequency time-series of O₂, NO₃, DIC, pCO₂, TA or pH. More specifically, a Kalman filter has been implemented in a very simplified biogeochemical model describing the dynamics of O₂, NO₃, DIC and TA and linking the concentration data to biogeochemical fluxes. Different sets of concentration data are assimilated and biogeochemical rates are estimated. The frequency of assimilation required to get acceptable results is investigated and is compared with the frequency of sampling in the field or in controlled experimental settings.Smoothing of the data to remove data noise before assimilation improves the estimation of the biogeochemical rates. The best estimated rates are obtained when assimilating O₂, NO₃ and TA although the assimilation of DIC instead of TA also gives satisfactory results. In case pH or pCO₂ is assimilated rather than DIC or TA, the linearization of the (now nonlinear) observation equation introduces perturbations and the Kalman filter behaves suboptimal. We conclude that, given the resolution of data required, the tool has potential to estimate biogeochemical rates of the carbonate system under controlled settings.     

Bimodality of the combined removal and signs-of-activities estimator for sampling closed animal populations

The possibility of a bimodal log-likelihood function arises with certain data when the combined removal and signs-of-activities estimator is used. Bimodal log-likelihoods may, in turn, yield disjoint confidence intervals for certain confidence levels. The hypothesis that bimodality is caused by the violation of the equal catchability assumption of the removal model, leading to the combination of contradictory data/models in the combined estimator is set forth. Simulations exploring the effect of the violation of removal model assumptions on estimation and inference showed that the assumption of unequal capture probability influenced the frequency of bimodal likelihoods; similarly, extreme parameter values for probability of capture influenced the number of excessively large confidence intervals produced. A sex-specific combined estimator is developed as a remedial model tailored to the problem. The simulations suggest that both the signs-of-activities estimator and the sex-specific estimator perform equally well over the range of simulations presented, though the signs-of-activities estimator is easier to implement.     

Assessing the effects of disease and bleaching on Florida Keys corals by fitting population models to data

Coral diseases have increased in frequency over the past few decades and have important influences on the structure and composition of coral reef communities. However, there is limited information on the etiologies of many coral diseases, and pathways through which coral diseases are acquired and transmitted are still in question. Furthermore, it is difficult to assess the impacts of disease on coral populations because outbreaks often co-occur with temperature-induced bleaching and anthropogenic stressors. We developed spatially explicit population models of coral disease and bleaching dynamics to quantify the impact of six common diseases on Florida Keys corals, including aspergillosis, dark spots, white band, white plague, white patch, and Caribbean yellow band. Models were fit to an 8-year data set of coral abundance, disease prevalence, and bleaching prevalence. Model selection was used to assess alternative pathways for disease transmission, and the influence of environmental stressors, including sea temperature and human population density, on disease prevalence and coral mortality. Classic disease transmission from contagious to susceptible colonies provided the best-fit model only for aspergillosis. For other diseases, external disease forcing, such as through a vector or directly from pathogens in the environment, provided the best fit to observed data. Estimates of disease reproductive ratio values (R₀) were less than one for each disease, indicating coral colonies were below densities required for diseases to become established through contagious spread alone. Incidences of white band and white patch disease were associated with greater susceptibility or slower recovery of bleached colonies, and no disease outbreaks were associated with periods of elevated sea temperatures alone. Projections of best-fit models indicated that, atleast during the period of this study, disease and bleaching did not have substantial impacts on populations and impaired rates of population growth appeared to be attributable to other stressors. By applying epidemiological models to field data, our study gives qualitative insights into the dynamics of coral diseases, relative stressor impacts, and directions for future research.     

Application of an evolutionary algorithm to the inverse parameter estimation of an individual-based model

Inverse parameter estimation of individual-based models (IBMs) is a research area which is still in its infancy, in a context where conventional statistical methods are not well suited to confront this type of models with data. In this paper, we propose an original evolutionary algorithm which is designed for the calibration of complex IBMs, i.e. characterized by high stochasticity, parameter uncertainty and numerous non-linear interactions between parameters and model output. Our algorithm corresponds to a variant of the population-based incremental learning (PBIL) genetic algorithm, with a specific “optimal individual” operator. The method is presented in detail and applied to the individual-based model OSMOSE. The performance of the algorithm is evaluated and estimated parameters are compared with an independent manual calibration. The results show that automated and convergent methods for inverse parameter estimation are a significant improvement to existing ad hoc methods for the calibration of IBMs.     

Assessing the effects of pseudo-absences on predictive distribution model performance

Modelling species distributions with presence data from atlases, museum collections and databases is challenging. In this paper, we compare seven procedures to generate pseudo-absence data, which in turn are used to generate GLM-logistic regressed models when reliable absence data are not available. We use pseudo-absences selected randomly or by means of presence-only methods (ENFA and MDE) to model the distribution of a threatened endemic Iberian moth species (Graellsia isabelae). The results show that the pseudo-absence selection method greatly influences the percentage of explained variability, the scores of the accuracy measures and, most importantly, the degree of constraint in the distribution estimated. As we extract pseudo-absences from environmental regions further from the optimum established by presence data, the models generated obtain better accuracy scores, and over-prediction increases. When variables other than environmental ones influence the distribution of the species (i.e., non-equilibrium state) and precise information on absences is non-existent, the random selection of pseudo-absences or their selection from environmental localities similar to those of species presence data generates the most constrained predictive distribution maps, because pseudo-absences can be located within environmentally suitable areas. This study shows that if we do not have reliable absence data, the method of pseudo-absence selection strongly conditions the obtained model, generating different model predictions in the gradient between potential and realized distributions.     

Dynamics of a low‐density tiger population in Southeast Asia in the context of improved law enforcement

Recovering small populations of threatened species is an important global conservation strategy. Monitoring the anticipated recovery, however, often relies on uncertain abundance indices rather than on rigorous demographic estimates. To counter the severe threat from poaching of wild tigers (Panthera tigris), the Government of Thailand established an intensive patrolling system in 2005 to protect and recover its largest source population in Huai Kha Khaeng Wildlife Sanctuary. Concurrently, we assessed the dynamics of this tiger population over the next 8 years with rigorous photographic capture‐recapture methods. From 2006 to 2012, we sampled across 624–1026 km² with 137–200 camera traps. Cameras deployed for 21,359 trap days yielded photographic records of 90 distinct individuals. We used closed model Bayesian spatial capture‐recapture methods to estimate tiger abundances annually. Abundance estimates were integrated with likelihood‐based open model analyses to estimate rates of annual and overall rates of survival, recruitment, and changes in abundance. Estimates of demographic parameters fluctuated widely: annual density ranged from 1.25 to 2.01 tigers/100 km², abundance from 35 to 58 tigers, survival from 79.6% to 95.5%, and annual recruitment from 0 to 25 tigers. The number of distinct individuals photographed demonstrates the value of photographic capture–recapture methods for assessments of population dynamics in rare and elusive species that are identifiable from natural markings. Possibly because of poaching pressure, overall tiger densities at Huai Kha Khaeng were 82–90% lower than in ecologically comparable sites in India. However, intensified patrolling after 2006 appeared to reduce poaching and was correlated with marginal improvement in tiger survival and recruitment. Our results suggest that population recovery of low‐density tiger populations may be slower than anticipated by current global strategies aimed at doubling the number of wild tigers in a decade.     

Predicting Recovery Criteria for Threatened and Endangered Plant Species on the Basis of Past Abundances and Biological Traits

Recovery plans for species listed under the U.S. Endangered Species Act are required to specify measurable criteria that can be used to determine when the species can be delisted. For the 642 listed endangered and threatened plant species that have recovery plans, we applied recursive partitioning methods to test whether the number of individuals or populations required for delisting can be predicted on the basis of distributional and biological traits, previous abundance at multiple time steps, or a combination of traits and previous abundances. We also tested listing status (threatened or endangered) and the year the recovery plan was written as predictors of recovery criteria. We analyzed separately recovery criteria that were stated as number of populations and as number of individuals (population‐based and individual‐based criteria, respectively). Previous abundances alone were relatively good predictors of population‐based recovery criteria. Fewer populations, but a greater proportion of historically known populations, were required to delist species that had few populations at listing compared with species that had more populations at listing. Previous abundances were also good predictors of individual‐based delisting criteria when models included both abundances and traits. The physiographic division in which the species occur was also a good predictor of individual‐based criteria. Our results suggest managers are relying on previous abundances and patterns of decline as guidelines for setting recovery criteria. This may be justifiable in that previous abundances inform managers of the effects of both intrinsic traits and extrinsic threats that interact and determine extinction risk. Predicción de Criterios de Recuperación para Especies de Plantas en Peligro y Amenazadas con Base en Abundancias Pasadas y Atributos Biológicos     

Use of Community‐Composition Data to Predict the Fecundity and Abundance of Species

Abstract: Species distribution models are critical tools for the prediction of invasive species spread and conservation of biodiversity. The majority of species distribution models have been built with environmental data. Community ecology theory suggests that species co‐occurrence data could also be used to predict current and potential distributions of species. Species assemblages are the products of biotic and environmental constraints on the distribution of individual species and as a result may contain valuable information for niche modeling. We compared the predictive ability of distribution models of annual grassland plants derived from either environmental or community‐composition data. Composition‐based models were built with the presence or absence of species at a site as predictors of site quality, whereas environment‐based models were built with soil chemistry, moisture content, above‐ground biomass, and solar radiation as predictors. The reproductive output of experimentally seeded individuals of 4 species and the abundance of 100 species were used to evaluate the resulting models. Community‐composition data were the best predictors of both the site‐specific reproductive output of sown individuals and the site‐specific abundance of existing populations. Successful community‐based models were robust to omission of data on the occurrence of rare species, which suggests that even very basic survey data on the occurrence of common species may be adequate for generating such models. Our results highlight the need for increased public availability of ecological survey data to facilitate community‐based modeling at scales relevant to conservation.     

Use of Multispecies Occupancy Models to Evaluate the Response of Bird Communities to Forest Degradation Associated with Logging

Forest degradation is arguably the greatest threat to biodiversity, ecosystem services, and rural livelihoods. Therefore, increasing understanding of how organisms respond to degradation is essential for management and conservation planning. We were motivated by the need for rapid and practical analytical tools to assess the influence of management and degradation on biodiversity and system state in areas subject to rapid environmental change. We compared bird community composition and size in managed (ejido, i.e., communally owned lands) and unmanaged (national park) forests in the Sierra Tarahumara region, Mexico, using multispecies occupancy models and data from a 2‐year breeding bird survey. Unmanaged sites had on average higher species occupancy and richness than managed sites. Most species were present in low numbers as indicated by lower values of detection and occupancy associated with logging‐induced degradation. Less than 10% of species had occupancy probabilities >0.5, and degradation had no positive effects on occupancy. The estimated metacommunity size of 125 exceeded previous estimates for the region, and sites with mature trees and uneven‐aged forest stand characteristics contained the highest species richness. Higher estimation uncertainty and decreases in richness and occupancy for all species, including habitat generalists, were associated with degraded young, even‐aged stands. Our findings show that multispecies occupancy methods provide tractable measures of biodiversity and system state and valuable decision support for landholders and managers. These techniques can be used to rapidly address gaps in biodiversity information, threats to biodiversity, and vulnerabilities of species of interest on a landscape level, even in degraded or fast‐changing environments. Moreover, such tools may be particularly relevant in the assessment of species richness and distribution in a wide array of habitats. Uso de Modelos de Ocupación para Múltiples Especies para Evaluar la Respuesta de las Comunidades de Aves a la Degradación de Bosques Asociada con la Tala