DAS field dataset to compare technologies and deployment scenarios – Antarctica Dataset

This report describes a Distributed Acoustic Sensing (DAS) dataset acquired by the British Antarctic Survey (BAS) and the University of Oxford in Antarctic during 2020. The field dataset contributes to the Deliverable D1.1 of the DigiMon project (DAS field dataset to compare technologies and deployment scenarios), which is associated with tasks 1.2 and 1.3 of the project

Constraints on North Anatolian Fault Zone Width in the Crust and Upper Mantle From S Wave Teleseismic Tomography

We present high-resolution S wave teleseismic tomography images of the western segment of the North Anatolian Fault Zone (NAFZ) in Turkey using teleseismic data recorded during the deployment period of the Dense Array for Northern Anatolia array. The array comprised 66 stations with a nominal station spacing of 7 km, thus permitting a horizontal and vertical resolution of approximately 15 km. We use the current S wave results with previously published P wave teleseismic tomography to produce maps of relative VP/VSanomalies, which we use to highlight the difference in overall composition of the three terranes separated by the northern (NNAF) and southern branches of the NAFZ. Our results show a narrow S wave low-velocity anomaly beneath the northern branch of the NAFZ extending from the upper crust, where it has a width of ∼10 km, to the lower crust, where it widens to ∼30 km. This low-velocity zone most likely extends into the upper mantle, where we constrain its width to be ≤ 50 km and interpret it as indicative of localized shear beneath the NNAF; this structure is similar to what has been observed for the NAFZ west of 32°, and therefore, we propose that the structure of the NNAF is similar to that of the NAFZ in the east. The southern branch of the NAFZ does not show a very strong signature in our images, and we conclude that it is most likely rooted in the crust, possibly accommodating deformation related to rotation of the Armutlu/Almacik Blocks situated between the two NAFZ branches

Oceanic and atmospheric forcing of Larsen C Ice-Shelf thinning

The catastrophic collapses of Larsen A and B ice shelves on the eastern Antarctic Peninsula have caused their tributary glaciers to accelerate, contributing to sea-level rise and freshening the Antarctic Bottom Water formed nearby. The surface of Larsen C Ice Shelf (LCIS), the largest ice shelf on the peninsula, is lowering. This could be caused by unbalanced ocean melting (ice loss) or enhanced firn melting and compaction (englacial air loss). Using a novel method to analyse eight radar surveys, this study derives separate estimates of ice and air thickness changes during a 15-year period. The uncertainties are considerable, but the primary estimate is that the surveyed lowering (0.066 ± 0.017 m yr−1) is caused by both ice loss (0.28 ± 0.18 m yr−1) and firn-air loss (0.037 ± 0.026 m yr−1). The ice loss is much larger than the air loss, but both contribute approximately equally to the lowering because the ice is floating. The ice loss could be explained by high basal melting and/or ice divergence, and the air loss by low surface accumulation or high surface melting and/or compaction. The primary estimate therefore requires that at least two forcings caused the surveyed lowering. Mechanisms are discussed by which LCIS stability could be compromised in the future. The most rapid pathways to collapse are offered by the ungrounding of LCIS from Bawden Ice Rise or ice-front retreat past a "compressive arch" in strain rates. Recent evidence suggests that either mechanism could pose an imminent risk

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

Differentiating flow, melt, or fossil seismic anisotropy beneath Ethiopia

Ethiopia is a region where continental rifting gives way to oceanic spreading. Yet the role that pre-existing lithospheric structure, melt, mantle flow, or active upwellings may play in this process is debated. Measurements of seismic anisotropy are often used to attempt to understand the contribution that these mechanisms may play. In this study, we use new data in Afar, Ethiopia along with legacy data across Ethiopia, Djibouti, and Yemen to obtain estimates of mantle anisotropy using SKS-wave splitting. We show that two layers of anisotropy exist, and we directly invert for these. We show that fossil anisotropy with fast directions oriented northeast-southwest may be preserved in the lithosphere away from the rift. Beneath the Main Ethiopian Rift and parts of Afar, anisotropy due to shear segregated melt along sharp changes in lithospheric thickness dominates the shear-wave splitting signal in the mantle. Beneath Afar, away from regions with significant lithospheric topography, melt pockets associated with the crustal and uppermost mantle magma storage dominate the signal in localized regions. In general, little anisotropy is seen in the uppermost mantle beneath Afar suggesting melt retains no preferential alignment. These results show the important role melt plays in weakening the lithosphere and imply that as rifting evolves passive upwelling sustains extension. A dominant northeast-southwest anisotropic fast direction is observed in a deeper layer across all of Ethiopia. This suggests that a conduit like plume is lacking beneath Afar today, rather a broad flow from the southwest dominates flow in the upper mantle

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

Subglacial lakes and hydrology across the Ellsworth Subglacial Highlands, West Antarctica

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

Measuring seismic attenuation in polar firn: method and application to Korff Ice Rise, West Antarctica

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

Distributed Acoustic Sensing (DAS) for natural microseismicity studies: A case study from Antarctica

Icequakes, microseismic earthquakes at glaciers, offer insights into the dynamics of ice sheets. For the first time in the Antarctic, we explore the use of fiber optic cables as Distributed Acoustic Sensors (DAS) as a new approach for monitoring basal icequakes. We present the use of DAS for studying icequakes as a case study for the application of DAS to microseismic datasets in other geological settings. Fiber was deployed on the ice surface at Rutford Ice Stream in two different configurations. We compare the performance of DAS with a conventional geophone network for: microseismic detection and location; resolving source and noise spectra; source mechanism inversion; and measuring anisotropic shear-wave splitting parameters. Both DAS array geometries detect fewer events than the geophone array. However, DAS is superior to geophones for recording the microseism signal, suggesting the applicability of DAS for ambient noise interferometry. We also present the first full-waveform source mechanism inversions using DAS anywhere, successfully showing the horizontal stick-slip nature of the icequakes. In addition, we develop an approach to use a 2D DAS array geometry as an effective multi-component sensor capable of accurately characterising shear-wave splitting due to the anisotropic ice fabric. Although our observations originate from a glacial environment, the methodology and implications of this work are relevant for employing DAS in other microseismic environments

Deep crustal melt plumbing of Bárðarbunga volcano, Iceland

Understanding magmatic plumbing within the Earth's crust is important for understanding volcanic systems and improving eruption forecasting. We discuss magma plumbing under Bárðarbunga volcano, Iceland, over a 4 year period encompassing the largest Icelandic eruption in 230 years. Microseismicity extends through the usually ductile region of the Earth's crust, from 7 to 22 km depth in a subvertical column. Moment tensor solutions for an example earthquake exhibits opening tensile crack behavior. This is consistent with the deep (>7 km) seismicity being caused by the movement of melt in the normally aseismic crust. The seismically inferred melt path from the mantle source is offset laterally from the center of the Bárðarbunga caldera by ~12 km, rather than lying directly beneath it. It is likely that an aseismic melt feed also exists directly beneath the caldera and is aseismic due to elevated temperatures and pervasive partial melt under the caldera