A practical sampling design for acoustic surveys of bats

Acoustic surveys are widely used for describing bat occurrence and activity patterns and are increasingly important for addressing concerns for habitat management, wind energy, and disease on bat populations. Designing these surveys presents unique challenges, particularly when a probabilistic sample is required for drawing inference to unsampled areas. Sampling frame errors and other logistical constraints often require survey sites to be dropped from the sample and new sites added. Maintaining spatial balance and representativeness of the sample when these changes are made can be problematic. Spatially balanced sampling designs recently developed to support aquatic surveys along rivers provide solutions to a number of practical challenges faced by bat researchers and allow for sample site additions and deletions, support unequal-probability selection of sites, and provide an approximately unbiased local neighborhood-weighted variance estimator that is efficient for spatially structured populations such as is typical for bats. We implemented a spatially balanced design to survey canyon bat (Parastrellus hesperus) activity along a stream network. The spatially balanced design accommodated typical logistical challenges and yielded a 25% smaller estimated standard error for the mean activity level than the usual simple random sampling estimator. Spatially balanced designs have broad application to bat research and monitoring programs and will improve studies relying on model-based inference (e.g., occupancy models) by providing flexibility and protection against violations of the independence assumption, even if design-based estimators are not used. Our approach is scalable and can be used for pre- and post-construction surveys along wind turbine arrays and for regional monitoring programs. © 2011 The Wildlife Society.     

Monitoring temporal trends in spatially structured populations: how should sampling effort be allocated between space and time?

Estimating temporal trends in spatially structured populations has a critical role to play in understanding regional changes in biological populations and developing management strategies. Designing effective monitoring programmes to estimate these trends requires important decisions to be made about how to allocate sampling effort among spatial replicates (i.e. number of sites) and temporal replicates (i.e. how often to survey) to minimise uncertainty in trend estimates. In particular, the optimal mix of spatial and temporal replicates is likely to depend upon the spatial and temporal correlations in population dynamics. Although there has been considerable interest in the ecological literature on understanding spatial and temporal correlations in species’ population dynamics, little attention has been paid to its consequences for monitoring design. We address this issue using model‐based survey design to identify the optimal allocation of sampling effort among spatial and temporal replicates for estimating population trends under different levels of spatial and temporal correlation. Based on linear trends, we show that how we should allocate sampling effort among spatial and temporal replicates depends crucially on the spatial and temporal correlations in population dynamics, environmental variation, observation error and the spatial variation in temporal trends. When spatial correlation is low and temporal correlation is high, the best option is likely to be to sample many sites infrequently, particularly when observation error and/or spatial variation in temporal trends are high. When spatial correlation is high and temporal correlation is low, the best option is likely to be to sample few sites frequently, particularly when observation error and/or spatial variation in temporal trends are low. When abundances are spatially independent, it is always preferable to maximise spatial replication. This provides important insights into how spatio‐temporal monitoring programmes should be designed to estimate temporal trends in spatially structured populations.     

Additional Circular Systematic Sampling Methods

Several systematic sampling methods have been used to estimate the population mean when size of the population is a multiple of sample size. Among these, only few methods have been extended and used to estimate mean of the population when its size is not a multiple of sample size. In this paper, new methods called balanced circular systematic sampling and centered circular systematic sampling are introduced by extending balanced systematic sampling method and centered systematic sampling method respectively. These methods are compared with circular systematic sampling using average variance of corrected sample means for populations exhibiting approximate linear and parabolic trends. The suggested methods are found suitable to estimate the population mean.     

Optimizing methods to estimate zooplankton concentration based on generalized patterns of patchiness inside ballast tanks and ballast water discharges

