On the computation of mass-change trends from GRACE gravity field time-series

Similar documents
QR factorization. Let P 1, P 2, P n-1, be matrices such that Pn 1Pn 2... PPA

Midterm Exam. Tuesday, September hour, 15 minutes

Survival Prediction Based on Compound Covariate under Cox Proportional Hazard Models

Least Squares Fitting (LSQF) with a complicated function Theexampleswehavelookedatsofarhavebeenlinearintheparameters

Determination of Antoine Equation Parameters. December 4, 2012 PreFEED Corporation Yoshio Kumagae. Introduction

Cyclone. Anti-cyclone

The Linear Regression Of Weighted Segments

As evident from the full-sample-model, we continue to assume that individual errors are identically and

Efficient Estimators for Population Variance using Auxiliary Information

Chapter 8. Simple Linear Regression

The Bernstein Operational Matrix of Integration

IMPROVED PORTFOLIO OPTIMIZATION MODEL WITH TRANSACTION COST AND MINIMAL TRANSACTION LOTS

Other Topics in Kernel Method Statistical Inference with Reproducing Kernel Hilbert Space

Spike-and-Slab Dirichlet Process Mixture Models

Pricing Asian Options with Fourier Convolution

Comparison of the Bayesian and Maximum Likelihood Estimation for Weibull Distribution

14. Poisson Processes

The Optimal Combination Forecasting Based on ARIMA,VAR and SSM

AML710 CAD LECTURE 12 CUBIC SPLINE CURVES. Cubic Splines Matrix formulation Normalised cubic splines Alternate end conditions Parabolic blending

Linear Regression Linear Regression with Shrinkage

VARIATIONAL ITERATION METHOD FOR DELAY DIFFERENTIAL-ALGEBRAIC EQUATIONS. Hunan , China,

Use of Non-Conventional Measures of Dispersion for Improved Estimation of Population Mean

Least squares and motion. Nuno Vasconcelos ECE Department, UCSD

Key words: Fractional difference equation, oscillatory solutions,

AN INCREMENTAL QUASI-NEWTON METHOD WITH A LOCAL SUPERLINEAR CONVERGENCE RATE. Aryan Mokhtari Mark Eisen Alejandro Ribeiro

COMPARISON OF ESTIMATORS OF PARAMETERS FOR THE RAYLEIGH DISTRIBUTION

Final Exam Applied Econometrics

Real-Time Systems. Example: scheduling using EDF. Feasibility analysis for EDF. Example: scheduling using EDF

SYRIAN SEISMIC CODE :

Nature and Science, 5(1), 2007, Han and Xu, Multi-variable Grey Model based on Genetic Algorithm and its Application in Urban Water Consumption

International Journal Of Engineering And Computer Science ISSN: Volume 5 Issue 12 Dec. 2016, Page No.

θ = θ Π Π Parametric counting process models θ θ θ Log-likelihood: Consider counting processes: Score functions:

The textbook expresses the stock price as the present discounted value of the dividend paid and the price of the stock next period.

The textbook expresses the stock price as the present discounted value of the dividend paid and the price of the stock next period.

The Signal, Variable System, and Transformation: A Personal Perspective

Fault Tolerant Computing. Fault Tolerant Computing CS 530 Probabilistic methods: overview

To Estimate or to Predict

APPLICATION OF ADAPTIVE NEURO-FUZZY INFERENCE SYSTEM IN INTEREST RATES EFFECTS ON STOCK RETURNS

Regression Approach to Parameter Estimation of an Exponential Software Reliability Model

Quantitative Portfolio Theory & Performance Analysis

Lecture 3 Topic 2: Distributions, hypothesis testing, and sample size determination

Solution. The straightforward approach is surprisingly difficult because one has to be careful about the limits.

Fully Fuzzy Linear Systems Solving Using MOLP

An Efficient Dual to Ratio and Product Estimator of Population Variance in Sample Surveys

Real-time Classification of Large Data Sets using Binary Knapsack

