Forecasting of Carbon Emission in China Based on Gradient Boosting Decision Tree Optimized by Modiﬁed Whale Optimization Algorithm

: As the global temperature continues to rise, people have become increasingly concerned about global climate change. In order to help China to effectively develop a carbon peak target completion plan, this paper proposes a carbon emission prediction model based on the improved whale algorithm-optimized gradient boosting decision tree, which combines four optimization methods and signiﬁcantly improves the prediction accuracy. This paper uses historical data to verify the superiority of the gradient boosting tree prediction model optimized by the improved whale algorithm. In addition, this study also predicted the carbon emission values of China from 2020 to 2035 and compared them with the target values, concluding that China can accomplish the relevant target values, which suggests that this research has practical implications for China’s future carbon emission reduction policies.


Introduction
As the world's largest developing country, China is also the world's largest carbon emitter [1]. Since entering the 21st century, China's carbon emissions have been growing rapidly, and its share in the world is also on the rise. In 2021, the terms carbon peak and carbon neutral were first written in the government's work report and became hot topics discussed by delegates. Xi Jinping, President of the People's Republic of China, stressed that the "14th Five-Year Plan" is a key period and a window period for a peak in carbon emissions, and an important plan needs to be made for the next five years to reduce carbon emissions, laying out a clear "construction plan" for the carbon emission peak in the 14th Five-Year Plan. In order to smoothly reduce carbon emissions, the accurate prediction of carbon emissions has become an important issue.
There are many kinds of prediction methods, such as the Bayes method [2][3][4], rule induction method [5,6] and KNN (K nearest neighbor) algorithm [7][8][9][10]. At present, the most widely used prediction methods are various algorithms in the field of artificial intelligence (AI), which are widely used because of their high computing speed and high computing precision.
Decision tree representation is one of the most widely used logic methods. It deduces the classification rules of decision tree representation from a group of unordered and irregular cases. Tso et al. used decision trees to predict electrical energy consumption [11]; Peng et al. used gradient driven regression trees to predict household energy consumption [12]. Decision tree analysis has been used to determine the influence factors of residents' carbon dioxide emissions under different travel modes [13]. One of the greatest advantages of the Many improved algorithms have been developed based on WOA. To overcome the shortcomings of low convergence accuracy and convergence speed, Yan et al. used logistic mapping to initialize population localization and improved the algorithm using inertial weighting [57]. To address the problems of the WOA algorithm's uneven population distribution, slow convergence speed and the tendency to fall into local optima, He et al. used Sobol sequences to initialize the population, a stochastic learning strategy to increase population diversity and a Corsi variance boosting algorithm to jump out of local optima to optimize the algorithm [58]. Similarly, Liu et al. introduced an optimal neighborhood perturbation strategy to solve the phenomenon of the premature algorithm [59]. Lin et al. used chaotic Tent mapping to determine the initial population position of a random generation algorithm in order to make the population distribution more uniform and speed up the convergence of the algorithm [60].
We found that various scholars use more than one optimization approach to optimize the WOA algorithm. Based on this, this paper combines multiple optimization methods to optimize the WOA algorithm in order to better improve the algorithm's optimization finding ability.
At present, the GBDT algorithm has not been widely used in carbon emission prediction. Based on this, this paper establishes the improved whale algorithm-optimized GBDT to predict China's carbon emissions.
Some of the innovations of the current study are as follows: (1) To improve the global search ability and local search ability, this paper uses four methods: the compound chaotic map, nonlinear convergence factor, local domain perturbation and reverse learning. In addition, this paper compares the optimization performance of the algorithm and proves that the optimized whale algorithm has strong optimization ability; (2) Based on the impact of the policy, this paper evaluates China's 2030 Carbon Peak Target, provides a scientific basis for carbon reduction policy making and puts forward relevant suggestions.

Standard Whale Optimization Algorithm
The specific steps of whale optimization algorithm (WOA) are discussed below.

Prey Encirclement Stage
Humpback whales can identify the location of prey and surround them, and WOA assumes that the current optimal candidate solution is the target prey, while the other search agents update their positions by "target prey". The formula is as follows: where i is the current number of iterations; X * (i) is the optimal candidate solution; X(i) is the current position vector; || is the absolute value sign; is the element-by-element multiplication; A and C are random vectors between [0, 1]; A is the convergence factor, which decreases linearly from 2 to 0.