Zooplankton populations are spatially heterogeneous in nature and inside ship ballast tanks. Sampling methods should take heterogeneity into account, particularly when estimating quantitative variables such as abundance or concentration. It is particularly important to generate unbiased estimates of zooplankton concentration in ballast water when assessing compliance with new international ballast water discharge standards. We measured spatial heterogeneity of zooplankton within ballast water using three sampling methodologies. In‐tank pump samples were collected at fixed depths within the vertical part of the ballast tank (side tank). Vertical net‐haul samples were collected from the upper portion of the tank as a depth‐integrated and historically relevant method. In‐line, time‐integrated samples were collected during ballast discharge by an isokinetic sample probe, likely representing the double bottom part of the ballast tank. The bias and precision associated with each sampling method were evaluated in reference to the estimated average abundance of the entire ballast tank, which was modeled from the data collected by all methods. In‐tank pump samples provided robust evidence for vertical stratification of zooplankton concentration in the side tank. A consistent trend was also observed for in‐line discharge samples, with zooplankton concentration decreasing through time as the ballast tank is being discharged. Sample representativeness, as compared to the tank average, varied depending on the depth or tank volume discharged. In‐line discharge samples provided the least biased and most precise estimate of average tank abundance (having lowest mean squared error) when collected during the time frame of 20%–60% of the tank volume being discharged. Results were consistent across five trips despite differences in ballast water source, season, and age.     

Viability analyses with habitat-based metapopulation models

Population viability analysis (PVA) models incorporate spatial dynamics in different ways. At one extreme are the occupancy models that are based on the number of occupied populations. The simplest occupancy models ignore the location of populations. At the other extreme are individual-based models, which describe the spatial structure with the location of each individual in the population, or the location of territories or home ranges. In between these are spatially structured metapopulation models that describe the dynamics of each population with structured demographic models and incorporate spatial dynamics by modeling dispersal and temporal correlation among populations. Both dispersal and correlation between each pair of populations depend on the location of the populations, making these models spatially structured. In this article, I describe a method that expands spatially structured metapopulation models by incorporating information about habitat relationships of the species and the characteristics of the landscape in which the metapopulation exists. This method uses a habitat suitability map to determine the spatial structure of the metapopulation, including the number, size, and location of habitat patches in which subpopulations of the metapopulation live. The habitat suitability map can be calculated in a number of different ways, including statistical analyses (such as logistic regression) that find the relationship between the occurrence (or, density) of the species and independent variables which describe its habitat requirements. The habitat suitability map is then used to calculate the spatial structure of the metapopulation, based on species-specific characteristics such as the home range size, dispersal distance, and minimum habitat suitability for reproduction. Received: April 1, 1999 / Accepted: October 29, 1999     

Environmental colour intensifies the Moran effect when population dynamics are spatially heterogeneous

Evidence for synchronous fluctuations of spatially separated populations is ubiquitous in the literature, including accounts within and across taxa. Among the few mechanisms explaining this phenomenon is the Moran effect, whereby independent populations are synchronized by spatially correlated environmental disturbances. The body of research on the Moran effect predominantly assumes that environmental disturbances within a local site are serially uncorrelated; that is, successive observations in time at a particular local site are independent. Yet, many environmental variables are known to possess strong temporal autocorrelation – a character which has often been described as 'colour'. The omission of environmental colour from research on the Moran effect may be due in part to the lack of methods capable of generating sets of time series with a desired colour and spatial correlation. Here I present a novel and simple method designated as 'phase partnering' to generate such sets of time series and I investigate the combined impact of spatial correlation and environmental colour on population synchrony in two common models of population dynamics. For linear population dynamics, and for a subset of nonlinear population dynamics, coloured environments intensify the Moran effect when population dynamics are spatially heterogeneous; in coloured environments the spatial correlation between populations more closely mimics the spatial correlation between their respective environments. Given that most environmental variables are coloured, these results imply that the Moran effect may be a far more significant driver of regional-scale population and interspecific synchrony than is currently believed.     

Emergence of long‐term balanced polymorphism under cyclic selection of spatially variable magnitude

A fundamental question in evolutionary biology is what promotes genetic variation at nonneutral loci, a major precursor to adaptation in changing environments. In particular, balanced polymorphism under realistic evolutionary models of temporally varying environments in finite natural populations remains to be demonstrated. Here, we propose a novel mechanism of balancing selection under temporally varying fitnesses. Using forward‐in‐time computer simulations and mathematical analysis, we show that cyclic selection that spatially varies in magnitude, such as along an environmental gradient, can lead to elevated levels of nonneutral genetic polymorphism in finite populations. Balanced polymorphism is more likely with an increase in gene flow, magnitude and period of fitness oscillations, and spatial heterogeneity. This polymorphism‐promoting effect is robust to small systematic fitness differences between competing alleles or to random environmental perturbation. Furthermore, we demonstrate analytically that protected polymorphism arises as spatially heterogeneous cyclic fitness oscillations generate a type of storage effect that leads to negative frequency dependent selection. Our findings imply that spatially variable cyclic environments can promote elevated levels of nonneutral genetic variation in natural populations.     

