Next Article in Journal
Waste Management, Waste Indicators and the Relationship with Sustainable Development Goals (SDGs): A Systematic Literature Review
Previous Article in Journal
Composite Intake Ports Interference Mechanism and Application in Dual Intake Channels of the Opposed-Piston Diesel Engine
Previous Article in Special Issue
Dynamic Calculation Method for Zonal Carbon Emissions in Power Systems Based on the Theory of Production Simulation and Carbon Emission Flow Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Factors Influencing Carbon Emission and Low-Carbon Development Levels in Shandong Province: Method Analysis Based on Improved Random Forest Partial Least Squares Structural Equation Model and Entropy Weight Method

1
School of Mathematics and Statistics, Changchun University, Changchun 130022, China
2
Key Laboratory of Intelligent Rehabilitation and Barrier-Free for the Disabled, Ministry of Education, Changchun University, Changchun 130022, China
3
Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
4
School of Landscape Architecture, Changchun University, Changchun 130022, China
*
Authors to whom correspondence should be addressed.
Sustainability 2024, 16(19), 8488; https://doi.org/10.3390/su16198488
Submission received: 2 September 2024 / Revised: 21 September 2024 / Accepted: 23 September 2024 / Published: 29 September 2024
(This article belongs to the Special Issue Energy Sources, Carbon Emissions and Economic Growth)

Abstract

:
Comprehensively clarifying the influencing factors of carbon emissions is crucial to realizing carbon emission reduction targets in China. To address this issue, this paper develops a four-level carbon emission influencing factor system from six perspectives: population, economy, energy, water resources, main pollutants, and afforestation. To analyze how these factors affect carbon emissions, we propose an improved partial least squares structural equation model (PLS-SEM) based on a random forest (RF), named RF-PLS-SEM. In addition, the entropy weight method (EWM) is employed to evaluate the low-carbon development level according to the results of the RF-PLS-SEM. This paper takes Shandong Province as an example for empirical analysis. The results demonstrate that the improved model significantly improves accuracy from 0.8141 to 0.9220. Moreover, water resources and afforestation have relatively small impacts on carbon emissions. Primary and tertiary industries are negative influencing factors that inhibit the growth of carbon emissions, whereas total energy consumption, the volume of wastewater discharged and of common industrial solid waste are positive and direct influencing factors, and population density is indirect. In particular, this paper explores the important role of fisheries in reducing carbon emissions and discusses the relationship between population aging and carbon emissions. In terms of the level of low-carbon development, the assessment system of carbon emission is constructed from four dimensions, namely, population, economy, energy, and main pollutants, showing weak, basic, and sustainable stages of low-carbon development during the 1997–2012, 2013–2020, and 2021–2022 periods, respectively.

Graphical Abstract

1. Introduction

According to the United Nations Intergovernmental Panel on Climate Change (IPCC), human social and economic activities have accelerated global warming, resulting in higher temperatures, rising sea levels, reduced biodiversity, and more frequent extreme climate events [1,2]. As a major contributor to global warming, CO 2 accounts for approximately two-thirds of global greenhouse gas (GHG) emissions [3]. The concentration of carbon dioxide in the atmosphere has increased from approximately 280 ppm in preindustrial times to 407.8 ppm in 2018 [4], and reached 420 ppm in 2021 [5]. Without adopting any effective action, the world’s average temperature could rise by 1.4–5.8 degrees centigrade over the next 100 years, posing a major threat to global sustainable development [6]. Therefore, the continued increase in carbon emissions has become a global concern. Many countries have adopted mitigation strategies to reduce carbon emissions, including financial incentives or tax subsidies, support funds, insurance premiums, and noneconomic incentives such as regulations, standards, and bans [7]. It is essential for developing countries to improve their carbon reduction efficiency, which was proposed by the Paris Agreement [8].
China is the largest carbon emitter in the world, accounting for nearly 30% of global carbon emissions; thus, China has an obligation to play the leading role in reducing carbon emissions [9]. Moreover, China has pledged to peak its carbon dioxide emissions by approximately 2030 and is striving to achieve carbon neutrality by 2060. Hence, exploring scientific and comprehensive carbon reduction mechanism measures is urgently needed in China.
Since China’s provinces differ greatly in terms of population, economy, industry, etc., it is important to study carbon emission reduction according to the actual situation of China’s regions [10]. The carbon emissions of 34 provinces in China in 2021 are shown in Figure 1.
Figure 1 shows that the Shandong, Hebei, and Jiangsu Provinces and the Inner Mongolia Autonomous Region have the highest carbon emissions among all provinces and regions. Moreover, the Shandong and Jiangsu Provinces are the main sources of carbon emissions in eastern China, and Shandong Province ranks first in China, far exceeding other provinces and regions. In particular, Shandong is a developing province with high energy consumption and rapid urbanization. Owing to its heavy industry, the energy structure is biased toward coal, and many industries have high energy consumption and emissions; hence, research on the carbon peak and carbon neutrality in Shandong Province has attracted much attention. Moreover, establishing an influencing factor system for carbon emissions is crucial to investigate the development of low carbon emissions in Shandong Province.

2. Literature Review

Systematically identifying the influencing factors of carbon emissions and formulating effective strategies are critical to achieving the “30–60” dual-carbon target and realizing the goal of low-carbon sustainable development, such as economic growth, energy structure, and population distribution [11,12,13]. Considering the importance of ecological factors in influencing carbon emissions, He et al. [14] analyzed the relationship between water and energy from a life cycle perspective. In addition, Xian et al. [15] investigated the emissions of SO 2 , NO X , CO, PM 2.5 , and other pollutants via the multiresolution emission inventory model for climate (MEIC), concluding that the pollution abatement policy had a greater inhibitory effect on controlling air pollution and carbon emissions. In addition, trees absorb carbon dioxide through photosynthesis, and forests soak up carbon dioxide in the form of biomass and soil carbon. Teng et al. [16] introduced carbon sequestration by afforestation based on the relaxation measures to study the change in carbon emission efficiency in China. Their findings about utilizing natural ecosystems to sequester carbon and combat climate change are important. Cabon et al. [17] highlighted the importance of considering processes other than photosynthesis to estimate how much carbon trees can sequester, and their findings are important for utilizing natural ecosystems to sequester carbon and combating climate change.
Currently, the main models for studying the factors affecting carbon emissions include the input–output [18,19], Durbin [20,21,22,23], index decomposition analysis [24,25,26], panel threshold [27,28,29], and dynamic spatial econometric models [30]. However, these models do not systematically consider the influencing factors, and lack internal feedback.
The structural equation method (SEM), a popular method for exploring causal relationships, is applied to address barriers [31,32]. The SEM can extract all of the features of hidden information from the observed data while considering the structure and influence within the indicators, and can directly analyze the unmeasured attributes. In recent years, the SEM has played an important role in low-carbon research, such as low-carbon agriculture [33], travel [34], and building industries [35]. These studies are based on questionnaires or visiting survey methods. However, for the issue of carbon emissions, panel data have the advantage of the time dimension, which better reveals the changing patterns between carbon emissions and their influencing factors and avoids subjective judgments. For the first time, Wei et al. [36] utilized panel data and partial least squares structural equation modeling (PLS-SEM) to construct the path relationship between carbon emissions and their multiple influencing factors, assessing the degree of influence of each factor on different regions in China. They considered a two-level indicator system for carbon emission influencing factors and adopted traditional methods to eliminate observable variables with load values less than 0.7, achieving a goodness-of-fit (GOF) of 0.862 for the model. To comprehensively establish a carbon emission influencing factor indicator system and improve model fit, this paper constructs a four-level carbon emission influencing factor indicator system encompassing population, economy, energy, major pollutants, water resources, and afforestation. The random forest (RF) model is utilized to select variables highly correlated with carbon emissions, which are then combined with PLS-SEM, referred to as RF-PLS-SEM. The empirical analysis is conducted in Shandong Province, with the results indicating a significant improvement in the GOF of the refined model.
The innovations of this paper are as follows:
(a)
This paper presents a four-level carbon emission influencing factor system, including six qualitative indicators, such as population, economy, energy, main pollutants, water resources, and afforestation, and 40 quantitative indicators, which is more comprehensive and systematic.
(b)
Compared with traditional PLS-SEM, the improved RF-PLS-SEM substantially enhances the GOF from 0.8141 to 0.9220, and the loading exceeds 0.8. To reveal concealed information within the data, we investigate the mediating variables for the indirect influencing factors via RF-PLS-SEM. In particular, as a negative primary factor, the economic variable is quadratically decomposed via RF-PLS-SEM to explore which factors are important in inhibiting carbon emissions.
(c)
Combining RF-PLS-SEM with the EWM, the carbon emission indicator system is used to calculate the low-carbon development score in Shandong Province. After feature selection and causal analysis by RF-PLS-SEM, the influencing factors are highly coupled in the low-carbon development evaluation model, and the direction of the indicators is determined according to the relationships between the data, to ensure high credibility of the evaluation results.

