From ferland at geod.nrcan.gc.ca Thu Jul 24 18:42:33 2003 From: ferland at geod.nrcan.gc.ca (Remi Ferland) Date: Thu, 24 Jul 2003 21:42:33 EDT Subject: [IGSREPORT-10456] Wk 1226 IGS SINEX Combination Message-ID: <200307250142.VAA27750@newton.geod.nrcan.gc.ca> ****************************************************************************** IGS Electronic Report 24 Jul 18:42:36 PDT 2003 Message Number 10456 ****************************************************************************** Author: R. Ferland , D. Hutchison =================================================================== Combined AC solution for GPS week 1226 - 03/07/06 to 03/07/12 =================================================================== 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 igs03P1226.snx Combined weekly solution. igs03P1226.erp Daily Earth Rotation Parameters. igs00p02.erp Daily Earth Rotation Parameters since GPS week 1013. IGS03P28.snx Cumulative Combined solution. igs03P1226.sum This file. igs03P1226.itr Residuals between the AC and the IGS00 for the RF stations. igs03P1226.res Residuals between the AC and the Combined Weekly solution. IGS03P28.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 IGS00 stations at the current epoch. ftp://igscb.jpl.nasa.gov/igscb/station/coord/IGS01P37_RS54.SNX.Z 3-4) The AC are compared to IGS00, 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 IGS00 of the current epoch using a 7-parameters transformation. The singular cumulative combined solution is aligned to the IGS00 of the reference epoch using a 14-parameters transformation and the 54 ITRF station positions and velocities. The IGS00 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 = IGS01P37_RS54.SNX igs = igs03P1226.snx IGS = IGS03P28.snx BuA = IERS Bulletin A mit ncl are used to quality control the weekly combination. gfz excluded from the geocenter estimation (Sec. 5-4). Solutions with manual station deletion: After unconstraining: esa SANT ngs KOKB LPGS 4-1) Starting with GPS week 1143, the solutions are aligned to the IGS01P37_RS54.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 | 169 | 40 | 24.770 | 1.000 | 24.766 | emr | 43 | 24 | 5.040 | 1.000 | 5.040 | esa | 58 | 38 | 11.647 | 1.000 | 11.648 | gfz | 177 | 44 | 7.056 | 1.000 | 7.055 | jpl | 110 | 32 | 4.885 | 1.000 | 4.885 | ngs | 76 | 45 | 27.874 | 1.000 | 27.873 | sio | 145 | 45 | 1.755 | 1.000 | 1.755 | mit | 240 | 44 | 1.731 | 1.067 | 1.848 | ncl | 136 | 45 | 1.648 | 1.000 | 1.648 | igs | 218 | 45 | 1.366 | 1.001 | 1.367 | ---------------|-----|-----|---------------|---------------|---------------| 5-2) Solution Residuals weighted Average & STD : ===================================================== 5-2-1) With respect to the IGS00 stations: ============================================= Solution | Weighted Average (mm) | STD (mm) | | N E U | N E U | ---------------|---------------------------|---------------------------| cod 39 | 1.3 0.4 0.4 | 3.1 3.9 13.2 | emr 23 | 0.7 -0.1 0.1 | 3.2 3.1 8.8 | esa 36 | 1.4 0.2 0.3 | 5.7 5.5 10.4 | gfz 43 | 1.7 0.5 -0.4 | 3.0 3.5 9.4 | jpl 31 | 1.5 0.1 0.4 | 3.5 2.6 10.7 | ngs 44 | -0.5 1.2 -0.5 | 5.0 4.9 14.0 | sio 44 | 1.8 0.5 0.3 | 3.5 3.3 6.7 | mit 43 | 1.5 0.5 0.4 | 3.2 3.4 7.0 | ncl 44 | 1.0 0.4 0.4 | 3.1 3.3 7.1 | ---------------|---------------------------|---------------------------| 5-2-2) With respect to the Weekly Combination: ================================================= Solution | Weighted Average (mm) | STD (mm) | | N E U | N E U | ---------------|---------------------------|---------------------------| cod 165 | 0.1 0.5 -1.5 | 1.3 1.4 7.3 | emr 42 | -0.3 -0.5 -3.3 | 2.0 1.8 6.4 | esa 57 | -1.2 0.1 0.1 | 5.8 4.6 7.1 | gfz 173 | 0.0 0.2 -0.4 | 1.2 1.4 4.9 | jpl 109 | 0.6 0.1 -0.2 | 1.9 1.3 6.4 | ngs 74 | -2.5 0.8 1.1 | 3.6 4.3 15.0 | sio 140 | 0.6 0.5 -0.0 | 1.8 2.7 4.7 | mit 211 | 0.1 0.1 0.1 | 1.5 1.1 2.9 | ncl 134 | -0.6 0.2 -0.1 | 1.0 0.6 3.0 | ITR 44 | -1.3 -0.8 0.6 | 2.9 2.8 7.3 | ---------------|---------------------------|---------------------------| 5-2-3) With respect to the Cumulative Combination: ===================================================== Solution | Weighted Average (mm) | STD (mm) | | N E U | N E U | ---------------|---------------------------|---------------------------| cod 133 | 1.1 0.4 -2.7 | 2.1 2.1 11.3 | emr 41 | 1.0 -0.5 -2.9 | 2.4 2.8 9.8 | esa 56 | -0.0 0.1 0.4 | 6.1 4.7 9.0 | gfz 144 | 1.0 0.1 -1.0 | 2.2 2.2 7.9 | jpl 87 | 1.6 -0.2 -1.0 | 2.6 2.2 10.1 | ngs 71 | -1.3 0.5 1.1 | 4.3 4.0 13.9 | sio 127 | 1.6 0.3 -1.0 | 3.1 3.4 6.0 | mit 155 | 1.1 0.1 -0.7 | 2.6 2.0 6.0 | ncl 125 | 0.5 0.1 -0.9 | 2.1 1.8 6.9 | igs 160 | 0.9 -0.0 -0.9 | 1.9 1.8 6.3 | ITR 54 | -0.6 -0.7 0.2 | 3.2 2.4 4.3 | ---------------|---------------------------|---------------------------| 5-3) 7-Parameters Transformation: AC ----> IGS01P37_RS54.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.008 0.005 -0.014 4.1 11.3 8.9 1.86 | emr | 0.025 -0.001 -0.001 -0.6 3.0 3.2 3.16 | esa | 0.158 -0.050 -0.035 -8.7 17.9 13.9 0.22 | gfz | 0.018 -0.048 0.050 -1.7 1.2 -0.6 0.86 | jpl | 0.007 -0.001 0.019 -2.2 3.5 6.2 0.55 | ngs | -0.046 -0.000 -0.052 -9.9 0.1 0.5 1.06 | sio | -0.020 -0.046 -0.057 5.1 2.2 24.3 1.40 | mit | 0.012 -0.001 -0.017 -1.0 -0.0 -2.4 1.33 | ncl | -0.003 -0.026 -0.037 -0.6 0.3 -2.1 1.25 | ---------------|---------------------------------------------------------------| 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.062 0.065 0.070 1.6 1.6 1.6 0.25 | emr | 0.084 0.081 0.093 2.1 2.1 2.1 0.33 | esa | 0.066 0.070 0.073 1.8 1.8 1.7 0.27 | gfz | 0.059 0.061 0.068 1.6 1.6 1.6 0.24 | jpl | 0.072 0.071 0.080 1.9 1.8 1.8 0.29 | ngs | 0.059 0.061 0.068 1.5 1.5 1.5 0.24 | sio | 0.058 0.060 0.068 1.5 1.5 1.5 0.24 | mit | 0.059 0.061 0.069 1.6 1.6 1.6 0.24 | ncl | 0.058 0.060 0.068 1.5 1.5 1.5 0.24 | ---------------|---------------------------------------------------------------| 5-4) Geocenter: ==================== | X Y Z | s X s Y s Z | | (mm) (mm) (mm) | (mm) (mm) (mm) | ---------------|---------------------------|---------------------------| igs | -1.4 4.2 6.4 | 2.1 2.1 3.9 | ---------------|---------------------------|---------------------------| 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 | --- 4.6 -3.9 --- -4.3 --- | emr | --- -4.1 3.7 2.0 -0.8 10.7 | esa | --- -5.8 -24.4 --- -22.7 --- | gfz | --- 7.2 2.2 0.5 -0.1 -2.7 | jpl | --- 10.6 2.4 -8.5 -0.1 5.5 | ngs | --- -13.1 15.9 -21.8 18.7 -18.4 | sio | --- -19.0 -3.8 1.9 2.1 -8.6 | mit | --- 15.9 -7.3 1.4 8.7 -9.7 | ncl | --- -3.7 -4.8 --- 6.2 --- | ITR | --- --- --- --- --- --- | BuA | --- -0.1 -6.4 1.7 26.8 -8.2 | ---------------|------------------------------------------------------------| 5-5-2) STD : =============== Solution | UT LOD XPO XPOR YPO YPOR | | us us .01mas .01ma/d .01mas .01ma/d | ---------------|------------------------------------------------------------| cod | --- 10.5 1.9 --- 3.1 --- | emr | --- 25.1 4.1 5.9 3.0 19.6 | esa | --- 7.3 5.5 --- 5.4 --- | gfz | --- 28.3 2.8 5.3 3.0 4.7 | jpl | --- 18.9 3.5 5.5 2.9 10.7 | ngs | --- 42.0 11.7 28.7 10.4 22.5 | sio | --- 11.4 1.7 7.1 1.5 5.1 | mit | --- 12.0 3.3 6.4 4.6 9.1 | ncl | --- 13.2 1.7 --- 1.3 --- | ITR | --- --- --- --- --- --- | BuA | --- 13.3 5.8 10.8 5.0 10.1 | ---------------|------------------------------------------------------------| 5-6) Outliers: =================== Thresholds set at 5 sigmas and 50 mm. Solution Solution | Station | N E U | Minus |code id soln| (mm) (mm) (mm)| -----------------------------------------|------------|-----------------------| | | | cod12267.snx jpl12267.snx |SUNM A ----| -1.5 1.7 -54.0| | | | IGS03P27.snx esa12267.snx |KOUR A ----| 17.3 -1.7 7.3| | | | IGS03P27.snx gfz12267.snx |OBET A ----| 0.8 3.9 168.1| | | | IGS03P27.snx gfz12267.snx |MIZU A ----| -6.1 13.1 -3.0| | | | IGS03P27.snx jpl12267.snx |GUAO A ----| -0.8 -6.4 -107.8| | | | IGS03P27.snx mit1226g.snx |GUAO A ----| -0.7 -6.6 -110.5| IGS03P27.snx mit1226g.snx |OBET A ----| 0.5 3.3 166.2| IGS03P27.snx mit1226g.snx |SUNM A ----| 2.0 -59.7 -159.0| | | | IGS03P27.snx mit1226g.snx |CASC A ----| 10.4 -6.3 -40.6| | | | IGS03P27.snx ncl1226g.snx |KELY A ----| -4.9 2.3 64.3| | | | igs03P1226.snx ncl1226g.snx |QAQ1 A ----| -2.9 -0.5 -31.4| | | | IGS03P27.snx igs03P1226.snx |LHAS A ----| 6.4 -4.4 -8.7| -----------------------------------------|------------|-----------------------| 5-7) Conflicts: ==================== IERS DOMES NO cod emr esa gfz jpl ngs sio IGS igs ============= ====== ====== ====== ====== ====== ====== ====== ====== ====== 50134M001 DARR A DARW A DARW A DARW A DARW A DARW A DARR A DARW A DARW A 66010M001 DAVR A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A 40405S031 GOLD A GOLD A GOLD A GOL2 A GOL2 A GOLD A GOLD A 22306M002 IISC A IISC A IISC B IISC A IISC A IISC B IISC A IISC A IISC A 13407S012 MADR A MADR A MAD2 A MAD2 A MADR A MADR A 10402M004 ONSA A ONSA A ONSA A ONSA A ONSA A ONSA A ONSA A ONSA A ONSA A OS0G A OS0G A 30314M002 SUTH A SUTH B 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 THU2 A THU3 A 50103M108 TIDB A TIDB A TIDB A TIDB A TID1 A TID2 A TIDB A TIDB A TIDB A TID1 A TID2 A 40508M002 CIC1 A CIC1 A CIC1 A CIC1 A MISSING IERS DOMES NOS AT gfz: ----------------------------- BISH A PALM A MISSING IERS DOMES NOS AT sio: ----------------------------- CAVA A SFEL A VOLT A BRAS A IERS DOMES NO mit ncl IGS igs ============= ====== ====== ====== ====== 13407S012 MADR A MADR A MADR A MAD2 A 43001M002 THU2 A THU3 A THU2 A THU3 A THU3 A 40405S031 GOLD A GOLD A GOLD A GOLD A GOL2 A 50103M108 TIDB A TIDB A TIDB A TIDB A TID1 A TID2 A 50134M001 DARR A DARW A DARW A DARR A DARW A DARW A 50107M004 YAR2 A YAR1 A YAR1 A 66010M001 DAVR A DAV1 A DAV1 A DAV1 A DAV1 A 14201M010 WTZR A WTZR A WTZR A WTZR A MISSING IERS DOMES NOS AT mit: ----------------------------- BISH A --