AI-Driven Prediction and Mapping of Soil Liquefaction Risks for Enhancing Earthquake Resilience in Smart Cities

1. Introduction

In an era of rapid urbanization and escalating environmental challenges, the concept of ‘smart cities’ as shown in Figure 1 has emerged as a pivotal framework for sustainable development [1,2,3,4]. These cities harness advanced technologies, such as IoT and big data analytics, to enhance quality of life, economic efficiency, and environmental sustainability [5]. A key element of this initiative is the integration of artificial intelligence (AI) into urban planning processes [6]. This integration not only revolutionizes resource management in urban areas but also enhances their resilience against natural hazards like earthquakes, torrential rain, and so on [7]. This study leads in applying smart technologies in geotechnical engineering, enhancing urban resilience through AI-driven risk prediction and zoning for soil liquefaction risk. By leveraging sophisticated AI algorithms to analyze extensive geotechnical data, this methodology facilitates the creation of detailed hazard maps. These maps are indispensable for effective urban planning and hazard mitigation, equipping urban planners and engineers with essential tools to make informed decisions that prioritize safety and sustainability.

The increasing frequency of urban hazards poses significant risks to densely populated areas, highlighting the critical importance of this study.

By leveraging predictive analytics, smart cities can proactively address these challenges, ensuring preparedness and adaptability within a dynamic urban environment. This study significantly enhances conventional methods by delivering a nuanced, precise, and dynamic evaluation of soil liquefaction risks, particularly in regions sensitive to geological and topographical variations. Moreover, the smart city design strategy outlined in this study encompasses comprehensive urban development, incorporating strategic infrastructure placement, optimized resource management, and improved citizen engagement through transparent data practices. By integrating AI-driven tools into the urban fabric, we ensure that cities are not only functional but also resilient and responsive to the needs of their inhabitants. The introduction of such advanced technological solutions tackles long-standing urban challenges and aligns with the broader objectives of smart cities. The predictive capabilities developed in this study represent a significant advancement in protecting and enhancing urban livability, marking a critical contribution to the conceptualization and design of smart cities worldwide.

This study focuses on the significant natural hazard of soil liquefaction [8,9,10,11,12], with examples illustrated in Figure 2, in Japan, where frequent seismic activity often causes severe ground instabilities. Soil liquefaction primarily occurs due to the shaking of loosely packed, water-saturated soils during earthquakes. This shaking increases pore water pressure and reduces the effective stress between soil particles, leading to fluid-like soil behavior, as depicted in Figure 2. Notable events, such as the 2011 Great East Japan Earthquake and the 2024 Noto Peninsula Earthquake, highlight the devastating impacts of these phenomena and underscore the critical need for accurate prediction and effective mitigation strategies [7]. This study utilizes AI algorithms to predict risks associated with soil liquefaction by analyzing extensive datasets from geological surveys and seismic activity. The ultimate goal is to develop an AI-driven predictive model that integrates various databases, including geotechnical and geographic information, using ensemble machine learning techniques. This model specifically aims to predict potential risks of soil liquefaction, with a particular focus on areas like Yokohama, Japan, which is known for its extensive reclaimed lands and high susceptibility to these hazards. By creating and validating hazard maps through case studies in such regions, the study aims to enhance urban planning and significantly contribute to the development of smart cities.

In the integration and visualization of the smart city concept illustrated in Figure 3, an AI-driven predictive model that synthesizes data from diverse databases, including geotechnical and geographic information, enhances urban resilience and promotes the development of a safer and more sustainable society. This approach contributes to the sustainable growth of smart cities and ensures the safety of their inhabitants. In smart cities, prioritizing the reduction of natural hazard risks is crucial. The AI-driven predictive model leverages geo-information databases to assess liquefaction risks, minimizing potential damage to buildings and infrastructure during earthquakes. The model’s high accuracy enables the identification of specific areas at risk of liquefaction, allowing for the proactive planning of effective countermeasures. Such capabilities expedite and streamline disaster responses, diminishing disruptions during emergencies. Additionally, this method provides a more detailed and accurate risk analysis than traditional approaches reliant on human judgment, thereby enhancing scientific decision-making in the planning and development within smart cities.

