Table 3 Statistics of the final gravimetric geoid, quasi-geoid, and their validation [m] modification, which is performed for a specific d/o and then linearly tapered to another higher d/o, all pairs formed from d/o 60 to d/o 300 have been tested. Since FFT needs gridded residual gravity anomalies, the grid was generated based on the irregular residual gravity anomalies over the Kingdom and prediction on a grid with LSC. To evaluate the different gravimetric geoid models resulting from the combination of number of parallels and modification degrees, evalua- tion with a set of available, high-accuracy, GNSS/Leveling dataset by GEOSA was performed. The best results were achieved with a Wong-Gore modification between d/o 80 and 100 and a multiband solution with 3 bands. Figure 3 depicts the final geoid height differences between the final gravimet- ric geoid and the GEOSA GNSS/Leveling geoid heights. In the same processing line, the quasi-geoid over the Kingdom has been determined and from that the geoid was once again estimated using the analytical evaluation of the quasi-geoid to geoid separation by Flury and Rummel (2009). Table 3 tabulates the final gravimetric geoid (see Fig. 2), the quasi- geoid, the difference between the gravimetric geoid and that