Abstract

In theory, the industrial agglomeration is a double-edged sword as there are both positive and negative externalities. China’s cities, with great disparities on degrees of the industrial agglomeration, often face different energy and carbon dioxide emission problems, which raise the question whether the industrial agglomeration promotes or inhibits energy efficiency and carbon dioxide emission. This paper explored the effects of the industrial agglomeration on carbon efficiency in China. Spatial econometric methods were implemented using panel data (2007–2016) of 285 cities above the prefecture level. The results revealed that industrial agglomerations have significant impacts on the urban carbon efficiency with significant spatial spillover effects. The agglomerations of the manufacturing and high-end productive service industries take positive effects on carbon efficiency while the low-end productive and living service industries take negative effects. As a comparison, we found that the agglomeration effects at the level of the megalopolis are greater than those at the national level, especially for the living services industry, in which the higher levels of agglomeration make the effects on carbon efficiency change from negative to positive. The divisions of labor for the central and common cities in the megalopolises are integrated into the industrial agglomeration. Furthermore, the fractional-order grey forecasting model is used in this paper. By the virtue of its advantage in dealing with small sample data which lack statistical rules, this paper makes an out-of-sample prediction of carbon efficiency and industrial agglomeration degree of Chinese cities. By adding the predicted results to the spatial correlation test, new evidence on the spatial correlation of carbon efficiency and spatial division of labor between cities is obtained. Based on the empirical results of the present study, we have proposed some policy recommendations.

1. Introduction

Since its reform and opening up, China’s rapid urbanization process has brought about a significant increase in energy consumption and CO2 emissions. During 1980–2014, China’s average annual growth rate of CO2 emissions was as much as 9.2%. The International Energy Agency (IEA) reported that China has, since 2007, far exceeded the CO2 emissions of the European Union and the United States [1]. At present, according to statistics from the World Bank, China’s CO2 emissions had reached 10.29 billion tons, accounting for 28% of global CO2 emissions in 2014. In comparison, in 2014, the CO2 emissions of the European Union and the United States were only 3.24 and 5.25 billion tons, respectively [2]. There is no doubt that the Chinese government faces an increasing amount of international and domestic pressure while keeping economic development stable. The Chinese government has made a list of commitments for reducing carbon emissions and to curb global climate change. In 2014, the Chinese government announced, in the “Sino US joint statement on climate change,” that it plans to reach the peak of CO2 emissions by 2030. At the 2015 Paris Climate Conference, the Chinese government committed to CO2 emissions per unit of GDP being 60–65% lower in 2030 than the level in 2005. Accordingly, China has enhanced a series of restrictions on CO2 emissions. For instance, China has controlled the expansion of the production capacity of the manufacturing industry because this industry produces high levels of pollution, consumes a large amount of energy, and is responsible for high CO2 emissions. China has intentionally reduced its reliance on the manufacturing industry to achieve industrial upgrading. However, the implementation of restrictions on the manufacturing industry is bound to slow down economic growth.

In the process of economic development, cities are regarded as carriers of economic activities and also account for about 85% of China’s CO2 emissions [3]. The rapid urbanization increases not only massive CO2 emissions but also the level of total factor productivity (TFP) with the accumulation of the population, economic factors, and activities in the cities [4]. An improvement of the TFP plays an essential role in promoting economic growth and transformation [5, 6], as well as improving efforts to reduce energy consumption and CO2 emissions [7, 8]. Therefore, the economies of scale and technology spillover effects will inevitably affect the quantity and efficiency of CO2 emission. In this sense, the effective use of urban industrial agglomeration effects is conducive to not only reducing CO2 emissions but also promoting economic growth.

Many scholars studied how urbanization improves carbon efficiency from the perspective of agglomeration economies. Most scholars identified that the driving factors for improving carbon efficiency are the aggregation of the population, technological progress, and industrial structure. For example, Li et al. studied the efficiency of CO2 emissions of China’s 31 manufacturing industries, from 2012 to 2016, and found that the high-tech industry ranked at the top of the list of industries with different levels of technology [9]. Wang et al. discovered that the energy structure, technological level, and rate of urbanization all have significant impacts on carbon efficiency [10]. Sun et al. evaluated the carbon efficiency and the global value chain (GVC) position index. They highlighted that, compared with the labor-intensive and resource-intensive manufacturing industries, the technology-intensive manufacturing industry in the GVC is most effective in improving carbon efficiency [11]. Zhang et al. indicated that enterprise size is a crucial factor for improving carbon efficiency [12]. Wang and Ma determined that the level of urbanization, the structure of energy consumption, and the level of industrialization are the main driving factors of CO2 emissions [13]. Wang et al. found that the gradual improvement of the level of technology is the main driving force for the improvement of the efficiency of CO2 emissions [14, 15]. Wang et al. found that the production technology shows heterogeneities at the level of region [16] and sector [17].

