Multiscale effects and drivers of landscape heterogeneity for water‐related ecosystem services in urban agglomerations

Abstract Understanding the multiscale impacts and drivers of urban agglomeration landscape patterns for ecosystem services (ESs), especially water‐related ecosystem services (WESs), is essential for the development of regional ecological management. However, the multiscale impacts and driving mechanisms of urban agglomeration landscape patterns for ESs have not been adequately explained. In this study, multivariate data were employed, and the InVEST model, trend test method, coupled GeoDetector and geographically and temporally weighted regression (GTWR) method were utilized to comprehensively explore the spatial and temporal changes in landscape patterns and WESs in the Pearl River Delta urban agglomeration (PRDUA) at various grid and administrative scales from 1990 to 2020 and to determine the driving mechanisms affecting WESs. The results indicated that the variation characteristics of landscape patterns and WESs in the PRDUA were consistent, forming a binary spatial structure of core and peripheral areas in an inverted “U” shape around the estuary of the Pearl River. The relationship between landscape patterns and WESs weakened with the increase of scale, and the correlation coefficient decreased by approximately 0.10 from 5 km to 10 km grid scale. Additionally, precipitation (PRE) was the main factor controlling WESs changes in the PRDUA, explaining more than 50% of the changes in WESs, and the regression coefficients ranged from 0.0825 to 0.1584. Changes in WESs were the result of the combined effects of natural factors, including PRE, landscape pattern, elevation, slope, and socioeconomic factors, such as population and gross domestic product (GDP). Overall, these findings could contribute to optimizing regional landscape patterns and fostering sustainable development of the ecological environment in urban agglomerations.