Search for: All records

SUMMARY The sensitivity of Rayleigh wave amplitude to Earth structure has applications to seismic tomography, both in cases where amplitude information is used to supplement phase velocity data to improve images of elastic parameters, and to correct amplitudes for local Earth structure in attenuation tomography. We review the theoretical basis of the ray theoretical approximation, in which the wave amplitudes are controlled by a combination of geometrical spreading and local changes in energy density due to Earth structure. We focus mainly on the latter effect, which we term the constant energy flux approximation. We investigate the ray theoretical basis for this approximation, test it against a full waveform simulation that verifies its accuracy and show how it can be used to compute the sensitivity of amplitude to elastic moduli and density. We investigate how perturbing these parameters in a set of simple Earth models affects Rayleigh wave amplitudes, and demonstrate that a slow velocity heterogeneity can cause either increased or reduced amplitudes, depending upon the depth of the heterogeneity and the observation frequency. Consequently, amplitude sensitivity can be either positive or negative, and its magnitude can vary significantly with frequency. Although an added complication, the very different behaviour of phase velocity and amplitudes to changes in Earth structure implies that the two types of data are complementary and suggest the effectiveness of using both in Rayleigh wave tomography. 

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 On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical Extreme Events Reconnaissance Association and the National Institute of Standards and Technology documented structural and nonstructural damage, including substantial damage to a historic masonry building in Lebanon, New Jersey. The USGS National Earthquake Information Center reported a focal depth of about 5 km, consistent with a lack of signal in Interferometric Synthetic Aperture Radar data. The focal mechanism solution is strike slip with a substantial thrust component. Neither mechanism’s nodal plane is parallel to the primary northeast trend of geologic discontinuities and mapped faults in the region, including the Ramapo fault. However, many of the relocated aftershocks, for which locations were augmented by temporary seismic deployments, form a cluster that parallels the general northeast trend of the faults. The aftershocks lie near the Tewksbury fault, north of the Ramapo fault. 

SUMMARY A key initial step in geophysical imaging is to devise an effective means of mapping the sensitivity of an observation to the model parameters, that is to compute its Fréchet derivatives or sensitivity kernel. In the absence of any simplifying assumptions and when faced with a large number of free parameters, the adjoint method can be an effective and efficient approach to calculating Fréchet derivatives and requires just two numerical simulations. In the Glacial Isostatic Adjustment problem, these consist of a forward simulation driven by changes in ice mass and an adjoint simulation driven by fictitious loads that are applied at the observation sites. The theoretical basis for this approach has seen considerable development over the last decade. Here, we present the final elements needed to image 3-D mantle viscosity using a dataset of palaeo sea-level observations. Developments include the calculation of viscosity Fréchet derivatives (i.e. sensitivity kernels) for relative sea-level observations, a modification to the numerical implementation of the forward and adjoint problem that permits application to 3-D viscosity structure, and a recalibration of initial sea level that ensures the forward simulation honours present-day topography. In the process of addressing these items, we build intuition concerning how absolute sea-level and relative sea-level observations sense Earth’s viscosity structure and the physical processes involved. We discuss examples for potential observations located in the near field (Andenes, Norway), far field (Seychelles), and edge of the forebulge of the Laurentide ice sheet (Barbados). Examination of these kernels: (1) reveals why 1-D estimates of mantle viscosity from far-field relative sea-level observations can be biased; (2) hints at why an appropriate differential relative sea-level observation can provide a better constraint on local mantle viscosity and (3) demonstrates that sea-level observations have non-negligible 3-D sensitivity to deep mantle viscosity structure, which is counter to the intuition gained from 1-D radial viscosity Fréchet derivatives. Finally, we explore the influence of lateral variations in viscosity on relative sea-level observations in the Amundsen Sea Embayment and at Barbados. These predictions are based on a new global 3-D viscosity inference derived from the shear-wave speeds of GLAD-M25 and an inverse calibration scheme that ensures compatibility with certain fundamental geophysical observations. Use of the 3-D viscosity inference leads to: (1) generally greater complexity within the kernel; (2) an increase in sensitivity and presence of shorter length-scale features within lower viscosity regions; (3) a zeroing out of the sensitivity kernel within high-viscosity regions where elastic deformation dominates and (4) shifting of sensitivity at a given depth towards distal regions of weaker viscosity. The tools and intuition built here provide the necessary framework to explore inversions for 3-D mantle viscosity based on palaeo sea-level data. 

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 The geothermal heat flux (GHF) is an important boundary condition for modeling the movement of the Antarctic ice sheet but is difficult to measure systematically at a continental scale. Earlier GHF maps suffer from low resolution and possibly biased assumptions in tectonism and crustal heat generation, resulting in significant uncertainty. We present a new GHF map for Antarctica constructed by empirically relating the upper mantle structure to known GHF in the continental United States. The new map, compared with previously seismologically determined one, has improved resolution and lower uncertainties. New features in this map include high GHF in the southern Transantarctic Mountains where warmer uppermost mantle is introduced by lithospheric removal and in the Thwaites Glacier region. Additionally, a modest GHF in the central West Antarctic Rift system near the Siple Coast and an absence of large‐scale regions with GHF greater than 90 mW/m²are found. 

