Spatial–Temporal Evolution and Driving Factors of Coupling Coordination between High-Quality Urban Development and Carbon Emissions Intensity in Guangdong Province

: Promoting urban green growth necessitates a dual focus on reducing carbon emissions and fostering high-quality development. However, the body of research on the coupling relationship between high-quality urban development and carbon emissions remains remarkably limited. Taking Guangdong province as an example, this study strived to establish a comprehensive evaluation sys-tem for high-quality urban development encompassing economic, societal, and ecological dimensions and further conducted an in-depth examination of the spatiotemporal pattern and driving forces of coupling coordination degree between high-quality urban development and carbon emission intensity during 2000–2017. The coupling coordination degree in Guangdong province has shown continuous growth, transitioning from moderate incoordination to moderate coordination. The coupling coordination degree showed the overall spatial distribution characteristics of “high in southeast and low in others”, with Zhuhai, Zhongshan, Foshan, Guangzhou, Dongguan, and Shenzhen as the core. Notably, technological advancement, environmental governance, and economic development emerge as pivotal factors that positively affect carbon emission intensity reduction, environmental quality improvement, and coupling coordination enhancement. This research provides valuable insights for achieving harmonized high-quality development in Guangdong province involving policies of regional differences, industrial competitiveness, and new-type urbanization.


Introduction
China's economic growth has shifted from a high-speed growth phase to a high-quality development phase.China needs to realize high-quality development by achieving carbon peak and carbon neutrality and adhering to the concept of green development.Cities are the spatial units and action carriers for energy conservation, emissions reduction, and ecological transformation [1].Their high-quality development contributes to economic transformation and upgrading, expanding consumption and investment, and achieving green development and sustainability, which holds profound significance for China's pursuit of high-quality development.As the largest emitter worldwide [2], China has the responsibility and capability of fulfilling the energy-conservation and emission-reduction task.Cities are the primary sources of carbon emissions and the critical regional units for achieving carbon reduction [3,4].Thus, promoting high-quality development while reducing carbon emissions is critical for China's government to achieve carbon goals and sustainable development.Guangdong province has been at the forefront of China's rapid economic growth.In 2022, Guangdong's GDP exceeded 12.9 trillion CNY, ranking first in China for 34 consecutive years.Meanwhile, the urbanization rate of permanent residents in Guangdong reached 71.3%, far higher than the national average (64.8%).In light of this, Guangdong had the largest volume of carbon emissions, with 266 million tons in 2022, but its carbon emission intensity was the third lowest nationwide.As a leading economic region in China, Guangdong performed well in economic growth and environmental performance, but it still faced various challenges to achieving sustainable and green development, such as pronounced urban-rural development disparities, lack of an effective mechanism of emissions reduction, and stubborn industrial structure.Therefore, it is imperative for Guangdong to incorporate the dual carbon strategy into high-quality urban development.The precondition is to explore the coupling relationship between urban development and carbon emission intensity over time and space.Moreover, the practice of Guangdong offers insights for other cities to achieve the balance of rapid economic growth and improved environmental efficiency.
As a concept of physical science, coupling refers to a phenomenon in which two or more than two systems influence each other via different interactions [5].The coupling degree measures the strength of interaction and the degree of correlation between the systems [6], but in some cases, it cannot explain the synergies between the systems [7].Hence, the concept of coupling coordination is proposed to reflect the coupling and synergies between the systems.Coordination, the process of organizing and synchronizing various elements or activities to work together effectively, can explain the sustainability of the system [8].The coupling coordination degree model is an adequate model used to measure both the coupling and coordination relationship among systems, which has been widely used in studies of high-quality development [9,10].Given the nature of these concepts, the analysis of the coupling coordination between high-quality urban development and carbon emissions would provide effective support for facilitating the coordinated development of the economy, society, and environment in a region.
Existing literature measured high-quality urban development based on single indicators.In particular, the scholars pay more attention to high-quality economic development [11,12].Yang et al. (2021) proposed high-quality economic development measurements regarding the ecological environment, economic efficiency, and economic structure [13].Guo et al. (2023) established an evaluation index system of high-quality economic development concerning the industrial structure, inclusive total factor productivity, technology innovation, ecological environment, and residents' living standards [14].Some studies focused on high-quality development concerning society and ecology [15,16].Artmann et al. (2019) proposed a systemic conceptual framework for sustainable cities by combining the ecological concepts of smart growth and green infrastructure [17].Wang et al. (2023) focused on a series of environmental indicators, such as average rainfall and PM2.5, to improve the framework of high-quality development [18].Moreover, some studies incorporated sustainability into urban high-quality development [19][20][21].For instance, Artmann et al. (2019) utilized multi-dimensional-ecological, economic, and social, as well as multi-scale-regional, city, and neighborhood-indicators to form a framework of urban sustainable development [17].Although the current analysis investigated the level of high-quality development based on a single system, it provides a reference for establishing a comprehensive evaluation index system of high-quality development.
Many scholars have studied the relationship between carbon emission and urban high-quality development with regard to a unilateral relationship [22][23][24] or a bilateral relationship [25][26][27].On the one hand, some literature analyzed the unidirectional relationship between carbon emission and specific aspects of high-quality development [28,29].Rapid industrial transformation requires massive energy use, resulting in a large-scale of carbon emissions, especially in India, Egypt, and other developing countries [30,31].In particular, Pan et al. (2020) took Stockholm County, Sweden, as a case study, analyzing the temporal change of greenhouse gas emissions under urban development [32].Wang et al. (2016) investigated the effect of urbanization on carbon emission in BRICS countries, esteeming that there is a long-run causality from urbanization to carbon emission [33].On the other hand, some scholars investigated the coupling coordination between carbon emission and specific aspects of high-quality development [34,35].Zhang et al. (2023) explored the coupling coordination degree between socio-economic development and carbon emissions of 76 contracting countries in The Paris Agreement, concluding that the coupling coordination degree at the global level is not symmetrical in spatial distribution.That is, European countries were generally rated higher than their African counterparts [36].Wang et al. (2015) applied a semi-parametric panel fixed-effects regression model to investigate the coupling between urbanization and carbon emissions in OECD countries, verifying the Environmental Kuznets Curve hypothesis [37].Huang et al. (2020) analyzed the temporal change and driving factors of the coupling coordination degree of urbanization and ecological environment in Kazakhstan [38].
Moreover, some scholars attempted to probe into the influencing factors of coupling coordination between urban development and carbon emissions [39][40][41].Yang et al. (2023) used coupled coordination models and panel regression models in 166 countries to identify the influencing factors of coupling coordination of sustainable development, concluding that carbon emissions have a significant spatial spill-over effect on the sustainable development coupling coordination degree [42].Taking Kazakhstan as an example, Huang et al. (2020) discovered that the GDP per capita and social fixed asset investment per capita, employment in industry, and services positively affected the coupling coordination degree of urbanization and ecological environment [38].With the IPAT/Kaya approach and variance analysis technology, Chontanawat (2018) investigated the driving factors of the relationship between increasing carbon emission and economic growth in the Association of Southeast Asian Nations during 1971-2013, concluding that population growth and per-capita income increase were the most important factors increasing carbon emission [43].However, only a few studies focused on the driving factors of the coupling coordination between urban development and carbon emission intensity [44,45].
This research contributes to the literature in the following aspects.First, this paper attempted to establish a comprehensive evaluation index system of high-quality development concerning innovation, coordination, greenness, openness, and sharing, whereas previous studies focused more on one or two aspects of high-quality development, especially the economy and ecology.Second, in contrast to previous literature that explored the coupling relationship between urban development and carbon emission, this research paid attention to carbon emission intensity, measuring the environmental performance and efficiency.Third, multiple aspects of high-quality urban development are addressed in order to investigate the driving forces of the coordinated development of carbon emission intensity and high-quality development from different perspectives.This study endeavors to provide some insightful implications for harmonized and effective high-quality development in Guangdong province and even in other regions.
The rest of this research is organized as follows: Section 2 explains the methodology and materials of the research; Section 3 reports the estimated results; Section 4 discusses the empirical results.Section 5 provides the policy suggestions and limitations.

