Search

Seismicity at glaciers, ice sheets, and ice shelves provides observational constraint on a number of glaciological processes. Detecting and locating this seismicity, specifically icequakes, is a necessary first step in studying processes such as basal slip, crevassing, imaging ice fabric, and iceberg calving, for example. Most glacier deployments to date use conventional seismic networks, comprised of seismometers distributed over the entire area of interest. However, smaller-aperture seismic arrays can also be used, which are typically sensitive to seismicity distal from the array footprint and require a smaller number of instruments. Here, we investigate the potential of arrays and array-processing methods to detect and locate subsurface microseismicity at glaciers, benchmarking performance against conventional seismic-network-based methods for an example at an Antarctic ice stream. We also provide an array-processing recipe for body-wave cryoseismology applications. Results from an array and a network deployed at Rutford Ice Stream, Antarctica, show that arrays and networks both have strengths and weaknesses. Arrays can detect icequakes from further distances, whereas networks outperform arrays in more comprehensive studies of a particular process due to greater hypocentral constraint within the network extent. We also gain new insights into seismic behaviour at the Rutford Ice Stream. The array detects basal icequakes in what was previously interpreted to be an aseismic region of the bed, as well as new icequake observations downstream and at the ice stream shear margins, where it would be challenging to deploy instruments. Finally, we make some practical recommendations for future array deployments at glaciers.

We present seismic measurements of the firn column at Korff Ice Rise, West Antarctica, including measurements of compressional-wave velocity and attenuation. We describe a modified spectral-ratio method of measuring the seismic quality factor (Q) based on analysis of diving waves, which, combined with a stochastic method of error propagation, enables us to characterise the attenuative structure of firn in greater detail than has previously been possible. Q increases from 56 ± 23 in the uppermost 12 m to 570 ± 450 between 55 and 77 m depth. We corroborate our method with consistent measurements obtained via primary reflection, multiple, source ghost, and critically refracted waves. Using the primary reflection and its ghost, we find Q = 53 ± 20 in the uppermost 20 m of firn. From the critical refraction, we find Q = 640 ± 400 at 90 m depth. Our method aids the understanding of the seismic structure of firn and benefits characterisation of deeper glaciological targets, providing an alternative means of correcting seismic reflection amplitudes in cases where conventional methods of Q correction may be impossible.

Observations of surface lowering on Larsen C Ice Shelf (LCIS), Antarctica, have prompted concern about its stability. In this study, an ocean model is used to investigate the extent to which changes in ocean forcing may have influenced ice loss and the distribution of stabilising marine ice beneath LCIS. The model uses a new bathymetry, containing a southern seabed trough discovered using seismic observations. The modelled extent of marine ice, thought to stabilise LCIS, is in good agreement with observations. Experiments applying idealised ocean warming yield an increase in melting over the southern trough. This is inconsistent with lowering observed in northern LCIS, suggesting oceanic forcing is not responsible for that signal. The marine ice extent and thickness reduces significantly under ocean warming, implying a high sensitivity of LCIS stability to changes in ocean forcing. This result could have wide implications for other cold-water ice shelves around Antarctica.

Basal properties beneath ice streams and glaciers are known to be a control for ice flow dynamics, hence knowledge of them is crucial for predicting sea level due to changes in glacial dynamics. Basal properties, processes and topography also drive the formation of subglacial landforms. Bed properties beneath Rutford Ice Stream (West Antarctica) have previously been described using seismic acoustic impedance measurements at a sparse spatial coverage. Here, we derive bed properties in a 15 x 17 km grid of surface radar data with coverage and sampling much higher than previous seismic studies. Bed reflection amplitudes in surface radar data were calibrated using sediment porosities (ranging from 0.4 – 0.5) derived from seismic acoustic impedance. We find the bed properties are spatially variable, consisting of low porosity material in some areas and soft sediment in other areas. Comparison of seismic and surface radar data imply the low porosity material to be a consolidated sediment or sedimentary rock. Mega-scale glacial lineations (MSGLs) are ubiquitous on the bed and consist of soft, high porosity, probably deforming sediment, consistent with previous interpretations of MSGLs. We find some MSGLs have high reflectivity on their crest, interpreted as water bodies overlying high porosity sediment, whereas the trough around and the upstream end of some landforms consist of low porosity material. Integrating these different observations, we place constraints on possible explanations for the occurrence of water on the crest of landforms.