3. Data and Research Method

3.1. Data Source and Processing

In this study, we initially select the influencing factors that are highly related to carbon emissions from the six dimensions of economy, population, energy, water resources, major pollutants, and afforestation, which include six qualitative indicators and 40 quantitative indicators. The panel data are selected from Shandong Province from 1997 to 2022, and the information and sources of the variables are shown in Table 1.

3.2. Method

3.2.1. Random Forest (RF)

Feature selection by RF is important in building a classification system; it not only filters out important indicators but also reduces the dimensionality of the data [53]. According to the idea of integrated learning, the RF collects multiple decision trees and averages the output of each decision tree to obtain the final output result, which measures the relative importance of each feature for feature selection. The G I N I index is taken as a measure of feature importance in the model, and the following formula calculates the G I N I index of node m:
G I N I m = 1 k = 1 | K | p m k 2 ,
where K represents the number of categories and p m k denotes the proportion of category k in node m. The G I N I index is the purity of the node; the larger the G I N I index, the lower the purity of the node, and the average change in the G I N I index serves as the level of importance of the feature.
The RF is adopted to select a subset of features retained in light of a preset feature threshold after the important features are obtained [54]. The following methods are commonly used for the selection of thresholds:
(1)
Retention of feature indicators with importance scores greater than 0;
(2)
Preservation of the first K feature indicators according to demand;
(3)
Screening out feature indicators with less than 10% feature significance.

3.2.2. Partial Least Squares Structural Equation Modeling (PLS-SEM)

PLS-SEM is a composite-based approach to SEM that uses linear combinations of variables to explain the variance of the target constructs in the structural model [55]. It is one of the most powerful techniques for accounting for the correlation between many measurable and nonmeasurable factors. Suppose that there are J groups of observable variables, and that each group contains p j variables. Thus, each group of observable variables can be expressed as X j = ( x j 1 , x j 2 , , x j p j ) , ( j = 1 , 2 , , J ) . Each group corresponds to a latent variable ξ j , which is assumed to be normalized, i.e., the mean value is 0, and the variance is 1. Hence, each group of observable variables X j and the corresponding latent variables ξ j constitute a measurement model, also known as an external model. The structure model is a description of the different latent variables of causality, usually expressed as Formula (2) [36]
ξ j = i j β i j ξ i + ζ j ,
where ζ j is the random error term, and its mean value of the residuals is 0, which is uncorrelated with ξ j , and where β i j is the coefficient. Formula (2) shows that the interdependence among the latent variables can be considered a causal association model, which is a causal chain with no loops. The causal association model can be represented as a correlation matrix whose dimension is the number of latent variables. If latent variable j explains latent variable i, the element in the matrix takes a value of 1; otherwise, it takes a value of 0; thus, the matrix is also called the internal design matrix.
PLS regression can estimate latent variables in two ways. The first is to calculate latent variables based on the correlation between observable and latent variables, also known as an external estimation. The second way is to evaluate a specific latent variable using other latent variables [32], which is an internal estimation; the result is denoted as Z j , which is shown in Formula (3):
Z j = i , β i j e i j Y i * ,
where “*” denotes the standardization of the estimate, Y i is the external estimate of the other potential variables, and e i j is the internal weight, which is calculated by the following Formula (4):
e i j = s i g n ( r ( Y i , Y j ) ) = 1 ( r ( Y i , Y j ) > 0 ) 0 ( r ( Y i , Y j ) = 0 ) 1 ( r ( Y i , Y j ) < 0 ) ,
where “ S i g n ” refers to the sign function, and r ( Y i , Y j ) denotes the correlation coefficient between the external weight estimates Y i and Y j . The weights for the internal estimates are calculated as follows:
W j = X j T X j 1 X j T Z j ,
where W j is the coefficient after the normal least squares regression of Z j .
The bootstrap method is used for statistical significance testing, but PLS-SEM has no underlying assumption that the data should be normally distributed, meaning parametric significance tests are not applicable for testing whether the coefficients are significant. Instead, the nonparametric bootstrap procedure is applied to PLS-SEM to test the statistical significance of the results, and parameter estimates (e.g., external weights, loadings, and path coefficients) are used to derive standard errors of the estimates. The GOF index can explain the model quality of the measurements and structural models, which are calculated as the geometric mean of the average communality and average R 2 . The GOF is used to evaluate the quality of the model, and a GOF value greater than 0.7 indicates that the model performs well. Loadings reflect the correlation between latent and observed variables, and the observed variables with loadings below 0.7 should be deleted [56].
Cronbach’s alpha (C. alpha) is the coefficient of reliability. When the data have a multidimensional structure, C. alpha is usually low. Dillon-Goldstein’s rho (DG. rho) is used to evaluate the measurement effectiveness of the set of indicators for their corresponding underlying construct. Acceptable values above 0.7 for C. alpha and DG. rho indicate high reliability in the block of interest. Furthermore, the average variance extracted (AVE) is used to assess the convergent validity of the latent variables. The AVE values of all the constructs are greater than 0.5, confirming significant reliability and validity [57].

3.2.3. Entropy Weight Method (EWM)

The EWM is an unbiased weighting method that determines the weights of indicators according to the information provided by the indicators themselves, avoiding the negative impact of subjective factors and making the results more credible [58]. The entropy weight is the parameter that describes the differences in the evaluation objectives. The lower the entropy value, the more information is provided, and the higher the weight.
The first step is data normalization, as follows:
r i j = x i j min { x i j } max { x i j } min { x i j } ,
r i j = max { x i j } x i j max { x i j } min { x i j } ,
where Formula (6) is used to normalize the positive index. The larger the values, the better the model’s efficacy. In contrast, Formula (7) is applied to the negative indicators, and smaller values are more desirable. In addition, r i j is the normalized value, x i j is the original value, and max { x i j } and min { x i j } represent the maximum and the minimum values in the dataset, respectively. The weight of each indicator H i is determined via the entropy method [59], which is defined in Formula (8), and the coefficients f i j and k are calculated by Formulas (9) and (10):
H i = k j = 1 n f i j ln f i j ,
f i j = r i j / j = 1 n r i j ,
k = 1 / ln n ,
where n is the number of evaluated objects. When f i j = 0 , we suppose that f i j ln f i j is also equal to 0.
The entropy weight [60] of each indicator w i is then calculated using Formula (11):
w i = ( 1 H i ) / ( m i = 1 m H i ) ,
where m is the total number of indicators. The carbon emissions score effectively measures low-carbon development by multiplying the entropy weight w i by the dimensionless value r i j . The indicators are summarized by the category layer, obtaining the score S c for each, as shown in Formula (12):
S c = i = 1 k ( w i × r i j ) ,
where S c is the low-carbon development level score of the c-th category, and k is the number of the outermost indicators included in class c.

4. Results and Discussion

4.1. ARIMA Projection of Carbon Emissions in 2022

Because the carbon emission data of Shandong Province in 2022 have not yet been published, this study employs an autoregressive integrated moving average model (ARIMA) to forecast the value based on the time series data in Shandong Province, which shares the national carbon emissions from 1997 to 2021. This approach enhances accuracy compared with directly extrapolating the 2022 data from the historical carbon emissions of Shandong Province. Moreover, this paper applies the grid search method to find the optimal parameters, which are p = 1, q = 0 and d = 0, determined as ARIMA (1, 0, 0) with a prediction value of 0.0819, which is multiplied by the actual national carbon emission data from 2022; thus, we obtain the carbon emission data of Shandong Province in 2022, which are 939.5187 MtCO 2 .

4.2. Analytical Results of PLS-SEM

This paper presents a four-level initial empirical indicator system of carbon emission influencing factors, which is shown in Figure 2.
The indicator system in Figure 2 includes 6 primary indicators of carbon emissions (CE), labeled in orange; 14 secondary indicators, labeled in yellow; 18 tertiary indicators, labeled in green; and 8 quadruple indicators, labeled in blue. The population (POP), economy (ECO), energy (ENE), main pollutants (MP), water resources (WR), and afforestation (AFF) are selected as latent variables for the CE. Initially, it is essential to ascertain the intrinsic relationships between the CE and the six latent variables. In addition, constructing paths between every two variables is crucial, where the path coefficients are calculated according to causality and correlation matrices. The literature shows that the POP influences the ECO [41,61], ENE [62], MP [63,64], WR [65] and AFF [66], and the ENE influences the ECO [67] and MP [68,69]. A PLS-SEM is constructed using the relationship between these latent and observable variables, the loadings of which are shown in Table 2.
Table 2 shows certain instances where the factor loadings between observable variables and their associated latent variables are less than 0.7, for which observable variables should be removed and the model rebuilt. After adjusting for variables, the GOF value is 0.8141, indicating that the indicators selected effectively construct the PLS-SEM, but the model precision can be further optimized.

