scholarly journals Altimetry, gravimetry, GPS and viscoelastic modelling data for the joint inversion for glacial isostatic adjustment in Antarctica (ESA STSE Project REGINA)

2017 ◽  
Author(s):  
Ingo Sasgen ◽  
Alba Martín-Español ◽  
Alexander Horvath ◽  
Volker Klemann ◽  
Elizabeth J. Petrie ◽  
...  
Keyword(s):  
Solid Earth ◽  
Joint Inversion ◽  
Ice Sheet ◽  
Data Sets ◽  
Glacial Isostatic Adjustment ◽  
Response Functions ◽  
Forward Modelling ◽  
Satellite Gravimetry ◽  
Isostatic Adjustment ◽  
The Antarctic

Abstract. A major uncertainty in determining the mass balance of the Antarctic ice sheet from measurements of satellite gravimetry, and to a lesser extent satellite altimetry, is the poorly known correction for the ongoing deformation of the solid Earth caused by glacial isostatic adjustment (GIA). In the past decade, much progress has been made in consistently modelling the ice sheet and solid Earth interactions; however, forward-modelling solutions of GIA in Antarctica remain uncertain due to the sparsity of constraints on the ice sheet evolution, as well as the Earth's rheological properties. An alternative approach towards estimating GIA is the joint inversion of multiple satellite data – namely, satellite gravimetry, satellite altimetry and GPS, which reflect, with different sensitivities, trends of recent glacial changes and GIA. Crucial to the success of this approach is the accuracy of the space-geodetic data sets. Here, we present reprocessed rates of surface-ice elevation change (Envisat/ICESat; 2003–2009), gravity field change (GRACE; 2003–2009) and bedrock uplift (GPS; 1995–2013.7). The data analysis is complemented by the forward-modelling of viscoelastic response functions to disc load forcing, allowing us to relate GIA-induced surface displacements with gravity changes for different rheological parameters of the solid Earth. The data and modelling results presented here are available in the Pangea archive; https://doi.pangaea.de/10.1594/PANGAEA.875745. The data sets are the input streams for the joint inversion estimate of present-day ice-mass change and GIA, focusing on Antarctica. However, the methods, code and data provided in this paper are applicable to solve other problems, such as volume balances of the Antarctic ice sheet, or to other geographical regions, in the case of the viscoelastic response functions. This paper presents the first of two contributions summarizing the work carried out within a European Space Agency funded study, REGINA.

scholarly journals Altimetry, gravimetry, GPS and viscoelastic modeling data for the joint inversion for glacial isostatic adjustment in Antarctica (ESA STSE Project REGINA)

2018 ◽  
Vol 10 (1) ◽  
pp. 493-523 ◽  
Author(s):  
Ingo Sasgen ◽  
Alba Martín-Español ◽  
Alexander Horvath ◽  
Volker Klemann ◽  
Elizabeth J. Petrie ◽  
...  
Keyword(s):  
Solid Earth ◽  
Joint Inversion ◽  
Ice Sheet ◽  
Forward Modeling ◽  
Data Sets ◽  
Glacial Isostatic Adjustment ◽  
Response Functions ◽  
Satellite Gravimetry ◽  
Isostatic Adjustment ◽  
The Antarctic

Abstract. The poorly known correction for the ongoing deformation of the solid Earth caused by glacial isostatic adjustment (GIA) is a major uncertainty in determining the mass balance of the Antarctic ice sheet from measurements of satellite gravimetry and to a lesser extent satellite altimetry. In the past decade, much progress has been made in consistently modeling ice sheet and solid Earth interactions; however, forward-modeling solutions of GIA in Antarctica remain uncertain due to the sparsity of constraints on the ice sheet evolution, as well as the Earth's rheological properties. An alternative approach towards estimating GIA is the joint inversion of multiple satellite data – namely, satellite gravimetry, satellite altimetry and GPS, which reflect, with different sensitivities, trends in recent glacial changes and GIA. Crucial to the success of this approach is the accuracy of the space-geodetic data sets. Here, we present reprocessed rates of surface-ice elevation change (Envisat/Ice, Cloud,and land Elevation Satellite, ICESat; 2003–2009), gravity field change (Gravity Recovery and Climate Experiment, GRACE; 2003–2009) and bedrock uplift (GPS; 1995–2013). The data analysis is complemented by the forward modeling of viscoelastic response functions to disc load forcing, allowing us to relate GIA-induced surface displacements with gravity changes for different rheological parameters of the solid Earth. The data and modeling results presented here are available in the PANGAEA database (https://doi.org/10.1594/PANGAEA.875745). The data sets are the input streams for the joint inversion estimate of present-day ice-mass change and GIA, focusing on Antarctica. However, the methods, code and data provided in this paper can be used to solve other problems, such as volume balances of the Antarctic ice sheet, or can be applied to other geographical regions in the case of the viscoelastic response functions. This paper presents the first of two contributions summarizing the work carried out within a European Space Agency funded study: Regional glacial isostatic adjustment and CryoSat elevation rate corrections in Antarctica (REGINA).

