Distribution and Change Characteristics of Ecosystem Services in Highly Urbanized Areas along Gradients of Human Activity Intensity: A Case Study of Shenzhen City, China

Comprehensively and objectively revealing the spatial relationship between ecosystem services (ESs) and human activity intensity (HAI) is vital for achieving sustainable development goals. However, existing studies still have an incomplete understanding of how ESs change with HAI gradients. Therefore, taking Shenzhen City, China, as an example region experiencing rapid urbanization, the distribution of ESs and HAI in 2010 and 2020 were quantified using the InVEST model and the human footprint index method; at the same time, the gradient perspective was introduced and the spatial and temporal correlation characteristics of ESs along 10 HAI gradient bands, from weak to strong, were captured by applying multiscale geographically weighted regression (MGWR) and the bivariate spatial autocorrelation model. The findings showed that (1) the HAI demonstrated an increasing trend (20.63 (2010) and 23.36 (2020)), and the area with high values of HAI (the 10th gradient band) was distributed in the western part of the study area; meanwhile, the area with low HAI values (the 1st gradient band) was more distributed in the eastern part of the study area. (2) On the whole, the average levels of water conservation, soil conservation, carbon storage, and habitat quality decreased from 2010 to 2020; the spatial distribution characteristics of these parameters were similar. (3) In general, ESs were negatively correlated with HAI, and the negative correlation ratio was more than 65%. At the same time, the spatial and temporal correlations between ES patterns and HAI under different gradient bands were significant. These findings can effectively alleviate the pressure on the ecosystem caused by human activities, which is of great significance for the sustainable development of highly urbanized regions.