We construct a new seismic model for central and West Antarctica by jointly inverting Rayleigh wave phase and group velocities along with P wave receiver functions. Ambient noise tomography exploiting data from more than 200 seismic stations deployed over the past 18 years is used to construct Rayleigh wave phase and group velocity dispersion maps. Comparison between the ambient noise phase velocity maps with those constructed using teleseismic earthquakes confirms the accuracy of both results. These maps, together with P receiver function waveforms, are used to construct a new 3‐D shear velocity (Vs) model for the crust and uppermost mantle using a Bayesian Monte Carlo algorithm. The new 3‐D seismic model shows the dichotomy of the tectonically active West Antarctica (WANT) and the stable and ancient East Antarctica (EANT). In WANT, the model exhibits a slow uppermost mantle along the Transantarctic Mountains (TAMs) front, interpreted as the thermal effect from Cenozoic rifting. Beneath the southern TAMs, the slow uppermost mantle extends horizontally beneath the traditionally recognized EANT, hypothesized to be associated with lithospheric delamination. Thin crust and lithosphere observed along the Amundsen Sea coast and extending into the interior suggest involvement of these areas in Cenozoic rifting. EANT, with its relatively thick and cold crust and lithosphere marked by high Vs, displays a slower Vs anomaly beneath the Gamburtsev Subglacial Mountains in the uppermost mantle, which we hypothesize may be the signature of a compositionally anomalous body, perhaps remnant from a continental collision.

Long-period seismic sources associated with glacier motion have been recently discovered, and an increase in ice flow over the past decade has been suggested on the basis of secular changes in such measurements. Their significance, however, remains uncertain, as a relationship to ice flow has not been confirmed by direct observation. Here we combine long-period surface-wave observations with simultaneous Global Positioning System measurements of ice displacement to study the tidally modulated stick-slip motion of the Whillans Ice Stream in West Antarctica. The seismic origin time corresponds to slip nucleation at a region of the bed of the Whillans Ice Stream that is likely stronger than in surrounding regions and, thus, acts like an 'asperity' in traditional fault models. In addition to the initial pulse, two seismic arrivals occurring 10-23 minutes later represent stopping phases as the slip terminates at the ice stream edge and the grounding line. Seismic amplitude and average rupture velocity are correlated with tidal amplitude for the different slip events during the spring-to-neap tidal cycle. Although the total seismic moment calculated from ice rigidity, slip displacement, and rupture area is equivalent to an earthquake of moment magnitude seven (M(w) 7), seismic amplitudes are modest (M(s) 3.6-4.2), owing to the source duration of 20-30 minutes. Seismic radiation from ice movement is proportional to the derivative of the moment rate function at periods of 25-100 seconds and very long-period radiation is not detected, owing to the source geometry. Long-period seismic waves are thus useful for detecting and studying sudden ice movements but are insensitive to the total amount of slip.

The downstream portion of Whillans Ice Stream, West Antarctica, moves primarily by stick‐slip motion. The observation of stick‐slip motion suggests that the bed is governed by velocity‐weakening physics and that the basal physics is more unstable than suggested by laboratory studies. The stick‐slip cycle of Whillans Ice Plain exhibits substantial variability in both the duration of sticky periods and in slip magnitude. To understand this variability, we modeled the forces acting on the ice stream during the stick phase of the stick‐slip cycle. The ocean tides introduce changes in the rate at which stress is applied to the ice plain. Increased loading rates promote earlier failure and vice versa. Results show that the bed of Whillans Ice Stream strengthens over time (healing) during the quiescent intervals in the stick‐slip cycle, with the bed weakening during slip events. The time‐dependent strengthening of the ice plain bed following termination of slip events indicates that the strength of the bed may vary by up to 0.35 kPa during the course of a single day.

P‐to‐S receiver functions (PRFs) from the Polar Earth Observing Network (POLENET) GPS and seismic leg of POLENET spanning West Antarctica and the Transantarctic Mountains deployment of seismographic stations provide new estimates of crustal thickness across West Antarctica, including the West Antarctic Rift System (WARS), Marie Byrd Land (MBL) dome, and the Transantarctic Mountains (TAM) margin. We show that complications arising from ice sheet multiples can be effectively managed and further information concerning low‐velocity subglacial sediment thickness may be determined, via top‐down utilization of synthetic receiver function models. We combine shallow structure constraints with the response of deeper layers using a regularized Markov chain Monte Carlo methodology to constrain bulk crustal properties. Crustal thickness estimates range from 17.0±4 km at Fishtail Point in the western WARS to 45±5 km at Lonewolf Nunataks in the TAM. Symmetric regions of crustal thinning observed in a transect deployment across the West Antarctic Ice Sheet correlate with deep subice basins, consistent with pure shear crustal necking under past localized extension. Subglacial sediment deposit thicknesses generally correlate with trough/dome expectations, with the thickest inferred subice low‐velocity sediment estimated as ∼0.4 km within the Bentley Subglacial Trench. Inverted PRFs from this study and other published crustal estimates are combined with ambient noise surface wave constraints to generate a crustal thickness map for West Antarctica south of 75°S. Observations are consistent with isostatic crustal compensation across the central WARS but indicate significant mantle compensation across the TAM, Ellsworth Block, MBL dome, and eastern and western sectors of thinnest WARS crust, consistent with low density and likely dynamic, low‐viscosity high‐temperature mantle.

The seismic velocity structure of Antarctica is important, both as a constraint on the tectonic history of the continent and for understanding solid Earth interactions with the ice sheet. We use Rayleigh wave array analysis methods applied to teleseismic data from recent temporary broadband seismograph deployments to image the upper mantle structure of central and West Antarctica. Phase velocity maps are determined using a two–plane wave tomography method and are inverted for shear velocity using a Monte Carlo approach to estimate three‐dimensional velocity structure. Results illuminate the structural dichotomy between the East Antarctic Craton and West Antarctica, with West Antarctica showing thinner crust and slower upper mantle velocity. West Antarctica is characterized by a 70–100 km thick lithosphere, underlain by a low‐velocity zone to depths of at least 200 km. The slowest anomalies are beneath Ross Island and the Marie Byrd Land dome and are interpreted as upper mantle thermal anomalies possibly due to mantle plumes. The central Transantarctic Mountains are marked by an uppermost mantle slow‐velocity anomaly, suggesting that the topography is thermally supported. The presence of thin, higher‐velocity lithosphere to depths of about 70 km beneath the West Antarctic Rift System limits estimates of the regionally averaged heat flow to less than 90 mW/m2. The Ellsworth‐Whitmore block is underlain by mantle with velocities that are intermediate between those of the West Antarctic Rift System and the East Antarctic Craton. We interpret this province as Precambrian continental lithosphere that has been altered by Phanerozoic tectonic and magmatic activity.