Applying Biomineralization Technology to Study the Effects of Rainfall Induced Soil Erosion

The rainless days and drought seasons reveal a tendency to lengthen the wet and dry period in recent years in Taiwan. In the bare riverbeds in central Taiwan, such as the Dajia and Zhuoshui rivers, fugitive dust is the common problem during the winter’s dry period with northeast monsoon. The study aims to use the biological method, Microbial-induced carbonate precipitation (MICP), to solidify the soil and implement a series of rainfall simulator experiments to reduce the Aeolian dust emission problems. Accordingly, the relationship between rainfall-induced soil erosion and its soil specimens are also discussed. The soil conducted the MICP for seven days as the curing age, and then be analyzed the degree of soil solidification under different conditions, rainfall intensity of 41 mm/h and 61 mm/h by the soil erosion experiment. The effect of soil solidification with various relative density of soils of 60%, 70%, and 80%, and hillside slopes of 5°, 10°, 20° were tested. The result indicated that, the higher the relative density of soils, the better the effect of soil solidification would be. The relative density of soils from 60% to 80% all kept the effect of soil solidification as applying to MICP. Therefore, it was important to select the curing age with the matching relative density of soils. Moreover, the most appropriate condition for the effects of soil solidification by MICP was the slopes below 10° and the curing age for seven days. The effect of soil solidification was still preserved in that with the high intensity rainfall (66 mm/h) due to the Aeolian dust emission commonly happening in the bare riverbeds with mild slopes.


