Influences of different factors on temporal correlations of GNSS observations

Size: px
Start display at page:

Download "Influences of different factors on temporal correlations of GNSS observations"

Transcription

1 2 T : oscillation period (unknown) C : scaling factor (unknown) ND : determined empirically a priori C,T: least-squares regression on sample ACF ND : first zero point of sample ACF h= ACF h π h ND ND T h : lag (epoch distance) h 2 C exp cos for h ND = C for h = 0 ACF E.g. analytical autocorrelation function (ACF) (Howind 2005) Modelling approach using residual time series Signal processing techniques Site-specific influences (e.g. multipath effects) Physical processes in the atmosphere Impacting factors GNSS observations are temporally correlated! Introduction KIT The Cooperation of Forschungszentrum Karlsruhe GmbH and Universität Karlsruhe (TH) Lag ACF X. Luo, M. Mayer, B. Heck Geodetic Institute Influences of different factors on temporal correlations of GNSS observations Geodetic Week 2009 September 22-24, Karlsruhe S4: Applied Geodesy and GNSS

2 Sample ACF Definition: x, K, x n are observations of a time series { X t }. The sample autocovariance function is γˆ( h ) : = n n h t = ( x The sample autocorrelation function is n t + h x )( x t x ) x t ˆ( γ h ) ˆ ρ ( h ) =, n < h < n. ˆ(0) γ n < h < n, where x = n t =. { X t } (weakly) stationary ˆ ρ ( h ) representative for the dependence in the data Distribution approximation of ˆ ρ ( h ) For large n, the ˆ ρ ( h ) of an iid* sequence with finite variance are approx. iid with distribution ˆ( ρ h ) ~ N (0, ). n * iid: independent and identically distributed 3 Case study: data base Multipath (MP): strong Tab. : GPS processing strategies and data base MP: weak TAAF HEDA Observations Processing interval Obs. weighting Tropospheric model (time span: 5 min) GPS -Hz double differences (phase) DOY2007: 6-8, UT: 5-8 h SNR*-based: w=f(snr) (cut-off: 3 ) Niell (dry) (a priori model) Niell (wet) (mapping function) Ambiguity resolution SIGMA strategy (L5, L3) Antenna correction Individual absolute calibration RATA Residual data base (: 285) Studentised double difference residuals (SDDR) in sidereal time Identical length=3600 epochs ( h) Atmos. conditions Meteorological surface data (DWD**) * SNR: Signal-to-Noise Ratio, ** DWD: Deutscher Wetterdienst ( SIBI Tab. 2: Considered factors impacting temporal correlations Factor Comparison between Multipath HEDA (strong) TAAF (weak) 48/5 Baseline length SIBI (42.5 km) RATA(203.7 km) 48/33 Rel. humidity wet days (83 %) dry days (46 %) 3/50 Fig. : SAPOS network in the area of the state of Baden-Württemberg (Southwest Germany) Wind force windy (4-5 m/s) calm (-2 m/s) 26/38 4

3 5 6 Modelling residual trends (a): Conventional approach (b): Improved approach with additional correlation analysis Fig. 2: Comparison of different approaches for modelling residual trends (conventional vs. improved) Modelling residual trends SDDR: SIBI097(all SDDR) Trend-eliminated SDDR: SIBI097(selected SDDR) Sample ACF Fig. 3: Selected representative example of comparison of trend modelling (conventional vs. improved)

4 7 8 Distribution of ND Fig. 4: Empirical distribution of ND (ND : first zero-point of sample ACF) Tab. 3: Statistical characteristics and hypothesis tests for normal distribution of ND Statistical characteristics Test for normality p-value ** 95% quantile Chi-square Jarque-Bera 0.46 *STD: standard deviation, **: interquartile range Multipath effects (MP) outlier HEDA (MP: strong) TAAF (MP: weak) Fig. 5: Influences of multipath effects on correlation length represented by ND Baseline HEDA Tab. 4: Influences of multipath effects on statistical characteristics of ND MP strong Length [km] Statistical characteristics of ND ** TAAF weak

5 9 0 Baseline length RATA (203.7 km) SIBI (42.5 km) Fig. 6: Influences of baseline length on correlation length represented by ND Baseline RATA Tab. 5: Influences of baseline length on statistical characteristics of ND MP weak Length [km] Statistical characteristics of ND ** SIBI weak Relative humidity (rh) Wet days 66, 73, 76 Dry days 67, 68, 70 Fig. 7: Influences of atmospheric relative humidity on correlation length represented by ND Days Wet Tab. 6: Influences of atmospheric relative humidity on statistical characteristics of ND relative humidity 83% 3 88 Statistical characteristics of ND ** Dry 46%

6 2 Wind force Windy days 72, 73, 77 Calm days 65, 68, 69 Fig. 8: Influences of wind force on correlation length represented by ND Days Windy Tab. 7: Influences of wind force on statistical characteristics of ND wind velocity 4-5 m/s Statistical characteristics of ND ** Calm -2 m/s Conclusions and outlook Temporal correlation analysis Residual-based, using sample ACF Assumption of stationary time series Improved approach for modelling residual trends Correlation length Represented by the first zero-point of sample ACF Asymptotic normal distribution (hypothesis tests) correlation length within this study: approx. 00 s Decreasing with stronger multipath effects and wind force Insignificantly influenced by baseline length and relative humidity Further research work Validating the statements related to atmospheric conditions Analysing the influences of satellite geometry on temporal correlations

7 Questions & comments Thank you very much for your attention! The project Improving the stochastic model of GPS observations by modelling physical correlations (HE433/6-/2) is supported by the Deutsche Forschungsgemeinschaft (DFG). GPS system modernisations + advanced mathematical modelling = accurate and reliable positioning results 3 Geodetic Institute Geodetic Universität Institute: X. Karlsruhe Luo, M. Mayer, (TH) B. Englerstraße Heck luo@gik.uni-karlsruhe.de 7, 763 Karlsruhe, Germany Tel.: +49 (0)72 Influences , of different Fax: +49 factors (0)72 on temporal , correlations home page: of GNSS observations

Atmospheric Water Vapor Effect on GNSS Signals and InSAR Data

Atmospheric Water Vapor Effect on GNSS Signals and InSAR Data Atmospheric Water Vapor Effect on GNSS Signals and InSAR Data Basic Concept and Preliminary Results F. Alshawaf 1, S. Hinz 1, A. Thiele 1 T. Fuhrmann 2, B. Heck 2, A. Knöpfler 2, X. Luo 2, M. Mayer 2,

More information

Constructing high-resolution, absolute maps of atmospheric water vapor by combining InSAR and GNSS observations

Constructing high-resolution, absolute maps of atmospheric water vapor by combining InSAR and GNSS observations Constructing high-resolution, absolute maps of atmospheric water vapor by combining InSAR and GNSS observations Fadwa Alshawaf, Stefan Hinz, Michael Mayer, Franz J. Meyer fadwa.alshawaf@kit.edu INSTITUTE

More information

ESTIMATING THE RESIDUAL TROPOSPHERIC DELAY FOR AIRBORNE DIFFERENTIAL GPS POSITIONING (A SUMMARY)

ESTIMATING THE RESIDUAL TROPOSPHERIC DELAY FOR AIRBORNE DIFFERENTIAL GPS POSITIONING (A SUMMARY) ESTIMATING THE RESIDUAL TROPOSPHERIC DELAY FOR AIRBORNE DIFFERENTIAL GPS POSITIONING (A SUMMARY) J. Paul Collins and Richard B. Langley Geodetic Research Laboratory Department of Geodesy and Geomatics

More information

Atmospheric delay. X, Y, Z : satellite cartesian coordinates. Z : receiver cartesian coordinates. In the vacuum the signal speed c is constant

Atmospheric delay. X, Y, Z : satellite cartesian coordinates. Z : receiver cartesian coordinates. In the vacuum the signal speed c is constant Atmospheric delay In the vacuum the signal speed c is constant c τ = ρ = ( X X ) + ( Y Y ) + ( Z Z ) S S S 2 S 2 S 2 X, Y, Z : receiver cartesian coordinates S S S X, Y, Z : satellite cartesian coordinates

More information

A correction model for zenith dry delay of GPS signals using regional meteorological sites. GPS-based determination of atmospheric water vapour

A correction model for zenith dry delay of GPS signals using regional meteorological sites. GPS-based determination of atmospheric water vapour Geodetc Week 00 October 05-07, Cologne S4: Appled Geodesy and GNSS A correcton model for zenth dry delay of GPS sgnals usng regonal meteorologcal stes Xaoguang Luo Geodetc Insttute, Department of Cvl Engneerng,

More information

Financial Econometrics and Quantitative Risk Managenent Return Properties

Financial Econometrics and Quantitative Risk Managenent Return Properties Financial Econometrics and Quantitative Risk Managenent Return Properties Eric Zivot Updated: April 1, 2013 Lecture Outline Course introduction Return definitions Empirical properties of returns Reading

More information

NANYANG TECHNOLOGICAL UNIVERSITY SEMESTER II EXAMINATION MAS451/MTH451 Time Series Analysis TIME ALLOWED: 2 HOURS

NANYANG TECHNOLOGICAL UNIVERSITY SEMESTER II EXAMINATION MAS451/MTH451 Time Series Analysis TIME ALLOWED: 2 HOURS NANYANG TECHNOLOGICAL UNIVERSITY SEMESTER II EXAMINATION 2012-2013 MAS451/MTH451 Time Series Analysis May 2013 TIME ALLOWED: 2 HOURS INSTRUCTIONS TO CANDIDATES 1. This examination paper contains FOUR (4)

More information

Econometría 2: Análisis de series de Tiempo

Econometría 2: Análisis de series de Tiempo Econometría 2: Análisis de series de Tiempo Karoll GOMEZ kgomezp@unal.edu.co http://karollgomez.wordpress.com Segundo semestre 2016 II. Basic definitions A time series is a set of observations X t, each

More information

On the realistic stochastic model of GPS observables: Implementation and Performance

On the realistic stochastic model of GPS observables: Implementation and Performance he International Archives of the Photogrammetry, Remote Sensing Spatial Information Sciences, Volume XL-/W5, 05 International Conference on Sensors & Models in Remote Sensing & Photogrammetry, 3 5 Nov

More information

Noise Characteristics in High Precision GPS Positioning

Noise Characteristics in High Precision GPS Positioning Noise Characteristics in High Precision GPS Positioning A.R. Amiri-Simkooei, C.C.J.M. Tiberius, P.J.G. Teunissen, Delft Institute of Earth Observation and Space systems (DEOS), Delft University of Technology,

More information

Minitab Project Report Assignment 3

Minitab Project Report Assignment 3 3.1.1 Simulation of Gaussian White Noise Minitab Project Report Assignment 3 Time Series Plot of zt Function zt 1 0. 0. zt 0-1 0. 0. -0. -0. - -3 1 0 30 0 50 Index 0 70 0 90 0 1 1 1 1 0 marks The series

More information

A. Barbu, J. Laurent-Varin, F. Perosanz, F. Mercier and J. Marty. AVENUE project. June, 20

A. Barbu, J. Laurent-Varin, F. Perosanz, F. Mercier and J. Marty. AVENUE project. June, 20 Efficient QR Sequential Least Square algorithm for high frequency GNSS Precise Point Positioning A. Barbu, J. Laurent-Varin, F. Perosanz, F. Mercier and J. Marty AVENUE project June, 20 A. Barbu, J. Laurent-Varin,

More information

MATRAG Measurement of Alpine Tropospheric Delay by Radiometer and GPS

MATRAG Measurement of Alpine Tropospheric Delay by Radiometer and GPS MATRAG Measurement of Alpine Tropospheric Delay by Radiometer and GPS Petra Häfele 1, Matthias Becker, Elmar Brockmann, Lorenz Martin, Michael Kirchner 1 University of the Bundeswehr Munich, 85577 Neubiberg,

More information

Analysis of the Accuracy of GMF, NMF, and VMF1 Mapping Functions with GPT 50 a Priori Zenith Constraint in Tropospheric Delay Modelling

Analysis of the Accuracy of GMF, NMF, and VMF1 Mapping Functions with GPT 50 a Priori Zenith Constraint in Tropospheric Delay Modelling Analysis of the Accuracy of GMF, NMF, and VMF1 Mapping Functions with GPT 50 a Priori Zenith Constraint in Tropospheric Delay Modelling Brian Makabayi 1 Addisu Hunegnaw 2 1 Assistant Lecturer, Department

More information

Use of ground-based GNSS measurements in data assimilation. Reima Eresmaa Finnish Meteorological Institute

Use of ground-based GNSS measurements in data assimilation. Reima Eresmaa Finnish Meteorological Institute Use of ground-based GNSS measurements in data assimilation Reima Eresmaa Finnish Meteorological Institute 16 June 2006 Outline 1) Introduction GNSS * positioning Tropospheric delay 2) GNSS as a meteorological