Firn densification profiles are an important parameter for ice-sheet mass balance and palaeoclimate studies. One conventional method of investigating firn profiles is using seismic refraction surveys, but these are difficult to upscale to large-area measurements. Distributed acoustic sensing (DAS) presents an opportunity for large-scale seismic measurements of firn with dense spatial sampling and easy deployment, especially when seismic noise is used. We study the feasibility of seismic noise interferometry (SI) on DAS data for characterizing the firn layer at the Rutford Ice Stream, West Antarctica. Dominant seismic energy appears to come from anthropogenic noise and shear-margin crevasses. The DAS cross-correlation interferometry yields noisy Rayleigh wave signals. To overcome this, we present two strategies for cross-correlations: (a) hybrid instruments—correlating a geophone with DAS, and (b) stacking of selected cross-correlation panels picked in the tau-p domain. These approaches are validated with results derived from an active survey. Using the retrieved Rayleigh wave dispersion curve, we inverted for a high-resolution 1D S-wave velocity profile down to a depth of 100 m. The profile shows a “kink” (velocity gradient inflection) at ∼12 m depth, resulting from a change of compaction mechanism. A triangular DAS array is used to investigate directional variation in velocity, which shows no evident variations thus suggesting a lack of azimuthal anisotropy in the firn. Our results demonstrate the potential of using DAS and SI to image the near-surface and present a new approach to derive S-velocity profiles from surface wave inversion in firn studies., Applied Geophysics and Petrophysics

Microseismicity, induced by the sliding of a glacier over its bed, can be used to characterize frictional properties of the ice‐bed interface, which are a key parameter controlling ice stream flow. We use naturally occurring seismicity to monitor spatiotemporally varying bed properties at Rutford Ice Stream, West Antarctica. We locate 230000 micro‐earthquakes with local magnitudes from –2.0 to –0.3 using 90 days of recordings from a 35‐station seismic network located ∼40 km upstream of the grounding line. Events exclusively occur near the ice‐bed interface and indicate predominantly flow‐parallel stick‐slip. They mostly lie within a region of interpreted stiff till and along the likely stiffer part of mega‐scale glacial landforms. Within these regions, micro‐earthquakes occur in spatially (<100 m radius) and temporally (mostly 1‐5 days activity) restricted event‐clusters (up to 4000 events), which exhibit an increase, followed by a decrease, in event magnitude with time. This may indicate event triggering once activity is initiated. Although ocean tides modulate the surface ice flow velocity, we observe little periodic variation in overall event frequency over time and conclude that water content, bed topography and stiffness are the major factors controlling microseismicity. Based on variable rupture mechanisms and spatiotemporal characteristics, we suggest the event‐clusters relate to three end‐member types of bed deformation: (1) continuous creation and seismogenic destruction of small‐scale bed‐roughness, (2) ploughed clasts and (3) flow‐oblique deformation during landform‐formation or along bedrock outcrops. This indicates that multiple processes, simultaneously active during glacial sliding, can accommodate stick‐slip behaviour and that the bed continuously reorganizes.

Antarctic ice sheet history is imprinted in the structure and fabric of the ice column. At ice rises, the signature of ice flow history is preserved due to the low strain rates inherent at these independent ice flow centres. We present results from a distributed acoustic sensing (DAS) experiment at Skytrain Ice Rise in the Weddell Sea sector of West Antarctica, aimed at delineating the englacial fabric to improve our understanding of ice sheet history in the region. This pilot experiment demonstrates the feasibility of an innovative technique to delineate ice rise structure. Both direct and reflected P- and S-wave energy, as well as surface wave energy, are observed using a range of source offsets, i.e. a walkaway vertical seismic profile, recorded using fibre optic cable. Significant noise, which results from the cable hanging untethered in the borehole, is modelled and suppressed at the processing stage. At greater depth where the cable is suspended in drilling fluid, seismic interval velocities and attenuation are measured. Vertical P-wave velocities are high (VINT=3984±218 m s−1) and consistent with a strong vertical cluster fabric. Seismic attenuation is high (QINT=75±12) and inconsistent with previous observations in ice sheets over this temperature range. The signal level is too low, and the noise level too high, to undertake analysis of englacial fabric variability. However, modelling of P- and S-wave travel times and amplitudes with a range of fabric geometries, combined with these measurements, demonstrates the capacity of the DAS method to discriminate englacial fabric distribution. From this pilot study we make a number of recommendations for future experiments aimed at quantifying englacial fabric to improve our understanding of recent ice sheet history.

