Experimental and Empirical Studies to Evaluate Liquefaction Resistance of Partially Saturated Sands

: Partially saturated sands with occluded gas bubbles may occur naturally in the ﬁeld or can be created artiﬁcially as a way of mitigating liquefaction effects. This study aimed to investigate the main parameters that inﬂuence the liquefaction resistance of partially saturated sands containing occluded gas bubbles. It also adopted a simpliﬁed approach and proposed empirical relationships that linked the liquefaction-resistance ratio to the bulk modulus, maximum volumetric strain, and excess pore-pressure ratio. Considering these objectives, a total of 60 stress-controlled dynamic triaxial tests were conducted. Nearly/fully saturated and partially saturated sand specimens prepared using sodium percarbonate were subjected to undrained cyclic loading with a loading frequency of 0.1 Hz. The loosely or densely packed sand specimens at different degrees of saturation (79–100%) were examined under different initial effective conﬁning stresses (25–100 kPa) and back pressures (100–400 kPa). The analysis of the test results indicated that the liquefaction resistance of the partially saturated sands constantly increased as the degree of saturation decreased, and this effect was more prominent under higher initial effective conﬁning stresses and lower back pressures. The adopted method satisfactorily captured this trend and showed reasonable matches between the empirical predictions and experimentally observed results from this study and previous research.


Introduction
Earthquake-induced liquefaction is an important phenomenon of natural hazards that causes significant deformation of soil deposits and failure of geotechnical structures.The catastrophic effects of liquefaction have been recurrently observed in many seismic events such as the 1964 Alaska (USA) and Niigata (Japan) earthquakes [1,2], the 1999 Kocaeli and Duzce (Turkey) earthquakes [3], the 2010 Haiti and Maule (Chile) earthquakes [4,5], the 2011 Tohoku (Japan) and Christchurch (New Zealand) earthquakes [6,7], and the 2016 Ecuador earthquake [8,9].
Liquefaction occurs in granular deposits when excess pore pressures develop under seismic loading and cause a large reduction in the effective stress and significant degradation of the soil strength and stiffness.In aftermath of the 1964 earthquakes, comprehensive research has been conducted to elucidate the important mechanisms and fundamental understanding of earthquake-induced liquefaction.These studies have mostly benefited from the analysis of case histories and/or laboratory tests on undisturbed or reconstituted samples and have produced several semi-empirical principles and frameworks that are widely used in current design practice for the evaluation of liquefaction susceptibility and its effects [10][11][12][13].
Natural soil deposits in seismically active zones may be found in a nearly/fully saturated, unsaturated, or partially saturated state depending on the groundwater level, degree of saturation (S r ), and distribution of gas bubbles in the soil [14,15].Granular soil deposits tend to contract when subjected to cyclic shear strains, which causes the rearrangement and densification of soil grains.For the nearly/fully saturated state, which generally occurs well below the groundwater level, interparticle pores are filled with practically incompressible pore fluid (e.g., water).During an earthquake loading in an undrained condition, the pore fluid cannot move or contract because the rate of loading and the bulk stiffness of the pore fluid are high, thereby causing the build-up of large excess pore pressures.For the unsaturated state, which usually occurs well above the groundwater level, interparticle pores are filled with a continuous air and water meniscus, and the effective stress and hence liquefaction resistance of the soil deposit increase due to the presence of a large matric suction.For the partially saturated state occurring below the groundwater level, pores are filled with water and gas that may exist in a continuous or an occluded bubble form depending on the level of saturation.The gas takes an occluded bubble form at high degrees of saturation (typically S r ≥ 80%), whereas the continuous gas phase is prevalent below this level [16].In the former case, the presence of entrapped gas bubbles in pore spaces of submerged soil significantly reduces the bulk modulus and increases the compressibility of pore fluid, namely, the gas-water mixture.The gas bubbles tend to reduce in volume and absorb the excess pore pressures during an earthquake loading, thereby enhancing the resistance of the soil to liquefaction [17].
Partially saturated soils can be formed by naturally occurring gas bubbles via biological activities [18,19] or via human activities.In the latter case, the degree of saturation of soil deposits is reduced by using various induced partial saturation (IPS) methods to enhance the liquefaction resistance.Air injection [20][21][22], electrolysis and drainage-recharge [23], use of denitrifying bacteria [24,25], and use of chemical composition [26,27] are some of the most widely used IPS techniques.Previous studies on IPS methods showed that gas bubbles created in partially saturated soils can stay in the pore spaces and that the partially saturated state can last for a sufficiently long time [23,28,29].
Previous element and physical modeling tests showed that several parameters such as the degree of saturation (S r ), initial effective confining stress (σ c ), and back pressure (u 0 ) affect the co-and post-liquefaction behavior of partially saturated soils.The excess pore pressure generation, ground surface, and structural settlement are significantly reduced as the degree of saturation of soils is lowered [30][31][32][33].Through element tests, the liquefaction resistance of partially saturated soils was shown to increase with a decreasing degree of saturation [34][35][36][37][38][39][40].The observations based on cyclic element tests also indicated that not only the degree of saturation is of importance, but the effects of initial effective confining or vertical stress and back pressure also must be accounted for [17,41].Partially saturated soils under a higher confining stress and lower back pressure appear to show greater resistance to liquefaction.
As mentioned previously, the presence of gas bubbles in the pore spaces plays an important role in soil liquefaction, and a true assessment of this aspect is important in geotechnical earthquake engineering practice.The current study discusses the most important parameters that govern the liquefaction behavior of partially saturated sands with occluded gas bubbles (e.g., S r ≥ 80%).Moreover, an empirical framework was adopted to obtain a quick and reasonable estimation of the liquefaction-resistance ratio in terms of the bulk modulus, maximum volumetric strain, and excess pore-pressure ratio.To achieve these goals, a series of 60 stress-controlled undrained dynamic triaxial tests were conducted on nearly/fully saturated and partially saturated sand specimens prepared with a chemical composition of sodium percarbonate.The test parameters (e.g., degree of saturation, initial effective confining stress, back pressure, and relative density) were changed.Subsequently, the predictions of the proposed empirical methodology were compared with the experimental data obtained from this study and collected from the published literature.This study is expected to provide significant insights that may enhance the use of this particular induced-partial-saturation method in design practice.

