Subject: [IGSREPORT-6139] Wk 1009 Preliminary IGS SINEX Combination
From: R. Ferland , D. Hutchison
Date: Thu, 3 Jun 1999 14:20:00 PDT

******************************************************************************
IGS Electronic Report    Thu Jun  3 14:20:00 PDT 1999      Message Number 6139
******************************************************************************

Author: R. Ferland , D. Hutchison
Subject: Wk 1009 Preliminary IGS SINEX Combination

===================================================================
 Combined AC  solution for GPS week 1009 -  99/05/09 to 99/05/15
===================================================================

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

igs99P1009.snx  Combined weekly solution. 
igs99P1009.erp  Daily Earth Rotation Parameters.
IGS99P19.snx    Cumulative Combined solution. 
igs99P1009.sum  This file.
igs99P1009.itr  Residuals between the AC and the ITRF96 for the RF stations.
igs99P1009.res  Residuals between the AC and the Combined Weekly solution.
IGS99P19.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
     ITRF96 Reference Frame (RF) stations at the current epoch.
     ftp://igscb.jpl.nasa.gov/igscb/station/coord/ITRF96_IGS_RS47.SNX.Z
3-4) The AC are compared to ITRF, 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 ITRF96 of the current
     epoch using a 7-parameters transformation. The singular cumulative
     combined solution is aligned to the ITRF96 of the reference epoch
     using a 14-parameters transformation and the 47 ITRF station positions
     and velocities.

     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:
===========

     igs = igs99P1009.snx
     IGS = IGS99P19.snx

     JPL = jpl1009g.snx
     jpl = jpl10097.snx

     ITR = ITRF96_IGS_RS47.SNX
     itr = ITRF97 with the current 47 Reference Frame stations plus
           AUCK BAHR CAS1 CHAT MCM4.
     JPL mit ncl itr are used to quality control the weekly combination.
     emr gfz ngs sio excluded from the geocenter estimation (Sec. 5-4).

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            | 112 |  40 |        37.068 |         1.005 |        37.255 |
emr            |  43 |  27 |         5.537 |         1.000 |         5.536 |
esa            |  58 |  38 |         7.502 |         0.998 |         7.489 |
gfz            |  60 |  33 |        19.662 |         1.002 |        19.698 |
jpl            |  57 |  25 |         3.956 |         1.003 |         3.966 |
ngs            |  71 |  43 |        46.582 |         1.006 |        46.884 |
sio            | 121 |  40 |         2.334 |         1.001 |         2.337 |
JPL            | 140 |  41 |         1.880 |         1.010 |         1.899 |
mit            | 142 |  42 |         3.611 |         0.999 |         3.606 |
ncl            |  82 |  40 |         2.902 |         1.000 |         2.902 |
itr            |  52 |  47 |        10.000 |         0.341 |         3.411 |
---------------|-----|-----|---------------|---------------|---------------|

5-2)      Solution Residuals weighted Average & R.M.S. :
========================================================

5-2-1)    With respect to ITRF96 Reference Frame stations:
=========================================================

Solution       |  Weighted Average (mm)    |        R.M.S. (mm)        |
               |      N        E        U  |      N        E        U  |
---------------|---------------------------|---------------------------|
cod            |     0.6     -1.6     -0.1 |     5.5      8.7     12.6 |
emr            |    -4.9     -0.7     -0.1 |     9.3     13.3     17.3 |
esa            |    -1.3     -2.9      1.0 |     4.7     11.7     19.7 |
gfz            |    -1.1     -1.2      0.5 |     5.9      6.4     11.9 |
jpl            |     0.5     -0.5      1.8 |     5.7      7.1     14.6 |
ngs            |    -0.7     -0.4      1.2 |    10.6     12.5     15.8 |
sio            |    -3.3     -1.2      0.0 |     6.2      8.4     12.4 |
JPL            |    -1.8     -0.9      2.3 |     4.4      5.0     13.8 |
mit            |    -0.4     -0.8      0.5 |     3.7      5.9      6.7 |
ncl            |     0.0     -1.0      0.2 |     3.9      5.0      6.1 |
itr            |    -0.3      0.0     -0.4 |     4.2      6.1      9.1 |
---------------|---------------------------|---------------------------|

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

Solution       |  Weighted Average (mm)    |        R.M.S. (mm)        |
               |      N        E        U  |      N        E        U  |