Methods and Data
First, this paper calculated carbon emission intensity and constructed a comprehensive evaluation index system of high-quality urban development levels based on existing literature.Second, Hausman test was employed to select the key influencing factors of coupling coordination degree between high-quality urban development and carbon emission intensity.Third, the entropy weight method was utilized to calculate a comprehensive score of high-quality urban development, and coupling coordination degree model was used to analyze the coupling coordination degree between high-quality urban development and carbon emission intensity.Last, Moran I's test and Spatial Lag Model were exploited to analyze the spatial correlation of driving factors of coupling coordination degree.Figure 1 illustrates the conceptual framework of this study.The bold in Figure 1 illustrates the four steps in the methodology.

Accounting of Carbon Emission Intensity
Carbon emission intensity refers to carbon dioxide emissions per GDP value, indicating how clean a country or region's economy is.Carbon emission intensity enjoys some advantages over carbon dioxide emissions.On the one hand, with the characteristic of normalization, it is possible to compare the carbon emission efficiency among regions, industries, and activities of different scales and production levels, offering a fair basis for assessing carbon emissions.On the other hand, as a valuable concept in understanding how efficiently carbon emissions are associated with specific activities, carbon emission intensity plays a critical role in environmental assessments, climate change mitigation strategies, and making informed decisions to reduce the carbon footprint of various processes and sectors.Official statistics lack data on the carbon emissions of prefectural-level cities, let alone the carbon emission intensity of prefectural-level cities.Given data availability, this paper applied the carbon emissions from Carbon Emission Accounts and Datasets (CEADs) 1 , covering from 2000 to 2017.Moreover, the GDP data of prefectural-level cities stemmed from the China City Statistical Yearbook.In this paper, the GDP deflator method is used to make the calculation results of different years comparable; that is, the