Volatility Change Point Detection Using Stochastic Differential Equations and Time Series Control Charts

(1) Cov(, ) E[( E( ))( E( ))]

Decompression diagram sampler_src (source files and makefiles) bin (binary files) --- sh (sample shells) --- input (sample input files)

Chapter 3: Maximum-Likelihood & Bayesian Parameter Estimation (part 1)

The ray paths and travel times for multiple layers can be computed using ray-tracing, as demonstrated in Lab 3.

RATIO ESTIMATORS USING CHARACTERISTICS OF POISSON DISTRIBUTION WITH APPLICATION TO EARTHQUAKE DATA

Correlation and Regression Analysis

EE 6885 Statistical Pattern Recognition

Moments of Order Statistics from Nonidentically Distributed Three Parameters Beta typei and Erlang Truncated Exponential Variables

Brownian Motion and Stochastic Calculus. Brownian Motion and Stochastic Calculus

BACKTESTING VAR ESTIMATION UNDER GARCH AND GJR-GARCH MODELS

Moving Least Squares Coordinates. Josiah Manson and Scott Schaefer Texas A&M University

Lecture Notes Forecasting the process of estimating or predicting unknown situations

Computational Functional Anatomy

ESTIMATION AND TESTING

ASYMPTOTIC EQUIVALENCE OF NONPARAMETRIC REGRESSION AND WHITE NOISE. BY LAWRENCE D. BROWN 1 AND MARK G. LOW 2 University of Pennsylvania

Continuous Time Markov Chains

Econ107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)

Research on portfolio model based on information entropy theory

Iterated Bernstein polynomial approximations

Towards combined global monthly gravity field solutions

Supplement Material for Inverse Probability Weighted Estimation of Local Average Treatment Effects: A Higher Order MSE Expansion

NUMERICAL EVALUATION of DYNAMIC RESPONSE

D 2 : Decentralized Training over Decentralized Data

GIA Requirements for Future Gravity Satellite Missions: Implications After More Than 7 Years of GRACE Data

Stability Criterion for BAM Neural Networks of Neutral- Type with Interval Time-Varying Delays

Solution set Stat 471/Spring 06. Homework 2

Multiple Choice Test. Chapter Adequacy of Models for Regression

Fundamentals of Speech Recognition Suggested Project The Hidden Markov Model

Optimized Calculation of Hourly Price Forward Curve (HPFC)

Solving Fuzzy Equations Using Neural Nets with a New Learning Algorithm

The Poisson Process Properties of the Poisson Process

Study of Correlation using Bayes Approach under bivariate Distributions

Modeling of the linear time-variant channel. Sven-Gustav Häggman

Data Collection Definitions of Variables - Conceptualize vs Operationalize Sample Selection Criteria Source of Data Consistency of Data

Mathematical Formulation

Study on one-dimensional consolidation of soil under cyclic loading and with varied compressibility *

Stochastic Petri Nets with Low Variation Matrix Exponentially Distributed Firing Time

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

Practice Final Exam (corrected formulas, 12/10 11AM)

Probability Bracket Notation and Probability Modeling. Xing M. Wang Sherman Visual Lab, Sunnyvale, CA 94087, USA. Abstract

Neural Network Global Sliding Mode PID Control for Robot Manipulators

EDUCATION COMMITTEE OF THE SOCIETY OF ACTUARIES ADVANCED TOPICS IN GENERAL INSURANCE STUDY NOTE CREDIBILITY WITH SHIFTING RISK PARAMETERS

Spatial-Temporal Separation Based on the Dynamic Recurrent Wavelet Neural Network Modelling for ASP Flooding

Journal of Mathematical Psychology

EMD Based on Independent Component Analysis and Its Application in Machinery Fault Diagnosis

Comparison of Dual to Ratio-Cum-Product Estimators of Population Mean

Expert-Statistical Processing of Data and the Method of Analogs in Solution of Applied Problems in Control Theory

FALL HOMEWORK NO. 6 - SOLUTION Problem 1.: Use the Storage-Indication Method to route the Input hydrograph tabulated below.

