Shrub encroachment represents a widespread shift in vegetation structure, yet its influence on the relationship between vegetation greening and ecosystem productivity across global aridity gradients remains poorly understood. Focusing on global grasslands as the study domain, we systematically examined shrub-encroached areas across diverse climatic zones spanning aridity gradients from arid to humid regions. Here, we integrated multi-source remote sensing datasets, including MODIS land cover, leaf area index (LAI), and gross primary productivity (GPP), with a global aridity index to systematically detect shrub encroachment. By combining trend analysis, Pettitt change-point detection, and a moving-window pairwise comparison approach, we characterized the spatiotemporal dynamics of encroachment and quantified its differential effects on vegetation dynamics. Our results showed that the global extent of shrub encroachment expanded continuously from 2002 to 2022, with an abrupt change detected between 2010 and 2014. Relative to comparable non-encroached coverage pixels (NSEC), comparable shrub-encroached coverage pixels (CSEC) exhibited generally increasing trends in LAI, fractional vegetation cover (FVC), and GPP, with enhancement magnitudes showing strong aridity gradient dependence, peaking in semi-arid and dry sub-humid regions, indicating that shrub encroachment is substantially regulated by water limitation. Furthermore, joint trend analysis revealed that the concurrent increase in both LAI and GPP (LAI+GPP+) was the dominant pattern, observed in approximately 70% of encroached areas, although decoupling persisted at high latitudes and in certain humid regions (LAI+GPP−). These findings demonstrate that shrub encroachment is fundamentally regulated by moisture gradients and exhibits pronounced spatial heterogeneity, providing new evidence for understanding carbon cycling dynamics and informing grassland ecosystem management under increasing aridity.
Dong et al. (Fri,) studied this question.