Comparing environmental influences on coral bleaching across and within species using clustered binomial regression

Differential susceptibility among reef-building coral species can lead to community shifts and loss of diversity as a result of temperature-induced mass bleaching events. We evaluate environmental influences on coral colony bleaching over an 8-year period in the Florida Keys, USA. Clustered binomial regression is used to develop models incorporating taxon-specific responses to the environment in order to identify conditions and species for which bleaching is likely to be severe. By building three separate models incorporating environment, community composition, and taxon-specific responses to environment, we show observed prevalence of bleaching reflects an interaction between community composition and local environmental conditions. Environmental variables, including elevated sea temperature, solar radiation, and reef depth, explained 90% and 78% of variability in colony bleaching across space and time, respectively. The effects of environmental variables were only partially explained (33% of variability) by corresponding differences in community composition. Taxon-specific models indicated individual coral species responded differently to local environmental conditions and had different sensitivities to temperature-induced bleaching. For many coral species, but not all, bleaching was exacerbated by high solar radiation. A 25% reduction in the probability of bleaching in shallow locations for one species may reflect an ability to acclimatize to local conditions. Overall, model results indicate predictions of coral bleaching require knowledge of not just the environmental conditions or community composition, but the responses of individual species to the environment. Model development provides a useful tool for coral reef management by quantifying the influence of the local environment on individual species bleaching sensitivities, identifying susceptible species, and predicting the likelihood of mass bleaching events with changing environmental conditions.     

Prospective spatial prediction of infectious disease: experience of New York State (USA) with West Nile Virus and proposed directions for improved surveillance

Infectious disease surveillance has become an international top priority due to the perceived risk of bioterrorism. This is driving the improvement of real-time geo-spatial surveillance systems for monitoring disease indicators, which is expected to have many benefits beyond detecting a bioterror event. West Nile Virus surveillance in New York State (USA) is highlighted as a working system that uses dead American Crows (Corvus brachyrhynchos) to prospectively indicate viral activity prior to human onset. A cross-disciplinary review is then presented to argue that this system, and infectious disease surveillance in general, can be improved by complementing spatial cluster detection of an outcome variable with predictive “risk mapping” that incorporates spatiotemporal data on the environment, climate and human population through the flexible class of generalized linear mixed models.

Designing environmental monitoring networks to measure extremes

This paper discusses challenges arising in the design of networks for monitoring extreme values over the domain of a random environmental space-time field {X _ij } i = 1, . . . , I denoting site and j = 1, . . . denoting time (e.g. hour). The field of extremes for time span r over site domain i = 1, . . . ,I is given by \(\{Y_{i(r+1)}=\max_{j=k}^{k+n-1} X_{ij}\}\) for k = 1 + rn, r = 0, . . . ,. Such networks must not only measure extremes at the monitored sites but also enable their prediction at the non-monitored ones. Designing such a network poses special challenges that do not seem to have been generally recognized. One of these problems is the loss of spatial dependence between site responses in going from the environmental process to the field of extremes it generates. In particular we show empirically that the intersite covariance Cov(Y _i(r+1),Y _i′(r+1)) can generally decline toward zero as r increases, for site pairs i ≠ i′. Thus the measured extreme values may not predict the unmeasured ones very precisely. Consequently high levels of pollution exposure of a sensitive group (e.g. school children) located between monitored sites may be overlooked. This potential deficiency raises concerns about the adequacy of air pollution monitoring networks whose primary role is the detection of noncompliance with air quality standards based on extremes designed to protect human health. The need to monitor for noncompliance and thereby protect human health, points to other issues. How well do networks designed to monitor the field monitor their fields of extremes? What criterion should be used to select prospective monitoring sites when setting up or adding to a network? As the paper demonstrates by assessing an existing network, the answer to the first question is not well, at least in the case considered. To the second, the paper suggests a variety of plausible answers but shows through a simulation study, that they can lead to different optimum designs. The paper offers an approach that circumvents the dilemma posed by the answer to the second question. That approach models the field of extremes (suitably transformed) by a multivariate Gaussian-Inverse Wishart hierarchical Bayesian distribution. The adequacy of this model is empirically assessed in an application by finding the relative coverage frequency of the predictive credibility ellipsoid implied by its posterior distribution. The favorable results obtained suggest this posterior adequately describes that (transformed) field. Hence it can form the basis for designing an appropriate network. Its use is demonstrated by a hypothetical extension of an existing monitoring network. That foundation in turn enables a network to be designed of sufficient density (relative to cost) to serve its regulatory purpose.     

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     

Integrated models to support multiobjective ecological restoration decisions

Many objectives motivate ecological restoration, including improving vegetation condition, increasing the range and abundance of threatened species, and improving species richness and diversity. Although models have been used to examine the outcomes of ecological restoration, few researchers have attempted to develop models to account for multiple, potentially competing objectives. We developed a combined state‐and‐transition, species‐distribution model to predict the effects of restoration actions on vegetation condition and extent, bird diversity, and the distribution of several bird species in southeastern Australian woodlands. The actions reflected several management objectives. We then validated the models against an independent data set and investigated how the best management decision might change when objectives were valued differently. We also used model results to identify effective restoration options for vegetation and bird species under a constrained budget. In the examples we evaluated, no one action (improving vegetation condition and extent, increasing bird diversity, or increasing the probability of occurrence for threatened species) provided the best outcome across all objectives. In agricultural lands, the optimal management actions for promoting the occurrence of the Brown Treecreeper (Climacteris picumnus), an iconic threatened species, resulted in little improvement in the extent of the vegetation and a high probability of decreased vegetation condition. This result highlights that the best management action in any situation depends on how much the different objectives are valued. In our example scenario, no management or weed control were most likely to be the best management options to satisfy multiple restoration objectives. Our approach to exploring trade‐offs in management outcomes through integrated modeling and structured decision‐support approaches has wide application for situations in which trade‐offs exist between competing conservation objectives.     

