Subject: [IGSREPORT-15005] Wk 1421 IGS SINEX Combination
From: "Martin Bourassa" <[email-masked]>
To: [email-masked]
Cc: [email-masked], [email-masked]
Date: Thu, 19 Apr 2007 00:36:40 -0400

******************************************************************************
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: [email-masked]
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