Bubble-Net Attack
Humpback whales spit bubbles to attack their prey, spiraling upwards to gradually reduce the envelope to obtain food, divided into two main mechanisms of contracting the envelope and spiraling to update the position to achieve mathematical visibility. The mathematical model can be expressed as follows: (a) Contraction of the envelope: encirclement shrinkage of the prey by means of the control parameter a of the coefficient variable A; (b) Spiral position renewal: after encircling the prey, the whale captures it in a spiral movement, which is mathematically modelled as follows: where D p = |X * (t) − X(t)| is the distance between the whale and its prey, b is the logarithmic spiral-shaped constant and l is a random number between [−1, 1]. Humpback whales attack their prey in bubble nets in a synchronous manner with both constricted encirclement and spiral position updates. In the mathematical model, the position information is updated using the same probability and can be expressed as where p is a random number on [0, 1].

Searching for Prey
In addition to the bubble-net hunting strategy, humpback whales need to search for prey randomly, and the process is mathematically modeled as follows: where X rand represents a randomly selected position vector in the current population.

Composite Chaotic Mapping
The initial populations of most current intelligent optimization algorithms are randomly generated in the search space, and the quality of the initialized populations has a great impact on the efficiency of the optimization algorithm. WOA has the disadvantages of reducing the diversity of individuals in the late iterations and easily falling into local optima when solving problems. Chaotic operators can enhance the diversity of individuals in a population by not repeatedly traversing all states within a certain range. This paper uses Tent and Chebyshev mappings to form a composite chaotic mapping as a new search method [61].
Tent mapping is defined as: Chebyshev mapping is defined as: When the order k ≥ 2, the interval of Chebyshev mapping is x ∈ [−1, 1], and the non-linear factor that generates chaos is the square operation.
By embedding the standard Tent mapping into Chebyshev mapping and changing its chaotic attractor, an improved composite chaotic mapping is obtained, whose equation is x n+1 = (cos(k · (arccos(1 − |2x n |)))) (9) When the order k ≥ 2, the complex chaotic mapping interval is x ∈ [−1, 1]. The non-linear factors that generate chaos are taken as absolute values and squared operations.

Non-Linear Convergence Factor
In order to more effectively balance the global exploration and local exploitation performance of the algorithm, a combined convergence factor strategy is proposed in this paper. This is specified as follows: α * decreases nonlinearly in the [0, 2] range. The convergence factor changes more in the first period for fast convergence and less in the later period for an improvement in the local search efficiency.

Local Neighbourhood Perturbation
There are two modes of search in the whale algorithm, namely, search around optimal individuals and search around random individuals, both of which make it difficult for the whale algorithm to operate outside the local optimum. If it converges quickly, it is likely to fall into a local optimum.
We set up a local perturbation optimization [62] to prevent the algorithm from falling into a local optimum.
If the optimal particle in the perturbation group is better adapted than the existing global optimal particle, then the best particle in the perturbation group is better than the existing global best particle, then global best particle is assigned to the global best particle and the global best particle is used to replace the worst particle, so as to retain the existing global optimum while skipping the local optimum. The global optimal particle is assigned to the global optimal particle, and the global optimal particle replaces the worst particle.
If no better particle appears, the local perturbation is considered, the local perturbation fails and the existing particles are not changed.

Reverse Learning Strategy
In order to enhance the population diversity of individual whales in the search optimization process, this paper introduces a backward learning mechanism in the iterative search process of the whale algorithm. The backward learning mechanism is used to help the algorithm escape the influence of the local optimum, which is conducive to the global convergence of the algorithm.
Reverse learning is defined as follows: Supposing Q = (x 1 , x 2 , x 3 , · · ·, x n ) is a feasible solution in an N-dimensional search space and satisfies x 1 , . . , n . Based on this, we can see, with respect to the definition of backward learning, that there is only one corresponding backward solution for each solution in the search space. We introduce a backward learning mechanism into the whale algorithm to achieve the goal of maintaining population diversity while avoiding becoming trapped in local optimization.

