From mbourass at nrcan.gc.ca Wed Apr 18 21:36:40 2007 From: mbourass at nrcan.gc.ca (Martin Bourassa) Date: Thu, 19 Apr 2007 00:36:40 -0400 Subject: [IGSREPORT-15005] Wk 1421 IGS SINEX Combination Message-ID: <200704190436.l3J4aeF1008715@newton.geod.nrcan.gc.ca> ****************************************************************************** IGS Electronic Report 18 Apr 21:36:55 PDT 2007 Message Number 15005 ****************************************************************************** Author: R. Ferland, M. Bourassa =============================================================== Combined AC solution for GPS week 1421 - 07/04/01 to 07/04/07 =============================================================== 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) igs07P1421.snx Combined weekly solution. igs07P1421.erp Daily Earth Rotation Parameters. igs00p03.erp Daily Earth Rotation Parameters since GPS week 860. IGS07P14.snx Cumulative Combined solution. igs07P1421.sum This file. igs07P1421.itr Residuals between the AC and the IGS05 for the RF stations. igs07P1421.res Residuals between the AC and the Combined Weekly solution. IGS07P14.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 IGS05 stations at the current epoch. ftp://igscb.jpl.nasa.gov/igscb/station/coord/IGS05.snx 3-4) The AC are compared to IGS05, 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 IGS05 of the current epoch using a 7-parameters transformation. The singular cumulative combined solution is aligned to the IGS05 of the reference epoch using a 14-parameters transformation and the 132 ITRF station positions and velocities. The IGS05 is a realization of ITRF2005. 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 = IGS05.snx igs = igs07P1421.snx IGS = IGS07P14.snx MIT = mit14217.snx mit = mit1421g.snx BuA = IERS Bulletin A jpl 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: From the original solution: cod SATA X Y Z emr SATA X Y Z gfz SATA X Y Z MIT SATA X Y Z sio SATA X Y Z mit SATA X Y Z After unconstraining: jpl AUCK GMSD ncl GMSD 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: ================================================== 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 | 162 | 95 | 23.567 | 1.000 | 23.558 | emr | 56 | 43 | 6.089 | 1.000 | 6.086 | esa | 68 | 57 | 18.989 | 1.000 | 18.983 | gfz | 180 | 102 | 5.284 | 1.000 | 5.282 | MIT | 148 | 90 | 2.685 | 0.999 | 2.683 | ngs | 162 | 107 | 36.258 | 1.000 | 36.250 | sio | 158 | 97 | 2.515 | 1.000 | 2.515 | jpl | 121 | 70 | 7.692 | 1.000 | 7.690 | mit | 247 | 110 | 1.882 | 1.000 | 1.881 | ncl | 164 | 106 | 4.139 | 1.000 | 4.138 | igs | 249 | 114 | 1.491 | 1.000 | 1.490 | ---------------|-----|-----|---------------|---------------|---------------| 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 94 | -1.0 0.0 0.3 | 3.7 2.8 8.1 | emr 43 | -0.3 0.3 0.4 | 3.5 2.5 9.7 | esa 57 | -1.3 0.2 -0.3 | 4.3 2.7 12.0 | gfz 102 | -0.4 0.2 0.1 | 2.8 2.3 7.1 | MIT 90 | -0.4 0.2 0.4 | 2.6 2.0 7.9 | ngs 107 | 0.3 0.0 0.0 | 2.5 2.3 5.8 | sio 96 | 0.6 0.1 0.5 | 2.5 2.4 7.6 | jpl 70 | -0.3 0.3 0.6 | 3.2 3.0 14.2 | mit 110 | -0.1 0.1 0.2 | 2.7 2.0 5.9 | ncl 106 | 0.0 0.0 0.3 | 0.4 0.4 1.2 | ---------------|---------------------------|---------------------------| 5-2-2) With respect to the Weekly Combination: ================================================= Solution | Weighted Average (mm) | STD (mm) | | N E U | N E U | ---------------|---------------------------|---------------------------| cod 159 | -0.6 -0.6 -0.5 | 1.9 2.2 4.1 | emr 56 | 0.0 0.5 0.2 | 1.7 1.6 7.3 | esa 68 | -0.9 0.6 1.1 | 3.4 2.5 12.2 | gfz 178 | -0.1 0.2 -0.3 | 1.0 1.5 3.6 | MIT 139 | -0.1 0.2 0.0 | 0.9 1.2 3.3 | ngs 159 | 0.6 -0.1 -0.5 | 1.7 1.9 4.2 | sio 153 | 0.9 0.1 0.0 | 0.8 1.7 2.9 | jpl 114 | 0.4 0.3 1.2 | 2.9 2.7 13.8 | mit 225 | 0.1 -0.1 -0.2 | 1.1 1.0 2.5 | ncl 159 | 0.3 0.0 -0.1 | 2.3 1.7 5.9 | ITR 113 | 0.2 -0.1 -0.4 | 2.5 1.7 6.0 | ---------------|---------------------------|---------------------------| 5-2-3) With respect to the Cumulative Combination: ===================================================== Solution | Weighted Average (mm) | STD (mm) | | N E U | N E U | ---------------|---------------------------|---------------------------| cod 152 | -1.0 -0.1 -0.2 | 2.6 2.9 6.3 | emr 55 | -0.4 0.7 -0.2 | 3.0 2.9 8.4 | esa 68 | -1.0 0.5 1.2 | 3.4 3.4 12.9 | gfz 169 | -0.4 0.4 -0.1 | 2.0 2.6 6.3 | MIT 134 | -0.4 0.4 0.6 | 1.9 2.3 5.9 | ngs 154 | 0.3 0.3 -0.3 | 2.1 2.5 5.2 | sio 145 | 0.5 0.0 0.7 | 1.5 2.7 6.6 | jpl 109 | 0.2 0.6 2.1 | 3.0 3.4 14.3 | mit 210 | -0.2 0.2 0.2 | 1.9 2.1 4.8 | ncl 154 | 0.0 0.2 0.6 | 1.9 1.8 5.1 | igs 229 | -0.4 0.3 0.4 | 1.6 1.9 4.9 | ITR 129 | -0.1 0.0 -0.1 | 1.8 1.3 5.2 | ---------------|---------------------------|---------------------------| 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.004 0.006 -0.001 0.0 0.1 0.4 -0.05 | emr | -0.013 0.016 0.012 8.5 -16.7 5.9 0.27 | esa | -0.015 -0.002 -0.003 4.7 -2.0 -3.4 -0.87 | gfz | -0.006 -0.030 0.000 0.5 -0.2 0.5 -0.11 | MIT | -0.018 0.004 -0.008 3.7 -2.8 15.4 0.34 | ngs | -0.030 -0.052 0.027 5.5 -1.5 6.0 -0.91 | sio | -0.079 -0.036 0.025 -0.3 -8.0 -9.6 0.44 | jpl | 0.163 -0.093 0.165 6.5 -5.7 22.7 0.51 | mit | -0.009 0.016 -0.007 -0.4 -0.4 0.6 -0.48 | ncl | -0.004 -0.000 0.001 0.4 -0.3 0.3 -0.18 | igs | 0.000 0.000 0.000 0.0 0.0 0.0 -0.10 | ---------------|---------------------------------------------------------------| 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.041 0.044 0.046 1.1 1.1 1.1 0.17 | emr | 0.060 0.062 0.069 1.6 1.6 1.6 0.25 | esa | 0.050 0.054 0.058 1.4 1.4 1.4 0.21 | gfz | 0.040 0.041 0.045 1.1 1.1 1.0 0.16 | MIT | 0.042 0.044 0.047 1.1 1.1 1.1 0.17 | ngs | 0.039 0.040 0.044 1.0 1.0 1.0 0.16 | sio | 0.040 0.043 0.046 1.1 1.1 1.1 0.17 | jpl | 0.047 0.050 0.055 1.3 1.3 1.3 0.20 | mit | 0.038 0.040 0.043 1.0 1.0 1.0 0.16 | ncl | 0.039 0.041 0.044 1.1 1.0 1.0 0.16 | igs | 0.000 0.000 0.000 0.0 0.0 0.0 0.56 | ---------------|---------------------------------------------------------------| 5-4) Geocenter: ==================== | X Y Z | s X s Y s Z | | (mm) (mm) (mm) | (mm) (mm) (mm) | ---------------|---------------------------|---------------------------| igs | 4.3 -1.7 5.3 | 3.1 3.2 4.6 | ---------------|---------------------------|---------------------------| 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 | --- --- -0.4 --- -7.1 --- | emr | --- 13.9 -0.4 22.5 -3.7 -13.9 | esa | --- -19.2 -1.6 --- -3.2 --- | gfz | --- -6.9 -3.4 -15.9 3.9 -20.8 | MIT | --- -10.0 0.9 -0.1 -0.5 -1.2 | ngs | --- 18.2 7.1 -7.0 1.3 5.5 | sio | --- 33.6 -6.4 -0.4 0.5 0.3 | jpl | --- 28.3 0.5 -0.7 -2.7 -19.2 | mit | --- 21.5 -7.4 0.0 0.3 0.5 | ncl | --- 35.9 2.6 --- -2.4 --- | ITR | --- --- --- --- --- --- | BuA | --- 17.6 -5.9 -0.5 29.6 0.0 | ---------------|------------------------------------------------------------| 5-5-2) STD : =============== Solution | UT LOD XPO XPOR YPO YPOR | | us us .01mas .01ma/d .01mas .01ma/d | ---------------|------------------------------------------------------------| cod | --- --- 1.9 --- 1.1 --- | emr | --- 5.4 0.7 13.1 2.4 11.8 | esa | --- 33.2 2.6 --- 2.6 --- | gfz | --- 12.9 1.9 4.9 1.8 10.9 | MIT | --- 5.3 1.1 1.9 1.9 4.9 | ngs | --- 8.5 4.7 16.1 9.5 17.3 | sio | --- 5.1 1.5 6.0 2.2 4.7 | jpl | --- 19.9 1.9 18.5 2.1 14.5 | mit | --- 7.6 1.5 9.4 2.4 5.2 | ncl | --- 4.1 1.6 --- 1.3 --- | ITR | --- --- --- --- --- --- | BuA | --- 21.4 2.8 6.8 3.2 12.4 | ---------------|------------------------------------------------------------| 5-6) Satellite Antenna Offset Residuals weighted Average & STD : ===================================================================== No Satellite phase center ... 5-7) Outliers: =================== Thresholds set at 5 sigmas and 50 mm. Solution Solution | Station | N E U | Minus |code id soln| (mm) (mm) (mm)| -----------------------------------------|------------|-----------------------| | | | ITRF_IGS05.snx MIT14217.snx |MAC1 A 4| 12.5 1.5 -5.5| | | | ITRF_IGS05.snx mit1421g.snx |TSKB A 2| -2.2 -4.2 -13.1| | | | cod14217.snx sio14217.snx |WUHN A 4| 15.0 1.4 -53.5| | | | cod14217.snx ncl1421g.snx |YSSK A 3| -0.4 21.3 2.5| | | | ngs14217.snx ncl1421g.snx |YSSK A 3| -5.2 17.5 -5.0| | | | igs07P1420.snx ncl1421g.snx |OUS2 A 2| 0.6 3.2 51.2| | | | IGS07P13.snx MIT14217.snx |ENG1 A 1| 26.3 -13.6 3.3| | | | IGS07P13.snx igs07P1421.snx |RESO A 1| 1.3 5.9 0.3| -----------------------------------------|------------|-----------------------| 5-8) Conflicts: ==================== IERS DOMES NO cod emr esa gfz jpl MIT ngs sio igs IGS ============= ====== ====== ====== ====== ====== ====== ====== ====== ====== ====== 43401M001 BDOS A BDOS A BDOS A BDOS A BDOS A BDOS A BDOS A 66010M001 DAVR A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A 40451M123 GODE A GODE A GODE A GODE A GODE A GODE A GODE A GODE A GODZ A 40405S031 GOLD A GOLD A GOLD A GOL2 A GOL2 A GOLD A GOL2 A GOLD A GOLD A 22306M002 IISC A IISC A IISC B IISC A IISC A IISC A IISC B IISC A IISC A IISC A 43005M002 KELY B KELY B KELY B KELY B KELY B KELY A KELY A KELY A 10422M001 KIR0 A KIR0 A KIR0 A KR0G A KR0G A 22003M001 PIMO A PIMO A PIMO B PIMO A PIMO A PIMO A PIMO A PIMO A PIMO A PIMO A 30314M002 SUTH A SUTH B SUTH A SUTH A SUTH A SUTH A SUTH A SUTH A 43001M002 THU2 A THU3 A THU3 A THU3 A THU3 A THU2 A THU3 A THU3 A 50107M004 YAR2 A YAR2 A YAR2 A YAR2 A YAR1 A YAR1 A 50103M108 TIDB A TIDB A TIDB A TIDB A TIDB A TID2 A TIDB A TIDB A TID1 A TIDB A TID2 A 13407S012 MADR A MAD2 A MADR A MAD2 A MADR A MADR A MADR A MADR A 40420M101 HARV A HARV A HARV A HARV A MISSING IERS DOMES NOS AT MIT: ----------------------------- IFRN A HLFX A PALM A MIA3 A KYW1 A ELEN A ENG1 A NETP A ADKS A AZCN A UPO1 A MTJO A CNCL A QUAR A KARA A MAR6 A MISSING IERS DOMES NOS AT sio: ----------------------------- KULU A SG21 A CNCL A KARA A IERS DOMES NO mit ncl igs IGS ============= ====== ====== ====== ====== 13407S012 MAD2 A MAD2 A MADR A MADR A MADR A MADR A 43005M002 KELY A 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 GODE A GODE A GODE A GODE A GODZ A 40405S031 GOLD A GOL2 A GOLD A GOLD A GOLD A 50103M108 TIDB A TID2 A TIDB A TIDB A TIDB A 50107M004 YAR2 A YAR2 A YAR1 A YAR1 A 20703S001 RAMO A RAMO A RAMO A ELRO A RAMO A 10422M001 KIR0 A KIR0 A KR0G A KR0G A 31101M001 WIND A WIND A WIND A WIND A MISSING IERS DOMES NOS AT mit: ----------------------------- IFRN A KULU A PALM A MIA3 A KYW1 A ENG1 A NETP A ADKS A SG21 A AZCN A UPO1 A MTJO A CNCL A QUAR A KARA A MAR6 A