Imaging the poro-elastic properties of glacier beds using ambient seismic noise monitoring : application to Whillans ice stream, Antarctica

Author(s):  
Aurélien Mordret ◽  
Gauthier Guerin ◽  
Diane Rivet ◽  
Brad Lipovsky ◽  
Brent Minchew

<p>Part of the movement that occurs on all glaciers in Antarctica is a continuous and stable movement that unloads the ice into the sea. The Whillans Ice Plain (WIP) is a portion of the Whillans ice stream that measures 8000 km² for an ice thickness of 800 meters. This glacier has a unique characteristic of moving thanks to tidally modulated stick-slip events twice a day. The slip speed varies laterally across the glacier.  We measured surface wave velocity variations computed from ambient seismic noise cross-correlation. The cross-correlations make it possible to monitor temporally and spatially the seismic velocities at the bed of the glacier, associated with changes in poro-elastic parameters and frictional properties of the glacial till. We averaged our observations for the 78 stick-slip events of our dataset and managed to achieve a 5 min temporal resolution along the 45 min long slip events. The results show a decrease in velocity of about 9% of the S-wave velocity in the subglacial sediment layer about 30 minutes after the initiation of the slip. This velocity drop mainly affects the central part of the glacier. A 10% increase in porosity could induce this velocity decrease due to dilatancy. Dilatant strengthening results from this porosity increase, which in turn keeps the glacier in a slow-sliding regime. The high rate of seismic cycles on such a large scale makes the Whillans ice stream a unique laboratory to study transient aseismic slips in glacial context but also in active tectonic faults one. </p>

2011 ◽  
Vol 305 (3-4) ◽  
pp. 283-289 ◽  
Author(s):  
J. Paul Winberry ◽  
Sridhar Anandakrishnan ◽  
Douglas A. Wiens ◽  
Richard B. Alley ◽  
Knut Christianson

Author(s):  
J. Paul Winberry ◽  
Sridhar Anandakrishnan ◽  
Richard B. Alley ◽  
Robert A. Bindschadler ◽  
Matt A. King

2019 ◽  
Vol 124 (2) ◽  
pp. 1601-1625 ◽  
Author(s):  
Paul M. Bremner ◽  
Mark P. Panning ◽  
R. M. Russo ◽  
Victor Mocanu ◽  
A. Christian Stanciu ◽  
...  

2021 ◽  
Author(s):  
◽  
Rachel Heckels

<p>Ambient seismic noise is used to examine the spatial and temporal surface wave velocity structures and ambient seismic noise fields in the vicinity of different fault zone environments. This study focuses on two distinct regions of central South Island, New Zealand. The Canterbury Plains is a sedimentary basin with many minor faults, which was considered to have low seismic hazard prior to the 2010 – 2011 Canterbury earthquake sequence. We focus on the time period immediately following the 2010 Darfield earthquake, which ruptured the previously unmapped Greendale Fault. The second region of interest is the central Southern Alps. The locked portion of the Alpine Fault currently poses one of the largest seismic hazards for New Zealand. The wealth of data from both permanent and temporary seismic deployments in these regions make them ideal areas in which to assess the effectiveness of ambient noise for velocity modelling in regions surrounding faults at different stages of their seismic cycles.  Temporal velocity changes are measured following the Mw 7.1 Darfield earthquake of 4 September 2010 in the Canterbury Plains. Nine-component cross-correlations are computed from temporary and permanent seismic stations lying on and surrounding the Greendale Fault. Using the Moving-Window Cross-Spectral method, surface wave velocity changes are calculated for the four months immediately following the earthquake until 10 January 2011, for 0.1 — 1.0 Hz. An average increase in seismic velocity of 0.14 ± 0.04 % is determined throughout the region, providing the first such estimate of postseismic relaxation rates in Canterbury. Depth analyses further showed that velocity changes are confined to the uppermost 5 km of the subsurface and we attribute this to postseismic relaxation via crack-healing of the Greendale Fault and throughout the surrounding region.  Rayleigh and Love wave dispersion is examined throughout the Canterbury region. Multi-component cross-correlation functions are analysed for group and phase dispersion curves. These are inverted using frequency-time analysis for 2-D phase and group velocity maps of Rayleigh and Love waves. A high-velocity zone to the southeast of the region coincides with volcanic rocks of Banks Peninsula. Dispersion curves generated from the surface wave tomography are further inverted for one-dimensional shear velocity profiles. These models show a thin, low-velocity near surface layer consistent with the basin sediments, which thins towards the foothills of the Southern Alps. A near-surface damage zone is identified along the length of the Greendale Fault, with consistent reduced Vs velocities to depth of up to 5 km.  Surface and shear wave velocity maps are computed for the central Southern Alps to image the seismic structure of the region. Tomographic surface maps at periods of 5 – 12 s are produced from dispersion measurements of three-component cross-correlation functions. At periods of 5 – 8 s a strong NE-SW trending velocity contrast highlights the Alpine Fault. One-dimensional shear velocity models, computed from the surface wave maps, are in agreement with previous models produced by other conventional methods. An analysis of surface wave amplitudes through signal-to-noise ratios of cross-correlations reveals strong directional effects. Calculated signal-to-noise ratios are up to eight times higher for surface waves travelling north-west than for waves travelling to the south or east. We attribute this to a combination of more energetic ocean wave signals from the Southern Ocean compared to the Tasman Sea.</p>