Introduction
Taiwan is surrounded by the volcano area of Pacific Ocean and it is located in the boundary polymerization of the Philippine Sea plate and Eurasian plate. The high frequency of seismic activity, mountainous island with steep landform result in frequent geological hazards that affect people's living and other security issues [1]. Besides, Taiwan is also located in subtropical zone of the western Pacific Ocean with a typical island climate. During the raining season, flash floods cause severe fluvial erosion. A large quantity of soil is accumulated to downstream regions, or even to the estuaries. In fact, serious riverbed erosion and soil loss happen due to the fragile shale in Taiwan. However, a large area of bare riverbed is formed due to insufficient river flows later in the drought season. Consequently, fine particles (Aeolian deposit) are lifted when strong winds are blowing, and transported by saltation and suspended in the northeast monsoon. The air quality is getting worse, especially in western central Taiwan, as a result of Aeolian dust [2]. Another reason that caused severe Aeolian dust is the 921 Earthquake. 921 Earthquake, also known as the Chi-Chi earthquake, was the most serious natural disaster in Taiwan. It was a 7.3 Richter magnitude scale earthquake that occurred in Chi-Chi, Nantou County, Taiwan on Tuesday, 21 September 1999 at 01:47:12 local time. 921 Earthquake resulted in serious landslides in the upstream of hill stream in Taiwan. A great number of soils accumulated to downstream regions, or even to the estuaries, which led to the Aeolian Dust [1,2]. Subsequently, since the Chi-Chi Weir has operated in 2002, the downstream of weir has serious sediment and bareness problems. The accumulation of mud and sand in the downstream is getting serious, and then Aeolian and fugitive dust directly affect the living conditions of residents. Overgrazing is also triggering and acceleration soil erosion processes, which leads to significant land degradation for agricultural lands, sedimentation of adjacent rivers and lakes, loss of livelihoods, and property damage [3]. Therefore, the innovative techniques and strategies to prevent soil erosion or reduce erosion rates is definitely a crucial issue [4].
The analysis of potential river dust events was conducted from 1994 to 2015 over western Taiwan, and it indicted that the significant river dust event with unprecedented PM10 concentration was examined [5]. Additionally, a similar study pointed the PM2.5 concentrations during the Aeolian dust episodes (ADEs) were 3-3.6 fold higher than those after the ADEs [6]. In the present situation, the prevention strategies for fugitive dust emission in Taiwan includes vegetation planting, bare surfaces watering, or soil adhesive spraying, etc., for example, planting is suitable for soil consolidation and it might prevent the Aeolian dust lifting and achieve the effect of sand-gathering. Additionally, construction of windbreaks along the river tends to reduce wind speed, so that it might prevent Aeolian dust to transport, and then be blocked by the windbreak. However, the windbreaks that are planted in impoverished and bare soil are often washed away by the floods. Thus, the bare surfaces watering must be carried out with high frequencies under strong winds condition. This study is adopting a full-scale rainfall simulator to directly simulate the rainfall, which has reached a maximum speed. In addition, the on-site soils are collected to conduct the experiment of erosion resistance. The methodology that is developed in this study is expected to apply to main estuary downstream in order that the serious eolian dust can be mitigated in Taiwan.
In the 1960s, American geologist Boquet pointed that bacteria could result in soil solidifying during the metabolism process. The bacteria that were used in this study were Bacillus Pasteurii DSM33 (renamed Sporosarcina Pasterurii DSM33), which belong to Gram-positive bacteria, and were commonly applied to research into Microbial-Induced Carbonate Precipitation (MICP). The bacteria were ubiquitous in the soil, and they usually produced a large amount of urease. With a series of chemical and environmental reactions, the bacteria combined with Ca 2+ to become CaCO3 or other cementing materials [7,8]. The microbial induced carbonate precipitation can enhance the resistance abilities to shear stress for moist sandy soil [9]. As mentioned above, the application of MICP is a natural process that is aided by living organisms and an eco-friendly, economical treatment. Therefore, it is an environmentally friendly engineering method [10]. An improved ecological method, carbonate precipitation by reaction of biological species (Microbial-Induced Carbonate Precipitation, MICP), is applied in this study. As the advanced technique, the MICP adopts the Microbial Metabolism to solidify the soil and make changes in its physical properties, such as shear strength, consolidation, porosity, erosion resistance, etc. [11][12][13]. The principle of MICP is that the urea hydrolysis is produced by microbial catalyzing to carbonate. The study pointed out that used native bacteria present in the soil to precipitate calcium carbonate can result in improvements in the mechanical properties of the soil, which allowed for the mitigation of the erosion processes that characterize the soil profile [14]. This technique can improve and enhance the solidification, as well as the physical properties of soil. Solidifying loose soil by MICP in the bare surface reduces the soil erosion as well as Aeolian dust emission. In short, the biotechnological technique is extensively applied in many engineering projects since environmental awareness has inreased in recent years [15][16][17][18]. The biological techniques were confirmed with more effective than other techniques [19], and is indeed a potential technology that has improved the engineering efficiency [15].
This study conducted experiments of rainfall simulator in a laboratory in order to apply biomineralization technology to the research of rainfall-runoff erosion. Approximately forty different kinds of rainfall simulators were used widely in a lot of erosion research due to the higher precisely experimental controllability and the retesting convenience. The major differences among 80% of simulated rainfall were the characteristics of nozzles. The Kentucky rainfall simulator was commonly used in the experiments. It is a small and simple facility with a diameter of 15.24 cm rainfall area [20].
A single-drop Rainfall simulator was used to evaluate the process of specific droplet, which induced the erosion [21,22]. A multi-drop rainfall simulator was applied to a laboratory and field scale study [23]. These above categories can be further classified into pressurized rainfall simulator and rainfall simulator with dripper-type needle. The pressurized rainfall simulator was used to simulate the rainfall intensity by pump changing the pressure. A lot of related research, such as rainfall properties comparison, in situ testing, and the effect of sediment concentration, was undertaken [23][24][25]. Rainfall simulator with dripper-type needle was used with the needles which control the falling droplet by gravity. Thus, the rainfall intensity was controlled by the flow discharge. It could be used to estimate the flow velocity which affects sediment [26], evaluates the impact of soil erosion [27,28], and analyzes the correlation between infiltration erosion and rainfall-runoff [29].
In this study, the experiment aimed to evaluate the effects of solidification by MICP and compared solidification and non-solidification by the means of destructive degree of erosion, the amount of soil loss, and the erosion resistance rate. Firstly, we adopted rainfall simulator to produce two extreme rainfall intensities: 41 and 66 mm/h to evaluate the effect of MICP technique that was applied on soil erosion prevention. Secondly, the analysis of soil erosion degree, amount of soil loss, and indexes of erosion resistance on both solidification and non-solidification were conducted. The paper introduced the research motivation and examined some lectures in the first section. Subsequently, the study presented the experimental instruments and some works of calibrations.