4.3. Screening Results for RF

To further increase the GOF value and optimize its pathways, this study proposes an improved model named RF-PLS-SEM. First, the RF is used to select indicators in the carbon emission indicator system, where we set the screening threshold at 0.1, meaning indicators with feature importance less than 0.1 are excluded; otherwise, they are retained. The screening results at all levels are shown in Figure 3.
In Figure 3, the horizontal axis represents feature importance, and the vertical axis represents the screened variable. Since the primary indicators are qualitative and feature screening requires actual data, thus CE is directly used to screen the secondary indicators. As illustrated in Figure 3a, the feature importance thresholds of the secondary indicators associated with the AFF and WR do not exceed 0.1. Consequently, the AFF and WR are excluded as primary indicators. In the existing literature, most studies were focused on cumulative afforestation area, which plays a significant role in carbon absorption. However, in this paper, we mainly consider the newly added afforestation area for the current year, which is sourced from the Shandong Statistical Yearbook. Since the photosynthetic capacity of new saplings is relatively low, their carbon sink effect is not significant in the initial stages. Therefore, the annual increase in afforestation area has a relatively low impact on carbon emissions. Furthermore, the rapid urbanization in Shandong Province has led to a decrease in afforestation areas. In a word, afforestation has a minimal impact on carbon emissions in Shandong Province. Moreover, the secondary indicators corresponding to the POP, ECO, ENE, and MP have not been completely eliminated; thus, these four primary indicators are retained. Similarly, the secondary indicators are utilized to screen the tertiary indicators, which are adopted to screen the quadruple indicators. Figure 3b,c show the screening results of the tertiary indicators. Specifically, Figure 3b show that the primary, secondary and primary industries (PI, SI, TI) are retained, and Figure 3c shows that coal (COAC), coke (COKC), crude oil (COC), natural gas (NGC), and electricity consumption (EC) are obtained, whereas diesel oil (DC), gasoline (GC), fuel oil (FC) and kerosene consumption (KC) are rejected. Figure 3d,e show the screening of the indicators at the fourth level, demonstrating the industry screening results for the secondary and tertiary industries, respectively, and none of the variables are excluded. The final indicator system is reconstructed per the above results as shown in Figure 4.
Figure 4 shows that the main factors influencing carbon emission in Shandong Province can be categorized into four aspects: population, energy, economy and main pollutants. After screening, the analysis includes four primary indicators, six secondary indicators, seven tertiary indicators, and eight quaternary indicators.
This paper adopts the explained variance score (EVS), mean squared error (MSE), mean absolute error (MAE) and root mean square error (RMSE) to compare the accuracy of the initial indicator system with that of the renewable system based on the RF. A comparison of the model accuracy is shown in Table 3.
EVS is the number of [0, 1] values; the larger the value, the better its prediction effect. The MSE, MAE, and RMSE are applied to calculate the prediction error. The smaller the value, the better the prediction will be. As shown in Table 3, the accuracy of both the energy and the total indicator systems are improved significantly compared with those of the prescreening indicators.

4.4. Analytical Results of the RF-PLS-SEM

After the variables are screened according to the RF, the GOF value of the model is significantly improved to 0.9220. The structural equation path diagram between the CE and each latent variable constructed is shown in Figure 5.
As shown in Figure 5, POP, ENE, and MP positively impact CE, and the ECO negatively impacts CE. The POP, whose influence coefficient is 0.3712, has some influence on CE, suggesting that the density of population (DP) somewhat contributes to CE. The path coefficient between the ECO and CE indicators is −0.7647, meaning that carbon emissions decrease as the economy grows. The following study addresses the inhibitory effect of the ECO indicator on CE by decomposing the ECO indicator and constructing a new structural equation model in detail.
The ENE, measured by total consumption production (TCP), has the highest positive impact on CE, and the path coefficient is 0.9987, meaning that Shandong is the province that consumes the most energy in China, and high energy consumption is associated with its rapid economic development. In addition, in 2020, coal consumption accounted for 64% of Shandong Province’s energy consumption structure, whereas clean energy consumption accounted for only 7.4%; thus, the issue of carbon emissions caused by energy consumption must be taken seriously.
In addition, the path coefficient of the MP and CE is 0.3666, indicating that the main pollutants have a positive influence on carbon emissions and that the volume of wastewater discharged (VWW) and of common industrial solid waste generated (VCISWG) have positive effects on the CE. The Ministry of Ecology and Environment noted that, in the context of “double carbon”, water environmental management has once again become a key area for reducing carbon emissions [70]. Most recent estimates have shown that global WWTPs directly emit approximately 650 Gg CO2e annually [14]. Moreover, owing to rapid industrialization, large amounts of industrial solid waste have been generated, resulting in different degrees of environmental pollution, such as carbon emissions. Shandong Province is a major industrial and chemical province, with more than 7200 enterprises producing chemical substances, making pollutant management more burdensome. The rational utilization of industrial solid waste can reduce carbon emissions during industrial production. Therefore, the resource utilization of industrial solid waste, such as fly ash, should receive increased attention from government departments [71].
Latent variables are linear combinations of observable variables. The 95% confidence interval is employed to determine whether the path coefficients are significant, and the test results are shown in Table 4.
As shown in Table 4, the P value between POP and CE is 0.1325, indicating that it fails to meet the threshold for statistical significance. However, this means that, as an indirect influencing factor, POP affects carbon emissions through the mediating variables. The rest of the variables pass the significance test, showing that the ENE measured by the TCP, the ECO by the GDP and the per capita GDP (PCG), and the MP measured by the VWW and the VCISWG are direct influencing factors. The loadings between the observed variables and their latent variables are shown in Table 5. The path loadings of the observable variables and their corresponding latent variables are all greater than 0.7, indicating that the model is effective.
Moreover, the model effectiveness test results are displayed in Table 6, showing that C.alpha and DG.rho are greater than 0.7, and that the AVE is greater than 0.5, indicating that the model has better reliability and validity. To summarize, the improved model shows a significant increase in the accuracy with the GOF value from 0.8141 to 0.9220 compared to the traditional PLS-SEM model. At the same time, some of the variable loadings in the original model are below 0.7 in Table 2, whereas all the loadings in the newly established model have been increased to above 0.7 in Table 5.

4.5. Analytical Results of the RF-PLS-SEM after Decomposing Economic Indicators

