Selecting the optimal intensity measure (IM) is essential for accurately assessing the seismic performance of the submarine shield tunnels in the layered liquefiable seabed. However, current research relies on simplistic ranking or filtering methods that neglect the different contributions of each evaluation criterion on IM's overall performance. To address this, this study begins by developing a numerical simulation method for nonlinear dynamic analysis, considering joint deformation, ocean environmental loads, and soil liquefaction, which is validated by experimental and theoretical methods. Subsequently, a fuzzy multiple criteria decision-making (FMCDM) method based on fuzzy probabilistic seismic demand models (FPSDM) is proposed, which integrates the fuzzy analytical hierarchical process (FAHP) for calculating weights and the fuzzy technique for order preference by similarity to ideal solution (FTOPSIS) for ranking IM alternatives. Finally, tunnel damage is classified into four states employing joint opening as the index for measuring damage, then the seismic fragility analysis is conducted. The results indicate that the optimal IM of a submarine shield tunnel situated in layered liquefiable seabed is sustained maximum velocity (SMV). Furthermore, the comparison between the fragility curves established using SMV and peak ground acceleration (PGA) reveals PGA, a frequently employed IM, notably undervaluing the seismic hazard.

This paper has attempted to determine the weighting levels of the soil and ground motion parameters (engineering bedrock depth (EBd), average shear wave velocity (Vs30), fundamental frequency (f0), peak ground acceleration (PGA), Joyner-Boore distance (Rjb), and epicenter distance (Repi)) in reflecting the actual damage status after the 2023 Kahramanmara & scedil; earthquakes, which have a wide impact area of 11 provinces. The analytical hierarchy method (AHP), a multi-criteria decision-making (MCDM) process, was used to analyze these parameter data sets obtained from 44 Disaster and Emergency Management Presidency of T & uuml;rkiye (AFAD) stations (Gaziantep, Hatay, Kahramanmara & scedil;, and Osmaniye). The priority order of the parameters before the analysis was systematically collected. These parameters were categorized into soil, ground motion and earthquake source-path properties. Considering the literature, these characteristics and their combined effects were systematically weighted with AHP under five groups. According to the weighted groups in the scope of the study, the actual damage data can be determined with a minimum accuracy rate of 70% (Group 1). In comparison, the best performance evaluation was 82% (Group 5). The parameter order and weights in the actual damage data evaluation are suggested as EBd-%28, PGA-%24, Vs30-%19, Rjb-%14, f0-%10, and Repi-%5 considering the very high accuracy rate of Group 5. This suggested weighting allows the rapid and effective estimation of the damage distribution after a possible earthquake only with soil, ground motion and earthquake source-path characteristics, even in cases where reliable structure data cannot be obtained.