More information

GNSS-specific local effects at the Geodetic Observatory Wettzell

GNSS-specific local effects at the Geodetic Observatory Wettzell GNSS-specific local effects at the Geodetic Observatory Wettzell Peter Steigenberger, Urs Hugentobler, Ralf Schmid Technische Universität München (TUM) Uwe Hessels, Thomas Klügel Bundesamt für Kartographie

More information

Module 3. Descriptive Time Series Statistics and Introduction to Time Series Models

Module 3. Descriptive Time Series Statistics and Introduction to Time Series Models Module 3 Descriptive Time Series Statistics and Introduction to Time Series Models Class notes for Statistics 451: Applied Time Series Iowa State University Copyright 2015 W Q Meeker November 11, 2015

More information

Prof. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis

Prof. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation

More information

The increasing intensity of the strongest tropical cyclones

The increasing intensity of the strongest tropical cyclones The increasing intensity of the strongest tropical cyclones James B. Elsner Department of Geography, Florida State University Tallahassee, Florida Corresponding author address: Dept. of Geography, The

More information

Satellite Navigation error sources and position estimation

Satellite Navigation error sources and position estimation Satellite Navigation error sources and position estimation Picture: ESA AE4E08 Sandra Verhagen Course 2010 2011, lecture 6 1 Today s topics Recap: GPS measurements and error sources Signal propagation

