Search for: All records

Abstract The firn layer covers 98% of Antarctica's ice sheets, protecting underlying glacial ice from the external environment. Accurate measurement of firn properties is essential for assessing cryosphere mass balance and climate change impacts. Characterizing firn structure through core sampling is expensive and logistically challenging. Seismic surveys, which translate seismic velocities into firn densities, offer an efficient alternative. This study employs Distributed Acoustic Sensing technology to transform an existing fiber‐optic cable near the South Pole into a multichannel, low‐maintenance, continuously interrogated seismic array. The data resolve 16 seismic wave propagation modes at frequencies up to 100 Hz that constrain P and S wave velocities as functions of depth. Using co‐located geophones for ambient noise interferometry, we resolve very weak radial anisotropy. Leveraging nearby SPICEcore firn density data, we find prior empirical density‐velocity relationships underestimate firn air content by over 15%. We present a new empirical relationship for the South Pole region. 

The response of the Antarctic Ice Sheet (AIS) to climate change is the largest uncertainty in projecting future sea level. The impact of three-dimensional (3D) Earth structure on the AIS and future global sea levels is assessed here by coupling a global glacial isostatic adjustment model incorporating 3D Earth structure to a dynamic ice-sheet model. We show that including 3D viscous effects produces rapid uplift in marine sectors and reduces projected ice loss for low greenhouse gas emission scenarios, lowering Antarctica’s contribution to global sea level in the coming centuries by up to ~40%. Under high-emission scenarios, ice retreat outpaces uplift, and sea-level rise is amplified by water expulsion from Antarctic marine areas. 

Abstract We develop a 3‐D isotropic shear velocity model for the Alaska subduction zone using data from seafloor and land‐based seismographs to investigate along‐strike variations in structure. By applying ambient noise and teleseismic Helmholtz tomography, we derive Rayleigh wave group and phase velocity dispersion maps, then invert them for shear velocity structure using a Bayesian Monte Carlo algorithm. For land‐based stations, we perform a joint inversion of receiver functions and dispersion curves. The forearc crust is relatively thick (35–42 km) and has reduced lower crustal velocities beneath the Kodiak and Semidi segments, which may promote higher seismic coupling. Bristol Bay Basin crust is relatively thin and has a high‐velocity lower layer, suggesting a dense mafic lower crust emplaced by the rifting processes. The incoming plate shows low uppermost mantle velocities, indicating serpentinization. This hydration is more pronounced in the Shumagin segment, with greater velocity reduction extending to 18 ± 3 km depth, compared to the Semidi segment, showing smaller reductions extending to 14 ± 3 km depth. Our estimates of percent serpentinization from V_Sreduction and V_P/V_Sare larger than those determined using V_Preduction in prior studies, likely due to water in cracks affecting V_Smore than V_P. Revised estimates of serpentinization show that more water subducts than previous studies, and that twice as much mantle water is subducted in the Shumagin segment compared to the Semidi segment. Together with estimates from other subduction zones, the results indicate a wide variation in subducted mantle water between different subduction segments. 

Abstract Oceanic plates experience extensive normal faulting as they bend and subduct, enabling fracturing of the incoming lithosphere. Debate remains about the relative importance of pre‐existing faults, plate curvature and other factors controlling the extent and style of bending‐related faulting. The subduction zone off the Alaska Peninsula is an ideal place to investigate controls on bending faulting as the orientation of the abyssal‐hill fabric with respect to the trench and plate curvature vary along the margin. Here, we characterize faulting between longitudes 161°W and 155°W using newly collected multibeam bathymetry data. We also use a compilation of seismic reflection data to constrain patterns of sediment thickness on the incoming plate. Although sediment thickness increases over 1 km from 156°W to 160°W, most sediments were deposited prior to the onset of bending faulting and thus should have limited impact on the expression of bend‐related fault strikes and throws in bathymetry data. Where magnetic anomalies trend subparallel to the trench (<30°) west of ∼156°W, bending faults parallel magnetic anomalies, implying that bending faults reactivate pre‐existing structures. Where magnetic anomalies are highly oblique (>30°) to the trench east of 156°W, no bending faults are observed. Summed fault throws increase to the west, including where pre‐existing structure orientations are constant (between 157 and 161°W), suggesting that another factor such as the increase in slab curvature must influence bending faulting. However, the westward increase in summed fault throws is more abrupt than expected for gradual changes in slab bending alone, suggesting potential feedbacks between pre‐existing structures, slab dip, and faulting. 

