From rferland at nrcan.gc.ca Thu Jan 13 19:48:30 2005 From: rferland at nrcan.gc.ca (Remi Ferland) Date: Thu, 13 Jan 2005 22:48:30 -0500 (EST) Subject: [IGSREPORT-12143] Wk 1303 IGS SINEX Combination Message-ID: <200501140348.j0E3mUF11544@cassini.geod.nrcan.gc.ca> ****************************************************************************** IGS Electronic Report 13 Jan 19:49:54 PST 2005 Message Number 12143 ****************************************************************************** Author: R. Ferland , D. Hutchison =================================================================== Combined AC solution for GPS week 1303 - 04/12/26 to 05/01/01 =================================================================== 1) Contacts: ============ e-mail: ferland at geod.nrcan.gc.ca Tel: (613)995-4002 (Remi Ferland) Fax: (613)995-3215 2) Products: ============ All products are available at cddisa.gsfc.nasa.gov and at garner.ucsd.edu igs04P1303.snx Combined weekly solution. igs04P1303.erp Daily Earth Rotation Parameters. igs00p02.erp Daily Earth Rotation Parameters since GPS week 1013. IGS04P52.snx Cumulative Combined solution. igs04P1303.sum This file. igs04P1303.itr Residuals between the AC and the IGb00 for the RF stations. igs04P1303.res Residuals between the AC and the Combined Weekly solution. IGS04P52.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 = igs04P1303.snx IGS = IGS04P52.snx MIT = mit13037.snx mit = mit1303g.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 BAKO jpl GMAS New stations in the cumulative solution IGS04P52.snx MAC1 A 2 NTUS A 2 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 | 163 | 74 | 35.232 | 0.828 | 29.172 | emr | 49 | 35 | 7.646 | 0.788 | 6.025 | esa | 77 | 59 | 10.524 | 0.855 | 9.001 | gfz | 179 | 81 | 10.135 | 0.998 | 10.110 | jpl | 121 | 55 | 5.816 | 0.993 | 5.773 | MIT | 140 | 64 | 3.438 | 0.980 | 3.369 | ngs | 130 | 71 | 33.375 | 0.946 | 31.566 | sio | 170 | 79 | 3.059 | 0.929 | 2.843 | mit | 268 | 81 | 1.080 | 0.971 | 1.048 | ncl | 164 | 82 | 3.726 | 0.813 | 3.030 | igs | 237 | 82 | 1.345 | 0.973 | 1.309 | ---------------|-----|-----|---------------|---------------|---------------| 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 74 | 0.2 0.2 0.4 | 3.1 3.6 12.5 | emr 35 | -1.1 0.4 1.9 | 3.6 3.1 13.3 | esa 58 | -1.4 0.3 1.7 | 3.2 3.1 13.8 | gfz 81 | -0.7 -0.0 0.1 | 3.4 3.2 12.1 | jpl 55 | -0.9 0.3 0.8 | 3.5 3.4 11.9 | MIT 63 | 0.0 0.3 0.7 | 3.2 3.7 11.5 | ngs 71 | -2.7 0.1 0.5 | 4.6 4.6 13.3 | sio 79 | 0.3 0.5 -0.3 | 3.5 3.7 9.1 | mit 81 | -0.3 0.2 0.8 | 2.8 3.0 10.1 | ncl 82 | 0.3 0.1 0.4 | 2.8 2.8 9.3 | ---------------|---------------------------|---------------------------| 5-2-2) With respect to the Weekly Combination: ================================================= Solution | Weighted Average (mm) | STD (mm) | | N E U | N E U | ---------------|---------------------------|---------------------------| cod 157 | -0.3 0.1 -1.6 | 1.4 1.4 5.1 | emr 47 | -1.3 -0.1 -0.2 | 2.1 1.3 7.7 | esa 73 | -1.5 -0.2 -0.4 | 2.8 3.0 12.0 | gfz 171 | -1.3 -0.2 -0.9 | 1.5 1.6 4.5 | jpl 119 | -0.9 -0.1 -0.8 | 2.0 2.4 8.3 | MIT 121 | 0.0 -0.1 -1.5 | 1.2 1.7 4.7 | ngs 122 | -2.7 -0.2 -1.0 | 2.9 3.8 8.1 | sio 150 | 0.5 0.4 -0.4 | 2.4 2.3 6.3 | mit 229 | -0.3 0.0 -0.8 | 0.6 0.5 2.6 | ncl 155 | -0.0 -0.1 -0.4 | 0.9 0.5 2.8 | ITR 81 | -0.0 -0.1 -0.8 | 2.5 2.6 10.4 | ---------------|---------------------------|---------------------------| 5-2-3) With respect to the Cumulative Combination: ===================================================== Solution | Weighted Average (mm) | STD (mm) | | N E U | N E U | ---------------|---------------------------|---------------------------| cod 130 | 0.1 0.6 -0.2 | 2.9 2.9 9.1 | emr 45 | -1.8 -0.6 2.3 | 3.0 3.8 10.8 | esa 74 | -1.6 -0.5 0.3 | 4.3 4.8 14.6 | gfz 156 | -1.1 0.1 0.3 | 2.7 2.7 8.5 | jpl 100 | -1.1 -0.0 0.9 | 3.1 3.3 9.9 | MIT 109 | 0.0 -0.0 -0.4 | 2.3 3.5 7.7 | ngs 119 | -2.4 -0.4 -0.1 | 4.2 6.0 10.6 | sio 142 | 0.9 0.6 -0.0 | 3.8 4.3 7.8 | mit 173 | -0.2 0.4 0.4 | 2.3 2.5 7.3 | ncl 145 | 0.3 0.2 0.5 | 2.5 2.4 6.4 | igs 173 | -0.1 0.3 0.8 | 1.9 1.9 6.6 | ITR 89 | -0.4 0.0 -0.5 | 1.2 1.2 3.8 | ---------------|---------------------------|---------------------------| 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.033 0.031 -0.010 14.3 11.8 4.3 2.96 | emr | -0.002 -0.078 0.029 6.4 10.4 -9.3 4.09 | esa | -0.014 -0.031 0.013 -4.5 9.4 2.1 3.04 | gfz | 0.050 -0.071 0.020 0.3 1.1 -4.8 2.18 | jpl | 0.144 -0.071 0.123 1.7 14.4 -5.1 1.49 | MIT | 0.184 -0.061 -0.006 8.3 10.6 36.6 3.73 | ngs | -0.072 0.026 0.007 2.4 3.4 53.7 3.19 | sio | -0.075 -0.026 0.019 19.0 8.9 49.5 3.20 | mit | -0.010 0.014 -0.034 -1.1 -0.5 -1.3 3.11 | ncl | 0.019 -0.054 -0.014 1.4 1.4 -3.5 2.78 | ---------------|---------------------------------------------------------------| 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.046 0.048 0.052 1.2 1.2 1.2 0.19 | emr | 0.064 0.066 0.078 1.7 1.7 1.7 0.27 | esa | 0.051 0.053 0.057 1.3 1.3 1.3 0.21 | gfz | 0.044 0.044 0.050 1.2 1.1 1.1 0.18 | jpl | 0.052 0.053 0.058 1.4 1.4 1.4 0.21 | MIT | 0.049 0.051 0.055 1.3 1.3 1.3 0.20 | ngs | 0.046 0.047 0.053 1.2 1.2 1.2 0.19 | sio | 0.045 0.046 0.051 1.2 1.2 1.2 0.18 | mit | 0.044 0.045 0.051 1.2 1.2 1.2 0.18 | ncl | 0.043 0.044 0.050 1.2 1.1 1.1 0.18 | ---------------|---------------------------------------------------------------| 5-4) Geocenter: ==================== | X Y Z | s X s Y s Z | | (mm) (mm) (mm) | (mm) (mm) (mm) | ---------------|---------------------------|---------------------------| igs | 8.8 11.4 7.4 | 3.9 4.0 5.4 | ---------------|---------------------------|---------------------------| 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.6 -6.0 --- 0.0 --- | emr | --- 24.2 -0.5 11.1 -4.1 12.6 | esa | --- 34.5 1.3 --- -0.0 --- | gfz | --- -9.5 6.1 0.1 5.9 16.7 | jpl | --- 28.9 -2.7 17.9 -0.6 -1.5 | MIT | --- -8.3 -0.9 1.0 1.6 -2.7 | ngs | --- 3.3 -4.7 -12.1 10.2 0.8 | sio | --- 3.3 -0.4 0.9 5.9 -0.2 | mit | --- 18.3 -2.7 0.4 1.8 -1.9 | ncl | --- 2.8 5.4 --- 2.7 --- | ITR | --- --- --- --- --- --- | BuA | --- 3.5 -5.8 -4.7 32.7 0.3 | ---------------|------------------------------------------------------------| 5-5-2) STD : =============== Solution | UT LOD XPO XPOR YPO YPOR | | us us .01mas .01ma/d .01mas .01ma/d | ---------------|------------------------------------------------------------| cod | --- 6.3 0.9 --- 1.8 --- | emr | --- 10.5 4.0 8.3 2.3 10.5 | esa | --- 15.1 2.8 --- 5.1 --- | gfz | --- 6.7 1.7 9.4 2.7 6.9 | jpl | --- 18.6 1.5 44.4 3.9 8.0 | MIT | --- 4.0 2.7 3.7 3.6 8.4 | ngs | --- 16.7 6.0 46.5 11.0 50.8 | sio | --- 3.9 1.6 10.4 2.6 14.0 | mit | --- 7.4 3.4 5.9 4.6 10.4 | ncl | --- 1.4 1.0 --- 2.0 --- | ITR | --- --- --- --- --- --- | BuA | --- 11.8 5.0 13.6 8.7 17.4 | ---------------|------------------------------------------------------------| 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 cod13037.snx |GOUG A ----| -0.2 4.4 61.6| | | | esa13037.snx MIT13037.snx |MCM4 A ----| 4.1 -4.0 51.8| | | | IGS04P51.snx cod13037.snx |IISC A ----| -1.9 -19.3 -1.5| IGS04P51.snx cod13037.snx |YAKT A ----| 23.6 -25.0 -23.7| | | | IGS04P51.snx gfz13037.snx |IISC A ----| -3.0 -19.2 5.6| IGS04P51.snx gfz13037.snx |YAKT A ----| 25.1 -27.2 -29.6| | | | IGS04P51.snx gfz13037.snx |KIRU A ----| -9.3 3.4 -52.9| IGS04P51.snx gfz13037.snx |RESO A ----| 14.1 2.9 -7.6| | | | IGS04P51.snx jpl13037.snx |IISC A ----| -1.9 -19.5 -6.3| IGS04P51.snx jpl13037.snx |YAKT A ----| 16.4 -21.5 -17.8| | | | IGS04P51.snx MIT13037.snx |TVST A ----| -18.9 4.6 -38.3| | | | IGS04P51.snx sio13037.snx |ELAT A ----| -3.8 -12.7 50.2| IGS04P51.snx sio13037.snx |KIRU A ----| -12.6 -1.6 -53.3| | | | IGS04P51.snx mit1303g.snx |IISC A ----| -0.6 -18.9 -0.3| IGS04P51.snx mit1303g.snx |KIRU A ----| -10.0 -0.5 -52.6| IGS04P51.snx mit1303g.snx |USUD A ----| 12.6 3.4 4.8| IGS04P51.snx mit1303g.snx |YAKT A ----| 22.8 -25.4 -25.5| | | | IGS04P51.snx ncl1303g.snx |IISC A ----| -0.3 -18.4 -0.6| IGS04P51.snx ncl1303g.snx |KIRU A ----| -9.5 -0.3 -54.7| IGS04P51.snx ncl1303g.snx |YAKT A ----| 26.1 -27.3 -20.4| | | | IGS04P51.snx igs04P1303.snx |ALGO A ----| -0.8 -3.1 20.8| IGS04P51.snx igs04P1303.snx |BAN2 A ----| -0.8 -13.1 -1.0| IGS04P51.snx igs04P1303.snx |CHUR A ----| -1.8 -6.2 6.1| IGS04P51.snx igs04P1303.snx |IISC A ----| -0.6 -18.3 -4.9| IGS04P51.snx igs04P1303.snx |MAG0 A ----| 8.2 9.7 -7.7| IGS04P51.snx igs04P1303.snx |RESO A ----| 14.1 5.1 -6.0| IGS04P51.snx igs04P1303.snx |SCH2 A ----| 6.7 9.0 4.5| IGS04P51.snx igs04P1303.snx |USUD A ----| 12.1 2.4 2.9| IGS04P51.snx igs04P1303.snx |YAKT A ----| 28.3 -33.3 -16.6| -----------------------------------------|------------|-----------------------| 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 THU3 A THU3 A THU3 A THU3 A THU2 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 TID1 A TIDB A TID2 A TIDB A TIDB A TID1 A TIDB A TID2 A 50134M001 DARW A DARR 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 66010M001 DAV1 A DAVR A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A 22306M002 IISC A IISC A IISC A IISC B IISC A IISC A IISC B IISC A IISC A 10422M001 KIR0 A KIR0 A KR0G A 24201M001 ULAB A ULAB A ULAB A ULAB A ULAB A ULAB A MISSING IERS DOMES NOS AT MIT: ----------------------------- HLFX A PALM A PUR3 A MIA3 A KYW1 A MCD1 A ARP3 A GUAX A UPO1 A OUSD A MTJO A HORN A CNCL A QUAR A KARA A MAC1 B TVST A KAYT A MAR6 A IENG A MISSING IERS DOMES NOS AT gfz: ----------------------------- BISH A PALM A MISSING IERS DOMES NOS AT ngs: ----------------------------- CHUM A MISSING IERS DOMES NOS AT sio: ----------------------------- KULU A SG21 A DOAR A OUSD A CNCL A KARA A MERS A CAVA A SFEL A VOLT A BRAS A IERS DOMES NO mit ncl IGS igs ============= ====== ====== ====== ====== 13407S014 MAD2 A 13407S012 MADR A MAD2 A MADR A MADR A MADR 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 11101M002 SOFI A SOFI A SOFI A MISSING IERS DOMES NOS AT mit: ----------------------------- MAC1 B BISH A