Geostatistical approaches and optimal additional sampling schemes for spatial patterns and future sampling of bird diversity

Aim To evaluate geostatistical approaches, namely kriging, co‐kriging and geostatistical simulation, and to develop an optimal sampling design for mapping the spatial patterns of bird diversity, estimating their spatial autocorrelations and selecting additional samples of bird diversity in a 2450 km² basin. Location Taiwan. Methods Kriging, co‐kriging and simulated annealing are applied to estimate and simulate the spatial patterns of bird diversity. In addition, kriging and co‐kriging with a genetic algorithm are used to optimally select further samples to improve the kriging and co‐kriging estimations. The association between bird diversity and elevation, and bird diversity and land cover, is analysed with estimated and simulated maps. Results The Simpson index correlates spatially with the normalized difference vegetation index (NDVI) within the micro‐scale and the macro‐scale in the study basin, but the Shannon diversity index only correlates spatially with NDVI within the micro‐scale. Co‐kriging and simulated annealing simulation accurately simulate the statistical and spatial patterns of bird diversity. The mean estimated diversity and the simulated diversity increase with elevation and decrease with increasing urbanization. The proposed optimal sampling approach selects 43 additional sampling sites with a high spatial estimation variance in bird diversity. Main conclusions Small‐scale variations dominate the total spatial variation of the observed diversity due to a lack of spatial information and insufficient sampling. However, simulations of bird diversity consistently capture the sampling statistics and spatial patterns of the observed bird diversity. The data thus accumulated can be used to understand the spatial patterns of bird diversity associated with different types of land cover and elevation, and to optimize sample selection. Co‐kriging combined with a genetic algorithm yields additional optimal sampling sites, which can be used to augment existing sampling points in future studies of bird diversity.     

How noise and coupling influence leading indicators of population extinction in a spatially extended ecological system

Anticipating critical transitions in spatially extended systems is a key topic of interest to ecologists. Gradually declining metapopulations are an important example of a spatially extended biological system that may exhibit a critical transition. Theory for spatially extended systems approaching extinction that accounts for environmental stochasticity and coupling is currently lacking. Here, we develop spatially implicit two-patch models with additive and multiplicative forms of environmental stochasticity that are slowly forced through population collapse, through changing environmental conditions. We derive patch-specific expressions for candidate indicators of extinction and test their performance via a simulation study. Coupling and spatial heterogeneities decrease the magnitude of the proposed indicators in coupled populations relative to isolated populations, and the noise regime and the degree of coupling together determine trends in summary statistics. This theory may be readily applied to other spatially extended ecological systems, such as coupled infectious disease systems on the verge of elimination.     

An improved survey method for monitoring population trends of Golden‐winged Warblers and other patchily distributed birds

Conventional surveys designed to monitor common and widespread species may fail to adequately track population changes of rare or patchily distributed species that are often of high conservation concern. We evaluated the performance of a new monitoring approach that employs both a spatially balanced sampling design and a targeted survey protocol designed to estimate population trends of one such patchily distributed species, the Golden‐winged Warbler (Vermivora chrysoptera), in the Appalachian Mountains Bird Conservation Region (BCR 28), USA. Our spatially balanced survey consisted of 105 sample quads (one‐quarter Delorme Atlas pages) across the current range of Golden‐winged Warblers within BCR 28, each with five sample points located in early successional habitat. From 2009 to 2013, collaborators visited each sample point once per year during the peak breeding season and conducted a 17‐min survey consisting of passive observation and playback of conspecific songs and mobbing vocalizations. We used multi‐season, single‐species occupancy models to estimate probability of quad occupancy, detection probability, and occupancy dynamics for Golden‐winged Warblers and closely related Blue‐winged Warblers (Vermivora cyanoptera). Our survey protocol resulted in high estimates of detection probability for Golden‐winged (92%) and Blue‐winged (79%) warblers, with 47% and 56% of quads estimated to be initially occupied, respectively. Derived population trend estimates (λ) indicated an average decline in population of 6% for Golden‐winged Warblers and 7% for Blue‐winged Warblers, resulting in estimated 21% and 22% declines, respectively, in quad occupancy after 5 yr. Our results demonstrate that coupling a spatially balanced survey design in appropriate habitat with a playback protocol to increase detection rates is a viable strategy for tracking populations of Golden‐winged Warblers in the Appalachian Mountains BCR. Similar survey methods should be considered for other rare, declining, or patchily distributed bird species that require targeted monitoring.     

