Impact of Urbanization and Industrial Structure on Carbon Emissions: Evidence from Huaihe River Eco-Economic Zone

: Rapid urbanization in China greatly contributes to carbon emissions, while the industrial structure greatly contributes to changes in the variation of carbon emissions. This research addresses the impact of urbanization and industrial structure on carbon emissions from 2010 to 2018, by focusing on the Huaihe River Eco-economic Zone, which is an important economic corridor along the north–south division of China. Although many studies have focused on investigating the impact of urbanization or industrial structure on carbon emissions, few studies further addressed an analysis of the impact of both on carbon emissions, using multiple measurement models. This paper reveals the holistic and local impact of industrial structure and urbanization on carbon emissions, by integrating a threshold regression model with geographically weighted regression. The results are as follows: (1) From a holistic point of view, industrial structure and urbanization had both, single threshold, and double threshold effects on carbon emissions in the Huaihe River Eco-economic Zone. (2) From a regional perspective, the coefﬁcients of industrial structure on carbon emissions were all positive, but the rate of increase gradually slowed down. The coefﬁcients of urbanization on carbon emissions were all negative, reaching a maximum value of negative effect in 2013. Understanding the holistic and local impact of urbanization and industrial structure on carbon emissions provides governments with differentiated and forward-looking suggestions for mitigating carbon emissions in the Huaihe River Eco-economic Zone.


