Introduction

The erosion of soils caused serious engineering problems of structures (i.e., landslide [17], bank failures [15], ground subsidence [16]) resulting to catastrophic and irreversible failures. Piping failure and rainfall erosion of the embankment are the most typical structural failures caused by soil erosion [13]. In 2005, the piping failure of a metro tunnel constructed in Taiwan was reported (Fig. 1a). They found that the main reason for the failure was the very high dispersive of the silty sand [14]. More recently, in 2018, an earthen and rockfill dam crumbled as a result of heavy summer rain alongside devastating consequences including loss of lives (13 persons) and numerous properties (300 persons homeless) [23].

Fig. 1
figure 1

a The cave-in in Kaohsiung, Taiwan b a sinkhole found in Seoul, South Korea

Another typical failure caused by erosion of the soil is a sinkhole (Fig. 1b). In recent years, many sinkholes have been discovered in Korea; about 70 sinkholes were reported between January 2012 and July 2014. Most of them were caused by the erosion of underground soils owing to the ruptured sewer pipes and heavy rainfall [12]. The most important factor among the common causes of these problems is the erosion characteristic of soils, but so far, most designs and constructions have not considered this part much in Korea.

To evaluate the erodibility of soils, many studies were conducted with different devices and testing methods: crumb test [8], jet erosion test [7, 10, 19], flow pump tests, and flume erosion tests [1, 18], pinhole test [2, 20, 21], hole erosion test [24].

Among these tests, the pinhole test is well-known as the most reliable for classification and identification of the erodibility of the soil since it directly measures the erodibility and dispersibility of the soil caused by the water flow in a small hole in the center of the soil sample [13, 21]. The advantage of this test is the small sample required, easy sample preparation, and testing procedure. According to ASTM (2006), a pinhole test is used to identify and classify the erodibility of the soil based on the flow rate of discharge water, cloudiness, the size of the hole, and the applied pressure. There are 6 levels of the dispersive ranging from very high dispersive soils to non-dispersive soils (D1, D2, ND4, ND3, ND2, ND1) [2].

This study investigates the erodibility of different soil types (clay, SM soil, SP-SM soil) using the pinhole test. The sieve analysis and hydrometer analysis were performed to evaluate the effect of the grain size distribution on the erodibility of the soil. Furthermore, the erosion coefficients and the critical shear stress were calculated as additional parameters for the evaluation of the erodibility of the soils.

Materials

In this study, soil samples were taken nearby the Nakdong river (Fig. 2). The grain size distribution curves, soil properties, and soil classification parameters are shown in Fig. 3 and Table 1, respectively. Based on the soil classification results, the tested soils are divided into three groups: SP-SM (namely RNB-2 and RNB-3), SM (namely RNB-10, RNB-11, and RNB-13), and undisturbed clay (namely RNB-23).

Fig. 2
figure 2

Soil sampling area

Fig. 3
figure 3

The grain size distribution of soils

Table 1 Basic soil properties and soil classification parameters

Experimental program

Pinhole test apparatus

Figure 4 shows the schematic drawing of the typical pinhole test apparatus. The specimen with the dimension of 32.8 mm in diameter and 38.0 mm in height was sandwiched with coarse sand layers at the top and bottom. The wire screen is placed between the specimen and coarse sand to prevent the penetration of the coarse sand into the specimen. A hole (1.2 mm in diameter) centered with a truncated cone is prefabricated by a No. 19 veterinarian hypodermic needle. At the top plate, the inlet valve is connected to the constant head tank in order to supply distilled water with a pH of 5.5–7.0. To measure the hydraulic head at the inlet, a pressure transducer with an accuracy of 0.01 kPa is employed. A stopwatch with a reading of 0.1 s is used to record the erosion time. The volume of outlet water is measured using graduated cylinders with capacities ranging from 10 to 100 mL.

Fig. 4
figure 4

Schematic diagram of the pinhole test equipment

Sample preparation

Before the commencement of the test, the disturbed soil (RNB-2, RNB-3, RNB-10, RNB-11, RNB-13) was dried with the oven and particles larger than 2 mm (No. 10 sieves) were removed. Afterward, the three layers of soils were compacted into the pinhole test cylinder until a total length of 38 mm is reached. It is important to recognize that the specimen should be compacted at natural water content to simulate its in-situ condition. A small rammer (350 g) with a dimension of 20 mm diameter, 120 mm length was used to compact the sample (Fig. 5). The samples were divided into three layers and one layer was compacted 20 times at natural water content to simulate its in-situ condition. Note that after compacted 20 times, the density of the compacted sample should be determined and compared with that of in-situ density. If the density of the compacted sample is higher than that of in-situ density, the compacted times are reduced until the difference between these densities lower than 5% and inverse. The density of the compacted sample can be determined by dividing the mass of the compacted soil by its volume. The mass of the compacted soil is determined by a weight scale while the volume of the soil is 32.09 cm3 (cylinder specimens with a diameter of 3.28 cm and a length of 3.8 cm). The truncated cone centering guide with a 1.5 mm diameter hole was inserted into the center of the top of the specimen using finger pressure (Fig. 5). The needle was inserted into the centering guide and forced through the soil sample. Then, the wire screen was placed on top of the specimen and the remaining void in the top of the test cylinder was filled with coarse sand. Finally, the top plate was assembled, connected with the distilled water tank, and the pressure transducer (Fig. 6).

