Description: Accurate quantification of the millennial-scale mass balance of the Greenland ice sheet (GrIS) and its contribution to global sea-level rise remain challenging because of sparse in situ observations in key regions. Glacial isostatic adjustment (GIA) is the ongoing response of the solid Earth to ice and ocean load changes occurring since the Last Glacial Maximum (LGM; ~21 thousand years ago) and may be used to constrain the GrIS deglaciation history. We use data from the Greenland Global Positioning System network to directly measure GIA and estimate basin-wide mass changes since the LGM. Unpredicted, large GIA uplift rates of +12 mm/year are found in southeast Greenland. These rates are due to low upper mantle viscosity in the region, from when Greenland passed over the Iceland hot spot about 40 million years ago. This region of concentrated soft rheology has a profound influence on reconstructing the deglaciation history of Greenland. We reevaluate the evolution of the GrIS since LGM and obtain a loss of 1.5-m sea-level equivalent from the northwest and southeast. These same sectors are dominating modern mass loss. We suggest that the present destabilization of these marine-based sectors may increase sea level for centuries to come. Our new deglaciation history and GIA uplift estimates suggest that studies that use the Gravity Recovery and Climate Experiment satellite mission to infer present-day changes in the GrIS may have erroneously corrected for GIA and underestimated the mass loss by about 20 gigatons/year.

Description: Accurate quantification of the millennial-scale mass balance of the Greenland ice sheet (GrIS) and its contribution to global sea-level rise remain challenging because of sparse in situ observations in key regions. Glacial isostatic adjustment (GIA) is the ongoing response of the solid Earth to ice and ocean load changes occurring since the Last Glacial Maximum (LGM; ~21 thousand years ago) and may be used to constrain the GrIS deglaciation history. We use data from the Greenland Global Positioning System network to directly measure GIA and estimate basin-wide mass changes since the LGM. Unpredicted, large GIA uplift rates of +12 mm/year are found in southeast Greenland. These rates are due to low upper mantle viscosity in the region, from when Greenland passed over the Iceland hot spot about 40 million years ago. This region of concentrated soft rheology has a profound influence on reconstructing the deglaciation history of Greenland. We reevaluate the evolution of the GrIS since LGM and obtain a loss of 1.5-m sea-level equivalent from the northwest and southeast. These same sectors are dominating modern mass loss. We suggest that the present destabilization of these marine-based sectors may increase sea level for centuries to come. Our new deglaciation history and GIA uplift estimates suggest that studies that use the Gravity Recovery and Climate Experiment satellite mission to infer present-day changes in the GrIS may have erroneously corrected for GIA and underestimated the mass loss by about 20 gigatons/year.

Description: In this study, a new method for computing the sensitivity of the glacial isostatic adjustment (GIA) forward solution with respect to the Earth's mantle viscosity, the so-called the forward sensitivity method (FSM), and a method for computing the gradient of data misfit with respect to viscosity parameters, the so-called adjoint-state method (ASM), are presented. These advanced formal methods complement each other in the inverse modelling of GIA-related observations. When solving this inverse problem, the first step is to calculate the forward sensitivities by the FSM and use them to fix the model parameters that do not affect the forward model solution, as well as identifying and removing redundant parts of the inferred viscosity structure. Once the viscosity model is optimized in view of the forward sensitivities, the minimization of the data misfit with respect to the viscosity parameters can be carried out by a gradient technique which makes use of the ASM. The aim is this paper is to derive the FSM and ASM in the forms that are closely associated with the forward solver of GIA developed by Martinec. Since this method is based on a continuous form of the forward model equations, which are then discretized by spectral and finite elements, we first derive the continuous forms of the FSM and ASM and then discretize them by the spectral and finite elements used in the discretization of the forward model equations. The advantage of this approach is that all three methods (forward, FSM and ASM) have the same matrix of equations and use the same methodology for the implementation of the time evolution of stresses. The only difference between the forward method and the FSM and ASM is that the different numerical differencing schemes for the time evolution of the Maxwell and generalized Maxwell viscous stresses are applied in the respective methods. However, it requires only a little extra computational time for carrying out the FSM and ASM numerically. An straightforward approach to compute the gradient of the data misfit is the brute-force method, whereby the partial derivatives of the misfit with respect to model parameters are approximated by the centred difference of two forward model runs. Although the brute-force method is useful for computing the gradient of the data misfit with respect to a small number of model parameters, it becomes expensive for a viscosity model with a large number of parameters. The ASM offers an efficient alternative for computing the gradient of the misfit since the computational time of the ASM is independent of the number of viscosity parameters. The ASM is thus highly efficient for calculating the gradient of the misfit for models with large numbers of parameters. However, the forward-model solution for each time step must be stored, hence the memory demands scale linearly with the number of time steps. This is the main drawback of the ASM.