---------------|---------------------------|---------------------------|
cod            |     3.2     -1.0     -0.1 |     4.6      3.5      8.9 |
emr            |    -2.9      0.1      0.8 |     5.4     10.0     14.1 |
esa            |     1.8     -4.2     -0.9 |     5.4     10.1     17.1 |
gfz            |     1.6     -1.0      0.2 |     2.7      3.3      6.6 |
jpl            |     4.9     -1.1      2.4 |     5.6      3.4      8.8 |
ngs            |     1.3     -1.1      0.7 |     8.1     13.1     15.2 |
sio            |    -0.6     -1.3     -0.9 |     2.4      3.9      7.1 |
JPL            |     0.9     -1.4      1.9 |     1.8      2.2      6.2 |
mit            |     2.2     -0.7      1.5 |     3.2      2.4      8.9 |
ncl            |     2.7     -0.8      0.4 |     3.0      1.7      5.7 |
itr            |     2.5      1.2     -1.0 |     4.7      7.3      9.6 |
ITR            |     2.3      1.8     -1.4 |     6.0      9.5     14.2 |
---------------|---------------------------|---------------------------|

5-2-3)    With respect to Cumulative Combination:
================================================

Solution       |  Weighted Average (mm)    |        R.M.S. (mm)        |
               |      N        E        U  |      N        E        U  |
---------------|---------------------------|---------------------------|
cod            |    -0.2     -0.9     -1.1 |     3.6      4.2     12.7 |
emr            |    -5.4      0.7     -1.7 |     7.1     10.7     18.8 |
esa            |    -1.2     -3.8     -2.7 |     5.6      9.0     19.2 |
gfz            |    -1.1     -0.4     -1.9 |     3.2      4.3     12.9 |
jpl            |     1.8     -0.4      0.0 |     3.6      4.4     12.7 |
ngs            |    -1.4     -1.4     -1.0 |     8.9     13.1     17.4 |
sio            |    -3.7     -1.2     -2.3 |     4.8      5.4     10.3 |
JPL            |    -2.4     -1.0      0.6 |     3.8      3.3     13.0 |
mit            |    -1.0     -0.4      0.5 |     2.7      3.3     12.9 |
ncl            |    -0.1     -0.5     -1.0 |     2.2      2.8     11.2 |
itr            |     0.3      0.4     -1.0 |     2.8      3.0      9.4 |
igs            |    -3.1      0.3     -1.1 |     4.2      4.2      9.8 |
ITR            |     0.4      0.9     -0.8 |     5.5      7.6     13.7 |
---------------|---------------------------|---------------------------|

5-3)      7-Parameters Transformation: AC ----> ITRF96_IGS_RS47.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.060   -0.009    0.002     13.7     -2.4     35.5     1.39 |
emr            |   0.111    0.005    0.001     -2.9    -54.1    312.2    -0.20 |
esa            |   0.008   -0.017    0.062     11.5     14.5      7.1     0.11 |
gfz            |   0.019    0.076   -0.043     -0.0      1.1      0.1    -0.06 |
jpl            |   0.066    0.017    0.060      6.3    -13.6     47.9     0.78 |
ngs            |  -0.140    0.248    0.141     -5.1     -7.7     48.1    -0.41 |
sio            |  -0.034    0.552    0.380     -4.9      5.3    145.4     0.97 |
JPL            |   0.076    0.139    0.058     -1.9      4.3      0.9     0.38 |
mit            |  -0.015    0.025    0.025     -0.1      0.4     -0.7     0.42 |
ncl            |   0.035   -0.018    0.054      1.4      1.3     -1.0     0.41 |
itr            |  -0.160    0.261    0.057     -0.5     -0.6     14.5    -1.43 |
---------------|---------------------------------------------------------------|

5-3-2)    Sigmas : 
==================

Solution       |    R X      R Y      R Z      T X      T Y      T Z      SCL  |
               |   (mas)    (mas)    (mas)     (mm)     (mm)     (mm)    (ppb) |
