Abstract
We develop a method for multiscale estimation of pollutant concentrations, based on a nonparametric spatial statistical model. We apply this method to estimate nitrate concentrations in groundwater over the mid-Atlantic states, using measurements gathered during a period of 10 years. A map of the fine-scale estimated nitrate concentration is obtained, as well as maps of the estimated county-level average nitrate concentration and similar maps at the level of watersheds and other geographic regions. The fine-scale and coarse-scale estimates arise naturally from a single model, without refitting or ad hoc aggregation. As a result, the uncertainty associated with each estimate is available, without approximations relying on high spatial density of measurements or parametric distributional assumptions.
Several risk measures are also obtained, including the probability of the pollutant concentration exceeding a particular threshold. These risk measures can be obtained at the fine scale, or at the level of counties or other regions.
The nonparametric Bayesian statistical model allows for this flexibility in estimation while avoiding strong assumptions. This method can be applied directly to estimate ozone concentrations in air, pesticide concentrations in groundwater, or any other quantity that varies over a geographic region, based on approximate measurements at some locations and perhaps of associated covariates. An S-PLUS package with this capability is provided as supplemental material.
Similar content being viewed by others
References
Abramowitz, M., and Stegun, I. A. (eds.) (1964), Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series, Vol. 55, Washington, DC: National Bureau of Standards.
Ator, S. W. (1998), “Nitrate and Pesticide Data for Waters of the Mid-Atlantic Region,” USGS Open File Report 98-158, U.S. Geological Survey, Reston, VA.
Ator, S. W., and Denis, J. M. (1997), “Relation of Nitrogen and Phosphorus in Ground Water to Land Use in Four Subunits of the Potomac River Basin,” USGS Water-Resources Investigations Report 97-4268, U.S. Geological Survey, Reston, VA.
Ator, S. W., and Ferrari, M. J. (1997), “Nitrate and Selected Pesticides in Ground Water of the Mid-Atlantic Region,” USGS Water-Resources Investigations Report 97-4139, U.S. Geological Survey, Reston, VA.
Besag, J., York, J., and Mollié, A. (1991), “Bayesian Image Restoration, With Two Applications in Spatial Statistics” (with comments), Annals of the Institute of Statistical Mathematics, 43, 1–59.
Best, N. G., Ickstadt, K., and Wolpert, R. L. (2000), “Spatial Poisson Regression for Health and Exposure Data Measured at Disparate Resolutions,” Journal of the American Statistical Association, 95, 1076–1088.
Chilès, J.-P., and Delfiner, P. (eds.) (1999), Geostatistics, Modeling Spatial Uncertainty, New York: Wiley.
Clyde, M. A., House, L. L., and Wolpert, R. L. (2006), “Nonparametric Models for Proteomic Peak Identification and Quantification,” in Bayesian Inference for Gene Expression and Proteomics, eds. K. A. Do, P. Muller, and M. Vannucci, Cambridge, U.K.: Cambridge University Press, pp. 293–308.
Cressie, N. (1993), Statistics for Spatial Data, New York: Wiley.
Cressie, N., and Chan, N. H. (1989), “Spatial Modeling of Regional Variables,” Journal of the American Statistical Association, 84, 393–401.
Faulkner, B. R. (2003), “Confronting the Modifiable Areal Unit Problem for Inference on Nitrate in Regional Shallow Ground Water,” in Groundwater Quality Modeling and Management Under Uncertainty, ed. S. Mishra, Reston, VA: American Society of Civil Engineers, pp. 248–259.
Fields Development Team (2004), Fields: Tools for Spatial Data, Boulder, CO: National Center for Atmospheric Research. Available at http://www.cgd.ucar.edu/stats/Software/Fields/ .
Gilks, W. R., Richardson, S., and Spiegelhalter, D. J. (eds.) (1996), Markov Chain Monte Carlo in Practice, New York: Chapman & Hall.
Green, P. J. (1995), “Reversible Jump Markov Chain Monte Carlo Computation and Bayesian Model Determination,” Biometrika, 82, 711–732.
Hamilton, P. A., Denver, J. M., Phillips, P. J., and Shedlock, R. J. (1993), “Water-Quality Assessment of the Delmarva Peninsula, Delaware, Maryland, and Virginia—Effects of Agricultural Activities on, and Distribution of, Nitrate and Other Inorganic Constituents in the Surficial Aquifer,” USGS Open File Report 93-40, U.S. Geological Survey, Reston, VA.
House, L. L., Clyde, M. A., and Wolpert, R. L. (2006), “Nonparametric Models for Peak Identification and Quantification in Mass Spectroscopy, With Application to MALDI-TOF,” Discussion Paper 2006-24, Duke University, Dept. of Statistical Science, available at ftp://ftp.isds.duke.edu/pub/WorkingPapers/06-24.html .
Ickstadt, K., and Wolpert, R. L. (1997), “Multiresolution Assessment of Forest Inhomogeneity,” in Case Studies in Bayesian Statistics, Vol. III, eds. C. Gatsonis, J. S. Hodges, R. E. Kass, R. E. McCulloch, P. Rossi, and N. D. Singpurwalla, New York: Springer-Verlag, pp. 371–386.
— (1999), “Spatial Regression for Marked Point Processes” (with comments), in Bayesian Statistics, Vol. 6, eds. J. M. Bernardo, J. O. Berger, A. P. Dawid, and A. F. M. Smith, Oxford: Oxford University Press, pp. 323–341.
Jones, G. L., Haran, M., Caffo, B. S., and Neath, R. (2006), “Fixed-Width Output Analysis for Markov Chain Monte Carlo,” Journal of the American Statistical Association, 101, 1537–1547.
LaMotte, A. E., and Greene, E. A. (2007), “Spatial Analysis of Land Use and Shallow Groundwater Vulnerability in the Watershed Adjacent to Assateague Island National Seashore, Maryland and Virginia, USA,” Environmental Geology, 52, 1413–1421.
Lophaven, S. N., Nielsen, H. B., and Søndergaard, J. (2002), “DACE: A Matlab Kriging Toolbox, Version 2.0,” Technical Report IMM-TR-2002-12, Technical University of Denmark, available at http://www.imm.dtu.dk/~hbn/dace/dace.pdf .
Nakaya, T. (2000), “An Information Statistical Approach to the Modifiable Areal Unit Problem in Incidence Rate Maps,” Environment and Planning A, 32, 91–109.
Nolan, B. T., Hitt, K. J., and Ruddy, B. C. (2002), “Probability of Nitrate Contamination of Recently Recharged Groundwaters in the Conterminous United States,” Environmental Science and Technology, 36, 2138–2145.
Roberts, G. O., Gelman, A., and Gilks, W. R. (1997), “Weak Convergence and Optimal Scaling of Random Walk Metropolis Algorithms,” The Annals of Applied Probability, 7, 110–120.
S-Plus (2007), S-PLUS 8 Programmer’s Guide, Seattle, WA: Insightful Corporation.
Stein, M. L. (1999), Interpolation of Spatial Data: Some Theory for Kriging, New York: Springer-Verlag.
Tierney, L. (1994), “Markov Chains for Exploring Posterior Distributions” (with discussion), The Annals of Statistics, 22, 1701–1762.
Tu, C. (2006), “Bayesian Nonparametric Modeling Using Lévy Process Priors With Applications for Function Estimation, Time Series Modeling, and Spatio-Temporal Modeling,” Ph.D. thesis, Duke University, Dept. of Statistical Science.
U. S. Environmental Protection Agency (1991), “Fact Sheet: National Primary Drinking Water Standards,” U.S. Government Printing Office, Washington, DC.
Wolpert, R. L., and Ickstadt, K. (1998a), “Poisson/Gamma Random Field Models for Spatial Statistics,” Biometrika, 85, 251–267.
— (1998b), “Simulation of Lévy Random Fields,” New York: Springer-Verlag.
Wolpert, R. L., Clyde, M. A., and Tu, C. (2006), “Lévy Adaptive Regression Kernels,” Discussion Paper 2006-08, Duke University, Dept. of Statistical Science, available at http://ftp.stat.duke.edu/WorkingPapers/06-08.html .
Author information
Authors and Affiliations
Corresponding author
Electronic Supplementary Material
Rights and permissions
About this article
Cite this article
Woodard, D.B., Wolpert, R.L. & O’Connell, M.A. Spatial Inference of Nitrate Concentrations in Groundwater. JABES 15, 209–227 (2010). https://doi.org/10.1007/s13253-009-0006-x
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s13253-009-0006-x