Variability of Bimodal Soil-Water Characteristic Curves under Different Confining Pressures

Soils with two subcurves of Soil-Water Characteristic Curve (SWCC) (dual porosity soils) might be found within various residual soils. Soils located in different depths have different confining pressure. Residual soils are found in the unsaturated zones due to the deep groundwater table. 'ere is a linear correlation between the hydraulic properties of the soil in the unsaturated area and that of its unsaturated properties. 'is study aims to examine the influence of the confining pressure towards the SWCC of dual porosity soil. 'e scope of this study involves measurements of the drying and wetting SWCC using Tempe cells, pressure plates, and an advanced triaxial apparatus. In this study, the mathematical equations were developed to explain the effect of confining pressure on SWCC. 'e experimental results indicated that the dual porosity soil exhibits bimodal characteristics for the drying curve of SWCC and it exhibits unimodal characteristics for the wetting curve of SWCC. As the confining pressure increases, the air entry values, the inflection points, and the standard deviation of drying SWCC increase. In addition, the hysteresis of SWCC is becoming smaller with the increasing confining pressure.


Introduction
Saturated soil consists solely of soil solids and water as the air spaces within have been filled with water. e principles of classical soil mechanics can be used to correlate the properties and behavior of saturated soil. Unsaturated soil has voids within that contain water partially comprising soil entities, water, and a combination of water and air in its fourth phase termed contractile skin [1,2]. As a result of the various phases, the characteristics of unsaturated soil deviate from those of typical soil mechanisms. e layers of soil on top of the groundwater table are examples of soil that is being formed normally by nature [3,4], while those in backfills and layered soil barrier systems are created artificially.
Since current geotechnical engineering operations mainly employ the usage of unsaturated soil over saturated soil, experimental studies on the parameters of unsaturated soil can reap many benefits. In unsaturated soil, the porewater pressure (u w ) is lower than the pore-air pressure (u a ); this accounts for its positive matric suction (u a −u w ) [1,5,6]. e relative humidity of the pore-water vapour found in soil is defined by the matric and osmotic suction which makes up the whole suction structure [7]. Since the changes in osmotic suction are minuscule in relation to deviations in matric suction, the full suction may consider matric suction as the sole value for ease of computation; the reverse is also possible [1,8,9].
Numerous residual and colluvial soils are gap-graded and demonstrate bimodal grain size distribution (GSD) [10][11][12][13][14]. Mercer et al. [15] recorded that gap-graded soil is mainly correlated with dual porosity structure. Dual porosity soil (soil exhibiting bimodal soil-water characteristic curve or SWCC) comprises substantially large amounts of coarse soil particles and great amounts of small soil particles.
is makes up the large pores and small pores, respectively [16]. e rationale behind the formation of the two summits within the pore-size distribution function was being researched upon by [17]. ey inferred that the occurrence of two different porous systems within the capillary domain created the bimodal pattern of the SWCC. e porous systems are matrix (intra-aggregate, textural, micro) pores inside soil aggregates and structural (interaggregate, macro) pores between the aggregates. e soil-water hydrostatics in the matrix domain is directly related to the arrangement of the soil skeleton. e shape, size, and stability of the aggregates are also factors affecting the linkage between structural pores. e air entry value of the matrix pores is the separating screen between matrix and structural pores. It also corresponds to the amount of residual water in the structural pores [18]. Typical structure and pore size distribution of a dual porosity soil (soil with bimodal SWCC) are presented in Figure 1. Feuerharmel et al. [19] recorded that SWCCs of most colluvium soils exhibit clear bimodal properties that correspond to the bimodal features of GSD. Soil with bimodal characteristics of SWCC has two air entry values, two inflection points, and two standard deviations due to the presence of large pores and small pores within the soil matrix. Figure 1 presents the typical shape of dual porosity soil (soil with bimodal SWCC) with its variables.
e matric suction has a function linked to the water content present in soil. Volumetric water content, gravimetric water content, or level of saturation can all be associated with the water content here [1]. As such, the water content can be drawn against matric suction values to obtain the SWCC. e SWCC can be used as a forecast for the absorption rate of unsaturated soils [20].
Typical models of SWCCs consist mainly of the unimodal form (one curve) or bimodal form (two subcurves). Correlations can also be drawn between the outlying form of the grain size distribution (GSD) curve and the SWCC. Historical researches pointed towards the indication that bimodal GSD soils are more probable to have bimodal pore size distributions [21], where the dual porosity is created by macropores and micropores. Saturated water content and dry density can act as indicators for the type of bimodal and unimodal characteristics that are demonstrated in the SWCC for soils with bimodal GSD [16,22,23]. e fine grains are unable to completely patch the empty voids between the larger grains; therefore, a bimodal SWCC would give rise to a bimodal pore size distribution [13]. e confining pressure applied on soil is one of the factors affecting its SWCC [24,25]. Ng and Pang [25] determined that the greater the net confining pressure applied on the soil, the greater the air entry value, and the slope of the SWCC is more gradual. Lee et al. [24] illustrate the test results showing that there is an inverse relationship between applied net confining pressure, at any given matric suction, and the saturated volumetric water content. However, it has a direct relationship with AEV in the soil.
Residual soils (dual-porosity soil) with two subcurves of SWCC can be found in many places in Singapore [26]. e soils located in different depths have different confining pressures. Residual soils are commonly located in the unsaturated zones due to the deep groundwater tables. Tropical regions like Singapore experience many rainfalls throughout the year [4]. e water infiltration into the soil layer results in pore-water pressure variations during rainy period. Hence, it will affect the behavior of dual porosity soil [27]. Other studies have also been conducted to investigate the SWCC of cohesive and noncohesive soils by Fattah et al. [28], Abd et al. [29], Fattah et al., [30], Abd et al. [31], and Fattah et al. [32].
Past studies only focused on the investigation of the variations of SWCC for unimodal soil due to different confining pressure.
ere are no studies on the effect of confining pressure on SWCC with bimodal characteristics or dual porosity soil. Examining the various possible outcomes on the SWCC of dual porosity soil with different confining pressure is the primary aim of this study. e scope of this study involves measurements of drying and wetting SWCC using Tempe cells, pressure plates, and an advanced triaxial apparatus. In addition, mathematical equations were developed to understand the effect of confining pressure on SWCC of dual porosity soils.