---------------|---------------------------------------------------------------|
cod            |   0.067    0.071    0.069      1.8      1.7      1.7     0.27 |
emr            |   0.079    0.078    0.087      2.1      2.1      2.0     0.32 |
esa            |   0.066    0.069    0.071      1.7      1.7      1.7     0.27 |
gfz            |   0.076    0.075    0.077      2.0      2.0      1.9     0.30 |
jpl            |   0.083    0.087    0.083      2.1      2.1      2.1     0.33 |
ngs            |   0.066    0.068    0.068      1.7      1.7      1.7     0.26 |
sio            |   0.067    0.071    0.069      1.8      1.7      1.7     0.27 |
JPL            |   0.067    0.071    0.068      1.8      1.8      1.7     0.27 |
mit            |   0.065    0.067    0.069      1.7      1.7      1.7     0.26 |
ncl            |   0.068    0.073    0.069      1.8      1.8      1.7     0.27 |
itr            |   0.060    0.062    0.064      1.6      1.6      1.6     0.24 |
---------------|---------------------------------------------------------------|

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

               |      X        Y        Z  |    s X      s Y      s Z  |
               |    (mm)     (mm)     (mm) |    (mm)     (mm)     (mm) |
---------------|---------------------------|---------------------------|
igs            |     8.8      2.5     34.5 |     5.1      5.0      7.4 |
---------------|---------------------------|---------------------------|

5-5)      ERP Residuals Weighted Average & R.M.S. :
===================================================

5-5-1)    Weighted Average :
============================

Solution       |      UT       LOD       XPO      XPOR       YPO      YPOR  |
               |      us        us    .01mas    .01ma/d   .01mas    .01ma/d |
---------------|------------------------------------------------------------|
cod            |     ---        2.3       9.0       2.9      21.0       -.7 |
emr            |     ---       24.5      23.6     -19.5       5.6      17.3 |
esa            |     ---       -9.9      -6.1      -1.6     -14.2       4.4 |
gfz            |     ---        5.5      -4.5      -1.6       -.4      -2.8 |
jpl            |     ---        1.3       3.1       8.4      11.9      12.9 |
ngs            |     ---     -331.9     -26.6     -48.1      54.0       3.6 |
sio            |     ---        ---       ---       ---       ---       --- |
JPL            |     ---        ---       ---       ---       ---       --- |
mit            |     ---       -9.7     -20.6       1.9        .5       2.1 |
ncl            |     ---        ---       ---       ---       ---       --- |
itr            |     ---        ---       ---       ---       ---       --- |
ITR            |     ---        ---       ---       ---       ---       --- |
---------------|------------------------------------------------------------|

5-5-2)    R.M.S. :
==================

Solution       |      UT       LOD       XPO      XPOR       YPO      YPOR  |
               |      us        us    .01mas    .01ma/d   .01mas    .01ma/d |
---------------|------------------------------------------------------------|
cod            |     ---        6.7      10.1      13.9      22.2      23.5 |
emr            |     ---       27.3      24.1      34.3      10.0      27.0 |
esa            |     ---       14.8      11.1      18.7      15.6      15.0 |
gfz            |     ---       15.4       5.4       9.6       2.9       7.1 |
jpl            |     ---        5.9       7.4      18.6      14.4      29.2 |
ngs            |     ---      371.2      27.3      81.5      55.0      37.1 |
sio            |     ---        ---       ---       ---       ---       --- |
JPL            |     ---        ---       ---       ---       ---       --- |
mit            |     ---       10.7      20.7       8.8       3.4      12.6 |
ncl            |     ---        ---       ---       ---       ---       --- |
itr            |     ---        ---       ---       ---       ---       --- |
ITR            |     ---        ---       ---       ---       ---       --- |
---------------|------------------------------------------------------------|

5-6)      Outliers:
===================

Threshold set at 5 sigmas

 Solution             Solution           |  Station   |    N       E       U  |
              Minus                      |code id soln|   (mm)    (mm)    (mm)|
-----------------------------------------|------------|-----------------------|
                                         |            |                       |
ITRF96_IGS_RS47.snx  esa10097.snx        |MADR  A ----|   10.6   -43.2    33.8|
                                         |            |                       |
ITRF96_IGS_RS47.snx  jpl10097.snx        |MADR  A ----|   19.6   -43.2    22.5|
                                         |            |                       |
ITRF96_IGS_RS47.snx  jpl10097.snx        |MAC1  A ----|   -2.6    22.0     9.2|
                                         |            |                       |
ITRF96_IGS_RS47.snx  sio10097.snx        |MADR  A ----|   10.6   -38.1    19.0|
                                         |            |                       |
ITRF96_IGS_RS47.snx  JPL1009g.snx        |MAC1  A ----|   -2.4    25.0     6.6|
ITRF96_IGS_RS47.snx  JPL1009g.snx        |MADR  A ----|   12.0   -42.5    27.7|
                                         |            |                       |