Assessing Normality. Assessing Normality. Assessing Normality. Assessing Normality. Normal Probability Plot for Normal Distribution.

Chapter 13 Student Lecture Notes 13-1

Big O Notation for Time Complexity of Algorithms

Research on SVM Prediction Model Based on Chaos Theory

A Novel Algorithm for Criminal Statistical Processing

Reliability Analysis of Sparsely Connected Consecutive-k Systems: GERT Approach

Transcription:

Geodäsche Woche 0 Nürberg, 7. 9. Sepember 0 O he compuao of mass-chage reds from GRACE gravy feld me-seres Olver Baur Isu für Welraumforschug, Öserrechsche Aademe der Wsseschafe

Movao Greelad lear mass-chage reds (jus a few examples... Velcoga ad Wahr (GRL, 005 Ramlle e al. (GPC, 006 Averagg 04/00 07/004 CSR erel 07/00 03/005 Averagg CNES/GRGS erel Che e al. (Scece, 006 Luhce e al. (Scece, 006 07/003 07/005 Mascos KBRR 04/00 /005 Forward CSR modelg Wouers e al. (GRL, 008 Baur e al. (JGR, 009 Velcoga (GRL, 009 Che e al. (JGR, 0 0/003 0/008 Forward CSR, modelg GFZ Leaage 08/00 07/008 CSR, quafcao GFZ, JPL 04/00 0/009 Averagg CSR erel 04/00 /009 Forward CSR modelg -50-00 -50-00 -50-300 G/yr Sources of dscrepacy mpac geoceer of GIA geoceer moo sgal from moo SLR Iverso mehod / pre-processg Perod of vesgao Daa used Daa mapulao (c 0, degree-

Movao Greelad mass-chage accelerao (aga, jus a few examples... Velcoga ad Wahr (Naure, 006 04/00 04/004: -95 ± 50 G/yr 05/004 04/006: -34 ± 6 G/yr Che e al. (JGR, 0 04/00 03/005: -44 ± 5 G/yr 04/005 /009: -48 ± 43 G/yr Velcoga (GRL, 009 Schrama e al. (JGR, 0 Rgo e al. (GRL, 0-5 -0-5 -30-35 -40 G/yr Sources of dscrepacy Iverso mehod / pre-processg Perod of vesgao Daa used Daa mapulao (c 0, degree- Compuao of mass-chage reds 3

Coe Lear red modelg No-lear red modelg Model seleco Coclusos Greelad Iverso mehod: Baur e al. (JGR, 009 Perod of vesgao: 05/00 04/0 (ad subses hereof Daa used: CSR (RL04 Daa mapulao : c 0 : SLR, geoceer moo: SLR, GIA: Paulso e al. (007 Compuao of mass-chage reds: varable Amazo bas 4

Lear red modelg Model fuco Tred, approxmaed by frs-order polyomal (regresso le y ( 0,,..., Addoal aual sgal y ( A cos( f B s( f 0,,..., F chage rae (G/yr Regresso le -3 ± 5 + aual sgal -3 ± 4 + sem-aual sgal -3 ± 3 + S dal alas -3 ± 3 + K dal alas -30 ± 3 + K dal alas -30 ± 3 5

Lear red modelg Model fuco Tred, approxmaed by frs-order polyomal (regresso le y ( 0,,..., Addoal aual ad sem-aual sgal ad dal alases y ( cos( s(,,..., ; 0 A f B f j j j j j j,... F chage rae (G/yr Regresso le -3 ± 5 + aual sgal -3 ± 4 + sem-aual sgal -3 ± 3 + S dal alas -3 ± 3 + K dal alas -30 ± 3 + K dal alas -30 ± 3 6

Lear red modelg Model fuco Tred, approxmaed by frs-order polyomal (regresso le y ( 0,,..., Addoal aual ad sem-aual sgal ad dal alases y ( cos( s(,,..., ; 0 A f B f j j j j j j,... F chage rae (G/yr Regresso le -3 ± 5 + aual sgal -3 ± 4 + sem-aual sgal -3 ± 3 + S dal alas -3 ± 3 K + K dal alas -30 ± 3 + K dal alas -30 ± 3 K S 7

Lear red modelg Smoohg of me seres Tred, approxmaed by frs-order polyomal (regresso le Addoal aual ad sem-aual sgal ad dal alases Ier-aual varaos 3-moh wdow Lowess fler low-pass fler lowess 0.5-33 ± G/yr lowess 0. -33 ± G/yr Velcoga (GRL, 009 8

Lear red modelg F chage rae (G/yr Regresso le 50 ± 8 + aual sgal 43 ± 6 + sem-aual sgal 43 ± 6 + S dal alas 43 ± 6 + K dal alas 36 ± 7 + K dal alas 36 ± 7 lowess 0.5 0. 39 4 ± 3 G/yr K K S 9

Lear red modelg Mass-chage accelerao Che e al. (JGR, 0 04/00 03/005: -9 ± 38 G/yr 04/00 03/005: -44 ± 5 G/yr 04/005 /009: -48 ± 43 G/yr Re-compuao - ± 4 G/yr -70 ± 3 G/yr -50 ± 8 G/yr 0

Lear red modelg Mass-chage accelerao Che e al. (JGR, 0 04/00 03/005: -44 ± 5 G/yr 04/005 /009: -48 ± 43 G/yr Velcoga ad Wahr (Naure, 006 04/00 04/004: -95 ± 50 G/yr 05/004 04/006: -34 ± 6 G/yr Velcoga (GRL, 009 04/00??/003: -37 G/yr??/007 0/009: -86 G/yr Choce of accelerao po?

Lear red modelg Subse perod: 5 years Subse shf: 6 mohs Subse perod: 6 years Subse shf: 6 mohs

3,... ;,...,,, (... ( 0 f y j j p p Tred approxmao by p-h order polyomal (regresso polyomal Regresso le Regresso parabola No-lear red modelg Model fuco

4 Tred approxmao by p-h order polyomal (regresso polyomal Regresso le Regresso parabola No-lear red modelg Model fuco,... ;,...,,, (... ( 0 f y j j p p

Model seleco Whch model s relable / o be preferred?? 5

6 (, / 0 ~ / ˆ, 0 :, 0 : p T H H Hypoheses esg of polyomal coeffces (absolue measure Iformao crera (relave measure Leave-oe-ou cross valdao (relave measure CV ˆ CV ( y y Whch model s relable / o be preferred? l( ˆ l( BIC ( AIC AIC, ˆ ˆ, ˆ l( AIC R ( R ( ˆ R c p e - p - - y y σ - R-square Aae IC Bayesa IC Model seleco

Model seleco Whch model s relable / o be preferred? Greelad Order ΔR ΔAIC c ΔBIC ΔCV Sgfca parameers 0.00 0.7 8. 69 all 0 0 0 0 all 3 0. 3.6 5945 ercep 4 0 3.4 8. 8865 ercep Secod-order polyomal s superor o oher red fucos accelerao: -5 ± 3 G/yr Amazo bas Order ΔR ΔAIC c ΔBIC ΔCV Sgfca parameers 0.043 0 0 0 all 0 0.7 3.3 373 oe Frs-order polyomal s superor o oher red fucos 7

Coclusos I bref I s a mess! More scefcally Treds umbers (ad her error bouds are hghly subjec o he uderlyg processg scheme Ierpreao of mass-chage reds ad mass-chage red varaos over shor perods should be avoded Never ever valdae chage raes agas reds umbers based o dfferg perods Tred model seleco should rouely be accompaed by hypoheses esg ad formao crera evaluao The compuao of mass varao from GRACE has bee becomg creasgly heerogeeous; coveos would mprove cossecy 8

Tred model Whch model s relable / o be preferred? CSR Nle bas ITG DMT 9

Tred model Subse perod: 5 years Subse shf: 6 mohs Subse perod: 6 years Subse shf: 6 mohs 0