Presentación realizada en: IAG-IASPEI 39th Joint Scientific Assembly celebrada en Kobe, Japón, del 30 de julio al 4 de agosto de 2017

Bock, Olivier

Bogusz, J.

Chimani, B.

Elias, M.

Gruszczynska, M.

Guijarro, José Antonio

Klos, A.

Malderen, Roeland van

Ning, T.

Pottiaux, E.

Zengin Kazancı, S.

English

Klos, Anna

Van Malderen, Roeland

Pottiaux, Eric

Bogusz, Janusz

Chimani, Barbara

Elias, Michal

Gruszczynska, Marta

Guijarro Pastor, José Antonio

Zengin Kazancı, Selma

Ning, Tong

Agencia Estatal de Meteorología

Performance of various homogenization tools on a synthetic benchmark dataset of GPS and ERA-interim IWV differencesAnna Klos1), Roeland Van Malderen2), Eric Pottiaux3), Olivier Bock4), Janusz Bogusz1), Barbara Chimani5), Michal Elias6), Marta Gruszczynska1), José Guijarro7), Selma Zengin Kazancı8), Tong Ning9)1) Military University of Technology, Warsaw, Poland,2) Royal Meteorological Institute of Belgium (RMI), Brussels, Belgium,3) Royal Observatory of Belgium (ROB), Brussels, Belgium,4) IGN LAREG, University Paris Diderot, Sorbonne Paris, France,5) Central Institute for Meteorology and Geodynamics, Austria,6) Research Institute of Geodesy, Topography and Cartography, Czech Republic,7) AEMET (Spanish Meteorological Agency), Spain,8) Karadeniz Technical University, Turkey,9) Lantmäteriet, Sweden.Motivation Data Synthetic benchmark   Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan1IAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan2Motivation & Introduction:1. COST action GNSS4SWEC „Advanced Global Navigation SatelliteSystems tropospheric products for monitoring severe weatherevents and climate”, WG3: Use of GNSS tropospheric productsfor climate monitoring.2. A proper homogenization of tropospheric dataset isindispensable, as the parameters of deterministic part,e.g. trend will be influenced by undetected breaks.3. Different groups / different methods / different estimates– the truth is not known.4. A synthetic benchmark dataset: a way to quantify results givenby various algorithms.Motivation Data Synthetic benchmark   Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan3Motivation & Introduction:How does it look in practice?A change in trend possible and very likely!What we aim at? Only real breakpoints not regime-like shifts should be corrected!Can anyone see more offsets?Motivation Data Synthetic benchmark   Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan4Data:1. IGS „repro1” troposphere products screened and converted to Integrated Water Vapor (IWV) by O. Bock. 2. 120 stations, daily observations, a period of 1995-2010.Dataset available at: https://doi.org/10.14768/06337394-73a9-407c-9997-0e380dac5590.Motivation Data Synthetic benchmark   Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan5Data:1. IGS „repro1” troposphere products screened and converted to Integrated Water Vapor (IWV) by O. Bock. 2. 120 stations, daily observations, a period of 1995-2010. 3. The IWV differences: ERAI-GPS were employed.Motivation Data Synthetic benchmark   Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan6Real data analysis:1. Manual homogenization + IGS log files: 221 epochs.2. Analysis of significant frequencies: Power Spectral Densities.3. Maximum Likelihood Estimation (MLE) employed.4. Noise analysis: AR(1)+WN chosen as the preferred noise model.Motivation Data Synthetic benchmark   Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan7Generation of the benchmark:3 variants of synthetic time series were generated:1. EASY dataset: seasonal signals + offsets + white noise (WN),2. LESS COMPLICATED dataset: same as 1. + autoregressive process of the first order (noise model = AR(1)+WN),3. FULLY COMPLICATED dataset: same as 2. + trend + gaps (up to 20% of missing data).Motivation Data Synthetic benchmark Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan8Generation of the benchmark:1. 120 series in each synthetic dataset simulated.2. Deterministic model of data taken directly from real differences: trend, seasonal signals, noise.3. Offsets simulated randomly.4. Number of offsets and exact epochs are blinded.Motivation Data Synthetic benchmark Tools   Results SummaryCan anyone see more offsets?IAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan9Verfication of the benchmark:FULLY-COMPLICATEDCoefficients of AR agree within 0.05.Amplitudes & trends agree with – see histograms.Motivation Data Synthetic benchmark Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan10Expected trend uncertainty:Preliminary estimation of which trend uncertainty can be expected from real differences, based on synthetic differences.Motivation Data Synthetic benchmark Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan11Algorithms:1. Sensitivity analysis: the identification of the epochs of the inserted breakpoints.2. Estimates of the trends of the 3 sets of synthetic IWV differences.Motivation Data Synthetic benchmark   Tools Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan12How to classify breaks?Defining a proper timewindow – 2 monthsMotivation Data Synthetic benchmark   Tools Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan13Offsets:Motivation Data Synthetic benchmark   Tools   Results SummaryAmplitudes of reported offsets:EASY, DAILY (SIM: 291):• method 1: 211,• method 2: 252,• method 3: 377,• method 5: 216.IAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan14Offsets:Motivation Data Synthetic benchmark   Tools   Results SummaryAmplitudes of reported offsets:FULLY-COMPLICATED, DAILY (SIM: 317):• method 1: 295,• method 2: 386,• method 3: 622,• method 5: 264.IAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan15Toolsperfor-mance:Motivation Data Synthetic benchmark   Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan16Trends:Motivation Data Synthetic benchmark   Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan17What’s next?1. A detailed assessment of tools – their sensitivity.2. Now – epochs were given to the participants to fine-tune their methods.3. Still looking for other contributions. Interested? Contact Roeland: roeland@meteo.be, Eric: eric.pottiaux@observatoire.be or me: anna.klos@wat.edu.pl4. A next generation of synthetic benchmark is ongoing.5. New results of blind homogenization by the September/October.6. Next homogenization workshop following two previous ones in Brussels (in 2016) and Warsaw (in 2017). If interested in partcipating – please, contact us Motivation Data Synthetic benchmark   Tools   Results SummaryIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan18And then… 1. The best performing tools are going to be employed to homogenize the IGS repro1.2. A need to define the reliable strategy for homogenization.Motivation Data Synthetic benchmark   Tools   Results Summaryhttp://iag-gnssclimate.oma.be/index.phpIAG-IASPEI Joint Scientific Meeting, July 30 – August 4, 2017, Kobe, Japan19AcknowledgmentsAnna Klos was supported by COST Action ES1206 GNSS4SWEC (www.gnss4swec.knmi.nl) during her stay at the Royal Observatory of Belgium (ROB) and the Royal Meteorological Institute of Belgium (RMI). Janusz Bogusz was supported by the Polish National Science Centre, grant No. UMO-2016/21/B/ST10/02353.The map was drawn in the Generic Mapping Tool (Wessel et al., 2013).Thank you!Motivation Data Synthetic benchmark   Tools   Results SummaryView publication stats

Performance of various homogenization tools on a synthetic benchmark dataset of GPS and ERA-interim IWV differences

https://repositorio.aemet.es/bitstream/20.500.11765/11990/1/IAGHomogenizationAK.pdf

Performance of various homogenization tools on a synthetic benchmark dataset of GPS and ERA-interim IWV differences

Abstract

Similar works

Full text

Available Versions

Agencia Estatal de Meteorología