More information

Benefits of State Space Modeling in GNSS Multi-Station Adjustment

Benefits of State Space Modeling in GNSS Multi-Station Adjustment Benefits of State Space Modeling in GNSS Multi-Station Adjustment Andreas Bagge, Gerhard Wübbena Martin Schmitz Geo++ GmbH D-30827 Garbsen, Germany www.geopp.de GeoInformation Workshop 2004, Istanbul Kultur

More information

LECTURE 10: MORE ON RANDOM PROCESSES

LECTURE 10: MORE ON RANDOM PROCESSES LECTURE 10: MORE ON RANDOM PROCESSES AND SERIAL CORRELATION 2 Classification of random processes (cont d) stationary vs. non-stationary processes stationary = distribution does not change over time more

More information

Humidity 3D field comparisons between GNSS tomography, IASI satellite observations and ALARO model. Belgian Institute for Space Aeronomy BIRA 3

Humidity 3D field comparisons between GNSS tomography, IASI satellite observations and ALARO model. Belgian Institute for Space Aeronomy BIRA 3 Oral Presentation, EGU0-85 Humidity D field comparisons between, H. Brenot, C. Champollion, A. Deckmyn, R. van Malderen, N. Kumps, R. Warnant, E. Goudenhoofdt, L. Delobbe and M. De Mazière contact: Belgian

More information

Econ 424 Time Series Concepts

Econ 424 Time Series Concepts Econ 424 Time Series Concepts Eric Zivot January 20 2015 Time Series Processes Stochastic (Random) Process { 1 2 +1 } = { } = sequence of random variables indexed by time Observed time series of length

More information

Modern Navigation. Thomas Herring

