scholarly journals Glacial‐isostatic adjustment models using geodynamically constrained 3D Earth structures

Author(s):  
M. Bagge ◽  
V. Klemann ◽  
B. Steinberger ◽  
M. Latinović ◽  
M. Thomas
2021 ◽  
Author(s):  
Meike Bagge ◽  
Volker Klemann ◽  
Bernhard Steinberger ◽  
Milena Latinović ◽  
Maik Thomas

<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>


2020 ◽  
Author(s):  
Robert Hartmann ◽  
Jörg Ebbing ◽  
Clinton P. Conrad

<p>The pseudo-spectral form of the sea level equation (SLE) requires the approximation of a radially-symmetric visco-elastic Earth. Thus, the resulting predictions of sea level change (SLC) and glacial isostatic adjustment (GIA) often ignore lateral variations in the Earth structure. Here, we assess the capabilities of a Multiple 1D Earth Approach (M1DEA) applied to large-scale ice load components with different Earth structures to account for these variations. In this approach the total SLC and GIA responses result from the superposition of individual responses from each load component, each computed globally assuming locally-appropriate 1D Earth structures. We apply the M1DEA to three separate regions (East Antarctica, West Antarctica, and outside Antarctica) to analyze uplift rates for a range of Earth structures and different ice loads at various distances. We find that the uplift response is mostly sensitive to the local Earth structure, which supports the usefulness of the M1DEA. However, stresses transmitted across rheological boundaries (e.g., producing peripheral bulges) present challenges for the M1DEA, but can be minimized under two conditions: (1) If the considered time period of ice loading for each component is consistent with the relaxation time of the local Earth structure. (2) If the load components can be subdivided according to the scale of the lateral variations in Earth structure. Overall, our results indicate that M1DEA could be a computationally much cheaper alternative to 3D finite element models, but further work is needed to quantify the relative accuracy of both methods for different resolutions, loads, and Earth structure variations.</p>


1992 ◽  
Vol 29 (11) ◽  
pp. 2418-2425 ◽  
Author(s):  
A. Mark Tushingham

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.


2020 ◽  
Author(s):  
Tanghua Li ◽  
Nicole Khan ◽  
Alisa Baranskaya ◽  
Timothy Shaw ◽  
W Richard Peltier ◽  
...  

2013 ◽  
Vol 5 (2) ◽  
pp. 2345-2388 ◽  
Author(s):  
P. Schmidt ◽  
B. Lund ◽  
J-O. Näslund

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