Applicable Theories
e GSD equation for the best fitting bimodal GSD used in this study was developed by [16] as shown in Soil suction, u a -u w (kPa) Figure 1: Soil with bimodal SWCC.

Applied and Environmental Soil Science
Here, we have the following: d � diameter of soil particle (mm) d r � parameter representing the maximum particle size of fine particles d m � parameter representing the geometric mean of soil particle diameter (mm) d min � parameter representing the minimum particle size that can be measured using hydrometer analysis (typical value � 0.0001 mm) d max � parameter representing the maximum diameter of soil particle (mm) s d � parameter representing the geometric standard deviation of GSD curve W � parameter representing the proportion of soil particle within coarse-or fine-grained particles d m � geometric mean of soil particle %d 0.001 � percentage of soil particle with diameter smaller than 0.001 mm Subscripts 1 and 2 represent subcurves 1 and 2, respectively.
For SWCC curves that display bimodal properties, equation (4) proposed by [16] was used: where we have the following: θ w � determined volumetric water content θ s � volumetric water content in saturated condition θ r � volumetric water content in residual condition ψ r � (u a −u w ) at residual condition (kPa) ψ m � (u a −u w ) associated with dominant pore size (kPa) ψ a � air-entry value of soil (kPa) s � SWCC geometric standard deviation erfc � the complimentary error function Subscripts 1 and 2 represent subcurves 1 and 2, respectively

Research Program
Laboratory tests were conducted on compacted soils in order to investigate the influence of confining pressure on bimodal SWCC. e compacted soils were produced by mixing coarse kaolin and Ottawa sand. e compacted soils were used in the research to avoid the heterogeneity of the soils for better analyses. ASTM graded sand, Ottawa sand, which is furnished by U.S. Silica Company, was selected to be mixed with the coarse kaolin which is produced by Kaolin Malaysia SDN BHD (Malaysia) to produce the sand-kaolin mixture. e index properties of the soil mixture of 50% kaolin and 50% Ottawa sand was obtained from the Atterberg limits test, grain size distribution test, hydrometer test, and specific gravity test. Index properties tests were carried out following the ASTM standards [33][34][35][36][37].
Following the steps recorded in [38], static compaction was performed. Even mixing of the soil was done manually. e water content started off at 12.2%, which correlates with the 95% maximum dry density at the dry optimum conditions of the compaction curve. ree separate soil samples with the same weight were obtained from the initial sample and put inside three circular moulds 50 mm wide. Uniform compaction using a static compaction machine then occurs to create folds of 10 mm with constant loading rate at 1 mm per minute. Lastly, physical data (i.e., weight and dimensions) for the 30 mm × 50 mm specimen were determined after it was extruded.
e experimental works to obtain drying and wetting SWCC were conducted using Tempe cell, pressure plate [39], and modified triaxial cell. e relevant procedures and equipments for each tests are explained in Subsections 3.1-3.3.