Gradient Boosting Tree
GBDT was proposed by Friedman, and it mainly solves the optimization problem of general loss functions. The core idea is to fit the residuals of the previous round of base learners by the negative gradient of the loss function, so that the residual estimate of each round gradually decreases. Thus, the output of each round of the base learner gradually approximates the true value. Fitting in the negative gradient direction ensures that the loss function decreases as quickly as possible in each round of training, accelerating convergence to a local or global optimum solution. The corresponding equation is as follows: where x is the input sample, h is the categorical regression tree, w is the parameters of the regression tree and a denotes the weights. Solving for the optimal model (16) is equivalent to minimizing the loss function, which is: Overall, GBDT is a process of fitting the residuals and superimposing them on F. In this process, the residuals become smaller and the loss (loss function) approaches the minimum.
As GBDT is modelled with a large number of hyperparameters, the training model is not guaranteed to be optimal. In this paper, MWOA is used to optimise the four parameters of min_samples_split, min_samples_leaf, min_weight and max_depth in the GBDT model.
Explanations of these parameters are given below: Min_samples_split specifies the minimum number of samples needed to split an internal node (not a leaf node).
Min_samples_leaf specifies the minimum number of samples to be included in each leaf node.
Min_weight_fraction_leaf specifies the minimum weight factor of the samples in a leaf node.
Max_depth specifies the maximum depth of each regression estimate. The maximum depth limits the number of nodes in the tree and this parameter is adjusted for the best performance.

MWOA-Based GBDT Prediction Model
Data are collected, and dimensionless processing is performed on the data. The gray correlation degree is applied to screen the influencing factors.
(1) Input data selection Collect the data and do dimensionless processing on the data. Apply the gray correlation degree to screen the influencing factors.
(2) GBDT prediction model based on MWOA On the basis of considering the influencing factors of carbon emission prediction, the MWOA-GBDT model is used to derive the energy consumption prediction results. The specific steps are as follows: Step 1: The parameters are set and calculated to initialize the algorithm population based on Equation (9); Step 2: The population fitness is calculated and ranked, and the current best position is recorded; Step 3: The parameters are updated. If p < 0.5, one can move to Step 4, and, conversely, the searching individual attacks the prey in a spiral motion; Step 4: If |A| ≤ 1, the whale swims toward the optimal individual. Conversely, the whale swims toward the random individual; Step 5: The magnitude of the fitness value between the current individual and the candidate optimal solution is used to determine whether to replace the candidate solution; Step 6: After solving the reverse solution of the updated optimal position, the solution with the larger fitness value between them is selected as the final updated position; Step 7: Partial disturbance is performed. The newly generated position is compared with the current optimal position and the individual fitness is compared to decide whether to replace the current optimal individual; Sustainability 2021, 13, 12302 7 of 18 Step 10: It is determined whether the end condition is reached. If it has not been reached, one must return to Step 2 to continue the iteration.
Finally, the parameters obtained from the optimization in MWOA as previously described are assigned to GBDT to complete the prediction.
The whole model is shown in Figure 1.
whale swims toward the random individual; Step 5: The magnitude of the fitness value between the current individual and the candidate optimal solution is used to determine whether to replace the candidate solution; Step 6: After solving the reverse solution of the updated optimal position, the solution with the larger fitness value between them is selected as the final updated position; Step 7: Partial disturbance is performed. The newly generated position is compared with the current optimal position and the individual fitness is compared to decide whether to replace the current optimal individual; Step 10: It is determined whether the end condition is reached. If it has not been reached, one must return to Step 2 to continue the iteration.
Finally, the parameters obtained from the optimization in MWOA as previously described are assigned to GBDT to complete the prediction.
The whole model is shown in Figure 1.

Confirmation of Input Values for the Prediction Model
There are many factors influencing carbon emissions, and we chose to select China's total GDP, residential consumption level, total import and export, industrial structure, urbanization rate and total energy consumption from 1990 to 2019 as initial values based on some past studies and experience (data from the World Bank, China Carbon Accounting Database and BP Energy Statistical Yearbook). In order to determine the main influencing factors more accurately and exclude the influence of secondary factors, we used gray correlation analysis to determine the final influencing factors.
The analysis results of this paper are shown in the following table.
In Table 1, we selected influences with a correlation of 0.7 or higher as input values for the prediction model. Based on the above, we chose GDP, population consumption level, total import and export, urbanization rate and total energy consumption from 1990 to 2019 as the data set for this experiment.

