report of the wut lac · 2016-04-07 · 8th epn lac workshop, brussels, may 15–16 report of the...

8th EPN LAC Workshop, Brussels, May 15–16

Report of the WUT LAC

T. Liwosz, J. Rogowski

Warsaw University of TechnologyDepartment of Geodesy and Geodetic Astronomy

Pl. Politechniki 1, 00-661 Warsaw

8th EPN LAC WorkshopBrussels, May 15–16, 2013

Liwosz and Rogowski Report of the WUT LAC


Contents

This report contains the following topics:1 Routine analysis:

ProductsStrategy descriptionNetwork characteristics

2 Analysis of Poland’s national GPS reference network

3 GNSS analysis with GLONASS-specific PCVs

4 Outlook



WUT Products

Currently WUT contributes to EPN with the following products:Final (operational)

weekly coordinate estimates (wutWWWW7.snx)daily coordinate estimates (wutWWWWD.snx)daily tropospheric delays (wutWWWWD.tro)

Rapid (since week 1565, January 2010)daily coordinate estimates (wutWWWWDr.snx)rapid multi-day time series (since February 2011); dedicatedweb page created and maintained by EPN CB

Reprocessed results from a period 1996–2006



Processing strategy

Summary of current strategy:

Software: Bernese GPS Software version 5.0

processing performed according to EPN LAC Guidelines

GPS and GLONASS (since week 1609) observations analyzed

CODE orbits and ERPs

antenna calibration tables: epn08 wwww.atx

Ambiguity resolution for GPS satellites: L1&L2, L5/L3, QIFdepending on baseline length

Reference Frame: IGS08

Sites used for NNT minimum constraints: usable IGb08stations (currently 16)



WUT Subnetwork

−40˚ −20˚ 0˚ 20˚ 40˚30˚

40˚

50˚

60˚

70˚

GPS

GNSS

IGb08

KELY, NYA1 and THU3 outside the map

Network status:79 stations:

61 GNSS (77%)16 IGb085 proposed

Since 2011:

4 new

4 former



Poland’s national GNSS reference network

In the end of 2011, for the Polish NMA – Head Office for Geodesyand Cartography (GUGiK), we processed Poland’s national GNSSreference network and created a new national realization ofETRS89 for Poland.



Poland’s national GNSS reference network

15˚ 20˚ 25˚

50˚

55˚

POLREF

EUREF−POL

1st class

EUVN

ASG−PL/CZEPOS/LITPOS/SAPOS/SKPOS

EPN/IGS

Two campaigns:1 20082 2010/2011

66 days

∼700 sites processed



National network: Data analysis and results

Software: Bernese GPS Software version 5.0

Data processed with two elevation masks: 5◦ and 10◦

EPN cumulative solution (C1600) as reference (35 stations)

Statistics of combined solutions:

Solution type Repeatability (mm) RMS # observationsN E U (mm)

WUT ASG-PL 5◦ 1.20 1.23 3.11 1.29 38 617 779WUT ASG-PL 10◦ 1.38 1.30 3.47 1.32 33 532 898

Final coordinates (solution WUT ASG-PL 5◦) expressed in:

ITRF2005 at epoch 2011.0

ETRF2000(R05) at epoch 2011.0



GNSS analysis with GLONASS-specific PCVs

Since igs08.atx, IGS supports the receiver GLONASS specific PCC.Tests performed to check the effect on coordinates:

three solutions with same baseline definition:1 GNSS: GNSS obs. with GPS PCVs for both systems2 GNSSPCV: GNSS obs.; GPS and GLONASS specific PCVs

for GPS and GLONASS obs. respectivelyour own implementation into Bernese GPS Software ver. 5.0constant values, for k = 0, used for all other channels kGLONASS freqs. used for GLONASS L3 PCV

3 GPS: only GPS obervations used

igs08 1682.atx model (with spikes)

GPS-only IGS08 sites as reference

data span: April 17, 2011 – December 17, 2011 (245 days)



GNSS Analysis – test network

−20˚ 0˚ 20˚ 40˚

30˚

40˚

50˚

60˚

AJAC

BELF

BOR1

BORJ

BPDL

BRUS

BUCU

BUTE

BYDG

COMO

DAREDRES

DUBR

EGLT

ENIS

EUSK

FOYL

GANP

GLSV

GOPE

GWWL

HERT

HOE2

HOFN

IGMI

INVR

