Search

This study presents the development of two new sedimentary velocity models for the San Francisco Bay Area (SFBA) to improve the near-surface representation of shear-wave velocity ($V_S$) for large-scale, broadband numerical simulations, with the ultimate goal of enhancing the representation of the sedimentary layers in the Bay Area community velocity model. The first velocity model is stationary and is based solely on $V_{S30}$; the second velocity model is spatially varying and has location-specific adjustments. They were developed using a dataset of 200 measured $V_S$ profiles. Both models were formulated within a hierarchical Bayesian framework, using a parameterization that ensures robust scaling. The spatially varying model includes a slope adjustment term modeled as a Gaussian process to capture site-specific effects based on location. Residual analysis shows that both models are unbiased for $V_S$ values up to 1000 m/sec. Along-depth variability models were also developed using within-profile residuals. The proposed models show higher $V_S$ in the San Jose area and Livermore Valley compared to the USGS Bay Area community velocity model by a factor of two or more in some cases. Goodness-of-fit (GOF) comparisons using one-dimensional linear site-response analysis at selected sites demonstrate that the proposed models outperform the USGS model in capturing near-surface amplification across a broad frequency range. Incorporating along-depth variability further improves the GOF scores by reducing over-amplification at high frequencies. These results underscore the importance of integrating data-driven models of the shallow crust, like the ones presented here, in coarser regional community velocity models to enhance regional seismic hazard assessments., Comment: Manuscript: 18 pages, 12 figures. Supplement: 10 pages, 17 figures

This article extends a recently proposed heterogeneous data assimilation technique for site characterization to estimate compression and shear wave velocity (Vp and Vs, respectively) and damping at Treasure Island and Delaney Park downhole arrays. The adopted method is based on the joint inversion of earthquake acceleration time series and experimental surface wave dispersion data, and including physical constraints to improve the inverse problem's well-posedness. We first use synthetic data at these two sites to refine the proposed approach and then apply the refined algorithm to real data sets available at the Treasure Island and Delaney Park downhole arrays. The joint inversion results show that the estimated Vs and Vp profiles are in very good agreement with measured profiles at these two sites. Our synthetic and real data experiment results suggest that Vp estimation from inversion at downhole arrays can be improved by integrating the water table depth information or the higher modes of the Rayleigh wave dispersion data. Depending on the site complexity, water table information can also help reduce uncertainties associated with damping estimation. In the last part of this article, we compare the performance of the inverted profiles to other methods used to incorporate spatial variability and wave scattering effects in 1D ground response analysis (GRA). The comparisons show that the joint inversion-based Vs and Vp profiles and damping ratios estimated in this article can effectively integrate the effects of spatial variability and wave scattering into 1D GRAs, especially at the Delaney Park downhole array, which is classified as a poorly modeled site using traditional 1D GRA.

We present a sequential data assimilation algorithm based on the ensemble Kalman inversion to estimate the near-surface shear wave velocity profile and damping when heterogeneous data and a priori information that can be represented in forms of (physical) equality and inequality constraints in the inverse problem are available. Although non-invasive methods, such as surface wave testing, are efficient and cost effective methods for inferring Vs profile, one should acknowledge that site characterization using inverse analyses can yield erroneous results associated with the inverse problem non-uniqueness. One viable solution to alleviate the inverse problem ill-posedness is to enrich the prior knowledge and/or the data space with complementary observations. In the case of non-invasive methods, the pertinent data are the dispersion curve of surface waves, typically resolved by means of active source methods at high frequencies and passive methods at low frequencies. To improve the inverse problem well-posedness, horizontal to vertical spectral ratio (HVSR) data are commonly used jointly with the dispersion data in the inversion. In this paper, we show that the joint inversion of dispersion and strong motion downhole array data can also reduce the margins of uncertainty in the Vs profile estimation. This is because acceleration time-series recorded at downhole arrays include both body and surface waves and therefore can enrich the observational data space in the inverse problem setting. We also show how the proposed algorithm can be modified to systematically incorporate physical constraints that further enhance its well-posedness. We use both synthetic and real data to examine the performance of the proposed framework in estimation of Vs profile and damping at the Garner Valley downhole array, and compare them against the Vs estimations in previous studies.

Ensemble Kalman methods constitute an increasingly important tool in both state and parameter estimation problems. Their popularity stems from the derivative-free nature of the methodology which may be readily applied when computer code is available for the underlying state-space dynamics (for state estimation) or for the parameter-to-observable map (for parameter estimation). There are many applications in which it is desirable to enforce prior information in the form of equality or inequality constraints on the state or parameter. This paper establishes a general framework for doing so, describing a widely applicable methodology, a theory which justifies the methodology, and a set of numerical experiments exemplifying it.