Space-time modeling of environmental monitoring data

This study describes and applies statistical methods for space-time modeling of data from environmental monitoring programs, e.g., within areas such as climate change, air pollution and aquatic environment. Such data are often characterized by sparse sampling in both the temporal and spatial dimensions. In order to improve the amount of information on the physical system in question we suggest using statistical modeling methods for monitoring data. Model predictions combined with observations could be analyzed directly to assess the environmental state or as forcing functions for time series models and deterministic, hydrodynamic models. To illustrate the approach we applied the proposed modeling methods to data from the Danish and Swedish marine monitoring programs. Time series with a weekly resolution were predicted from observations of dissolved inorganic nitrogen (DIN) from the Kattegat basin (1993–1997). DIN observations were sparse, irregularly distributed and comprised approximately 10% of the generated time series.     

A process-convolution approach to modelling temperatures in the North Atlantic Ocean

This paper develops a process-convolution approach for space-time modelling. With this approach, a dependent process is constructed by convolving a simple, perhaps independent, process. Since the convolution kernel may evolve over space and time, this approach lends itself to specifying models with non-stationary dependence structure. The model is motivated by an application from oceanography: estimation of the mean temperature field in the North Atlantic Ocean as a function of spatial location and time. The large amount of this data poses some difficulties; hence computational considerations weigh heavily in some modelling aspects. A Bayesian approach is taken here which relies on Markov chain Monte Carlo for exploring the posterior distribution.     

Hierarchical Bayesian Spatial Models for Multispecies Conservation Planning and Monitoring

Abstract: Biologists who develop and apply habitat models are often familiar with the statistical challenges posed by their data's spatial structure but are unsure of whether the use of complex spatial models will increase the utility of model results in planning. We compared the relative performance of nonspatial and hierarchical Bayesian spatial models for three vertebrate and invertebrate taxa of conservation concern (Church's sideband snails [Monadenia churchi], red tree voles [Arborimus longicaudus], and Pacific fishers [Martes pennanti pacifica]) that provide examples of a range of distributional extents and dispersal abilities. We used presence–absence data derived from regional monitoring programs to develop models with both landscape and site‐level environmental covariates. We used Markov chain Monte Carlo algorithms and a conditional autoregressive or intrinsic conditional autoregressive model framework to fit spatial models. The fit of Bayesian spatial models was between 35 and 55% better than the fit of nonspatial analogue models. Bayesian spatial models outperformed analogous models developed with maximum entropy (Maxent) methods. Although the best spatial and nonspatial models included similar environmental variables, spatial models provided estimates of residual spatial effects that suggested how ecological processes might structure distribution patterns. Spatial models built from presence–absence data improved fit most for localized endemic species with ranges constrained by poorly known biogeographic factors and for widely distributed species suspected to be strongly affected by unmeasured environmental variables or population processes. By treating spatial effects as a variable of interest rather than a nuisance, hierarchical Bayesian spatial models, especially when they are based on a common broad‐scale spatial lattice (here the national Forest Inventory and Analysis grid of 24 km² hexagons), can increase the relevance of habitat models to multispecies conservation planning.     

Hierarchical modeling for extreme values observed over space and time

We propose a hierarchical modeling approach for explaining a collection of spatially referenced time series of extreme values. We assume that the observations follow generalized extreme value (GEV) distributions whose locations and scales are jointly spatially dependent where the dependence is captured using multivariate Markov random field models specified through coregionalization. In addition, there is temporal dependence in the locations. There are various ways to provide appropriate specifications; we consider four choices. The models can be fitted using a Markov Chain Monte Carlo (MCMC) algorithm to enable inference for parameters and to provide spatio–temporal predictions. We fit the models to a set of gridded interpolated precipitation data collected over a 50-year period for the Cape Floristic Region in South Africa, summarizing results for what appears to be the best choice of model.

State space models with spatial deformation

Space deformation has been proposed to model space-time varying observation processes with non-stationary spatial covariance structure under the hypothesis of temporal stationarity. In real applications, however, the temporal stationarity assumption is inappropriate and unrealistic. In this work we propose a spatial-temporal model whose temporal trend is modeled through state space models and a spatially varying anisotropy is modeled through spatial deformation, under the Bayesian approach. A distinctive feature of our approach is the consideration of model uncertainty in an unified framework. Our model has a clear advantage over the ones proposed so far in the literature when the main objective of the study is to perform spatial interpolation for fixed points in time. Approximations of the posterior distributions of the model parameters are obtained via Markov chain Monte Carlo methods. This allows for prediction of the process values in space and time as well as handling of missing values. Two applications are presented: the first one to model concentrations of sulfur dioxide in the eastern United States and the second one to model monthly minimum temperatures in the State of Rio de Janeiro.     

