Fitting a triaxial ellipsoid to a geoid model

The aim of this work is the determination of the parameters of Earth’s triaxiality through a geometric fitting of a triaxial ellipsoid to a set of given points in space, as they are derived from a geoid model. Starting from a Cartesian equation of an ellipsoid in an arbitrary reference system, we develop a transformation of its coefficients into the coordinates of the ellipsoid center, of the three rotation angles and the three ellipsoid semi-axes. Furthermore, we present different mathematical models for some special and degenerate cases of the triaxial ellipsoid. We also present the required mathematical background of the theory of least-squares, under the condition of minimization of the sum of squares of geoid heights. Also, we describe a method for the determination of the foot points of the set of given space points. Then, we prepare suitable data sets and we derive results for various geoid models, which were proposed in the last fifty years. Among the results, we report the semi-axes of the triaxial ellipsoid of geometric fitting with four unknowns to be 6378171.92 m, 6378102.06 m and 6356752.17 m and the equatorial longitude of the major semi-axis –14.9367 degrees. Also, the parameters of Earth’s triaxiality are directly estimated from the spherical harmonic coefficients of degree and order two. Finally, the results indicate that the geoid heights with reference to the triaxial ellipsoid are smaller than those with reference to the oblate spheroid and the improvement in the corresponding rms value is about 20 per cent