Research Method and Design
Different experimental variables, such as rainfall intensities, relative densities of soil, and hillside slopes, were designed to discuss the effect of MICP conducting on surface soil. The bacteria this study used were Bacillus Pasteurii DSM33 and they were commonly applied to research MICP. The bacteria were ubiquitous in the soil, and thy usually produced a large amount of urease. Urease catalyzes the hydrolysis of urea and then it would produce CO2, NH4 + , and OH − . The pH value would gradually increase, as OH − was released. The process of MICP was that CO2 converted to CO3 2− , and the CO3 2− combined with Ca 2+ to become CaCO3. The hydroxyl ions (OH−) that were generated from the conversion of ammonia to ammonium result in the increase in local pH value, which leads to the decomposition of bicarbonate to carbonate ions. The bacteria that were used in this study were characterized as positively aerobic. It could induce calcium carbonate deposition under the condition of coexistence with urea and CaCl2. It could repair soil cracks and form an acid-base as well as a corrosion-resistant protective layer on the surface of the soils. The urease was still active, even if the bacteria died. Meanwhile, the reaction of inducing calcium carbonate deposition continued. Table 1 shows the experiment settings. Firstly, solidifying soil was added with both Bacillus Pasteurii DSM33 and calcium carbonate; the curing age of soil was seven days. In contrast, nothing was added to the non-solidifying soil. Subsequently, the slopes that the study presents take account of the sand drift issues at the estuaries in the future; thus, the moderate slopes, 5°, 10°, and 20°, are adopted. The relative density ranges (60%, 70%, 80%) are chosen because the comprehensively evaluation and the reference of previous study, , are adopted [30]. The experiment was conducted with a self-designed box for more accurate collecting erosion soils. This was a soil erosion box, which was 30 cm × 23 cm × 40 cm, with a mesh rack that was fixed inside, which was 15 cm × 15 cm × 15 cm, as shown in Figure 1a,b. The experiment was conducted with equipment which was produced by ourselves, a soil erosion box with a mesh rack fixed inside. They were both designed to collect the soil form eroding so that it could be more precise and convenient to measure the total mass of washout soil. In general, the Standard Proctor Compaction Test could obtain the values of the maximum dry density, minimum dry density, optimum moisture content, and all of the above values could determine the relative density of soil. After obtaining the maximum dry density values, the necessary amount of experimental soil with the relative density of 60%, 70%, and 80% were produced, respectively. The way of conducting the experiment was that two cameras were used to record the experiment processes and the experimental rainfall lasted for one hour. By different experimental settings, firstly, the photos were taken every minute to observe the soil erosion. Secondly, the ultrasonic distance measuring device was applied to measure the elevations of surface soil specimen. Thirdly, the destructive degree of soils erosion were analyzed and classified. Fourthly, the amounts of soil specimen loss after erosion were then measured and analyzed. Finally, the destructive degree of soil specimens with the MICP and non-MICP were also analyzed. The ultrasonic distance measurement device is one of the best instruments for sensing proximity and detecting distances with high reliability. The experiment utilizes it to scan the whole sample surfaces to obtain the distance variations and estimate the surface erosion levels. The washout soils are collected by the designed soil erosion box to obtain the weight (g) after the soils are dried.
The major reason for soil erosion was taken as rainfall intensity and its kinetic energy, while it cannot fully explain the cause of soil erosion. Therefore, it is difficult to clearly define the amount of soil erosion [31]. Each specimen was set a distance of 5 cm from one another, which is shown on Figure 2, to efficiently observe the erosion and calculate the soil loss by camera. After the rainfall simulator was ceased, top clean water in the erosion box was removed.
The mass of washout soil was gathered and it was put into the oven with 105°C heat for 24 h. Lastly, the soil after heat (dried soil) was taken as the amount of soil loss and it was compared with the original amount of soil; therefore, the Erosion Resistance rate could be calculated, which is shown as Equation (1).
where LSn is the loss of non-solidification soils and LSs means the loss of solidification soils.

