An official website of the United States government
Abstract The Amundsen Sea Embayment of the West Antarctic Ice Sheet contains Thwaites and Pine Island Glaciers, two of the most rapidly changing glaciers in Antarctica. To date, Pine Island and Thwaites Glaciers have only been observed by independent airborne radar sounding surveys, but a combined cross‐basin analysis that investigates the basal conditions across the Pine Island‐Thwaites Glaciers boundary has not been performed. Here, we combine two radar surveys and correct for their differences in system parameters to produce unified englacial attenuation and basal relative reflectivity maps spanning both Pine Island and Thwaites Glaciers. Relative reflectivities range from −24.8 to +37.4 dB with the highest values beneath fast‐flowing ice at the ice sheet margin. By comparing our reflectivity results with previously derived radar specularity and trailing bed echoes at Thwaites Glacier, we find a highly diverse subglacial landscape and hydrologic conditions that evolve along‐flow. Together, these findings highlight the potential for joint airborne radar analysis with ground‐based seismic and geomorphological observations to understand variations in the bed properties and cross‐catchment interactions of ice streams and outlet glaciers.
more »
« less
Signal characteristics of surface seismic explosive sources near the West Antarctic Ice Sheet divide
Abstract Seismic imaging in 3-D holds great potential for improving our understanding of ice sheet structure and dynamics. Conducting 3-D imaging in remote areas is simplified by using lightweight and logistically straightforward sources. We report results from controlled seismic source tests carried out near the West Antarctic Ice Sheet Divide investigating the characteristics of two types of surface seismic sources, Poulter shots and detonating cord, for use in both 2-D and 3-D seismic surveys on glaciers. Both source types produced strong basal P-wave and S-wave reflections and multiples recorded in three components. The Poulter shots had a higher amplitude for low frequencies (<10 Hz) and comparable amplitude at high frequencies (>50 Hz) relative to the detonating cord. Amplitudes, frequencies, speed of source set-up, and cost all suggested Poulter shots to be the preferred surface source compared to detonating cord for future 2-D and 3-D seismic surveys on glaciers.
more »
« less
- Award ID(s):
- 1739027
- PAR ID:
- 10579312
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Glaciology
- Volume:
- 70
- ISSN:
- 0022-1430
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
SUMMARY Elastodynamic Green’s functions are an essential ingredient in seismology as they form the connection between direct observations of seismic waves and the earthquake source. They are also fundamental to various seismological techniques including physics-based ground motion prediction and kinematic or dynamic source inversions. In regions with established 3-D models of the Earth’s elastic structure, such as southern California, 3-D Green’s functions can be computed using numerical simulations of seismic wave propagation. However, such simulations are computationally expensive, which poses challenges for real-time ground motion prediction and uncertainty quantification in source inversions. In this study, we address these challenges by using a reduced-order model (ROM) approach that enables the rapid evaluation of approximate Green’s functions. The ROM technique developed approximates three-component time-dependent surface velocity wavefields obtained from numerical simulations of seismic wave propagation. We apply our ROM approach to a 50 km $$\times$$ 40 km area in greater Los Angeles accounting for topography, site effects, 3-D subsurface velocity structure, and viscoelastic attenuation. The ROM constructed for this region enables rapid computation ($$\approx 0.0001$$ CPU hr) of complete, high-resolution (500 m spacing), 0.5 Hz surface velocity wavefields that are accurate for a shortest wavelength of 1.0 km for a single elementary moment tensor source. Using leave-one-out cross validation, we measure the accuracy of our Green’s functions for the CVM-S velocity model in both the time domain and frequency domain. Averaged across all sources, receivers, and time steps, the error in the rapid seismograms is less than 0.01 cm s−1. We demonstrate that the ROM can accurately and rapidly reproduce simulated seismograms for generalized moment tensor sources in our region, as well as kinematic sources by using a finite fault model of the 1987 $$M_\mathrm{ W}$$ 5.9 Whittier Narrows earthquake as an example. We envision that rapid, accurate Green’s functions from reduced-order modelling for complex 3-D seismic wave propagation simulations will be useful for constructing real-time ground motion synthetics and source inversions with high spatial resolution.more » « less
-
Abstract Debris‐covered glaciers (DCG) and rock glaciers have been increasingly studied in recent years because of the role they play within local watersheds, glacial ablation models due to climate change, and as analogs for buried ice features on planetary bodies such as Mars. Characterizing the supraglacial debris layer is a large part of these efforts. Geophysical exploration of DCG has mostly excluded active seismic methods, with the exception of refraction studies, due to the attenuating properties of the debris cover and field survey efficiency. We evaluate the accuracy, field efficiency, and effectiveness of seismic refraction, reflection, and surface‐wave surveys for determining the elastic properties of the debris layer and any underlying layers on DCG using the Sourdough Rock Glacier in Southcentral Alaska as a test site. We provide evidence for imaging an ultra‐shallow seismic reflection from the bottom of the loose debris layer using ultra‐dense receiver arrays and compare it to ground‐penetrating radar (GPR) images taken along the same profiles. We also detail how reliable dispersion curve images can be extracted from the surface wave package of the seismic data using the multi‐channel analysis of surface waves technique, which allows for the (s)‐wave profile to be inverted for. We find this could be a valuable addition to the toolbox of future geophysical investigations on DCG.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jog.2024.41