scholarly journals Estimation of present-day glacial isostatic adjustment, ice mass change and elastic vertical crustal deformation over the Antarctic ice sheet

2017 ◽  
Vol 63 (240) ◽  
pp. 703-715 ◽  
Author(s):  
BAOJUN ZHANG ◽  
ZEMIN WANG ◽  
FEI LI ◽  
JIACHUN AN ◽  
YUANDE YANG ◽  
...  
Keyword(s):  
Antarctic Peninsula ◽  
Crustal Deformation ◽  
Ice Sheet ◽  
East Antarctica ◽  
Mass Change ◽  
Glacial Isostatic Adjustment ◽  
Antarctic Ice Sheet ◽  
Isostatic Adjustment ◽  
The Antarctic ◽  
Vertical Crustal Deformation

ABSTRACTThis study explores an iterative method for simultaneously estimating the present-day glacial isostatic adjustment (GIA), ice mass change and elastic vertical crustal deformation of the Antarctic ice sheet (AIS) for the period October 2003–October 2009. The estimations are derived by combining mass measurements of the GRACE mission and surface height observations of the ICESat mission under the constraint of GPS vertical crustal deformation rates in the spatial domain. The influence of active subglacial lakes on GIA estimates are mitigated for the first time through additional processing of ICESat data. The inferred GIA shows that the strongest uplift is found in the Amundsen Sea Embayment (ASE) sector and subsidence mostly occurs in Adelie Terre and the East Antarctica inland. The total GIA-related mass change estimates for the entire AIS, West Antarctica Ice Sheet (WAIS), East Antarctica Ice Sheet (EAIS), and Antarctic Peninsula Ice Sheet (APIS) are 43 ± 38, 53 ± 24, −23 ± 29 and 13 ± 6 Gt a−1, respectively. The overall ice mass change of the AIS is −46 ± 43 Gt a−1 (WAIS: −104 ± 25, EAIS: 77 ± 35, APIS: −20 ± 6). The most significant ice mass loss and most significant elastic vertical crustal deformations are concentrated in the ASE and northern Antarctic Peninsula.

scholarly journals Sensitivity of inverse glacial isostatic adjustment estimates over Antarctica

2020 ◽  
Vol 14 (1) ◽  
pp. 349-366
Author(s):  
Matthias O. Willen ◽  
Martin Horwath ◽  
Ludwig Schröder ◽  
Andreas Groh ◽  
Stefan R. M. Ligtenberg ◽  
...  
Keyword(s):  
Mass Balance ◽  
Bias Correction ◽  
Total Mass ◽  
Process Models ◽  
Mass Change ◽  
Quality Data ◽  
Data Sets ◽  
Glacial Isostatic Adjustment ◽  
Satellite Gravimetry ◽  
Isostatic Adjustment

