Partitioning vegetation response to anthropogenic stress to develop multi-taxa wetland indicators.

Emergent plants can be suitable indicators of anthropogenic stress in coastal wetlands if their responses to natural environmental variation can be parsed from their responses to human activities in and around wetlands. We used hierarchical partitioning to evaluate the independent influence of geomorphology, geography, and anthropogenic stress on common wetland plants of the U.S. Great Lakes coast and developed multi-taxa models indicating wetland condition. A seven-taxon model predicted condition relative to watershed-derived anthropogenic stress, and a four-taxon model predicted condition relative to within-wetland anthropogenic stressors that modified hydrology. The Great Lake on which the wetlands occurred explained an average of about half the variation in species cover, and subdividing the data by lake allowed us to remove that source of variation. We developed lake-specific multi-taxa models for all of the Great Lakes except Lake Ontario, which had no plant species with significant independent effects of anthropogenic stress. Plant responses were both positive (increasing cover with stress) and negative (decreasing cover with stress), and plant taxa incorporated into the lake-specific models differed by Great Lake. The resulting models require information on only a few taxa, rather than all plant species within a wetland, making them easier to implement than existing indicators.     

Estimation of population parameters from aerial counts of North American mallards: a cautionary tale.

The accuracy of population estimates strongly interferes with our ability to obtain unbiased estimates of population parameters based on analyses of time series of population fluctuations. Here we use long-term data on fluctuations in the size of Mallard populations collected as part of the May Breeding Waterfowl Survey covering a large section of North America. We assume a log-linear model of density dependence and use a hierarchical Bayesian state-space approach in which all parameters are assumed to be realizations from a common underlying distribution. Thus, parameters for different populations are not allowed to vary independently of each other. We then simulated independent time series of aerial counts, using the estimated parameters and adding various levels of observation error. These simulations showed that the estimates of stochastic population growth rate and strength of density dependence were biased even when moderate sampling errors were present. In contrast, the estimates of the environmental stochasticity and the carrying capacity were unbiased even for short time series and large observation error. Our results underline the importance of reducing the magnitude of sampling error in the design of large-scale monitoring programs of population fluctuations.     

Hierarchical linear mixed models in multi-stage sampling soil studies

The issue of variances of different soil variables prevailing at different sampling scales is addressed. This topic is relevant for soil science, agronomy and landscape ecology. In multi-stage sampling there are randomness components in each stage of sampling which can be taken into account by introducing random effects in analysis through the use of hierarchical linear mixed models (HLMM). Due to the nested sampling scheme, there are several hierarchical sub-models. The selection of the best model can be carried out through likelihood ratio tests (LRTs) or Wald tests, which are asymptotically equivalent under standard conditions. However, when the comparison leads to a restricted hypothesis of variance components, standard conditions are not maintained, which leads to more elaborated versions of LRTs. These versions are not disseminated among environmental scientists. The present study shows the modeling of soil data from a sampling where sites, fields within sites, transects within fields, and sampling points within transects were selected in order to take samples from different vegetation types (open and shade). For soil data, several sub-models were compared using Wald tests, classic LRTs and adjusted LRTs where the distribution of the test statistic under the null hypothesis is the Chi-square mixture of Chi-square distributions. The inclusion of random effects via HLMM and suggested by the latest version of LRT allowed us to detect effects of vegetation type on soil properties that were not detected under a classical ANOVA.     

Using the negative binomial distribution to model overdispersion in ecological count data

A Poisson process is a commonly used starting point for modeling stochastic variation of ecological count data around a theoretical expectation. However, data typically show more variation than implied by the Poisson distribution. Such overdispersion is often accounted for by using models with different assumptions about how the variance changes with the expectation. The choice of these assumptions can naturally have apparent consequences for statistical inference. We propose a parameterization of the negative binomial distribution, where two overdispersion parameters are introduced to allow for various quadratic mean-variance relationships, including the ones assumed in the most commonly used approaches. Using bird migration as an example, we present hypothetical scenarios on how overdispersion can arise due to sampling, flocking behavior or aggregation, environmental variability, or combinations of these factors. For all considered scenarios, mean-variance relationships can be appropriately described by the negative binomial distribution with two overdispersion parameters. To illustrate, we apply the model to empirical migration data with a high level of overdispersion, gaining clearly different model fits with different assumptions about mean-variance relationships. The proposed framework can be a useful approximation for modeling marginal distributions of independent count data in likelihood-based analyses.     

Seasonal components of avian population change: joint analysis of two large-scale monitoring programs

