We have conducted three-dimensional (3D) 0–7.5 Hz physics-based wave propagation simulations to model the seismic response of the Long Valley Dam (LVD), which has formed Lake Crowley in Central California, to estimate peak ground motions and settlement of the dam expected during maximum credible earthquake (MCE) scenarios on the nearby Hilton Creek Fault (HCF). We calibrated the velocity structure, anelastic attenuation model, and the overall elastic properties of the dam via linear simulations of a Mw3.7 event as well as the Mw6.2 Chalfant Valley earthquake of 1986, constrained by observed ground motions on and nearby the LVD. The Statewide California Earthquake Center (SCEC) Community Velocity Model CVM-S4.26.M01 superimposed with a geotechnical layer using [Formula: see text] information tapered from the surface to a 700-m depth was used in the simulations. We found optimal fit of simulated and observed ground motions at the LVD using frequency-independent attenuation of [Formula: see text] ([Formula: see text] in m/s). Using the calibrated model, we simulated 3D nonlinear ground motions at the LVD for Mw6.6 rupture scenarios on the HCF using an Iwan-type, multi-yield-surface technique. We use a two-step method where the computationally expensive nonlinear calculations were carried out in a small domain with the plane wave excitation along the bottom boundary obtained from a full-domain 3D linear finite-fault simulation. Our nonlinear MCE simulation results show that peak ground velocities (PGVs) and peak ground accelerations (PGAs) as high as 72 cm/s and 0.55 g, respectively, can be expected at the crest of the LVD. Compared with linear ground motion simulation results, our results show that Iwan nonlinear damping reduces PGAs on the dam crest by up to a factor of 8 and increasingly depletes the high-frequency content of the waves toward the dam crest. We find horizontal relative displacements of the material inside the dam of up to [Formula: see text] and up to [Formula: see text] of vertical subsidence, equivalent to 1% of the dam height.
This content will become publicly available on July 1, 2024
- NSF-PAR ID:
- 10457303
- Date Published:
- Journal Name:
- GEOPHYSICS
- Volume:
- 88
- Issue:
- 4
- ISSN:
- 0016-8033
- Page Range / eLocation ID:
- KS101 to KS112
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
SUMMARY This paper deals with in situ characterization of the small-strain shear-wave velocity VS and damping ratio DS from an advanced interpretation of Multi-channel Analysis of Surface Waves (MASW) surveys. A new approach based on extracting Rayleigh wave data using the CFDBFa method has been discussed in the companion paper. This paper focuses on mapping the experimental Rayleigh wave phase velocity and attenuation into profiles of VS and DS versus depth, which is achieved through a joint inversion procedure. The joint inversion of phase velocity and attenuation data utilizes a newly developed Monte Carlo global search algorithm, which implements a smart sampling procedure. This scheme exploits the scaling properties of the solution of the Rayleigh eigenvalue problem to modify the trial earth models and improve the matching with the experimental data. Thus, a reliable result can be achieved with a limited number of trial ground models. The proposed algorithm is applied to the inversion of synthetic data and of experimental data collected at the Garner Valley Downhole Array site, as described in the companion paper. In general, inverted soil models exhibit well-defined VS profiles, whereas DS profiles are affected by larger uncertainties. Greater uncertainty in the inverted DS profiles is a direct result of higher variability in the experimental attenuation data, the limited wavelength range at which reliable values of attenuation parameters can be retrieved, and the sensitivity of attenuation data to both DS and VS. Nonetheless, the resulting inverted earth models agree well with alternative in situ estimates and geological data. The results stress the feasibility of retrieving both stiffness and attenuation parameters from active-source MASW testing and the effectiveness of extracting in situ damping ratio estimates from surface wave data.
-
more » « less
SUMMARY A robust, in situ estimate of shear-wave velocity VS and the small-strain damping ratio DS (or equivalently, the quality factor QS) is crucial for the design of buildings and geotechnical systems subjected to vibrations or earthquake ground shaking. A promising technique for simultaneously obtaining both VS and DS relies on the Multichannel Analysis of Surface Waves (MASW) method. MASW can be used to extract the Rayleigh wave phase velocity and phase attenuation data from active-source seismic traces recorded along linear arrays. Then, these data can be inverted to obtain VS and DS profiles. This paper introduces two novel methodologies for extracting the phase velocity and attenuation data. These new approaches are based on an extension of the beamforming technique which can be combined with a modal filter to isolate different Rayleigh propagation modes. Thus, the techniques return reliable phase velocity and attenuation estimates even in the presence of a multimode wavefield, which is typical of complex stratigraphic conditions. The reliability and effectiveness of the proposed approaches are assessed on a suite of synthetic wavefields and on experimental data collected at the Garner Valley Downhole Array and Mirandola sites. The results reveal that, under proper modelling of wavefield conditions, accurate estimates of Rayleigh wave phase velocity and attenuation can be extracted from active-source MASW wavefields over a broad frequency range. Eventually, the estimation of soil mechanical parameters also requires a robust inversion procedure to map the experimental Rayleigh wave parameters into soil models describing VS and DS with depth. The simultaneous inversion of phase velocity and attenuation data is discussed in detail in the companion paper.
-
more » « less
Seismic imaging and monitoring of the near-surface structure are crucial for the sustainable development of urban areas. However, standard seismic surveys based on cabled or autonomous geophone arrays are expensive and hard to adapt to noisy metropolitan environments. Distributed acoustic sensing (DAS) with pre-existing telecom fiber optic cables, together with seismic ambient noise interferometry, have the potential to fulfill this gap. However, a detailed noise wavefield characterization is needed before retrievingcoherent waves from chaotic noise sources. We analyze local seismic ambient noise by tracking five-month changes in signal-to-noise ratio (SNR) of Rayleigh surface wave estimated from traffic noise recorded by DAS along the straight university campus busy road. We apply the seismic interferometry method to the 800 m long part of the Penn State Fiber-Optic For Environment Sensing (FORESEE) array. We evaluate the 160 virtual shot gathers (VSGs) by determining the SNR using the slant-stack technique. We observe strong SNR variations in time and space. We notice higher SNR for virtual source points close to road obstacles. The spatial noise distribution confirms that noise energy focuses mainly on bumps and utility holes. We also see the destructive impact of precipitation, pedestrian traffic, and traffic along main intersections on VSGs. A similar processing workflow can be applied to various straight roadside fiber optic arrays in metropolitan areas.
-
more » « less
There is a growing need to characterize the engineering material properties of the shallow subsurface in three dimensions for advanced engineering analyses. However, imaging the near-surface in three dimensions at spatial resolutions required for such purposes remains in its infancy and requires further study before it can be adopted into practice. To enable and accelerate research in this area, we present a large subsurface imaging data set acquired using a dense network of three-component (3C) nodal stations acquired in 2019 at the Garner Valley Downhole Array (GVDA) site. Acquisition of this data set involved the deployment of 196 stations positioned on a 14 × 14 grid with a 5 m spacing. The array was used to acquire active-source data generated by a vibroseis truck and an instrumented sledgehammer, and passive-wavefield data containing ambient noise. The active-source acquisition included 66 vibroseis and 209 instrumented sledgehammer source locations. Multiple source impacts were recorded at each source location to enable stacking of the recorded signals. The active-source recordings are provided in terms of both raw, uncorrected units of counts and corrected engineering units of meters per second. For each source impact, the force output from the vibroseis or instrumented sledgehammer was recorded and is provided in both raw counts and engineering units of kilonewtons. The passive-wavefield data include 28 h of ambient noise recorded over two nighttime deployments. The data set is shown to be useful for active-source and passive-wavefield three-dimensional imaging and other subsurface characterization techniques, which include horizontal-to-vertical spectral ratios (HVSRs), multichannel analysis of surface waves (MASW), and microtremor array measurements (MAM).