Abstract. Glacial isostatic adjustment (GIA) is a major source of uncertainty for ice and ocean mass balance estimates derived from satellite gravimetry. In Antarctica the gravimetric effect of cryospheric mass change and GIA are of the same order of magnitude. Inverse estimates from geodetic observations hold some promise for mass signal separation. Here, we investigate the combination of satellite gravimetry and altimetry and demonstrate that the choice of input data sets and processing methods will influence the resultant GIA inverse estimate. This includes the combination that spans the full GRACE record (April 2002–August 2016). Additionally, we show the variations that arise from combining the actual time series of the differing data sets. Using the inferred trends, we assess the spread of GIA solutions owing to (1) the choice of different degree-1 and C20 products, (2) viable candidate surface-elevation-change products derived from different altimetry missions corresponding to different time intervals, and (3) the uncertainties associated with firn process models. Decomposing the total-mass signal into the ice mass and the GIA components is strongly dependent on properly correcting for an apparent bias in regions of small signal. Here our ab initio solutions force the mean GIA and GRACE trend over the low precipitation zone of East Antarctica to be zero. Without applying this bias correction, the overall spread of total-mass change and GIA-related mass change using differing degree-1 and C20 products is 68 and 72 Gt a−1, respectively, for the same time period (March 2003–October 2009). The bias correction method collapses this spread to 6 and 5 Gt a−1, respectively. We characterize the firn process model uncertainty empirically by analysing differences between two alternative surface mass balance products. The differences propagate to a 10 Gt a−1 spread in debiased GIA-related mass change estimates. The choice of the altimetry product poses the largest uncertainty on debiased mass change estimates. The spread of debiased GIA-related mass change amounts to 15 Gt a−1 for the period from March 2003 to October 2009. We found a spread of 49 Gt a−1 comparing results for the periods April 2002–August 2016 and July 2010–August 2016. Our findings point out limitations associated with data quality, data processing, and correction for apparent biases.

3D glacial-isostatic adjustment models using geodynamically constrained Earth structures

2021 ◽  
Author(s):  
Meike Bagge ◽  
Volker Klemann ◽  
Bernhard Steinberger ◽  
Milena Latinović ◽  
Maik Thomas
Keyword(s):  
Solid Earth ◽  
Sea Level ◽  
Ice Sheet ◽  
Sea Level Change ◽  
Glacial Isostatic Adjustment ◽  
Level Change ◽  
Isostatic Adjustment ◽  
Mantle Viscosity ◽  
The Earth ◽  
Earth Structures

<p>The interaction between ice sheets and the solid Earth plays an important role for ice-sheet stability and sea-level change and hence for global climate models. Glacial-isostatic adjustment (GIA) models enable simulation of the solid Earth response due to variations in ice-sheet and ocean loading and prediction of the relative sea-level change. Because the viscoelastic response of the solid Earth depends on both ice-sheet distribution and the Earth’s rheology, independent constraints for the Earth structure in GIA models are beneficial. Seismic tomography models facilitate insights into the Earth’s interior, revealing lateral variability of the mantle viscosity that allows studying its relevance in GIA modeling. Especially, in regions of low mantle viscosity, the predicted surface deformations generated with such 3D GIA models differ considerably from those generated by traditional GIA models with radially symmetric structures. But also, the conversion from seismic velocity variations to viscosity is affected by a set of uncertainties. Here, we apply geodynamically constrained 3D Earth structures. We analyze the impact of conversion parameters (reduction factor in Arrhenius law and radial viscosity profile) on relative sea-level predictions. Furthermore, we focus on exemplary low-viscosity regions like the Cascadian subduction zone and southern Patagonia, which coincide with significant ice-mass changes.</p>

scholarly journals Dynamic regimes of the Greenland Ice Sheet emerging from interacting melt-elevation and glacial isostatic adjustment feedbacks

2021 ◽  
Author(s):  
Maria Zeitz ◽  
Jan M. Haacker ◽  
Jonathan F. Donges ◽  
Torsten Albrecht ◽  
Ricarda Winkelmann
Keyword(s):  
Global Warming ◽  
Solid Earth ◽  
Ice Sheet ◽  
Greenland Ice Sheet ◽  
Lapse Rate ◽  
Glacial Isostatic Adjustment ◽  
Partial Recovery ◽  
Ice Volume ◽  
Isostatic Adjustment