Previous Literature
Liquefaction is directly associated with the build-up of excess pore pressures under cyclic loading, and therefore its prediction is of engineering significance.For this purpose, different stress-based [42][43][44], and energy-based predictive models [45] have been proposed in the literature.These models are intended for nearly or fully saturated sands with or without fines, and their robustness relies heavily on the soil and test-dependent empirical coefficients.
The gas bubbles in pore spaces reduce the bulk modulus of the pore fluid (gas-water mixture) and therefore affect the generation of excess pore pressures in partially saturated sands.Previous studies [46,47] suggested that excess pore pressures during one loading cycle (∆u) can be correlated with the volumetric strain increment (∆ε vd ), the rebound modulus of the soil skeleton (E r ), the porosity of the soil deposit (n), and the bulk modulus of the pore fluid (K f ): This model requires the determination of several soil parameters and experimental constants that define the magnitude of volume change and elastic rebound of soils under drained conditions.
For partially saturated sands with occluded gas bubbles (S r ≥ 80%), the surface tension between the gas and water can be ignored, and the gas bubble pressure is considered to be almost equal to pore water pressure.The bulk modulus of the pore fluid (K f ) for this particular case can be estimated using the relationship given in Equation ( 2) proposed by Koning [48]: where K w and K g are the bulk modulus of water and gas, respectively.The occluded gas bubbles in pore spaces of partially saturated sand initially exist in equilibrium with the surrounding pore water.During cyclic loading in undrained conditions, the change in pore water (∆u w ) and gas pressure (∆u g ) will be the same as the excess pore pressure (∆u).In this condition, the volumetric strains (ε v ) suffered by partially saturated sands can be estimated using the theoretical formulation in Equation (3): where e is the void ratio of the soil and P 0 is the absolute pore water pressure.The sand suffers the largest (potential or maximum) volumetric strain (ε * v ) when the excess pore pressure reaches the initial effective confining stress (σ c ).
Okamura and Soga [17] showed that a unique relationship exists between the maximum volumetric strain (ε * v ) and liquefaction resistance ratio (LRR) that is defined as the liquefaction resistance of partially saturated soils normalized with that of fully saturated soils.Few strain-based or stress-based empirical models were proposed for excess pore pressure generation [26,49] and for settlement in partially saturated soils [50,51].

Adopted Empirical Methodology
In the field of geotechnical earthquake engineering, studies on partially saturated soils are markedly less common than on fully saturated soils, which is probably due to difficulties and problems encountered during laboratory research or field investigation.A review of the literature indicated that only a few models are available to predict the important aspects of the liquefaction mechanism and response of partially saturated sands with occluded gas bubbles, so this topic requires further attention.This study used a potentially simple empirical methodology to evaluate the liquefaction-resistance ratio (LRR) in terms of pore fluid bulk modulus (K f ), maximum volumetric strain (ε * v ), and excess pore-pressure ratio (r u ).This method was built on the basis of experimental observations, analysis of the database compiled from the literature, and theoretical justification.
The liquefaction-resistance ratio (LRR) can be estimated as a function of the maximum volumetric strain (ε * v ) or the bulk modulus of the pore fluid (K f ): As Equation (3) indicates, the maximum volumetric strain (ε * v ) is a function of the degree of saturation (S r ), initial effective confining stress (σ c ), back or pore fluid pressure (u 0 ), and void ratio (e) or relative density (D r ): The normalized excess pore-pressure ratio ( r u−part r u−sat ) can be inversely correlated to the liquefaction-resistance ratio (LRR): where r u−sat and r u−part are the excess pore-pressure ratios generated in fully and partially saturated sands, respectively.

Testing Device and Test Material
The dynamic triaxial experiments were conducted at Mus Alparslan University using a fully automated test system (Dynatriax EmS).The triaxial apparatus was manufactured by Wykeham Farrance-Controls Groups (Milan, Italy); its main components are shown in Figure 1. with occluded gas bubbles, so this topic requires further attention.This study used a potentially simple empirical methodology to evaluate the liquefaction-resistance ratio () in terms of pore fluid bulk modulus ( ), maximum volumetric strain ( * ), and excess pore-pressure ratio ( ).This method was built on the basis of experimental observations, analysis of the database compiled from the literature, and theoretical justification.
The liquefaction-resistance ratio () can be estimated as a function of the maximum volumetric strain ( * ) or the bulk modulus of the pore fluid ( ): As Equation (3) indicates, the maximum volumetric strain ( * ) is a function of the degree of saturation ( ), initial effective confining stress ( ), back or pore fluid pressure ( ), and void ratio () or relative density ( ): The normalized excess pore-pressure ratio ( ) can be inversely correlated to the liquefaction-resistance ratio (): where  and  are the excess pore-pressure ratios generated in fully and partially saturated sands, respectively.

Testing Device and Test Material
The dynamic triaxial experiments were conducted at Mus Alparslan University using a fully automated test system (Dynatriax EmS).The triaxial apparatus was manufactured by Wykeham Farrance-Controls Groups (Milan, Italy); its main components are shown in Figure 1.The electromechanical actuator of the Dynatriax EmS system has a 15 kN dynamic load capacity, and the system is capable of producing loading cycles with various wave The electromechanical actuator of the Dynatriax EmS system has a 15 kN dynamic load capacity, and the system is capable of producing loading cycles with various wave shapes (e.g., sinusoidal) with frequencies ranging from 0.01 to 10 Hz.The software and data-acquisition system of the equipment provides automatic and effective control of the saturation, consolidation, and cyclic-shearing stages of the triaxial tests.The basic instrumentation of the system consists of a submersible load cell (with a 0-25 kN measurement range and accuracy of ±0.05% FS); pressure transducers that measure the pore water, cell, and back pressure (0-1000 kPa measurement range and accuracy of ±0.25% FS); external displacement transducers (±25 mm measurement range and accuracy of 0.2% FS); and a 100 cc automatic volume-change device (precision of ±0.1 cc).The instruments were periodically calibrated to minimize the uncertainties of the measurements made during the test program.
The triaxial specimens were prepared using Sile sand (AFS 55/60).Figure 2 presents the particle size distribution (PSD) curve of this sand and its comparison with different types of sand that are widely used in geotechnical laboratories around the world.Based on the liquefaction boundaries of Tsuchida [52], it is a highly liquefiable soil.The index properties of Sile sand are given in Table 1.
Appl.Sci.2023, 13, x FOR PEER REVIEW 5 of 23 shapes (e.g., sinusoidal) with frequencies ranging from 0.01 to 10 Hz.The software and data-acquisition system of the equipment provides automatic and effective control of the saturation, consolidation, and cyclic-shearing stages of the triaxial tests.The basic instrumentation of the system consists of a submersible load cell (with a 0-25 kN measurement range and accuracy of ±0.05% FS); pressure transducers that measure the pore water, cell, and back pressure (0-1000 kPa measurement range and accuracy of ±0.25% FS); external displacement transducers (±25 mm measurement range and accuracy of 0.2% FS); and a 100 cc automatic volume-change device (precision of ±0.1 cc).The instruments were periodically calibrated to minimize the uncertainties of the measurements made during the test program.The triaxial specimens were prepared using Sile sand (AFS 55/60).Figure 2 presents the particle size distribution (PSD) curve of this sand and its comparison with different types of sand that are widely used in geotechnical laboratories around the world.Based on the liquefaction boundaries of Tsuchida [52], it is a highly liquefiable soil.The index properties of Sile sand are given in Table 1.

