The global navigation satellite system (GNSS) techniques are effective for monitoring landslides and earth movement in high-rise and coastal regions. The GNSS techniques support construction activities against earthquake effects, leading to an early warning for the catastrophic crisis and risk management. The US National Oceanic and Atmospheric Administration’s (NOAA’s) national geodetic survey (NGS) introduced the concept of continuously operating reference stations (CORS) two decades ago. This system is designed to provide users with the geodetic latitude, longitude, height, orthometric height, geopotential, acceleration of gravity, and vertical deflection at any point in the USA (Snay & Soler 2008).
The international GNSS service (IGS) provides a worldwide collection of CORS network stations. Walpersdorf et al. (2007) reported that the IGS stations in Africa have many characteristics. The network is mainly situated in a coastal area with long baselines, and this location affects the accuracy of network precision. In addition, the network is very active with ionosphere activities.
Egypt is rapidly undertaking many infrastructure projects such as highways, railway projects, and the construction of many new cities. One of the main obstacles in obtaining reliable continuous GNSS solutions is the lack of coverage for IGS stations. Consequently, the Egyptian Surveying Authority (ESA) established the first permanent Egyptian CORS network in January 2012. The network consisted of 40 stations covering mainly the Nile valley and its delta (Figure 1). Furthermore, according to the ESA (2012), the network was adjusted relative to the International Terrestrial Reference Frame (ITRF) system at the epoch October 23, 2011 (ITRF 2020). Moreover, the ESA older national passive GNSS reference network, called the high-accuracy reference network (HARN) (Figure 1), established in 1995, was referenced to the ITRF1994 (epoch 1996).
Zumberge et al. (1997) introduced the precise point positioning (PPP) estimation model two decades ago. The PPP technique that requires just one GNSS receiver faces many error sources, and these errors need to be modeled or eliminated. Four types of errors are considered: satellite-dependent, receiver-dependent, atmospheric (troposphere and ionosphere), and geophysical errors, including solid earth tides, polar tides, ocean tidal loading, and atmospheric tidal loading. Several software packages and online services provide the PPP solution (Table 1).
List of available PPP software and free online services
Software package | Provider | |
---|---|---|
PPP software package | gLAB | Research group of Astronomy and Geomatics (gAGE) from the Universitat Politecnica de Catalunya (UPC) (Sanz et al. 2012) |
Net_Diff | GNSS Analysis Center, Shanghai Astronomical Observatory, Chinese Academy of Sciences ( |
|
PRIDE PPP-AR | GNSS Research Center of Wuhan University, China ( |
|
GAMP | Research group at the German Research Centre for Geosciences GFZ (Zhou et al. 2018) | |
Bernese GNSS V. 5.2 | Astronomical Institute of the University of Bern (AIUB), Switzerland (Dach & Walser 2015) | |
GIPSY-OASIS II | Jet Propulsion Laboratory (JPL), California Institute of Technology, USA ( |
|
RTKLIB | Tomoja Takasu at Tokyo University Marine Science and Technology, Japan ( |
|
NAPEOS | European Space Agency ( |
|
GAMIT/GLOBK | Department of Earth, Atmospheric, and Planetary Sciences at the Massachusetts Institute of Technology (MIT) ( |
|
PPP free online service | CSRS-PPP | Canadian Geodetic Survey (CGS) ( |
APPS-PPP | Jet Propulsion Laboratory (JPL), California Institute of Technology, USA ( |
|
GAPS-PPP | GNSS Analysis and Positioning Software, Department of Geodesy and Geomatics Engineering University of New Brunswick ( |
|
magicGNSS | GMV Aerospace and Defence S.A., Spain ( |
Researchers have recently shown increasing interest in evaluating the PPP accuracy of CORS network in Africa. Abdallah (2015) investigated different static PPP solutions for four CORS in Germany and found that solutions reached centimeter-level accuracy using the static technique for a convergence time of 24 h. This result matched the one obtained by Ayhan and Almuslmani (2021). In addition, Abdallah and Schwieger (2016) introduced a study of the PPP solution for the IGS-CORS in Africa. They concluded that the PPP solution obtained using (Canadian Spatial Reference System) CSRS-PPP online service produced three-dimensional (3D) root mean square error (RMSE) between 4 and 6 mm. Similar results based on GNSS datasets in Egypt were reported by El Shouny and Miky (2019) and El Manaily et al. (2017). Rabah et al. (2016) proposed the utilization of PPP for geodetic datum maintenance and update in Egypt. Jamieson and Gillins (2018) determined the accuracy of static PPP solution for a static campaign, which consisted of six stations with 10 h of observation data. The study showed less than 1 and 1.5 cm accuracies for the horizontal and vertical directions, respectively. To evaluate the accuracy of the new Nigerian GNSS network (NIGNET), Isioye et al. (2019) studied different PPP online services with different observation durations (1, 2, and 6 h, and up to 24 h). The study used free online PPP services (e.g., GNSS Analysis and Positioning Software (GAPS-PPP) and magicGNSS). They found that the online GNSS PPP processing services can provide users with reliable results. Further, the results showed that the 24-h observation files produced results with a millimeter to centimeter level of accuracy by processing with magicGNSS and GAPS-PPP services.
This study aims to contribute to this growing area of research by exploring the accuracy of static PPP solutions for the Egyptian permeant CORS (ESA-CORS). The observation data of ESA-CORS are accessible only for the ESA and unavailable for users. Therefore, this paper assesses the significance and reliability of using CSRS-PPP as a solution to obtain a geodetic coordinate, such as a static PPP solution. The overall structure of the study comprises five sections, including this introductory section. Section 2 begins by laying out the theoretical approach of PPP solution and explaining the processing parameters using CSRS-PPP online service. Section 3 discusses the dataset and evaluation methodology. Section 4 presents the analysis and discusses the results obtained. Finally, the conclusion provides a summary of our findings.
As shown in equations (1)–(4), the PPP position is determined based on the processing of the following ionosphere-free combination of the undifferenced code Δiono is the ionospheric delay, Δtrop is the tropospheric delay, Δsol is the solid earth tides error, Δpol is the pole tides error, Δocn is the ocean loading effect, Δatm is the atmospheric loading effect, Δmul is the multipath effect, Δpcv is the antenna phase center variation, and ∈
The tropsphereic zenith delay (Δtrop) consi sts of dry and wet components. The dry component (zenith hydrostatic delay) is about 90% of the delay; this part can be modeled using mathematical models. The wet component (zenith wet delay) is unpredictable; it depends on water vapor and is typically less than 30 cm (Hofmann-Wellenhof, Lichtenegger & Wasle 2008; Bar-Sever, Kroger & Borjesson 1998). The total tropospheric zenith delay is expressed in equation (5). Due to the consideration of the arbitrary zenith angle of the signal, the estimation of the delay is expressed in the function of the elevation angle of the satellite
Table 2 presents the processing parameters of CSRS-PPP software. The reference system for the software is based on ITRF 2014 (epoch 2019 according to measurement date); the obtained coordinates have Cartesians XYZ format and Ellipsoidal/(Universal Transverse Mercator) UTM system. The IGS final ephemerides are used during processing with a satellite orbit of a 15-min interval and a satellite clock of a 30-s interval. The ionospheric delay is eliminated for dual-frequency data using the ionospheric-free linear combination, and the second-order errors are equally considered. For single-frequency data, the ionospheric delay is modeled using the final global ionospheric maps (GIM) from IGS. Regarding the tropospheric delay, the processing tool that is based on the global pressure and temperature (GPT) model uses the Davis model for the hydrostatic delay. For the wet part, the processing mechanism uses the Hopfield model that is based on the GPT model. In addition, the Vienna mapping function (VMF1) is utilized (Boehm, Werl & Schuh 2006). The satellite and receiver antenna phase variations are based on the IGS-ANTEX (NGS 2021) format (CSRS-PPP 2021). Furthermore, the PPP-based coordinates in ITRF 2014 (epoch 2019) were compared with the recent network adjustment in the same frame and epoch.
CSRS-PPP processing parameters
Item | Processing parameter |
---|---|
GNSS system | GPS and GLONASS |
Reference system | ITRF2014 (epoch 2019) |
Coordinate format | XYZ/Ellipsoidal/UTM |
Satellite orbit and clock ephemeris | IGS final (Orbit: 15-minute interval & clock: 30-second interval) |
Satellite phase center offsets | IGS-ANTEX (NGS 2021) |
Receiver phase center offsets | IGS-ANTEX (NGS 2021) |
Elevation angle | 7.5° |
Sampling rate | 30 s |
Ionospheric model | Linear ionosphere-free combination + second-order parameters |
Tropospheric model | Dry: Davis (GPT) (Sakurai 1985; Kouba 2009) Wet: Hopfield model (GPT) (Hopfield 1969; Kouba 2009) VMF1 (Böhm & Schuh 2004) |
Software | CSRS-PPP |
Observation data | Single/dual frequency and static/kinematic |
Data transfer | |
Ocean tide loading | FES 2004 (Le Provost & Lyard 1997) |
Overall, 32 ESA-CORS were analyzed for three consecutive days (DOY: 201–203/2019) with an observation interval of 30 s (Figure 2a) using CSRS-PPP online service (see Table 2). For global impact, six IGS-CORS were also processed for the same observation days (Figure 2b); these stations are surrounding Egypt. The complete list of processed stations is presented in Table 3. The height of ESA-CORS varied between 30 and 149 m. Moreover, the height of IGS-CORS (shadowed stations in Table 3) varied between 31 and 886 m. Figure 2a shows the ESA-CORS and Figure 2b shows the IGS-CORS.
List of CORS
Station ID | Latitude | Longitude | H | Station ID | Latitude | Longitude | H |
---|---|---|---|---|---|---|---|
ADFO | N24°58′46.241″ | E32°52′49.913″ | 105.97 | QANA | N26°09′22.293″ | E32°43′18.173″ | 101.20 |
ADWH | N28°41′41.549″ | E30°46′01.677″ | 64.61 | QANT | N30°51′27.105″ | E32°18′36.257″ | 34.11 |
ASHM | N30°17′40.001″ | E30°59′04.720″ | 39.67 | RMDN | N30°17′26.508″ | E31°44′23.437″ | 149.50 |
AYAT | N29°37′06.841″ | E31°15′22.315″ | 58.82 | RSHD | N31°24′19.102″ | E30°25′17.957″ | 36.85 |
BADR | N30°34′57.049″ | E30°42′37.032″ | 47.94 | SDAT | N30°22′31.934″ | E30°29′34.676″ | 71.50 |
BLTM | N31°33′11.924″ | E31°05′40.480″ | 30.37 | SHKH | N31°06′24.609″ | E30°56′26.923″ | 30.55 |
BNHA | N30°28′15.304″ | E31°10′45.118″ | 46.02 | SMLT | N28°18′38.917″ | E30°42′39.376″ | 71.19 |
CARO | N30°01′55.725″ | E31°12′55.058″ | 75.57 | SUEF | N29°04′53.755″ | E31°05′30.069″ | 60.47 |
DMNH | N31°02′27.912″ | E30°27′57.358″ | 41.63 | SUZE | N29°59′18.836″ | E32°31′19.498″ | 47.13 |
ETSA | N29°14′15.413″ | E30°47′49.999″ | 41.64 | TANT | N30°47′36.413″ | E30°59′35.224″ | 40.81 |
HMOL | N31°18′35.684″ | E31°09′02.195″ | 44.08 | THAT | N26°46′09.649″ | E31°29′50.601″ | 87.69 |
ISML | N30°36′15.307″ | E32°18′00.940″ | 43.27 | TSHT | N26°07′06.090″ | E32°05′40.335″ | 88.37 |
ISNA | N25°17′20.274″ | E32°33′28.518″ | 102.49 | WAKF | N26°03′41.760″ | E32°25′52.706″ | 109.83 |
KBER | N30°43′56.659″ | E31°40′21.508″ | 33.15 | DYNG | N38°04′42.788″ | E23°55′56.765″ | 510.56 |
LXOR | N25°42′29.159″ | E32°39′09.500″ | 110.58 | NICO | N35°08′27.558″ | E33°23′47.216″ | 190.02 |
MNSH | N26°29′00.621″ | E31°47′39.230″ | 96.24 | NKLG | N00°21′14.072″ | E09°40′19.660″ | 31.49 |
MNZL | N31°09′36.216″ | E31°56′04.274″ | 36.95 | NOT1 | N36°52′33.049″ | E14°59′23.246″ | 126.35 |
MOUS | N27°38′36.827″ | E30°51′27.876″ | 82.44 | RAMO | N30°35′51.388″ | E34°45′47.314″ | 886.84 |
MTMR | N30°54′45.358″ | E30°09′31.184″ | 33.16 | YKRO | N06°52′14.023″ | W05°14′24.327″ | 270.26 |
To obtain the reference solution of ESA-CORS, Trimble Business Center (TBC) V 5.0 was used to process our geodetic network. The final precise satellite orbits “.SP3” and the satellite clock corrections “.clk” were downloaded from the IGS website and included in data processing. To achieve such a high-precision GNSS solution, ionospheric and tropospheric delays as well as ocean tide loading were considered. The processing strategy consisted of two steps: firstly, six ESA-CORS (ADFO, ALEX, CARO, MOUS, QANT, and SUZE) were first processed with six tied IGS-CORS (YKRO, DYNG, NICO, NKLG, NOT1, and RAMO). The reference solution for IGS-stations was obtained from
As shown in equation (7), the evaluation methodology is presented using estimation of the error values δ
To assess the results, the consistency of results obtained from each case study should be evaluated separately.
In Figure 5, the accuracy obtained in east, north, and height directions for DOY 201–202 and 203/2019 is presented. Furthermore, Table 4 shows the estimated statistics for all stations for the three consecutive days and the mean values. This table also shows the calculated values of maximum, mean, minimum error values, RMSE, and SD. SD95% is estimated by excluding the outliers.
Statistics of static-PPP solution for ESA-CORS in millimeters
DOY201 | DOY202 | DOY203 | Overall mean values | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
E | N | h | E | N | h | E | N | h | E | N | h | |
Max. | 20.0 | 15.0 | 44.0 | 10.0 | 15.0 | 48.0 | 9.0 | 10.0 | 48.0 | 13.0 | 13.3 | 46.7 |
6.5 | 4.2 | 18.5 | 6.0 | 4.6 | 17.4 | 4.3 | 4.2 | 17.6 | 5.6 | 4.3 | 17.8 | |
Min. | 0.0 | 1.0 | 0.0 | 1.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.3 | 0.3 | 0.0 |
RMSE | 7.7 | 5.6 | 22.3 | 6.4 | 5.8 | 22.7 | 4.8 | 5.0 | 22.2 | 6.3 | 5.5 | 22.4 |
SD | 4.2 | 3.8 | 12.8 | 2.3 | 3.6 | 14.8 | 2.1 | 3.8 | 13.7 | 2.9 | 3.7 | 13.8 |
SD95% | 2.1 | 1.8 | 5.8 | 2.0 | 1.8 | 6.3 | 1.9 | 1.8 | 7.0 | 2.0 | 1.8 | 6.4 |
For DOY 201, the horizontal components showed an error of up to 20 mm with an average value of 4–6 mm using an RMSE of 6–7 mm (SD = 4 mm). Only ADWH, BADR, QANA, ADFO, and WAKF stations showed an error between 1 and 1.5 cm. It is clear from the result obtained that the error presented a value of up to 44 mm (
The error for DOY203 was up to 10 mm (
Three stations continuously lacked observation data for 1 day or three successive days (Table 5). In addition to the tabulated stations, other stations also had some losses in the observation data for various epochs. These continuous lacks or small losses affected the quality of PPP solutions. ADFO station lacked data for 4 h at the beginning of DOY 201, and the data also had some other losses (see Figure 6). However, some high errors were reported in the north direction. ADWH station lacked observation data for 1 h on DOY 201, which yielded a high error in the east direction. Furthermore, ASHM and AYAT stations had continuous observation data that reflected the high error in height direction for DOY 201. BADR station had an extreme lack of observation data for DOY 201, which led to high errors. For stations BLTM, DMNH, KBER, QANA, and TSHT, the continuous lack of observation data decreased PPP accuracy in the height direction. Several issues were identified for ISNA, LXOR, MNSH, MOUS, MTMR, RSHD, SMLT, and THAT stations; the solution reported a high error for height direction due to several data losses during the convergence time for the 3 days. Finally, WAKF station had a high error for the height direction for DOY 202/203. The possible reason for the losses is bad Global System for Mobile Communications (GSM) transfer for the observation data to the center in Cairo, especially for the stations in Upper Egypt (e.g., WAKF, ISNA, LXOR, and MNSH).
Statistics of continuous lack of observation data for ESA-CORS
Station ID | DOY 201 | DOY 202 | DOY 203 | Station ID | DOY 201 | DOY 202 | DOY 203 |
---|---|---|---|---|---|---|---|
ADFO | 4 h | - | - | QANA | 9 h | 9 h | 9 h |
ADWH | 1 h | - | - | QANT | 1 h | - | - |
BLTM | 6 h | 12 h | 8 h | RMDN | 3 h | - | - |
DMNH | 13 h | 6 h | 5 h | THAT | 1 h | - | - |
ISNA | 2 h | - | - | TSHT | 8 h | - | - |
KBER | 16 h | 16 h | 16 h | WAKF | 22 h | - | - |
MNSH | 1 h | - | - |
Figure 7 presents the PPP errors in the east, north, and height directions; in addition, Table 6 shows the detailed statistics of the PPP solution (max., min., average, RMS, SD, and SD95%). For DOY201, the horizontal error showed up to 3 mm (
Statistics of static-PPP solution for IGS-CORS in millimeters
DOY201 | DOY202 | DOY203 | Mean | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
E | N | h | E | N | h | E | N | h | E | N | h | |
Max. | 3.0 | 3.0 | 8.0 | 2.0 | 2.0 | 12.0 | 4.0 | 2.0 | 3.0 | 2.3 | 8.3 | |
1.3 | 1.2 | 3.3 | 1.0 | 1.3 | 5.5 | 1.2 | 1.2 | 3.0 | 1.2 | 1.2 | 3.9 | |
Min. | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 1.0 | 0.0 | 1.0 | 2.0 | 0.0 | 0.3 | 1.0 |
RMSE | 1.7 | 1.6 | 4.5 | 1.3 | 1.5 | 6.7 | 1.8 | 1.2 | 3.2 | 1.6 | 1.4 | 4.8 |
SD | 1.2 | 1.2 | 3.3 | 0.9 | 0.8 | 4.2 | 1.5 | 1.2 | 1.3 | 1.2 | 1.1 | 2.9 |
SD95% | 1.0 | 0.8 | 1.0 | 0.9 | 0.8 | 1.0 | 0.5 | 0.8 | 1.3 | 0.8 | 0.8 | 1.1 |
Finally, from Table 6, the error obtained in DOY203 for the east direction was up to 4 mm (μ = 1.2 mm) with an RMSE of 1.8 mm (SD = 1.5 mm). For north direction, the obtained error was up to 2 mm (
Establishing GNSS CORS in Egypt is considered one of the main means to provide precise positioning and navigation for all development applications. This study provided a framework for investigating the accuracy of 32 ESA-CORS in addition to six IGS-CORS. Three consecutive observation days were analyzed using CSRS-PPP’s free PPP online service.
The overall PPP solution for ESA-CORS showed up to 13 mm errors with a mean value
Regarding the PPP solution for six IGS-CORS, this study has shown that the errors for East and North were up to 3 mm with a mean value of 1.2 mm (RMSE was circa 1.5 mm). In the height direction, the PPP solution was up to 8 mm (