Rainfall Simulator and Calibrations
Generally, rainfall erosivity represents the primary driver for particle detachment in splash soil erosion [27]. The study conducted the experiments to evaluate the rainfall-runoff by rainfall simulator. The rainfall simulator was used in this experiment due to its higher precisely experimental controllability and retesting convenience. A rainfall simulator and dripper-type needle were applied in this study, as shown in Figure 3. This instrument was placed in the laboratory of National Chung Hsing University. The available distances between the rainfall simulator and the soil specimens was 8.5 m. The rainfall reaches terminal velocity with sufficient distances. The calibration demonstrates that the distance of 8.5 m is available for obtaining the rainfall's terminal velocity. Additionally, the available experimental area for rainfall was about one by four meters square. Figure 4 shows a sketch of the configuration of this rainfall induced erosion experiment. It can be seen that the rainfall simulator connected with the motor to control the droplet size and conduct the rainfall uniformity test. All of the above strategies could be adjusted by turning the motor.  Rainfall simulators have often been used to mimic natural rainfall for studies of various landsurface and water interaction processes [28]. However, the simulator could not fully evaluate the realistic soil erosion [26], thus the result that was evaluated by the rainfall simulator should be retested to obtain the more accurate data [32]. There were three sets of experiments to ensure the accuracy and availability of experiments, rainfall uniformity, rain droplet size, and rain drop terminal velocity. Therefore, the very set is tested twice and then an average is taken. For the Rainfall Simulator, its rainfall intensity was controlled by the gate of discharge (liter per millimeter, L/mm), the rainfall uniformity was regulated by the frequencies of turning motor (revolutions per minute, rpm), and the rain drop terminal velocity was calculated by high-speed photography method which would elaborate in the following section.
The calibrations were conducted by giving different amount of discharges that were controlled by the outlet valve to obtain their corresponding rainfall intensity in order to obtain the most appropriate rainfall intensity and the uniformity by the rainfall simulator. Additionally, the minimum relative error of rainfall uniformity was measured by different frequencies of turning motor (rpm). The minimum relative error could be treated as the most suitable configuration for the specific inputs. Figure 5 shows the experimental results, which presents the relative errors versus the revolution of motor with various outflow.
. It could be clearly seen that the minimum error was 600 (rpm) as the outflow rose from 0.5 to 2.0 (L/min.). On the top of that, the best performance was as the outflow was 2 (L/min.), the minimum error would be 800 (rpm). Table 2 lists the configuration results of most suitable revolution of motor per minute and rainfall intensity. In general, using suitable motor rotation to generate rainfall intensity could reach a more uniform state of rainfall. According to the definition of Central Weather Bureau in Taiwan, heavy rain in low-lying areas is 40 mm/hr. Thus, the study adopts the 41mm and 66mm (higher than the heavy rain) to process the experiment based on the distinctive rainfall pattern in Taiwan. Rainfall droplet sizes and its terminal velocity was required to verify before experiments. The study used the high-speed photography method to measure the mentioned rainfall properties. The moving tracks of droplets can be captured by the high-speed camera and particle sizes, and terminal velocity could then be calculated. The stain method that used the dying filter paper was applied to determine the droplet size.
The approach was based on the assumption that the droplet was proportional to stain diffusion area. The study used the droplet mass to calculate its size and the corresponding stain diffusion area. The experimental results could estimate a regression between the droplet size and the stain diffusion area. This experiment used the filter paper of Advantec NO.1 and then soaked it in the dye. The dried filter paper was put on the experimental zone that was under the rainfall simulator. As the rain droplet dripped down the filer paper, the area of stain diffusion was calculated. The droplet size could be obtained via substituting the area into the regressions, which was shown as the third column of Table 3. Theoretically, the terminal velocity of droplet can be estimated via the following equation [33]. The study calculated the theoretic terminal velocity by substituting droplet size from the third column of Table 3, as listed on the forth column of Table 3.
where D is rain drop diameter in which a sphere is assumed and d C is the dimensionless drag coefficient. The w ρ and a ρ are the densities of water and air, respectively.
High-speed photography method was adopted by Particle Image Velocimetry (PIV) in order to calibrate experimental terminal velocity in this study. Using a high-speed camera, Phantom Miro ex4 which is manufactured by AMETEK Company at Wayne, N.J., USA, as shown as Figure 6, in PIV results in accurate, detailed, and fast images for motion analysis. It has specially designed CMOS sensors that are capable of varying frame rates, resolutions (800 × 600 pixel), and exposure times (2-788 μs). It can be adjusted according to the different experiments settings to present the vivid and clear experimental observation and data. The PVI is utilized and set for 600 photos in one second in the experiment, with 883 photos in total.
The trajectory of raindrop is presented as bright dot without turbulence and elongation. The terminal velocity is obtained from the function of landing elevation and droplet size. In addition, the 8.2 m is determined as the landing elevation conducted in the experiment. Thus, five pictures are used to find a same droplet for the displacement by comparison and measurement to achieve the terminal velocity. Lastly, the terminal velocity is the displacement divided by the time interval. The experiments were conducted while using Phantom Miro ex4 high-speed camera combined the zoom lens of Nikon 50 mm. The pictures were set 600 photos per second. The study used a UTOPIA laser, manufactured by UTOPIA Instruments Company at Taipei, Taiwan, to bright several spots on the picture for taking their trajectory. The results used five successive photos to identify the same raindrop and the displacement of the same rain drop was measured, as shown as Figure 7. Subsequently, the value of terminal velocity was calculated by the displacement amount dividing the time interval. These values are experimental terminal velocity, which are listed in the fifth column of Table 3, and its errors are calculated in the sixth column. The data indicated that theoretical velocities were between 7.88 and 8.18 (m/s), and the experimental results were between 7.50 and 8.25 (m/s). That revealed that the errors between the theoretical and experimental values were less than 10%, and some sets of experiments were even less than 1%. Therefore, it can be assured that experimental terminal velocity can be treated as a theoretical value.

