Experimental Study on Mechanism of Wave-Induced Liquefaction of Sand-Clay Seabed

: In this study, a series of laboratory experiments for the response of wave induced clay-sand seabed were carried out to clarify the mechanism of liquefaction of clayey seabed. The experiments were conducted in an 80 m long wave ﬂume. In the tests, the sand-clay beds were mixed with various clay contents (CC) from 0.5% to 15% and were tested for given wave conditions. The pore water pressure and the water elevation were measured in each test. Soil properties tests and scanning electron microscope (SEM) experiments on di ﬀ erent seabed samples were carried out to further explore the mechanism of liquefaction. The experimental results indicated that the amplitude and accumulation of the excess pore water pressure (EPP) varied with di ﬀ erent CC in the sand-clay bed. With the introduction of CC, micro-structure and properties (such as permeability and compressibility) of bed soils changed. Sand-clay bed presented more susceptibility to liquefy compared with pure sand bed. CC promoted seabed liquefaction, even if the added amount was very small (CC is 0.5%), however when CC exceeded a certain value (10% in this study), the mixed bed will not be liqueﬁed. This phenomenon can be well explained by the micro-structure of sand-clay bed. CC within a sandy seabed, does not only a ﬀ ect the permeability, but also change the compressibility of seabed soils. For example, the microfabric of seabed vulnerable to liquefaction is loose. Clay aggregations generally gathered at the sand particle contact points. This microfabric is easily compressed under wave loads and allowed pore water to ﬂow, resulting in the accumulation of pore water pressure. On the other hand, the microfabric of seabed that was resistant to liquefaction appeared to be more compact. Due to clay-ﬁlled gaps between the sand particles, the pore water is more di ﬃ cult to ﬂow when seabed was compressed. Furthermore, the tendency of seabed liquefaction is closely related to CC.


