Bo Li. National Center for Atmospheric Research. Based on joint work with:
|
|
- Eileen Georgina Nash
- 5 years ago
- Views:
Transcription
1 Nonparametric Assessment of Properties of Space-time Covariance Functions and its Application in Paleoclimate Reconstruction Bo Li National Center for Atmospheric Research Based on joint work with: Marc Genton: University of Geneva & TAMU Mike Sherman: TAMU Doug Nychka and Caspar Ammann: NCAR Supported by the National Science Foundation
2
3 Outline of this talk Testing properties of space-time covariance function 1. Introduction 2. Development of testing methodology 3. Data analysis and simulations Reconstructing paleoclimate Apply the testing procedures to solve the problem 1. Significance and data 2. Why use BHM framework? 3. Challenges and recipe
4 Introduction Spatio-temporal data sets are common in geophysical and environmental sciences strictly stationary space-time random field {Z(s, t) : s R d, t R} C(h, u) = cov{z(s, t), Z(s + h, t + u)} Modeling such data often relies on parametric covariance models (Cressie and Huang, 1999; Gneiting, 2002; Stein, 2005) Make assumptions to simplify the covariance function in large and complex space-time data sets
5 Ways of simplifying covariance function Full symmetry: C(h, u) = C(h, u), or if C(h, u) = C( h, u) t 2 wind direction t 1 z(s 1, t 2 ) z(s 2, t 2 ) z(s 1, t 1 ) h z(s 2, t 1 ) Z: CO 2 concentration h = s 2 s 1, u = t 2 t 1 Cov(Z(s 1, t 1 ), Z(s 2, t 2 )) = C(h, u) = C( h, u) Cov(Z(s 2, t 1 ), Z(s 1, t 2 )) = C(h, u) = C( h, u)? C(h, u) = C(h, u)
6 Ways of simplifying covariance function Separability: C(h, u) = C(h, 0)C(0, u) C(0, 0) C(h, 0): pure spatial covariance C(0, u): pure temporal covariance Separability = Full symmetry Taylor Hypothesis: There exists a velocity vector v such that C(vu, 0) = C(0, u) for all u...
7 Ways of simplifying covariance function Important because they simplify the structure of the model and its inference, and ease possibly extensive computations Random field: Z(s i, t j ), i = 1,..., 100, j = 1,..., 100 Size of space-time covariance Σ: 10,000 10,000 If Σ is separable: Σ = Σ 1 Σ 2 Σ 1 = Σ 1 1 Σ 1 2 size of spatial covariance Σ 1 : size of temporal covariance Σ 2 :
8 Ways of simplifying covariance function Not always appropriate: prevailing winds, water flows, and atmosphere circulation in geoscience, meteorology and ecology Pacific Ocean wind data (Cressie and Huang,1999) Irish wind speed data (Gneiting et al., 2007) Sulfate concentration levels (Jun and Stein, 2007) Their graphic-based diagnostics are useful but difficult to assess and open to interpretation Need a formal and unified approach...
9 Existing approaches to test specific properties Developed the test in a particular situation or with the assumption of distribution Shitan and Brockwell (1995), Guo and Billard (1998), Genton and Koul (2007): test for separability in SAR/QAR Lu and Zimmerman (2001) and Guan, Sherman and Calvin (2004): nonparametric tests for spatial isotropy Scaccia and Martin (2002, 2005): spectral method to test the sym. and sep. for spatial lattice processes Fuentes (2006): nonparametric test for sep. of a spatio-temporal process based on a spectral method Mitchell, Genton and Gumpertz (2005, 2006): LRT for sep. in the context of multivariate repeated measures assuming multinormality, with extension to space-time
10 Outline of this talk Testing properties of space-time covariance function 1. Introduction 2. Development of testing methodology 3. Data analysis and simulations Reconstructing paleoclimate Apply the testing procedures to solve the problem 1. Significance and data 2. Why use BHM framework? 3. Challenges and recipe
11 Hypotheses Λ: a set of space-time lags (h, u) Ĉ(h, u): an estimator of C(h, u) G = {C(h, u), (h, u) Λ}, Ĝ = {Ĉ(h, u), (h, u) Λ} H 0 : Af(G) = 0 where A: contrast matrix of row rank q, and f = (f 1,..., f r ) T : real-valued functions, differentiable at G (Lu and Zimmerman, 2001) Full symmetry: H 1 0 : C(h, u) C(h, u) = 0, (h, u) Λ Separability: H 2 0: C(h, u)/c(h, 0) C(0, u)/c(0, 0) = 0, (h, u) Λ
12 Examples A 1 G = 0 and A 2 f(g) = 0 Λ = {(0, 0), (h 1, u 1 ), (h 1, u 1 ), (h 2, u 2 ), (h 2, u 2 ), (h 1, 0), (0, u 1 ), (h 2, 0), (0, u 2 )} G = (C(0, 0), C(h 1, u 1 ), C(h 1, u 1 ), C(h 2, u 2 ), C(h 2, u 2 ), f(g) = C(h 1, 0), C(0, u 1 ), C(h 2, 0), C(0, u 2 )) T A 1 = [ ( ) T C(h 1,u 1 ) C(h, C(h 1, u 1 ) 1,0) C(h, C(h 2,u 2 ) 1,0) C(h, C(h 2, u 2 ) 2,0) C(h, C(0,u 1), C(0,u 2) 2,0) C(0,0) C(0,0) A 2 = ]
13 Asymptotics of sample space-time covariances In many situations, the observations are taken from a fixed space S R d at regularly spaced times T n = {1,..., n} Increasing index sets D n = S T n Ĉn(h, u): sample based estimator of C(h, u) based on obs. in a sequence of increasing index sets D n Ĝn = {Ĉn(h, u) : (h, u) Λ}: estimator of G computed over D n. We make no assumptions on the marginal or the joint distribution of observations other than mild moment and mixing conditions on the random field
14 Asymptotics of sample space-time covariances Under zero mean assumption for Z: Ĉ n (h, u) = 1 n S(h) n u S(h) t=1 Z(s, t)z(s + h, t + u) Strong mixing condition: α(u) = O(u ɛ ) for some ɛ > 0 Moment condition for Ĉn(h, u): sup n E{ n{ĉn(h, u) C(h, u)} 2+δ } < for some δ > 0
15 Asymptotics of sample space-time covariances Theorem: Let {Z(s, t), s R d, t Z} be a strictly stationary spatio-temporal random field observed in D n = S T n, where S R d and T n = {1,..., n}. Assume t Z cov{z(0, 0)Z(h 1, u 1 ), Z(s, t)z(s + h 2, t + u 2 )} < for all h 1, h 2 S, s S(h 2 ) and all finite u 1, u 2, then Σ = lim n ncov(ĝn, Ĝn) exists. Further with mixing and moment conditions we have n(ĝn G) d N m (0, Σ) as n.
16 Test statistics Li, Genton and Sherman (2008): a n (Ĝn G) d N m (0, Σ) under mild moment and mixing conditions, and for a wide variety of data structures d a n {f(ĝn) f(g)} N r (0, B T ΣB) where B ij i = 1,..., m, j = 1,..., r = f j / G i Test Statistic for H 0 : Af(G) = 0: TS = a 2 n{af(ĝn)} T (AB T ΣBA T ) 1 {Af(Ĝn)} Examples: full symmetry and separability d χ 2 q Estimate the matrix B empirically (replacing G with Ĝ n ) and estimate Σ using subsampling techniques
17 Estimate Σ using Subsampling Data observed over fixed spatial domain S and increasing time domain T n Form overlapping S l(n) subblocks using moving subblock window along time Choice of optimal block length, l(n), in the sense of minimizing the MSE of estimators in a variety of contexts (e.g., Lahiri, 2003) For simplicity we follow Carlstein (1986): l(n) = ( 2γ 1 γ 2 ) 2/3 ( 3n 2 )1/3, where we estimate γ by γ n = Ĉ n (0, 1)/Ĉn(0, 0) but often works well in practice, see Lahiri (2003)
18 Choice of contrast matrices Choice of A is not unique For example: pick only the first two rows of A to form a new test statistic whose asymptotic distribution follows χ 2 2 Although these tests will have approximately the same size, the power depends on the specific choice Generally: preferable to use lags combining small spatial and temporal lags Take features of physical process into consideration while choosing testing lags
19 Outline of this talk Testing properties of space-time covariance function 1. Introduction 2. Development of testing methodology 3. Data analysis and simulations Reconstructing paleoclimate Apply the testing procedures to solve the problem 1. Significance and data 2. Why use BHM framework? 3. Challenges and recipe
20 Irish wind speed data Time series of daily average wind speed at 11 synoptic meteorological stations in Ireland during the period Prevailing westerly wind Gneiting et al.(2007) fit spatio-temporal model C(h, u) = [ { } (1 ν)(1 λ) c h exp 1 + a u 2α (1 + a u 2α ) β/2 + ν 1 ν δ h=0 ] +λ ( 1 1 2ν h 1 vu symmetry parameter λ, space-time interaction parameter β both take values in [0,1] λ = 0: fully symmetric; λ = 0, β = 0: separable Estimates from training data: ˆλ = and ˆβ = ) +
21 Testing full symmetry and separability 5 pairs of stations among the 55 pairs and time lags u = 1 and 2 days 1) prevailing westerly wind suggests choosing the 5 pairs of stations with the smallest ratio of h 2 /h 1 (h 1 WE; h 2 NS) 2) the 5 pairs of stations with the largest ratio of h 2 /h 1 3) the 5 pairs with the shortest spatial distance h among the 55 pairs
22 Malin Head smallest h 2 h 1 largest h 2 h 1 shortest h Belmullet Clones Claremorris Mullingar Birr Dublin Irish wind stations and the five pairs selected for the test Shannon Kilkenny Valentia Roche s Point
23 Testing full symmetry and separability TS1: test statistic for full symmetry TS2: test statistic for separability 1) TS1 = 262.7, p-value = 0; TS2 = 445.2, p-value = 0 reject assumptions of full symmetry and separability 2) TS1 = 20.2, TS2 = ) TS1 = 132.9, TS2 = Results in 1) detect further departure from full symmetry and separability
24 Assessing size and power Empirical sizes and powers for testing full symmetry and separability in the Irish wind speed data λ β Full Symmetry Separability 0 0 size=0.074 [0.067] (0.083) size=0.084 [0.073] (0.091) 0 1 size=0.070 [0.081] (0.092) pow=0.945 [1.000] (0.503) pow=1.000 [0.141] (1.000) pow=1.000 [0.510] (1.000) pow=1.000 [0.309] (1.000) pow=1.000 [0.959] (1.000) pow=0.900[0.095](0.529) pow=0.985 [0.125] (0.445) size=0.095 [0.075] (0.078) pow=0.717 [0.929] (0.306) pow=0.876[0.083] (0.478) pow=1.000[0.941] (0.867) Nominal level is Sizes/powers in [ ] are obtained by using the five station pairs with the largest ratio of h 2 /h 1, and sizes/powers in ( ) are obtained by using the five station pairs with the shortest h.
25 Features of the testing methodology Unified framework for testing a variety of assumptions for covariance functions Can be applied to a wide range of data sets Easy to implement Detection of asymmetric and nonseparable features in Irish wind speed data Simulation experiments show reliability and accuracy
26 Extensions for climate data Test Taylor s hypothesis in precipitation Extend to multivariate random fields: e.g., temperature, precipitation and geopotential heights in climatology Test properties of cross-covariances C ij (k), i, j: index of variates, k: lag vector C ij (k) behaves differently than C(k)
27 Extensions for climate data New definitions of simplications: Multivariate symmetry: C ij (k) = C ji (k) or C ij (k) = C ij ( k) Multivariate Separability: C ij (k) = ρ(k)t for a spatial correlation ρ(k) and covariance matrix T between variates Linear model of coregionalization (LMC): C(k) = r g=1 ρ g (k)t g
28 Outline of this talk Testing properties of space-time covariance function 1. Introduction 2. Development of testing methodology 3. Data analysis and simulations Reconstructing paleoclimate Apply the testing procedures to solve the problem 1. Significance and data 2. Why use BHM framework? 3. Challenges and recipe
29 Our paleoclimate team Caspar Ammann, Doug Nychka: NCAR Gene Wahl: Alfred University Jack Williams: University of Wisconsin-Madison Malcolm Hughes, Michael Evans: University of Arizona Mark Berliner: Ohio State University Robert Harris: Oregon State University
30 Climate change
31 Climate change
32 Climate change
33 The problem and the data Important to our life, understand climate change New research problem for statisticians Observed: Meterological variables (temperature, precipitation, 500hPa geopotential height fields) 1850-present Proxies: Tree rings, Pollens, Ice core, Borehole... past millennium Unobserved: Meteorological variables before 1849
34 Tree rings
35 Ice core
36 Pollen
37 Borehole
38 Proxies Each proxy may be skillful at different time scales Tree ring: decadal Pollen: bi-decadal to semi-centennial Borehole: centennial and onward Goal: Reconstruct the past NH temperature, and NA temperature, precipitation and 500 hpa heights field. Exploit all three sources of proxies Incorporate energy balance and external drivers
39 Bayesian Hierarchical Model (BHM) A natural framework: Three hierarchies: Data Stage: [Proxies Geophysical Process, Parameters] Likelihood of Proxies given underlying process Process Stage: [Geophysical Process Parameters] Physical model of Geophysical Process Parameter Stage: [Parameters] Specify the prior of parameters
40
41 Illustration of using BHM One type of proxy, say, 14 tree rings Construct past temperature only Target is northern hemisphere mean temperature, i.e. single time series of temperature No spatial locations of proxy and temperature are considered Data Model: P t T t = β 0 + β 1 T t + e t, e t AR(2)
42 Process Model Radiative forcings: Explosive volcanism, Solar Activity Changes and Anthropogenic forcings Force
43 Process Model µ: Radiative forcings T t µ t = a + bµ t + e t, e t AR(2) Process model is essentially an AR(2) time series as the target is NH mean temperature Bring the knowledge of radiative forcings into the reconstruction
44 Prior β 0(i) N(0, σ 2 0), β 1(i) N(0, σ 2 1) a N(0, σ 2 a), b N(0, σ 2 b ) φ 2ɛ unif( 1, 1), φ 1ɛ φ 2ɛ (1 φ 2ɛ ) unif( 1, 1) φ 2 unif( 1, 1), φ 1 φ 2 (1 φ 2 ) unif( 1, 1) σ 2 P IG(q P, r P ), σ 2 T IG(q T, r T )
45 An example: simulated numerical data Output from global coupled climate model Degrees C pseudo-proxy series sampled from 14 individual grid boxes in the climate model
46 An example of BHM: simulated numerical data Input: 14 pseudo-proxies and full model temperature Reconstruct the past temperature ( ) Degree C
47 Conclusion from the example The posterior mean matches the trend of the numerical data The numerical data is within the 95% prediction band The posterior mean of parameters are close to those directly estimated from numerical data BHM works well in reconstructing the past temperature
48 Apply this simple BHM to 14 real proxies Input: 14 proxies (mostly are tree rings) and the instrumental temperature Reconstruct the past temperature ( ) Degree C
49 Reconstruct real world NH temperature Uncertainty of decadal average and decadal maximum Degree C
50 Current work Data Stage: [Proxies Temperature field, Parameters] Some issues: Proxies {}}{ TreeRing, Pollen, Borehole 1. dating errors of pollen and borehole 2. forward algorithm to convert the surface temperature to pollen assemblage and borehole depth temperature Challenge in process stage: spatio-temporal model of temperature random field
51 Current work Take advantage from long numerical simulations of climate model Apply the unified testing procedures to determine an appropriate space-time covariance function Spatially isotropic? Space-time symmetric? Space-time separable?... Hope all the answers are Yes!
52 Future work Data Stage: [Proxies Geophysical process, Parameters] Geophysical process {}}{ Temperature, Precipitation, 500hPa heights Need in process stage: A multivariate space-time geophysical process model Guide to choose this model: Multivariate symmetric? Multivariate separable? If unfortunately LMC is needed, how many components are necessary?
53 Summary Unified framework for testing a variety of assumptions for covariance functions of stationary spatiotemporal random fields Bayesian hierarchical models (BHM) to reconstruct the past Northern Hemisphere temperatures Building BHM for reconstruction of multivariate spacetime climate random field based on testing procedures
54
55 References for this talk Li, B., Genton, M. G. and Sherman, M., A nonparametric assesment of properties of space-time covariance functions. (2007) JASA, Vol. 102, Li, B., Nychka, W. D. and Ammann, C. M., The Hockey Stick and the 1990s: A statistical perspective on reconstructing hemispheric temperatures. (2007), Tellus, Vol. 59, Li, B., Genton, M. G. and Sherman, M., On the asymptotic joint distribution of sample space-time covariance estimators. (2008) Bernoulli, in press Li, B., Genton, M. G. and Sherman, M., Testing the covariance structure of multivariate random fields, under review by Biometrika Li, B., Murthi, A., Bowman, K. P., North, G. R., Genton, M. G.,and Sherman, M., Statistical tests of Taylor s hypothesis: An application to precipitation fields. (2007) under review by Journal of Hydrometeorology
Predictive spatio-temporal models for spatially sparse environmental data. Umeå University
Seminar p.1/28 Predictive spatio-temporal models for spatially sparse environmental data Xavier de Luna and Marc G. Genton xavier.deluna@stat.umu.se and genton@stat.ncsu.edu http://www.stat.umu.se/egna/xdl/index.html
More informationA TEST FOR STATIONARITY OF SPATIO-TEMPORAL RANDOM FIELDS ON PLANAR AND SPHERICAL DOMAINS
Statistica Sinica 22 (2012), 1737-1764 doi:http://dx.doi.org/10.5705/ss.2010.251 A TEST FOR STATIONARITY OF SPATIO-TEMPORAL RANDOM FIELDS ON PLANAR AND SPHERICAL DOMAINS Mikyoung Jun and Marc G. Genton
More informationA test for stationarity of spatio-temporal random fields on planar and spherical domains
A test for stationarity of spatio-temporal random fields on planar and spherical domains Mikyoung Jun and Marc G. Genton 1 June 13, 2010 ABSTRACT: A formal test for weak stationarity of spatial and spatio-temporal
More informationGeostatistical Space-Time Models, Stationarity, Separability and Full Symmetry. Technical Report no. 475
Geostatistical Space-Time Models, Stationarity, Separability and Full Symmetry Tilmann Gneiting, 1,2 Marc G. Genton 3 and Peter Guttorp 1,4 1 Department of Statistics, University of Washington Seattle,
More informationSymmetry and Separability In Spatial-Temporal Processes
Symmetry and Separability In Spatial-Temporal Processes Man Sik Park, Montserrat Fuentes Symmetry and Separability In Spatial-Temporal Processes 1 Motivation In general, environmental data have very complex
More informationWhat s for today. Introduction to Space-time models. c Mikyoung Jun (Texas A&M) Stat647 Lecture 14 October 16, / 19
What s for today Introduction to Space-time models c Mikyoung Jun (Texas A&M) Stat647 Lecture 14 October 16, 2012 1 / 19 Space-time Data So far we looked at the data that vary over space Now we add another
More informationNonstationary cross-covariance models for multivariate processes on a globe
Nonstationary cross-covariance models for multivariate processes on a globe Mikyoung Jun 1 April 15, 2011 Abstract: In geophysical and environmental problems, it is common to have multiple variables of
More informationSeparable approximations of space-time covariance matrices
ENVIRONMETRICS Environmetrics 2007; 18: 681 695 Published online 8 May 2007 in Wiley InterScience (www.interscience.wiley.com).854 Separable approximations of space-time covariance matrices Marc G. Genton
More informationNon-stationary Cross-Covariance Models for Multivariate Processes on a Globe
Scandinavian Journal of Statistics, Vol. 38: 726 747, 2011 doi: 10.1111/j.1467-9469.2011.00751.x Published by Blackwell Publishing Ltd. Non-stationary Cross-Covariance Models for Multivariate Processes
More informationGeostatistical Space-Time Models, Stationarity, Separability, and Full Symmetry
CHAPTER 4 Geostatistical Space-Time Models, Stationarity, Separability, and Full Symmetry Tilmann Gneiting, Marc G. Genton, and Peter Guttorp Contents 4.1 Introduction...151 4.2 Geostatistical space-time
More informationA General Class of Nonseparable Space-time Covariance
A General Class of Nonseparable Space-time Covariance Models Thaís C.O. Fonseca and Mark F.J. Steel Department of Statistics, University of Warwick, U.K. Summary. The aim of this work is to construct nonseparable,
More informationProbabilistic Models of Past Climate Change
Probabilistic Models of Past Climate Change Julien Emile-Geay USC College, Department of Earth Sciences SIAM International Conference on Data Mining Anaheim, California April 27, 2012 Probabilistic Models
More informationTHE HOCKEY STICK CONTROVERSY AND PCA
THE HOCKEY STICK CONTROVERSY AND PCA 1 In 1998/1999, Mann, Bradley and Hughes (MBH) wrote two papers reconstructing temperatures over the last 1000 years from proxy data (tree rings, corals, ice cores
More informationModels for models. Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research
Models for models Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research Outline Statistical models and tools Spatial fields (Wavelets) Climate regimes (Regression and clustering)
More informationNon-gaussian spatiotemporal modeling
Dec, 2008 1/ 37 Non-gaussian spatiotemporal modeling Thais C O da Fonseca Joint work with Prof Mark F J Steel Department of Statistics University of Warwick Dec, 2008 Dec, 2008 2/ 37 1 Introduction Motivation
More informationA Framework for Daily Spatio-Temporal Stochastic Weather Simulation
A Framework for Daily Spatio-Temporal Stochastic Weather Simulation, Rick Katz, Balaji Rajagopalan Geophysical Statistics Project Institute for Mathematics Applied to Geosciences National Center for Atmospheric
More informationPREDICTIVE SPATIO-TEMPORAL MODELS FOR SPATIALLY SPARSE ENVIRONMENTAL DATA
Statistica Sinica 15(25), 547-568 PREDICTIVE SPATIO-TEMPORAL MODELS FOR SPATIALLY SPARSE ENVIRONMENTAL DATA Xavier de Luna and Marc G. Genton Umeå University and Texas A&M University Abstract: We present
More informationReply MICHAEL E. MANN
15 JULY 2007 N O T E S A N D C O R R E S P O N D E N C E 3699 Reply MICHAEL E. MANN Department of Meteorology, and Earth and Environmental Systems Institute, The Pennsylvania State University, University
More informationClimate Change: the Uncertainty of Certainty
Climate Change: the Uncertainty of Certainty Reinhard Furrer, UZH JSS, Geneva Oct. 30, 2009 Collaboration with: Stephan Sain - NCAR Reto Knutti - ETHZ Claudia Tebaldi - Climate Central Ryan Ford, Doug
More informationStatistical Science: Contributions to the Administration s Research Priority on Climate Change
Statistical Science: Contributions to the Administration s Research Priority on Climate Change April 2014 A White Paper of the American Statistical Association s Advisory Committee for Climate Change Policy
More informationThe Temperature Proxy Controversy
School of Mathematics February 8, 2012 Overview Introduction 1 Introduction 2 3 4 A provocative statement Figures often beguile me, particularly when I have the arranging of them myself; in which case
More informationThis is a repository copy of Estimation and testing for covariance-spectral spatial-temporal models.
This is a repository copy of Estimation and testing for covariance-spectral spatial-temporal models. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/89424/ Version: Accepted
More informationSurface Temperature Reconstructions for the Last 2,000 Years. Statement of
Surface Temperature Reconstructions for the Last 2,000 Years Statement of Gerald R. North, Ph.D. Chairman, Committee on Surface Temperature Reconstructions for the Last 2,000 Years National Research Council
More informationWeather generators for studying climate change
Weather generators for studying climate change Assessing climate impacts Generating Weather (WGEN) Conditional models for precip Douglas Nychka, Sarah Streett Geophysical Statistics Project, National Center
More informationSpace-time analysis using a general product-sum model
Space-time analysis using a general product-sum model De Iaco S., Myers D. E. 2 and Posa D. 3,4 Università di Chieti, Pescara - ITALY; sdeiaco@tiscalinet.it 2 University of Arizona, Tucson AZ - USA; myers@math.arizona.edu
More informationAn Introduction to Spatial Statistics. Chunfeng Huang Department of Statistics, Indiana University
An Introduction to Spatial Statistics Chunfeng Huang Department of Statistics, Indiana University Microwave Sounding Unit (MSU) Anomalies (Monthly): 1979-2006. Iron Ore (Cressie, 1986) Raw percent data
More informationNon-Fully Symmetric Space-Time Matérn Covariance Functions
Non-Fully Symmetric Space-Time Matérn Covariance Functions Tonglin Zhang and Hao Zhang Purdue University Abstract The problem of nonseparable space-time covariance functions is essential and difficult
More informationTesting Lack of Symmetry in Spatial-Temporal Processes. Man Sik Park and Montserrat Fuentes 1
Testing Lack of Symmetry in Spatial-Temporal Processes Man Sik Park and Montserrat Fuentes 1 Institute of Statistics Mimeo Series# 2585 SUMMAY Symmetry is one of the main assumptions that are frequently
More informationThe Matrix Reloaded: Computations for large spatial data sets
The Matrix Reloaded: Computations for large spatial data sets The spatial model Solving linear systems Matrix multiplication Creating sparsity Doug Nychka National Center for Atmospheric Research Sparsity,
More informationCross-covariance Functions for Tangent Vector Fields on the Sphere
Cross-covariance Functions for Tangent Vector Fields on the Sphere Minjie Fan 1 Tomoko Matsuo 2 1 Department of Statistics University of California, Davis 2 Cooperative Institute for Research in Environmental
More informationHalf-Spectral Space-Time Covariance Models
Half-Spectral Space-Time Covariance Models Michael T. Horrell & Michael L. Stein University of Chicago arxiv:1505.01243v1 [stat.me] 6 May 2015 September 18, 2018 Abstract We develop two new classes of
More informationHierarchical Modeling for Multivariate Spatial Data
Hierarchical Modeling for Multivariate Spatial Data Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department
More informationSpatio-temporal precipitation modeling based on time-varying regressions
Spatio-temporal precipitation modeling based on time-varying regressions Oleg Makhnin Department of Mathematics New Mexico Tech Socorro, NM 87801 January 19, 2007 1 Abstract: A time-varying regression
More informationA full scale, non stationary approach for the kriging of large spatio(-temporal) datasets
A full scale, non stationary approach for the kriging of large spatio(-temporal) datasets Thomas Romary, Nicolas Desassis & Francky Fouedjio Mines ParisTech Centre de Géosciences, Equipe Géostatistique
More informationOn Sampling Errors in Empirical Orthogonal Functions
3704 J O U R N A L O F C L I M A T E VOLUME 18 On Sampling Errors in Empirical Orthogonal Functions ROBERTA QUADRELLI, CHRISTOPHER S. BRETHERTON, AND JOHN M. WALLACE University of Washington, Seattle,
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota,
More informationClimate Change and Global Warming
Climate Change and Global Warming Michael E. Mann Department of Environmental Sciences University of Virginia ETP 401 Guest Lecture University of Virginia April 12, 2005 `The balance of evidence suggests
More informationPaper Review: NONSTATIONARY COVARIANCE MODELS FOR GLOBAL DATA
Paper Review: NONSTATIONARY COVARIANCE MODELS FOR GLOBAL DATA BY MIKYOUNG JUN AND MICHAEL L. STEIN Presented by Sungkyu Jung April, 2009 Outline 1 Introduction 2 Covariance Models 3 Application: Level
More informationROBUST MEASUREMENT OF THE DURATION OF
ROBUST MEASUREMENT OF THE DURATION OF THE GLOBAL WARMING HIATUS Ross McKitrick Department of Economics University of Guelph Revised version, July 3, 2014 Abstract: The IPCC has drawn attention to an apparent
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 informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Alan Gelfand 1 and Andrew O. Finley 2 1 Department of Statistical Science, Duke University, Durham, North
More informationParameter Estimation in the Spatio-Temporal Mixed Effects Model Analysis of Massive Spatio-Temporal Data Sets
Parameter Estimation in the Spatio-Temporal Mixed Effects Model Analysis of Massive Spatio-Temporal Data Sets Matthias Katzfuß Advisor: Dr. Noel Cressie Department of Statistics The Ohio State University
More informationImperfect Data in an Uncertain World
Imperfect Data in an Uncertain World James B. Elsner Department of Geography, Florida State University Tallahassee, Florida Corresponding author address: Dept. of Geography, Florida State University Tallahassee,
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Andrew O. Finley Department of Forestry & Department of Geography, Michigan State University, Lansing
More informationA spatio-temporal model for extreme precipitation simulated by a climate model
A spatio-temporal model for extreme precipitation simulated by a climate model Jonathan Jalbert Postdoctoral fellow at McGill University, Montréal Anne-Catherine Favre, Claude Bélisle and Jean-François
More informationProxy reconstructions of Pacific decadal variability. March 8, 2006
Proxy reconstructions of Pacific decadal variability Michael N. Evans Laboratory of Tree-Ring Research/Geosciences/Atmospheric Sciences ATMO/GEOS 513 web pages: http://ic.ltrr.arizona.edu/ic/enso/ March
More informationMulti-resolution models for large data sets
Multi-resolution models for large data sets Douglas Nychka, National Center for Atmospheric Research National Science Foundation NORDSTAT, Umeå, June, 2012 Credits Steve Sain, NCAR Tia LeRud, UC Davis
More informationMultiple timescale coupled atmosphere-ocean data assimilation
Multiple timescale coupled atmosphere-ocean data assimilation (for climate prediction & reanalysis) Robert Tardif Gregory J. Hakim H L H University of Washington w/ contributions from: Chris Snyder NCAR
More informationA New Bayesian Variable Selection Method: The Bayesian Lasso with Pseudo Variables
A New Bayesian Variable Selection Method: The Bayesian Lasso with Pseudo Variables Qi Tang (Joint work with Kam-Wah Tsui and Sijian Wang) Department of Statistics University of Wisconsin-Madison Feb. 8,
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 informationSpatial Statistics for Energy Challenges
Spatial Statistics for Energy Challenges Tilmann Gneiting Heidelberg Institute for Theoretical Studies (HITS) Karlsruhe Institute of Technology (KIT) tilmann.gneiting@h-its.org http://www.h-its.org/english/research/cst/
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Andrew O. Finley 1 and Sudipto Banerjee 2 1 Department of Forestry & Department of Geography, Michigan
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 informationConstraining simulated atmospheric states by sparse empirical information
Constraining simulated atmospheric states by sparse empirical information Martin Widmann University of Birmingham Hugues Goosse University of Louvain-la-Neuve NCAR workshop Data assimilation and climate
More informationHierarchical Modelling for Multivariate Spatial Data
Hierarchical Modelling for Multivariate Spatial Data Geography 890, Hierarchical Bayesian Models for Environmental Spatial Data Analysis February 15, 2011 1 Point-referenced spatial data often come as
More informationThe Matrix Reloaded: Computations for large spatial data sets
The Matrix Reloaded: Computations for large spatial data sets Doug Nychka National Center for Atmospheric Research The spatial model Solving linear systems Matrix multiplication Creating sparsity Sparsity,
More informationCBMS Lecture 1. Alan E. Gelfand Duke University
CBMS Lecture 1 Alan E. Gelfand Duke University Introduction to spatial data and models Researchers in diverse areas such as climatology, ecology, environmental exposure, public health, and real estate
More informationFusing point and areal level space-time data. data with application to wet deposition
Fusing point and areal level space-time data with application to wet deposition Alan Gelfand Duke University Joint work with Sujit Sahu and David Holland Chemical Deposition Combustion of fossil fuel produces
More informationPossible causes of data model discrepancy in the temperature history of the last Millennium
This is a post-peer-review, pre-copyedit version of an article published in Scientific Reports. The final authenticated version is available online at: http://dx.doi.org/10.1038/s41598-018-25862-2 Possible
More informationthe long tau-path for detecting monotone association in an unspecified subpopulation
the long tau-path for detecting monotone association in an unspecified subpopulation Joe Verducci Current Challenges in Statistical Learning Workshop Banff International Research Station Tuesday, December
More informationWrapped Gaussian processes: a short review and some new results
Wrapped Gaussian processes: a short review and some new results Giovanna Jona Lasinio 1, Gianluca Mastrantonio 2 and Alan Gelfand 3 1-Università Sapienza di Roma 2- Università RomaTRE 3- Duke University
More informationRecent Climate History - The Instrumental Era.
2002 Recent Climate History - The Instrumental Era. Figure 1. Reconstructed surface temperature record. Strong warming in the first and late part of the century. El Ninos and major volcanic eruptions are
More informationTests for separability in nonparametric covariance operators of random surfaces
Tests for separability in nonparametric covariance operators of random surfaces Shahin Tavakoli (joint with John Aston and Davide Pigoli) April 19, 2016 Analysis of Multidimensional Functional Data Shahin
More informationModelling Ireland s exchange rates: from EMS to EMU
From the SelectedWorks of Derek Bond November, 2007 Modelling Ireland s exchange rates: from EMS to EMU Derek Bond, University of Ulster Available at: https://works.bepress.com/derek_bond/15/ Background
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 informationPart 5: Spatial-Temporal Modeling
Part 5: Spatial-Temporal Modeling 1 Introductory Model Formulation For a broad spectrum of applications within the realm of spatial statistics, commonly the data vary across both space and time giving
More informationModels for spatial data (cont d) Types of spatial data. Types of spatial data (cont d) Hierarchical models for spatial data
Hierarchical models for spatial data Based on the book by Banerjee, Carlin and Gelfand Hierarchical Modeling and Analysis for Spatial Data, 2004. We focus on Chapters 1, 2 and 5. Geo-referenced data arise
More informationThe Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations
The Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations John R. Michael, Significance, Inc. and William R. Schucany, Southern Methodist University The mixture
More informationBayesian SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling
Bayesian SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling Jon Wakefield Departments of Statistics and Biostatistics University of Washington 1 / 37 Lecture Content Motivation
More informationExtratropical transition of North Atlantic tropical cyclones in variable-resolution CAM5
Extratropical transition of North Atlantic tropical cyclones in variable-resolution CAM5 Diana Thatcher, Christiane Jablonowski University of Michigan Colin Zarzycki National Center for Atmospheric Research
More informationClimate Change Insights from the. Centuries
Climate Change Insights from the Paleoclimate Record of Past Centuries Michael E. Mann Penn State University with contributions from: Caspar Ammann, Raymond Bradley, Elizabeth Crespin, Jeff Donnelly, Greg
More informationModeling daily precipitation in Space and Time
Space and Time SWGen - Hydro Berlin 20 September 2017 temporal - dependence Outline temporal - dependence temporal - dependence Stochastic Weather Generator Stochastic Weather Generator (SWG) is a stochastic
More informationClimate variations in the eastern Mediterranean over the last millennium; New evidence from proxy reconstructions and high resolution climate models
Climate variations in the eastern Mediterranean over the last millennium; New evidence from proxy reconstructions and high resolution climate models Elena Xoplaki, Jürg Luterbacher Juan José Gómez-Navarro,
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationBayesian hierarchical modelling for data assimilation of past observations and numerical model forecasts
Bayesian hierarchical modelling for data assimilation of past observations and numerical model forecasts Stan Yip Exeter Climate Systems, University of Exeter c.y.yip@ex.ac.uk Joint work with Sujit Sahu
More informationStatistical analysis of a spatio-temporal model with location dependent parameters and a test for spatial stationarity
Statistical analysis of a spatio-temporal model with location dependent parameters and a test for spatial stationarity Suhasini Subba Rao Department of Statistics, Texas A&M University, College Station,
More informationBayesian inference for multivariate extreme value distributions
Bayesian inference for multivariate extreme value distributions Sebastian Engelke Clément Dombry, Marco Oesting Toronto, Fields Institute, May 4th, 2016 Main motivation For a parametric model Z F θ of
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 informationSTATS 200: Introduction to Statistical Inference. Lecture 29: Course review
STATS 200: Introduction to Statistical Inference Lecture 29: Course review Course review We started in Lecture 1 with a fundamental assumption: Data is a realization of a random process. The goal throughout
More informationSemi-parametric estimation of non-stationary Pickands functions
Semi-parametric estimation of non-stationary Pickands functions Linda Mhalla 1 Joint work with: Valérie Chavez-Demoulin 2 and Philippe Naveau 3 1 Geneva School of Economics and Management, University of
More informationEstimation of non-stationary spatial covariance structure
Estimation of non-stationary spatial covariance structure DAVID J NOTT Department of Statistics, University of New South Wales, Sydney 2052, Australia djn@mathsunsweduau WILLIAM T M DUNSMUIR Division of
More informationStatistical Inference
Statistical Inference Liu Yang Florida State University October 27, 2016 Liu Yang, Libo Wang (Florida State University) Statistical Inference October 27, 2016 1 / 27 Outline The Bayesian Lasso Trevor Park
More information2. Outline of the MRI-EPS
2. Outline of the MRI-EPS The MRI-EPS includes BGM cycle system running on the MRI supercomputer system, which is developed by using the operational one-month forecasting system by the Climate Prediction
More informationFundamentals of Data Assimila1on
014 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO July 14-16, 014 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
More informationLocal Likelihood Bayesian Cluster Modeling for small area health data. Andrew Lawson Arnold School of Public Health University of South Carolina
Local Likelihood Bayesian Cluster Modeling for small area health data Andrew Lawson Arnold School of Public Health University of South Carolina Local Likelihood Bayesian Cluster Modelling for Small Area
More informationSimple example of analysis on spatial-temporal data set
Simple example of analysis on spatial-temporal data set I used the ground level ozone data in North Carolina (from Suhasini Subba Rao s website) The original data consists of 920 days of data over 72 locations
More informationSpatial Statistics with Image Analysis. Outline. A Statistical Approach. Johan Lindström 1. Lund October 6, 2016
Spatial Statistics Spatial Examples More Spatial Statistics with Image Analysis Johan Lindström 1 1 Mathematical Statistics Centre for Mathematical Sciences Lund University Lund October 6, 2016 Johan Lindström
More informationSecond-Order Analysis of Spatial Point Processes
Title Second-Order Analysis of Spatial Point Process Tonglin Zhang Outline Outline Spatial Point Processes Intensity Functions Mean and Variance Pair Correlation Functions Stationarity K-functions Some
More informationDETECTION AND ATTRIBUTION
DETECTION AND ATTRIBUTION Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260, USA rls@email.unc.edu 41st Winter Conference in Statistics
More informationImpact of METOP ASCAT Ocean Surface Winds in the NCEP GDAS/GFS and NRL NAVDAS
Impact of METOP ASCAT Ocean Surface Winds in the NCEP GDAS/GFS and NRL NAVDAS COAMPS @ Li Bi 1,2 James Jung 3,4 Michael Morgan 5 John F. Le Marshall 6 Nancy Baker 2 Dave Santek 3 1 University Corporation
More informationDetection and Attribution of Climate Change. ... in Indices of Extremes?
Detection and Attribution of Climate Change... in Indices of Extremes? Reiner Schnur Max Planck Institute for Meteorology MPI-M Workshop Climate Change Scenarios and Their Use for Impact Studies September
More informationToday s Climate in Perspective: Hendrick Avercamp ( ) ~1608; Rijksmuseum, Amsterdam
Today s Climate in Perspective: Paleoclimate Evidence Hendrick Avercamp (1585-1634) ~1608; Rijksmuseum, Amsterdam Observations Instrumental surface temperature records? (Le Treut et al., 2007 IPCC AR4
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 informationNOTES AND CORRESPONDENCE
5666 J O U R N A L O F C L I M A T E VOLUME 20 NOTES AND CORRESPONDENCE Comments on Testing the Fidelity of Methods Used in Proxy-Based Reconstructions of Past Climate : The Role of the Standardization
More informationCopulas, a novel approach to model spatial and spatio-temporal dependence
Copulas, a novel approach to model spatial and spatio-temporal dependence Benedikt Gräler 1, Hannes Kazianka 2, Giovana Mira de Espindola 3 1 Institute for Geoinformatics, University of Münster, Germany
More informationSTATISTICAL INTERPOLATION METHODS RICHARD SMITH AND NOEL CRESSIE Statistical methods of interpolation are all based on assuming that the process
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 STATISTICAL INTERPOLATION METHODS
More informationStatistical Sciences Symposium, 2014
College of Letters and Science Department of Statistics UC Davis Statistical Sciences Symposium, 2014 Spatial-Temporal Statistics: Methods and Applications SATURDAY April 26 8:30 a.m. 5:15 p.m. Mathematical
More informationAn integrated speleothem proxy and climate modeling study of the last deglacial climate in the Pacific Northwest
An integrated speleothem proxy and climate modeling study of the last deglacial climate in the Pacific Northwest Jerry Potter UC Davis, University of Michigan, NASA GSFC Isabel Montanez, University of
More informationExtending the use of surface-sensitive microwave channels in the ECMWF system
Extending the use of surface-sensitive microwave channels in the ECMWF system Enza Di Tomaso and Niels Bormann European Centre for Medium-range Weather Forecasts Shinfield Park, Reading, RG2 9AX, United
More informationAsymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands
Asymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands Elizabeth C. Mannshardt-Shamseldin Advisor: Richard L. Smith Duke University Department
More information