Laboratory Characterization and Discrete Element Modeling of Shrinkage and Cracking Behavior of Soil in Farmland

Soil desiccation cracks are common in farmland under dry conditions, which can alter soil water movement by providing preferential flow paths and thus affect water and fertilizer use efficiency. Understanding the mechanism of soil shrinkage and cracking is of great significance for optimizing field management by crack utilization or prevention. The behavior of soil shrinkage and cracking was monitored during drying experiments and analyzed with the help of a digital image processing method. The results showed that during shrinkage, the changes in soil height and equivalent diameter with water content differed significantly. The height change consisted of a rapid decline stage and a residual stage, while the equivalent diameter had a stable stage before the rapid decline stage. The VG-Peng model was suitable to fit the soil shrinkage characteristic curves, and the curves revealed that the soil shrinkage contained structural shrinkage, proportional shrinkage, residual shrinkage, and zero shrinkage stages. According to the changes in evaporation intensity, soil water evaporation could be divided into three stages: stable stage, declining stage, and residual stage. Cracks first formed in the defect areas and edge areas of the soil, and they mainly propagated in the stable evaporation stage. Crack development was dominated by an increase in crack length during the early cracking stage, while the propagation of crack width played a major role during the later stage. At the end of drying, the contribution ratio of crack length and width to the crack area was approximately 30% and 70%, respectively. The box-counting fractal dimension of the stabilized cracks was approximately 1.65, indicating that the crack network had significant self-similarity. The experimental results were used to implement the discrete element method to model the process of soil shrinkage and cracking. The models could effectively simulate the variation characteristics of soil height and equivalent diameter during shrinkage, as well as the variation characteristics of crack ratio and length density during cracking, with acceptable relative errors. In particular, the modeled morphology of the crack network was highly similar to the experimental observation. Our results provide new insights into the characterization and simulation of soil desiccation cracks, which will be conducive to understanding crack evolution and soil water movement in farmland.