Hybrid Geoid Model : Theory and Application in Brazil

Determination of the ellipsoidal height by Global Navigation Satellite Systems (GNSS) is becoming better known and used for purposes of leveling with the aid of geoid models. However, the disadvantage of this method is the quality of the geoid models, which degrade heights and limit the application of the method. In order to provide better quality in transforming height using GNSS leveling, this research aims to develop a hybridization methodology of gravimetric geoid models EGM08, MAPGEO2015 and GEOIDSP2014 for the State of São Paulo, providing more consistent models with GNSS technology. Radial Basis Function (RBF) neural networks were used to obtain the corrector surface, based on differences between geoid model undulations and the undulations obtained by GNSS tracking in benchmarks. The experiments showed that the most suitable interpolation for correction modeling is the linear RBF. Checkpoints indicate that the geoid hybrid models feature root mean square deviation ± 0.107, ± 0.104 and ± 0.098 m, respectively. The results shows an improvement of 30 to 40% in consistencies compared with the gravimetric geoids, providing users with better quality in transformation of geometric to orthometric heights.


INTRODUCTION
The altimetric component is fundamental for a variety of scientific and engineering applications, such as topographic and geodesic surveying, coastal studies, engineering projects, dams, irrigation and other things.However, the GNSS (Global Navigation Satellite Systems) reference surface (ellipsoid) has no connection with the Earth's gravity field, thus equal values of geometric height can occur in different levels or uneven height values may occur on the same equipotential surface of the gravity field (Freitas and Blitzkow 1999).For this reason, many altimetric applications cannot directly use the height component from GNSS.In the first place the geometric height must be reduced to normal, orthometric or dynamic height, and usually quasi-geoid or geoid models are used to aid this conversion ( Featherstone 2000, Kotsakis and Sideris 2001, Torge 2001, Seeber 2003).
Some users such as oceanographers and geophysicists require decimeter quality for geoid model (Kuroishi 2009), but certain geodetic and cartographic applications require a quality close to 1 cm, such as global unification of vertical datum and deglaciation studies (Plag and Pearlman 2009).
The geometric heights obtained by GNSS/leveling in benchmarks allow users to have the distance between the reference ellipsoid and the geoid surface, i.e., the geoid undulation or geoid height.This value is considered closer to the real static geoid than the gravimetric geoid models, for this reason benchmarks are used for evaluation of gravimetric geoid models.However, we must bear in mind that even the geoid from benchmarks is not the "true" geoid.Since different types of measurements and techniques lead to the benchmark geoid, systematic and aleatory errors are present, such as leveling errors, GNSS positioning errors and permanent tidal systems errors.
Another application of GNSS/leveling in benchmarks is the refinement of the gravimetric geoid which has been done by many countries around the world, here we call hybridization process.The integration of GNSS/leveling and gravimetric geoid models gives rise to hybrid geoid models.The principle method is fitting a corrector surface obtained by interpolation in a chosen gravimetric geoid.The generation of this surface is discussed in depth in subsequent sections.Studies have shown that these models are compatible with geodetic reference frames and have provided better quality results in the transformation of geometric to orthometric heights than gravimetric geoid models.
The improvement potential of hybrid geoid models is to the centimeter level, e.g., the case of Japan, with an root mean square difference (RMSD) of ± 0.02 m (Miyahara et al. 2014); cities of China (Hong Kong and Shenzhen) showed an RMS of ± 0.01 m (Chen and Yang 2001); the United States (CONUScontiguous United States) presented an RMS of ± 0.01 m (Roman et al. 2010); Turkey achieved an RMS of ± 0.03 m (Erol and Erol 2013); and South Africa observed an RMS of ± 0.07 m (Chandler and Merry 2010).Other countries have also used the hybridization, in the case of Canada (NRCan 2014), South Korea (Lee et al. 2012, 2013) and Australia (Fotopoulos et al. 2002, Featherstone et al. 2005).
This paper describes a case study showing the methodology applied to obtain a hybrid geoid model in a region very important economically in Brazil, the State of São Paulo.The corrector surfaces were generated by neural networks RBF and the Kriging method.

