### Refine

#### Year of publication

- 2004 (2) (remove)

#### Document Type

- Doctoral Thesis (1)
- Preprint (1)

#### Language

- English (2) (remove)

#### Keywords

- Gravitationsfeld (2) (remove)

Nowadays one of the major objectives in geosciences is the determination of the gravitational field of our planet, the Earth. A precise knowledge of this quantity is not just interesting on its own but it is indeed a key point for a vast number of applications. The important question is how to obtain a good model for the gravitational field on a global scale. The only applicable solution - both in costs and data coverage - is the usage of satellite data. We concentrate on highly precise measurements which will be obtained by GOCE (Gravity Field and Steady State Ocean Circulation Explorer, launch expected 2006). This satellite has a gradiometer onboard which returns the second derivatives of the gravitational potential. Mathematically seen we have to deal with several obstacles. The first one is that the noise in the different components of these second derivatives differs over several orders of magnitude, i.e. a straightforward solution of this outer boundary value problem will not work properly. Furthermore we are not interested in the data at satellite height but we want to know the field at the Earth's surface, thus we need a regularization (downward-continuation) of the data. These two problems are tackled in the thesis and are now described briefly. Split Operators: We have to solve an outer boundary value problem at the height of the satellite track. Classically one can handle first order side conditions which are not tangential to the surface and second derivatives pointing in the radial direction employing integral and pseudo differential equation methods. We present a different approach: We classify all first and purely second order operators which fulfill that a harmonic function stays harmonic under their application. This task is done by using modern algebraic methods for solving systems of partial differential equations symbolically. Now we can look at the problem with oblique side conditions as if we had ordinary i.e. non-derived side conditions. The only additional work which has to be done is an inversion of the differential operator, i.e. integration. In particular we are capable to deal with derivatives which are tangential to the boundary. Auto-Regularization: The second obstacle is finding a proper regularization procedure. This is complicated by the fact that we are facing stochastic rather than deterministic noise. The main question is how to find an optimal regularization parameter which is impossible without any additional knowledge. However we could show that with a very limited number of additional information, which are obtainable also in practice, we can regularize in an asymptotically optimal way. In particular we showed that the knowledge of two input data sets allows an order optimal regularization procedure even under the hard conditions of Gaussian white noise and an exponentially ill-posed problem. A last but rather simple task is combining data from different derivatives which can be done by a weighted least squares approach using the information we obtained out of the regularization procedure. A practical application to the downward-continuation problem for simulated gravitational data is shown.

The following two papers present recent developments in multiscale ocean circulation modeling and multiscale gravitational field modeling that have been presented at the 2nd International GOCE User Workshop 2004 in Frascati. Part A - Multiscale Modeling of Ocean Circulation In this paper the applicability of multiscale methods to oceanography is demonstrated. More precisely, we use convolutions with certain locally supported kernels to approximate the dynamic topography and the geostrophic flow. As data sets the French CLS01 data are used for the mean sea surface topography and are compared to the EGM96 geoid. Since those two data sets have very different levels of spatial resolutions the necessity of an interpolating or approximating tool is evident. Compared to the standard spherical harmonics approach, the strongly space localizing kernels improve the possibilities of local data analysis here. Part B - Multiscale Modeling from EIGEN-1S, EIGEN-2, EIGEN-GRACE01S, GGM01, UCPH2002_0.5, EGM96 Spherical wavelets have been developed by the Geomathematics Group Kaiserslautern for several years and have been successfully applied to georelevant problems. Wavelets can be considered as consecutive band-pass filters and allow local approximations. The wavelet transform can also be applied to spherical harmonic models of the Earth's gravitational field like the most up-to-date EIGEN-1S, EIGEN-2, EIGEN-GRACE01S, GGM01, UCPH2002_0.5, and the well-known EGM96. Thereby, wavelet coefficients arise. In this paper it is the aim of the Geomathematics Group to make these data available to other interested groups. These wavelet coefficients allow not only the reconstruction of the wavelet approximations of the gravitational potential but also of the geoid, of the gravity anomalies and other important functionals of the gravitational field. Different types of wavelets are considered: bandlimited wavelets (here: Shannon and Cubic Polynomial (CuP)) as well as non-bandlimited ones (in our case: Abel-Poisson). For these types wavelet coefficients are computed and wavelet variances are given. The data format of the wavelet coefficients is also included.