Abstract. The stability of the Greenland Ice Sheet under global warming is governed by a number of dynamic processes and interacting feedback mechanisms in the ice sheet, atmosphere and solid Earth. Here we study the long-term effects due to the interplay of the competing melt-elevation and glacial isostatic adjustment (GIA) feedbacks for different temperature step forcing experiments with a coupled ice-sheet and solid-Earth model. Our model results show that for warming levels above 2 °C, Greenland could become essentially ice-free on the long-term, mainly as a result of surface melting and acceleration of ice flow. These ice losses can be mitigated, however, in some cases with strong GIA feedback even promoting the partial recovery of the Greenland ice volume. We further explore the full-factorial parameter space determining the relative strengths of the two feedbacks: Our findings suggest distinct dynamic regimes of the Greenland Ice Sheets on the route to destabilization under global warming – from recovery, via quasi-periodic oscillations in ice volume to ice-sheet collapse. In the recovery regime, the initial ice loss due to warming is essentially reversed within 50,000 years and the ice volume stabilizes at 61–93 % of the present-day volume. For certain combinations of temperature increase, atmospheric lapse rate and mantle viscosity, the interaction of the GIA feedback and the melt-elevation feedback leads to self-sustained, long-term oscillations in ice-sheet volume with oscillation periods of tens to hundreds of thousands of years and oscillation amplitudes between 15–70 % of present-day ice volume. This oscillatory regime reveals a possible mode of internal climatic variability in the Earth system on time scales on the order of 100,000 years that may be excited by or synchronized with orbital forcing or interact with glacial cycles and other slow modes of variability. Our findings are not meant as scenario-based near-term projections of ice losses but rather providing insight into of the feedback loops governing the "deep future" and, thus, long-term resilience of the Greenland Ice Sheet.

scholarly journals Sensitivity of inverse glacial isostatic adjustment estimates over Antarctica

2019 ◽  
Author(s):  
Matthias O. Willen ◽  
Martin Horwath ◽  
Ludwig Schröder ◽  
Andreas Groh ◽  
Stefan R. M. Ligtenberg ◽  
...  
Keyword(s):  
Mass Balance ◽  
Bias Correction ◽  
Total Mass ◽  
Mass Change ◽  
Data Sets ◽  
Glacial Isostatic Adjustment ◽  
Satellite Gravimetry ◽  
Isostatic Adjustment ◽  
Time Period ◽  
Combining Data

Abstract. Glacial isostatic adjustment (GIA) is a major source of uncertainty in estimated ice and ocean mass balance that are based on satellite gravimetry. In particular over Antarctica the gravimetric effect of cryospheric mass change and GIA are of the same order of magnitude. Inverse estimates from geodetic observations are promising for separating the two superimposed mass signals. Here, we investigate the combination of satellite gravimetry and altimetry and how the choice of input data sets and processing details affect the inverse GIA estimates. This includes the combination for almost full GRACE lifespan (2002-04/2016-08). Further we show results from combining data sets on time-series level. Specifically on trend level, we assess the spread of GIA solutions that arises from (1) the choice of different degree-1 and C20 products, (2) different surface elevation change products derived from different altimetry missions and associated to different time intervals, and (3) the uncertainty of firn-process models. The decomposition of the total-mass signal into the ice-mass signal and the apparent GIA-mass signal depends strongly on correcting for apparent biases in initial solutions by forcing the mean GIA and GRACE trend over the low precipitation zone of East Antarctica to be zero. Prior to bias correction, the overall spread of total-mass change and apparent GIA-mass change using differing degree-1 and C20 products is 68 and 72 Gt a−1, respectively, for the same time period (2003-03/2009-10). The bias correction suppresses this spread to 6 and 5 Gt a−1, respectively. We characterise the firn-process model uncertainty empirically by analysing differences between two alternative surface-mass-balance products. The differences propagate to a 21 Gt a−1 spread in apparent GIA-mass-change estimates. The choice of the altimetry product poses the largest uncertainty on debiased mass-change estimates. The overall spread of debiased GIA-mass change amounts to 18 and 49 Gt a−1 for a fixed time period (2003-03/2009-10) and various time periods, respectively. Our findings point out limitations associated with data processing, correction for apparent biases, and time dependency.

Antarctic glacial isostatic adjustment: a new assessment

2005 ◽  
Vol 17 (4) ◽  
pp. 541-553 ◽  
Author(s):  
ERIK R. IVINS ◽  
THOMAS S. JAMES
Keyword(s):  
Mass Balance ◽  
Antarctic Peninsula ◽  
Ross Sea ◽  
Ice Sheet ◽  
Gravity Change ◽  
Glacial Isostatic Adjustment ◽  
Crustal Motion ◽  
Isostatic Adjustment ◽  
Time Space ◽  
The Antarctic

