Local and Regional Geoid Determination
Total Page:16
File Type:pdf, Size:1020Kb
LOCAL AND REGIONAL GEOID DETERMINATION FROM VECTOR AIRBORNE GRAVIMETRY DISSERTATION Presented in Partial Fulfillment of the Requirements for The Degree Doctor of Philosophy in the Graduate School of The Ohio State University By Juan Gilberto Serpas, M.S. ***** The Ohio State University 2003 Dissertation Committee: Approved by Professor Christopher Jekeli, Adviser Professor C. K. Shum ____________________________ Adviser Professor Burkhard Schaffrin Geodetic Science and Surveying ABSTRACT The local geoid in a test area in the Canadian Rocky Mountains is computed using airborne gravimetry data. The geoid is computed by the use of the vertical and horizontal components (VC and HC) of the gravity disturbance vector. In addition, an attempt to combine the three components by the use of least squares collocation is done. The technique of using crossovers to estimate for biases and trends in the gravity signals and the use of minimal control in the form of constraints in the crossover adjustment are studied. Moreover, the downward continuation as well as the direct and indirect effects due to removal and restoration of the masses are investigated. An expression for the effect of the masses applied directly to disturbing potential is provided. Comparison of the predicted components of the gravity disturbing vector with control data indicates that the vertical component is better determined than the horizontal component. The estimated accuracy for the vertical components is on the order of 4 mGal, whereas for the horizontal components it is on the order of 8 to 12 mGal. Both geoid estimates coming from the vertical and horizontal components of the gravity disturbance vector, computed using Hotine’s and line integral, show the same level of accuracy when compared to the Canadian geoid. Relative geoid accuracies on the order of 3 to 7 cm for the VC geoid, and on the order of 4 to 12 cm for the HC geoid are achieved. ii The VC geoid suffers from edge effects on the results, while the HC geoid is highly dependent on ground control. In order to alleviate the use of full ground control for the HC geoid, the computation of the geoid at two crossing tracks is explored. Regarding the estimation of the geoid using least squares collocation to combine the three components of the gravity disturbance vector (3C-LSC), we observe differences in the range of 4 to 6 cm, without including edge effects, with respect to the Canadian geoid. Comparing the 3C-LSC results with those from the VC geoid using Hotine’s integral, the 3C-LSC are comparable and improved for some lines, in terms of standard deviation. In general the result from the 3C-LSC are better than those from HC, by line integral. On the other hand, the use of only the vertical component by least squares collocation (VC-LSC) provides, in general, better results than those from 3C-LSC, and those from the VC and HC by the use of Hotine’s and line integral. We could expect better results for the case of 3C-LSC if we are able to improve the quality of the measured horizontal components of the gravity disturbance vector. The application of a wave correlation filter to both HC and VC component geoid is also explored, and promising results for the improvement of the accuracy of the combined geoid are observed. iii Dedicated to my family iv ACKNOWLEDGMENTS I wish to express my sincere gratitude to my advisor Professor Christopher Jekeli for proposing the topic for my dissertation, for his intellectual support and guidance, and for his patience throughout the development of this research. I also want to acknowledge Professor C.K. Shum and Professor Burkhard Schaffrin for reviewing this study and being part of the evaluation committee. Finally, I want to thank my family for their support and patience. Thank you Maricela, Gilberto Antonio, Juan Marcelo and José Joaquín. Without you and your support and love this work would have not been possible. This work was done with support under the contract No. NMA202-98-1- 1110(NURI98). v VITA October 29, 1959 …………………………… San Miguel, El Salvador 1996 …………………………………………. M.S. Geodetic Science and Surveying, The Ohio State University 1996 – 1998 …………………………………. Researcher and Instructor Universidad Nacional, Costa Rica 1998 – present ………………………………. Graduate Teaching and Research Associate, The Ohio State University PUBLICATIONS Research publications 1. Jekeli, C. and Serpas J. G. (2003) Review and numerical assessment of the direct topographical reduction in geoid determination. Journal of Geodesy. Vol 77, 3-4, pp. 226-239. 2. Serpas Juan G. (2002). Helmert condensation method and the geoid. Proceedings of the Symposium of Geodesy. Dept. of Geodetic Science and Surveying, The Ohio State University, Columbus Ohio, USA. 3. Amador J.A., Burgos A.S., Soley F.G., Serpas G. and Vargas A. (1993) Remote sensing at the Centre for Geophysical Research, University of Costa Rica. Int. J. Remote Sensing. Vol. 14, No.12, 2391-2401. 4. Serpas, J (1993) Ajuste de redes geodesicas. Revista del Colegio federado de ingenieros y arquitectos de Costa Rica . Año 36 No. 6/93, p.46-49. vi FIELDS OF STUDY Major Field: Geodetic Science and Surveying vii TABLE OF CONTENTS PAGE Abstract............................................................................................................................... ii Dedication.......................................................................................................................... iv Acknowledgments............................................................................................................... v Vita..................................................................................................................................... vi List of tables....................................................................................................................... xi List of figures................................................................................................................... xiv Chapters 1. INTRODUCTION .......................................................................................................... 1 1.1 Geoid determination.................................................................................................. 1 1.2 Motivation and background ...................................................................................... 3 1.3 Chapter descriptions.................................................................................................. 6 2. LOCAL AND REGIONAL GEOID DETERMINATION FROM AIRBORNE GRAVIMETRY DATA...................................................................................................... 8 2.1 Introduction............................................................................................................... 8 2.2 The gravity disturbances........................................................................................... 9 2.3 Hotine’s integral for geoid determination............................................................... 10 2.4 Relative geoid by profile integration ...................................................................... 13 2.5 Least Squares Collocation (LSC)............................................................................ 17 2.5.1 Basic formulas ................................................................................................. 17 2.5.2 The covariance matrices .................................................................................. 19 viii 2.5.2.1 Global covariance function....................................................................... 19 2.5.2.2 Local covariance function......................................................................... 21 2.6 Continuation of harmonic functions ....................................................................... 24 2.6.1 The upward continuation ................................................................................. 25 2.6.2 The downward continuation: ........................................................................... 27 2.6.2.1 The gradient method ................................................................................. 27 2.6.2.2 Iterative solution of Poisson’s integral ..................................................... 28 2.6.2.3 Inversion of Poisson’s integral.................................................................. 30 2.7 Reduction of observations for geoid determination................................................ 32 2.7.1 Traditional approach for geoid computation.................................................... 32 2.7.2 Helmert’s condensation method ...................................................................... 35 2.7.3 Helmert’s condensation method as a function of potential.............................. 38 3. NUMERICAL ASPECTS OF GEOID COMPUTATION........................................... 42 3.1 Introduction............................................................................................................. 42 3.2 The Fourier Transform and the FFT ....................................................................... 42 3.3 Numerical implementation of Hotine’s integral ..................................................... 47 3.4 Numerical implementation of the line integral for relative geoid determination ... 52 3.5 Numerical aspect for LSC....................................................................................... 53 3.6 Numerical treatment for the downward continuation ............................................. 58 3.6.1 The gradient