ITRF96_IGS_RS47.snx  JPL1009g.snx        |DAV1  A ----|  -17.7   -16.7    -7.0|
ITRF96_IGS_RS47.snx  JPL1009g.snx        |GODE  A ----|    2.3   -14.6   -17.3|
                                         |            |                       |
ITRF96_IGS_RS47.snx  mit1009g.snx        |MADR  A ----|   12.3   -41.7    22.7|
                                         |            |                       |
ITRF96_IGS_RS47.snx  ncl1009g.snx        |MAC1  A ----|    0.0    20.4     1.7|
                                         |            |                       |
IGS99P18.snx         sio10097.snx        |PRDS  A ----|   -4.9     1.0   -61.6|
IGS99P18.snx         sio10097.snx        |YAKA  A ----|  -28.5     3.4   -10.6|
                                         |            |                       |
IGS99P18.snx         JPL1009g.snx        |WES2  A ----|    9.1   -11.0     7.6|
IGS99P18.snx         JPL1009g.snx        |YAKA  A ----|  -25.4     5.4   -17.6|
                                         |            |                       |
IGS99P18.snx         mit1009g.snx        |HERS  A ----|   -4.7   -16.9    -2.6|
IGS99P18.snx         mit1009g.snx        |YAKA  A ----|  -24.8     5.7   -13.3|
                                         |            |                       |
IGS99P18.snx         igs99P1009.snx      |BOGT  A ----|   -4.2    -1.2   -46.6|
IGS99P18.snx         igs99P1009.snx      |HERS  A ----|   -6.0    -9.9    -9.3|
IGS99P18.snx         igs99P1009.snx      |WES2  A ----|    7.7   -10.0     2.6|
IGS99P18.snx         igs99P1009.snx      |YAKA  A ----|  -25.9     5.7   -21.7|
-----------------------------------------|------------|-----------------------|

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

STATIONS WITH INCONSISTENT ANTENNA ECCENTRICITY INFORMATION:
(FILE: / STATION / L1 (UNE) / L2 (UNE) / ECCENTRICITY (UNE))
(UNE = Up North East (metres))

jpl1009g.snx: BAKO    .078   .000  -.003   .074  -.003  -.001  1.676   .000   .000
mit1009g.snx: BAKO    .078   .000  -.003   .074  -.003  -.001   .000   .000   .000
   - Unresolvable. No log file. Antenna parameters not listed in igs.snx.

jpl1009g.snx: BSHM    .000   .000   .000   .000   .000   .000   .000   .000   .000
mit1009g.snx: BSHM    .081   .000  -.001   .078  -.002   .004   .000   .000   .000
   - Unresolvable. No log file. Antenna parameters not listed in igs.snx.

jpl1009g.snx: SUTH    .000   .000   .000   .000   .000   .000   .119   .000   .000
mit1009g.snx: SUTH    .110   .000   .000   .128   .000   .000   .084   .000   .000
   - Both wrong. Should be L1/L2 offset, eccentricity 0.110, 0.128, 0.119 M. respectively.

jpl1009g.snx: TROM    .078   .000   .000   .096   .000   .000  2.369   .000   .000
mit1009g.snx: TROM    .110   .000   .000   .128   .000   .000   .000   .000   .000
   - Both wrong. Should be L1/L2 offset, eccentricity 0.078, 0.096, 2.475 M. respectively.

jpl1009g.snx: WUHN    .110   .000   .000   .128   .000   .000  2.355   .000   .000
mit1009g.snx: WUHN    .110   .000   .000   .128   .000   .000  2.320   .000   .000
   - Uncertain. Most recent log file supports an antenna height of 2.32 M., contrary to 2.355 M.
     given in igs.snx. This log file gives reference point as BCR instead of customary ARP.

jpl1009g.snx: BARB    .000   .000   .000   .000   .000   .000   .000   .000   .000
mit1009g.snx: GOUG    .110   .000   .000   .128   .000   .000   .003   .000   .000
jpl1009g.snx: MDVO    .000   .000   .000   .000   .000   .000   .146   .000   .000
jpl1009g.snx: MONP    .000   .000   .000   .000   .000   .000   .119   .000   .000
mit1009g.snx: VESL    .110   .000   .000   .128   .000   .000  -.011   .000   .000

