The coupling coordination degree (CCD) of the production-living-ecological space (PLES) functional index is an indicator of regional sustainable development potential. However, previous studies have failed to reveal the driving mechanisms of the CCD of PLES functional index on a global scale. Therefore, this study employed the CCD model to evaluate the CCD of the PLES functional index and spatial regression models to measure the heterogeneous drivers using multi-source data in 2000, 2010, and 2020. The results demonstrated that ecological spaces dominated (85%) globally, while living spaces comprised the smallest share (3%). The ecological functional index was higher than the production and living functional indices. Further analyses revealed that topographic factors were the main restricting factors for PLES; the proportion of production and living spaces decreased with increasing altitude and slope, whereas ecological spaces showed the opposite trend. The global CCD values of the PLES functional index in 2000, 2010, and 2020 were 0.186, 0.189, and 0.198, respectively, showing an increasing trend. High CCD was generally observed in areas with dense population and industry where human activity systems interact with natural ecosystems. The formation of the CCD of the global PLES functional pattern results from the joint action of natural and socioeconomic factors, with pronounced spatial heterogeneity. Our findings can help optimize global territorial space utilization, improve territorial space utilization efficiency, and realize global sustainable development goals.

Urbanization is inevitably accompanied by drastic changes in regional land use and therefore presents an evident influence on ecosystems. Habitat quality (HQ) reflects the ability of a habitat to provide suitable conditions for the survival of an individual or population, and clarifying the interrelationships between urbanization level (UL) and HQ provides insights into sustainable urbanzation and ecosystems conservation. Much attention has focused on how urbanization related to HQ at multi-scales, but few studies have analyzed the interactive coercive relationships between UL and HQ on the global county level. To address this gap, we adopted the bivariate spatial autocorrelation and coupling coordination degree (CCD) model to determine their interactive coercive relationships at the county level globally. Results showed that the global average UL was 0.0807, 0.0838, and 0.0857 in 2000, 2010, and 2020, respectively, with a continuously increasing trend. The global average HQ was 0.6186, 0.6133, and 0.6111 for 2000, 2010, and 2020, reflecting opposite declining trends. The Moran’s I of population urbanization and HQ in 2000, 2010, and 2020 globally remained negative but showed an increasing trend, with values of −0.189, −0.228, and −0.254, respectively, while those of economic urbanization and HQ and land urbanization and HQ also remained negative and exhibited a similar increasing trend. The spatial autocorrelations for UL and HQ in different dimensions indicated that deserts and plateaus functioned as catchment areas with low UL and HQ. The CCD between UL and HQ in each county globally ranged from 0 to 0.70, while CCD between land urbanization and HQ was the most optimal at 0–0.90. The CCD between population urbanization and HQ was similar to that of economic urbanization in terms of spatial distribution. Nevertheless, the CCD of land urbanization in HQ exhibited the greatest degree of coordination among the three urbanization dimensions during the study period. These findings provided important support for sustainable urbanization and ecosystem protection globally.

Urban agglomerations (UAs), which play a significant role in socioeconomic development and urbanization, are confronted with mounting ecological stress and a profound global imbalance in ecosystem services (ES). Understanding the conflict and coordination of knowledge about ES supply and demand (ESSD) can facilitate more efficacious guidance on the ecological sustainable development of UAs. Nevertheless, the characteristics of the conflict and coordination relationship between ESSD in Chinese UAs remain unclear, and further investigation into the interactive coercive relationship between ESSD is warranted. Consequently, we employed spatial regression and coupled coordination models to elucidate the conflict and coordination relationship between ESSD, utilizing multi-source data on Chinese UAs from 2000 to 2020. We found that ES supply in the UAs decreased, while ES demand increased. Furthermore, the coupling coordination degree between ESSD demonstrated an increase trend. The overall coupling coordination degrees between ESSD in UAs were 0.260, 0.285, and 0.311 in 2000, 2010, and 2020, respectively. The central UAs were identified as stress areas, whereas the peripheral areas were classified as non-stress areas. This study offered valuable insights into the interactive relationship between ESSD in UAs and provided a basis for formulating differentiated policies for the sustainable development of ecosystems and human activities.