^{1}

^{2}

^{3}

Since 1990, positioning technology has undergone a dramatic improvement in terms of accuracy and accessibility. Prior to 1990, national geodetic datums were considered to be static with fixed coordinates assigned to the datum origin by convention. Datum coordinates were dynamic only as a consequence of re-observation and re-adjustment, or localized disturbance. This strategy has been suitable for terrestrial surveys within stable continental areas. The accessibility of Precise Point Positioning “PPP” is now widespread with the provision of free online processing facilities by various geodetic agencies. Such ready accessibility to a global datum has serious implication for this technology with static geocentric datums. The Geocentric Datum of Egypt (ITRF1994 Epoch 1996 based on GPS Observation Campaign 1996), for example, is now offset by -42.0 to +32.0cm from the ITRF2008 Epoch 2015 due to the inexorable tectonic movements of the Nubian plate since the datum realization in 1996. Unless the local geocentric datums are kept in step with global terrestrial reference frames, discrepancies between the two will increase in magnitude over time and will become discernable, even using PPP. The PPP users nowadays are unaware of this datum offset and incorrectly assume that a coordinated datum monument is in error if there is found to be a disagreement greater than the level of accuracy of the positioning technique. In this study, to avoid the discrepancies between the applied datum and On-Line GNSS Services, a semi-kinematic datum and related deformation mode are proposed for Egypt. The result also confirmed that the PMM model was giving good performance for the other parts out of Egypt. The evaluation study shows the best performances for the Egyptian deformation model over the ITRF2008 PMM.