The ECO is the negative indicator of CE, and the observed variables of the ECO are divided into PCG and GDP, while GDP includes primary industry (PI), secondary industry (SI), and tertiary industry (TI). To study which economic indicators hinder carbon emissions, this paper further analyses the pathways of carbon emission influencing factors, including decomposed economic variables. The RF-PLS-SEM model, which includes detailed economic variables and other influencing factors, has a GOF value of 0.9392, with all the variable loadings surpassing 0.8. The model demonstrates strong reliability and validity, as evidenced in Table 7. The results of the optimized path adjustment model are depicted in Figure 6. In addition, the PI and TI are negative indicators of CE, and both ENE and POP slightly hinder the TI.
The path coefficient of the PI and CE is −0.4603, indicating that primary industry inhibits the growth of carbon emissions. In China, the value added of primary industry mainly comes from agriculture, forestry, fishery, and animal husbandry. Urbanization has reduced the number of fishing grounds, thereby increasing CO 2 emissions [72], and long-term empirical statistics show that for every 1% increase in the negative impact of fishing vessels on aquatic organisms, the carbon intensity increases by 0.55%. The latest data prove that fishery factors play a role in reducing carbon emissions. Because fisheries not only emit fewer carbon emissions than animal husbandry, but also reduce carbon intensity by farming seaweed [73], they play an important role in reducing carbon emissions. Shandong Province is also very rich in fishery resources, with more than 2300 km of coastline and a sea area of over 500,000 square kilometers, and the marine fishery industry is one of the traditionally advantageous industries. Marine pasture is an important factor in absorbing carbon dioxide from the atmosphere [74] and has successfully addressed carbon emissions in China. Shandong Province is rich in various marine biological resources, such as seaweed and marine microorganisms, which can absorb and fix carbon dioxide. Shandong Province has been pushing forward the construction of modern marine ranches and currently has 67 national marine ranching demonstration zones, the most in the country.
The path coefficient of TI and CE is −0.8281, indicating that tertiary industries suppress the growth of carbon emissions. With the rapid development of e-commerce, the technological innovation capacity of modern service industries, such as information technology, has also significantly increased. Improvements to energy efficiency and innovation levels are the main driving force for upgrading production technology [75]. Tertiary industry gradually eliminates the excess production capacity of secondary industries, which helps reduce unnecessary industrial pollution emissions and has a significant inhibitory effect on carbon emission intensity. Research results show that tertiary industries, especially those that involve technological innovation activities, can stimulate more low-carbon and clean technology innovations [76].
The path coefficient of the POP and TI is −0.1629, meaning that when the population increases, the tertiary value added also decreases. Moreover, the tertiary industry has a negative effect on carbon emissions, while population has an indirect effect on carbon emissions, confirming that the population in Shandong Province slightly hinders the development of the tertiary industry, thereby indirectly affecting carbon emissions through the technological innovation capacity of the TI. According to the data obtained from the seventh population census in 2020, 11 provinces in China had an elderly population of 10 million, while Shandong’s exceeded 20 million, making its elderly population the largest in the country, and indicating aging will continue to accelerate. The population of older adults in Shandong Province is growing rapidly, while the proportion of young and middle-aged people is declining annually, which poses challenges to traditional service industries that rely heavily on labor. Many emerging sectors within the tertiary industry require continuous technological innovation, but the population aging trend somewhat restricts enhancing innovation capabilities. This study revealed that the increase in population density in Shandong Province has a slight inhibitory effect on the development of tertiary industries. Therefore, actively developing the “silver economy” and establishing a sound elderly care service system (including health care, tourism, education, etc.) are crucial for promoting economic development within the tertiary industry. This approach transforms the adverse aspects of the increasing elderly population density into a favorable factor that drives economic growth. Furthermore, reducing talent loss in Shandong Province is crucial for achieving the optimal labor allocation in an aging society.
The path coefficient between ENE and TI is −0.1277, indicating that energy consumption not only directly affects carbon emissions, but also slightly inhibits the development of tertiary industries, which is unfavorable for innovating low-carbon clean technologies. The path coefficient between SI and CE is 0.9109, indicating it has a significant effect on CE. China’s economy has long been overly dependent on secondary industries, which have become the economy’s leading industrial sector. The development of secondary industries, especially the manufacturing industry, tends to consume more resources and emit more carbon dioxide than tertiary industries. Therefore, to achieve sustainable economic development, Shandong Province must change its economic development mode as soon as possible and increase the proportion of tertiary industries.

4.6. Low-Carbon Development Level Score

Low-carbon development is a multiobjective issue that seeks not only to reduce greenhouse gas emissions but also to ensure economic growth. Low-carbon development is a complex dynamic system, and the subsystems are interconnected, interact, and constrain one another [77]. The carbon emission indicator system covers the population, economy, energy, and ecology, is comprehensive and systematic, and can be used to evaluate the current status of the low-carbon development level in Shandong Province. The indicators screened are of greater fitness by quantitative analysis, and they are positive or negative indicators according to the path coefficients; thus, the evaluation results are highly convincing. Moreover, there is an inverse relationship between the level of low-carbon development and carbon emissions; that is, if a factor has a positive effect on carbon emissions, it will have a negative effect on the sustainable development score, and vice versa. Influential factors obtained using the RF-PLS-SEM are applied to the EWM to calculate the low-carbon development score in Shandong Province; the results are shown in Table 8.
The weights of the indicators in Table 8 are used to calculate the low-carbon development level score. The plus and minus signs represent the direction of the effect of each indicator on low-carbon development based on the quantitative analysis of the RF-PLS-SEM, respectively. The higher the score, the greater the effect of a low-carbon development level from 1997 to 2021; the results are shown in Figure 7.
In Figure 7, the POP_score, ECO_score, ENE_score, and MP_score represent the scores for population, economy, energy, and main pollutants, respectively; the larger the area of the bar, the higher the indicator score for that year. The ENE_score in Shandong Province was higher and more stable from 1997 to 2000, with a decline to approximately 2001 and a significant drop after 2005. The proportion of POP_score and MP_score of the total score gradually decrease for the low-carbon development level. The ECO_score also shows the opposite state of change compared with the ENE_score. In addition, the level of economic development barely increased before 2005 and then slowly increased. After 2011, the growth rate increased dramatically, and the economy reached the level of low-carbon development in Shandong Province. To adapt to the new development situation in recent years, Shandong Province has been promoting industrial upgrading, continuously eliminating backward production capacity and high energy-consuming production capacity and promoting the economy to better embrace the new era. The overall score trend of the low-carbon development level in Shandong Province is shown in Figure 8.
As shown in the bar chart in Figure 8, the low-carbon development level score of Shandong Province remained stable in the early stage from 1997 to 2000, then began to decline in 2001, and reached its lowest value in 2006. However, it began to increase in 2007, and continued increasing until 2022. The line graph represents the growth rate of the low-carbon development level score, which reveals an initial downward trend, with a negative growth rate from 2000 to 2006 and a negative peak in 2005. Drawing on previous research results, the level of low-carbon development is divided into five stages on average: poor, weak, basic, sustainable and strong, with an overall threshold value of [0, 1] [60]. According to the calculations, the stages of low-carbon development were weak, basic, and sustainable during the 1997–2012, 2013–2020, and 2021–2022 periods, respectively. Moreover, 2021 is an important turning point toward sustainable development, which is the opening year of the “14th Five-Year Plan”. Shandong Province focused on transforming and upgrading traditional energy to achieve positive results. Specifically, investment in high-tech industries grew by 11.6%, crude steel production increased by 12.59 million tonnes, and new and renewable energy generation capacity reached 100.6 billion kilowatt-hours. The quality development of Shandong Province was on the hoof in 2022, and innovation-driven results were significant.

5. Conclusions and Policy Implications

This paper presents a four-level indicator system for the influencing factors of carbon emissions, and proposes an improved RF-PLS-SEM model, combining the EWM to analyze the influencing factors of carbon emissions in Shandong Province. The results demonstrate that population, economy, energy, and main pollutants have greater impacts on carbon emissions, whereas water resources and afforestation have relatively lower impacts. Furthermore, the improved model shows a significant increase in accuracy, with the GOF value increasing from 0.8141 to 0.9220 compared with the traditional PLS-SEM. The specific conclusions and policy recommendations of this paper are as follows:
(1)
The primary and tertiary industries in Shandong Province negatively influence carbon emissions, and the secondary industry significantly contributes to carbon emissions. As one of the larger marine provinces, Shandong Province has abundant marine and fishery resources, and has pushed forward the development of modernized marine pastures, which can help reduce carbon emissions; thus, primary industry plays an important role in reducing carbon emissions. It is essential to develop its advantages in marine and fisheries, promote green development and the upgrading of fisheries, conserve aquatic biological resources, etc., which can also provide some references for other regions with marine resources in Shandong Province. In addition, tertiary industries restrict carbon emissions, but secondary industries significantly promote carbon emissions; thus, accelerating the upgrading of the industrial structure and increasing the proportion of tertiary industries in the national economic system are crucial. Moreover, expanding fiscal subsidies and investment channels for green finance can have a positive effect on realizing the “dual carbon” goals in Shandong Province.
(2)
The population density in Shandong Province has an indirect influence on carbon emissions and slightly inhibits the development of tertiary industries. Shandong Province’s aging population presents both challenges and opportunities for the economic development of tertiary industries. By formulating scientific response strategies and measures, we can fully leverage the market opportunities presented by population aging and actively develop a “silver economy” tailored to the demand characteristics of the elderly population, covering multiple sectors, such as pensions, health care, tourism, and education. By providing diversified products and services to meet the consumption needs of elderly individuals, we can drive the tertiary industry’s transformation, upgrading, and high-quality development.
(3)
Energy, measured by total consumption production, has the greatest positive impact on CE; this not only directly accelerates carbon emissions, but also slightly inhibits the development of tertiary industries, which is adverse for the innovation of low-carbon clean technologies. Shandong Province has vigorously developed its economy with increasing energy consumption, which has led to environmental deterioration; however, this issue has improved with the proposal of sustainable development policies. Shandong Province can adhere to the concept of green development in resource recycling. In addition, promoting a cleaner and low-carbon energy transition and curbing the development of high-energy consumption and high-emission projects can result in a win–win outcome of reducing energy consumption and developing the economy.
(4)
Of the main pollution-influencing factors in Shandong Province, the volume of discharged wastewater and common industrial solid waste generated has a direct and positive influence on carbon emissions. The wastewater treatment industry accounts for the largest share of the environmental protection industry; therefore, policies should be formulated to promote energy-efficient products and equipment and accelerate eliminating old and inefficient equipment. Moreover, the output and trend of solid waste must be assessed by local departments, policies related to solid waste should be formulated, and facilities for solid waste treatment should be constructed in Shandong Province.
(5)
The scores at different stages of low-carbon development in Figure 8 show a trend of steady progress, decline, and then growth. Specifically, the trend was weak, basic, and sustainable in the 1997–2012, 2013–2020, and 2021–2022 periods, respectively, Moreover, the growth rate of the scores from 2000 to 2006 is negative, whereas those of the other years are positive. According to the results of the low-carbon development level assessment, the inhibitory effect of economic factors on carbon emissions must be enhanced in Shandong Province. Moreover, reducing energy consumption and pollution, accelerating industrial upgrading, and promoting green technological innovation are pivotal to achieving the target of reducing carbon emissions as soon as possible.
This paper mainly focuses on the influencing factors of carbon emission, carbon emission mechanisms and the low carbon development level in Shandong Province, and the results will not only provide a scientific basis for Shandong Province to draw up carbon emission policies, but also help Shandong Province to achieve the peak carbon goal by 2030. Furthermore, it will provide some references for other provinces and regions. As we know, the prediction of carbon emission intensity is one of the hottest issues. In the future, we will focus on the prediction of carbon emission intensity according to the improved neural network algorithm. Additionally, we will also consider the impact of extreme weather and policies on carbon reduction strategies.