Modern Navigation. Thomas Herring 12.215 Modern Navigation Thomas Herring Basic Statistics Summary of last class Statistical description and parameters Probability distributions Descriptions: expectations, variances, moments Covariances

More information

Impact of Tropospheric Delay Gradients on Total Tropospheric Delay and Precise Point Positioning

Impact of Tropospheric Delay Gradients on Total Tropospheric Delay and Precise Point Positioning International Journal of Geosciences, 016, 7, 645-654 Published Online May 016 in SciRes. http://www.scirp.org/journal/ijg http://dx.doi.org/10.436/ijg.016.75050 Impact of Tropospheric Delay Gradients

More information

Reliability and Risk Analysis. Time Series, Types of Trend Functions and Estimates of Trends

Reliability and Risk Analysis. Time Series, Types of Trend Functions and Estimates of Trends Reliability and Risk Analysis Stochastic process The sequence of random variables {Y t, t = 0, ±1, ±2 } is called the stochastic process The mean function of a stochastic process {Y t} is the function

More information

Satellite Navigation PVT estimation and integrity

Satellite Navigation PVT estimation and integrity Satellite Navigation PVT estimation and integrity Picture: ESA AE4E8 Sandra Verhagen Course 1 11, lecture 7 1 Satellite Navigation (AE4E8 Lecture 7 Today s topics Position, Velocity and Time estimation

More information

On the statistical significance of climatic trends estimated from GRUAN tropospheric time series

On the statistical significance of climatic trends estimated from GRUAN tropospheric time series On the statistical significance of climatic trends estimated from GRUAN tropospheric time series Fadwa Alshawaf, Galina Dick, Jens Wickert 1 On the statistical significance of climatic trends estimated

More information

Econometrics of Panel Data

Econometrics of Panel Data Econometrics of Panel Data Jakub Mućk Meeting # 9 Jakub Mućk Econometrics of Panel Data Meeting # 9 1 / 22 Outline 1 Time series analysis Stationarity Unit Root Tests for Nonstationarity 2 Panel Unit Root

More information

EESC Geodesy with the Global Positioning System. Class 7: Relative Positioning using Carrier-Beat Phase

EESC Geodesy with the Global Positioning System. Class 7: Relative Positioning using Carrier-Beat Phase EESC 9945 Geodesy with the Global Positioning System Class 7: Relative Positioning using Carrier-Beat Phase GPS Carrier Phase The model for the carrier-beat phase observable for receiver p and satellite

More information

LECTURES 2-3 : Stochastic Processes, Autocorrelation function. Stationarity.

LECTURES 2-3 : Stochastic Processes, Autocorrelation function. Stationarity. LECTURES 2-3 : Stochastic Processes, Autocorrelation function. Stationarity. Important points of Lecture 1: A time series {X t } is a series of observations taken sequentially over time: x t is an observation

More information

Non-Stationary Time Series and Unit Root Testing

Non-Stationary Time Series and Unit Root Testing Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity

More information

Problem Set 1 Solution Sketches Time Series Analysis Spring 2010

Problem Set 1 Solution Sketches Time Series Analysis Spring 2010 Problem Set 1 Solution Sketches Time Series Analysis Spring 2010 1. Construct a martingale difference process that is not weakly stationary. Simplest e.g.: Let Y t be a sequence of independent, non-identically

More information

The Power of the KPSS Test for Cointegration when Residuals are Fractionally Integrated 1

The Power of the KPSS Test for Cointegration when Residuals are Fractionally Integrated 1 The Power of the KPSS Test for Cointegration when Residuals are Fractionally Integrated 1 by Philipp Sibbertsen 2 and Walter Krämer Fachbereich Statistik, Universität Dortmund, D-44221 Dortmund, Germany

More information

Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data

Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data Johannes Böhm, Arthur Niell, Paul Tregoning, and Harald Schuh Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data Geophysical Research Letters Vol. 33,

More information

Nonlinear time series

Nonlinear time series Based on the book by Fan/Yao: Nonlinear Time Series Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 27, 2009 Outline Characteristics of

More information

The Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data

The Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data The Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data J. Boehm, A. Niell, P. Tregoning, H. Schuh Troposphere mapping functions are used in the analyses

More information

Introduction to Economic Time Series

Introduction to Economic Time Series Econometrics II Introduction to Economic Time Series Morten Nyboe Tabor Learning Goals 1 Give an account for the important differences between (independent) cross-sectional data and time series data. 2

More information

REFINED AND SITE-AUGMENTED TROPOSPHERIC DELAY MODELS FOR GNSS

REFINED AND SITE-AUGMENTED TROPOSPHERIC DELAY MODELS FOR GNSS REFINED AND SITE-AUGMENTED TROPOSPHERIC DELAY MODELS FOR GNSS Daniel Landskron, Gregor Möller, Armin Hofmeister, Johannes Böhm, and Robert Weber Technische Universität Wien, Austria Gußhausstraße 27-29,

More information

APPLIED ECONOMETRIC TIME SERIES 4TH EDITION

APPLIED ECONOMETRIC TIME SERIES 4TH EDITION APPLIED ECONOMETRIC TIME SERIES 4TH EDITION Chapter 2: STATIONARY TIME-SERIES MODELS WALTER ENDERS, UNIVERSITY OF ALABAMA Copyright 2015 John Wiley & Sons, Inc. Section 1 STOCHASTIC DIFFERENCE EQUATION

More information

Empirical Market Microstructure Analysis (EMMA)

Empirical Market Microstructure Analysis (EMMA) Empirical Market Microstructure Analysis (EMMA) Lecture 3: Statistical Building Blocks and Econometric Basics Prof. Dr. Michael Stein michael.stein@vwl.uni-freiburg.de Albert-Ludwigs-University of Freiburg

More information

STAT 520: Forecasting and Time Series. David B. Hitchcock University of South Carolina Department of Statistics

STAT 520: Forecasting and Time Series. David B. Hitchcock University of South Carolina Department of Statistics David B. University of South Carolina Department of Statistics What are Time Series Data? Time series data are collected sequentially over time. Some common examples include: 1. Meteorological data (temperatures,

More information

Near Real Time atmosphere model based on GNSS and meteorological data from ASG-EUPOS reference stations

Near Real Time atmosphere model based on GNSS and meteorological data from ASG-EUPOS reference stations 10th Czech-Polish Workshop: Szklarska Porba / Poland - November 5-7, 2009 1/20 Near Real Time atmosphere model based on GNSS and meteorological data from ASG-EUPOS reference stations Bosy J. (1), Rohm

More information

Climate Monitoring with Radio Occultation Data

Climate Monitoring with Radio Occultation Data Climate Monitoring with Radio Occultation Data Systematic Error Sources C. Rocken, S. Sokolovskiy, B. Schreiner, D. Hunt, B. Ho, B. Kuo, U. Foelsche Radio Occultation Claims Most stable Global Thermometer

More information

Econometrics I: Univariate Time Series Econometrics (1)

Econometrics I: Univariate Time Series Econometrics (1) Econometrics I: Dipartimento di Economia Politica e Metodi Quantitativi University of Pavia Overview of the Lecture 1 st EViews Session VI: Some Theoretical Premises 2 Overview of the Lecture 1 st EViews

More information

Deutscher Wetterdienst

Deutscher Wetterdienst Deutscher Wetterdienst Modelling the Volcanic Ash Episode: Experiences with COSMO-ART Detlev Majewski (FE1) Bernhard Vogel, Heike Vogel (KIT) Thomas Hanisch, Jochen Förstner (FE13), Ulrich Pflüger (FE15)

More information

Romanian Economic and Business Review Vol. 3, No. 3 THE EVOLUTION OF SNP PETROM STOCK LIST - STUDY THROUGH AUTOREGRESSIVE MODELS

Romanian Economic and Business Review Vol. 3, No. 3 THE EVOLUTION OF SNP PETROM STOCK LIST - STUDY THROUGH AUTOREGRESSIVE MODELS THE EVOLUTION OF SNP PETROM STOCK LIST - STUDY THROUGH AUTOREGRESSIVE MODELS Marian Zaharia, Ioana Zaheu, and Elena Roxana Stan Abstract Stock exchange market is one of the most dynamic and unpredictable

More information

Stationary and nonstationary variables

Stationary and nonstationary variables Stationary and nonstationary variables Stationary variable: 1. Finite and constant in time expected value: E (y t ) = µ < 2. Finite and constant in time variance: Var (y t ) = σ 2 < 3. Covariance dependent

More information

High Rate GPS Solutions

High Rate GPS Solutions High Rate GPS Solutions High rate GPS data (1 Hz or higher) Network solution Fixed a local reference clock Bias fixed Sub daily position estimates solutions Position becomes stochastic parameter Fairly

More information

Lecture 3 Stationary Processes and the Ergodic LLN (Reference Section 2.2, Hayashi)

Lecture 3 Stationary Processes and the Ergodic LLN (Reference Section 2.2, Hayashi) Lecture 3 Stationary Processes and the Ergodic LLN (Reference Section 2.2, Hayashi) Our immediate goal is to formulate an LLN and a CLT which can be applied to establish sufficient conditions for the consistency

More information

Problem set 1 - Solutions

Problem set 1 - Solutions EMPIRICAL FINANCE AND FINANCIAL ECONOMETRICS - MODULE (8448) Problem set 1 - Solutions Exercise 1 -Solutions 1. The correct answer is (a). In fact, the process generating daily prices is usually assumed

More information

Investigations Concerning the Reliability and the External Accuracy of GPS Real-Time Measurements

Investigations Concerning the Reliability and the External Accuracy of GPS Real-Time Measurements Investigations Concerning the Reliability and the External Accuracy of GPS Real-Time Measurements Dr. Michael ILLNER, Germany Key words: real time, GPS, network adjustment, horizontal position, height,

More information

Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation

Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation Z. W. LI, X. L. DING Department of Land Surveying and Geo-Informatics, Hong Kong Polytechnic University, Hung

More information

STOCHASTIC MODELING OF MONTHLY RAINFALL AT KOTA REGION

STOCHASTIC MODELING OF MONTHLY RAINFALL AT KOTA REGION STOCHASTIC MODELIG OF MOTHLY RAIFALL AT KOTA REGIO S. R. Bhakar, Raj Vir Singh, eeraj Chhajed and Anil Kumar Bansal Department of Soil and Water Engineering, CTAE, Udaipur, Rajasthan, India E-mail: srbhakar@rediffmail.com

More information

An application of the GAM-PCA-VAR model to respiratory disease and air pollution data

An application of the GAM-PCA-VAR model to respiratory disease and air pollution data An application of the GAM-PCA-VAR model to respiratory disease and air pollution data Márton Ispány 1 Faculty of Informatics, University of Debrecen Hungary Joint work with Juliana Bottoni de Souza, Valdério

More information

Sebastian Strasser, Torsten Mayer-Gürr

Sebastian Strasser, Torsten Mayer-Gürr Sebastian Strasser, Torsten Mayer-Gürr Institute of Geodesy, Graz University of Technology WG Theoretical Geodesy and Satellite Geodesy Geodetic Week 2015, Stuttgart, Germany Sebastian Strasser 16.09.2015

More information

EESC Geodesy with the Global Positioning System. Class 6: Point Positioning using Pseuduorange

EESC Geodesy with the Global Positioning System. Class 6: Point Positioning using Pseuduorange EESC 9945 Geodesy with the Global Positioning System Class 6: Point Positioning using Pseuduorange GPS Positioning Solutions Point Positioning: Determination of the coordinates of a point with respect

More information

Carrier Phase Integer Ambiguity Resolution Recent results and open issues

Carrier Phase Integer Ambiguity Resolution Recent results and open issues Carrier Phase Integer Ambiguity Resolution Recent results and open issues Sandra Verhagen DEOS, Delft University of Technology, The Netherlands Kolloquium Satellitennavigation, TU München, 9 February 2010

More information

Approximation of ambiguity covariance matrix for integer de-correlation procedure in single-epoch GNSS positioning

Approximation of ambiguity covariance matrix for integer de-correlation procedure in single-epoch GNSS positioning he 9 th International Conference ENVIRONMENAL ENGINEERING 22 23 May 24, Vilnius, Lithuania SELECED PAPERS eissn 229-792 / eisbn 978-69-457-64-9 Available online at http://enviro.vgtu.lt Section: echnologies

More information

Time Series Analysis -- An Introduction -- AMS 586

Time Series Analysis -- An Introduction -- AMS 586 Time Series Analysis -- An Introduction -- AMS 586 1 Objectives of time series analysis Data description Data interpretation Modeling Control Prediction & Forecasting 2 Time-Series Data Numerical data

More information

On 1.9, you will need to use the facts that, for any x and y, sin(x+y) = sin(x) cos(y) + cos(x) sin(y). cos(x+y) = cos(x) cos(y) - sin(x) sin(y).

On 1.9, you will need to use the facts that, for any x and y, sin(x+y) = sin(x) cos(y) + cos(x) sin(y). cos(x+y) = cos(x) cos(y) - sin(x) sin(y). On 1.9, you will need to use the facts that, for any x and y, sin(x+y) = sin(x) cos(y) + cos(x) sin(y). cos(x+y) = cos(x) cos(y) - sin(x) sin(y). (sin(x)) 2 + (cos(x)) 2 = 1. 28 1 Characteristics of Time

More information

GPS Geodesy - LAB 7. Neglecting the propagation, multipath, and receiver errors, eq.(1) becomes:

GPS Geodesy - LAB 7. Neglecting the propagation, multipath, and receiver errors, eq.(1) becomes: GPS Geodesy - LAB 7 GPS pseudorange position solution The pseudorange measurements j R i can be modeled as: j R i = j ρ i + c( j δ δ i + ΔI + ΔT + MP + ε (1 t = time of epoch j R i = pseudorange measurement

More information

Usage of a Correction Model to Enhance the Evaluation of the Zenith Tropospheric Delay

Usage of a Correction Model to Enhance the Evaluation of the Zenith Tropospheric Delay International Journal of Applied Engineering Research ISSN 097-562 Volume 11, Number 6 (2016) pp 68-65 Usage of a Correction Model to Enhance the Evaluation of the Zenith Tropospheric Delay Meriem Jgouta,

More information

Homework 2. For the homework, be sure to give full explanations where required and to turn in any relevant plots.

Homework 2. For the homework, be sure to give full explanations where required and to turn in any relevant plots. Homework 2 1 Data analysis problems For the homework, be sure to give full explanations where required and to turn in any relevant plots. 1. The file berkeley.dat contains average yearly temperatures for

More information

Testing for IID Noise/White Noise: I

Testing for IID Noise/White Noise: I Testing for IID Noise/White Noise: I want to be able to test null hypothesis time series {x t } or set of residuals {r t } is IID(0, 2 ) or WN(0, 2 ) there are many such tests, including informal test

More information

Introduction to Least Squares Adjustment for geodetic VLBI

Introduction to Least Squares Adjustment for geodetic VLBI Introduction to Least Squares Adjustment for geodetic VLBI Matthias Schartner a, David Mayer a a TU Wien, Department of Geodesy and Geoinformation Least Squares Adjustment why? observation is τ (baseline)

More information

Evaluation of the impact of atmospheric pressure loading modeling on GNSS data analysis

Evaluation of the impact of atmospheric pressure loading modeling on GNSS data analysis Evaluation of the impact of atmospheric pressure loading modeling on GNSS data analysis R. Dach a, J. Böhm b, S. Lutz a, and P. Steigenberger c a Astronomical Institute, University of Bern, Bern, Switzerland

More information

Chapter 2: Unit Roots

Chapter 2: Unit Roots Chapter 2: Unit Roots 1 Contents: Lehrstuhl für Department Empirische of Wirtschaftsforschung Empirical Research and undeconometrics II. Unit Roots... 3 II.1 Integration Level... 3 II.2 Nonstationarity

More information

A time series is called strictly stationary if the joint distribution of every collection (Y t

A time series is called strictly stationary if the joint distribution of every collection (Y t 5 Time series A time series is a set of observations recorded over time. You can think for example at the GDP of a country over the years (or quarters) or the hourly measurements of temperature over a

More information

ERAD Water vapor observations with SAR, microwave radiometer and GPS: comparison of scaling characteristics

ERAD Water vapor observations with SAR, microwave radiometer and GPS: comparison of scaling characteristics Proceedings of ERAD (2002): 190 194 c Copernicus GmbH 2002 ERAD 2002 Water vapor observations with SAR, microwave radiometer and GPS: comparison of scaling characteristics D. N. Moisseev 1, R. F. Hanssen

More information

SIGMA-F: Variances of GPS Observations Determined by a Fuzzy System

SIGMA-F: Variances of GPS Observations Determined by a Fuzzy System SIGMA-F: Variances of GPS Observations Determined by a Fuzzy System A. Wieser and F.K. Brunner Engineering Surveying and Metrology, Graz University of Technology, Steyrergasse 3, A-8 Graz, Austria Keywords.

More information

WG1 Overview. PP KENDA for km-scale EPS: LETKF. current DA method: nudging. radar reflectivity (precip): latent heat nudging 1DVar (comparison)

WG1 Overview. PP KENDA for km-scale EPS: LETKF. current DA method: nudging. radar reflectivity (precip): latent heat nudging 1DVar (comparison) WG1 Overview Deutscher Wetterdienst, D-63067 Offenbach, Germany current DA method: nudging PP KENDA for km-scale EPS: LETKF radar reflectivity (precip): latent heat nudging 1DVar (comparison) radar radial

More information

THE EFFECT OF PHYSICAL CORRELATIONS ON THE AMBIGUITY RESOLUTION AND ACCURACY ESTIMATION IN GPS DIFFERENTIAL POSITIONING

THE EFFECT OF PHYSICAL CORRELATIONS ON THE AMBIGUITY RESOLUTION AND ACCURACY ESTIMATION IN GPS DIFFERENTIAL POSITIONING THE EFFECT OF PHYSICAL CORRELATIONS ON THE AMBIGUITY RESOLUTION AND ACCURACY ESTIMATION IN GPS DIFFERENTIAL POSITIONING A. E-S. EL RABBANY May 1994 TECHNICAL REPORT NO. 170 PREFACE In order to make our

More information

Atmospheric phase correction for ALMA with water-vapour radiometers

Atmospheric phase correction for ALMA with water-vapour radiometers Atmospheric phase correction for ALMA with water-vapour radiometers B. Nikolic Cavendish Laboratory, University of Cambridge January 29 NA URSI, Boulder, CO B. Nikolic (University of Cambridge) WVR phase

More information

PAPER 206 APPLIED STATISTICS

PAPER 206 APPLIED STATISTICS MATHEMATICAL TRIPOS Part III Thursday, 1 June, 2017 9:00 am to 12:00 pm PAPER 206 APPLIED STATISTICS Attempt no more than FOUR questions. There are SIX questions in total. The questions carry equal weight.

More information

Applied time-series analysis

Applied time-series analysis Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 18, 2011 Outline Introduction and overview Econometric Time-Series Analysis In principle,

More information

Cross-validation methods for quality control, cloud screening, etc.

Cross-validation methods for quality control, cloud screening, etc. Cross-validation methods for quality control, cloud screening, etc. Olaf Stiller, Deutscher Wetterdienst Are observations consistent Sensitivity functions with the other observations? given the background

More information

P. Cipollini, H. Snaith - A short course on Altimetry. Altimetry 2 - Data processing (from satellite height to sea surface height)

P. Cipollini, H. Snaith - A short course on Altimetry. Altimetry 2 - Data processing (from satellite height to sea surface height) P. Cipollini, H. Snaith - A short course on Altimetry Altimetry 2 - Data processing (from satellite height to sea surface height) 1 2 Satellite height to sea surface height The altimeter measures the altitude

More information

Developments at DWD: Integrated water vapour (IWV) from ground-based GPS

Developments at DWD: Integrated water vapour (IWV) from ground-based GPS 1 Working Group on Data Assimilation 2 Developments at DWD: Integrated water vapour (IWV) from ground-based Christoph Schraff, Maria Tomassini, and Klaus Stephan Deutscher Wetterdienst, Frankfurter Strasse

More information

Økonomisk Kandidateksamen 2004 (I) Econometrics 2

Økonomisk Kandidateksamen 2004 (I) Econometrics 2 Økonomisk Kandidateksamen 2004 (I) Econometrics 2 This is a closed-book exam (uden hjælpemidler). Answer all questions! The group of questions 1 to 4 have equal weight. Within each group, part (a) represents

More information

NGA GNSS Division Precise Ephemeris Parameters

NGA GNSS Division Precise Ephemeris Parameters NGA GNSS Division Precise Ephemeris Parameters Precise Ephemeris Units. Earth-centered, Earth-fixed Coordinate system Position Velocity GPS time Trajectory interval Standard Trajectory Optional Trajectory

More information

Definition of a Stochastic Process

Definition of a Stochastic Process Definition of a Stochastic Process Balu Santhanam Dept. of E.C.E., University of New Mexico Fax: 505 277 8298 bsanthan@unm.edu August 26, 2018 Balu Santhanam (UNM) August 26, 2018 1 / 20 Overview 1 Stochastic

More information

Part II. Time Series

Part II. Time Series Part II Time Series 12 Introduction This Part is mainly a summary of the book of Brockwell and Davis (2002). Additionally the textbook Shumway and Stoffer (2010) can be recommended. 1 Our purpose is to

More information

Lecture 2: ARMA(p,q) models (part 2)

Lecture 2: ARMA(p,q) models (part 2) Lecture 2: ARMA(p,q) models (part 2) Florian Pelgrin University of Lausanne, École des HEC Department of mathematics (IMEA-Nice) Sept. 2011 - Jan. 2012 Florian Pelgrin (HEC) Univariate time series Sept.

More information

Time Series Analysis

Time Series Analysis Time Series Analysis Christopher Ting http://mysmu.edu.sg/faculty/christophert/ christopherting@smu.edu.sg Quantitative Finance Singapore Management University March 3, 2017 Christopher Ting Week 9 March

More information

Non-Stationary Time Series and Unit Root Testing

Non-Stationary Time Series and Unit Root Testing Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity

More information

Autoregressive Moving Average (ARMA) Models and their Practical Applications

Autoregressive Moving Average (ARMA) Models and their Practical Applications Autoregressive Moving Average (ARMA) Models and their Practical Applications Massimo Guidolin February 2018 1 Essential Concepts in Time Series Analysis 1.1 Time Series and Their Properties Time series:

More information

Review Session: Econometrics - CLEFIN (20192)

Review Session: Econometrics - CLEFIN (20192) Review Session: Econometrics - CLEFIN (20192) Part II: Univariate time series analysis Daniele Bianchi March 20, 2013 Fundamentals Stationarity A time series is a sequence of random variables x t, t =

More information

Topic 1: Atmosphere and Climate

Topic 1: Atmosphere and Climate Topic 1: Atmosphere and Climate Peter Braesicke Fügen Sie auf der Masterfolie ein frei wählbares Bild ein (z.b. passend zum Vortrag) KIT Universität des Landes Baden-Württemberg und nationales Forschungszentrum

More information

Impact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames

Impact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames Impact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames J. Böhm, H. Spicakova, L. Urquhart, P. Steigenberger, H. Schuh Abstract We compare the influence of two different a priori gradient

More information

Torsten Mayer-Gürr Institute of Geodesy, NAWI Graz Technische Universität Graz

Torsten Mayer-Gürr Institute of Geodesy, NAWI Graz Technische Universität Graz GGOS and Reference Systems Introduction 2015-10-12 Torsten Mayer-Gürr Institute of Geodesy, NAWI Graz Technische Universität Graz Torsten Mayer-Gürr 1 Course and exam Lecture Monday 14:00 16:00, A111 (ST01044)

More information

Nonlinear Time Series Modeling

Nonlinear Time Series Modeling Nonlinear Time Series Modeling Part II: Time Series Models in Finance Richard A. Davis Colorado State University (http://www.stat.colostate.edu/~rdavis/lectures) MaPhySto Workshop Copenhagen September

More information

Impact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames

Impact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames Impact of A Priori Gradients on VLBI-Derived Terrestrial Reference Frames J. Böhm, H. Spicakova, L. Urquhart, P. Steigenberger, H. Schuh Abstract We compare the influence of two different a priori gradient

More information

Large sample distribution for fully functional periodicity tests

Large sample distribution for fully functional periodicity tests Large sample distribution for fully functional periodicity tests Siegfried Hörmann Institute for Statistics Graz University of Technology Based on joint work with Piotr Kokoszka (Colorado State) and Gilles

More information

Dirk Schlabing and András Bárdossy. Comparing Five Weather Generators in Terms of Entropy

Dirk Schlabing and András Bárdossy. Comparing Five Weather Generators in Terms of Entropy Dirk Schlabing and András Bárdossy Comparing Five Weather Generators in Terms of Entropy Motivation 1 Motivation What properties of weather should be reproduced [...]? Dirk Schlabing & András Bárdossy,

More information

MAT3379 (Winter 2016)

MAT3379 (Winter 2016) MAT3379 (Winter 2016) Assignment 4 - SOLUTIONS The following questions will be marked: 1a), 2, 4, 6, 7a Total number of points for Assignment 4: 20 Q1. (Theoretical Question, 2 points). Yule-Walker estimation

More information

Downscaling in Time. Andrew W. Robertson, IRI. Advanced Training Institute on Climate Variability and Food Security, 12 July 2002

Downscaling in Time. Andrew W. Robertson, IRI. Advanced Training Institute on Climate Variability and Food Security, 12 July 2002 Downscaling in Time Andrew W. Robertson, IRI Advanced Training Institute on Climate Variability and Food Security, 12 July 2002 Preliminaries Crop yields are driven by daily weather variations! Current

More information

FIN822 project 2 Project 2 contains part I and part II. (Due on November 10, 2008)

FIN822 project 2 Project 2 contains part I and part II. (Due on November 10, 2008) FIN822 project 2 Project 2 contains part I and part II. (Due on November 10, 2008) Part I Logit Model in Bankruptcy Prediction You do not believe in Altman and you decide to estimate the bankruptcy prediction

More information

Statistics of stochastic processes

Statistics of stochastic processes Introduction Statistics of stochastic processes Generally statistics is performed on observations y 1,..., y n assumed to be realizations of independent random variables Y 1,..., Y n. 14 settembre 2014

More information