An official website of the United States government
Abstract The structure of fault zones and the ruptures they host are inextricably linked. Fault zones are narrow, which has made imaging their structure at seismogenic depths a persistent problem. Fiber‐optic seismology allows for low‐maintenance, long‐term deployments of dense seismic arrays, which present new opportunities to address this problem. We use a fiber array that crosses the Garlock Fault to explore its structure. With a multifaceted imaging approach, we peel back the shallow structure around the fault to see how the fault changes with depth in the crust. We first generate a shallow velocity model across the fault with a joint inversion of active source and ambient noise data. Subsequently, we investigate the fault at deeper depths using travel‐time observations from local earthquakes. By comparing the shallow velocity model and the earthquake travel‐time observations, we find that the fault's low‐velocity zone below the top few hundred meters is at most unexpectedly narrow, potentially indicating fault zone healing. Using differential travel‐time measurements from earthquake pairs, we resolve a sharp bimaterial contrast at depth that suggests preferred westward rupture directivity.
more »
« less
Bayesian Seismic Refraction Inversion for Critical Zone Science and Near‐Surface Applications
Abstract The critical zone (CZ) is the region of the Earth’s surface that extends from the bottom of the weathered bedrock to the tree canopy and is important because of its ability to store water and support ecosystems. A growing number of studies use active source shallow seismic refraction to explore and define the size and structure of the CZ across landscapes. However, measurement uncertainty and model resolution at depth are generally not evaluated, which makes the identification and interpretation of CZ features inconclusive. To reliably resolve seismic velocity with depth, we implement a Transdimensional Hierarchical Bayesian (THB) framework with reversible‐jump Markov Chain Monte Carlo to generate samples from the posterior distribution of velocity structures. We also perform 2D synthetic tests to explore how well THB traveltime inversion can resolve different subsurface velocity structures. We find that THB recovers both sharp changes in velocity as well as gradual velocity increases with depth. Furthermore, we explore the velocity structure in a series of ridge‐valley systems in northern California. The posterior velocity model shows an increasing thickness of low velocity material from channels to ridgetops along a transect parallel to bedding strike, implying a deeper weathering zone below ridgetops and hillslopes than below channels. The THB method enhances the ability to reliably image CZ structure, and the model uncertainty estimates it yields provides an objective way to interpret deep CZ structure. The method can be applied across other near‐surface studies, especially in the presence of significant surface topography.
more »
« less
- Award ID(s):
- 2012616
- PAR ID:
- 10451091
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 22
- Issue:
- 5
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This dataset contains the codes and data used in the manuscript “Influence of Subsurface Critical Zone Structure on Hydrological Partitioning in Mountainous Headwater Catchments” submitted to Geophysical Research Letters. The software requirement are summarized in requirement.txt; hydrologic modeling input data are in the folder TLnewtest2sfb2; the observation data used in the simulation are indicated as comments in the python scripts. Note that the hydrologic modeling was run in HPC (Linux system) with parallel computing. Below are the abstract of the manuscript: “Headwater catchments play a vital role in regional water supply and ecohydrology, and a quantitative understanding of the hydrological partitioning in these catchments is critically needed, particularly under a changing climate. Recent studies have highlighted the importance of subsurface critical zone (CZ) structure in modulating the partitioning of precipitation in mountainous catchments; however, few existing studies have explicitly taken into account the 3D subsurface CZ structure. In this study, we designed realistic synthetic catchment models based on seismic velocity-estimated 3D subsurface CZ structures. Integrated hydrologic modeling is then used to study the effect of the shape of the weathered bedrock bottom on various hydrologic fluxes and storages in mountainous headwater catchments. Numerical results show that the shape of the weathered bedrock bottom not only affects the magnitude but also the peak time of both streamflow and subsurface dynamic storage.”more » « less
-
Abstract Understanding how soil thickness and bedrock weathering vary across ridge and valley topography is needed to constrain the flowpaths of water and sediment production within a landscape. Here, we investigate saprolite and weathered bedrock properties across a ridge‐valley system in the Northern California Coast Ranges, USA, where topography varies with slope aspect such that north‐facing slopes have thicker soils and are more densely vegetated than south‐facing slopes. We use active source seismic refraction surveys to extend observations made in boreholes to the hillslope scale. Seismic velocity models across several ridges capture a high velocity gradient zone (from 1,000 to 2,500 m/s) located ∼4–13 m below ridgetops that coincides with transitions in material strength and chemical depletion observed in boreholes. Comparing this transition depth across multiple north‐ and south‐facing slopes, we find that the thickness of saprolite does not vary with slope aspects. Additionally, seismic survey lines perpendicular and parallel to bedding planes reveal weathering profiles that thicken upslope and taper downslope to channels. Using a rock physics model incorporating seismic velocity, we estimate the total porosity of the saprolite and find that inherited fractures contribute a substantial amount of pore space in the upper 6 m, and the lateral porosity structure varies strongly with hillslope position. The aspect‐independent weathering structure suggests that the contemporary critical zone structure at Rancho Venada is a legacy of past climate and vegetation conditions.more » « less
-
Abstract Headwater catchments play a vital role in regional water supply and ecohydrology, and a quantitative understanding of the hydrological partitioning in these catchments is critically needed, particularly under a changing climate. Recent studies have highlighted the importance of subsurface critical zone (CZ) structure in modulating the partitioning of precipitation in mountainous catchments; however, few existing studies have explicitly taken into account the 3D subsurface CZ structure. In this study, we designed realistic synthetic catchment models based on seismic velocity‐estimated 3D subsurface CZ structures. Integrated hydrologic modeling is then used to study the effects of the shape of the weathered bedrock and the associated storage capacity on various hydrologic fluxes and storages in mountainous headwater catchments. Numerical results show that the weathered bedrock affects not only the magnitude but also the peak time of both streamflow and subsurface dynamic storage.more » « less
-
SUMMARY The near-surface seismic structure (to a depth of about 1000 m), particularly the shear wave velocity (VS), can strongly affect the propagation of seismic waves and, therefore, must be accurately calibrated for ground motion simulations and seismic hazard assessment. The VS of the top (<300 m) crust is often well characterized from borehole studies, geotechnical measurements, and water and oil wells, while the velocities of the material deeper than about 1000 m are typically determined by tomography studies. However, in depth ranges lacking information on shallow lithological stratification, typically rock sites outside the sedimentary basins, the material parameters between these two regions are typically poorly characterized due to resolution limits of seismic tomography. When the alluded geological constraints are not available, models, such as the Southern California Earthquake Center (SCEC) Community Velocity Models (CVMs), default to regional tomographic estimates that do not resolve the uppermost VS values, and therefore deliver unrealistically high shallow VS estimates. The SCEC Unified Community Velocity Model (UCVM) software includes a method to incorporate the near-surface earth structure by applying a generic overlay based on measurements of time-averaged VS in top 30 m (VS30) to taper the upper part of the model to merge with tomography at a depth of 350 m, which can be applied to any of the velocity models accessible through UCVM. However, our 3-D simulations of the 2014 Mw 5.1 La Habra earthquake in the Los Angeles area using the CVM-S4.26.M01 model significantly underpredict low-frequency (<1 Hz) ground motions at sites where the material properties in the top 350 m are significantly modified by the generic overlay (‘taper’). On the other hand, extending the VS30-based taper of the shallow velocities down to a depth of about 1000 m improves the fit between our synthetics and seismic data at those sites, without compromising the fit at well-constrained sites. We explore various tapering depths, demonstrating increasing amplification as the tapering depth increases, and the model with 1000 m tapering depth yields overall favourable results. Effects of varying anelastic attenuation are small compared to effects of velocity tapering and do not significantly bias the estimated tapering depth. Although a uniform tapering depth is adopted in the models, we observe some spatial variabilities that may further improve our method.more » « less
