Glacial Isostatic Adjustment with 3D Earth models: A comparison of case studies of deglacial relative sea level records of North America and Russian Arctic

The Canadian landmass of North America and the Russian Arctic were covered by large ice sheets during the Last Glacial Maximum, and have been key areas for Glacial Isostatic Adjustment (GIA) studies. Previous GIA studies have applied 1D models of Earth’s interior viscoelastic structure; however, seismic tomography, field geology and recent studies reveal the potential importance of 3D models of this structure. Here, using the latest quality-controlled deglacial sea-level databases from North America and the Russian Arctic, we investigate the effects of 3D structure on GIA predictions. We explore scaling factors in the upper mantle (β_UM) and lower mantle (β_LM) and the 1D background viscosity model (η_o) with predictions of of the ICE-6G_C (VM5a) glaciation/deglaciation model of Peltier et al (2015, JGR) in these two regions, and compare with the best fit 3D viscosity structures.

We compute gravitationally self-consistent relative sea-level histories with time dependent coastlines and rotational feedback using both the Normal Mode Method and Coupled Laplace-Finite Element Method. A subset of 3D GIA models is found that can fit the deglacial sea-level databases for both regions. These databases cover both the near and intermediate field regions. However, North America and Russian Arctic prefer different 3D structures (i.e., combinations of (η_o, β_UM, β_LM)) to provide the best fits. The Russian Arctic database prefers a softer background viscosity model (η_o), but larger scaling factors (β_UM, β_LM) than those preferred by the North America database.

Outstanding issues include the uncertainty of the history of local glaciation history. For example, preliminary modifications of the ice model in Russian Arctic reveal that the misfits of 1D models can be significantly reduced, but still fit less well than the best fit 3D GIA model.An additional issue concerns the extent to which the 3D models are able to improve both fits in North America and Russian Arctic when compared with 1D internal structure (ICE-6G_C VM5a & ICE-7G VM7), will be assessed in a preliminary fashion.