Using glacier seismicity for phase velocity measurements and Green's function retrieval

High-melt areas of glaciers and ice sheets foster a rich spectrum of ambient seismicity. These signals not only shed light on source mechanisms (e.g. englacial fracturing, water flow, iceberg detachment, basal motion) but also carry information about seismic wave propagation within glacier ice. Here, we present two approaches to measure and potentially monitor phase velocities of high-frequency seismic waves (≥1 Hz) using naturally occurring glacier seismicity. These two approaches were developed for data recorded by on-ice seasonal seismic networks on the Greenland Ice Sheet and a Swiss Alpine glacier. The Greenland data set consists of continuous seismograms, dominated by long-term tremor-like signals of englacial water flow, whereas the Alpine data were collected in triggered mode producing 1-2 s long records that include fracture events within the ice (‘icequakes'). We use a matched-field processing technique to retrieve frequency-dependent phase velocity measurements for the Greenland data. In principle, this phase dispersion relationship can be inverted for ice sheet thickness and bed properties. For these Greenland data, inversion of the dispersion curve yields a bedrock depth of 541 m, which may be too small by as much as 35 per cent. We suggest that the discrepancy is due to lateral changes in ice sheet depth and bed properties beneath the network, which may cause unaccounted mixing of surface wave modes in the dispersion curve. The Swiss Alpine icequake records, on the other hand, allow for reconstruction of the impulse response between two seismometers. The direct and scattered wave fields from the vast numbers of icequake records (tens of thousands per month) can be used to measure small changes in englacial velocities and thus monitor structural changes within the ic

Seismic noise interferometry reveals transverse drainage configuration beneath the surging Bering Glacier

Subglacial drainage systems are known to critically control ice flows, but their spatial configuration and temporal evolution are poorly constrained due to inaccessibility. Here we report a 12‐year‐long monitoring of the drainage underneath Bering Glacier, Alaska, by correlating ambient noise recorded at two seismic stations on the sides of the glacier. We find that the seismic surface waves traveling across Bering Glacier slowed down by 1–2% during its latest 2008–2011 surge, likely due to the switch of the subglacial drainage from a channelized system to a distributed system. In contrast to current models, the relative amplitude of velocity reductions for Rayleigh and Love waves requires the distributed drainage to be highly anisotropic and aligned perpendicular to the ice flow direction. We infer that the subglacial water flow is mainly through a network of transverse basal crevasses during surges and thus can sustain the high water pressure and ice flow speed

Applications of SAR Interferometry in Earth and Environmental Science Research

This paper provides a review of the progress in regard to the InSAR remote sensing technique and its applications in earth and environmental sciences, especially in the past decade. Basic principles, factors, limits, InSAR sensors, available software packages for the generation of InSAR interferograms were summarized to support future applications. Emphasis was placed on the applications of InSAR in seismology, volcanology, land subsidence/uplift, landslide, glaciology, hydrology, and forestry sciences. It ends with a discussion of future research directions

Mapping crustal shear wave velocity structure and radial anisotropy beneath West Antarctica using seismic ambient noise

Using 8‐25s period Rayleigh and Love wave phase velocity dispersion data extracted from seismic ambient noise, we (i) model the 3D shear wave velocity structure of the West Antarctic crust and (ii) map variations in crustal radial anisotropy. Enhanced regional resolution is offered by the UK Antarctic Seismic Network. In the West Antarctic Rift System (WARS), a ridge of crust ~26‐30km thick extending south from Marie Byrd Land separates domains of more extended crust (~22km thick) in the Ross and Amundsen Sea Embayments, suggesting along‐strike variability in the Cenozoic evolution of the WARS. The southern margin of the WARS is defined along the southern Transantarctic Mountains (TAM) and Haag Nunataks‐Ellsworth Whitmore Mountains (HEW) block by a sharp crustal thickness gradient. Crust ~35‐40km is modelled beneath the Haag Nunataks‐Ellsworth Mountains, decreasing to ~30‐32km km thick beneath the Whitmore Mountains, reflecting distinct structural domains within the composite HEW block. Our analysis suggests that the lower crust and potentially the mid crust is positively radially anisotropic (VSH > VSV) across West Antarctica. The strongest anisotropic signature is observed in the HEW block, emphasising its unique provenance amongst West Antarctica's crustal units, and conceivably reflects a ~13km thick metasedimentary succession atop Precambrian metamorphic basement. Positive radial anisotropy in the WARS crust is consistent with observations in extensional settings, and likely reflects the lattice‐preferred orientation of minerals such as mica and amphibole by extensional deformation. Our observations support a contention that anisotropy may be ubiquitous in continental crust

Mapping Crustal Shear Wave Velocity Structure and Radial Anisotropy Beneath West Antarctica Using Seismic Ambient Noise

Using 8- to 25-s-period Rayleigh and Love wave phase velocity dispersion data extracted from seismic ambient noise, we (i) model the 3-D shear wave velocity structure of the West Antarctic crust and (ii) map variations in crustal radial anisotropy. Enhanced regional resolution is offered by the UK Antarctic Seismic Network. In the West Antarctic Rift System (WARS), a ridge of crust ∼26–30km thick extending south from Marie Byrd Land separates domains of more extended crust (∼22km thick) in the Ross and Amundsen Sea Embayments, suggesting along-strike variability in the Cenozoic evolution of the WARS. The southern margin of the WARS is defined along the southern Transantarctic Mountains and Haag-Ellsworth Whitmore Mountains (HEW) block by a sharp crustal thickness gradient. Crust ∼35–40km is modeled beneath the Haag Nunataks-Ellsworth Mountains, decreasing to ∼30–32km thick beneath the Whitmore Mountains, reflecting distinct structural domains within the composite HEW block. Our analysis suggests that the lower crust and potentially the middle crust is positively radially anisotropic (VSH \u3e VSV) across West Antarctica. The strongest anisotropic signature is observed in the HEW block, emphasizing its unique provenance among West Antarctica\u27s crustal units, and conceivably reflects a ∼13-km-thick metasedimentary succession atop Precambrian metamorphic basement. Positive radial anisotropy in the WARS crust is consistent with observations in extensional settings and likely reflects the lattice-preferred orientation of minerals such as mica and amphibole by extensional deformation. Our observations support a contention that anisotropy may be ubiquitous in the continental crust

Seismic Array Studies of Antarctica and Madagascar

The scope of this dissertation is broad, involving seismic array studies from Antarctica and Madagascar, and includes aspects of glaciology and oceanography as well as solid Earth geophysics. Chapter 2 focuses on the study of stickslip motion of the Whillans Ice Stream, West Antarctica. It includes methods combining seismic array and GPS time series, from ice stream based-sensors, to determine source dynamics in the framework of an earthquake source. The source characteristics are then analyzed to explain far-field seismic observations of ice stream- sourced surface waves detected throughout West Antarctica. Locations of asperities, or sticky- spots, that cause the Whillans Ice Stream to accelerate and generate seismic energy are found. Some of these asperities are in close proximity to the grounding line, where properties of the bed are altered through tidal flexure of the ice shelf and the influx of water into the subglacial till. Chapter 3 explores ocean generated microseismic noise that is also detected on these ice stream seismometers, with the geometry of the array providing excellent azimuthal resolution. Stacked cross-correlations of seismograms enhance microseismic energy generated by the Southern Ocean in the form of both surface and body waves. The frequency spectra of these waves is analyzed as well as applying seismic array techniques, such as beamforming. Each frequency band provides different information on the source regions of that particular microseism suggesting multiple source mechanisms. Microseisms are modeled using ocean state hindcasts to compare with observations and identify microseism source regions and improve understanding of the effect of sea ice. It is shown that single-frequency microseisms are heavily damped by the presence of sea ice over the continental shelf. Long-period double frequency microseisms are observed and modeled to be sourced in the deep ocean. Short-period double frequency microseisms are also influenced by sea ice seasonality; however, this chapter provides evidence that shows that a component of this band may be sourced in the deep ocean. The focus of Chapter 4 moves away from Antarctica, to Madagascar and the analysis of the first island-wide deployment of broadband seismometers. The priorities of this project are to better understand the crustal and upper mantle structure of Madagascar, and to assess the intraplate volcanism on the island from a seismological point of view for the first time. This chapter presents a surface wave tomography study producing the first shear velocity model of the crust and upper mantle of the island. A range of commonly employed surface wave methods is used to calculate phase velocities across the island. These are then amalgamated and inverted for shear velocity in the crust and the upper mantle. Low velocity regions are shown to extend to upper mantle depths beneath the center and north of the island above which lie intraplate volcanic provinces. This suggests that the mantle lithosphere has been significantly thinned, explaining the relatively high topography observed for a fragment of continental crust

Characterisation of the subglacial environment using geophysical constrained Bayesian inversion techniques

An accurate characterization of the inaccessible subglacial environment is key to accurately modelling the dynamic behaviour of ice sheets and glaciers, crucial for predicting sea-level rise. The composition and water content of subglacial material can be inferred from measurements of shear wave velocity (Vs) and bulk electrical resistivity (R), themselves derived from Rayleigh wave dispersion curves and transient electromagnetic (TEM) soundings. Conventional Rayleigh wave and TEM inversions can suffer from poor resolution and non-uniqueness. In this thesis, I present a novel constrained inversion methodology which applies a Markov chain Monte Carlo implementation of Bayesian inversion to produce probability distributions of geophysical parameters. MuLTI (Multimodal Layered Transdimensional Inversion) is used to derive Vs from Rayleigh wave dispersion curves, and its TEM variant, MuLTI-TEM, for evaluating bulk electrical resistivity. The methodologies can include independent depth constraints, drawn from external data sources (e.g., boreholes or other geophysical data), which significantly improves the resolution compared to conventional unconstrained inversions. Compared to such inversions, synthetic studies suggested that MuLTI reduces the error between the true and best-fit models by a factor of 10, and reduces the vertically averaged spread of the Vs distribution twofold, based on the 95% credible intervals. MuLTI and MuLTI-TEM were applied to derive Vs and R profiles from seismic and TEM electromagnetic data acquired on the terminus of the Norwegian glacier Midtdalsbreen. Three subglacial material classifications were determined: sediment (Vs 1600 m/s, R > 500 Ωm) and weathered/fractured bedrock containing saline water (Vs > 1900 m/s, R < 50 Ωm). These algorithms offer a step-change in our ability to resolve and quantify the uncertainties in subsurface inversions, and show promise for constraining the properties of subglacial aquifers beneath Antarctic ice masses. MuLTI and MuLTITEM have both been made publicly available via GitHub to motivate users, in the cryosphere and other environmental settings, for continued advancement

Estimation of full moment tensors, including uncertainties, for earthquakes, volcanic events, and nuclear explosions

Thesis (Ph.D.) University of Alaska Fairbanks, 2016We present a catalog of full seismic moment tensors for 63 events from Uturuncu volcano in Bolivia. The events were recorded during 2011-2012 in the PLUTONS seismic array of 24 broadband stations. Most events had magnitudes between 0.5 and 2.0 and did not generate discernible surface waves; the largest event was Mw 2.8. For each event we computed the misfit between observed and synthetic waveforms, and we used first-motion polarity measurements to reduce the number of possible solutions. Each moment tensor solution was obtained using a grid search over the six-dimensional space of moment tensors. For each event we show the misfit function in eigenvalue space, represented by a lune. We identify three subsets of the catalog: (1) 6 isotropic events, (2) 5 tensional crack events, and (3) a swarm of 14 events southeast of the volcanic center that appear to be double couples. The occurrence of positively isotropic events is consistent with other published results from volcanic and geothermal regions. Several of these previous results, as well as our results, cannot be interpreted within the context of either an oblique opening crack or a crack-plus-double-couple model. Proper characterization of uncertainties for full moment tensors is critical for distinguishing among physical models of source processes. A seismic moment tensor is a 3×3 symmetric matrix that provides a compact representation of a seismic source. We develop an algorithm to estimate moment tensors and their uncertainties from observed seismic data. For a given event, the algorithm performs a grid search over the six-dimensional space of moment tensors by generating synthetic waveforms for each moment tensor and then evaluating a misfit function between the observed and synthetic waveforms. 'The' moment tensor M₀ for the event is then the moment tensor with minimum misfit. To describe the uncertainty associated with M₀, we first convert the misfit function to a probability function. The uncertainty, or rather the confidence, is then given by the 'confidence curve' P(V), where P(V) is the probability that the true moment tensor for the event lies within the neighborhood of M that has fractional volume V. The area under the confidence curve provides a single, abbreviated 'confidence parameter' for M0. We apply the method to data from events in different regions and tectonic settings: 63 small (Mw 4) earthquakes in the southern Alaska subduction zone, and 12 earthquakes and 17 nuclear explosions at the Nevada Test Site. Characterization of moment tensor uncertainties puts us in better position to discriminate among moment tensor source types and to assign physical processes to the events

Scaling full seismic waveform inversions

The main goal of this research study is to scale full seismic waveform inversions using the adjoint-state method to the data volumes that are nowadays available in seismology. Practical issues hinder the routine application of this, to a certain extent theoretically well understood, method. To a large part this comes down to outdated or flat out missing tools and ways to automate the highly iterative procedure in a reliable way. This thesis tackles these issues in three successive stages. It first introduces a modern and properly designed data processing framework sitting at the very core of all the consecutive developments. The ObsPy toolkit is a Python library providing a bridge for seismology into the scientific Python ecosystem and bestowing seismologists with effortless I/O and a powerful signal processing library, amongst other things. The following chapter deals with a framework designed to handle the specific data management and organization issues arising in full seismic waveform inversions, the Large-scale Seismic Inversion Framework. It has been created to orchestrate the various pieces of data accruing in the course of an iterative waveform inversion. Then, the Adaptable Seismic Data Format, a new, self-describing, and scalable data format for seismology is introduced along with the rationale why it is needed for full waveform inversions in particular and seismology in general. Finally, these developments are put into service to construct a novel full seismic waveform inversion model for elastic subsurface structure beneath the North American continent and the Northern Atlantic well into Europe. The spectral element method is used for the forward and adjoint simulations coupled with windowed time-frequency phase misfit measurements. Later iterations use 72 events, all happening after the USArray project has commenced, resulting in approximately 150`000 three components recordings that are inverted for. 20 L-BFGS iterations yield a model that can produce complete seismograms at a period range between 30 and 120 seconds while comparing favorably to observed data

3C Seismic Interferometry at the Polymetallic Kylylahti Deposit, Outokumpu District, Finland

Non peer reviewe