From igs-rf at ign.fr Wed Sep 29 06:03:50 2010 From: igs-rf at ign.fr (igs-rf at ign.fr) Date: Wed, 29 Sep 2010 15:03:50 +0200 (CEST) Subject: [IGSREPORT-19011] Wk 1601 IGS SINEX Combination Message-ID: <20100929130350.5F0CDD4E7C8@newsmtp.ign.fr> Authors: B. Garayt, P. Rebischung =============================================================== Combined AC solution for GPS week 1601 - 10/09/12 to 10/09/18 =============================================================== 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 igs10P1601.snx Combined weekly solution igs10P1601.erp Daily Earth Rotation Parameters igs10P1601.sum This file igs10P1601.itr Residuals between AC solutions and IGS05 for the RF stations igs10P1601.res Residuals between AC solutions and combined weekly solution IGS10P37.snx Cumulative combined solution IGS10P37.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 = mit16017.snx mig = mit1601g.snx ITR = IGS05.snx igs = igs10P1601.snx IGS = IGS10P37.snx The following solutions are not used in the combination but only for comparison: mig ncl Solutions with manual station deletion: jpl NURK 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. 4-4) Starting with week 1598, the first iterative combinations are performed using AC solutions minimally constrained over the 7 Helmert parameters. (Until week 1597, they were constrained over the 3 rotations only.) The final combination is however still performed with AC solutions minimally constrained in rotations only. The combination is hardly affected by this change except that it is now more sensitive to outliers. 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 | 230 | 60 | 6.634 | 25.728 | emr | 82 | 42 | 2.338 | 7.303 | esa | 131 | 35 | 2.839 | 6.089 | gfz | 198 | 57 | 5.845 | 21.086 | grg | 123 | 45 | 5.706 | 9.249 | jpl | 123 | 44 | 5.146 | 12.967 | mit | 255 | 63 | 0.594 | 2.945 | ngs | 215 | 63 | 13.667 | 41.500 | sio | 266 | 61 | 0.730 | 2.653 | mig | 422 | 67 | ------ | 2.347 | ncl | 222 | 65 | ------ | 3.842 | igs | 397 | 67 | ------ | ------ | ---------------|-----|-----|---------------|---------------| 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 60 | -0.0 0.1 -0.1 | 3.7 3.9 9.4 | emr 42 | 0.3 0.3 0.0 | 3.7 3.7 8.9 | esa 35 | 0.1 -0.1 0.0 | 4.1 4.2 10.0 | gfz 57 | 0.8 0.2 -0.3 | 3.9 3.9 8.8 | grg 45 | 1.5 0.2 0.1 | 3.4 4.2 8.8 | jpl 44 | 1.2 0.0 -0.1 | 3.0 3.6 9.6 | mit 63 | 0.9 0.4 0.2 | 3.5 3.7 9.4 | ngs 63 | 0.3 0.2 0.1 | 3.5 3.4 9.0 | sio 61 | 0.8 0.2 0.2 | 3.4 3.7 8.8 | mig 67 | 0.0 0.1 0.2 | 3.4 3.6 9.2 | ncl 65 | 0.6 0.2 -0.0 | 3.5 3.4 9.0 | ---------------|---------------------------|---------------------------| 5-2-2) With respect to the Weekly Combination: ================================================= Solution | Weighted Average (mm) | STD (mm) | | N E U | N E U | ---------------|---------------------------|---------------------------| cod 220 | 0.0 -0.0 0.0 | 0.8 0.8 2.5 | emr 82 | 0.0 0.0 0.0 | 1.8 1.0 3.8 | esa 127 | 0.1 0.0 0.2 | 2.2 1.2 4.0 | gfz 193 | 0.1 0.0 -0.0 | 0.9 1.0 3.0 | grg 121 | 0.7 -0.0 0.1 | 2.7 2.0 5.6 | jpl 120 | -0.0 0.0 0.2 | 1.4 2.3 4.5 | mit 238 | 0.1 0.1 -0.0 | 0.5 0.6 1.5 | ngs 214 | -0.2 -0.0 -0.0 | 1.0 1.2 3.0 | sio 256 | -0.0 -0.1 0.0 | 1.2 0.8 2.1 | mig 396 | -0.2 0.0 -0.0 | 0.3 0.4 0.8 | ncl 222 | -0.1 0.0 0.0 | 0.7 0.7 1.0 | ITR 67 | -0.1 -0.2 -0.4 | 3.4 3.7 9.2 | ---------------|---------------------------|---------------------------| 5-2-3) With respect to the Cumulative Combination: ===================================================== Solution | Weighted Average (mm) | STD (mm) | | N E U | N E U | ---------------|---------------------------|---------------------------| cod 214 | 0.1 -0.0 -0.2 | 2.2 2.0 6.1 | emr 82 | 0.3 0.1 -0.1 | 2.9 1.9 5.5 | esa 119 | -0.0 -0.0 0.0 | 3.3 2.2 7.3 | gfz 193 | 0.3 0.0 -0.2 | 2.3 2.4 6.1 | grg 118 | 0.7 0.0 -0.0 | 3.2 2.4 6.5 | jpl 116 | -0.1 -0.0 0.0 | 3.7 3.6 6.4 | mit 232 | 0.1 0.1 -0.1 | 2.2 2.1 5.4 | ngs 207 | -0.1 0.0 -0.2 | 2.3 2.4 5.7 | sio 243 | 0.2 -0.0 -0.2 | 2.1 2.1 4.9 | mig 359 | -0.2 -0.0 -0.2 | 2.6 2.5 5.8 | ncl 221 | 0.1 -0.0 -0.2 | 2.2 2.0 5.3 | igs 360 | 0.0 -0.0 -0.2 | 2.5 2.4 5.7 | ITR 115 | -0.9 -0.1 -0.1 | 3.1 3.5 12.6 | ---------------|---------------------------|---------------------------| 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.016 -0.010 -0.033 0.7 1.5 0.7 -0.27 | emr | -0.000 -0.010 0.004 6.1 7.5 -3.2 0.05 | esa | 0.055 0.011 -0.044 3.9 3.5 8.4 -0.17 | gfz | 0.032 0.010 -0.051 -0.1 1.8 0.5 -0.14 | grg | 0.035 0.060 -0.015 3.5 12.4 4.8 -0.18 | jpl | -0.002 -0.043 0.001 3.8 9.2 -4.5 -0.00 | mit | 0.017 0.020 -0.031 4.2 3.8 9.5 -0.28 | ngs | 0.081 -0.099 0.021 6.2 4.8 5.5 -0.16 | sio | 0.103 -0.110 -0.007 5.3 6.5 -5.9 -0.14 | mig | 0.019 0.019 0.010 -0.8 2.0 0.9 -0.25 | ncl | 0.027 0.206 -0.118 -4.6 0.6 0.8 -0.21 | igs | 0.000 0.000 0.000 0.0 0.0 0.0 -0.15 | ---------------|---------------------------------------------------------------| 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.033 0.035 0.039 0.9 0.9 0.9 0.14 | emr | 0.035 0.038 0.043 1.0 1.0 0.9 0.15 | esa | 0.046 0.049 0.056 1.3 1.2 1.2 0.19 | gfz | 0.033 0.035 0.037 0.9 0.9 0.9 0.14 | grg | 0.038 0.040 0.044 1.0 1.0 1.0 0.16 | jpl | 0.040 0.044 0.045 1.1 1.1 1.1 0.17 | mit | 0.033 0.034 0.037 0.9 0.9 0.9 0.13 | ngs | 0.032 0.033 0.035 0.9 0.8 0.8 0.13 | sio | 0.031 0.032 0.035 0.8 0.8 0.8 0.13 | mig | 0.031 0.032 0.035 0.8 0.8 0.8 0.13 | ncl | 0.031 0.032 0.035 0.8 0.8 0.8 0.13 | igs | 0.000 0.000 0.000 0.0 0.0 0.0 0.10 | ---------------|---------------------------------------------------------------| 5-4) Geocenter: ==================== | X Y Z | s X s Y s Z | | (mm) (mm) (mm) | (mm) (mm) (mm) | ---------------|---------------------------|---------------------------| igs | 4.5 6.4 4.5 | 1.3 1.1 3.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 | --- --- --- --- --- --- | emr | --- -7.5 0.8 18.7 5.1 28.8 | esa | --- 12.1 -3.4 16.9 -2.2 18.2 | gfz | --- -3.4 -1.5 -0.5 -1.4 -5.9 | grg | --- --- -5.8 1.8 0.4 -2.0 | jpl | --- 0.0 0.0 11.2 -1.8 24.2 | mit | --- 0.5 0.7 -9.1 0.1 -2.4 | ngs | --- -10.5 -1.9 13.3 2.3 0.8 | sio | --- 5.9 0.6 18.8 -0.1 -5.0 | mig | --- 3.6 0.2 17.7 -1.3 -3.1 | ncl | --- --- --- --- --- --- | BuA | --- -8.4 -0.4 14.0 2.5 9.4 | ---------------|------------------------------------------------------------| 5-5-2) STD: ============== Solution | UT LOD XPO XPOR YPO YPOR | | us us .01mas .01ma/d .01mas .01ma/d | ---------------|------------------------------------------------------------| cod | --- --- --- --- --- --- | emr | --- 11.3 2.9 8.5 2.4 19.3 | esa | --- 11.6 2.2 6.8 1.4 12.5 | gfz | --- 9.7 1.3 7.8 1.3 13.3 | grg | --- --- 3.4 11.7 2.6 8.1 | jpl | --- 14.8 2.0 9.7 1.8 20.6 | mit | --- 0.8 0.7 1.5 0.7 2.3 | ngs | --- 18.3 4.2 6.1 5.7 20.4 | sio | --- 5.9 1.0 6.5 0.5 8.2 | mig | --- 6.0 1.1 6.2 0.6 7.6 | ncl | --- --- --- --- --- --- | BuA | --- 12.6 4.6 9.5 7.3 12.2 | ---------------|------------------------------------------------------------| 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 8 sigma. Solution Stage of rejection | Station | N E U | |code id soln| (mm) (mm) (mm)| -----------------------------------------|------------|-----------------------| | | | cod16017.snx Combination: Iter. 1 |HOB2 A 2| -2.3 -1.0 -7.9| esa16017.snx Combination: Iter. 1 |HOB2 A 2| -7.1 0.5 -4.5| grg16017.snx Combination: Iter. 1 |AREQ A 6| -11.0 2.6 -9.8| grg16017.snx Combination: Iter. 1 |COCO A 4| -12.0 -3.2 -4.5| jpl16017.snx Combination: Iter. 1 |GUAT A 1| 1.2 9.9 -13.2| jpl16017.snx Combination: Iter. 1 |PETS A 3| -1.5 -8.1 8.6| mit16017.snx Combination: Iter. 1 |DAV1 A 2| -0.4 2.1 0.2| mit16017.snx Combination: Iter. 1 |MAW1 A 1| -0.8 2.0 1.4| mit16017.snx Combination: Iter. 1 |SYOG A 1| -1.7 2.4 2.2| sio16017.snx Combination: Iter. 1 |UNSA A 2| 4.3 -0.8 0.2| sio16017.snx Combination: Iter. 1 |SANT A 4| 4.8 -1.2 -5.2| sio16017.snx Combination: Iter. 1 |AREQ A 6| 3.5 -1.5 0.8| sio16017.snx Combination: Iter. 1 |OUS2 A 3| 2.1 -2.5 0.4| sio16017.snx Combination: Iter. 1 |MAC1 A 7| 3.9 -1.0 2.0| sio16017.snx Combination: Iter. 1 |HOB2 A 2| 4.4 -0.3 7.0| esa16017.snx Combination: Iter. 2 |TIXJ A 2| -6.1 -1.3 -11.1| sio16017.snx Combination: Iter. 2 |STR1 A 1| 2.7 -0.7 3.7| | | | igs10P1601.snx Stacking: Iter. 1 |ALRT A 1| -12.8 -3.3 -5.4| igs10P1601.snx Stacking: Iter. 1 |SIO3 A 12| 5.1 -13.1 -2.3| igs10P1601.snx Stacking: Iter. 1 |TASH A 1| -9.9 5.0 67.0| igs10P1601.snx Stacking: Iter. 2 |METS A 1| 12.1 -3.7 -17.7| | | | -----------------------------------------|------------|-----------------------|