GRAVIMETRIC GEOID MODELS
Among the several geoid models already published, the EGM (Earth Gravitational Model) series deserves mention.The EGM96 was the first model of the series.This model was conducted by the NGA (National Geo-spatial Intelligence Agency) in partnership with the GSFC (Goddard Space Flight Center) in mid of 90s.The latest model in the series, the EGM08, reached degrees 2190 and order 2159 of spherical harmonic function coefficients, with a resolution of 5' and accuracy around ± 0.15 m (Pavlis et al. 2012).A global analysis by Pavlis et al. (2008) of EGM08 undulations compared to the undulations obtained by GNSS tracking in 12,305 benchmarks discloses an RMS ± 0.103 m.Discrepant results were obtained by Barthelmes and Köhler (2012), using 12,036 points in comparison.Then, the RMS obtained was ± 0.239 m.For Brazil, the statistics obtained by Matos et al. (2012) also diverge from ICGEM (2015).The first, using 844 points, obtained an RMS ± 0.330 m, while the ICGEM, with 1,112 points, showed an RMS ± 0.460 m.
In the case of Brazil, MAPGEO2010 (Mapa Geoidal -2010) was the official geoid model for the country until 2015, when the newest model came out, MAPGEO2015.This series of national-scale models was developed in partnership between IBGE (Instituto Brasileiro de Geografia e Estatística) and EPUSP (Escola Politécnica da Universidade de São Paulo).The methodology used to create these models was Fast Fourrier Transform (FFT) and the modified Stokes integral, with use of the Removed-Restore technique, where the EIGEN-6C4 (Förste et al. 2014) to degree and order 200 was used to obtain the medium and long wavelengths.The last geoid model has an improvement of 20% in comparison with its antecessor in consistence with GNSS positioning, providing to users an RMS of ± 0.17 m, comparing 592 points tracked by GNSS (IBGE 2015).
Since the area of Brazil is 8,515,767 km 2 , the common practice which in a lot of States has been to develop its own local geoid model for scientific and technical purposes.The State of São Paulo, as an important economic region, has many different geoid models.Since 2008 the GEOIDSP (Geoide de São Paulo) series have been elaborated, with the same methodology of the MAPGEO series.The last local geoid model developed for the State of São Paulo was the GEOIDSP2014 (Guimarães et al. 2015), which used the EIGEN-6C3STAT (Förste et al. 2013) up to degree and order 200 for obtaining medium and long wavelengths.According to Guimarães et al. (2015), the validation of GEOIDSP2014 through 168 points tracked by GNSS showed an RMS ± 0.21 m.
The schedule to determine the reginal geoid model using FFT can be described in 5 steps (Blitzkow et al. 2009): 1. Calculation of point free air gravity anomalies through terrestrial gravimetric data (coordinates, orthometric height and gravity acceleration); 2. Calculation of complete Bouguer anomalies in order to derive mean free air gravity anomalies.The 5' x 5' grid of these anomalies was computed from point gravity data.Over the ocean, DTU10 was used; 3. Calculation of Helmert gravity anomalies referred to the surface of the Earth.These are obtained from the mean free air anomaly by adding Direct Topographical Effect (DTE), Direct Atmospheric Effect (DAE) and Secondary Indirect Topographical Effect (SITE) (Ellmann and Vaníček 2007); 4. Modified Stokes' kernel integration with the use of the spectral decomposition to calculate the cogeoid; 5. Primary Indirect Topographical Effect (PITE) was added to co-geoid heights to obtain geoid heights (Martinec andVaníček 1994, Martinec 1998).

