Mapping crustal shear wave velocity structure and radial anisotropy beneath West Antarctica using seismic ambient noise

Using 8‐25s period Rayleigh and Love wave phase velocity dispersion data extracted from seismic ambient noise, we (i) model the 3D 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 ~26‐30km thick extending south from Marie Byrd Land separates domains of more extended crust (~22km 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 (TAM) and Haag Nunataks‐Ellsworth Whitmore Mountains (HEW) block by a sharp crustal thickness gradient. Crust ~35‐40km is modelled beneath the Haag Nunataks‐Ellsworth Mountains, decreasing to ~30‐32km 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 mid crust is positively radially anisotropic (VSH > VSV) across West Antarctica. The strongest anisotropic signature is observed in the HEW block, emphasising its unique provenance amongst West Antarctica's crustal units, and conceivably reflects a ~13km 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 continental crust

Glacial Isostatic Adjustment in response to changing Late Holocene behaviour of ice streams on the Siple Coast, West Antarctica

The Siple Coast region of Antarctica contains a number of fast-flowing ice streams, which control the dynamics and mass balance of the region. These ice streams are known to undergo stagnation and reactivation cycles, which lead to ice thickness changes that may be sufficient to excite a viscous solid Earth response (glacial isostatic adjustment; GIA). This study aims to quantify Siple Coast ice thickness changes during the last 2000 yr in order to determine the degree to which they might contribute to GIA and associated present-day bedrock uplift rates. This is important because accurate modelling of GIA is necessary to determine the rate of present-day ice-mass change from satellite gravimetry. Recently-published reconstructions of ice-stream variability were used to create a suite of kinematic models for the stagnation-related thickening of Kamb Ice Stream since ∼1850 AD, and a GIA model was used to predict present-day deformation rates in response to this thickening. A number of longer-term loading scenarios, which include the stagnation and reactivation of ice streams across the Siple Coast over the past 2000 yr, were also constructed, and used to investigate the longer term GIA signal in the region. Uplift rates for each of the ice loading histories, based on a range of earth models, were compared with regional GPS-observed uplift rates and an empirical GIA estimate. We estimate Kamb Ice Stream to have thickened by 70–130 m since stagnation ∼165 years ago. Modelled present-day vertical motion in response to this load increase peaks at −17 mm yr–1 (i.e. 17 mm yr–1 subsidence) for the weakest earth models tested here. Comparison of the solid Earth response to ice load changes throughout the last glacial cycle, including ice stream stagnation and reactivation across the Siple Coast during the last 2000 yr, with an empirical GIA estimate suggests that the upper mantle viscosity of the region is greater than 1 × 1020 Pa s. When upper mantle viscosity values of 1 × 1020 Pa s or smaller are combined with our suite of ice-load scenarios we predict uplift rates across Siple Coast that are at least 4 mm yr–1 smaller than those predicted by the empirical GIA estimate. GPS data are unable to further constrain model parameters due to the distance of the GPS sites from the study area. Our results demonstrate that Late Holocene ice load changes related to the stagnation and reactivation of ice streams on the Siple Coast may play a dominant role in defining the present-day uplift signal. However, both the detailed Earth structure and deglacial history of the region need to be better constrained in order to reduce uncertainties associated with the GIA signal of this region

The impact of lateral variations in lithospheric thickness on glacial isostatic adjustment in West Antarctica

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 1D Earth models). Over-simplifying 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 3D 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 3D model with results derived from a 1D Earth model that has a uniform lithospheric thickness defined as the average of the 3D lithospheric thickness. Irrespective of deglacial history and sub-lithospheric mantle viscosity, we find higher gradients of present-day uplift rates (i.e. higher amplitude and shorter wavelength) in West Antarctica when using the 3D models, due to the thinner-than-1D-average lithosphere prevalent in this region. The second experiment uses seismically-inferred temperature as 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 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 3D 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 GPS observations as well as an improved GIA correction for GRACE data

Rapid bedrock uplift in the Antarctic Peninsula explained by viscoelastic response to recent ice unloading

Since 1995 several ice shelves in the Northern Antarctic Peninsula have collapsed and triggered ice-mass unloading, invoking a solid Earth response that has been recorded at continuous GPS (cGPS) stations. A previous attempt to model the observation of rapid uplift following the 2002 breakup of Larsen B Ice Shelf was limited by incomplete knowledge of the pattern of ice unloading and possibly the assumption of an elastic-only mechanism. We make use of a new high resolution dataset of ice elevation change that captures ice-mass loss north of 66°S to first show that non-linear uplift of the Palmer cGPS station since 2002 cannot be explained by elastic deformation alone. We apply a viscoelastic model with linear Maxwell rheology to predict uplift since 1995 and test the fit to the Palmer cGPS time series, finding a well constrained upper mantle viscosity but less sensitivity to lithospheric thickness. We further constrain the best fitting Earth model by including six cGPS stations deployed after 2009 (the LARISSA network), with vertical velocities in the range 1.7 to 14.9 mm/yr. This results in a best fitting Earth model with lithospheric thickness of 100–140 km and upper mantle viscosity of View the MathML source6×1017–2×1018 Pas – much lower than previously suggested for this region. Combining the LARISSA time series with the Palmer cGPS time series offers a rare opportunity to study the time-evolution of the low-viscosity solid Earth response to a well-captured ice unloading event

Uplift rates from a new high-density GPS network in Palmer Land indicate significant late Holocene ice loss in the southwestern Weddell Sea

The measurement of ongoing ice-mass loss and associated melt water contribution to sea-level change from regions such as West Antarctica is dependent on a combination of remote sensing methods. A key method, the measurement of changes in Earth's gravity via the GRACE satellite mission, requires a potentially large correction to account for the isostatic response of the solid Earth to ice-load changes since the Last Glacial Maximum. In this study, we combine glacial isostatic adjustment modelling with a new GPS dataset of solid Earth deformation for the southern Antarctic Peninsula to test the current understanding of ice history in this region. A sufficiently complete history of past ice-load change is required for glacial isostatic adjustment models to accurately predict the spatial variation of ongoing solid Earth deformation, once the independently-constrained effects of present-day ice mass loss have been accounted for. Comparisons between the GPS data and glacial isostatic adjustment model predictions reveal a substantial misfit. The misfit is localized on the southwestern Weddell Sea, where current ice models under-predict uplift rates by approximately 2 mm yr−1. This under-prediction suggests that either the retreat of the ice sheet grounding line in this region occurred significantly later in the Holocene than currently assumed, or that the region previously hosted more ice than currently assumed. This finding demonstrates the need for further fieldwork to obtain direct constraints on the timing of Holocene grounding line retreat in the southwestern Weddell Sea and that GRACE estimates of ice sheet mass balance will be unreliable in this region until this is resolved

The uppermost mantle seismic velocity structure of West Antarctica from Rayleigh wave tomography: insights into tectonic structure and geothermal heat flow

We present a shear wave model of the West Antarctic upper mantle 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 ∼5% between this anomaly and the inner WARS translates to a temperature difference of ∼125-200 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 ∼85 km for the WARS, Marie Byrd Land, and Thurston Island block. This increases to ∼96 km in the Ellsworth Mountains. A surface heat flow of ∼60mW/m2 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 ∼50mW/m2 geotherm best explains the velocities in the Ellsworth Mountains, and a ∼60mW/m2 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