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. [ABSTRACT FROM AUTHOR]

The accuracy of sea-level rise predictions is highly dependent on reliably understanding the subglacial environment beneath ice streams. Bedforms result from the interaction between ice and its substrate, and therefore have the potential to shed light on poorly understood basal conditions influencing ice dynamics. However, theoretical models of bedform formation are often based on observations from deglaciated areas or from sparsely-sampled geophysical surveys over glaciated regions. Here, we use high-resolution three- and two-dimensional radar and seismic data to reveal details of the initiation and evolution of a subglacial bedform beneath Rutford Ice Stream, West Antarctica. Radar surveys at 20 m and 50 m line spacing allow detailed imaging of bed topography, including a moat up to 55 m deep, surrounding the upstream end of a 50-m high and >18-km long bedform. Many models rely on either a topographical or a locally resistant seed point to initiate bedform formation. The bedform described here is mostly composed of soft sediment (porosity >0.3) and lacks a boulder or outcrop, suggesting the bedform initiated without a topographic seed point. Sediment at the upstream end of bedforms appears stiffer over a distance of 2.3 km. We suggest sediment inhomogeneities in the initially flat bed cause the deposition of sediment, which, assuming resistant enough, acts as a seed point for bedform extension and moat erosion. The moat's geometry and its truncation of other bedforms suggest that it was eroded after the deposition of surrounding bedforms. These observations from a modern ice stream deliver information of subglacial processes involved in the initiation as well as in situ high-resolution topography and properties of bedforms and moats. Using these observations numerical models can be tested and developed accordingly. • Evidence for a bedform initiating without pre-existing topographic perturbations. • A moat around an elongated bedform revealed from high-resolution radar data. • Constraints of bedform formation mechanism from a modern ice stream. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In recent decades, rapid ice shelf disintegration along the Antarctic Peninsula has had a global impact through enhancing outlet glacier flow and hence sea level rise and the freshening of Antarctic Bottom Water. Ice shelf thinning due to basal melting results from the circulation of relatively warm water in the underlying ocean cavity. However, the effect of sub-shelf circulation on future ice shelf stability cannot be predicted accurately with computer simulations if the geometry of the ice shelf cavity is unknown. To address this deficit for Larsen C Ice Shelf, West Antarctica, we integrate new water column thickness measurements from recent seismic campaigns with existing observations. We present these new data here along with an updated bathymetry grid of the ocean cavity. Key findings include a relatively deep seabed to the southeast of the Kenyon Peninsula, along the grounding line and around the key ice shelf pinning-point of Bawden Ice Rise. In addition, we can confirm that the cavity's southern trough stretches from Mobiloil Inlet to the open ocean. These areas of deep seabed will influence ocean circulation and tidal mixing and will therefore affect the basal-melt distribution. These results will help constrain models of ice shelf cavity circulation with the aim of improving our understanding of sub-shelf processes and their potential influence on ice shelf stability. The datasets are comprised of all the new point measurements of seabed depth. We present the new depth measurements here, as well as a compilation of previously published measurements. To demonstrate the improvements to the sub-shelf bathymetry map that these new data provide we include a gridded data product in the Supplement of this paper, derived using the additional measurements of both offshore seabed depth and the thickness of grounded ice. The underlying seismic datasets that were used to determine bed depth and ice thickness are available at 10.5285/315740B1-A7B9-4CF0-9521-86F046E33E9A (Brisbourne et al., 2019), 10.5285/5D63777D-B375-4791-918F-9A5527093298 (Booth, 2019), 10.5285/FFF8AFEE-4978-495E-9210-120872983A8D (Kulessa and Bevan, 2019) and 10.5285/147BAF64-B9AF-4A97-8091-26AEC0D3C0BB (Booth et al., 2019). [ABSTRACT FROM AUTHOR]

In recent decades, rapid ice-shelf disintegration along the Antarctic Peninsula has had a global impact through enhancing outlet glacier flow, and hence sea level rise, and the freshening of Antarctic Bottom Water. Ice shelf thinning due to basal melting results from the circulation of relatively warm water in the underlying ocean cavity. However, the effect of sub-shelf circulation on future ice-shelf stability cannot be predicted accurately with computer simulations if the geometry of the ice-shelf cavity is unknown. To address this deficit for Larsen C Ice Shelf, west Antarctica, we integrate new water-column thickness measurements with existing observations. We present these new data here along with an updated bathymetry grid of the ocean cavity. Key findings include relatively deep seabed to the south-east of the Kenyon Peninsula, along the grounding line and around the key ice shelf pinning point of Bawden Ice Rise. In addition, we can confirm that the cavity’s southern trough stretches from Mobiloil Inlet to the open ocean. These areas of deep seabed will influence ocean circulation and tidal mixing, and will therefore affect the basal-melt distribution. These results will help constrain models of ice-shelf cavity circulation with the aim of improving our understanding of sub-shelf processes and their potential influence on ice shelf stability. The data set comprises all point measurements of seabed depth and a gridded data product, derived using additional measurements of both offshore seabed depth and the thickness of grounded ice. We present all new depth measurements here as well as a compilation of previously published measurements used in the gridding process. The gridded data product is included in the supplementary material. [ABSTRACT FROM AUTHOR]

