Search for: All records

Abstract Relative seismic velocity changes are being increasingly used to monitor changes in groundwater. However, it remains challenging to verify its implementation in watersheds without direct groundwater well measurements. In this study, we conduct a 12‐year observation in a watershed of the Yellowstone National Park (YNP). We find that the seasonal fluctuations and long‐term trend of the measured are highly correlated with the estimated baseflow, which serves as a constraint for groundwater changes. We integrate the estimated baseflow into a poroelastic mechanism and conduct two simulations based on pressure diffusion. These simulations closely match with our observed variations. In addition, our analysis suggest that the measured is primarily influenced by hydrologic pressure diffusion rather than surface air temperature. We conclude that the baseflow analysis can further enhance the seismic monitoring of groundwater changes. 

Abstract To explore seismic structures beneath the Australian continents and subduction zone geometry around the Australian plate, we introduce a new radially‐anisotropic shear‐wavespeed model, AU21. By employing full‐waveform inversion on data from 248 regional earthquakes and 1,102 seismographic stations, we iteratively refine AU21, resulting in 32,655 body‐wave and 35,897 surface wave measurements. AU21 reveals distinct shear‐wavespeed contrasts between the Phanerozoic eastern continental margin and the Precambrian western and central Australia, with the lithosphere‐asthenosphere boundary estimated at 250–300 km beneath central and western Australia. Notably, a unique weak radial anisotropy layer at 80–150 km is identified beneath the western Australian craton, possibly due to alignments of dipping layers or tilted symmetry axes of anisotropic minerals. Furthermore, slow anomalies extending to the uppermost lower mantle beneath the east of New Guinea, Tasmania, and the Tasman Sea indicate deep thermal activities, likely contributing to the formation of a low wavespeed band along the eastern Australian margin. In addition, our findings demonstrate the stagnant Tonga slab within the mantle transition zone and the Kermadec slab's penetration through the 660‐km discontinuity into the lower mantle. 

Abstract The architecture of magma plumbing systems plays a fundamental role in volcano eruption and evolution. However, the precise configuration of crustal magma reservoirs and conduits responsible for supplying eruptions are difficult to explore across most active volcanic systems. Consequently, our understanding of their correlation with eruption dynamics is limited. Axial Seamount is an active submarine volcano located along the Juan de Fuca Ridge, with known eruptions in 1998, 2011, and 2015. Here we present high-resolution images of P-wave velocity, attenuation, and estimates of temperature and partial melt beneath the summit of Axial Seamount, derived from multi-parameter full waveform inversion of a 2D multi-channel seismic line. Multiple magma reservoirs, including a newly discovered western magma reservoir, are identified in the upper crust, with the maximum melt fraction of ~15–32% in the upper main magma reservoir (MMR) and lower fractions of 10% to 26% in other satellite reservoirs. In addition, a feeding conduit below the MMR with a melt fraction of ~4–11% and a low-velocity throat beneath the eastern caldera wall connecting the MMR roof with eruptive fissures are imaged. These findings delineate an asymmetric shallow plumbing system beneath Axial Seamount, providing insights into the magma pathways that fed recent eruptions. 

SUMMARY For seismographic stations with short acquisition duration, the signal-to-noise ratios (SNRs) of ambient noise cross-correlation functions (CCFs) are typically low, preventing us from accurately measuring surface wave dispersion curves or waveform characteristics. In addition, with noisy CCFs, it is difficult to extract relatively weak signals such as body waves. In this study, we propose to use local attributes to improve the SNRs of ambient noise CCFs, which allows us to enhance the quality of CCFs for stations with limited acquisition duration. Two local attributes: local cross-correlation and local similarity, are used in this study. The local cross-correlation allows us to extend the dimensionality of daily CCFs with computational costs similar to global cross-correlation. Taking advantage of this extended dimensionality, the local similarity is then used to measure non-stationary similarity between the extended daily CCFs with a reference trace, which enables us to design better stacking weights to enhance coherent features and attenuate incoherent background noises. Ambient noise recorded by several broad-band stations from the USArray in North Texas and Oklahoma, the Superior Province Rifting EarthScope Experiment in Minnesota and Wisconsin and a high-frequency nodal array deployed in the northern Los Angeles basin are used to demonstrate the performance of the proposed approach for improving the SNR of CCFs. 

Abstract Increasing deployment of dense arrays has facilitated detailed structure imaging for tectonic investigation, hazard assessment and resource exploration. Strong velocity heterogeneity and topographic changes have to be considered during passive source imaging. However, it is quite challenging for ray‐based methods, such as Kirchhoff migration or the widely used teleseismic receiver function, to handle these problems. In this study, we propose a 3‐D passive source reverse time migration strategy based on the spectral element method. It is realized by decomposing the time reversal full elastic wavefield into amplitude‐preserved vector P and S wavefields by solving the corresponding weak‐form solutions, followed by a dot‐product imaging condition to get images for the subsurface structures. It enables us to use regional 3‐D migration velocity models and take topographic variations into account, helping us to locate reflectors at more accurate positions than traditional 1‐D model‐based methods, like teleseismic receiver functions. Two synthetic tests are used to demonstrate the advantages of the proposed method to handle topographic variations and complex velocity heterogeneities. Furthermore, applications to the Laramie array data using both teleseismic P and S waves enable us to identify several south‐dipping structures beneath the Laramie basin in southeast Wyoming, which are interpreted as the Cheyenne Belt suture zone and agree with, and improve upon previous geological interpretations. 

