Present‐day ice‐mass changes in Antarctica typically deform the solid Earth elastically, and this signal needs removing from Global Positioning System (GPS) observations of displacement before they can be used to constrain models of glacial isostatic adjustment. However, much of West Antarctica is underlain by weak upper mantle, meaning these short‐term fluctuations may also cause viscous (transient or steady‐state) deformation of the Earth. We model the viscoelastic response of the Earth to surface mass balance (SMB) variability in the Antarctic Peninsula and find an improved fit to GPS data at most sites compared to elastic only. Viscoelastic modeling constrains upper mantle steady‐state viscosity in the northern Peninsula to 5 × 10 17 to 2 × 10 18 Pa s, and >1 × 10 18 Pa s for the mid to southern Peninsula. In the northern Peninsula, removing viscoelastic displacement caused by SMB variability from GPS time series increases estimated uplift rates by up to 3 mm/yr compared with using an elastic‐only correction.

Raised shorelines and associated lacustrine sediments in the central Estrecho de Magallanes (Strait of Magellan) have been interpreted as products of cordilleran glaciers impounding a large proglacial lake and preventing drainage to the South Pacific and Southern Ocean during the Late glacial between c. 15.0 and 12.0 cal ka BP. However, a growing body of glacial geological evidence points towards an earlier retreat of the Magellan cordilleran ice dome, insufficient to dam lakes at that time. We critically re-evaluate the extant evidence for the c. 15.0-12.0 cal ka BP lake, here named 'Lago Kawesqar', and provide further sedimentological and chronological evidence for its existence. We also provide new cosmogenic surface nuclide dating of erratic and bedrock samples collected from extensive field campaigns that confirm the rapid and widespread retreat of the Magellan ice fields to the inner fjords of the Fuegian archipelago by c. 16.0 ka. To resolve the apparent paradox between these two lines of evidence we propose that glacial isostatic adjustment led to a topographic barrier to lake drainage rather than an ice dam. We use Glacial Isostatic Adjustment modelling to demonstrate that rapid isostatic recovery following the early deglaciation after c. 17.0 cal ka BP likely led to elevation of the present shallow south-western coastal margin of the Fuegian archipelago. Final drainage of Lago Kawesqar was probably caused by neotectonic subsidence of the same margin along the boundary of the South American - Scotia tectonic plates at c. 12.0 cal ka BP.

ANET-POLENET (Antarctic Network of the Polar Earth Observing Network) bedrock GNSS sites in the Ross Sea region of Antarctica surround an LGM load center in the Siple region of the Ross Embayment and record crustal motion due to GIA. Rather than a radial pattern of horizontal motion away from the former load, we instead observe three primary patterns of deformation; 1) motions are reversed towards the load in the southern region of the Transantarctic Mountains (TAM), 2) motions are radially away from the load in the Marie Byrd Land (MBL) region, and 3) an overall gradient in motion is present, with magnitudes progressively increasing from East to West Antarctica. We investigate the effects of alternative Earth model and ice loading scenarios, with the goal of understanding these distinct patterns of horizontal bedrock motion and their drivers. Using GIA models with a range of 1D Earth models, alternative ice loading scenarios for the Wilkes Subglacial Basin (LGM time scale) and the Siple Coast (centennial and millennial time scales) are explored. We find that no 1D model, regardless of the Earth model and ice loading scenario used, reproduces all three distinct patterns of observed motion at the same time. For select ice loading scenarios we also examine the influence of more complex rheology by invoking a boundary in Earth properties beneath the Transantarctic Mountains. This approach accounts for the strong lateral gradient in Earth properties across the continent by effectively separating East and West Antarctica into two different Earth model profiles. Some of our GIA models utilizing 3D Earth structure reproduce predicted motions that match all three observed patterns of deformation, and we find that a multiple order magnitude of change in upper mantle viscosity between East and West Antarctica is required to fit the observations.

