Subject: [IGSREPORT-11162] Wk 1261 IGS SINEX Combination (Resubmitted)
From: Remi Ferland <[email-masked]>
To: [email-masked] (IGS Report Mailbox)
Cc: [email-masked] (gendt), [email-masked] (nisn)
Date: Fri, 26 Mar 2004 10:26:25 EST

******************************************************************************
IGS Electronic Report   26 Mar 08:02:01 PST 2004      Message Number 11162
******************************************************************************

Author:  R. Ferland , D. Hutchison

===================================================================
 Combined AC  solution for GPS week 1261 -  04/03/07 to 04/03/13
===================================================================

1) Contacts:
============
e-mail: [email-masked]
Tel:    (613)995-4002 (Remi Ferland)
Fax:    (613)995-3215

2) Products:
============
All products are available at cddisa.gsfc.nasa.gov

igs04P1261.snx  Combined weekly solution. 
igs04P1261.erp  Daily Earth Rotation Parameters.
igs00p02.erp    Daily Earth Rotation Parameters since GPS week 1013.
IGS04P10.snx    Cumulative Combined solution. 
igs04P1261.sum  This file.
igs04P1261.itr  Residuals between the AC and the IGb00 for the RF stations. 
igs04P1261.res  Residuals between the AC and the Combined Weekly solution.
IGS04P10.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 = igs04P1261.snx
     IGS = IGS04P10.snx

     MIT = mit12617.snx
     mit = mit1261g.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:
             cod  MANZ
             sio  CONZ  MANZ
             gfz  BAKO

     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.

     4-3) This message was resubmitted, because the original one
     could not be delivered by the local mail system.

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 |  76 |        26.365 |         1.000 |        26.366 |
emr            |  44 |  31 |         4.412 |         1.000 |         4.412 |
gfz            | 173 |  85 |         8.958 |         1.000 |         8.959 |
jpl            | 116 |  60 |         5.156 |         1.000 |         5.156 |
MIT            | 142 |  70 |         2.777 |         1.000 |         2.778 |
ngs            |  79 |  79 |        26.335 |         1.000 |        26.336 |
sio            | 175 |  87 |         2.128 |         1.000 |         2.128 |
mit            | 263 |  88 |         1.701 |         1.000 |         1.701 |
ncl            | 158 |  88 |         2.852 |         1.034 |         2.950 |
igs            | 243 |  92 |         1.400 |         0.998 |         1.398 |
---------------|-----|-----|---------------|---------------|---------------|

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         76 |     0.4      0.0     -1.0 |     2.4      2.4      8.3 |
emr         31 |    -0.7     -0.1      0.2 |     2.5      2.9      9.9 |
gfz         85 |    -0.1      0.0     -0.5 |     2.6      2.6      8.7 |
jpl         60 |    -0.7     -0.1      0.5 |     2.8      2.8     10.4 |
MIT         70 |     0.1      0.2      0.1 |     2.2      2.5      9.3 |
ngs         78 |    -0.1      0.8     -0.6 |     7.4      5.8     12.1 |
sio         87 |    -0.5      0.3      0.0 |     3.0      3.5      7.8 |
mit         88 |    -0.0      0.2      0.0 |     2.4      2.7      8.6 |
ncl         88 |     0.3     -0.1      0.1 |     2.5      2.9      7.3 |
---------------|---------------------------|---------------------------|

5-2-2)    With respect to the Weekly Combination:
=================================================

Solution       |  Weighted Average (mm)    |          STD (mm)         |
               |      N        E        U  |      N        E        U  |
