skip to main content


Title: Ambient high-frequency seismic surface waves in the firn column of central west Antarctica
Abstract Firn is the pervasive surface material across Antarctica, and its structures reflect its formation and history in response to environmental perturbations. In addition to the role of firn in thermally isolating underlying glacial ice, it defines near-surface elastic and density structure and strongly influences high-frequency (> 5 Hz) seismic phenomena observed near the surface. We investigate high-frequency seismic data collected with an array of seismographs deployed on the West Antarctic Ice Sheet (WAIS) near WAIS Divide camp in January 2019. Cross-correlations of anthropogenic noise originating from the approximately 5 km-distant camp were constructed using a 1 km-diameter circular array of 22 seismographs. We distinguish three Rayleigh (elastic surface) wave modes at frequencies up to 50 Hz that exhibit systematic spatially varying particle motion characteristics. The horizontal-to-vertical ratio for the second mode shows a spatial pattern of peak frequencies that matches particle motion transitions for both the fundamental and second Rayleigh modes. This pattern is further evident in the appearance of narrow band spectral peaks. We find that shallow lateral structural variations are consistent with these observations, and model spectral peaks as Rayleigh wave amplifications within similarly scaled shallow basin-like structures delineated by the strong velocity and density gradients typical of Antarctic firn.  more » « less
Award ID(s):
1739027
NSF-PAR ID:
10318913
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Glaciology
ISSN:
0022-1430
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Ambient seismic recordings taken at broad locations across Ross Ice Shelf and a dense array near West Antarctic Ice Sheet (WAIS) Divide, Antarctica, show pervasive temporally variable resonance peaks associated with trapped seismic waves in near-surface firn layers. These resonance peaks feature splitting on the horizontal components, here interpreted as frequency-dependent anisotropy in the firn and underlying ice due to several overlapping mechanisms driven by ice flow. Frequency peak splitting magnitudes and fast/slow axes were systematically estimated at single stations using a novel algorithm and compared with good agreement with active source anisotropy measurements at WAIS Divide determined via active sources recorded on a 1 km circular array. The approach was further applied to the broad Ross Ice Shelf (RIS) array, where anisotropy axes were directly compared with visible surface features and ice shelf flow lines. The near-surface firn, depicted by anisotropy above 30 Hz, was shown to exhibit a novel plastic stretching mechanism of anisotropy, whereby the fast direction in snow aligns with accelerating ice shelf flow. 
    more » « less
  2. Abstract Continuous seismic observations across the Ross Ice Shelf reveal ubiquitous ambient resonances at frequencies >5 Hz. These firn-trapped surface wave signals arise through wind and snow bedform interactions coupled with very low velocity structures. Progressive and long-term spectral changes are associated with surface snow redistribution by wind and with a January 2016 regional melt event. Modeling demonstrates high spectral sensitivity to near-surface (top several meters) elastic parameters. We propose that spectral peak changes arise from surface snow redistribution in wind events and to velocity drops reflecting snow lattice weakening near 0∘C for the melt event. Percolation-related refrozen layers and layer thinning may also contribute to long-term spectral changes after the melt event. Single-station observations are inverted for elastic structure for multiple stations across the ice shelf. High-frequency ambient noise seismology presents opportunities for continuous assessment of near-surface ice shelf or other firn environments. 
    more » « less
  3. Abstract

    Explosions and earthquakes are effectively discriminated byP/Samplitude ratios for moderate magnitude events (M ≥ 4) observed at regional to teleseismic distances (≥200 km). It is less clear ifP/Sratios are effective explosion discriminants for lower magnitudes observed at shorter distances. We report new tests ofP/Sdiscrimination using a dense seismic array in a continental volcanic arc setting near Mount St. Helens, with 23 single‐fired borehole explosions (ML0.9–2.3) and 406 earthquakes (ML1–3.3). The array provides up to 95 three‐component broadband seismographs, and most source‐receiver distances are <120 km. Additional insight is provided by ~3,000 vertical component geophone recordings of each explosion. Potential controls on local distanceP/Sratios are investigated, including frequency range, distance, magnitude, source depth, number of seismographs, and site effects. A frequency band of about 10–18 Hz performs better than lower or narrower bands because explosion‐inducedSwave amplitudes diminish relative toPfor higher frequencies. Source depth and magnitude exhibited weak influences onP/Sratios. Site responses for earthquakes and explosions are correlated with each other and with shallow crustalVpandVsfrom traveltime tomography. Overall, the results indicate high potential for local distanceP/Sexplosion discrimination in a continental volcanic arc setting, with ≥98% true positives and ≤6.3% false positives when using the array median from ≥16 stations. Performance is reduced for smaller arrays, especially those with ≤4 stations, thereby emphasizing the importance of array data for discrimination of low magnitude explosions.

     
    more » « less
  4. SUMMARY

    We have simulated 0–5 Hz deterministic wave propagation for a suite of 17 models of the 2014 Mw 5.1 La Habra, CA, earthquake with the Southern California Earthquake Center Community Velocity Model Version S4.26-M01 using a finite-fault source. Strong motion data at 259 sites within a 148 km × 140 km area are used to validate our simulations. Our simulations quantify the effects of statistical distributions of small-scale crustal heterogeneities (SSHs), frequency-dependent attenuation Q(f), surface topography and near-surface low-velocity material (via a 1-D approximation) on the resulting ground motion synthetics. The shear wave quality factor QS(f) is parametrized as QS, 0 and QS, 0fγ for frequencies less than and higher than 1 Hz, respectively. We find the most favourable fit to data for models using ratios of QS, 0 to shear wave velocity VS of 0.075–1.0 and γ values less than 0.6, with the best-fitting amplitude drop-off for the higher frequencies obtained for γ values of 0.2–0.4. Models including topography and a realistic near-surface weathering layer tend to increase peak velocities at mountain peaks and ridges, with a corresponding decrease behind the peaks and ridges in the direction of wave propagation. We find a clear negative correlation between the effects on peak ground velocity amplification and duration lengthening, suggesting that topography redistributes seismic energy from the large-amplitude first arrivals to the adjacent coda waves. A weathering layer with realistic near-surface low velocities is found to enhance the amplification at mountain peaks and ridges, and may partly explain the underprediction of the effects of topography on ground motions found in models. Our models including topography tend to improve the fit to data, as compared to models with a flat free surface, while our distributions of SSHs with constraints from borehole data fail to significantly improve the fit. Accuracy of the velocity model, particularly the near-surface low velocities, as well as the source description, controls the resolution with which the anelastic attenuation can be determined. Our results demonstrate that it is feasible to use fully deterministic physics-based simulations to estimate ground motions for seismic hazard analysis up to 5 Hz. Here, the effects of, and trade-offs with, near-surface low-velocity material, topography, SSHs and Q(f) become increasingly important as frequencies increase towards 5 Hz, and should be included in the calculations. Future improvement in community velocity models, wider access to computational resources, more efficient numerical codes and guidance from this study are bound to further constrain the ground motion models, leading to more accurate seismic hazard analysis.

     
    more » « less
  5. null (Ed.)
    Abstract Ocean swell interacting with Antarctic ice shelves produces sustained (approximately, 2×106 cycles per year) gravity-elastic perturbations with deformation amplitudes near the ice front as large as tens to hundreds of nanostrain. This process is the most energetically excited during the austral summer, when sea ice-induced swell attenuation is at a minimum. A 2014–2017 deployment of broadband seismographs on the Ross Ice shelf, which included three stations sited, approximately, 2 km from the ice front, reveals prolific swell-associated triggering of discrete near-ice-front (magnitude≲0) seismic subevents, for which we identify three generic types. During some strong swell episodes, subevent timing becomes sufficiently phase-locked with swell excitation, to create prominent harmonic features in spectra calculated across sufficiently lengthy time windows via a Dirac comb effect, for which we articulate a theoretical development for randomized interevent times. These events are observable at near-front stations, have dominant frequency content between 0.5 and 20 Hz, and, in many cases, show highly repetitive waveforms. Matched filtering detection and analysis shows that events occur at a low-background rate during all swell states, but become particularly strongly excited during large amplitude swell at rates of up to many thousands per day. The superimposed elastic energy from swell-triggered sources illuminates the shelf interior as extensional (elastic plate) Lamb waves that are observable more than 100 km from the ice edge. Seismic swarms show threshold excitation and hysteresis with respect to rising and falling swell excitation. This behavior is consistent with repeated seismogenic fracture excitation and growth within a near-ice-front damage zone, encompassing fracture features seen in satellite imagery. A much smaller population of distinctly larger near-front seismic events, previously noted to be weakly associated with extended periods of swell perturbation, likely indicate calving or other larger-scale ice failures near the shelf front. 
    more » « less