Spatial uncertainty analysis: Propagation of interpolation errors in spatially distributed models

被引:59
作者
Phillips, DL
Marks, DG
机构
[1] U.S. Environmental Protection Agency, USEPA Natl. Hlth./Environ. E., Corvallis, OR 97333
[2] U.S. Geological Survey, USEPA Natl. Hlth./Environ. E., Corvallis, OR 97333
关键词
error; evapotranspiration; kriging; spatial patterns; uncertainty analysis;
D O I
10.1016/0304-3800(95)00191-3
中图分类号
Q14 [生态学(生物生态学)];
学科分类号
071012 ; 0713 ;
摘要
In simulation modelling, it is desirable to quantify model uncertainties and provide not only point estimates for output variables but confidence intervals as well. Spatially distributed physical and ecological process models are becoming widely used, with runs being made over a grid of points that represent the landscape. This requires input values at each grid point, which often have to be interpolated from irregularly scattered measurement sites, e.g., weather stations. Interpolation introduces spatially varying errors which propagate through the model. We extended established uncertainty analysis methods to a spatial domain for quantifying spatial patterns of input variable interpolation errors and how they propagate through a model to affect the uncertainty of the model output. We applied this to a model of potential evapotranspiration (PET) as a demonstration. We modelled PET for three time periods in 1990 as a function of temperature, humidity, and wind on a 10-km grid across the U.S. portion of the Columbia River Basin. Temperature, humidity, and wind speed were interpolated using kriging from 700-1000 supporting data points. Kriging standard deviations (SD) were used to quantify the spatially varying interpolation uncertainties. For each of 5693 grid points, 100 Monte Carlo Simulations were done, using the kriged values of temperature, humidity, and wind, plus random error terms determined by the kriging SDs and the correlations of interpolation errors among the three variables. For the spring season example, kriging SDs averaged 2.6 degrees C for temperature, 8.7% for relative humidity, and 0.38 m s(-1) for wind. The resultant PET estimates had coefficients of variation (CVs) ranging from 14% to 27% for the 10-km grid cells. Maps of PET means and CVs showed the spatial patterns of PET with a measure of its uncertainty due to interpolation of the input variables. This methodology should be applicable to a variety of spatially distributed models using interpolated inputs.
引用
收藏
页码:213 / 229
页数:17
相关论文
共 26 条
[1]  
[Anonymous], 1987, Wind energy resource atlas of the United States
[2]  
Brutsaert W., 2013, Evaporation into the Atmosphere: Theory, History and Applications
[3]  
DEUTSCH CV, 1992, GSLIB GEOSTATISTICAL, P117
[4]  
GARDNER GH, 1983, NW83555 STUDSV EN AB
[5]  
GARDNER RH, 1990, PROCESS MODELING OF FOREST GROWTH RESPONSES TO ENVIRONMENTAL STRESS, P208
[6]   A DISTRIBUTION-FREE APPROACH TO INDUCING RANK CORRELATION AMONG INPUT VARIABLES [J].
IMAN, RL ;
CONOVER, WJ .
COMMUNICATIONS IN STATISTICS PART B-SIMULATION AND COMPUTATION, 1982, 11 (03) :311-334
[7]  
Isaaks EH., 1989, INTRO APPL GEOSTATIS, P278
[8]   WHEN DO WE NEED A TREND MODEL IN KRIGING [J].
JOURNEL, AG ;
ROSSI, ME .
MATHEMATICAL GEOLOGY, 1989, 21 (07) :715-739
[9]  
King AW., 1991, QUANTITATIVE METHODS, P479
[10]   CLIMATE AND ENERGY EXCHANGE AT THE SNOW SURFACE IN THE ALPINE REGION OF THE SIERRA-NEVADA .2. SNOW COVER ENERGY-BALANCE [J].
MARKS, D ;
DOZIER, J .
WATER RESOURCES RESEARCH, 1992, 28 (11) :3043-3054