---------------|---------------------------|---------------------------|
cod        159 |     0.2     -0.1     -0.2 |     1.5      1.1      5.6 |
emr         44 |    -0.5      0.1      0.3 |     1.7      1.5      5.3 |
gfz        169 |    -0.1     -0.2      0.2 |     1.3      1.4      5.2 |
jpl        113 |    -1.1     -0.2      0.4 |     1.6      1.6      8.1 |
MIT        130 |     0.5      0.1     -0.3 |     1.3      1.6      4.2 |
ngs         78 |    -0.0      0.7     -0.3 |     6.7      5.6     10.4 |
sio        162 |    -0.1      0.3      0.3 |     1.8      2.4      3.8 |
mit        234 |    -0.0      0.0     -0.2 |     0.6      0.8      3.6 |
ncl        155 |     0.3     -0.1      0.5 |     0.5      0.7      3.0 |
ITR         91 |     0.2      0.1     -0.5 |     2.2      2.2      7.1 |
---------------|---------------------------|---------------------------|

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.2      0.1 |     2.4      2.5      7.8 |
emr         41 |    -0.9     -0.2     -0.5 |     2.3      2.8      8.6 |
gfz        157 |    -0.6     -0.2      0.5 |     2.5      2.4      8.4 |
jpl        100 |    -1.4     -0.4      1.2 |     2.7      2.9      9.3 |
MIT        113 |     0.2     -0.1      0.1 |     2.1      2.7      7.5 |
ngs         79 |    -0.3      0.7     -0.3 |     7.1      5.7     11.8 |
sio        150 |    -0.4      0.2      0.9 |     2.6      3.2      6.8 |
mit        181 |    -0.3     -0.0      0.4 |     2.2      2.7      7.4 |
ncl        145 |     0.1     -0.3      0.6 |     2.3      2.6      6.4 |
igs        185 |    -0.4      0.0      0.4 |     1.8      2.0      6.0 |
ITR         96 |    -0.2      0.1     -0.2 |     0.6      0.5      1.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.017    0.006    0.024      8.1     14.8      9.8     2.65 |
emr            |  -0.011   -0.021    0.004      1.9      9.8      3.1     3.58 |
gfz            |   0.036    0.030    0.032     -1.8      2.8     -1.9     1.69 |
jpl            |   0.106   -0.033    0.063     -7.2      7.0      9.6     1.00 |
MIT            |   0.141   -0.008   -0.017      0.4      2.4      8.1     4.07 |
ngs            |   0.127   -0.215    0.036     13.0      6.2    -32.6     2.23 |
sio            |   0.021    0.004    0.084      6.0     -2.1     43.1     2.56 |
mit            |  -0.010   -0.000    0.001      0.1     -0.3     -0.9     2.81 |
ncl            |   0.003   -0.038   -0.042      2.2      1.6     -2.5     2.60 |
---------------|---------------------------------------------------------------|

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.046    0.051      1.2      1.2      1.2     0.19 |
emr            |   0.068    0.070    0.082      1.8      1.8      1.8     0.28 |
gfz            |   0.043    0.043    0.048      1.1      1.1      1.1     0.17 |
jpl            |   0.050    0.050    0.056      1.3      1.3      1.3     0.20 |
MIT            |   0.047    0.047    0.052      1.2      1.2      1.2     0.19 |
ngs            |   0.044    0.044    0.051      1.2      1.2      1.2     0.18 |
sio            |   0.042    0.043    0.048      1.1      1.1      1.1     0.17 |
mit            |   0.042    0.043    0.048      1.1      1.1      1.1     0.17 |
ncl            |   0.042    0.042    0.048      1.1      1.1      1.1     0.17 |
---------------|---------------------------------------------------------------|

5-4)      Geocenter:
====================

               |      X        Y        Z  |    s X      s Y      s Z  |
               |    (mm)     (mm)     (mm) |    (mm)     (mm)     (mm) |
---------------|---------------------------|---------------------------|
igs            |     9.0      8.4      8.1 |     1.8      1.6      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            |     ---      -33.1       0.3       ---      -3.9       --- |
emr            |     ---        2.6      -0.1      -8.0      -2.4       1.9 |
gfz            |     ---       24.4       3.0       5.2       4.3      -2.5 |
jpl            |     ---       36.1       1.5     -21.5      -2.0      15.7 |
MIT            |     ---       -1.2       2.5       1.1       0.7      -0.4 |
ngs            |     ---      -97.3      12.7       ---      25.5       --- |
sio            |     ---       -4.1      -3.2      -5.1      -2.6      -2.7 |
mit            |     ---       -9.5      -0.6       0.4       3.2       0.3 |
ncl            |     ---       -4.7       1.2       ---      -5.5       --- |
ITR            |     ---        ---       ---       ---       ---       --- |
BuA            |     ---      -25.0       0.1      -0.4      31.7      -1.7 |
---------------|------------------------------------------------------------|

5-5-2)    STD :
===============

Solution       |      UT       LOD       XPO      XPOR       YPO      YPOR  |
               |      us        us    .01mas    .01ma/d   .01mas    .01ma/d |
---------------|------------------------------------------------------------|
cod            |     ---        9.3       1.6       ---       1.6       --- |
emr            |     ---       13.7       2.8      17.5       3.6      10.5 |
gfz            |     ---        3.7       3.3       9.0       6.1       9.1 |
jpl            |     ---       27.7       3.0      15.3       3.5      11.7 |
MIT            |     ---        6.9       2.3       5.3       3.1      12.6 |
ngs            |     ---       36.2       9.0       ---      19.5       --- |
sio            |     ---       11.5       4.2      10.0       2.9       9.7 |
mit            |     ---        6.7       2.4       4.9       3.1       7.2 |
ncl            |     ---        0.9       1.5       ---       1.5       --- |
ITR            |     ---        ---       ---       ---       ---       --- |
BuA            |     ---       10.8       3.4      12.4       6.3      14.5 |
---------------|------------------------------------------------------------|

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    cod12617.snx        |YAKT  A ----|   11.5   -19.7     0.6|
                                         |            |                       |
IGS03P33_RS99.snx    gfz12617.snx        |YAKT  A ----|   14.6   -22.2   -11.7|
                                         |            |                       |
