A latent Gaussian process model for the spatial distribution of liquefaction manifestation

This paper presents a model for distributing zones of liquefaction and nonliquefaction for use in regional liquefaction risk analysis. There are two broad methodologies that have been used to evaluate liquefaction risk on the regional scale: (a) application of site-specific procedures using soil properties inferred from geology, or (b) application of geospatial proxies for liquefaction. The first approach will tend to predict similar liquefaction probabilities across broad areas with similar geology, water table depths, and shaking intensities. The second approach yields the probability of liquefaction, which can be interpreted as the portion of the area affected by liquefaction (%Aliq). Neither approach, however, gives an informed prediction of the spatial distribution of liquefaction and the resulting displacements, which are particularly important for assessments of seismic risk for spatially distributed infrastructure systems. We propose a methodology for incorporating spatial correlation into a geospatial proxy for liquefaction to create maps of liquefaction and nonliquefaction for a given earthquake scenario. First, we describe a latent Gaussian process that is assumed to govern the spatial distribution of liquefaction. Next, a database of empirical observations of liquefaction is used to obtain the coefficients that describe that latent Gaussian process. The proposed model yields random realizations of maps of liquefaction and nonliquefaction conditioned on a map of (%Aliq). Such maps can be used to constrain the area over which displacements are estimated using soil properties inferred from geology and are therefore a critical component in reducing bias in assessments of liquefaction risk at the regional scale.