An official website of the United States government
Fractures in Earth's critical zone influence groundwater flow and storage and promote chemical weathering. Fractured materials are difficult to characterize on large spatial scales because they contain fractures that span a range of sizes, have complex spatial distributions, and are often inaccessible. Therefore, geophysical characterizations of the critical zone depend on the scale of measurements and on the response of the medium to impulses at that scale. Using P-wave velocities collected at two scales, we show that seismic velocities in the fractured bedrock layer of the critical zone are scale-dependent. The smaller-scale velocities, derived from sonic logs with a dominant wavelength of ~0.3 m, show substantial vertical and lateral heterogeneity in the fractured rock, with sonic velocities varying by 2,000 m/s over short lateral distances (~20 m), indicating strong spatial variations in fracture density. In contrast, the larger-scale velocities, derived from seismic refraction surveys with a dominant wavelength of ~50 m, are notably slower than the sonic velocities (a difference of ~3,000 m/s) and lack lateral heterogeneity. We show that this discrepancy is a consequence of contrasting measurement scales between the two methods; in other words, the contrast is not an artifact but rather information—the signature of a fractured medium (weathered/fractured bedrock) when probed at vastly different scales. We explore the sample volumes of each measurement and show that surface refraction velocities provide reliable estimates of critical zone thickness but are relatively insensitive to lateral changes in fracture density at scales of a few tens of meters. At depth, converging refraction and sonic velocities likely indicate the top of unweathered bedrock, indicative of material with similar fracture density across scales.
more »
« less
2D Near‐Surface Full‐Waveform Tomography Reveals Bedrock Controls on Critical Zone Architecture
Abstract For decades, seismic imaging methods have been used to study the critical zone, Earth's thin, life‐supporting skin. The vast majority of critical zone seismic studies use traveltime tomography, which poorly resolves heterogeneity at many scales relevant to near‐surface processes, therefore limiting progress in critical zone science. Full‐waveform tomography can overcome this limitation by leveraging more seismic data and enhancing the resolution of geophysical imaging. In this study, we apply 2D full‐waveform tomography to match the phases of observed seismograms and elucidate previously undetected heterogeneity in the critical zone at a well‐studied catchment in the Laramie Range, Wyoming. In contrast to traveltime tomograms from the same data set, our results show variations in depth to bedrock ranging from 5 to 60 m over lateral scales of just tens of meters and image steep low‐velocity anomalies suggesting hydrologic pathways into the deep critical zone. Our results also show that areas with thick fractured bedrock layers correspond to zones of slightly lower velocities in the deep bedrock, while zones of high bedrock velocity correspond to sharp vertical transitions from bedrock to saprolite. By corroborating these findings with borehole imagery, we hypothesize that lateral changes in bedrock fracture density majorly impact critical zone architecture. Borehole data also show that our full‐waveform tomography results agree significantly better with velocity logs than previously published traveltime tomography models. Full‐waveform tomography thus appears unprecedentedly capable of imaging the spatially complex porosity structure crucial to critical zone hydrology and processes.
more »
« less
- Award ID(s):
- 2012353
- PAR ID:
- 10493243
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Earth and Space Science
- Volume:
- 11
- Issue:
- 2
- ISSN:
- 2333-5084
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Sinkhole collapse may result in significant property damage and even loss of life. Early detection of sinkhole attributes (buried voids, raveling zones) is critical to limit the cost of remediation. One of the most promising ways to obtain subsurface imaging is 3D seismic full-waveform inversion. For demonstration, a recently developed 3D Gauss-Newton full-waveform inversion (3D GN-FWI) method is used to detect buried voids, raveling soils, and characterize variable subsurface soil/rock layering. It is based on a finite-difference solution of 3D elastic wave equations and Gauss-Newton optimization. The method is tested first on a data set constructed from the numerical simulation of a challenging synthetic model and subsequently on field data collected from two separate test sites in Florida. For the field tests, receivers and sources are placed in uniform 2D surface grids to acquire the seismic wavefields, which then are inverted to extract the 3D subsurface velocity structures. The inverted synthetic results suggest that the approach is viable for detecting voids and characterizing layering. The field seismic results reveal that the 3D waveform analysis identified a known manmade void (plastic culvert), unknown natural voids, raveling, as well as laterally variable soil/rock layering including rock pinnacles. The results are confirmed later by standard penetration tests, including depth to bedrock, two buried voids, and a raveling soil zone. Our study provides insight into the application of the 3D seismic FWI technique as a powerful tool in detecting shallow voids and other localized subsurface features.more » « less
-
SUMMARY Long-period (T > 10 s) shear wave reflections between the surface and reflecting boundaries below seismic stations are useful for studying phase transitions in the mantle transition zone (MTZ) but shear-velocity heterogeneity and finite-frequency effects complicate the interpretation of waveform stacks. We follow up on a recent study by Shearer & Buehler (hereafter SB19) of the top-side shear wave reflection Ssds as a probe for mapping the depths of the 410-km and 660-km discontinuities beneath the USArray. Like SB19, we observe that the recorded Ss410s-S and Ss660s-S traveltime differences are longer at stations in the western United States than in the central-eastern United States. The 410-km and 660-km discontinuities are about 40–50 km deeper beneath the western United States than the central-eastern United States if Ss410s-S and Ss660s-S traveltime differences are transformed to depth using a common-reflection point (CRP) mapping approach based on a 1-D seismic model (PREM in our case). However, the east-to-west deepening of the MTZ disappears in the CRP image if we account for 3-D shear wave velocity variations in the mantle according to global tomography. In addition, from spectral-element method synthetics, we find that ray theory overpredicts the traveltime delays of the reverberations. Undulations of the 410-km and 660-km discontinuities are underestimated when their wavelengths are smaller than the Fresnel zones of the wave reverberations in the MTZ. Therefore, modelling of layering in the upper mantle must be based on 3-D reference structures and accurate calculations of reverberation traveltimes.more » « less
-
Abstract Within Earth's critical zone, weathering processes influence landscape evolution and hillslope hydrology by creating porosity in bedrock, transforming it into saprolite and eventually soil. In situ weathering processes drive much of this transformation while preserving the rock fabric of the parent material. Inherited rock fabric in regolith makes the critical zone anisotropic, affecting its mechanical and hydrological properties. Therefore, quantifying and studying anisotropy is an important part of characterising the critical zone, yet doing so remains challenging. Seismic methods can be used to detect rock fabric and infer mechanical and hydrologic conductivity anisotropy across landscapes. We present a novel way of measuring seismic anisotropy in the critical zone using Rayleigh and Love surface waves. This method leverages multi‐component surface seismic data to create a high‐resolution model of seismic anisotropy, which we compare with a nuclear magnetic resonance log measured in a nearby borehole. The two geophysical data sets show that seismic anisotropy and porosity develop at similar depths in weathered bedrock and both reach their maximum values in saprolite, implying that in situ weathering enhances anisotropy while concurrently generating porosity in the critical zone. We bolster our findings with in situ measurements of seismic and hydrologic conductivity anisotropy made in a 3 m deep soil excavation. Our study offers a fresh perspective on the importance of rock fabric in the development and function of the critical zone and sheds new insights into how weathering processes operate.more » « less
-
Despite progress in tomographic imaging of Earth’s interior, a number of critical questions regarding the large-scale structure and dynamics of the mantle remain outstanding. One of those questions is the impact of phase-boundary undulations on global imaging of mantle heterogeneity and on geodynamic (i.e. convection-related) observables. To address this issue, we developed a joint seismic-geodynamic-mineral physical tomographic inversion procedure that incorporates lateral variations in the depths of the 410- and 660-km discontinuities. This inversion includes S-wave traveltimes, SS precursors that are sensitive to transition-zone topography, geodynamic observables/data (free-air gravity, dynamic surface topography, horizontal divergence of tectonic plates and excess core-mantle boundary ellipticity) and mineral physical constraints on thermal heterogeneity. Compared to joint tomography models that do not include data sensitivity to phase-boundary undulations in the transition zone, the inclusion of 410- and 660-km topography strongly influences the inference of volumetric anomalies in a depth interval that encompasses the transition zone and mid-mantle. It is notable that joint tomography inversions, which include constraints on transition-zone discontinuity topography by seismic and geodynamic data, yield more pronounced density anomalies associated with subduction zones and hotspots. We also find that the inclusion of 410- and 660-km topography may improve the fit to the geodynamic observables, depending on the weights applied to seismic and geodynamic data in the inversions. As a consequence, we find that the amplitude of non-thermal density anomalies required to explain the geodynamic data decreases in most of the mantle. These findings underline the sensitivity of the joint inversions to the inclusion of transition-zone complexity (e.g. phase-boundary topography) and the implications for the inferred non-thermal density anomalies in these depth regions. Finally, we underline that our inferences of 410- and 660-km topography avoid a commonly employed approximation that represents the contribution of volumetric heterogeneity to SS-wave precursor data. Our results suggest that this previously employed correction, based on a priori estimates of uppermantle heterogeneity, might be a significant source of error in estimating the 410- and 660-km topography.more » « less
