## Harmonic Spline-Wavelets on the 3-dimensional Ball and their Application to the Reconstruction of the Earth´s Density Distribution from Gravitational Data at Arbitrarily Shaped Satellite Orbits

- We introduce splines for the approximation of harmonic functions on a 3-dimensional ball. Those splines are combined with a multiresolution concept. More precisely, at each step of improving the approximation we add more data and, at the same time, reduce the hat-width of the used spline basis functions. Finally, a convergence theorem is proved. One possible application, that is discussed in detail, is the reconstruction of the Earth´s density distribution from gravitational data obtained at a satellite orbit. This is an exponentially ill-posed problem where only the harmonic part of the density can be recovered since its orthogonal complement has the potential 0. Whereas classical approaches use a truncated singular value decomposition (TSVD) with the well-known disadvantages like the non-localizing character of the used spherical harmonics and the bandlimitedness of the solution, modern regularization techniques use wavelets allowing a localized reconstruction via convolutions with kernels that are only essentially large in the region of interest. The essential remaining drawback of a TSVD and the wavelet approaches is that the integrals (i.e. the inner product in case of a TSVD and the convolution in case of wavelets) are calculated on a spherical orbit, which is not given in reality. Thus, simplifying modelling assumptions, that certainly include a modelling error, have to be made. The splines introduced here have the important advantage, that the given data need not be located on a sphere but may be (almost) arbitrarily distributed in the outer space of the Earth. This includes, in particular, the possibility to mix data from different satellite missions (different orbits, different derivatives of the gravitational potential) in the calculation of the Earth´s density distribution. Moreover, the approximating splines can be calculated at varying resolution scales, where the differences for increasing the resolution can be computed with the introduced spline-wavelet technique.

Author: | Martin J. Fengler, Dominik Michel, Volker Michel |
---|---|

URN (permanent link): | urn:nbn:de:hbz:386-kluedo-13714 |

Serie (Series number): | Schriften zur Funktionalanalysis und Geomathematik (16) |

Document Type: | Preprint |

Language of publication: | English |

Year of Completion: | 2005 |

Year of Publication: | 2005 |

Publishing Institute: | Technische Universität Kaiserslautern |

Tag: | GOCE <Satellitenmission>; GRACE <Satellitenmission> ; Kugel ; Spline-Wavelets ; reguläre Fläche GOCE <satellite mission>; GRACE <satellite mission> ; ball ; regular surface ; spline-wavelets |

GND-Keyword: | CHAMP <Satellitenmission>; Harmonische Spline-Funktion ; Inverses Problem ; Kugelfunktion ; Mehrskalenanalyse ; Regularisierung ; Sobolev-Raum |

Faculties / Organisational entities: | Fachbereich Mathematik |

DDC-Cassification: | 510 Mathematik |

MSC-Classification (mathematics): | 41A15 Spline approximation |

42C40 Wavelets and other special systems | |

45K05 Integro-partial differential equations [See also 34K30, 35R09, 35R10, 47G20] | |

86A22 Inverse problems [See also 35R30] |