SWCC Testing Using Tempe Cell.
Tempe cell is mainly utilized at low pressures of one to hundred kilopascals, whereas the pressure plate was used at higher pressures of hundred to five hundred kilopascals [40]. Daily recording of the mass of the specimen was then to be taken. Following the formation of the soil specimen, with equal proportions of kaolin and sand, it was then tailored so that it goes precisely into the brass cylinder. Application of deoxygenated purified water by operating the vacuum desiccator onto the porous ceramic plate will make it completely saturated. e brass cylinder is then rested onto the plate. e volume of the cylindrical specimen can be obtained using the data of its height and diameter. Figure 2 depicts the setup of the Tempe cell. e main part of the Tempe cell and pressure plate is the porous ceramic plate. e procedures of SWCC testing were carried out following the study in [42]. Prior to linking up the Tempe cell to the pressure outlet of the manifold, saturation of the soil specimen was performed. e saturation process was ceased when the water content of the soil specimen reached its saturated value. After the saturation stage, the weight of the saturated specimen and the Tempe cell was recorded. e air pressure system was connected from the top of the Tempe cell. Air pressure was then applied to the specified value to create respective matric suction in the soil specimen. Air would not flow through the porous ceramic plate when higher air pressure was introduced into the cell, as long as the air pressure did not exceed the air-entry value of the ceramic plate and the ceramic plate was kept saturated. For this reason, the bottom of the Tempe cell was connected to a container filled with water to maintain the degree of saturation of the ceramic plate. Weighing of the specimen at regular intervals was necessary to obtain sufficient data for the plot of the water volume change in the specimen at various matric suctions. e ceramic disk and water compartment were flushed every time after the readings were taken to ensure that the ceramic disk and the water compartment are in saturated condition.
Monitor and increase the air pressure inside the cell to the required amount using pressure gauges. By increasing the pore-air pressure (u a ), it meant that the matric suction (u a −u w ) would increase because the pore-water pressure (u w ) remained at zero throughout the SWCC test. When measurements of the mass of the specimen need to be taken at the scheduled time intervals, the Tempe cell will be unplugged from the pressure outlet. After no further increase in mass is observed, the air pressure was turned up and this will cause the matric suction (u a −u w ) in the Tempe cell to be increased. Rerun these steps until the final matric suction (u a −u w ) is obtained. e drying SWCC measurement followed the matric suction incremental steps of 0.5, 1, 3, 5, 10, 20, 30, 40, 60, and 80 kPa. en, the specimen was moved into pressure plate. Once the drying process in pressure plate is completed, the wetting process was started. e wetting process was conducted by reducing the matric suction within the soil specimen. e wetting SWCC measurement at Tempe cell followed the matric suction decremental steps of 50, 20, 10, 5, and 1 kPa.

SWCC Testing Using Pressure Plate.
Operations at elevated pressure scale higher (u a −u w ), requiring a pressure plate to be employed. Pressure plate consists of a pressure chamber that contains a saturated ceramic plate. e pressure plate was connected to a burette for measurement of water volume change and maintenance of degree of saturation of the ceramic plate. ere must be a good contact between the soil specimen and the ceramic plate to ensure the water flow between the specimen and the plate is continuous. e highest air pressure applied on the pressure plates in this study is 400 kPa. Figure 3 illustrates the setup of a pressure plate which comprises a chamber with a pressure plate cell and a burette being held vertically upright using a stand. roughout this experiment, the soil specimen was taken out of the pressure plate at designated time intervals to take measurements of its mass. After no significant change in mass was being observed, the air pressure was raised to increase the matric suction (u a −u w ) in the pressure plate. ese steps were repeated until the final matric suction (u a −u w ). e drying SWCC measurement at pressure plate followed the matric suction incremental steps of 100, 200, and 400 kPa. Once the drying process in pressure plate is completed, the wetting process was started. e wetting SWCC measurement at pressure plate followed the matric suction decremental steps of 200 and 100 kPa.