A passive seismology experiment was conducted across the main overdeepening of Storglaciaren in the Tarfala valley, northern Sweden, to investigate the spatial and temporal distribution of basal microseismic waveforms in relation to known dynamics of this small polythermal sub-arctic glacier. The high ablation rate made it difficult to keep geophones buried and well coupled to the glacier during the experiment and reduced the number of days of good-quality data collection. The characterization of typical and atypical waveforms showed that the dominant waveforms were from near-surface events such as crevassing. Waveforms resembling basal microseismic signals were very rare, and seldom observed on more than two seismic stations simultaneously. The analysis of waveforms, amplitudes and particle motions suggested a near-field origin for most events. Even though basal sliding is known to occur in the overdeepening, no convincing examples of basal waveforms were detected, suggesting basal microseismic signals are rare or difficult to detect beneath polythermal glaciers like Storglaciaren. We discuss the reasons for failing to locate basal signals, consider the origin of common waveforms and make recommendations for setting up passive seismology experiments on glaciers with high ablation rates. [ABSTRACT FROM AUTHOR]

Presents the anisotropic structure of the Hikurangi subduction zone in New Zealand. Assessment of the shear-wave anisotropy of surface-wave polarizations; Elaboration of symmetry axis of the anisotropic structure; Anisotropic characteristics in the subducted crust underlying the mantle.

Investigates Rayleigh-wave dispersion to determine the shear-wave velocity structure in North Island, New Zealand. Implication of surface-wave study for the thickness of the upper-mantle; Use of anisotropy to explain the presence of high velocities; Indication of the result of the mantle velocities from the isotropic inversions of the S-wave velocity.

Ice is one of the most anisotropic natural materials but it is common knowledge in our community that anisotropy in ice is: a) a second order factor that can be ignored, b) easily assimilated into isotropic models by tuning certain parameters, or c) so complicated that it is better to ignore. Interestingly, any measurement of anisotropy in polar ice shows that ice develops crystalline preferred-orientation fabrics that produce strong mechanical, optical and dielectric anisotropic properties. Here, we look in detail at polarimetric radar measurements in different regions of East and West Antarctica. We use a ground-based phase-sensitive frequency-modulated continuous-wave radar (ApRES) and a matrix-based model to study the radio-waves depolarization and anisotropic scattering (Fujita et al, 2006). We find that anisotropy is widespread and an excellent archive of past ice flow conditions. We use it here to extract details of the recent deglaciation of the Ross and Ronny ice shelves over the last thousands of years. Also, because the length of the record is related to the advection time of the ice flow, we use our method to show that a candidate ice core location near Dome C, for the Beyond Epica Oldest Ice Core project, has been under steady ice-flow conditions over the last million years or so. [ABSTRACT FROM AUTHOR]

Acoustic impedance is a measurable quantity that can be recovered using active seismic surveys. Its value depends on the speed of sound in the sediment and the sediment density. For granular sediments there are theories that relate the speed of sound to the effective pressure, the difference between ice overburden and subglacial water pressure. This provides a link between the seismic observable and a key control on basal sliding. For sediments that follow a Coulomb friction law, the higher the effective pressure the higher the shear stress supported by the sediments. In this presentation we investigate the links between acoustic impedance, effective pressure and basal drag that are predicted by theories of the speed of sound in granular materials. We use results for basal drag from an ice flow model inversion in Antarctica to predict how acoustic impedance would vary under the hypothesis that sliding is controlled by a Coulomb sliding law. [ABSTRACT FROM AUTHOR]