We present a combined analysis of data from two large-scale surveys of bird populations. The North American Breeding Bird Survey is conducted each summer; the Christmas Bird Count is conducted in early winter. The temporal staggering of these surveys allows investigation of seasonal components of population change, which we illustrate with an examination of the effects of severe winters on the Carolina Wren (Thryothorus ludovicianus). Our analysis uses a hierarchical log-linear model with controls for survey-specific sampling covariates. Temporal change in population size is modeled seasonally, with covariates for winter severity. Overall, the winter-spring seasons are associated with 82% of the total population variation for Carolina Wrens, and an additional day of snow cover during winter-spring is associated with an incremental decline of 1.1% of the population.     

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.     

Species distribution modelling—Effect of design and sample size of pseudo-absence observations

We explored the effect of varying pseudo-absence data in species distribution modelling using empirical data for four real species and simulated data for two imaginary species. In all analyses we used a fixed study area, a fixed set of environmental predictors and a fixed set of presence observations. Next, we added pseudo-absence data generated by different sampling designs and in different numbers to assess their relative importance for the output from the species distribution model. The sampling design strongly influenced the predictive performance of the models while the number of pseudo-absences had minimal effect on the predictive performance. We attribute much of these results to the relationship between the environmental range of the pseudo-absences (i.e. the extent of the environmental space being considered) and the environmental range of the presence observations (i.e. under which environmental conditions the species occurs). The number of generated pseudo-absences had a direct effect on the predicted probability, which translated to different distribution areas. Pseudo-absence observations that fell within grid cells with presence observations were purposely included in our analyses. We discourage the practice of excluding certain pseudo-absence data because it involves arbitrary assumptions about what are (un)suitable environments for the species being modelled.     

Efficiency of composite sampling for estimating a lognormal distribution

In many environmental studies measuring the amount of a contaminant in a sampling unit is expensive. In such cases, composite sampling is often used to reduce data collection cost. However, composite sampling is known to be beneficial for estimating the mean of a population, but not necessarily for estimating the variance or other parameters. As some applications, for example, Monte Carlo risk assessment, require an estimate of the entire distribution, and as the lognormal model is commonly used in environmental risk assessment, in this paper we investigate efficiency of composite sampling for estimating a lognormal distribution. In particular, we examine the magnitude of savings in the number of measurements over simple random sampling, and the nature of its dependence on composite size and the parameters of the distribution utilizing simulation and asymptotic calculations.     

Spatial heterogeneity influences native and nonnative plant species richness

Spatial heterogeneity may have differential effects on the distribution of native and nonnative plant species richness. We examined the effects of spatial heterogeneity on native and nonnative plant species richness distributions in the central part of Rocky Mountain National Park, Colorado, USA. Spatial heterogeneity around vegetation plots was characterized using landscape metrics, environmental/topographic variables (slope, aspect, elevation, and distance from stream or river), and soil variables (nitrogen, clay, and sand). The landscape metrics represented five components of landscape heterogeneity and were measured at four spatial extents (within varying radii of 120, 240, 480, and 960 m) using the FRAGSTATS landscape pattern analysis program. Akaike's Information Criterion adjusted for small sample size (AICc) was used to select the best models from a set of multiple linear regression models developed for native and nonnative plant species richness at four spatial extents and three levels of ecological hierarchy (i.e., landscape, land cover, and community). Both native and nonnative plant species richness were positively correlated with edge density, Simpson's diversity index and interspersion/juxtaposition index, and were negatively correlated with mean patch size. The amount of variation explained at four spatial extents and three hierarchical levels ranged from 30% to 70%. At the landscape level, the best models explained 43% of the variation in native plant species richness and 70% of the variation in nonnative plant species richness (240-m extent). In general, the amount of variation explained was always higher for nonnative plant species richness, and the inclusion of landscape metrics always significantly improved the models. The best models explained 66% of the variation in nonnative plant species richness for both the conifer land cover type and lodgepole pine community. The relative influence of the components of spatial heterogeneity differed for native and nonnative plant species richness and varied with the spatial extent of analysis and levels of ecological hierarchy. The study offers an approach to quantify spatial heterogeneity to improve models of plant biodiversity. The results demonstrate that ecologists must recognize the importance of spatial heterogeneity in managing native and nonnative plant species.     

Rapid assessment of postfire plant invasions in coniferous forests of the western United States.