Specimen Preparation
In this study, nearly/fully saturated and partially saturated sand specimens with a height/diameter ratio of about 2 were prepared in loose or dense conditions.For the cylindrical triaxial specimens, a porous stone and filter paper was positioned in the bottom part of the cell (base pedestal), and a rubber latex membrane was installed using a pair of O-rings.A rigid split mold (140 mm in height and 70 mm in diameter) was used to support the specimens.Suction was applied with a vacuum pump to remove air between the membrane and mold and to hold them together.The reconstituted, uniform, dry sand specimens were prepared inside the mold.The air (dry)-pluviation method was used for the loose specimens, while the dry-tamping method was chosen for the densely packed sand specimens.A filter paper, porous stone, and top adaptor (cap) were placed over the specimens, and membrane ends were sealed with another pair of O-rings.Vacuum pressure (e.g., 20 kPa) was applied to the specimens so that they could stand up by themselves and remain undeformed following the removal of the split mold.The triaxial cell was carefully lowered over the specimens and filled with water.The specimens were kept under suction until an effective confining stress of 20 kPa was applied by increasing the cell pressure.Figure 3 shows examples of test images recorded during the specimen preparation.

Specimen Preparation
In this study, nearly/fully saturated and partially saturated sand specimens with a height/diameter ratio of about 2 were prepared in loose or dense conditions.For the cylindrical triaxial specimens, a porous stone and filter paper was positioned in the bottom part of the cell (base pedestal), and a rubber latex membrane was installed using a pair of O-rings.A rigid split mold (140 mm in height and 70 mm in diameter) was used to support the specimens.Suction was applied with a vacuum pump to remove air between the membrane and mold and to hold them together.The reconstituted, uniform, dry sand specimens were prepared inside the mold.The air (dry)-pluviation method was used for the loose specimens, while the dry-tamping method was chosen for the densely packed sand specimens.A filter paper, porous stone, and top adaptor (cap) were placed over the specimens, and membrane ends were sealed with another pair of O-rings.Vacuum pressure (e.g., 20 kPa) was applied to the specimens so that they could stand up by themselves and remain undeformed following the removal of the split mold.The triaxial cell was carefully lowered over the specimens and filled with water.The specimens were kept under suction until an effective confining stress of 20 kPa was applied by increasing the cell pressure.Figure 3 shows examples of test images recorded during the specimen preparation.For the nearly/fully saturated sand tests, the dry specimens set in the triaxial cell and subjected to effective confining stress of 20 kPa were flushed with carbon dioxide gas.This process was repeated with de-aired water.For the partially saturated sand tests, the dry specimens under an effective confining stress of 20 kPa were flushed with a known volume of sodium percarbonate-water solution.This solution was prepared by mixing finely powdered sodium percarbonate with de-aired water in specified ratios.The chemical compound of sodium percarbonate reacted with water (for 1 day) and generated oxygen bubbles inside the pore spaces of specimens, which produced a partially saturated condition.The degree of saturation ( ) of the specimens was adjusted by changing the initial concentration of the solution.Details of the specimen-preparation procedure in partially saturated sands were elaborated in Zeybek [57].For the nearly/fully saturated sand tests, the dry specimens set in the triaxial cell and subjected to effective confining stress of 20 kPa were flushed with carbon dioxide gas.This process was repeated with de-aired water.For the partially saturated sand tests, the dry specimens under an effective confining stress of 20 kPa were flushed with a known volume of sodium percarbonate-water solution.This solution was prepared by mixing finely powdered sodium percarbonate with de-aired water in specified ratios.The chemical compound of sodium percarbonate reacted with water (for 1 day) and generated oxygen bubbles inside the pore spaces of specimens, which produced a partially saturated condition.The degree of saturation (S r ) of the specimens was adjusted by changing the initial concentration of the solution.Details of the specimen-preparation procedure in partially saturated sands were elaborated in Zeybek [57].