The density structure of firn has implications for hydrological and climate modelling, and ice-shelf stability. The structure of firn can be evaluated from depth models of seismic velocity, widely obtained with Herglotz–Wiechert inversion (HWI), an approach that considers the slowness of refracted seismic arrivals. However, HWI is strictly appropriate only for steady-state firn profiles and the inversion accuracy can be compromised where firn contains ice layers. In these cases, full waveform inversion (FWI) may yield more success than HWI. FWI extends HWI capabilities by considering the full seismic waveform and incorporates reflected arrivals. Using synthetic firn density profiles, assuming both steady- and non-steady-state accumulation, we show that FWI outperforms HWI for detecting ice slab boundaries (5–80 m thick, 5–80 m deep) and velocity anomalies within firn. FWI can detect slabs thicker than one wavelength (here, 20 m, assuming a maximum frequency of 60 Hz) but requires the starting velocity model to be accurate to ±2.5%. We recommend for field practice that the shallowest layers of velocity models are constrained with ground-truth data. Nonetheless, FWI shows advantages over established methods, and should be considered when the characterisation of firn ice slabs is the goal of the seismic survey. [ABSTRACT FROM AUTHOR]

The density structure of firn has implications for hydrological and climate modelling and for ice shelf stability. The firn structure can be evaluated from depth models of seismic velocity, widely obtained with Herglotz-Wiechert inversion (HWI), an approach that considers the slowness of refracted seismic arrivals. However, HWI is appropriate only for steady-state firn profiles and the inversion accuracy can be compromised where firn contains ice layers. In these cases, Full Waveform Inversion (FWI) can be more successful than HWI. FWI extends HWI capabilities by considering the full seismic waveform and incorporates reflected arrivals, thus offering a more accurate estimate of a velocity profile. We show the FWI characterisation of the velocity model has an error of only 1.7% for regions (vs. 4.2% with HWI) with an ice slab (20 m thick, 40 m deep) in an otherwise steady-state firn profile. [ABSTRACT FROM AUTHOR]

Alex Brisbourne reports on a one-day discussion meeting organized by SEIS-UK for seismologists to discuss current seismological research undertaken in the UK and to consider how the provision of new seismic equipment and training facility may enhance that work. [ABSTRACT FROM AUTHOR]

The density structure of firn has implications for hydrological and climate modelling and for ice shelf stability. The firn structure can be evaluated from depth models of seismic velocity, widely obtained with Herglotz-Wiechert inversion (HWI), an approach that considers the slowness of refracted seismic arrivals. However, HWI is appropriate only for steady-state firn profiles and the inversion accuracy can be compromised where firn contains ice layers. In these cases, Full Waveform Inversion (FWI) can be more successful than HWI. FWI extends HWI capabilities by considering the full seismic waveform and incorporates reflected arrivals, thus offering a more accurate estimate of a velocity profile. We show the FWI characterisation of the velocity model has an error of only 1.7% for regions (vs. 4.2% with HWI) with an ice slab (20 m thick, 40 m deep) in an otherwise steady-state firn profile. [ABSTRACT FROM AUTHOR]

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 stabilizing marine ice beneath LCIS. The model uses a new bathymetry, containing a southern seabed trough discovered using seismic observations. The modeled extent of marine ice, thought to stabilize LCIS, is in good agreement with observations. Experiments applying idealized 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. Plain Language Summary: Satellite observations have revealed a lowering in recent decades of the surface of Larsen C Ice Shelf (LCIS), Antarctica, which has led to concern about its stability. By modeling ocean conditions under LCIS, we investigate the extent to which ocean melting may have caused the ice to thin, leading to the observed lowering, or altered the pattern of marine ice beneath LCIS. Marine ice forms when seawater freezes to the base of the ice shelf, and is thought to stabilize LCIS. The model uses a new seabed data set that contains a wide, deep seabed trough in the south, found by a seismic survey. In modeled ocean warming experiments, an increase in melting is concentrated in this southern region. However, greater lowering has been observed in the north, suggesting that changes in ocean conditions are not responsible for the lowering. The calculated pattern of marine ice at the base of LCIS looks similar to the observed pattern. With a warmer ocean, marine ice is significantly reduced in crucial regions of the ice shelf. This shows that the stability of LCIS is sensitive to changes in ocean conditions and other ice shelves around Antarctica are likely sensitive to these changes too. Key Points: An ocean model with new Larsen C bathymetry shows greatest melting and melt sensitivity where seismic data indicate a deep southern troughThe calculated marine ice distribution under Larsen C Ice Shelf, based on model melt/freeze rates, is in good agreement with observationsA reduction in marine ice with ocean warming implies a threat to Larsen C stability, with wide implications for cold‐water ice shelves [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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, which may be important for subglacial traction, till continuity and landform development. 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 shows that 90% of the 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 radarsounded basal interface. Given the precision of our measurements, the upper limit of subglacial erosion observed here is 500mma-1, far exceeding erosion rates reported for glacial settings from proglacial sediment yields, but substantially below subglacial erosion rates of 1.0ma-1 previously reported from repeat geophysical surveys in West Antarctica. [ABSTRACT FROM AUTHOR]