Description: The GFZ German Research Centre for Geosciences as part of the GRACE Science Data System (SDS) is currently reprocessing the complete GRACE mission data. This new Level-2 data release (RL06 in the SDS nomenclature) will be based on reprocessed Level-1B instrument data (RL03), updated processing standards and background models and will take care of limitations known from previous RL05. Examples are the application of the latest RL06 Atmosphere and Ocean Dealiasing Model, update of the ocean tide model, implementation of the most recent IERS conventions or improvements in GFZ´s GPS data processing. This 15+ year time series of monthly Level-2 spherical harmonics and underlying processing standards will then serve for the continuation with GRACE-FO (Follow-on) data expected for early 2018. In parallel a team of GFZ, the Alfred-Wegener-Institute Bremerhaven and TU Dresden has developed and implemented a portal at GFZ where users can download dedicated Level-3 products for hydrological, oceanic and polar research activities. This portal is expected to be made public by the end of 2017. The presentation will show the status and examples of these new RL06 Level-2 products and prototype Level-3 products based on GFZ’s RL05a Level-2 monthly solutions.

Description: Time-resolved satellite gravimetry has revolutionized understanding of mass transport in the Earth system. Since 2002, the Gravity Recovery and Climate Experiment (GRACE) has enabled monitoring of the terrestrial water cycle, ice sheet and glacier mass balance, sea level change and ocean bottom pressure variations, as well as understanding responses to changes in the global climate system. Initially a pioneering experiment of geodesy, the time-variable observations have matured into reliable mass transport products, allowing assessment and forecast of a number of important climate trends, and improvements in service applications such as the United States Drought Monitor. With the successful launch of the GRACE Follow-On mission, a multi-decadal record of mass variability in the Earth system is within reach.

Description: We employ a coupled model for ice-sheet dynamics and Maxwell viscoelastic solid-Earth dynamics, including a gravitationally consistent description of sea level. With this model, we study the influence of the solid Earth on the future evolution of the West Antarctic Ice Sheet (WAIS). Starting from steady-state conditions close to the present-day configuration of the Antarctic Ice Sheet, we apply different atmospheric and oceanic forcings and solid-Earth rheologies in order to analyse the retreat of the WAIS. Climate forcing is the primary control on the occurrence of WAIS collapse. For moderate climate forcing and weak solid-Earth rheologies, however, we find that the relative sea level (RSL) fall associated with the viscoelastic solid-Earth response due to unloading by WAIS retreat limits the retreat to the Amundsen Sea embayment on time scales of several millennia, whereas stiffer Earth structures yield a collapse under these conditions. Under stronger climate forcing, weak Earth structures associated with the West Antarctic rift system produce a delay of up to 5000 years in comparison to stiffer, Antarctic-average solid-Earth rheologies. Furthermore, we find that sea-level rise from an assumed fast deglaciation of the Greenland Ice Sheet induces WAIS collapse in the presence of higher asthenosphere viscosities in cases when the climatic forcing is too weak to force WAIS collapse alone.