Fig. 5
figure 5

Preparation of disturbed soil samples

Fig. 6
figure 6

Pinhole test set-up

For clay, two samples were taken from different depths [RNB-23 (16.5 m) and RNB-23 (18 m)], trimmed to a 38 mm long to fit snugly into the test cylinder. The test procedure for undisturbed samples is the same as that of disturbed samples (Fig. 7). It should be noted that the perimeter of the top of the clay sample was sealed with moulding clay to avoid the leach of water between the soil sample and the wall of the test cylinder. Besides, for some stiff clay samples, the penetration of the truncated cone guide nipple may be very difficult and results in the damage of the sample. Accordingly, the truncated cone guide nipple should be omitted in this case.

Fig. 7
figure 7

Preparation of clay sample

Experimental procedure

The experimental procedure is shown in Fig. 8. To measure the flow rate of the discharge, the outlet water was collected at an interval of 60 s. Note that if there is no outlet water at the start, the test should be stopped and then re-puncture the hole on the specimen. The cloudiness of the effluent, the flow rate, and the hole sizes are recorded as dispersive classification parameters. Based on these parameters, the soil can be classified as nondispersive (ND1 and ND2), slightly to moderately dispersive (ND3), moderately dispersive (ND4), dispersive (D2), and highly dispersive (D1) (Fig. 8).

Fig. 8
figure 8

Pinhole test procedure and classification

Results

Group 1: clay

Figure 9 shows the flow rate and dispersive grade of clay samples (RNB-23 (16.5 m) and RNB-23 (18 m)). The hole size of the samples increased from 1.2 at the initial to 2.2 mm at the final stage of the test under the final applied head pressure of 380 mm. However, the water collected from the test is completely clear (Fig. 10). Based on the results, the clay samples are classified as slightly to moderately dispersive soil and slow erosion (ND3).

Fig. 9
figure 9

Flow rate and dispersive grade of clay

Fig. 10
figure 10

a Hole size after test and b water collected from test

Group 2: SP-SM soil (RNB-2 and RNB-3)

The testing results of group 2 soils are presented in Fig. 11. Under the head of 50 mm, the final flow rate through the specimen is 1.33 mL/s after the testing time of 5 min. The diameters of the hole after the test were 3.73 mm (RNB-2) and 3.33 mm (RNB-3), which is approximately three times larger than those of the initial hole diameter (1.2 mm). Furthermore, the water collected from the test is very dark from the side and top of the cylinder (Fig. 12). As a result, RNB-2 and RNB-3 are classified as highly dispersive soil and rapid erosion (D1).

Fig. 11
figure 11

Flow rate and dispersive grade of SP-SM soils

Fig. 12
figure 12

a Hole size after test and b water collected from the test of SP-SM soils

Group 3: SM soil (RNB-10, RNB-11, and RNB-13)

The flow rate and dispersive grade of group 3 soils (RNB-10, RNB-11, RNB-13) are presented in Fig. 13. After 10 min under the head of 50 mm, the final flow rate through the specimens ranged from 1.17 to 1.25 mL/s. The hole size after the test varies from 2.40 to 2.56 mm, which is slightly smaller than that of group 2 soils (SP-SM). The flow at the end of the test from side and top is dark (Fig. 14). Thus, the SM soils are classified as rapid erosion soil. Since both the final flow rate and hole size after the test of group 3 size is lower than that of group 2 soils; therefore, as presented in Fig. 13, the group 3 soils (SM soils) is classified as dispersive soil (D2).

Fig. 13
figure 13

Flow rate and dispersive grade of SP-SM soils

Fig. 14
figure 14

a Hole size after test and b water collected from test

Discussion

Effect of grain size distribution on the erosion of the soils

Regarding the aforementioned results, SP-SM soils (RNB-2 and RNB-3) have the highest erosion level among three soil types. It is probably due to the fact that RNB-2 and RNB-3 have a low percent of the small particle size between 7.4 and 7.0% (Passing No. 200) respectively (Table 1). Furthermore, RNB-2 and RNB-3 curves are too uniform, i.e., Cu < 2. According to Holtz and Kovacs (2011), very poorly graded soils have a Cu of 2 or 3 such as beach sands, whereas very well graded soils may have Cu of 15 or greater [11]. Thus, the void between the large particles is not filled with small particles as shown in Fig. 15 [26]. Therefore, RNB-2 and RNB-3 are more susceptible to erosion based on the piping failure mechanism [13]. The erosion result of SP-SM soils also agreed very well with the backward erosion piping mechanism presented in previous studies [9, 22, 25]. The backward erosion initially progresses in small channels rather than a single large channel. They found that poor grade soils with low uniformity coefficients (Cu < 3) are subject to backward erosion piping [9, 22, 25]. In this case, the large void between uniform soil particles may consider as small channels, where the small soil particles are transported at the initial process of backward erosion piping.

