An official website of the United States government
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
Using a Dense Seismic Array to Determine Structure and Site Effects of the Two Towers Earthflow in Northern California
Abstract We deployed a network of 68 three-component geophones on the slow-moving Two Towers earthflow in northern California. We compute horizontal-to-vertical spectral ratios (HVSRs) from the ambient seismic field. The HVSRs have two prominent peaks, one near 1.23 Hz and another between 4 and 8 Hz at most stations. The 1.23 Hz resonance is a property of the background noise field and may be due to a velocity contrast at a few hundred meters depth. We interpret the higher frequency peaks as being related to slide deposits and invert the spectral ratios for shallow velocity structure using in situ thickness measurements as a priori constraints on the inversion. The thickness of the shallowest, low-velocity layer is systematically larger than landslide thicknesses inferred from inclinometer data acquired since 2013. Given constraints from field observations and boreholes, the inversion may reflect the thickness of deposits of an older slide that is larger in spatial extent and depth than the currently active slide. Because the HVSR peaks measured at Two Towers are caused by shallow slide deposits and represent frequencies that will experience amplification during earthquakes, the depth of the actively sliding mass may be less relevant for assessing potential slide volume and associated hazard than the thicknesses determined by our inversions. More generally, our results underscore the utility of combining both geotechnical measurements and subsurface imaging for landslide characterization and hazard assessment.
more »
« less
- Award ID(s):
- 1848302
- PAR ID:
- 10183745
- Date Published:
- Journal Name:
- Seismological Research Letters
- Volume:
- 91
- Issue:
- 2A
- ISSN:
- 0895-0695
- Page Range / eLocation ID:
- 913 to 920
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY We have simulated 0–5 Hz deterministic wave propagation for a suite of 17 models of the 2014 Mw 5.1 La Habra, CA, earthquake with the Southern California Earthquake Center Community Velocity Model Version S4.26-M01 using a finite-fault source. Strong motion data at 259 sites within a 148 km × 140 km area are used to validate our simulations. Our simulations quantify the effects of statistical distributions of small-scale crustal heterogeneities (SSHs), frequency-dependent attenuation Q(f), surface topography and near-surface low-velocity material (via a 1-D approximation) on the resulting ground motion synthetics. The shear wave quality factor QS(f) is parametrized as QS, 0 and QS, 0fγ for frequencies less than and higher than 1 Hz, respectively. We find the most favourable fit to data for models using ratios of QS, 0 to shear wave velocity VS of 0.075–1.0 and γ values less than 0.6, with the best-fitting amplitude drop-off for the higher frequencies obtained for γ values of 0.2–0.4. Models including topography and a realistic near-surface weathering layer tend to increase peak velocities at mountain peaks and ridges, with a corresponding decrease behind the peaks and ridges in the direction of wave propagation. We find a clear negative correlation between the effects on peak ground velocity amplification and duration lengthening, suggesting that topography redistributes seismic energy from the large-amplitude first arrivals to the adjacent coda waves. A weathering layer with realistic near-surface low velocities is found to enhance the amplification at mountain peaks and ridges, and may partly explain the underprediction of the effects of topography on ground motions found in models. Our models including topography tend to improve the fit to data, as compared to models with a flat free surface, while our distributions of SSHs with constraints from borehole data fail to significantly improve the fit. Accuracy of the velocity model, particularly the near-surface low velocities, as well as the source description, controls the resolution with which the anelastic attenuation can be determined. Our results demonstrate that it is feasible to use fully deterministic physics-based simulations to estimate ground motions for seismic hazard analysis up to 5 Hz. Here, the effects of, and trade-offs with, near-surface low-velocity material, topography, SSHs and Q(f) become increasingly important as frequencies increase towards 5 Hz, and should be included in the calculations. Future improvement in community velocity models, wider access to computational resources, more efficient numerical codes and guidance from this study are bound to further constrain the ground motion models, leading to more accurate seismic hazard analysis.more » « less
-
Abstract Shallow bedrock strength controls both landslide hazard and the rate and form of erosion, yet regional patterns in near‐surface mechanical properties are rarely known quantitatively due to the challenge in collectingin situmeasurements. Here we present seismic and geomechanical characterizations of the shallow subsurface across the central Himalayan Range in Nepal. By pairing widely distributed 1D shear wave velocity surveys and engineering outcrop descriptions per the Geological Strength Index classification system, we evaluate landscape‐scale patterns in near‐surface mechanical characteristics and their relation to environmental factors known to affect rock strength. We find that shallow bedrock strength is more dependent on the degree of chemical and physical weathering, rather than the mineral and textural differences between the metamorphic lithologies found in the central Himalaya. Furthermore, weathering varies systematically with topography. Bedrock ridge top sites are highly weathered and have S‐wave seismic velocities and shear strength characteristics that are more typical of soils, whereas sites near valley bottoms tend to be less weathered and characterized by high S‐wave velocities and shear strength estimates typical of rock. Weathering on hillslopes is significantly more variable, resulting in S‐wave velocities that range between the ridge and channel endmembers. We speculate that variability in the hillslope environment may be partly explained by the episodic nature of mass wasting, which clears away weathered material where landslide scars are recent. These results underscore the mechanical heterogeneity in the shallow subsurface and highlight the need to account for variable bedrock weathering when estimating strength parameters for regional landslide hazard analysis.more » « less
-
Abstract Predicting rainfall‐induced landslide motion is challenging because shallow groundwater flow is extremely sensitive to the preexisting moisture content in the ground. Here, we use groundwater hydrology theory and numerical modeling combined with five years of field monitoring to illustrate how unsaturated groundwater flow processes modulate the seasonal pore water pressure rise and therefore the onset of motion for slow‐moving landslides. The onset of landslide motion at Oak Ridge earthflow in California’s Diablo Range occurs after an abrupt water table rise to near the landslide surface 52–129 days after seasonal rainfall commences. Model results and theory suggest that this abrupt rise occurs from the advection of a nearly saturated wetting front, which marks the leading edge of the integrated downward flux of seasonal rainfall, to the water table. Prior to this abrupt rise, we observe little measured pore water pressure response within the landslide due to rainfall. However, once the wetting front reaches the water table, we observe nearly instantaneous pore water pressure transmission within the landslide body that is accompanied by landslide acceleration. We cast the timescale to reach a critical pore water pressure threshold using a simple mass balance model that considers variable moisture storage with depth and explains the onset of seasonal landslide motion with a rainfall intensity‐duration threshold. Our model shows that the seasonal response time of slow‐moving landslides is controlled by the dry season vadose zone depth rather than the total landslide thickness.more » « less
-
We present a real-time spectral-scanning frequency-modulated continuous wave (FMCW) 3D imaging and velocimetry system that can produce 3D depth maps at 33 Hz, with 48° × 68° field of view (FOV) and 32.8-cm depth range. Each depth map consists of 507 × 500 pixels, with 0.095° × 0.14° angular resolution and 2.82-mm depth resolution. The system employs a grating for beam steering and a telescope for angular FOV magnification. Quantitative depth, reflectivity, and axial velocity measurements of a static 3D printed depth variation target and a moving robotic arm are demonstrated.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0220190206
