Please wait a minute...
Journal of Arid Land  2020, Vol. 12 Issue (6): 937-949    DOI: 10.1007/s40333-020-0106-6
Research article     
Freeze-thaw effects on erosion process in loess slope under simulated rainfall
SU Yuanyi1,2, LI Peng1,2,*(), REN Zongping1,2, XIAO Lie1,2, ZHANG Hui3
1State Key Laboratory of Eco-hydraulics in Northwest Arid Region of China, Xi'an University of Technology, Xi'an 710048, China
2Key Laboratory of National Forestry Administration on Ecological Hydrology and Disaster Prevention in Arid Regions, Xi'an 710048, China
3China JIKAN Research Institute of Engineering Investigations and Design Co. Ltd., Xi'an 710048, China
Download: HTML     PDF(633KB)
Export: BibTeX | EndNote (RIS)      


Seasonal freeze-thaw processes have led to severe soil erosion in the middle and high latitudes. The area affected by freeze-thaw erosion in China exceeds 13% of the national territory. So understanding the effect of freeze-thaw on erosion process is of great significance for soil and water conservation as well as for ecological engineering. In this study, we designed simulated rainfall experiments to investigate soil erosion processes under two soil conditions, unfrozen slope (UFS) and frozen slope (FS), and three rainfall intensities of 0.6, 0.9 and 1.2 mm/min. The results showed that the initial runoff time of FS occurred much earlier than that of the UFS. Under the same rainfall intensity, the runoff of FS is 1.17-1.26 times that of UFS; and the sediment yield of FS is 6.48-10.49 times that of UFS. With increasing rainfall time, rills were produced on the slope. After the appearance of the rills, the sediment yield on the FS accounts for 74%-86% of the total sediment yield. Rill erosion was the main reason for the increase in soil erosion rate on FS, and the reduction in water percolation resulting from frozen layers was one of the important factors leading to the advancement of rills on slope. A linear relationship existed between the cumulative runoff and the sediment yield of UFS and FS (R2>0.97, P<0.01). The average mean weight diameter (MWD) on the slope erosion particles was as follows: UFS0.9 (73.84 μm)>FS0.6 (72.30 μm)>UFS1.2 (72.23 μm)>substrate (71.23 μm)>FS1.2 (71.06 μm)>FS0.9 (70.72 μm). During the early stage of the rainfall, the MWD of the FS was relatively large. However, during the middle to late rainfall, the particle composition gradually approached that of the soil substrate. Under different rainfall intensities, the mean soil erodibility (MK) of the FS was 7.22 times that of the UFS. The ratio of the mean regression coefficient C2 (MC2) between FS and UFS was roughly correspondent with MK. Therefore, the parameter C2 can be used to evaluate soil erodibility after the appearance of the rills. This article explored the influence mechanism of freeze-thaw effects on loess soil erosion and provided a theoretical basis for further studies on soil erosion in the loess hilly regions.