A programmable borehole measurement system was deployed in hot water drilled ice holes during the ‘Bed Access and Monitoring of Ice Sheet History’ (BEAMISH) project to drill to the bed of the Rutford Ice Stream in West Antarctica. This system operates autonomously (no live data) after deployment, and records borehole diameter (non-contact measurement), water column pressure, heading and inclination. Three cameras, two sideways looking and one vertical, are also included for visual inspection of hole integrity and sediments. The system is small, lightweight (~35.5 kg) and low power using only 6 ‘D’ cell sized lithium batteries, making it ideal for transport and use in remote field sites. The system is 2.81 m long and 165 mm in diameter, and can be deployed attached to the drill hose for measurements during drilling or on its own deployment line afterwards. The full system is discussed in detail, highlighting design strengths and weaknesses. Data from the BEAMISH project are also presented in the form of camera images showing hole integrity, and sensor data used to calculate borehole diameter through the full length of the hole. These data are used to show confidence in hole verticality and subsurface cavity development and connection.

Subglacial water plays an important role in ice sheet dynamics and stability. Subglacial lakes are often located at the onset of ice streams and have been hypothesised to enhance ice flow downstream by lubricating the ice–bed interface. The most recent subglacial-lake inventory of Antarctica mapped nearly 400 lakes, of which ∼ 14 % are found in West Antarctica. Despite the potential importance of subglacial water for ice dynamics, there is a lack of detailed subglacial-water characterisation in West Antarctica. Using radio-echo sounding data, we analyse the ice–bed interface to detect subglacial lakes. We report 33 previously uncharted subglacial lakes and present a systematic analysis of their physical properties. This represents a ∼ 40 % increase in subglacial lakes in West Antarctica. Additionally, a new digital elevation model of basal topography of the Ellsworth Subglacial Highlands was built and used to create a hydropotential model to simulate the subglacial hydrological network. This allows us to characterise basal hydrology, determine subglacial water catchments and assess their connectivity. We show that the simulated subglacial hydrological catchments of the Rutford Ice Stream, Pine Island Glacier and Thwaites Glacier do not correspond to their ice surface catchments.

Icequakes at or near the bed of a glacier have the potential to allow us to investigate theinteraction of ice with the underlying till or bedrock. Understanding this interaction is important forstudying basal sliding of glaciers and ice streams, a critical process in ice dynamics models used to con-strain future sea-level rise projections. However, seismic observations on glaciers can be dominated byseismic energy from surface crevassing. We present a method of automatically detecting basal icequakesand discriminating them from surface crevassing, comparing this method to a commonly used spectrum-based method of detecting icequakes. We use data from Skeidararjökull, an outlet glacier of theVatnajökull Ice Cap, South-East Iceland, to demonstrate that our method outperforms the commonlyused spectrum-based method. Our method detects a higher number of basal icequakes, has a lowerrate of incorrectly identifying crevassing as basal icequakes and detects an additional, spatially inde-pendent basal icequake cluster. We also show independently that the icequakes do not originate fromnear the glacier surface. We conclude that the method described here is more effective than currentlyimplemented methods for detecting and discriminating basal icequakes from surface crevassing.

The crystal orientation fabric of ice reflects its flow history, information which is required to better constrain projections of future ice sheet behavior. Here we present a novel combination of polarimetric phase‐sensitive radar and seismic anisotropy measurements to provide independent and consistent constraints on ice fabric at Korff Ice Rise, within the Weddell Sea sector of West Antarctica. The nature and depth distribution of fabric in the ice column is constrained using the azimuthal variation in (1) the received power anomaly and phase difference of polarimetric vertical radar soundings and (2) seismic velocities and shear wave splitting measurements. Radar and seismic observations are modeled separately to determine the nature and strength of fabric within the ice column. Both methods indicate ice fabric above 200‐m depth which is consistent with present‐day ice‐divide flow. However, both measurements also indicate an oblique girdle fabric below 230‐m depth within the ice column, inconsistent with steady state divide flow. Our interpretation is that this deeper fabric is a remnant fabric from a previous episode of flow, which is currently being overwritten by ongoing fabric development associated with the present‐day flow regime. The preexisting fabric is consistent with ice flow from the south prior to ice‐divide formation, in agreement with models of Holocene ice sheet evolution. These findings apply new constraints to the flow history at Korff Ice Rise prior to divide formation and demonstrate the capacity of radar and seismic measurements to map fabric and thus constrain past ice flow.