. Accounting of Carbon Emission Intensity
Carbon emission intensity refers to carbon dioxide emissions per GDP value, indicating how clean a country or region's economy is.Carbon emission intensity enjoys some advantages over carbon dioxide emissions.On the one hand, with the characteristic of normalization, it is possible to compare the carbon emission efficiency among regions, industries, and activities of different scales and production levels, offering a fair basis for assessing carbon emissions.On the other hand, as a valuable concept in understanding how efficiently carbon emissions are associated with specific activities, carbon emission intensity plays a critical role in environmental assessments, climate change mitigation strategies, and making informed decisions to reduce the carbon footprint of various processes and sectors.Official statistics lack data on the carbon emissions of prefectural-level cities, let alone the carbon emission intensity of prefectural-level cities.Given data availability, this paper applied the carbon emissions from Carbon Emission Accounts and Datasets (CEADs) 1 , covering from 2000 to 2017.Moreover, the GDP data of prefectural-level cities stemmed from the China City Statistical Yearbook.In this paper, the GDP deflator method is used to make the calculation results of different years comparable; that is, the appropriate GDP index (CNY at 2000 constant prices) was used to adjust all of the monetary data.The calculation formula is as follows: where i and t represent the city and the year, respectively.Carbon emission intensity was assumed to be a negative indicator, normalized and standardized.
2.1.2.Evaluation Index of High-Quality Urban Development Level High-quality development is a highly intricate economic, social, and ecological system [46,47].Research on high-quality development that helps build a strong foundation for the definition and features of high-quality development has contemporary implications [48].From the perspective of "double carbon" goals, this paper referred to the Five Principles for Development proposed at the Fifth Plenary Session of the 18th Central Committee of the Communist Party of China 2 : innovation, coordination, greenness, openness, and sharing.Based on the studies of Gan et al. (2022) [49] and Bei (2018) [50], we established an evaluation index system for high-quality urban development levels (Table 1).The original data were from the China City Statistical Yearbook and Guangdong Statistics Yearbook.Some missing data were supplemented and improved using the interpolation method.* It refers to the comprehensive utilization of generated solid waste, rather than simple disposal or landfilling.This includes activities such as waste sorting, recycling, reuse, and energy recovery, aiming to minimize the negative impact on the environment.Sources: Adapted from Ref. [49].2022, Gan, W.; Yao, W.; Huang, S. and Ref. [50].2018, Bei J.

Hausman Test
A large and growing body of literature has shown that economic development, environmental governance, and technological innovation have contributed significantly to improving coupling coordination between carbon emissions and high-quality development [34,36,51].According to He and Liu (2022) [52] and Jiang et al. (2022) [40], we built a system of influencing factors of coupling coordination between carbon emissions and high-quality urban development, involving technology development, environmental governance, and economic development (Table 2).Technology development consists of human capital and technology expenditure intensity.Environmental governance incorporates a qualified amount of industrial wastewater and a percentage of municipal sewage disposal.Proportion of fixed asset investment, number of employed persons at the year-end, financial efficiency, and population density constituted economic development.The data originated from the China City Statistical Yearbook, the Guangdong Statistical Yearbook, and statistical bulletins of prefectural-level cities.This paper applied the entropy weight method to process the indicators.The entropy weight method is an objective weighting approach avoiding the biases induced by human subjective judgment [53].It determines the weight of each factor based on the variation in their values.The entropy weight method is widely accepted in environmental assessment, corporate performance, investment decisions, and energy use.One of the most outstanding advantages of this method is that it is assumed to be suitable for all decision-making processes requiring weight determination [54].The uncertainty of variables and the influences of controlling elements on outcomes can be measured, so this approach is highly acceptable in modeling and mapping natural hazards [55,56].Therefore, the entropy weight method was utilized to measure the high-quality urban development and carbon emission intensity in Guangdong province during 2000-2017.
(1) Data standardization.Suppose there are m evaluated objects, and each object has n indicators.Then, m × n judgment matrix can be constructed.We further standardized the indicators and obtained a new matrix R = r ij m×n .r ij ∈ [0, 1].The formula is as follows: If the indicator has a positive trend: If the indicator has a negative trend: (2) Calculation of entropy.We can calculate the contribution of the i year under the j index with the following equation: Next, the entropy value can be calculated by where Land 2023, 12, 2082 8 of 20 (3) Definition of the weight of entropy.The weight of the j indicator is determined with the following formula: (4) Calculation of indices.The comprehensive score in the year i can be calculated as

Coupling Coordination Degree Model
According to the definition of coupling coordination mentioned above, this paper constructed a coupling coordination degree model between carbon emission intensity and high-quality urban development.The formulas are as follows: where C represents the degree of coupling between carbon emission intensity and urban high-quality development; D represents the degree of coupling coordination, ranging from 0 to 1; and T reflects the overall scores of these two systems.In Guangdong province, carbon emission intensity is as equally important as high-quality development, so the values of α and β are equivalent, i.e., α = β = 0.5.
Based on the degree of coupling coordination (D), the coupling of carbon emission intensity and high-quality urban development was divided into eight stages, ranging from extreme incoordination to high-quality coordination (Table 3).The higher degree of coupling coordination indicates a strong relationship between high-quality urban development and carbon emission intensity, and high-quality development and carbon emission tend to make sound progress in coordination.When the D value is less than 0.5, it is coordinated at a lower level and called incoordination, indicating lower-quality urban development and worse carbon emission [46].

Moran's I Test
Since another main objective of this research is to analyze the driving factors of coupling coordination between carbon emission intensity and high-quality urban development, we decided to apply the panel models with more sample information.Before we conducted the panel analysis, we applied the spatial autocorrelation model to detect the correlation of observed values between neighboring locations within a geographic space.Moran's I is a correlation coefficient that measures the overall spatial autocorrelation of carbon emission Land 2023, 12, 2082 9 of 20 intensity and high-quality urban development in this study, which ranges from −1 to 1.The formula is as follows: In the above formula, x i and x j represent the data of city i and j, respectively.w ij presents spatial adjacency matrix.The results can be classified into three categories.First, when Moran's I coefficient exceeds 0, it suggests a positive spatial correlation among the subjects.A coefficient closer to 1 signifies a more pronounced positive correlation, showing increased spatial clustering of these subjects.Second, conversely, a Moran's I coefficient below 0 implies a negative spatial correlation.A coefficient nearing -1 reveals a heightened negative correlation, denoting a more prominent spatial disparity among the subjects.Third, a Moran's I coefficient close to 0 indicates that the subjects are distributed randomly, suggesting no spatial autocorrelation [58].

