Spatial Variation in Responses of Plant Spring Phenology to Climate Warming in Grasslands of Inner Mongolia: Drivers and Application.

Plant spring phenology in grasslands distributed in the Northern Hemisphere is highly responsive to climate warming. The growth of plants is intricately influenced by not only air temperature but also precipitation and soil factors, both of which exhibit spatial variation. Given the critical impact of the plant growth season on the livelihood of husbandry communities in grasslands, it becomes imperative to comprehend regional-scale spatial variation in the response of plant spring phenology to climate warming and the effects of precipitation and soil factors on such variation. This understanding is beneficial for region-specific phenology predictions in husbandry communities. In this study, we analyzed the spatial pattern of the correlation coefficient between the start date of the plant growth season (SOS) and the average winter-spring air temperature (WST) of Inner Mongolia grassland from 2003 to 2019. Subsequently, we analyzed the importance of 13 precipitation and soil factors for the correlation between SOS and average WST using a random forest model and analyzed the interactive effect of the important factors on the SOS using linear mixing models (LMMs). Based on these, we established SOS models using data from pastoral areas within different types of grassland. The percentage of areas with a negative correlation between SOS and average WST in meadow and typical grasslands was higher than that in desert grasslands. Results from the random forest model highlighted the significance of snow cover days (SCD), soil organic carbon (SOC), and soil nitrogen content (SNC) as influential factors affecting the correlation between SOS and average WST. Meadow grasslands exhibited significantly higher levels of SCD, SOC, and SNC compared to typical and desert grasslands. The LMMs indicated that the interaction of grassland type and the average WST and SCD can effectively explain the variation in SOS. The multiple linear models that incorporated both average WST and SCD proved to be better than models utilizing WST or SCD alone in predicting SOS. These findings indicate that the spatial patterns of precipitation and soil factors are closely associated with the spatial variation in the response of SOS to climate warming in Inner Mongolia grassland. Moreover, the average WST and SCD, when considered jointly, can be used to predict plant spring phenology in husbandry communities.