The Filchner‐Ronne Ice Shelf, the ocean cavity beneath it and the Weddell Sea that bounds it, form an important part of the global climate system by modulating ice discharge from the Antarctic Ice Sheet and producing cold dense water masses that feed the global thermohaline circulation. A prerequisite for modeling the ice sheet and oceanographic processes within the cavity is an accurate knowledge of the sub‐ice‐sheet bedrock elevation, but beneath the ice shelf where airborne radar cannot penetrate, bathymetric data are sparse. This paper presents new seismic point measurements of cavity geometry from a particularly poorly sampled region south of Berkner Island that connects the Filchner and Ronne ice shelves. An updated bathymetric grid formed by combining the new data with existing datasets reveals several new features. In particular, a sill running between Berkner Island and the mainland could alter ocean circulation within the cavity and change our understanding of paleo‐ice‐stream flow in the region. Also revealed are deep troughs near the grounding lines of Foundation and Support Force ice streams, which provide access for seawater with melting potential. Running an ocean tidal model with the new bathymetry reveals large differences in tidal current velocities, both within the new gridded region and further afield, potentially affecting sub‐ice‐shelf melt rates.

Projections of sea-level rise contributions from West Antarctica's dynamically thinning ice streams contain high uncertainty because some of the key processes involved are extremely challenging to observe. An especially poorly observed parameter is sub-decadal stability of ice-stream beds. Only two previous studies have made repeated geophysical measurements of ice-stream beds at the same locations in different years, but both studies were limited in spatial extent. Here, we present the results from repeat radar measurements of the bed of Pine Island Glacier, West Antarctica, conducted 3–6 years apart, along a cumulative ~ 60 km of profiles. Analysis of the correlation of bed picks between repeat surveys show that 90 % of the ice-stream bed displays no significant change despite the glacier increasing in speed by up to 40 % over the last decade. We attribute the negligible detection of morphological change at the bed of Pine Island Glacier to the ubiquitous presence of a deforming till layer, wherein sediment transport is in steady state, such that sediment is transported along the basal interface without inducing morphological change to the radar-sounded bed. Significant change was only detected in one 500 m section of the bed where a change in bed morphology occurs with a difference in vertical amplitude of 3–5 m. Given the precision of our measurements, the maximum possible erosion rate that could go undetected along our profiles is 500 mm a-1, far exceeding erosion rates reported for glacial settings from proglacial sediment yields, but substantially below subglacial erosion rates of 1000 mm a-1 previously reported from repeat geophysical surveys in West Antarctica.

Here we present new measurements of an anisotropic ice fabric in a fast moving (377 ma−1) ice stream in West Antarctica. We use ∼6000 measurements of shear wave splitting observed in microseismic signals from the bed of Rutford Ice Stream, to show that in contrast to large-scale ice flow models, which assume that ice is isotropic, the ice in Rutford Ice Stream is dominated by a previously unobserved type of partial girdle fabric. This fabric has a strong directional contrast in mechanical properties, shearing 9.1 times more easily along the ice flow direction than across flow. This observed fabric is likely to be widespread and representative of fabrics in other ice streams and large glaciers, suggesting it is essential to consider anisotropy in data-driven models to correctly predict ice loss and future flow in these regions. We show how passive microseismic monitoring can be effectively used to provide these data.

Although 90% of Antarctica's discharge occurs via its fast-flowing ice streams, our ability to project future ice sheet response has been limited by poor observational constraints on the ice-bed conditions used in numerical models to determine basal slip. We have helped address this observational deficit by acquiring and analyzing a series of seismic reflection profiles to determine basal conditions beneath the main trunk and tributaries of Pine Island Glacier (PIG), West Antarctica. Seismic profiles indicate large-scale sedimentary deposits. Combined with seismic reflection images, measured acoustic impedance values indicate relatively uniform bed conditions directly beneath the main trunk and tributaries, comprising a widespread reworked sediment layer with a dilated sediment lid of minimum thickness 1.5 ± 0.4 m. Beneath a slow-moving intertributary region, a discrete low-porosity sediment layer of 7 ± 3 m thickness is imaged. Despite considerable basal topography, seismic observations indicate that a till layer at the ice base is ubiquitous beneath PIG, which requires a highly mobile sediment body to maintain an abundant supply. These results are compatible with existing ice sheet models used to invert for basal shear stress: existing basal conditions upstream will not inhibit further rapid retreat of PIG if the high-friction region currently restraining flow, directly upstream of the grounding line, is breached. However, small changes in the pressure regime at the bed, as a result of stress reorganization following retreat, may result in a less-readily deformable bed and conditions which are less likely to maintain high ice-flow rates.