MWOA-GBDT-Based Carbon Emission Forecasting in China
In this paper, the K-fold cross-validation method was used to prevent the overfitting of the model on the test set. The data from 1990 to 2010 were divided into training and validation sets, the verification set at a ratio of 4:1 and the remaining sample data from 2010 to 2019 were used as the test set. Each time, we used the five groups of sample data as the training data and the remaining group of sample data as the verification set for a total of five times, in order to use the average accuracy to verify the model. As the sample data in this study are small, we adopted five-fold cross-validation for every model, in which each group had four samples. The training set had 16 samples and the validation set has 4 samples.
The simulations were run on an Inter(R) Core(TM) CUP: i7-9750H, 16 GB RAM, 2.60 GHz main frequency computer, and the program was implemented using Python programming. The population size was set to 30.
The algorithm parameters are set as shown in Table 2, where the first decimal of each parameter is rounded to the nearest whole number.  [1,10] The predicted results are shown in the graph below. Relative error (RE) refers to the absolute error caused by the measurement and the ratio of the measured (agreed) true value multiplied by 100% of the resulting value, expressed as a percentage. The closer the RE value is to 0, the better the variables of the equation explain y and the better the model fits the data. The closer the RE value is to 1, the worse the model fits.
In general, the relative error is more reflective of the degree of confidence in the measurement. The actual values for the period of 1990-2019 are compared with the predicted values predicted using MWOA-GBDT, as shown in Table 3 and Figure 2.  As can be seen from the above graphs, the prediction curves for energy consumption formed using the MWOA-GBDT model show an excellent fit to the actual curves and the predictions are quite significant, with the relative error at all prediction points not even exceeding 1%. This proves that MWOA-GBDT has superior prediction performance. As can be seen from the above graphs, the prediction curves for energy consumption formed using the MWOA-GBDT model show an excellent fit to the actual curves and the predictions are quite significant, with the relative error at all prediction points not even exceeding 1%. This proves that MWOA-GBDT has superior prediction performance.

Comparison of Prediction Results between Models
In order to further demonstrate the reliability of MWOA-GBDT prediction models, we selected five models (WOA-GBDT, PSO-GBDT, GWO-GBDT, BA-GBDT and GBDT) to predict the data via the same samples, and the prediction results and prediction errors of each model are analyzed and studied in this paper. Figure 3 shows the actual carbon emission values and the predicted carbon emission values of each model.    Figure 4, the MWOA-GBDT model has the smallest prediction error, followed by PSO-GBDT.
In addition to visually comparing the magnitude and distribution of errors, MAPE (mean absolute percentage error), RMSE (root mean square error) and MAE (mean absolute error) were used to compare the forecasting accuracy. The metrics selected by us are formulated as follows: The prediction results of the MWOA-GBDT model fit the actual carbon emission values best, and the fitting effect of PSO-GBDT is also very significant. The prediction results of GWO-GBDT have large fluctuation, and the prediction results of GBDT have an obvious error. The box plots in the figure show the minimum, first quartile, median, third quartile and maximum values of the relative errors in the six models. Figures 4 and 5 compare the relative error magnitudes of the models. As seen in Figure 4, the MWOA-GBDT model has the smallest prediction error, followed by PSO-GBDT.  In addition to visually comparing the magnitude and distribution of errors, MAPE (mean absolute percentage error), RMSE (root mean square error) and MAE (mean absolute error) were used to compare the forecasting accuracy. The metrics selected by us are formulated as follows:  MAPE is an indicator that represents the average deviation of the predicted outcome from the true outcome. MAPE takes into account not only the error between the predicted and true values, but also the ratio between the error and the true value. In addition, MAPE can also evaluate the goodness of the model.
The indicator RMSE represents the deviation of the observed value from the true value.
The indicator MAR represents the average of the absolute values of the deviations of the predicted values from the true values.
The smaller the value of these three indicators, the better the accuracy of the prediction model in describing the experimental data.
As we can see from Table 4, among the three indicators MAPE, RMSE and MAR, MWOA-GBDT has the smallest value, which proves that the model accuracy is relatively optimal. It is closely followed by PSO-GBDT, and the unoptimized GBDT has the relatively worst prediction accuracy.

China's Carbon Emission Projections for 2020-2035
China has proposed to achieve the "carbon peak" target by 2030, and, therefore, this paper uses MWOA-GBDT to project carbon emission values for 2020-2035 to help China formulate relevant policies. In addition to the carbon peak target, at the UN Climate Ambition Summit, President Xi announced that, by 2030, China's carbon emissions intensity will drop by more than 65% compared to 2005.
Based on this, this paper uses MWOA-GBDT to project China's carbon emissions and carbon intensity data for 2020-2035 in conjunction with existing policies and data to determine whether China will be able to efficiently meet its carbon peaking targets.

