A Methodology for Susceptibility Assessment of Wave-Induced Seabed Liquefaction in Silt-Dominated Nearshore Environments

: Wave-induced seabed liquefaction significantly jeopardizes the stability of marine structures and the safety of human life. Susceptibility assessment is key to enabling spatial predictions and establishing a solid foundation for effective risk analysis and management. However, the current research encounters various challenges, involving an incomplete evaluation system, poor applicability of methods, and insufficient databases. These issues collectively hinder the accuracy of susceptibility assessments, undermining their utility in engineering projects. To address these challenges, a susceptibility assessment method with the safety factor was developed as the key assessment parameter, allowing for a comprehensive susceptibility assessment across the silt-dominated nearshore environment using Empirical Bayesian Kriging (EBK). The safety factor is determined by combining the cyclic stress ratio (CSR) and the cyclic resistance ratio (CRR), which characterize wave loadings and sediment properties in the study area, respectively. This method was applied in the Chengdao region of the Yellow River Estuary, China, a typical silt-dominated nearshore environment where wave-induced liquefaction events have been reported as being responsible for multiple oil platform and pipeline accidents. By collecting the regional wave and seabed sediment data from cores spanning from 1998 to 2017, the safety factors were calculated, and a zonal map depicting the susceptibility assessment of wave-induced seabed liquefaction was created. This study can serve as a valuable reference for the construction and maintenance of marine engineering in liquefaction-prone areas.


Introduction
Seabed liquefaction is a phenomenon associated with sediment instability.Under the influence of external forces (e.g., wave action, earthquakes, and human activities), the internal stress state of the sediment changes.Specifically, when the excess pore pressure is equivalent to the effective stress of the overlying sediment, the state of the sediment changes from solid to fluid, and the sediment becomes liquefied [1-3].Wave-induced seabed liquefaction is caused by an external load cyclically applied to the seabed, generating differential loading on the seafloor through pressure waves that induce a series of cyclic shear stresses in the underlying sediment.If the induced shear stress exceeds the strength, it may lead to substantial deformation or liquefaction failure [1].This sequence of this event can lead to various engineering safety problems and compromise the stability of Figure 1.Effects of seabed liquefaction on coastal environments.Because of seabed liquefaction, differential settlement of coastal areas causes cracks in nearshore houses and walls (a); sediment strength is reduced, causing destabilization of offshore platforms (b); and the seabed soil around the trench turned into a dense liquid that flooded the trench and caused the pipeline to fail (c) (images from Di Fiore et al. [6], Chávez et al. [7], and Damgaard et al. [8]).
To predict and evaluate the occurrence probability and severity of seabed liquefaction effectively, numerous researchers have started applying risk assessment and management techniques to this field [9][10][11][12].These studies contribute to the quantitative analysis of the spatial extent of disasters, geological disaster mechanisms, and triggering factors [13,14].The insights derived from these analyses provide recommendations for predisaster planning, disaster response, and post-disaster recovery efforts, ultimately mitigating losses.This field has garnered significant attention from governments, industries, and the scientific research community [15].Nevertheless, the complex nature of waveinduced seabed liquefaction mechanisms makes it challenging to identify assessment factors.This difficulty hinders the advancement of risk assessment and management in this specific context.Among these challenges, the susceptibility assessment of seabed liquefaction emerges as a pivotal aspect in the risk assessment of seabed liquefaction.It directly impacts the accuracy of risk assessment and requires the advancement of reliable theories and practical methodologies.Susceptibility assessment of wave-induced seabed liquefaction primarily involves the creation of a disaster inventory, selection of conditioning factors, and identification of assessment methods [16].Because of the complexity of waveinduced seabed liquefaction mechanisms, the choice of different monitoring tools for liquefaction under different sediment conditions directly affects the establishment of the disaster inventory and the selection of condition factors, and ultimately determines the research methods [17].Therefore, priority should be given to addressing the accuracy of susceptibility assessment methods.The existing research methods lack the capability to adequately evaluate susceptibility across various sediment types, highlighting the necessity for the development of new methods to enhance evaluation accuracy [17,18].
Currently, susceptibility assessment methods for seabed liquefaction predominantly revolve around statistical approaches and deterministic analysis.The statistical approach enables both qualitative and quantitative evaluations of geological disasters by employing statistical principles to describe relationships among conditioning factors [19].Consequently, this method, including bivariate statistical analysis, multivariate statistical analysis, and machine learning (such as neural network models and maximum entropy models), finds extensive application in the susceptibility analysis of marine geological disasters [16].However, the statistical approach primarily focuses on the probability of disasters and allocates less consideration to the mechanism of disasters [19].In contrast, the deterministic approach focuses on the mechanics of disaster occurrence.It has gained Figure 1.Effects of seabed liquefaction on coastal environments.Because of seabed liquefaction, differential settlement of coastal areas causes cracks in nearshore houses and walls (a); sediment strength is reduced, causing destabilization of offshore platforms (b); and the seabed soil around the trench turned into a dense liquid that flooded the trench and caused the pipeline to fail (c) (images from Di Fiore et al. [6], Chávez et al. [7], and Damgaard et al. [8]).
To predict and evaluate the occurrence probability and severity of seabed liquefaction effectively, numerous researchers have started applying risk assessment and management techniques to this field [9][10][11][12].These studies contribute to the quantitative analysis of the spatial extent of disasters, geological disaster mechanisms, and triggering factors [13,14].The insights derived from these analyses provide recommendations for pre-disaster planning, disaster response, and post-disaster recovery efforts, ultimately mitigating losses.This field has garnered significant attention from governments, industries, and the scientific research community [15].Nevertheless, the complex nature of wave-induced seabed liquefaction mechanisms makes it challenging to identify assessment factors.This difficulty hinders the advancement of risk assessment and management in this specific context.Among these challenges, the susceptibility assessment of seabed liquefaction emerges as a pivotal aspect in the risk assessment of seabed liquefaction.It directly impacts the accuracy of risk assessment and requires the advancement of reliable theories and practical methodologies.Susceptibility assessment of wave-induced seabed liquefaction primarily involves the creation of a disaster inventory, selection of conditioning factors, and identification of assessment methods [16].Because of the complexity of wave-induced seabed liquefaction mechanisms, the choice of different monitoring tools for liquefaction under different sediment conditions directly affects the establishment of the disaster inventory and the selection of condition factors, and ultimately determines the research methods [17].Therefore, priority should be given to addressing the accuracy of susceptibility assessment methods.The existing research methods lack the capability to adequately evaluate susceptibility across various sediment types, highlighting the necessity for the development of new methods to enhance evaluation accuracy [17,18].
Currently, susceptibility assessment methods for seabed liquefaction predominantly revolve around statistical approaches and deterministic analysis.The statistical approach enables both qualitative and quantitative evaluations of geological disasters by employing statistical principles to describe relationships among conditioning factors [19].Consequently, this method, including bivariate statistical analysis, multivariate statistical analysis, and machine learning (such as neural network models and maximum entropy models), finds extensive application in the susceptibility analysis of marine geological disasters [16].However, the statistical approach primarily focuses on the probability of disasters and allocates less consideration to the mechanism of disasters [19].In contrast, the deterministic approach focuses on the mechanics of disaster occurrence.It has gained significant attention from researchers for its ability to quantitatively assess susceptibility to wave-induced seabed liquefaction [16].A 3-D assessment methodology for seabed liquefaction has been proposed, based on cyclic triaxial tests, liquefaction potential evaluation criteria, and the nearshore spectral windwave (NSW) model.Nonetheless, this method did not take into account the effects of pore water pressure dissipation that could take place during a storm period, rendering it unsuitable for completing the susceptibility assessment of shallow sediment liquefaction [20].Di Fiore et al. explored wave-induced seabed liquefaction through cyclic triaxial tests and evaluated the liquefaction potential of the Magoodhoo coralline island and coastal regions [6].However, their calculation method for CSR (the cyclic stress ratio) was only applicable to wave heights ranging from 0.2 to 3.0 m.Furthermore, the evaluation factor was solely wave height, resulting in an oversimplified assessment model that compromises result accuracy [6].In summary, the current evaluation method lacks a comprehensive consideration of evaluation factors for generating zoning maps based on evaluation results, and it falls short in assessing different sediment types.Therefore, it is crucial to introduce a more precise methodology to overcome these challenges.
In this study, a method for the susceptibility assessment of wave-induced seabed liquefaction is introduced, with a particular focus on safety factors, which have been widely used in susceptibility assessment for various natural processes.This method characterizes wave loading through CSR and employs CRR (the cyclic resistance ratio) to represent sediment properties to calculate safety factors.The calculation of CSR follows the method proposed by Ishihara et al. [21], while CRR analysis integrates considerations of sediment characteristics.Utilizing this approach in the Chengdao region involved calculating CSR values for diverse wave recurrence periods (5-year, 10-year, and 25-year).Concurrently, sediment properties determined CRR at varying depths, assessed through cyclic triaxial torsion shear tests at sampling points.Establishing the smallest safety factor among these points facilitated the identification of the wave-induced seabed liquefaction susceptibility assessment.The susceptibility zoning map was generated using Empirical Bayesian Kriging (EBK) in the Arc Geographic Information System 10.7 (ArcGIS 10.7).This method is applicable to seabed areas with silt-dominated environments and can maximize the accuracy of small-scale susceptibility assessments of wave-induced seabed liquefaction, even when data are limited.