Abstract Seismic anisotropy provides essential information for characterizing the orientation of deformation and flow in the crust and mantle. The isotropic structure of the Antarctic crust and upper mantle has been determined by previous studies, but the azimuthal anisotropy structure has only been constrained by mantle core phase (SKS) splitting observations. This study determines the azimuthal anisotropic structure of the crust and mantle beneath the central and West Antarctica based on 8—55 s Rayleigh wave phase velocities from ambient noise cross‐correlation. An anisotropic Rayleigh wave phase velocity map was created using a ray—based tomography method. These data are inverted using a Bayesian Monte Carlo method to obtain an azimuthal anisotropy model with uncertainties. The azimuthal anisotropy structure in most of the study region can be fit by a two‐layer structure, with one layer at depths of 0–15 km in the shallow crust and the other layer in the uppermost mantle. The azimuthal anisotropic layer in the shallow crust of West Antarctica, where it coincides with strong positive radial anisotropy quantified by the previous study, shows a fast direction that is subparallel to the inferred extension direction of the West Antarctic Rift System. Fast directions of upper mantle azimuthal anisotropy generally align with teleseismic shear wave splitting fast directions, suggesting a thin lithosphere or similar lithosphere‐asthenosphere deformation. However, inconsistencies in this exist in Marie Byrd Land, indicating differing ancient deformation patterns in the shallow mantle lithosphere sampled by the surface waves and deformation in the deeper mantle and asthenosphere sampled more strongly by splitting measurements. 

Abstract The deployment of seismic stations and the development of ambient noise tomography and new analysis methods provide an opportunity for higher resolution imaging of Antarctica. Here we review recent seismic structure models and describe their implications for the dynamics and history of the Antarctic upper mantle. Results show that most of East Antarctica is underlain by continental lithosphere to depths of ∼ 200 km. The thickest lithosphere is found in a band 500-1000 km west of the Transantarctic Mountains, representing the continuation of cratonic lithosphere with Australian affinity beneath the ice. Dronning Maud Land and the Lambert Graben show much thinner lithosphere, consistent with Phanerozoic lithospheric disruption. The Transantarctic Mountains mark a sharp boundary between cratonic lithosphere and the warmer upper mantle of West Antarctica. In the Southern Transantarctic Mountains, cratonic lithosphere has been replaced by warm asthenosphere, giving rise to Cenozoic volcanism and an elevated mountainous region. The Marie Byrd Land volcanic dome is underlain by slow seismic velocities extending through the transition zone, consistent with a mantle plume. Slow velocity anomalies beneath the coast from the Amundsen Sea Embayment to the Antarctic Peninsula likely result from upwelling of warm asthenosphere during subduction of the Antarctic-Phoenix spreading center. 

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. 

Abstract Observations of teleseismic earthquakes using broadband seismometers on the Ross Ice Shelf (RIS) must contend with environmental and structural processes that do not exist for land-sited seismometers. Important considerations are: (1) a broadband, multi-mode ambient wavefield excited by ocean gravity wave interactions with the ice shelf; (2) body wave reverberations produced by seismic impedance contrasts at the ice/water and water/seafloor interfaces and (3) decoupling of the solid Earth horizontal wavefield by the sub-shelf water column. We analyze seasonal and geographic variations in signal-to-noise ratios for teleseismic P-wave (0.5–2.0 s), S-wave (10–15 s) and surface wave (13–25 s) arrivals relative to the RIS noise field. We use ice and water layer reverberations generated by teleseismic P-waves to accurately estimate the sub-station thicknesses of these layers. We present observations consistent with the theoretically predicted transition of the water column from compressible to incompressible mechanics, relevant for vertically incident solid Earth waves with periods longer than 3 s. Finally, we observe symmetric-mode Lamb waves generated by teleseismic S-waves incident on the grounding zones. Despite their complexity, we conclude that teleseismic coda can be utilized for passive imaging of sub-shelf Earth structure, although longer deployments relative to conventional land-sited seismometers will be necessary to acquire adequate data.