Influence of 3D Earth Structure on Glacial Isostatic Adjustment in the Russian Arctic

Analyses of glacial isostatic adjustment (GIA) and deglacial relative sea-level (RSL) change in the Russian Arctic deliver important insights into the Earth's viscosity structure and the deglaciation history of the Eurasian ice sheet complex. Here, we validate the 1D GIA models ICE-6G_C (VM5a) and ICE-7G_NA (VM7) and select new 3D GIA models in the Russian Arctic against a quality-controlled deglacial RSL database of >500 sea-level data points from 24 regions. Both 1D models correspond to the RSL data along the southern coast of the Barents Sea and Franz Josef Land from ∼11 ka BP to present but show notable misfits (>50 m at 10 ka BP) with the White Sea data. We find 3D model predictions of deglacial RSL resolve most of the misfits with the observed data for the White Sea while retaining comparable fits in other regions of the Russian Arctic. Our results further reveal: (a) RSL in the western Russian Arctic is sensitive to elastic lithosphere with lateral thickness variation and 3D viscosity structure in the upper mantle; and (b) RSL in the whole Russian Arctic is less sensitive to 3D viscosity structure in the lower mantle compared to the upper mantle. The 3D models reveal a compromise in the upper mantle between the background viscosity and scaling factor to best fit the RSL data, which needs to be considered in future 3D GIA studies.