The prediction of crustal motions and gravity change driven by glacial isostatic adjustment (GIA) in Antarctica is critically dependent on the reconstruction of the configuration and thickness of the ice sheet during the Late Pleistocene and Holocene. The collection and analysis of field data to improve the reconstruction has occurred at an accelerated pace during the past decade. At the same time, space-based imaging and altimetry, combined with on-ice velocity measurements using Global Positioning System (GPS) geodesy, has provided better assessments of the present-day mass balance of the Antarctic ice sheet. Present-day mass change appears to be dominated by deglaciation that is, in large part, a continuation of late-Holocene evolution. Here a new ice load model is constructed, based on a synthesis of the current constraints on past ice history and present-day mass balance. The load is used to predict GIA crustal motion and geoid change. Compared to existing glacioisostatic models, the new ice history model is significantly improved in four aspects: (i) the timing of volume losses in the region ranging from the Ross Sea sector to the Antarctic Peninsula, (ii) the maximum ice heights in parts of the Ellsworth and Transantarctic Mountains, (iii) maximum grounding line position in Pine Island Bay, the Antarctic Peninsula, and in the Ross Sea, (iv) incorporation of present-day net mass balance estimates. The predicted present-day GIA uplift rates peak at 14–18 mm yr−1 and geoid rates peak at 4–5 mm yr−1 for two contrasting viscosity models. If the asthenosphere underlying West Antarctica has a low viscosity then the predictions could change substantially due to the extreme sensitivity to recent (past two millennia) ice mass variability. Future observations of crustal motion and gravity change will substantially improve the understanding of sub-Antarctic lithospheric and mantle rheology.

scholarly journals A Joint Inversion Estimate of Antarctic Ice Sheet Mass Balance Using Multi-Geodetic Data Sets

2019 ◽  
Vol 11 (6) ◽  
pp. 653 ◽  
Author(s):  
Chunchun Gao ◽  
Yang Lu ◽  
Zizhan Zhang ◽  
Hongling Shi
Keyword(s):  
Mass Balance ◽  
Joint Inversion ◽  
Ice Sheet ◽  
Mass Change ◽  
West Antarctica ◽  
Amundsen Sea ◽  
Antarctic Ice Sheet ◽  
Forward Models ◽  
Uplift Rates ◽  
The Antarctic

Many recent mass balance estimates using the Gravity Recovery and Climate Experiment (GRACE) and satellite altimetry (including two kinds of sensors of radar and laser) show that the ice mass of the Antarctic ice sheet (AIS) is in overall decline. However, there are still large differences among previously published estimates of the total mass change, even in the same observed periods. The considerable error sources mainly arise from the forward models (e.g., glacial isostatic adjustment [GIA] and firn compaction) that may be uncertain but indispensable to simulate some processes not directly measured or obtained by these observations. To minimize the use of these forward models, we estimate the mass change of ice sheet and present-day GIA using multi-geodetic observations, including GRACE and Ice, Cloud and land Elevation Satellite (ICESat), as well as Global Positioning System (GPS), by an improved method of joint inversion estimate (JIE), which enables us to solve simultaneously for the Antarctic GIA and ice mass trends. The GIA uplift rates generated from our JIE method show a good agreement with the elastic-corrected GPS uplift rates, and the total GIA-induced mass change estimate for the AIS is 54 ± 27 Gt/yr, which is in line with many recent GPS calibrated GIA estimates. Our GIA result displays the presence of significant uplift rates in the Amundsen Sea Embayment of West Antarctica, where strong uplift has been observed by GPS. Over the period February 2003 to October 2009, the entire AIS changed in mass by −84 ± 31 Gt/yr (West Antarctica: −69 ± 24, East Antarctica: 12 ± 16 and the Antarctic Peninsula: −27 ± 8), greater than the GRACE-only estimates obtained from three Mascon solutions (CSR: −50 ± 30, JPL: −71 ± 30, and GSFC: −51 ± 33 Gt/yr) for the same period. This may imply that single GRACE data tend to underestimate ice mass loss due to the signal leakage and attenuation errors of ice discharge are often worse than that of surface mass balance over the AIS.

Observations of postglacial uplift at Churchill, Manitoba