Bayesian hierarchical spatio-temporal modelling of trends and future projections in the ocean wave climate with a \text{ CO }_2 regression component

Bad weather and rough seas continue to be a major cause for ship losses and is thus a significant contributor to the risk to maritime transportation. This stresses the importance of taking severe sea state conditions adequately into account, with due treatment of the uncertainties involved, in ship design and operation in order to enhance safety. Hence, there is a need for appropriate stochastic models describing the variability of sea states. These should also incorporate realistic projections of future return levels of extreme sea states, taking into account long-term trends related to climate change and inherent uncertainties. The stochastic ocean wave model presented in this paper exploits the flexible framework of Bayesian hierarchical space-time models. It allows modelling of complex dependence structures in space and time and incorporation of physical features and prior knowledge, yet at the same time remains intuitive and easily interpreted. Furthermore, by taking a Bayesian approach, the uncertainties of the model parameters are also taken into account. A regression component with $\text{ CO }_2$ as an explanatory variable has been introduced in order to extract long-term trends in the data. The model has been fitted by monthly maximum significant wave height data for an area in the North Atlantic ocean. The different components of the model will be outlined in the paper, and the results will be discussed. Furthermore, a discussion of possible extensions to the model will be given.     

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.     

Spatio-temporal log-Gaussian Cox processes for modelling wildfire occurrence: the case of Catalonia, 1994–2008

Wildfires have become one of the principal environmental problems in the Mediterranean basin. While fire plays an important role in most terrestrial plant ecosystems, the potential hazard that it represents for human lives and property has led to the application of fire exclusion policies that, in the long term, have caused severe damage, mainly due to the increase of fuel loadings in forested areas, in some forest systems. The lack of an easy solution to forest fire management highlights the importance of preventive tasks. The observed spatio-temporal pattern of wildfire occurrences may be idealized as a realization of some stochastic process. In particular, we may use a space–time point pattern approach for the analysis and inference process. We studied wildfires in Catalonia, a region in the north-east of the Iberian Peninsula, and we analyzed the spatio-temporal patterns produced by those wildfire incidences by considering the influence of covariates on trends in the intensity of wildfire locations. A total of 3,166 wildfires from 1994–2008 have been recorded. We specified spatio-temporal log-Gaussian Cox process models. Models were estimated using Bayesian inference for Gaussian Markov Random Field through the integrated nested Laplace approximation algorithm. The results of our analysis have provided statistical evidence that areas closer to humans have more human induced wildfires, areas farther have more naturally occurring wildfires. We believe the methods presented in this paper may contribute to the prevention and management of those wildfires which are not random in space or time.     

Sampling variability and estimates of density dependence: a composite-likelihood approach

It is well known that sampling variability, if not properly taken into account, affects various ecologically important analyses. Statistical inference for stochastic population dynamics models is difficult when, in addition to the process error, there is also sampling error. The standard maximum-likelihood approach suffers from large computational burden. In this paper, I discuss an application of the composite-likelihood method for estimation of the parameters of the Gompertz model in the presence of sampling variability. The main advantage of the method of composite likelihood is that it reduces the computational burden substantially with little loss of statistical efficiency. Missing observations are a common problem with many ecological time series. The method of composite likelihood can accommodate missing observations in a straightforward fashion. Environmental conditions also affect the parameters of stochastic population dynamics models. This method is shown to handle such nonstationary population dynamics processes as well. Many ecological time series are short, and statistical inferences based on such short time series tend to be less precise. However, spatial replications of short time series provide an opportunity to increase the effective sample size. Application of likelihood-based methods for spatial time-series data for population dynamics models is computationally prohibitive. The method of composite likelihood is shown to have significantly less computational burden, making it possible to analyze large spatial time-series data. After discussing the methodology in general terms, I illustrate its use by analyzing a time series of counts of American Redstart (Setophaga ruticilla) from the Breeding Bird Survey data, San Joaquin kit fox (Vulpes macrotis mutica) population abundance data, and spatial time series of Bull trout (Salvelinus confluentus) redds count data.     