Testing Stages and Test Program
The specimens were saturated using the back-pressure method in which the cell and back pressures were steadily increased with a 10 kPa difference (Figure 4a).The back pressures (u 0 ) applied to the nearly fully saturated sand specimens were high (300-400 kPa), while they were relatively low (100-200 kPa) for the partially saturated sand specimens.The final B values (Skempton's pore pressure coefficient) were recorded; this was followed by the isotropic consolidation of the specimens (Figure 4b).The effective confining stress (σ c ) was in the range of 25-100 kPa in the tests.The consolidated specimens were cyclically sheared with a loading frequency of 0.1 Hz (Figure 4c).The applied cyclic-stress ratios (CSR) ranged from 0.1 to 0.45 depending on the test requirements.The cyclic loading was continued until the excess pore-pressure ratios (r u ) reached 1. Finally, excess pore pressures were allowed to dissipate (Figure 4d).The experiments were performed following the principles in ASTM, D5311/D5311M-13 [58].

Testing Stages and Test Program
The specimens were saturated using the back-pressure method in which the cell and back pressures were steadily increased with a 10 kPa difference (Figure 4a).The back pressures ( ) applied to the nearly fully saturated sand specimens were high (300-400 kPa), while they were relatively low (100-200 kPa) for the partially saturated sand specimens.The final B values (Skempton's pore pressure coefficient) were recorded; this was followed by the isotropic consolidation of the specimens (Figure 4b).The effective confining stress ( ) was in the range of 25-100 kPa in the tests.The consolidated specimens were cyclically sheared with a loading frequency of 0.1 Hz (Figure 4c).The applied cyclic-stress ratios () ranged from 0.1 to 0.45 depending on the test requirements.The cyclic loading was continued until the excess pore-pressure ratios ( ) reached 1. Finally, excess pore pressures were allowed to dissipate (Figure 4d).The experiments were performed following the principles in ASTM, D5311/D5311M-13 [58].Tables 2 and 3 present the testing parameters and typical experimental results obtained in the nearly/fully saturated and partially saturated tests.A total of 60 stress-controlled dynamic triaxial tests were conducted in undrained conditions; 16 of these tests were on nearly or fully saturated sand specimens, which provided a benchmark for the assessment of the influential parameters in partially saturated sand behavior.The rest of the 44 tests were carried out on partially saturated sand specimens.The target relative densities ( ) were 40% and 85% for the loose and dense specimens, respectively; these were recalculated upon completion of the consolidation of the specimens.The degrees of saturation ( ) of the partially saturated specimens were determined using the soil phase relationships while considering the amount of pore fluid leaving and entering the specimens during the saturation and/or desaturation process.Following the back-and cellpressure ramp, the final degree of saturation was in the range of 99-100% and 79-96% for the nearly/fully saturated and partially saturated sand specimens, respectively.Tables 2 and 3 present the testing parameters and typical experimental results obtained in the nearly/fully saturated and partially saturated tests.A total of 60 stress-controlled dynamic triaxial tests were conducted in undrained conditions; 16 of these tests were on nearly or fully saturated sand specimens, which provided a benchmark for the assessment of the influential parameters in partially saturated sand behavior.The rest of the 44 tests were carried out on partially saturated sand specimens.The target relative densities (D r ) were 40% and 85% for the loose and dense specimens, respectively; these were recalculated upon completion of the consolidation of the specimens.The degrees of saturation (S r ) of the partially saturated specimens were determined using the soil phase relationships while considering the amount of pore fluid leaving and entering the specimens during the saturation and/or desaturation process.Following the back-and cell-pressure ramp, the final degree of saturation was in the range of 99-100% and 79-96% for the nearly/fully saturated and partially saturated sand specimens, respectively.

Relationship between Degree of Saturation and B Parameter
Zeybek [57] conducted a series of beaker tests on partially saturated sand specimens prepared with a sodium percarbonate-water solution and established a correlation between the weight of the sodium percarbonate normalized with that of the de-aired water (R s/w ) and the degree of saturation (S r ).The prediction of his empirical equation was compared with the experimental data accomplished in the current study.Figure 5 shows that the degree of saturation (S r ) was reduced as the sodium percarbonate concentration in the aqueous solution was increased, and this trend was captured well by the empirical relation proposed by Zeybek [57]. Figure 5 also contains a plot of the degree of saturation values against the B measurements that shows a drop in the degree of saturation from nearly 1 to 0.8 as the B value decreased from almost 1 to 0.08; the power function in Equation ( 7) appeared to provide the best fit line for these datasets.

Relationship between Degree of Saturation and B Parameter
Zeybek [57] conducted a series of beaker tests on partially saturated sand specimens prepared with a sodium percarbonate-water solution and established a correlation between the weight of the sodium percarbonate normalized with that of the de-aired water ( / ) and the degree of saturation ( ).The prediction of his empirical equation was compared with the experimental data accomplished in the current study.Figure 5 shows that the degree of saturation ( ) was reduced as the sodium percarbonate concentration in the aqueous solution was increased, and this trend was captured well by the empirical relation proposed by Zeybek [57]. Figure 5 also contains a plot of the degree of saturation values against the B measurements that shows a drop in the degree of saturation from nearly 1 to 0.8 as the B value decreased from almost 1 to 0.08; the power function in Equation (7) appeared to provide the best fit line for these datasets.

Excess Pore-Pressure Generation and Stress Degradation
Typical experimental results in terms of excess pore pressure (u e ) generation and effective stress paths (q − p ) are presented in Figure 6 to show the liquefaction behavior of the nearly/fully and partially saturated sand specimens.The plots on the left-hand side of the figure (Figure 6a) show the triaxial test data for loose, fully saturated sand (S r ≈ 100%) tested under a back pressure of u 0 = 300 kPa and an initial effective confining stress of σ c = 100 kPa.The results on the right-hand side (Figure 6b) correspond to test data for loose, partially saturated sand (S r ≈ 95%) examined under a back pressure of u 0 = 100 kPa and an initial effective confining stress of σ c = 100 kPa.In both cases, the specimens were subjected to a cyclic-stress ratio of CSR = 0.2.
the figure (Figure 6a) show the triaxial test data for loose, fully saturated sand ( ≈ 100%) tested under a back pressure of  = 300 kPa and an initial effective confining stress of  = 100 kPa.The results on the right-hand side (Figure 6b) correspond to test data for loose, partially saturated sand ( ≈ 95%) examined under a back pressure of  = 100 kPa and an initial effective confining stress of  = 100 kPa.In both cases, the specimens were subjected to a cyclic-stress ratio of  = 0.2.
When both sand specimens were cyclically sheared in undrained conditions, the pore pressures accumulated and become almost equal to the initial total stress.The mean effective stresses ( ) were reduced and approached zero as the excess pore pressures ( ) increased and converged up to the initial effective confining stresses, which corresponded to liquefaction ( = 1).Under similar loading conditions, the degradation of the mean effective stress was consistently slower, and the number of cycles required for liquefaction ( ) was consistently larger in the partially saturated sands than their fully saturated counterparts.

Parameters Affecting Liquefaction Behavior of Partially Saturated Sands
This section presents the liquefaction resistance curves defined by the relationship between the cyclic-stress ratio () and the number of cycles required for liquefaction ( ) in order to discuss the effects of the initial effective confining stress ( ), degree of saturation ( ), back pressure ( ), and relative density ( ) on the sand behavior when subjected to undrained cyclic loading.When both sand specimens were cyclically sheared in undrained conditions, the pore pressures accumulated and become almost equal to the initial total stress.The mean effective stresses (p ) were reduced and approached zero as the excess pore pressures (u e ) increased and converged up to the initial effective confining stresses, which corresponded to liquefaction (r u = 1).Under similar loading conditions, the degradation of the mean effective stress was consistently slower, and the number of cycles required for liquefaction (N liq ) was consistently larger in the partially saturated sands than their fully saturated counterparts.

Parameters Affecting Liquefaction Behavior of Partially Saturated Sands
This section presents the liquefaction resistance curves defined by the relationship between the cyclic-stress ratio (CSR) and the number of cycles required for liquefaction (N liq ) in order to discuss the effects of the initial effective confining stress (σ c ), degree of saturation (S r ), back pressure (u 0 ), and relative density (D r ) on the sand behavior when subjected to undrained cyclic loading.

Effective Confining Stress
The effect of the initial effective confining stress on the liquefaction behavior of the sands was examined using the results from the tests in which the sand specimens were isotropically consolidated under different stress levels (σ c = 25, 50, or 100 kPa) while the other testing parameters (e.g., S r , u 0 , and D r ) were analogous.
Figure 7 shows the liquefaction-resistance curves for the nearly/fully saturated sand specimens.It can be seen that the CSR at a given N liq decreased as the σ c increased from 50 to 100 kPa.This implied that under comparable loading and experimental conditions, the resistance of the nearly or fully saturated sands to liquefaction was reduced with an increase in the initial effective confining stress.This result agreed with many research findings in the field of soil liquefaction [59,60].

Effective Confining Stress
The effect of the initial effective confining stress on the liquefaction behavior of the sands was examined using the results from the tests in which the sand specimens were isotropically consolidated under different stress levels ( = 25, 50, or 100 kPa) while the other testing parameters (e.g.,  ,  , and  ) were analogous.
Figure 7 shows the liquefaction-resistance curves for the nearly/fully saturated sand specimens.It can be seen that the CSR at a given  decreased as the  increased from 50 to 100 kPa.This implied that under comparable loading and experimental conditions, the resistance of the nearly or fully saturated sands to liquefaction was reduced with an increase in the initial effective confining stress.This result agreed with many research findings in the field of soil liquefaction [59,60].
Figure 8 exhibits the cyclic liquefaction-resistance curves obtained from two sets of tests conducted at two different degrees of saturation (e.g.,  = 94-96% and  = 84-86%).It appeared that for a constant relative density, back pressure, and degree of saturation, the CSR at a specific  increased as the initial effective confining stress was increased from 25 kPa to 100 kPa.This indicated that the partially saturated sand specimens were more resistant to liquefaction at high effective confining stresses.Similar observations were made in previous studies [17,41].It appeared that for a constant relative density, back pressure, and degree of saturation, the CSR at a specific N liq increased as the initial effective confining stress was increased from 25 kPa to 100 kPa.This indicated that the partially saturated sand specimens were more resistant to liquefaction at high effective confining stresses.Similar observations were made in previous studies [17,41].

Degree of Saturation
Figure 9 depicts the variation in the CSR with N liq at different degrees of saturation in the range of 79-100% under a constant relative density, effective confining stress, and back pressure.The cyclic-stress ratios at the same number of cycles were found to increase as the nearly/fully saturated sand specimens were desaturated using the IPS method.Irrespective of the level of saturation, the liquefaction resistance of the partially saturated specimens was larger than that of the nearly/fully saturated specimens.Further improvement in the liquefaction resistance was obtained via a further reduction in the degree of saturation, which was in agreement with previous experimental observations [41,[61][62][63].

Degree of Saturation
Figure 9 depicts the variation in the  with  at different degrees of saturation in the range of 79-100% under a constant relative density, effective confining stress, and back pressure.The cyclic-stress ratios at the same number of cycles were found to increase as the nearly/fully saturated sand specimens were desaturated using the IPS method.Irrespective of the level of saturation, the liquefaction resistance of the partially saturated specimens was larger than that of the nearly/fully saturated specimens.Further improvement in the liquefaction resistance was obtained via a further reduction in the degree of saturation, which was in agreement with previous experimental observations [41,[61][62][63].

Degree of Saturation
Figure 9 depicts the variation in the  with  at different degrees of saturation in the range of 79-100% under a constant relative density, effective confining stress, and back pressure.The cyclic-stress ratios at the same number of cycles were found to increase as the nearly/fully saturated sand specimens were desaturated using the IPS method.Irrespective of the level of saturation, the liquefaction resistance of the partially saturated specimens was larger than that of the nearly/fully saturated specimens.Further improvement in the liquefaction resistance was obtained via a further reduction in the degree of saturation, which was in agreement with previous experimental observations [41,[61][62][63].

Back Pressure
The CSR-N liq curves recorded at varying back pressures are presented in Figure 10 to show the effect of the back pressure.The nearly/fully saturated sand specimens tested at back pressures of 300 and 400 kPa showed an almost identical trend, which suggested that the back pressure had an insignificant effect on the liquefaction resistance of such soils.The results of the partially saturated sand tests under back pressures of 100, 150, and 200 kPa showed that back pressure affected the liquefaction behavior of the partially saturated sand specimens, which were observed to exhibit a lower liquefaction resistance at larger back pressures.This again correlated well with the previous findings of researchers [17,41].
at back pressures of 300 and 400 kPa showed an almost identical trend, which suggested that the back pressure had an insignificant effect on the liquefaction resistance of such soils.The results of the partially saturated sand tests under back pressures of 100, 150, and 200 kPa showed that back pressure affected the liquefaction behavior of the partially saturated sand specimens, which were observed to exhibit a lower liquefaction resistance at larger back pressures.This again correlated well with the previous findings of researchers [17,41].

Relative Density
Figure 11 shows the test results obtained in the nearly/fully saturated and partially saturated sand tests conducted at two different relative densities and with the same effective confining stress.Whilst an increase in the relative density caused an increase in the liquefaction resistance of both the nearly/fully saturated and partially saturated sand specimens, this effect was much more pronounced in the former case.The increasing liquefaction resistance of the sand specimens with an increasing relative density can be ascribed to their contraction/dilation tendency.Loose sands have a higher tendency to contract, whereas dense sands show less contractive or dilative behavior.Under comparable test conditions, dense sands tend to generate lower excess pore pressures during undrained cyclic loading, and they exhibit greater liquefaction resistance than loose sands.

Relative Density
Figure 11 shows the test results obtained in the nearly/fully saturated and partially saturated sand tests conducted at two different relative densities and with the same effective confining stress.Whilst an increase in the relative density caused an increase in the liquefaction resistance of both the nearly/fully saturated and partially saturated sand specimens, this effect was much more pronounced in the former case.The increasing liquefaction resistance of the sand specimens with an increasing relative density can be ascribed to their contraction/dilation tendency.Loose sands have a higher tendency to contract, whereas dense sands show less contractive or dilative behavior.Under comparable test conditions, dense sands tend to generate lower excess pore pressures during undrained cyclic loading, and they exhibit greater liquefaction resistance than loose sands.

Liquefaction-Resistance Ratio
The aforementioned test results were presented herein in terms of the normalized liquefaction resistance-or simply the liquefaction resistance ratio ( )-which will form the basis of the empirical relations developed in the subsequent sections.The lique-

Liquefaction-Resistance Ratio
The aforementioned test results were presented herein in terms of the normalized liquefaction resistance-or simply the liquefaction resistance ratio (LRR)-which will form the basis of the empirical relations developed in the subsequent sections.The liquefactionresistance ratio was calculated as below: where (CSR N=20 ) sat and (CSR N=20 ) part are the cyclic-stress ratios that correspond to 20 cycles for the nearly/fully saturated and partially saturated sands, respectively.Figure 12 shows that for a given loading cycle (N = 20), the liquefaction-resistance ratio (LRR) showed a tendency to increase with a decreasing degree of saturation.Similarly, as the effective confining stress increased, the LRR markedly increased in the partially saturated sand specimens and slightly decreased in the nearly/fully saturated sand specimens.

Proposed Empirical Functions
This section introduces empirical correlations for the practical assessment of the normalized liquefaction resistance of partially saturated soils in relation to the bulk modulus, maximum (potential) volumetric strain, and excess pore-pressure ratio.The empirical formulas were produced using the curve-fitting toolbox in MATLAB, and their estimations were validated using experimental data from the literature.

Bulk Modulus
As stated earlier, for partially saturated sands with occluded bubbles, the bulk modulus of the pore fluid (gas-water mixture) is an essential parameter that controls the liq- The effective confining stress plays a significant role in the stress-strain-strength behavior of saturated sands.Depending on the effective confining stress value, saturated sands may exhibit strain-softening or strain-hardening behavior during shear, and the sand behavior may change from contractive to dilative or vice versa [64].Loose sand specimens at high effective confining stresses are expected to show a highly contractive behavior and suffer a larger build-up of excess pore pressures during undrained cyclic loading.The slight decrease observed in the liquefaction resistance of fully/nearly saturated sand specimens with increasing confining stress could be ascribed to the enhanced contractive tendency of the sand.On the other hand, the liquefaction behavior of the partially saturated sand specimens with occluded gas/air bubbles was controlled by the pore-fluid compressibility or potential volumetric strains.The occluded gas/air bubbles in the pore spaces played a role in absorbing the excess pore pressures and eventually increased the resistance of the sand to liquefaction.The compressibility of the pore fluid and hence the liquefaction resistance were enhanced by lowering the degree of saturation and/or increasing the effective confining stress.

Proposed Empirical Functions
This section introduces empirical correlations for the practical assessment of the normalized liquefaction resistance of partially saturated soils in relation to the bulk modulus, maximum (potential) volumetric strain, and excess pore-pressure ratio.The empirical formulas were produced using the curve-fitting toolbox in MATLAB, and their estimations were validated using experimental data from the literature.

Bulk Modulus
As stated earlier, for partially saturated sands with occluded bubbles, the bulk modulus of the pore fluid (gas-water mixture) is an essential parameter that controls the liquefaction behavior.The bulk modulus of the pore fluid (K f ) was calculated using Equation (2) for a range of experiments conducted at varying degrees of saturation, confining stress, and back pressures.The bulk modulus of the partially saturated sands (K f −part ) was normalized with that of their saturated counterparts (K f −sat ) and plotted against the degree of saturation (S r ) and liquefaction-resistance ratio (LRR) as shown in Figure 13.A trend line was fitted to the data points that produced the relationship below: Appl.Sci.2023, 13, x FOR PEER REVIEW 16 of 23 ratio, and the constant of 0.9983 corresponds to the effects of other factors.The negative exponent suggests the decay of the liquefaction-resistance ratio with the increasing normalized bulk modulus of the pore fluid.It can be observed in Figure 13 that the bulk modulus of the pore fluid significantly dropped as gas bubbles were generated in the nearly/fully saturated sand specimens.The reduction in the degree of saturation of the partially saturated sands was followed by a reduction in the bulk modulus of the pore fluid.This figure also reveals that the liquefaction-resistance ratio increased with a decreasing pore-fluid bulk modulus or an increasing-pore fluid compressibility, which was the reciprocal of the bulk modulus ( = ).

Maximum Volumetric Strains
The maximum volumetric strain ( * ) incorporates the effects of several parameters that govern the liquefaction response of partially saturated sands.Establishing a direct correlation between the liquefaction-resistance ratio ( ) and maximum volumetric strain is of great importance in engineering.The maximum volumetric strains were calculated using Equation (3) for different initial conditions that were considered during the experimental work in the present study.The calculated  * values are plotted against the liquefaction-resistance ratios () on the left-hand side of Figure 14; a power function describing the general trend was proposed as follows: This function has a dependent variable of the liquefaction-resistance ratio (LRR) and an independent variable of the normalized bulk modulus ( ) at its base.The exponent of −0.9 determines the growth/decay rate and shape of the function.The scaling factor of 0.0017 indicates the amplitude of the decrease/increase in the liquefaction resistance ratio, and the constant of 0.9983 corresponds to the effects of other factors.The negative exponent suggests the decay of the liquefaction-resistance ratio with the increasing normalized bulk modulus of the pore fluid.
It can be observed in Figure 13 that the bulk modulus of the pore fluid significantly dropped as gas bubbles were generated in the nearly/fully saturated sand specimens.The reduction in the degree of saturation of the partially saturated sands was followed by a reduction in the bulk modulus of the pore fluid.This figure also reveals that the liquefactionresistance ratio increased with a decreasing pore-fluid bulk modulus or an increasing-pore fluid compressibility, which was the reciprocal of the bulk modulus (C f = 1 K f ).

Maximum Volumetric Strains
The maximum volumetric strain (ε * v ) incorporates the effects of several parameters that govern the liquefaction response of partially saturated sands.Establishing a direct correlation between the liquefaction-resistance ratio (LRR) and maximum volumetric strain is of great importance in engineering.The maximum volumetric strains were calculated using Equation (3) for different initial conditions that were considered during the experimental work in the present study.The calculated ε * v values are plotted against the liquefactionresistance ratios (LRR) on the left-hand side of Figure 14; a power function describing the general trend was proposed as follows:  By considering a fully submerged sand element, the liquefaction-resistance ratios were estimated using Equation (10) and plotted against the degrees of saturation for five different field cases (Figure 15).In this analysis, back pressures corresponding to the pore water pressure in the field and the effective confining stresses used during the experimental studies were considered.For a buoyant unit weight of 10 kN/m , Case 1 represents a condition in which the water table was at the ground surface, whereas the groundwater table was located at 5-10 m above the ground surface for other cases.
The evidence in Figure 15 indicates that irrespective of the sand's depth and the position of the water table, the liquefaction-resistance ratios increased with decreasing degrees of saturation.In the right plot in the figure, it is apparent that the liquefaction-resistance ratio of sand under an equivalent depth (e.g., 10 m) slightly decreased as the water table rose, which caused an increase in the pore water pressure that acted on the sand element.This may be ascribed to the decrease in the pore-fluid compressibility with increasing back or pore water pressure [41].It is also evident in the left plot that under a constant initial pore water pressure, the liquefaction resistance of the sand increased quite  [17,35,36].
In this equation, the positive exponent (namely, 0.45) indicates the growth of the liquefaction-resistance ratio (LRR) with an increase in the normalized potential volumetric strains (ε * v ).The amplitude of the increase in the liquefaction-resistance ratio with increasing normalized potential volumetric strains is characterized by a scaling factor of 0.9.As the exponent is smaller than 1, the slope of the curve is expected to continuously decrease.
The right-hand side of Figure 14 compares the predictions of the power law relationship in Equation (10) with previously published experimental data obtained for different types of sand tested under a wide range of initial conditions (e.g., S r and D r ) and stress levels (e.g., σ c and u 0 ).It appeared that the normalized liquefaction resistance of partially saturated sands increased with the maximum volumetric strains and followed a non-linear trend.Most of the experimental data points were successfully captured by the proposed empirical power function.
By considering a fully submerged sand element, the liquefaction-resistance ratios were estimated using Equation (10) and plotted against the degrees of saturation for five different field cases (Figure 15).In this analysis, back pressures corresponding to the pore water pressure in the field and the effective confining stresses used during the experimental studies were considered.For a buoyant unit weight of 10 kN/m 3 , Case 1 represents a condition in which the water table was at the ground surface, whereas the groundwater table was located at 5-10 m above the ground surface for other cases.

Excess Pore-Pressure Ratios
The excess pore-pressure ratio ( ), which accommodates the effects of excess pore pressure and the initial effective confining stress, is an important design parameter.This study assumed that a reasonable assessment of excess pore-pressure ratios can be achieved by establishing an inverse correlation between the liquefaction-resistance ratio () and the normalized excess pore-pressure ratio ( ) as described in Equation (11).
A similar approach was undertaken in previous research [50,51].
A group of experimental data was gathered from the literature through a comprehensive review of the elemental (e.g., undrained cyclic triaxial or simple shear) tests conducted on different types of sand with various degrees of saturation, effective confining stresses, and back pressures.The left window in Figure 16 shows the relationship between the liquefaction-resistance ratio and degree of saturation; this was obtained via an analysis of the test data from the present study and previous research.The right window in Figure 16 shows the normalized excess pore-pressure ratios ( ) calculated using Equation (11) as a function of the degree of saturation.Figure 16 also compares the predictions of Equations ( 10) and (11) with the experimental data.It is clear in Figure 16 that the proposed method successfully captured the decreasing trend in normalized excess pore-pressure ratios with decreasing degrees of saturation.A relative density of 40%; a back pressure of 100 kPa; and effective confining stresses of 20 kPa (shallower depth), 50 kPa (middle depth), and 100 kPa (deeper depth) provided an upper, an average, and lower boundaries for the normalized excess pore-pressure ratios, respectively.There were some discrepancies between the experimental and empirical re- The evidence in Figure 15 indicates that irrespective of the sand's depth and the position of the water table, the liquefaction-resistance ratios increased with decreasing degrees of saturation.In the right plot in the figure, it is apparent that the liquefactionresistance ratio of sand under an equivalent depth (e.g., 10 m) slightly decreased as the water table rose, which caused an increase in the pore water pressure that acted on the sand element.This may be ascribed to the decrease in the pore-fluid compressibility with increasing back or pore water pressure [41].It is also evident in the left plot that under a constant initial pore water pressure, the liquefaction resistance of the sand increased quite significantly with an increasing depth or effective confining stress.These results may imply that in comparison with the initial pore water pressure, the degree of saturation and effective confining stress played a more important role and were the controlling factors in the liquefaction resistance of the partially saturated sands.

Excess Pore-Pressure Ratios
The excess pore-pressure ratio (r u ), which accommodates the effects of excess pore pressure and the initial effective confining stress, is an important design parameter.This study assumed that a reasonable assessment of excess pore-pressure ratios can be achieved by establishing an inverse correlation between the liquefaction-resistance ratio (LRR) and the normalized excess pore-pressure ratio ( r u−part r u−sat ) as described in Equation (11).A similar approach was undertaken in previous research [50,51].
A group of experimental data was gathered from the literature through a comprehensive review of the elemental (e.g., undrained cyclic triaxial or simple shear) tests conducted on different types of sand with various degrees of saturation, effective confining stresses, and back pressures.The left window in Figure 16 shows the relationship between the liquefaction-resistance ratio and degree of saturation; this was obtained via an analysis of the test data from the present study and previous research.The right window in Figure 16 shows the normalized excess pore-pressure ratios ( r u−part r u−sat ) calculated using Equation ( 11) as a function of the degree of saturation.Figure 16 also compares the predictions of Equations ( 10) and ( 11) with the experimental data.The predictions of the empirical method were also compared with the physical modeling (e.g., centrifuge and shaking table) test data, which provided further validation of the results with different test and loading conditions.The excess pore-pressure data gathered from centrifuge tests [65] and 1 g shaking table tests [66] were processed and presented in the form of a relationship between the normalized excess pore-pressure ratio and degree of saturation (Figure 17).It appeared that the proposed methodology successfully captured the decreasing trend in the experimentally determined excess pore-pressure ratios with decreasing degrees of saturation.There were some differences in the measured and predicted values that can be mainly ascribed to a limited number of test data and experimental uncertainties such as the non-uniform distribution of the air/gas bubbles observed in physical model tests [65].From an overall response perspective, this difference is expected to be acceptable in design practice when considering the complex nature of the problem under consideration.It is clear in Figure 16 that the proposed method successfully captured the decreasing trend in normalized excess pore-pressure ratios with decreasing degrees of saturation.A relative density of 40%; a back pressure of 100 kPa; and effective confining stresses of 20 kPa (shallower depth), 50 kPa (middle depth), and 100 kPa (deeper depth) provided an upper, an average, and lower boundaries for the normalized excess pore-pressure ratios, respectively.There were some discrepancies between the experimental and empirical results, particularly at an S r of 95%.From a design perspective, this variation was, however, understandable when considering the experimental uncertainties and complexities inherent in assessing the cyclic behavior of partially saturated sands.
The predictions of the empirical method were also compared with the physical modeling (e.g., centrifuge and shaking table) test data, which provided further validation of the results with different test and loading conditions.The excess pore-pressure data gathered from centrifuge tests [65] and 1 g shaking table tests [66] were processed and presented in the form of a relationship between the normalized excess pore-pressure ratio and degree of saturation (Figure 17).It appeared that the proposed methodology successfully captured the decreasing trend in the experimentally determined excess pore-pressure ratios with decreasing degrees of saturation.There were some differences in the measured and predicted values that can be mainly ascribed to a limited number of test data and experimental uncertainties such as the non-uniform distribution of the air/gas bubbles observed in physical model tests [65].From an overall response perspective, this difference is expected to be acceptable in design practice when considering the complex nature of the problem under consideration.
sure ratios with decreasing degrees of saturation.There were some differences in the measured and predicted values that can be mainly ascribed to a limited number of test data and experimental uncertainties such as the non-uniform distribution of the air/gas bubbles observed in physical model tests [65].From an overall response perspective, this difference is expected to be acceptable in design practice when considering the complex nature of the problem under consideration.[65] and 1 g shaking table [66] test data.
The above results may suggest that when following the proposed methodology, the liquefaction resistance of partially saturated sand can be estimated if the liquefaction resistance of the nearly/fully saturated soils is known; the same applies to excess porepressure ratios.

Conclusions and Recommendations
Studies on the cyclic behavior of induced partially saturated sands are limited.Current design practices are mostly based on the methods developed for nearly/fully saturated soils, but these are unable to consider the effects of the presence of gas bubbles in soils.This study was intended to provide experimental insights into the testing parameters that affect the liquefaction behavior of partially saturated sands with occluded gas bubbles.Moreover, this study proposed an alternative and potentially simple empirical methodology in which the liquefaction-resistance ratios were related in a unique way to the bulk modulus, maximum volumetric strains, and excess pore-pressure ratios.
In this regard, a series of undrained dynamic triaxial tests were conducted on nearly/fully saturated and partially saturated sand specimens.The loose and dense specimens were reconstituted using the air-pluviation and dry-tamping methods.The partially saturated sands were prepared by percolating a sodium percarbonate-water solution through the specimens.The effects of four different testing parameters (namely, the degree of saturation, effective confining stress, back pressure, and relative density) were discussed with the aid of liquefaction-resistance curves.The validity and applicability of the proposed empirical approach were evaluated using the test results from the current study and the existing literature.The following conclusions were drawn: (1) The liquefaction behavior of sand was observed to change as the triaxial specimens were turned from nearly/fully saturated to a partially saturated state.The liquefaction resistance of the partially saturated sand specimens appeared to increase as the degree of saturation was reduced.(2) The undrained cyclic behavior of the partially saturated sand was not only sensitive to the degree of saturation but also to the level of the initial effective confining stress and back pressure (or initial pore water pressure) acting on the sand.At the same degrees of saturation, the partially saturated sands showed a greater resistance to liquefaction under larger initial effective confining stresses and lower back pressures.
Amongst these parameters, the effect of the back pressure was less significant than those of the degrees of saturation and confining stress.(3) The positive effect of the confining stress and the adverse effect of the back pressure require consideration in design practice.During the practical application of IPS methods for the treatment of liquefiable deposits in the field, the shallower layers should be desaturated more than the deeper layers.(4) The proposed empirical approach was shown to enable a practical assessment of the normalized liquefaction resistance of partially saturated sands as a function of the pore-fluid bulk modulus, maximum volumetric strains, and normalized excess pore-pressure ratios.The model seemed to successfully capture the trends followed by the experimental data and produced physically sensible estimations for a complex topic of geotechnical earthquake engineering.
This study provided additional data and valuable insights that may enhance the scientific understanding of partially saturated soil behavior and the practical use of this particular IPS method in engineering projects.The experimental observations and empirical analysis were valid only for submerged, partially saturated clean sand specimens with occluded gas bubbles in which the effect of the matric suction was insignificant.Future studies should focus on sand deposits that contain some amount of silts and/or clays, and the robustness of the proposed methodology should be tested for different types of soils tested under a wider range of relative densities, degrees of saturation, and stress levels.

Figure 2 .
Figure 2. Particle size distribution curve for Sile sand and its comparison with different types of sand and liquefaction boundaries suggested by Tsuchida [52].

Figure 2 .
Figure 2. Particle size distribution curve for Sile sand and its comparison with different types of sand and liquefaction boundaries suggested by Tsuchida [52].

Figure 3 .
Figure 3. Examples of test images recorded during the specimen preparation.

Figure 3 .
Figure 3. Examples of test images recorded during the specimen preparation.

Figure 5 .
Figure 5. Relationship between the degree of saturation S r and Skempton's pore-pressure parameter B.

Figure 7 .
Figure 7. Effect of confining stress on liquefaction resistance of nearly/fully saturated sand specimens.

Figure 7 .
Figure 7. Effect of confining stress on liquefaction resistance of nearly/fully saturated sand specimens.

Figure 8
Figure8exhibits the cyclic liquefaction-resistance curves obtained from two sets of tests conducted at two different degrees of saturation (e.g., S r = 94-96% and S r = 84-86%).It appeared that for a constant relative density, back pressure, and degree of saturation, the CSR at a specific N liq increased as the initial effective confining stress was increased from 25 kPa to 100 kPa.This indicated that the partially saturated sand specimens were more resistant to liquefaction at high effective confining stresses.Similar observations were made in previous studies[17,41].

Figure 8 .
Figure 8.Effect of confining stress on liquefaction resistance of partially saturated sand specimens.

Figure 9 .
Figure 9.Effect of degree of saturation on liquefaction resistance of partially saturated specimens.

Figure 8 .
Figure 8.Effect of confining stress on liquefaction resistance of partially saturated sand specimens.

Figure 8 .
Figure 8.Effect of confining stress on liquefaction resistance of partially saturated sand specimens.

Figure 9 .
Figure 9.Effect of degree of saturation on liquefaction resistance of partially saturated specimens.Figure 9. Effect of degree of saturation on liquefaction resistance of partially saturated specimens.

Figure 9 .
Figure 9.Effect of degree of saturation on liquefaction resistance of partially saturated specimens.Figure 9. Effect of degree of saturation on liquefaction resistance of partially saturated specimens.

Figure 10 .
Figure 10.Effect of back pressure on liquefaction resistance of sand: (a) nearly/fully saturated sand specimens; (b) partially saturated sand specimens.

Figure 10 .
Figure 10.Effect of back pressure on liquefaction resistance of sand: (a) nearly/fully saturated sand specimens; (b) partially saturated sand specimens.

Figure 11 .
Figure 11.Effect of relative density on liquefaction resistance of sand: (a) nearly/fully saturated sand specimens; (b) partially saturated sand specimens.

Figure 11 .
Figure 11.Effect of relative density on liquefaction resistance of sand: (a) nearly/fully saturated sand specimens; (b) partially saturated sand specimens.

Figure 12 .
Figure 12.Variation in liquefaction resistance and liquefaction-resistance ratio with the degree of saturation and effective confining stress.

Figure 12 .
Figure 12.Variation in liquefaction resistance and liquefaction-resistance ratio with the degree of saturation and effective confining stress.

Figure 13 .
Figure 13.Variation in liquefaction-resistance ratio with normalized bulk modulus of pore fluid.

Figure 13 .
Figure 13.Variation in liquefaction-resistance ratio with normalized bulk modulus of pore fluid.

23 Figure 15 .
Figure 15.Variation in liquefaction-resistance ratio with the degree of saturation under different field conditions.

Figure 15 .
Figure 15.Variation in liquefaction-resistance ratio with the degree of saturation under different field conditions.

Table 1 .
The index properties of Sile sand.

Table 1 .
The index properties of Sile sand.

Table 2 .
Test parameters and typical results for nearly or fully saturated sand specimens.

Table 3 .
Test parameters and typical results for partially saturated sand specimens.