Introduction
The increase in carbon emissions has become a serious social problem, because it hinders the development of society and the global economy [1]. The upgrading of industrial structure and urbanization lead to an increase in carbon emissions. According to the World Bank, China's carbon emissions have exceeded the sum of the European Union and the United States; the three territories with the highest carbon emissions in the world [2]. Their total carbon emissions accounted for more than half of the world. China's carbon emissions reached 9.429 Gt in 2018, accounting for 27.8% of global carbon emissions [3]. As the largest developing country and the largest carbon emitter, China strives to achieve its promise to reach the peak of its total carbon dioxide emissions around 2030 and become carbon neutral by 2060 [4].
The Huaihe River Eco-economic Zone (HREZ) covers the provinces of Jiangsu, Henan, and Anhui. Accordingly, tertiary industry in the region accounts for about 10%, and the average annual growth rate of urbanization is about 3.56%. Specifically, the backward production pattern leads to increased carbon emissions in the zone. Furthermore, the carbon emissions of the HREZ are significantly higher than that of other regions. This study is conducive to a reduction of carbon emissions in the region, and even in China as a whole, as well as the early achievement of the carbon neutrality target. Therefore, it is imperative for the HREZ to focus on the impacts of urbanization and industrial structure on carbon emissions. At the same time, it is important to understand the internal mechanism of carbon emissions in the zone, which would be helpful to achieve China's carbon emissions reduction commitments.
Urbanization and industrialization are the main characteristics of the economic growth in China and promote rapid increases in carbon emissions [5]. Many methods for urbanization and industrial structure, such as the STIRPAT model, the Kaya index method, the LDMI method, and the threshold model, have been used to study the influencing factors on carbon emissions [6]. The panel threshold regression model aims to incorporate a certain threshold value as an unknown variable in the regression model, and to empirically test and estimate the corresponding threshold and the effect of the threshold [7][8][9]. In addition, the STIRPAT equation helps to disaggregate the factors influencing carbon emissions. Many researchers have proposed the STIRPAT equation, to explore the impact of carbon emissions [10][11][12][13][14]. With the development of spatial measurement in recent years, several scholars have included spatial factors into the framework for consideration of carbon emission driving factors. According to Tobler's first law of geography, all attribute values on a geographic surface are related to each other, but closer values are more strongly related than more distant ones [15]. Geographically weighted regression (GWR) is a local spatial regression technique for modeling the varying relationships in many application domains. Regression is a statistical method used to find the relationship between two or more variables that are interdependent and affected by one another.
Although the economy [16], technology [17], urbanization, industrial structure, trade openness [18][19][20], energy consumption [21][22][23], and financial development [11] have been considered in various previous studies, the analysis of both urbanization and industrial structure has not been addressed. Moreover, the optimization of urbanization and industrial structure can help achieve national targets for reducing carbon intensity.
Therefore, although urbanization and industrial structure are two important factors influencing carbon emissions, the combined impact of the two on carbon emissions is not fully understood. Second, the threshold model and GWR were not adequately combined to study carbon emissions in previous studies. The threshold model and GWR are beneficial for studying the carbon emissions of the whole region and each part of the region, respectively. Third, the HREZ has not been adequately addressed. It is one of the most promising areas in the east-central part of China, and the research of carbon emissions is of great value for its long-term development. This paper explores the holistic and local impact of industrial structure and urbanization on the carbon emissions of the HREZ from 2010 to 2018. We aimed to understand how they function, regarding carbon emissions. The main contributions are as follows: (1) This paper combines a threshold model with GWR, to study the impact of urbanization and industrial structure on carbon emissions, from both holistic and local perspectives. (2) Through research on the HREZ, relevant policies can be extended to other regions to mitigate environmental pressure.
With the development of a low carbon economy, the driving factors of carbon emissions have consistently been a research hotspot. Urbanization and industrialization are the main characteristics of economic growth in China and promote rapid increases in carbon emissions. At present, research on carbon emissions in combination with urbanization and industrial structure is mainly concentrated in three strands. The first strand is whether an environmental Kuznets curve (EKC) of carbon emissions exists. Namely, that carbon emissions will eventually appear as an inverted U-shaped curve with the development of urbanization and industrial structure [24][25][26][27][28]. For instance, [10] investigated the heterogeneous dynamic impacts of rural-urban migration and industrial agglomeration on environmental degradation. The combined impacts of both industrial agglomeration and rural-urban migration confirmed EKC in models of the tertiary sector. Reference [29] used the EKC to discuss the relationship between urbanization and carbon emissions in three typical models. The second strand is to use various models or methods to explore the causal relationship and driving factors among urbanization, industrial structure, and carbon emissions, such as cointegration analysis, structural decomposition models, decou- pling theory, and LMDI [30][31][32][33][34], which showed that urbanization level is the main driving force promoting decoupling elasticity, both in the long-term and short-term, contributing to a significant increase in carbon emissions. Reference [35] explored the threshold effects of industrial structure and urbanization on carbon emissions. With an increase in the level of population aging, industrial structure and urbanization will lead to an increase in carbon emissions, but the relationship fluctuates in three intervals. The last strand is the exploration of carbon emission rights and carbon compensation to establish efficient and stable operations in the carbon market to achieve emissions reduction targets [36,37]. For example, [37] improved the calculation method of the interregional carbon compensation cost, based on the theoretical carbon deficit. While, reference [38] examined how the announcement of an organizational carbon offset program affected consumption in a CO 2 -intensive everyday activity. As for the impact of other factors on carbon emissions, foreign direct investment (FDI) has both a promoting and an inhibiting effect on carbon emissions performance [39,40]. There is a relationship among economic growth, energy consumption structure, and carbon emissions [41,42]. In terms of technology, low-carbon technology can increase the potential for carbon emissions reduction [43][44][45].
In terms of regional measurement, scholars have focused more on studies at the national provincial, national strategic regional, and municipal levels. Among the national strategic regions, scholars have measured the carbon efficiency and influencing factors of the Yangtze River Economic Belt [46][47][48], Beijing-Tianjin-Hebei Urban Agglomeration [29,49,50], Belt and Road Initiative [40,51], and Pearl River Delta [52]. Although the HREZ is very important to China's economic development and for achieving goals for carbon abatement, there is currently a lack of research on such aspects.
Scholars have studied the influencing factors of carbon emissions more widely, and have made many achievements regarding the influencing factors on carbon emissions. However, there has been very little research on the HREZ. The government has not yet realized the importance of the study of carbon emissions in the HREZ for the national carbon abatement. The HREZ is located between the Yangtze River Economic Belt and the Huanghe River Basin, but its relative overall economic development is lagging. It is one of the regions with the greatest development potentials in the central and eastern part of China. For example, the GDP per capita was 49121 CNY/person in 2018, which is a fraction of the Yangtze River Economic Belt and the Huanghe River Basin. The HREZ, with a planned area of 243,000 square kilometers, has great development potential. This has important strategic significance for promoting the balanced development of China, such as being a bridge connecting the south regions to the north regions, and the east regions to the west regions. It includes several key development areas of heavy industry in China. Therefore, it is valuable to research the factors influencing carbon emissions in the HREZ. The HREZ overlaps with other important strategic regions in China, so its successful experience with carbon abatement could be applied to other geographical regions in China.
As the goal of carbon emissions reduction is to set up joint efforts, with emissions mitigation suggestions for each region, it is thus necessary to identify changes of carbon emissions under the consideration of regional variations; however, to our knowledge, corresponding research on the HREZ from a differentiation perspective is absent. Thus, this paper has the following novel features: (1) From multiple perspectives, the article combines industrial structure and urbanization to study their impacts on carbon emissions. (2) The combination of a threshold model and GWR is conducive to analysis and research from both planar and spatial perspectives. (3) The HREZ is the research object, and currently lacks theoretical and empirical studies on improving carbon emissions efficiency in China. This is essential for China's carbon emission reduction and regional development.

