Description: The ability of any satellite gravity mission concept to monitor mass transport processes in the Earth system is typically tested well ahead of its implementation by means of various simulation studies. Those studies often extend from the simulation of realistic orbits and instrumental data all the way down to the retrieval of global gravity field solution time-series. Basic requirement for all these simulations are realistic representations of the spatio-temporal mass variability in the different sub-systems of the Earth, as a source model for the orbit computations. For such simulations, a suitable source model is required to represent (i) high-frequency (i.e., subdaily to weekly) mass variability in the atmosphere and oceans, in order to realistically include the effects of temporal aliasing due to non-tidal high-frequency mass variability into the retrieved gravity fields. In parallel, (ii) low-frequency (i.e., monthly to interannual) variability needs to be modelled with realistic amplitudes, particularly at small spatial scales, in order to assess to what extent a new mission concept might provide further insight into physical processes currently not observable. The new source model documented here attempts to fulfil both requirements: Based on ECMWF’s recent atmospheric reanalysis ERA-Interim and corresponding simulations from numerical models of the other Earth system components, it offers spherical harmonic coefficients of the time-variable global gravity field due to mass variability in atmosphere, oceans, the terrestrial hydrosphere including the ice-sheets and glaciers, as well as the solid Earth. Simulated features range from sub-daily to multiyear periods with a spatial resolution of spherical harmonics degree and order 180 over a period of 12 years. In addition to the source model, a de-aliasing model for atmospheric and oceanic high-frequency variability with augmented systematic and random noise is required for a realistic simulation of the gravity field retrieval process, whose necessary error characteristics are discussed. The documentation of the updated ESA Earth System Model (updated ESM) for gravity mission simulation studies is organized as follows: The characteristics of the updated ESM along with some basic validation is presented in Volume 1. A detailed comparison to the original ESA ESM (Gruber et al., 2011) is provided in Volume 2, while Volume 3 contains the description of a strategy to derive realistic errors for the de-aliasing model of high-frequency mass variability in atmosphere and ocean.

Description: Der glazial-isostatische Ausgleich in Island infolge des rezenten Abschmelzens der Vatnajökull-Eiskappe wird durch die Viskositätsverteilung im Erdinnern und durch die Details der Abschmelzgeschichte kontrolliert. Interpretationen der Ergebnisse von GPS- und Schweremeßkampagnen im Zeitintervall 1991–2000 bzw. 1992–1999 mit Hilfe lateral homogener Erdmodelle zur Bestimmung der Lithosphärenmächtigkeit, Asthenosphärenmächtigkeit und Asthenosphärenviskosität sind bislang nicht voll zufriedenstellend gewesen. Insbesondere nahe des Eisrandes war die Anpassung der berechneten Landhebung und Schwereänderung an die Beobachtungsdaten nur unzureichend, was mit der Nichtberücksichtigung des Island-Plumes in den lateral homogenen Erdmodellen zusammenhängen kann. In der vorliegenden Arbeit wird für die Modellierung der Landhebung und Schwereänderung ein Programmpaket verwendet, daß die Berechnung auflastinduzierter Störungen eines Maxwellviskoelastischen, inkompressiblen, selbstgravitierenden, sphärischen Erdmodells gestattet. Um das Vorhandensein des Plumes unter dem Vatnajökull zu simulieren, wird eine axialsymmetrische Viskositätsverteilung verwendet, wobei der Plumeradius und die Plumeviskosität freie Parameter sind. Basierend auf seismischen Ergebnissen wird über dem Plume eine 6 km mächtige Lithosphäre angenommen, die sich im peripheren Bereich des Plumes auf 35 km verdickt. Die Abschmelzgeschichte des Vatnajökulls beruht auf Interpretationen geomorphologischer und klimatologischer Untersuchungen und wird durch eine mit dem Plume koaxiale Last mit parabolischem Profil und zeitabhängigem Radius simuliert. Die Ergebnisse der Modellierung favorisieren einen Plumeradius von ~ 80 km und eine Plumeviskosität von (0.3–1.0) × 1018 Pa s.

Description: We review the historical, geological, tide-gauge, GPS and gravimetric evidence advanced in favour of or against continuing land uplift around Hudson Bay, Canada. After this, we reanalyse the tide-gauge and GPS data for Churchill using longer time series than those available to previous investigators. The dependence of the mean rate of relative sea-level change obtained on the length and mid-epoch of the observation interval considered is investigated by means of the newly developed linear-trend analysis diagram. For studying the shorter-period variability of the tide-gauge record, the continuous-wavelet transform is used. The mean rate of land uplift obtained from GPS is based on a new analysis using IGS solutions of GFZ. Furthermore, sea-level indicators from the Churchill region representing the relative sea-level history during the past 8000 a are included. Finally, the four types of observable are jointly inverted in terms of mantle viscosity. The optimum values are 3×10^20 Pa s and 1.6 × 10^22 Pa s for the upper- and lower-mantle viscosities, respectively.