2016 ◽  
Vol 10 (1) ◽  
pp. 385-399 ◽  
Author(s):  
B. P. Lipovsky ◽  
E. M. Dunham

Abstract. During the 200 km-scale stick slip of the Whillans Ice Plain (WIP), West Antarctica, seismic tremor episodes occur at the ice–bed interface. We interpret these tremor episodes as swarms of small repeating earthquakes. The earthquakes are evenly spaced in time, and this even spacing gives rise to spectral peaks at integer multiples of the recurrence frequency ∼ 10–20 Hz. We conduct numerical simulations of the tremor episodes that include the balance of forces acting on the fault, the evolution of rate- and state-dependent fault friction, and wave propagation from the fault patch to a seismometer located on the ice. The ice slides as an elastic block loaded by the push of the upstream ice, and so the simulated basal fault patch experiences a loading velocity equal to the velocity observed by GPS receivers on the surface of the WIP. By matching synthetic seismograms to observed seismograms, we infer fault patch area ∼ 10 m2, bed shear modulus ∼ 20 MPa, effective pressure ∼ 10 kPa, and frictional state evolution distance ∼ 1 μm. Large-scale slip events often occur twice daily, although skipped events have been increasing in frequency over the last decade. The amplitude of tremor (recorded by seismometers on the ice surface) is greater during the double wait time events that follow skipped events. The physical mechanism responsible for these elevated amplitudes may provide a window into near-future subglacial conditions and the processes that occur during ice-stream stagnation.


Author(s):  
J Salomón ◽  
C Pastén ◽  
S Ruiz ◽  
F Leyton ◽  
M Sáez ◽  
...  

Summary The seismic response of the Santiago City, the capital of Chile with more than 5.5 million inhabitants, is controlled by the properties of the shallower quaternary deposits and the impedance contrast with the underlying Abanico formation, among other factors. In this study, we process continuous records of ambient seismic noise to perform an ambient seismic noise tomography with the aim of defining the shallower structure of the Abanico formation underneath the densely populated metropolitan area of Santiago, Chile. The seismic signals were recorded by a network consisting of 29 broadband seismological stations and 12 accelerograph stations, located in a 35 × 35 km2 quadrant. We used the average coherency of the vertical components to calculate dispersion curves from 0.1 to 5 Hz and Bootstrap resampling to estimate the variance of the travel times. The reliable frequency band of the dispersion curves was defined by an empirical method based on sign normalization of the coherency real part. The ambient noise tomography was solved on a domain discretized into 256 2 × 2 km2 cells. Using a regularized weighted least squares inversion, we inverted the observed travel-times between stations, assuming straight ray paths, in order to obtain 2D phase velocity maps from 0.2 Hz to 1.1 Hz, linearly spaced every 0.05 Hz, in 157 of the 256 square cells of the domain. In each square cell with information, dispersion curves were assembled and used to invert shear wave velocity profiles, which were interpolated using the ordinary Kriging method to obtain a 3D shear wave velocity model valid from 0.6 to 5 km depth. The 3D velocity model shows that the Abanico formation is stiffer in the south of the study area with larger velocity anomalies towards the shallower part of the model. The value of the shear wave velocity narrows with depth, reaching an average value of 3.5 km/s from 3 to 5 km depth.


2016 ◽  
Vol 121 (7) ◽  
pp. 1295-1309 ◽  
Author(s):  
J. R. Leeman ◽  
R. D. Valdez ◽  
R. B. Alley ◽  
S. Anandakrishnan ◽  
D. M. Saffer

2014 ◽  
Vol 60 (222) ◽  
pp. 795-807 ◽  
Author(s):  
J. Paul Winberry ◽  
Sridhar Anandakrishnan ◽  
Richard B. Alley ◽  
Douglas A. Wiens ◽  
Martin J. Pratt

AbstractWe summarize new observations of the deceleration and stick–slip motion of Whillans Ice Stream (WIS), Antarctica. We refine the location of the large sticky spots that resist motion between slip events, the locations of which are controlled by the patterns of subglacial water flow. Our examination of the long-term velocity time series for the ice stream reveals that the decadal-scale deceleration is not occurring at a steady rate, but varies at the sub-decadal timescale. This unsteady deceleration modulates the temporal evolution of a broad (~50 km across) surface-elevation bulge forming at the junction between the relatively narrow upstream portion of the ice stream and broad ice plain that constitutes the downstream end of WIS. Comparison of observations from April 2003 and November 2010 reveals significant changes in the tidally modulated stick–slip cycle that regulates motion on the ice plain. We observe that the timing of slip events has become less regular in response to decreased flow speed in the upstream portions of the ice stream. The decreased regularity of slip events has reduced the release of stored elastic strain during slip events, increasing the rate of deceleration.


Sign in / Sign up

Export Citation Format

Share Document