Accounting for uncertainty in ecological analysis: the strengths and limitations of hierarchical statistical modeling

Analyses of ecological data should account for the uncertainty in the process(es) that generated the data. However, accounting for these uncertainties is a difficult task, since ecology is known for its complexity. Measurement and/or process errors are often the only sources of uncertainty modeled when addressing complex ecological problems, yet analyses should also account for uncertainty in sampling design, in model specification, in parameters governing the specified model, and in initial and boundary conditions. Only then can we be confident in the scientific inferences and forecasts made from an analysis. Probability and statistics provide a framework that accounts for multiple sources of uncertainty. Given the complexities of ecological studies, the hierarchical statistical model is an invaluable tool. This approach is not new in ecology, and there are many examples (both Bayesian and non-Bayesian) in the literature illustrating the benefits of this approach. In this article, we provide a baseline for concepts, notation, and methods, from which discussion on hierarchical statistical modeling in ecology can proceed. We have also planted some seeds for discussion and tried to show where the practical difficulties lie. Our thesis is that hierarchical statistical modeling is a powerful way of approaching ecological analysis in the presence of inevitable but quantifiable uncertainties, even if practical issues sometimes require pragmatic compromises.

Space–time zero-inflated count models of Harbor seals

Environmental data are spatial, temporal, and often come with many zeros. In this paper, we included space–time random effects in zero-inflated Poisson (ZIP) and ‘hurdle’ models to investigate haulout patterns of harbor seals on glacial ice. The data consisted of counts, for 18 dates on a lattice grid of samples, of harbor seals hauled out on glacial ice in Disenchantment Bay, near Yakutat, Alaska. A hurdle model is similar to a ZIP model except it does not mix zeros from the binary and count processes. Both models can be used for zero-inflated data, and we compared space–time ZIP and hurdle models in a Bayesian hierarchical model. Space–time ZIP and hurdle models were constructed by using spatial conditional autoregressive (CAR) models and temporal first-order autoregressive (AR(1)) models as random effects in ZIP and hurdle regression models. We created maps of smoothed predictions for harbor seal counts based on ice density, other covariates, and spatio-temporal random effects. For both models predictions around the edges appeared to be positively biased. The linex loss function is an asymmetric loss function that penalizes overprediction more than underprediction, and we used it to correct for prediction bias to get the best map for space–time ZIP and hurdle models.

Dynamic Routing in Translucent WDM Optical Networks

Translucent WDM optical networks use sparse placement of regenerators to overcome the impairments and wavelength contention introduced by fully transparent networks, and achieve a performance close to fully opaque networks with much less cost. Our previous study proved the feasibility of translucent networks using sparse regeneration technique. We addressed the placement of regenerators based on static schemes allowing only fixed number of regenerators at fixed locations. This paper furthers the study by proposing a suite of dynamical routing schemes. Dynamic allocation, advertisement and discovery of regeneration resources are proposed to support sharing transmitters and receivers between regeneration and access functions. This study follows the current trend in optical networking industry by utilizing extension of IP control protocols. Dynamic routing algorithms, aware of current regeneration resources and link states, are designed to smartly route the connection requests under quality constraints. A hierarchical network model, supported by the MPLS-based control plane, is also proposed to provide scalability. Experiments show that network performance is improved without placement of extra regenerators.

Cetaceans of Southeast Alaska: Distribution and Seasonal Occurrence