FILE ncl1009g.snx MISSING THE SITE/GPS_PHASE_CENTER BLOCK.
DEFAULTING TO KNOWN L1/L2 OFFSETS AND ECCENTRICITIES.

 
ANTENNA PHASE CENTRE CONFLICTS IN FILE jpl1009g.snx
(MODEL / SER. NO. / L1 (UNE) / L2 (UNE))
(UNE = Up North East (metres))
   DORNE MARGOLIN B     -----    .078   .000   .000   .096   .000   .000
   DORNE MARGOLIN B     -----    .000   .000   .000   .000   .000   .000
 
SITES AFFECTED BY ABOVE CONFLICTS:
   Stations: NYAL 

 STATION NAME/POINT CODES CHANGED IN mit1009g.snx 
 ENFORCED FROM TEMPLATE FILE igs.snx .

    Old     New      IERS
Code/Pt Code/Pt  DOMES NO
------- -------  --------
YAKZ  A YAKA  A 12353M001

      1 Change  in all.

 STATION NAME/POINT CODES CHANGED IN jpl1009g.snx 
 ENFORCED FROM TEMPLATE FILE igs.snx .

    Old     New      IERS
Code/Pt Code/Pt  DOMES NO
------- -------  --------
YAKZ  A YAKA  A 12353M001

      1 Change  in all.

THE FOLLOWING ITRF96 FIDUCIAL STATIONS ARE MISSING FROM jpl1009g.snx :
   Stations: HART KERG MAS1

THE FOLLOWING ITRF96 FIDUCIAL STATIONS ARE MISSING FROM mit1009g.snx :
   Stations: HART KERG MAS1

THE FOLLOWING ITRF96 FIDUCIAL STATIONS ARE MISSING FROM ncl1009g.snx :
   Stations: BRAZ HART KERG MAD2 MAS1 OHIG TROM

SITE/ID BLOCK IERS DOMES NOS. VS STATION CODE/PT FIELDS IN COMPARED GNAAC FILES:
(ENTRIES APPEAR ONLY WHERE GNAACS DIFFER.)

IERS DOMES NO    JPL    MIT    NCL
============= ====== ====== ======
    40472S003 AMCT A              
              AMC2 A              
    42202M005 AREQ A AREQ A AREQ A
              ARE1 A              
    40405S031 GOL2 A GOLD A GOL2 A
              GOLD A GOL2 A GOLD A
              GOL1 A              
    30302M007 HARK A HARK A HARK A
              HAR1 A              
    13407S012 MAD2 A MAD2 A       
              MADR A MADR A       
    12309M002 MDVO A              
    12309M001        MDVO A       
    50103M108 TID2 A TIDB A TID2 A
              TIDB A TID2 A TIDB A
 
MISSING IERS DOMES NOS AT JPL:
-----------------------------
  ALIC A  BAKO A  BSHM A  ELAT A  KABR A  KARR A  KATZ A  MAW1 A  RAMO A  RIOG A
  TELA A  TOW2 A  URUM A
 
MISSING IERS DOMES NOS AT MIT:
-----------------------------
  KATZ A  KABR A

SITE/ID BLOCK IERS DOMES NOS. VS STATION CODE/PT FIELDS IN ABBREVIATED .ssc FILES:
(ENTRIES APPEAR ONLY WHERE ACS DIFFER.)

IERS DOMES NO    COD    NRC    ESA    GFZ    JPL    NOA    SIO
============= ====== ====== ====== ====== ====== ====== ======
    30302M007        HART B HART B HARK A HARK A HART B       
    10317M003 NYAL A               NYA1 A NYA1 A        NYA1 A
    10317M001        NYAL A NYAL A               NYAL A       
    10302M006 TROM A               TRO1 A               TRO1 A
    10302M003               TROM A               TROM A       
    12355M002 PETR A               PETP A PETP A PETP A PETP A
    40472S004 AMCT A                                    AMC2 A
    40405S031 GOLD A GOLD A GOLD A GOL2 A GOL2 A GOLD B GOL2 A
    42202M005 AREQ A AREQ B AREQ B AREQ A        AREQ B AREQ A
    50103M108 TIDB A TIDB A TIDB A TID2 A TID2 A TIDB A TID2 A
    13407S012               MADR A        MAD2 A MADR A MAD2 A
    22601M001               NTUS A                      NTUS A
    22601M002                                    NTUS A       
 
MISSING IERS DOMES NOS AT NOA:
-----------------------------
  BAKO A
 
MISSING IERS DOMES NOS AT SIO:
-----------------------------
  KATZ A  KABR A

[Mailed From: Remi Ferland <[email-masked]>]