No abstract

Accurate GIA models are required for correcting measurements of mass change in Antarctica and for improving knowledge of the sub-surface, especially in areas of large current ice loss such as the Amundsen Sea Embayment(ASE). Regionally, seismic and gravity data suggests lateral differences in viscosity (3D). Furthermore, mantle flow laws allow for a stress-dependent effective viscosity which changes over time (3D-s). In this study we investigate whether models with 3D/3D-s have significant effects on the uplift in the region. We use a finite element model with composite rheology consisting of diffusion and dislocation creep, forced by an ice deglaciation model starting in 1900. We use its uplift predictions as synthetic observations to test the performance of 1D model inversion in the presence of viscosity variations. Stress-dependent rheology results in lower viscosity beneath the load and a more localized uplift pattern. We demonstrate that the background stress from earlier ice load changes can both increase or decrease the influence of stress-induced effective viscosity changes. For the ASE, fitting 1D models to 3D model uplift results in a best fitting model with viscosity that represents the average of a large area, while for 3D-s rheology, local viscosity is more influential. 1D models are statistically indistinguishable from 3D/3D-s viscosity with current GPS stations. However, 3D and 3D-

Abstract. We present a finite-element model of post-seismic solid Earth deformation built in the software package Abaqus (version 2018). The model is global and spherical, includes self-gravitation and is built for the purpose of calculating post-seismic deformation in the far field (>∼300 km) of major earthquakes. An earthquake is simulated by prescribing slip on a fault plane in the mesh and the model relaxes under the resulting change in stress. Both linear Maxwell and biviscous (Burgers) rheological models have been implemented and the model can be easily adapted to include different rheological models and lateral variations in Earth structure, a particular advantage over existing models. We benchmark the model against an analytical coseismic solution and an existing open-source post-seismic model code, demonstrating good agreement for all fault geometries tested. Due to the inclusion of self-gravity, the model has the potential for predicting deformation in response to multiple sources of stress change, for example, changing ice thickness in tectonically active regions.

Abstract. We present a finite-element model of postseismic solid Earth deformation built in the software package ABAQUS (version 2018). The model is global and spherical, and includes self-gravitation and is built for the purpose of calculating postseismic deformation in the far-field (> ~ 300 km) of major earthquakes. An earthquake is simulated by prescribing slip on a fault plane in the mesh and the model relaxes under the resulting change in stress. Both linear Maxwell and biviscous (Burgers) rheological models have been implemented and the model can be easily adapted to include different rheological models and lateral variations in Earth structure, a particular advantage over existing models. We benchmark the model against an analytical coseismic solution and an existing open-source postseismic model, demonstrating good agreement for all fault geometries tested. Due to the inclusion of self-gravity the model has the potential for predicting deformation in response to multiple sources of stress change, for example, changing ice thickness in tectonically active regions.

No abstract