1. Introduction
Internal erosion is defined as the erosion of soil particles due to water flowing through the soil [
1,
2]. It occurs when the seepage forces of water are larger than the resistance forces of soil, thus forcing the soil particles to start moving due to the water flow. These seepage forces cause erosion in both natural soil deposits and hydraulic and embankment-fill structures. Additionally, erosion does not occur uniformly through a structure; however, it increases in places with a concentration of seepage and water velocity [
3]. Various modifications happen in the soil through the erosion process that led to alterations in the void ratio, hydraulic conductivity, and shear strength of the soil. Many factors affect the erosion of the soil. Factors include soil water content; unit weight; plasticity index; void ratio; soil swell; and percent passing sieve #200, soil pH, and clay minerals [
4]. Other studies [
5,
6] showed that internal soil erosion is affected by soil type, water and material properties, and soil–piping resistance relationships. Generally, these factors resemble the physical and mineral composition of the soil [
4].
Various types of internal erosion are recognized in hydraulic structures, such as piping, suffusion, blowout (heave, uplift), Suffusion, and backward erosion [
5,
6,
7]. Piping involves the formation of a continuous tunnel, and as a result, water will flow between the upstream and downstream of the earth structure. This eases the erosion of the surrounding soil. The heave resulted in an imbalance of seepage uplift force and buoyant weight of soil particles, causing the lifting of soil particles. Suffusion erosion occurs when the small particles are transmitted and pass through the pores of the larger cohesionless granular soils. Backward erosion takes place when soil particles are separated from each other by seepage forces, causing them to move backward from downstream to upstream [
7].
The main factors that affect the erosion phenomenon can be listed as the erodibility of soil, water velocity inside the soil mass, and the geometry of the earth-structure. Therefore, erodibility can be defined as the relationship between water velocity and the corresponding erosion rate experienced by the soil. However, there is a problem with the definition as water velocity is a vector quantity and varies along the flow, and it is theoretically zero at the soil–water intercede. Therefore, the best definition of erodibility can be designated as the relationship between the erosion rate (Ż) and the shear stress (τ), and it is noticed that there is a direct relationship between the erosion rate, shear stress, and the type of soil [
4]. Briaud [
4] stated that grain size controls the coarse-grained soil erosion and plasticity controls the fine soil erosion. Sherard et al. [
8] showed that soil piping resistance is higher in well-compacted high-plasticity clays more than well-graded coarse sand and sand gravel mixture [
9]. Another study by Ke and Takahashi [
9] demonstrated that non-cohesive gap-graded soil erodes more because of the reduction in the strength and efficiency of shearing resistant forces due to the absence of some grain sizes.
Other physical soil factors affecting the erodibility of the soil are soil laying at the site, the degree of compaction, and the homogeneity of the soil [
10]. Khoshghalb et al. [
11] investigated the low-compaction level on the internal erosion behavior of the soil. The study revealed that at low compaction level the soil does not necessarily erode more if the suitable water content is used. Skempton and Brogon [
12] and Tomlison and Vaid [
13] introduced new factors affecting the erosion of the soil such as water flow velocity, flow direction, hydraulic gradient, and possible chemical reactions in the soil. Furthermore, soil erodibility is affected by other factors such as the unit weight of soil, water content, void ratio, plasticity index, soil swell, pH of soil, clay minerals, and percent of soil passing sieve #200. Moreover, soil erosion depends on the physical and mineral composition of soil [
12].
Studies by Luthi [
14] and Banu and Attom [
5,
6] revealed that approximately 30–50% of embankment failures and are caused by the internal erosion of the soil. These percentages are considered a very high percentage and indicated that the majority of the earth fill dams and embankments are failed due to internal erosion. Therefore, the need to stabilize the soil against internal erosion is a vital step to prevent the catastrophic consequences of hydraulic embankment structure failure and save lives and properties.
Several materials and techniques have been used in soil stabilization such as cement stabilization [
15,
16], lime stabilization [
5,
6,
17], bituminous stabilization [
18,
19], electrical stabilization [
20,
21], grouting [
22,
23], fly ash stabilization [
24,
25], burned sludge stabilization [
26,
27,
28], sisal fiber stabilization [
29,
30], polypropylene stabilization [
31,
32], and compaction [
33]. These materials and techniques were used to increase the shear strength properties, as well as to reduce settlement, compressibility, and expansion of the soils. However, few studies were found in the literature about soil stabilization against internal erosion. A recent study by Rodríguez and Cardoso [
34] used bio-cementation through enzymatically induced carbonate precipitation treatment to prevent soil erosion. The study concluded that the use of such material will prevent the scour erosion. On the other hand, researchers have recently developed several tests to evaluate soil internal erosion such as the Pinhole Test [
8], Jet Erosion Test [
35,
36], Hole Erosion Test, and Slot Erosion Test [
37]. The development of these tests allowed for the investigation of the internal erosion phenomena in more depth and to understand the behavior of failure mechanism throughout the failure process. In addition, interested researchers can use these tests to study the effect of different materials on the internal erosion and to stabilize both cohesive and non-cohesive soils against internal erosion.
Clayey soils are the core materials used in earth fill embankments and pond projects. Such projects are mainly used for water retention and water diversion in addition to other hydraulic purposes. Since the main and common reasons behind failure of earth fill embankment is piping and internal erosion, it is essential to stabilize the clay against this phenomenon using a well-known and efficient stabilizer.
Lime is used in stabilizing soils against swelling and to improve the mechanical properties of weak soils [
5,
6]. Lime stabilizes the soil due to the soil–lime reaction that generates long-term cementitious pozzolanic bonding. This reaction takes place when the calcium in the lime reacts with the aluminates and silicates in the clay, resulting in a stable and stronger soil. However, using of lime at a large scale and long term has also been investigated. Most of the investigations revealed that lime could improve the chemical and biological properties of soils, and that the long term use of lime will increase the soil chemical properties [
38]. Another study showed that liming acid soils will keep the pH value suitable to grow crops [
39]. Additionally, the effect of freezing and thawing on soil stabilized with lime has been studied. The study revealed that the strength of highly plastic clay at 28 days curing time increased by 15 times and 3 times for low-plastic soil but decreased by about 10–15% after three freeze–thaw cycles [
40]. The microstructure was also influenced by the addition of lime and resulted in decreasing the pores [
41]. Recently, many researchers investigated the use of lime and lime mixed with different materials to stabilize the clay soil and improve its’ mechanical properties [
41,
42,
43,
44].
This research aims to assess the effectiveness of lime as a stabilizing agent against internal erosion in clay soils and to advance the understanding of internal erosion behavior. Two types of clayey soils, namely, CH and ML, were selected and mixed with Quick Lime at varying percentages by dry weight of the clay. The Hole Erosion Test (HET) was employed to evaluate the internal erosion indices of the soils. The findings of this study are expected to contribute practical insights towards constructing safer dams and more efficient earth fill embankments while mitigating the risks associated with soil internal erosion failure. Additionally, this research will pave the way for exploring the utilization of alternative stabilization materials, such as solid wastes and other environmentally hazardous substances, in civil engineering practices.
3. Materials and Testing Program
Two types of clay soils were selected, and the initial physical properties such as specific gravity, compaction parameters, Atterberg’s limits, maximum dry unit weight, and optimum moisture content were determined in accordance with ASTM standard procedures. The initial physical properties are shown in
Table 2.
In this research, an adequate number of specimens from the two soils were prepared in the hole erosion test cylindrical mold with 10 cm diameter and 10.16 cm height. The soil was passed through sieve no. 40 and then the two types of soil were mixed with different lime content by dry weight of the soil at five different percentages, namely, 1%, 2%, 3%, 4%, and 5%. Six identical samples were prepared and tested for each lime content at three different curing times (24 h, 48 h, and seven days). Two identical samples were tested, and the average value of the results was taken. Two samples from each soil with zero lime content were prepared and tested as reference samples. All specimens were prepared in the standard mold at 95% relative compaction and optimum water content conditions. A 6 mm diameter hole was drilled throughout the longitudinal direction of the samples to simulate a concentrated leak pipe of 10.16 cm in length.
Figure 1 shows a prepared specimen with the 6 mm drilled hole at the center of the specimen.
The HET was employed to evaluate the erosion rate of the compacted samples at each percentage for the two types of soils. In this test, the sample was placed between two different heads of water: inlet (upper-stream) and outlet (downstream). The water was allowed to flow throughout the hole with a difference in head between inlet and outlet of a value in the range of 500 to 1200 mm. The head value was chosen within this range to initiate erosion. The chosen head between the inlet and outlet remained constant throughout the test but might be changed after a minimum of 45 min of running test time if no failure in the sample takes place. The hydraulic gradient, the flow rate at any time, and the diameter of the hole at the end of the test were measured.
Figure 2 shows the schematic diagram of the Hole Erosion Test.
The diameter of hole at any time throughout the test can be determined using the recorded data. Additionally, the erosion rate per unit surface area and wall shear stress along the hole was calculated at any time of the test. The I
HET and the critical shear stress (τ
c) values were determined from the above measured and calculated parameters.
Figure 3 shows the steps of calculating the erosion rate index and the critical shear stress; all symbols are defined in
Table 3.
Figure 4 shows the friction factor (FT.
t) versus time. The relationship was linear. The friction factors throughout the test timing were determined from the linear equation shown in
Figure 5.
To obtain the diameter at any time,
Figure 4 was used. At a specific time, FT can be obtained, and the corresponding Q will already be measured during the test. All other variables were constant and known throughout the test. By applying Equation (1), the diameter at any time can be calculated.
is the diameter at time t.
To obtain the erosion rate index I
HET and the critical shear τ
c, the best fit line between the points relating the rate of erosion per unit area (ɛ
o) and the shear stress (τ
t) at different test timings was drawn. Additionally, the equation of the best-fit line was obtained. The results are shown in
Figure 5.
4. Results and Discussion
4.1. Effect of Adding Lime on Soil Erosion Parameters
Two types of clay soils: Soil #1 (CH) and Soil #2 (ML), with different percentages of lime, were tested. Soil 1 was mixed with 0.5%, 1.0%, 1.5%, and 2.0% of lime. Soil #2 was mixed with 1.0%, 2.0%, 3.0%, 4.0%, and 5.0% of lime. The percentage of added lime was defined by the dry weight of the soil.
For Soil #1, a sample containing 5% lime was tested after 1 day, and no erosion was noticed after 2 h of continuous water flow. The flow rate of the water path was almost constant throughout the test time, and the final diameter exhibited no change from the initial diameter. Hence, it was excluded from calculations, and the percentages for Soil #1 were reduced to be in the range between 0.5% and 2.0% at 0.5% increments. Moreover, both the 1.5% sample and 2.0% sample were stable when tested after 24 h of curing time; therefore, no further testing at 48 h and 7 days was conducted.
4.1.1. Effect of Adding Lime on the Diameter of the Hole
To visualize the effect of adding lime on the final diameter of the water path, the relationship between the final diameters of the water path versus the lime percentage on both soil types was plotted and is shown in
Figure 6. As shown in
Figure 5, the relation between the percentage lime added and the final diameter of the water path was inversely related. An increase in the lime percentage added to the soil would achieve a reduction in the final diameter of the water path and be closer to the initial diameter.
Figure 6 implies that 2% of lime by dry weight of Soil #1 reduced the change in diameter by 97%. This behavior was also noticed when 5% lime was added to Soil #2. For Soil #1, the diameter at 0% lime was 22 mm in comparison of 6.5 mm at 2.0% lime. Similarly, in Soil #2, the diameter was equal to 45 mm at 0% lime, while at 5%, the diameter was 9 mm.
4.1.2. Effect of Adding Lime on Critical Shear Stress
Critical shear stress (τ
c) was plotted against percentage of lime added to both soils, as shown in
Figure 7. It is noticed that an increase in lime percentage would lead to an increase in the critical shear stress of the soil. For Soil #1 at 2% lime, the critical shear stress reached 100 N/m
2 at 24 h of curing time. This was a 50% increase in critical shear stress at the stabilizing lime percentage of 2% compared to 50 N/m
2 critical shear stress at 0% lime.
Similarly, with a greater percentage, a 203% increase in critical shear stress was achieved for Soil #2 with the addition of 5% lime by dry weight of soil. The critical shear stress increased from 33 N/m2 at 0% lime to 100 N/m2 at 5% lime addition and 48 h curing time.
4.1.3. Effect of Adding Lime on IHET
The main output of the HET analysis is to determine the I
HET. I
HET is used to classify the soil based on its resistance to erosion. Wan and Fell [
37] classified soil based on the value of the I
HET as per
Table 1. I
HET value ranges from less than 2 (very rapid erodibility) to more than 6 (extremely slow erodibility).
Figure 8 below shows the relationship between lime percentage added and I
HET value at 24 h curing time. It is obvious, for both soils, that an increase in lime percentage led to an increase in I
HET value and hence an increase in the stability of the soil.
Figure 9 shows the overall I
HET value increase for Soil #1. The I
HET, with 0% lime, was 4.5, which is classified as moderately slow erodible soil as per Wan and Fell [
37]. With the addition of 2% lime by dry weight of the soil, the I
HET value increased to 5.7 at 24 h of curing time. The new value of the I
HET of 5.7 classifies Soil #1 as a very slow erodible soil.
Similarly,
Figure 10 shows the overall I
HET value increase for Soil #2. The I
HET of Soil #2, with 0% lime, was 2.9, which is classified as very rapidly erodible soil. With the addition of 5% lime by dry weight of soil, the I
HET value increased to 5.7 at 7 days curing time. The new value of the I
HET of 5.7 classifies Soil #2 as a very slow erodible soil.
Table 4 and
Table 5 summarize the results obtained in
Figure 7,
Figure 8 and
Figure 9 as the final results for Soil #1 and Soil #2, respectively. The tables present the erosion parameters including τ
c, I
HET, and Φ
f, noting that the initial diameter (Φ
0) was 6 mm.
Figure 11 shows the relationship between the final diameter of the water path and the I
HET in both soils. This shows that the larger the I
HET, the lower the final diameter of the water path in the soil sample and the less erodible the soil.
Table 6 and
Table 7 summarize the change in erosion parameters such as final diameter of water path (Φ
f), τ
c, and the I
HET for Soil #1 and Soil #2, respectively. A description of the erodibility of soil is given based on the I
HET value in the last column of the same tables.
The findings from
Table 6 and
Table 7 indicate that the addition of 2% lime by dry weight effectively stabilized Soil #1 against internal erosion. This conclusion is supported by the obtained I
HET of 5.7, which classifies the soil as very slow erosion. Furthermore, the HET conducted on a soil sample treated with 2% lime showed insignificant change in diameter or flow over a two-hour period, indicating sustained stability.
Similarly, for Soil #2, a 5% lime addition by dry weight proved sufficient for internal erosion stabilization. The resultant IHET of 5.7 underscores the soil’s classification as very slow erosion. Additionally, the HET carried out on the 5% lime-treated soil demonstrated consistent diameter and flow characteristics over a continuous four-hour period, reaffirming the soil’s enduring stability.
4.2. Effect of Lime Percentage and Curing Time on Soil Erosion Parameters
Careful consideration of curing time is essential in engineering projects involving materials, particularly those subject to chemical interactions when mixed. This period is crucial for enhancing the strength and stability of the materials utilized. In this study, Quick Lime served as a stabilizing agent for two types of clayey soils, known to achieve near-complete curing within the initial 48 h.
The testing of samples, incorporating varying percentages of lime, occurred at three distinct curing intervals: 24 h, 48 h, and 7 days. Subsequent sections will delve into the analysis of how curing time impacts erosion parameters. Moreover, the combined influence of increased lime percentage and extended curing duration will be examined.
The erosion parameters under scrutiny include the diameter of the Φf, τc, and IHET. These variables will be examined concerning both lime percentage and curing time to gain deeper insights into their interplay.
4.2.1. Effect of Lime Percentage and Curing Time on Diameter of the Water Path
Figure 12 and
Figure 13 below show the relationship between the final diameter of the water path versus the percentage of lime added at three different curing times (24 h, 48 h, and 7 days) for both Soil #1 and Soil #2. The diameter curve for both soils at 7 days is the lowest among the three curing curves. This indicates that the change in final diameter becomes lower with the increase in curing time.
4.2.2. Effect of Lime Percentage and Curing Time on Critical Shear Stress
Figure 14 and
Figure 15 present the relationship between the critical shear stress and the percentage of lime added at three different curing times (24 h, 48 h, and 7 days) for both soils, Soil #1 and Soil #2, respectively. The critical shear stress for both soils at 7 days was the highest among the three curing bars at each lime percentage. However, the increase in the percentages as shown in the figures implies that the lime gained most of its strength at 24 h curing time.
An increase in lime percentage will continue increasing the critical shear stress. However, if lime is added further than the lime fixation point, the soil will be stabilized, while the extra lime will work on increasing the soil strength. That is why at higher lime percentages, the critical shear stress value exhibited a larger jump between different curing times, as shown in the bar charts below.
4.2.3. Effect of Lime Percentage and Curing Time on Erosion Rate Index (IHET)
Figure 16 and
Figure 17 present the relationship between the Erosion Rate Index and the percentage of lime added at three different curing times (24 h, 48 h, and 7 days) for Soil #1 and Soil #2, respectively. The I
HET for both soils at 7 days was the highest amongst the three curing bars at every lime percentage. This indicated that the I
HET of the soil increased as the curing time increased.
Moreover, with a higher lime percentage, the escalation in the IHET diminished. This can be attributed to two primary reasons. Firstly, lime fixation was attained, implying that a further rise in lime content will not enhance soil stability against internal erosion; instead, it will bolster the soil’s strength. Secondly, quick lime typically undergoes full curing at 48 h, and extending the curing period beyond this duration will not yield a significant effect.
4.3. Effect of Soil Type and Gradation of the Soil on the Percentage of Lime Needed to Stabilize the Soil against Internal Erosion
As noted in
Section 4.1 previously, Soil #1 is categorized as CH soil, while Soil #2 falls under the ML soil classification according to the Unified Soil Classification System. CH soil is characterized by its clayey nature with low plasticity and a high clay content, evident from the grain size distribution presented in
Table 2, which indicates that Soil #1 contains 65% clay. In contrast, Soil #2, designated as ML soil, is described as low-plastic silt with only 17% clay content. This discrepancy elucidates why Soil #1 required only 2% lime of dry mass of soil for stabilization, whereas Soil #2 necessitated 5% lime. The lime primarily reacts with the clayey constituents within the soil, accounting for the disparity in required percentages. Consequently, the soil type and the proportion of clayey material therein significantly influence the amount of lime needed to mitigate internal erosion.