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.