510 likes | 669 Vues
Tutorial: Analysis Tools for GRACE time-variable harmonic coefficients Jürgen Kusche , Annette Eicker, Ehsan Forootan Bonn University, Germany IGCP 565 Workshop November 21-22, 2011 University of the Witwatersrand, Johannesburg, South Africa. Analysis Tools.
E N D
Tutorial: Analysis Tools for GRACE time-variable harmonic coefficientsJürgen Kusche, Annette Eicker, Ehsan ForootanBonn University, GermanyIGCP 565 WorkshopNovember 21-22, 2011University of the Witwatersrand, Johannesburg, South Africa
Analysis Tools For everybody who has worked with GRACE data, it has become clear that GRACE solutions (say SH coefficients converted to TWS, total water storage) require some post-processing, • to suppress correlated noise, remove ‚stripes‘ ( filtering) • to extract the dominating ‚modes‘ of temporal variability ( PCA, …) Being in general use, these analysis tools always remove signal content together with ‚noise‘. For any comparison of GRACE data with geophysical modelling, it is imperative therefore that the same tool is applied to both. For getting ‚absolute‘ amplitudes, rates, etc., it is imperative to consider the ‚bias‘ of an analysis technique. Note: If you download GRACE gridded products from GRACE Tellus website, GFZ ICGEM or others, some filtering has been applied already. Try to understand what has been done to the data, and what has to be done to be fully compatible with model output.
Analysis Chain: Filtering GRACE SHC SHA geophysical model SHC Globally defined? Treatment degree 1 & 2, temporal anomalies wrt epoch t, restore AOD Temporal alignement, possibly remove deg. 1 consistent to GRACE Filtering (in spectral domain) Filtering (in spectral domain) Multiplication in spectral domain Mapping to space domain or basin averaging Mapping to space domain or basin averaging OR Filtering (in space domain) OR Filtering (in space domain) Convolution on the sphere
Analysis Chain: Filtering - quite often GRACE SHC Geophysical model Treatment degree 1 & 2, temporal anomalies wrt epoch t, restore AOD Temporal alignement, possibly remove deg. 1 consistent to GRACE Filtering (in spectral domain) Filtering (in spectral domain) Mapping to space domain or basin averaging Mapping to space domain or basin averaging Filtering (in space domain)
Analysis Tools: PCA Filtering (in spectral domain) Filtering (in spectral domain) Mapping to space domain or basin averaging Mapping to space domain or basin averaging Computation of EOFs from GRACE or model Project gridded data onto EOFs: PCs for GRACE and model Comparison Truncated (filtered) reconstruction
Analysis Chain: PCA - also possible… GRACE SHC SHA geophysical model SHC Treatment degree 1 & 2, temporal anomalies wrt epoch t, restore AOD Temporal alignement, possibly remove deg. 1 consistent to GRACE PCA & truncated reconstruction (in spectral domain) Mapping to space domain or basin averaging
Part I: Filtering techniques and their application to GRACE data What is the purpose of filtering a GRACE solution? Unfiltered GRACE solution
Part I: Filtering techniques and their application to GRACE data Filtering • Attempts to suppress ‚noise‘ in data (here: SH coefficients) • Requires that we have an a-priori knowledge of expected ‚noise‘ ( characterize spectral behaviour of noise) • Filters take this into account either implicitly (‚deterministic filters‘ or explicitly ‚stochastic filters‘) • Also regularizing or ‚constraining‘ GRACE normal equations corresponds to some kind of filtering (Kusche 2007, Klees et al 2008, Swenson & Wahr 2011) Data sets e.g. • GRACE-derived SH coefficients or maps of geoid, gravity anomalies,TWS • SH coefficients or maps of TWS and other data from model output
Part I: Filtering techniques and their application to GRACE data What is the purpose of filtering a GRACE solution? Boxcar-filtered GRACE solution: Remove stripes by averaging, convolution
Part I: Filtering techniques and their application to GRACE data What is the purpose of filtering a GRACE solution? Boxcar-filtered GRACE solution: Remove stripes by averaging, convolution
Part I: Filtering techniques and their application to GRACE data What is the purpose of filtering a GRACE solution? Boxcar-filtered GRACE solution: Remove stripes by averaging, convolution
Part I: Filtering techniques and their application to GRACE data What is the purpose of filtering a GRACE solution? Gaussian-filtered GRACE solution: convolution with a smooth kernel
Part I: Filtering techniques and their application to GRACE data Filtering a GRACE solution in spatial and in spectral domain
Part I: Filtering techniques and their application to GRACE data Notations I use • Positive and negative order • Fully normalized • All factors are put into the coefficients Potential Surface Mass (involves Love numbers)
Part I: Filtering techniques and their application to GRACE data Isotropic filters
Part I: Filtering techniques and their application to GRACE data Isotropic filters: Gaussian filter harmonic degree n
Part I: Filtering techniques and their application to GRACE data Destriping filters Isotropic Non-isotropic
Part I: Filtering techniques and their application to GRACE data Destriping filters Swenson & Wahr 2006
Part I: Filtering techniques and their application to GRACE data Destriping filters Kusche 2007
Part I: Filtering techniques and their application to GRACE data Destriping filters and WRMS reduction GRACE TWS WRMS [cm], 3 filters stronger filtering less
Part I: Filtering techniques and their application to GRACE data Basin averaging … Ocean mass t t
Part I: Filtering techniques and their application to GRACE data Basin averaging (spatial domain) (spectral domain)
Part I: Filtering techniques and their application to GRACE data Basin averaging Note: Equality requires the integral is discretized at an error small enough the SH truncation degree is big enough OR BOTH F and O are band-limited (spectral domain)
Part I: Filtering techniques and their application to GRACE data Smoothed basin averaging and bias
Part I: Filtering techniques and their application to GRACE data Leakage problem and filter bias: basin averaging, what happens? GRACE (truncated SH) Model
Part I: Filtering techniques and their application to GRACE data Leakage problem and filter bias: basin averaging, what happens? GRACE (truncated SH) Model + noise
Part I: Filtering techniques and their application to GRACE data Leakage problem and filter bias: basin averaging, what happens? GRACE (truncated SH) Model Exact averaging Spectral leakage Lost signal
Part I: Filtering techniques and their application to GRACE data Leakage problem and filter bias: basin averaging, what happens? GRACE (truncated SH) Model Lost signal (leaking out) Leaking out Bias can be computed from model
Part I: Filtering techniques and their application to GRACE data Leakage problem and filter bias: basin averaging, what happens? Lost signal (leaking out) Signal added by surrounding region (leaking in)
Part II: Principle component analysis and related ideas Empirical orthogonal function analysis Principle components Related concepts
Part II: Principle component analysis and related ideas PCA • attempts to find a relatively small number of independent modes in a data set that convey as much as possible information without redundancy • can be used to explore the structure of the data variability in an objective way, i.e. without assumptions on periodic behaviour etc. • and EOF analysis are the same. Data sets e.g. • GRACE-derived maps of TWS, TWS and other maps from model output • Sea level anomalies • Other related spatial fields (SST, SLP, …)
Part II: Principle component analysis and related ideas What does PCA do? PCA uses a set of orthogonal functions (EOFs) to represent a spatio-temporal data field in the following way EOF ej =e (, ) show spatial patterns of the major factors (‚modes‘) that account for temporal variations. PC dj;i = d(t) tells how the amplitude of EOF varies with time. epochs (time) PCs (expansion coefficients) grid points EOFs (new basis)
Part II: Principle component analysis and related ideas How do we obtain the PCs? Lets assume the EOFs are orthogonal (why are they called EOF, after all) and normalized. The PCs are found as an orthogonal projection of the data onto the new basis functions (the EOFs). We can try to reconstruct the original data using only the ‚major‘ EOFs
Part II: Principle component analysis and related ideas PC1 What do we get from PCA? 62% t EOF1 PC2 24% t EOF2 … PCn <1% t t EOFn
Part II: Principle component analysis and related ideas Example (I) GRACE TWS 400km Gaussian filtered (note: EOF PC has unit of data)
Part II: Principle component analysis and related ideas How do we choose the EOFs? Total variance Then, the maximum variance is concentrated in a single EOF if I.e., maximize subject to Or, solve an eigenvalue problem max. variance
Part II: Principle component analysis and related ideas How do we choose the EOFs? The EOFs are found as the eigenvectors of the data covariance matrix. Some remarks (I) • The covariance as above is temporal, i.e. it considers auto- and covariances of time series per grid point (for n grid points) (we‘ll come back to this point)
Part II: Principle component analysis and related ideas Representing data in an eigenvector basis
Part II: Principle component analysis and related ideas The EOFs are found as the eigenvectors of the data covariance matrix. Some remarks (II) • It is an empirical realization. Should we know the true covariance, we might better use this one instead of the empirical one • All data has been considered as (perfectly) centered • Eigenvectors require normalization and a further convention for uniqueness, e.g.
Part II: Principle component analysis and related ideas The EOFs are found as the eigenvectors of the data covariance matrix. Some remarks (III) • Alternatively, we could use SVD applied to the data matrix • Each EOF explains a fraction of the total variance, given by the ratio of the EV vs. TV (total variance) • It is common to choose the number of EOFs as to ‚explain‘’ 90% (or …) of the TV
Part II: Principle component analysis and related ideas Example GRACE global analysis (left: EV, right: cumulative percentage of TV) (courtesy E. Forootan)
Part II: Principle component analysis and related ideas The EOFs are found as the eigenvectors of the data covariance matrix. Some remarks (IV) • Instead of computing the temporal data covariance, we may compute the spatial covariance (spatial variance and covariances for the p epochs) • Requires less memory for p < n • Temporal and spatial covariance matrices share the same p EVs • k-th EOF (spatial) ~ k-th PC (temporal) and vice versa
Part II: Principle component analysis and related ideas Some remarks (V) • Linear transformations of the data lead to new eigenvectors and –values Therefore: • EOFs and PCs computed on a regional grid look different from EOFs/PCs computed on a global grid • GRACE: EOFs of geoid change look different from EOF of TWS change • PCA applied to GRACE SH coefficients (EOF filter of Schrama & Wouters) looks different from PCA applied to grids.
Part II: Principle component analysis and related ideas Example (I) revisited Trend (decrease)+ annual Semi-annual (modulation of annual) Annual signals out of phase Trend (increase)+ annual Trend ? 2007: unusually strong rainfall in Congo
Part II: Principle component analysis and related ideas Examples (II) Wouters & Schrama (2007) Direct EOF filtering of GRACE SH coefficients Top: Unfiltered/Gaussian, Middle: EOF filtered, Bottom: difference. Unit [cm]
Part II: Principle component analysis and related ideas How many modes (EOFs) should we retain? In other words, how many % of the data TV should we reconstruct? North et al. 1982, Month. Weath. Rev.: Consider the spatio-temporal data as stochastic, i.e. perturbed by e.g. Gaussian noise. Then, the covariance and the eigenvalues / -vectors will be stochastic as well. In first order… If the sampling error in the eigenvalue is comparable to the spacing of the eigenvalues, then the sampling error of the EOF will be comparable to the nearby EOF. And then it is time to truncate.
Part II: Principle component analysis and related ideas North et al (1982), reprinted in Von Storch & Zwiers
Part II: Principle component analysis and related ideas Related concepts: (Orthogonal) EOF rotation REOFs are still orthogonal - RPCs are correlated now. How can we use this degree of freedom?
Part II: Principle component analysis and related ideas Related concepts: REOF, how can we use this degree of freedom? I.e. how do we find the rotation matrix? E.g. VARIMAX, maximize the variance of the square of the REOFs (i.e. the spreading of the total variability of the modes) E.g. ICA, minimize 3th / 4th statistical moments of the REOFs to maximize the independence of the RPCs
Part II: Principle component analysis and related ideas simulated PCA VARIMAX Dommenget & Latif (2001)