GNSS POSITIONING AND GEOID MODEL INTEGRATION
The integration of satellite positioning technologies with the various existing geoid models has different purposes, one of which gives rise to GNSS leveling technique and other results in the geoid hybrid models.
For the development of these techniques, knowledge of the methods and concepts of GNSS positioning is fundamental to obtain the ellipsoidal height and geoid models to calculate the geoid undulations.
Regarding the gravimetric geoid models, many computation techniques can be used, such as numerical integration of Stokes formula directly or using Fast Fourier Transform (FFT), Least Square Collocation (LSC), spherical harmonic functions developed in series, direct calculation by the difference between the ellipsoidal height and orthometric height (from spirit leveling) and radar altimetry (in the oceans) (Featherstone andSideris 1998, Sideris 2013).But a combination or integration between techniques can provide better results in terms of compatibility.
Although there are many techniques and methods to obtain gravimetric geoid models, these are commonly validated by GNSS survey in benchmarks (Seeber 2003).For this reason the integration between these geoid models and GNSS leveling data can make it more consistent with satellite tracking.
Three surfaces are considered in geometric to orthometric transformation height: the physical surface of the Earth, where geodesic operations are carried out; the geoid surface, adopted as altimetric reference; and the ellipsoidal surface, limiting the geometric model adopted for the Earth, where geodesic calculations are performed.This assumption can be seen in Figure 1.A simplified mathematical model of the relationship between geometric (h) and orthometric (H) height is considered from the coincidence between the vertical (v) and normal (n) lines.
According to Jekeli (2000) and Sjöberg ( 2006), the consideration of coincidence between the lines can lead to an error of 0.00001% to 0.00002%.This approach limits error to no more 1.6 mm in an extreme case where the elevation reaches 8,000 m.
Following these assumptions, the methodology summarized to obtain hybrid geoid models is shown in Figure 2. The input data for hybrid geoid generation are GNSS/leveling geoid heights (N GN SS ) and a gravimetric geoid grid (N model ).The hybridization is processed in three steps.First, the calculation of offsets for each GNSS/leveling point; second, the generation of corrector surface; and third, fitting the obtained surface to the gravimetric geoid grid.
Initially the hybrid models depend on offsets obtained between GNSS data from benchmarks (N GN SS ) and a gravimetric geoid model (N model ).These offsets, calculated according equation 2, are regarded as corrections (∆N ) to the gravimetric geoid model.However, before treating the offset as corrections it is necessary to transform them in an equal grid resolution of gravimetric geoid model, i.e., a corrector surface.Figure 3   Once the offset interpolation is finished, the corrector surface generation is complete and can be added in the gravimetric geoid model grid (equation 3).The hybrid geoid model is obtained from the sum of the grids.
The refinement of the conventional geoid is performed by the corrector surface (∆N ).For this reason, adequate surface representation is crucial for a good quality of the hybrid model.There are several interpolation methods found in the literature for modeling geoid model offsets, such as the multiquadric (Hardy 1975), the weighted average (Chen and Luo 2004), the bilinear and bicubic (Featherstone 2000), the poly-nomial surfaces (Fotopoulos et al. 2002); Shepard (Chen and Luo 2004); and Kriging (Chandler and Merry 2010).However, the most widely used method of interpolation is LSC (Fukuda et al. 1997, Featherstone 2000, Fotopoulos et al. 2002, Featherstone et al. 2005, Kuroishi et al. 2002, Chen and Luo 2004, Roman et al. 2004, 2010, Miyahara et al. 2014).An interesting point of this interpolating is that it considers the spatial distribution of observational errors, thus the method requires a covariance matrix of data.Some authors have been using Artificial Neural Networks (ANN) to obtain geoid undulations (Maia and Segatine 2003, Kavzoglu and Saka 2005, Veronez et al. 2011).With the results, this type of interpolator has been found to have great potential in modeling offsets.Moreover, Fazio and Roisenberg (2013) prove mathematically that these interpolators deliver results equivalent to the Kriging method and can produce identical results if properly configured.In particular, RBF (Radial Basis Function) neural networks are much faster than the Kriging method because the RBF interpolation weights need not be recomputed for each interpolated point.
However, the selection of the interpolation method to represent the corrector surface must be performed empirically as specific interpolation can better represent the spatial reality of the data and the study site (Chen and Luo 2004).

