Abstract:In arid and semiarid regions, shrubs are often spatially distributed in patchy patterns, allowing shrubs to self-regulate and adapt to dry conditions. Shrub encroachment (also known as shrub invasion) is a typical pattern of patchy shrub distribution. This phenomenon represents one of the major environmental problems faced by the world's grassland ecosystems, and widely occurs in Africa, the USA, Australia, Asia, Europe. Shrub encroachment is exacerbated by overgrazing, prairie fires, climate change, atmospheric carbon dioxide concentrations, and biotic and abiotic environmental factors. Precipitation is considered to be an influential factor of patchy shrub distribution in arid and semiarid regions, which may regulate shrub coverage and the patch size. Theoretically, a steady shrub ecosystem would have a balance between water supply and demand, while certain levels of precipitation should correspond with optimal shrub coverage, i.e., the optimal ratio of shrub patch size to the total area in the steady patchy shrub ecosystem. The maximum entropy production principle (MEP) to living systems may be used to quantitatively explain optimization in the non-equilibrium process, such as the evolution of biological macromolecules and the plant optimization theory at different scales. Shrub patches and interspace may be considered as a water gathering system, in which entropy production is mainly caused by two processes: (1) soil water mixing in shrub patches and (2) either rainfall or outside runoff. For entropy flow, the main contributing factors are precipitation, evaporation, and transpiration processes in shrub patches. In this study, a modified model for a shrub patch ecosystem was built. It was applied to explain the moisture gathering in shrub patches. The model was based on the entropy change theory and Klausmier moisture gathering model, consisting of a pair of partial differential equations on soil moisture and plant biomass in the two-dimension plane. Based on the principle of maximum entropy production (MEP), the optimal ratio of shrub patches responding to different annual precipitation may be simulated. The simulation results showed a negative correlation between biomass and soil water content, which was verified by field experiments. When the ecosystem was steady (and the MEP principle was satisfied), the annual precipitation and the optimal ratio of shrub patches had a certain linear relationship. The existing data for Inner Mongolia was used to set up the model parameters, from which the optimal ratio of shrub patches responding to different annual precipitation values varying from 50 mm to 450 mm was simulated. Then the optimal shrub area ratios for practical application were determined in arid and semiarid regions. For example a region with 300 mm precipitation produced a shrub patch ratio of around 30%. This result provided theoretical support and practical guidance for ecosystem protection and plant recovery in arid and semiarid regions. The simulated values varied considerably, which may be because the effects of surface runoff were simplified to a linear relationship with precipitation in the modified model. Therefore, surface runoff sensitivity to patchy shrub distribution patterns requires further research.