Effect of Initial Conditions on the Pore Structure and Bimodal Soil–Water Characteristic Curve of Compacted Granite Residual Soil

: Granite residual soil typically forms complex pore structures and exhibits high water sensitivity due to physical and chemical weathering processes. Changes in initial compaction conditions significantly affect the mechanical and hydraulic properties of in situ granite residual soil sub-grades, with these variations fundamentally related to changes in pore structure and soil–water characteristics. This study investigates the pore structure and bimodal soil–water characteristic curve (SWCC) of a compacted granite residual soil through laboratory tests and mercury intrusion porosimetry tests. Nine initial conditions were selected based on potential in situ compaction conditions of subgrades, and their effects on the pore size distribution (PSD) and SWCC were thoroughly analyzed. The results show strong correlations between bimodal pore structure and SWCC. The size and volume of inter-aggregate pores exhibit noticeable changes with initial conditions, affecting SWCC within the low and middle suction range. Conversely, the intra-aggregate pores, which constitute over 60% of the pore structures, remain nearly intact across different initial conditions, resulting in similar SWCCs within the high suction range. As the compaction energy increases, the inter-aggregate pores are compressed and lead to a higher water retention capacity. In addition, as the compaction water content increases, the SWCC becomes less sensitive to compaction energy after the aggregates in the pore structure are fully saturated. Additionally, a three-dimensional bi-modal SWCC equation is proposed and validated using test data with an R 2 value above 0.98. These findings offer valuable insights for the design and quality control of granite residual soil subgrades.


Introduction
Granite is widely distributed in two-thirds of the Earth's crust.Due to the development of joints in granite, physical and chemical weathering can spread inward along the joints, forming thick weathering crusts [1,2].The upper soil layer of the weathering crust on granite is known as the granite residual soil, which partly preserves the structure of The in situ granite residual soil subgrades are typically in an unsaturated state, where the pore water pressure within the soil is lower than the pore air pressure.When the pores are connected to the atmosphere, the pore air pressure equals the atmospheric pressure, resulting in a negative pore water pressure.The pressure difference between pore air and pore water, known as matric suction, significantly influences the properties of the subgrade [5,6].From the perspective of unsaturated soil mechanics, the long-term performance of granite residual soil subgrade is directly related to its unsaturated shear strength and deformation characteristics, which are directly affected by the initial compaction condition, moisture variations, and fill height.Fundamentally, these properties are closely related to the pore structure and soil-water characteristics of the unsaturated soil [6,7].Variations in compaction water content and initial dry density (i.e., the initial condition) of compacted granite residual soil subgrades result in differences in initial pore structure and initial matric suction, and thus affect hydraulic and mechanical properties [8][9][10].The determination of SWCC equations can be further utilized to calculate the unsaturated effective stress, and quantitatively predict the mechanical behaviors of granite residual soil subgrades.Therefore, investigating the pore structure and soil-water characteristic curve (SWCC) of granite residual soil under different initial conditions can capture these effects and will be beneficial for designing and ensuring the long-term performance of granite residual soil subgrades [11][12][13].
In the past few decades, several researchers have investigated the effect of the initial conditions on the SWCC.Romero et al. [8] and Romero and Vaunat [14] found that SWCCs change obviously under different initial conditions.Ng and Pang [15,16] and Tarantino [17] investigated the effects of compaction water contents and dry densities on the SWCC and found that the influence of compaction water contents on the SWCC was more significant than the initial dry densities.Birle et al. [9] found that the effect of initial dry density was obvious when the SWCC was expressed by the degree of saturation-matric suction relationship.As the initial dry density increases, the air-entry value of the SWCC increases.In addition, the effects of the initial conditions on SWCCs can be different in different suction ranges.
With the development of microscopic structure detection techniques, it was found that the pore structure of unsaturated soils could be influenced by the initial conditions, and thus affected the SWCC [18].To obtain the pore structure information, two methods are commonly used: the scanning electronic microscope (SEM) and mercury intrusion porosimetry (MIP) [19,20].Through MIP tests, the pore structure characteristic can be quantitatively obtained.The pore size distribution of different soils can follow unimodal or bimodal patterns [21][22][23][24][25].For the soil with a bimodal pore structure, the pores are commonly separated into two categories: the inter-aggregate pores correspond to a larger dominant pore size, and the intra-aggregate pores correspond to a smaller dominant pore size [8,24,25].The inter-aggregate pores mainly change due to soil compaction, and the intra-aggregate pores change during the soil wetting-drying process [26].The SWCC of the soil with a bimodal pore size distribution (PSD) is shown as a bimodal "s" curve [27][28][29].Compared to a unimodal SWCC, the bimodal SWCC was more easily affected by the initial conditions and more complicated to determine [30][31][32].
Although the SWCCs for different initial conditions can be determined, respectively, its application was restricted by the lack of quantitative correlation among these SWCC equations.For this reason, researchers were devoted to propose a single equation to cover the change in SWCCs under different initial conditions.The void ratio (e) is one of the more commonly used indexes to quantify the effects of density on the SWCC, and thus the three-dimensional (3D) SWCC could be presented in the degree of saturation/matric suction/void ratio (Sr-s-e) space.Based on the Van Genuchten model [33] and Fredlund and Xing model [34], several 3D SWCC equations were proposed, as shown in Table 1 [35][36][37][38].In addition, Wijaya and Leong [39] established a 3D SWCC equation using a transition function and the parameters can be obtained from graphic determination method.
However, the equations are mainly established and validated for the unimodal SWCC, the determination of the 3D bimodal SWCC still requires further investigation.