IGS03P33_RS99.snx    jpl12617.snx        |IISC  A ----|   -2.6     2.2  -166.3|
                                         |            |                       |
IGS03P33_RS99.snx    jpl12617.snx        |YAKT  A ----|    2.5   -15.2    15.3|
                                         |            |                       |
gfz12617.snx         jpl12617.snx        |ASPA  A ----|   -3.9    -0.2   -71.3|
                                         |            |                       |
gfz12617.snx         ngs12617.snx        |KOKB  A    2|   -1.5    13.7    56.7|
                                         |            |                       |
jpl12617.snx         ncl1261g.snx        |ASPA  A ----|    4.8    -2.6    52.4|
                                         |            |                       |
ngs12617.snx         ncl1261g.snx        |KOKB  A    2|    2.5   -11.5   -52.0|
                                         |            |                       |
IGS04P09.snx         gfz12617.snx        |OBET  A ----|   -0.6    -3.5   -69.3|
                                         |            |                       |
IGS04P09.snx         gfz12617.snx        |BAY2  A ----|    1.2    -4.0    59.4|
IGS04P09.snx         gfz12617.snx        |ZWEN  A ----|  -11.1    -6.1     5.6|
                                         |            |                       |
IGS04P09.snx         MIT12617.snx        |NETT  A ----|    6.6   -16.3   -23.9|
                                         |            |                       |
IGS04P09.snx         mit1261g.snx        |MANZ  A    2| -198.4   388.1    11.1|
                                         |            |                       |
IGS04P09.snx         mit1261g.snx        |BAY2  A ----|    0.8    -1.4    54.2|
IGS04P09.snx         mit1261g.snx        |OBET  A ----|    0.5    -3.5   -74.1|
                                         |            |                       |
IGS04P09.snx         igs04P1261.snx      |YSSK  A ----|    8.7    -3.9     3.6|
IGS04P09.snx         igs04P1261.snx      |ZWEN  A ----|   -9.0    -6.7     4.2|
-----------------------------------------|------------|-----------------------|

5-7)      Conflicts:
====================

IERS DOMES NO    MIT    cod    emr    gfz    jpl    ngs    sio    IGS    igs
============= ====== ====== ====== ====== ====== ====== ====== ====== ======
    13407S012 MAD2 A               MADR A MAD2 A        MAD2 A MADR A MADR A
                                                        MADR A
    43001M002 THU3 A THU2 A        THU3 A        THU3 A THU2 A THU3 A THU3 A
    40405S031 GOL2 A GOLD A GOL2 A GOLD A        GOLD B GOL2 A GOLD A GOLD A
    50103M108 TID1 A        TID2 A TIDB A TIDB A TIDB A TID2 A TIDB A TIDB A
              TID2 A               TID2 A               TIDB A
              TIDB A
    50134M001 DARR A DARR A        DARW A DARR A DARW A DARW A DARW A DARR A
              DARW A                      DARW A                      DARW A
    50107M004 YAR2 A               YAR2 A               YAR2 A YAR1 A YAR2 A
                                                               YAR2 A 
    66010M001 DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A DAV1 A
              DAVR A DAVR A                                           DAVR A
    22306M002        IISC A IISC A IISC A IISC A IISC B IISC A IISC A IISC A
    12313M002 MDVJ A IRKJ A               IRKJ A
    10422M001        KIR0 A                                           KR0G A
    10023M001                      LROC A                      LROC A LROC A
 
MISSING IERS DOMES NOS AT MIT:
-----------------------------
  HLFX A  PALM A  PUR3 A  PUR2 A  ELEN A  UPO1 A  OUSD A  HORN A  NETT A  CNCL A
  QUAR A  KARA A  SAMP A  MAR6 A  IENG A
 
MISSING IERS DOMES NOS AT gfz:
-----------------------------
  BISH A  PALM A
 
MISSING IERS DOMES NOS AT sio:
-----------------------------
  KULU A  PUR2 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
============= ====== ====== ====== ======
    13407S012 MAD2 A MAD2 A MADR A MADR A
              MADR A 
    43001M002 THU2 A THU3 A THU3 A THU3 A
              THU3 A 
    40405S031 GOLD A GOL2 A GOLD A GOLD A
              GOL2 A 
    50103M108 TID2 A TID2 A TIDB A TIDB A
              TIDB A TIDB A
              TID1 A 
    50134M001 DARR A DARR A DARW A DARR A
              DARW A DARW A        DARW A
    50107M004 YAR2 A YAR2 A YAR1 A YAR2 A
                            YAR2 A 
    12313M002 IRKJ A 
              MDVJ A 
    66010M001 DAV1 A DAV1 A DAV1 A DAV1 A
              DAVR A               DAVR A
    10422M001 KIR0 A               KR0G A
 
MISSING IERS DOMES NOS AT mit:
-----------------------------
  BISH A
--