1992 ◽  
Vol 29 (11) ◽  
pp. 2418-2425 ◽  
Author(s):  
A. Mark Tushingham
Keyword(s):  
Sea Level ◽  
Tide Gauge ◽  
Ice Sheet ◽  
Laurentide Ice Sheet ◽  
Glacial Isostatic Adjustment ◽  
Tide Gauge Record ◽  
Isostatic Adjustment ◽  
Absolute Gravimetry ◽  
The Earth ◽  
Level Data

Churchill, Manitoba, is located near the centre of postglacial uplift caused by the Earth's recovery from the melting of the Laurentide Ice Sheet. The value of present-day uplift at Churchill has important implications in the study of postglacial uplift in that it can aid in constraining the thickness of the ice sheet and the rheology of the Earth. The tide-gauge record at Churchill since 1940 is examined, along with nearby Holocene relative sea-level data, geodetic measurements, and recent absolute gravimetry measurements, and a present-day rate of uplift of 8–9 mm/a is estimated. Glacial isostatic adjustment models yield similar estimates for the rate of uplift at Churchill. The effects of the tide-gauge record of the diversion of the Churchill River during the mid-1970's are discussed.

scholarly journals Comparing a thermo-mechanical Weichselian ice sheet reconstruction to GIA driven reconstructions: aspects of earth response and ice configuration

2013 ◽  
Vol 5 (2) ◽  
pp. 2345-2388 ◽  
Author(s):  
P. Schmidt ◽  
B. Lund ◽  
J-O. Näslund
Keyword(s):  
Sea Level ◽  
Ice Sheet ◽  
Younger Dryas ◽  
Slow Phase ◽  
Glacial Isostatic Adjustment ◽  
Relative Sea Level ◽  
Isostatic Adjustment ◽  
Earth Models ◽  
Uplift Rates ◽  
Best Fit

Abstract. In this study we compare a recent reconstruction of the Weichselian ice-sheet as simulated by the University of Main ice-sheet model (UMISM) to two reconstructions commonly used in glacial isostatic adjustment (GIA) modeling: ICE-5G and ANU (also known as RSES). The UMISM reconstruction is carried out on a regional scale based on thermo-mechanical modelling whereas ANU and ICE-5G are global models based on the sea-level equation. The Weichselian ice-sheet in the three models are compared directly in terms of ice volume, extent and thickness, as well as in terms of predicted glacial isostatic adjustment in Fennoscandia. The three reconstructions display significant differences. UMISM and ANU includes phases of pronounced advance and retreat prior to the last glacial maximum (LGM), whereas the thickness and areal extent of the ICE-5G ice-sheet is more or less constant up until LGM. The final retreat of the ice-sheet initiates at earliest time in ICE-5G and latest in UMISM, while ice free conditions are reached earliest in UMISM and latest in ICE-5G. The post-LGM deglaciation style also differs notably between the ice models. While the UMISM simulation includes two temporary halts in the deglaciation, the later during the Younger Dryas, ANU only includes a decreased deglaciation rate during Younger Dryas and ICE-5G retreats at a relatively constant pace after an initial slow phase. Moreover, ANU and ICE-5G melt relatively uniformly over the entire ice-sheet in contrast to UMISM which melts preferentially from the edges. We find that all three reconstructions fit the present day uplift rates over Fennoscandia and the observed relative sea-level curve along the Ångerman river equally well, albeit with different optimal earth model parameters. Given identical earth models, ICE-5G predicts the fastest present day uplift rates and ANU the slowest, ANU also prefers the thinnest lithosphere. Moreover, only for ANU can a unique best fit model be determined. For UMISM and ICE-5G there is a range of earth models that can reproduce the present day uplift rates equally well. This is understood from the higher present day uplift rates predicted by ICE-5G and UMISM, which results in a bifurcation in the best fit mantle viscosity. Comparison of the uplift histories predicted by the ice-sheets indicate that inclusion of relative sea-level data in the data fit can reduce the observed ambiguity. We study the areal distributions of present day residual surface velocities in Fennoscandia and show that all three reconstructions generally over-predict velocities in southwestern Fennoscandia and that there are large differences in the fit to the observational data in Finland and northernmost Sweden and Norway. These difference may provide input to further enhancements of the ice-sheet reconstructions.

Sign in / Sign up

Export Citation Format

Share Document