In general, technological progress, economies of scale, and the aggregation of the population will significantly affect the carbon efficiency. However, these studies did not discuss the agglomeration of the segmented service industries. Thus, it is unknown that how the impact of the spatial agglomeration of economic activity on carbon efficiency is derived from the agglomeration of services. Besides, according to the China Megalopolises Development Report 2010, megalopolises cover 21.13% of the national land area, accounting for 48.99% of the total population, 51.41% of the urban population, and 46.7% of the urban members. In addition, they create 81.94% of the added value of the second industry and 83.5% of the added value of the third industry. The megalopolis is an advanced form of spatial organization with a concentration of economic actions resulting from highly developed industrialization and urbanization [18]. The spatial characteristics of a megalopolis contribute to providing transportation, communication, information, and other infrastructures, as well as market conditions, which further stimulate the agglomeration of industries natively. However, few scholars have studied the relationship between industrial agglomeration and carbon efficiency in megalopolises.

Compared to the extant research, the present study makes three major contributions: (1) the present study built a spatial Durbin model to study the effects of the heterogeneity of the industrial agglomeration on carbon efficiency at the national level and at the level of the megalopolis; (2) we explored the different impacts from various service industry agglomerations on urban carbon efficiency and identified specific impacts in megalopolises; (3) comparisons between the manufacturing and segmented service industries were made, and finally, we sought the optimal development path of cities for the dual goals of economic development and the reduction of CO2 emissions.

2. Theoretical Frameworks and Research Hypotheses

At present, most urban-related studies focused on the relationship between different types of industrial agglomeration and the labor productivity of a single city, but no consistent conclusions could be drawn. First, some studies showed that industrial agglomeration has caused serious environmental pollution along with the expansion of enterprise-scale and large-scale use of energy from fossil fuels. Virkanen found that industrial agglomeration in southern Finland caused heavy metal pollution of the local water resources and atmosphere [19]. Verhoef and Nijkamp found that industrial agglomeration harms the quality of the environment, leading to the aggravation of environmental pollution [20]. Frank et al. showed that industrial agglomeration worsens air quality in sampled cities in the European Union [21]. Duc et al. surveyed industrial companies on both sides of the Vietnamese river and found that industrial agglomeration caused river pollution [22]. In contrast to mentioned above, some scholars further discovered that industrial agglomeration is conducive to the progress and diffusion of environmental protection technology, promoting production and energy efficiencies, and reducing environmental pollution [2327]. Han et al. analyzed the impact of Marshall’s agglomeration, such as labor concentration effects, economies of scale of intermediate inputs, and space spillover effects, on industrial energy efficiency. They pointed out that the space availability of intermediate inputs and space technology spillover effects significantly increase industrial energy efficiency [28]. Furthermore, Han et al. constructed a spatial Durbin model to explore the effects of industrial specialization and diversified agglomerations of manufacturing on CO2 emissions. They provided an evidence that both of them can reduce CO2 emissions locally and in surrounding cities through the agglomeration of economic externalities [29]. Zhang et al. analyzed data from 18 cities in Henan Province and proposed that industrial agglomeration reduces CO2 emissions in the surrounding areas and that there is a negative impact on the intensity of CO2 emissions [30]. Shen et al. used a dynamic spatial panel model and demonstrated that the agglomeration of the productive service industries has a significant positive space spillover effect on energy efficiency in surrounding areas [31]. Jacobs et al. found that diversified agglomeration is available for the high-end productive service industry in central cities, such as large and international cities. They found that the industries of these cities have the characteristics of high added value within products, a large service radius, and a small transaction frequency, thus affecting the surrounding areas with significant spatial spillover effects of improvements in energy efficiency [32]. However, in the existing literature, there is an absence of research on the low-end productive services and living services. In addition, there are no comparative studies between the manufacturing and service industries, or with the inter-heterogeneity in the subdivided service industries, regarding the effects of industrial agglomeration on urban carbon efficiency. Therefore, we propose Hypothesis 1.

Hypothesis 1. The agglomeration effects of various industries show heterogeneity, and their impacts on urban carbon efficiency have diversity spatially.
The urban area provides transportation, communication, information, and other infrastructures as well as market conditions and stimulates the agglomeration of industries. The development of industrial agglomeration is part of the trend of urbanization, resulting in the industries being close to the source of innovation of the central city to obtain advanced technology. When the urbanization is in different stages, the internal industrial agglomeration pattern has different characteristics. The expansion of cities promotes the diffusion and re-agglomeration of industries. The second industry cluster is dominant in the developing stage of urbanization, while the third industry cluster is evident in the developed stage. The coordinated development of regional industrial agglomeration and urbanization promotes the overall development of the regional economy. Therefore, cities of different scales and grades are adjacent to each other or cities of the same scale and grade in space. This will form a megalopolis through strengthening the division of labor and cooperation between cities. Huang et al. argued that the development of a megalopolis is beneficial to the improvement of urban ecological efficiency and that the degree of improvement is affected by the city level [33]. Yue et al. found that the impact of an industrial structural adjustment on the energy intensity of cities in the Pearl River Delta is not apparent [34]. Chen et al. found that the transfer of high energy-consuming industries, which are led by capital-intensive industries, from the developed region to the developing region, can effectively reduce CO2 emissions in the developed regions [35]. Li et al. indicated that industrial transfer between any two cities within a megalopolis would bring about a decline in energy intensity [36]. Xiao et al. revealed that the performance of the total-factor CO2 emissions of the tertiary industry is better than that of the secondary industry. However, prospects for the low-carbon development of the tertiary industry, in the Yangtze River Delta, are not optimistic [37]. Therefore, we propose Hypothesis 2.