Simulation of Influencing Factors
Before predicting China's carbon emissions, it is necessary to forecast the influencing factors in advance. Since China has only proposed carbon neutrality and carbon peaking targets come 2020, a considerable number of important policies have been or are being introduced. The traditional gray forecasting method is based on the past situation, and it does not match with the impact factors that we need. Therefore, it is necessary to combine the factors with China's policy planning to determine the trend of some influencing factors, and then combine them with the gray forecasting method to obtain more accurate forecasting results. In this paper, all of the indicators are set according to China's current plans and targets.

GDP
Due to the COVID-19 pandemic, the GDP growth rate for 2020 is 2.3% per year. At China's two sessions, China set its expected GDP growth target for 2021 at 6% or more. The Fifth Plenary Session of the 19th Central Committee proposed that GDP per capita should reach the level of medium developed countries by 2035, implying that China's GDP growth rate needs to be maintained in a reasonable range over the next 15 years. In the process of shifting from high growth to high quality development, China's total economic base is increasing in size, while the potential economic growth rate is expected to slow down due to declining labor productivity growth and demographic changes, and then the GDP growth rate will gradually decline over time.

Consumption Level of the Population
Under the influence of the COVID-19 pandemic, residents' consumption levels have dropped; however, under the joint promotion of policy support factors, such as a stable employment situation and recovery of business activities due to economic recovery and continuous increase in livelihood protection, residents' disposable income has improved in tandem with the economic situation. In the long run, the 14th Five-Year Plan period will be a critical period for China to move towards becoming a high-income country. The Fifth Plenary Session of the 19th CPC Central Committee clearly put forward that raising people's income level is a very important issue. We will strive to reach a new level of livelihood and well-being, with the growth of people's income and economic growth essentially synchronized and the distribution structure significantly improved. In the future, China's economic upward growth pulling force will mainly rely on consumer demand, which mainly requires raising residents' income, and better income distribution can improve overall consumer demand and raise residents' consumption level.

Total Imports and Exports
China's Ministry of Commerce in July in the "14th Five-Year Plan" for business development proposed that, by 2025, China's total imports and exports of goods should increase from USD 4.65 trillion in 2020 to USD 5.1 trillion, with an average annual growth rate of 2%. Officials have indicated that China will expand imports of high-quality consumer goods, advanced technology, important equipment and energy resources, and improve the quality of exports. The global target of "double carbon" will also greatly enhance the global competitiveness of China's manufacturing industry and improve the country's export capacity. However, the continued evolution of the overseas COVID-19 pandemic has exacerbated the complexity of the foreign trade situation and put pressure on future import and export trade.

Industrial Structure
At present, the industrial base of China's service industry trade development is still relatively weak, and the share of the service industry in the value added of GDP is 54.5%, which is lower than the share of world service industry in the value added of the world GDP of 67%. Under the double carbon target, China will prevent the continued development of high energy-consuming and high-emission industries and accelerate the development of strategic emerging industries, high-tech industries and modern service industries. As a result, the share of the service sector will increase until it reaches the world average.

Urbanization Rate
The integration of peak carbon and carbon neutrality with urbanization is proceeding in an orderly manner. Globally, the countries and regions that have already achieved peak carbon and that are promoting carbon neutrality are generally developed countries and regions. The common feature of these developed countries and regions is that they have already achieved urbanization. In a certain sense, they are building low-carbon and zerocarbon cities in the post-urbanization and post-industrialization period. The difference in China is that urbanization has not yet been completed and is still in the process of shaping changes, which provides a scenario where carbon peaking and carbon neutrality can be fully utilized, and each city needs to target its urbanization process according to its own carbon pressure. During the 2030-2035 period, China's resident population urbanization rate is also expected to approach or reach its peak and enter a relatively stable state.

Total Energy Consumption
In the future, because China is still a developing country, China's total energy consumption will remain on the increase for a longer period of time, and the growth rate of energy consumption is expected to show a gradual slowdown. The government's work report proposes that energy consumption per unit of GDP will be reduced by 13.5% during the 14th Five-Year Plan period. The National Development and Reform Commission (NDRC) has proposed a target of controlling total energy consumption to within six billion tonnes of standard coal by 2030. With the continuous development of high-tech carbon reduction technologies and policy subsidies, the total energy consumption will continue to decrease in later years. At present, China is entering the late stages of industrialization and the overall energy consumption demand has seen a post-peak decline. The World and China Energy Outlook 2050, published in 2016, states that China's energy consumption will peak at around 2035 under the influence of accelerated economic restructuring and policies to control total energy consumption. Based on the current new mission targets and new energy policies, total energy consumption is expected to peak by 2030.
Based on the above, we combined policy regulations and the gray forecasting method to obtain the data on influencing factors in China for the period of 2020-2035.

