Variance–covariance analysis of two high-resolution regional least-squares quasi-geoid models

This paper investigates the full variance–covariance (VC) matrix of two high-resolution regional quasi-geoid models, utilizing a spherical radial basis function parameterization. Model parameters were estimated using weighted least-squares techniques and variance component estimation (VCE) for data weighting. The first model, known as the “RCR model,” is computed through the remove–compute–restore method, incorporating various local gravity and radar altimeter datasets. The second model, the “combined model,” includes the GOCO05s satellite-only global geopotential model as an additional dataset with a full-noise VC matrix. Validation of the noise VC matrix scaling for each quasi-geoid model is achieved by comparing observed and formal noise standard deviations of differences between geometric and gravimetric height anomalies at GPS height markers in the Netherlands. Analysis of the noise VC matrix of height anomalies at grid nodes reveals significantly smaller formal noise standard deviations for the RCR model compared to the combined model. This difference is attributed to VCE assigning larger weights to the GOCO05s dataset, which exhibits greater noise standard deviations for the specific spatial scales used. Additionally, the formal noise standard deviations of height anomaly differences, relevant for GNSS-heighting, favor the RCR model. However, the disparity between the two models is smaller than implied by the height anomaly noise standard deviations. This is due to the combined model’s noise autocorrelation function displaying a longer correlation length (67 km) in contrast to the RCR model’s (17 km). Consequently, the combined model exhibits a greater reduction in noise variance for height anomaly differences relative to white noise compared to the RCR model.