Threshold Regression Model
Following [35], we apply the static panel threshold regression model proposed by [53] to attain the goal of this paper. According to the STIRPAT model, this paper selected variables of carbon emissions (C, tones), urbanization (U,%), industrial structure (M,%), structure of energy consumption (S,%), energy intensity (Q, tones/million CNY), and GDP per capita (Y, CNY/person) to construct the equation. In addition, industrial structure and GDP per capita reflect economic growth effects, energy intensity and energy consumption structure reflect technological advances, and urbanization reflects social effects. The static panel threshold regression models, in a form of STIRPAT, are as follows: where C i,t denotes carbon emissions (tones), M i,t denotes industrial structure (%), U i,t denotes urbanization (%), S i,t denotes energy consumption structure (%), Q i,t denotes energy intensity (tones/million CNY), Y i,t denotes GDP per capita (CNY/person), µ is a constant term, ε is a random disturbance term, subscripts i, t are prefecture-level administrative unit and year, LnX i,t is a threshold variable, γ is a threshold value, and I(·) is an indicator function; when LnX i,t ≤ γ, I(·) = 1 , when LnX i,t > γ, I(·) = 0 .

GWR Method
GWR relies on standard multiple parameter regression. Unlike classical regression, the coefficients are not constant in this technique. It is an effective technique for exploring spatial nonstationarity. GWR model is specified in the following form: where y i is the value at location i, (u i , v i ) are the geographic coordinates of position i, β 0 (u i , v i ) and β j (u i , v i ) are the intercept and local coefficients at location i, respectively, k is the number of factors, x ik is the explanatory variable at location i, and ε i accounts for the random error term at location i.

Data Sources and Processing
Referring to the relevant literature, it was found that scholars in various countries had different emphases in the selection of influencing factors of carbon emissions. This paper selected the variables of carbon emissions, urbanization, industrial structure, structure of energy consumption, energy intensity, and GDP per capita. These were all derived from the statistical yearbooks of cities in the HREZ. This paper's analysis covers 5 provinces (excluding Henan, Jiangsu, Anhui, Hubei, Shandong) in China. The related data of all variables, spanning 2010 to 2018, were sourced from the China Statistical Yearbook and China Energy Statistical Yearbook [54,55]. To eliminate the impact of price inflation, GDP per capita was converted into constant prices from 2010 using the corresponding indices. A long period can clearly see show an increasing or decreasing trend in carbon emissions. In 2018, the state compiled the "Development Plan of Huaihe River Eco-economic Zone". Therefore, a summary of the development over a long period before 2018 is more conducive for the targeted measures for future carbon abatement.
There are two common methods to calculate carbon emissions. One is based on the reference method provided by the Intergovernmental Panel on Climate Change (IPCC) [56]. The IPCC is the United Nations body for assessing the science related to climate change. IPCC scientists provide a comprehensive summary of the drivers of climate change, its impacts, and future risks, and how adaptation and mitigation can reduce those risks [57]. The other is the input-output (I-O) method. The I-O method represents the economic dependencies between different sectors and uses matrix representations to describe the inter-industry relationships in the economy [58]. While, the input-output method is updated every five years, the method provided by the IPCC can better capture the impact of factors on regional carbon emissions on a year-by-year basis. The carbon emissions were calculated as follows: where i is the fossil fuel type; E i is the terminal consumption of the i th fossil fuel (tones); NCV i is the average low calorific value (kj/kg); CEF i is the respective carbon emissions coefficient for the 8 fossil fuels (kg/tj); CC i is the carbon content (t/tj); COF i is the carbon oxidation factor; and 44 and 12 are the molecular weights of CO 2 and C, respectively.