Large earthquakes in the vicinity of Antarctica have the potential to cause postseismic viscoelastic deformation affecting measurements of displacement that are used to constrain models of glacial isostatic adjustment (GIA). In November 2013, a Mw 7.7 strike‐slip earthquake occurred in the Scotia Sea, 650 km from the Antarctic Peninsula. GPS time series from the northern Peninsula show a change in rate after this event, indicating a far‐field postseismic deformation signal is present. In this study, we use a finite element model with a suite of 1D and 3D Earth structures to investigate the extent of postseismic deformation in the Antarctic Peninsula. Model output is compared with GPS time series to place constraints on the Earth structure in this region. The preferred Earth structure has a thin lithosphere combined with a Burgers rheology with steady‐state viscosity of 4 × 10 18 Pa s and transient viscosity one order of magnitude lower. Our study shows that including 3D Earth structure does not improve the fit. Using the best fitting Earth structure, we run a forward model of the nearby 2003 Mw 7.6 strike‐slip earthquake and combine the predictions for both earthquakes. We show that postseismic deformation is widespread across the northern Peninsula with rates of horizontal deformation up to 1.65 mm/yr for the period 2015–2020, a signal that persists for decades. These results suggest that much of Antarctica may be deforming due to recent postseismic deformation and this signal needs to be accounted for when using GPS observations to constrain geophysical models.

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 This chapter reviews glacial isostatic adjustment (GIA) and post-seismic deformation in Antarctica. It discusses numerical models and their inputs, and observations and inferences that have been made from them. Both processes are controlled by mantle viscosity but their forcings are different. Ongoing GIA induced by the loss of ice since the last glacial maximum (LGM) could have amounted to 5–15 m of global sea-level rise. However, mantle viscosity is so low in parts of West Antarctica ( c. 10 18 Pa s) that changes in ice thickness over the last centuries and decades have controlled the current uplift rates there. The uplift due to GIA has promoted ice-sheet stability since the LGM, and in West Antarctica GIA is a significant negative feedback on the current decline of the ice sheet. Post-seismic deformation following the 1998 earthquake near the Balleny Islands south of New Zealand has been detected in global navigation satellite system (GNSS) data and compared to model outputs. The best-fitting viscosity for this area is c. 10 19 Pa s, similar to GIA-based estimates for the Antarctic Peninsula. Future work should focus on unifying descriptions of viscosity across geodynamic models, and integrating information from seismic, gravity, experimental and geological data.

In Antarctica, Global Positioning System (GPS) vertical time series exhibit non‐linear signals over a wide range of temporal scales. To explain these non‐linearities, a number of hypotheses have been proposed, among them the short‐term rapid solid Earth response to contemporaneous ice mass change. Here we use GPS vertical time series to reveal the solid Earth response to variations in surface mass balance (SMB) in the Southern Antarctic Peninsula (SAP). At four locations in the SAP we show that interannual variations of SMB anomalies cause measurable elastic deformation. We use regional climate model SMB products to calculate the induced displacement assuming a perfectly elastic Earth. Our results show a reduction of the misfit when fitting a linear trend to GPS time series corrected for the elastic response to SMB variations. Our results imply that, for a better understanding of the glacial isostatic adjustment signal in Antarctica, SMB variability must be considered.

Large earthquakes in the vicinity of Antarctica have the potential to cause post-seismic deformation on the continent, affecting measurements of displacement and gravity field change from GRACE or those attempting to constrain models of glacial isostatic adjustment. In November 2013 a magnitude 7.7 strike-slip earthquake occurred in the Scotia Sea around 650 km from the northern tip of the Antarctic Peninsula. GPS coordinate time series from the Peninsula region show a change in rate after this event indicating a far-field post-seismic deformation signal is present. At these far-field locations, the effects of fault after-slip are likely negligible and hence we consider the deformation to be due to post-seismic viscoelastic deformation. Here we use a global spherical finite element model to investigate the extent of post-seismic viscoelastic deformation in the northern Antarctic Peninsula. We investigate possible 1D earth models that can fit the GPS data and consider the effect of including a simple 3D earth structure in the region. These results, combined with previous results showing East Antarctica is still deforming following 1998 Mw 8.2 intraplate earthquake, suggest that much of Antarctica is deforming due to recent post-seismic deformation.

We present a new, open-source viscoelastic solid earth deformation model, Elmer/Earth. Using the multi-physics finite-element package Elmer, a model to compute viscoelastic material deformation has been implemented into the existing linear elasticity solver routine. Unlike approaches often implemented in engineering codes, our solver accounts for the restoring force of buoyancy within a system of layers with depth-varying density. It does this by directly integrating the solution of the system rather than by applying stress-jump conditions in the form of Winkler foundations on inter-layer boundaries, as is usually needed when solving the minimization problem given by the stress divergence in commercial codes. We benchmarked the new model with results from a commercial finite-element engineering package (ABAQUS, v2018) and another open-source code that uses viscoelastic normal mode theory, TABOO, using a flat-earth setup loaded by a cylindrical disc of 100 km in diameter and 100 m in height at the density of ice. Evaluating the differences in predicted surface deformation at the centre of the load and two distinctive distances (100 and 200 km), average deviations of 7 and 2.7 cm of Elmer/Earth results to ABAQUS and TABOO, respectively, were observed. In view of more than 100 cm maximum vertical deformation and the different numerical methods and parameters, these are very encouraging results. Elmer is set up as a highly scalable parallel code and distributed under the (L)GPL license, meaning that large-scale computations can be made without any licensing restrictions. Scaling figures presented in this paper show good parallel performance of the new model. Additionally, the high-fidelity ice-sheet code Elmer/Ice utilizes the same source base as Elmer and thereby the new model opens the way to undertaking high-resolution coupled ice-flow–solid-earth deformation simulations, which are required for robust projections of future sea-level rise and glacial isostatic adjustment.