In this study, the authors develop a novel AI-driven predictive model that substantially advances the state of the art in assessing soil liquefaction risks, particularly in urban areas prone to seismic activity. Unlike previous studies [13,14,15,16], which primarily employed empirical or semi-empirical methods, this AI-driven model utilizes advanced ensemble machine learning techniques to integrate comprehensive geotechnical and geographic data. This integration yields a more precise, dynamic, and detailed assessment of soil liquefaction risks. The model is distinct in its ability to synthesize data from multiple sources, including real-time geological surveys and seismic records, thereby enhancing the accuracy of risk mapping and the efficacy of urban planning strategies. Moreover, the AI-driven predictive model addresses a critical gap in existing methodologies [13,14,15,16] by providing high-resolution hazard maps that are crucial for enhancing urban resilience in smart cities. By conducting a comparative analysis with models described in recent studies [13,14,15,16], this study not only fills a vital knowledge gap by improving the granularity and accuracy of predictions but also broadens the application of AI in geotechnical engineering, setting a new benchmark for future study and practical applications in earthquake-prone urban environments.

2. AI-Driven Predictive Models

In statistics and machine learning, ensemble machine learning trains multiple machine learning models (algorithms) and combines their predictions to obtain improved performance and better generalization ability than individual models. The fundamental idea behind ensemble machine learning is the recognition that each machine learning model has limitations and can make errors. Hence, ensemble machine learning aims to improve classification performance by harnessing the strengths of multiple models. Limitations of some machine learning models include high variance, high bias, and low accuracy [17,18]. However, several studies have shown that ensemble machine learning often achieves higher accuracy than a single machine learning model [19]. Ensemble machine learning can limit the variance and bias errors associated with a single model; for example, ensemble machine learning with a bagging algorithm reduces variance without increasing bias, while ensemble machine learning with a boosting algorithm reduces bias [20,21,22]. Overall, ensemble machine learning is more robust than an individual machine learning model [23]. Algorithms used in ensemble machine learning are mainly divided into three categories: bagging, boosting, and stacking. In this study, the boosting algorithm was selected for the ensemble machine learning.

2.1. Ensemble Machine Learning

Boosting is an algorithm that sequentially fits “weak” classifiers to different weightings of the observations in a dataset. Those observations, that the previous classifier poorly predicts, receive greater weight in the next iteration. The final classifier is a weighted average of all the weak classifiers. In a wide variety of problems, the weighting scheme and final classifier of boosting algorithm are merged, and the result has proven to be an effective method for reducing bias and variance.

The ensemble machine learning model is established using a boosting algorithm, and the procedure is illustrated in Figure 4, where sample data are randomly extracted from the training dataset. In this study, ensemble machine learning with the boosting algorithm was adopted, and an AI-driven predictive model was created by integrating a neural network with a gradient boosting decision trees model. By combining the strengths of different machine learning models, it is possible to achieve higher prediction accuracy than with a single model.

A neural network model shown in Figure 5a is a soft computing technique used today by engineers around the world for a wide variety of nonlinear problems that can occur in civil engineering [24]. And it is an interconnected collection of simple processing elements, units or nodes, whose functionality is loosely based on the animal neuron. The processing capability of the network is stored in the inter-unit connection strengths or weights, which are obtained through a process of adaptation to, or learning from, a set of training patterns [25,26]. Recently, neural networks have achieved outstanding results in learning image representations for various vision tasks, including image classification [27] and object detection [28]. These neural networks typically consist of several components, including convolutional, fully connected, and pooling layers. By stacking several of these layers, convolutional neural networks are able to learn complex features that are highly invariant and discriminative [29]. Previous studies have confirmed it as the most suitable and widely accepted statistical approach for soil liquefaction assessment [30,31].

