Spatial heterogeneity of soil nutrients and plant species in herb-dominated communities of contrasting land use

Recent interest in spatial pattern in terrestrial ecosystems has come from an awareness of theintimate relationship between spatial heterogeneity of soil resources and maintenance of plant species diversity. Soil and vegetation can vary spatially inresponse to several state factors of the system. In this study, we examined fine-scale spatial variability of soil nutrients and vascular plant species in contrasting herb-dominated communities (a pasture and an oldfield) to determine degree of spatial dependenceamong soil variables and plant community characteristics within these communities by sampling at 1-m intervals. Each site was divided into 25 1-m 2 plots. Mineral soil was sampled (2-cm diameter, 5-cm depth) from each of four 0.25-m2 quarters and combined into a single composite sample per plot. Soil organic matter was measured as loss-on-ignition. Extractable NH4 and NO3 were determined before and after laboratory incubation to determine potential net N mineralization and nitrification. Cations were analyzed using inductively coupled plasma emission spectrometry. Vegetation was assessed using estimated percent cover. Most soiland plant variables exhibited sharp contrasts betweenpasture and old-field sites, with the old field having significantly higher net N mineralization/nitrification, pH, Ca, Mg, Al, plant cover, and species diversity, richness, and evenness. Multiple regressions revealedthat all plant variables (species diversity, richness,evenness, and cover) were significantly related to soil characteristics (available nitrogen, organic matter,moisture, pH, Ca, and Mg) in the pasture; in the old field only cover was significantly related to soil characteristics (organic matter and moisture). Both sites contrasted sharply with respect to spatial pattern of soil variables, with the old field exhibiting a higher degree of spatial dependence. These results demonstrate that land-use practices can exert profound influence on spatial heterogeneity of both soil properties and vegetation in herb-dominated communities.

Computational Techniques for Spatial Logistic Regression with Large Datasets

In epidemiological work, outcomes are frequently non-normal, sample sizes may be large, and effects are often small. To relate health outcomes to geographic risk factors, fast and powerful methods for fitting spatial models, particularly for non-normal data, are required. We focus on binary outcomes, with the risk surface a smooth function of space. We compare penalized likelihood models, including the penalized quasi-likelihood (PQL) approach, and Bayesian models based on fit, speed, and ease of implementation. A Bayesian model using a spectral basis representation of the spatial surface provides the best tradeoff of sensitivity and specificity in simulations, detecting real spatial features while limiting overfitting and being more efficient computationally than other Bayesian approaches. One of the contributions of this work is further development of this underused representation. The spectral basis model outperforms the penalized likelihood methods, which are prone to overfitting, but is slower to fit and not as easily implemented. Conclusions based on a real dataset of cancer cases in Taiwan are similar albeit less conclusive with respect to comparing the approaches. The success of the spectral basis with binary data and similar results with count data suggest that it may be generally useful in spatial models and more complicated hierarchical models.

GEOSTATISTICAL INFERENCE UNDER PREFERENTIAL SAMPLING

Geostatistics involves the fitting of spatially continuous models to spatially discrete data (Chil`es and Delfiner, 1999). Preferential sampling arises when the process that determines the data-locations and the process being modelled are stochastically dependent. Conventional geostatistical methods assume, if only implicitly, that sampling is non-preferential. However, these methods are often used in situations where sampling is likely to be preferential. For example, in mineral exploration samples may be concentrated in areas thought likely to yield high-grade ore. We give a general expression for the likelihood function of preferentially sampled geostatistical data and describe how this can be evaluated approximately using Monte Carlo methods. We present a model for preferential sampling, and demonstrate through simulated examples that ignoring preferential sampling can lead to seriously misleading inferences. We describe an application of the model to a set of bio-monitoring data from Galicia, northern Spain, in which making allowance for preferential sampling materially changes the inferences.

Bayesian Geostatistical Design

This paper describes the use of model-based geostatistics for choosing the optimal set of sampling locations, collectively called the design, for a geostatistical analysis. Two types of design situations are considered. These are retrospective design, which concerns the addition of sampling locations to, or deletion of locations from, an existing design, and prospective design, which consists of choosing optimal positions for a new set of sampling locations. We propose a Bayesian design criterion which focuses on the goal of efficient spatial prediction whilst allowing for the fact that model parameter values are unknown. The results show that in this situation a wide range of inter-point distances should be included in the design, and the widely used regular design is therefore not the optimal choice.

BIVARIATE BINOMIAL SPATIAL MODELLING LOA loa PREVALENCE IN TROPICAL AFRICA

We present a state-of-the-art application of smoothing for dependent bivariate binomial spatial data to Loa loa prevalence mapping in West Africa. This application is special because it starts with the non-spatial calibration of survey instruments, continues with the spatial model building and assessment and ends with robust, tested software that will be used by the field scientists of the World Health Organization for online prevalence map updating. From a statistical perspective several important methodological issues were addressed: (a) building spatial models that are complex enough to capture the structure of the data but remain computationally usable; (b)reducing the computational burden in the handling of very large covariate data sets; (c) devising methods for comparing spatial prediction methods for a given exceedance policy threshold.