SWCC Testing Using Modified Triaxial Cell.
Modified triaxial cell was utilized to perform SWCC tests with the application of confining pressures. A modified triaxial apparatus, as described by [1], was used for conducting the SWCC tests. A 5-bar high air-entry ceramic disk with a thickness of 6.35 mm and a diameter of 50 mm was sealed on the pedestal by applying slow settling epoxy glue along its circumference.
e specimen was placed directly on the saturated ceramic disk while filter paper, coarse corundum stone, and top cap were placed on top of the specimen. Rubber membrane was put on the specimen and secured with O-rings at the top cap as well as the pedestal. All specimens were saturated at the beginning of the test to have a uniform initial condition and to reduce the matric suction to a low value. Saturation was performed by applying cell pressure, σ3, and back pressure, uw, from digital pressure and volume controller (DPVC). A net confining pressure of 10 kPa was maintained to prevent significant swelling of the specimen until the pore-water pressure parameter, B, was larger than 0.95 as suggested by [43].
Isotropic consolidation of the specimen was applied directly after saturation. Consolidation ceases only after the equilibrium phase in the pore-water volume is achieved. Once consolidation was completed, drying process of the  SWCC was conducted by applying the designated matric suction using the axis translation technique [44]. Pore-water and cell pressures and volume changes were applied and controlled by DPVCs. Drying process of SWCC was achieved by following the designated incremental steps of matric suction while wetting process of SWCC was achieved by following the designated decremental steps of matric suction. After the matric suction and water volume change reached equilibrium, the matric suction was increased to next step by decreasing the pore-water pressure while maintaining constant confining and pore-air pressures. After the specimen reached the maximum matric suction 440 kPa in this research, wetting stage of the specimen was started by increasing the pore-water pressure while maintaining constant confining and pore-air pressures. Before incrementing or decrementing matric suction, the water compartment was flushed to remove diffused air from the water compartment.

Presentations and Results
e results of index property tests are shown in Table 1. e grain size distribution curve of the soil mixture is presented in Figure 4. Figure 5 illustrates the compaction curve derived, which is established by the dry density and water content determined for the respective soil samples. It also indicates the 95% maximum dry density and the maximum dry density. Figure 5 shows that the optimum dry density is 1.84 Mg/m 3 at 14.2% water content. e 95% maximum dry density at the dry optimum condition is 1.75 Mg/m 3 at 12.20% water content. Figure 6 illustrates the SWCC obtained from Tempe cell tests. e SWCC data from laboratory tests were best fitted using (4). en, the fitting parameter was compared with the variables of SWCC (i.e., AEV1, AEV2, inflection point 1, and inflection point 2) that were manually determined from the SWCC curve graph. Hence, the fitting parameters are close to the variables of SWCC. erefore, the fitting parameters from (4) were used to determine SWCC variables. e fitting parameters of SWCC in the Tempe cell tests are collated in Table 2. e air entry value 1 is 10 kPa, the air entry value 2 is 40 kPa, the inflection point 1 is 15 kPa, the inflection point 2 is 105.6 kPa, the standard of deviation 1 is 0.5, and the standard of deviation 2 is 1.202 for the drying curve. e inflection point is 5 kPa for the wetting curve. Figure 7 shows the SWCC from Triaxial Test with 25 kPa of net confining pressure. e determination of SWCC variables follow the same procedures as performed for SWCC data in Figure 6. e fitting parameters of SWCC from triaxial cell tests with 25 kPa net confining pressure are      Table 2. e air entry value 1 is 12 kPa, the air entry value 2 is 60 kPa, the inflection point 1 is 20 kPa, the inflection point 2 is 200 kPa, the standard of deviation 1 is 0.6, and the standard of deviation 2 is 1.35 for the drying curve. e inflection point is 15 kPa for the wetting curve. Figure 8 shows the SWCC for Triaxial Test with 50 kPa confining pressure. e determination of SWCC variables follows the same procedures as performed for SWCC data in Figure 6. e fitting parameters of SWCC from triaxial cell tests with 50 kPa net confining pressure are summarized in Table 2. e air entry value 1 is 20 kPa, the air entry value 2 is 90 kPa, the inflection point 1 is 35 kPa, the inflection point 2 is 300 kPa, the standard of deviation 1 is 0.65, and the standard of deviation 2 is 1.35 for the drying curve. e inflection point is 30 kPa for the wetting curve. Figure 9 shows the SWCC for Triaxial Test with 75 kPa net confining pressure. e determination of SWCC variables follows the same procedures as performed for SWCC data in Figure 6. e fitting parameters of SWCC from triaxial cell tests with 75 kPa net confining pressure are summarized in Table 2. e air entry value 1 is 30 kPa, the air entry value 2 is 200 kPa, the inflection point 1 is 50 kPa, the