Abstract Over the past decade, the seismicity rate in the state of Oklahoma has increased significantly, which has been linked to industrial operations, such as saltwater injection and hydraulic fracturing. Taking advantage of induced earthquakes and recently deployed seismometers, we construct a 3‐D radially anisotropic seismic velocity model for the crust of Oklahoma by using full waveform inversion. To mitigate the well‐known cycle‐skipping problem, we use misfit functions based on phase and waveform differences in several frequency bands. Relative velocity perturbations in the inverted model allow us to delineate major geological provinces in Oklahoma, such as the Anadarko Basin and the Cherokee Platform/Shelf. In addition, radial anisotropy in the inverted model reflects deformation within the crust of Oklahoma, which might correlate with sedimentary layering, microcracks/fractures, as well as dominant orientation of anisotropic minerals. The crystalline basement beneath Oklahoma can be inferred from the new velocity model, which enables us to better classify induced seismicity in current earthquake catalogs. Furthermore, synthetic experiments suggest that the new velocity model enables us to better constrain earthquake locations in Oklahoma, especially for determining their depths, which are important for investigating induced seismicity. 

Abstract Significant imbalances in terrestrial water storage (TWS) and severe drought have been observed around the world as a consequence of climate changes. Improving our ability to monitor TWS and drought is critical for water‐resource management and water‐deficit estimation. We use continuous seismic ambient noise to monitor temporal evolution of near‐surface seismic velocity,dv/v, in central Oklahoma from 2013 to 2022. The deriveddv/vis found to be negatively correlated with gravitational measurements and groundwater depths, showing the impact of groundwater storage on seismic velocities. The hydrological effects involving droughts and recharge of groundwater occur on a multi‐year time scale and dominate the overall derived velocity changes. The thermoelastic response to atmospheric temperature variations occurs primarily on a yearly timescale and dominates the superposed seasonal velocity changes in this study. The occurrences of droughts appear simultaneously with local peaks ofdv/v, demonstrating the sensitivity of near‐surface seismic velocities to droughts. 

Abstract One important feature of the Greenland Ice Sheet (GrIS) change is its strong seasonal fluctuation. Taking advantage of deployed seismographic stations in Greenland, we apply cross‐component auto‐correlation of seismic ambient noise to measure in‐situ near surface relative velocity change (dv/v) in different regions of Greenland. Our results demonstrate thatdv/vmeasurements for most stations have less than 3 months lag times in comparison to the surface mass change. These various lag times may provide us constraints for the thickness of the subglacial till layer over different regions in Greenland. Moreover, in southwest Greenland, we observe a change in the long‐term trend ofdv/vfor three stations, which might be consistent with the mass change rate (dM/dt) due to the “2012–2013 warm‐cold transition.” These observations suggest that seismic noise auto‐correlation technique may be used to monitor both seasonal and long‐term changes of the GrIS. 

Abstract By fitting observed data with predicted seismograms, least‐squares migration (LSM) computes a generalized inverse for a subsurface reflectivity model, which can improve image resolution and reduce artifacts caused by incomplete acquisition. However, the large computational cost of LSM required for simulations and migrations limits its wide applications for large‐scale imaging problems. Using point‐spread function (PSF) deconvolution, we present an efficient and stable high‐resolution imaging method. The PSFs are first computed on a coarse grid using local ray‐based Gaussian beam Born modeling and migration. Then, we interpolate the PSFs onto a fine‐image grid and apply a high‐dimensional Gaussian function to attenuate artifacts far away from the PSF centers. With 2D/3D partition of unity, we decompose the traditional adjoint migration results into local images with the same window size as the PSFs. Then, these local images are deconvolved by the PSFs in the wavenumber domain to reduce the effects of the band‐limited source function and compensate for irregular subsurface illumination. The final assembled image is obtained by applying the inverse of the partitions for the deconvolved local images. Numerical examples for both synthetic and field data demonstrate that the proposed PSF deconvolution can significantly improve image resolution and amplitudes for deep structures, while not being sensitive to velocity errors as the data‐domain LSM. 

Abstract Laboratory experiments and geodynamic simulations demonstrate that poloidal- and toroidal-mode mantle flows develop around subduction zones. Here, we use a new 3-D azimuthal anisotropy model constructed by full waveform inversion, to infer deep subduction-induced mantle flows underneath Middle America. At depths shallower than 150 km, poloidal-mode flow is perpendicular to the trajectory of the Middle American Trench. From 300 to 450 km depth, return flows surround the edges of the Rivera and Atlantic slabs, while escape flows are inferred through slab windows beneath Panama and central Mexico. Furthermore, at 700 km depth, the study region is dominated by the Farallon anomaly, with fast axes perpendicular to its strike, suggesting the development of lattice-preferred orientations by substantial stress. These observations provide depth-dependent seismic anisotropy for future mantle flow simulations, and call for further investigations about the deformation mechanisms and elasticity of minerals in the transition zone and uppermost lower mantle.