The ability to monitor transient motion along faults is critical to improving our ability to understand many natural phenomena such as landslides and earthquakes. Here, we usedata from a GPS and ...seismometer network that were deployed to monitor the regularly repeating glacial earthquakes of Whillans Ice Stream, West Antarctica to show that a unique pattern of precursory slip precedes complete rupture along the bed of the ice stream. Additionally, we show that rupture can be independently tracked by increased levels of microseismic activity, including harmonic tremor, that are coincident with the onset of slip at any location, thus providing a remote means of monitoring stress and rupture propagation during the glacial earthquakes.
Key Points
Whillans Ice Stream glacial earthquakes begin with a nucleation phase
Sliding may be tracked by microseismic emissions
Harmonic tremor accompanies glacial earthquakes
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
Antarctica's fast-flowing ice streams drain the ice sheet, with their velocity modulated by subglacial water systems. Current knowledge of these water systems is limited to the shallow portions near ...the ice-bed interface, but hypothesized deeper groundwater could also influence ice streaming. Here, we use magnetotelluric and passive seismic data from Whillans Ice Stream, West Antarctica, to provide the first observations of deep sub-ice stream groundwater. Our data reveal a volume of groundwater within a >1-kilometer-thick sedimentary basin that is more than an order of magnitude larger than the known subglacial system. A vertical salinity gradient indicates exchange between paleo seawater at depth and contemporary basal meltwater above. Our results provide new constraints for subglacial water systems that affect ice streaming and subglacial biogeochemical processes.
Totten Glacier is a fast‐moving East Antarctic outlet with the potential for significant future sea‐level contributions. We deployed four autonomous phase‐sensitive radars on its ice shelf to monitor ...ice‐ocean interactions near its grounding zone and made active source seismic observations to constrain gravity‐derived bathymetry models. We observe an asymmetry in basal melting with mean melt rates along the grounding zone differing by up to 20 m/a. Our new bathymetry model reveals that this melt rate asymmetry coincides with an asymmetry in water column thickness and that the low‐melting ice‐shelf portion is shielded from the main cavity circulation. A 2‐year record yields year‐to‐year melt rate variability of 7–9 m/a with no seasonal cycle. Our results highlight the key role of bathymetry near grounding lines for accurate modeling of ice‐shelf melt, and the importance of sustained multi‐year monitoring, especially at ice‐shelf cavities where the dominant melt rate drivers vary primarily inter‐annually.
Plain Language Summary
The point were the Antarctic Ice Sheet goes afloat on the ocean represents a critical region, where minor variations in melt rates can impact glacier flow and influence the rate of sea‐level rise. East Antarctica's Totten Glacier holds the potential to raise global sea level by several meters. Therefore, to understand the conditions it is exposed to, we measured melt rates for 2 years in several key locations near the point where the ice first touches the ocean. Our new measurements of the shape of the Totten Ice Shelf cavity help explain an observed spatial pattern of basal melting and together with local melt rate data resolve a disagreement between existing melt rate estimates from remote‐sensing methods.
Key Points
Totten Glacier melt rates vary spatially between 0 and over 20 m/a; differences are explained by water column thickness variations from updated bathymetry
Temporal melt rate variability is primarily inter‐annual; melt rates differ by 7–9 m/a over two observed years and there is no clear seasonal cycle
Contrary to previous findings, we find no topographic barriers to the intrusion of warm water to the Totten Glacier grounding zone
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
Crary Ice Rise formed after the Ross Ice Shelf re-grounded ~1 kyr BP. We present new ice-penetrating radar data from two systems operating at center frequencies of 7 and 750 MHz that confirm the ice ...rise is composed of a former ice shelf buried by subsequent accumulation. Stacks of englacial diffraction hyperbolas are present almost everywhere across the central ice rise and extend up to ~350 m above the bed. In many cases, bed reflections beneath the diffraction hyperbolas are obscured for distances up to 1 km. Waveform modeling indicates that the diffraction hyperbolas are likely caused by marine ice deposits in former basal crevasses and rifts. The in-filling of rifts and basal crevasses may have strengthened the connection between the ice rise and the surrounding ice shelf, which could have influenced local and regional ice dynamics. Three internal reflection horizons mark the upper limit of disturbed ice and diffraction hyperbolas in different sections of the ice rise, indicating at least three stages of flow stabilization across the ice rise. A surface lineation visible in MODIS imagery corresponds spatially to deepening and strong deformation of these layers, consistent with the characteristics of former grounding lines observed elsewhere in Antarctica.
Cryoseismology is a powerful toolset for progressing the understanding of the structure and dynamics of glaciers and ice sheets. It can enable the detection of hidden processes such as brittle ...fracture, basal sliding, transient hydrological processes, and calving. Addressing the challenge of detecting signals from many different processes, we present a novel approach for the semi-automated detection of events and event-like noise, which is well-suited for use as Part 1 of a workflow where unsupervised machine learning will be used as Part 2 (Latto et al., 2024) to facilitate the main reconnaissance of diverse detected event types. Implemented in the open-source and widely used ObsPy Python package, the multi-STA/LTA algorithm constructs a hybrid characteristic function from a set of short-term average (sta)–long-term average (lta) pairs (refer to Sect. 2 in the main text for an explanation of how uppercase and lowercase STA/sta and LTA/lta abbreviations are differentiated). We apply the algorithm to data from a seismic array deployed on the Whillans Ice Stream (WIS) in West Antarctica (austral summer 2010–2011) to form a “catch-all” catalogue of events and event-like noise. The new algorithm compares favorably with standard approaches, yielding a diversity of seismic events, including all previously identified stick-slip events (Pratt et al., 2014), teleseisms, and other noise-type signals. In terms of a catalogue overview, we investigate a partial association of seismicity with the tidal cycle and a slight association with ice temperature changes of the Antarctic summer. The new algorithm and workflow will assist in the comparison of different glacier environments using seismology, the identification of process change over time, and the targeting of possible subsequent high-resolution studies.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
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 (Mw 7), seismic amplitudes are modest (Ms 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.
Full text
Available for:
DOBA, IJS, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
We have located 117 previously undetected seismic events mainly occurring between 2015 and 2017 that originated from glacial, tectonic, and volcanic processes in central West Antarctica using data ...recorded on Polar Earth Observing Network (POLENET/ANET) and UK Antarctic Network (UKANET) seismic stations. The seismic events, with local magnitudes (ML) ranging from 1.1 to 3.5, are predominantly clustered in four geographic regions; the Ellsworth Mountains, Thwaites Glacier, Pine Island Glacier, and Mount Takahe. Eighteen of the events are in the Ellsworth Mountains and can be attributed to a mixture of glacial and tectonic processes. The largest event noted in this study was a mid‐crustal (∼19 km focal depth; ML 3.5) normal mechanism earthquake beneath Thwaites Glacier. We also located 91 glacial events near the grounding zones of Thwaites Glacier and Pine Island Glacier that are predominantly associated with time periods of significant calving activity. Eight events, likely arising from volcano‐tectonic processes, occurred beneath Mount Takahe. Using Pn travel times from the seismic events, we find laterally variable uppermost mantle structure in central West Antarctica. On average, the Ellsworth Mountains are underlain by a faster mantle lid (VPn = ∼8.4 km/s) compared to the Amundsen Sea Embayment region (VPn = ∼8.1 km/s). Within the Amundsen Sea Embayment itself, we find mantle lid velocities ranging from ∼8.05 to 8.18 km/s. Laterally heterogeneous uppermost mantle structure, indicative of variable thermal and rheological structure, likely influences both geothermal heat flux and glacial isostatic adjustment spatial patterns and rates within central West Antarctica.
Plain Language Summary
In this study, we have identified and located 117 seismic events in central West Antarctica. The seismic events originate from both solid‐Earth and glacial processes. While the seismic events originating from solid‐Earth processes are associated with tectonic and volcano‐tectonic activity in central West Antarctica, the seismic events arising from glacial processes are predominantly associated with time periods of significant calving activity at Thwaites Glacier and Pine Island Glacier. Using select P‐waves produced by seismic events located in this study, we investigate the uppermost portion of the mantle beneath central West Antarctica and find that the temperature of the uppermost mantle likely varies by at least 160 K across central West Antarctica. Variable uppermost mantle temperatures have numerous implications for the dynamics and evolution of the West Antarctic Ice Sheet.
Key Points
Located 117 previously undetected seismic events of glacial, tectonic, and volcano‐tectonic origin
Laterally variable Pn velocity structure across central West Antarctica
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
Given the high number and diversity of events in a typical cryoseismic dataset, in particular those recorded on ice sheet margins, it is desirable to use a semi-automated method of grouping similar ...events for reconnaissance and ongoing analysis. We present a workflow for employing semi-unsupervised cluster analysis to inform investigations of the processes occurring in glaciers and ice sheets. In this demonstration study, we make use of a seismic event catalogue previously compiled for the Whillans Ice Stream, for the 2010–2011 austral summer (outlined in Part 1, Latto et al., 2024). We address the challenges of seismic event analysis for a complex wave field by clustering similar seismic events into groups using characteristic temporal, spectral, and polarization attributes of seismic time series with the k-means++ algorithm. This provides the basis for a reconnaissance analysis of a seismic wave field that contains local events (from the ice stream) set in an ambient wave field that itself contains a diversity of signals (mostly from the Ross Ice Shelf). As one result, we find that two clusters include stick-slip events that diverge in terms of length and initiation locality (i.e., central sticky spot and/or the grounding line). We also identify a swarm of high-frequency signals on 16–17 January 2011 that are potentially associated with a surface melt event from the Ross Ice Shelf. Used together with the event detection presented in Part 1, the semi-automated workflow could readily be generalized to other locations and, as a possible benchmark procedure, could enable the monitoring of remote glaciers over time and comparisons between locations.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Changing rates of water input can affect both the flow of glaciers and ice sheets and their propensity to crevasse. Here we examine geodetic and seismic observations during two substantial ...(10–18-times background velocity) rain-induced glacier accelerations at Haupapa/Tasman Glacier, New Zealand. Changes in rain rate result in glacier acceleration and associated uplift, which propagate down-glacier. This pattern of acceleration results in a change to the strain rate field, which correlates with an order of magnitude increase in the apparent seismicity rate and an overall down-glacier migration in located seismicity. After each acceleration event the apparent seismicity rate decreases to below the pre-acceleration rate for 3 days. This suggests that seismic events associated with surface crevasse growth occur early during phases of glacier acceleration due to elevated extensional stresses, and then do not occur again until stresses recover.
Understanding the processes that affect streaming ice flow and the mass balance of glaciers and ice sheets requires sound knowledge of their subglacial environments. Previous studies have shown that ...an extensive deformable subglacial sediment layer favors fast ice-stream flow. However, areas of high basal drag, termed sticky spots, are of particular interest because they inhibit the fast flow of the overriding ice. The stick-slip behavior of Whillans Ice Stream (WIS) is perhaps the most conspicuous manifestation of a subglacial sticky spot. We present new ice-thickness and seismic-reflection measurements collected over the main sticky spot in the ice plain of WIS, allowing us to elucidate its role in the behavior of the ice stream. Ice-thickness and surface-elevation data show that the sticky spot occupies a subglacial topographic high. Water flow in response to the hydrological potential gradient will be routed around the sticky spot if effective pressures are similar on the sticky spot and elsewhere. The seismic experiment imaged a laterally continuous basal layer approximately 6 m thick, having compressional wave velocities of greater than 1800 m s−1 and density greater than 1800 kg m−3, indicative of a till layer that is stiffer than corresponding till beneath well-lubricated parts of the ice stream. This layer likely continues to deform under the higher shear stress of the sticky spot, and some water may be pumped up onto the sticky spot during motion events.
•A laterally continuous basal till layer of approximately 6 m thickness is imaged.•Stiff till beneath the sticky spot on Whillans Ice Stream (WIS) influences its flow.•High bed elevation and low hydropotential are the likely cause of the sticky spot.•Effective pressure at the sticky spot is higher than the streaming part of the WIS.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK, ZRSKP