Patterns of Spatial Autocorrelation in Stream Water Chemistry

Geostatistical models are typically based on symmetric straight-line distance, which fails to represent the spatial configuration, connectivity, directionality, and relative position of sites in a stream network. Freshwater ecologists have explored spatial patterns in stream networks using hydrologic distance measures and new geostatistical methodologies have recently been developed that enable directional hydrologic distance measures to be considered. The purpose of this study was to quantify patterns of spatial correlation in stream water chemistry using three distance measures: straight-line distance, symmetric hydrologic distance, and weighted asymmetric hydrologic distance. We used a dataset collected in Maryland, USA to develop both general linear models and geostatistical models (based on the three distance measures) for acid neutralizing capacity, conductivity, pH, nitrate, sulfate, temperature, dissolved oxygen, and dissolved organic carbon. The spatial AICC methodology allowed us to fit the autocorrelation and covariate parameters simultaneously and to select the model with the most support in the data. We used the universal kriging algorithm to generate geostatistical model predictions. We found that spatial correlation exists in stream chemistry data at a relatively coarse scale and that geostatistical models consistently improved the accuracy of model predictions. More than one distance measure performed well for most chemical response variables, but straight-line distance appears to be the most suitable distance measure for regional geostatistical modeling. It may be necessary to develop new survey designs that more fully capture spatial correlation at a variety of scales to improve the use of weighted asymmetric hydrologic distance measures in regional geostatistical models.