HYBRID GEOID MODELS
Since nowadays GNSS positioning is widely used for leveling purposes, users must perform height transformation.For this, it is essential to focus on more accurate geoid models.For this research, methods and techniques on the geoid calculation have been improved.The most accurate results for geoid models are centimeter or decimeter-level, according to region extension, methodology, geopotential model, gravimetric data and GNSS/leveling data for validation.
More gravimetric data can make the geoid model more realistic, but in some cases it does not show any improvement in comparison with GNSS/leveling points.This happens not only because of uncertainties of gravimetric and GNSS/ leveling data but also because the two techniques are not completely compatible.The cause of this imperfect compatibly is source of geoid undulation quantities, while geoid models uses gravimetric data with only-physical meaning the GNSS/leveling in benchmarks provides geoid height with parcial-physical or geometric-physical effect.In this context, the hybridization process of traditional gravimetric geoid takes place.Some authors may use a different nomenclature for geoids resulting from this methodology.Generally, improvement in local contexts or "refinement" of the geoid generates no surface correction separately and simpler surfaces are used to describe the hybrid model.Featherstone et al. (1998) refer to these models as "combined geoids", since an adjustment by least squares of the data source is proposed with gravimetric GNSS leveling data.This type of combination in areas between 1,000 and 3,000 km 2 have generated RMS between ± 0.01 and ± 0.05 m.
The hybridization methodology shows interesting points since hybrid models treats offsets between geoid models and GNSS/leveling in benchmarks as corrections.Therefore, this method can be seen as a way of making geoid models more consistent with GNSS positioning.The residual inconsistency present in hybrid geoid models are due to uncertainty of GNSS positioning, datum, leveling networks, gravimetry and geopotential models.
Several countries have adopted this hybrid geoids proposal getting more consistent results with the positioning methods.The importance and the global concern is evident for the users of GNSS in converting heights with better quality.Among the countries that have adopted this methodology are: The United States -CONUS - (Roman et al. 2004, 2010, NGS 2015), Canada (NRCan 2014), Australia (Fotopoulos et al. 2002, Featherstone and Sideris 1998, Featherstone et al. 2005), South Korea (Lee et al. 2012, 2013), Japan (Kuroishi et al. 2002, Miyahara et al. 2014); China (Chen and Yang 2001, Chen and Luo 2004), Turkey (Erol et al. 2008, Erol andErol 2013) and South Africa (Chandler and Merry 2010).
Brazil still does not have a hybrid geoid available on a regional scale for civil uses.Goldani (2006) has presented a geoid model developed only with benchmark data tracked by GNSS for the State of Paraná, the GEOIDERAAP, where RMS obtained was ± 0.10 m, evidencing the model improvement over the national geoid when using this type of data in the model.The State of São Paulo has been the subject of several studies to improve the geoid through alternative methods with integration of GNSS leveling (Maia and Segatine 2003, Kavzoglu and Saka 2005, Veronez et al. 2011).
The analysis of discrepancies between undulations of benchmarks, tracked by GNSS, and the hybrid geoid model at control points indicate an interpolating noise related with its own mathematical characteristics, while the quality of the hybrid geoid model is shown by discrepancies at checkpoints considering all errors involved.
Each country can adopt a specific methodology for generation of corrections and validation.For example, the geoids in South Africa and China were evaluated by the conventional method of quality analysis.Other countries, such as Australia and the USA, do not use the conventional method of validation because the LSC methodology provides realistic estimates of the hybrid model precision.
In the case of the Japanese models, an empirical method validation is done considering the independence of other data and the values obtained confirm that the SD is compatible with the estimated model RMS.More details can be seen in Kuroishi et al. (2002) and Miyahara et al. (2014).

