Geo 479579 Geostatistics Ch 17 Cokriging Introduction Data

Geo 479/579: Geostatistics Ch 17. Cokriging

Introduction Data sets often contain more than one variable of interest These variables are usually spatially crosscorrelated

The Cokriging System A method for estimation that minimizes the variance of the estimation error by exploiting the crosscorrelation between several variables Cross-correlated information contained in the secondary variable should help reduce the variance of the estimation errors

The Cokriging System When is the secondary variable useful in estimates? Primary variable of interest is under sampled then the only information we have is the cross correlated information

The Cokriging System The cokriging estimate is a linear combination of both primary and secondary data values This is the Equation used in Ordinary Kriging pg 279

The Cokriging System The development of the cokriging system is identical to the development of ordinary kriging system Estimation Error R can be defined as This is a modification of the error estimation in Ordinary Kriging( pg 279)

The Cokriging System Using matrix notation we can write w = { a 1, a 2, a 3, …an, b 1, b 2, b 3, …bm} Z = { U 1, U 2…. . Ui, V 1, …. Vj}

The Cokriging System Using Equation 9. 14 (p 216), 12. 6 (p 283), 16. 3 (p 372) we can write

This is similar to Chapter 16

The Cokriging System 1) Unbiasedness condition v Note: Other nonbias conditions are also possible

It is similar to unbiasedness in Ordinary Kriging (p 281) § We set error at as 0:

The Cokriging System 2) Minimizing error variance Lagrangean Relaxation:

The Cokriging System Lagrangean Relaxation: Original Lagrange parameter: (12. 9)

The Cokriging System Equating n+m+2 partial derivatives of Var{R} to zero, we get the following system of equations

§ This is similar to minimizing the varianves of error in Ordinary Kriging § The set of weights that minimize the error variance under the unbiasedness condition satisfies the following n+1 equations - ordinary kriging system: (12. 11) (12. 12)

Minimization of the Error Variance § The ordinary kriging system expressed in matrix (12. 13) (12. 14)

The Cokriging System Positive definiteness must hold for the set of auto- and cross-variograms (Eq 16. 44, p 391).

The Cokriging System If the primary and secondary variables both exist at all data locations and the auto- and crossvariograms are proportional to the same basic model then the cokriging estimates will be identical to the ordinary kriging estimates

A Cokriging Example

A Case Study v Compares cokriging and ordinary kriging v Undersampled variable U is estimated using 275 U & 470 V sample data for cokriging and only the 275 U data for ordinary kriging

Ordinary kriging 275 U values Using eq 17. 11 for the variogram model

Cokriging 275 U and 470 V values Using eq 17. 11 for the variogram model Two non-bias conditions 1) uses the initial conditions 2) uses only one nonbias condition

In the alternate unbiased condition, the unknown U value is now estimated as a weighted linear combination of nearby V values adjusted by a constant so that their mean is equal to the mean of the U values

Negative estimates occur because of the nonbias condition

Cokriging with two nonbias conditions is less than satisfactory A physical process with both negative and positive weighting scheme is difficult to imagine Cokriging with one nonbias condition considerably improved the spread of errors and bias Though we had to calculate global means of U and V

If the spatial continuity is modeled using semivariograms then they can be converted to covariance values for cokriging matrices by following equation:

The Cokriging System If we want an estimate over a local area A, there are two options: 1) Average of point estimations within A

The Cokriging System 2) Replace all the covariance terms and in point cokriging system, with average covariance values and

The Cokriging System With the unbiasedness conditions, we can calculate error variance as follows