JOZE/2

KHAR

KLOP

KURE

LAMA

LODZ

MAN2

MARS

MATE

METS

NOT1

ONSA

PENC

POLV

POUS

REDZ

REYK

SASS

SKE0

SMLA

SVTL

SWKI

TORI

USDL

VEN1

VILL

WROC

WSRT

WTZRZIMM

ZYWI

Orange: GNSS site; Grey: GPS-only reference site



GNSS Analysis – examples of coordinate differences

−3

−2

−1

0

1

2

3

No

rth

(m

m)

GNSSPCV−GNSS

−3

−2

−1

0

1

2

3

Ea

st

(mm

)

−9

−6

−3

0

3

6

9

Up

(m

m)

1640 1650 1660

GPS week

−3

−2

−1

0

1

2

3

No

rth

(m

m)

GNSSPCV−GNSS

−3

−2

−1

0

1

2

3

Ea

st

(mm

)

−9

−6

−3

0

3

6

9

Up

(m

m)

1640 1650 1660

GPS week

−3

−2

−1

0

1

2

3

No

rth

(m

m)

GNSSPCV−GNSS

−3

−2

−1

0

1

2

3

Ea

st

(mm

)

−9

−6

−3

0

3

6

9

Up

(m

m)

1640 1650 1660

GPS week

DARE WTZR DRES



GNSS analysis – comparison

Mean differences over 245 days between GNSS solution andGNSSPCV solution and their standard deviations.

−2

−1

0

1

2

Nort

h (

mm

)

−2

−1

0

1

2

East (m

m)

−4

−2

0

2

4

Up (

mm

)

DARE

INVR

BORJ

DRES

HOE2

WTZR

TORI

BELF

BUCU

EGLT

ENIS

EUSK

FOYL

KURE

LAMA

MARS

PENC

WROC

AJAC

HERT

JOZ2

MAN2

MATE

VEN1

GLSV

KHAR

SMLA

COMO

HOFN

POUS

SASS

GOPE

IGMI

REYK

SVTL

DUBR

GANP

BPDL

BUTE

BYDG

GWWL

KLOP

LODZ

REDZ

SWKI

USDL

ZYWI

LEIAR25_LEIT

LEIAR25.R3_LEIT

LEIAR25.R3_NONE

LEIAT504GG_LEIS

LEIAT504GG_NONE

NOV702GG_NONE

TPSCR3_GGD_CONE

TPSCR.G3_TPSH

TRM55971.00_NONE

TRM55971.00_TZGD



GNSS analysis cont.

After introducing receiver GLONASS-specific PCVs in our GNSSanalysis:

no repeatability change of station positions

in general solutions did not get closer to GPS-only solutions

better correlation between height change and ZTD changewrt. GPS-only solution:

−6

−4

−2

0

2

4

6

∆ZT

D (

mm

)

−6 −4 −2 0 2 4 6

∆Height (mm)

GPS minus GNSS

R=−0.74

(a)

−6

−4

−2

0

2

4

6

∆ZT

D (

mm

)

−6 −4 −2 0 2 4 6

∆Height (mm)

GPS minus GNSSPCV

R=−0.88

(b)

Difference of correlation coefficients (R) significant (p < 0.05)



GNSS Analysis – difference of coordinates: REYK

−6

−4

−2

0

2

4

6

No

rth

(m

m)

GPS−GNSS

−6

−4

−2

0

2

4

6

GPS−GNSSPCV

−6

−4

−2

0

2

4

6

GNSSPCV−GNSS

−6

−4

−2

0

2

4

6

Ea

st

(mm

)

−6

−4

−2

0

2

4

6

−6

−4

−2

0

2

4

6

−12

−8

−4

0

4

8

12

Up

(m

m)

1640 1650 1660

GPS week

−12

−8

−4

0

4

8

12

1640 1650 1660

GPS week

−12

−8

−4

0

4

8

12

1640 1650 1660

GPS week



Outlook

Currently WUT contributes to the EPN with final and rapidproducts:

plans to prepare and submit NRT results

WUT can participate in repro2

We could possibly generate solutions with global sites (initialanalysis already done)we would use the new version of Bernese GNSS Software (5.2)for GNSS data analysis for EUREF

software available at WUT (∼2 weeks); no tests performed sofar


report of the wut lac · 2016-04-07 · 8th epn lac workshop, brussels, may 15–16 report of the...

Documents