A spatial model of bird abundance as adjusted for detection probability

Modeling the spatial distribution of animals can be complicated by spatial and temporal effects (i.e. spatial autocorrelation and trends in abundance over time) and other factors such as imperfect detection probabilities and observation-related nuisance variables. Recent advances in modeling have demonstrated various approaches that handle most of these factors but which require a degree of sampling effort (e.g. replication) not available to many field studies. We present a two-step approach that addresses these challenges to spatially model species abundance. Habitat, spatial and temporal variables were handled with a Bayesian approach which facilitated modeling hierarchically structured data. Predicted abundance was subsequently adjusted to account for imperfect detection and the area effectively sampled for each species. We provide examples of our modeling approach for two endemic Hawaiian nectarivorous honeycreepers: 'i'iwi Vestiaria coccinea and 'apapane Himatione sanguinea .     

降水量和温度对植被覆盖指数影响的空间非平稳性特征——以新疆伊犁河谷地区为例

多变量空间相关分析多基于时间序列数据,对数据时长与统计要求严格,空间非平稳性特征分析可以利用单期数据分析多变量之间的相关性。通过空间变系数回归模型分析了2006年和2011年的新疆伊犁地区降水量和温度对植被覆盖度指数影响的空间变化特征,利用局部线性地理加权回归(GWR)方法估计得到了回归系数曲面,揭示出变量间相互影响的空间异质性,同时利用线性回归最小二乘估计进行了对比。结果表明:(1)空间变系数回归模型可以用于变量间的空间相关分析;(2)局部线性GWR估计方法明显优于线性回归最小二乘估计;(3)拟合结果表明,伊犁地区降水量和温度对植被覆盖指数的影响具有显著的空间非平稳性特征;(4)模型估计误差是降水、气温之外的地形、地貌及人类活动等多种因素造成的,需进一步研究。方法可为具有空间非平稳性特征变量间空间相关性分析以及植被覆盖指数的空间模拟分布提供思路和方法。     

有色环境噪音对空间异质种群动态同步性的影响

随着种群动态和空间结构研究兴趣的增加,激发了大量的有关空间同步性的理论和实验的研究工作。空间种群的同步波动现象在自然界广泛存在,它的影响和原因引起了很多生态学家的兴趣。Moran定理是一个非常重要的解释。但以往的研究大多假设环境变化为空间相关的白噪音。越来越多的研究表明很多环境变化的时间序列具有正的时间自相关性,也就是说用红噪音来描述更加合理。因此,推广经典的Moran效应来处理空间相关红噪音的情形很有必要。利用线性的二阶自回归过程的种群模型,推导了两种群空间同步性与种群动态异质性和环境变化的时间相关性(即环境噪音的颜色)之间的关系。深入分析了种群异质性和噪音颜色对空间同步性的影响。结果表明种群动态异质性不利于空间同步性,但详细的关系比较复杂。而红色噪音的同步能力体现在两方面:一方面,本身的相关性对同步性有贡献;另一方面,环境变化时间相关性可以通过改变种群密度依赖来影响同步性,但对同步性的影响并无一致性的结论,依赖于种群的平均动态等因素。这些结果对理解同步性的机理、利用同步机理来制定物种保护策略和害虫防治都有重要的意义。     

Comparative effects of Aroclor 1254 (polychlorinated biphenyls) and phenanthrene on glucose uptake by freshwater microbial populations.