Results of the Threshold Regression Model
The results of the fixed effects model are as follows. The R 2 value was 0.81 and the p-value also passed the test. The p-value of urbanization, industrial structure, energy intensity, and energy consumption structure were ≤1% or lower, and the p-value of GDP per capita was 0.028. Each variable passed the significance test in the model. The details are shown in the following, Table 1.  Note: values in parentheses are the T-statistics, **, *** denote significance at 10%, 5%, and 1% levels, respectively. Then, we tested the threshold effect. We modified the threshold variables and the threshold explanatory variables. The results indicated the existence of single or multiple threshold effects with GDP per capita and energy intensity as threshold variables, respectively. First, the results of GDP per capita as a threshold variable are presented in Table 2. The p-value of urbanization, energy intensity, and energy consumption structure were all ≤1% or lower. From the significance of the model, each variable passed the significance test.
From the perspective of industrial structure, the threshold of LnY i,t is −2.725. Before and after the threshold, each unit increase in industrial structure leads to a reduction of 0.35 units followed by an increase of 0.075 units in carbon emissions.
(1) In the early stage of economic development, the HREZ was committed to developing strategic emerging industries. The region made use of new knowledge and technologies to develop the industrial structure in the direction of low consumption, low pollution, and high efficiency. Meanwhile, the HREZ emphasized the development of an ecological industry. Through the introduction of clean energy and advanced technology, the government aimed at transitioning the HREZ from resource-intensive industry to technology-intensive industry. (2) With the increase of GDP per capita, tertiary industry tended to be saturated. The zone relied on secondary industry, which is the traditionally preponderant industry. The development of secondary industry was more difficult to separate from the From the perspective of urbanization, the thresholds of LnY i,t are −2.725 and −1.389. When the GDP per capita exceeds the first and second thresholds, each unit of urbanization, first makes carbon emissions decrease by 0.08, then increase by 0.08, and finally decline by 0.001.
(3) Before the first threshold, the GDP per capita of the zone was low and the urbanization process was just beginning. The workforce moved from the countryside to the cities. The gathering of urban population reduced the cost of public goods. Along with the change of lifestyle in the zone, the quality of residents continued to improve. Residents have become more aware of environmental protection. The change of residents has greatly reduced the negative impact of urbanization on the environment. (4) Between the first and second thresholds, the urbanization process entered the middle stage and became the main driving factor affecting carbon emissions. The consumption structure and level of residents improved rapidly. The increasing use of energy-intensive products, such as cars and appliances led to high carbon emissions in the HREZ. (5) After the second threshold, the impact of urbanization on carbon emissions slowed down significantly. Despite the advanced economy, the spillover from industrial consumption was compensated by urbanization to a certain extent. The corresponding likelihood ratio function graphs, with GDP per capita as the threshold variable, are shown in Figure 1. Moreover, the results of energy intensity as a threshold variable are shown in Table 3. The p-value of energy intensity and energy consumption structure are ≤1% or lower, and the p-value of urbanization is 0.042. However, the p-value of GDP per capita does not pass the significance test.
From the perspective of industrial structure, the thresholds of LnQ i,t are −2.013 and −2.410. As the energy intensity exceeded the first and the second thresholds, each unit of industrial structure made carbon emissions first decrease by 0.16, then increase by 0.13, and finally decrease by 0.01.
(1) Originally, secondary industry was not the main force of economic development and was less dependent on energy consumption. Due to the low energy intensity, the upgrading of industrial structure meant more production resources were given to tertiary industry. The clean nature of tertiary industry led to a reduction in carbon emissions. (2) After the energy intensity exceeded the first threshold, the energy dependence of the zone was significantly higher than before. Industries with a high energy consumption and high emissions accounted for a large proportion of economic development. Some regions were relatively backward in technology, so the development of tertiary industry could not offset the high carbon emissions caused by the secondary industry. (3) With the further increase of energy intensity up to the second threshold, tertiary industry reached its heyday. The regional industrial structure was continuously optimized under the requirements of the new development concept. Green development of secondary industry and the low carbon consumption of tertiary industry compensated for the spillover from the high energy consumption of industry.
From the perspective of urbanization, the thresholds of LnQ i,t are −1.321 and −1.709. As the energy intensity exceeded the first and second thresholds, each unit of urbanization made carbon emissions, first decrease by 0.01, then increase by 0.08, and finally decrease by 0.03. Moreover, the results of energy intensity as a threshold variable are shown in Table  3. The p-value of energy intensity and energy consumption structure are ≤1% or lower, and the p-value of urbanization is 0.042. However, the p-value of GDP per capita does not pass the significance test. Table 3. Estimation of panel threshold model with energy intensity as the threshold variable.
(1) At the outset of the urbanization process, the increase of urbanization meant more rural populations moving to the cities. The industrial development of the zone was in the primary stage. A low energy intensity meant low carbon emissions. The positive impact was optimal when the first threshold was reached. (2) With the increase of energy intensity, energy-intensive industries in the HREZ developed rapidly. The lack of industrial innovation capacity and the pressure of economic transformation were the main problems of the development. With the migration of a large number of people, the economic development relied on high-consumption and high-pollution secondary industry. The high energy intensity led to a significant increase in carbon emissions in the HREZ. (3) After the energy intensity exceeded the second threshold, the industrial development reached a certain scale. The focus of development shifted to tertiary industry, which is dominated by the service industry. The tertiary sector emits fewer pollutants and has a suppressive effect on carbon emissions. The low carbon consumption of the tertiary sector compensated for the high carbon consumption of the secondary sector.
Corresponding likelihood ratio function graphs, with energy intensity as the threshold variable, are shown in Figure 2.

