Quasigeoid modelling in New Zealand to unify multiple local vertical datums

Abstract

One goal of modern geodesy is the global unification of vertical datums so that height data from them can be properly integrated. This thesis studies the unification of the 13 disparate levelling- and tide-gauge-based vertical datums in New Zealand (NZ). It proposes a new NZ-wide single vertical datum based on a gravimetric quasigeoid model to unify the existing local vertical datums. This will also include methods to transform height data in terms of the existing datums to the new datum and vice versa. After defining and comparing the main types of height system and vertical datum used around to world, the system of heights used in NZ was shown to be normal-orthometric. Consequently, datum unification was achieved using a quasigeoid model, as opposed to a geoid model. The quasigeoid was computed by combining the GRACE-based GGM02 and EGM96 global geopotential models with land gravity data (40,737 observations) and a 56-m resolution digital elevation model (DEM). Marine gravity data came from a least-squares collocation combination of 1,300,266 crossover-adjusted ship track observations and gravity anomalies derived from multi-mission satellite altimetry.To ensure that the best quasigeoid was computed for the NZ datasets, a number of computation processes were compared and contrasted. The Hammer chart, fast Fourier transform (FFT) and prism integration methods of computing terrain corrections (TCs) were compared. This showed that the prism integration TC is the most realistic in NZ. The mean Helmert gravity anomalies, required for numerical integration of Stokes’s formula, were computed via refined Bouguer anomalies with the prism TCs, and reconstruction with heights from the DEM used to ‘reconstruct’ more representative mean anomalies. In addition, five deterministic modifications to Stokes’s formula were compared. There was little difference between three of them, so the Featherstone et al. (1998) modification ( 0 y = 1.5°, M = 40) was chosen because it is theoretically better than its predecessors. The global geopotential, gravimetric geoid, sea surface topography and geodetic boundary-value problem approaches to vertical datum unification were then contrasted. As none was likely to be effective in NZ, a new iterative quasigeoid approach was adopted. This procedure computes an initial quasigeoid from existing data on the various local vertical datums to estimate the vertical datum offsets from co-located GPS-levelling data. These offsets were then subsequently applied to the gravity observations by way of additional reductions to the initially computed quasigeoid to reduce gravity anomaly biases caused by the vertically offset datums. These adjusted gravity anomalies were then used to compute a new quasigeoid model, and the process repeated until the computed offsets between the local vertical datums (or equivalently two quasigeoid solutions) converged, which took only two iterations.The computed offsets were then compared with spirit-levelled height differences among adjoining datums, where these were available, giving an average agreement of 7 cm. Since the local vertical datums are effectively unified, the new national vertical datum for NZ will comprise the iteratively computed gravimetric quasigeoid model, accompanied by local vertical datums. This approach is implemented to give a new national vertical datum for NZ. When used with the appropriate offset, this enables the transformation of heights in terms of the national vertical datum to the 13 precise-levelling datums and the ellipsoidal national geodetic datum, NZ Geodetic Datum 2000.