[IGSREPORT-14388] Wk 1396 IGS SINEX Combination
Martin Bourassa
mbourass at nrcan.gc.ca
Thu Oct 26 16:27:43 PDT 2006
******************************************************************************
IGS Electronic Report 26 Oct 16:27:48 PDT 2006 Message Number 14388
******************************************************************************
Author: R. Ferland , M. Bourassa
===================================================================
Combined AC solution for GPS week 1396 - 06/10/08 to 06/10/14
===================================================================
1) Contacts:
============
e-mail: RFerland at NRCan.gc.ca
Tel: (613)995-4002 (Remi Ferland)
Fax: (613)995-3215
2) Products:
============
All products are available at
CDDIS (cddis.gsfc.nasa.gov) and SOPAC (garner.ucsd.edu)
igs06P1396.snx Combined weekly solution.
igs06P1396.erp Daily Earth Rotation Parameters.
igs00p02.erp Daily Earth Rotation Parameters since GPS week 1013.
IGS06P41.snx Cumulative Combined solution.
igs06P1396.sum This file.
igs06P1396.itr Residuals between the AC and the IGb00 for the RF stations.
igs06P1396.res Residuals between the AC and the Combined Weekly solution.
IGS06P41.res Residuals between the AC and the Cumulative Combined solution.
Where IGSyyPww & igsyyPwwww represent:
yy Last 2 digits of the year
ww week in year yy
wwww GPS week
P GPS technique used
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) The AC solutions are unconstrained according to the information provided
in the SOLUTION/APRIORI Block.
3-3) The AC solutions are aligned (using a 7-parameters) and compared to the
IGb00 stations at the current epoch.
ftp://igscb.jpl.nasa.gov/igscb/station/coord/IGS03P33_RS99.SNX.Z
3-4) The AC are compared to IGb00, to each other, to the previous week
solution and to the cumulative solution in order to identify and remove
outliers.
3-5) The weekly AC solutions are combined and added to the cumulative
combined solution. The cumulative combination is done to ensures 14
parameters singularity.
The weekly combined solution is aligned to the IGb00 of the current
epoch using a 7-parameters transformation. The singular cumulative
combined solution is aligned to the IGb00 of the reference epoch
using a 14-parameters transformation and the 99 ITRF station positions
and velocities. The IGb00 is a realization of ITRF2000.
For more details, see:
Blewitt, G., C. Boucher and J. Kouba, 1997, Procedure for ITRF continuous
realization, the report of the IGS GB, presented to the IGS Governing
Board Meeting held in San Francisco, Dec. 11.
Kouba, J., J. Ray and M.M. Watkins, IGS Reference Frame Realization, 1998
IGS Analysis Center Workshop Proceedings, European Space Operations
Centre, Darmstadt, Germany.
The combination software was developed at NRCan.
4) Remarks:
===========
ITR = IGS03P33_RS99.SNX
igs = igs06P1396.snx
IGS = IGS06P41.snx
MIT = mit13967.snx
mit = mit1396g.snx
BuA = IERS Bulletin A
mit ncl are used to quality control the weekly combination.
gfz ngs excluded from the geocenter estimation (Sec. 5-4).
Solutions with manual station deletion:
After unconstraining:
esa CEBR NNOR
mit CEBR NNOR
4-1) Starting with GPS week 1143, the solutions are aligned to the
IGS03P33_RS99.SNX realization of ITRF2000.
4-2) Starting with GPS week 1078, UT1-UTC from IERS Bulletin A
is included in the erp files.
5) Results:
===========
5-1) AC (Variance Factor)^0.5 estimates:
==================================================
The residuals used in the estimation of the Variance Factor
are the differences between the cumulative combined solution
and each AC weekly solution.
SF = (Variance Factor)^0.5 = (Chi^2/dof)^0.5
Final(SF) = Apriori(SF) * Aposteriori(SF)
Solution |# sta| #RF | Apriori(SF)|Aposteriori(SF)| Final(SF)|
---------------|-----|-----|---------------|---------------|---------------|
cod | 173 | 67 | 21.124 | 1.000 | 21.124 |
emr | 41 | 28 | 5.385 | 1.000 | 5.385 |
esa | 66 | 48 | 16.001 | 1.000 | 16.001 |
gfz | 182 | 75 | 4.696 | 1.000 | 4.696 |
jpl | 110 | 43 | 6.109 | 1.000 | 6.109 |
MIT | 148 | 61 | 2.554 | 0.998 | 2.548 |
ngs | 186 | 69 | 36.816 | 1.000 | 36.817 |
sio | 162 | 68 | 2.197 | 0.997 | 2.191 |
mit | 293 | 75 | 0.913 | 1.000 | 0.913 |
ncl | 174 | 75 | 2.196 | 0.998 | 2.191 |
igs | 271 | 77 | 1.399 | 0.998 | 1.396 |
---------------|-----|-----|---------------|---------------|---------------|
5-2) Solution Residuals weighted Average & STD :
=====================================================
5-2-1) With respect to the IGb00 stations:
=============================================
Solution | Weighted Average (mm) | STD (mm) |
| N E U | N E U |
---------------|---------------------------|---------------------------|
cod 67 | 0.0 0.2 -0.5 | 4.1 3.0 11.5 |
emr 27 | -1.5 0.0 0.4 | 4.5 3.9 13.0 |
esa 47 | -1.5 0.3 0.2 | 4.4 3.2 11.2 |
gfz 75 | -1.9 0.0 -0.3 | 4.3 3.7 12.9 |
jpl 43 | 0.1 0.1 0.5 | 4.0 4.0 12.5 |
MIT 61 | 0.2 0.1 0.3 | 4.0 4.2 10.7 |
ngs 69 | -1.6 -0.4 -0.0 | 4.6 5.5 14.0 |
sio 68 | 0.0 0.3 0.4 | 4.2 4.4 11.2 |
mit 74 | 0.7 0.1 -0.2 | 3.7 3.7 11.1 |
ncl 75 | 0.1 0.2 -0.0 | 3.7 3.2 9.7 |
---------------|---------------------------|---------------------------|
5-2-2) With respect to the Weekly Combination:
=================================================
Solution | Weighted Average (mm) | STD (mm) |
| N E U | N E U |
---------------|---------------------------|---------------------------|
cod 169 | 0.6 0.1 -1.7 | 1.5 1.8 4.8 |
emr 39 | 0.2 0.8 -0.8 | 2.2 1.9 7.1 |
esa 63 | -0.5 -0.3 0.5 | 2.6 2.2 11.8 |
gfz 178 | -1.1 -0.1 -0.6 | 1.0 1.1 3.6 |
jpl 107 | 0.9 -0.3 -0.5 | 1.6 3.0 10.1 |
MIT 137 | 0.9 0.3 -2.1 | 1.6 2.3 4.6 |
ngs 178 | -0.1 -0.9 -1.5 | 1.8 3.4 6.7 |
sio 154 | 0.7 0.2 -2.3 | 1.8 2.6 4.3 |
mit 262 | 1.5 0.2 -1.7 | 0.8 0.9 2.2 |
ncl 168 | 0.9 -0.1 -1.4 | 0.9 0.8 2.1 |
ITR 76 | 0.8 0.0 -1.6 | 3.6 3.6 10.8 |
---------------|---------------------------|---------------------------|
5-2-3) With respect to the Cumulative Combination:
=====================================================
Solution | Weighted Average (mm) | STD (mm) |
| N E U | N E U |
---------------|---------------------------|---------------------------|
cod 157 | -0.7 0.3 -0.6 | 2.6 2.2 7.0 |
emr 39 | -1.4 1.1 0.1 | 2.9 2.3 8.4 |
esa 64 | -2.1 -0.4 0.1 | 3.2 2.8 12.1 |
gfz 167 | -2.3 -0.0 -0.0 | 2.3 2.2 6.6 |
jpl 98 | -0.5 -0.1 0.2 | 2.6 3.6 10.7 |
MIT 130 | -0.4 0.3 -1.4 | 2.3 2.8 5.9 |
ngs 152 | -1.4 -0.6 -0.8 | 2.9 4.7 9.2 |
sio 145 | -0.5 0.1 -2.0 | 2.7 3.2 6.7 |
mit 214 | 0.4 0.3 -0.8 | 2.1 2.2 5.4 |
ncl 161 | -0.3 -0.1 -0.5 | 2.2 2.1 4.9 |
igs 220 | -1.1 0.1 0.6 | 1.9 1.9 5.2 |
ITR 84 | -0.4 -0.1 -0.8 | 2.6 2.8 8.1 |
---------------|---------------------------|---------------------------|
5-3) 7-Parameters Transformation: AC ----> IGS03P33_RS99.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.017 -0.016 0.008 13.5 2.4 -1.1 3.30 |
emr | 0.013 -0.012 0.005 0.5 7.1 -13.7 4.21 |
esa | 0.005 0.019 0.023 0.0 6.4 -13.6 1.92 |
gfz | 0.037 -0.011 -0.011 -1.8 2.0 -4.5 2.54 |
jpl | 0.180 -0.017 0.135 -6.0 6.0 -5.4 1.56 |
MIT | -0.019 -0.018 -0.004 8.8 14.0 -6.4 3.86 |
ngs | 0.053 -0.110 0.015 11.1 4.5 1.3 3.23 |
sio | 0.099 -0.026 -0.015 1.2 4.8 -6.6 3.61 |
mit | -0.007 0.009 -0.017 -0.8 -1.1 -1.2 3.14 |
ncl | 0.150 -0.108 0.044 2.0 1.7 -5.4 3.07 |
igs | 0.000 0.000 0.000 0.0 0.0 0.0 3.04 |
---------------|---------------------------------------------------------------|
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.049 0.049 0.055 1.3 1.3 1.3 0.20 |
emr | 0.073 0.078 0.086 2.0 2.0 2.0 0.31 |
esa | 0.055 0.057 0.065 1.5 1.5 1.5 0.23 |
gfz | 0.046 0.045 0.053 1.2 1.2 1.2 0.19 |
jpl | 0.060 0.059 0.068 1.5 1.5 1.5 0.24 |
MIT | 0.050 0.050 0.057 1.3 1.3 1.3 0.20 |
ngs | 0.047 0.047 0.055 1.2 1.2 1.2 0.19 |
sio | 0.047 0.048 0.054 1.2 1.2 1.2 0.19 |
mit | 0.045 0.046 0.053 1.2 1.2 1.2 0.19 |
ncl | 0.046 0.046 0.052 1.2 1.2 1.2 0.19 |
igs | 0.000 0.000 0.000 0.0 0.0 0.0 0.94 |
---------------|---------------------------------------------------------------|
5-4) Geocenter:
====================
| X Y Z | s X s Y s Z |
| (mm) (mm) (mm) | (mm) (mm) (mm) |
---------------|---------------------------|---------------------------|
igs | 5.9 7.9 -6.2 | 3.0 3.1 4.1 |
---------------|---------------------------|---------------------------|
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 | --- 6.5 0.8 --- 0.6 --- |
emr | --- 11.3 -4.0 13.0 -1.4 -18.8 |
esa | --- 3.9 4.3 --- 3.0 --- |
gfz | --- -4.5 4.0 2.2 -3.7 3.8 |
jpl | --- 4.4 -3.0 7.6 0.8 -10.1 |
MIT | --- -14.9 3.2 -0.6 -3.7 -0.5 |
ngs | --- 17.3 0.7 40.2 11.9 -2.9 |
sio | --- -18.0 -0.8 0.0 -7.0 2.7 |
mit | --- 9.0 1.8 -0.4 -1.8 0.4 |
ncl | --- 13.6 3.0 --- 2.0 --- |
ITR | --- --- --- --- --- --- |
BuA | --- 11.5 -1.3 -0.5 34.8 0.8 |
---------------|------------------------------------------------------------|
5-5-2) STD :
===============
Solution | UT LOD XPO XPOR YPO YPOR |
| us us .01mas .01ma/d .01mas .01ma/d |
---------------|------------------------------------------------------------|
cod | --- 9.2 3.2 --- 1.3 --- |
emr | --- 6.2 2.2 26.8 6.6 10.7 |
esa | --- 20.9 3.4 --- 4.0 --- |
gfz | --- 10.9 2.5 6.4 1.4 14.4 |
jpl | --- 16.4 1.8 26.7 3.3 15.7 |
MIT | --- 3.3 3.7 7.4 2.5 4.7 |
ngs | --- 30.7 6.5 10.1 7.3 19.0 |
sio | --- 6.6 4.3 9.0 4.6 5.3 |
mit | --- 8.0 4.2 7.9 2.6 9.6 |
ncl | --- 3.7 1.3 --- 1.3 --- |
ITR | --- --- --- --- --- --- |
BuA | --- 14.2 9.2 14.0 8.6 13.6 |
---------------|------------------------------------------------------------|
5-6) Outliers:
===================
Thresholds set at 5 sigmas and 50 mm.
Solution Solution | Station | N E U |
Minus |code id soln| (mm) (mm) (mm)|
-----------------------------------------|------------|-----------------------|
| | |
IGS03P33_RS99.snx cod13967.snx |GOUG A ----| 7.5 -2.1 74.2|
IGS03P33_RS99.snx cod13967.snx |HOFN A 2| 7.2 -9.2 -37.9|
IGS03P33_RS99.snx cod13967.snx |URUM A ----| 2.8 -10.1 20.2|
IGS03P33_RS99.snx cod13967.snx |ZIMM A ----| -2.9 7.9 -29.6|
| | |
IGS03P33_RS99.snx emr13967.snx |GOUG A ----| 7.3 -2.1 53.6|
| | |
IGS03P33_RS99.snx jpl13967.snx |GOUG A ----| 15.4 -12.8 66.3|
IGS03P33_RS99.snx jpl13967.snx |OHI2 A ----| -16.0 -1.5 -56.9|
| | |
IGS03P33_RS99.snx MIT13967.snx |DARW A ----| -11.0 0.5 4.4|
IGS03P33_RS99.snx MIT13967.snx |OHI2 A ----| -13.9 0.0 -44.9|
| | |
IGS03P33_RS99.snx sio13967.snx |GOUG A ----| 3.9 -4.6 50.8|
IGS03P33_RS99.snx sio13967.snx |HOFN A 2| 7.5 -10.0 -40.8|
| | |
IGS03P33_RS99.snx ncl1396g.snx |HOFN A 2| 6.5 -11.7 -32.5|
IGS03P33_RS99.snx ncl1396g.snx |OHI2 A ----| -12.1 0.1 -45.4|
| | |
emr13967.snx esa13967.snx |VILL A ----| 0.3 -0.3 -50.2|
| | |
jpl13967.snx MIT13967.snx |SUNM A ----| 8.2 -7.2 63.2|
| | |
IGS06P40.snx MIT13967.snx |SAMP A 2| -20.2 -19.3 4.8|
| | |
IGS06P40.snx mit1396g.snx |OHI2 A ----| -2.7 -6.0 -35.4|
| | |
IGS06P40.snx igs06P1396.snx |ISPA A ----| -1.3 8.4 8.2|
IGS06P40.snx igs06P1396.snx |KOUR A ----| -9.6 -2.4 6.2|
IGS06P40.snx igs06P1396.snx |NTUS A 2| -3.4 -7.2 5.1|
-----------------------------------------|------------|-----------------------|
5-7) Conflicts:
====================
IERS DOMES NO MIT cod emr esa gfz jpl ngs sio IGS igs
============= ====== ====== ====== ====== ====== ====== ====== ====== ====== ======
43001M002 THU2 A THU3 A THU3 A THU3 A THU3 A THU3 A THU2 A THU3 A THU3 A
THU3 A
42202M005 AREQ A AREQ A AREQ A AREQ A AREQ A AREQ A AREQ B AREQ A AREQ A AREQ A
40451M123 GODE A GODZ A GODE A GODE A GODE A GODE A GODE A GODE A
40405S038 GOL2 A
40405S031 GOLD A GOLD A GOLD A GOLD A GOL2 A GOLD A GOLD A
50103M108 TIDB A TIDB A TIDB A TIDB A TIDB A TID1 A TID2 A TIDB A TIDB A
TID1 A TIDB A
TID2 A
50134M002 DARR A
50134M001 DARW A DARR A DARW A DARW A DARW A DARR A DARW A DARW A DARW A DARW A
DARW A
22003M001 PIMO A PIMO A PIMO A PIMO B PIMO A PIMO A PIMO A PIMO A PIMO A PIMO A
22306M002 IISC A IISC A IISC B IISC A IISC A IISC B IISC A IISC A IISC A
43401M001 BDOS A BDOS A BDOS A BDOS A BDOS A
BDOS A BDOS A
43005M002 KELY B KELY B KELY B KELY B KELY B KELY B KELY A KELY A
10422M001 KIR0 A KIR0 A KR0G A KR0G A
20703S001 RAMO A RAMO A RAMO A RAMO A ELRO A RAMO A
RAMO A
20808M001 TRAB A TRAB A TRAB A TRAB A TRAB A
MISSING IERS DOMES NOS AT MIT:
-----------------------------
HLFX A PALM A PUR3 A MIA3 A KYW1 A MCD1 A ELEN A AZCN A OUSD A MTJO A
HORN A CNCL A QUAR A KARA A SAMP A MAR6 A IENG A
MISSING IERS DOMES NOS AT gfz:
-----------------------------
PALM A
MISSING IERS DOMES NOS AT ngs:
-----------------------------
SVTL A HAMM A MCN1 A MLF1 A SUTB A
MISSING IERS DOMES NOS AT sio:
-----------------------------
KULU A SG21 A OUSD A CNCL A KARA A
MISSING IERS DOMES NOS AT igs:
-----------------------------
AZCN A FAAA A
IERS DOMES NO mit ncl IGS igs
============= ====== ====== ====== ======
13407S014 MAD2 A
13407S012 MADR A MAD2 A MADR A MADR A
MADR A
43005M002 KELY B KELY B KELY A KELY A
43401M001 BDOS A BDOS A BDOS A BDOS A
BDOS A BDOS A
43001M002 THU2 A THU3 A THU3 A THU3 A
40451M123 GODZ A GODE A GODE A GODE A
GODE A
50103M108 TIDB A TID1 A TIDB A TIDB A
TID2 A
TIDB A
50134M002 DARR A
50134M001 DARW A DARR A DARW A DARW A
DARW A
12309M005 MDVJ A MDVJ A MDVJ A MDVJ A
MISSING IERS DOMES NOS AT mit:
-----------------------------
SUTB A AIS1 A
MISSING IERS DOMES NOS AT igs:
-----------------------------
AZCN A FAAA A
More information about the IGSREPORT
mailing list