Description: Pine Island Glacier (PIG) in the Amundsen Sea sector of the West Antarctic Ice Sheet (WAIS) is losing mass and contributing to global sea-level rise at an accelerating rate. Although recent observations and modeling have identified the incursion of relatively warm Circumpolar Deep Water (CDW) beneath the PIG ice shelf (PIGIS) as the main driver of this ice-mass loss, the lack of precise bathymetry limits furthering our understanding of the ice-ocean interactions and improving the accuracy of modeling. Here we present updated bathymetry and sediment distribution beneath the PIGIS, modeled by the inversion of aerogravity data with constraints from active-source seismic data, observations from an autonomous underwater vehicle, and the regional gravity-anomaly field derived from satellite gravity observations. Modeled bathymetry shows a submarine ridge beneath the middle of PIGIS that rises ∼350 to 400 m above the surrounding sea floor, with a minimum water-column thickness of ∼200 m above it. This submarine ridge continues across the whole width of the 45-km wide ice shelf, with no deep troughs crossing it, confirming the general features of the previously predicted sub-ice-shelf ocean circulation. However, the relatively low resolution of the aerogravity data and limitations in our inversion method leave a possibility that there is an undetected, few-kilometers-wide or narrower trough that may alter the predicted sub-ice-shelf ocean circulation. Modeled sediment distribution indicates a sedimentary basin of up to ∼800 m thick near the current grounding zone of the main PIG trunk and extending farther inland, and a region seaward of the submarine ridge where sediments are thin or absent with exposed crystalline basement that extends seaward into Pine Island Bay. Therefore, the submarine ridge marks the transition from a thick sedimentary basin providing a smooth interface over which ice could flow easily by sliding or sediment deformation, to a region with no to little sediments and instead a rough interface over which ice flows mainly by deformation. We hypothesize that the post-Last Glacial Maximum retreat of PIG stabilized at this location because of the spatial transition in basal conditions. This in turn supports the hypothesis that the recent retreat of PIG was strongly forced, probably by changes in ocean circulation, rather than occurring because of ongoing response to the end of the ice age or other changes inland of or beneath PIG.

Description: The GFZ German Research Centre for Geosciences, in collaboration with the Alfred-Wegener-Institute (AWI) Bremerhaven and Technische Universität Dresden, has developed and implemented a new web portal called GravIS (Gravity Information Service). At GravIS, user-friendly ("Level-3") products based on the most recent GFZ GRACE (and in near future also GRACE Follow-on) Level-2 release of monthly global gravity field solutions are visualized and described. These Level-3 products comprise globally gridded mass anomalies as well as basin average time series and are available for three dedicated applications: (1) terrestrial water storage over non-glaciated regions; (2) bottom pressure variations in ocean basins, and; (3) ice mass changes in both Antarctica and Greenland. Several postprocessing steps have been applied to the Level-2 spherical harmonic coefficients before mass anomaly inversion to achieve the highest possible accuracy of the Level-3 products. The poster gives an overview of the GravIS web portal and examples of the latest Level-3 products are shown. All Level-3 products will also be available for download via GFZ’s data archive ISDC.

Description: A new synthetic model of the time-variable global gravity field is now available based on realistic mass variability in atmosphere, oceans, terrestrial water storage, continental ice-sheets, and the solid Earth. The updated ESA Earth System Model is provided in Stokes coefficients up to degree and order 180 with a temporal resolution of 6 h covering the time period 1995–2006, and can be readily applied as a source model in future gravity mission simulation studies. The model contains plausible variability and trends in both low-degree coefficients and the global mean eustatic sea level. It depicts reasonable mass variability all over the globe at a wide range of frequencies including multi-year trends, year-to-year variability, and seasonal variability even at very fine spatial scales, which is important for a realistic representation of spatial aliasing and leakage. In particular on these small spatial scales between 50 and 250 km, the model contains a range of signals that have not been reliably observed yet by satellite gravimetry. In addition, the updated Earth System Model provides substantial high-frequency variability at periods down to a few hours only, thereby allowing to critically test strategies for the minimization of temporal aliasing.