ARTIFICIAL NEURAL NETWORK RBF
The Artificial Neural Networks have been chosen as one of the focuses in this research study because of the facility of application, ability to learn from examples, robustness and speed of processing (Tubb 1993).This technique presents a mathematical model inspired by the neural structure of intelligent organisms, acquiring knowledge through experiences such as pattern recognition, interpolations and ratings (Wang 2003).
The RBF neural networks have a similar architecture to the MLP (multilayer perceptron) networks model.However RBF networks use distances between points for data interpolation and have only one hidden layer.The functions use distances as input, from this variable a number of neurons are activated (Wright 2003).A very interesting feature of these networks is their speed and efficiency.
The RBF networks use radial basis functions as neurons, distributed spatially.The functions use distances as input, from this variable a number of neurons are activated (Wright 2003, Haykin 2009).Some examples of radial basis functions are shown below: where r is a distance; and µ is an additional smoothing parameter that needs to be estimated for each function according to the data.
The distance (φ (i,j) ) between point (i) and the center of the radial basis function of each neuron (j) during the interpolation process of the RBF networks can activate more or fewer neurons, resulting in a higher or lower contribution from neurons to the final result.
The parameters of radial basis function, the position of the center of each neuron of radial basis function, the number of neurons used and the equation of weights are defined by training the RBF Networks.During this process the RBF Networks establish a single function to the data set from the center of the radial basis function.The definition of the center of the radial basis function can be made mainly by three methods.The exact RBF kernel was used in this work.In this method, center of the radial basis function is placed on the sample data, data for forming a neuron sample.In this case, the calculation of the weights is simpler.It is necessary to solve a system of linear equations, seen that the amount of sample data and variables to be estimated is the same: The result of the training process is neurons of radial basis functions with certain degrees of activation to a network of neurons.From this mesh it is possible to perform a check of the quality of the interpolation, at other sample points at which the values of the interpolated variable are known within the mesh neurons.From the differences between the generated surface and the checkpoints the average and dispersion of data reflecting the quality of the interpolator to mesh interest can be obtained.This quality, expressed in RMS shows the adjustment or the ability of the RBF network to represent real distortions.