Spatial Lag Model
The spatial correlation should be considered in delving into the driving factors of the coupling coordination relationship between carbon emission intensity and high-quality development.This research employed the spatial econometric model for examining the spatial effects of these driving factors in Guangdong province.The Spatial Lag Model (SLM) is a regression model incorporating spatial dependence into standard regression models like linear regression.It is commonly used in spatial econometrics to analyze geospatial data where observations are influenced by spatial proximity to other observations.The SLM enjoys several advantages.On the one hand, the SLM helps reduce model bias.Ignoring spatial dependence can lead to biased and inefficient estimates.SLM helps in reducing or eliminating this bias.Ignoring spatial dependence can lead to biased and inefficient estimates.On the other hand, the SLM is adaptable and can be combined with other modeling methods, such as incorporating spatially lagged independent variables or integrating time dynamics for spatiotemporal analyses.According to the results of the LM test, Wald test, and LR test, the Spatial Lag Model is more suitable for exploring the spatial effects of high-quality urban development than the Spatial Errors Model.The formula is as follows: where ρ denotes coefficient of endogenous interaction effects, reflecting the degree of spatial spillover or diffusion effects between adjacent cities; α stands for linear correlation coefficients; ε indicates random errors.

Study Area
Situated on the southern coast of mainland China, Guangdong province is a pioneer of reform and opening up (Figure 2).From 1978-2021, Guangdong's economy has grown more than 400 times and replaced Jiangsu province as China's most prominent economic province.In 2022, the gross domestic production (GDP) of Guangdong province amounted to 12.91 trillion CNY, with a year-on-year growth of 1.9%.As a significant economic province nationwide, Guangdong province accounted for about 6.7% of the national energy consumption and 5% of national carbon emissions, supporting about 9% of the country's permanent resident population and 10.7% of the national economy.There are 21 prefecture-level cities in Guangdong province, which is divided into Pearl River Delta (PRD), East Guangdong, West Guangdong, and North Guangdong.In general, the economic development of Guangdong follows the principle of core-periphery, with the PRD as the most developed region.In 2022, with a population of 78.29 million (61.6% of the province), the PRD region produced 10.47 trillion CNY (81% of the province).There is a significant regional disparity in the province.In 2020, Guangdong province's urbanization rate reached 74.2%, entering a mature period of urbanization development, and

Trends in Land Urbanization
As shown in Figure 3, the carbon emission intensity in Guangdong province decreased from 1.93 CNY/ton to 0.71 CNY/ton by 63.2% during 2000-2017.The carbon emission intensity score increased from 10.4 to 17.4 during the study period, with a growth of 67.7%.According to Equation (1), carbon emission intensity is decided by the volume of carbon emission and GDP.From 2000 to 2017, carbon emissions experienced a double increase in Guangdong province, while GDP grew sixfold.In other words, Guangdong province has been likely to witness an improvement in energy efficiency and economic structure during the study period.Meanwhile, similar to previous studies that investigated the trend of high-quality urban development in Guangdong [59,60], there was a marked increase in high-quality urban development from 6.3 to 10.7, indicating that Guangdong tended to possess a more vital ability to achieve higher-quality development.
More importantly, there was a significant difference in environmental quality and high-quality development at different stages of economic development (Figure 3).There was also a striking variation in the coupling coordination degree between carbon emission intensity and high-quality development.The results in Figure 3 revealed that the coupling coordination degree experienced stable and slow growth, ranging from 3.5 to 5.2 during the study period.During the initial period (2000)(2001)(2002)(2003)(2004)(2005)(2006), the coupling coordination degree of Guangdong province stayed at 3.6.from 2007 to 2017; however, there has been a significant rise of approximately 30% in the coupling coordination degree.This result demonstrates that there has been a strong integration between carbon emission intensity and high-quality urban development in Guangdong province driven by rapid economic growth and improving environmental quality, which is also supported by Chontanawat (2018) [43].

