The Gravity field and steady-state Ocean Circulation Explorer (GOCE) gravity satellite mission was a part of ESA’s (European Space Agency) Living Planet program  launched in the period between March 2009 and October 2013. Since begin of 2010, consecutive highly accurate GOCE-based (GOCE-only and combined) static Global Geopotential Models (GGMs) were developed using GOCE observations. The latest releases of GOCE-based GGMs have been made available to the user community in the mid of 2014, which incorporate the complete mission data set exhibiting a significant higher spatial resolution  . Up to now, GOCE-based GGMs represent the best estimate of the gravity field anomalies and mean quasi geoid heights observed from the space with a precision of about 1 mGal and 1 - 2 cm, respectively.
Three types of GOCE-based GGMs are developed through the European Space Agency (ESA) project GOCE High-level Processing Facility (HPF)   , processed mainly by three approaches: direct (DIR)  , time-wise (TIM)  and space-wise (SPW)  . For more description, the readers are referred to  . The maximum achi- eved spatial resolution of GOCE-based GGMs of types SPW, TIM and DIR ranges from 83.5 km, 71.5 km and 66.8 km, respectively based on the formula:
where is the maximum spherical harmonics (SH) degree; is the mean radius of the Earth at the reference ellipsoid (≈6378.137 km); D is the spatial resolution (i.e. half of wavelength) in km.
Moreover, other solutions of GOCE-based GGMs have been developed by different research groups such as the Gravity Observation Combination (GOCO) models  , the European Improved Gravity models of the Earth by New techniques (EIGEN) model series    , the Delft Gravity Model (DGM) solution  and the Bonn (Institute of Geodesy and Geoinformation, University of Bonn) solution (ITG-GOCE02)  .
To ensure the best accuracy of GOCE-based GGMs, especially for the local/regional geoid modeling issue, evaluation with the use of ground-based gravity data is quite important. Thus, our goal here is to examine which of GOCE-based GGMs approximate closely the gravity field over Egypt and in which spectral bands they do deliver improved information. We have to mention here that several studies concerning the validation as well as the accuracy assessment of GOCE-based GGMs have been performed in different regions of the world by different research teams, e.g.    world- wide;   in Germany;  in Central Europe:  in Poland;  in Norway;  in Brazil;   in Sudan and  in Saudi Arabia.
In the following, the gravity field functionals in terms of quasi geoid heights and free-air anomalies as determined from different of GOCE-based GGMs compared with the corresponding ones from the Earth Global Model 2008 (EGM2008,  ) at both global and regional scales over Egypt are examined. In addition, a comparison between the ground-based (terrestrial) free-air anomalies and the corresponding one obtained from GOCE-based anomalies is performed. For consistent validation in the later comparison, maximum d/o of GOCE-based GGMs will be completed using EGM2008 up to d/o 2190. In addition, the gravity signals induced beyond 2190 are compensated from local topography information which is considered as a remaining omission error.
The datasets used in this study consist of: 1) GOCE-based GGMs and EGM2008; 2) terrestrial free-air gravity anomalies and GPS/levelling data collected over the Egyptian region; and 3) High-resolution topographic data from the SRTM30_PLUS (Shuttle Radar Topography Mission) digital terrain model.
2.1. GOCE-Based GGMs and EGM2008
Over the past years, numerous GOCE-based GGMs have been developed differing in their observations time period and their processing strategies. The time span of GOCE observations range from release 1 to 5 of about 0.16, 0.66, 1.5, 2.75, 4 years, respectively. In the current investigation, eight GOCE-based GGMs based on mainly the three approaches of direct (DIR), time-wise (TIM) and space-wise (SPW) are used. These models are available via the International Centre for Global Earth Models (ICGEMs) web-service icgem.gfz-potsdam.de/ICGEM/ in terms of geopotential spherical harmonic coefficients (SH). The characteristics of those models are summarized in Table 1.
Moreover, beside GOCE-based GGMs, EGM2008 Global Geopotential Model is used in this study as state-of-the-art high resolution GGM for evaluating GOCE-based GGMs as well as for estimating higher gravity signal, i.e. the gravity signal beyond the applied maximum resolution of GOCE-based GGMs. When developing EGM2008, over Africa, the values of 5\ × 5\ mean gravity anomalies were synthesized using GGM02S spherical harmonic coefficients  for degrees 2 - 60; coefficients for degrees 61 - 360 were augmented with those of EGM96  , and for further degrees, i.e. 361 - 2159 coefficients were augmented from the analysis of the residual terrain modelling  . This model has been extensively evaluated worldwide. Its fit to the African geoid measured with standard deviation of geoid height differences was evaluated as equal to 0.73 m  .
2.2. Terrestrial Data
2.2.1. Free-Air Terrestrial Gravity Data
The available point gravity anomalies dataset in this study consists of 930 stations. Figure 1 shows the distribution of the free-air gravity anomalies for Egypt used for the current study indicated from two different sources. The date of these observations and
Table 1. The main characteristics of the investigated GOCE-based GGMs.
Figure 1. Distribution of the 784 and 146 gravity stations from older gravity network and SRI network, respectively and 17 GPS/levelling benchmarks in Egypt.
their accuracy vary greatly between the most recent Egyptian National Gravity Standardization Network of 1997 (ENGSN97), which was established by the Survey Research Institute (SRI) (146 stations), and older gravity (784 stations) surveys carried out by many private organizations several decades ago. The accuracy of ENGSN97 gravity values is 0.022 mGal, while the accuracy estimate for older gravity data is 0.5 mGal on average  .
2.2.2. GPS/Levelling Data
GPS/levelling dataset available here consists of 17 stations of the Egyptian National High Accuracy Reference Network (HARN) observed by the Egyptian Survey Authority to form the New Egyptian Datum 1995 (NED-95), distributed over Egypt as shown in Figure 1. In this network, the GPS observations were tied to the International Geodetic Stations (IGS) reference system. It should be mentioned that, the HARN network consists of 30 stations, but 13 stations (located in remote areas) have no observed orthometric heights and consequently, no undulations could be obtained for these stations. The precision of the HARN network is 0.1 parts per million (ppm), which can be written in another form as 1:10,000,000. Although the HARN network consists of 30 stations, only 17 points among them have first order levels. The precision of geoid undulations at these stations has been suggested from the provider of about 1 cm or lower. Despite the fact that the distribution of the GPS/levelling dataset of 17 stations is rather sparse, one may use them for the assessment purpose since these stations cover almost the whole Egyptian territory in different areas  .
2.3. High-Resolution Terrain Data
The SRTM30_PLUS (Shuttle Radar Topography Mission of spatial resolution about ~900 m, 30 arc-sec) data is used in this study to compute the topographic potential effect on the geoid over Egypt region. The SRTM30_PLUS has been provided by the institution of Oceanography, University of California. The detailed information concerning the development of the SRTM30_PLUS are given in the website http://topex.ucsd.edu/WWW_html/srtm30_plus.html. Figure 2 depicts the topography of Egypt from the SRTM30_PLUS.
3. Difference Amplitude of EGM2008 vs. GOCE-Based GGMs
Quasi geoid heights computed from GOCE-based GGMs and the corresponding ones obtained from EGM2008 on both global and regional scales are now compared. On the global scale, degree variances are often used to quantify the power of signal and error in the gravity field estimates at various spectral wavelengths to quantify the powers of signal and error in the gravity field estimates at various spatial wavelengths as (  , pp. 98)
with and are the differences between SH coefficients of GOCE-based
Figure 2. Topographic map showing heights in [m] as deduced from the SRTM30_PLUS data.
GGMs and the corresponding ones from the EGM2008 and are degree and order of spherical harmonics, respectively. They read as follows
Accordingly, the difference degree variances of geoid heights as shown in Figure 3 were computed here as:
With is the mean radius of the Earth at the reference ellipsoid (≈6378.137 km), the fully normalized spherical harmonics, and, are a preferred tool for many theoretical and practical applications in geodesy, especially for the representation of the Earth’s gravitational potential V 
where G is the geocentric gravitational constant; M is the mass of the Earth; are the spherical coordinates: geocentric radius, co-latitude and longitude of the computation point, respectively; is the fully normalized associated Legendre function; is the normal gravity at point; is the maximum SH degree of geopotential model applied. Therefore, most of the existing gravity field models
Figure 3. Difference degree variances in terms of geoid heights [m] between GOCE-based GGMs and EGM2008.
are formulated in terms of SH coefficients. In the geodesy representation, the related quantities such as gravity anomalies and quasi geoid heights computed with the following formulae (  , p. 272):
As seen from Figure 3, some points can be addressed here: 1) the GOCE-based direct solutions (DIR) outperform the time-wise (TIM) and space-wise (SPW) ones, particularly at the long spectrum of the gravity wavelength, where the observations of the Gravity Recovery And Climate Experiment (GRACE) mission contribute the DIR solutions. 2) Higher errors of the TIM models are provided at lower harmonics. The reason may be back to that no a-priori gravity information, neither as reference model, were considered within the processing strategies of the TIM releases. However, the TIM releases provide outperformed accuracy than the DIR solutions in such cases (cf. releases 4 and 5 of TIM and DIR) at the short wavelength of the gravity spectrum due to the contribution of the satellite gravity gradiometer (SGG) only measurements. 3) The space-wise solution, SPW_R4 model, provides high variances at lower harmonics because EGM2008 was used to model the degree variances and to calibrate the error of the estimated gravitational potential along track. 4) The Eigen-6C4 solution has the least variances especially at higher spherical harmonic degrees because of the contribution of EGM2008 altimetry/terrestrial data into Eigen-6C series (Eigen-6C3stat and -6C4).
On a regional scale, the spectral consistency between GOCE-based GGMs and EGM- 2008 should be taken in to consideration toward a reliable comparison. This is due to the fact that GOCE-based GGMs differ in resolution, besides their processing strategies with respect to EGM2008. Accordingly, the solution of each model has been restricted from SH d/o 2 up to (10 d/o step) on grids from 1.8˚ × 1.8˚, 1.64˚ × 1.64˚・・・ 0.66˚ × 0.66˚, 0.64˚ × 0.64˚, respectively. We have focused in our study only on SH degrees from 100 to 280 for two main reasons; first to examine which of GOCE-based GGMs approximate closely the gravity field over Egypt especially at the medium to short wavelength of the gravity spectral bands, where the contribution of GOCE is significant. Second, to investigate at which spectral bands they deliver improved information.
4. Comparison of GOCE-Based GGMs with the EGM2008 over Egypt
The GOCE-based GGMs specified in section 3 have been chosen as a subject for the evaluation over the study area of Egypt using EGM2008. Since the investigated GGMs differ in their data processing methodology and spatial resolution as well (see section 1), the spectral consistency between them should be taken into account. In other words, the GOCE-based GGMs and EGM2008 are examined at the same SH spectral bands. Thus, each model has been restricted from SH d/o 2 up to (10 d/o step), on grids from 1.8˚ × 1.8˚, 1.64˚ × 1.64˚・・・ 0.62˚ × 0.62˚, 0.6˚ × 0.6˚, respectively.
Figure 4 and Figure 5 show the standard deviations (std.) of differences between the free-air gravity anomalies δ∆g (Equation 9) and quasi geoid heights (Equation 10) obtained from GOCE-based GGMs and the corresponding ones from EGM2008 at the same maximum d/o of their spherical harmonics:
The results, represented in Figure 4 and Figure 5, indicate that beside GOCE data the contribution of terrestrial/altimetry data to the Eigen-6C3 and Eigen-6C4 are clearly visible from d/o 170 onward. This is because the reason mentioned in section 3 that Eigen-6C series have been developed with the use of altimetry/terrestrial data obtained from EGM2008 and DUT12 GGM. From d/o 200 onward, the std. have increased rapidly for all the GOCE-based GGMs solutions, since the coefficients beyond d/o 180 - 200 were estimated with the use of Kaula’s rule    and because above d/o 200 noise starts to dominate signals  .
Except that, the contribution of GOCE at higher harmonics is clearly visible. For instance, comparing the DIR_R4 with DIR_R5 of Figure 4 at the same spectral domain (e.g. SH d/o 210), we find that the former model estimates free-air gravity anomalies as
Figure 4. Standard deviations of the differences between gravity anomalies in [mGal] from the investigated GOCE-based GGMs () and the corresponding ones from EGM2008 truncated at the same.
Figure 5. Standard deviations of the differences between quasi geoid heights in [m] from the investigated GOCE- based GGMs () and the corresponding ones from EGM2008 truncated at the same.
well as quasi geoid heights better than the latter one. From d/o 200 onward, the std. for each of the relatively recent GOCE-based GGM releases outperform the former release of the same type. This can be obviously seen from Table 2, which shows the statistics of
Table 2. Statistics of the differences between free-air gravity anomalies in mGal (left side) and quasi geoid heights in cm (right side) obtained from GOCE-based GGMs () and the corresponding ones calculated from EGM2008.
the differences in terms of standard deviations between free-air gravity anomalies δ∆g and quasi geoid heights. The std. differences of quasi geoid heights and free-air gravity anomalies between GOCE-based GGMs and EGM2008 are at the level of about 29 - 32 cm and 6 - 7 mGal for d/o 200, respectively, of about 28 - 37 cm and 6 - 9 mGal for d/o 240 and of about 38 cm and 10.7 mGal for d/o 280.
5. Comparison of GOCE-Based GGMs with Ground-Based Data over Egypt
To provide a reliable external validation and accuracy assessment of the investigated GOCE-based GGMs, the use of ground-truth data (terrestrial gravity and/or GPS/le- velling) is quite important. In the following, the validation strategies of both free-air gravity anomalies and quasi geoid heights obtained from GOCE-based GGMs using the corresponding ones obtained from terrestrial gravity and GPS/levelling data over Egypt region are discussed.
5.1. Comparison of GOCE-Based GGMs with Terrestrial Gravity Data
Regarding the comparison between gravity anomalies as derived from the GOCE-based GGMs and those obtained from terrestrial data, the different spectral consistency is an important issue that should be taken into consideration. Spectral consistency means that we add the contribution of the short and very-short wavelength gravity signals beyond the applied maximum degree of GOCE-based GGMs. For this purpose, we add the contribution of short and very-short gravity signals using EGM2008 and the residual terrain modeling procedure as computed from the Terrain Correction TC program  , respectively.
This is due to the fact that gravity anomalies provided by the GOCE-based GGMs have a finite wavelength range of the gravity spectrum, whereas gravity anomalies obtained from terrestrial measurements contain “theoretically” the full spectral information of the gravity field, particularly the high- and very-high-frequency ranges, which are not included in the GGMs as described in Figure 6.
Numerous approaches were developed to remove the spectral inconsistency between both datasets in order to get a reliable comparison, e.g. the spectral enhancement method   , the use of Gaussian “low-pass” filtering  , and by means of orbit residuals and geoid comparisons  . Two approaches have been implemented in this study to compare between the GOCE- and ground-based gravity anomalies and geoid heights. The first method uses the spectral enhancement method (SEM) to compensate the high frequency gravity signal that is not included in GOCE-based GGMs. In the second method, a simple low-pass filtering is applied. In the following, both methods are discussed.
In  , the spectral gap between GOCE-based GGMs and terrestrial observations is bridged partially by a combination of: 1) the high-degree spectral bands of EGM2008 up to d/o 2190, and 2) the omitted gravity signal from d/o 2190 onward estimated from the topography data, i.e. from SRTM30_PLUS (Shuttle Radar Topography Mission of spatial resolution about ~900 m, 30 arc-sec), using the RTM method (see  ). Accordingly, free-air anomalies from GOCE-based GGMs, ∆gGOCE, have been computed over the 930 stations in the medium and short spectrum of gravity wavelength of d/o 100, 110, 120・・・ 300 (with 10 step). We should mention here that based on the results indicated in section 4 (Table 2), we focused on the last release of GOCE-based GGMS of type DIR_R5, TIM_R5, SPW_R4, where the contribution of GOCE observations is clearly visible in addition to the Eigen-6C4.
The remaining gravity signals in the spectral range from 101, 111, 121・・・ 301 up to 2190 were compensated from EGM2008 (∆gEGM2008). The RTM omission error estimates (i.e. from 2190 onward) as implied by the residual topography (∆gTOPO) has been determined using the TC program of  . The summation ∆gGGM reads then
The term ∆gGGM possess almost the full spectral power and can be now compared with terrestrial free-air gravity.
Figure 7 shows the standard deviations of differences (d∆gSEM) between the terrestrial
Figure 6. Spectral range (SH d/o) and spatial resolution (D in km) for satellite and terrestrial data.
Figure 7. Standard deviations of the differences between terrestrial free-air gravity anomalies [mGal] and the corresponding ones from GOCE-based GGMs extended by the EGM2008 and RTM over Egypt.
free-air gravity anomalies and the corresponding ones from GOCE-based GGMs, while the std. are indicated in Table 3. One can clearly see that the gravity anomalies as derived from the GOCE-based GGMs indicates slightly better fit to the terrestrial data than EGM2008 over Egypt. The improvement is about 0.61 mGal in terms of standard deviation of the differences. Moreover, the contribution of the SPW_R4 extended with EGM2008 and RTM based on Equation 11 is obviously seen from SH d/o 120 up to d/o 200. We can infer from Tables 1-3 that the GOCE-based GGM of type SPW_R4 are the appropriated satellite model that approximate the gravity field over Egypt region than those of type DIR_R5 and TIM_R5. Accordingly, when modeling the local quasi-geoid heights over Egypt, the SPW_R4 is strongly suggested by subtracting its wavelength spectrum up to SH d/o 200, from the terrestrial gravity anomalies. On a spatial domain, the outcome of Equations (11)-(12) is represented in Figure 8 as follow: Figure 8(a) free-air anomalies from the SPW_R4 GOCE-based GGM (∆gGOCE) up to SH d/o 200, Figure 8(b) the remaining short gravity signal from 201 up to 2190 from EGM2008 (∆gEGM2008), Figure 8(c) the omission error estimates based on RTM data (∆gTOPO), Figure 8(d) the summation of Figures 8(a)-(c) (i.e. ∆gGGM), Figure 8(e) terrestrial free air anomalies for the 930 sites (∆gTerr) and Figure 8(f) the difference between 8(d) and 8(e) (d∆gSEM).
In the second approach, a simple low pass filter (LPF) based on reciprocal distances from the computation point was applied, following  , for generating a grid of terrestrial free-air gravity anomalies in order to mimic the spatial resolution of the investigated GOCE-based GGMs at SH d/o 200 (≈100 km). The weighting function W of this low pass filter is expressed as follows.
Table 3. Statistics of the std. differences between terrestrial free-air gravity anomalies and GOCE- based GGMs at SH d/o 200 and 240.
where xi, yi are the components of the distance from the computation point in x, and y directions, respectively. In this study, the central weight Wc and the power p were both fixed to 1; the maximum distance from the computation point was fixed at a spatial scale of 100 km.
The resulting terrestrial gravity signals after applying the low pass filter are presented in Figure 9.
The terrestrial free-air gravity anomalies ∆gLPF at the location of gravity stations (Figure 1) were then extrapolated from the filtered grid as shown in Figure 9. The differences δ∆gLPF between filtered terrestrial free-air gravity anomalies (∆gLPF) and the corresponding ones (∆gGGM) computed from GOCE-based GGMs were obtained as follows
Table 4 shows the std. of the differences between free-air gravity anomalies derived from GOCE-based GGMs and the filtered ones obtained from terrestrial gravity data. The results indicate that the std. of the differences between free-air gravity anomalies from GOCE-based GGMs and filtered terrestrial gravity data are within the range of 11.96 - 12.17 mGal. It should be pointed out here that the fit of all GOCE-based GGMs and EGM08 to the terrestrial gravity anomalies is similar when using the SEM approach, whereas in case of using the LPF approach, GOCE-based GGMs exhibit improvement with respect to EGM08 by about 0.69 mGal in terms of the std. of their differences. The use of LPF approach results in about 14% reduction in std. of the differences for the investigated GOCE-based GGMs.
We can conclude here that EGM08 is not accurate enough for modelling higher frequency gravity signals (e.g. d/o 201 - 2190) over EGYPT, since it is based on fill-in observations over that area. The means of the differences slightly increased (by about 8.90 mGal) when applying low pass filter approach.
Figure 8. Free-air gravity anomalies in [mGal] of (a) the SPW_R4 GOCE-based GGM up to SH d/o 200, (b) the remained gravity signal from 201 up to 2190 from EGM2008, (c) the omission error estimates based on RTM, (d) the summation of (a) to (c), (e) terrestrial free air anomalies and (f) the difference between (d) and (e).
Figure 9. The terrestrial gravity signals after applying the low pass filter.
Table 4. Statistics of the std. of the differences between terrestrial free-air gravity anomalies after LPF approach and GOCE-based GGMs at SH d/o 200.
5.2. Comparison of GOCE-Based GGMs with GPS/Levelling Data
Here the quasi geoid heights derived from GOCE-based GGMs have been validated with the corresponding ones of ground-based GPS/levelling date on different spectral bands following the same procedures described in section 5.1. Consequently, the term was compensated by adding the term that represents the high frequency information from EGM2008 up to d/o 2190. Moreover, the term which depicts the very high frequency gravity signal beyond SH d/o 2190 was also added. It is expressed in planar approximation as (Forsberg 1984).
The complete formulae can be written as follows spectral enhancement method   :
Now, the term possess almost the full spectral domain and can be compared with the quasi geoid heights from GPS/levelling data ():
The differences () in terms of standard deviations are shown in Figure 10 and Table 5. The least differences are provided by the SPW_R4 from SH d/o 190 up to d/o ~280 (from 43.2 cm to 37.9 cm, respectively). Whereas, the std. of geoid height differences between the GPS/levelling and EGM2008 () are about 53.8 cm.
Obviously, it can be concluded that, implementing the spectral enhancement method helped much in improving the quality of GOCE-based GGMs. A considerable reduction of about 29.55% (i.e. improvement factor of about one third) in terms of standard deviations of quasi geoid heights differences is observed after compensating the medium/ short and very-short wavelength gravity signal using EGM2008 and RTM.
6. Fitting of Geoid Models from GOCE-Based GGMs to GPS/Levelling
The results indicated in section 5.2 are represented in terms of std. of the differences between GPS/levelling data and GOCE-based GGMs. However, it should be noted that the mean values of these differences are at the level of 0.30 m. This can be interpreted as a bias which may be attributed to the offset between two equipotential surfaces formed by the gravimetric geoid model derived by GOCE and Egypt vertical datum (NED-95, see section 2.2.2).
In order to minimize this bias and to determine a better fit of geoid models determined from GOCE-based GGMs and the GPS/levelling data over Egypt, we have applied the parameter transformation models of corrector surface similar to  and  . Based on this, three types of parameter transformation models have been investigated; the four-, five- and seven-parameter transformation models. The four-parameter transformation model reads  :
While the five-parameter transformation model reads  :
And the seven-parameter transformation model reads e.g.  :
Figure 10. Standard deviations of the differences between quasi geoid heights [meters] from the terrestrial GPS/ levelling and GOCE-based GGMs extended by EGM2008 and RTM over Egypt.
where are the transformation parameters between two datums; (φi, λi) are the latitude and longitude of the point; vi is the residual of random noise term and e is the ellipsoid eccentricity.
Table 6 gives the statistics of the differences between quasi geoid heights [m] as obtained from GPS/levelling data and the corresponding of GOCE-based GGMs of types DIR_ R5, TIM_R5 and SPW_R4 before and after removing the biases effect applying Equations 17-19. The results presented in Table 6 indicate that applying the 4-, 5-, and 7-parameter transformation models removes the systematic error (mean values in Table 6). The assessed accuracy of quasi geoid heights differences in terms of std. after fitting gravimetric geoid models obtained from the GGMs investigated to the corresponding ones from GPS/levelling data is at the level of about 32 cm. Table 6 indicates that the difference among the three applied transformation models in terms of their resulting statistics provide consistent results. This may indicate that practically all transformation models implemented in this study are suitable for fitting the gravimetric geoid model to the geometrical one obtained from GPS/levelling data.
This paper investigated the gravity field solutions in terms of free-air anomalies and
Table 5. Statistics of the differences between GPS/levelling and GOCE-based GGMs at SH d/o 200 and 240.
Table 6. Statistics of the differences between quasi geoid heights [m] as derived from GPS/level- ling data over Egypt and GOCE-based GGMs at SH d/o 200 extended by EGM2008 and RTM before and after applying 4-, 5-, and 7-parameter transformation models.
quasi geoid heights as derived from the ground-based data and EGM2008 over the area of Egypt as a validation strategy of the GOCE-based GGMs. To validate the GOCE- based GGMs with ground-based data, two strategies have been applied: the spectral enhancement method and low-pass filter approach.
First, gravity functionals in term of free-air gravity anomalies and quasi geoid heights derived from GOCE-based GGMs were compared with the corresponding ones calculated from EGM2008 within the spectral range of the gravitational wavelength from d/o 100 to 300. Our findings indicate some deviations of the quasi geoid heights and the free-air gravity anomalies obtained from GOCE-based GGMs and the corresponding ones from EGM2008 within the level of 29 - 32 cm and 6 - 7 mGal, respectively, at SH d/o 200, at the level of 28 - 37 cm and 6 - 9 mGal at d/o 240 and at the level of 38 cm and 10.7 mGal at d/o 280.
Second, free-air gravity anomalies obtained from GOCE-based GGMs were compared with the corresponding ones from terrestrial gravity data using the spectral enhancement method and low pass filter approach. In the former method, the free-air gravity anomalies obtained from GOCE-based GGMs agree with the corresponding terrestrial ones within the range of 14.1 mGal at SH d/o 200 providing slight improvement than EGM2008. This may back to poor coverage of gravity data in the area of Egypt and neighboring countries when developing EGM2008. In case of quasi geoid heights, a considerable reduction of about 29% in terms of standard deviations of quasi geoid heights differences is observed after implementing the spectral enhancement method. This means that the SEM strategy helps much in improving the quality of GOCE-based GGMs.
Regarding the low pass filter strategy, GOCE-based GGMs exhibit improvement with respect to EGM08 by about 0.69 mGal in terms of the std. of their differences. The use of LPF approach results in about 14% reduction in std. of the differences for the investigated GOCE-based GGMs and the terrestrial gravity data.
Finally, we have applied three types of transformation models, namely four-, five- and seven-parameters transformation model, in order to deal with the problem of datum inconsistencies between different datums of both datasets (i.e. geoid models obtained from GOCE-based GGMs and from GPS/levelling data). Our findings showed that the tree applied models exhibit all most consistent results. They suggest that the standard deviations of vertical datum over the Egyptian territory are at level of about 32 cm.
To sum up, the GOCE-based SPW_R4 (up to SH d/o 200) model is suggested to be used as a reference model when modeling the local quasi-geoid heights from the terrestrial gravity data, since it approximates the gravity field well over Egypt. This is a subject of our future consideration.
The spatial representations of the results have been plotted using the GMT5 (Generic Mapping Tools) software  . The computations of the topography potential effect in the geoid were carried out using the TC software and RTM procedure  .