Susceptibility Assessment Model of Wave-Induced Seabed Liquefaction
Currently, various methods exist for the susceptibility assessment of wave-induced seabed liquefaction, commonly categorized as a statistical approach or deterministic analysis [16].Deterministic analysis, as discussed earlier, allows for a more comprehensive exploration of the mechanism behind wave-induced seabed liquefaction, aligning closely with its characteristics and thus finding extensive application [16].The safety factor employed in this study falls within a deterministic analysis, offering quantitative insights into liquefaction susceptibility.This method was initially proposed by Seed et al. in 1971 [22].Its versatility allows for application across study areas with diverse sediment types and various scenarios, promising a more precise susceptibility of wave-induced seabed liquefaction [23].This method employs CSR to signify the extent of external load impact on sediment liquefaction and uses CRR to indicate sediment resistance properties [22,24].The safety factor can be determined by comparing CSR and CRR at varying depths according to Equation (1): where FS is the safety factor; CSR is the cyclic stress ratio; and CRR is the cyclic resistance ratio.The seabed liquefaction susceptibility of each area is further classified according to the characteristics of the study area and the probability of liquefaction, combined with the safety factor [25].
It is usually considered that the liquefaction will not occur in the region when FS > 1 and that the liquefaction will occur in the region when FS < 1.The interrelation between the safety factor and reliability had been explored, proposing that the safety factor can be specified based on reliability levels, and the allocation of the safety factor should be based on having the same reliability associated with each failure mode [26].The influence of geotechnical parameter variability on the liquefaction potential of tailing dams has been explored [27].The safety factor for different liquefaction probabilities has been compared, highlighting the use of probabilistic methods for failure probability and the safety factor against liquefaction analysis [27].Therefore, the division intervals of safety factors need to be determined based on the actual conditions of the research area, considering the characteristics of disaster events, and based on different liquefaction probabilities.

Cyclic Stress Ratio
CSR represents the ratio of the amplitude of shear stress to the effective confining stress, which is used to quantify the deformability or strength of soils under cyclic loading conditions [22].According to linear wave theory, a wave event can be described by fundamental wave characteristics [28].Thus, various authors have suggested simplifying the CSR calculation based on this theory [21].When analyzing wave conditions, the following criteria must be considered: (1) to reduce computational complexity, the irregular loads, induced by actual waves on the seabed surface, are simplified to uniform loads and (2) because this method is designed for shallow sea areas, it assumes that all waves are linear [23,29,30].
Regarding the characteristics of waves, they can be considered to consist of an infinite number of wave trains with the same amplitude and wavelength.The passage of such an array of waves over the ocean creates harmonic pressure waves on the seafloor, increasing the pressure under the crest and reducing it under the trough, as shown in Figure 2 [21].Considering the rectangular coordinate system, the sea surface elevation of a wave field η may be represented by Ishihara and Yamazaki's Equation (2) [21]: where H is the wave height, k = 2π L is the wave number, L is the wave length, ω = 2π T is the wave frequence, T is the wave period, x is the spatial coordinate in the x direction, and t is the time coordinate.For linear waves, k and ω are related by the dispersion relationship according to Equation (3) [1]: where g is the gravitational acceleration.Equation (3) also can be rewritten as Equation ( 4): where L is the wave length, T is the wave period, and d is the water depth.Regarding the characteristics of wave-induced cyclic stress, the seabed deposit i sumed to consist of a homogeneous elastic material extending to an infinite depth.B sinesq's classical solution for the two-dimensional plane strain problem can be emplo to calculate the stresses acting on it [31].The stress induced in the seabed is, theref analyzed by applying a sinusoidally changing load on the horizontal surface from m to plus infinity according to Equation (5): Regarding the characteristics of wave-induced cyclic stress, the seabed deposit is assumed to consist of a homogeneous elastic material extending to an infinite depth.Boussinesq's classical solution for the two-dimensional plane strain problem can be employed to calculate the stresses acting on it [31].The stress induced in the seabed is, therefore, analyzed by applying a sinusoidally changing load on the horizontal surface from minus to plus infinity according to Equation (5): where L is the wave length, x is the spatial coordinate in the x direction, T is the wave period, and t is the time coordinate.According to the small amplitude wave, the amplitude of the pressure fluctuation exerted on the sea bottom, p 0 , by the travelling wave is given by Horikawa's Equation ( 6) [32]: where ρ w is the density of seawater, g is the gravitational acceleration, H is the wave height, d is the water depth, and L is the wave length.Considering the uniform distribution and infinite depth of the sediments, the amplitude of maximum shear stress can be determined by analyzing the relationship between the horizontal stress, vertical stress, and shear stress of the harmonic pressure wave.The analysis utilizes the classical Boussinesq solution for two-dimensional plane strain problems according to Equation (7) [21]: where τ vh is the shear stress, p 0 is the amplitude of the pressure fluctuation exerted on the sea bottom, L is the wave length, and z is the sampling depth.Since the magnitude of σ ′ v is expressed by σ ′ v = ρ ′ gz, CSR is given by Equation (8) [21]: where τ vh is the shear stress, σ ′ v is the vertical effective stress, p 0 is the amplitude of the pressure fluctuation exerted on the sea bottom, ρ ′ is the submerged unit mass of soils in the seabed, g is the gravitational acceleration, L is the wave length, and z is the sampling depth.According to the formula, comprehensive data on fundamental wave characteristics like wave height, wavelength, wave period, and water depth is essential.This information not only allows for the calculation of CSR but also facilitates the simulation of CSR under extreme wave conditions.

