Changes in the hydrological processes in alpine soil constitute one of the several key problems encountered with studying watershed hydrology and ecosystem stability against the background of global warming. A typically developing thermokarst lake was chosen as a subject for a study using model simulation based on observations of soil physical properties, infiltration processes, and soil moisture. The results showed that the selected thermokarst lake imposed certain changes on the soil infiltration processes and, with the degree of impact intensifying, the initial infiltration rate decreased. The greatest reduction was achieved in the area of moderate impact. However, the stable infiltration rate and cumulative infiltration gradually increased in the surface layer at a depth of 10 and 20 cm, both decreasing initially and then increasing, which is correlated significantly with soil textures. Moreover, the cumulative infiltration changed in line with steady infiltration rate. Based on a comparative analysis, the Horton model helps better understand the effect on the soil infiltration processes of the cold alpine meadow close to the chosen thermokarst lake. In conclusion, the formation of the thermokarst lake reduced the water holding capacity of the alpine meadow soil and caused the hydraulic conductivity to increase, resulting in the reduction of runoff capacity in the area of the thermokarst lake. 相似文献
Discrete element method has been widely adopted to simulate processes that are challenging to continuum-based approaches. However, its computational efficiency can be greatly compromised when large number of particles are required to model regions of less interest to researchers. Due to this, the application of DEM to boundary value problems has been limited. This paper introduces a three-dimensional discrete element–finite difference coupling method, in which the discrete–continuum interactions are modeled in local coordinate systems where the force and displacement compatibilities between the coupled subdomains are considered. The method is validated using a model dynamic compaction test on sand. The comparison between the numerical and physical test results shows that the coupling method can effectively simulate the dynamic compaction process. The responses of the DEM model show that dynamic stress propagation (compaction mechanism) and tamper penetration (bearing capacity mechanism) play very different roles in soil deformations. Under impact loading, the soil undergoes a transient weakening process induced by dynamic stress propagation, which makes the soil easier to densify under bearing capacity mechanism. The distribution of tamping energy between the two mechanisms can influence the compaction efficiency, and allocating higher compaction energy to bearing capacity mechanism could improve the efficiency of dynamic compaction.
