A method for the interpolation of nonnegative functions with an application to contaminant load estimation |
| |
Authors: | A M Michalak P K Kitanidis |
| |
Institution: | (1) Department of Civil and Environmental Engineering, Stanford University, Stanford, CA 94305-4020, USA;(2) Climate Monitoring and Diagnostics Laboratory (CMDL), National Oceanic and Atmospheric Administration (NOAA), Mailcode R/CMDL1, 325 Broadway, Boulder, Colorado 80305-3328, USA |
| |
Abstract: | The objective of this work is to extend kriging, a geostatistical interpolation method, to honor parameter nonnegativity. The new method uses a prior probability distribution based on reflected Brownian motion that enforces this constraint. The work presented in this paper focuses on interpolation problems where the unknown is a function of a single variable (e.g. time), and is developed both for the case with and without measurement error in the available data. The algorithms presented for conditional simulations are computationally efficient, particularly in the case with no measurement error. We present an application to the interpolation of dissolved arsenic concentration data from the North Fork of the Humboldt River, Nevada. |
| |
Keywords: | Inference under constraints Nonnegativity Reflected Brownian motion Geostatistics Kriging Gibbs sampler |
本文献已被 SpringerLink 等数据库收录! |
|