Hypothesis 2. The diversified characteristics of the industrial agglomeration of a megalopolis in space result in its effects on carbon efficiency being different.
The structural optimization effect is the driving force for improving carbon efficiency in a megalopolis. First, the development of megalopolises improves the proportion of production service industries and promotes the diversification of production in central cities. Second, the development of megalopolises promotes their industrialization progress and the specialization of production in common cities. Henderson pointed out that there is an inverted U-shaped relationship between urban scale and industrial agglomeration. With the expansion of the urban scale, the competition for urban public resources, the rise of factor costs, the congestion of transportation, and other factors will accelerate the industrial diffusion. Simultaneously, some industries leapfrog to the same level or to a low gradient city, namely, Friedman’s “centre-periphery” mode, which leads to the initial formation of megalopolises [38, 39]. In recent years, megalopolises have played an essential role in China’s new urbanization process [40]. In the Yangtze River Delta, the development of the service industry in Shanghai far exceeds other industries, but the manufacturing industry still leads and plays a decisive role in Jiangsu and Zhejiang [41]. The Pearl River Delta has entered the post-industrialization period and has developed with the promotion of the agglomeration of the service industry and manufacturing diffusion [35]. In the Circum-Bohai Sea Economic Zone, the heavy industry accounts for too large a share of the overall economy, which restricts the division of labor between the cities [42]. Yu et al. found that there is a negative correlation between the size of a city and carbon emissions and further demonstrated that a megalopolis is more efficient with regard to carbon emissions [43]. Liu et al. deemed that the social, land, population, and economic urbanization are all related to urbanization and conducted an assessment of the multiple effects of these components of urbanization on CO2 emissions in the Pearl River Delta [4446]. Zheng et al. predicted that the energy consumption of the Yangtze River Delta would continue to rise, and the main source of energy consumption will be manufacturing. They also suggested that the government should update its methods and concepts related to sustainable development and limit the development of energy-intensive industries [41]. Ouyang et al. found that GDP per capita, environmental regulation strength, factor input structure, foreign direct investment, local government spending, and the degree of openness have positive impacts on the industrial total-factor energy efficiency in the Pearl River Delta [47]. Li et al. argued that socioeconomic factors, such as the industrial structure, economic level, carbon intensity, and other spatial factors, play significant roles in decreasing carbon intensity in the Yangtze River Delta [48]. Therefore, the development of megalopolises improves the proportion of the tertiary industry in central cities and also promotes the progress of both industrialization and specialized production in common cities. The development of megalopolises decentralizes the diseconomy of agglomeration of central cities and promotes the flow of the elements of production to common cities. However, although the aforementioned studies have contributed to the theoretical discussion on the mechanism by which industrial agglomeration influences economic development of megalopolis, they still lack empirical support. On the basis of these studies, we propose the following.

Hypothesis 3. On the internal side of a megalopolis, the division of labor between central and common cities is coupled with industrial agglomeration for improving carbon efficiency.

3. Data and Methodology