Estimation Results of GWR Model
The Moran's Index values obtained for all the years are greater than zero, which suggests that there is positive autocorrelation or a highly clustered pattern. The same is seen by looking at the p-values. The obtained p value is less than 0.05 for all the years, ruling out the basic assumption of randomness and independence in the data values. The results of Moran's I are shown in Table 4. To show the spatial heterogeneity of the impact of industrial structure on carbon emissions, the spatial impact of industrial structure from 2010 to 2018 is visualized. The

Estimation Results of GWR Model
The Moran's Index values obtained for all the years are greater than zero, which suggests that there is positive autocorrelation or a highly clustered pattern. The same is seen by looking at the p-values. The obtained p value is less than 0.05 for all the years, ruling out the basic assumption of randomness and independence in the data values. The results of Moran's I are shown in Table 4. To show the spatial heterogeneity of the impact of industrial structure on carbon emissions, the spatial impact of industrial structure from 2010 to 2018 is visualized. The results are shown in Figure 3a-d. industries, such as labor and land. In order to maintain the current scale of other industries, the production relied more on large mechanical production equipment and energy consumption. The coefficient was lower in Xiaogan, Suizhou, Nanyang, Xinyang, Zhumadian, Pingdingshan, Luohe, Zhoukou, Shangqiu, and Bozhou. These cities are deep inland and had a weak foundation in secondary industry. Due to the backwardness of their production, the reliance on mechanical equipment and energy was low. They relied more on labor for production, so the impact on carbon emissions was lower. From 2012 to 2015, the trend was roughly the same as before. The coefficient increased significantly, especially in Jiangsu and Shandong Provinces, which were significantly higher than those in southwestern cities. Due to the development of transportation and high-tech industries, the rapid growth of the transportation sector made it the largest source of carbon emissions outside of industry. The negative effect of the transportation sector on carbon emissions was less significant in the central region, because of less convenient transportation and a lower population.
However, the cities in the southwestern zone underwent a rapid economic transformation after 2016. The government promoted the introduction of new industries. All cities reshaped industrial structures and promoted industrial clustering. Industrial restructuring drove down the growth rate of carbon emissions. The gap between the southwestern zone and other cities became smaller. The positive influence of industrial structure on carbon emissions in the HREZ from 2010 to 2018 become more and more obvious, but the trend of growth slowed down. In general, it showed a decreasing trend from the northeast to the southwest, and a more obvious east-west gap.
The cities with high coefficients in 2010 were Suqian, Lianyungang, Huai'an, Yancheng, Yangzhou, and Taizhou. The positive influence on carbon emissions was significantly higher in the eastern cities than in the western cities. With their excellent location and resources, these cities seized the opportunity for development. The accelerated development of tertiary industry led to the accelerated flow of production factors to the tertiary industries, such as labor and land. In order to maintain the current scale of other industries, the production relied more on large mechanical production equipment and energy consumption. The coefficient was lower in Xiaogan, Suizhou, Nanyang, Xinyang, Zhumadian, Pingdingshan, Luohe, Zhoukou, Shangqiu, and Bozhou. These cities are deep inland and had a weak foundation in secondary industry. Due to the backwardness of their production, the reliance on mechanical equipment and energy was low. They relied more on labor for production, so the impact on carbon emissions was lower.
From 2012 to 2015, the trend was roughly the same as before. The coefficient increased significantly, especially in Jiangsu and Shandong Provinces, which were significantly higher than those in southwestern cities. Due to the development of transportation and high-tech industries, the rapid growth of the transportation sector made it the largest source of carbon emissions outside of industry. The negative effect of the transportation sector on carbon emissions was less significant in the central region, because of less convenient transportation and a lower population.
However, the cities in the southwestern zone underwent a rapid economic transformation after 2016. The government promoted the introduction of new industries. All cities reshaped industrial structures and promoted industrial clustering. Industrial restructuring drove down the growth rate of carbon emissions. The gap between the southwestern zone and other cities became smaller.
In summary, the government needs to avoid the negative influence of transportation and high-tech industries on carbon emissions. At the same time, traditional industries should be reformed and upgraded to avoid a reliance on large-scale equipment with high energy consumption.
Similarly, the spatial impacts of urbanization are shown in Figure 4a-d. The negative impact of urbanization on carbon emissions increased at first and then decreased. The growth of carbon emissions was inhibited by urbanization in the early stage. When urban development strategies and spatial layouts were adjusted, the energy-saving effects of urbanization became apparent.