Hierarchical modeling of count data with application to nuclear fall-out

We propose a method for a Bayesian hierarchical analysis of count data that are observed at irregular locations in a bounded domain of R². We model the data as having been observed on a fine regular lattice, where we do not have observations at all the sites. The counts are assumed to be independent Poisson random variables whose means are given by a log Gaussian process. In this article, the Gaussian process is assumed to be either a Markov random field (MRF) or a geostatistical model, and we compare the two models on an environmental data set. To make the comparison, we calibrate priors for the parameters in the geostatistical model to priors for the parameters in the MRF. The calibration is obtained empirically. The main goal is to predict the hidden Poisson-mean process at all sites on the lattice, given the spatially irregular count data; to do this we use an efficient MCMC. The spatial Bayesian methods are illustrated on radioactivity counts analyzed by Diggle et al. (1998).     

Modeling freshwater fish distributions using multiscale landscape data: A case study of six narrow range endemics

Species distribution models (SDMs) have become integral tools in scientific research and conservation planning. Despite progress in the assessment of various statistical models for use in SDMs, little has been done in way of evaluating appropriate ecological models. In this paper, we evaluate the multiscale filter framework as a suitable theoretical model for predicting freshwater fish distributions in the upper Green River system (Ohio River drainage), USA. The spatial distributions of six fishes with contrasting biogeographies were modeled using boosted regression trees and multiscale landscape data. Species biogeography did not appear to affect predictive performance and all models performed well statistically with receiver operating characteristic area under the curve (AUC) ranging from 0.87 to 0.98. Predictive maps show accurate estimations of ranges for five of six species based on historical collections. The relative influence of each type of environmental feature and spatial scale varied markedly with between species. A hierarchical effect was detected for narrowly distributed species. These species were highly influenced by soil composition at larger spatial scales and land use/land cover (LULC) patterns at more proximal scales. Conversely, LULC pattern was the most influential feature for widely distributed at all spatial scales. Using multiscale data capable of capturing hierarchical landscape influences allowed production of accurate predictive models and provided further insight into factors controlling freshwater fish distributions.     

Space-time statistics for environmental and agricultural related phenomena

This paper presents an overview of space-time statistical procedures to analyse agricultural and environmental related phenomena. It starts with an application on root-rot development in cotton. Dependence modelling in space and time is done with the space-time variogram. Various kriging interpolators are presented for making predictions in space and time. Simulated annealing is used to design an optimal monitoring network for estimation of space-time variograms. In the application no clear indication was found for anisotropy, although strong evidence exists that the disease not only proceeds within rows but also jumps between rows. The optimal sampling scheme showed a spatial clustering of observations at the first and the last monitoring day and less observations at intermediate times.     

A single-chain-based multidimensional Markov chain model for subsurface characterization

Multidimensional Markov chain models in geosciences were often built on multiple chains, one in each direction, and assumed these 1-D chains to be independent of each other. Thus, unwanted transitions (i.e., transitions of multiple chains to the same location with unequal states) inevitably occur and have to be excluded in estimating the states at unobserved locations. This consequently may result in unreliable estimates, such as underestimation of small classes (i.e., classes with smaller than average areas) in simulated realizations. This paper presents a single-chain-based multidimensional Markov chain model for estimation (i.e., prediction and conditional stochastic simulation) of spatial distribution of subsurface formations with borehole data. The model assumes that a single Markov chain moves in a lattice space, interacting with its nearest known neighbors through different transition probability rules in different cardinal directions. The conditional probability distribution of the Markov chain at the location to be estimated is formulated in an explicit form by following the Bayes’ Theorem and the conditional independence of sparse data in cardinal directions. Since no unwanted transitions are involved, the model can estimate all classes fairly. Transiogram models (i.e., 1-D continuous Markov transition probability diagrams) are used to provide transition probability input with needed lags to generalize the model. Therefore, conditional simulation can be conducted directly and efficiently. The model provides an alternative for heterogeneity characterization of subsurface formations.

Hierarchical modeling of abundance in closed population capture–recapture models under heterogeneity