The effects of polychlorinated biphenyl (PCB) and phenanthrene stress on glucose uptake by natural microbial populations were examined by the heterotrophic potential technique. Temporal and spatial distributions in glucose uptake velocities were examined for natural samples as well as PCB- and phenanthrene-stressed samples. Statistical analysis indicated significant variability among the various samples. It was demonstrated that the environmental variables contributed significantly to the variability in uptake kinetics. Although general trends indicated a PCB-induced stimulation in uptake velocities, these trends were in part masked by sample variability. Data analysis indicated no statistically significant PCB or phenanthrene effect on either total glucose uptake velocities or the proportion of 14CO2 evolved, as compared to natural unstressed samples.     

Comparative effects of Aroclor 1254 (polychlorinated biphenyls) and phenanthrene on glucose uptake by freshwater microbial populations.

The effects of polychlorinated biphenyl (PCB) and phenanthrene stress on glucose uptake by natural microbial populations were examined by the heterotrophic potential technique. Temporal and spatial distributions in glucose uptake velocities were examined for natural samples as well as PCB- and phenanthrene-stressed samples. Statistical analysis indicated significant variability among the various samples. It was demonstrated that the environmental variables contributed significantly to the variability in uptake kinetics. Although general trends indicated a PCB-induced stimulation in uptake velocities, these trends were in part masked by sample variability. Data analysis indicated no statistically significant PCB or phenanthrene effect on either total glucose uptake velocities or the proportion of 14CO2 evolved, as compared to natural unstressed samples.     

Metapopulation dynamics in a changing climate: Increasing spatial synchrony in weather conditions drives metapopulation synchrony of a butterfly inhabiting a fragmented landscape

Habitat fragmentation and climate change are both prominent manifestations of global change, but there is little knowledge on the specific mechanisms of how climate change may modify the effects of habitat fragmentation, for example, by altering dynamics of spatially structured populations. The long‐term viability of metapopulations is dependent on independent dynamics of local populations, because it mitigates fluctuations in the size of the metapopulation as a whole. Metapopulation viability will be compromised if climate change increases spatial synchrony in weather conditions associated with population growth rates. We studied a recently reported increase in metapopulation synchrony of the Glanville fritillary butterfly (Melitaea cinxia) in the Finnish archipelago, to see if it could be explained by an increase in synchrony of weather conditions. For this, we used 23 years of butterfly survey data together with monthly weather records for the same period. We first examined the associations between population growth rates within different regions of the metapopulation and weather conditions during different life‐history stages of the butterfly. We then examined the association between the trends in the synchrony of the weather conditions and the synchrony of the butterfly metapopulation dynamics. We found that precipitation from spring to late summer are associated with the M. cinxia per capita growth rate, with early summer conditions being most important. We further found that the increase in metapopulation synchrony is paralleled by an increase in the synchrony of weather conditions. Alternative explanations for spatial synchrony, such as increased dispersal or trophic interactions with a specialist parasitoid, did not show paralleled trends and are not supported. The climate driven increase in M. cinxia metapopulation synchrony suggests that climate change can increase extinction risk of spatially structured populations living in fragmented landscapes by altering their dynamics.     

Spatial heterogeneity in climate change effects decouples the long‐term dynamics of wild reindeer populations in the high Arctic

The ‘Moran effect’ predicts that dynamics of populations of a species are synchronized over similar distances as their environmental drivers. Strong population synchrony reduces species viability, but spatial heterogeneity in density dependence, the environment, or its ecological responses may decouple dynamics in space, preventing extinctions. How such heterogeneity buffers impacts of global change on large‐scale population dynamics is not well studied. Here, we show that spatially autocorrelated fluctuations in annual winter weather synchronize wild reindeer dynamics across high‐Arctic Svalbard, while, paradoxically, spatial variation in winter climate trends contribute to diverging local population trajectories. Warmer summers have improved the carrying capacity and apparently led to increased total reindeer abundance. However, fluctuations in population size seem mainly driven by negative effects of stochastic winter rain‐on‐snow (ROS) events causing icing, with strongest effects at high densities. Count data for 10 reindeer populations 8–324 km apart suggested that density‐dependent ROS effects contributed to synchrony in population dynamics, mainly through spatially autocorrelated mortality. By comparing one coastal and one ‘continental’ reindeer population over four decades, we show that locally contrasting abundance trends can arise from spatial differences in climate change and responses to weather. The coastal population experienced a larger increase in ROS, and a stronger density‐dependent ROS effect on population growth rates, than the continental population. In contrast, the latter experienced stronger summer warming and showed the strongest positive response to summer temperatures. Accordingly, contrasting net effects of a recent climate regime shift—with increased ROS and harsher winters, yet higher summer temperatures and improved carrying capacity—led to negative and positive abundance trends in the coastal and continental population respectively. Thus, synchronized population fluctuations by climatic drivers can be buffered by spatial heterogeneity in the same drivers, as well as in the ecological responses, averaging out climate change effects at larger spatial scales.     