Discussion
In this paper, we analyzed the factors and spatial patterns influencing carbon emissions in the HREZ from 2010 to 2018. Our findings are summarized as follows: (1) Urbanization and industrial structure have double and single threshold effects on the The cities with low coefficients in 2010 were Xuzhou, Suqian, Lianyungang, Huai'an, Yancheng, Yangzhou, Taizhou, Linyi, Suzhou, and Chuzhou, all of which are located in the coastal region. With the input of production factors from abroad, the increase in urbanization rate made the industries gather in cities. Manufacturing industries with serious environmental pollution were gathered around cities. The increase in urbanization increased carbon emissions. While the northwestern cities, such as Zaozhuang, Huaibei, Fuyang, and Xinyang, had a less negative impact on carbon emissions. The urbanization of these cities started later and the urbanization process was much smaller than the eastern cities. Their development mostly relied on secondary industry, so urbanization could not offset the negative impact of secondary sector on carbon emissions.
The coefficient reached the minimum value in 2013. The coefficients of Jining, Zaozhuang, Linyi, Lianyungang, Xuzhou, Suizhou, Suqian, Huai'an, and Yancheng were significantly lower than the western cities. The distribution in the following two years was consistent with that in 2013. The boundary of Heze, Huaibei, and Suzhou showed a clear east-west gap, in which the correlation of urbanization with carbon emissions in the eastern zone was significantly smaller than in the western zone. The rapid urbanization of eastern cities led to a disorderly development and blind expansion of urban scale, which brought a series of problems in urban land and led to a significant increase in carbon emissions.
The coefficient continued to increase from 2016 to 2018. The zone attached importance to the transportation and tertiary industries. The government solved the problems of employment and resource allocation. However, the coefficient in 2018 was larger than in 2010. The negative effect of urbanization was significantly reduced. The rapid development of urbanization led to the concentration of the urban population and the growth of urban areas, which was a great test for the urban environmental situation. At the same time, urbanization brought an improvement in the living standard of residents and the quality of the population.
In a word, the zone was in the stage of urbanization development, from high speed to high quality. Developed cities should play a radiating role, to avoid the waste of resources and pollution caused by blind expansion. The economically backward regions should accelerate their urbanization process.

Discussion
In this paper, we analyzed the factors and spatial patterns influencing carbon emissions in the HREZ from 2010 to 2018. Our findings are summarized as follows: (1) Urbanization and industrial structure have double and single threshold effects on the carbon emissions, respectively, with GDP per capita as the threshold variable. Urbanization and industrial structure have double threshold effects on the carbon emissions, with energy intensity as the threshold variable. (2) The coefficients of industrial structure are all positive, but the rate of increase gradually decreases. They show a decreasing trend from northeast to southwest. Before 2016, the positive impact on carbon emissions in Jiangsu and Shandong provinces was significantly higher than in the southwestern cities. Cities bounded by Shangqiu, Huabei, Bengbu, and Huainan had developed tertiary industries and other industries relying on large machinery production equipment in order to compensate for the development of tertiary industry. After 2016, the zone began to move in a favorable direction, towards reducing carbon emissions, gradually meeting the low carbon development targets of the thirteenth Five-Year Plan of China. In order to improve the efficiency of carbon emissions and reduce the regional differences in the HREZ, some suggestions are proposed as follows: (1) The industrial structure and urbanization development speed of the HREZ need to be adjusted according to the thresholds. Local governments in the HREZ should adjust regional measures and policies according to the marginal benefits of industrial structure and urbanization. (2) Eastern cities in the zone need to avoid the adverse effects of transportation and high-tech industries on carbon emissions, due to the unreasonable development of industrial structure. In addition, it also needs to reduce the dependence of traditional industries on large equipment with a high energy consumption.