Trends in Land Urbanization
As shown in Figure 3, the carbon emission intensity in Guangdong province decreased 1.93 CNY/ton to 0.71 CNY/ton by 63.2% during 2000-2017.The carbon emission intensity score increased from 10.4 to 17.4 during the study period, with a growth of 67.7%.According to Equation (1), carbon emission intensity is decided by the volume of carbon emission and GDP.From 2000 to 2017, carbon emissions experienced a double increase in Guangdong province, while GDP grew sixfold.In other words, Guangdong province has been likely to witness an improvement in energy efficiency and economic structure during the study period.Meanwhile, similar to previous studies that investigated the trend of high-quality urban development in Guangdong [59,60], there was a marked increase in high-quality urban development from 6.3 to 10.7, indicating that Guangdong tended to possess a more vital ability to achieve higher-quality development.
More importantly, there was a significant difference in environmental quality and high-quality development at different stages of economic development (Figure 3).There was also a striking variation in the coupling coordination degree between carbon emission intensity and high-quality development.The results in Figure 3 revealed that the coupling coordination degree experienced stable and slow growth, ranging from 3.5 to 5.2 during the study period.During the initial period (2000)(2001)(2002)(2003)(2004)(2005)(2006), the coupling coordination degree of Guangdong province stayed at 3.6 from 2007 to 2017; however, there has been a significant rise of approximately 30% in the coupling coordination degree.This result demonstrates that there has been a strong integration between carbon emission intensity and highquality urban development in Guangdong province driven by rapid economic growth and improving environmental quality, which is also supported by Chontanawat (2018) [43].Compared to a coupling degree, a coupling coordination degree is more conducive to accurately and reasonably assessing the coordination between environmental quality and high-quality urban development [61].Figure 4 illustrates the coupling coordination degree of 21 prefectural-level cities in Guangdong province.During the study period, the mean coupling coordination degree fluctuated between 0.4, experiencing a stable growth from 0.3 to 0.5.Moreover, before 2007, Guangdong province was in the stage of mild disorder, while it has entered the stage of mild incoordination since 2008.In 2015, Guangdong province was in the stage of mild coordination, which might be due to energy conservation and environmental protection measures that were reinforced with the enhancement of high-quality urban development.This research investigated a similar trend of the coupling coordination degree between urban development and carbon emission in other studies [43,62].Compared to a coupling degree, a coupling coordination degree is more conducive to accurately and reasonably assessing the coordination between environmental quality and high-quality urban development [61].Figure 4 illustrates the coupling coordination degree of 21 prefectural-level cities in Guangdong province.During the study period, the mean coupling coordination degree fluctuated between 0.4, experiencing a stable growth from 0.3 to 0.5.Moreover, before 2007, Guangdong province was in the stage of mild disorder, while it has entered the stage of mild incoordination since 2008.In 2015, Guangdong province was in the stage of mild coordination, which might be due to energy conservation and environmental protection measures that were reinforced with the enhancement of high-quality urban development.This research investigated a similar trend of the coupling coordination degree between urban development and carbon emission in other studies [43,62].Compared to a coupling degree, a coupling coordination degree is more conducive to accurately and reasonably assessing the coordination between environmental quality and high-quality urban development [61].Figure 4 illustrates the coupling coordination degree of 21 prefectural-level cities in Guangdong province.During the study period, the mean coupling coordination degree fluctuated between 0.4, experiencing a stable growth from 0.3 to 0.5.Moreover, before 2007, Guangdong province was in the stage of mild disorder, while it has entered the stage of mild incoordination since 2008.In 2015, Guangdong province was in the stage of mild coordination, which might be due to energy conservation and environmental protection measures that were reinforced with the enhancement of high-quality urban development.This research investigated a similar trend of the coupling coordination degree between urban development and carbon emission in other studies [43,62].Likewise, the coupling coordination degree witnessed slow growth at the city level during the study period (Figure 5).More importantly, there has been a decrease in the gap in the coupling coordination degree between prefectural-level cities over time.For instance, the difference between the highest and lowest level of coupling coordination degree was 0.40 in 2000, whereas this difference was only 0.36 in 2017.During 2000-2017, Qingyuan, Zhongshan, Huizhou, and Shaoguan experienced the fastest growth of the coupling coordination degree, with rates of 208%, 133%, 76%, and 65%, demonstrating that these cities performed well in coordination development between ecological systems and high-quality development.This result might be related to the terrain of these cities.For instance, the terrain of Shaoguan is mainly mountainous and hilly, and the plain and platform area accounts for about 20%.More importantly, Shaoguan has made great progress in environmental protection, such as pollution control and environmental governance.Likewise, the coupling coordination degree witnessed slow growth at the city level during the study period (Figure 5).More importantly, there has been a decrease in the gap in the coupling coordination degree between prefectural-level cities over time.For instance, the difference between the highest and lowest level of coupling coordination degree was 0.40 in 2000, whereas this difference was only 0.36 in 2017.During 2000-2017, Qingyuan, Zhongshan, Huizhou, and Shaoguan experienced the fastest growth of the coupling coordination degree, with rates of 208%, 133%, 76%, and 65%, demonstrating that these cities performed well in coordination development between ecological systems and high-quality development.This result might be related to the terrain of these cities.For instance, the terrain of Shaoguan is mainly mountainous and hilly, and the plain and platform area accounts for about 20%.More importantly, Shaoguan has made great progress in environmental protection, such as pollution control and environmental governance.Only Guangzhou, Zhuhai, and Foshan were in the stage of barely coordination, and Shenzhen was in the stage of mild coordination.Most cities were in a low degree of coupling coordination development at the initial stages.In 2004, the number of cities in the period of serious incoordination declined, and Shenzhen and Guangzhou began to enter the stage of moderate coordination and mild coordination.When it came to 2008, most cities were in the period of moderate incoordination, whereas Meizhou was still in the stage of serious incoordination.What stands out in Figure 6 is that Qingyuan improved its level of coupling coordination degree, switching from extreme incoordination to moderate incoordination during 2000-2008.This can be mainly attributed to the strong awareness of environmental protection.