Author Contributions

Conceptualization, Y.Z. and Y.G.; Data curation, Y.G.; Formal analysis, Y.C.; Investigation, D.Z.; Methodology, Y.G.; Project administration, Y.Z.; Software, Y.G.; Validation, Y.Z. and Y.G.; Visualization, J.M.; Writing—original draft, Y.Z. and Y.G.; Writing—review and editing, Y.C. and J.M. All the authors were informed about each step of manuscript processing, including submission, revision, revision reminders, etc., via email. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the Jilin Provincial Department of Science and Technology (No. 20230101232JC), in part by the National Natural Science Foundation of China (No. 41701054), and, in part, by the Project Grant for Teaching Cases of Graduate Students in Jilin Province (No. JJKH20230100YJG).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors are grateful to the anonymous referees for their careful reading and various corrections, which greatly improved the exposition of the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Brown, P.T.; Hanley, H.; Mahesh, A.; Reed, C.; Strenfel, S.J.; Davis, S.J.; Kochanski, A.K.; Clements, C.B. Climate warming increases extreme daily wildfire growth risk in California. Nature 2023, 621, 760–766. [Google Scholar] [CrossRef] [PubMed]
  2. Wernberg, T.; Thomsen, M.S.; Baum, J.K.; Bishop, M.J.; Bruno, J.F.; Coleman, M.A.; Filbee-Dexter, K.; Gagnon, K.; He, Q.; Murdiyarso, D.; et al. Impacts of climate change on marine foundation species. Annu. Rev. Mar. Sci. 2024, 16, 247–282. [Google Scholar] [CrossRef] [PubMed]
  3. Aslam, B.; Hu, J.; Shahab, S.; Ahmad, A.; Saleem, M.; Shah, S.S.A.; Javed, M.S.; Aslam, M.K.; Hussain, S.; Hassan, M. The nexus of industrialization, GDP per capita and CO2 emission in China. Environ. Technol. Innov. 2021, 23, 101674. [Google Scholar] [CrossRef]
  4. Tian, M.H.; Shen, L.D.; Liu, X.; Bai, Y.N.; Hu, Z.H.; Jin, J.H.; Feng, Y.F.; Liu, Y.; Yang, W.T.; Yang, Y.L.; et al. Response of nitrite-dependent anaerobic methanotrophs to elevated atmospheric CO2 concentration in paddy fields. Sci. Total Environ. 2021, 801, 149785. [Google Scholar] [CrossRef]
  5. Misra, A.; Jha, A. How to combat atmospheric carbon dioxide along with development activities? A mathematical model. Phys. D Nonlinear Phenom. 2023, 454, 133861. [Google Scholar] [CrossRef]
  6. Wang, Y.S.; Gu, J.D. Ecological responses, adaptation and mechanisms of mangrove wetland ecosystem to global climate change and anthropogenic activities. Int. Biodeterior. Biodegrad. 2021, 162, 105248. [Google Scholar] [CrossRef]
  7. Fekete, H.; Kuramochi, T.; Roelfsema, M.; den Elzen, M.; Forsell, N.; Höhne, N.; Luna, L.; Hans, F.; Sterl, S.; Olivier, J.; et al. A review of successful climate change mitigation policies in major emitting economies and the potential of global replication. Renew. Sustain. Energy Rev. 2021, 137, 110602. [Google Scholar] [CrossRef]
  8. Gu, G.; Zhang, W.; Cheng, C. Mitigation effects of global low carbon technology financing and its technological and economic impacts in the context of climate cooperation. J. Clean. Prod. 2022, 381, 135182. [Google Scholar] [CrossRef]
  9. Xu, G.; Dong, H.; Xu, Z.; Bhattarai, N. China can reach carbon neutrality before 2050 by improving economic development quality. Energy 2022, 243, 123087. [Google Scholar] [CrossRef]
  10. Zheng, J.; Mi, Z.; Coffman, D.; Milcheva, S.; Shan, Y.; Guan, D.; Wang, S. Regional development and carbon emissions in China. Energy Econ. 2019, 81, 25–36. [Google Scholar] [CrossRef]
  11. Shen, L.; Wu, Y.; Lou, Y.; Zeng, D.; Shuai, C.; Song, X. What drives the carbon emission in the Chinese cities?—A case of pilot low carbon city of Beijing. J. Clean. Prod. 2018, 174, 343–354. [Google Scholar] [CrossRef]
  12. Chang, K.; Du, Z.; Chen, G.; Zhang, Y.; Sui, L. Panel estimation for the impact factors on carbon dioxide emissions: A new regional classification perspective in China. J. Clean. Prod. 2021, 279, 123637. [Google Scholar] [CrossRef]
  13. Hao, J.; Gao, F.; Fang, X.; Nong, X.; Zhang, Y.; Hong, F. Multi-factor decomposition and multi-scenario prediction decoupling analysis of China’s carbon emission under dual carbon goal. Sci. Total Environ. 2022, 841, 156788. [Google Scholar] [CrossRef]
  14. He, X.; Li, Z.; Xing, C.; Li, Y.; Liu, M.; Gao, X.; Ding, Y.; Lu, L.; Liu, C.; Li, C.; et al. Carbon footprint of a conventional wastewater treatment plant: An analysis of water-energy nexus from life cycle perspective for emission reduction. J. Clean. Prod. 2023, 429, 139562. [Google Scholar] [CrossRef]
  15. Xian, B.; Xu, Y.; Chen, W.; Wang, Y.; Qiu, L. Co-benefits of policies to reduce air pollution and carbon emissions in China. Environ. Impact Assess. Rev. 2024, 104, 107301. [Google Scholar] [CrossRef]
  16. Teng, X.; Liu, F.p.; Chiu, Y.H. The change in energy and carbon emissions efficiency after afforestation in China by applying a modified dynamic SBM model. Energy 2021, 216, 119301. [Google Scholar] [CrossRef]
  17. Green, J.K.; Keenan, T.F. The limits of forest carbon sequestration. Science 2022, 376, 692–693. [Google Scholar] [CrossRef] [PubMed]
  18. Jiang, T.; Yu, Y.; Jahanger, A.; Balsalobre-Lorente, D. Structural emissions reduction of China’s power and heating industry under the goal of “double carbon”: A perspective from input-output analysis. Sustain. Prod. Consum. 2022, 31, 346–356. [Google Scholar] [CrossRef]
  19. Zhou, Z.; Zeng, C.; Li, K.; Yang, Y.; Zhao, K.; Wang, Z. Decomposition of the decoupling between electricity CO2 emissions and economic growth: A production and consumption perspective. Energy 2024, 293, 130644. [Google Scholar] [CrossRef]
  20. Chen, H.; Yi, J.; Chen, A.; Peng, D.; Yang, J. Green technology innovation and CO2 emission in China: Evidence from a spatial-temporal analysis and a nonlinear spatial durbin model. Energy Policy 2023, 172, 113338. [Google Scholar] [CrossRef]
  21. Dabuo, F.T.; Du, J.; Madzikanda, B.; Coulibaly, P.T. Influence of research and development, environmental regulation, and consumption of energy on CO2 emissions in China—Novel spatial Durbin model perspective. Environ. Sci. Pollut. Res. 2023, 30, 29065–29085. [Google Scholar] [CrossRef]
  22. Liu, X.; Niu, Q.; Dong, S.; Zhong, S. How does renewable energy consumption affect carbon emission intensity? Temporal-spatial impact analysis in China. Energy 2023, 284, 128690. [Google Scholar] [CrossRef]
  23. Xu, J.; Wang, J.; Wang, T.; Li, C. Impact of industrial agglomeration on carbon emissions from dairy farming—Empirical analysis based on life cycle assessmsent method and spatial durbin model. J. Clean. Prod. 2023, 406, 137081. [Google Scholar] [CrossRef]
  24. Quan, C.; Cheng, X.; Yu, S.; Ye, X. Analysis on the influencing factors of carbon emission in China’s logistics industry based on LMDI method. Sci. Total Environ. 2020, 734, 138473. [Google Scholar] [CrossRef]
  25. He, Y.; Xing, Y.; Zeng, X.; Ji, Y.; Hou, H.; Zhang, Y.; Zhu, Z. Factors influencing carbon emissions from China’s electricity industry: Analysis using the combination of LMDI and K-means clustering. Environ. Impact Assess. Rev. 2022, 93, 106724. [Google Scholar] [CrossRef]
  26. Wang, M.; Zhu, C.; Cheng, Y.; Du, W.; Dong, S. The influencing factors of carbon emissions in the railway transportation industry based on extended LMDI decomposition method: Evidence from the BRIC countries. Environ. Sci. Pollut. Res. 2023, 30, 15490–15504. [Google Scholar] [CrossRef]
  27. Aye, G.C.; Edoja, P.E. Effect of economic growth on CO2 emission in developing countries: Evidence from a dynamic panel threshold model. Cogent Econ. Financ. 2017, 5, 1379239. [Google Scholar] [CrossRef]
  28. Dong, F.; Li, Y.; Gao, Y.; Zhu, J.; Qin, C.; Zhang, X. Energy transition and carbon neutrality: Exploring the non-linear impact of renewable energy development on carbon emission efficiency in developed countries. Resour. Conserv. Recycl. 2022, 177, 106002. [Google Scholar] [CrossRef]
  29. Bai, L.; Guo, T.; Xu, W.; Liu, Y.; Kuang, M.; Jiang, L. Effects of digital economy on carbon emission intensity in Chinese cities: A life-cycle theory and the application of non-linear spatial panel smooth transition threshold model. Energy Policy 2023, 183, 113792. [Google Scholar] [CrossRef]
  30. Wang, Y.; Liu, J.; Zhao, Z.; Ren, J.; Chen, X. Research on carbon emission reduction effect of China’s regional digital trade under the “double carbon” target—Combination of the regulatory role of industrial agglomeration and carbon emissions trading mechanism. J. Clean. Prod. 2023, 405, 137049. [Google Scholar] [CrossRef]
  31. Mohandes, S.R.; Kineber, A.F.; Abdelkhalek, S.; Kaddoura, K.; Elsayed, M.; Hosseini, M.R.; Zayed, T. Evaluation of the critical factors causing sewer overflows through modeling of structural equations and system dynamics. J. Clean. Prod. 2022, 375, 134035. [Google Scholar] [CrossRef]
  32. Yin, Q.; Wang, Y.; Xu, Z.; Wan, K.; Wang, D. Factors influencing green transformation efficiency in China’s mineral resource-based cities: Method analysis based on IPAT-E and PLS-SEM. J. Clean. Prod. 2022, 330, 129783. [Google Scholar] [CrossRef]
  33. Hou, J.; Hou, B. Farmers’ adoption of low-carbon agriculture in China: An extended theory of the planned behavior model. Sustainability 2019, 11, 1399. [Google Scholar] [CrossRef]
  34. Yin, J.; Shi, S. Analysis of the mediating role of social network embeddedness on low-carbon household behaviour: Evidence from China. J. Clean. Prod. 2019, 234, 858–866. [Google Scholar] [CrossRef]
  35. Zhu, R.; Li, L. SEM-Based Analysis of Carbon Emission Reduction Pathway Study during the Materialization Stage of Prefabricated Buildings: Evidence from Shenyang and Guiyang, China. J. Environ. Public Health 2022, 2022, 9721446. [Google Scholar] [CrossRef]
  36. Wei, Y.; Zhu, X.; Li, Y.; Yao, T.; Tao, Y. Influential factors of national and regional CO2 emission in China based on combined model of DPSIR and PLS-SEM. J. Clean. Prod. 2019, 212, 698–712. [Google Scholar] [CrossRef]
  37. Li, S.; Diao, H.; Wang, L.; Li, L. A complete total-factor CO2 emissions efficiency measure and “2030• 60 CO2 emissions targets” for Shandong Province, China. J. Clean. Prod. 2022, 360, 132230. [Google Scholar] [CrossRef]
  38. Shan, Y.; Guan, D.; Zheng, H.; Ou, J.; Li, Y.; Meng, J.; Mi, Z.; Liu, Z.; Zhang, Q. China CO2 emission accounts 1997–2015. Sci. Data 2018, 5, 170201. [Google Scholar] [CrossRef]
  39. Shan, Y.; Huang, Q.; Guan, D.; Hubacek, K. China CO2 emission accounts 2016–2017. Sci. Data 2020, 7, 54. [Google Scholar] [CrossRef] [PubMed]
  40. Liu, J.; Li, S.; Ji, Q. Regional differences and driving factors analysis of carbon emission intensity from transport sector in China. Energy 2021, 224, 120178. [Google Scholar] [CrossRef]
  41. Cui, Y.; Khan, S.U.; Sauer, J.; Zhao, M. Exploring the spatiotemporal heterogeneity and influencing factors of agricultural carbon footprint and carbon footprint intensity: Embodying carbon sink effect. Sci. Total Environ. 2022, 846, 157507. [Google Scholar] [CrossRef] [PubMed]
  42. Yang, H.; Lu, Z.; Shi, X.; Muhammad, S.; Cao, Y. How well has economic strategy changed CO2 emissions? Evidence from China’s largest emission province. Sci. Total Environ. 2021, 774, 146575. [Google Scholar] [CrossRef] [PubMed]
  43. Lan, B.; Dong, K.; Li, L.; Lei, Y.; Wu, S.; Hua, E.; Sun, R. CO2 emission reduction pathways of iron and steel industry in Shandong based on CO2 emission equity and efficiency. Resour. Policy 2023, 81, 103406. [Google Scholar] [CrossRef]
  44. Murshed, M.; Rashid, S.; Ulucak, R.; Dagar, V.; Rehman, A.; Alvarado, R.; Nathaniel, S.P. Mitigating energy production-based carbon dioxide emissions in Argentina: The roles of renewable energy and economic globalization. Environ. Sci. Pollut. Res. 2022, 29, 16939–16958. [Google Scholar] [CrossRef]
  45. Kemfert, C.; Präger, F.; Braunger, I.; Hoffart, F.M.; Brauers, H. The expansion of natural gas infrastructure puts energy transitions at risk. Nat. Energy 2022, 7, 582–587. [Google Scholar] [CrossRef]
  46. Jiang, Q.; Khattak, S.I.; Rahman, Z.U. Measuring the simultaneous effects of electricity consumption and production on carbon dioxide emissions (CO2e) in China: New evidence from an EKC-based assessment. Energy 2021, 229, 120616. [Google Scholar] [CrossRef]
  47. Afshar, A.; Khosravi, M.; Molajou, A. Assessing adaptability of cyclic and non-cyclic approach to conjunctive use of groundwater and surface water for sustainable management plans under climate change. Water Resour. Manag. 2021, 35, 3463–3479. [Google Scholar] [CrossRef]
  48. Li, M.; Cao, X.; Liu, D.; Fu, Q.; Li, T.; Shang, R. Sustainable management of agricultural water and land resources under changing climate and socio-economic conditions: A multi-dimensional optimization approach. Agric. Water Manag. 2022, 259, 107235. [Google Scholar] [CrossRef]
  49. Wang, X.C.; Klemeš, J.J.; Wang, Y.; Dong, X.; Wei, H.; Xu, Z.; Varbanov, P.S. Water-Energy-Carbon Emissions nexus analysis of China: An environmental input-output model-based approach. Appl. Energy 2020, 261, 114431. [Google Scholar] [CrossRef]
  50. Pahunang, R.R.; Buonerba, A.; Senatore, V.; Oliva, G.; Ouda, M.; Zarra, T.; Muñoz, R.; Puig, S.; Ballesteros, F.C., Jr.; Li, C.W.; et al. Advances in technological control of greenhouse gas emissions from wastewater in the context of circular economy. Sci. Total Environ. 2021, 792, 148479. [Google Scholar] [CrossRef]
  51. Guo, W.; Xi, B.; Huang, C.; Li, J.; Tang, Z.; Li, W.; Ma, C.; Wu, W. Solid waste management in China: Policy and driving factors in 2004–2019. Resour. Conserv. Recycl. 2021, 173, 105727. [Google Scholar] [CrossRef]
  52. Brilli, L.; Carotenuto, F.; Chiesi, M.; Fiorillo, E.; Genesio, L.; Magno, R.; Morabito, M.; Nardino, M.; Zaldei, A.; Gioli, B. An integrated approach to estimate how much urban afforestation can contribute to move towards carbon neutrality. Sci. Total Environ. 2022, 842, 156843. [Google Scholar] [CrossRef] [PubMed]
  53. Kashef, S.; Nezamabadi-pour, H.; Nikpour, B. Multilabel feature selection: A comprehensive review and guiding experiments. Wiley Interdiscip. Rev. Data Min. Knowl. Discov. 2018, 8, e1240. [Google Scholar] [CrossRef]
  54. Yao, R.; Li, J.; Hui, M.; Bai, L.; Wu, Q. Feature selection based on random forest for partial discharges characteristic set. IEEE Access 2020, 8, 159151–159161. [Google Scholar] [CrossRef]
  55. Ali, F.; Rasoolimanesh, S.M.; Sarstedt, M.; Ringle, C.M.; Ryu, K. An assessment of the use of partial least squares structural equation modeling (PLS-SEM) in hospitality research. Int. J. Contemp. Hosp. Manag. 2018, 30, 514–538. [Google Scholar] [CrossRef]
  56. Talebzadehhosseini, S.; Garibay, I. The interaction effects of technological innovation and path-dependent economic growth on countries overall green growth performance. J. Clean. Prod. 2022, 333, 130134. [Google Scholar] [CrossRef]
  57. Dash, G.; Paul, J. CB-SEM vs PLS-SEM methods for research in social sciences and technology forecasting. Technol. Forecast. Soc. Chang. 2021, 173, 121092. [Google Scholar] [CrossRef]
  58. Zou, Z.H.; Yi, Y.; Sun, J.N. Entropy method for determination of weight of evaluating indicators in fuzzy synthetic evaluation for water quality assessment. J. Environ. Sci. 2006, 18, 1020–1023. [Google Scholar] [CrossRef]
  59. Tan, S.; Yang, J.; Yan, J.; Lee, C.; Hashim, H.; Chen, B. A holistic low carbon city indicator framework for sustainable development. Appl. Energy 2017, 185, 1919–1930. [Google Scholar] [CrossRef]
  60. Cunha-Zeri, G.; Guidolini, J.F.; Branco, E.A.; Ometto, J.P. How sustainable is the nitrogen management in Brazil? A sustainability assessment using the Entropy Weight Method. J. Environ. Manag. 2022, 316, 115330. [Google Scholar] [CrossRef]
  61. Khan, I.; Hou, F.; Irfan, M.; Zakari, A.; Le, H.P. Does energy trilemma a driver of economic growth? The roles of energy use, population growth, and financial development. Renew. Sustain. Energy Rev. 2021, 146, 111157. [Google Scholar] [CrossRef]
  62. Rahman, M.M. Exploring the effects of economic growth, population density and international trade on energy consumption and environmental quality in India. Int. J. Energy Sect. Manag. 2020, 14, 1177–1203. [Google Scholar] [CrossRef]
  63. Shi, G.; Lu, X.; Deng, Y.; Urpelainen, J.; Liu, L.C.; Zhang, Z.; Wei, W.; Wang, H. Air pollutant emissions induced by population migration in China. Environ. Sci. Technol. 2020, 54, 6308–6318. [Google Scholar] [CrossRef] [PubMed]
  64. Cao, L.; Li, L.; Wu, Y. How does population structure affect pollutant discharge in China? Evidence from an improved STIRPAT model. Environ. Sci. Pollut. Res. 2021, 28, 2765–2778. [Google Scholar] [CrossRef]
  65. Ma, B.; Tian, G.; Kong, L. Spatial-temporal characteristics of China’s industrial wastewater discharge at different scales. Environ. Sci. Pollut. Res. 2020, 27, 8103–8118. [Google Scholar] [CrossRef]
  66. Wu, X.; Wang, S.; Fu, B.; Liu, J. Spatial variation and influencing factors of the effectiveness of afforestation in China’s Loess Plateau. Sci. Total Environ. 2021, 771, 144904. [Google Scholar] [CrossRef] [PubMed]
  67. Abbasi, K.R.; Shahbaz, M.; Jiao, Z.; Tufail, M. How energy consumption, industrial growth, urbanization, and CO2 emissions affect economic growth in Pakistan? A novel dynamic ARDL simulations approach. Energy 2021, 221, 119793. [Google Scholar] [CrossRef]
  68. Khan, I.; Hou, F.; Le, H.P. The impact of natural resources, energy consumption, and population growth on environmental quality: Fresh evidence from the United States of America. Sci. Total Environ. 2021, 754, 142222. [Google Scholar] [CrossRef]
  69. Sharma, R.; Shahbaz, M.; Kautish, P.; Vo, X.V. Does energy consumption reinforce environmental pollution? Evidence from emerging Asian economies. J. Environ. Manag. 2021, 297, 113272. [Google Scholar] [CrossRef]
  70. Cosgrove, W.J.; Loucks, D.P. Water management: Current and future challenges and research directions. Water Resour. Res. 2015, 51, 4823–4839. [Google Scholar] [CrossRef]
  71. Hou, H.; Su, L.; Guo, D.; Xu, H. Resource utilization of solid waste for the collaborative reduction of pollution and carbon emissions: Case study of fly ash. J. Clean. Prod. 2023, 383, 135449. [Google Scholar] [CrossRef]
  72. Mahmood, H.; Alkhateeb, T.T.Y.; Furqan, M. Oil sector and CO2 emissions in Saudi Arabia: Asymmetry analysis. Palgrave Commun. 2020, 6, 88. [Google Scholar] [CrossRef]
  73. Feng, C.; Ye, G.; Jiang, Q.; Zheng, Y.; Chen, G.; Wu, J.; Feng, X.; Si, Y.; Zeng, J.; Li, P.; et al. The contribution of ocean-based solutions to carbon reduction in China. Sci. Total Environ. 2021, 797, 149168. [Google Scholar] [CrossRef] [PubMed]
  74. Wan, X.; Xiao, S.; Li, Q.; Du, Y. Evolutionary policy of trading of blue carbon produced by marine ranching with media participation and government supervision. Mar. Policy 2021, 124, 104302. [Google Scholar] [CrossRef]
  75. Jin, B.; Han, Y.; Kou, P. Dynamically evaluating the comprehensive efficiency of technological innovation and low-carbon economy in China’s industrial sectors. Socio-Econ. Plan. Sci. 2023, 86, 101480. [Google Scholar] [CrossRef]
  76. Walheer, B. Labor productivity and technology heterogeneity. J. Macroecon. 2021, 68, 103290. [Google Scholar] [CrossRef]
  77. Wang, Y.; Fang, X.; Yin, S.; Chen, W. Low-carbon development quality of cities in China: Evaluation and obstacle analysis. Sustain. Cities Soc. 2021, 64, 102553. [Google Scholar] [CrossRef]
Figure 1. Carbon emissions of 34 provinces in China in 2021. Source of data: the carbon emission data for 30 provinces in China (excluding Hong Kong, Macao, Taiwan, and Tibet Autonomous Region) in 2021 were obtained from the China Emission Accounts and Datasets (CEADs), and the data for Hong Kong, Macao, Taiwan, and the Tibet Autonomous Region were obtained from the Environment and Ecology Bureau, the Environmental Protection Bureau of the Macao Special Administrative Region, the BIOSIS Previews database and the China Tibet News Network, respectively. The green line in the figure represents the Qinling-Huaihe River demarcation line, and the blue line, and the blue lines on the map represent islands or coastlines.
Figure 1. Carbon emissions of 34 provinces in China in 2021. Source of data: the carbon emission data for 30 provinces in China (excluding Hong Kong, Macao, Taiwan, and Tibet Autonomous Region) in 2021 were obtained from the China Emission Accounts and Datasets (CEADs), and the data for Hong Kong, Macao, Taiwan, and the Tibet Autonomous Region were obtained from the Environment and Ecology Bureau, the Environmental Protection Bureau of the Macao Special Administrative Region, the BIOSIS Previews database and the China Tibet News Network, respectively. The green line in the figure represents the Qinling-Huaihe River demarcation line, and the blue line, and the blue lines on the map represent islands or coastlines.
Sustainability 16 08488 g001
Figure 2. Carbon emission empirical indicator system.
Figure 2. Carbon emission empirical indicator system.
Sustainability 16 08488 g002
Figure 3. Feature screening results for RF.
Figure 3. Feature screening results for RF.
Sustainability 16 08488 g003
Figure 4. The system of terminal indicators.
Figure 4. The system of terminal indicators.
Sustainability 16 08488 g004
Figure 5. Latent variables path analysis of carbon emissions. Note: the paths representing positive influencing factors are depicted as solid black lines, whereas the paths for negative influencing factors are shown as red dashed lines.
Figure 5. Latent variables path analysis of carbon emissions. Note: the paths representing positive influencing factors are depicted as solid black lines, whereas the paths for negative influencing factors are shown as red dashed lines.
Sustainability 16 08488 g005
Figure 6. Path analysis of disaggregated economic indicators and carbon emissions and their influencing factors. Note: the paths representing positive influencing factors are depicted as solid black lines, whereas those for negative influencing factors are shown as red dashed lines.
Figure 6. Path analysis of disaggregated economic indicators and carbon emissions and their influencing factors. Note: the paths representing positive influencing factors are depicted as solid black lines, whereas those for negative influencing factors are shown as red dashed lines.
Sustainability 16 08488 g006
Figure 7. Low-carbon development level score by subsystem.
Figure 7. Low-carbon development level score by subsystem.
Sustainability 16 08488 g007
Figure 8. Overall score for the level of low-carbon development. Note: The blue line graph represents the growth rate of the low-carbon development level score.
Figure 8. Overall score for the level of low-carbon development. Note: The blue line graph represents the growth rate of the low-carbon development level score.
Sustainability 16 08488 g008
Table 1. Basic data information.
Table 1. Basic data information.
Indicator (Variable Name)UnitSourceLiterature
Carbon emissions (CE) Mt CO 2 CEADs[3,37,38,39]
Population (POP)   
total population (TP)10,000 personsShandong Statistical Yearbook[40]
density of population (DP)person/sqcmShandong Statistical Yearbook[41]
Economic (ECO)   
per capita GDP (PCG)yuanShandong Statistical Yearbook[3]
gross domestic product (GDP)
primary industry (PI)
100 million
yuan
Shandong Statistical 
Yearbook
[40,42]
secondary industry (SI)
tertiary industry (TI)
agriculture, forestry, animal husbandry and fishery (AFAHF)
industry (IND)
100 million
yuan
National Bureau of 
Statistics
[43]
construction (CON)
wholesale and retail trade (WRT)
hotels and catering services (HCS)
transport storage and postal services (TSPS)
financial intermediation (FI)
real estate (RE)
Energy (ENE)   
total energy production (TEP)
total consumption production (TCP)
10,000 tons of
SCE
Shandong Statistical 
Yearbook; China Energy
Statistical Yearbook
[44]
[42]
 
fuel oil consumption (FOC)
coal consumption (COAC)
coke consumption (COKC)
10,000 tonsShandong Statistical
Yearbook; China Energy
Statistical Yearbook
[42]
crude oil consumption (COC)
kerosene consumption (KC)
diesel oil consumption (DOC)
gasoline  consumption (GC)
natural gas consumption (NGC)10,000 million cu·mShandong Statistical Yearbook; China Energy Statistical Yearbook[45]
electricity consumption (EC)10,000 million kW·hShandong Statistical Yearbook; China Energy Statistical Yearbook[46]
water resources (WR)   
water supply (WS)
surface water resources (SWR)
10,000 million
cu·m
Shandong Statistical 
Yearbook
[47,48]
ground water resources (GWR)
total amount of water resources (TAWR)10,000 million
cu·m
China Water Resources 
Bulletin; Shandong
Water Resources Bullet
[49]
Main pollutants (MP)
volume of waste water discharged (VWW)10,000 tonsShandong Statistical Yearbook[50]
volume of common industrial solid waste generated (VCISWG)
volume of common industrial solid waste utilized (VCISWU)
10,000 tonsShandong Statistical 
Yearbook
[15,51]
volume of sulfur dioxide discharged (VSDD)
volume of particulate emissions (VPE)
Afforestation (AFF)   
total area of afforestation (TAA)
protection forests (PF)
hectareChina Forestry Statistical 
Yearbook
[52]
by-product forests (BF)
fuel forests (FF)
forests for special purpose (FSP)
Note:The primary industry includes AFAHF; the secondary industry includes IND, CON; the tertiary industry includes WRT, HCS, TSPS, FI, RE.
Table 2. The external loadings between latent and observable variables in the PLS-SEM.
Table 2. The external loadings between latent and observable variables in the PLS-SEM.
Variable Relation       LoadingVariable RelationLoading
POP-TP0.9980MP-VCISWU0.9830
POP-DP0.9980MP-VPE−0.4780
ENE-TCP0.6210MP-VSDD−0.8710
ENE-TEC0.9690MP-VWW0.6870
ECO-GDP1.0000WR-TAWR0.7230
ECO-PCG1.0000WR-WS−0.8790
MP-VCISWG0.9830AFF-TAA1.0000
Table 3. Error comparison between the original and the reconstructed indicator system.
Table 3. Error comparison between the original and the reconstructed indicator system.
EVSMSEMAERMSE
ENE-original0.9411 0.01010.07740.1003
ENE-reconstructed0.9891  0.00070.02330.0257
total-original0.8645 0.00990.05650.0997
total-reconstructed0.9902 0.00180.03790.0424
Table 4. Statistical significance test for latent variables.
Table 4. Statistical significance test for latent variables.
VariableEstimateStd. ErrorT-Valuep-Value
Population0.37100.23701.57000.1325
Energy0.99900.23104.33000.0003
Economy−0.76500.1670−4.58000.0002
Main pollutants0.36700.10003.67000.0014
Table 5. The external loadings between latent and observable variables in the RF-PLS-SEM.
Table 5. The external loadings between latent and observable variables in the RF-PLS-SEM.
Variable RelationLoadingVariable RelationLoading
POP-DP1.0000ECO-GDP1.0000
ENE-COEC0.9640ECO-PCG1.0000
ENE-COAEC0.9240MP-VCISW0.9260
ENE-COKEC0.9630MP-VWW0.8770
ENE-NGEC0.9130 
Table 6. Reliability and validity test of the RF-PLS-SEM model for latent variables.
Table 6. Reliability and validity test of the RF-PLS-SEM model for latent variables.
VariableC. alphaDG. rhoAVE
POP1.00001.00001.0000
ENE0.95700.96900.8860
ECO1.00001.00001.0000
MP0.77300.89800.8130
CE1.00001.00001.0000
Table 7. Effectiveness testing of the RF-PLS-SEM model for segmented variables.
Table 7. Effectiveness testing of the RF-PLS-SEM model for segmented variables.
VariableC. alphaDG. rhoAVE
DP1.00001.00001.0000
TEC0.95700.96900.8860
PCG1.00001.00001.0000
SI0.98100.99100.9810
PI1.00001.00001.0000
TI0.99400.99700.9950
POI0.77300.89800.8120
CE1.00001.00001.0000
Table 8. Weights and directions of subsystem indicators.
Table 8. Weights and directions of subsystem indicators.
SubsystemsIndicatorWeightDirection
POPDP0.0498   −
 PCG0.0759   +
