Array processing in cryoseismology: a comparison to network-based approaches at an Antarctic ice stream

Seismicity at glaciers, ice sheets, and ice shelves provides observational constraint on a number of glaciologi- cal processes. Detecting and locating this seismicity, specifi- cally 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 re- quire a smaller number of instruments. Here, we investigate the potential of arrays and array-processing methods to de- tect and locate subsurface microseismicity at glaciers, bench- marking 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 de- tects 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 mar- gins, where it would be challenging to deploy instruments. Finally, we make some practical recommendations for future array deployments at glaciers

Seismic Noise Interferometry and Distributed Acoustic Sensing (DAS): Inverting for the Firn Layer S ‐Velocity Structure on Rutford Ice Stream, Antarctica

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

How dynamic are ice-stream beds?

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 radar-sounded basal interface. Given the precision of our measurements, the upper limit of subglacial erosion observed here is 500 mm a‾¹, far exceeding erosion rates reported for glacial settings from proglacial sediment yields, but substantially below subglacial erosion rates of 1.0 m a‾¹ previously reported from repeat geophysical surveys in West Antarctica

Highly variable friction and slip observed at Antarctic ice stream bed

The slip of glaciers over the underlying bed is the dominant mechanism governing the migration of ice from land into the oceans, with accelerating slip contributing to sea-level rise. Yet glacier slip remains poorly understood, and observational constraints are sparse. Here we use passive seismic observations to measure both frictional shear stress and slip at the bed of the Rutford Ice Stream in Antarctica using 100,000 repetitive stick-slip icequakes. We find that basal shear stresses and slip rates vary from 10$^4$ to 10$^7$ Pa and 0.2 to 1.5 m per day, respectively. Friction and slip vary temporally over the order of hours, and spatially over 10s of metres, due to corresponding variations in effective normal stress and ice–bed interface material. Our findings suggest that the bed is substantially more complex than currently assumed in ice stream models and that basal effective normal stresses may be significantly higher than previously thought. Our observations can provide constraints on the basal boundary conditions for ice-dynamics models. This is critical for constraining the primary contribution of ice mass loss in Antarctica and hence for reducing uncertainty in sea-level rise projections

Constraining recent ice flow history at Korff Ice Rise, West Antarctica, using radar and seismic measurements of ice fabric

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

Automated detection of basal icequakes and discrimination from surface crevassing

Icequakes at or near the bed of a glacier have the potential to allow us to investigate the interaction of ice with the underlying till or bedrock. Understanding this interaction is important for studying basal sliding of glaciers and ice streams, a critical process in ice dynamics models used to constrain future sea level rise projections. However, seismic observations on glaciers can be dominated by seismic energy from surface crevassing. We present a method of automatically detecting basal icequakes and discriminating them from surface crevassing, comparing this method to a commonly used spectrum-based method of detecting icequakes. We use data from Skeidararjo ̈kull, an outlet glacier of the Vatnaj ̈okull Ice Cap, South-East Iceland, to demonstrate that our method outperforms the commonly used spectrum-based method. Our method detects a higher number of basal icequakes, has a lower rate of incorrectly identifying crevassing as basal icequakes and detects an additional, spatially independent basal icequake cluster. We also show independently that the icequakes do not originate from near the glacier surface. We conclude that the method described here is more effective than currently implemented methods for detecting and discriminating basal icequakes from surface crevassing

A multiphase seismic investigation of the shallow subduction zone, southern North Island, New Zealand

The shallow structure of the Hikurangi margin, in particular the interface between the Australian Plate and the subducting Pacific Plate, is investigated using the traveltimes of direct and converted seismic phases from local earthquakes. Mode conversions take place as upgoing energy from earthquakes in the subducted slab crosses the plate interface. These PS and SP converted arrivals are observed as intermediate phases between the direct P and S waves. They place an additional constraint on the depth of the interface and enable the topography of the subducted plate to be mapped across the region. 301 suitable earthquakes were recorded by the Leeds (Tararua) broad-band seismic array, a temporary line of three-component short-period stations, and the permanent stations of the New Zealand national network. This provided coverage across the land area of southern North Island, New Zealand, at a total of 17 stations. Rays are traced through a structure parametrized using layered B-splines and the traveltime residuals inverted, simultaneously, for hypocentre relocation, interface depth and seismic velocity. The results are consistent with sediment in the northeast of the study region and gentle topography on the subducting plate. This study and recent tectonic reconstructions of the southwest Pacific suggest that the subducting plate consists of captured, oceanic crust. The anomalous nature of this crust partly accounts for the unusual features of the Hikurangi margin, e.g. the shallow trench, in comparison with the subducting margin further north

Bathymetry and bed conditions of Lago Subglacial CECs, West Antarctica

Although over 600 Antarctic subglacial lakes have been identified using radar and satellite observations, the bathymetry and bed properties, which are key to understanding conditions within the lake, have been determined in very few localities. We present measurements of water column thickness and lakebed properties from Lago Subglacial CECs (SLC), located beneath 2653 m of ice at the Rutford-Institute-Minnesota divide in Antarctica. Seismic profiles indicate a maximum water column thickness of 301.3 ± 1.5 m, at the widest part of the lake, with an estimated lake volume of 2.5 ± 0.3 km3. Seismic imaging and measurements of the reflection strength at the ice base and lakebed indicate >15 m of high-porosity fine-grained sediment in the central section of the lakebed, consistent with a depositional sequence with an age of up to 0.5 Ma. These observations, along with previous radar measurements and modelling, indicate a low-energy sedimentary environment with a long water-residence time. As such, SLC is a suitable target for exploration via direct access to recover sediment records of ice sheet and climate history and investigate microbial life with long periods of isolation

Downhole distributed acoustic seismic profiling at Skytrain Ice Rise, West Antarctica

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