Cyclic Resistance Ratio
CRR is a parameter that measures the resistance of sediments to liquefaction.It represents the ratio between the average cyclic shear stress and the vertical effective stress of sediments [22].Different sediment types require specific calculation methods, and the operational ease varies among various test methods.The field test methods of CRR under seismic liquefaction were summarized by referring to the research results of Youd et al., as shown in Table 1 [24].The primary field test includes SPT (Standard Penetration Test), CPT (Cone Penetration Test), measurement of V S (Shear Wave Velocity), and BPT.SPT is a test conducted during a test boring in the field to measure the approximate soil resistance to penetration of a split-spoon sampler at various depths below the ground surface [33].This test was intended to measure the number of blows required to drive a standard sampling tube into the ground to a certain depth [34].As for CPT, an elongated metal probe (typically conical in shape) is pushed into the ground, and the resistance and soil displacement during the penetration process are measured.It is performed by pushing an instrumented probe with a specific diameter into the earth at a constant speed while simultaneously measuring the cone resistance and sleeve friction resistance [35].The V S test is a geoengineering test used to determine the shear wave velocity in soil or rock.A shear wave is a wave that propagates through shear deformation in a material, and its speed is related to the shear modulus of the material.Shear wave velocity tests are often used to study the elastic properties and mechanical behavior of soil or rock [36].When conducting liquefaction assessments, SASWs are useful for V S profile surveys of liquefaction sites [37].Becker penetration resistance is defined as the number of blows required to drive the casing through an increment of 300 mm [24].The formula for calculating CRR corresponding to the above method is mentioned in Table 2.It is essential to choose the appropriate method based on the actual conditions of CRR calculation, significantly enhancing the accuracy of CRR and FS and consequently improving the precision of susceptibility assessment results.

SASW
Note: N 1,60 is the standard penetration test blowcount value corrected for overburden, energy, equipment, and procedural factors; θ i is the set of unknown model coefficients; FC is the fines content; M w is the earthquake moment magnitude; σ ′ v is the vertical effective stress; P a is the atmospheric pressure (1 atm); Φ is the standard cumulative normal distribution; P L is the probability of triggering of liquefaction; q c,1 is the normalized cone tip resistance; R f is the friction ratio; V s1 is the effective stress normalized shear-wave velocity.
The above method is suitable for the calculation of CRR in the terrestrial environment.However, the monitoring difficulty in the marine environment is greater than that in the terrestrial environment, and the in situ test is more difficult and costly [40,41].Therefore, the laboratory test has attracted more attention in the field of submarine geological disaster susceptibility assessment.The cyclic triaxial torsion shear test is a widely used sediment testing method suitable for saturated sand and saturated silt in nearshore areas [31].It offers a streamlined computation process with fewer evaluation factors, facilitating a rapid assessment of sediment resistance to liquefaction.This test reflects the cyclic loading characteristics of the sediment, a characteristic that is one of the determining factors in the susceptibility assessment of wave-induced seabed liquefaction [31].Introduced by Ishihara et al., this method selects the dynamic stress ratio when the double amplitude deviator strain reaches 5% as CRR for the sediments at the given point [21].Based on this standard, they conducted a large number of cyclic triaxial tests, summarized the relevant rules, and proposed the calculation according to the relative density of sediments and the coefficient of earth pressure at rest, as shown in Equation (9) [21]: where D r is the relative density and K 0 is the coefficient of earth pressure at rest.However, this method is applicable to sandy sediments, but not to silt, so some scholars have found another universal rule, which indicates that the cyclic stress ratio at 100 load cycles is determined to be closest to this criterion, as illustrated in Equation (10) [21,23].Since there exists an exponential relationship between the number of cyclic loadings and the cyclic dynamic stress in the cyclic triaxial test, as depicted in Equation ( 11), the CRR can be calculated by fitting this exponential relationship with convincing data obtained through the tests [23].
where τ L is the horizontal shear stress, σ ′ v is the vertical effective stress, σ d is the cyclic dynamic stress, σ C is the average consolidation pressures, σ 1 is the axial consolidation pressures, and σ 3 is the lateral consolidation pressures.
where τ L is the horizontal shear stress, σ ′ v is the vertical effective stress, N is the number of cyclic loadings, and a and b are the fitting coefficients.

Geospatial Interpolation
The previous process allows us to calculate the safety factor for a specific location.However, to complete a susceptibility assessment of marine geological disasters for the entire study area, spatial interpolation tools within the Arc Geographic Information System 10.7 (ArcGIS 10.7) need to be employed.These tools maximize the utilization of existing data points to comprehensively delineate the susceptibility assessment of liquefaction across the entire area [9].To reconstruct a continuous attribute distribution over the entire study area, the spatial interpolation method utilizes a finite set of sampling points S = {(x i , f i ), i = 1, 2, 3, . . ., n} to establish an interpolation function ( f : x → x) that estimates attribute values at any given data point within the study area [10].Spatial interpolation can be categorized into geostatistical interpolation and deterministic interpolation.A widely used geostatistical interpolation method is EBK, which is valued for its simplicity and accurate results [42,43].This method accounts for introduced error by estimating the underlying semivariance function, whereas other kriging methods calculate the semivariance function from a known data location and apply this single semivariance function to make predictions at an unknown location [9,44].By not considering the uncertainty in the estimation of the semivariance function, other kriging methods underestimate the standard error of the prediction, while the EBK is more accurate [44].Additionally, EBK requires minimal interactive modeling, excels at predicting unstable data, and is particularly well-suited for smaller datasets [9].

Susceptibility Assessment System of Wave-Induced Seabed Liquefaction
Based on the above discussion of the susceptibility assessment and spatial interpolation methods, the specific steps for partitioning the susceptibility assessment of wave-induced seabed liquefaction are shown in Figure 3.

Susceptibility Assessment System of Wave-Induced Seabed Liquefaction
Based on the above discussion of the susceptibility assessment and spatial inte tion methods, the specific steps for partitioning the susceptibility assessment of wa duced seabed liquefaction are shown in Figure 3.  [14,22,43,44]).
This study proceeds through the following steps: (1) Identify the study area through a comprehensive review of historical data.Su rize the fundamental wave characteristics within this region and establish a co hensive disaster inventory.(2) Calculate the CSR of all measurement points using Equation ( 8), considerin measured water depth and wave parameters of different wave recurrence pe Utilize the cyclic triaxial torsion shear test data from sampling points to cal Figure 3. Process for the susceptibility assessment of wave-induced seabed liquefaction (modified from [14,22,43,44]).
This study proceeds through the following steps: (1) Identify the study area through a comprehensive review of historical data.Summarize the fundamental wave characteristics within this region and establish a comprehensive disaster inventory.(2) Calculate the CSR of all measurement points using Equation ( 8), considering the measured water depth and wave parameters of different wave recurrence periods.Utilize the cyclic triaxial torsion shear test data from sampling points to calculate CRR at the known depth of each sampling point.Select the minimum CRR value for each measurement point as the CRR value for susceptibility assessment.Calculate the safety factor of wave-induced seabed liquefaction at the known depth of each sampling point under the different wave conditions by applying Equation (1).Select the minimum safety factor value for each point under varying wave conditions.( 3) Analyze and process the data through ArcGIS 10.7, incorporating interpolation methods to generate a liquefaction susceptibility zoning map for the study area.Subsequently, conduct an in-depth analysis of the zoning results.