Hierarchical modeling of abundance in space or time using closed-population mark-recapture under heterogeneity (model \(\hbox {M}_{\text {h}}\) ) presents two challenges: (i) finding a flexible likelihood in which abundance appears as an explicit parameter and (ii) fitting the hierarchical model for abundance. The first challenge arises because abundance not only indexes the population size, it also determines the dimension of the capture probabilities in heterogeneity models. A common approach is to use data augmentation to include these capture probabilities directly into the likelihood and fit the model using Bayesian inference via Markov chain Monte Carlo (MCMC). Two such examples of this approach are (i) explicit trans-dimensional MCMC, and (ii) superpopulation data augmentation. The superpopulation approach has the advantage of simple specification that is easily implemented in BUGS and related software. However, it reparameterizes the model so that abundance is no longer included, except as a derived quantity. This is a drawback when hierarchical models for abundance, or related parameters, are desired. Here, we analytically compare the two approaches and show that they are more closely related than might appear superficially. We exploit this relationship to specify the model in a way that allows us to include abundance as a parameter and that facilitates hierarchical modeling using readily available software such as BUGS. We use this approach to model trends in grizzly bear abundance in Yellowstone National Park from 1986 to 1998.     

GIS-basierte Modellierung von Gewässerimmissionen

Goal and Scope

Environmental assessment of aquatic micro pollutants should consider the spatial and temporal variability of emission, transport and transformation. Simulation models coupled with Geographic Information Systems (GIS) provide digital maps of concentration patterns caused by the overlay of multipoint and diffuse emissions and natural attenuation processes in river basins. The paper gives an overview on GIS-based models for river basins and demonstrates the applicability by using some illustrating examples with GREAT-ER.

Main Features

Georeferenced models have several advantages: visualization of concentration patterns, investigation of spatial and temporal concentration profiles, analysis of exceedance of environmental quality standards, embedding in integrated river basin management systems.

Results and Discussion

GIS-based models allow a more realistic assessment. Monitoring programmes should be designed to deliver appropriate measured data for the evaluation and improvement of models.

Recommendation and Perspectives

The combination of digital maps, simulation models and environmental monitoring would provide better approaches for the risk assessment and water quality management of aquatic micro pollutants.     

Spatial structure effects on the detection of patches boundaries using local operators

Landscapes exhibit various degrees of spatial heterogeneity according to the differential intensity and interactions among processes and disturbances that they are subjected to. The management of these spatially dynamical landscapes requires that we can accurately map them and monitor the evolution of their spatial arrangement through time. Such a mapping requires first the delineation of various spatial features present in the landscape such as patches and their boundaries. However, there are several environmental (spatial variability) as well as technical (spatial resolution) factors that impair our ability to accurately delineate patches and their boundaries as polygons. Here, we investigate how the spatial structure and spatial resolution of the data affect the accuracy of detecting patches and their boundaries over simulated landscapes and real data. Simulated landscapes consisted of two patches with parameterized spatial properties (patches’ level of spatial autocorrelation, mean value and variance) separated by a boundary of known location. Real data allowed the investigation of a more complex landscape where there is a known transition between two forest domains with unknown spatial properties. Boundary locations are defined using the lattice-wombling edge detector at various aggregation levels and the degree of patch homogeneity is determined using Getis-Ord’s G*. Results show that boundary detection using a local edge detector is greatly affected by the spatial conditions of the data, namely variance, abruptness of the spatial gradient between two patches and patches’ level of spatial autocorrelation. They also suggest that data aggregation is not a panacea for bringing out the ecological process creating the patches and that indicators derived from local measures of spatial association can be complementary tools for analysing spatial structures affecting boundary delineation.

Fragmentation profiles for real and simulated landscapes

When a natural landscape is represented by a series of categorical raster maps of varying resolution, a multiresolution characterization of spatial pattern can be obtained in which entropy is computed at each resolution conditional on the next coarser resolution. The series of entropy values is plotted as a function of resolution, resulting in a multiresolution profile of fragmentation pattern in the landscape. If a categorical raster map is available at a single resolution only, a series of degraded maps at increasingly coarser resolutions is generated and the fragmentation profile is computed for this series. An algorithm has been developed for obtaining the profile directly from the single resolution map without having to generate and store the coarser resolution maps. A hierarchical stochastic model is described for simulating categorical raster maps and the fragmentation profile of the generating process is obtained in terms of the model parameters. These process profiles provide benchmarks for assessing empirical profiles obtained from raster maps of actual landscapes. Methods of the paper are applied to several watersheds of Pennsylvania using landcover maps derived from satellite imagery. These examples indicate that characteristic landscape types induce characteristic features in their fragmentation profiles.