STUDY CASE FOR THE STATE OF SÃO PAULO
A total of 338 benchmarks tracked by GNSS were made available by public organizations such as national agencies (IBGE), universities (Universidade Estadual Paulista "Júlio de Mesquita Filho"-UNESP -and the Universidade de São Paulo -USP), and private companies.All GNSS data were available in South America Geocentric Reference System frame (SIRGAS 2000(SIRGAS , epoch 2000.4).4).The geometric height SD provided by organizations is predominantly ± 0.10 m.While ortometric height SD in leveling references, adjusted according to IBGE (2011), is around ± 0.06 and ± 0.08 m. Figure 4 illustrates the standard deviations of both quantities.The statistics of distortions, i.e., the discrepancy between the GNSS leveling data, considering all 338 points, and the geoid models studied (EGM2008, MAPGEO2015 e GEOIDSP2014) are presented in Table I.SD of data between the models is similar, though the accuracy factor (RMS) of EGM2008 is impaired due to a higher trend of the mean in model.MAPGEO2010 and GEOIDSP2014 have a similar consistency through the use of the same method to obtain the geoid model and the long and medium wavelength, also be provided by a geopotential model of the same series, EIGEN.The skewness and kurtosis coefficients indicate proximity that distortions have in relation to a normal distribution.A skewness coefficient close by zero (null) values indicates a symmetric distribution of distortions.In the case of EGM2008 distortions, this coefficient showed less symmetrical distribution with respect to other models.Discrepancies obtained with EGM2008, MAPGEO2015 and GEOIDESP2014 models showed similar platykurtic behavior.
Data distortions are shown in Figure 5 (a) in histograms (left) and punctuality (right).The presence of high discrepancies, some of them bigger than 1 meter, is noticed in data.Since these points can indicate an existence of errors in orthometric and/or geometric heights, it is proposed a filtering for elimination of gross errors or outliers data in order to obtain a hybrid geoid model of the best possible quality.We assume that distortions between geoid models and GNSS data should present a normal distribution when filtering data.Considering the amount of sample points, using the interval [µ ± 3σ] statistically provides a discard of 0.27% of data without non-normal error, equivalent to approximately one point.From this idea, the filtering distortion data was carried out by adopting this range in each model.Also, two steps are proposed for the filtering.First the selection of data in relation to the total sample of data, the second step is to disregard points that have inconsistent distortion with a group of neighboring points, i.e. geometrically adjacent.
In the initial step, intervals are calculated based on the full set of data statistics (Table I), i.e., 338 points tracked by GNSS.The feature of this filter is the elimination of spaced distortions from the averaged data.As a result, the statistical distortions are considerably reduced.Respectively for EGM08, MAPGEO2010 and GEOIDSP2014 the RMS reduced 0.068 m, 0.065 m and 0.066 m; while the maximum distortion values decreased 0.972 m, 0.683 m and 0.700 m.Table II shows the statistics of distortions, after filtering the data.
The second filtering step also used averages and standard deviations, but, checking locally if each distortion is within the region range [µ ± 3σ].However, the selection of the next points requires the definition of some parameters: the number of minimum or maximum points that must be taken into account and/or the minimum or maximum distance to search for neighboring points.
An analysis shows that, considering 2 • distance as the search radius, 94% of data has at least 30 neighbor points, corresponding to approximately 10% of the total amount of data from the study area.However, borders are more isolated in the study area, with greater difficulty in finding neighboring points.Consequently these points need to overreach the search radius to meet a certain minimum number of points.The solution found was to reduce the number of points for the calculation of statistics.It was observed that the use of 10 neighboring points for a search does not exceed a maximum radius of 2 • and 86% of the data points was still found at a distance of at most 1 • .Therefore the additional parameter for step 2 was to consider the first 10 closer neighbors to check whether each distortion is within the region range [µ ± 3σ].
After the second filtering step, new statistical distortions were calculated (Table III).Again the statistics showed reductions.Respectively for EGM08, MAPGEO2010 and GEOIDSP2014, the RMS reduced 0.008 m, 0.015 m and 0.015 m; and the maximum values distortions decreased 0.168 m, 0.262 m and 0.232 m.In addition, comparing statistics before and after the filtering process, we noted higher attenuations.For EGM08, MAPGEO2010 and GEOIDSP2014, respectively, the RMS decreased 0.076 m, 0.080 m and 0.081 m; while the maximum values reduced 1.140 m, 0.945 m and 0.932 m.
The data distortions, after second filtering step, are shown in Figure 5 (b) in histograms (left) and punctuality (right).
A summary of the amount of remaining points after each filtering stage is shown in Table IV.

HYDEGM, HYDMAPGEO AND HYDGEOIDSP
After filtering data of benchmarks tracked by GNSS, three hybrid geoids were obtained from three different gravimetric geoid models where linear RBF was used for corrector surface generation.The first was generated from EGM2008 with degree and order 2160 (HYDEGM), the second model was modeled using MAPGEO2015 (HYDMAPGEO) and the third model using GEOIDSP2014 (HYDGEOIDSP).Step 1 330 331 331 Step 2 319 317 315 Discrepancies between the undulations of hybrid models and benchmarks tracked by GNSS, which did not participate in the model generation (10%), were used for validation.For comparative purposes, the statistics of validation points for gravimetric geoids are shown in Table V.The statistics obtained for the geoid hybrid models, evaluated with the same points used in the validation of their respective gravimetric geoids is shown in Table VI.
Comparing the statistics of gravimetric and hybrid models, it is possible to notice reduction trends: HYDEGM showed a reduction of 0.104 m; HYDMAPGEO decreased 0.035 m; and HYDGEOIDSP indicated an attenuation of 0.047 m.Regarding standard deviations, hybrid models showed reductions of 0.054, 0.043 and 0.062 m, respectively.Consequently, with the application of the proposed methodology, the RMS of hybrid models decreased between 30 and 40% compared to the gravimetric models, indicating model improvements.Another factor that deserves mention is the reduction of the maximum offset values.
In order to evaluate the use of RBF neural networks for modeling of distortions another interpolation method was used while maintaining the same methodology of choice points for generation of surface dis-  The differences between means of hybrid models HYDEGM, HYDMAPGEO and HYDGEOIDSP by linear RBF and Kriging reached 0.008 m; while RMS has differences of 0.024, 0.034 and 0.018 m; and the maximum value differences were between 0.024 and 0.104 m.The hybrid models by linear RBF therefore showed an RMS of 15 to 25% lower compared to hybrid models by Kriging.This comparison shows that the linear RBF neural networks were more suitable in the study case.