Field Sample
The topsoil samples (<30cm) were collected from Fongshan river in Hsinchu, as shown in Figure  8, which was a typical type in Taiwan. The soil is rich in iron oxide minerals, aluminum, and clay minerals. In addition, the greater plasticity, permeability (silt content: 16.69%), and physical properties, as well as clear configuration, are also its characteristic. However, the soil pH is revealed as acidic (5.1-5.5); therefore, it has poor fertility. Afterwards, the samples are dried at 105° for 24 h and then weighed so that the soil filter experiment can be conducted. Lastly, the soil specimen filtered by #8 (<2.36mm) was utilized to process the Standard Compaction Test. The specimen was made as a cubic of 15 cm by 15 cm by 5 cm and the relative density of specimen amounts that the study required were 60%, 70%, and 80% which could be calculated while using the following equation.
where max d γ represents the maximum dry density of soil (g/cm 3 ); is minimum dry density of soil (g/cm 3 ); Dr is relative density (%); and, V is volume of soil model (L 3 ). The process of using calcium to induce the deposition of calcium carbonate was a natural and eco-friendly process. The shape of Bacillus Pasteurii was close to a round rod or a spiral and its cell diameter was about 0.5 μm~3 μm. The diameter of some strains could be as small as 0.2μm. Therefore, the bacillus strain was more suitable for sandy soil with a hole size larger than 6μm. It is ideally suited to a cell diameter greater than 6μm environment in soils, especially in sandy soil. However, the clay soil is adequate for other smaller bacteria. The study used the soil from Fengshan creek basin, which was a suitable environment for bacillus.
According to previous studies, raindrop splash that enlarges the pore on the soil surface results in sever soil erosion. However, applying MICP to soil presents significant anti-erosion capacity which declines the permeability to prevent the erosion [34]. Moreover, on the basis of Maleki et al. (2016), the erosion of soil samples is experimentally investigated in a wind tunnel under the condition of wind velocity of 45 km/h [35]. The application of MICP on the soil surface can be an effective alternative for the soil wind erosion control. The precipitation of calcium carbonate on the surface is not only reduction in permeability to provide an impervious layer, but also adequate strength for the prevention of erosions. Therefore, the soil still has protective effects as the microbial induced carbonate is presented. The density is taken as a variable to examine whether the rainfall alters the density of soil or not for the purpose of understanding the effect of solidifying with different density of soil. The different slops are also taken as variables in the experiment, because the adequate gradient, splash erosion, runoff erosion, and variation of MICP application are crucial for the study. Therefore, the following experiments and discussions were conducted.