Using 8- to 25-s-period Rayleigh and Love wave phase velocity dispersion data extracted from seismic ambient noise, we (i) model the 3-D shear wave velocity structure of the West Antarctic crust and (ii) map variations in crustal radial anisotropy. Enhanced regional resolution is offered by the UK Antarctic Seismic Network. In the West Antarctic Rift System (WARS), a ridge of crust similar to 26-30 km thick extending south from Marie Byrd Land separates domains of more extended crust (similar to 22 km thick) in the Ross and Amundsen Sea Embayments, suggesting along-strike variability in the Cenozoic evolution of the WARS. The southern margin of the WARS is defined along the southern Transantarctic Mountains and Haag-Ellsworth Whitmore Mountains (HEW) block by a sharp crustal thickness gradient. Crust similar to 35-40 km is modeled beneath the Haag Nunataks-Ellsworth Mountains, decreasing to similar to 30-32 km thick beneath the Whitmore Mountains, reflecting distinct structural domains within the composite HEW block. Our analysis suggests that the lower crust and potentially the middle crust is positively radially anisotropic (V-SH>V-SV) across West Antarctica. The strongest anisotropic signature is observed in the HEW block, emphasizing its unique provenance among West Antarctica's crustal units, and conceivably reflects a similar to 13-km-thick metasedimentary succession atop Precambrian metamorphic basement. Positive radial anisotropy in the WARS crust is consistent with observations in extensional settings and likely reflects the lattice-preferred orientation of minerals such as mica and amphibole by extensional deformation. Our observations support a contention that anisotropy may be ubiquitous in the continental crust.

We present a shear wave model of the West Antarctic upper mantle to similar to 200 km depth with enhanced regional resolution from the 2016-2018 UK Antarctic Seismic Network. The model is constructed from the combination of fundamental mode Rayleigh wave phase velocities extracted from ambient noise (periods 8-25 s) and earthquake data by two-plane wave analysis (periods 20-143 s). We seek to (i) image and interpret structures against the tectonic evolution of West Antarctica, and (ii) extract information from the seismic model that can serve as boundary conditions in ice sheet and glacial isostatic adjustment modelling efforts. The distribution of low velocity anomalies in the uppermost mantle suggests that recent tectonism in the West Antarctic Rift System (WARS) is mainly concentrated beneath the rift margins and largely confined to the uppermost mantle (<180 km). On the northern margin of the WARS, a pronounced low velocity anomaly extends eastward from beneath the Marie Byrd Land dome toward Pine Island Bay, underlying Thwaites Glacier, but not Pine Island Glacier. If of plume-related thermal origin, the velocity contrast of similar to 5% between this anomaly and the inner WARS translates to a temperature difference of similar to 125-200 degrees C. However, the strike of the anomaly parallels the paleo-Pacific convergent margin of Gondwana, so it may reflect subduction-related melt and volatiles rather than anomalously elevated temperatures, or a combination thereof. Motivated by xenolith analyses, we speculate that high velocity zones imaged south of the Marie Byrd Land dome and in the eastern Ross Sea Embayment might reflect the compositional signature of ancient continental fragments. A pronounced low velocity anomaly underlying the southern Transantarctic Mountains (TAM) is consistent with a published lithospheric foundering hypothesis. Taken together with a magnetotelluric study advocating flexural support of the central TAM by thick, stable lithosphere, this points to along-strike variation in the tectonic history of the TAM. A high velocity anomaly located in the southern Weddell Sea Rift System might reflect depleted mantle lithosphere following the extraction of voluminous melt related to Gondwana fragmentation. Lithospheric thickness estimates extracted from 1D shear wave velocity profiles representative of tectonic domains in West Antarctica indicate an average lithospheric thickness of similar to 85 km for the WARS, Marie Byrd Land, and Thurston Island block. This increases to similar to 96km in the Ellsworth Mountains. A surface heat flow of similar to 60 mW/m(2) and attendant geotherm best explains lithospheric mantle shear wave velocities in the central WARS and in the Thurston Island block adjacent to Pine Island Glacier; a similar to 50 mW/m(2) geotherm best explains the velocities in the Ellsworth Mountains, and a similar to 60 mW/m(2) geotherm best explains a less well-constrained velocity profile on the southern Antarctic Peninsula. We emphasise that these are regional average (many hundreds of km) heat flow estimates constrained by seismic data with limited sensitivity to upper crustal composition. (C) 2019 Elsevier B.V. All rights reserved.