Spatial Patterns of Coupling Coordination Degree
Since entering 2010, Guangdong province made significant progress in improving the coupling coordination degree.Except for the cities in North Guangdong, other cities  Only Guangzhou, Zhuhai, and Foshan were in the stage of barely coordination, and Shenzhen was in the stage of mild coordination.Most cities were in a low degree of coupling coordination development at the initial stages.In 2004, the number of cities in the period of serious incoordination declined, and Shenzhen and Guangzhou began to enter the stage of moderate coordination and mild coordination.When it came to 2008, most cities were in the period of moderate incoordination, whereas Meizhou was still in the stage of serious incoordination.What stands out in Figure 6 is that Qingyuan improved its level of coupling coordination degree, switching from extreme incoordination to moderate incoordination during 2000-2008.This can be mainly attributed to the strong awareness of environmental protection.

Spatial Patterns of Coupling Coordination Degree
respectively.Apparently, Shenzhen and Guangzhou were the pillar cities of the coupling coordination degree of high-quality regional development.In 2017, spatially, the coupling coordination degree showed Zhuhai, Zhongshan, Foshan, Guangzhou, Dongguan, and Shenzhen as the core; Jiangmen, Zhaoqing, and Huizhou as the sub-core; and other cities as the edge.Obviously, the PRD region leads in the high-quality urban development in the province.The spatial agglomeration characteristics of coupling coordination degree between carbon emission intensity and high-quality development were measured using the global Moran's I index of spatial autocorrelation in Guangdong province.As shown in Table 4, Moran's I index of coupling coordination degree almost passed the significance level test, except for 2001-2004.In other words, no solid spatial autocorrelation existed between carbon emission intensity and high-quality development during 2001-2004.A possible explanation for this might be the small spatial scale, constraining the significance of Moran's I index.Moreover, during 2004-2017, the value of Moran's I index was positive, demonstrating that the coupling coordination degree of carbon emission intensity and high-quality economic development exhibits strong positive spatial autocorrelation, implying a polarization that higher-quality development usually accompanies lower emission intensity.There were high-value coupling coordination agglomeration and low-value agglomeration.Since entering 2010, Guangdong province made significant progress in improving the coupling coordination degree.Except for the cities in North Guangdong, other cities entered the stage of mild incoordination in 2012.More importantly, Zhongshan and Foshan were in a period of primary coordination.Shenzhen and Guangzhou improved the coupling coordination degree level, entering the moderate and benign coordination stages, respectively.Apparently, Shenzhen and Guangzhou were the pillar cities of the coupling coordination degree of high-quality regional development.In 2017, spatially, the coupling coordination degree showed Zhuhai, Zhongshan, Foshan, Guangzhou, Dongguan, and Shenzhen as the core; Jiangmen, Zhaoqing, and Huizhou as the sub-core; and other cities as the edge.Obviously, the PRD region leads in the high-quality urban development in the province.
The spatial agglomeration characteristics of coupling coordination degree between carbon emission intensity and high-quality development were measured using the global Moran's I index of spatial autocorrelation in Guangdong province.As shown in Table 4, Moran's I index of coupling coordination degree almost passed the significance level test, except for 2001-2004.In other words, no solid spatial autocorrelation existed between carbon emission intensity and high-quality development during 2001-2004.A possible explanation for this might be the small spatial scale, constraining the significance of Moran's I index.Moreover, during 2004-2017, the value of Moran's I index was positive, demonstrating that the coupling coordination degree of carbon emission intensity and high-quality economic development exhibits strong positive spatial autocorrelation, implying a polarization that higher-quality development usually accompanies lower emission intensity.There were high-value coupling coordination agglomeration and low-value agglomeration.
From 2004 to 2017, the values of Global Moran's I varied from 0.137 to 0.476 and were significant at the 90% confidence level.This indicated that the coupling coordination degree in the study area showed a stronger spatial correlation during the study period.Moreover, Global Moran's I experienced a slight decline during 2007-2008 and 2014-2015, demonstrating a steady decrease in the coupling coordination degree cluster.The former decrease was partly because of the financial crisis of 2008, reducing carbon emissions to some extent, while differentiated regional environmental protection entry policies influenced the latter decrease in Guangdong province since 2013 [63].

Analysis of Influencing Factors
As shown in Table 5, the fixed effect is more suitable for estimating the model because of the statistical significance of the P-value of the Hausman test.The value of the Hausman test was negative (−66.2) because the progressive assumption of the random-effects model cannot be satisfied.According to the results of the LM test and LR test, the Spatial Lag Model is more suitable for analyzing the driving factors of coupling coordination degree between carbon emission intensity and urban high-quality development.In other ways, a city's coupling coordination degree is relevant to other cities.Overall, technology development, environmental governance, and economic development significantly and positively affected the coupling coordination degree in Guangdong province.First, the enhancement of education levels was conducive to accelerating the cultivation of high-quality talents and improving the foundation of scientific research, which is in accordance with Wang et al. (2017) [64] and Zhang et al. (2021) [65].The increased number of individuals receiving higher education further facilitated technological innovation and innovative capabilities.This finding is likely to be related to the importance given to education in Guangdong.To some extent, the increase in technology expenditure indicated a promotion of technological input and enterprise reform, thereby accelerating high-quality regional development.
Second, environmental governance had a significant and positive effect on the coupling coordination degree.Strict environmental control and continuous improvement in pollutant treatment can reduce carbon emission intensity, improve environmental quality, and enhance the coupling coordination degree.Existing studies have confirmed that provinces in eastern China have conformed to the "Porter Hypothesis" [66,67], suggesting that appropriate environmental regulations also promote technological innovation and high-quality urban development.
Third, economic development had a positive correlation with the coupling coordination degree.Improvements in the employment environment, enhancement of financial efficiency, and effective control of population density made an outstanding contribution to laying a solid foundation for high-quality urban development, realizing urban transformation from "quantitative change" to "qualitative change".These results are consistent with Shi et al. (2020) [68], who found that economic development tended to have positive relationships with coupling coordination degrees in Eastern China.Based on the value of coefficients, economic development played the most significant role in the promotion of the coupling coordination degree, followed by environmental governance during the study period.