Limited by the availability of data, the present study selected panel data (2007–2016) for 285 cities in China. The original data sources were the China City Statistical Yearbook, China Energy Statistics Yearbook, and Beijing Shu Huitong Environmental Technology Research Institute Data Service Network (http://www.3edata.com). To ensure the consistency and comparability of the statistical data, the present study used the relatively stable statistical caliber of the city.

3.1. Measurement of CO2 Emissions at Prefecture-Level Cities

Given the availability of data sources, we selected four forms of energy consumption, electricity, natural gas, liquefied petroleum gas, and transportation industries, to measure CO2 emissions [49].

To better reflect the actual situation in China, the reference coefficient of coal for various energy conversions and the average low calorific value were adopted from the Chinese National Standard GB/T2589-2008. The data of the unit calorific value of the carbon content and the oxidation rate used to calculate the CO2 emission coefficient were obtained from the “2006 IPCC National Greenhouse Gas Emission Inventory Guide.” The formula used to calculate the CO2 emissions is as follows:where denotes the CO2 emissions, is the i-th energy, is the j-th city (region), denotes the CO2 emissions generated by the i-th energy consumption in city j, and is the consumption of the i-th energy in city j, and represents the CO2 emission coefficient of the i-th energy in city j.

The average low calorific value, carbon content per unit calorific value, carbon oxidation rates of the principal energies, and calculated CO2 emission factors are shown in Table 1.

Data on the electricity consumption of the whole society were collected from the China City Statistical Yearbook. Glaeser and Kahn argued that there is an independent CO2 emission coefficient in each regional power grid [50]. Accordingly, we divide China into six regions: North China, Northeast China, East China, Central China, Northwest China, and South China. The CO2 emission coefficient for electricity was measured according to the baseline emission coefficient in the regional grid, which was issued by the China Certified Emission Reduction Exchange Information Platform.

The data on the supplies of natural gas and liquefied petroleum gas in a region were obtained from the data service network of the Beijing Shu Huitong Environmental Technology Research Institute. The CO2 emission coefficient of natural gas is 0.2666 kg/m3, and the CO2 emission coefficient of liquefied petroleum gas is 3.1013 kg/m3.

Due to differences between the regions of China and incomplete statistical data, it was difficult to obtain data such as mileage travelled by vehicle, the fuel used per mile for different types of motor vehicles, etc. According to public data, it is less accurate to use the “bottom-up” method to calculate CO2 emissions from transportation. In the present study, we use a “top-down” algorithm to calculate the CO2 emissions by transportation in various cities in China. The CO2 emissions from transportation mainly come from direct emissions from the combustion of fossil fuels. In contrast, the indirect emissions of CO2 from electricity and heat consumption are less. Therefore, in the present study, the calculations of CO2 emissions by transportation included only the direct emissions from the consumption of fossil fuels.

At the provincial level, the Regional Energy Balance Sheet contains the energy consumed by transportation, warehousing, and postal industries. The present study excluded the energy consumption of the warehousing and postal industries. This was achieved by using the income ratios of transportation, warehousing, and postal industries in the “China Economic Census Yearbook” to estimate energy consumption and calculate CO2 emissions for the transportation industry only. Then, referring to Li et al. [51], we obtained the values of the passenger and freight turnovers from the China Economic Network Statistics Database and the energy intensity data in the Transportation Energy Data Book (TEDB). These were used to calculate the ratio of CO2 emissions between passenger and freight transport. According to this ratio, we divided the CO2 emissions by the transport activities for passenger and freight at the provincial level. At the city level, we calculated the proportion of the volume of the passenger transport of each city and that of the freight transport of each city in a province. We then divided the provincial CO2 emissions into the city level in terms of each calculated proportion. Finally, we aggregated the CO2 emissions data of the passenger and freight transports as the total emissions for the transportation industry of the prefecture-level cities.

According to Figure 1, from 2007 to 2016, the total CO2 emissions of China’s prefecture-level cities increased from 2022.77 to 3499.70 million tons, with an average annual growth rate of 6.31%. The CO2 emissions of cities in megalopolises increased from 1065.89 to 1741.97 million tons, with an average annual growth rate of 5.63%. The CO2 emissions of other cities increased from 956.89 to 1757.72 million tons, and the average growth rate was 7.04%. There has been an annual decrease in the proportion of CO2 emissions in megalopolises. During 2007–2016, the 70 cities in three megalopolises contributed an average of 50.37% of the total CO2 emissions among 285 above prefecture-level cities across the country. This means that the megalopolises are still the main sources of China’s CO2 emissions.

3.2. Spatial Distribution of CO2 Emissions

According to Tobler’s first law, everything is related to other things, and similar things are more closely related to each other. This means that spatial data can interact with each other due to geographical location. We explored the spatial distribution of CO2 emissions in 2007, 2011, and 2016 in above prefecture-level cities. As shown in Figure 2, CO2 emissions show distinct spatial agglomeration characteristics. CO2 emissions decreased gradually along the east, central, and west areas. The Pearl River Delta, Yangtze River Delta, Beijing-Tianjin-Hebei Economic Zone, and Chengdu-Chongqing Economic Zone, which are the largest and most developed areas in China, showed higher levels of emissions. Shandong and Liaoning Provinces showed significantly high concentrations of emissions, which combined with the Tianjin-Hebei Economic Zone form the Circum-Bohai Sea Economic Zone.

3.3. Measurement of the Agglomeration of the Service and Manufacturing

As for the measurement of industrial agglomeration, a series of measurement methods are available, such as market concentration, location entropy, Herfindahl index, spatial Gini coefficient, geographic concentration index, and MS index. Each of them has its own focus. The location entropy measures industrial agglomeration through industrial specialization, which is conducive to indicating the advantageous industries in the region. The Herfindahl index reflects the concentration of enterprises in the industry, and the spatial Gini coefficient reflects the regional agglomeration of the industry. The geographic concentration index comprehensively evaluates the degree and structure of industrial agglomeration. Accordingly, to find the dominant industries and leading industries in different cities, the present research concentrated on evaluating the impacts of different industrial agglomerations, in different regions, on the local and surrounding carbon efficiencies. The number of employees for each industry was obtained from public data. Therefore, we used the location entropy to measure the level of urban industrial agglomeration. The formula for calculating the industrial agglomeration level iswhere indicates the number of employees in industry i in city j, indicates the total number of employees in city j, indicates the total number of employees in industry i in the country, and indicates the total number of employees in the country. The larger the value of , the more obvious the agglomeration advantage. When , it indicates that the development intensity of the industry in the region is higher than the national average for similar industries; when , it indicates that the development intensity of the industry in the region has a significant comparative advantage.

In the manufacturing sector, due to the lack of complete classification statistics for the above prefecture-level city, the preset study analyzed it as a whole and could not discuss the segmented manufacturing industries any further. In 2011, China established classification standards for the 15 subindustries of the service sector. Complete employment statistics were available at the city level, and the preset study categorized them accordingly. In the classification of the service industries, referring to Li and Li [52] and Lu [53], we divided the service sector into three categories: (1) the productive service industries whose service targets are the manufacturers of industrial goods; (2) the living service industries who serve to the final consumer; and (3) the public service industries whose providers are governments. The public administrations determine the size of the public service industry, and it was not discussed in the present study. In terms of the intensity of knowledge accumulation in the productive service industry, there exist high-end productive service industries with strong knowledge accumulation and low-end productive service industries with low knowledge accumulation [54]. Accordingly, this paper focused on the specific industries within manufacturing, productive services (high-end and low-end), and living services (see Table 2).

Based on the original data of employees in 285 above prefecture-level cities in China, the location entropy calculation formula was used to measure the average annual value of agglomeration level of each industry in various cities (Figure 3). Figure 4 shows the average location entropy in the different areas.

As shown in Figure 3, manufacturing showed a significant aggregation in space, especially in the eastern coastal areas. However, there were no apparent spatial agglomerations in the distributions of the service industries. In Figure 4, the average degree of the manufacturing agglomeration of each megalopolis was higher than that at the national level. The levels of manufacturing agglomeration decreased in the order of the Pearl River Delta, Yangtze River Delta, and Circum-Bohai Sea Economic Zone, respectively. The average concentrations of all kinds of service industries in the Circum-Bohai Sea Economic Zone were higher than those in the Pearl River Delta and the Yangtze River Delta. The concentration of the productive service industry in the Pearl River Delta was the lowest of the three megalopolises.

3.4. Fractional-Order Grey Forecasting Model

At present, there are two popular methods to make forecasting: (1) the white box models, commonly used as the differential equations, and (2) the black box models, such as the machine learning models. The grey models take the merits of above models. Especially, the grey models are good at modeling time series forecasting in small sample [55]. The grey model is between the black model and the white model. It contains both unknown information and known information. At the same time, the factors in the system have uncertain relations, and the biggest characteristic is that the accumulation process can reflect the regularity of the data. This model has been applied well in many fields since its establishment [56]. The old grey models all are integer derivative models, which belong to ideal memory models and are not suitable for describing some irregular phenomena. In fact, the grey model is usually fractional. For example, the weather is the derivative of fractional order of climate, so it is the bridge between cause and effect in the complex system. Moreover, when the objects with fractional characteristics are described by fractional order, the essential characteristics and behaviors of the objects can be better revealed [57]. Therefore, this method is capable to describe dynamic behavior of the processes and systems with high accuracy and concise models. The economic phenomena are the activities of multiagent social system, so fractional-order grey forecasting is capable for addressing our topic.

Referring to Wu [55], this paper used the fractional-order grey forecasting model to forecast the carbon efficiency and industrial agglomeration degree of each prefecture-level city from 2017 to 2021. The modeling steps are as follows.

Step 1. The cumulative sequence of order was calculated:

Step 2. A non-negative sequence, , should be an accumulation operator of order . Stipulate , , and it is called a cumulative sequence of order .

Step 3. Substituting into Step 2, the least square method is used to estimate parameter .

Step 4. Using to predict .

Step 5. Decrease by order . That is,

4. Econometric Models and Main Results

4.1. Basic Model

Referring to Griliches [58], we assumed that the Cobb–Douglas function can manifest the production function of city i at period t:where is the desirable output, , , and are the labor, capital, and energy consumption for urban production, respectively, denotes the total factor productivity, which is related to economies of economy, technology change, technology spillover, and institutional arrangements, , , and are the production elasticities, and is the error term. Divide both sides of equation (5) by C (CO2 emission):where is the desirable output per unit of CO2 emissions, denoting the carbon efficiency. , , and are input factor-carbon ratios including labor, capital, and energy consumption, describing the efficiencies of input factors. Therefore, the efficiencies of input factors are the decompositions of carbon efficiency. Take the logarithm form of equation (6) as follows:where . Also, we take government intervention (GOV) [59] and foreign direct investment (FDI) [60] as control variables in equation (7). The controlled model iswhere and are the elasticities of the control variables. Transform equation (8) into a matrix form:

In equation (9), , , , , , and are the vectors of the logarithms of independent variables and is the vector of error terms.

Moreover, the improvement of total factor productivity brought about by technological progress is another important way to improve carbon efficiency, other than improving the efficiency of input factors. The technology spillover and economies of scale, based on the theory of agglomeration economies, are two main driving forces promoting the total factor productivity [61, 62], which also will further strengthen carbon efficiency. Moreover, the effects of agglomeration not only affect the carbon efficiency in the local area but also affect surrounding cities [61, 63]. Therefore, we composed total factor productivity with the industrial agglomerations of native and surrounding cities. Accordingly, can be set aswhere is the exogenous variable noting an advanced technology for all cities national-wide; , , , and are manufacturing, low-end productive service, high-end productive service, and living service agglomerations for city i at time t, respectively; ,, , and indicate the elasticities of the four kinds of industrial agglomeration; and , and are geometrically weighted values of the agglomerations of manufacturing, low-end productive service, high-end productive service, living service, and the total factor productivity of neighboring cities. The degree of effects of industrial agglomeration among cities is described by identical , , , , and , respectively. However, the net effects of these effects are dominated by the connectivity between each city, which can be described with the exogenous friction terms , where and . The closer the connectivity exists, the higher the value takes. We took the logarithm of equation (10) and transformed it into a matrix form:where , , , and are the vectors of the logarithms of industrial agglomeration and is the matrix of the friction term . If , , , , and and if , , and are not eigenvalues of , equation (11) can be transformed to

Substitute equation (12) into equation (9) and multiply equation (9) by on both sides:where , , , , , and . Describe equation (13) with the subscript i and obtain the carbon efficiency of city i:

Equation (14), which includes both local effects and spatial effects, is called a spatial Durbin model. The descriptive statistics for the above variables are shown in Table 3. The original data sources were the China City Statistical Yearbook, China Energy Statistics Yearbook, and Beijing Shu Huitong Environmental Technology Research Institute Data Service Network (http://www.3edata.com).

4.2. Spatial Weighted Matrix and Spatial Correlation Analysis

Based on the availability of data sources, we constructed and standardized the weighted spatial matrix with a geographical distance matrix, where the sum of elements of each row is equal to one:where is the distance between each pair of cities based on latitudinal and longitudinal data, where . The geographical attenuation parameter takes the value of 2. On the basis of , the present study used Moran’s I index to measure the autocorrelation in the geographical distribution of the explained variables and the core explanatory variables:where is the value of the variable for region i. Moran’s I index takes a value within [-1, 1]. A positive value indicates that the distribution of the spatial agglomeration of the variable of interest is in a specific region with similar values; otherwise, it will follow the discrete distribution. If Moran’s I value is close to zero, observed values are randomly distributed and have no relationship in space. As shown in Table 4, Moran’s I values of carbon efficiency, manufacturing agglomeration, and high-end productive service agglomeration all meet the significance test at the 1% level, indicating that these variables have significant spatial dependencies. The low-end productive service passed the significance test for most years, but a small part of the living service industrial agglomeration degree passed the significance test. Therefore, the results showed that it is necessary to consider the spatial effects of industrial agglomeration on carbon efficiency.

In order to make up the lag of the released data in the official statistics and also to make the division of functional space between cities more convincing, in this part, we use the data from 2007 to 2016 and the fractional-order grey prediction model to forecast the data of the urban carbon efficiency and industrial agglomeration degree from 2017 to 2021 and measure the Moran value according to the prediction results (i.e., all the predictions for Moran’s I index were calculated again based on the predicted data). The results show that the spatial correlation of carbon efficiency is expected to stabilize at a high level during 2017–2021 after experiencing the “high-low-high-low” fluctuation. For the manufacturing and high-end productive service, which are of representative significance in the spatial division of urban labor, their Moran’s I value continued to rise from 2007 to 2014. After reaching the maximum value in 2015, both of them showed a decreasing trend in the following year. According to the results of fractional-order grey prediction, their spatial correlation decreased year by year from 2017 to 2021, which also well confirmed the changing trend. It is worth noting that the spatial correlation of the living service industry in each city gradually emerged, and its spatial correlation showed a stable upward trend from 2017 to 2021. The results give us the following implications. (1) In terms of the spatial correlation of carbon efficiency, except for carbon dioxide due to its natural characteristics, the economic correlation between cities is increasingly close. The carbon efficiency presents a “high-high” agglomeration trend in space, and the spatial correlation degree expressed by Moran’s I values remains at a high level in the predicted years. (2) The essence of the spatial division of labor in urban agglomerations discussed above is that cities of different levels focus on the development of different industries. According to the Moran value of the predicted values, it can be found that the pattern of urban functional spatial division of labor still exists at the national level. Especially for manufacturing and high-end productive service, local government should choose industries with comparative advantages according to their development status and resource endowment, to form “staggered peak” development among cities. China’s overall economic development will continue to advance in this development direction. (3) As a living service industry whose main function is to meet a number of necessary needs in residents’ life, the development of this industry is highly correlated with urban economic conditions and residents’ wealth. The increasing degree of spatial correlation of this industry year by year also proves that economic development among cities is more closely correlated.

4.3. Determination of Spatial Econometric Model

In order to accurately reflect the degree of spatial effects and the reasons for the spatial dependence, it is necessary to select a suitable spatial econometric model before doing the regression analysis. Referring to Elhorst [64], the method of “specific to general” was used in the preset study to identify the spatial econometric model (Table 5). First, the Lagrange multiplier (LM) tests were used to choose the spatial autoregressive model (SAR) or the spatial error model (SEM). If both LM_lag and R-LM_lag fail the significance test, we choose the SAR. If both LM_err and R_LM_err fail the significance test, we choose the SEM. In the present study, we rejected the SEM because both LM_err and R_LM_err passed the significance tests. However, we could not choose the SAR on the condition that the LM_lag passed the significance test and the R-LM_lag failed the significance test. We need to determine whether the spatial Durbin model (SDM) could degenerate into SAR and SEM through performing the LR test and Wald test. Because both of these results significantly rejected the null hypothesis, we chose the SDM for the spatial econometric regression. Then, we performed a Hausman test to determine whether the model uses fixed or random effects. However, when calculating the variance inflation factor (VIF), we found that the addition of urban fixed effects would produce multicollinearity problem (VIF>10) among the explanatory variables (Table 6). Since the result of the Hausman test significantly rejected the null hypothesis, we used the time fixed spatial Durbin model.

4.4. Estimated Results of Spatial Econometric Model

According to Table 5, we chose the SDM with the time fixed effect to estimate equation (14). To take the robustness test, we also estimated OLS, SAR, SEM, and SLX with the time fixed effect (Table 7). Most explanatory variables passed the significance test, which indicated that the empirical data support the theoretical model and variables. Also, urban carbon efficiency is affected by both the endogenous and exogenous spatial interaction effects.

The spatial econometric model includes the interaction terms between the spatial weight matrix and the explained variables. It also includes the interaction terms between the spatial weight matrix and the explanatory variables. Consequently, the total effect includes two parts: the direct effect of the local area and the indirect effect caused by other areas (also known as spatial spillover effect). As the spatial spillover effect, the change of explanatory variables will cause the change of the carbon efficiency both in the local area and surrounding areas. Therefore, the change of carbon efficiency in surrounding areas will eventually cause the change of carbon efficiency in the local area owing to the feedback effect [65]. In other words, the change of the explanatory variables of city i at time t will not only have a direct impact on the explained variable of the city itself but also have an indirect impact on the explained variable of other cities (i.e., the adjacent areas of the city) and will eventually affect the area in turn (because for other areas, the area is also the “other area”). Therefore, considering the global effect in the spatial Durbin model, the coefficients of the explanatory variables only represent the direction and significance of their influences on the explained variable in the SDM model but do not reflect their actual marginal effects, which was raised as an issue by LeSage and Pace [66]. Accordingly, to solve the problem of the coefficients in the spatial econometric model being difficult to explain, we further separated the direct and indirect effects (Table 8). The direct effects contain the pure direct effects and the feedback effects on local carbon efficiency. The indirect effects contain the spatial impact of the explanatory variables from the surrounding areas on the local carbon efficiency, reflecting the spatial spillover effect.

As shown in Table 8, we found that the agglomeration effects of various industries showed significant heterogeneity. Thus, Hypothesis 1 is validated. The agglomeration of manufacturing industry and the agglomeration of high-end productive services have significant impacts on the urban carbon efficiency. In addition, the impact of the high-end productive service is greater than that of the manufacturing industry, while the low-end productive services and living services have negative impacts on the urban carbon efficiency. Li et al. indicated that the greater demand for energy from fossil fuels of the manufacturing industry results in more CO2 emissions than other industries [67]. However, in the process of manufacturing industry agglomeration, the knowledge spillover, economies of scale of intermediate inputs, the information exchange among enterprises, and the multiple infrastructure construction will further promote carbon efficiency [29, 30, 68]. The impacts of segmented service industries showed significant heterogeneity. In the high-end productive service industry, the learning effect of knowledge spillover among enterprises, with the high added value of products, has brought significant improvements in efficiency. On the contrary, the agglomeration of the low-end productive service industries, such as transportation, warehousing and postal industry, and leasing and business service industry, cause an increase in direct energy consumption and “crowding effect” as bad as urban congestion and resource waste to worsen the urban carbon efficiency. Therefore, in order to achieve low-carbon development, priority should be given to the high-end productive service industry with the characteristics of knowledge- and capital-intensiveness.

The impacts of industrial agglomerations on carbon efficiency from neighboring areas, that is, the spatial spillover effect brought by urban industrial agglomeration, showed that all industries will have positive impacts on carbon efficiency. The impact of the high-end productive service industry is not significant, which shows that China is still in the early stage of industrial transformation. The high-end productive service industry with a large service radius and small transaction frequency has not shown significant spatial spillover effects. The impact of the living service industry is not significant either, which indicated that the living service industry has the characteristic of a native service.

As for other explanatory variables, the direct effect coefficients of labor-carbon ratio, capital-carbon ratio, and energy-carbon ratio are significantly positive, noting that the improvement of efficiencies of input factors has significant positive effects on the improvement of carbon efficiency. The direct and indirect effects of FDI are both significantly positive, which shows that FDI plays a positive role in promoting China’s economic development and shows that the “pollution heaven hypothesis” is not currently available for China [67]. The coefficient of lnGOV is significantly negative, which indicates that the government should not interfere too much in the economy. Otherwise, it will lead to an inefficiency of resource allocation and the “race to the bottom” in energy utilization, which will not be conducive to improving urban carbon efficiency.

4.5. The Effects of Industrial Agglomeration on Carbon Efficiency in a Megalopolis

At the level of the megalopolis, the agglomerations of manufacturing, high-end productive service, and living service industries all play important roles in promoting the urban carbon efficiency (Table 9), while the low-end productive service industry mainly has a negative impact on the carbon efficiency. Compared with the estimated results at the national level (Table 8), the manufacturing industry and the high-end industry have the same impact direction, but the degree of industrial agglomeration in promoting carbon efficiency in the megalopolis is significantly higher than that at the national level. The reason is that the spatial concentration of industries and population in a megalopolis improve the frequencies of knowledge spillover between manufacturing and service enterprises and reduce the transportation cost of goods exchanged between enterprises. They ultimately improve the efficiency of the economic operation and the efficiency of the energy utilization of the whole region. As far as the control variables are concerned, the direct effects of labor-carbon ratio, capital-carbon ratio, and energy-carbon ratio are significantly positive, while the direct effect of lnGOV is negative. This indicates that at both the national or level of the megalopolis, too much policy intervention will reduce carbon efficiency.

Comparing three large megalopolises, the direct effect of manufacturing agglomeration in the Yangtze River Delta is the largest, and the positive spatial spillover effect of manufacturing agglomeration in the Pearl River Delta is significantly greater than that in other megalopolises. In terms of the segmented service industry, the low-end productive service industry has no significant impact in the Circum-Bohai Sea Economic Zone. The direct effect in the Yangtze River Delta city group is significantly negative, and the indirect effect in the Pearl River Delta city group is significantly negative. The high-end productive service industry has a significant positive impact in each megalopolis, and the impact degrees are arranged by the Pearl River Delta, the Yangtze River Delta, and the Circum-Bohai Sea Economic Zone from large to small. All of the above are consistent with Hypothesis 2. It is worth noting that the spatial effect of high-end productive service agglomeration in the Pearl River Delta megalopolis is far greater than other megalopolises, which provides a good example for the development of megalopolis.

We know that the manufacturing industry is inclined to be located near to the cities with productive services [69], while there also exists a “crowding out” effect between the manufacturing and service industries [70]. We further analyzed the trend of the ratio of the employment number in central cities (municipalities, subprovincial cities, and provincial capital cities) in megalopolises to measure the level of the division of labor. As shown in Figure 5, the employment proportion of high-end productive service increased in the central cities but that of manufacturing and low-end productive service decreased year by year, which partially confirm Hypothesis 3. Therefore, the central cities tend to cultivate the high-end productive service industry and gradually become “technology pools” and “market areas” for improving carbon efficiency of themselves and surrounding areas. The common cities are inclined to develop their manufacturing industry with complementary functions for neighboring areas. Therefore, the divisions of labor for central and common cities in the megalopolis are coupled with the industrial agglomeration for improving carbon efficiency.

5. Conclusion

There are two constraints on China’s economic development: economic growth and the control of CO2 emissions. In the present study, we used the spatial econometric method to explore the impact of industrial agglomeration on carbon efficiency with panel data of 285 cities above prefecture level and three megalopolises in China, for the period 2007–2016. The results showed the following. (1) At the national level, all industrial agglomerations, except the living service industry, take positive roles in improving the urban carbon efficiency, and the impact from low-end productive service agglomeration mainly comes from the spatial spillover effect. (2) Comparing the industrial agglomerations of the whole country and the megalopolises, we found that manufacturing industry, high-end productive service industry, and living service industry of each megalopolis have significant influences on the urban carbon efficiency, and the influences are greater than the average national level. (3) In the three megalopolises, manufacturing agglomeration, high-end productive service agglomeration, and living service agglomeration all significantly improve carbon efficiency. The impacts of manufacturing agglomeration and high-end productive service agglomeration on carbon efficiency are ranked from large to small in the order of the Pearl River Delta, the Yangtze River Delta, and the Circum-Bohai Sea Economic Zone. The low-end productive service agglomerations in the Yangtze River Delta and Pearl River Delta have negative effects on carbon efficiency with the crowding effect.

Accordingly, we propose several policy implications.

The empirical examination showed that manufacturing agglomeration and high-end productive service agglomeration improve carbon efficiency in the city itself and surrounding cities. With the purpose of keeping stable economic growth and reducing CO2 emissions, the Chinese government should “unbound” the factor market and establish a “common market” between cities and release the “institutional dividend” through breaking down administrative barriers and regional trade protectionism. According to the heterogeneity of the industrial agglomeration effects, the government should accord priority to specific industries from the perspective of regional comparative advantages. It should avoid the excessive pursuit of industrial advancement in the competition of cities, especially the “comparison and pursuit” for the development of high-end productive service industry among cities.

At the level of megalopolises, compared with the national level, the agglomerations of manufacturing and high-end productive service have greater impacts on carbon efficiency. The central cities have more chances to obtain institutional and policy priority and should cultivate diversified high-end productive service industries as well as become the “technology pool” and “market area” for improving carbon efficiency in native and surrounding areas. Since the manufacturing and low-end productive service industries are more dependent on savings on production costs and improvement of transaction efficiency, the common cities should form complementary functions for central cities and make the coordinated development of manufacturing and low-end productive services.

Data Availability

The original data sources were the China City Statistical Yearbook, China Energy Statistics Yearbook, and Beijing Shu Huitong Environmental Technology Research Institute Data Service Network (http://www.3edata.com).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (grant no. 71704095) (a study on regional air pollutant emission reduction target, step-by-step benchmark path and technical prediction based on environmental production technology) and Young Scholars Program of Shandong University, Weihai.