Fire is a natural part of most forest ecosystems in the western United States, but its effects on nonnative plant invasion have only recently been studied. Also, forest managers are engaging in fuel reduction projects to lessen fire severity, often without considering potential negative ecological consequences such as nonnative plant species introductions. Increased availability of light, nutrients, and bare ground have all been associated with high-severity fires and fuel treatments and are known to aid in the establishment of nonnative plant species. We use vegetation and environmental data collected after wildfires at seven sites in coniferous forests in the western United States to study responses of nonnative plants to wildfire. We compared burned vs. unburned plots and plots treated with mechanical thinning and/or prescribed burning vs. untreated plots for nonnative plant species richness and cover and used correlation analyses to infer the effect of abiotic site conditions on invasibility. Wildfire was responsible for significant increases in nonnative species richness and cover, and a significant decrease in native cover. Mechanical thinning and prescribed fire fuel treatments were associated with significant changes in plant species composition at some sites. Treatment effects across sites were minimal and inconclusive due to significant site and site x treatment interaction effects caused by variation between sites including differences in treatment and fire severities and initial conditions (e.g., nonnative species sources). We used canonical correspondence analysis (CCA) to determine what combinations of environmental variables best explained patterns of nonnative plant species richness and cover. Variables related to fire severity, soil nutrients, and elevation explained most of the variation in species composition. Nonnative species were generally associated with sites with higher fire severity, elevation, percentage of bare ground, and lower soil nutrient levels and lower canopy cover. Early assessments of postfire stand conditions can guide rapid responses to nonnative plant invasions.     

Building statistical models to analyze species distributions.

Models of the geographic distributions of species have wide application in ecology. But the nonspatial, single-level, regression models that ecologists have often employed do not deal with problems of irregular sampling intensity or spatial dependence, and do not adequately quantify uncertainty. We show here how to build statistical models that can handle these features of spatial prediction and provide richer, more powerful inference about species niche relations, distributions, and the effects of human disturbance. We begin with a familiar generalized linear model and build in additional features, including spatial random effects and hierarchical levels. Since these models are fully specified statistical models, we show that it is possible to add complexity without sacrificing interpretability. This step-by-step approach, together with attached code that implements a simple, spatially explicit, regression model, is structured to facilitate self-teaching. All models are developed in a Bayesian framework. We assess the performance of the models by using them to predict the distributions of two plant species (Proteaceae) from South Africa's Cape Floristic Region. We demonstrate that making distribution models spatially explicit can be essential for accurately characterizing the environmental response of species, predicting their probability of occurrence, and assessing uncertainty in the model results. Adding hierarchical levels to the models has further advantages in allowing human transformation of the landscape to be taken into account, as well as additional features of the sampling process.     

A hierarchical model for spatial capture-recapture data

Estimating density is a fundamental objective of many animal population studies. Application of methods for estimating population size from ostensibly closed populations is widespread, but ineffective for estimating absolute density because most populations are subject to short-term movements or so-called temporary emigration. This phenomenon invalidates the resulting estimates because the effective sample area is unknown. A number of methods involving the adjustment of estimates based on heuristic considerations are in widespread use. In this paper, a hierarchical model of spatially indexed capture-recapture data is proposed for sampling based on area searches of spatial sample units subject to uniform sampling intensity. The hierarchical model contains explicit models for the distribution of individuals and their movements, in addition to an observation model that is conditional on the location of individuals during sampling. Bayesian analysis of the hierarchical model is achieved by the use of data augmentation, which allows for a straightforward implementation in the freely available software WinBUGS. We present results of a simulation study that was carried out to evaluate the operating characteristics of the Bayesian estimator under variable densities and movement patterns of individuals. An application of the model is presented for survey data on the flat-tailed horned lizard (Phrynosoma mcallii) in Arizona, USA.     

Predicting species distributions from samples collected along roadsides

Predictive models of species distributions are typically developed with data collected along roads. Roadside sampling may provide a biased (nonrandom) sample; however, it is currently unknown whether roadside sampling limits the accuracy of predictions generated by species distribution models. We tested whether roadside sampling affects the accuracy of predictions generated by species distribution models by using a prospective sampling strategy designed specifically to address this issue. We built models from roadside data and validated model predictions at paired locations on unpaved roads and 200 m away from roads (off road), spatially and temporally independent from the data used for model building. We predicted species distributions of 15 bird species on the basis of point-count data from a landbird monitoring program in Montana and Idaho (U.S.A.). We used hierarchical occupancy models to account for imperfect detection. We expected predictions of species distributions derived from roadside-sampling data would be less accurate when validated with data from off-road sampling than when it was validated with data from roadside sampling and that model accuracy would be differentially affected by whether species were generalists, associated with edges, or associated with interior forest. Model performance measures (kappa, area under the curve of a receiver operating characteristic plot, and true skill statistic) did not differ between model predictions of roadside and off-road distributions of species. Furthermore, performance measures did not differ among edge, generalist, and interior species, despite a difference in vegetation structure along roadsides and off road and that 2 of the 15 species were more likely to occur along roadsides. If the range of environmental gradients is surveyed in roadside-sampling efforts, our results suggest that surveys along unpaved roads can be a valuable, unbiased source of information for species distribution models.     

