Nonstationary covariance structures II

1. NRCSE Nonstationary covariance structures II

2. Drawbacks with deformation approach Oversmoothing of large areas Local deformations not local part of global fits Covariance shape does not change over space Limited class of nonstationary covariances

3. Nonstationary spatial covariance: Basic idea: the parameters of a localvariogram model–nugget, range, sill, and anisotropy–vary spatially. For deformation, which do? Not nugget nor sill.

4. Major approaches: Haas, 1990, Moving window kriging Kim, Mallock & Holmes, 2005, Piecewise Gaussian modeling Nott & Dunsmuir, 2002, Biometrika, Average of locally stationary processes Fuentes, 2002, Kernel averaging of orthogonal, locally stationary processes. Pintore & Holmes, 2005, Fourier and Karhunen-Loeve expansions Higdon & Swall, 1998, 2000, Gaussian moving averages or “process convolution” model Nychka, Wikle & Royle, 2002. Wavelet expansion.

5. Thetford revisited:Spectrum of canopy heights Ridge in neg quadrat 46.4% of var Low freq ridge in pos quadrat 19.7% of var High freq ridge in pos quadrat 6.1% of var

6. Looking at several adjacent plots Features depend on spatial location HF q<0 LF q>0 LF q<0 HF q>0 other

7. Moving window approach For kriging at s, consider only sites “near” s Near defined small enough to make stationarity/isotropy reasonable

9. The moving window approach In purely spatial context, pick the ns nearest sites (for spacetime, use a window in space and time and grow until get ns sites) Fit a model based only on the chosen sites and perform appropriate kriging Difficulty: not a single model, so can be discontinuous Overall covariance may not be positive definite

10. Kernel averaging Fuentes (2000): Introduce uncorrelated stationary processes Zk(s), k=1,...,K, defined on subregions Sk and construct where wk(s) is a weight function related to dist(s,Sk). Then

11. Spectral version so where Hence

12. Estimating spectrum Asymptotically

13. Details K = 9; h = 687 km Mixture of Matérn spectra

14. An example Models-3 output, 81x87 grid, 36km x 36km. Hourly averaged nitric acid concentration week of 950711.

15. Models-3 fit

16. A spectral approach to nonstationary processes Spectral representation: s slowly varying square integrable, Y uncorrelated increments Hence is the space-varying spectral density Observe at grid; use FFT to estimate in nbd of s

17. Testing for nonstationarity U(s,w) = log has approximately mean f(s,w) = log fs(w) and constant variance in s and w. Taking s1,...,sn and w1,...,wm sufficently well separated, we have approximately Uij = U(si,wj) = fij+eij with the eij iid. We can test for interaction between location and frequency (nonstationarity) using ANOVA.

18. Details The general model has The hypothesis of no interaction (ij=0) corresponds to additivity on log scale: (uniformly modulated process: Z0 stationary) Stationarity corresponds to Tests based on approx 2-distribution (known variance)

19. Models-3 revisited

20. Moving averages A simple way of constructing stationary sequences is to average an iid sequence . A continuous time counterpart is , where x is a random measure which is stationary and assigns independent random variables to disjoint sets, i.e., has stationary and independent increments.

22. Covariance In the squared exponential kernel case:

23. Lévy-Khinchine General representation for psii: n is the Lévy measure, and xt is the Lévy process. We can construct it from a Poisson measure H(du,ds) on R2 with intensity E(H(du,ds))=n(du)ds and a standard Brownian motion Bt as

24. Examples Brownian motion with drift: xt~N(mt,s2t) n(du)=0. Poisson process: xt~Po(lt) m=s2=0, n(du)=ld{1}(du) Gamma process: xt~G(at,b) m=s2=0, n(du)=ae-bu1(u>0)du/u Cauchy process: m=s2=0, n(du)=bu-2du/p

25. Spatial moving averages We can replace R for t or s with R2 (or a general metric space) We can replace b(t-s) by bt(s) to relax stationarity We can let the intensity measure for H be an arbitrary measure n(ds,du)

26. Gaussian moving averages Higdon (1998), Swall (2000): Let x be a Brownian motion without drift, and . This is a Gaussian process with correlation Account for nonstationarity by letting the kernel b vary with location:

28. Details yields an explicit covariance function which is squared exponential with parameter changing with spatial location. The prior distribution describes “local ellipses,” i.e., smoothly changing random kernels.

29. Prior parametrization Bivariate normal covariance can be described by an ellipse (level curve), determined by area, center and one focus. Focus chosen independently Gaussian with isotropic squared exponential covariance. Another parameter describes the range of influence (scale) of a given ellipse. Prior gamma.

30. Local ellipses Focus Ellipse determined by focus, center, and fixed area Focus components independent Gaussian time series

31. Posterior Compute posterior on a grid Mostly Metropolis-Hastings algorithm Full model allows the size of the kernels to vary across locations

32. Example Piazza Road Superfund cleanup. Dioxin applied to road to prevent dust storms seeped into groundwater.

33. Posterior distribution