EXPLORING SPATIAL CORRELATION IN RIVERS

1. EXPLORING SPATIAL CORRELATION IN RIVERS by Joshua French

2. Introduction A city is required to extends its sewage pipelines farther in its bay to meet EPA requirements. How far should the pipelines be extended? The city doesn’t want to spend any more money than it needs to extend the pipelines. It needs to find a way to make predictions for the waste levels at different sites in the bay.

3. Usually we might try to interpolate the data using a linear model. Usually we assume observations are independent. For spatial data however, we intuitively know that response values for points close together should be more similar than points separated by a great distance. We can use the correlation between sampling sites to make better predictions with our model.

4. The Road Ahead • Methods • Introduction to the Variogram • Exploratory Analysis • Sample Variogram • Modeling the Variogram • Analysis • 3 types of results • Conclusions • Future Work

5. Introduction to the Variogram Spatial data is often viewed as a stochastic process. For each point x, a specific property Z(x) is viewed as a random variable with mean µ, variance σ2, higher-order moments, and a cumulative distribution function.

6. Each individual Z(xi) is assumed to have its own distribution, and the set {Z(x1),Z(x2),…} is a stochastic process. The data values in a given data set are simply a realization of the stochastic process. For a spatial process, second-order stationarity is often assumed.

7. Second-order stationarity implies that the mean is the same everywhere: i.e. E[Z(xj)]=µ for all points xj. It also implies that Cov(Z(xj),Z(xk)) becomes a function of the distance xj toxk.

8. Thus, Cov(Z(xj),Z(xk)) = Cov(Z(x),Z(x+h)) = Cov(h) where h measures the distance between two points.

9. Looking at the variance of differences Var[Z(x)-Z(x+h)] =E[ (Z(x)-Z(x+h))2 ] = 2 γ(h) Assuming second-order stationarity, γ(h)=Cov(0)-Cov(h). γ(h) is known as the semi-variogram. The plot of γ(h) on h is known as the variogram.

10. Things to know about variograms: • γ(h)= γ(-h). Because it is an even function, usually only positive lag distances are shown. • Nugget effect - by definition, γ(0)= 0. In practice however, sample variograms often have a positive value at lag 0. This is called the “nugget effect”.

11. Tend to increase monotonically • Sill – the maximum variance of the variogram • Range – the lag distance at which the sill is reached. Observations are not correlated past this distance. The following figure shows these features

12. Variogram Example

13. Exploratory Analysis The data studied is the longitudinal profile of the Ohio River. Instead of worrying about the river network with streams, tributaries, and other factors, we simply look at the Ohio River as a one-dimensional object.

14. The Ohio River

15. Longitudinal Profile of the Ohio River Sampling Sites

16. Before we model variograms, we should explore the data. We need to make sure that the data analyzed satisfies second-order stationarity If there is an obvious trend in the data, we should remove it and analyze the residuals. If the variance increases or decreases with lag distance, then we should transform the variable to correct this.

17. It is fairly easy to check for stationarity of this data set using a scatter plot.

18. If the data contains outliers, we should do analysis both with and without outliers present. If G1>1, then we should transform the data to approximate normality if possible.

19. 3.3 The Sample Variogram One of the previous definitions of semivariance is: The logical estimator is: where N(h) is the number of pairs of observations associated with that lag.

20. Sample Variogram Example

21. Modeling the Variogram Our goal is to estimate the true variogram of the data. There were four variogram models used to model the sample variogram: the spherical, Gaussian, exponential, and Matern models.

22. Variogram Models

23. Analysis The data analyzed is a set of particle size and biological variables for the Ohio River. The data was collected by “The Ohio River Valley Sanitation Commission. This is better known as ORSANCO. There were between 190 and 235 unique sampling sites, depending on the variable.

24. ORSANCO data collection

25. The results of the analysis fell into three main groups: • Able to fit the sample variogram well • Not able to fit the sample variogram well • Analysis not reasonable

26. Good Results: Number of Individuals at a site After correcting for skewness by doing a log transformation, there are a number of outliers. We analyze the data both with and without the outliers.

27. log(Num Individuals) Sample Variogramwith outliers

28. log(Num Individuals) Sample Variogramwithout outliers

29. We were not able to model the sample variogram perfectly, but we were able to detect some amount of spatial correlation in the data, especially when the outliers were removed. We are able to obtain reasonable estimates of the nugget, sill, and variance.

30. Poor Results: Percent Sand After doing exploratory spatial analysis and removing a trend, we fit the sample variogram of the percent sand residuals.

31. Sample Variogram of percent sand residuals

32. The sample variogram does not really increase monotonically with distance. Our variogram models cannot fit this very well. Though we can obtain estimates of the nugget, sill, and range, the estimates cannot be trusted.

33. No results: Percent Hardpan This variable was so badly skewed that analysis was not reasonable. The skewness coefficient is 12.38. This is extremely high.

34. QQplot of Percent Hardpan

35. Scatter plot of Percent Hardpan

36. The data is nearly all zeros! There is also an erroneous data value. A percentage cannot be greater than 100%. Data analysis does not seem reasonable. Our data does not meet the conditions necessary to use the spatial methods discussed.

37. Conclusions Able to fit sample variogram reasonably well – percent gravel, number of individuals, number of species Not able to fit sample variogram well – percent sand, percent detritivore, percent simple lithophilic individuals, percent invertivore No results – remaining variables

38. Summary of Results

39. Future Work Things to consider in future analysis: - The water flows in only one-direction. A point downstream cannot affect a point upstream - Natural features such as tributaries may impact spatial correlation - Manmade features such as dams may impact spatial correlation

40. Concluding Thought Before you criticize someone, you should walk a mile in their shoes. That way, when you criticize them, you’re a mile away and you have their shoes. - Jack Handey