Glacial isostatic adjustment of the Pacific Coast of North America: the influence of lateral Earth structure

Abstract

SUMMARY The Pacific Coast of Central North America is a geodynamically complex region which has been subject to various geophysical processes operating on different timescales. Glacial isostatic adjustment (GIA), the ongoing deformational response of the solid Earth to past deglaciation, is an important geodynamic process in this region. In this study, we apply earth models with 3-D structure to determine if the inclusion of lateral structure can explain the poor performance of 1-D models in this region. Three different approaches are used to construct 3-D models of the Earth structure. For the first approach, we adopt an optimal 1-D viscosity structure from previous work and add lateral variations based on four global seismic shear wave velocity anomalies and two global lithosphere thickness models. The results based on these models indicate that the addition of lateral structure significantly impacts modelled RSL changes, but the data-model fits are not improved. The global seismic models are limited in spatial resolution and so two other approaches were considered to produce higher resolution models of 3-D structure: inserting a regional seismic model into two of the global seismic models and, explicitly incorporating regional structure of the Cascadia subduction zone and vicinity, that is the subducting slab, the overlying mantle wedge and the plate boundary interface. The results associated with these higher resolution models do not reveal any clear improvement in satisfying the RSL observations, suggesting that our estimates of lateral structure are inaccurate and/or the data-model misfits are primarily due to limitations in the adopted ice-loading history. The different realizations of 3-D Earth structure gives useful insight to uncertainty associated with this aspect of the GIA model. Our results indicate that improving constraints on the deglacial history of the southwest sector of the Cordilleran ice sheet is an important step towards developing more accurate of GIA models for this region.