Research Equation Parameters
Gallipoli et al. [35]  Sr: degree of saturation S: matric suction e: void ratio e0: initial void ratio w: gravimetric water content wsat: saturated gravimetric water content sr: residual matric suction n0: initial porosity Gs: soil gravity Ψ,φ,ng, mg, At, Bt, af, nf, mf = fitting parameters m: slope of segment i si: matric suction at the intersection between segment i and segment i -1 ki: curvature parameter between segment i and segment i -1

Wijaya and
Leong [39] ( ) ) In this paper, a high-liquid limit granite residual soil was sampled from Jiangxi, China.To address the significant variability in the compaction water content and degree of compaction of granite residual soil subgrades, nine different initial conditions were selected based on the potential in situ compaction conditions.The pore size distributions of the granite residual soil were investigated using MIP tests, and the filter paper method and vapor equilibrium method were employed to determine the bimodal SWCC.The corrections between the changes in PSD and SWCC were analyzed in detail.Based on the results, the 3D bimodal SWCC equation was proposed and validated using the test data.These findings will be valuable for the design and quality control of granite residual soil subgrades.

Materials
The material utilized in this study was a granite residual soil sampled from Jiangxi, China, with the grain-size distribution curve presented in Figure 1a.The physical properties of the soil were measured following Chinese standard JTG 3430-2020 (2020) [40] and are shown in Table 2.According to the Unified Soil Classification System, the tested soil was classified as a fine-grained high-liquid limit clay (CH).Also, the chemical composition of the tested granite residual soil was analyzed using X-ray powder diffraction (XRD) to obtain X-ray diffraction patterns [41][42][43].Subsequently, semi-quantitative analysis was conducted using JADE software, as shown in Figure 2. Results indicate that the soil mainly comprises quartz and kaolinite, along with gibbsite, lizardite, and goethite.where, w0, natural moisture content; wp, plastic limit; wL, liquid limit; Ip, plasticity index; Gs, specific gravity.To describe the change in the initial conditions of the compacted granite residual soil, nine different initial conditions were selected based on the proctor curves under different compaction energies, as shown in Figure 1b.The dried soil was moistened to three compaction water contents (wc) around the optimum moisture content (OMC), namely OMC − 5%, OMC, and OMC + 5%.After that, a different compaction energy was selected according to JTG 3430-2020 (2020) [40].The soil specimen of the proctor test was compacted in five layers.For each layer, a compaction hammer with the compaction energy of 2677.2 kJ/m 3 was applied for 98, 50, and 30 blows, respectively, and the proctor types were named as modified proctor (MP), standard proctor (SP), and reduced proctor (RP).To prepare test specimens, the air-dried tested soil was passed through 2 mm sieves, mixed with water to reach the compaction water content, and then statically compacted to the target dry density.The compaction water content, initial dry density, and initial void ratio of all test conditions are shown in Table 3.

