Kurzfassung
Um die Auswirkungen von Eismassenvariationen auf das lokale Schwerefeld zu untersuchen, wird ein numerischer Ansatz zur Schwere-Vorwärtsmodellierung entwickelt und vorgestellt. Diese Untersuchungen bauen auf einem synthetisch generierten Gletschermodell für die Nordinsel der Novaya Zemlya Inselgruppe auf, das sowohl die geometrische Struktur als auch die 3D-Dichteverteilung beinhaltet. Durch Modifikationen der Modellparameter wie Eisdicke und Dichteverteilung im Eiskörper werden die zu erwartenden Veränderungen im Schweresignal untersucht. Die modellierte Topographie des Felsuntergrundes kann ebenfalls hinsichtlich unterschiedlicher Annahmen auf Differenzen im resultierenden Schweresignal betrachtet werden. Die Simulationen mit realistisch angenommenen Modellparametern ergeben Gravitationsunterschiede von wenigen mGal. Weiters wird mit Hilfe des Vorwärtsmodellierungsansatzes die Auswirkung der Eismassenveränderungen der letzten 60 Jahre untersucht, die in Form von zwei digitalen Geländemodellen gegeben sind. Der abgeschätzte Effekt auf das Schwerefeld erreicht eine maximale Amplitude von 6 mGal über den gesamten Zeitraum, bzw. eine durchschnittliche Veränderung von ca. 1 mGal pro Jahrzehnt. In weiterer Folge wird in diesem Beitrag ein Konzept vorgestellt, wie Gradientenbeobachtungen der ESA Satellitenmission GOCE für eine regionale Schwerefeldlösung verwendet werden können. Im Gegensatz zur offiziellen Zielsetzung, der Bestimmung eines globalen statischen Schwerefelds basierend auf der Gesamtheit aller Beobachtungen, werden hier die Messungen als Direktbeobachtungen über einem räumlich begrenzten Gebiet eingeführt und die Schwerefeldlösung über die Methode der Kollokation nach kleinsten Quadraten errechnet. Dazu werden die rauschbehafteten Gradientendaten nach dem Wiener-Ansatz gefiltert und die für die Kollokation notwendigen Kovarianzfunktionen abgeleitet.Weiters wird die Problematik des Koordinatenrahmens diskutiert und ein möglicher Lösungsansatz vorgestellt. Mit einem realen GOCE Gradienten Datensatz für November 2009 wird eine Schwerefeldlösung in Form von Schwereanomalien für das oben gente Untersuchungsgebiet berechnet. Mit der verwendeten Meth_x000C_F. Heuberger and D. Rieser: Impact of glacier changes on the local gravity field ... ode und Datenkonfiguration kann das Schwerefeld mit einer geschätzten Genauigkeit von 4 mGal bestimmt werden. Die schwierige Gegenüberstellung der beiden Ansätze (Gravitation aus Vorwärtsmodellierung und Schwereanomalien aus Satellitengradiometrie) wird diskutiert.

Abstract
A numerical approach to gravity forward modelling is developed and introduced in order to investigate the effects of ice mass changes on the local gravity field. These studies are based on a synthetic glacier model of the northern island of Novaya Zemlya, which incorporates geometrical as well as 3D-density information. By modifying the model parameters like ice thickness and the density distribution in the interior of the ice body, the changes that can be expected in the gravity signal are estimated. Furthermore, different assumptions on the underlying bedrock topography can also be evaluated with respect to the resulting gravity signal. Simulations with realistic model parameters yield to gravity attraction differences in the order of a few mGal. Based on given digital elevation models featuring ice mass changes within the last 60 years, the forward modelling approach allows the investigation of the impact of ice change on the gravity field. The estimated effect on the gravity field reaches a maximum amplitude of 6 mGal over the whole period, implying an average change of 1 mGal per decade. In addition, a concept for using gradient observations of ESAs satellite mission GOCE for regional gravity field determination is introduced in this paper. In contrast to the official objectives, i.e. the generation of a global static gravity field based on the entirety of observations, here the measurements are introduced as in situ observations over a spatially restricted area and the gravity field is determined by means of Least Squares Collocation. For this purpose the noisy gradient data are filtered using the Wiener approach and the covariance functions required for collocation are derived. Furthermore, the problematic issue of the coordinate frame is discussed and a possible solution is presented. Finally, a gravity field solution based on real GOCE gradient data for November 2009 is generated for the above mentioned study area in terms of gravity anomalies. With this method and the chosen data configuration it is possible to determine the gravity field with an estimated accuracy of 4 mGal. The difficult comparison of gravity attractions from numerical forward modelling and gravity anomalies from the space-borne gradiometry is discussed.