Forecast Results
Based on our previous settings of the corresponding influencing factors, we used the constructed MWOA-GBDT model to forecast the carbon emission values of China for the period of 2020-2035. The prediction results are shown in Figure 6 and Table 5 and the comparison results of the predictions are shown in Table 6. Based on the above, we combined policy regulations and the gray forecasting method to obtain the data on influencing factors in China for the period of 2020-2035.

Forecast Results
Based on our previous settings of the corresponding influencing factors, we used the constructed MWOA-GBDT model to forecast the carbon emission values of China for the period of 2020-2035. The prediction results are shown in Figure 6 and Table 5 and the comparison results of the predictions are shown in Table 6.     Carbon intensity refers to the amount of carbon dioxide emissions per unit of GNP growth. This indicator is used to measure the relationship between a country's economy and its carbon emissions. If a country's economic growth is accompanied by a decline in carbon dioxide emissions per unit of GNP, then the country has achieved a low carbon development model. The predicted results are shown in Table 6. To show more clearly the trend of carbon emission values and carbon emission intensity, a line graph is drawn in this paper, as shown in Figure 6.
As can be seen in Figure 6, the carbon emission values show an increasing and then decreasing trend, peaking in 2028, while the carbon emission intensity shows a decreasing trend year by year. The comparison of the results in Table 6 shows that both policy target values for China can be achieved under the existing policies, while the peak carbon target will be achieved earlier in 2028.

Conclusions and Recommendations
Based on the results in the previous subsection, the following conclusions can be drawn in this paper.
(1) Total energy consumption, urbanization level and industrial structure are the top three factors with the highest correlation. Therefore, relevant suggestions can be made in terms of low carbon consumption to further reduce carbon emissions; (2) By comparing the performance of the optimization algorithm and the prediction model, the MWOA-GBDT model constructed in this paper has an excellent prediction capability; (3) Due to the superiority of MWOA-GBDT in error comparisons, the prediction results of carbon emission consumption have practical significance. According to the prediction results of 2020-2035, China can achieve the carbon emission-related target in 2030 under the existing policies.
To ensure that China's carbon targets are met on time, this paper also makes the following observations. As China is currently in a period of high-quality development, many high-energy consumption and high-carbon projects are still being deployed, and the Fifth Plenary Session of the 19th CPC Central Committee has put forward a clear requirement to "formulate a carbon emission peaking action plan by 2030". It is necessary to make some relevant recommendations.
The recommendations in this paper are as follows: (1) Acceleration of the adjustment of industrial structure and implementation of industry transformation.
China's high-carbon industries account for more than 70% of total energy consumption. From the screening of the above-mentioned influencing factors, there is a close connection between industrial structure and energy consumption. Strictly controlling the new production capacity of high energy-consuming industries and promoting the transformation and upgrading of traditional high energy-consuming industries are the key measures to reduce the value of carbon emissions.
(2) Acceleration of the construction of a perfect carbon market system.
Building a national carbon emission trading market is an important institutional innovation to control greenhouse gas emissions and promote green and low-carbon de-velopment using market mechanisms, and it has become one of the core policy tools to achieve China's peak CO2 emission target and implement carbon-neutral vision.
(3) Building a clean, low-carbon, safe and efficient energy system. Carbon dioxide emissions from energy consumption account for nearly 90% of China's total carbon dioxide emissions, and total energy consumption is highly correlated with the value of carbon emissions. Promoting the effective allocation and utilization of energy resources and effectively improving the energy-saving technology and energy management capacity of key energy-using enterprises are important tasks.
(4) Promotion of a carbon-labeling system for products and low-carbon consumption by the public.
On the basis of setting mandatory product energy efficiency standards, the state can establish a product energy efficiency rating and certification system to compel enterprises to express the effectiveness and level of product energy efficiency levels on product packaging, which can effectively guide consumers' demand and purchase behavior.
In summary, we forecasted China's future carbon emissions for the sake of its future policy development and made some relevant recommendations based on the current situation.
Author Contributions: X.C.: conceptualization, methodology and writing-original draft preparation; S.E. and D.N.: data curation, analysis, and writing-original draft preparation; B.C. and J.F.: resources, investigation, and validation. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement:
The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to privacy reasons.

Conflicts of Interest:
The authors declare no conflict of interest.