Next Article in Journal
Research Progress and Hotspot Analysis of Urban Heat Island Effects Based on Cite Space Analysis
Previous Article in Journal
Impact of Land Value Tax on the Equity of Planning Outcomes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial–Temporal Pattern Characteristics and Impact Factors of Carbon Emissions in Production–Living–Ecological Spaces in Heilongjiang Province, China

1
School of Architecture, Harbin Institute of Technology, Harbin 150006, China
2
Key Laboratory of Cold Region Urban and Rural Human Settlement Environment Science and Technology, Ministry of Industry and Information Technology, Harbin 150006, China
*
Author to whom correspondence should be addressed.
Land 2023, 12(6), 1153; https://doi.org/10.3390/land12061153
Submission received: 8 May 2023 / Revised: 27 May 2023 / Accepted: 29 May 2023 / Published: 30 May 2023

Abstract

:
Under the threat of global climate change, China has proposed a dual carbon goal of peak carbon and carbon neutrality. As the vital carrier for territorial spatial planning, production–living–ecological (PLE) spaces drive carbon emissions and are important to the dual carbon goals. In this study, carbon emissions and sinks of PLE spaces in cities in Heilongjiang Province from 2005 to 2020 were calculated and spatial–temporal changes were analyzed. The carbon emission structure was analyzed in segmentation sectors. The land use changes and socioeconomic factors on carbon emissions were analyzed, and emission reduction strategies were implemented. The results show the following: (1) Carbon emissions from production and living spaces increased yearly. Carbon sinks were smaller than emissions, but capacity was stable. (2) Higher-emission cities were concentrated in southwest Heilongjiang, and carbon emission differences between regions gradually increased. (3) Among carbon emission sectors, agricultural and household made up smaller proportions, while animal husbandry, industrial, transportation, and traffic travel contributed most. Carbon emission structures were transformed by adjusting urban development and industrial structure. (4) For most cities, industrial space was the main emission space, but agricultural production and urban–rural living spaces dominated in some cities. (5) GDP, urbanization rate, and area of city paved roads suppressed emissions in cities with decreased carbon emission grades. The industrial structure and coal consumption inhibited emissions in cities with maintaining and increasing carbon emissions grades.

1. Introduction

Climate change has become a major challenge for society as a whole. Climate change caused by greenhouse gas emissions, especially carbon dioxide emissions, constantly reminds mankind to achieve carbon emission reduction commitments promptly [1,2,3]. In the past few decades, China has been the world’s largest CO2 emitter, accounting for approximately 29% of global CO2 emissions [4]. To alleviate the dual pressures of domestic environmental degradation and international climate negotiations, China announced greenhouse gas emission targets in the “U.S.–China Joint Announcement on Climate Change”, the core of which is a commitment to peak CO2 emissions by 2030 or sooner [5]. In the 75th session of the United Nations General Assembly, China reaffirmed the “dual carbon goals” of peaking CO2 emissions by 2030 and achieving carbon neutrality by 2060 [6]. China is in the development stage and faces the crises and challenges of tightening resource constraints, increasingly serious environmental pollution, and gradual degradation of ecosystems [7].
The planning concepts of “intensive and efficient production space, livable and moderate living space, and beautiful ecological space” proposed by China at the 18th National Congress are known as production–ecological–living (PLE) spaces and have become important means to implement the dual carbon goals [8]. PLE spaces not only indicate space, but also reflect their production, living, and ecological functions. A consensus has been reached regarding the classification of PLE spaces based on the multifunctionality of land use [9,10]. Many studies on urban functional space and multifunctionality have similarities with Chinese studies on land use classification [11,12,13], that is, they carry out studies related to the coordination and balance of regional production, ecological, and living space functions from the perspective of spatial planning and management. Recently, various studies on carbon emissions in PLE spaces have been conducted, including the construction of a theoretical framework [14], identification and classification of spatial functions [15,16,17], and evaluation and optimization strategies [18,19,20].
Since China put forward territorial spatial planning and the dual carbon goals, carbon emissions from territorial spaces have become a focal point of research. Xu et al. constructed a greenhouse gas accounting framework suitable for territorial spatial planning that solved the problem of carbon emission accounting in urban and central urban areas [21]; King et al. analyzed the distribution characteristics of carbon emissions in territorial space and compared the differences in each functional area [22,23,24]. Some scholars have also explored territorial spatial planning methods and optimization strategies of provincial, municipal, and county administrative units with the goal of carbon emission reduction [25,26,27]; some scholars have further studied the carbon emission reduction of territorial space through carbon balance zoning and carbon emission reduction capacity measurement [28]. For example, Feng et al. constructed a carbon metabolism network model in an urban PLE space, analyzed the changes in horizontal carbon flow, and evaluated the comprehensive effect of PLE space changes using an ecological network analysis method [29]. Guo et al. calculated the carbon emissions of PLE spaces in Chongqing and used the STIRPAT model to analyze influencing factors [30]. Commonly used carbon emission analysis methods include the LMDI model [31], STIRPAT model [32], and Kaya constant equation [33]; however, each method has a different research applicability. Duro et al. decomposed international inequalities in per capita CO2 emissions into the Kaya equation to analyze its inequality components [34]. Liu et al. adopted STIRPAT to assess the effect of urban household energy consumption [35]. Yang et al. examined the carbon emissions of the transportation sector of the Yangtze River Economic Belt based on the LMDI model [36].
The existing research provides solid theoretical support and rich methods for this study. However, most studies have focused on carbon emissions from territorial space from a macro perspective. Regarding provincial carbon emissions, few studies have calculated and analyzed municipal cities in the same provincial area from the perspective of PLE spaces. Furthermore, both land use and socioeconomic factors impact carbon emissions, but few scholars have discussed this impact. The study of carbon emissions from PLE spaces can indicate the contribution of carbon emissions in each city more deeply for clearer comparison. Research from the perspectives of PLE spaces, land use, and segmentation sectors can compensate for these deficiencies. The STIRPAT model could add impact factors according to the situation of the research object. It is more conducive for policymakers carrying out planning to promote urban carbon neutrality while simultaneously considering the overall characteristics of the province.
Heilongjiang Province is the main grain-producing area in China and is crucial to national food security. Moreover, Heilongjiang Province was an early industrial base. Upgrading energy and industrial structures is an important step in achieving the dual carbon goals. Population shrinkage has increased the challenge of achieving competitive development in Heilongjiang Province [37]. There are certain particularities in low-carbon development in other provinces and cities.
Therefore, this study takes Heilongjiang Province as an example, classifies the PLE spaces, constructs a calculation system for PLE space carbon emissions, and calculates emissions for each city. Subsequently, the spatial and temporal differentiation characteristics were analyzed. On this basis, the subdivision of carbon emission sectors was used to further analyze the main carbon sources and carbon emission structures of each city. This was followed by an analysis of the impact of land use change on carbon emissions. Finally, combined with the development characteristics of cities in Heilongjiang Province, a STIRPAT model was constructed to explore the factors influencing carbon emissions. The results provide a planning basis and scientific reference for territorial spatial planning and carbon emission reduction in Heilongjiang Province. The technical roadmap is shown in Figure 1.

2. Study Area and Methodologies

2.1. Study Area and Data Sources

2.1.1. Study Area

Heilongjiang Province is located in Northeast China, adjacent to eastern Russia, bordering eastern Mongolia, and is the northernmost and easternmost provincial administrative region (Figure 2), with an area of 473,000 km2, accounting for 4.9% of China’s total land area [38]. Heilongjiang Province is the primary energy industry base in China and contributes a large amount of carbon emissions. Additionally, the province is the main granary ensuring China’s food security, and protecting its cultivated land areas is of great significance to strategic national development. Heilongjiang Province governs 13 prefecture-level administrative regions, including 12 prefecture-level cities and 1 prefecture. Owing to data limitations, the study area included the 12 prefecture-level cities and excluded Daxing’anling prefecture.