Study Area
Numerous studies indicate that nearshore silty seabed are typical areas prone to liquefaction, posing a significant threat to the stability of nearshore marine engineering [7,[45][46][47].In this context, the Chengdao Island area, where silty sediments are widely distributed, has been selected as the study area.The Chengdao region is situated to the north of the Yellow River Estuary in China, exhibiting a southwest-to-northeast topographical gradient, as depicted in Figure 4.This area primarily consists of two subdeltas [47].The first subdelta was formed when the Yellow River was redirected to the Shenxiangou flow path from 1953 to 1964, while the second subdelta was established when the river was diverted to the Diaokou flow path from 1964 to 1976 [48].The underwater Yellow River Estuary spans an area of approximately 3000 km 2 .A distinctive feature of this region is its remarkably gentle gradient, typically less than 0.4 • , and an exceptionally low slope drop ranging from 1% to 5% [48].This makes the area highly susceptible to storm surge hazards [1].Between 1949 and 2005, more than 150 storm surge events occurred in the region, with an average of two to three events per year, of which 36 storm surges caused geological disasters, most of which occurred in the summer and fall [49].
sampling point under the different wave conditions by applying Equation (1).Select the minimum safety factor value for each point under varying wave conditions.
(3) Analyze and process the data through ArcGIS 10.7, incorporating interpolation methods to generate a liquefaction susceptibility zoning map for the study area.Subsequently, conduct an in-depth analysis of the zoning results.

Study Area
Numerous studies indicate that nearshore silty seabed are typical areas prone to liquefaction, posing a significant threat to the stability of nearshore marine engineering [7,[45][46][47].In this context, the Chengdao Island area, where silty sediments are widely distributed, has been selected as the study area.The Chengdao region is situated to the north of the Yellow River Estuary in China, exhibiting a southwest-to-northeast topographical gradient, as depicted in Figure 4.This area primarily consists of two subdeltas [47].The first subdelta was formed when the Yellow River was redirected to the Shenxiangou flow path from 1953 to 1964, while the second subdelta was established when the river was diverted to the Diaokou flow path from 1964 to 1976 [48].The underwater Yellow River Estuary spans an area of approximately 3000 km 2 .A distinctive feature of this region is its remarkably gentle gradient, typically less than 0.4°, and an exceptionally low slope drop ranging from 1% to 5% [48].This makes the area highly susceptible to storm surge hazards [1].Between 1949 and 2005, more than 150 storm surge events occurred in the region, with an average of two to three events per year, of which 36 storm surges caused geological disasters, most of which occurred in the summer and fall [49].The Chengdao region holds significant engineering importance as it encompasses the Chengdao oil field, which is a crucial oil-producing area.Numerous engineering facilities, including oil platforms, submarine pipelines, and breakwaters, have been established in this region [51][52][53].However, the marine engineering environment in the Chengdao The Chengdao region holds significant engineering importance as it encompasses the Chengdao oil field, which is a crucial oil-producing area.Numerous engineering facilities, including oil platforms, submarine pipelines, and breakwaters, have been established in this region [51][52][53].However, the marine engineering environment in the Chengdao region is notably complex, with frequent occurrences of seabed liquefaction leading to numerous destabilization accidents involving offshore engineering facilities.For instance, in November 2003, a seabed deformation sliding event occurred near the oil production platform CB12B, which led to the rupture and interruption of two submarine cables.In May 2009, the CB25A-CB25B subsea crude oil transmission pipeline fractured near the endpoint of CB25A [52].The investigation data indicated that the accident site was approximately 50-60 m from the CB25A platform, 80 m from the intersection of the Center II-CB1A subsea pipeline, and 200 m from the Center II-CB1A submarine cable.Furthermore, a capsizing incident occurred on the Shengli Operation No. 3 maintenance platform near the CB22C well group platform of Chengdu Island Oilfield in 2010, putting 36 people in danger [54].Unfortunately, two individuals lost their lives [7].The Shengli Operation No. 3 platform and the adjacent seabed were adversely affected by intense wind and waves induced by a typhoon.This led to the liquefaction of the soft substrate beneath the platform, and as the liquefaction zone gradually expanded, it triggered instability, ultimately causing the capsizing of the platform [7].

Cyclic Stress Ratio under Various Wave Return Periods
The Chengdao region is influenced by prevailing monsoon winds, experiencing dominant north winds in the winter and south winds in the summer [52].Wave patterns in this area are primarily controlled by surface winds, characterized by short wind zones, rapid wave growth, high wave heights, short wave periods, and limited attenuation distances.These characteristics are associated with the semi-enclosed nature of the Chengdao region, which is connected solely to the Yellow Sea through the narrow Bohai Strait.Furthermore, the engineering geological conditions of the underwater delta impact wave propagation, limiting the distance waves can travel [55].Under normal sea conditions, wave heights typically do not exceed 1.5 m.However, during extreme sea conditions, wave heights can surpass 5.8 m, and the measured maximum current velocity can reach 1.5 m/s [1].It is estimated that the extreme water level may increase by 1.8-3.2m per year, with an average increase of 2.2 m in the Chengdao region [46].These extreme conditions are primarily generated by weather phenomena such as typhoons, cyclones, and cold waves [56].Observational and numerical simulation data provide a breakdown of wave information for different wind conditions in the Chengdao region, as shown in Table 3.According to linear wave theory, an extreme wave event can be considered to consist of many different waves with characteristics, and the wave components can be characterized by wave period, wavelength, and wave height [31].The peak wave stress in the study area is calculated under the wave conditions of a 5-year, 10-year, and 25-year return period to analyze the distribution of cyclic stress within the study area across varying wave conditions and sampling depths.Considering that the maximum liquefaction depth in the Chengdao region is less than 15 m, CSRs for three different wave conditions at sampling depths of 5 m, 10 m, and 15 m were further calculated [13].The interpolation method was utilized to generate a zonal map of the peak wave pressures under different wave conditions (ρ w = 1025 kg/m 3 is the density of seawater and g = 9.8 m/s 2 is the gravitational acceleration), as shown in Figure 5.   Overall, a discernible spatial regularity characterizes the distribution of CSR, with contours closely aligned with the coastline.The lowest CSR value is observed under the wave conditions of a 5-year return period, situated at a depth of 15 m in the southeastern offshore area of the study region, with a value of 0.0004.Conversely, the highest CSR value is found under the wave conditions of a 25-year return period, located at a depth of 5 m in the northeastern offshore area, with a value of 0.305.A deeper analysis of the CSR distribution at different depths but the same wave conditions reveals a general increasing trend from nearshore to offshore areas.However, a distinct spatial pattern emerges in the northern nearshore part of the Chengdao region, characterized by a concentration of higher CSR values where numerous offshore platforms are located.Additionally, CSR displays a negative correlation with sampling depth, ranging from 0.03 to 0.31 at 5 m depth, 0.01 to 0.23 at 10 m depth, and 0.0004 to 0.16 at 15 m depth.This indicates that deeper sediment experiences a reduced impact from wave loading.By analyzing the distribution of CSR at the same depth but under different wave conditions, it can be found that higher wave conditions lead to an increased CSR, which, in turn, amplifies the likelihood of seabed liquefaction.The precision of the interpolation results is primarily attributed to the limited distribution of interpolation points.However, regarding the interpolation of CSR, the known data points are mostly distributed in nearshore areas, with fewer points in offshore areas.As a result, the accuracy of interpolation results in offshore areas will be lower than in nearshore areas.

Cyclic Resistance Ratio of Sediments
The sediment types in the Chengdao region can be categorized into four types, including clay, silty clay, silt, and silty sand, as shown in Figure 6 [52].In the Chengdao oil field beach and shallow marine stratigraphy, silt is the main sediment type [51].However, the nearshore area experiences strong hydrodynamic forces and sediment redeposition, leading to a more intricate sediment composition [49].This results in an alternation between coarser and finer sediment types, with silt predominating locally.Moreover, the distribution of silt in the far shore area slightly surpasses that of clay [46,57].
is found under the wave conditions of a 25-year return period, located at a depth of 5 m in the northeastern offshore area, with a value of 0.305.A deeper analysis of the CSR distribution at different depths but the same wave conditions reveals a general increasing trend from nearshore to offshore areas.However, a distinct spatial pattern emerges in the northern nearshore part of the Chengdao region, characterized by a concentration of higher CSR values where numerous offshore platforms are located.Additionally, CSR displays a negative correlation with sampling depth, ranging from 0.03 to 0.31 at 5 m depth, 0.01 to 0.23 at 10 m depth, and 0.0004 to 0.16 at 15 m depth.This indicates that deeper sediment experiences a reduced impact from wave loading.By analyzing the distribution of CSR at the same depth but under different wave conditions, it can be found that higher wave conditions lead to an increased CSR, which, in turn, amplifies the likelihood of seabed liquefaction.The precision of the interpolation results is primarily attributed to the limited distribution of interpolation points.However, regarding the interpolation of CSR, the known data points are mostly distributed in nearshore areas, with fewer points in offshore areas.As a result, the accuracy of interpolation results in offshore areas will be lower than in nearshore areas.

Cyclic Resistance Ratio of Sediments
The sediment types in the Chengdao region can be categorized into four types, including clay, silty clay, silt, and silty sand, as shown in Figure 6 [52].In the Chengdao oil field beach and shallow marine stratigraphy, silt is the main sediment type [51].However, the nearshore area experiences strong hydrodynamic forces and sediment redeposition, leading to a more intricate sediment composition [49].This results in an alternation between coarser and finer sediment types, with silt predominating locally.Moreover, the distribution of silt in the far shore area slightly surpasses that of clay [46,57].The sediment strength in the Chengdao region exhibits a notable pattern, which can be observed in Figure 7, depicting the distribution of undrained shear strength in depth ranges of 0-5 m, 5-10 m, and 10-15 m.The exact calculation results are presented in Appendix A. Because of more rapid sedimentation, the shallow sediments usually do not reach complete consolidation, so their undrained shear strengths are lower, and the range The sediment strength in the Chengdao region exhibits a notable pattern, which can be observed in Figure 7, depicting the distribution of undrained shear strength in depth ranges of 0-5 m, 5-10 m, and 10-15 m.The exact calculation results are presented in Appendix A. Because of more rapid sedimentation, the shallow sediments usually do not reach complete consolidation, so their undrained shear strengths are lower, and the range of variation is smaller [57].On the contrary, the undrained shear strength of the deeper sediments roughly shows a gradual increase along the delta and outward [57].This phenomenon results from a gradual reduction in sedimentation rates and a progressive thinning of sediment particles.Furthermore, regional trends in the sediment water content, the pore ratio, and the liquid-plastic limit follow a similar pattern, with these parameters gradually increasing with water depth and distance from the Diaokou inlet [48].These characteristics collectively create conditions characterized by low consolidation, a high water content, a high pore ratio, and low strength, which increases the risk of geological disasters such as sliding, thixotropy, and liquefaction.
ning of sediment particles.Furthermore, regional trends in the sediment water content, the pore ratio, and the liquid-plastic limit follow a similar pattern, with these parameters gradually increasing with water depth and distance from the Diaokou inlet [48].These characteristics collectively create conditions characterized by low consolidation, a high water content, a high pore ratio, and low strength, which increases the risk of geological disasters such as sliding, thixotropy, and liquefaction.To calculate CRR, this study used the cyclic triaxial shear test data.The method for this experiment was the consolidation of undrained shear, with an isotropic consolidation process (i.e., axial consolidation pressure ( ) equaled lateral consolidation pressure ( )).After consolidation was completed, the confining pressure and axial pressure remained unchanged, while a sinusoidal axial load with a frequency of 1 Hz was applied.This caused cyclic shear stress on a 45° inclined plane within the soil sample, simulating the shear forces induced by wave loading within the seabed.The same set of soil samples was subjected to at least three different amplitudes of dynamic loading to obtain a relationship curve between the number of cycles and the dynamic stress ratio.Data from 29 groups of cyclic triaxial tests conducted between 1998 and 2017 with sediment samples in the study area were collected and organized.The relationship between the dynamic stress ratio and the number of cycles at various depths was fitted to calculate the cyclic resistance ratio (CRR) for each sampling point.The exact calculation results are provided in Appendix B. The minimum CRR value at each point was selected as the final CRR. Figure 8 only presents the relationship between the dynamic stress ratio and the number of cycles for a portion of different sampling points, along with the method for determining CRR.However, because of the nonuniform distribution of 29 sampling points, it was necessary to utilize Empirical Bayesian Interpolation to extrapolate the CRR distribution across the continuous surface of the study area, as illustrated in Figure 9.To calculate CRR, this study used the cyclic triaxial shear test data.The method for this experiment was the consolidation of undrained shear, with an isotropic consolidation process (i.e., axial consolidation pressure (σ 1 ) equaled lateral consolidation pressure (σ 3 )).After consolidation was completed, the confining pressure and axial pressure remained unchanged, while a sinusoidal axial load with a frequency of 1 Hz was applied.This caused cyclic shear stress on a 45 • inclined plane within the soil sample, simulating the shear forces induced by wave loading within the seabed.The same set of soil samples was subjected to at least three different amplitudes of dynamic loading to obtain a relationship curve between the number of cycles and the dynamic stress ratio.Data from 29 groups of cyclic triaxial tests conducted between 1998 and 2017 with sediment samples in the study area were collected and organized.The relationship between the dynamic stress ratio and the number of cycles at various depths was fitted to calculate the cyclic resistance ratio (CRR) for each sampling point.The exact calculation results are provided in Appendix B. The minimum CRR value at each point was selected as the final CRR. Figure 8 only presents the relationship between the dynamic stress ratio and the number of cycles for a portion of different sampling points, along with the method for determining CRR.However, because of the nonuniform distribution of 29 sampling points, it was necessary to utilize Empirical Bayesian Interpolation to extrapolate the CRR distribution across the continuous surface of the study area, as illustrated in Figure 9.
CRR mainly reflects the ability of sediments in the study area to resist liquefaction.In general, variations in CRR within the study area are relatively minor and closely linked to sediment type.The northwestern region predominantly features silt sediments, with most of this area exhibiting CRR values exceeding 0.27, indicating greater resistance to liquefaction.Furthermore, CRR correlates with undrained shear strength, which is notably lower in the southeastern region at different depths, resulting in a CRR within the range of 0.2 to 0.22.Unlike the distribution of CSR, CRR does not follow a coastline-based distribution pattern.Instead, they display one high CRR center and three low CRR centers arranged in a concentric circle formation.These interpolation results are influenced by the distribution of sampling points and the selection of the interpolation method.
A comprehensive analysis of CSR, CRR, and safety factors at various sampling depths at the same location reveals interesting patterns.Under the same wave condition, CSR and CRR exhibit a gradual decrease as sampling depth increases.The reduction in CSR (Figure 5) is associated with the attenuation of stress as waves propagate downward through the sediment layers.Furthermore, the reduction in the CRR trend could be attributed to the robust hydrodynamic forces in the nearshore area, which generate intense shear stresses, inhibiting the settling of sediment particles toward the seabed and consequently resulting in a decline in CRR with greater depth.However, as sampling depth increases, the magnitude of reduction in CSR surpasses that of CRR, resulting in an increase in the safety factor and a decreased susceptibility to wave-induced seabed liquefaction.This implies that the depth of sediment may significantly influence the susceptibility to wave-induced seabed liquefaction, leading to a reduced susceptibility to disasters.CRR mainly reflects the ability of sediments in the study area to resist liquefaction.In general, variations in CRR within the study area are relatively minor and closely linked to sediment type.The northwestern region predominantly features silt sediments, with most of this area exhibiting CRR values exceeding 0.27, indicating greater resistance to liquefaction.Furthermore, CRR correlates with undrained shear strength, which is notably lower in the southeastern region at different depths, resulting in a CRR within the range of 0.2 to 0.22.Unlike the distribution of CSR, CRR does not follow a coastline-based distribution pattern.Instead, they display one high CRR center and three low CRR centers arranged in a concentric circle formation.These interpolation results are influenced by the distribution of sampling points and the selection of the interpolation method.
A comprehensive analysis of CSR, CRR, and safety factors at various sampling depths at the same location reveals interesting patterns.Under the same wave condition, CSR and CRR exhibit a gradual decrease as sampling depth increases.The reduction in CSR (Figure 5) is associated with the attenuation of stress as waves propagate downward through the sediment layers.Furthermore, the reduction in the CRR trend could be   CRR mainly reflects the ability of sediments in the study area to resist liquefaction.In general, variations in CRR within the study area are relatively minor and closely linked to sediment type.The northwestern region predominantly features silt sediments, with most of this area exhibiting CRR values exceeding 0.27, indicating greater resistance to liquefaction.Furthermore, CRR correlates with undrained shear strength, which is notably lower in the southeastern region at different depths, resulting in a CRR within the range of 0.2 to 0.22.Unlike the distribution of CSR, CRR does not follow a coastline-based distribution pattern.Instead, they display one high CRR center and three low CRR centers arranged in a concentric circle formation.These interpolation results are influenced by the distribution of sampling points and the selection of the interpolation method.
A comprehensive analysis of CSR, CRR, and safety factors at various sampling depths at the same location reveals interesting patterns.Under the same wave condition, CSR and CRR exhibit a gradual decrease as sampling depth increases.The reduction in CSR (Figure 5) is associated with the attenuation of stress as waves propagate downward through the sediment layers.Furthermore, the reduction in the CRR trend could be

Susceptibility Assessment
The safety factor for wave-induced seabed liquefaction at various known depths of sampling points was calculated following the previously outlined procedure.These calculations were conducted under different wave conditions corresponding to 5-year, 10-year, and 25-year return periods.We considered the smallest safety factor among those calculated at different known depths as the safety factor for wave-induced seabed liquefaction at that specific point, as detailed in Table 4. Subsequently, we utilized EBK in ArcGIS 10.7 to generate a distribution map of the safety factor of wave-induced seabed liquefaction, as shown in Figure 10.[48].The safety factor and liquefaction probability of sediments against liquefaction were calculated under different wave conditions, taking into account the variability in wave and soil parameters.Based on the Chengdao region, a division of the safety factor was proposed.If the probability is controlled within 5%, the safety factor must be greater than 1.8, i.e., the area is a non-liquefaction zone when FS > 1.8 and a liquefaction susceptibility zone when FS < 1.8 [48].Broadly, the safety factor in the western part of the study area surpasses that in the eastern part, showing substantial spatial variation.The extent of the susceptibility zone expands with more severe wave conditions.Under the wave condition in a 5-year return period, the safety factor in the study area ranges from 0.93 to 14.95, with the liquefaction susceptibility zone encompassing 64.46% of the total study area.Under the wave condition in a 10-year return period, the safety factor ranges from 0.839 to 14.48, covering 65.07% of the study area.Under the wave condition in a 25-year return period, the safety factor varies from 0.76 to 13.22, with the susceptibility zone encompassing 68.89% of the study area.The distribution of the safety factor shows a concentric circle structure centered on the sampling point, with a high safety factor center occurring in the northwestern area and a low safety factor center in the southeastern area.The high safety factor center was located in the northeast of the study area,   [48].The safety factor and liquefaction probability of sediments against liquefaction were calculated under different wave conditions, taking into account the variability in wave and soil parameters.Based on the Chengdao region, a division of the safety factor was proposed.If the probability is controlled within 5%, the safety factor must be greater than 1.8, i.e., the area is a non-liquefaction zone when FS > 1.8 and a liquefaction susceptibility zone when FS < 1.8 [48].Broadly, the safety factor in the western part of the study area surpasses that in the eastern part, showing substantial spatial variation.The extent of the susceptibility zone expands with more severe wave conditions.Under the wave condition in a 5-year return period, the safety factor in the study area ranges from 0.93 to 14.95, with the liquefaction susceptibility zone encompassing 64.46% of the total study area.Under the wave condition in a 10-year return period, the safety factor ranges from 0.839 to 14.48, covering 65.07% of the study area.Under the wave condition in a 25-year return period, the safety factor varies from 0.76 to 13.22, with the susceptibility zone encompassing 68.89% of the study area.The distribution of the safety factor shows a concentric circle structure centered on the sampling point, with a high safety factor center occurring in the northwestern area and a low safety factor center in the southeastern area.The high safety factor center was located in the northeast of the study area, and its range gradually expanded and developed significantly with an enhancement in wind and wave action, while the range of the high safety factor center in the northwest gradually decreased and the range of the low safety factor center in the southeast gradually expanded (Figure 10).Notably, the safety factor in the western part of the study area under the 25-year return period wave conditions ranges between 0.98 and 1.0.There is an area extending southwestward along the sampling points, and its range changes most significantly, from 18.73% for the wave condition in a 5-year return period to 43.11% for the wave condition in a 25-year return period.Combining the distribution of sampling points and susceptibility areas, it is observed that most sampling points are located in high-safety areas.However, there are still some sampling points with lower safety factors.This is of concern because some important submarine engineering facilities are in these zones.For example, submarine cables from KD48 to KD34C, near borehole K1-1, are buried in an area with the lowest safety factor.This site is particularly susceptible to wave-induced seabed liquefaction under different wave conditions.
The method can realize the preliminary susceptibility assessment of wave-induced seabed liquefaction, which provides a reference value for the construction of marine engineering and maintenance of facilities, but there remains room for its enhancement.Firstly, the accuracy of susceptibility assessments in areas with limited or sparse data points is compromised because of their spatial distribution and the lack of data.Secondly, the current method does not fully account for the influence of sampling depth changes on the safety factor, making it unable to provide a comprehensive assessment of the three-dimensional spatial susceptibility to wave-induced seabed liquefaction.Consequently, the evaluation results may deviate from the actual conditions.To enhance the accuracy of this method, future research can explore wave-induced seabed liquefaction principles through extensive indoor experiments, take into account the impact of sampling depth on susceptibility more comprehensively, and potentially employ machine learning-based models to achieve a more precise evaluation.These improvements could lead to more reliable results in evaluating susceptibility and better support marine engineering construction and facility maintenance decisions.

Conclusions
This study modifies the formula for calculating the safety factor, which makes the method suitable for studying the susceptibility assessment of seafloor sediment liquefaction under wave loads.This method employs the safety factor as a primary evaluation criterion and combines an analysis of the cyclic stress ratio (CSR) and the cyclic resistance ratio (CRR) to provide a comprehensive assessment.The Chengdao region in China, characterized by extensive silt distribution, was selected as the study area.To create a regional evaluation zoning map for the susceptibility of wave-induced seabed liquefaction, Empirical Bayesian Kriging (EBK) in the Arc Geographic Information System 10.7 (ArcGIS 10.7) was utilized.The following conclusions can be drawn from this study: (1) A method for the susceptibility assessment of wave-induced seabed liquefaction was modified using the safety factor.This safety factor was determined by employing CSR to account for the wave characteristics in the study area and CRR to reflect the sediment properties.This method was then applied to the Chengdao region, where a substantial amount of data was collected to evaluate the susceptibility of waveinduced seabed liquefaction.The EBK method was further utilized to create zoning maps for the susceptibility assessment of wave-induced seabed liquefaction in the Chengdao region under varying wave conditions.(2) CSR primarily characterizes the wave conditions in the study area and is a crucial factor for the susceptibility assessment of wave-induced seabed liquefaction.In this study, CSR values are calculated using basic wave elements based on linear wave theory.To gain a deeper insight into the impact of waves on the susceptibility assessment of wave-induced seabed liquefaction and enhance the accuracy of the evaluation, we believe that multiple extreme wave conditions should be calculated for CSR.(3) CRR primarily represents the sediment characteristics of the study area, and its calculation method should be tailored to the sediment types in the study area.This study presents several testing methods, such as SPT (Standard Penetration Test), CPT (Cone Penetration Test), measurement of V S (shear wave velocity), and BPT, with a specific focus on the cyclic triaxial test, which is convenient for rapid computation, and the data are easily obtainable, making it applicable to various types of sediments.(4) The safety factor serves as an effective indicator of wave-induced seabed liquefaction susceptibility assessment in the study area.The safety factor integrates the wave characteristics and sediment properties of the study area.Combined with spatial interpolation methods, it allows for a continuous assessment of wave-induced seabed liquefaction susceptibility across the study area.The specific categorization of susceptibility levels also needs to consider the characteristics of the study area and liquefaction probabilities.Through the calculation of safety factors, it was observed that CSR is more sensitive to changes in sampling depth compared with CRR, exerting a greater impact on the safety factor.Therefore, in future studies, it is crucial to consider the sampling depth as a key reference factor for a more comprehensive evaluation of wave-induced seabed liquefaction susceptibility in three-dimensional space.
Overall, the methodology of this study provides a tool for the susceptibility assessment of wave-induced seabed liquefaction, setting the stage for subsequent risk assessment of such disasters.This method can be directly applied to seabed areas in the study region where the sediment type is silt, providing a spatial reference for stability assessment and site selection of marine engineering facilities within the study area.When this method is applied to other types of sediment in the study area, the formula of the cyclic resistance ratio needs to be adjusted according to the characteristics of the sediment.However, we also need to recognize that there is still room for improvement in this method, especially in the amount and distribution of data, and more in-depth studies are needed to enhance the accuracy and reliability of the assessment.Marine Science Data Center for furnishing the 2018 annual wave element data for the Chengdao area.We acknowledge with gratitude the provision of cyclic triaxial test data from the sampling points, the comprehensive physical and mechanical property data of sediment in the Chengdao area by the Sinopec Petroleum Engineering Corporation, conducted in collaboration with Ocean University of China, the First Institute of Oceanography, Ministry of Natural Resources, Qingdao Huanhai Marine Engineering Survey Institute, and Sinopec Petroleum Engineering Corporation.Additionally, we would like to extend our appreciation to Mengjiao Liu for her invaluable assistance in the design and production of the figures.
Conflicts of Interest: Jinkun Liu and Fang Hou were employed by the Sinopec Petroleum Engineering Corporation.The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Figure 4 .
Figure 4. Overview of the Chengdao region and the study zone.The large image shows the location of the Chengdao region and the study zone, with arrows indicating the direction of ocean currents in the study area; the small image in the upper right corner represents the study zone, where yellow indicates the distribution of offshore platforms in existing publicly available data (images modified from Zhang et al. [47] and Wang et al. [50]).

Figure 4 .
Figure 4. Overview of the Chengdao region and the study zone.The large image shows the location of the Chengdao region and the study zone, with arrows indicating the direction of ocean currents in the study area; the small image in the upper right corner represents the study zone, where yellow indicates the distribution of offshore platforms in existing publicly available data (images modified from Zhang et al. [47] and Wang et al. [50]).
ditions and sampling depths.Considering that the maximum liquefaction depth in the Chengdao region is less than 15 m, CSRs for three different wave conditions at sampling depths of 5 m, 10 m, and 15 m were further calculated [13].The interpolation method was utilized to generate a zonal map of the peak wave pressures under different wave conditions ( 1025 kg/m is the density of seawater and  9.8 m/s is the gravitational acceleration), as shown in Figure 5.

Figure 5 .
Figure 5. CSR of the study area (the small image in Figure 4) under different wave conditions at different sampling depths: (a) CSR of the wave condition in a 5-year return period at 5 m depth; (b) CSR of the wave condition in a 5-year return period at 10 m depth; (c) CSR of the wave condition in a 5-year return period at 15 m depth; (d) CSR of the wave condition in a 10-year return period at 5 m depth; (e) CSR of the wave condition in a 10-year return period at 10 m depth; (f) CSR of the wave condition in a 10-year return period at 15 m depth; (g) CSR of the wave condition in a 25-year return period at 5 m depth; (h) CSR of the wave condition in a 25-year return period at 10 m depth; and (i) CSR of the wave condition in a 25-year return period at 15 m depth.

Figure 5 .
Figure 5. CSR of the study area (the small image in Figure 4) under different wave conditions at different sampling depths: (a) CSR of the wave condition in a 5-year return period at 5 m depth; (b) CSR of the wave condition in a 5-year return period at 10 m depth; (c) CSR of the wave condition in a 5-year return period at 15 m depth; (d) CSR of the wave condition in a 10-year return period at 5 m depth; (e) CSR of the wave condition in a 10-year return period at 10 m depth; (f) CSR of the wave condition in a 10-year return period at 15 m depth; (g) CSR of the wave condition in a 25-year return period at 5 m depth; (h) CSR of the wave condition in a 25-year return period at 10 m depth; and (i) CSR of the wave condition in a 25-year return period at 15 m depth.

Figure 6 .
Figure 6.Distribution of sediment types in the study area.Blue dots represent sample points (data from Liu et al. [51] and Liu et al. [52]).

Figure 6 .
Figure 6.Distribution of sediment types in the study area.Blue dots represent sample points (data from Liu et al. [51] and Liu et al. [52]).

J 29 Figure 8 .
Figure 8.The relationship between the dynamic stress ratio of the selected sampling points and the number of cycles.The dotted line represents the dynamic stress ratio corresponding to the cycle number of 100, which is considered the CRR of the sampling point.

Figure 9 .
Figure 9.The interpolated distribution represents the minimum values of known CRRs for the sampling points within the 15 m depth range.

Figure 8 .
Figure 8.The relationship between the dynamic stress ratio of the selected sampling points and the number of cycles.The dotted line represents the dynamic stress ratio corresponding to the cycle number of 100, which is considered the CRR of the sampling point.

J 29 Figure 8 .
Figure 8.The relationship between the dynamic stress ratio of the selected sampling points and the number of cycles.The dotted line represents the dynamic stress ratio corresponding to the cycle number of 100, which is considered the CRR of the sampling point.

Figure 9 .
Figure 9.The interpolated distribution represents the minimum values of known CRRs for the sampling points within the 15 m depth range.

Figure 9 .
Figure 9.The interpolated distribution represents the minimum values of known CRRs for the sampling points within the 15 m depth range.

Figure 10 .
Figure 10.Grading evaluation of wave-induced seabed liquefaction susceptibility in the Chengdao region under wave conditions in a 5-year return period (a), a 10-year return period (b), and a 25year return period (c), respectively.

Figure 10
Figure 10 reveals the susceptibility assessment result of wave-induced seabed liquefaction in the Chengdao region.The zoning of susceptibility refers to previous research by Chang et al.[48].The safety factor and liquefaction probability of sediments against liquefaction were calculated under different wave conditions, taking into account the variability in wave and soil parameters.Based on the Chengdao region, a division of the safety factor was proposed.If the probability is controlled within 5%, the safety factor must be greater than 1.8, i.e., the area is a non-liquefaction zone when FS > 1.8 and a liquefaction susceptibility zone when FS < 1.8[48].Broadly, the safety factor in the western part of the study area surpasses that in the eastern part, showing substantial spatial variation.The extent of the susceptibility zone expands with more severe wave conditions.Under the wave condition in a 5-year return period, the safety factor in the study area ranges from 0.93 to 14.95, with the liquefaction susceptibility zone encompassing 64.46% of the total study area.Under the wave condition in a 10-year return period, the safety factor ranges from 0.839 to 14.48, covering 65.07% of the study area.Under the wave condition in a 25-year return period, the safety factor varies from 0.76 to 13.22, with the susceptibility zone encompassing 68.89% of the study area.The distribution of the safety factor shows a concentric circle structure centered on the sampling point, with a high safety factor center occurring in the northwestern area and a low safety factor center in the southeastern area.The high safety factor center was located in the northeast of the study area,

Figure 10 .
Figure 10.Grading evaluation of wave-induced seabed liquefaction susceptibility in the Chengdao region under wave conditions in a 5-year return period (a), a 10-year return period (b), and a 25-year return period (c), respectively.

Figure 10
Figure 10 reveals the susceptibility assessment result of wave-induced seabed liquefaction in the Chengdao region.The zoning of susceptibility refers to previous research by Chang et al.[48].The safety factor and liquefaction probability of sediments against liquefaction were calculated under different wave conditions, taking into account the variability in wave and soil parameters.Based on the Chengdao region, a division of the safety factor was proposed.If the probability is controlled within 5%, the safety factor must be greater than 1.8, i.e., the area is a non-liquefaction zone when FS > 1.8 and a liquefaction susceptibility zone when FS < 1.8[48].Broadly, the safety factor in the western part of the study area surpasses that in the eastern part, showing substantial spatial variation.The extent of the susceptibility zone expands with more severe wave conditions.Under the wave condition in a 5-year return period, the safety factor in the study area ranges from 0.93

Table 1 .
Comparison of advantages and disadvantages of various field tests for CRR.

Table 2 .
Summary of methods for calculating CRR.
Note: H 1/10 is the effective wave height.

Table 4 .
CSRs, CRRs, and the safety factor minima calculated for 29 core sampling points.

Table A4 .
Sediment cyclic triaxial torsion shear test data in the research area and CRR.