Fig. 15
figure 15

Packing structures of a poorly graded soil (low Cu) and b well-graded soil (high Cu)

The frequency of each soil grain size is shown in Fig. 16. It is interesting to note that the grain size with the highest frequency of the rapid erosion soils (SM and SP-SM) ranged between 0.2 and 0.6 mm, while the slow erosion soil (CL clays) has the highest frequency of grain size ranging from 0.05 to 0.09 mm. Clearly, particles ranging between 0.2 and 0.6 mm of SP-SM soils (D1 class) are significantly larger than those of SM soils (D2 class). It is found that the grain size of 0.2–0.6 mm range is the most susceptible to soil erosion. The result of this study is in good agreement with the previous study [27].

Fig. 16
figure 16

Frequency of each grain size

Analysis of the critical shear stress of soils from the pinhole test data

The pinhole test uses the flow rate, diameter of the hole, cloudiness of the effluent, and applied pressure head to evaluate the erodilility of the soil. In this study, a process to analyze the erosion rate of soils by calculating the critical shear stress of soils from the pinhole test data is proposed. The empirical erosion law can be expressed as follows [5, 6, 24]:

$$ \dot{\varepsilon } = k_{er} (\tau - \tau_{c} ) $$
(1)

where \(\dot{\varepsilon }\) = erosion rate (kg m−2 s−1), τ = shear stress (Pa), τc = critical shear stress, and ker = erosion coefficient.

Firstly, the shear stress in Eq. (1) can be calculated from the force equilibrium equation as follows:

$$ 2\pi rL\tau = \pi r^{2} \Delta P $$
(2)

From Eq. (2), the shear stress τ can be calculated as follows:

$$ \tau = - \nabla P\frac{r(t)}{2} $$
(3)

where: L is the length of specimen (0.038 m), \(\Delta P\) is head pressure (0.05 m, 0.18 m, or 0.38 m); thus, \(\nabla P = - \Delta P/L\) is the applied pressure gradient in the direction of the flow, r(t) (m) is the hole size corresponding to the time t, and it can be calculated as follow [3,4,5]:

$$ r(t) = r_{i} \left( {\frac{{Q_{t} }}{{Q_{i} }}} \right)^{\frac{2}{5}} $$
(4)

where: ri (mm) is the initial hole size corresponding to the initial flow rate \(Q_{i}\)(m3/s), \(Q(t)\)(m3/s) is the flow rate corresponding to the time t.

Secondly, the erosion rate \(\dot{\varepsilon }(t)\) in Eq. (1) is calculated as follows:

$$ \dot{\varepsilon }(t) = \frac{{M_{So - er} (t)}}{A \times t} $$
(5)

where,

$$ M_{So - er} (t) = \rho_{s} L\pi (r_{(t)}^{2} - r_{i}^{2} ) $$
(6)

where: t = reference time (s), MSo-er = mass of soil eroded at time t (kg), r(t) = hole size corresponding to the time t (s), ri = initial hole size (m).

The critical shear stress analysis results from pinhole test data are shown in Fig. 17. The linear relationship between erosion rate and shear stress is observed regardless of soil types; thus, the proposed empirical erosion law (Eq. (1)) is reasonable for determining critical shear stress. The erosion coefficients and critical shear stress of soils calculated from proposed equations are presented in Table 2. The erosion coefficients of clay are approximately 13 times and 16 times smaller than those of the SM and SP-SM soils, respectively. Accordingly, the erosion coefficients can be used to distinguish between the low erosion (ND3 and ND4) and high erosion soils (D1 and D2). Besides, the lower dispersive class, the higher the critical shear stress is identified (Table 2). It can be carefully concluded that the critical shear stress might be used to classify the erosion of the soil such as τc > 3.5 Pa (ND3), 3.0 Pa < τc < 3.5 Pa (D2), and τc < 3.0 Pa (D1).

Fig. 17
figure 17

Critical shear stress of a clay and b SM and SP-SM soils

Table 2 Dispersive classification and critical shear stress of soils

Conclusions

  1. 1.

    This paper evaluates the erosion of the soil using a pinhole test. The clay, SP-SM, and SM soil were classified as slightly to moderately dispersive soil (ND3), highly dispersive soil (D1), and dispersive soil (D2).

  2. 2.

    The particle size distribution and coefficient of uniformity of soils are significant factors that affect the erosion of the soil. The well-graded soil has a low erodibility compared to that of the poorly graded one. The soil with a grain size ranging from 0.2 to 0.6 mm is most susceptible to soil erosion.

  3. 3.

    The erosion coefficients and the critical shear stress were calculated from the results of pinhole test as an additional parameter for classifying the erosion of the soil. The erosion coefficient is feasible to use to identify between the low erosion (ND3 and ND4) and high erosion soils (D1 and D2). Meanwhile, the critical shear stress might be used as a classification parameter of erosion of the soil as: τc > 3.5 Pa (ND3), 3.0 Pa < τc < 3.5 Pa (D2), and τc < 3.0 Pa (D1). However, this result should be further studied.