A Regional Gravimetric and Hybrid Geoid Model in Northern Greece from Dedicated Gravity Campaigns

The determination of physical heights is of key importance for a wide spectrum of geoscientific applications and, in particular, for engineering projects. The main scope of the present work is focused on the determination of a high-accuracy and high-resolution gravimetric and hybrid geoid model, to determine orthometric heights without the need of conventional leveling. Both historical and newly acquired gravity data have been collected during dedicated gravity campaigns, around the location of a dedicated GNSS network as well as areas where the existing land gravity database presented voids. Geoid determination was based on the classical remove–compute–restore (RCR) technique and spectral and stochastic approaches. The low frequencies have been modeled based on the XGM2019e global geopotential model (GGM) and the topographic effects have been evaluated with the residual terrain model (RTM) reduction. The evaluation of the final geoid model was performed over 462 GNSS/leveling benchmarks (BMs), where the newly determined gravimetric geoid has shown an improvement of 3.1 cm, in the std of the differences to the GNSS/leveling BMs, compared to the latest national geoid model. A deterministic and stochastic fit to the GNSS/leveling data has been performed, investigating various choices for the parametric models and analytical covariance functions. The scope was to determine a hybrid geoid model, tailored to the area and GNSS/leveling data, which will be the one used for the direct estimation of high-accuracy orthometric heights from GNSS observations. After the deterministic fit, a std to the GNSS/leveling data of 10.1 cm has been achieved, with 54.8% and 83.1% of the absolute height differences being below the 1 cm and 2 cm per square root km of baseline length. The final hybrid geoid model, i.e., after the stochastic treatment of the adjusted residuals, gave a std of the difference to the GNSS/leveling data of 1.1 cm, with 99.8% and 99.9% of the height difference being smaller than the 1 cm and 2 cm standard errors, thus achieving a 1 cm accuracy regional geoid.