Discussion
An accurate grasp of the coupling and coordination relationship between carbon emission intensity and high-quality urban development will help accelerate the process of Guangdong province achieving its double carbon goals (carbon peaking and carbon neutrality).It involves the coordination and development of economic, social, and ecological systems.However, existing research findings primarily focus on the unidirectional and coupling relationship between carbon emissions and specific indicators or dimensions of high-quality urban development.It is worth mentioning that a city is a complex system formed by the interaction, influence, and mutual constraints of various elements such as society, economy, and ecology.Therefore, combined with a new development philosophy, this research attempted to establish an evaluation index system of high-quality urban development based on innovation, coordination, greenness, openness, and sharing.Moreover, we analyzed the coupling coordination degree between carbon emission intensity and high-quality development in 21 prefectural-level cities in Guangdong province during 2000-2017 with the entropy weight method, coupling coordination degree model, spatial autocorrelation method, and Spatial Lag Model.This study further analyzed the spatiotemporal characteristics and the driving factors of the coupling coordination degree.
Consistent with previous studies concerning carbon emission intensity [69], this study concluded that the carbon emission intensity witnessed a substantial decrease of 63.2% in Guangdong province during 2000-2017.In general, Guangdong takes the lead in highquality urban development and carbon emission reduction attributed to technological innovation, industrial transformation, and government incentives [69].Therefore, Guangdong made progress in harmony between high-quality development and carbon emission intensity, with a slow and stable increase in coupling coordination degree from 2000 to 2017.However, there was a significant regional disparity in the province.Spatially, the coupling coordination degrees in the PRD region were far higher than that in other regions, and the lowest coupling coordination degree was in North Guangdong, which is supported by previous studies concerning the divergence of the coupling coordination degree in Guangdong [70,71].This was mainly attributed to an imbalance of regional development concerning the economy, society, and ecology in the province [72].In other words, future policy makers should consider regional differences instead of adopting a "one size fits all" approach.Moreover, during the study period, the coupling coordination degree of carbon emission intensity and high-quality urban development exhibited strong positive spatial autocorrelation in Guangdong province, demonstrating that cities are placing more emphasis on environmental protection and sustainable development through enacting more effective policy instruments.This result coincides with current studies that show that balancing rapid economic growth and environmental sustainability is a priority for many regions to achieve inclusive development [72].Such a trend aligns with the global requirements for low-carbon, green, and sustainable development.It also indicates that cities have recognized the critical balance between economic development and environmental protection and are progressively achieving it.
According to the results of influencing factors, technology development, environmental governance, and economic development had a significant positive impact on the coupling coordination degree in Guangdong.This finding was also supported by Fengyin et al. [73], who stated that these factors were the key drivers of promoting high-quality development.These factors offer innovation elements and material bases to support a higher coupling coordination degree between high-quality urban development and carbon emission intensity.More importantly, compared with other factors, economic development had the highest coefficient, positively affecting the promotion of the coupling coordination degree.Over the study period, the economic growth reached over 10% annually, and the urbanization growth has been above the national level in Guangdong province.This result implies that economic development is a priority for achieving higher-quality development.
Furthermore, previous studies prove that strict environmental control is conducive to improving environmental performance, but the impact on the growth of productivity will become weaker with stricter environmental control [74].In other words, stricter environmental governance will help improve the coupling coordination degree between high-quality urban development and carbon emission intensity.This is in accordance with China's goal of upgrade and reconstruction as well as industrial structure adjustment.Recently, Guangdong has pioneered clean production, recycling transformation, and a new energy industry, promoting comprehensive green transformation of economic and social development.