Abstract We examine upper mantle anisotropy across the Antarctic continent using 102 new shear wave splitting measurements obtained from teleseismic SKS, SKKS, and PKS phases combined with 107 previously published results. For the new measurements, an eigenvalue technique is used to estimate the fast polarization direction and delay time for each phase arrival, and high‐quality measurements are stacked to determine the best‐fit splitting parameters at each seismic station. The ensemble of splitting measurements shows largely NE‐SW‐oriented fast polarization directions across Antarctica, with a broadly clockwise rotation in polarization directions evident moving from west to east across the continent. Although the first‐order pattern of NE‐SW‐oriented polarization directions is suggestive of a single plate‐wide source of anisotropy, we argue the observed pattern of anisotropy more likely arises from regionally variable contributions of both lithospheric and sub‐lithospheric mantle sources. Anisotropy observed in the interior of East Antarctica, a region underlain by thick lithosphere, can be attributed to relict fabrics associated with Precambrian tectonism. In contrast, anisotropy observed in coastal East Antarctica, the Transantarctic Mountains (TAM), and across much of West Antarctica likely reflects both lithospheric and sub‐lithospheric mantle fabrics. While sub‐lithospheric mantle fabrics are best associated with either plate motion‐induced asthenospheric flow or small‐scale convection, lithospheric mantle fabrics in coastal East Antarctica, the TAM, and West Antarctica generally reflect Jurassic—Cenozoic tectonic activity. 

Abstract We have located 117 previously undetected seismic events mainly occurring between 2015 and 2017 that originated from glacial, tectonic, and volcanic processes in central West Antarctica using data recorded on Polar Earth Observing Network (POLENET/ANET) and UK Antarctic Network (UKANET) seismic stations. The seismic events, with local magnitudes (M_L) ranging from 1.1 to 3.5, are predominantly clustered in four geographic regions; the Ellsworth Mountains, Thwaites Glacier, Pine Island Glacier, and Mount Takahe. Eighteen of the events are in the Ellsworth Mountains and can be attributed to a mixture of glacial and tectonic processes. The largest event noted in this study was a mid‐crustal (∼19 km focal depth;M_L3.5) normal mechanism earthquake beneath Thwaites Glacier. We also located 91 glacial events near the grounding zones of Thwaites Glacier and Pine Island Glacier that are predominantly associated with time periods of significant calving activity. Eight events, likely arising from volcano‐tectonic processes, occurred beneath Mount Takahe. Using Pn travel times from the seismic events, we find laterally variable uppermost mantle structure in central West Antarctica. On average, the Ellsworth Mountains are underlain by a faster mantle lid (V_Pn = ∼8.4 km/s) compared to the Amundsen Sea Embayment region (V_Pn = ∼8.1 km/s). Within the Amundsen Sea Embayment itself, we find mantle lid velocities ranging from ∼8.05 to 8.18 km/s. Laterally heterogeneous uppermost mantle structure, indicative of variable thermal and rheological structure, likely influences both geothermal heat flux and glacial isostatic adjustment spatial patterns and rates within central West Antarctica.