Spatially balanced sampling and ground‐level imagery for vegetation monitoring on reclaimed well pads

Land reclamation associated with natural gas development has become increasingly important to mitigate land surface disturbance in western North America. Since well pads occur on sites with multiple land use and ownership, the progress and outcomes of these efforts are of interest to multiple stakeholders including industry, practitioners and consultants, regulatory agents, private landowners, and the scientific community. Reclamation success criteria often vary within, and among, government agencies and across land ownership type. Typically, reclamation success of a well pad is judged by comparing vegetation cover from a single transect on the pad to a single transect in an adjacent reference site and data are collected by a large number of technicians with various field monitoring skills. We utilized “SamplePoint” image analysis software and a spatially balanced sampling design, called balanced acceptance sampling, to demonstrate how spatially explicit quantitative data can be used to determine if sites are meeting various reclamation success criteria and used chi‐square tests to show how sites in vegetation percent cover differ from a statistical standpoint. This method collects field data faster than traditional methods. We demonstrate how quantitative and spatially explicit data can be utilized by multiple stakeholders, how it can improve upon current reference site selection, how it can satisfy reclamation monitoring requirements for multiple regulatory agencies, how it may help improve future seed mix selection, and discuss how it may reduce costs for operations responsible for reclamation and how it may reduce observer bias.     

Modeling unobserved sources of heterogeneity in animal abundance using a Dirichlet process prior

Summary .   In surveys of natural populations of animals, a sampling protocol is often spatially replicated to collect a representative sample of the population. In these surveys, differences in abundance of animals among sample locations may induce spatial heterogeneity in the counts associated with a particular sampling protocol. For some species, the sources of heterogeneity in abundance may be unknown or unmeasurable, leading one to specify the variation in abundance among sample locations stochastically. However, choosing a parametric model for the distribution of unmeasured heterogeneity is potentially subject to error and can have profound effects on predictions of abundance at unsampled locations. In this article, we develop an alternative approach wherein a Dirichlet process prior is assumed for the distribution of latent abundances. This approach allows for uncertainty in model specification and for natural clustering in the distribution of abundances in a data-adaptive way. We apply this approach in an analysis of counts based on removal samples of an endangered fish species, the Okaloosa darter. Results of our data analysis and simulation studies suggest that our implementation of the Dirichlet process prior has several attractive features not shared by conventional, fully parametric alternatives.     

Revealing cryptic spatial patterns in genetic variability by a new multivariate method

Increasing attention is being devoted to taking landscape information into account in genetic studies. Among landscape variables, space is often considered as one of the most important. To reveal spatial patterns, a statistical method should be spatially explicit, that is, it should directly take spatial information into account as a component of the adjusted model or of the optimized criterion. In this paper we propose a new spatially explicit multivariate method, spatial principal component analysis (sPCA), to investigate the spatial pattern of genetic variability using allelic frequency data of individuals or populations. This analysis does not require data to meet Hardy-Weinberg expectations or linkage equilibrium to exist between loci. The sPCA yields scores summarizing both the genetic variability and the spatial structure among individuals (or populations). Global structures (patches, clines and intermediates) are disentangled from local ones (strong genetic differences between neighbors) and from random noise. Two statistical tests are proposed to detect the existence of both types of patterns. As an illustration, the results of principal component analysis (PCA) and sPCA are compared using simulated datasets and real georeferenced microsatellite data of Scandinavian brown bear individuals (Ursus arctos). sPCA performed better than PCA to reveal spatial genetic patterns. The proposed methodology is implemented in the adegenet package of the free software R.