[IGSREPORT-20733] Wk 1681 IGS SINEX Combination
igs-rf at ign.fr
igs-rf at ign.fr
Thu Apr 12 07:17:39 PDT 2012
Authors: A. Duret, B. Garayt, P. Rebischung
================================================================================
Combination of AC solutions for GPS week 1681 - 12/03/25 to 12/03/31
================================================================================
1) Contacts:
============
e-mail: igs-rf at ign.fr
Tel: 0033143988197 (Bruno Garayt)
Fax: 0033143988450
2) Products:
============
Products are available at:
CDDIS: ftp://cddis.gsfc.nasa.gov/gps/products
SOPAC: ftp://garner.ucsd.edu/pub/products
ENSG: ftp://igs.ensg.eu/pub/igs/products
IGN: ftp://igs.ign.fr/pub/igs/products
igs12P1681.snx Combined weekly solution
igs12P1681.erp Daily Earth Rotation Parameters
igs12P1681.sum This file
igs12P1681.itr Residuals between AC solutions and IGS08 for the RF stations
igs12P1681.res Residuals between AC solutions and combined weekly solution
3) Combination Strategy:
========================
The AC solutions are processed using the following procedure:
3-1) The AC solutions are checked/corrected for possible format problems,
station name inconsistencies, etc...
3-2) Constrained AC solutions are unconstrained according to the information
given in their SOLUTION/APRIORI blocks and then minimally constrained.
AC solutions given in the form of a normal matrix are inverted by adding
minimal constraints.
3-4) The weekly AC solutions are combined. The 7-parameter singularity is
avoided by using internal constraints. Iterations are performed until no
residuals exceed 2 cm or 5 sigma.
3-5) The combined solution is compared to IGS08 propagated at current week
(datum). Some outliers may be removed from the datum. Finally, the
combined solution is aligned to the datum using a 7-parameter
transformation estimated using the IGS08 core network.
3-6) The AC solutions NOT used in the combination are compared to the combined
solution in order to identify and remove outliers.
4) Remarks:
===========
mit = mit16817.snx
mig = mit1681g.snx
ITR = IGS08.snx
igs = igs12P1681.snx (weekly combined solution)
IGS = IGS12P13.snx (cumulative solution)
The following solutions are not used in the combination but only for
comparison:
jpl mig
Solutions with manual station deletion:
esa LPGS
gfz CRAR
mit CORB LPGS
ngs CHIN
sio CRAR CANT LPGS NEIA
jpl LPGS
mig NEIA SAMO TONG LAUT CRAR KATZ KABR ELAT TELA
5) Results:
============
5-1) AC (Variance Factor)^0.5 estimates:
=============================================
At each iteration, the variance factor for solution s is estimated by :
VF = (vs'*Ps*vs)/fs with fs = nobss - tr(As*Q*As'*Ps)
vs : residuals of solution s
Ps : weight matrix of solution s
nobss : number of observations of solution s
As : design matrix of solution s
Q : total covariance matrix
This variance factor is used to rescale the covariance matrix of solution s
for the next iteration.
Solution |# sta| #RF |#core| Final VF |
---------------|-----|-----|-----|---------------|
cod | 261 | 99 | 54 | 7.776 |
emr | 83 | 56 | 50 | 2.021 |
esa | 135 | 57 | 43 | 1.794 |
gfz | 191 | 99 | 55 | 6.440 |
grg | 123 | 77 | 52 | 2.890 |
mit | 269 | 102 | 53 | 0.951 |
ngs | 229 | 104 | 57 | 13.871 |
sio | 282 | 114 | 50 | 1.241 |
jpl | 149 | 79 | 51 | ------ |
mig | 526 | 136 | 58 | ------ |
igs | 448 | 135 | 58 | ------ |
---------------|-----|-----|-----|---------------|
5-2) Solution Residuals weighted Average & STD:
====================================================
5-2-1) With respect to the IGS08 stations:
=============================================
Solution | Weighted Average (mm) | STD (mm) |
| N E U | N E U |
---------------|---------------------------|---------------------------|
cod 99 | 0.1 -0.1 0.1 | 3.4 3.3 8.0 |
emr 56 | 0.7 -0.1 0.1 | 3.6 3.7 8.2 |
esa 57 | 0.6 -0.2 -0.2 | 4.6 4.0 7.8 |
gfz 99 | 0.2 0.0 -0.3 | 3.7 3.0 8.5 |
grg 77 | 0.4 -0.2 0.0 | 3.4 3.4 7.2 |
mit 102 | -0.3 0.1 -0.3 | 3.9 3.8 7.6 |
ngs 104 | 0.5 0.1 -0.1 | 3.5 3.4 7.9 |
sio 114 | 0.1 0.0 0.4 | 4.0 4.1 9.6 |
jpl 79 | 0.3 -0.1 0.1 | 3.8 3.0 9.3 |
mig 136 | -0.1 0.1 -0.3 | 3.4 3.4 7.3 |
---------------|---------------------------|---------------------------|
5-2-2) With respect to the Weekly Combination:
=================================================
Solution | Weighted Average (mm) | STD (mm) |
| N E U | N E U |
---------------|---------------------------|---------------------------|
cod 255 | -0.3 -0.1 0.0 | 1.0 0.7 3.0 |
emr 83 | 0.4 -0.1 -0.0 | 1.1 0.8 3.2 |
esa 133 | -0.2 0.0 0.1 | 1.0 1.1 3.1 |
gfz 189 | -0.2 0.1 -0.0 | 0.9 1.0 2.7 |
grg 121 | 0.4 -0.0 -0.1 | 1.1 1.0 4.0 |
mit 244 | 0.0 0.0 -0.2 | 0.9 0.8 2.3 |
ngs 223 | 0.5 0.0 -0.1 | 1.1 1.5 3.1 |
sio 264 | -0.4 -0.1 0.6 | 1.9 1.7 5.7 |
jpl 126 | 0.1 0.0 0.0 | 1.3 1.2 5.1 |
mig 439 | -0.2 0.0 -0.1 | 0.2 0.2 0.8 |
ITR 135 | -0.1 -0.1 -0.1 | 3.4 3.3 7.2 |
---------------|---------------------------|---------------------------|
5-2-3) With respect to the Cumulative Combination:
=====================================================
Solution | Weighted Average (mm) | STD (mm) |
| N E U | N E U |
---------------|---------------------------|---------------------------|
cod 232 | 0.0 -0.0 0.0 | 2.0 1.9 5.3 |
emr 82 | 0.5 -0.0 0.0 | 1.9 2.3 5.9 |
esa 130 | 0.2 0.1 0.2 | 2.3 2.4 6.0 |
gfz 188 | 0.2 0.2 0.0 | 2.4 2.2 5.2 |
grg 121 | 0.8 0.1 -0.0 | 2.1 2.3 6.1 |
mit 236 | 0.2 0.0 -0.2 | 2.1 2.4 5.6 |
ngs 218 | 0.8 0.1 -0.0 | 2.0 2.2 5.6 |
sio 254 | -0.1 -0.0 0.6 | 3.3 3.0 8.0 |
jpl 123 | 0.4 0.0 0.1 | 2.6 2.1 7.1 |
mig 402 | -0.0 0.1 -0.1 | 2.2 2.3 5.0 |
igs 408 | 0.1 0.0 -0.1 | 2.2 2.3 4.9 |
ITR 173 | -0.0 -0.1 0.0 | 2.4 2.7 4.8 |
---------------|---------------------------|---------------------------|
5-3) 7-Parameters Transformation: AC ----> IGS08
=====================================================
These transformation parameters are estimated using the IGS08 core stations
available in the AC solutions.
5-3-1) Estimates:
====================
Solution | R X R Y R Z T X T Y T Z SCL |
| (mas) (mas) (mas) (mm) (mm) (mm) (ppb) |
---------------|---------------------------------------------------------------|
cod | -0.056 0.000 -0.001 -0.4 -0.3 -1.0 0.17 |
emr | -0.004 0.008 0.013 2.9 -3.8 9.1 0.70 |
esa | 0.035 -0.137 0.127 3.3 0.9 0.3 0.42 |
gfz | -0.045 0.027 -0.027 -1.0 0.4 -0.7 0.23 |
grg | -0.059 0.043 -0.014 -1.1 -0.9 0.4 0.32 |
mit | -0.047 0.031 -0.051 6.1 -5.1 11.7 0.50 |
ngs | -0.065 0.007 0.014 3.2 -2.1 13.8 0.58 |
sio | -0.024 -0.054 0.125 -2.4 -4.1 -5.1 0.85 |
jpl | -0.005 -0.028 0.026 -5.1 -0.9 9.5 0.49 |
mig | -0.029 0.022 -0.046 -0.7 -0.5 0.3 0.44 |
igs | 0.000 0.000 0.000 0.0 0.0 0.0 0.46 |
---------------|---------------------------------------------------------------|
5-3-2) STD:
==============
Solution | R X R Y R Z T X T Y T Z SCL |
| (mas) (mas) (mas) (mm) (mm) (mm) (ppb) |
---------------|---------------------------------------------------------------|
cod | 0.032 0.035 0.036 0.9 0.9 0.9 0.13 |
emr | 0.031 0.034 0.036 0.8 0.8 0.8 0.13 |
esa | 0.034 0.038 0.038 0.9 0.9 0.9 0.14 |
gfz | 0.031 0.034 0.035 0.8 0.8 0.8 0.13 |
grg | 0.028 0.030 0.032 0.8 0.8 0.8 0.12 |
mit | 0.032 0.034 0.035 0.8 0.8 0.8 0.13 |
ngs | 0.029 0.031 0.033 0.8 0.8 0.8 0.12 |
sio | 0.041 0.044 0.046 1.1 1.1 1.1 0.17 |
jpl | 0.032 0.035 0.038 0.9 0.9 0.9 0.14 |
mig | 0.028 0.030 0.032 0.8 0.8 0.8 0.12 |
igs | 0.000 0.000 0.000 0.0 0.0 0.0 0.22 |
---------------|---------------------------------------------------------------|
5-4) Geocenter:
====================
| X Y Z | s X s Y s Z |
| (mm) (mm) (mm) | (mm) (mm) (mm) |
---------------|---------------------------|---------------------------|
igs | 1.7 -1.3 5.1 | 0.4 0.4 0.5 |
---------------|---------------------------|---------------------------|
5-5) ERP Residuals Weighted Average & STD:
===============================================
5-5-1) Weighted Average:
===========================
Solution | UT LOD XPO XPOR YPO YPOR |
| us us .01mas .01ma/d .01mas .01ma/d |
---------------|------------------------------------------------------------|
cod | --- --- --- --- --- --- |
emr | --- -6.2 2.1 -11.4 2.8 20.8 |
esa | --- 18.7 -1.2 -13.4 -0.4 -8.5 |
gfz | --- -2.1 -0.9 11.8 -1.6 -5.6 |
grg | --- --- 1.5 3.9 0.5 -1.9 |
mit | --- 6.3 1.6 -8.9 -0.8 4.9 |
ngs | --- -0.1 -3.6 13.4 -0.1 -1.8 |
sio | --- -11.8 0.0 -27.8 4.8 -2.9 |
jpl | --- 4.4 0.9 10.5 1.2 7.4 |
mig | --- 46.1 0.7 -15.4 2.9 -4.0 |
BuA | --- 4.8 -4.2 8.2 -0.4 8.3 |
---------------|------------------------------------------------------------|
5-5-2) STD:
==============
Solution | UT LOD XPO XPOR YPO YPOR |
| us us .01mas .01ma/d .01mas .01ma/d |
---------------|------------------------------------------------------------|
cod | --- --- --- --- --- --- |
emr | --- 4.3 2.3 9.2 1.5 18.1 |
esa | --- 3.4 1.0 6.7 1.9 5.5 |
gfz | --- 6.2 1.5 9.7 2.2 13.0 |
grg | --- --- 1.6 5.9 1.8 5.7 |
mit | --- 4.8 1.2 7.9 1.0 3.5 |
ngs | --- 0.1 4.3 6.1 4.8 9.5 |
sio | --- 3.2 2.3 6.8 1.8 5.4 |
jpl | --- 6.8 1.1 9.8 2.9 5.4 |
mig | --- 2.9 2.5 7.0 1.7 4.1 |
BuA | --- 10.1 4.7 9.0 5.5 10.5 |
---------------|------------------------------------------------------------|
5-6) Satellite Antenna Offset Residuals weighted Average & STD:
====================================================================
No Satellite phase center ...
5-7) Outliers:
===================
Outliers may be rejected at different stages of the processing:
- During the combination (stage 3-4)), outliers may be removed from the
weekly AC solutions. Thresholds are set to 2 cm and 5 sigma.
- During stage 3-6), outliers may be removed from the AC solutions NOT
included in the combination. Thresholds are set to 2 cm and 5 sigma.
- During the stacking of the weekly combined solution on the cumulative
solution (stage 3-7)), outliers may be removed from the weekly
solution. Thresholds are set to 3 cm and 8 sigma.
Solution Stage of rejection | Station | N E U |
|code id soln| (mm) (mm) (mm)|
-----------------------------------------|------------|-----------------------|
| | |
grg16817.snx Combination: Iter. 1 |NOVM A 1| 42.4 -2.3 35.2|
grg16817.snx Combination: Iter. 1 |POTS A 3| 7.8 3.0 -3.0|
sio16817.snx Combination: Iter. 1 |SANT A 5| 2.5 -0.6 21.5|
sio16817.snx Combination: Iter. 2 |GUAO A 1| 1.5 0.5 14.7|
sio16817.snx Combination: Iter. 3 |PBRI A 1| 5.9 -0.6 0.1|
| | |
mit1681g.snx Comp. to igs12P1681.snx |NOVM A 1| 0.1 5.9 3.8|
| | |
igs12P1681.snx Stacking: Iter. 1 |ALRT A 1| 4.6 8.6 8.8|
igs12P1681.snx Stacking: Iter. 1 |CTWN A 2| 1.6 -2.6 49.3|
igs12P1681.snx Stacking: Iter. 1 |DARW A 3| 7.9 5.0 -11.0|
igs12P1681.snx Stacking: Iter. 1 |SMST A 2| -5.7 4.4 85.0|
igs12P1681.snx Stacking: Iter. 1 |TWTF A 2| 8.7 -0.1 -3.7|
| | |
-----------------------------------------|------------|-----------------------|
6) Check of metadata in AC solutions:
=====================================
Starting with week 1585, receivers, antennas and eccentricities reported in
AC solutions are checked versus station logs. Station logs are gathered
from the following sources:
igs ftp://igscb.jpl.nasa.gov/pub/station/log
igsf ftp://igscb.jpl.nasa.gov/pub/station/oldlog
epn ftp://epncb.oma.be/pub/station/log
rgp ftp://rgpdata.ign.fr/pub/logsheet
cors ftp://www.ngs.noaa.gov/cors/station_log
argn ftp://ftp.ga.gov.au/geodesy-outgoing/gnss/logs
nrcan ftp://rtopsdata.geod.nrcan.gc.ca/pub/gps/stations_log
sirgas ftp://ftp.sirgas.org/pub/gps/DGF/station/log
rbmc ftp://geoftp.ibge.gov.br/RBMC/dados/log
unavco ftp://data-out.unavco.org/pub/logs
ucsd ftp://garner.ucsd.edu/pub/docs/site_logs
ign ftp://igs.ign.fr/pub/igs/igscb/station/log
6-1) Receivers:
====================
AC | Station | Used Receiver | Receiver from log | Source |
----|---------|----------------------|----------------------|--------|
esa | AREQ A | JAVAD TRE_G3TH | ASHTECH UZ-12 | igs |
gfz | BIS2 A | AOA BENCHMARK ACT | ASHTECH Z-XII3 | cors |
sio | DGAV A | JPS EGGDT | JAVAD TRE_G3TH DELTA | igs |
mit | ELRO A | TPS NET-G3A | LEICA RS500 | ucsd |
sio | MANA A | TRIMBLE 4000SSI | TRIMBLE NETR9 | igs |
mit | NETT A | TRIMBLE NETRS | TRIMBLE 4000SSI | ucsd |
sio | OHI3 A | LEICA GRX1200GGPRO | LEICA GRX1200+GNSS | igs |
mit | SIO5 A | TPS NET-G3A | ASHTECH Z-XII3 | cors |
sio | SPT0 A | JPS E_GGD | JAVAD TRE_G3TH DELTA | igs |
cod | TABV A | JAVAD TRE_G3TH DELTA | JPS EGGDT | jpl |
mit | YRCM A | LEICA SR530 | LEICA RS500 | ucsd |
mit | ZAMB A | ROGUE SNR-8000 | ASHTECH Z-XII3 | igs |
sio | ZAMB A | ROGUE SNR-8000 | ASHTECH Z-XII3 | igs |
----|---------|----------------------|----------------------|--------|
6-2) Antennas:
===================
AC | Station | Used Antenna | Antenna from log | Source |
----|---------|----------------------|----------------------|--------|
ngs | AIRA A | TRM29659.00 NONE | TRM29659.00 DOME | igs |
ngs | ALBH A | AOAD/M_T NONE | AOAD/M_T SCIS | igs |
mit | ALON A | LEIAT504 SNOW | LEIAT504 LEIS | ucsd |
grg | AREQ A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
mit | AREQ A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
ngs | AREQ A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
sio | AREQ A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
grg | ARTU A | ASH700936D_M NONE | ASH700936D_M DOME | igs |
ngs | ARTU A | ASH700936D_M NONE | ASH700936D_M DOME | igs |
gfz | BIS2 A | AOAD/M_T NONE | ASH700829.3 SNOW | cors |
ngs | BLYT A | ASH701945B_M NONE | ASH701945B_M ENCL | igs |
grg | CAS1 A | AOAD/M_T NONE | AOAD/M_T AUST | igs |
ngs | CAS1 A | AOAD/M_T NONE | AOAD/M_T AUST | igs |
grg | CEDU A | AOAD/M_T NONE | AOAD/M_T AUST | igs |
ngs | CEDU A | AOAD/M_T NONE | AOAD/M_T AUST | igs |
ngs | CLAR A | ASH700936A_M NONE | ASH700936A_M SCPL | igs |
mit | CONO A | TRM33429.20+GP UNKN | TRM33429.20+GP NONE | cors |
mit | COVX A | TRM55971.00 NONE | TRM57971.00 NONE | cors |
grg | CRO1 A | ASH701945G_M NONE | ASH701945G_M JPLA | igs |
ngs | CRO1 A | ASH701945G_M NONE | ASH701945G_M JPLA | igs |
mit | CSAR A | TRM29659.00 SNOW | TRM29659.00 NONE | ucsd |
grg | FAIR A | ASH701945G_M NONE | ASH701945G_M JPLA | igs |
ngs | FAIR A | ASH701945G_M NONE | ASH701945G_M JPLA | igs |
ngs | GODE A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
ngs | GUAM A | ASH701945B_M NONE | ASH701945B_M JPLA | igs |
mit | HILB A | TRM33429.00+GP UNKN | TRM33429.00+GP NONE | cors |
grg | KIR0 A | AOAD/M_T NONE | AOAD/M_T OSOD | igs |
grg | KOUC A | TRM57971.00 NONE | TRM57971.00 TZGD | igs |
ngs | KOUC A | TRM57971.00 NONE | TRM57971.00 TZGD | igs |
grg | LHAZ A | ASH701941.B NONE | ASH701941.B SNOW | igs |
ngs | LHAZ A | ASH701941.B NONE | ASH701941.B SNOW | igs |
ngs | MAC1 A | AOAD/M_T NONE | AOAD/M_T AUST | igs |
ngs | MAG0 A | ASH700936D_M NONE | ASH700936D_M JPLA | igs |
grg | MAR6 A | AOAD/M_T NONE | AOAD/M_T OSOD | igs |
grg | MAW1 A | AOAD/M_T NONE | AOAD/M_T AUST | igs |
ngs | MAW1 A | AOAD/M_T NONE | AOAD/M_T AUST | igs |
grg | MCIL A | TRM29659.00 NONE | TRM29659.00 DOME | igs |
grg | MCM4 A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
ngs | MCM4 A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
grg | MDO1 A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
ngs | MDO1 A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
grg | MTKA A | ASH701073.1 NONE | ASH701073.1 SNOW | igs |
mit | NETT A | TRM55971.00 NONE | TRM29659.00 SCIS | ucsd |
ngs | NVSK A | TRM33429.00+GP NONE | TRM33429.00+GP DOME | igs |
ngs | NYA1 A | ASH701073.1 NONE | ASH701073.1 SNOW | igs |
grg | NYAL A | AOAD/M_B NONE | AOAD/M_B DOME | igs |
ngs | NYAL A | AOAD/M_B NONE | AOAD/M_B DOME | igs |
grg | ONSA A | AOAD/M_B NONE | AOAD/M_B OSOD | igs |
ngs | ONSA A | AOAD/M_B NONE | AOAD/M_B OSOD | igs |
ngs | PETS A | ASH701933B_M NONE | ASH701933B_M SCIS | igs |
sio | RWSN A | ASH700936D_M UNKN | ASH700936D_M NONE | sirgas |
grg | SANT A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
ngs | SANT A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
ngs | SHAO A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
mit | SNFD A | TRM33429.20+GP UNKN | TRM33429.20+GP NONE | cors |
grg | SPT0 A | AOAD/M_T NONE | AOAD/M_T OSOD | igs |
ngs | SUTB A | ASH700936C_M NONE | ASH700936C_M SCIS | cors |
grg | SYOG A | AOAD/M_T NONE | AOAD/M_T DOME | igs |
ngs | SYOG A | AOAD/M_T NONE | AOAD/M_T DOME | igs |
cod | TABV A | TPSCR.G3 SCIS | ASH700936A_M SNOW | jpl |
grg | THU3 A | ASH701073.1 NONE | ASH701073.1 SCIS | igs |
ngs | THU3 A | ASH701073.1 NONE | ASH701073.1 SCIS | igs |
ngs | TIDB A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
sio | VBCA A | LEIAT303 UNKN | LEIAT303 NONE | sirgas |
grg | VIS0 A | AOAD/M_T NONE | AOAD/M_T OSOD | igs |
ngs | WAB2 A | ASH700936F_C NONE | ASH700936F_C SNOW | igs |
ngs | WILL A | AOAD/M_T NONE | AOAD/M_T EMRA | igs |
grg | WSRT A | AOAD/M_T NONE | AOAD/M_T DUTD | igs |
ngs | WSRT A | AOAD/M_T NONE | AOAD/M_T DUTD | igs |
ngs | YAKT A | ASH701933B_M NONE | ASH701933B_M SCIS | igs |
ngs | YAR2 A | AOAD/M_T NONE | AOAD/M_T JPLA | igs |
mit | YRCM A | LEIAT504 SNOW | LEIAT504 LEIS | ucsd |
ngs | YSSK A | ASH701933B_M NONE | ASH701933B_M DOME | igs |
----|---------|----------------------|----------------------|--------|
6-3) Eccentricities:
=========================
AC | Station | Used Eccentricity | Eccentricity from log | Source |
----|---------|-------------------------|-------------------------|--------|
----|---------|-------------------------|-------------------------|--------|
6-4) No station logs found for:
====================================
CRAR A 66001M004
CTWN A M
MCCJ A M
OBE3 A 14208M005
TIXG A 12360M001
UNB3 A 40146M002
More information about the IGSREPORT
mailing list