[IGSREPORT-18568] Wk 1581 IGS SINEX Combination
igs-rf at ign.fr
igs-rf at ign.fr
Wed May 12 10:05:14 PDT 2010
******************************************************************************
IGS Electronic Report 12 May 10:05:38 PDT 2010 Message Number 18568
******************************************************************************
Authors: B. Garayt, P. Rebischung
===============================================================
Combined AC solution for GPS week 1581 - 10/04/25 to 10/05/01
===============================================================
1) Contacts:
============
e-mail: igs-rf at ign.fr
Tel: 0033143988197 (Bruno Garayt)
Fax: 0033143988450
2) Products:
============
All products are available at:
CDDIS: ftp://cddis.gsfc.nasa.gov/gps/products
SOPAC: ftp://garner.ucsd.edu/pub/products
IGN: ftp://igs.ensg.eu/pub/igs/products
ftp://igs.ign.fr/pub/igs/products
igs10P1581.snx Combined weekly solution
igs10P1581.erp Daily Earth Rotation Parameters
igs10P1581.sum This file
igs10P1581.itr Residuals between AC solutions and IGS05 for the RF stations
igs10P1581.res Residuals between AC solutions and combined weekly solution
IGS10P17.snx Cumulative combined solution
IGS10P17.res Residuals between AC solutions and cumulative solution
igs00p03.erp Daily Earth Rotation Parameters since GPS week 860
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 names 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 (esa) are inverted by
adding minimal constraints.
3-4) The weekly AC solutions are combined. The 7-parameters singularity is
avoided by using internal constraints. Iterations are performed until no
residuals exceed 5 cm or 5 sigma.
3-5) The combined solution is compared to IGS05 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-parameters
transformation.
3-6) The AC solutions NOT used in the combination are compared to the combined
solution in order to identify and remove outliers.
3-7) The weekly combined solution is added to the cumulative solution.
Iterations are performed until no residuals exceed 3 cm or 3 sigma in the
weekly solution.
3-8) The updated cumulative solution is compared to IGS05 (datum). Some outliers
may be removed from the datum. Finally, the cumulative solution is aligned
to the datum using a 14-parameters transformation.
4) Remarks:
===========
mit = mit15817.snx
mig = mit1581g.snx
ITR = IGS05.snx
igs = igs10P1581.snx
IGS = IGS10P17.snx
grg, mig, ncl are not used in the combination.
They are only used for comparison.
Solutions with manual station deletion:
mit SOFI
ngs BHR1 OSN1 SOFI
sio SOFI
grg SOFI
4-1) Starting with GPS week 1400, the solutions are aligned to the
IGS05.snx realization of ITRF2005.
4-2) Starting with GPS week 1078, UT1-UTC from IERS Bulletin A
is included in the erp files.
4-3) Starting with GPS week 1412, the accumulated ERP file has been
changed from "igs00p02.erp" to "igs00p03.erp". The main changes
are:
1) Correction of the pole position from IGb00 to IGS05 using the
transformation in:
http://igscb.jpl.nasa.gov/mail/igsmail/2006/msg00178.html
2) Marginal ERP's in GPS week 1161 were replaced by reprocessed
values.
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.
Another scale factor is computed for each solution s, by comparing it to
the official cumulative solution of previous week (IGS) :
SF = (vs'*Ps*vs)/(3*nstas) with vs : residuals of the comparison
Ps : weight matrix of solution s
nstas : number of common stations
Solution |# sta| #RF | Final VF | SF |
---------------|-----|-----|---------------|---------------|
cod | 246 | 66 | 6.721 | 24.147 |
emr | 89 | 50 | 2.541 | 5.935 |
esa | 149 | 44 | 2.751 | 4.738 |
gfz | 200 | 66 | 3.790 | 15.993 |
jpl | 124 | 44 | 5.042 | 10.566 |
mit | 259 | 71 | 0.627 | 2.700 |
ngs | 223 | 71 | 14.817 | 39.754 |
sio | 269 | 61 | 1.133 | 2.544 |
grg | 128 | 50 | ------ | 8.862 |
mig | 429 | 70 | ------ | 1.719 |
ncl | 230 | 71 | ------ | 3.076 |
igs | 397 | 71 | ------ | ------ |
---------------|-----|-----|---------------|---------------|
5-2) Solution Residuals weighted Average & STD:
====================================================
5-2-1) With respect to the IGS05 stations:
=============================================
Solution | Weighted Average (mm) | STD (mm) |
| N E U | N E U |
---------------|---------------------------|---------------------------|
cod 66 | 0.0 0.0 0.2 | 3.0 3.4 8.8 |
emr 50 | 0.1 0.3 0.1 | 3.2 3.0 8.6 |
esa 44 | -0.1 -0.0 -0.4 | 3.2 3.5 10.3 |
gfz 66 | 0.1 -0.1 0.0 | 2.8 3.3 9.5 |
jpl 44 | 0.3 -0.3 -0.2 | 2.6 3.1 8.4 |
mit 71 | -0.1 0.1 0.0 | 2.8 3.5 9.1 |
ngs 71 | -0.0 0.0 -0.2 | 2.7 4.0 9.5 |
sio 61 | 0.3 0.1 -0.1 | 2.7 2.9 7.0 |
grg 50 | 1.0 -0.1 -0.2 | 3.2 3.6 8.1 |
mig 70 | 0.1 0.1 -0.1 | 2.6 3.2 8.1 |
ncl 71 | -0.1 0.1 -0.2 | 2.7 3.4 8.6 |
---------------|---------------------------|---------------------------|
5-2-2) With respect to the Weekly Combination:
=================================================
Solution | Weighted Average (mm) | STD (mm) |
| N E U | N E U |
---------------|---------------------------|---------------------------|
cod 232 | -0.0 -0.1 0.0 | 1.0 0.9 2.5 |
emr 89 | 0.5 0.1 -0.1 | 1.8 0.8 4.3 |
esa 145 | 0.3 -0.0 0.2 | 2.2 1.4 4.9 |
gfz 200 | 0.1 -0.1 0.0 | 0.5 0.6 2.0 |
jpl 124 | 0.4 -0.0 -0.2 | 1.8 2.4 5.0 |
mit 245 | -0.0 0.0 -0.1 | 0.5 0.6 1.6 |
ngs 219 | -0.1 0.1 0.0 | 1.0 1.6 3.1 |
sio 244 | 0.0 0.1 0.0 | 1.4 1.2 2.8 |
grg 123 | 0.6 -0.2 -0.2 | 3.0 1.8 4.9 |
mig 393 | -0.1 0.0 -0.0 | 0.3 0.2 0.7 |
ncl 229 | 0.0 0.0 0.0 | 0.6 0.7 0.7 |
ITR 71 | 0.0 -0.2 -0.1 | 2.7 3.2 8.6 |
---------------|---------------------------|---------------------------|
5-2-3) With respect to the Cumulative Combination:
=====================================================
Solution | Weighted Average (mm) | STD (mm) |
| N E U | N E U |
---------------|---------------------------|---------------------------|
cod 214 | -0.2 -0.2 -0.1 | 1.9 2.1 7.4 |
emr 88 | -0.1 0.1 -0.1 | 2.3 1.9 6.7 |
esa 136 | -0.3 -0.1 0.1 | 2.7 2.3 6.7 |
gfz 197 | -0.1 -0.2 -0.1 | 1.6 1.9 8.3 |
jpl 114 | -0.1 -0.1 -0.2 | 2.6 3.3 8.0 |
mit 234 | -0.5 -0.2 -0.0 | 1.7 2.2 5.0 |
ngs 209 | -0.2 -0.1 -0.3 | 1.8 2.7 8.6 |
sio 230 | -0.3 -0.1 0.1 | 2.1 1.9 4.7 |
grg 122 | -0.0 -0.3 -0.1 | 3.0 2.4 6.1 |
mig 344 | -0.4 -0.1 -0.2 | 1.8 2.1 6.6 |
ncl 224 | -0.3 -0.1 -0.4 | 1.5 1.9 7.1 |
igs 346 | -0.4 -0.1 -0.2 | 1.8 2.2 6.8 |
ITR 116 | -0.8 -0.1 -0.1 | 3.0 3.3 12.1 |
---------------|---------------------------|---------------------------|
5-3) 7-Parameters Transformation: AC ----> IGS05.snx
=========================================================
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.010 -0.002 -0.006 0.2 0.5 0.4 0.00 |
emr | 0.005 0.005 -0.012 1.4 0.5 0.4 0.21 |
esa | 0.047 0.054 -0.058 1.8 2.8 2.7 0.08 |
gfz | -0.001 0.013 -0.046 -0.1 1.2 0.5 -0.02 |
jpl | 0.031 -0.033 0.004 1.8 1.6 0.7 0.43 |
mit | 0.022 -0.002 -0.017 2.8 -2.1 -4.8 0.15 |
ngs | -0.003 -0.064 0.010 4.6 -0.1 2.9 0.09 |
sio | 0.028 -0.003 -0.009 0.5 1.5 0.4 0.20 |
grg | 0.004 0.003 0.006 -2.0 3.3 7.7 0.45 |
mig | 0.024 0.008 -0.006 -0.2 2.2 0.2 0.15 |
ncl | 0.031 0.187 -0.119 -4.1 1.2 0.0 0.13 |
igs | 0.000 0.000 0.000 0.0 0.0 0.0 0.14 |
---------------|---------------------------------------------------------------|
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.030 0.030 0.033 0.8 0.8 0.8 0.12 |
emr | 0.030 0.033 0.036 0.8 0.8 0.8 0.13 |
esa | 0.040 0.043 0.046 1.1 1.1 1.1 0.17 |
gfz | 0.032 0.032 0.035 0.8 0.8 0.8 0.13 |
jpl | 0.033 0.035 0.039 0.9 0.9 0.9 0.14 |
mit | 0.030 0.030 0.032 0.8 0.8 0.8 0.12 |
ngs | 0.031 0.032 0.034 0.8 0.8 0.8 0.13 |
sio | 0.024 0.025 0.026 0.6 0.6 0.6 0.10 |
grg | 0.033 0.034 0.036 0.9 0.9 0.9 0.13 |
mig | 0.027 0.028 0.030 0.7 0.7 0.7 0.11 |
ncl | 0.028 0.029 0.031 0.8 0.7 0.7 0.11 |
igs | 0.000 0.000 0.000 0.0 0.0 0.0 0.13 |
---------------|---------------------------------------------------------------|
5-4) Geocenter:
====================
| X Y Z | s X s Y s Z |
| (mm) (mm) (mm) | (mm) (mm) (mm) |
---------------|---------------------------|---------------------------|
igs | 3.2 -0.1 0.1 | 1.5 1.2 3.2 |
---------------|---------------------------|---------------------------|
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 | --- -0.5 4.4 0.9 1.4 30.9 |
esa | --- 12.6 -1.5 -5.9 -0.6 3.0 |
gfz | --- -9.6 1.1 -4.1 -0.3 -2.5 |
jpl | --- -9.1 3.7 -3.1 -1.2 11.7 |
mit | --- 2.9 -2.4 12.6 -0.5 -7.0 |
ngs | --- -10.0 0.8 -20.8 -3.2 0.9 |
sio | --- 11.4 0.9 -46.4 5.8 8.1 |
grg | --- -13.0 5.2 17.3 -10.5 13.1 |
mig | --- 14.2 2.0 -36.0 4.4 8.1 |
ncl | --- --- --- --- --- --- |
BuA | --- -6.2 -1.6 -5.3 4.1 1.9 |
---------------|------------------------------------------------------------|
5-5-2) STD:
==============
Solution | UT LOD XPO XPOR YPO YPOR |
| us us .01mas .01ma/d .01mas .01ma/d |
---------------|------------------------------------------------------------|
cod | --- --- --- --- --- --- |
emr | --- 3.5 3.0 11.7 3.3 11.6 |
esa | --- 4.1 1.3 7.7 2.4 8.4 |
gfz | --- 6.4 1.9 4.0 2.0 5.9 |
jpl | --- 15.4 1.5 7.7 2.2 15.9 |
mit | --- 2.5 0.9 3.5 1.1 5.6 |
ngs | --- 6.1 5.5 11.3 7.7 13.1 |
sio | --- 2.8 1.3 6.4 0.9 4.0 |
grg | --- 19.1 5.2 17.0 3.7 15.4 |
mig | --- 3.0 1.3 6.8 0.9 3.1 |
ncl | --- --- --- --- --- --- |
BuA | --- 14.7 4.8 10.1 4.2 9.8 |
---------------|------------------------------------------------------------|
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 5 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 5 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 3 sigma.
Solution Stage of rejection | Station | N E U |
|code id soln| (mm) (mm) (mm)|
-----------------------------------------|------------|-----------------------|
| | |
sio15817.snx Combination: Iter. 1 |SFER A 5| -1.3 -4.6 3.5|
sio15817.snx Combination: Iter. 1 |OHI3 A 2| 3.4 -0.1 -0.7|
sio15817.snx Combination: Iter. 1 |KOKB A 5| -1.5 -3.5 -5.4|
sio15817.snx Combination: Iter. 1 |CHAT A 1| 4.5 -3.0 1.0|
sio15817.snx Combination: Iter. 1 |WGTT A 1| 4.6 -2.4 1.2|
sio15817.snx Combination: Iter. 1 |OUS2 A 3| 4.4 -2.5 1.3|
sio15817.snx Combination: Iter. 1 |MCM4 A 3| 5.2 -1.5 -9.6|
sio15817.snx Combination: Iter. 1 |MAC1 A 7| 7.2 -1.2 -2.6|
sio15817.snx Combination: Iter. 1 |STR1 A 1| 5.3 -0.7 4.6|
sio15817.snx Combination: Iter. 1 |TIDB A 2| 5.8 -0.1 4.8|
sio15817.snx Combination: Iter. 1 |HOB2 A 2| 7.7 -0.5 7.0|
sio15817.snx Combination: Iter. 1 |CEDU A 1| 4.5 1.4 3.7|
sio15817.snx Combination: Iter. 1 |PERT A 5| 5.6 2.3 -1.7|
sio15817.snx Combination: Iter. 1 |CAS1 A 3| 6.2 2.4 -7.2|
sio15817.snx Combination: Iter. 1 |POL2 A 1| -0.4 3.7 -1.6|
sio15817.snx Combination: Iter. 1 |SEY1 A 3| 3.7 8.2 -2.5|
sio15817.snx Combination: Iter. 2 |STR2 A 1| 5.6 -1.5 4.9|
sio15817.snx Combination: Iter. 2 |PARK A 1| 6.4 -1.0 3.1|
sio15817.snx Combination: Iter. 2 |NNOR A 1| 5.8 2.9 4.4|
sio15817.snx Combination: Iter. 3 |DAV1 A 2| 3.9 4.2 -8.3|
sio15817.snx Combination: Iter. 4 |AUCK A 2| 5.0 -2.8 1.1|
sio15817.snx Combination: Iter. 5 |MAW1 A 1| 3.5 4.7 -7.5|
| | |
grg15817.snx Comp. to igs10P1581.snx |CONZ A 3| -3.4 27.8 5.1|
grg15817.snx Comp. to igs10P1581.snx |AREQ A 6| -15.0 -26.4 -6.2|
grg15817.snx Comp. to igs10P1581.snx |KIT3 A 2| -6.6 -4.7 20.7|
mit1581g.snx Comp. to igs10P1581.snx |SOFI A 1| -1.3 -0.7 78.6|
mit1581g.snx Comp. to igs10P1581.snx |MCM4 A 3| 0.9 -0.2 -5.2|
ncl1581g.snx Comp. to igs10P1581.snx |SOFI A 1| -1.6 -0.5 123.4|
| | |
igs10P1581.snx Stacking: Iter. 1 |BOGT A 4| 0.7 -9.5 30.6|
igs10P1581.snx Stacking: Iter. 1 |DUBO A 4| -12.4 -1.8 34.4|
igs10P1581.snx Stacking: Iter. 1 |NETP A 1| -3.8 -3.7 33.0|
igs10P1581.snx Stacking: Iter. 1 |QUI1 A 1| 6.5 -9.8 -56.3|
igs10P1581.snx Stacking: Iter. 1 |SSIA A 4| 3.5 -4.5 30.9|
igs10P1581.snx Stacking: Iter. 1 |XIAN A 1| 1.4 28.1 -33.6|
| | |
-----------------------------------------|------------|-----------------------|
More information about the IGSREPORT
mailing list