ALBERT

Description: Glacial isostatic adjustment is dominated by Earth rheology resulting in a variability of relative sea-level (RSL) predictions of more than 100 meters during the last glacial cycle. Seismic tomography models reveal significant lateral variations in seismic wavespeed, most likely corresponding to variations in temperature and hence viscosity. Therefore, the replacement of 1D Earth structures by a 3D Earth structure is an essential part of recent research to reveal the impact of lateral viscosity contrasts and to achieve a more consistent view on solid-Earth dynamics. Here, we apply the VIscoelastic Lithosphere and MAntle model VILMA to predict RSL during the last deglaciation. We create an ensemble of geodynamically constrained 3D Earth structures which is based on seismic tomography models while considering a range of conversion factors to transfer seismic velocity variations into viscosity variations. For a number of globally distributed sites, we discuss the resulting variability in RSL predictions, compare this with regionally optimized 1D Earth structures, and validate the model results with relative sea-level data (sea-level indicators). This study is part of the German Climate Modeling initiative PalMod aiming the modeling of the last glacial cycle under consideration of a coupled Earth system model, i.e. including feedbacks between ice-sheets and the solid Earth.

Description: The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3-D lithosphere–asthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 (WSM2016) and the observation-based residual topography. We derive the upper mantle thermal structure from either a heat flow model combined with a seafloor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. The modeled stress field directions, using only the mantle heterogeneities below 300 km, are not perturbed much when the effects of lithosphere and crust above 300 km are added. In contrast, modeled stress magnitudes and dynamic topography are to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a good correlation of 0.51 in continents, but this correction leads to no significant improvement of the fit between the WSM2016 and the resulting lithosphere stresses. In continental regions with abundant heat flow data, TM1 results in relatively small angular misfits. For example, in western Europe the misfit between the modeled and observation-based stress is 18.3°. Our findings emphasize that the relative contributions coming from shallow and deep mantle dynamic forces are quite different for the lithospheric stress field and dynamic topography.

Description: 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.

Description: 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.

Description: Ice-penetrating radar1,2,3 and ice core drilling4 have shown that large parts of the north-central Greenland ice sheet are melting from below. It has been argued that basal ice melt is due to the anomalously high geothermal flux1,4 that has also influenced the development of the longest ice stream in Greenland1. Here we estimate the geothermal flux beneath the Greenland ice sheet and identify a 1,200-km-long and 400-km-wide geothermal anomaly beneath the thick ice cover. We suggest that this anomaly explains the observed melting of the ice sheet’s base, which drives the vigorous subglacial hydrology3 and controls the position of the head of the enigmatic 750-km-long northeastern Greenland ice stream5. Our combined analysis of independent seismic, gravity and tectonic data6,7,8,9 implies that the geothermal anomaly, which crosses Greenland from west to east, was formed by Greenland’s passage over the Iceland mantle plume between roughly 80 and 35 million years ago. We conclude that the complexity of the present-day subglacial hydrology and dynamic features of the north-central Greenland ice sheet originated in tectonic events that pre-date the onset of glaciation in Greenland by many tens of millions of years.

Description: Glacial-isostatic adjustment (GIA) is the key process controlling relative sea-level (RSL) and paleo-topography. The viscoelastic response of the solid Earth is controlled by its viscosity structure. Therefore, the appropriate choice of Earth structure for GIA models is still an important area of research in geodynamics. We construct 18 3D Earth structures that are derived from seismic tomography models and are geodynamically constrained. We consider uncertainties in 3D viscosity structures that arise from variations in the conversion from seismic velocity to temperature variations (factor r) and radial viscosity profiles (RVP). We apply these Earth models to a 3D GIA model, VILMA, to investigate the influence of such structure on RSL predictions. The variabilities in 3D Earth structures and RSL predictions are investigated for globally distributed sites and applied for comparisons with regional 1D models for ice center (North America, Antarctica) and peripheral regions (Central Oregon Coast, San Jorge Gulf). The results from 1D and 3D models reveal substantial influence of lateral viscosity variations on RSL. Depending on time and location, the influence of factor r and/or RVP can be reverse, for example, the same RVP causes lowest RSL in Churchill and largest RSL in Oregon. Regional 1D models representing the structure beneath the ice and 3D models show similar influence of factor r and RVP on RSL prediction. This is not the case for regional 1D models representing the structure beneath peripheral regions indicating the dependence on the 3D Earth structure. The 3D Earth structures of this study are made available.