In 1992, an Egyptian Survey Authority “ESA” steering committee developed a plan for the creation of new datum for Egypt, with the following approach [

・ First, observing approximately 30 stations at approximately 200 km interval, covering all of Egypt, creating a High Accuracy Reference Network (HARN). Both high absolute and relative accuracies are required for these stations.

・ Second, establishing the Notational Agricultural Cadastral Network (NACN) relative to these 30 stations, covering the green area of Egypt (Nile Valley and the Delta) at 30 - 40 km intervals. This station spacing was selected to allow for further densification with single frequency receivers, see

・ Third, densifying this network at a station spacing of approximately 5 km for use as cadastral control at the governorate level.

・ Finally, replacing the existing Egyptian Mercator grid with a new modified UTM coordinate system.

Each station was observed for six sessions, every session was 6 hours with 30 seconds epoch interval. The observation time was planned to produce 1:10,000,000 (Order A) for HARN and 1:1,000,000 (Order B) for NACN relative network accuracy standard between stations. The results of analyzing both of them were defined in ITRF1994 Epoch 1996.

Many users of the geodetic system are non-technical users, for whom managing the dynamics of the datum represents a potential annoyance and complexity. ESA does not manage the dynamics of the geodetic system which enables other spatial datasets connected to it to be updated. The implementation of a “Continuously Operating Reference Stations (CORS)” network in Egypt in both ESA & NRIAG also represents a set of issues that need to be part of a semi-dynamic datum. This paper represents some of the implications and limitation for the users of geodetic and related datasets when implementing a semi-dynamic datum.

The realization of International Terrestrial Reference System is named (conventional) Terrestrial Reference Frame (TRF) and it is carried out through the coordinates of a set of points on the earth serving as reference points. A conventional TRF is defined as a set of physical points with precisely determined coordinates in a specific coordinate system that is the realization of an ideal Terrestrial Reference System [

necessary information of the datum definition of the Terrestrial Reference Frame (TRF) (origin, scale, and orientation), combination of reference frames provided by independent techniques has long been the standard method to implement global terrestrial reference frames.

In the classical sense, a geodetic datum is a reference surface, generally an ellipsoid of revolution of adopted size and shape, with origin, orientation, and scale defined by a geocentric terrestrial frame. Once an ellipsoid is selected, coordinates of a point in space can be given in Cartesian or geodetic (curvilinear) coordinates (geodetic longitude, latitude, and ellipsoid height). Furthermore, the classical concept of geodetic datum implies that a datum’s coordinates are fixed and does not change with time except for the effect of local tectonic motion. Thus, the coordinate frame of a geodetic datum should be somewhat attached to the plate and move with it in such a way that the coordinates of the points will not change as a consequence of plate rotation.

Two types of geodetic datum can be defined namely a static and kinematic geodetic datum. A static datum is thought of as a traditional geodetic datum where all sites are assumed to have coordinates which are fixed or unchanging with time. This is an incorrect assumption since the surface of the earth is constantly changing because of tectonic motion. The major limitation with a static geodetic datum arises from the ongoing divergence between ITRF and the fixed coordinates of a static geodetic network due to rigid plate motion. While the network may not necessarily be deforming internally to any significant degree, the lithospheric plate on which the network sits is moving as a rigid body over the Earth’s asthenospheric mantle. The impact of this deformation is noticeable where precision GNSS techniques are used to compute ITRF coordinates. Another limitation of a static datum arises from the processing of long GNSS baselines. If the static coordinates of a reference station are held fixed, rigid plate rotation of a long baseline will degrade the precision of the point computation as a function of time [

To show the drawback of the currently used Egypt static datum, namely ITRF1994 Epoch 1996, the Helwan IGS Station with coordinate values computed by IGS computing center in ITRF2008 Epoch 2015.274, was transferred by the published transformation parameters to ITRF1994 Epoch 1996 (http://itrf.ign.fr/trans_para.php). The related Reference Frame coordinate values, related velocity parameters and the differences are tabulated in

A kinematic datum (also referred to as a dynamic datum in literature and practice) such as the latest realization of the International Terrestrial Reference Frame ITRF2008 [

Reference Frame | HELWAN | ||
---|---|---|---|

IGS-ITRF2008 Epoch 2015.274 | 4,728,141.098 | 2,879,662.549 | 3,157,147.092 |

Published ITRF2008 Epoch 2005 | 4,728,141.309 | 2,879,662.406 | 3,157,146.932 |

Published Vx, Vy, Vz (m/year) | −0.0211 | 0.0143 | 0.016 |

ITRF1994 Epoch 1996 | 4,728,141.516 | 2,879,662.289 | 3,157,146.776 |

Differences | −0.418 | 0.260 | 0.316 |

absolute geodetic reference frame such as ITRF to be propagated to a static reference epoch, a process described in this paper as “epochwise”, in order to distinguish it from classical kinematic conformal transformation methods such as the 14 parameter (7 parameters and their rates of change) Bursa-Wolf transformation [

Global Navigation Satellite Systems (GNSS) inherently use an ITRF aligned reference frame for estimation of orbit ephemerides and derived positions of GNSS sensor locations. GNSS devices which use ITRF or closely aligned systems position users in agreement with the underlying kinematic frame. However, in practice there are a number of very significant drawbacks to a kinematic datum. Surveys undertaken at different epochs cannot be combined or integrated unless a deformation model is applied rigorously, or is embedded within the data, and the data are correctly time-tagged.

The limitations of a kinematic datum soon become apparent as kinematic and ground-fixed coordinates become misaligned. Kinematic coordinates increase the risk of litigation where legal definition of “fixed” coordinates is contested. They also increase the likelihood of errors in engineering surveys and damage to sub-surface infrastructure arising from misinterpretation or absence of epoch metadata.

A semi-kinematic datum is one where a deformation is modeled as an implicit component of the datum, thus enabling kinematic coordinates acquired by space geodetic techniques to be transformed consistently and accurately to a fixed reference epoch (static datum) over time. Many countries and regions which straddle major plate boundaries have adopted a semi-kinematic (or semi-dynamic) geodetic datum in order to prevent degradation of the datum as a function of time due to ongoing crustal deformation that is occurring within the country. The two major drivers for the adoption of a semi-kinematic datum in these countries have been:

・ The widespread adoption of precision GNSS techniques for surveying and positioning

・ The need to maintain consistency of coordinates to support combination and integration of spatial information acquired at different epochs [

High precision GNSS positioning and navigation is very rapidly highlighting the disparity between global kinematic reference frames such as ITRF and WGS84, and traditional static geodetic datum. The disparity is brought about by the increasingly widespread use of PPP and the sensitivity of these techniques to deformation of the Earth due to plate tectonics.

In order for surveys undertaken at different epochs to be combined or integrated within a kinematic datum, a deformation model has to be applied rigorously, or be embedded within the data. For the model to be applied correctly, all data has to be correctly time-tagged with the epoch of acquisition. A deformation model implies that the deformation is modeled with respect to an Earth-Centered Earth-Fixed (ECEF) reference frame consistent with the latest IERS Conventions, such as ITRF. A deformation model can be used to predict ITRF (absolute) site velocities at any location. The ITRF site velocity is then used to compute the displacement of the location between the epoch of measurement or acquisition and a reference epoch.

A deformation model can also be used for processing of long GNSS baselines by forward propagation of the ITRF coordinates of a CORS reference station from a reference epoch. This approach enables processing and network analysis to be accomplished wholly within ITRF thereby eliminating the effects of un-modeled plate rotation and deformation between the measurement and reference epochs. The deformation model can then be used in reverse propagation mode to compute the rover station coordinates at the reference epoch. A deformation model can be become an integral component of any semi-kinematic datum definition, with reference epoch (static) coordinates of the datum realized by epoch projection.

In order to estimate a plate motion model consistent with ITRF2008, detailed analyses of the ITRF2008 velocity field were undertaken by [_{i}, located on plate p:

As it is shown in

Plate | N.S. | (mas/a) | WRMS | ||||
---|---|---|---|---|---|---|---|

E | N | ||||||

ARAB | 4 | 1.202 ± 0.082 | −0.054 ±0.100 | 1.485 ± 0.063 | 0.531±0.027 | 0.23 | 0.15 |

EURA | 69 | −0.083 ± 0.008 | −0.534 ± 0.007 | 0.750 ± 0.008 | 0.257 ± 0.002 | 0.34 | 0.28 |

NUBI | 11 | 0.095 ± 0.009 | −0.598 ± 0.007 | 0.723 ± 0.009 | 0.262 ± 0.003 | 0.26 | 0.35 |

SOMA | 3 | −0.080 ± 0.028 | −0.745 ± 0.030 | 0.897± 0.012 | 0.325 ± 0.007 | 0.28 | 0.21 |

In 2006, NRIAG started the establishment of the Egyptian Permanent GPS Network (EPGN). Basically, the site selection was aimed to cover geographically all the Egyptian territory but also considering the tectonic setting of Egypt. As shown in

To calculate the velocity field referred to ITRF2008, in addition to the 16 EPGN stations [

To convert the aforementioned plate rotation rates, given by PMM in

The Pole Rotation Rate can be computed by:

The Pole Latitude:

The Pole Longitude:

The velocity Estimation by PMM can be obtained by:

Parameter | Egypt DM | ITRF2008 PMM |
---|---|---|

Euler pole Latitude (˚) | 50.30 | 50.055 |

Euler pole Longitude (˚) | −80.51 | −80.973 |

Angular velocity (˚/myr) | 0.267 | 0.262 |

PPP, global RTK and GNSS post-processing services are now used extensively to provide realizations of ITRF and WGS84 globally with a precision of a few centimeters. Unless these instantaneous realizations are transformed to a static or semi-kinematic datum using a suitable deformation model, repeat surveys using these techniques will result in datum divergence as a function of time arising from the effects of un-modeled tectonic plate motion. Transformations from kinematic ITRF to a static datum are conventionally done by either using the site velocity (computed either from ITRF definition or from a plate motion model).

a) Using Published IGS Displacement Values

The solution of CSRS-PPP SW is now used extensively to provide realizations of ITRF globally with a precision of a few centimeters. The CSRS-PPP’s service computes the absolute position for Helwan and Rabat, where both of them is considered as one of IGS stations that were used in computing the ITRF2008m IGS stations, at Epoch 2014.907 and Epoch 2015.324 respectively and the results are:

PPP solution for Helwan at ITRF2008 Epoch 2014.907

The published coordinates/standard deviations (95%) for the PHLW3310 RINEX file are as follow:

Latitude (ITRF08): 29 51 41.5732 dms/0.002 m, Longitude (ITRF08): 31 20 36.2409 dms/0.004 m, Ellipsoidal Height (ITRF08): 148.729 m/0.009 m Cartesian coordinates (X Y Z in meters) à 4,728,141.193, 2,879,662.605, 3,157,147.146

PPP solution for Rabat at ITRF2008 Epoch 2015.324

The published coordinates/standard deviations (95%) for the rabt1000 RINEX file are as follow:

Latitude (ITRF08): 33 59 53.1790 dms/0.003 m, Longitude (ITRF08): −6 51 15.4350 dms/0.006 m, Ellipsoidal Height (ITRF08): 90.086 m/0.013 m Cartesian coordinates (X Y Z in metres) à 5,255,617.590, −631,745.508, 3,546,322.700

The coordinates of the HELWAN and RABAT stations at ITRF2008 reference epoch (t_{0} = 2005) were published by IGS at http://itrf.ensg.ign.fr/site_info_and_select/solutions_extraction.php, and given in

To transfer the Helwan and Rabat coordinates to be in ITRF2008 at the Epoch 2014.907 and 2015.324, the approach represented by Equation (5) can be applied:

where:

The transferred coordinates of the updated Helwan and Rabat stations are tabulated in

b) Using Egyptian DM Values in the Displacement Approach

To see how the Egyptian DM can be used in the displacement approach, the following comments should be considered:

§ In the displacement approach, the geocentric cartesian (XYZ) coordinates is always used, so the need for

ITRF2008 Station Positions for Helwan at Epoch 2005.0 and Velocities | |||||||||
---|---|---|---|---|---|---|---|---|---|

30101M001 | HELWAN | PHLW | GNSS | 4,728,141.309 | 2,879,662.406 | 3,157,146.932 | 0.002 | 0.002 | 0.002 |

−0.0211 | 0.0143 | 0.016 | 0.0007 | 0.0005 | 0.0005 | ||||

ITRF2008 Station Positions for Rabat at Epoch 2005.0 and Velocities | |||||||||

35001M002 Rabat RABAT GNSS | 5,255,617.683 | −631,745.687 | 3,546,322.552 | 0.001 | 0.001 | 0.001 | |||

−.0088 | 0.0174 | 0.0142 | 0.0001 | 0.0000 | 0.0001 |

Helwan | Sol. Type | X | Y | Z |
---|---|---|---|---|

PPP-Solution | 4,728,141.193 | 2,879,662.605 | 3,157,147.146 | |

Transformed IGS-Sol. | 4,728,141.100 | 2,879,662.548 | 3,157,147.091 | |

Differences | 0.093 | 0.057 | 0.055 | |

Rabat | Sol. Type | X | Y | Z |

PPP-Solution | 5,255,617.590 | −631,745.508 | 3,546,322.700 | |

Transformed IGS-Sol. | 5,255,617.592 | −631,745.508 | 3,546,322.699 | |

Differences | −0.0026 | 0.0002 | 0.0021 |

XYZ velocities is a must, namely V_{XYZ};

§ The Egyptian deformation model, the velocity is given in topocentric (ENU), namely V_{enu};

§ To convert between the geocentric to topocentric velocities, rotation matrix, for the point’s latitude (φ) and longitude (λ) can be used:

By applying the above equation for Helwan and Rabat we get the following the results using Egyptian DM in the displacement approach.

For Helwan

For Rabat

To evaluate the computed velocity V_{xyz} of both Deformation Models, we transfer the coordinate values of Helwan and Rabat defined in ITRF2008 Epoch 2014.901 & 2015.324 to be compared with the published values of IGS defined in ITRF2008 Epoch 2005, given in

As it is depicted in

In case of the transformation between two epochs defined in two different reference frames, the following model can be used in the transformation process:

Helwan | Egypt V_{xyz} (m) | −0.018,889,752 | 0.013,798,388 | 0.015,601,729 |
---|---|---|---|---|

PPP EGY-P(XYZ)_{Epoch}_{ }_{2005} | 4,728,141.380 | 2,879,662.468 | 3,157,146.991 | |

Publ. IGS-P(XYZ)_{Epoch}_{ }_{2005} | 4,728,141.309 | 2,879,662.406 | 3,157,146.932 | |

Differences | 0.071 | 0.062 | 0.059 | |

Rabat | Egypt V_{xyz} (m) | −0.011,826 | 0.015,601 | 0.014,649 |

PPP EGY-P(XYZ)_{Epoch}_{ }_{2005} | 5,255,617.712 | −631,745.669 | 3,546,322.549 | |

Pub. IGS-P(XYZ)_{Epoch}_{ }_{2005} | 5,255,617.683 | −631,745.687 | 3,546,322.552 | |

Differences | 0.0290 | 0.018 | −0.003 |

and

where:

Rigid Plate movement is conventionally defined by a rotation rate about an Euler Pole

Equations (10)-(12) were used to compute (Ω_{X}, Ω_{Y}, Ω_{Z}) for the Nubian Plate using the Euler parameters determined by ITRF2008 PMM & Egypt’s Deformation Model, the results were depicted in

A site velocity in Cartesian format _{X}, Ω_{Y}, Ω_{Z }in radians per million years) using:

By introducing a reference epoch

Bu using Equation (13) for Helwan & Rabat and the computed angular velocity for Egypt DM & PMM, given in

Deformation Model | Absolute Pole Cartesian angular Velocity for Nubian Plate | ||
---|---|---|---|

Egypt DM | 0.000,419 | −0.002,930 | 0.003,580 |

PMM | 0.000,461 | −0.002,899 | 0.003,505 |

Egypt DM at (t_{o}) | PMM at (t_{o}) | IGS Sol. at (t_{o}) | Egypt Diff. | PMM-Diff | |
---|---|---|---|---|---|

Helwan | 4,728,141.387 | 4,728,141.384 | 4,728,141.309 | 0.078 | 0.075 |

2,879,662.450 | 2,879,662.455 | 2,879,662.406 | 0.044 | 0.049 | |

3,157,146.997 | 3,157,146.997 | 3,157,146.932 | 0.055 | 0.065 | |

Rabat | 5,255,617.674 | 5,255,617.673 | 5,255,617.683 | −0.009 | −0.010 |

−631,745.687 | −631,745.681 | −631,745.687 | 0.000 | 0.006 | |

3,546,322.544 | 3,546,322.546 | 3,546,322.552 | −0.008 | −0.006 |

For Egypt Euler Parameters:

For PMM Euler Parameters:

As it is demonstrated in

The above equation can also be used to realize a static geocentric datum aligned with ITRF at a specific reference epoch. Instantaneous ITRF positions measured at different locations and at different epochs on the same rigid plate can be related to the static datum at the reference epoch by using the same parameters (

The above equations are accurate within rigid plate zones, however within deforming zones additional parameters derived from deformation models (e.g. Finite Element Model and Fault Locking models) are required in order to maintain consistency between different epochs [

where:

PPP, global RTK and GNSS post-processing services are now used extensively to provide realizations of ITRF and WGS84 globally with a precision of a few centimeters. Unless these instantaneous realizations are transformed to a static or semi-kinematic datum by using a suitable deformation model, repeat surveys using these techniques will result in datum divergence as a function of time arising from the effects of un-modeled tectonic plate motion. Transformations from kinematic ITRF to a static datum are conventionally done by either using the site velocity (computed either from ITRF definition or from a plate motion model).

To avoid the discrepancies between the currently used datum for Egypt, namely ITRF1994 Epoch 1996, and On-Line GNSS Services, a semi-kinematic datum and related deformation models were applied and tested for Egypt. The results of some point samples show that the local Deformation Model for Egypt is giving less discrepancy than the most recent published velocities of ITRF2008 and ITRF2008 PMM. The results also confirm that the PMM model is giving good performance for the other parts out of Egypt.

MostafaRabah,AhmedShaker,MagdaFarhan, (2015) Towards a Semi-Kinematic Datum for Egypt. Positioning,06,49-60. doi: 10.4236/pos.2015.63006