Publisher's copyright statement: c © The Author(s) 2018. Published by Oxford University Press on behalf of The Royal Astronomical Society. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Differences in predictions of Glacial Isostatic Adjustment (GIA) for Antarctica persist due to uncertainties in deglacial history and Earth rheology. The Earth models adopted in many GIA studies are defined by parameters that vary in the radial direction only and represent a global average Earth structure (referred to as 1-D Earth models). Oversimplifying the actual Earth structure leads to bias in model predictions in regions where Earth parameters differ significantly from the global average, such as West Antarctica. We investigate the impact of lateral variations in lithospheric thickness on GIA in Antarctica by carrying out two experiments that use different rheological approaches to define 3-D Earth models that include spatial variations in lithospheric thickness. The first experiment defines an elastic lithosphere with spatial variations in thickness inferred from seismic studies. We compare the results from this 3-D model with results derived from a 1-D Earth model that has a uniform lithospheric thickness defined as the average of the 3-D lithospheric thickness. Irrespective of the deglacial history and sublithospheric mantle viscosity, we find higher gradients of present-day uplift rates (i.e. higher amplitude and shorter wavelength) in West Antarctica when using the 3-D models, due to the thinner-than-1-D-average lithosphere prevalent in this region. The second experiment uses seismically inferred temperature as an input to a power-law rheology, thereby allowing the lithosphere to have a viscosity structure. Modelling the lithosphere with a power-law rheology results in a behaviour that is equivalent to a thinner lithosphere model, and it leads to higher amplitude and shorter wavelength deformation compared with the first experiment. We conclude that neglecting spatial variations in lithospheric thickness in GIA models will result in predictions of peak uplift and subsidence that are biased low in West Antarctica. This has important implications for ice-sheet modelling studies as the steeper gradients of uplift predicted from the more realistic 3-D model may promote stability in marine-grounded regions of West Antarctica. Including lateral variations in lithospheric thickness, at least to the level of considering West and East Antarctica separately, is important for capturing short-wavelength deformation and it has the potential to provide a better fit to Global Positioning System observations as well as an improved GIA correction for the Gravity Recovery and Climate Experiment data.

Many studies of Holocene relative sea-level (RSL) changes across Antarctica assume that their reconstructions record uplift from glacial isostatic adjustment caused by the demise of the Last Glacial Maximum (LGM) ice sheets. However, recent analysis of GPS observations suggests that mantle viscosity beneath the Antarctic Peninsula is weaker than previously thought, which would imply that solid Earth motion is not controlled by post-LGM ice-sheet retreat but instead by late Holocene ice-mass changes. If this hypothesis is correct, one might expect to find Holocene RSL records that do not reflect a monotonic decrease in the rate of RSL fall but show variations in the rate of RSL change through the Holocene. We present a new record of late Holocene RSL change from Torgersen Island near Palmer Station in the western Antarctic Peninsula that shows an increase in the rate of relative sea-level fall from 3.0 +/- 1.2 mm/yr to 5.1 +/- 1.8 mm/yr during the late Holocene. Independent studies of the glacial history of the region provide evidence of ice-sheet changes over similar time scales that may be driving this change. When our RSL records are corrected for sea-surface height changes associated with glacial isostatic adjustment (GIA), the rate of post-0.79 ka land uplift at Torgersen Island, 53 1.8 mm/yr, is much higher than the rate of uplift recorded at a nearby GPS site at Palmer Station prior to the Larsen B breakup in 2002 AD (1998-2002 AD; <0.1 mm/yr), but similar to the rates observed after 2002 AD (2002-2013 AD; 6-9 mm/yr). This substantial variation in uplift rates further supports the hypothesis that Holocene RSL rates of change are recording responses to late Holocene and recent changes in local ice loading rather than a post-LGM signal across portions of the Antarctic Peninsula. Thus middle-to-late Holocene RSL data may not be an effective tool for constraining the size of the LGM ice sheet across portions of the Antarctic Peninsula underlain by weaker mantle. In addition, current global-scale GIA models are unable to predict our observed changes in late Holocene RSL. Complexities in Earth structure and neoglacial history need to be taken into consideration in GIA models used for correcting modern satellite-based observations of ice-mass loss. (C) 2018 Elsevier Ltd. All rights reserved.