A gradient boosting decision trees model shown in Figure 5b is a powerful machine learning tool that is widely used in many applications, including multi-class classification. This model uses decision trees as the base algorithm and sums the predictions of a series of trees. At each step, a new decision tree is trained to fit the residual errors between the ground truth and the current prediction. Gradient boosting decision trees are popular due to their accuracy, efficiency, and interpretability [32].

In ensemble machine learning model, the tuning of hyperparameters plays a critical role in optimizing model performance. Hyperparameters in ensemble algorithms, such as the number of base algorithms, the learning rate, and the depth of decision trees, require careful calibration to balance the tradeoff between bias and variance, thereby improving the generalizability of the model. The process of hyperparameter tuning involves a systematic search for the optimal settings that maximize the predictive accuracy of the ensemble. Techniques such as grid search, random search, and Bayesian optimization are commonly used for this purpose. These methods effectively explore the hyperparameter space and identify configurations that lead to the best performance on validation data sets. In this study, random search was employed to optimize model performance within the ensemble machine learning framework.

2.2. Reliability of Predictive Results by AI-Driven Predictive Model

In order to evaluate the performance of the AI-driven predictive model, “Accuracy” and “Coefficient of determination ($R^2$)” are selected as evaluation indicators.

“Accuracy” is the most basic index and is an evaluation index that indicates the percentage of correctly predicted predictions. This study calculated the “Accuracy” using Equation (1) [33,34] and the “Accuracy” was used as an evaluation index for predicting soil classification in this study.

$$Accuracy=\frac{TP+TN}{TP+FN+FP+TN}$$

(1)

where, $TP$ is true positives, $TN$ is true negatives, $FP$ is false positives, and $FN$ is false negatives.

“Coefficient of determination ($R^2$)” [35] is a key evaluation metric in predictive modeling. It quantifies the proportion of the variance in the dependent variable that is predictable from the independent variables. $R^2$ is bounded above by 1, which indicates a perfect fit. Although there is no theoretical lower bound for $R^2$, a value of 0 suggests that the model predicts as well as the mean of the dependent variable, represented by the horizontal line $y=K$, where $K$ is the mean of the target values of all training points. Values of $R^2$ less than 0 indicate a fit worse than this trivial model. Consequently, $R^2$ values are typically considered within the range of 0 to 1. The coefficient of determination remains unchanged under linear transformations of the scale of measurement for the independent variables. A value close to 1 indicates a model with high predictive accuracy. $R^2$ was calculated using Equation (2).

$$R^2=1-\frac{\sum _i{(y_i-x_i)}^2}{\sum _i({y_i-\overline{y})}^2}$$

(2)

where, $y_i$ is the measured value, $x_i$ is the predicted value, and $\overline{y}$ is the average of measured values.

3. Soil Liquefaction Risk Prediction with AI-Driven Predictive Model

3.1. Procedures on Soil Liquefaction Risk Prediction

In this study, data from the geological columnar section provided by the Ministry of Land, Infrastructure, Transport and Tourism, accessible through their online database, was utilized for the target location, Yokohama, Japan.

Data pertaining to “N-value,” “soil classification,” and “groundwater levels” were collected from the geological columnar section for a total of 214 locations, as depicted in Figure 6. Additionally, “elevation” data was obtained from the Geospatial Information Authority of Japan, using location information. An example of these collection records is shown in Figure 7. The criteria for determining soil liquefaction, according to the Highway Bridge Specifications and Commentary, are as follows:

(1)

The groundwater level must be within 10 m of the ground surface, and saturation must occur within 20 m of the ground surface.

(2)

The fine particle content ($FC$) must exceed 35%, with a plasticity index of 15 or less.

(3)

The average grain size (D₅₀) must be 10 mm or less, and the 10% grain size (D₁₀) must be 1 mm or less.