Effects of surrounding land use and water depth on seagrass dynamics relative to a catastrophic algal bloom

Seagrasses are the foundation of many coastal ecosystems and are in global decline because of anthropogenic impacts. For the Indian River Lagoon (Florida, U.S.A.), we developed competing multistate statistical models to quantify how environmental factors (surrounding land use, water depth, and time [year]) influenced the variability of seagrass state dynamics from 2003 to 2014 while accounting for time‐specific detection probabilities that quantified our ability to determine seagrass state at particular locations and times. We classified seagrass states (presence or absence) at 764 points with geographic information system maps for years when seagrass maps were available and with aerial photographs when seagrass maps were not available. We used 4 categories (all conservation, mostly conservation, mostly urban, urban) to describe surrounding land use within sections of lagoonal waters, usually demarcated by land features that constricted these waters. The best models predicted that surrounding land use, depth, and year would affect transition and detection probabilities. Sections of the lagoon bordered by urban areas had the least stable seagrass beds and lowest detection probabilities, especially after a catastrophic seagrass die‐off linked to an algal bloom. Sections of the lagoon bordered by conservation lands had the most stable seagrass beds, which supports watershed conservation efforts. Our results show that a multistate approach can empirically estimate state‐transition probabilities as functions of environmental factors while accounting for state‐dependent differences in seagrass detection probabilities as part of the overall statistical inference procedure.     

Optimizing conservation in species-specific agricultural landscapes

Recovery of grassland birds in agricultural landscapes is a global imperative. Agricultural landscapes are complex, and the value of resource patches may vary substantially among species. The spatial extent at which landscape features affect populations (i.e., scale of effect) may also differ among species. There is a need for regional-scale conservation planning that considers landscape-scale and species-specific responses of grassland birds to environmental change. We developed a spatially explicit approach to optimizing grassland conservation in the context of species-specific landscapes and prioritization of species recovery and applied it to a conservation program in Kentucky (USA). We used a hierarchical distance-sampling model with an embedded scale of effect predictor to estimate the relationship between landscape structure and abundance of eastern meadowlarks (Sturnella magna), field sparrows (Spizella pusilla), and northern bobwhites (Colinus virginianus). We used a novel spatially explicit optimization procedure rooted in multi-attribute utility theory to design alternative conservation strategies (e.g., prioritize only northern bobwhite recovery or assign equal weight to each species’ recovery). Eastern meadowlarks and field sparrows were more likely to respond to landscape-scale resource patch adjacencies than landscape-scale patch densities. Northern bobwhite responded to both landscape-scale resource patch adjacencies and densities and responded strongly to increased grassland density. Effects of landscape features on local abundance decreased as distance increased and had negligible influence at 0.8 km for eastern meadowlarks (0.7–1.2 km 95% Bayesian credibility intervals [BCI]), 2.5 km for field sparrows (1.5–5.8 km 95% BCI), and 8.4 km for bobwhite (6.4–26 km 95% BCI). Northern bobwhites were predicted to benefit greatly from future grassland conservation regardless of conservation priorities, but eastern meadowlark and field sparrow were not. Our results suggest similar species can respond differently to broad-scale conservation practices because of species-specific, distance-dependent relationships with landscape structure. Our framework is quantitative, conceptually simple, customizable, and predictive and can be used to optimize conservation in heterogeneous ecosystems while considering landscape-scale processes and explicit prioritization of species recovery.     

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.     

Accident modelling and analysis in process industries

Accident modelling is a methodology used to relate the causes and effects of events that lead to accidents. This modelling effectively seeks to answer two main questions: (i) Why does an accident occur, and (ii) How does it occur. This paper presents a review of accident models that have been developed for the chemical process industry with in-depth analyses of a class of models known as dynamic sequential accident models (DSAMs). DSAMs are sequential models with a systematic procedure to utilise precursor data to estimate the posterior risk profile quantitatively. DSAM also offers updates on the failure probabilities of accident barriers and the prediction of future end states. Following a close scrutiny of these methodologies, several limitations are noted and discussed, and based on these insights, future work is suggested to enhance and improve this category of models further.     

Modelling fertility: an application of count regression models

Often the lifecycle data occur as count of the vital events and are recorded as integers. The purpose of this article is to model the fertility behavior based on religious, educational, economic, and occupational characteristics. The responses of classified groups according to these determinants are examined for significant influence on fertility using Poisson regression model (PRM) based on the National Family Health Survey -3 dataset. The observed and predicted probabilities under PRM indicate modal value of two children for the Poisson distribution modeled data. Presence of dominance of two child in the data motivates the authors to adopt multinomial regression model (MRM) in order to link fertility with various socioeconomic indicators responsible for fertility variation. Choice of the explanatory factors is limited to the availability of data. Trends and patterns of preference for birth counts suggest that religion, caste, wealth, female education, and occupation are the dominant factors shaping the observed birth process. Empirical analysis suggests that both the models used in the study perform similarly on the sample data. However, fitting of MRM by taking birth count of two as comparison category shows improved Akaike information criterion and consistent Akaike information criterion values. Current work contributes to the existing literature as it attempts to provide more insight into the determinants of Indian fertility using Poisson and MRM.