Discussion
Based on the SWCC data presented in Section 5, the Air Entry Value increases with the increase of the net confining pressure as illustrated in Figures 10 and 11. SWCC tests with 0 kPa net confining pressure generated the lowest Air Entry Value, whereas SWCC tests with 75 kPa net confining pressure produced the highest Air Entry Value. is could happen since the higher the net confining pressure, the higher the density of the soil and hence the pores become smaller and it makes water more difficult to flow out of the pores. Figure 10 shows that AEV1 � 0.272 (confining pressure) + 7.8. Figure 10 shows that AEV2 � 2.04 (confining pressure) + 21.
Based on the SWCC data, the inflection point within subcurve 1 and subcurve 2 of the SWCC increases with the increase in the net confining pressure as illustrated in Figures 12 and 13. SWCC tests with 0 kPa net confining pressure produced the lowest inflection point, whereas SWCC tests with 75 kPa net confining pressure generated the highest inflection point. is might happen since the higher net confining pressure results in slower water flows out from the pores of the soil for suctions beyond Air Entry Values of the soil. Figure 12 indicates that inflection point 1 � 0.48 (confining pressure) + 12. Figure 13 shows that inflection point 2 � 6.93 (confining pressure) + 54.
Based on the SWCC data, the standard deviation of SWCC increases with the increase in the net confining pressure as illustrated in Figures 14 and 15. SWCC tests with 0 kPa net confining pressure generated the lowest standard deviation value, whereas SWCC tests with 75 kPa net confining pressure produced the highest standard deviation of SWCC. is might happen since the higher the standard deviation of SWCC, the gentler the slope of SWCC graphs and hence water flows slower from the pores of the soil for suctions beyond the air entry values. Figure 14 indicates that standard deviation 1 � 0.0032 (confining pressure) + 0.505. Figure 15 shows that standard deviation 2 � 0.0036 (confining pressure) + 1.216.
Based on the observation of all SWCC, the drying curve exhibits bimodal characteristics, but the wetting curve exhibits unimodal characteristics. e difference in the drying curve and wetting curve was because of the rearrangement of  Applied and Environmental Soil Science 7 soil particles and the entrapped air. e rearrangement of the soil particles might happen since the soil experienced a volume change (shrink) during the drying process. Hence, there was no more entrapped pores during wetting process and SWCC became unimodal. e drying process also resulted in the entrapped air within the pores of the soil. erefore, during the wetting phase, the soil specimen is unable to absorb back the same amount of water.
Based on the SWCC data, the water-entry value increases with the increase in the net confining pressure, as presented in Table 2. Wetting SWCC tests with 0 kPa net confining pressure generated the lowest water-entry value, whereas wetting SWCC tests with 75 kPa net confining pressure produced the highest water-entry value. is could happen since the higher the net confining pressure, the higher the density of soil. As a result, the soil pore is smaller and the water is more difficult to flow into the soil pore.

Conclusions
In conclusion, for the dual porosity soil, it exhibits a bimodal SWCC for the drying curve and it exhibits a unimodal SWCC for the wetting curve. e rearrangement of the soil particles and entrapped air brought about the difference between the drying and wetting curves.
As the net confining pressure increases, the air entry value increases, which is due to the increase in the density of the soil. When the net confining pressure rises, the inflection points on the drying curve increase, which shows that there is a decrease in the rate of water flowing out after the air entry values. As the net confining pressure increases, the standard deviation increases.
For the wetting curve of the SWCC, as the net confining pressure increases, the water-entry value increases. As the net confining pressure increases, the soil pore is smaller and the water is more difficult to flow into the soil pore.
is study is limited to sandy soil with no significant volume change. Further studies must be conducted to investigate the effect of confining pressure on fine-grained soil with high volume change.
Data Availability e experimental data used to support the findings of this study are included within the article.