Abstract
The steady state two dimensional groundwater flow equation with constant transmissivities was studied by Whittle in 1954 as a stochastic Laplace equation. He showed that the correlation function consisted of a modified Bessel function of the second kind, order 1, multiplied by its argument. This paper uses this pioneering work of Whittle to fit an aquifer head field to unequally spaced observations by maximum likelihood. Observational error is also included in the model. Both the isotropic and anisotropic cases are considered. The fitted field is then calculated on a two dimensional grid together with its standard deviation. The method is closely related to the use of two-dimensional splines for fitting surfaces to irregularly spaced observations.
Similar content being viewed by others
References
Ansley, C. F. 1979: An algorithm for the exact likelihood of a mixed autoregressive-moving average process. Biometrika 66, 59–65
Bodo, B.A.; Thompson, M.E.; Unny, T.E. 1978: A review on stochastic differential equations for applications in hydrology. Stochastic Hydrol. Hydraul. 1, 81–100
Dennis, J. E.; Schnabel, R. B. 1983: Numerical Methods for Unconstrained Optimization and Nonlinear Equations. Englewood Cliffs, New Jersey: Prentice-Hall
Graybill, F. A. 1976: Theory and Application of the Linear Model. North Scituate, Massachusetts: Duxbury Press
Jones, R. H. 1980: Maximum likelihood fitting of ARMA models to time series with missing observations. Technometrics 22, 389–395
Jones, R. H. 1981: Fitting continuous time autoregressions to discrete data. In: Findley, D. F. (ed.) Applied Time Series Analysis II, pp. 651–632 New York: Academic Press
Jones, R.H. 1987: Serial correlation in unbalanced mixed models. Bulletin of the International Statistical Institute, Proceedings of the 46th Session, Tokyo, 8–16 September, 1987, Book 4, 105–122.
Jones, R. H.; Kon. J. 1985: Two dimensional linear filtering and prediction applied to aquifer height data. In: Preprint, Rocky Mountain Conference on Medical Applications of Statistics, University of Colorado School of Medicine, Box B-119, Denver, CO 890262, 31–38
Jones, R. H.; Tryon, P. V. 1987: Continuous time series models for unequally spaced data applied to modeling atomic clocks. SIAM J. Sci. Stat. Comput. 8, 71–81
Kalman, R. E. 1960: A new approach to linear filtering and prediction problems.. Tran. ASME (J. Basic Engineering) 82D, 35–45
Kalman, R. E.; Bucy, R. S. 1961. New results in linear filtering and prediction theory. Tran. ASME (J. Basic Engineering) 83D, 95–108
Lenfest, L. W., Jr. 1986: Ground-water levels and use of water for irrigation in the Saratoga Valley, southcentral Wyoming, 1980–81. Water-Resources Investigation Report 84-4040. Cheyenne: U.S. Geological Survey.
Mardia, K. V.; Marshall, R. J. 1984: Maximum likelihood estimation for models for residual covariance in spatial regression. Biometrika 71, 135–146
Serrano, S.E. 1988a: General solution to a random advective-dispersive equation in porous media. Part 1: Stochasticity in the sources and in the boundaries. Stochastic Hydrol. Hydraul. 2, 79–98
Serrano, S.E. 1988b: General solution to a random advective-dispersive equation in porous media. Part 2: Stochasticity in the parameters. Stochastic Hydrol. Hydraul. 2, 99–112
Serrano, S.E.; Unny, T.E. 1987: Predicting groundwater flow in a phreatic aquifer. J. of Hydrol. 95, 241–268
Serrano, S.E.; Unny, T.E.; Lennox, W.C. 1985a: Analysis of stochastic groundwater problems. Part I: Deterministic partial differential equations in groundwater flow. A functional-analytical approach. J. of Hydrol. 82, 247–263
Serrano, S.E.; Unny, T.E.; Lennox, W.C. 1985b: Analysis of stochastic groundwater problems. Part II: Stochastic partial differential equations in groundwater flow. A functional-analytical approach. J. of Hydrol. 82, 265–284
Serrano, S.E.; Unny, T.E.; Lennox, W.C. 1985c: Analysis of stochastic groundwater problems. Part III: Approximate solution of stochastic partial differential equations. J. of Hydrol. 82, 285–306
Unny, T.E. 1989: Stochastic partial differential equations in groundwater hydrology: Part 1—Theory. This journal (in press), to appear in Vol 3 (3)
Vecchia, A. V. 1985: A general class of models for stationary two-dimensional random processes. Biometrika 72, 281–291
Vecchia, A. V. 1988: Estimation and identification for continuous spatial processes. J. R. Statist. Soc. B 50, 297–312
Wahba, G. 1984: Surface fitting with scattered noisy data on Euclidean d-space and on the sphere. Rocky Mountain J. of Math. 14, 281–299
Whittle, P. 1954: On stationary processes in the plane. Biometrika 41, 434–449
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Jones, R.H. Fitting a stochastic partial differential equation to aquifer data. Stochastic Hydrol Hydraul 3, 85–96 (1989). https://doi.org/10.1007/BF01544074
Accepted:
Issue Date:
DOI: https://doi.org/10.1007/BF01544074