The AI-driven predictive model identifies layers with an N-value of 50 or more extending over at least 2 m as indicative of a bearing layer. Generally, soil or rock with an N-value of 20 or more is considered suitable as a load-bearing layer. If the N-value is between 30 and 50, the layer is deemed appropriate for the load-bearing layer of civil engineering and building structures. An N-value of 50 or more suggests a very solid layer, suitable as a load-bearing layer capable of supporting large structures such as high-rise condominiums. Therefore, in this study, continuous strata with an N-value exceeding 50 over a length of more than 3 m are designated as load-bearing layers [36].

Major soil classification systems in Japan have evolved in conjunction with public works projects [37,38]. In Japan, soils with similar properties are grouped based on their particle size distribution. This study categorizes “soil type classifications,” derived from geological columnar sections, into five categories: sand and gravel with fine content of 35% or less (prone to soil liquefaction), and silt, clay, and bedrock with fine content of 50% or more. Additionally, when multiple soil types, such as sandy silt or sandy clay mixed with gravel, are encountered, they are classified according to the component with the largest mass proportion, following the engineering classification method for soil materials.

Concerning “groundwater levels,” the geological columnar section contained sparse data, making it impractical to gather the necessary information for predictive analysis. Consequently, the borehole water level was assumed to represent the groundwater level, which was uniformly set at GL-1.0m.

Soil liquefaction risk is assessed for layers that satisfy the previously mentioned criteria. The AI-driven predictive model was trained using the data on “N-value” and “soil classification” to predict these parameters at unknown locations. The predicted data was then incorporated into the soil liquefaction potential index (LPI) to assess the risk of soil liquefaction.

3.2. Assumed Earthquake Motions for Soil Liquefaction Risk Prediction

Level 1 earthquake motion refers to lower-intensity seismic events that structures are expected to withstand without significant damage, typically maintaining operational functionality post-event. In contrast, Level 2 earthquake motion pertains to higher-intensity events. Buildings are designed to endure these, possibly sustaining some structural damage, to ensure occupant safety.

Regulations require consideration of two specific types of Level 2 earthquake ground motions: Type I for large plate-boundary earthquakes and Type II for inland earthquakes, such as the 1995 Hyogoken-Nanbu Earthquake. Type I earthquake motion is characterized by shorter duration and higher frequency vibrations, commonly experienced in regions close to tectonic plate boundaries. Conversely, Type II earthquake motion involves longer duration and lower frequency vibrations, typically occurring in areas more distant from plate boundaries. This classification into Types I and II aids in designing structures tailored to the specific seismic risks associated with the frequency and duration of ground shaking typical of a region.

3.3. Soil Liquefaction Potential Index

The soil liquefaction potential index ($LPI$) was utilized for analysis of soil liquefaction risk in this study. $LPI$, introduced in 1982, serves as a metric to assess the potential for soil liquefaction at a site, with higher values indicating greater risk. This index aggregates the factor of safety against liquefaction ($F_L$ value), which is weighted by depth, to evaluate the overall likelihood of soil liquefaction occurring at various depths within a site. A series of case studies have applied the LPI to assess liquefaction risks at targeted locations [38,39,40].

The $F_L$ value, essential for predicting soil liquefaction, quantifies the resistance to liquefaction at each depth. When the $F_L$ value is 1.0 or less, the occurrence of soil liquefaction is deemed likely. In this study, the $LPI$ was calculated using Equations (3)-(16) to comprehensively determine the risk of soil liquefaction at the study site.

$$LPI={\int }_0^{20}F\times (1-0.5x)d_x$$

(3)