Test Method
To obtain the SWCC of tested soil, a wide range suction measurement is needed.Therefore, the SWCCs for the low-medium and high suction ranges were measured using the filter paper method and the vapor equilibrium method, respectively.The two tests were conducted under controlled temperature conditions at 20 °C.Additionally, the MIP test was conducted to investigate the pore structure of the soil under different initial conditions.
(1) Filter paper method In this paper, the Whatman No. 42 filter paper with a 4.5 cm diameter was utilized.To ensure consistency in the initial conditions, the tested granite residual soil was compacted under the selected initial condition in a cutting ring measuring 6.18 cm in diameter and 2 cm in height, and was vapor-moistened or air-dried to eight different water contents around OMC, specifically OMC − 11.5%, −9.5%, −7.5%, −5%, −2.5%, +0%, +2.5%, +5%.Subsequently, a filter paper was placed between two soil specimens, and the specimens were placed inside a sealed box for 14 days to attain suction equilibrium, as depicted in Figure 3a.After that, the filter paper was weighed using an analytical balance with an accuracy of 0.0001 g.The calibration curve between the water content of the filter paper (wfp) and the matric suction (s) of soils is represented by Equation (1), following ASTM D5298-16 testing protocols [44] The saturated water content was determined individually using the vacuum saturation chamber (Figure 3b) to acquire the SWCC data point corresponding to zero suction.(2) Vapor equilibrium method In this paper, supersaturated salt solutions were utilized to provide different relative humidity levels, and different suction values could be obtained via the vapor equilibrium process [45,46].Five salt solutions were selected and summarized in Table 4.The test procedure is shown in Figure 4.The tested soil was compacted into a 1 cm height cylinder and was placed into a glass beaker.To create a sealed environment, a desiccator fitted with a porous disk was employed, with the specimen placed in a glass beaker suspended above the disk.The specimens were placed in the desiccators for a month to achieve vapor equilibrium.Subsequently, the change in the water content of soil specimens was measured to obtain the SWCC points.(3) Mercury intrusion porosimetry test For investigating the pore structure of the granite residual soil under various initial conditions, the widely employed mercury intrusion porosimetry (MIP) test was conducted [24,25,40].The Micromeritics AutoPore IV 9510 device(Micromeritics, Atlanta, GA, USA) with a maximum intrusion pressure of 414 MPa was used in this study.As depicted in Figure 5, the soil specimens were compacted and subsequently cut into 1 cm 3 cubes to match the dilatometer dimensions.Then, the cubes underwent freeze-drying at −50 °C for 24 h to exclude moisture without altering the pore structure [24].The dried cube was placed into a sealed dilatometer and was surrounded by mercury at a low pressure value.Then, the mercury pressure was increased step by step until it reached 414 MPa.The pore diameters (d) were recorded continuously as the mercury pressure increased, thus microstructure indexes of the tested soil, such as the cumulative intruded volume (ΣiΔVi), the cumulative intruded void ratio (ΣiΔVi/Vs), and the log-differential pore volume (−ΔVi/Δlog(di), could be obtained.The relationship between log-differential pore volume and pore diameter was named pore size distribution following related research [24,25].

Pore Size Distribution Curves
The pore size distributions of the tested soil under the nine initial conditions are depicted in Figure 6.For all the initial conditions, the soil specimens present bimodal pore size distributions.The intra-aggregate pores occupy a major part of pore structures and are nearly intact among different initial conditions.In contrast, the peak PSD value and the corresponding peak pore size of inter-aggregate pores change obviously, which is consistent with the results in related studies [47][48][49].Figure 6 shows that, for the same proctor type, the peak PSD value of inter-aggregate pores decreases as the compaction water content increases, for most cases.Compared with the specimen compacted at the dry side (i.e., wc lower than OMC) and wet side (i.e., wc higher than OMC), the specimens compacted at the optimum moisture content have the smallest peak inter-aggregate pore size and the narrowest inter-aggregate pore size distribution range.As the compaction energy decreases, the PSD curves are more easily affected by the compaction water content.
The change in inter-aggregate pores can be explained by differences in the initial void ratio.With higher compaction energy or less deviation between the compaction water content and OMC, the initial void ratio of soil specimens decreases, and thus the interaggregate pores were compressed due to the compaction effects [26].
Figure 7 shows that the change in PSD is closely correlated with the compaction water content.For specimens compacted at the dry side, an increased compaction energy reduces the peak pore size and pore size distribution range.For specimens compacted at OMC, the peak PSD values change slightly, and the peak pore sizes remain unchanged.For specimens compacted at the wet side, the peak pore size decreases as the compaction energy increases, while the peak PSD values remain unchanged.Compared with specimens compacted at the dry side, the specimens compacted at the wet side are less affected by the compaction energy.

Pore Size Distribution Indexes
To quantitatively analyze the pore size distribution characteristics, the bimodal pore structure is divided into two categories, as depicted in Figure 8.The diameter corresponding to the point of zero slope on the PSD curve is chosen to separate the two pore categories, and, thus, the total void ratio measured by MIP tests (et), inter-aggregate void ratio (einter), intra-aggregate void ratio (eintra), the proportion of intra-aggregate pores (pintra) [48], and the water ratio (ew) [47] can be calculated, as expressed in Equation (2) through Equation (6).Inter-aggregate pores where Vinter is the volume of inter-aggregate pores, Vintra is the volume of intra-aggregate pores, Vs is the volume of soil particles, and Vw is the volume of water.
The variation of the pore size distribution among the different initial conditions can be quantitatively captured by the pintra, which can capture the proportion of nearly intact intra-aggregate pores and changeable inter-aggregate pores.To correlate PSD with the initial conditions, the pintra is compared with the compaction water contents, dry density, and degree of saturation, respectively.Figure 9a shows that the pintra values for the specimens compacted at the dry side are separated from each other, while for the OMC and OMC + 5% conditions, the pintra values are nearly independent of proctor types.Figure 9b shows that an increased compaction energy reduces the variation in pintra.In addition, soil specimens with similar dry densities have different pintra values, which verify the results in related research that the soil compacted at the dry side or wet side show different pore structures [50].Figure 9c shows that the relationship between the degree of saturation and pintra can be described by the single-valued function, and the exponential fitting curve is depicted in Figure 9d.Thus, the pore size distribution parameter and initial condition parameter can be quantitatively correlated.To further investigate the correlation between the degree of saturation and pore structure, the filling degree of adsorbed water in intra-aggregate pores and the filling degree of capillary water in inter-aggregate pores can be quantitatively represented by the adsorbed saturation (Sr, a) and capillary saturation (Sr, c), respectively [48,51].The two indexes are calculated by Equations ( 7) and ( 8) and are compared in Figure 9.
It can be seen from Figure 10 that the intra-aggregate pores are not saturated under OMC − 5%, while for the OMC and OMC + 5% conditions, the intra-aggregate pores are fully saturated, and the capillary water fills in the inter-aggregate pores at different degrees.As the compaction water content increases from OMC − 5% to OMC, the filling of adsorbed water affects the volume and size of the aggregates significantly, thus affecting the PSD curve and corresponding pintra significantly.Once filled into the pores, the adsorbed water forms strong contact with the soil particles and is hard to dewater by the compaction effect.Thus, the pore structure becomes stable after the Sr, a reach 1.0.As the compaction water content increases from OMC to OMC + 5%, the pore structure is also affected by the Sr, c.Since the inter-aggregate pores occupy a relatively small proportion in the pore structure of the test soil, the PSD and corresponding pintra are less sensitive to the change in compaction water content and compaction effort, compared with the dry side.

Determination of the SWCC
The SWCC data obtained from the filter paper method, vapor equilibrium method, and saturated water content tests were summarized in Figure 11.To supplement the data of the low-suction section, the PSD curve is used to verify the air-entry value and the first inflection point of the SWCC [24,25].The measured SWCCs for all nine initial conditions show a typical bimodal pattern.Thus, the bimodal SWCCs are determined following the method proposed by Li et al. [52], as expressed in Equation ( 9).The parameters in Equation (9) are obtained graphically, as conceptionally depicted in Figure 12.In this equation, Sr, sat represents the degree of saturation under the fully saturated condition, Sr, R represents the residual degree of saturation, sa1 and sa2 are the air-entry values of inter-aggregate and intra-aggregate pores, and sR1 and sR2 represent the residual suction of inter-aggregate and intra-aggregate pores, respectively.

Effects of the Initial Conditions on the SWCC
The variations of bimodal SWCC among the different initial conditions are depicted in Figures 13 and 14. Figure 13 shows that, for the same proctor type, the SWCC of soil specimens compacted at OMC is higher than the specimens compacted at the wet side and Drainage of water in intra-aggregate pores dry side.The higher SWCC suggests a better water retention capacity, which means the water in the pores is harder to be drained under the same matric suction.Also, specimens compacted at the wet side generally show a better water retention capacity than those at the dry side.This phenomenon can be explained from the perspective of effective stress.The soil compacted at the wet side has a higher degree of saturation, leading to increased effective stress, hence improving the water retention capacity.Furthermore, as the compaction energy increases, the difference in SWCC decreases, which indicates that a better compaction condition can reduce the variation of water retention capacity with the moisture.Figure 14 shows that, for the same compaction water content type (i.e., OMC − 5%), an increased compaction energy improves the water retention capacity of compacted soil.For specimens compacted at the dry side, an increased compaction energy improves the water retention capacity significantly.In contrast, for the soil specimens compacted at OMC and OMC + 5%, the SWCC is less affected by the compaction energy.
To quantitatively compare the SWCC equations, the fitting parameters of the SWCCs are summarized in Table 5.Among all nine initial conditions, the air-entry values of interaggregate pores sa1 are less than 20 kPa, and the residual suction values of intra-aggregate pores sR2 are greater than 13,000 kPa, which suggests a wide effective suction range.Under the same compaction energy, the fitting parameters of the soil specimens compacted at OMC are generally higher than the corresponding values under OMC − 5% and OMC + 5%.This indicates that the granite residual soil compacted under the optimum moisture content has the strongest water retention capacity within the full suction range.Additionally, for the soil specimens compacted at OMC − 5% and OMC + 5% under a certain compaction energy, the relation of the first air-entry value sa1 and the first residual suction value sR1 are affected by the compaction energy, while their second air-entry value sa2 and the second residual suction sR2 are similar to each other.
Under the same compaction water content type, when the compaction energy changes, the first and second air-entry values sa1 and sa2 change obviously with the compaction energy.The average values under the modified proctor are 1.80 and 1.54 times higher compared with the specimens under the reduced proctor, respectively.In addition, the first and second residual suction sR1 and sR2 change significantly for the specimens compacted at OMC − 5%, and only change slightly for the specimens compacted at OMC and its wet side.

Correlation between Pore Structure and SWCC
Among the nine initial conditions, the difference in SWCCs of the tested compacted granite residual soil is obvious within the suction range of 10-1000 kPa, changes less within the suction range of 1000-10,000 kPa, and tends to be similar in the high suction range above 10,000 kPa.This phenomenon is consistent with the change of PSD curve among the different initial conditions.The size and volume of inter-aggregate pores change obviously with the initial conditions, which correspond to the change of the SWCC within the low-and middle-suction ranges.In contrast, the intra-aggregate pores are almost not affected by the initial conditions, thus the difference in SWCCs among the different initial conditions within the high-suction range is unobvious.
As the compaction energy increases, the size and volume of inter-aggregate pores decreases significantly (Figure 7).According to the Young-Laplace equation expressed by Equation (10), the matric suction is inversely proportional to the pore diameter.Thus, as the inter-aggregate pores are compressed, the water retention capacity of the tested soil is improved.

cos
where s is the matric suction; d is the pore diameter; Tw is the surface tension of water; and θw is the contact angle between the water and soil interface.
As the compaction water content increases from OMC − 5% to OMC, the aggregates in the pore structure are not fully filled by the adsorbed water, and the pintra values change obviously.This suggests an unstable pore structure, and, thus, the SWCC is affected by the compaction energy significantly.In contrast, as the compaction water content increases to OMC or above, the Sr, a reach 1.0 and the main proportion of the pore structure tends to be stable, thus the SWCC is less sensitive to the change in the compaction energy.

Determination of the 3D Bimodal SWCC Equation
The analysis above suggests that the pore structure and the SWCC are closely correlated with the initial conditions.To facilitate the application of the SWCC, a single equation can be established to describe the water retention capacity of the granite residual soil under the various initial conditions by extending the equation from 2D to 3D.Considering the formation of existing 3D SWCC functions [35,36,39], the initial void ratio is selected as the third variable to characterize the influence of the initial conditions on the SWCC. Figure 15 shows the correlation of the initial void ratio e0 and all the five SWCC fitting parameters.It can be seen that all relationships can be described using linear fitting, and the equations are also listed in Figure 15.The 3D SWCC equation and all the test data are depicted in Figure 16.It can be seen that the proposed 3D SWCC equation can capture the main characteristics of the bimodal SWCC and its variations among the different initial conditions.Figure 17 shows that the 3D bimodal SWCC equation fits the measured data well for the entire suction range with a high R2 value, which proves the effectiveness of the method used in this paper.
As the initial void ratio increases, all the fitting parameters of Equation ( 9) decrease, and the SWCCs tend to be lower in the degree of saturation/matric suction/initial void ratio (Sr-s-e0) space.The range of matric suction and degree of saturation covered by the first "s" shape increases with the increase in the initial void ratio.The starting suction value of the second "s" shape increases with the increase in the initial void ratio, accompanied by a gradual decrease in the slope of the descending section.These findings suggest that an increase in the density of granite residual soil results in a greater compression of inter-aggregate pores, leading to an overall enhancement of water retention capacity, especially within the middle-suction range.
Given the linear relationship between each fitting parameter and the initial void ratio, a recommended approach for utilizing the proposed 3D bimodal SWCC equation is presented.First, compute the initial void ratios based on physical parameters obtained from soil proctor curves.Subsequently, select several initial conditions with significantly different initial void ratios for SWCC testing, and then determine the proposed 3D bimodal SWCC equation.Following this, the bimodal SWCC for the untested initial conditions can be calculated.This approach can significantly reduce the number and duration of SWCC tests, facilitating the application of bimodal SWCCs in the design and research of unsaturated soil-water characteristic curves in practical engineering.

Conclusions
This paper investigated the bimodal SWCC and pore structure of a compacted granite residual soil using several laboratory tests.Nine initial conditions were selected based on the proctor curves, and the effect of the initial conditions on the PSD and SWCCs was thoroughly examined.Based on the results and findings, the following conclusions can be derived:

•
The tested granite residual soil presents a bimodal pore size distribution and bimodal SWCC for all the initial conditions, demonstrating a clear correlation between pore structure and SWCC variations.In contrast, intra-aggregate pores remain relatively unchanged, resulting in less pronounced differences in SWCCs within the high-suction range.

•
Incomplete filling of aggregates by adsorbed water leads to noticeable changes in pintra values, and the SWCC is significantly affected by the compaction energy.After the saturation of aggregates, the stabilized pore structure reduces sensitivity of the SWCC to the compaction energy.

•
The proposed 3D bimodal SWCC model captures the main characteristics and variations of SWCCs among initial conditions in the Sr-s-e0 space and exhibits strong agreement with laboratory test data (R 2 > 0.98).

Figure 1 .
Figure 1.(a) The grain-size distribution of tested soil and (b) the proctor curves under different compaction energies.

Figure 2 .
Figure 2. The X-ray powder diffraction results of the granite residual soil.

Figure 3 .
Figure 3. (a) The placement of filter paper method and (b) the vacuum saturation fixture.

Figure 4 .
Figure 4. (a) The compacted soil specimens, (b) the glass beaker and (c) the sealed containers containing selected supersaturated salt solutions.

Figure 5 .
Figure 5. (a) The soil specimens for the MIP test, (b) the freeze-drying device and (c) the MIP test device.

Figure 6 .
Figure 6.Comparison of the pore size distribution under different proctor types: (a) modified proctor, (b) standard proctor and (c) reduced proctor.

Figure 8 .
Figure 8. Division of the pore structure based on the PSD curve.

Figure 9 .
Figure 9. Relationship between pintra and initial condition parameters: (a) compaction water content, (b) dry density, (c) degree of saturation and (d) fitting curve of pintra and degree of saturation.

Figure 10 .
Figure 10.Comparison of the Sr, a and Sr, c for all the initial conditions.

Figure 12 .
Figure 12.Determination of the parameters for the bimodal SWCC fitting curve.

Figure 16 .
Figure 16.Three-dimensional SWCC based on Equation (9) and all the test data.

Figure 17 .
Figure 17.Comparison of measured and predicted degree of saturation.

•
The size and volume of inter-aggregate pores change noticeably with the initial conditions, leading to changes of the SWCCs within the low-and middle-suction range.

Table 1 .
Three-dimensional SWCC equations in related research.

Table 2 .
Selected physical properties of the soil.

Table 3 .
The compaction water contents and dry densities of the tested initial conditions. :

Table 4 .
Relative humidity and suction of selected salt solutions.

Table 5 .
Summary of the bimodal SWCC fitting parameters.