In this paper, two techniques for calculating the geoid-to-quasigeoid separation are employed. One of them is GPS/Levelling customary method as a criterion where the geoid undulation and height anomaly are computed by subtracting the ellipsoid height attained via GPS from the orthometric height and normal height, respectively. Another approach is Sjöberg’s equation. We have used of the ICGEM website for definition of the variables of the Sjöberg’s equation, as the applied reference model is the EGM2008 global geopotential model and WGS84 reference ellipsoid. The investigations are performed over the stations of the GPS/Leveling network related to three selected areas in desert, mountain and flatland namely the Lout, Zagros and Khuzestan in Iran and afterward the correlation coefficient between the geoid-to-quasigeoid separation calculated using the satellite data in Sjöberg’s equation and GPS/Levelling method is estimated. The results indicate a straight correlation between the estimated separations from the two methods as its value for the Lout is 0.754, for the Zagros is 0.497 and for the Khuzestan is 0.659. consequently, using the satellite data in Sjöberg’s equation for the regions where there are not the GPS/Levelling and land gravity data, specially for the even areas, yield a satisfactory response of the geoidto-quasigeoid separation.

An increased use of global navigation techniques for positioning, and in particular for height determination, led to a growing need for precise models of height reference surface, i.e. geoid or quasigeoid. Geoid or quasigeoid heights at a cm accuracy level, provided on growing number ofGPS/levelling sites, can not only be used for quality control of gravimetric geoid but they also can be integrated with gravity data for geoid/quasigeoid modelling. Such a model is of particular use for surveying practice. A method of quasigeoid modelling based on GPS/levelling data with support of geopotential model and gravity data was developed. The components of height anomaly are modelled with the deterministic part that consists of height anomaly based on EGM96 geopotential model and Molodensky's integral, as well as the polynomial representing trend, and from the stochastic part represented by the isotropic covariance function. Model parameters, i.e. polynomial coefficients and covariance function parameters are determined in a single process of robust estimation, resistant to the outlying measurements. The method was verified using almost a thousand height anomalies from the sites of the EUREF-POL, POLREF, EUVN'97 and WSSG (Military Satellite Geodetic Network) networks in Poland as well as geopotential model refined with gravity data in l' x l' grid. The estimated average mean square error of quasigeoid height is at the level of O.Ol m. The outlying measurements were efficiently detected.