$$F=\left\{\begin{array}{c}1-F_L(F_L<1.0)\\ 0 (F_L\ge 1.0)\end{array}\right.$$

(4)

$$F_L=\frac{R}{L}$$

(5)

where, $x$ represents the depth (m) from the ground surface; $F$ is the soil liquefaction resistance factor, which is governed by the factor of safety against soil liquefaction ($F_L$), $R$ is the dynamic shear strength ratio, and $L$ represents the earthquake shear stress ratio. $R$ is determined based on the standards for Level 1 and Level 2 earthquake motions.

$$R=C_W\times R_L$$

(6)

In the case of Level 1 and Type I of Level 2 earthquake motions:

$$C_W=1.0$$

(7)

In the case of Type II of Level 2 earthquake motion:

$$C_W=\left\{\begin{array}{c} 1.0 (R_L\le 0.1)\\ 3.3R_L+0.67 ({0.1<R}_L\le 0.4)\\ 2.0 \left({0.4<R}_L\right)\end{array}\right.$$

(8)

$$R_L=\left\{ \begin{array}{c} 0.0882\sqrt{(0.85N_a+2.1)/1.7} (N_a<14)\\ 0.0882\sqrt{N_a/1.7}+1.6\times {10}^{-6}\times {\left(N_a-14\right)}^{4.5 }(N_a\ge 14)\end{array}\right.$$

(9)

In the case of sand:

$$N_a=cFC\times \left(N_1+2.47\right)-2.47 (D_{50}<2mm)$$

(10)

$$N_1=170N/({{\sigma }_{vb}}^{\text{'}}+70)$$

(11)

$$C_{FC}=\left\{\begin{array}{c} 1 (0\%\le FC<10\%)\\ (FC+20)/30 (10\%\le FC\le 40\%)\\ (FC-16)/12 (40\%\le FC)\end{array}\right.$$

(12)

In the case of gravel:

$$N_a=\left\{1-\right.\left.0.36{\mathrm{l}\mathrm{o}\mathrm{g}}_{10}\times \left(\frac{D_{50}}{2}\right)\right\}{N}_1(D_{50}\ge 2mm)$$

(13)

$${L=r}_d{k}_{hgL}{{\sigma }_v{/\sigma }_v}^{\text{'}}$$

(14)

$$r_d=1.0-0.015x$$

(15)

$$k_{hgL}=C_Z{k}_{hgL0}$$

(16)

where, the explanation and definition of each coefficient and variable symbol used in Equations (3)–(16) are summarized in Table 1.

Table 2 presents the correction coefficients based on the characteristics of earthquake motion. Table 3 details the physical properties of different soil classifications related to existing groundwater.

When the $LPI$ equals 0, the risk or potential for soil liquefaction is very low. For 0 < $LPI$ ≤ 5, the risk or potential is low; for 5 < $LPI$ ≤ 15, it is high; and for $LPI$ > 15, it is very high. Table 4 summarizes the correlation between $LPI$ and liquefaction risk levels.

4. Results and Discussion

This study concentrates on geotechnical data derived from the geological columnar sections of Yokohama. Here, an AI-driven predictive model is employed to forecast the distribution of N-values and soil classifications. Case 1 involves predicting N-values, while Case 2 focuses on predicting soil classifications.

4.1. N-Values and Soil Classification Predictions

4.1.1. Case 1: Predictive Results for N-Value

Using geotechnical data collected from the geological columnar sections of Yokohama as training data, the AI-driven predictive model estimated the N-value at sites previously assumed to be unknown, though they were actually known.

To compare prediction accuracy based on data volume, the model’s performance was evaluated at 214 sites and at half that number. Model 1 used data from 117 sites, while Model 2 used data from all 214 sites.

Figure 8 shows the coefficient of determination ($R^2$) for predicting the N-value in Model 1 and Model 2, respectively. From the evaluation results for the N-value, it was observed that the $R^2$ tends to increase with increasing depth. This may be attributed to the increasing number of points (load-bearing layer) where the N-value exhibits a consistent range of fluctuation at greater depths, and the range of N-values at these points becomes narrower. Therefore, it can be inferred that predictions in deeper ground can be made with relatively high accuracy. Comparing prediction accuracy based on the number of data points reveals that the accuracy of the N-value predictions improves as the number of points increases. By gathering more data points to create training data, more reliable and precise predictions can be expected.

Therefore, Model 2 was selected for prediction. The N-value was predicted using data from 214 sites in Yokohama. Given the differing geologies of reclaimed land and plateaus, Yokohama was categorized into two types: ‘Landfill area (red)’ and ‘Plateau area (blue)’. Areas assumed to be unknown points were distinguished from known points in each category. Point 1 was selected from the landfill area and Point 2 from the plateau area, with their locations illustrated in Figure 6. Additionally, this study predicted the N-value, assumed to be an unknown point, using following three prediction procedures. The method that achieved the highest accuracy was selected:

(1)

A prediction procedure involving learning at every 1 m interval, without using the prediction results as data.

(2)

A prediction procedure that makes predictions at every 1 m interval from the ground surface to the subsurface, incorporating prediction results shallower than the predicted depth into the learning process.

(3)

A prediction procedure that makes predictions at every 1 m interval from 20 m below ground to the surface, incorporating prediction results deeper than the predicted depth into the learning process.

Figure 9 shows the results of predicting the N-value using each of the prediction procedures (1) to (3). Looking at the results for Point 1, while the results show that the N-value is low, prediction procedures (1) and (2) show that the N-value is 30 or more from points 10 m onwards. This is because the number of collected geotechnical data for a total of 214 points is more in the plateau area than in the landfill area, and there is much data that reaches the load-supporting layer from points around 10 m or higher and converges to an N-value of 30 or more. It is thought that this is having an influence. On the other hand, with prediction procedure (3), although the numerical value exceeds the survey results, it appears that similar behavior is predicted at points beyond 10 m while maintaining a relatively low N-value. As a result, at Point 1, prediction procedure (3) is more accurate than (1) and (2). Next, looking at the results for Point 2, all prediction procedures showed behavior similar to the survey results, and it can be said that predictions were made with relatively high accuracy. Among them, prediction procedure (3) can also predict the tendency of the survey N-value at a depth of around 5 m, and it is clear that the overall prediction accuracy is higher than (1) and (2). Based on the results of Points 1 and 2, it is thought that the prediction procedure (3) is better at making predictions at large depths because it learns and makes predictions from underground. Based on the above results, the following steps in this study. Therefore, the prediction procedure (3) was chosen for prediction.

4.1.2. Case 2: Predictive Results of Soil Classification

Similar to Case 1, geotechnical data collected from the geological columnar sections of Yokohama were utilized as training data to develop the AI-driven predictive model for soil classification. The results of the model’s prediction accuracy are displayed in Figure 10. Accuracy metrics were employed to assess the model’s performance. The results indicate high accuracy starting from a depth of 1.2 m, suggesting that predictions are generally reliable. Observations show that surface-near points are predominantly sandy, and as depth increases, the load-bearing layer is encountered more frequently, and the presence of bedrock points increases, which correlates with improved accuracy. Additionally, the prediction accuracy improves as the number of data points increases.

In addition, three prediction procedures of (1) to (3) used for N-value estimation were also applied to predict soil classification at unknown points. The results are detailed in Figure 11. According to these tables, the results at Point 1 show that the soil type, as identified from geotechnical data gathered from geological columnar sections, is characterized as sand up to about 3 m deep and silt below this depth. Furthermore, the geotechnical data reveal intermittent clay layers from 8 m to 15 m deep, which prediction procedures (1) and (2) failed to detect. In contrast, the prediction procedure (3), which analyzes data from deeper layers upwards, successfully predicted the presence of clay with relative accuracy.

For Point 2, mixed soil types are observed up to 10 m depth, beyond which the soil reaches the load-bearing layer and is identified as bedrock. The predictions accurately reflected the mixed soil types up to 10 m; however, they deviated from the geotechnical data at certain depths, resulting in less accurate classifications. At depths greater than 10 m, Prediction procedures (1) and (2) predicted only silt and failed to identify the bedrock. This discrepancy is likely due to the dataset predominantly classifying the ground that reaches the load-bearing layer as silt rather than bedrock. However, akin to its performance at Point 1, the prediction procedure (3) successfully predicted the presence of bedrock layers and mixed soil types below 10 m with greater accuracy. Consequently, the prediction procedure (3) was chosen for further soil classification predictions in this study.

4.2. Creation of Soil Liquefaction Risk Mapping

4.2.1. Comparison of Self-Created and Existing Soil Liquefaction Risk Maps

During the site selection and planning stages for engineering structures and human settlements, the soil liquefaction potential of an area is a critical factor considered in geotechnical engineering. Accurately evaluating the soil liquefaction potential is the first step toward mitigating soil liquefaction hazards. Soil liquefaction risk maps have been prepared by many autonomous bodies in Japan to assess the impact of soil liquefaction in each target area [41,42,43,44]. These maps are essential tools in the development of smart cities, where resilient urban planning and infrastructure are paramount. By integrating these hazard maps into urban planning processes, smart cities can optimize their development strategies to enhance safety and sustainability in earthquake-prone areas.

The creation of these maps can generally be divided into three steps: The first step involves calculating the spatial distributions of the instrumental seismic intensity for ten recent earthquakes that resulted in soil liquefaction. In the second step, the results of these seismic intensity distributions are compared with liquefied sites from past earthquakes to estimate the ratio of soil liquefaction occurrence. The third step involves merging geomorphological classes with similar soil liquefaction susceptibilities and applying a regression method to estimate the final probability of soil liquefaction occurrence based on instrumental seismic intensity.

Unlike previous approaches [13,14,15,16], this study first makes predictions for unknown locations and then uses these results to create hazard maps. Soil liquefaction risk maps were developed for three types of earthquake motions. Level 1 earthquake motion refers to seismic activity that is likely to occur once or twice during the lifespan of a structure; however, the intensity of such shaking is relatively low and is not expected to damage the structure. Figure 12a displays a self-created soil liquefaction risk map for Level 1 earthquake motion, while Figure 12b depicts the official soil liquefaction risk map for the same level published by Yokohama. These maps feature a 100 m x 100 m grid and are color-coded according to the soil liquefaction risk. Conversely, Level 2 earthquake motion represents a seismic event that is less likely to occur but could be extremely severe, potentially causing significant damage. This level is subdivided into Type I, which anticipates a major earthquake at the plate boundary, and Type II, which considers an inland earthquake. Figure 13a shows a self-created soil liquefaction risk map for Level 2 earthquake motion with Type I, while Figure 13b shows the official soil liquefaction risk map for the same level with type published by Yokohama. Figure 14a shows a self-created soil liquefaction risk map for Level 2 earthquake motion with Type II, while Figure 14b shows the official soil liquefaction risk map for the same level with type published by Yokohama.

The soil liquefaction risk maps created in this study reveal a relatively larger area at high risk of soil liquefaction compared to the maps published by Yokohama. Another major advantage of using AI models is the ability to quickly update maps as new geotechnical information and seismic data become available. This is a time-consuming task using standard methods, but can be done efficiently using AI. The maps created in this study were divided into detailed risk levels, making it easier to develop countermeasures for each level of risk. The highly accurate maps created by AI will provide important information to local residents and governments, helping them to formulate specific action plans to improve the safety of local communities. This is expected to minimize damage in the event of a liquefaction disaster. It can be confirmed that the liquefaction risk maps created in this study have many advantages over the existing maps provided by Yokohama. This will enable safer and more accurate disaster planning and further improve preparedness for future disasters.

4.2.2. Advantage of Soil Liquefaction Risk Maps Created by AI-Driven Predictive Model

The integration of AI into soil liquefaction risk prediction and mapping represents a transformative development in geotechnical engineering, particularly within the framework of smart cities. This paper explores the advantages of using AI-driven predictive models to produce detailed and dynamic soil liquefaction risk maps, enhancing urban resilience to earthquakes and other seismic activities. AI-driven predictive models represent a significant advancement over traditional liquefaction risk assessment methods, which typically rely on historical data and empirical formulas. While useful, traditional methods often fail to capture the complex, dynamic interactions within urban environments or adapt to rapid changes in urban landscapes. In contrast, AI models can continuously learn and adapt by integrating multiple data sources, including real-time geotechnical, geographic, and seismic data. This integration allows for a more nuanced and accurate assessment of liquefaction risk, providing urban planners and engineers with powerful tools for proactive urban development and disaster mitigation.

As shown in Table 5, one of the key advantages of AI-driven liquefaction risk maps is their ability to dynamically update and refine predictions as new data becomes available. This capability is critical in urban areas where conditions can change rapidly due to construction, land use changes, and environmental factors. AI models can adapt to these changes, maintaining the accuracy and relevance of risk maps. Additionally, AI’s ability to process large datasets enables the analysis of complex variables and interactions that would be too time-consuming or impossible for traditional methods, resulting in more detailed and comprehensive risk assessments. Another significant advantage of AI-driven hazard maps is their contribution to emergency preparedness and response strategies. Detailed hazard maps created from AI predictions can guide emergency planning and response efforts, helping to allocate resources more efficiently and improve response times during crises. For example, during an earthquake, AI-driven maps can quickly identify areas most at risk of soil liquefaction, allowing emergency services to prioritize these areas for immediate attention. In addition, these AI-driven maps facilitate better urban planning and infrastructure development. By predicting potential liquefaction zones, city planners can avoid building critical infrastructure in high-risk areas or implement special construction techniques to mitigate the risks. This proactive approach not only saves costs associated with repairing damage after a disaster but also protects lives by reducing the risk of catastrophic failure during seismic events.

The advantages of soil liquefaction risk maps created by AI-driven predictive models are numerous and significant. They provide a more accurate, detailed, and adaptive tool for managing geotechnical risk in urban environments. These maps enhance urban planning, improve emergency preparedness, and contribute to the overall resilience and sustainability of cities. As urban areas continue to grow and face complex challenges, the role of AI in geotechnical engineering will become increasingly important, marking a pivotal shift toward more intelligent, responsive, and resilient urban environments.

5. Conclusions

In this study, a self-developed AI-driven predictive model was trained to analyze geotechnical data collected from geological columnar sections at known points to predict the ‘N-value’ and ‘soil classification’ at unknown points. The predicted data were then used in the Liquefaction Potential Index ($LPI$) to assess soil liquefaction risk and to create maps for various levels of seismic motion. The findings include:

(1)

It was confirmed that the larger the training dataset used in the AI-driven predictive model, the higher the accuracy of the predictions.

(2)

The prediction procedure, which estimates the N-value and soil classification from 20 meters below ground to 1 meter above ground and incorporates learning from results deeper than the predicted depth, was found to be the most accurate.

(3)

The AI-driven predictive model provided more detailed soil liquefaction risk mapping by seismic motion level compared to existing mappings.

Future challenges include enhancing prediction accuracy. Despite utilizing the most accurate procedure available in this study, the predictions were not entirely consistent, with some demonstrating low accuracy due to biases in the training data. Further data collection, exploration of data processing methods, combinations of explanatory variables, and the development of new prediction techniques are necessary. In terms of assessing trigger risks, further improvements in prediction accuracy and the exploration of new evaluation methods may facilitate the creation of more precise hazard maps. Beyond assessing risks associated with soil liquefaction, this methodology could feasibly be applied to comprehensively evaluate other hazard risks and to generate maps predicting potential damage from various hazards.