Introduction
Wave-induced liquefaction of seabed is an important phenomenon in coastal engineering, because it will cause reduction of effective stress between the soil grains, and then uneven settlement of foundation caused by the reduction will result in the damage of structures. As reported in the literature, there are two main mechanisms of seabed liquefaction under waves: residual vs. momentary liquefaction. Residual liquefaction is associated with the accumulated pore-water pressure gradient, whereas the momentary liquefaction is associated with the amplitude of pressure oscillations [1]. For most of the coastal zones, seabed soils are generally classified into three categories, i.e., sand, silt Based on the aforementioned researches of seabed samples by SEM, a hypothesis was proposed to explain the mechanism of clayey seabed liquefaction.
Natural seabed commonly contains various types of soil particles, especially clay, and a small amount of clay in the seabed could lead to a significant impact on its dynamic behavior. To date, most of the researches mainly focused on the critical value of CC that liquefaction would not occur when the CC exceeded it [7,[20][21][22]. However, there are a few aspects concerning whether CC will affect seabed when it is very small. Both Gratchev et al. [7] and Kirca et al. [20] tested the seabed with 5% CC (sand-clay and silt-clay seabed), but they didn't consider the seabed with less than 5% CC. It is helpful to understand the liquefaction mechanism of clayey soils with studying seabed with very low CC. Therefore, this study endeavored to understand the influence of CC on the dynamic response of sand-clay bed. Laboratory experiments were carried out by changing of the concentration of clay from 0.5% to 15% in the mixed bed. The pore water pressure and the corresponding wave height were measured and analyzed.

Wave Flume
The experiments were conducted in a wave flume, 80.0 m long, 0.5 m wide, and 1.5 m high, with a pit filled with mixed clay-sand seabed for the investigation of seabed response ( Figure 1). Waves were generated by a piston-type wave-generator located at upstream end of the flume. At the near side of the wave paddle, a fixed bed with 1:10 slope adjacent to a 7.5 m long platform was set in front of a soil trench. Behind the sand trench, a 6 m long platform was set to retain the trench. The soil trench, 0.25 m deep, 1.0 m long and 0.5 m wide, is used to contain the mixed clay-sand seabed. A sloping wave absorber was located at downstream end to eliminate wave reflection. In total, 3 wave gauges (Yu Fan Technology Co, Chengdu) were set up in the flume, one is in front of the trench to measure the incident wave height, other one is behind the trench to measure the wave height after passing through the soil bed and the third gauge is located above sensors. J. Mar. Sci. Eng. 2020, 8,66 3 of 16 seabed samples by SEM, a hypothesis was proposed to explain the mechanism of clayey seabed liquefaction. Natural seabed commonly contains various types of soil particles, especially clay, and a small amount of clay in the seabed could lead to a significant impact on its dynamic behavior. To date, most of the researches mainly focused on the critical value of CC that liquefaction would not occur when the CC exceeded it [7,[20][21][22]. However, there are a few aspects concerning whether CC will affect seabed when it is very small. Both Gratchev et al. [7] and Kirca et al. [20] tested the seabed with 5% CC (sand-clay and silt-clay seabed), but they didn't consider the seabed with less than 5% CC. It is helpful to understand the liquefaction mechanism of clayey soils with studying seabed with very low CC. Therefore, this study endeavored to understand the influence of CC on the dynamic response of sand-clay bed. Laboratory experiments were carried out by changing of the concentration of clay from 0.5% to 15% in the mixed bed. The pore water pressure and the corresponding wave height were measured and analyzed.

Wave Flume
The experiments were conducted in a wave flume, 80.0 m long, 0.5 m wide, and 1.5 m high, with a pit filled with mixed clay-sand seabed for the investigation of seabed response ( Figure 1). Waves were generated by a piston-type wave-generator located at upstream end of the flume. At the near side of the wave paddle, a fixed bed with 1:10 slope adjacent to a 7.5 m long platform was set in front of a soil trench. Behind the sand trench, a 6 m long platform was set to retain the trench. The soil trench, 0.25 m deep, 1.0 m long and 0.5 m wide, is used to contain the mixed clay-sand seabed. A sloping wave absorber was located at downstream end to eliminate wave reflection. In total, 3 wave gauges (Yu Fan Technology Co, Chengdu) were set up in the flume, one is in front of the trench to measure the incident wave height, other one is behind the trench to measure the wave height after passing through the soil bed and the third gauge is located above sensors.

Soil Trench
Pore-water pressure sensors (HC-25, RuiHeng Chang Tai Technology Co, Beijing) were located in the middle of the trench at 3 different depths (Figure 2), z = 0.09, 0.16 and 0.23 m, where z was the vertical distance measured downward from the soil surface. Soil samples used for experiments were made indoor, by mixing natural sand and commercial kaolin with the determined content ratios.
Preparation of the sand-clay mixtures was as follow: sand-clay mixtures were made by using a handheld mixer that allowed mixing without any significant grinding of the mixture. First, fine sand was poured into the mixing tank in a nearly dry condition. Then, clay was separately poured into a small tank and mixed with appropriate water until it reached the liquid state. Subsequently, the liquid clay was added into the mixing tank and running the mixer to form a uniformed soil sample. Finally, the prepared soil sample was placed in the trench of the flume.

Soil Trench
Pore-water pressure sensors (HC-25, RuiHeng Chang Tai Technology Co, Beijing) were located in the middle of the trench at 3 different depths (Figure 2), z = 0.09, 0.16 and 0.23 m, where z was the vertical distance measured downward from the soil surface. Soil samples used for experiments were made indoor, by mixing natural sand and commercial kaolin with the determined content ratios.
Preparation of the sand-clay mixtures was as follow: sand-clay mixtures were made by using a handheld mixer that allowed mixing without any significant grinding of the mixture. First, fine sand was poured into the mixing tank in a nearly dry condition. Then, clay was separately poured into a small tank and mixed with appropriate water until it reached the liquid state. Subsequently, the liquid clay was added into the mixing tank and running the mixer to form a uniformed soil sample. Finally, the prepared soil sample was placed in the trench of the flume.

Soil Mechanics Test
The particle size distribution of materials was obtained by the Mastersizer 3000 laser particle size analyzer ( Figure 3). The of two kinds of soils are 0.163 mm (fine sand) and 0.009 mm (clay). The physical properties of seabeds are summarized in Table 1, CC was defined as the ratio of the weight of clay to the weight of sand-clay mixture, where is the density of soil, k is the permeability coefficient, is the rate of water content, e is the void ratio of the soil, is dry density of soil, is maximum dry density and is minimum dry density, is the relative density and is given by Equation (1).
is the coefficient of lateral earth pressure at rest and is given by = 1 − , is the friction angle. In this study, of all soil samples are 31 and is 0.48. It should be mentioned that k is obtained by variable head experiments because the seepage pipe would be blocked by CC in constant head tests.

Soil Mechanics Test
The particle size distribution of materials was obtained by the Mastersizer 3000 laser particle size analyzer ( Figure 3). The d 50 of two kinds of soils are 0.163 mm (fine sand) and 0.009 mm (clay). The physical properties of seabeds are summarized in Table 1, CC was defined as the ratio of the weight of clay to the weight of sand-clay mixture, where ρ s is the density of soil, k is the permeability coefficient, ω is the rate of water content, e is the void ratio of the soil, ρ d is dry density of soil, ρ dmax is maximum dry density and ρ dmin is minimum dry density, D r is the relative density and is given by Equation (1). k 0 is the coefficient of lateral earth pressure at rest and is given by k 0 = 1 − sinϕ, ϕ is the friction angle. In this study, ϕ of all soil samples are 31 and k 0 is 0.48. It should be mentioned that k is obtained by variable head experiments because the seepage pipe would be blocked by CC in constant head tests.

Soil Mechanics Test
The particle size distribution of materials was obtained by the Mastersizer 3000 laser particle size analyzer ( Figure 3). The of two kinds of soils are 0.163 mm (fine sand) and 0.009 mm (clay). The physical properties of seabeds are summarized in Table 1, CC was defined as the ratio of the weight of clay to the weight of sand-clay mixture, where is the density of soil, k is the permeability coefficient, is the rate of water content, e is the void ratio of the soil, is dry density of soil, is maximum dry density and is minimum dry density, is the relative density and is given by Equation (1).
is the coefficient of lateral earth pressure at rest and is given by = 1 − , is the friction angle. In this study, of all soil samples are 31 and is 0.48. It should be mentioned that k is obtained by variable head experiments because the seepage pipe would be blocked by CC in constant head tests.   Note: "CC0.5" represents a mixture of 0.5% clay content and 99.5% fine sand, and so on.
Compression modulus, E s , was also measured by compression test. The compression coefficient, a v , is obtained by the height change, ∆h i , of the confined soil sample after loading different loads, p. a v is calculated by Equation (2), e i is void ratio of soil under loading and is given by Equation (3). E s is given by E s = (1 + e)/a v and results show in Table 2.

Test Condition
The water depth above the soil surface was kept as 0.45 m for all tests. The percentage of clay content in the experiments are 0, 0.5, 1, 2.5, 5, 7.7, 9.9 and 14.2. The wave height in these tests are all 0.14 m and the wave period are all 1.6 s.

Procedure
The procedures in these experiments were described below: (1) Place the well-mixed soil in the trench; (2) Fill up the flume with water gradually to a depth of 0.45 m; (3) Wait 20-22 hours to allow the soil to consolidate; (4) Calibrate water gauges and preheating pressure transmitters to record the pore pressure; (5) Switch on the waves 40 minutes; (6) Switch off the waves, and wait until the suspended sediment settles; (7) Empty the flume slowly; (8) Empty the sediment pit, and re mix soil sample if the next experiment is for the clay-sand soil with the same clay content. Prepare for the next test.

Pore Pressure Response of Seabed
The relationship between clay and seabed liquefaction had been investigated by Gratchev et al. [7], Liu and Jeng [21,22], Kirca et al. [20] and Chavez et al. [23] through cyclic load simulation experiments and wave flume experiments. It was found that susceptibility of seabed liquefaction was increased with the addition of CC and there is a critical value of CC which liquefaction would not occur when the CC exceeded the critical value. However, it was neglected what would happen with low CC. To get further insight on mechanism of clayey seabed liquefaction, the effects of low CC on seabed liquefaction was investigated in this study. Table 3 gives an overview of experiments. The parameters of excess pore water pressure (EPP) are listed in Table 3. D is the distance form seabed surface to the EPP sensor. P ACC is the maximum of accumulated EPP. P AMP is the amplitude of EPP. N W is the number of waves required for EPP reaching P ACC . σ 0 is the critical mean normal effective stress that used to discriminate seabed liquefaction. It is given by Equation (4) in which γ is the submerged specific weight of the soil and k 0 is the coefficient of lateral earth pressure. σ v,0 is the initial vertical effective stress that used to discriminate seabed liquefaction in engineering practice [27]. It is given by Equation (5). T d is the time of pressure dissipation, the time point when EPP starts to drop under wave action. T c is the time of soil consolidation, the time point when EPP drop to the lowest point, which represents the completion of the consolidation. T d is the time of pressure dissipation. T c is the time of soil consolidation. From the Table 3, it can be concluded that the addition of CC caused liquefaction of seabed which would not liquefy without it. Even if only 0.5% clay was added to the seabed, a large accumulation of EPP occurred under the wave action. The accumulated EPP dissipated when the clay content is low, 0.5%-1%. Figure 4 shows EPP time series at z = 0.23 m in sand-clay cases. With the introduction of waves, EPP begins to build up. Sumer et al. [6] described the mechanism of buildup of EPP with following words. The waves generate cyclic shear strains in the soil. If the soil is loose, the cyclic shear strains will gradually rearrange the soil grains at the expense of the pore volume of the soil. The latter effect will pressurize the water in the pores, and presumably this will lead to the buildup of EPP in the case of nearly undrained soil.  With the buildup of EPP, an upward-directed pressure gradient was generated. Because of the low permeability of the seabed, pore water cannot be squeezed from the seabed. When the seabed was liquefied and in a fluid state, the pores between particles gradually became larger, the drainage and wave-induced compaction occurred. Due to the limitation of experimental facilities, The measurements show that seabed is not liquefied in the sand-alone case, where P ACC is much smaller than σ 0 . On the other hand, in the well mixed soil bed containing certain amount of CC, it is found that the EPP accumulated obviously. As can be seen from Figure 4, the lower clay content, the faster EPP accumulated. When CC is small, 0.5-4.9, the amplitude of pore pressure is large and the change of EPP is irregular. When CC is large, 7.7-9.9, the amplitude of pore pressure is smaller and the EPP is regular. This is because the liquefaction of low CC cases is more intense, resulting in severe deformation of the seabed. The change of seabed elevation leads to the variation of pressure. In addition, the time required for EPP accumulation to reach its maximum increases with the increase of CC. This is due to the increase of CC changes the permeability and plasticity of soil. In summary, with the increase of CC, the accumulation of EPP slowed down. As a result, the value of final EPP accumulation become smaller and amplitude of EPP as well.
With the buildup of EPP, an upward-directed pressure gradient was generated. Because of the low permeability of the seabed, pore water cannot be squeezed from the seabed. When the seabed was liquefied and in a fluid state, the pores between particles gradually became larger, the drainage and wave-induced compaction occurred. Due to the limitation of experimental facilities, wave-making time is limited. Therefore, wave-induced compaction process occurred only in low CC cases, CC0.5 and CC1.0. The compaction phenomena for the cases with larger CC that need comparatively longer wave acting time was not observed. EPP oscillated sinusoidally because of the wave action. In order to obtain a clearer tendency of EPP, the curve of EPP was smoothed by averaging the value over a wave period. Figure 5 shows period-averaged EPP in low CC cases at depths z = 0.09, 0.16 and 0.23 m. Period-averaged EPP is defined by Equation (6). It is seen that the observed EPP response in Figure 6 is fully consistent with the conceptual description of the liquefaction and compaction process in sand-bed by Sumer et al. [28] and Miyamoto et al. [29]. The process of liquefaction and compaction is classified into four stages, Pressure builds up, Liquefaction, Dissipation starts and Consolidation completed in Figure 6. Dissipation is shown clearly in Figure 5. It can be concluded that, with the decrease of CC, duration of Liquefaction decreases, time point of dissipation starts advances, and duration of Dissipation decreases. liquefaction and compaction process in sand-bed by Sumer et al. [28] and Miyamoto et al. [29]. The process of liquefaction and compaction is classified into four stages, Pressure builds up, Liquefaction, Dissipation starts and Consolidation completed in Figure 6. Dissipation is shown clearly in Figure 5. It can be concluded that, with the decrease of CC, duration of Liquefaction decreases, time point of dissipation starts advances, and duration of Dissipation decreases. = 1 (6) Figure 5. Period-averaged EPP in low clay contents (CC) bed (the wave height is 0.14 m and the wave period is 1.6 s).       Figure 8 shows the number of waves, N w , for period-averaged EPP to reach P ACC plotted against the CC at different depth. It can be seen that N w increases with the CC when CC is more than 2.5 regardless of the depth, and N w increases with the depth regardless of the CC except CC14.2. It is due to the decrease of permeability of mixed soil bed with the increase of CC, and more waves are needed to transmit wave energy to deeper soil layer under the same wave condition. When CC is 14.2%, wave energy cannot be transmitted to z = 0.23 m, thus, N w at z = 0.16m is bigger than N w at z = 0.23 m.  , for period-averaged EPP to reach plotted against the CC at different depth. It can be seen that increases with the CC when CC is more than 2.5 regardless of the depth, and increases with the depth regardless of the CC except CC14.2. It is due to the decrease of permeability of mixed soil bed with the increase of CC, and more waves are needed to transmit wave energy to deeper soil layer under the same wave condition. When CC is 14.2%, wave energy cannot be transmitted to z=0.23 m, thus, at z=0.16m is bigger than at z=0.23 m.

Mechanism of Clayey Seabed Liquefaction
As discussed previously, sand-clay mixture with CC less than approximately 10% is susceptible to liquefaction, the liquefaction susceptibility decreasing with increasing CC. The reason for the influence of CC on liquefaction of clayey seabed is linked to cohesive strength of clayey soil, seabed permeability and seabed compressibility.
Up to now, it is hard to directly measure the cohesive strength between particles in mixed clayey bed. Kirca et al. [20] measured undrained shear strength of pure clay. The cohesive strength of pure clay is rather small, 1.75 kPa, and the value will be smaller for mixture with low CC. Therefore, cohesion is not expected to play a major role in wave-induced liquefaction of sand-clay mixture, unless the CC is much higher than the critical value for the sand grains to be completely encapsulated by the clay matrices.
Seabed permeability is an important factor affecting seabed liquefaction [30]. It was studied in the past in conjunction with laboratory simulations of seismic-induced liquefaction by Prakash [31]. His results show that the smaller the permeability, the more susceptible the soil is to liquefaction. Seabed cannot liquefy with large permeability because all EPP in the soil would dissipate as rapidly as they develop. Gratchev et al. [7] studied permeability through observing the microfabric of artificial mixtures by using SEM. Mixtures were divided into four microfabric according to CC, non-plastic, low plasticity clayey sand, medium plasticity clayey sand and high plasticity clayey sand. It pointed out that CC reduced the number of contact points between the sand particles and clay bridge in low plasticity clayey sand provided a structure that was easier to rearrange under cyclic shear strains. Pore of high plasticity clay sand were filled with clay, so the permeability was very low and pore water cannot flow. Kirca et al. [20] also explained the liquefaction of silt-clay mixture with the theory of Gratchev et al. [7]. Wiebicke et al. [32] studied granular fabric under shearing by x-ray computed tomography. However, aforementioned theory cannot explain the liquefaction of very low CC cases. The coefficient of permeability for sand and low CC mixture are almost same, as shown in Table 1. Figure 9 shows microfabric of seabed obtained by SEM (Hitachi SU8100). It can be seen that clay aggregations generally gathered at the sand particle contact points when CC is low (CC0.5, CC1.0, CC2.5 and CC4.9) , and clay produced a matrix that made microfabric of seabed more compact when CC is high (CC7.7, CC9.9 and CC 14.2). With the increase of CC, clay aggregations and matrix increase, and porosity between particles decreases. Gratchev et al. [7] suggested that the aggregation

Mechanism of Clayey Seabed Liquefaction
As discussed previously, sand-clay mixture with CC less than approximately 10% is susceptible to liquefaction, the liquefaction susceptibility decreasing with increasing CC. The reason for the influence of CC on liquefaction of clayey seabed is linked to cohesive strength of clayey soil, seabed permeability and seabed compressibility.
Up to now, it is hard to directly measure the cohesive strength between particles in mixed clayey bed. Kirca et al. [20] measured undrained shear strength of pure clay. The cohesive strength of pure clay is rather small, 1.75 kPa, and the value will be smaller for mixture with low CC. Therefore, cohesion is not expected to play a major role in wave-induced liquefaction of sand-clay mixture, unless the CC is much higher than the critical value for the sand grains to be completely encapsulated by the clay matrices.
Seabed permeability is an important factor affecting seabed liquefaction [30]. It was studied in the past in conjunction with laboratory simulations of seismic-induced liquefaction by Prakash [31]. His results show that the smaller the permeability, the more susceptible the soil is to liquefaction. Seabed cannot liquefy with large permeability because all EPP in the soil would dissipate as rapidly as they develop. Gratchev et al. [7] studied permeability through observing the microfabric of artificial mixtures by using SEM. Mixtures were divided into four microfabric according to CC, non-plastic, low plasticity clayey sand, medium plasticity clayey sand and high plasticity clayey sand. It pointed out that CC reduced the number of contact points between the sand particles and clay bridge in low plasticity clayey sand provided a structure that was easier to rearrange under cyclic shear strains. Pore of high plasticity clay sand were filled with clay, so the permeability was very low and pore water cannot flow. Kirca et al. [20] also explained the liquefaction of silt-clay mixture with the theory of Gratchev et al. [7]. Wiebicke et al. [32] studied granular fabric under shearing by x-ray computed tomography. However, aforementioned theory cannot explain the liquefaction of very low CC cases. The coefficient of permeability for sand and low CC mixture are almost same, as shown in Table 1. Figure 9 shows microfabric of seabed obtained by SEM (Hitachi SU8100). It can be seen that clay aggregations generally gathered at the sand particle contact points when CC is low (CC0.5, CC1.0, CC2.5 and CC4.9), and clay produced a matrix that made microfabric of seabed more compact when CC is high (CC7.7, CC9.9 and CC 14.2). With the increase of CC, clay aggregations and matrix increase, and porosity between particles decreases. Gratchev et al. [7] suggested that the aggregation of clay at contact points of sand produced significant pore spaces and the contact of sand grains were easily destroyed during cyclic loading. The destruction of grains connection was the trigger for seismic-induced liquefaction. Different from seismic load, wave load comes from the fluctuation of sea water above the seabed. Figure 10 shows the mechanism of liquefaction. Seabed is compressed and seepage force is downward when wave is at the crest. The pore water pressure increases greatly due to the compression of soil skeleton and seepage effect. Seabed would rebound and seepage force is upward when wave is at the trough. The pore water pressure dissipates through drainage and rebound of soil skeleton. Due to the decrease of permeability and plastic deformation of the seabed after compression of the soil bed, neither drainage nor rebound of soil skeleton can be balanced with the compression and seepage effect of the seabed at wave crest, so the pore water pressure accumulates. As described by Li [33], excess pore water pressure is related to the change of soil volume. When the compression coefficient of soil skeleton is 0, any seepage action will not produce excess pore water pressure. So, the compressibility of seabed has a great influence on liquefaction. Compression modulus of seabed, E s , was measured by compression test and shown in Table 2. It can be seen that with the increase of CC, E s decreases and compressibility of seabed increases. Figure 11 shows the compression progress of clayey seabed. It is obviously that there will be clay aggregations at the sand grains contact points (from [7] and Figure 9) when CC is low. The contact friction force between particles that maintained the shape of soil skeleton becomes smaller due to the aggregations, so grains are easy to slip and rotate when seabed subjected to wave load. In a word, clay content improves the compressibility of clay-sand seabed and thus increases the possibility of seabed liquefaction.
In summary, the cohesive strength has little influence on the liquefaction of clayey seabed, while the permeability and compressibility play important roles in the liquefaction. As shown in the experimental results, CC does not only affect the permeability of seabed, but also change the structure of seabed, making it easier to compress. Under the influence of permeability and compressibility, the pore water pressure develops easily when the wave load acts on the seabed. In this experiment, the seabed with permeability of 0.000363-0.000613 cm/s and compression modulus of 10.53-7.60 MPa will liquefy under laboratory conditions. of clay at contact points of sand produced significant pore spaces and the contact of sand grains were easily destroyed during cyclic loading. The destruction of grains connection was the trigger for seismic-induced liquefaction. Different from seismic load, wave load comes from the fluctuation of sea water above the seabed. Figure 10 shows the mechanism of liquefaction. Seabed is compressed and seepage force is downward when wave is at the crest. The pore water pressure increases greatly due to the compression of soil skeleton and seepage effect. Seabed would rebound and seepage force is upward when wave is at the trough. The pore water pressure dissipates through drainage and rebound of soil skeleton. Due to the decrease of permeability and plastic deformation of the seabed after compression of the soil bed, neither drainage nor rebound of soil skeleton can be balanced with the compression and seepage effect of the seabed at wave crest, so the pore water pressure accumulates. As described by Li [33], excess pore water pressure is related to the change of soil volume. When the compression coefficient of soil skeleton is 0, any seepage action will not produce excess pore water pressure. So, the compressibility of seabed has a great influence on liquefaction.
Compression modulus of seabed, , was measured by compression test and shown in Table 2. It can be seen that with the increase of CC, decreases and compressibility of seabed increases. Figure 11 shows the compression progress of clayey seabed. It is obviously that there will be clay aggregations at the sand grains contact points (from [7] and Figure 9) when CC is low. The contact friction force between particles that maintained the shape of soil skeleton becomes smaller due to the aggregations, so grains are easy to slip and rotate when seabed subjected to wave load. In a word, clay content improves the compressibility of clay-sand seabed and thus increases the possibility of seabed liquefaction.  In summary, the cohesive strength has little influence on the liquefaction of clayey seabed, while the permeability and compressibility play important roles in the liquefaction. As shown in the experimental results, CC does not only affect the permeability of seabed, but also change the structure of seabed, making it easier to compress. Under the influence of permeability and compressibility, the pore water pressure develops easily when the wave load acts on the seabed. In this experiment, the seabed with permeability of 0.000363-0.000613 cm/s and compression modulus of 10.53-7.60 MPa will liquefy under laboratory conditions.   In summary, the cohesive strength has little influence on the liquefaction of clayey seabed, while the permeability and compressibility play important roles in the liquefaction. As shown in the experimental results, CC does not only affect the permeability of seabed, but also change the structure of seabed, making it easier to compress. Under the influence of permeability and compressibility, the pore water pressure develops easily when the wave load acts on the seabed. In this experiment, the seabed with permeability of 0.000363-0.000613 cm/s and compression modulus of 10.53-7.60 MPa will liquefy under laboratory conditions.

Conflicts of Interest:
The authors declare no conflict of interest.