Geostatistics in Hydrology: Kriging interpolation
|
|
- Corey Wilson
- 6 years ago
- Views:
Transcription
1 Chapter Geostatistics in Hydrology: Kriging interpolation Hydrologic properties, such as rainfall, aquifer characteristics (porosity, hydraulic conductivity, transmissivity, storage coefficient, etc.), effective recharge and so on, are all functions of space (and time) and often display a high spatial variability, also called heterogeneity. This variability is not in general random. It is a general rule that these properties display a so called scale effect, i.e., if we take measurements at two different points the difference in the measured values dicreases as the two points come closer to each other. It is convenient in certain cases to consider these properties as random functions having a given spatial structure, or in other word having a given spatial correlation, which can be conveniently described using appropriate statistical quantities. These variables are called regionalized variables []. The study of regionalized variables starts from the ability to interpolate a given field starting from a limited number of observation, but preserving the theoretical spatial correlation. This is accomplished by means of a technique called kriging developed by [, ] and largely applied in hydrology and other earth science disciplines [,, 7, ] for the spatial interpolation of various physical quantities given a number of spatially distributed measurments. Although theoreticallly kriging cannot be considered superior to other surface fitting techniques [7] and the use of a few arbitrary parameters may lead absurd results, this method is capable of obtaining objective interpolations evaluating at the same time the quality of the results. In the next sections we discuss first the statistical hypothesis that are needed to develop the theory of kriging and then we proceed at describing the method under different assumptions and finally report a few applications.. Statistical assumptions Let Z( x, ξ) be a random function (RF) simulating our hydrological quantity of interest. We recall that with this notation we imply that x denotes a point in space (two or three dimensional
2 CAPTER. GEOSTATISTICS IN HYDROLOGY space), while ξ denotes the state variable in the space of the realizations. In other words: Z( x, ξ) is a random variable at point x representing the entire set of realizations of the RF at point x ; Z( x, ξ ) is a particular realization of the RF Z; Z( x, ξ ) is a measurement point. Given a set of sampled values Z( x i, ξ ), i =,,... we want to reconstruct Z( x, ξ ), i.e. a possible realization of Z( x, ξ).. Kriging for Weak or Second Order stationary RF An RF is said to be second order stationary if:. Constant mean: E[Z( x, ξ)] = m (.). the autocovariance (another name for the covariance) is a function of the distance between the reference points x and x : cov[ x, x ] = E[(Z( x, ξ) m)(z( x, ξ) m)] = C( h) (.) In practice second order stationarity implies that the first two statistical moments (expected value and covariance) be translation invariant. Note that by saying that E[Z( x, ξ)] = m we imply that effectively the expected value taken on all the possible realizations ξ does not vary with space. However, for a given realization Z( x, ξ ) is a function of x. Because of (.), the covariance (.) can be written as: cov[ x, x ] = E[(Z( x, ξ) m)(z( x, ξ) m)] = E[Z( x, ξ)z( x, ξ)] me[z( x, ξ)] E[Z( x, ξ)]m + m = E[Z( x, ξ)z( x, ξ)] m (.) Obviously if h = we have the definition of variance, also called the dispersion variance: C() = var[z] = σ Z For simplicity of notation, from now on, we will denote x = x, h = h and we will drop the variable ξ in the RF.
3 . SECOND ORDER STATIONARY RF Case with m and C( h) known. new variable Y (x) with zero mean: If E[Z(x)] = m and C(h) are known, then we can define a Given the observed values: Y (x) = Z(x) m E[Y (x)] = x x... x n Y Y... Y n with Y i = Y (x i ) being the observation at point x i, we look for a linear estimator Y (x ) of Y (x ) at point x using the observed values. The form of the estimator is: Y (x ) = λ i Y i (.) i= Note that the estimator (.) is a realization of the RF Y (x, ξ) = λ i Y (x i, ξ) i= The weights λ i are calculated by imposing that the statistical error ɛ(x ) = Y (x ) Y (x ) has zero expecetd value and minimum variance: E[ɛ(x )] = var[ɛ(x )] = E[(Y (x ) Y (x )) ] = minimum Substituting eq. (.) in the expression of the variance we have: E[(Y (x ) Y (x )) ] = E[( λ i Y i Y ) ] i= = E[( λ i Y i Y )( λ i Y i Y )] i= i= i= i= = E[( λ i Y i )( λ i Y i )] E[ λ i Y i Y ] + E[Y ] = i= λ i λ j E[Y i Y j ] j= i= i= λ i E[Y i Y ] + E[Y ] but, since m = E[Y i Y j ] = C(x i x j ) + m = C(x i x j )
4 CAPTER. GEOSTATISTICS IN HYDROLOGY and is the dispersion variance of Y. Then: E[Y ] = C() = var[y ] E[(Y (x ) Y (x )) ] = λ i λ j C(x i x j ) λ i C(x i x ) i= j= i= + C() The minimum is found by setting to zero the first partial derivatives: λ i ( E[(Y (x ) Y (x )) ] ) = λ j C(x i x j ) C(x i x ) = j= j =,..., n This yields a linear system of equations: where matrix C is given by: C = Cλ = b (.) C() C(x x )... C(x x n ) C(x n x ) and the right hand side vector b is given by: C(x x ) b = C(x n x ) Matrix C is the spatial covariance matrix and does not depend upon x. It can be shown that if all the x j s are distinct then C is positive definite, and thus the linear system (.) can be solved with either direct or iterative methods. Once the solution vector λ is obtained, equation (.) yields the estimation of our regionalized variable at point x. Thus the calculated value for λ is actually function of the estimation point x. If we want to change the estimation point x, for example if we need to obtain a spatial distribution of our regionalized variable, we need to solve the linear system (.) for different values of x. In this case it is convenient to factorize matrix C using Cholesky decomposition and then proceed to the solution for the different right hand side vectors. C()
5 .. KRIGING WITH THE INTRINSIC HYPOTHESIS Evaluation of the estimation variance The estimation variance is defined as the variance of the error: ɛ = Y Y Hence: var[y Y ] = E[(Y Y ) ] E[(Y Y )] but: E[(Y Y )] = E[Y ] E[Y ] = λ i E[Y i ] E[Y ] = (E[Y ] = m = ) i and since we finally obtain: λ j C(x i x j ) = C(x i x ) j var[y Y ] = E[(Y Y ) ] = λ i λ j C(x i x j ) i j i = i λ i C(x i x ) + C() λ i C(x i x ) + C() = var[y ] i λ i C(x i x ) which, since i λ ic(x i x ) >, shows that the estimation variance of Y is smaller than the dispersion variance of Y (the real variance of the RF). In statistical terms, we can interpret this result by saying that since we have observed Y at some points x i then the uncertainty on Y decreases. It is important to remark the difference between estimation and dispersion variance. The latter is representative of the variation interval of the RF Y within the interpolation domain, while the estimation variance represents the residual uncertainty in the estimation of the realization Y of Z when n observations are available. The dispersion variance is a constant, while the estimation variance varies from point to point and is zero at the observation points. Our original variable was Z = Y + m and its estimate is thus: Z = m + i λ i (Z i m) var[z Z ] = var[z ] i λ i C(x i x ). Kriging with the intrinsic hypothesis The hypothesis of second order stationarity of the RF is not always satisfied, for example C() increases with the distance, violating hypothesis (.). In this case the Intrinsic hypothesis must be used, in which we assume that the first order increments δ = Y (x + h) Y (x)
6 CAPTER. GEOSTATISTICS IN HYDROLOGY range High spatial correlation Low spatial correlation range C(h) γ(h) sill h h Figure.: Behavior of the covariance as a function of distance (left) and the corresponding variogram (right). are second order RF: E[Y (x + h) Y (x)] = m(h) = (.) var[y (x + h) Y (x)] = γ(h) (.7) where the function γ(h) is called the variogram. If the mean m(h) is not zero an obvious change of variable is required. The variogram is defined as the mean quadratic increment of Y (x) (divided by ) for any two points x i and x j separated by a distance h: γ(h) = var[y (x + h) Y (x)] = E[(Y (x + h) Y (x)) ] (.8) and is related to the covariance function by: γ(h) = E[(Y (x + h) Y (x)) ] = E[Y (x + h)] E[Y (x + h)y (x)] + E[Y (x)] = C() C(h) The intrinsic hypothesis requires a finite value for the mean of Y (x) but not for its variance. In fact, hypothesis (.), as changed into (.), implies (.8), but not viceversa. The covariance C(h) has a decreasing behavior as shown in Fig... When C(h) is known then the variogram can be directly calculated. When C() is finite, the variogram γ(h) is bounded asymptotically by this value. The value of h at which the asymptot can be considered achieved is called the range, while C() is called the sill (see Fig..).
7 .. KRIGING WITH THE INTRINSIC HYPOTHESIS 7.. The variogram The variogram is usually calculated from the experimental observations, and describes the spatial structure of the RF. It can be shown [] that given x,..., x n are n points belonging to the domain of interpolation, for all the coefficients η,..., η n R satisfing: then i= η i = i= η i η j γ(x i x j ) j= i.e. γ(h) is a positive semidefinite function, and γ(h) lim = h h i.e. γ(h) is tends to infinity slowlier than h as h. In principle there γ(h) could assume different behaviors also with the direction of vector h ( anisotropy ), but this is in general not easily verifiable due to the limited number of data points usually available for hydrologic variables. If the experimental variogram displays anisotropy, then the intrinsic hypothesis is not verified and one has to use the so called universal kriging [, ]. The most commonly used isotropic variograms are shown in Fig.. and are of the form:. polinomial variogram:. exponential variogram:. gaussian variogram: γ(h) = ωh α < α < γ(h) = ω [ e αh] γ(h) = ω [ e (αh)]. spherical variogram: γ(h) = { [ ω h ( ) ] h α α ω h α h > α where ω and α are real constants. The variogram is estimated from the available observations in the following manner. The data points are subdivided into a prefixed number of classes based on the distances between the measurement locations. For each pair i and j of points and for each class calculate:. the number M of points that fall within the class;
8 8 CAPTER. GEOSTATISTICS IN HYDROLOGY 8 α=. α= α=.. α α α....8 α/ α Figure.: Behavior of the most commmonly used variograms: polinomial (top left); exponential (top right); gaussian (bottom left); spherical (bottom right)
9 .. KRIGING WITH THE INTRINSIC HYPOTHESIS 9. the average distance of the class;. the half of the mean quadratic increment (Yi Y j ) /M In general the pairs are not uniformly distributed among the different classes as usually there are more pairs for the smaller distances. Thus the experimental variogram will be less meaningful as h increases. A best fit procedure together with visual inspection is then used to select the most appropriate variogram and evaluate its optimal parameters. Remark. An experimental variogram that is not bounded above (for example the polinomial variogram with α ) implies an infinite variance, and thus the covariance does not exist. Only the intrinsic hypothesis is acceptable. If the variogram achieves a sill then the phenomenon has a finite variance and the covariance exists. Remark. For every variogram, γ() =. However sometimes the data may display a jump at the origin. This apparent discontinuity is called the nugget effect and can be due to measurement errors or to microregionalization effects that are not evidenced at the scale of the actual data points. If the nugget effect is present the variogram needs to be changed into: γ(h) = δ + γ (h) where δ is the jump at the origin and γ (h) the variogram without the jump. In Fig.. we report a variogram with the nugget effect together with some other special cases that may be encountered. Remark. When we look for an estimator of the type Y (x ) = λ i Y i i= from the intrinsic hypothesis, i.e. eq. (.), we see that our random variable Y (x) must have a constant mean m, so that: E[Y (x )] = E[ λ i Y i ] = m i= and, since E[Y (x i )] = E[Y i ] = m, we readily obtain: λ i = i=
10 CAPTER. GEOSTATISTICS IN HYDROLOGY γ(h) γ(h) δ δ measured values nugget effect random structure (no correlation) h h γ(h) γ(h) nested structure periodic structure h h γ(h) hole effect Figure.: Example of typical spatial correletion structures that may be encountered in analyzing measured data h
11 .. REMARKS ABOUT KRIGING The condition of minimum variance becomes now a constrained minimization problem which can be solve by introducing a lagrange multiplier µ. The system (.) must be rewritten as: γ(x x )... γ(x x n ) λ γ(x x )... γ(x n x ) λ n =. γ(x n x )... µ. Remarks about Kriging a. Kriging is a BLUE (Best Linear Unbiased Estimator) interpolator. In other word it is a Linear estimator that matches the correct expected value of the population (Unbiased) and that minimizes the variance of the observations (Best). b. Kriging is an exact interpolator if no errors are present. In fact, if we set x = x i in (.) we obtain immediately λ i =, λ j =, j =,..., n, j i. c. If we assume that the the error ɛ is Gaussian, then we can associate to the estimate Y (x ) a confidence interval. For example, the 9% confidence interval is ±σ where: Then the kriging estimator (.) becomes: σ = var[y (x ) Y (X )] Y (x ) = λ i Y i ± σ i= d. The solution of the linear system does not depend on the observed value but only on x i and x. e. A map of the estimated regionalized variable, and possibly its confidence intervals, can be obtained be defining a grid and solving the linear system for each point in the grid.. Kriging with uncertainties We now assume that the observations Y i are affected by measurement errors ɛ i, and that:. the errors ɛ i have zero mean: E[ɛ i ] = i =,..., n. the errors are uncorrelated: cov[ɛ i, ɛ j ] = i j
12 CAPTER. GEOSTATISTICS IN HYDROLOGY. the errors are not correlated with the RF: cov[ɛ i, Y i ] =. the variance σ i of the errors is a known quantity and can vary from point to point. The new coefficient matrix C of the linear system (.) is changed by adding to the main diagonal the quantity σi : σ C := C σn and everything proceeds as in the standard case.. Validation of the interpolation model The chosen model (in practice the variogram) can be validated by interpolating observed values. If n observations Y (x i ), i =,..., n are available, the validation process proceeds as follows: For each j, j =,..., n: discard point (x j, Y (x j )); estimate the Y (x j ) by solving the kriging system having set x = x j and using the remaining points x i, i j for the interpolation; evaluate the estimation error ɛ j = Y j Y j, The model can be considered theoretically valid if the error distribution is approximately gaussian with zero mean and unit variance (N(, )), i.e. satisfies the following:. there is no bias: n ɛ i i=. the estimation variance σ i is coherent with the error standard deviation: ( ) Y i Y i n i= One can also look at the behavior of the interpolation error at each point looking at the mean square error of the vector ɛ: Q = ɛ i n σ i i=
13 .7. COMPUTATIONAL ASPECTS The uncertainties connected to the choice of the theoretical variogram from the experimental data can be minimized by anaylizing the validation test. In fact, among all the possible variograms γ(h), that close to the origin display a slope compatible with the observations and gives rise to a theoretically coherent model, one can choose the variogram with the smallest value of Q..7 Computational aspects In the validation phase n linear systems of dimension n need to be solved. The system matrices are obtained by dropping one row and one column of the complete kriging matrix. This can be efficiently accomplished by means of intersections of n -dimensional lines with appropriate coordinate n-dimensional planes. Note that the kriging matrix C is symmetric, and thus its eigenvalues µ i are real. However, since µ i = Tr(C) = c ii = i= where Tr(C) is the trace of matrix C, it follows that some of the eigenvalues must be negative and thus C is not positive definite. For this reason, the solution of the linear systems is usually obtained by means of direct methods, such as Gaussian elimination or Choleski decomposition. Full Pivoting is often necessary to maintain stability of the algorithm..8 Kriging with moving neighborhoods Generally the experimental variogram is most accurate for small values of h, with uncertainties growing rapidly when h is large. The influence of this problem may be decreased by using moving neighborhoods. With this variant, only the points that lie within a prefixed radius R from point x are considered, provided that we are left with an adequate number of data (Fig..). The radius R is selected so that the lag h will remain within the range of maximum certainty for γ(h). This approach leads also to high saving in the computational cost of the procedure because each linear system is now much smaller (Gaussian elimination has a computational cost proportional to n, where n is the number of equations)..9 Detrending In certain cases the observations display a definite trend that needs to be taken into account. This is the case, for example, when one needs to interpolate piezometric heads observed from wells in an aquifer system where a regional gradient is present. To remove the trend from the data it is possible to work with residuals (= measurements - trend) that have a constant mean. However this detrending procedure may be dangerous as it may introduce a bias in the results. For this reason it is important that the trend be recognized not only from the raw data but i=
14 CAPTER. GEOSTATISTICS IN HYDROLOGY x R Figure.: Example of interpolation with moving neighborhood also from the physical behavior of the system from which the data come. If the trend cannot be removed from the observations by simple subtraction, the universal kriging approach [, ] can be used.
15 .. EXAMPLE OF APPLICATION OF KRIGING FOR THE RECONSTRUCTION OF AN ELEVATION. Example of application of kriging for the reconstruction of an elevation map ([DPSOHRI. U L J L Q J The range of the observations is - 8
16 Variogramma (m) Experimental variogram using classes n. coppie Lag Distance (km) Experimental variogram using classes Lag Distance (km) The shape is similar: robustness with respect to different class subdivisions Variogramma (m) CAPTER. GEOSTATISTICS IN HYDROLOGY
17 7 First try Linear variogram γ = h elevation Kriging reconstruction is robust evenwhen the variogram is not accurate Estimation variance meaningless interpolation error ε EXAMPLE OF APPLICATION 7
18 Second try Linear variogram with nugget effect Variogram (m) nugget effect = Lag Distance (km) 8 CAPTER. GEOSTATISTICS IN HYDROLOGY
19 Elevation Reconstructed values are smoothed because of nugget effect: large estimation variance Estimation variance ε EXAMPLE OF APPLICATION 9
20 Third try Variogram without nugget effect and with large slope at the origin Variogram (m) polinomial variogram: ω = = Lag Distance (km) CAPTER. GEOSTATISTICS IN HYDROLOGY
21 Elevation Estimation variance is large far away from observation points Estimation variance ε E E- -9E- -E E- -E-. E E E-7 E EXAMPLE OF APPLICATION
22 Fourth try Gaussian variogram: small slope at the origin Variogram (m) ω = 7 = Lag Distance (km) CAPTER. GEOSTATISTICS IN HYDROLOGY
23 Elevation small and smooth estimation variance; small errors Estimation variance ε 7 -E E E E E- 7E- -E E- 8E EXAMPLE OF APPLICATION
24 CAPTER. GEOSTATISTICS IN HYDROLOGY
25 Bibliography [] de Marsily, G. (98). Spatial variability of properties in porous media: A stochastic approach. In J. Bear and M. Corapcioglu, editors, Fundemantals of Transport Phenomena in Porous Media, pages 7 79, Dordrecht. Martinus Nijhoff Publ. [] Delhomme, J. P. (97). Applications de la thèorie des varaibles règionalisèes dans les sciences de l eau. Ph.D. thesis, Ecoles des Mines, Fontainebleaus (France). [] Delhomme, J. P. (978). Kriging in the hydrosciences. Adv. Water Resources,,. [] Matheron, G. (99). Le krigeage universel. Technical Report, Paris School of Mines. Cah. Cent. Morphol. Math., Fontainbleau. [] Matheron, G. (97). The theory of regionalized variables and its applications. Technical Report, Paris School of Mines. Cah. Cent. Morphol. Math., Fontainbleau. [] Volpi, G. and Gambolati, G. (978). On the use of a main trend for the kriging technique in hydrology. Adv. Water Resources,, 9. [7] Volpi, G., Gambolati, G., Carbognin, L., Gatto, P., and Mozzi, G. (979). Groundwater contour mapping in venice by stochastic interpolators.. results. Water Resour. Res.,, 9 97.
Basics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods. Hans Wackernagel. MINES ParisTech.
Basics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods Hans Wackernagel MINES ParisTech NERSC April 2013 http://hans.wackernagel.free.fr Basic concepts Geostatistics Hans Wackernagel
More information7 Geostatistics. Figure 7.1 Focus of geostatistics
7 Geostatistics 7.1 Introduction Geostatistics is the part of statistics that is concerned with geo-referenced data, i.e. data that are linked to spatial coordinates. To describe the spatial variation
More informationLecture 9: Introduction to Kriging
Lecture 9: Introduction to Kriging Math 586 Beginning remarks Kriging is a commonly used method of interpolation (prediction) for spatial data. The data are a set of observations of some variable(s) of
More informationPRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH
PRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH SURESH TRIPATHI Geostatistical Society of India Assumptions and Geostatistical Variogram
More informationPoint-Referenced Data Models
Point-Referenced Data Models Jamie Monogan University of Georgia Spring 2013 Jamie Monogan (UGA) Point-Referenced Data Models Spring 2013 1 / 19 Objectives By the end of these meetings, participants should
More informationBasics of Point-Referenced Data Models
Basics of Point-Referenced Data Models Basic tool is a spatial process, {Y (s), s D}, where D R r Chapter 2: Basics of Point-Referenced Data Models p. 1/45 Basics of Point-Referenced Data Models Basic
More informationInfluence of parameter estimation uncertainty in Kriging: Part 2 Test and case study applications
Hydrology and Earth System Influence Sciences, of 5(), parameter 5 3 estimation (1) uncertainty EGS in Kriging: Part Test and case study applications Influence of parameter estimation uncertainty in Kriging:
More informationGeostatistics for Gaussian processes
Introduction Geostatistical Model Covariance structure Cokriging Conclusion Geostatistics for Gaussian processes Hans Wackernagel Geostatistics group MINES ParisTech http://hans.wackernagel.free.fr Kernels
More informationHandbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp
Handbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp Marcela Alfaro Córdoba August 25, 2016 NCSU Department of Statistics Continuous Parameter
More informationEstimation of the transmissmty of the Santiago aquifer, Chile, using different geostatistical methods
Groundwater Management: Quantity and Quality (Proceedings of the Benidorm Symposium, October 1989). IAHS Publ. no. 188,1989. Estimation of the transmissmty of the Santiago aquifer, Chile, using different
More informationAssessment of Hydraulic Conductivity Upscaling Techniques and. Associated Uncertainty
CMWRXVI Assessment of Hydraulic Conductivity Upscaling Techniques and Associated Uncertainty FARAG BOTROS,, 4, AHMED HASSAN 3, 4, AND GREG POHLL Division of Hydrologic Sciences, University of Nevada, Reno
More informationIs there still room for new developments in geostatistics?
Is there still room for new developments in geostatistics? Jean-Paul Chilès MINES ParisTech, Fontainebleau, France, IAMG 34th IGC, Brisbane, 8 August 2012 Matheron: books and monographs 1962-1963: Treatise
More informationHighly Robust Variogram Estimation 1. Marc G. Genton 2
Mathematical Geology, Vol. 30, No. 2, 1998 Highly Robust Variogram Estimation 1 Marc G. Genton 2 The classical variogram estimator proposed by Matheron is not robust against outliers in the data, nor is
More informationChapter 4 - Fundamentals of spatial processes Lecture notes
TK4150 - Intro 1 Chapter 4 - Fundamentals of spatial processes Lecture notes Odd Kolbjørnsen and Geir Storvik January 30, 2017 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites
More informationIndex. Geostatistics for Environmental Scientists, 2nd Edition R. Webster and M. A. Oliver 2007 John Wiley & Sons, Ltd. ISBN:
Index Akaike information criterion (AIC) 105, 290 analysis of variance 35, 44, 127 132 angular transformation 22 anisotropy 59, 99 affine or geometric 59, 100 101 anisotropy ratio 101 exploring and displaying
More informationKriging models with Gaussian processes - covariance function estimation and impact of spatial sampling
Kriging models with Gaussian processes - covariance function estimation and impact of spatial sampling François Bachoc former PhD advisor: Josselin Garnier former CEA advisor: Jean-Marc Martinez Department
More informationA STUDY OF MEAN AREAL PRECIPITATION AND SPATIAL STRUCTURE OF RAINFALL DISTRIBUTION IN THE TSEN-WEN RIVER BASIN
Journal of the Chinese Institute of Engineers, Vol. 24, No. 5, pp. 649-658 (2001) 649 A STUDY OF MEAN AREAL PRECIPITATION AND SPATIAL STRUCTURE OF RAINFALL DISTRIBUTION IN THE TSEN-WEN RIVER BASIN Jet-Chau
More informationSpatial Backfitting of Roller Measurement Values from a Florida Test Bed
Spatial Backfitting of Roller Measurement Values from a Florida Test Bed Daniel K. Heersink 1, Reinhard Furrer 1, and Mike A. Mooney 2 1 Institute of Mathematics, University of Zurich, CH-8057 Zurich 2
More informationIntroduction. Semivariogram Cloud
Introduction Data: set of n attribute measurements {z(s i ), i = 1,, n}, available at n sample locations {s i, i = 1,, n} Objectives: Slide 1 quantify spatial auto-correlation, or attribute dissimilarity
More informationChapter XXII The Covariance
Chapter XXII The Covariance Geometrical Covariogram Random Version : -Sets -Functions - Properties of the covariance Intrinsic Theory : - Variogram J. Serra Ecole des Mines de Paris ( 2000 ) Course on
More informationDrill-Holes and Blast-Holes
Drill-Holes and Blast-Holes Serge Antoine Séguret 1 and Sebastian De La Fuente 2 Abstract The following is a geostatistical study of copper measurements on samples from diamond drill-holes and blast-holes.
More informationStatistícal Methods for Spatial Data Analysis
Texts in Statistícal Science Statistícal Methods for Spatial Data Analysis V- Oliver Schabenberger Carol A. Gotway PCT CHAPMAN & K Contents Preface xv 1 Introduction 1 1.1 The Need for Spatial Analysis
More informationCovariance function estimation in Gaussian process regression
Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian
More information11/8/2018. Spatial Interpolation & Geostatistics. Kriging Step 1
(Z i Z j ) 2 / 2 (Z i Zj) 2 / 2 Semivariance y 11/8/2018 Spatial Interpolation & Geostatistics Kriging Step 1 Describe spatial variation with Semivariogram Lag Distance between pairs of points Lag Mean
More informationMultivariate Geostatistics
Hans Wackernagel Multivariate Geostatistics An Introduction with Applications Third, completely revised edition with 117 Figures and 7 Tables Springer Contents 1 Introduction A From Statistics to Geostatistics
More informationCoregionalization by Linear Combination of Nonorthogonal Components 1
Mathematical Geology, Vol 34, No 4, May 2002 ( C 2002) Coregionalization by Linear Combination of Nonorthogonal Components 1 J A Vargas-Guzmán, 2,3 A W Warrick, 3 and D E Myers 4 This paper applies the
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 y Kriging Step 1 Describe spatial variation with Semivariogram (Z i Z j ) 2 / 2 Point cloud Map 3
More informationDeterministic solution of stochastic groundwater flow equations by nonlocalfiniteelements A. Guadagnini^ & S.P. Neuman^
Deterministic solution of stochastic groundwater flow equations by nonlocalfiniteelements A. Guadagnini^ & S.P. Neuman^ Milano, Italy; ^Department of. Hydrology & Water Resources, The University ofarizona,
More informationThe presentation is entitled ANISOTROPY OF THE ROCK QUALITY DESIGNATION & ITS GEOSTATISTICAL EVALUATION».
It is my pleasure to present the work done jointly by my collaborator Doctor Serge Séguret from France, and myself, on data provided by my team and by my colleague Claudio Rojas The presentation is entitled
More informationKriging Luc Anselin, All Rights Reserved
Kriging Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline Principles Kriging Models Spatial Interpolation
More informationInvestigation of Monthly Pan Evaporation in Turkey with Geostatistical Technique
Investigation of Monthly Pan Evaporation in Turkey with Geostatistical Technique Hatice Çitakoğlu 1, Murat Çobaner 1, Tefaruk Haktanir 1, 1 Department of Civil Engineering, Erciyes University, Kayseri,
More informationL26: Advanced dimensionality reduction
L26: Advanced dimensionality reduction The snapshot CA approach Oriented rincipal Components Analysis Non-linear dimensionality reduction (manifold learning) ISOMA Locally Linear Embedding CSCE 666 attern
More informationAcceptable Ergodic Fluctuations and Simulation of Skewed Distributions
Acceptable Ergodic Fluctuations and Simulation of Skewed Distributions Oy Leuangthong, Jason McLennan and Clayton V. Deutsch Centre for Computational Geostatistics Department of Civil & Environmental Engineering
More informationGeostatistics: Kriging
Geostatistics: Kriging 8.10.2015 Konetekniikka 1, Otakaari 4, 150 10-12 Rangsima Sunila, D.Sc. Background What is Geostatitics Concepts Variogram: experimental, theoretical Anisotropy, Isotropy Lag, Sill,
More informationBuilding Blocks for Direct Sequential Simulation on Unstructured Grids
Building Blocks for Direct Sequential Simulation on Unstructured Grids Abstract M. J. Pyrcz (mpyrcz@ualberta.ca) and C. V. Deutsch (cdeutsch@ualberta.ca) University of Alberta, Edmonton, Alberta, CANADA
More informationENGRG Introduction to GIS
ENGRG 59910 Introduction to GIS Michael Piasecki October 13, 2017 Lecture 06: Spatial Analysis Outline Today Concepts What is spatial interpolation Why is necessary Sample of interpolation (size and pattern)
More informationIntroduction. Spatial Processes & Spatial Patterns
Introduction Spatial data: set of geo-referenced attribute measurements: each measurement is associated with a location (point) or an entity (area/region/object) in geographical (or other) space; the domain
More informationSpatial Data Mining. Regression and Classification Techniques
Spatial Data Mining Regression and Classification Techniques 1 Spatial Regression and Classisfication Discrete class labels (left) vs. continues quantities (right) measured at locations (2D for geographic
More information7. MULTIVARATE STATIONARY PROCESSES
7. MULTIVARATE STATIONARY PROCESSES 1 1 Some Preliminary Definitions and Concepts Random Vector: A vector X = (X 1,..., X n ) whose components are scalar-valued random variables on the same probability
More informationA hybrid Marquardt-Simulated Annealing method for solving the groundwater inverse problem
Calibration and Reliability in Groundwater Modelling (Proceedings of the ModelCARE 99 Conference held at Zurich, Switzerland, September 1999). IAHS Publ. no. 265, 2000. 157 A hybrid Marquardt-Simulated
More informationBusiness Economics BUSINESS ECONOMICS. PAPER No. : 8, FUNDAMENTALS OF ECONOMETRICS MODULE No. : 3, GAUSS MARKOV THEOREM
Subject Business Economics Paper No and Title Module No and Title Module Tag 8, Fundamentals of Econometrics 3, The gauss Markov theorem BSE_P8_M3 1 TABLE OF CONTENTS 1. INTRODUCTION 2. ASSUMPTIONS OF
More informationSpatial Statistics with Image Analysis. Lecture L02. Computer exercise 0 Daily Temperature. Lecture 2. Johan Lindström.
C Stochastic fields Covariance Spatial Statistics with Image Analysis Lecture 2 Johan Lindström November 4, 26 Lecture L2 Johan Lindström - johanl@maths.lth.se FMSN2/MASM2 L /2 C Stochastic fields Covariance
More informationProblems in Space-Time Kriging of Geohydrological Data 1. Shahrokh Rouhani 2 and Donald E. Myers 3
Mathematical Geology, Vol. 22, No. 5, 1990 Problems in Space-Time Kriging of Geohydrological Data 1 Shahrokh Rouhani 2 and Donald E. Myers 3 Spatiotemporal variables constitute a large class of geohydrological
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Department of Forestry & Department of Geography, Michigan State University, Lansing Michigan, U.S.A. 2 Biostatistics,
More informationAdaptive Sampling of Clouds with a Fleet of UAVs: Improving Gaussian Process Regression by Including Prior Knowledge
Master s Thesis Presentation Adaptive Sampling of Clouds with a Fleet of UAVs: Improving Gaussian Process Regression by Including Prior Knowledge Diego Selle (RIS @ LAAS-CNRS, RT-TUM) Master s Thesis Presentation
More informationOptimization of Sampling Locations for Variogram Calculations
WATER RESOURCES RESEARCH, VOL. 23, N 3, PAGES 496-500, MARCH 1987 Optimization of Sampling Locations for Variogram Calculations A. W. WARRICK Department of Soil and Water Science, The University of Arizona,
More informationSimple Linear Regression
Simple Linear Regression ST 430/514 Recall: A regression model describes how a dependent variable (or response) Y is affected, on average, by one or more independent variables (or factors, or covariates)
More informationMultivariate Time Series: VAR(p) Processes and Models
Multivariate Time Series: VAR(p) Processes and Models A VAR(p) model, for p > 0 is X t = φ 0 + Φ 1 X t 1 + + Φ p X t p + A t, where X t, φ 0, and X t i are k-vectors, Φ 1,..., Φ p are k k matrices, with
More informationInterpolation and 3D Visualization of Geodata
Marek KULCZYCKI and Marcin LIGAS, Poland Key words: interpolation, kriging, real estate market analysis, property price index ABSRAC Data derived from property markets have spatial character, no doubt
More informationStatistical signal processing
Statistical signal processing Short overview of the fundamentals Outline Random variables Random processes Stationarity Ergodicity Spectral analysis Random variable and processes Intuition: A random variable
More informationOn dealing with spatially correlated residuals in remote sensing and GIS
On dealing with spatially correlated residuals in remote sensing and GIS Nicholas A. S. Hamm 1, Peter M. Atkinson and Edward J. Milton 3 School of Geography University of Southampton Southampton SO17 3AT
More informationBayesian Transgaussian Kriging
1 Bayesian Transgaussian Kriging Hannes Müller Institut für Statistik University of Klagenfurt 9020 Austria Keywords: Kriging, Bayesian statistics AMS: 62H11,60G90 Abstract In geostatistics a widely used
More informationSampling of variables with proportional effect
Sampling of variables with proportional effect Samuel Canchaya, R&D Senior Geologist, Cía. de Minas Buenaventura SAA, Lima, Peru, samuel.canchaya@buenaventura.pe ABSTRACT According to the theory of regionalized
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry
More informationDerivatives of Spatial Variances of Growing Windows and the Variogram 1
Mathematical Geology, Vol. 3, No. 7, 000 Derivatives of Spatial Variances of Growing Windows and the Variogram 1 J. A. Vargas-Guzmán, D. E. Myers, 3 and A. W. Warrick 4 Geostatistical analysis of spatial
More informationSimple closed form formulas for predicting groundwater flow model uncertainty in complex, heterogeneous trending media
WATER RESOURCES RESEARCH, VOL. 4,, doi:0.029/2005wr00443, 2005 Simple closed form formulas for predicting groundwater flow model uncertainty in complex, heterogeneous trending media Chuen-Fa Ni and Shu-Guang
More informationMUCK - A ovel Approach to Co-Kriging
Clark, W.V. Harper and K.L. Basinger, MUCK -- A novel approach to co-kriging in Proc. DOE/AECL 87 Conference on Geostatistical, Sensitivity and Uncertainty Methods for Ground-Water Flow and Radionuclide
More informationData Mining. Dimensionality reduction. Hamid Beigy. Sharif University of Technology. Fall 1395
Data Mining Dimensionality reduction Hamid Beigy Sharif University of Technology Fall 1395 Hamid Beigy (Sharif University of Technology) Data Mining Fall 1395 1 / 42 Outline 1 Introduction 2 Feature selection
More informationSpatial and Environmental Statistics
Spatial and Environmental Statistics Dale Zimmerman Department of Statistics and Actuarial Science University of Iowa January 17, 2019 Dale Zimmerman (UIOWA) Spatial and Environmental Statistics January
More informationComparison of rainfall distribution method
Team 6 Comparison of rainfall distribution method In this section different methods of rainfall distribution are compared. METEO-France is the French meteorological agency, a public administrative institution
More informationAdvances in Locally Varying Anisotropy With MDS
Paper 102, CCG Annual Report 11, 2009 ( 2009) Advances in Locally Varying Anisotropy With MDS J.B. Boisvert and C. V. Deutsch Often, geology displays non-linear features such as veins, channels or folds/faults
More informationChapter 2. Some basic tools. 2.1 Time series: Theory Stochastic processes
Chapter 2 Some basic tools 2.1 Time series: Theory 2.1.1 Stochastic processes A stochastic process is a sequence of random variables..., x 0, x 1, x 2,.... In this class, the subscript always means time.
More informationSPATIAL ELECTRICAL LOADS MODELING USING THE GEOSTATISTICAL METHODS
19 th International CODATA Conference THE INFORMATION SOCIETY: NEW HORIZONS FOR SCIENCE Berlin, Germany 7-1 November 24 SPATIAL ELECTRICAL LOADS MODELING USING THE GEOSTATISTICAL METHODS Barbara Namysłowska-Wilczyńska
More informationPrincipal Component Analysis
Principal Component Analysis Laurenz Wiskott Institute for Theoretical Biology Humboldt-University Berlin Invalidenstraße 43 D-10115 Berlin, Germany 11 March 2004 1 Intuition Problem Statement Experimental
More informationDefect Detection using Nonparametric Regression
Defect Detection using Nonparametric Regression Siana Halim Industrial Engineering Department-Petra Christian University Siwalankerto 121-131 Surabaya- Indonesia halim@petra.ac.id Abstract: To compare
More informationBME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAW
7 VIII. BME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAW A wide variety of natural processes are described using physical laws. A physical law may be expressed by means of an
More informationReliability Theory of Dynamically Loaded Structures (cont.)
Outline of Reliability Theory of Dynamically Loaded Structures (cont.) Probability Density Function of Local Maxima in a Stationary Gaussian Process. Distribution of Extreme Values. Monte Carlo Simulation
More informationLECTURE 10 LINEAR PROCESSES II: SPECTRAL DENSITY, LAG OPERATOR, ARMA. In this lecture, we continue to discuss covariance stationary processes.
MAY, 0 LECTURE 0 LINEAR PROCESSES II: SPECTRAL DENSITY, LAG OPERATOR, ARMA In this lecture, we continue to discuss covariance stationary processes. Spectral density Gourieroux and Monfort 990), Ch. 5;
More informationREML Estimation and Linear Mixed Models 4. Geostatistics and linear mixed models for spatial data
REML Estimation and Linear Mixed Models 4. Geostatistics and linear mixed models for spatial data Sue Welham Rothamsted Research Harpenden UK AL5 2JQ December 1, 2008 1 We will start by reviewing the principles
More informationTime-lapse filtering and improved repeatability with automatic factorial co-kriging. Thierry Coléou CGG Reservoir Services Massy
Time-lapse filtering and improved repeatability with automatic factorial co-kriging. Thierry Coléou CGG Reservoir Services Massy 1 Outline Introduction Variogram and Autocorrelation Factorial Kriging Factorial
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationLimit Kriging. Abstract
Limit Kriging V. Roshan Joseph School of Industrial and Systems Engineering Georgia Institute of Technology Atlanta, GA 30332-0205, USA roshan@isye.gatech.edu Abstract A new kriging predictor is proposed.
More informationDimension Reduction Techniques. Presented by Jie (Jerry) Yu
Dimension Reduction Techniques Presented by Jie (Jerry) Yu Outline Problem Modeling Review of PCA and MDS Isomap Local Linear Embedding (LLE) Charting Background Advances in data collection and storage
More informationEnsemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations
Ensemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations Alberto Guadagnini (1,), Marco Panzeri (1), Monica Riva (1,), Shlomo P. Neuman () (1) Department of
More informationSIO 221B, Rudnick adapted from Davis 1. 1 x lim. N x 2 n = 1 N. { x} 1 N. N x = 1 N. N x = 1 ( N N x ) x = 0 (3) = 1 x N 2
SIO B, Rudnick adapted from Davis VII. Sampling errors We do not have access to the true statistics, so we must compute sample statistics. By this we mean that the number of realizations we average over
More informationLecture 5 Geostatistics
Lecture 5 Geostatistics Lecture Outline Spatial Estimation Spatial Interpolation Spatial Prediction Sampling Spatial Interpolation Methods Spatial Prediction Methods Interpolating Raster Surfaces with
More information4th HR-HU and 15th HU geomathematical congress Geomathematics as Geoscience Reliability enhancement of groundwater estimations
Reliability enhancement of groundwater estimations Zoltán Zsolt Fehér 1,2, János Rakonczai 1, 1 Institute of Geoscience, University of Szeged, H-6722 Szeged, Hungary, 2 e-mail: zzfeher@geo.u-szeged.hu
More informationConditional Distribution Fitting of High Dimensional Stationary Data
Conditional Distribution Fitting of High Dimensional Stationary Data Miguel Cuba and Oy Leuangthong The second order stationary assumption implies the spatial variability defined by the variogram is constant
More informationA robust statistically based approach to estimating the probability of contamination occurring between sampling locations
A robust statistically based approach to estimating the probability of contamination occurring between sampling locations Peter Beck Principal Environmental Scientist Image placeholder Image placeholder
More informationAnomaly Density Estimation from Strip Transect Data: Pueblo of Isleta Example
Anomaly Density Estimation from Strip Transect Data: Pueblo of Isleta Example Sean A. McKenna, Sandia National Laboratories Brent Pulsipher, Pacific Northwest National Laboratory May 5 Distribution Statement
More informationWhat s for today. Random Fields Autocovariance Stationarity, Isotropy. c Mikyoung Jun (Texas A&M) stat647 Lecture 2 August 30, / 13
What s for today Random Fields Autocovariance Stationarity, Isotropy c Mikyoung Jun (Texas A&M) stat647 Lecture 2 August 30, 2012 1 / 13 Stochastic Process and Random Fields A stochastic process is a family
More informationNonstationary spatial process modeling Part II Paul D. Sampson --- Catherine Calder Univ of Washington --- Ohio State University
Nonstationary spatial process modeling Part II Paul D. Sampson --- Catherine Calder Univ of Washington --- Ohio State University this presentation derived from that presented at the Pan-American Advanced
More informationKriging for processes solving partial differential equations
Kriging for processes solving partial differential equations Karl Gerald van den Boogaart 1st July 2001 Abstract Physical laws are often expressed in terms of partial differential equations. However these
More information5: MULTIVARATE STATIONARY PROCESSES
5: MULTIVARATE STATIONARY PROCESSES 1 1 Some Preliminary Definitions and Concepts Random Vector: A vector X = (X 1,..., X n ) whose components are scalarvalued random variables on the same probability
More informationIntroduction to Geostatistics
Introduction to Geostatistics Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, 2017 1 Department of Biostatistics, Bloomberg School of Public Health, Johns Hopkins University, Baltimore,
More informationLinear Models 1. Isfahan University of Technology Fall Semester, 2014
Linear Models 1 Isfahan University of Technology Fall Semester, 2014 References: [1] G. A. F., Seber and A. J. Lee (2003). Linear Regression Analysis (2nd ed.). Hoboken, NJ: Wiley. [2] A. C. Rencher and
More informationImproving Spatial Data Interoperability
Improving Spatial Data Interoperability A Framework for Geostatistical Support-To To-Support Interpolation Michael F. Goodchild, Phaedon C. Kyriakidis, Philipp Schneider, Matt Rice, Qingfeng Guan, Jordan
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2
More informationarxiv: v1 [stat.me] 3 Dec 2014
A Generalized Convolution Model and Estimation for Non-stationary Random Functions F. Fouedjio, N. Desassis, J. Rivoirard Équipe Géostatistique - Centre de Géosciences, MINES ParisTech 35, rue Saint-Honoré,
More informationLet x be an approximate solution for Ax = b, e.g., obtained by Gaussian elimination. Let x denote the exact solution. Call. r := b A x.
ESTIMATION OF ERROR Let x be an approximate solution for Ax = b, e.g., obtained by Gaussian elimination. Let x denote the exact solution. Call the residual for x. Then r := b A x r = b A x = Ax A x = A
More information1 Introduction to Generalized Least Squares
ECONOMICS 7344, Spring 2017 Bent E. Sørensen April 12, 2017 1 Introduction to Generalized Least Squares Consider the model Y = Xβ + ɛ, where the N K matrix of regressors X is fixed, independent of the
More informationA MultiGaussian Approach to Assess Block Grade Uncertainty
A MultiGaussian Approach to Assess Block Grade Uncertainty Julián M. Ortiz 1, Oy Leuangthong 2, and Clayton V. Deutsch 2 1 Department of Mining Engineering, University of Chile 2 Department of Civil &
More information1 Linear Difference Equations
ARMA Handout Jialin Yu 1 Linear Difference Equations First order systems Let {ε t } t=1 denote an input sequence and {y t} t=1 sequence generated by denote an output y t = φy t 1 + ε t t = 1, 2,... with
More informationSpatial Process Estimates as Smoothers: A Review
Spatial Process Estimates as Smoothers: A Review Soutir Bandyopadhyay 1 Basic Model The observational model considered here has the form Y i = f(x i ) + ɛ i, for 1 i n. (1.1) where Y i is the observed
More informationLecture 1: Systems of linear equations and their solutions
Lecture 1: Systems of linear equations and their solutions Course overview Topics to be covered this semester: Systems of linear equations and Gaussian elimination: Solving linear equations and applications
More informationBackward Error Estimation
Backward Error Estimation S. Chandrasekaran E. Gomez Y. Karant K. E. Schubert Abstract Estimation of unknowns in the presence of noise and uncertainty is an active area of study, because no method handles
More informationTypes of Spatial Data
Spatial Data Types of Spatial Data Point pattern Point referenced geostatistical Block referenced Raster / lattice / grid Vector / polygon Point Pattern Data Interested in the location of points, not their
More informationNotes on Time Series Modeling
Notes on Time Series Modeling Garey Ramey University of California, San Diego January 17 1 Stationary processes De nition A stochastic process is any set of random variables y t indexed by t T : fy t g
More informationLinear Regression. In this problem sheet, we consider the problem of linear regression with p predictors and one intercept,
Linear Regression In this problem sheet, we consider the problem of linear regression with p predictors and one intercept, y = Xβ + ɛ, where y t = (y 1,..., y n ) is the column vector of target values,
More informationAssessing the covariance function in geostatistics
Statistics & Probability Letters 52 (2001) 199 206 Assessing the covariance function in geostatistics Ana F. Militino, M. Dolores Ugarte Departamento de Estadstica e Investigacion Operativa, Universidad
More information