Key wordsunfrozen slope (UFS)      frozen slope (FS)      simulated rainfall      soil size selectivity      soil erodibility      loess hilly region     
Received: 21 March 2020      Published: 10 November 2020
Corresponding Authors:
About author: *LI Peng (E-mail:
Cite this article:

SU Yuanyi, LI Peng, REN Zongping, XIAO Lie, ZHANG Hui. Freeze-thaw effects on erosion process in loess slope under simulated rainfall. Journal of Arid Land, 2020, 12(6): 937-949.

URL:     OR

Fig. 1 Experimental device diagrams of test equipment real diagram (a) and schematic diagram of the simulated rainfall system (b)
Treatment Initial runoff time (min) Time of rill appearance (min) Average value
Runoff rate (mm/min) Sediment yield rate during the Ⅰ-stage (g/min) Sediment yield rate during the Ⅱ-stage (g/min)
UFS0.9 24.2 42 532.30 13.97 62.67
UFS1.2 9.0 33 878.62 19.67 121.46
FS0.6 10.1 27 233.45a 38.08a 125.72a
FS0.9 5.5 21 620.28b 134.24b 382.24b
FS1.2 2.6 12 1104.61c 139.65b 523.00c
Table 1 Slope runoff, sediment yield and its main moments
Fig. 2 Topography of the unfrozen slope (UFS; a) and frozen slope (FS; b) under 1.2 mm/min rainfall intensity
Fig. 3 Temporal variations in runoff rate (a) and sediment yield rate (b) under different rainfall intensities. The numbers, 0.6, 0.9 and 1.2, after UFS and FS represent the rainfall intensities (mm/min).
Treatment Ⅰ-stage MC1 Ⅱ-stage MC2
UFS0.9 y=0.033x-3.052 (R2=0.9975**) 0.027 y=0.049x-37.537 (R2=0.9986**) 0.066
UFS1.2 y=0.021x+11.128 (R2=0.9823**) y=0.083x-174.74 (R2=0.9786**)
FS0.6 y=0.269x+23.462 (R2=0.9865**) 0.247 y=0.424x-142.10 (R2=0.9968**) 0.487
FS0.9 y=0.297x-3.618 (R2=0.9976**) y=0.575x-469.78 (R2=0.9981**)
FS1.2 y=0.176x+67.248 (R2=0.9861**) y=0.463x-704.45 (R2=0.9972**)
Table 2 Cumulative runoff and cumulative sediment yield fitted equation under different erosion stages
Fig. 4 Temporal variations of the mean weight diameter (MWD) of particles (a) and the MWD (b). Different lowercase letters represent significant difference between different rainfall intensities at 0.05 level.
Fig. 5 Temporal variations of eroded sediment particle contents under different soil treatments
Treatment Clay Fine silt Coarse silt Fine sand
I-stage Ⅱ-stage I-stage Ⅱ-stage I-stage Ⅱ-stage I-stage Ⅱ-stage
UFS0.9 0.030 0.025 16.753 15.162 44.733 47.237 37.615 36.945
UFS1.2 0.037 0.033 19.113 16.410 43.639 47.176 37.303 35.439
FS0.6 0.045 0.072 16.272 17.375 45.065 46.465 38.395 35.990
FS0.9 0.057 0.163 17.299 21.176 46.507 45.611 36.857 33.473
FS1.2 0.027 0.110 12.960 19.543 45.340 45.160 40.569 34.512
Table 3 Average percentage (%) of the effective particle size of sediment under different treatments
Treatment A (kg/m2) MA (kg/m2) R (MJ·mm/(m2·h)) LS K (kg·h/(MJ·mm)) MK (kg·h/(MJ·mm))
UFS0.9 1.16 1.955 0.13 1.16 5.07 8.03
UFS1.2 2.75 0.14 1.16 10.99
FS0.6 3.45 11.150 0.07 1.16 27.15 58.00
FS0.9 12.15 0.10 1.16 68.22
FS1.2 17.85 0.13 1.16 78.64
Table 4 Soil erodibility factor (K) calculated from the rainfall simulation experiments for the UFS and FS
[1]   Asadi H, Ghadiri H, Rose C W, et al. 2007. An investigation of flow-driven soil erosion processes at low streampowers. Journal of Hydrology, 342(1-2):134-142.
[2]   Bissonnais Y L. 1996. Aggregate stability and assessment of soil crustability and erodibility: I. Theory and methodology. European Journal of Soil Science, 47(4):425-437.
[3]   Bochove E V, Danielle P, Pelletier F. 2000. Effects of freeze-thaw and soil structure on nitrous oxide produced in a clay soil. Soil Science Society of America Journal, 64(5):1638-1643.
[4]   Bullock M S, Nelson S D, Kemper W D. 1988. Soil cohesion as affected by freezing, water content, time and tillage. Soil Science Society of America Journal, 52(3):770-776.
[5]   Chang E H, Li P, Li Z B, et al. 2019a. Using water isotopes to analyze water uptake during vegetation succession on abandoned cropland on the Loess Plateau, China. Catena, 181:104095.
[6]   Chang E H, Li P, Li Z B, et al. 2019b. The impact of vegetation successional status on slope runoff erosion in the Loess Plateau of China. Water, 11(12):2614.
[7]   Cheng Y T, Li P, Xu G C. 2018. The effect of soil water content and erodibility on losses of available nitrogen and phosphorus in simulated freeze-thaw conditions. Catena, 166:21-33.
[8]   Cruse R M, Mier R, Mize C W. 2001. Surface residue effects on erosion of thawing soils. Soil Science Society of America Journal, 65(1):178-184.
[9]   Dagesse D F. 2010. Freezing-induced bulk soil volume changes. Canadian Journal of Soil Science, 90(3):389-401.
[10]   Edwards L M, Burney J R. 1989. The effect of antecedent freeze-thaw frequency on runoff and soil loss from frozen soil with and without subsoil compaction and ground cover. Canadian Journal of Soil Science, 69(4):799-811.
[11]   Farenhorst A, Bryan R B. 1995. Particle size distribution of sediment transported by shallow flow. Catena, 25(1-4):47-62.
[12]   Ferrick M G, Gatto L W. 2005. Quantifying the effect of a freeze-thaw cycle on soil erosion: laboratory experiments. Earth Surface Processes and Landforms, 30(10):1305-1326.
[13]   Gatto L W. 2000. Soil freeze-thaw-induced changes to a simulated rill: potential impacts on soil erosion. Geomorphology, 32(1):147-160.
[14]   Henry H A L. 2007. Soil freeze-thaw cycle experiments: Trends, methodological weaknesses and suggested improvements. Soil Biology and Biochemistry, 39(5):977-986.
[15]   Issa O M, Bissonnais Y L, Planchon O. 2006. Soil detachment and transport on field- and laboratory-scale interrill areas: Erosion processes and the size-selectivity of eroded sediment. Earth Surface Processes and Landforms, 31(8):929-939.
[16]   Kirkby M J. 1980. Modelling water erosion processes. In: Kirkby M J, Morgan R P C. Soil Erosion. Chichester: Wiley, 183-196.
[17]   Koiter A J, Owens P N, Petticrew E L. 2013. The behavioural characteristics of sediment properties and their implications for sediment fingerprinting as an approach for identifying sediment sources in river basins. Earth-Science Reviews, 125:24-42.
[18]   Layton J B, Skidmore E L, Thompson C A. 1993. Winter-associated changes in dry-soil aggregation as influenced by management. Soil Science Society of America Journal, 57(6):1568.
[19]   Li G Y, Fan H M. 2014. Effect of freeze-thaw on water stability of aggregates in a black soil of northeast China. Pedosphere, 24(2):285-290.
[20]   Li X, Jin R, Pan X D. 2012. Changes in the near-surface soil freeze-thaw cycle on the Qinghai-Tibetan Plateau. International Journal of Applied Earth Observation and Geoinformation, 17:33-42.
[21]   Li Z, Wu P T, Feng H, et al. 2009. Simulated experiment on effect of soil bulk density on soil infiltration capacity. Transactions of the CSAE, 25(6):40-44. (in Chinese)
[22]   Øygarden L. 2003. Rill and gully development during an extreme winter runoff event in Norway. Catena, 50(2-4):217-242.
[23]   Oztas T, Fayetorbay F. 2003. Effect of freezing and thawing processes on soil aggregate stability. Catena, 52(1):1-8.
[24]   Pawluk S. 1988. Freeze-thaw effects on granular structure reorganization for soil materials of varying texture and moisture content. Canadian Journal of Soil Science, 68(3):485-494.
[25]   Ren Z P, Su Y Y, Li P, et al. 2018. Runoff scouring experiment on sand-covered slope. Journal of Soil and Water Conservation, 32(3):29-35, 41. (in Chinese)
[26]   Sahin U, Anapali O. 2007. Short communication: The effect of freeze-thaw cycles on soil aggregate stability in different salinity and sodicity conditions. Spanish Journal of Agricultural Research, 5(3):431-434.
[27]   Sharratt B S, Lindstrom M J, Benoit G R, et al. 2000. Runoff and soil erosion during spring thaw in the Northern U.S. Corn Belt. Journal of Soil and Water Conservation, 55(4):487-494.
[28]   Shen H O, Zheng F L, Wen L L, et al. 2016. Impacts of rainfall intensity and slope gradient on rill erosion processes at loessial hillslope. Soil and Tillage Research, 155:429-436.
[29]   Shi P, Qin Y, Liu Q, et al. 2019. Soil respiration and response of carbon source changes to vegetation restoration in the Loess Plateau, China. Science of the Total Environment, 707:135507, doi: 10.1016/j.scitotenv.2019.135507.
[30]   Shi P, Feng Z H, Gao H D, et al. 2020. Has "Grain for Green" threaten food security on the Loess Plateau of China? Ecosystem Health and Sustainability, 6(1):1709560, doi: 10.1080/20964129.2019.1709560.
[31]   Shi Z H, Fang N F, Wu F Z, et al. 2012. Soil erosion processes and sediment sorting associated with transport mechanisms on steep slopes. Journal of Hydrology, 454- 455:123-130.
[32]   Su Y Y, Li P, Ren Z P, et al. 2020. Slope erosion and hydraulics during thawing of the sand-covered Loess Plateau. Water, 12(9):2461.
[33]   Sutherland R A, Wan Y, Ziegler A D, et al. 1996. Splash and wash dynamics: an experimental investigation using an oxisol. Geoderma, 69(1):85-103.
[34]   Ting J M, Torrence, Martin R, et al. 1983. Mechanisms of strength for frozen sand. Journal of Geotechnical Engineering, 109(10):1286-1302.
[35]   Wan Y, El-Swaify S A. 1998. Characterizing interrill sediment size by partitioning splash and wash processes. Soil Science Society of America Journal, 62(2):430.
[36]   Wang B, Zheng F, RöMkens M J M. 2013. Soil erodibility for water erosion: A perspective and Chinese experiences. Geomorphology, 187(1):1-10.
[37]   Wang L, Tang L L, Wang X, et al. 2010. Effects of alley crop planting on soil and nutrient losses in the citrus orchards of the three gorges region. Soil and Tillage Research, 110(2):243-250.
[38]   Wang L, Shi Z H, Wang J. 2014. Rainfall kinetic energy controlling erosion processes and sediment sorting on steep hillslopes: A case study of clay loam soil from the Loess Plateau, China. Journal of Hydrology, 512:168-176.
[39]   Wang L, Wang Y, Saskia K, et al. 2018. Effect of soil management on soil erosion on sloping farmland during crop growth stages under a large-scale rainfall simulation experiment, 10(6):921-931.
[40]   Wang S J. 2004. Characteristics of freeze and thaw weathering and its contribution to sediment yield in middle yellow river basin. Bulletin of Soil and Water Conservation, 24(6):1-5. (in Chinese)
[41]   Wang T, Li P, Ren Z P, et al. 2017. Effects of freeze-thaw on soil erosion processes and sediment selectivity under simulated rainfall. Journal of Arid Land, 9(2):234-243.
[42]   Wang T, Li P, Hou J M, et al. 2018. Response of the meltwater erosion to runoff energy consumption on loessal slopes. Water, 10(11):1522.
[43]   Wang T, Li P, Li Z B, et al. 2019. The effects of freeze-thaw process on soil water migration in dam and slope farmland on the Loess Plateau, China. Science of The Total Environment, 666:721-730.
[44]   Wang T, Li P, Liu Y, et al. 2020. Experimental investigation of freeze-thaw meltwater compound erosion and runoff energy consumption on loessal slopes. Catena, 185:104310.
[45]   Wang W B, Shu X, Zhang Q F. 2015. Effects of freeze-thaw cycles on the soil nutrient balances, infiltration, and stability of cyanobacterial soil crusts in northern China. Plant and Soil, 386(1-2):263-272.
[46]   Weindorf D C. 2006. Relationships between permeability and erodibility of cultivated acrisols and cambisols in subtropical China. Pedosphere, 16(3):34-41.
[47]   Wischmeier W H, Smith D D. 1965. Predicting rainfall erosion losses: a guide to conservation planning. Agriculture Handbook. Washington D C, United States Department of Agriculture, 537:5-8.
[48]   Wischmeier W H, Johnson C B, Cross B V. 1971. A soil erodibility nomograph for farmland and construction sites. Journal of Soil Water Conservation, 26:93-189.
[49]   Xiao L, Yao K H, Li P, et al. 2019a. Effects of freeze-thaw cycles and initial soil moisture content on soil aggregate stability in natural grassland and Chinese pine forest on the Loess Plateau of China. Journal of Soils and Sediments, 20(3):1222-1230.
[50]   Xiao L, Zhang Y, Li P, et al. 2019b. Effects of freeze-thaw cycles on aggregate-associated organic carbon and glomalin-related soil protein in natural-succession grassland and Chinese pine forest on the Loess Plateau. Geoderma, 334:1-8.
[51]   Xiao P Q, Yao W Y, Shen Z Z, et al. 2011. Experimental study on erosion process and hydrodynamics mechanism of alfalfa grassland. Journal of Hydraulic Engineering, 42(2):232-237. (in Chinese)
[52]   Zhang Y, Zhang H, Li Z B. 2018. Process of runoff and sediment yield and relationship between water and sand of frozen soil slope in loess area under different rainfall intensities. Transactions of the Chinese Society of Agricultural Engineering, 34(11):136-14. (in Chinese)
[53]   Zhang Y, Li P, Liu X J, et al. 2019. Effects of farmland conversion on the stoichiometry of carbon, nitrogen, and phosphorus in soil aggregates on the Loess Plateau of China. Geoderma, 351:188-196.
[54]   Zhao B H, Li Z B, Li P, et al. 2017. Spatial distribution of soil organic carbon and its influencing factors under the condition of ecological construction in a hilly-gully watershed of the Loess Plateau, China. Geoderma, 296:10-17.
[55]   Zhou L L, Wang T L, Fan H M. 2009. Effects of incompletely thawed layer on black soil slope rainfall erosion. Journal of Soil and Water Conservation, 23(6):1-4, 37. (in Chinese)
[1] GAO Zhiyong, WANG Xing. Spatial variability of leaf wetness under different soil water conditions in rainfed jujube (Ziziphus jujuba Mill.) in the loess hilly region, China[J]. Journal of Arid Land, 2022, 14(1): 70-81.
[2] LIU Weichao, FU Shuyue, YAN Shengji, REN Chengjie, WU Shaojun, DENG Jian, LI Boyong, HAN Xinhui, YANG Gaihe. Responses of plant community to the linkages in plant-soil C:N:P stoichiometry during secondary succession of abandoned farmlands, China[J]. Journal of Arid Land, 2020, 12(2): 215-226.
[3] Linhua WANG, Yafeng WANG, SASKIA Keesstra, ARTEMI Cerdà, Bo MA, Faqi WU. Effect of soil management on soil erosion on sloping farmland during crop growth stages under a large-scale rainfall simulation experiment[J]. Journal of Arid Land, 2018, 10(6): 921-931.
[4] Tian WANG, Peng LI, Zongping REN, Guoce XU, Zhanbin LI, Yuanyuan YANG, Shanshan TANG, Jingwei YAO. Effects of freeze-thaw on soil erosion processes and sediment selectivity under simulated rainfall[J]. Journal of Arid Land, 2017, 9(2): 234-243.
[5] Yi ZHOU, Bing GUO, ShiXin WANG, HePing TAO. An estimation method of soil wind erosion in Inner Mongolia of China based on geographic information system and remote sensing[J]. Journal of Arid Land, 2015, 7(3): 304-317.