Relative Density and Erosion
This section analyzes the surface erosion from different relative densities of specimen, both solidified and non-solidified. This group of experiments were examined with the curing age for seven days, relative density of 60%, 70%, 80%, one-hour constant rainfall, rainfall intensity of 41 (mm/h), and 66 (mm/h) to discuss the soil erosion differences between the control group (non-solidified) and the test group (solidified). The slope for these sets of experiment is zero.
The elevation difference and erosion resistance rate were required in order to evaluate the effectiveness of solidification. The experimental surface elevations of non-solidification subtracted from solidification were their differences as elevation. The degree of surface erosion will affect the elevation of the specimen. If the degree of surface erosion is higher, then the elevation of the specimen that is measured will be lower. Figures 9 and 10 show that, by different rainfall intensities, the erosion resistance experiment presents the results in different relative densities, which is also a reference of relative density for the further experiments in the future. Figures 9 and 10 plot the results of elevation differences. The positive result represented that the erosion of the non-solidifiedspecimen is more severe than the solidification one. Then erosion resistance rate was calculated by the loss of soil. The result was shown in Tables 4 and 5 by comparing the elevation differences with erosion resistance rate.   The relative density of soil of 80% revealed the best resistance among others to compare the difference in erosion between control group and testing group in Tables 4 and 5. The solidification effect is also the most significant on both rainfall scenarios. It also revealed that the more solid the soil was, the more difficult the soil to be dispersed. The experimental results then indicated, with the intensity of rainfall increase, the ability to resist erosion was decreased. However, the relative density of 60%~80% specimen still had the effect of solidification with high intensity rainfall (66 mm/h). The MICP could effectively improve the shear strength for soil [36]. The improvements varied with soil densities. With MICP treatment, the improvement in shear strength of the soil specimens are significantly higher, as the density of soil specimens are higher. In addition, the relative density of soils significantly affected the amount of calcite precipitation, which was measured to be approximately 40%-50% higher in 80% relative density soil than that in 60% relative density soil. Our experiments, furthered to discuss the rainfall intensity by using the rainfall simulator to simulate the real rainy scene, also showed similar results to Ng et al. (2013) and Kim et al. (2016), where the higher relative density samples revealed better resistance to soil erosion [36,37].
According to the destructive degree of surface soil, the values for comparison from the control group (non-solidification) were subtracted from the testing group (solidification), as shown on Figures 9 and 10. The red color represented the accumulation of soil. Blue color represented the erosion of soil. The relative density of soil 80% revealed that particularly significant solidification was also shown in the figure. The higher the relative density, the more complete the soil specimen. Additionally, the specimen was evenly destructed with less erosion holes. The distribution of red is presented significant in the figure. With comprehensive analyzing, in the case of rainfall intensity 41mm/h with the relative density of soil, 80% revealed particularly significant values in solidification.

Slope and Erosion
In this section, the relative density of Dr = 80% was fixed for two different rainfall intensities, 41(mm/h), 66(mm/h), and the specimens are placed on different slopes (5°, 10°, 20°) for rainfall strike experiments. Additionally, the destructive degree were classified and the effectiveness of solidification were analyzed.
Rainfall intensity I = 41 (mm/h), 66 (mm/h), and slope (5°, 10°, 20°) were used as the independent variables. With different rainfall intensities, while comparing the different slopes the specimen was placed with the destructive degree of solidification and non-solidification were analyzed. The experimental surface elevation of non-solidification that were subtracted from solidification with different slopes were their differences in elevation, and were then plotted as a graph. The positive value represented that the erosion of non-solidification specimen is more severe than the solidification one. Subsequently, the erosion resistance rate was calculated by the loss of soil. The result was compared the elevation differences with erosion resistance rate. This group of experiments were examined for the different hillside slopes of 5°, 10°, 20°, one-hour constant rainfall, rainfall intensity of 41 (mm/h) and 66 (mm/h) to discuss the soil erosion differences between the control group (non-solidification) and the test group (solidification), which Tables 6 and 7 show. The results with a bearing capacity of footing on slope of 1.5H:1V (33.69°) have significantly increased after MICP treatment based on Pusadkar et al (2017) [38]. Moreover, the results thatwere proposed by Shao et al (2017) for applying MICP technology on silt slope anti-erosion provide valuable enlightenments [34]. Under heavy rainfall, the microbial soil significantly demonstrates the anti-erosion effectiveness on the slope by the high strength of surface. Nevertheless, our study has to further consider various hillside slopes and rainfall intensities. By adding more variables in our experiment, it conforms to real rainfall scenario and even increase the accuracies of the data. The results showed that the hillside slopes of 5% that revealed the best erosion resistance for the sheet eroding were mainly a kind of rainfall splash erosion. The solidification was presented significant. When the hillside slopes increased with the runoff erosion effect, the ability of soil erosion resistance was decreased. The ability of soil erosion resistance was the weakest, as 20° of hillside slope. However, the value was still positive, with rainfall intensity I = 41 (mm/h), which indicated that it still had the effect of solidification. With the rainfall intensity increased, the value of soil erosion resistance with slopes 20° was negative. It indicated that the minimum of solidification was slopes 10°, with I = 66 mm/h. In short, as the slope elevates, the function of erosion resistance will alleviate. The phenomenon is found by the result. Thus, a negative erosion resistance rate represents that the soils have no effect of solidifying; for instance, the ability of soil erosion resistance is the weakest, with 20° of hillside slope, which has no effect of solidification, as shown in the result.
According to the destructive degree of surface soil, the control group (non-solidification) subtracted from testing group (solidification) were the values as comparison, as shown on Figures 11  and 12. The red color represented the accumulation of soil; Blue color represented the erosion of soil. The hillside slope of 5°revealed particularly significant in solidification was also shown on figure.
The higher the hillside slope, the less the effect of solidification. The rainfall simulator platform is wobbled by the motors to simulate the real, random rainfall. In general, the same trajectories of raindrops do not exist, which also the reason why the Figures 9-12 do not reveal density homogenous of surface erosion.

4.Conclusions
The study aims to discuss the effect of MICP conducting on surface soil with rainfall simulator experiments. Experiments with different rainfall intensities, relative densities of soil, and hillside slopes were designed. The results revealed that, firstly, the lower density the soil, the weaker the effect of solidification. Therefore, the soil with low density should be applied by the biological method, Microbial-induced carbonate precipitation (MICP) with conducting different curing ages to improve our understanding on the effect of solidification. However, the higher the density the soil is, the shorter the curing age can be. Thus, with the curing age for seven days and the relative density of 60%~80%, it still has a solidifying effect. even at the high intensity rainfall (I = 66 mm/h). Secondly, the steeper the hillside slope, the less the effect of solidification. The hillside slopes increased with the runoff erosion effect; the ability of soil erosion resistance was decreased. The ability of soil erosion resistance was the weakest, at 20° of hillside slope. With the curing age for seven days, the minimum of solidification was slopes 10° with the high rainfall intensity (I = 66 mm/h).
In short, according to those experiments, in the case of applying the MICP method, different densities of soil should correspond with the suitable curing age to demonstrate maximal efficiency on solidification. Meanwhile, the hillside slopes increased with the runoff erosion effect and the ability of soil erosion resistance was decreased. Therefore, the steeper the hillside slope, the less the effect of solidification.