This study proposes a methodology for assessing seismic liquefaction hazard by implementing high-resolution three-dimensional (3D) ground models with high-density/high-precision site investigation data acquired in an area of interest, which would be linked to geotechnical numerical analysis tools. It is possible to estimate the vulnerability of earthquake-induced geotechnical phenomena (ground motion amplification, liquefaction, landslide, etc.) and their triggering complex disasters across an area for urban development with several stages of high-density datasets. In this study, the spatial-ground models for city development were built with a 3D high-precision grid of 5 m x 5 m x 1 m by applying geostatistic methods. Finally, after comparing each prediction error, the geotechnical model from the Gaussian sequential simulation is selected to assess earthquake-induced geotechnical hazards. In particular, with seven independent input earthquake motions, liquefaction analysis with finite element analyses and hazard mappings with LPI and LSN are performed reliably based on the spatial geotechnical models in the study area. Furthermore, various phenomena and parameters, including settlement in the city planning area, are assessed in terms of geotechnical vulnerability also based on the high-resolution spatial-ground modeling. This case study on the high-precision 3D ground model-based zonations in the area of interest verifies the usefulness in assessing spatially earthquake-induced hazards and geotechnical vulnerability and their decision-making support.

Earthquake event keeps increasing every year, and the recent cases of earthquake hazards invoke the necessity of seismic study in Korea, as geotechnical earthquake hazards, such as strong ground motion, liquefaction and landslides, are a significant threat to structures in industrial hub areas including coastal facilities. In this study, systemized framework of integrated assessment of earthquake-induced geotechnical hazard was established using advanced geospatial database. And a visible simulation of the framework was specifically conducted at two coastal facility areas in Incheon. First, the geospatial-grid information in the 3D domain were constructed with geostatistical interpolation method composed of multiple geospatial coverage mapping and 3D integration of geo-layer construction considering spatial outliers and geotechnical uncertainty. Second, the behavior of site-specific seismic responses were assessed by incorporating the depth to bedrock, mean shear wave velocity of the upper 30 m, and characteristic site period based on the geospatial-grid. Third, the normalized correlations between rock-outcrop accelerations and the maximum accelerations of each grid were determined considering the site-specific seismic response characteristics. Fourth, the potential damage due to liquefaction was estimated by combining the geospatial-grid and accelerations correlation grid based on the simplified liquefaction potential index evaluation method.