2.1.2. Data Sources

The land use data for 2005, 2010, 2015, and 2020 used in this study were obtained from the Resource and Environmental Science Data Center of the Chinese Academy of Sciences (http://www.resdc.cn/) (accessed on 20 October 2022). The administrative boundary data for Heilongjiang Province were obtained from the National Geographic Information Resource Catalog Service System (https://www.Webmap.cn) (accessed on 20 October 2022). Statistically related data on energy consumption and economy were derived from or calculated based on the yearbooks of national, provincial, and local cities. The conversion coefficients of various energy sources were obtained from the China Energy Statistics Yearbook.

2.2. Classification System for Production–Living–Ecological Spaces

The construction of a scientific and reasonable PLE space classification system is the premise and basis for studying the characteristics and influencing factors of PLE space carbon emissions [39]. The territorial spatial pattern is a comprehensive reflection of the coupling of and interactions between natural ecological processes and human social systems, which perform a variety of functions [40,41,42,43]. Scholars have conducted extensive research on the classification of PLE spaces, mainly based on the dominant functions of different land use types [44]. Based on previous studies [45,46,47], this study considered the versatility of territorial space as the starting point, combined the land use classification system of the Chinese Academy of Sciences and the classification of land use status (GB/T21010-2007), classified the PLE spaces in the study area, and constructed a classification system (Table 1).

2.3. Methods

2.3.1. Calculation of Carbon Emissions in PLE Spaces

In PLE spaces, production and living spaces are the main carbon sources, whereas the ecological space is the carbon sink. The carbon emissions and sinks of each space were calculated separately, and the carbon emissions of Heilongjiang Province were summarized. The specific formula is as follows:
C E P L E   S p a c e = C E P S + C E L S C E E S
where C E P L E   S p a c e represents the total carbon emissions (CE) of the PLE space, C E P S represents the carbon emissions of the production space, C E L S represents the carbon emissions of the living space, and C E E S   represents the carbon sink in the ecological space.
  • Production spaces
The carbon emissions from the production space mainly originate from agricultural production spaces and industrial production spaces. Carbon emissions from APS include agricultural production activities and animal husbandry. Carbon emissions from agricultural production activities include agricultural fertilization, crop cultivation, agricultural machinery use, irrigation, and CH4 emissions from paddy fields. Carbon emissions from animal husbandry include CO2 emissions from livestock respiration and CH4 emissions from livestock and poultry intestinal fermentation and feces. The primary livestock and poultry breeding varieties in Heilongjiang Province are pigs, cattle, sheep, and poultry [48]. CH4 emissions must be converted into carbon emissions. According to the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report, the greenhouse effect induced by 1 ton of CH4 is equivalent to that induced by 6.8182 tons of carbon [49,50,51,52]. The formula for calculating agricultural carbon emissions is as follows:
C E A P S = i = 1 n T i × α i
where CE denotes the carbon emissions in APS, Ti is the carbon emissions from each different APS source, and αi is the coefficient for each APS carbon source.
Carbon emissions from IPS include industry and transportation. Industrial carbon emissions were calculated using the estimation method proposed by the IPCC in 2006 and 2019 and the carbon emissions of 15 end-use energy sources, including raw coal, coke, natural gas, diesel, fuel oil, liquefied petroleum gas, petroleum coke, gasoline, heat, and electricity [53]. The formula for the calculations is as follows:
C E I n d u s t r y = i = 1 n E i × N C V i × E F i × O i
where C E I n d u s t r y represents industrial carbon emissions; i represents energy type; Ei represents energy consumption; NCVi represents average low-calorific value, that is, the heat released per unit of fossil energy combustion; EFi represents emissions coefficient, that is, the total carbon emissions per unit net calorific value generated by burning fossil fuels; and Oi represents carbon oxidation rate, that is, the oxidation rate of burning fossil fuels.
According to production and life functions, transportation can be divided into two parts: transportation outside the living space and transportation inside the living space.
The former portion of carbon emissions belongs to the production space which includes passenger transportation and freight transportation [54]. The latter is residents’ travel using different modes of transportation. To distinguish this from the transportation of the production space, the inside part is called traffic travel. Typically, it is difficult to compare the two types of transportation [55]. Therefore, passenger turnover was converted into freight turnover, and transportation carbon emissions were calculated. The relevant parameters are listed in Table 2, and the specific formula is as follows:
C E T r a n s p o r t a t i o n = i = 1 n ( T P a s s e n g e r × α i + T F r e i g h t ) × β i
In the formula, C E T r a n s p o r t a t i o n represents transportation carbon emissions, T P a s s e n g e r represents passenger turnover, α i represents the passenger and cargo conversion coefficient, T F r e i g h t represents passenger turnover, β i   represents the carbon dioxide emission coefficient of the ith transportation mode, and i represents the type of transportation. Transportation in this study included railways, highways, water transportation, and civil aviation and did not involve pipeline or other transportation modes.
The formula of carbon emissions in production is as follows:
C E P S = C E A P S + C E I n d u s t r y + C E T r a n s p o r t a t i o n
2.
Living spaces
Carbon emissions from living spaces comprise resident energy consumption and traffic travel. The carbon emissions from residential energy consumption were calculated using the electricity, natural gas, and liquefied petroleum gas consumed by residents. The specific formulas and coefficients refer to industrial carbon emissions. Carbon emissions from traffic travel were calculated based on ownership of private vehicles, buses, and taxis; average annual mileage; and carbon dioxide emission coefficient [56,57]. The relevant parameters are listed in Table 3, and the specific formula is as follows:
C E L S = C E R E C + i = 1 3 Q i × L i × β i
where C E L S   represents the carbon emissions from living space, C E R E C   represents the carbon emissions from resident energy consumption, and i represents different traffic travel modes, namely, private vehicles, buses, and taxis. Qi represents vehicle ownership of the i th means of travel mode, Li represents the average annual mileage of the i th means of travel mode, and βi represents the carbon emission coefficients of the i th means of travel mode.
3.
Ecological spaces
As the main carbon sinks, ecological spaces play an important role in the absorption of carbon emissions. The carbon sink capacity of an ecological space was calculated using carbon sink coefficients and areas of different land use types. Ecological spaces include forest ecological spaces, grassland ecological spaces, water ecological spaces, and other ecological spaces. The study area comprises a temperate coniferous and broadleaf mixed forest and temperate meadow steppe. The carbon sink coefficients of the FES and GES were 0.585 tC/hm2/a and 0.021 tC/hm2/a, respectively. Considering previous research on water carbon sinks in China and its northeastern region, the adopted WES carbon sink coefficient was −0.0395 tC/hm2/a [58,59,60,61]. OES refers to barren mountain gravel land, sandy land, bare land, and other difficult-to-use land. Its carbon emissions and absorption capacity are relatively weak, and the carbon absorption rate changes little with geographical location. According to previous research [62], the carbon sink coefficient of OES is 0.005 tC/hm2/a. The ecological space carbon sink equation is as follows:
C E E S = i = 1 4 A i × θ i
where CEES represents the carbon sink of an ecological space; i represents different ecological spaces, namely, FES, GES, WES, and OES; A i   represents the area of i th ecological spaces; and   θ i represents the carbon sink coefficient of the i th ecological space.

2.3.2. Spatial Autocorrelation

The spatial autocorrelation index can be used to measure the spatial aggregation or dispersion of regional carbon emissions. Autocorrelation analysis includes global and local spatial autocorrelation [63,64,65]. A global index was used to detect the spatial pattern of the entire region, and a single value was used to reflect the spatial autocorrelation of carbon emissions in the region. The global autocorrelation index is calculated as follows:
I = i = 1 n j = 1 n W ij ( x i   x ¯ ) ( x j   x ¯ ) [ i = 1 n j = 1 n W ij ( x i   x ¯ ) 2 ]
where Wij is the spatial weight matrix, xi and xj are the attribute values of spatial units, x is the average value, and n is the total number of cities. The global Moran index I ranges from −1 to 1. When I < 0, the distribution of the spatial elements is negatively correlated, and when I is closer to −1, the negative correlation is stronger. When I > 0, the distribution of the spatial elements is positively correlated, and when I is closer to 1, the positive correlation is stronger. When I = 0, the spatial elements are randomly distributed in space and there is no spatial correlation.

2.3.3. STIRPAT Model

A comparison of the three commonly used models for the analysis of carbon emission impact factors reveals that the factors in the decomposition of Kaya identity need to be analyzed year by year or by time and are constrained by the necessary identity. The LMDI method achieved no residual decomposition and quantified the contribution rate of each influencing factor in a specific year. However, the elasticity of each factor could not be examined; that is, when other factors remained unchanged, the change in carbon emissions caused by the change in one factor could not be determined. However, the traditional STIRPAT model includes only three factors—population, economic development level, and technical variables—that cannot comprehensively describe the impact of various socioeconomic aspects on CO2 emissions. In contrast to the Kaya constant equation and LMDI model, STIRPAT allows the impact of each factor to be estimated as a parameter, and researchers can extend the model for their own research purposes [66,67,68,69].
The STIRPAT model is an extensible stochastic environmental impact assessment model based on the IPAT model. As a nonlinear random model with multiple independent variables, the STIRPAT model can incorporate appropriate independent variables according to conditions specific to each location when conducting carbon emissions research. The IPAT model is intuitive and concise, but the STIRPAT model overcomes the limitations of the proportional relationship between influencing factors and dependent variables and increases its flexibility and adaptability using the following formula:
I = a Pb Ac Td e
where I, P, A, and T represent environmental pressure, population, affluence, and technical level, respectively. a is the model coefficient; b, c, and d are the exponential parameters of P, A, and T, respectively, and e is the random error of the model.
The equation is converted to logarithmic form as follows:
ln I = ln a + b ln P + c ln A + d ln T + ln e
To comprehensively and systematically analyze the factors influencing carbon emissions, gross regional product (GDP), industrial structure (IS), urbanization rate (UR), coal consumption (CC), area of city paved roads (CPR), per capita disposable income of urban residents (DIU), highway mileage (HM), construction land area (CLA), and total power of agricultural machinery (PAM) were added and expanded into the model. Because of the population loss in Heilongjiang Province, the urbanization rate was used to represent the population factor. The STIRPAT model was as follows:
ln C E = ln a + b 1 ln G D P + b 2 ln I S + b 3 ln U R + b 4 ln C C + b 5 ln C P R + b 6 ln D I U + b 7 ln H M + b 8 ln C L A + b 9 ln P A M + ln e
where a is the model coefficient; b1, b2, b3, b4, b5, b6, b7, b8, and b9 are exponential parameters; and e is the random error of the model.

3. Results

3.1. Results of Carbon Emissions in PLE Spaces

3.1.1. Characteristics of Production Space Carbon Emissions

For the results of production space carbon emissions, the Jenks Natural Breaks method was used to divide that of 12 cities into five grades, namely, high, higher, medium, lower, and low grade (Figure 3). The same grading method was used in living space, ecological space, and PLE space.
The production space is the main carbon source in Heilongjiang Province. Except for Jixi and Mudanjiang, which showed a downward trend in carbon emissions, the other cities showed upward trends. The city with the highest production space carbon emissions was Daqing, which maintained a high-grade status through the entire study period, followed by Harbin and Qiqihar, which also maintained higher grades, and Heihe, which maintained a low grade. The carbon emission grades of the production spaces in other cities changed, among which Mudanjiang and Jixi changed from medium to lower grades in 2010 and 2020, respectively. Jiamusi was at a lower grade in 2005, changed to a medium grade in 2010 and 2015, and returned to a lower grade in 2020. Yichun was in a low grade from 2005 to 2015 and changed to a lower grade in 2020. Carbon emissions from the production space were gradually transferred from the north of Heilongjiang Province to the southwest.

3.1.2. Characteristics of Living Space Carbon Emissions

The carbon emission capacity of living spaces was slightly weaker than that of production spaces, but the production space of each city increased annually (Figure 4). Harbin has always been the city with the highest living space carbon emissions, high grade, and a high growth rate. With improvements in living standards and increased numbers of motor vehicles, the living space carbon emissions of Harbin may exceed those of the production space in the future, making it the first city in Heilongjiang Province whose living space is the main carbon emissions source. Daqing had the next highest grade. Hegang and Heihe remained in lower grades through the whole study period. Other cities changed to varying degrees. Qiqihar had a higher grade from 2005 to 2015 and transformed to a medium grade in 2020. While it had a higher grade in 2015, Mudanjiang remained at a medium grade in 2005, 2010, and 2020. Qitaihe changed from a lower grade to a low grade. Yichun, Jiamusi, Shuangyashan, Jixi, and Suihua changed between low, lower, and medium grades and finally stabilized at lower grades. Living space carbon emissions changed from a higher overall grade to a higher grade in the southwest and a lower grade in other regions.

3.1.3. Characteristics of Ecological Space Carbon Sinks

Carbon emissions from Heilongjiang Province were much higher than fixation in carbon sinks. The carbon sinks of cities other than Qiqihar, Jixi, Hegang, Daqing, and Mudanjiang showed downward trends, but the magnitude of variation was small (Figure 5). Overall, the carbon sink grade of the ecological space was relatively stable. Heihe contributed the highest carbon sink, followed by Yichun and Harbin, Jixi, and Shuangyashan. Cities with lower-grade sinks included Suihua, Hegang, and Jiamusi. Among these, Hegang was in a lower grade in 2005, transformed into a medium grade in 2010 and 2015, and returned to a lower grade in 2020. The cities with the weakest carbon sink capacities were Qiqihar, Daqing, and Qitaihe; a low grade was maintained during the study period.

3.1.4. Characteristics of PLE Space Total Carbon Emissions

There was little change in high- and low-carbon emission grades and more changes in medium-carbon emission grades (Figure 6). Daqing and Harbin had high total carbon emissions, Qiqihar had a higher grade, Heihe had a low grade, and Qitaihe had a medium grade. Hegang, Jiamusi, and Mudanjiang remained in the lower grades through the study period. The total carbon emissions grade of Jixi changed from medium to lower in 2020. Shuangyashan and Suihua changed from low to medium grades. Yichun changed from a low grade to a lower grade. Cities with higher total carbon emissions in Heilongjiang Province were more concentrated in the southwest, and the total carbon emissions of cities in other regions were relatively low.

3.2. Spatial–Temporal Characteristics of PLE Space Carbon Emissions

According to the results of the Moran’s I test in Table 4, the Moran’s I index of total carbon emissions from 2005 to 2020 fluctuates greatly, first showing a decreasing trend, then increasing, and finally decreasing, from 0.171999 in 2005, 0.148051 in 2010, 0.186807 in 2015 to 0.179023 in 2020. Moran’s I index was greater than 0. However, the Moran’s I index in 2005 and 2010 did not pass the test with the absolute value of Z greater than 1.65 and the p-value less than 0.1, which indicates that the spatial correlation of carbon emissions in PLE spaces in Heilongjiang Province was weak and did not show significant aggregation during 2005–2010. The Moran’s I index in 2015 and 2020 passed the test with the absolute value of Z greater than 1.65 and the p-value less than 0.1, indicating that the spatial correlation is slightly stronger in 2015–2020, and there is a positive spatial agglomeration phenomenon. The overall difference in carbon emissions between the regions gradually widened.

3.3. Structural Features of Carbon Emissions in Segmentation Sectors

To further study the structure of carbon emissions, they were subdivided into six sectors: agricultural carbon emissions (ACE), animal husbandry carbon emissions (AHCE), industrial carbon emissions (ICE), transportation carbon emissions (TCE), household carbon emissions (HCE), and traffic travel carbon emissions (TTCE). The structure of carbon emissions and the contribution rate of each subsector were determined.
ACE and HCE accounted for a relatively small share of the subdivided sectors (Figure 7). The carbon emission structures of the five cities changed. In Harbin, the proportions of ACE, AHCE, and ICE decreased annually, and the largest contributing sector changed from ICE to TTCE. Furthermore, the ICE sector became the second largest contributor to carbon emissions. The carbon emissions structure of Shuangyashan changed from being dominated by ICE and AHCE to ICE alone. In Mudanjiang, the largest contributor changed from ICE to both AHCE and ICE. However, the leading carbon emission sectors in the other eight cities did not change. In Qiqihar, Jixi, Hegang, Daqing, and Qitaihe, ICE was the sector with the largest proportion of carbon emissions. Yichun and Jiamusi maintained the joint dominance of AHCE and ICE. The main sector of carbon emissions in Heihe was AHCE.

3.4. Impact Factors of Carbon Emissions in PLE Spaces

3.4.1. Impact Factors of PLE Space Transfer on Carbon Emissions

From Figure 8, it can be observed that APS and FES occupied a large proportion, and GES, WES, URLS, IPS, and OES occupied relatively small proportions. GES and OES appeared to decrease in space, and IPS and URLS appeared to increase in proportion. In the process of space change, FES in cities with strong carbon sink capacity was reduced. Combined with the carbon emissions of the PLE spaces, APS, IPS, and URLS were the main carbon sources, and construction land was the main carbon source. Figure 9 presents the carbon emissions of APS, IPS, and URLS. APS was the main carbon source in Heihe, and the main industries were agriculture and animal husbandry. IPS and URLS were the main carbon sources in Harbin, which was the only city in Heilongjiang Province where URLS contributed a large proportion of carbon emissions. The differences in carbon emissions contributed by APS, URLS, and IPS in Mudanjiang were not significant. For the other nine cities, including Qiqihar, IPS was the main carbon source. Even in Daqing, where GES had a larger share, FES was the main carbon sink in the city.
From the perspective of land use space, seven cities had increased APS areas: Harbin, Qiqihar, Shuangyashan, Daqing, Yichun, Jiamusi, and Heihe. Harbin, Shuangyashan, Yichun, and Jiamusi had reduced APS carbon emissions, indicating that green agriculture achieved its initial results in these cities. The other five cities with reduced APS areas were Jixi, Hegang, Qitaihe, Mudanjiang, and Suihua, among which Jixi and Hegang had increased APS, indicating that the implemented effects of green agriculture in these two cities were not significant. The IPS area and IPS carbon emissions of the ten cities other than Qitaihe and Mudanjiang increased, indicating that the increase in IPS area promoted carbon emissions. The IPS areas of Qitaihe and Mudanjiang also increased, but their IPS carbon emissions decreased, indicating that industrial structure upgrades in the cities had a mitigating effect on carbon emissions. The URLS area and URLS carbon emissions of each city increased; thus, increasing URLS area will increase carbon emissions.

3.4.2. STIRPAT Model Results

From the results of the STIRPAT model in Table 5, it is clear that 12 models were all better fitted. GDP of each city was positively correlated with carbon emissions, indicating that urban development in Heilongjiang Province is costly in terms of carbon emissions, and that economic development has a relatively significant role in promoting carbon emissions in Shuangyashan and Qitaihe. The proportions of secondary production in Harbin, Daqing, Yichun, Jiamusi, and Suihua were inversely related to carbon emissions, and for these cities, industrial restructuring effectively suppressed carbon emissions. The suppression effect of industrial restructuring in several other cities is not yet clear; however, it had a more obvious promoting effect on carbon emissions in Jixi, Qitaihe, and Mudanjiang. The urbanization rates of Harbin, Qiqihar, Jixi, Daqing, Yichun, Jiamusi, Mudanjiang, and Heihe were inversely related to carbon emissions, indicating that urbanization rate plays a suppressive role in carbon emissions. The total population and rural population of each city decreased, the population gathered in the cities, and the increase in the urbanization rate suppressed the carbon emissions of these cities.
Coal consumption in the cities other than Yichun was positively correlated with carbon emissions. Although coal consumption in Yichun decreased, coal consumption there played a suppressive role in carbon emissions. As for other cities, despite the gradual optimization of the urban energy structure, energy consumption currently relies mainly on coal; coal consumption in each city increased annually, and coal consumption promoted urban carbon emissions. Its promoting effects were more evident in Qiqihar, Shuangyashan, and Qitaihe. Except for Jixi, Qitaihe, and Mudanjiang, increases in the urban road area and road mileage promoted urban carbon emissions. Improvements in road facilities promoted development of traffic travel and transportation, which in turn increased transportation carbon emissions. The negative coefficient of per capita disposable income for Jixi and Qitaihe indicated that the higher the disposable income, the lower the carbon emissions, showing that residents gradually develop low-carbon awareness and can actively choose a low-carbon lifestyle. For the other 10 cities, the increase in per capita disposable income did not awaken low-carbon or environmental protection awareness but rather promoted carbon emissions due to increased consumption. The per capita disposable income of Qiqihar played a relatively significant role in promoting carbon emissions. Except for Qiqihar, Jixi, Qitaihe, and Mudanjiang, the increase in construction land area promoted urban carbon emissions; however, the inhibitory effect in Qiqihar was not significant. The coefficient of the total agricultural machinery power was negative in Jixi, which was the only city with a negative coefficient for this factor, indicating that an increase in the total agricultural machinery power input suppressed carbon emissions in Jixi. The other 11 cities had positive coefficients for the total power of agricultural machinery, indicating that the greater the input of the total power of agricultural machinery, the greater the carbon emissions of the cities. The contribution of the total agricultural machinery power to carbon emissions was relatively significant in the cities of Daqing and Qitaihe.

4. Discussion

For Jixi and Qitaihe, their carbon emissions grades decreased. Combined with the results of the STIRPAT model, the urbanization rate, the area of city paved roads at year-end, the per capita disposable income of urban residents, highway mileage, construction land area, and the total power of agricultural machinery had an inhibitory effect on carbon emissions. These influencing factors can be adjusted in combination with urban development to maintain a state of urban carbon emission reduction. For the other 10 cities, whose carbon emissions grade was maintaining or increasing, the industrial structure, urbanization rate, and coal consumption were factors that inhibited urban emissions. Carbon emissions can be reduced by upgrading industrial structures, optimizing population layouts, and reducing coal consumption.
When analyzing the factors influencing carbon emissions, it was concluded that an increase in the urbanization rate suppressed carbon emissions, which differs from the findings of other scholars [70,71]. When most scholars consider more developed regions as research objects, population and urbanization rates show increasing trends and are positively related to carbon emissions. In contrast, the total, urban, and rural populations of Heilongjiang Province show a trend of loss; therefore, the result is contrary to that of other scholars. However, the same conclusion was obtained in other studies in Heilongjiang Province [72]. The energy structure of Heilongjiang Province has gradually been upgraded, and the energy structure is still dominated by coal, accompanied by a yearly increase in carbon emissions. The energy consumption structure of the city, which is overly dependent on coal, has led to a low degree of energy optimization. The inhibitory effect of the decreasing share of coal consumption on carbon emissions is not yet sufficient to balance the negative growth of carbon emissions [73], and the energy consumption structure of Heilongjiang Province has not yet been significantly improved [74]. In the 14th Five-year Plan of Heilongjiang Province, the proportion of industrial output value must increase to 30% by 2025. For most cities, it will be necessary to build a new industrial system and develop strategic new industries to increase the proportion of industry and reach the targets while ensuring green and low-carbon development to help achieve the “dual carbon goals”.
Heilongjiang Province is in a slower development stage than those of the other provinces, and urban construction and economic development should be synergistic with the dual carbon goals, neither one pursued at the expense of the other. Simultaneously, it is important to ensure the living standards of residents and create a livable environment. This study aims to investigate the grade of carbon emission changes and the influencing factors of each city, and to provide a theoretical basis for low-carbon development. Conducted from the perspective of the three living spaces, this study can be better connected with territorial spatial planning and provide reference for planners and policymakers to promote the “dual carbon goals”.
Cities in Heilongjiang Province have gradually introduced new energy vehicles as they develop. Owing to the limitations of the data sources, this study does not distinguish between the carbon emissions of new energy and traditional vehicles, and the carbon emissions of all motor vehicles are still calculated according to the carbon emissions of traditional vehicles. This may cause carbon emissions from traffic travel to decrease or slowly decrease, affecting the analysis of the carbon emission structure and influencing factors. In future research, the carbon emissions from new energy vehicles and traditional vehicles should be calculated separately. The impact of the introduction of new energy vehicles on carbon emissions from transportation, living space, and the total carbon emissions of the PLE spaces in Heilongjiang Province should be further analyzed. Moreover, owing to the limitations of the data sources, this study does not discuss the carbon emissions generated by waste treatment or the carbon emissions of industrial production processes. In future research, we will continue to investigate these areas, supplement our calculation of carbon emissions, and comprehensively analyze the factors influencing carbon emissions in PLE spaces. The obtained 12 STIRPAT models can reflect the impact and mechanism of socioeconomic factors on the carbon emissions of each city. In addition, carbon emissions in the future can be predicted through the model, to judge the time and amount of carbon peak and carbon neutrality. It can be an effective tool for strategic planning and goal formulation of urban carbon control and emission reduction.

5. Conclusions

In this study, we calculated the PLE space carbon emissions in the prefecture-level cities of Heilongjiang Province, analyzed the characteristics of spatial–temporal features and the structure of carbon emissions in segmentation sectors, and analyzed the influence of land use and socioeconomic factors on carbon emissions, reaching the following conclusions:
  • Carbon emissions in production and living spaces increased yearly, and Daqing was the city with the highest carbon emissions in production space, followed by Harbin and Qiqihar. Harbin was the city with the highest carbon emissions in living space, followed by Daqing. The carbon sinks of all cities were much smaller than the carbon emissions, and except for a small number of cities that showed an increase in carbon sinks, carbon sinks showed decreasing trends. However, the changes were not significant, and the overall pattern of carbon sink capacity was relatively stable.
  • Cities with higher total carbon emissions were concentrated in the southwestern part of Heilongjiang Province, whereas cities in other regions were at a relatively low grade. There was a positive spatial agglomeration phenomenon in 2015 and 2020, and the overall difference in carbon emissions between regions gradually widened.
  • Among the six carbon emission segmentation sectors, the proportions of ACE and HCE were smaller, and urban carbon emissions consisted mainly of AHCE, ICE, TCE, and TTCE. The carbon emission structure of each city was also transformed by the adjustment of urban development and industrial structure.
  • IPS was the main carbon source space for the other nine cities, including Qiqihar. APS was the main carbon source for Heihe. For Harbin, IPS and URLS were the main carbon sources, and Harbin was the only city in Heilongjiang Province where URLS contributed a large proportion of carbon emissions. For Mudanjiang, there were few differences in the carbon emissions contributed by APS, URLS, and IPS. Even in Daqing, where GES had a larger proportion, FES was the main carbon sink in the city.
  • The economy contributed to carbon emissions in all 12 cities. Furthermore, the carbon emissions grade of Jixi and Qitaihe decreased. Urbanization rate, area of city paved roads, per capita disposable income of urban residents, road mileage, construction land area, and total agricultural machinery power inhibited carbon emissions in these cities. For cities with maintaining and increasing carbon emissions grades, industrial structure and coal consumption are factors that could suppress carbon emissions.

Author Contributions

Original draft, methodology, software, visualization, X.W. and R.G.; conceptualization, R.G., X.W., T.W. and C.D.; supervision and visualization, X.W., T.W. and C.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key Research and Development Project, grant number 2018YFC0704705.

Data Availability Statement

Not applicable.

Acknowledgments

We are very grateful to the anonymous reviewers and editors for their helpful reviews and critical comments.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

PLE spacesProduction–living–ecological spacesCCCoal consumption
APSAgricultural Production Space CPRArea of city paved roads
IPSIndustrial Production SpaceDIUPer capita disposable income of urban residents
URLSUrban and Rural Living Space HMHighway mileage
FESForestland Ecological Space CLAConstruction land area
GESGrassland Ecological Space PAMTotal power of agricultural machinery
WESWater Ecological Space ACEAgricultural carbon emissions
OESOther Ecological Space AHCEAnimal husbandry carbon emissions
CECarbon emissionsICEIndustrial carbon emissions
GDPGross regional productTCETransportation carbon emissions
ISIndustrial structureHCEHousehold carbon emissions
URUrbanization rateTTCETraffic travel carbon emissions

References

  1. Sun, Z.R.; Liu, Y.D.; Yu, Y.N. China’s carbon emission peak pre–2030: Exploring multi-scenario optimal low-carbon behaviors for China’s regions. J. Clean. Prod. 2019, 231, 963–979. [Google Scholar] [CrossRef]
  2. Wang, Q.; Li, S.Y.; Li, R.R.; Jiang, F. Underestimated impact of the COVID–19 on carbon emission reduction in developing countries–a novel assessment based on scenario analysis. Environ. Res. 2022, 204, 111990. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, J.H. An empirical study on China’s energy supply-and-demand model considering carbon emission peak constraints in 2030. Engineering 2017, 3, 512–517. [Google Scholar] [CrossRef]
  4. Wu, L.P.; Chen, Y.; Feylizadeh, M.R.; Liu, W.J. Estimation of China’s macro-carbon rebound effect: Method of integrating data envelopment analysis production model and sequential Malmquist-Luenberger index. J. Cleaner Prod. 2018, 198, 1431–1442. [Google Scholar] [CrossRef]
  5. Wang, S.J.; Liu, X.P. China’s city-level energy-related CO2 emissions: Spatiotemporal patterns and driving forces. Appl. Energy 2017, 200, 204–214. [Google Scholar] [CrossRef]
  6. Bai, X.M.; Shi, P.J.; Liu, Y.S. Realizing China’s urban dream. Nature 2014, 509, 158–160. [Google Scholar] [CrossRef]
  7. Liu, J.L.; Liu, Y.S.; Li, Y.R. Classification evaluation and spatial-temporal pattern analysis of production-living-ecological space in China. Acta Geogr. Sin. 2017, 72, 1290–1304. [Google Scholar]
  8. Zhang, H.; He, J.; Yang, X.Y.; Peng, Q.R. Zoning and governance of county ecological space for the increase of carbon sinks: A case study of the Beijing-Tianjin-Hebei region. Planners 2022, 38, 32–40. [Google Scholar]
  9. Zhang, Z.; Li, J.M. Spatial suitability and multi-scenarios for land use: Simulation and policy insights from the production-living-ecological perspective. Land Use Policy 2022, 119, 106219. [Google Scholar] [CrossRef]
  10. Huang, J.C.; Lin, H.X.; Qi, X.X. A literature review on optimization of spatial development pattern based on ecological-production-living space. Proc. Geogr. 2017, 36, 378–391. [Google Scholar]
  11. Quezada, C.R.; Jorquera, F. Urban fabrics to eco-friendly blue–green for urban wetland development. Sustainability 2021, 13, 13745. [Google Scholar] [CrossRef]
  12. Paracchini, M.L.; Pacini, C.; Jones, M.L.M.; Pérez-Soba, M. An aggregation framework to link indicators associated with multifunctional land use to the stakeholder evaluation of policy options. Ecol. Indic. 2011, 11, 71–80. [Google Scholar] [CrossRef]
  13. Tomlinson, S.J.; Dragosits, U.; Levy, P.E.; Thomson, A.M.; Moxley, J. Quantifying gross vs. net agricultural land use change in Great Britain using the Integrated Administration and Control System. Sci. Total Environ. 2018, 628–629, 1234–1248. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, T.; Kazak, J.; Han, Q.; de Vries, B. A framework for path-dependent industrial land transition analysis using vector data. Eur. Plan. Stud. 2019, 27, 1391–1412. [Google Scholar] [CrossRef]
  15. Xie, X.T.; Li, X.S.; Fan, H.P.; He, W.K. Spatial analysis of production-living-ecological functions and zoning method under symbiosis theory of Henan, China. Environ. Sci. Pollut. Res. 2021, 28, 69093–69110. [Google Scholar] [CrossRef]
  16. Fu, C.; Tu, X.Q.; Huang, A. Identification and characterization of production–living–ecological space in a central urban area based on POI data: A case study for Wuhan, China. Sustainability 2021, 13, 7691. [Google Scholar] [CrossRef]
  17. Tu, X.Q.; Fu, C.; Huang, A.; Chen, H.L.; Ding, X. DBSCAN spatial clustering analysis of urban production-living-ecological space based on POI data: A case study of central urban Wuhan, China. Int. J. Environ. Res. Public Health 2022, 19, 5153. [Google Scholar] [CrossRef]
  18. Tian, F.H.; Li, M.Y.; Han, X.L.; Liu, H.; Mo, B.X. A production–living–ecological space model for land-use optimisation: A case study of the core Tumen River region in China. Ecol. Model 2020, 437, 109310. [Google Scholar] [CrossRef]
  19. Xiao, X.Y.; Huang, X.; Jiang, L.L.; Jin, C.X. Empirical study on comparative analysis of dynamic degree differences of land use based on the optimization model. Geocarto Int. 2022, 2025919, 1–18. [Google Scholar] [CrossRef]
  20. Wang, D.; Jiang, D.; Fu, J.Y.; Lin, G.; Zhang, J.L. Comprehensive assessment of production–living–ecological space based on the coupling coordination degree model. Sustainability 2020, 12, 2009. [Google Scholar] [CrossRef]
  21. Xu, Y.J.; Li, T.F.; Xu, L.L. Greenhouse gas inventory model for territorial spatial master plans. Clim. Chang. Res. 2022, 18, 355–365. [Google Scholar]
  22. Wang, Z.Y.; Hu, C.G.; Zhang, Y.; Liu, H.L. Characteristics of carbon emission pattern in Hubei territorial space and its improvement measures. Planners 2022, 38, 24–31. [Google Scholar]
  23. Xu, Y.; Guo, N.; Ru, K.L.; Fan, S.L. Characteristics and optimization strategies of territorial space zone in Fujian Province, China based on carbon neutrality. Chin. J. Appl. Ecol. 2022, 33, 500–508. [Google Scholar]
  24. Zhu, S.L. Low-carbon development of urbanization areas in the context of land spatial development pattern- Case study of Heilongjiang, Guizhou and Guangzhou Province. Clim. Chang. Res. 2016, 12, 132–138. [Google Scholar]
  25. Yan, F.Y.; Gong, Y.S.; Yang, Y.W. Exploration of ideas and methods of county land spatial planning under carbon emission reduction targets. Urban Dev. Stud. 2022, 29, 119–128. [Google Scholar]
  26. Xun, W.H. Territorial space planning path from the perspective of carbon neutrality, Shenyang. Planners 2022, 38, 88–92. [Google Scholar]
  27. Huang, X.J.; Zhang, A.L.; Zhao, R.Q.; Gao, F.; Zhang, H.H. Carbon emission peak, carbon neutrality and territorial spatial planning implementation mechanism. Mod. Urban Res. 2022, 1, 399. [Google Scholar]
  28. Liu, Y.S.; Lin, F.; Sun, Y.; Zhu, L.K. Measurement and obstacle analysis of carbon emission reduction capability of county-level territorial spatial planning: A case study of Yangtze River Delta Ecological Green Integration Demonstration Area. J. Hum. Settl. West China 2023, 38, 9–16. [Google Scholar]
  29. Feng, X.; Li, Y.; Zhang, L.; Xia, C.; Yu, E.; Yang, J. Carbon metabolism in urban production–living–ecological space based on ecological network analysis. Land 2022, 11, 1445. [Google Scholar] [CrossRef]
  30. Guo, L.; Bai, L.F.; Wang, C.Y.; Liu, Y.X.; Yang, Y.Z.; Guo, X.H. Growth mechanism and trend prediction of carbon emission in Chongqing based on production-living-ecological space. Environ. Pollut. Control 2022, 44, 816–823. [Google Scholar]
  31. Ang, B.W.; Su, B. Carbon emission intensity in electricity production: A global analysis. Energy Policy 2016, 94, 56–63. [Google Scholar] [CrossRef]
  32. Hashmi, R.; Alam, K. Dynamic relationship among environmental regulation, innovation, CO2 emissions, population, and economic growth in OECD countries: A panel investigation. J. Clean. Prod. 2019, 231, 1100–1109. [Google Scholar] [CrossRef]
  33. Kaya, Y. Impact of Carbon Dioxide Emission Control on GNP Growth: Interpretation of Proposed Scenarios; Paper Presented to the Energy and Industry Subgroup: Paris, France, 1990. [Google Scholar]
  34. Duro, J.A.; Padilla, E. International inequalities in per capita CO2 emissions: A decomposition methodology by Kaya factors. Energy Econ. 2006, 28, 170–187. [Google Scholar] [CrossRef]
  35. Liu, M.Z.; Liu, X.X. Study on the influencing factors of urban household energy consumption in China based on STIRPAT Model. Resour. Environ. Yangtze Basin 2017, 26, 1111–1122. [Google Scholar]
  36. Yang, S.H. An LMDI-based investigation of the changes in carbon emissions of the transportation sector in the Yangtze River Economic Belt. China Environ. Sci. 2022, 42, 4817–4826. [Google Scholar]
  37. Ma, Z.P.; Li, C.G.; Zhang, J. Understanding urban shrinkage from a regional perspective: Case study of northeast China. J. Urban Plann. Deve. 2020, 146, 05020025. [Google Scholar] [CrossRef]
  38. Guo, R.; Wu, T.; Wu, X.C.; Luigi, S.; Wang, Y.Q. Simulation of urban land expansion under ecological constraints in Harbin-Changchun urban agglomeration, China. Chin. Geogra. Sci. 2022, 32, 438–455. [Google Scholar] [CrossRef]
  39. Geng, B.S.; Zhu, W.R.; Shi, P.L. A functional land use classification for ecological, production and living spaces in the Taihang Mountains. J. Resour. Ecol. 2019, 10, 246–255. [Google Scholar]
  40. Kuang, W.H. Issues regarding on spatial pattern change of national land space and its overall implementation on beautiful vision in new era. Resour. Sci. 2019, 41, 23–32. [Google Scholar]
  41. Fu, J.C.; Zhang, S.L. Functional assessment and coordination characteristics of production, living, ecological function-A case study of Henan Province, China. Int. J. Environ. Res. Public Health 2021, 18, 8051. [Google Scholar] [CrossRef]
  42. Zou, L.L.; Liu, Y.S.; Wang, J.Y.; Yang, Y.Y. An analysis of land use conflict potentials based on ecological-production-living function in the southeast coastal area of China. Ecol. Indic. 2021, 122, 107297. [Google Scholar] [CrossRef]
  43. Tao, Y.Y.; Wang, Q.X. Quantitative recognition and characteristic analysis of production-living-ecological space evolution for five resource-based cities: Zululand, Xuzhou, Lota, Surf Coast and Ruhr. Remote Sens. 2021, 13, 1563. [Google Scholar] [CrossRef]
  44. Li, J.S.; Sun, W.; Li, M.Y.; Meng, L.L. Coupling coordination degree of production, living and ecological spaces and its influencing factors in the Yellow River Basin. J. Clean. Prod. 2021, 298, 126803. [Google Scholar]
  45. Yang, Q.K.; Duan, X.J.; Wang, L.; Jin, Z.F. Land use transformation based on ecological-production-living spaces and associated eco-environment effects: A case study in the Yangtze River Delta. Sci. Geogr. Sin. 2018, 38, 97–106. [Google Scholar]
  46. Kong, D.Y.; Chen, H.G.; Wu, K.S. The evolution of land, eco-environmental effects and its influencing factors in China. J. Nat. Resour. 2021, 36, 1116–1135. [Google Scholar]
  47. Deng, Y.X.; Yang, R. Influence mechanism of production-living-ecological space change in the urbanization process of Guangdong Province, China. Land 2022, 10, 1357. [Google Scholar] [CrossRef]
  48. Zhao, R.Q.; Qin, M.Z. Temporospatial variation of partial carbon source/sinks of farmland ecosystem in coastal, China. J. Ecol. Rural Environ. 2007, 91, 1–6, 11. [Google Scholar]
  49. Tian, Y.; Zhang, J.; He, Y.-Y. Research on spatial-temporal characteristics and driving factor of agricultural carbon emissions in China. J. Integr. Agric. 2014, 13, 1393–1403. [Google Scholar]
  50. Zafeiriou, E.; Azam, M. CO2 emissions and economic performance in EU agriculture: Some evidence from Mediterranean countries. Ecol. Indic. 2017, 81, 104–114. [Google Scholar] [CrossRef]
  51. West, T.O.; Marland, G. A synthesis of carbon sequestration, carbon emissions, and net carbon flux in agriculture: Comparing tillage practices in the United States. Agric. Ecosyst. Environ. 2002, 91, 217–232. [Google Scholar] [CrossRef]
  52. IPCC. 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Available online: www.ipcc.ch/about/preparingreports/ (accessed on 16 October 2020).
  53. Wu, Y.; Xu, B. When will China’s carbon emissions peak? Evidence from judgment criteria and emissions reduction paths. Energy Rep. 2022, 8, 8722–8735. [Google Scholar] [CrossRef]
  54. Zhang, Q.; Tao, X.M.; Yang, P. Research on carbon emissions from metropolis urban passenger transport and countermeasures. China Popul. Resour. Environ. 2012, 22, 35–42. [Google Scholar]
  55. Jia, S.P.; Mao, B.H.; Liu, S.; Sun, Q.P. Calculation and analysis of transportation energy consumption level in China. J. Transp. Syst. Eng. Inf. Technol. 2010, 10, 22–27. [Google Scholar] [CrossRef]
  56. Lin, J.Y.; Meng, F.X.; Cui, S.H.; Yu, Y.; Zhao, S.N. Carbon footprint analysis on urban energy use: A case study of Xiamen, China. Acta Ecol. Sin. 2012, 32, 3782–3794. [Google Scholar]
  57. Ma, J.; Zhou, S.H.; Mitchell, G.; Zhang, J.T. CO2 emission from passenger travel in Guangzhou, China: A small area simulation. Appl. Geogr. 2018, 98, 121–132. [Google Scholar] [CrossRef]
  58. Sun, H.; Liang, H.M.; Chang, X.L.; Cui, Q.C.; Tao, Y. Land use patterns on carbon emission and spatial association in China. Econ. Geogr. 2015, 35, 154–162. [Google Scholar]
  59. Duan, X.N.; Wang, X.K.; Lu, F.; Ouyang, Z.Y. Carbon sequestration and its potential by wetland ecosystems in China. Acta Ecol. Sin. 2008, 2, 463–469. [Google Scholar]
  60. Zhang, Y.; Jiang, B.; Zhao, Y.H.; Zhao, Y. Study on carbon emission effect of urban land use in Northeast China based on land use change. Environ. Sci. Technol. 2022, 45, 209–217. [Google Scholar]
  61. Fang, J.Y.; Guo, Z.D.; Park, S.L.; Chen, A.P. Estimation of terrestrial vegetation carbon sinks in China from 1981 to 2000. Sci. Sin.–Terrae 2007, 37, 804–812. [Google Scholar]
  62. Zhou, J.; Yang, L.; Dong, M.N.; Jiang, L.L.; Wang, Y.; Zang, S.Y. Carbon emission intensity change analysis based on land use in Harbin City from 2004 to 2012. Sci. Geogr. Sin. 2015, 35, 322–327. [Google Scholar]
  63. Zhao, X.F.; Huang, X.J.; Zhang, X.Y.; Zhu, D.M.; Lai, L.; Zhong, T.Y. Application of spatial autocorrelation analysis to COD, SO2 and TSP emissions in Jiangsu Province. Environ. Sci. 2009, 30, 1580–1587. [Google Scholar]
  64. Jaber, A.; Csonka, B. How do land use, built environment and transportation facilities affect bike-sharing trip destinations? Promet. Traffic. Transp. 2023, 35, 119–132. [Google Scholar] [CrossRef]
  65. Xu, B. A spatial analysis of an effective path for low-carbon development of energy-intensive industries. Sustain. Prod. Consum. 2023, 37, 227–241. [Google Scholar] [CrossRef]
  66. Ma, M.d.; Yan, R.; Cai, W.G. An extended STIRPAT model-based methodology for evaluating the driving forces affecting carbon emissions in existing public building sector: Evidence from China in 2000–2015. Nat. Hazard 2017, 89, 741–756. [Google Scholar] [CrossRef]
  67. Anselin, L. Local indicators of spatial association—ISA. Geogr. Anal. 2010, 27, 93–115. [Google Scholar] [CrossRef]
  68. Dietz, T.; Rosa, E.A. Effects of population and affluence on CO2 emissions. Proc. Natl. Acad. Sci. USA 1997, 94, 175–179. [Google Scholar] [CrossRef]
  69. York, R.; Rosa, E.A.; DIETZ, T. STIRPAT, IPAT and ImPACT: Analytic tools for unpacking the driving forces of environmental impacts. Ecol. Econ. 2003, 46, 351–365. [Google Scholar] [CrossRef]
  70. Dietz, T.; Frank, K.A.; Whitley, C.T.; Kelly, J.; Kelly, R. Political influences on greenhouse gas emissions from US states. Proc. Natl. Acad. Sci. USA 2015, 112, 8254–8259. [Google Scholar] [CrossRef]
  71. Qi, X.; Guo, X.N. A comparative study of the impact of urbanization on carbon emissions in emerging economies. Ecol. Econ. 2022, 38, 101–108. [Google Scholar]
  72. Wang, Y.N.; Sun, Y.X.; Li, W.X. Research on the impact of population–related factors on carbon emission: An empirical analysis based on the varying-intercept model. Ecol. Econ. 2017, 33, 19–23. [Google Scholar]
  73. Xing, J.; Song, N.Z.; Li, W.T.; Li, C.Q.Q. Study on carbon intensity reduction target in Heilongjiang Province during the 14th Five-Year Plan period based on STIRPAT model. Environ. Sci. Manag. 2021, 46, 77–81. [Google Scholar]
  74. Ren, J.Q.; Wu, J.H.; Xia, J.Y.; Yin, Y.; Zhou, Z. Primary energy consumption and its structure in Heilongjiang Province. J. Beijing Univ. Chem. Technol. Nat. Sci. Ed. 2018, 45, 89–94. [Google Scholar]
Figure 1. Technical roadmap.
Figure 1. Technical roadmap.
Land 12 01153 g001
Figure 2. Geographical distribution of 13 administrative regions in Heilongjiang.
Figure 2. Geographical distribution of 13 administrative regions in Heilongjiang.
Land 12 01153 g002
Figure 3. Production space carbon emissions in 2005, 2010, 2015, and 2020 (unit: 104 t).
Figure 3. Production space carbon emissions in 2005, 2010, 2015, and 2020 (unit: 104 t).
Land 12 01153 g003
Figure 4. Living space carbon emissions in 2005, 2010, 2015, and 2020 (unit: 104 t).
Figure 4. Living space carbon emissions in 2005, 2010, 2015, and 2020 (unit: 104 t).
Land 12 01153 g004
Figure 5. Ecological space carbon emissions in 2005, 2010, 2015, and 2020 (unit: 104 t).
Figure 5. Ecological space carbon emissions in 2005, 2010, 2015, and 2020 (unit: 104 t).
Land 12 01153 g005
Figure 6. PLE space total carbon emissions in 2005, 2010, 2015, and 2020 (unit: 104 t).
Figure 6. PLE space total carbon emissions in 2005, 2010, 2015, and 2020 (unit: 104 t).
Land 12 01153 g006
Figure 7. Radar diagram of carbon emissions in segmentation sectors (unit: %).
Figure 7. Radar diagram of carbon emissions in segmentation sectors (unit: %).
Land 12 01153 g007
Figure 8. Land use in Heilongjiang Province in 2005, 2010, 2015, and 2020.
Figure 8. Land use in Heilongjiang Province in 2005, 2010, 2015, and 2020.
Land 12 01153 g008
Figure 9. Carbon emissions in APS, IPS, and URLS (unit: 104 t).
Figure 9. Carbon emissions in APS, IPS, and URLS (unit: 104 t).
Land 12 01153 g009
Table 1. The PLE space classification system of Heilongjiang Province.
Table 1. The PLE space classification system of Heilongjiang Province.
Primary ClassificationSecondary ClassificationCorresponding Land Use Type
Production spaceAgricultural Production Space (APS)Paddy field, dry land
Industrial Production Space (IPS)Other construction land
Living space Urban and Rural Living Space (URLS)Urban land and rural settlement
Ecological spaceForestland Ecological Space (FES)Forestland, shrub area, wood land, other forest land
Grassland Ecological Space (GES)High coverage grassland, medium coverage grassland, low coverage grassland
Water Ecological Space (WES)River and canals; lakes; reservoir, pit, and ponds; bottom land
Other Ecological Space (OES)Swampland, bare soil, bare rock
Table 2. Transport passenger and freight conversion coefficients and carbon emission coefficients [54,55].
Table 2. Transport passenger and freight conversion coefficients and carbon emission coefficients [54,55].
Transportation ModeTransport Passenger and Freight Conversion Coefficients (ton/Person)Carbon Emission Coefficients
(kgCO2/t·km)
Railways10.327
Highways0.10.028
Water transportation0.30.053
Civil aviation0.0721.961
Table 3. Carbon emission coefficients of travel mode and average annual mileage [56,57].
Table 3. Carbon emission coefficients of travel mode and average annual mileage [56,57].
Travel Mode Carbon Emission Coefficients/kg/100 kmAverage Annual Mileage/104 km
Private vehicles22.31.5
Bus88.16.5
Taxi28.310
Table 4. Moran’s I index of PLE space carbon emissions in Heilongjiang Province from 2005 to 2020.
Table 4. Moran’s I index of PLE space carbon emissions in Heilongjiang Province from 2005 to 2020.
Moran’s IndexZ-Scorep-Value
20050.1719991.5970660.110251
20100.1480511.4157420.156851
20150.1868071.6886970.091277
20200.1790231.6881870.091375
Table 5. STIRPAT Model Results.
Table 5. STIRPAT Model Results.
CitySTIRPAT Model R2Sig.
Harbin ln C E = 0.130 ln G D P - 0.101 ln I S 0.089 ln U R + 0.114 ln C C + 0.117 ln C P R + 0.129 ln D I U + 0.076 ln H M + 0.131 ln C L A + 0.131 ln P A M 0.9780.013
Qiqihar ln C E = 0.091 ln G D P + 0.118 ln I S 0.039 ln U R + 0.254 ln C C + 0.062 ln C P R + 0.253 ln D I U + 0.206 ln H M 0.007 ln C L A + 0.083 ln P A M 0.9800.022
Jixi ln C E = 0.065 ln G D P + 0.211 ln I S - 0.191 ln U R + 0.137 ln C C 0.222 ln C P R 0.109 ln D I U 0.090 ln H M 0.108 ln C L A 0.125 ln P A M 0.9220.033
Hegang ln C E = 0.094 ln G D P + 0.041 ln I S + 0.041 ln U R + 0.043 ln C C + 0.192 ln C P R + 0.197 ln D I U + 0.191 ln H M + 0.204 ln C L A + 0.060 ln P A M 0.9850.018
Shuang-yashan ln C E = 0.284 ln G D P + 0.074 ln I S + 0.014 ln U R + 0.296 ln C C + 0.047 ln C P R + 0.154 ln D I U + 0.061 ln H M + 0.046 ln C L A + 0.064 ln P A M 0.9650.025
Daqing ln C E = 0.130 ln G D P 0.073 ln I S 0.056 ln U R + 0.163 ln C C + 0.136 ln C P R + 0.132 ln D I U + 0.083 ln H M + 0.129 ln C L A + 0.199 ln P A M 0.9750.018
Yichun ln C E = 0.167 ln G D P 0.017 ln I S 0.021 ln U R 0.038 ln C C + 0.105 ln C P R + 0.190 n D I U + 0.190 ln H M + 0.081 ln C L A + 0.187 ln P A M 0.9820.046
Jiamusi ln C E = 0.143 ln G D P 0.019 ln I S 0.089 ln U R + 0.142 ln C C + 0.147 ln C P R + 0.151 ln D I U + 0.042 ln H M + 0.138 ln C L A + 0.167 ln P A M 0.9640.036
Qitaihe ln C E = 0.300 ln G D P + 0.193 n I S + 0.120 ln U R + 0.321 ln C C 0.171 ln C P R 0.147 ln D I U 0.086 ln H M 0.139 ln C L A + 0.209 ln P A M 0.9680.032
Mudanjiang ln C E = 0.124 ln G D P + 0.382 ln I S 0.152 ln U R + 0.152 ln C C 0.090 ln C P R + 0.068 ln D I U 0.180 ln H M 0.104 ln C L A + 0.141 ln P A M 0.9900.010
Heihe ln C E = 0.113 ln G D P + 0.050 ln I S 0.065 ln U R + 0.192 ln C C + 0.107 ln C P R + 0.102 ln D I U + 0.123 ln H M + 0.125 ln C L A + 0.169 n P A M 0.9140.007
Suihua ln C E = 0.144 ln G D P 0.076 ln I S + 0.010 ln U R + 0.170 ln C C + 0.091 ln C P R + 0.121 ln D I U + 0.144 n H M + 0.143 ln C L A + 0.124 ln P A M 0.9750.025
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, R.; Wu, X.; Wu, T.; Dai, C. Spatial–Temporal Pattern Characteristics and Impact Factors of Carbon Emissions in Production–Living–Ecological Spaces in Heilongjiang Province, China. Land 2023, 12, 1153. https://doi.org/10.3390/land12061153

AMA Style

Guo R, Wu X, Wu T, Dai C. Spatial–Temporal Pattern Characteristics and Impact Factors of Carbon Emissions in Production–Living–Ecological Spaces in Heilongjiang Province, China. Land. 2023; 12(6):1153. https://doi.org/10.3390/land12061153

Chicago/Turabian Style

Guo, Rong, Xiaochen Wu, Tong Wu, and Chao Dai. 2023. "Spatial–Temporal Pattern Characteristics and Impact Factors of Carbon Emissions in Production–Living–Ecological Spaces in Heilongjiang Province, China" Land 12, no. 6: 1153. https://doi.org/10.3390/land12061153

APA Style

Guo, R., Wu, X., Wu, T., & Dai, C. (2023). Spatial–Temporal Pattern Characteristics and Impact Factors of Carbon Emissions in Production–Living–Ecological Spaces in Heilongjiang Province, China. Land, 12(6), 1153. https://doi.org/10.3390/land12061153

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop