A 3D Finite Element Model with Improved Spatial Resolution to Investigate the Effect of Varying Viscosity on Antarctica

B. Blank1, H. Hu1, W. van der Wal1, F Pappa2, J. Ebbing21Delft University of Technology 2Christian-Albrechts-University of KielSince the beginning of the 2000's time-variable gravity data from GRACE has proved to be an effective method for mapping ice mass loss in Antarctica. However, Glacial Isostatic Adjustment (GIA) models are required to correct for GIA induced mass changes. While most GIA models have adopted an Earth model that only varies radially in parameters, it has long been clear that the Earth structure also varies with longitude and latitude. For this study a new global 3D GIA model has been developed within the finite element software package ABAQUS, which can be modified to operate on a spatial resolution down to 50 km locally. The model is being benchmarked against normal model models for surface loading. It will be used to investigate the effects of a 3D varying lithosphere and upper asthenosphere in Antarctica. Viscosity which will be computed from temperature estimates with laboratory based flow laws. A new 3D temperature map of the Antarctic lithosphere has been developed within ESA's GOCE+ project based on seismic data as well as on GOCE and GRACE inferred gravity gradients. Output from the GIA model with this new temperature estimates will be compared to that of 1D viscosity profiles and other recent 3D viscosity models based on seismic data. From these side to side comparisons we want to investigate the influence of the viscosity map on uplift rates and horizontal movement. Finally the results can be compared to GPS measurement to investigate the validity of all models.