Description: The ability of any satellite gravity mission concept to monitor mass transport processes in the Earth system is typically tested well ahead of its implementation by means of various simulation studies. Those studies often extend from the simulation of realistic orbits and instrumental data all the way down to the retrieval of global gravity field solution time-series. Basic requirement for all these simulations are realistic representations of the spatio-temporal mass variability in the different sub-systems of the Earth, as a source model for the orbit computations. For such simulations, a suitable source model is required to represent (i) high-frequency (i.e., sub-daily to weekly) mass variability in the atmosphere and oceans, in order to realistically include the effects of temporal aliasing due to non-tidal high-frequency mass variability into the retrieved gravity fields. In parallel, (ii) low-frequency (i.e., monthly to interannual) variability needs to be modelled with realistic amplitudes, particularly at small spatial scales, in order to assess to what extent a new mission concept might provide further insight into physical processes currently not observable. The new source model documented here attempts to fulfil both requirements: Based on ECMWF’s recent atmospheric reanalysis ERA-Interim and corresponding simulations from numerical models of the other Earth system components, it offers spherical harmonic coefficients of the time-variable global gravity field due to mass variability in atmosphere, oceans, the terrestrial hydrosphere including the ice-sheets and glaciers, as well as the solid Earth. Simulated features range from sub-daily to multiyear periods with a spatial resolution of spherical harmonics degree and order 180 over a period of 12 years. In addition to the source model, a de-aliasing model for atmospheric and oceanic high-frequency variability with augmented systematic and random noise is required for a realistic simulation of the gravity field retrieval process, whose necessary error characteristics are discussed. The documentation is organized as follows: The characteristics of the updated ESM along with some basic validation are presented in Volume 1 of this report (Dobslaw et al., 2014). A detailed comparison to the original ESA ESM (Gruber et al., 2011) is provided in Volume 2 (Bergmann-Wolf et al., 2014), while Volume 3 (Forootan et al., 2014) contains a description of the strategy to derive a realistically noisy de-aliasing model for the high-frequency mass variability in atmosphere and oceans. The files of the updated ESA Earth System Model for gravity mission simulation studies are accessible at DOI:10.5880/GFZ.1.3.2014.001.

Description: The ability of any satellite gravity mission concept to monitor mass transport processes in the Earth system is typically tested well ahead of its implementation by means of various simulation studies. Those studies often extend from the simulation of realistic orbits and instrumental data all the way down to the retrieval of global gravity field solution time-series. Basic requirement for all these simulations are realistic representations of the spatio-temporal mass variability in the different sub-systems of the Earth, as a source model for the orbit computations. For such simulations, a suitable source model is required to represent (i) high-frequency (i.e., sub-daily to weekly) mass variability in the atmosphere and oceans, in order to realistically include the effects of temporal aliasing due to non-tidal high-frequency mass variability into the retrieved gravity fields. In parallel, (ii) low-frequency (i.e., monthly to interannual) variability needs to be modelled with realistic amplitudes, particularly at small spatial scales, in order to assess to what extent a new mission concept might provide further insight into physical processes currently not observable. The new source model documented here attempts to fulfil both requirements: Based on ECMWF’s recent atmospheric reanalysis ERA-Interim and corresponding simulations from numerical models of the other Earth system components, it offers spherical harmonic coefficients of the time-variable global gravity field due to mass variability in atmosphere, oceans, the terrestrial hydrosphere including the ice-sheets and glaciers, as well as the solid Earth. Simulated features range from sub-daily to multiyear periods with a spatial resolution of spherical harmonics degree and order 180 over a period of 12 years. In addition to the source model, a de-aliasing model for atmospheric and oceanic high-frequency variability with augmented systematic and random noise is required for a realistic simulation of the gravity field retrieval process, whose necessary error characteristics are discussed. The documentation is organized as follows: The characteristics of the updated ESM along with some basic validation are presented in Volume 1 of this report (Dobslaw et al., 2014). A detailed comparison to the original ESA ESM (Gruber et al., 2011) is provided in Volume 2 (Bergmann-Wolf et al., 2014), while Volume 3 (Forootan et al., 2014) contains a description of the strategy to derive a realistically noisy de-aliasing model for the high-frequency mass variability in atmosphere and oceans. The files of the updated ESA Earth System Model for gravity mission simulation studies are accessible at DOI:10.5880/GFZ.1.3.2014.001.

Description: The objective of this study is to calculate the RSL height at the location of the Berkner Island ice core during the last glacial cycle using a viscoelastic earth model and several glacial histories. Earlier glacial cycles are not considered, the basic assumption being that the sediments were deposited during the last glacial cycle. However, to some extent the results can be extrapolated further back into the past. The results impose a constraint on the time of deposition of the sediments. Conversely, if the time of deposition is determined, e. g. by luminescence dating, this study provides a constraint on the glacial history of the Ronne Ice Shelf.

Description: The work contained in this thesis considers deformations of the Earth, which are produced by the loads of the last ice-age glacial sheets. The forces the Earth sets against the surface loads are the buoyancy force of the Earth's mantle and the opposing force by the elastic flexure of the lithosphere. Because the time scale of the ice-age of some 100,000 years is short with respect to geological time scales, the viscoelastic behaviour of the Earth has to be considered. Viscoelasticity results in a retarded response of the Earth, which is observed as postglacial uplift in previously glaciated regions, 8,000 years after deglaciation. To model the buoyancy of the Earth's mantle, often a viscous incompressible fluid of homogeneous density is assumed. More recent studies consider also compressibility of the mantle material, but keep the homogeneous density. This results in an inconsistent reference state, because the self compression due to hydrostatic pressure is neglected. These models are discussed here, and the problems are shown, which arise from the description of the field equations for a viscoelastic compressible gravitating continuum in a half-space geometry. The opposing force by the elastic flexure of the lithosphere is determined by the flexural rigidity of the lithospheric plate. If we consider viscoelastic layers in the lithosphere, the flexural rigidity is reduced. Therefore, the overall thickness of a viscoelastic layered lithosphere is much larger than its effective elastic thickness deduced from assuming one elastic plate. Consequently the effective elastic thickness looses its merit for assessing the lithosphere thickness. We show, how strong effective elastic thickness and lithosphere thickness may differ, in which way the viscoelastic structure of the Earth influences this difference for glacial loads and which consequences arise for the lithospheric stress state.