Poor-data and data-poor species stock assessment using a Bayesian hierarchical approach

Appropriate inference for stocks or species with low-quality data (poor data) or limited data (data poor) is extremely important. Hierarchical Bayesian methods are especially applicable to small-area, small-sample-size estimation problems because they allow poor-data species to borrow strength from species with good-quality data. We used a hammerhead shark complex as an example to investigate the advantages of using hierarchical Bayesian models in assessing the status of poor-data and data-poor exploited species. The hammerhead shark complex (Sphyrna spp.) along the Atlantic and Gulf of Mexico coasts of the United States is composed of three species: the scalloped hammerhead (S. lewini), the great hammerhead (S. mokarran), and the smooth hammerhead (S. zygaena) sharks. The scalloped hammerhead comprises 70-80% of the catch and has catch and relative abundance data of good quality, whereas great and smooth hammerheads have relative abundance indices that are both limited and of low quality presumably because of low stock density and limited sampling. Four hierarchical Bayesian state-space surplus production models were developed to simulate variability in population growth rates, carrying capacity, and catchability of the species. The results from the hierarchical Bayesian models were considerably more robust than those of the nonhierarchical models. The hierarchical Bayesian approach represents an intermediate strategy between traditional models that assume different population parameters for each species and those that assume all species share identical parameters. Use of the hierarchical Bayesian approach is suggested for future hammerhead shark stock assessments and for modeling fish complexes with species-specific data, because the poor-data species can borrow strength from the species with good data, making the estimation more stable and robust.     

Climate, habitat, and species interactions at different scales determine the structure of a Neotropical bat community

Climate, habitat, and species interactions are factors that control community properties (e.g., species richness, abundance) across various spatial scales. Usually, researchers study how a few properties are affected by one factor in isolation and at one scale. Hence, there are few multi-scale studies testing how multiple controlling factors simultaneously affect community properties at different scales. We ask whether climate, habitat structure, or insect resources at each of three spatial scales explains most of the variation in six community properties and which theory best explains the distribution of selected community properties across a rainfall gradient. We studied a Neotropical insectivorous bat ensemble in the Isthmus of Panama with acoustic monitoring techniques. Using climatological data, habitat surveys, and insect captures in a hierarchical sampling design we determined how much variation of the community properties was explained by the three factors employing two approaches for variance partitioning. Our results revealed that most of the variation in species richness, total abundance, and feeding activity occurred at the smallest spatial scale and was explained by habitat structure. In contrast, climate at large scales explained most of the variation in individual species' abundances. Although each species had an idiosyncratic response to the gradient, species richness peaked at intermediate levels of precipitation, whereas total abundance was very similar across sites, suggesting density compensation. All community properties responded in a different manner to the factor and scale under consideration.     

若尔盖高寒沼泽植物群落与环境因子的关系

了解高寒沼泽植物群落与环境因子之间的关系,对揭示物种对生态环境的需求以及高寒沼泽的保护管理有重要意义,然而目前尚不清楚影响若尔盖高寒沼泽植物分布的主要环境因素。作者于2011—2012年间在若尔盖高寒沼泽随机调查32个样地,收集了植被样方数据、土壤理化性质数据、地表水位等实测数据,借助方差分析、CCA分析和相关分析等经典的统计方法,研究高寒沼泽植物群落物种组成和不同地表积水状况下多样性的差异,以及高寒沼泽植物群落与环境因子的关系。结果显示：若尔盖高寒沼泽植物共有151个种,隶属于39科98属;多样性指数随地表积水的减少呈现出增加的趋势。典范对应分析（CCA）结果表明,第一轴与水位、土壤水分含量、土壤有机碳、土壤总氮显著负相关,与土壤容重显著正相关,第二轴与裸斑率和鼠洞数显著相关。前2个轴一共解释了79．8％的物种与环境因子的关系,与环境因子的相关系数均在0．87以上;其中,水分和土壤养分条件是影响物种分布的主导因子,其次是裸斑面及啮齿动物活动。植物多样性指数与土壤水分负相关,与土壤有机碳和土壤全氮显著负相关,与土壤容重显著正相关。若尔盖沼泽植物群落物种多样性除了受土壤水分和养分条件的影响以外,啮齿动物的活动可能是影响若尔盖沼泽物种分布以及促进群落进一步演替的主要因素之一。     

