2D Geotechnical site-response analysis including soil heterogeneity and wave scattering

This study describes an approach for modeling wave scattering and the spatial variability of ground motion in geotechnical site-response analysis by modeling soil heterogeneity through 2D correlated random fields. Importantly, the required site-specific inputs to apply the proposed approach in a practical setting are the same as those associated with conventional 1D site-response analysis. The results, which are affected by wave scattering attenuation, are compared to those from conventional laterally homogeneous 1D site-response analyses and 1D analyses with randomized velocity profiles extracted from heterogeneous 2D velocity model realizations. A sensitivity study, involving 5400 2D model realizations, investigates the influence of random field input parameters on wave scattering and site response. The computed ground surface acceleration waveforms and transfer functions show that this method is capable of scattering seismic waves. Multiple ground-motion intensity measures are analyzed to quantify this influence and distinguish between the effects of 1D vertical heterogeneities and averaging across many nodes and realizations, from the effects of wave scattering and 2D ground-motion phenomena. The redistribution of ground-motion energy across wider frequency bands and scattering attenuation of high-frequency waves in the 2D analyses resemble features observed in empirical transfer functions computed in other studies. While analyses with 1D randomized velocity profiles are able to replicate median results from 2D analyses for some low-frequency intensity measures (e.g. transfer functions at [Formula: see text] Hz, and spectral acceleration at the fundamental period), medians and standard deviations of high-frequency intensity measures (e.g. transfer function at [Formula: see text] Hz, [Formula: see text], and Arias intensity), which are influenced by wave scattering, are not appropriately captured. Given the equivalent input information requirements as conventional 1D analysis, and the availability of large computational resources, we advocate that the proposed 2D (and eventually 3D) approach is a fruitful path forward to improve the modeling of site-response physics and realize improved predictive capabilities.