Conclusions
This study aimed to develop an evaluation index system for high-quality urban development focusing on innovation, coordination, green practices, openness, and sharing.Taking Guangdong as an example, we examined the relationship between carbon emission intensity and high-quality development across 21 prefectural-level cities from 2000 to 2017 using the entropy weight method, the coupling coordination degree model, the spatial autocorrelation method, and the Spatial Lag Model.The main results are as follows.
First, Guangdong province witnessed a gradual and consistent rise in its coupling coordination degree, transitioning from a phase of moderate incoordination to moderate coordination during the study period.Second, regarding spatial distribution, Zhuhai, Zhongshan, Foshan, Guangzhou, Dongguan, and Shenzhen formed the primary core, with Jiangmen, Zhaoqing, and Huizhou serving as secondary hubs and the remaining cities acting as peripheral areas.Third, in Guangdong province, the degree of coupling coordination between carbon emission intensity and high-quality urban development shows a pronounced positive spatial correlation.Moreover, technology development, environmental governance, and economic growth positively influenced the coupling coordination degree.
Based on the above findings, the following policy recommendations are provided.First, it is imperative to regulate differentiated economic development strategies at the city level.This study confirmed that there was significant regional disparity of coupling coordination degree in the province.Therefore, differentiated economic transformation strategies are suitable for optimal industrial transformation.Second, it is necessary to reinforce inter-city collaborations concerning technological innovation, skill training, and investment projects to uplift the sub-core cities like Jiangmen, Zhaoqing, and Huizhou.This research discovered that peripheral cities have the potential to achieve harmonized growth of high-quality urban development by fostering connections between them and the core cities.Third, a more mature policy system involving economic development, environmental governance, and technology development should be comprehensively enhanced to give full play to these factors.The results proved that these three factors played a significant positive role in improving the coupling coordination degree between high-quality urban development and carbon emission intensity.In other words, the integration of technology development, environmental governance, and economic development helps Guangdong achieve the goal of higher-quality development.
In summary, Guangdong's high-quality urban development process has significant impacts not only on China's economic, social, and environmental development but also serves as a model for global sustainable growth.China, especially Guangdong, provides valuable lessons and best practices in various aspects, including urban planning, infrastructure development, land utilization, and resource management.Nevertheless, there are still some limitations to this study.On the one hand, this study is limited by the absence of the latest data, failing to have a broader view of the relationship between carbon emission intensity and high-quality development in Guangdong province.Therefore, future research could investigate this relationship by collecting the latest data.On the other hand, this research only considered several second-level indicators for each influencing factor due to data accessibility.In contrast, many factors concerning society, economy, and ecology could influence the coupling coordination degree.These questions are left for future studies.

Land 2023 ,
12,  x FOR PEER REVIEW 4 of 20 development and carbon emission intensity.Last, Moran I's test and Spatial Lag Model were exploited to analyze the spatial correlation of driving factors of coupling coordination degree.Figure1illustrates the conceptual framework of this study.The bold in Figure1illustrates the four steps in the methodology.

Figure 3 .
Figure 3. Evaluation score in Guangdong province from 2000 to 2019.Overall score refers to the score of high-quality development level, calculated based on Equation (9).

Figure 3 .
Figure 3. Evaluation score in Guangdong province from 2000 to 2019.Overall score refers to the score of high-quality development level, calculated based on Equation (9).

Land 2023 , 20 Figure 3 .
Figure 3. Evaluation score in Guangdong province from 2000 to 2019.Overall score refers to the score of high-quality development level, calculated based on Equation (9).

Figure 4 .
Figure 4. Boxplot of coupling coordination degree in Guangdong province in 2000-2017.Figure 4. Boxplot of coupling coordination degree in Guangdong province in 2000-2017.

Figure 4 .
Figure 4. Boxplot of coupling coordination degree in Guangdong province in 2000-2017.Figure 4. Boxplot of coupling coordination degree in Guangdong province in 2000-2017.

Figure 6
Figure6illustrates the spatial characteristics of coupling coordination degree in2000,  2004, 2008, 2012, and 2017.In 2000, West Guangdong was in the stage of moderate incoordination, while East Guangdong was primarily in the stage of serious incoordination.Only Guangzhou, Zhuhai, and Foshan were in the stage of barely coordination, and Shenzhen was in the stage of mild coordination.Most cities were in a low degree of coupling coordination development at the initial stages.In 2004, the number of cities in the period of serious incoordination declined, and Shenzhen and Guangzhou began to enter the stage of moderate coordination and mild coordination.When it came to 2008, most cities were in the period of moderate incoordination, whereas Meizhou was still in the stage of serious incoordination.What stands out in Figure6is that Qingyuan improved its level of coupling coordination degree, switching from extreme incoordination to moderate incoordination during 2000-2008.This can be mainly attributed to the strong awareness of environmental protection.Since entering 2010, Guangdong province made significant progress in improving the coupling coordination degree.Except for the cities in North Guangdong, other cities

Figure 6
Figure 6 illustrates the spatial characteristics of coupling coordination degree in 2000, 2004, 2008, 2012, and 2017.In 2000, West Guangdong was in the stage of moderate incoordination, while East Guangdong was primarily in the stage of serious incoordination.Only Guangzhou, Zhuhai, and Foshan were in the stage of barely coordination, and Shenzhen was in the stage of mild coordination.Most cities were in a low degree of coupling coordination development at the initial stages.In 2004, the number of cities in the period of serious incoordination declined, and Shenzhen and Guangzhou began to enter the stage of moderate coordination and mild coordination.When it came to 2008, most cities were in the period of moderate incoordination, whereas Meizhou was still in the stage of serious incoordination.What stands out in Figure 6 is that Qingyuan improved its level of coupling coordination degree, switching from extreme incoordination to moderate incoordination during 2000-2008.This can be mainly attributed to the strong awareness of environmental protection.

Figure 6 .
Figure 6.Spatial evolution of coupling coordination degree between carbon emission intensity and high-quality urban development in 2000, 2004, 2008, 2012, and 2017.

Figure 6 .
Figure 6.Spatial evolution of coupling coordination degree between carbon emission intensity and high-quality urban development in 2000, 2004, 2008, 2012, and 2017.

Table 1 .
Evaluation index system of high-quality urban development level.

Table 2 .
Influencing factors of coupling coordination between carbon emissions and high-quality urban development.

Table 4 .
Spatial autocorrelation coefficient of coupling coordination degree between carbon emission intensity and high-quality development.

Table 4 .
Spatial autocorrelation coefficient of coupling coordination degree between carbon emission intensity and high-quality development.

Table 5 .
Drive mechanism measurement estimation results.