Surrogates for Macrofungi and Mosses in Reservation Planning

Abstract: Our knowledge of cryptogam taxonomy and species distributions is currently too poor to directly plan for their conservation. We used inventory data from four distinct vegetation types, near Hobart Tasmania, to address the proposition that vegetation type, vascular plant taxon composition, and environmental variables can act as surrogates for mosses and macrofungi in reservation planning. The four vegetation types proved distinct in their taxon composition for all macrofungi, mosses, and vascular plants. We tested the strength of the relationships between the composition of cryptogam taxonomic groups and vascular plant composition and between the environmental variables and canopy cover. Taxon composition of woody vascular plants and vascular plants was the best predictor of the taxon composition of mosses and macrofungi. Combinations of environmental variables and canopy cover were also strong predictors of the taxon composition of mosses and macrofungi. We used an optimization routine for vascular plant taxa and woody plant species and determined the representation of cryptogam taxa in these selections. We identified sites with approximately 10% and 30% of the greatest proportions of vascular plants and woody vascular plants and calculated representation of mosses and macrofungi at these sites. We compared the results of these site selections with random site selections and random selections stratified by vegetation type. Random selection of sites by vegetation type generally captured more cryptogams than site selection by vascular plants at the 10% level. Vascular plant and woody plant taxon composition, vegetation type, and environmental and structural characteristics, all showed promise as surrogates for capturing common cryptogams in reserve systems.     

Fitting complex population models by combining particle filters with Markov chain Monte Carlo

We show how a recent framework combining Markov chain Monte Carlo (MCMC) with particle filters (PFMCMC) may be used to estimate population state-space models. With the purpose of utilizing the strengths of each method, PFMCMC explores hidden states by particle filters, while process and observation parameters are estimated using an MCMC algorithm. PFMCMC is exemplified by analyzing time series data on a red kangaroo (Macropus rufus) population in New South Wales, Australia, using MCMC over model parameters based on an adaptive Metropolis-Hastings algorithm. We fit three population models to these data; a density-dependent logistic diffusion model with environmental variance, an unregulated stochastic exponential growth model, and a random-walk model. Bayes factors and posterior model probabilities show that there is little support for density dependence and that the random-walk model is the most parsimonious model. The particle filter Metropolis-Hastings algorithm is a brute-force method that may be used to fit a range of complex population models. Implementation is straightforward and less involved than standard MCMC for many models, and marginal densities for model selection can be obtained with little additional effort. The cost is mainly computational, resulting in long running times that may be improved by parallelizing the algorithm.     

Using multilevel spatial models to understand salamander site occupancy patterns after wildfire

Studies of the distribution of elusive forest wildlife have suffered from the confounding of true presence with the uncertainty of detection. Occupancy modeling, which incorporates probabilities of species detection conditional on presence, is an emerging approach for reducing observation bias. However, the current likelihood modeling framework is restrictive for handling unexplained sources of variation in the response that may occur when there are dependence structures such as smaller sampling units that are nested within larger sampling units. We used multilevel Bayesian occupancy modeling to handle dependence structures and to partition sources of variation in occupancy of sites by terrestrial salamanders (family Plethodontidae) within and surrounding an earlier wildfire in western Oregon, USA. Comparison of model fit favored a spatial N-mixture model that accounted for variation in salamander abundance over models that were based on binary detection/non-detection data. Though catch per unit effort was higher in burned areas than unburned, there was strong support that this pattern was due to a higher probability of capture for individuals in burned plots. Within the burn, the odds of capturing an individual given it was present were 2.06 times the odds outside the burn, reflecting reduced complexity of ground cover in the burn. Ther was weak support that true occupancy was lower within the burned area. While the odds of occupancy in the burn were 0.49 times the odds outside the burn among the five species, the magnitude of variation attributed to the burn was small in comparison to variation attributed to other landscape variables and to unexplained, spatially autocorrelated random variation. While ordinary occupancy models may separate the biological pattern of interest from variation in detection probability when all sources of variation are known, the addition of random effects structures for unexplained sources of variation in occupancy and detection probability may often more appropriately represent levels of uncertainty.