Aim To assess the distribution, group size, seasonal occurrence and annual trends of cetaceans. Location The study area included all major inland waters of Southeast Alaska. Methods Between 1991 and 2007, cetacean surveys were conducted by observers who kept a constant watch when the vessel was underway and recorded all cetaceans encountered. For each species, we examined distributional patterns, group size, seasonal occurrence and annual trends. Analysis of variance (anova F) was used to test for differences in group sizes between multiple means, and Student’s t-test was used to detect differences between pairwise means. Cetacean seasonal occurrence and annual trends were investigated using a generalized linear model framework. Results Humpback whales (Megaptera novaeangliae) were seen throughout the region, with numbers lowest in spring and highest in the fall. Fin whale (Balaenoptera physalus) and minke whale (Balaenoptera acutorostrata) distributions were more restricted than that reported for humpback whales, and the low number of sightings precluded evaluating seasonal trends. Three killer whale (Orcinus orca) eco-types were documented with distributions occurring throughout inland waters. Seasonal patterns were not detected or could not be evaluated for resident and offshore killer whales, respectively; however, the transient eco-type was more abundant in the summer. Dall’s porpoise (Phocoenoides dalli) were distributed throughout the region, with more sightings in spring and summer than in fall. Harbour porpoise (Phocoena phocoena) distribution was clumped, with concentrations occurring in the Icy Strait/Glacier Bay and Wrangell areas and with no evidence of seasonality. Pacific white-sided dolphins (Lagenorhynchus obliquidens) were observed only occasionally, with more sightings in the spring. For most species, group size varied on both an annual and seasonal basis. Main conclusions Seven cetacean species occupy the inland waters of Southeast Alaska, with distribution, group size, seasonal occurrence and annual trends varying by species. Future studies that compare spatial and temporal patterns with other features (e.g. oceanography, prey resources) may help in identifying the key factors that support the high density and biodiversity of cetaceans found in this region. An increased understanding of the region’s marine ecology is an essential step towards ensuring the long-term conservation of cetaceans in Southeast Alaska.

Spatial methods for plot-based sampling of wildlife populations

Classical sampling methods can be used to estimate the mean of a finite or infinite population. Block kriging also estimates the mean, but of an infinite population in a continuous spatial domain. In this paper, I consider a finite population version of block kriging (FPBK) for plot-based sampling. The data are assumed to come from a spatial stochastic process. Minimizing mean-squared-prediction errors yields best linear unbiased predictions that are a finite population version of block kriging. FPBK has versions comparable to simple random sampling and stratified sampling, and includes the general linear model. This method has been tested for several years for moose surveys in Alaska, and an example is given where results are compared to stratified random sampling. In general, assuming a spatial model gives three main advantages over classical sampling: (1) FPBK is usually more precise than simple or stratified random sampling, (2) FPBK allows small area estimation, and (3) FPBK allows nonrandom sampling designs.

Long-term recovery patterns of arctic tundra after winter seismic exploration

In response to the increasing global demand for energy, oil exploration and development are expanding into frontier areas of the Arctic, where slow-growing tundra vegetation and the underlying permafrost soils are very sensitive to disturbance. The creation of vehicle trails on the tundra from seismic exploration for oil has accelerated in the past decade, and the cumulative impact represents a geographic footprint that covers a greater extent of Alaska’s North Slope tundra than all other direct human impacts combined. Seismic exploration for oil and gas was conducted on the coastal plain of the Arctic National Wildlife Refuge, Alaska, USA, in the winters of 1984 and 1985. This study documents recovery of vegetation and permafrost soils over a two-decade period after vehicle traffic on snow-covered tundra. Paired permanent vegetation plots (disturbed vs. reference) were monitored six times from 1984 to 2002. Data were collected on percent vegetative cover by plant species and on soil and ground ice characteristics. We developed Bayesian hierarchical models, with temporally and spatially autocorrelated errors, to analyze the effects of vegetation type and initial disturbance levels on recovery patterns of the different plant growth forms as well as soil thaw depth. Plant community composition was altered on the trails by species-specific responses to initial disturbance and subsequent changes in substrate. Long-term changes included increased cover of graminoids and decreased cover of evergreen shrubs and mosses. Trails with low levels of initial disturbance usually improved well over time, whereas those with medium to high levels of initial disturbance recovered slowly. Trails on ice-poor, gravel substrates of riparian areas recovered better than those on ice-rich loamy soils of the uplands, even after severe initial damage. Recovery to pre-disturbance communities was not possible where trail subsidence occurred due to thawing of ground ice. Previous studies of disturbance from winter seismic vehicles in the Arctic predicted short-term and mostly aesthetic impacts, but we found that severe impacts to tundra vegetation persisted for two decades after disturbance under some conditions. We recommend management approaches that should be used to prevent persistent tundra damage.