ECOAFAHF0.0749   +
 CON0.0056   −
 IND0.0593   −
 WRT0.0978   +
 TSPS0.0774   +
 HCS0.0799   +
 FI0.1280   +
 RE0.0913   +
ENECOAC0.0867   −
 COKC0.0361   −
 COC0.0324   −
 NGC0.0210   −
 EC0.0367   −
MPVWW0.0199   −
VCISWG0.0274   −
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

Zhu, Y.; Guo, Y.; Chen, Y.; Ma, J.; Zhang, D. Factors Influencing Carbon Emission and Low-Carbon Development Levels in Shandong Province: Method Analysis Based on Improved Random Forest Partial Least Squares Structural Equation Model and Entropy Weight Method. Sustainability 2024, 16, 8488. https://doi.org/10.3390/su16198488

AMA Style

Zhu Y, Guo Y, Chen Y, Ma J, Zhang D. Factors Influencing Carbon Emission and Low-Carbon Development Levels in Shandong Province: Method Analysis Based on Improved Random Forest Partial Least Squares Structural Equation Model and Entropy Weight Method. Sustainability. 2024; 16(19):8488. https://doi.org/10.3390/su16198488

Chicago/Turabian Style

Zhu, Yingjie, Yinghui Guo, Yongfa Chen, Jiageng Ma, and Dan Zhang. 2024. "Factors Influencing Carbon Emission and Low-Carbon Development Levels in Shandong Province: Method Analysis Based on Improved Random Forest Partial Least Squares Structural Equation Model and Entropy Weight Method" Sustainability 16, no. 19: 8488. https://doi.org/10.3390/su16198488

APA Style

Zhu, Y., Guo, Y., Chen, Y., Ma, J., & Zhang, D. (2024). Factors Influencing Carbon Emission and Low-Carbon Development Levels in Shandong Province: Method Analysis Based on Improved Random Forest Partial Least Squares Structural Equation Model and Entropy Weight Method. Sustainability, 16(19), 8488. https://doi.org/10.3390/su16198488

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