Constraining subglacial geology using ambient noise Rayleigh wave ellipticity

Abstract

B asal slip is an important mechanism by which glaciers and ice-sheets flow, and is a major source of uncertainty in simulations of ice-mass loss and sea level rise from the Greenland Ice Sheet (GrIS). Sub-ice geology is a dominant control on ice flow velocity with fast flow often coinciding with the presence of deformable subglacial till eroded from underlying sedimentary rocks. The subglacial geology of Greenland has received relatively little attention thus far and its control on ice flow is poorly understood. Seismic studies of the crust beneath the GrIS have been limited due to a lack of seismic stations and the reliance on earthquake event data. However, in the past decade, there has been a rapid increase in the number of both permanent and temporary seismic stations deployed in Greenland as well developments in ambient noise methods, allowing for improved spatial resolution of crustal geology. Ellipticity measurements, the ratio of the horizontal to vertical component of a Rayleigh wave, have been shown to be particularly sensitive to the geological structure directly beneath the station. Ambient noise H/V measurements have been used for decades in geotechnical and civil engineering for site characterisation, making them a well-suited technique to determine the subglacial geology of the GrIS. Using all available broadband stations deployed on Greenland from 2012 to 2018 we extract Rayleigh wave ellipticity measurement from ambient noise data using the degree-of-polarization (DOP) method where meaningful signals are defined as a waveform with an arbitrary polarization which remains stable for a given time window. We invert these ellipticity measurements in the period range of 4 ¿ 9 s to generate Vs profiles of the first 5 km beneath each station. Our inversions indicate that: (1) off-ice stations along the margins of the GrIS produce a good agreement with the litho1.0 model to within error and (2) an additional subglacial layer 1.0 - 2.0km thick with a Vs < 3.0km is necessary to match the data recorded at several of the on-ice stations. We attribute these observations to the widespread presence of sedimentary rocks beneath the GrIS, potentially capable of sustaining extensive subglacial till layers that can support enhanced basal slip.