CONCLUSIONS
When GNSS positioning is used for leveling purposes, users need to take account of height transformation.So, an accurate geoid model becomes an essential topic and methods and techniques on the geoid calculation are being improved.Since the consistency can also represent data quality, some authors focus on correcting the inconsistency of the geoid model.This method was emphasized in this paper.
Given the current need for a geoid undulation model of better quality for converting GNSS positioning geometric heights in orthometric heights, three geoid models were developed from the EGM08, MAP-GEO2015 and GEOIDSP2014 using the hybridization method in the State of São Paulo, respectively: the HYDEGM the HYDMAPGEO and HYDGEOIDSP.In these models the use of benchmarks tracked by GNSS data is essential in the process of hybridization.338 level references, these adjusted by the IBGE (2011), were used in the study area.However, after carrying out data filtering remained 319, 317 and 315 points for generating the hybrid models, respectively to HYDEGM, HYDMAPGEO and HYDGEOIDSP.
Since a specific interpolation can better represent the spatial reality of the data and the study site, it was found essential the selection of the interpolation method to represent the corrector surface, which must be performed empirically.In this study two interpolators were tested: RBF neural network and Kriging.The RBF interpolation showed 15 to 25% lower RMS by comparison with Kriging.HYDEGM, HYDMAPGEO and HYDGEOIDSP showed, respectively, an RMS of 0.107 m, 0.104 and 0.098 m.
Results of the standard deviations of the distortions between hybrids geoids and leveling data tracked by GNSS reflect consistency between the GNSS positioning and leveling reference points used to validate the models.In conclusion, the errors of the methods of positioning and leveling limit the accuracy of hybrid geoids.

Figure 1 -
Figure 1 -Physical, geoid and ellipsoid surface (see the colors in the online version).

Figure 2 -
Figure 2 -Methodology of Hybrid Geoid Models (see the colors in the online version).
illustrates the offsets ∆N 1 to ∆N 10 (top) and the resultant corrector surface (down) obtained by interpolation.Once a benchmarks network is available, the modeling corrections can be performed by interpolation, estimating correction (∆N ) near the points with known offset.According to the expression 2, positive corrections indicate that the geoid model has a positive trend and negative values indicate a negative trend for the region.

Figure 3 -
Figure 3 -Corrector surface generation (see the colors in the online version).

Figure 5 -
Figure 5 -Distortions (in meters) of geoid models EGM2008, MAPGEO2015 and GEOIDSP2014 regarding benchmarks tracked by GNSS (a) before and (b) after filtering data (see the colors in the online version).
Surface correction of gravimetric geoids EGM2008, MAPGEO2015 and GEOIDSP2014 (top-down) are shown on the left side of Figure6, where dots represent data used during the training phase of linear RBF, and triangles are points used for validation statistics of hybrid geoid models.On the right side in the same figure discrepancies of validation points with benchmarks tracked by GNSS are shown.

Figure 6 -
Figure 6 -Geoid corrector surfaces (left) generated using 90% points randomly selected (dots), and 10% points for checking (triangles), and punctual discrepancy of validation points and benchmarks tracked by GNSS (right) (see the colors in the online version).

TABLE VI Consistency (in meters) of hybrid geoids HYDEGM, HYDMAPGEO and HYDGEOIDSP generated by linear RBF networks regarding validation points.
Kriging was chosen as the interpolation method since this method presents equivalent results to the RBF networks, according to Fazio and Roisenberg (2013).The statistics for validation points concerning hybrid geoids generated by the Kriging interpolation method with linear variogram is shown in TableVII.