1. Introduction
Unreasonable development has caused environmental destruction, therefore, it is the general trend to protect the environment and save energy. Global collaboration promotes the rapid development of the world; immoderate use of natural resources, however, has given rise to land degradation, deforestation and biodiversity loss, and so on. Three wastes in industrial production cause soil pollution, water pollution and air pollution, the rapid development of the global secondary industry promotes the burning of fossil fuels in large quantities, and the greenhouse gas released intensifies the greenhouse effect, causing global warming. The melting of the polar ice cap causes the sea level to rise, and some river deltas with low altitude and fertile land are submerged. At the same time, it also causes seawater to pour into the harbor, which pollutes underground water sources and aggravates the salinization of land. We can know from the notice issued by “The State of Global Climate 2020” that the global average temperature in 2020 was about 1.2 degrees Celsius higher than the pre-industrial level. In the face of natural disasters, human beings are extremely helpless. In order to slow down the trend of climate warming, the United Nations adopted the United Nations Framework Convention on Climate Change in New York on 9 May 1992. In 1997, the Kyoto Protocol of the United Nations Framework Convention on Climate Change was successfully formulated, and it provided legally binding quantitative emission reduction and emission limitation targets for developed countries. In December 2017, twenty-nine countries around the world had signed the Joint Statement on Carbon Neutrality. By September 2019, 66 countries had agreed at the United Nations Climate Action Summit that lucid waters and lush mountains are invaluable assets, and had formed the Climate Ambition Alliance. All these measures have accelerated global carbon neutrality. Britain, Sweden, France, Denmark, New Zealand, Hungary and other countries have written the goal of carbon neutrality into their laws. The EU announced that it will become the first “carbon neutral” land in the world in 2050.
China is a big carbon emitter, so it is imperative to promote peak carbon dioxide emissions and achieve carbon neutrality. According to the data of the seventh national census, China’s population has exceeded 1.4 billion, accounting for 21.5% of the world’s total population. Abundant human resources have promoted the rapid development of the secondary industry, which is dominated by manufacturing. China’s economy is developing steadily, among which the traditional manufacturing industry with high energy consumption and high carbon emission is still the main industry in China. In 2019, China’s total carbon emissions reached 10.17 billion tons, accounting for 28% of the global carbon emissions, and China’s industrial carbon emissions accounted for more than 50% of China’s total carbon emissions. Therefore, Zhang et al. [
1] suggested that adjusting the industrial structure and the energy structure have become two obstacles on the road of carbon neutrality in China. Cai et al. [
2] used standard methods to calculate urban carbon dioxide emissions, and established a data set of urban carbon dioxide emissions in China. As the largest developing country in the world and a responsible big country, China passed the Energy Conservation Law of the People’s Republic of China on 1 November 1997 in order to slow down the global warming trend and play a leading role among developing countries. China released the white paper “China’s Policies and Actions to Address Climate Change” in October 2008. In the meantime, China also actively participates in global climate change negotiations, strengthens communication, coordination and cooperation with other countries in the world, and makes contributions to jointly addressing the challenges of climate change and promoting global sustainable development. In September 2020, China proposed at the United Nations General Assembly that carbon dioxide emissions would peak before 2030, and that it strives to achieve carbon neutrality before 2060; In 2021, at the National People’s Congress, peak carbon dioxide emissions and carbon neutrality were written into the government work report for the first time. In the same year, the basic ideas and important measures to realize peak carbon dioxide emissions and carbon neutrality were put forward at the ninth meeting of the Central Committee of Finance and Economics. At the National People’s Congress in 2021, peak carbon dioxide emissions and carbon neutrality were written into the government work report for the first time, and China put forward the basic ideas and important measures to realize peak carbon dioxide emissions and carbon neutrality at the ninth meeting of the Central Committee of Finance and Economics in the same year.
China is a big carbon emitting country and a big energy consumption country. In 2010, the proportion of carbon emissions from coal in primary energy accounted for about 70%. For this reason, Zou et al. [
3] found that the research and development of new energy will greatly promote the realization of carbon neutralization in China. New energy has become the protagonist of the third energy transformation, and will lead the future of carbon neutrality. In [
4,
5,
6,
7], the authors suggest that developing low-carbon cities, optimizing industrial structure, reducing carbon emissions from steel industry, improving carbon emission reduction technology and reducing carbon sequestration cost are important measures for China to realize peak carbon dioxide emissions and carbon neutrality ahead of schedule.
In the field of applied mathematics, there are a few researches on China’s carbon neutrality. Wang et al. [
8] innovatively constructed traditional Markov probability transfer matrix and spatial Markov probability transfer matrix to explore the temporal and spatial evolution of China’s urban carbon emission performance and predict the long-term trend of carbon emission performance. Chen [
9] put forward the energy supply and demand model under two related carbon emission scenarios, namely, China’s planned peak and advanced peak scenarios, and suggested that low carbon would be a basic feature of the change of energy supply and demand structure, and non-fossil energy would replace oil as the second largest energy source. Industrial structure and energy consumption structure all have significant influence on carbon dioxide emissions, especially industrial energy intensity. In [
10], Guo used economic accounting methods to estimate the potential of China’s industrial carbon emission reduction from the perspective of structural emission reduction and intensity emission reduction, and further discussed the influence of industrial internal structure adjustment and energy structure optimization on industrial carbon emission peak and emission reduction potential. According to the data, during the 20 years from 2000 to 2019, the proportion of coal decreased from the original peak of 72.5% to 57.7%, and the natural gas resources increased from 2.2% to 8.4%. Li et al. [
11] used the generalized Weng model to predict the regional natural gas production in China, and the prediction results show that the peak natural gas production will reach 323 billion cubic meters per year in 2036. In [
12], the scholars investigate the relationship between energy consumption, economic growth and carbon dioxide emissions in Pakistan by using the annual time series data from 1965 to 2015. The estimation results of ARDL show that energy consumption and economic growth have both increased CO
2 emissions in Pakistan in the short and long term. In [
13], based on China’s provincial panel data from 2004 to 2016, Liu et al. empirically analyzed the impact of ecological civilization construction on carbon emission intensity by using spatial Durbin model based on STIRPAT model. The above-mentioned scholars only consider one of carbon emission and carbon absorption, but not both. We know that only by considering both of them can we accurately and reasonably forecast the time when China will achieve peak carbon dioxide emissions and carbon neutrality. In this paper, we selected industrial structure and energy consumption structure as the influencing factors of carbon emissions.
As the processes of carbon emission and carbon absorption are time-varying processes, we can describe them by continuous differential equations. Furthermore, considering that carbon emission and carbon absorption are not only related to the current time, but also to the past time, we can use the delayed differential equation model to describe the phenomenon of the dynamic system of carbon emission and carbon absorption more truly and accurately. There is a lot of research work on delayed differential equations in many fields, such as biology, medicine, physics, and so on [
14,
15,
16,
17,
18]. At present, there are few research achievements in describing carbon emission and carbon absorption model by using delayed differential equations, so the purpose of this paper is to use delayed differential equations to describe carbon emission and carbon absorption model.
The motivation of this paper is as follows. Firstly, according to the Chinese government’s goal of achieving peak carbon dioxide emissions by 2030 and carbon neutrality by 2060, we want to make some predictions and analyze whether China can achieve it under the current policy based on the carbon absorption and emission model. If there is some gap between the simulated results of the model and the ideal goal, we can put forward some policy suggestions to achieve China’s peak carbon dioxide emissions carbon neutrality goal by combining the model with the actual situation. Secondly, considering that many scholars have studied the peak carbon dioxide emissions and carbon neutrality in China from different fields, but there are few related studies on the use of delayed differential equations to describe this problem, and we try to establish a carbon absorption and emission model from the perspective of delayed differential equations to analyze the problem, and analyze this problem from different angles to see if we can get new results. Thirdly, this paper establishes a two-dimensional delayed differential equation model of carbon emission and carbon absorption, which is different from models models cited in the literature [
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13]. We focus on analyzing the existence and stability of equilibrium point, and the existence of system bifurcation, and studying the long-term change process of carbon emission and carbon absorption. The models cited in the literature [
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13] include traditional Markov probability transfer matrix and spatial Markov probability transfer matrix, energy supply and demand model, generalized Weng’s model, spatial Durbin model based on STIRPAT model, etc. The above models do not analyze the amount of carbon absorption, but the two-dimensional delayed differential equation proposed by us is not only related to carbon emission, but also to carbon absorption. We also use the knowledge of differential equations to analyze the stability of the model, focusing on the long-term stability of carbon emission and absorption, and so on.
The rest of the content is arranged as follows. In
Section 2, we establish a delayed carbon absorption-emission model based on the carbon emission and carbon absorption. In
Section 3, we analyze the existence and stability of equilibria and the existence of Hopf bifurcation for the model with time delay. In
Section 4, we derive the normal form of the Hopf bifurcation of the above model and analyze the stability of the bifurcating periodic solutions. In
Section 5, we present numerical simulations to verify the correctness of our analysis. Finally, the conclusion is drawn in
Section 6.
2. Mathematical Modeling
In this paper, we consider carbon emission and carbon absorption together, and analyze the problem of carbon neutrality under China’s industrial adjustment. With the rapid development of economy and technology, we assume that carbon emissions and absorption are in a competitive relationship as a whole; this is because in the early stage of China’s economic development, the proportion of traditional industries has increased year by year. In 2007, the added value of China’s secondary industry accounted for 47.6% of the total proportion. At the same time, China’s clean energy development technology is not mature enough, the coal consumption is large and the utilization rate is low. In order to achieve economic growth, traditional high-carbon emission industries are developed, and natural resources are over-exploited, resulting in a significant increase in carbon emissions, immature carbon storage technology, and a corresponding reduction in carbon absorption. As the global climate is gradually warming, the greenhouse effect is obvious year on year, and mankind is facing serious natural disasters. China has gradually realized this great development idea of lucid waters and lush mountains are invaluable assets. In order to implement this correct development idea, our government is actively committed to reducing the coal proportion, improving the energy utilization rate, developing clean energy, shifting from the traditional high-carbon secondary industry to the green and sustainable tertiary industry, reducing carbon emissions, increasing the vegetation coverage, and striving to build a green city. When we only consider the competitive relationship between carbon emission and absorption, we can obtain the following model,
where
represents the annual growth rate of carbon emission,
represents the annual growth rate of carbon absorption,
represents China’s carbon emission amount at time
t,
represents China’s carbon absorption amount at time
t,
means the maximum capacity of carbon emissions and
means the maximum capacity of carbon absorption.
means the competition coefficient coefficient of carbon emissions relative to carbon absorption and
represents the competition coefficient of carbon absorption relative to carbon emission.
We think that adding carbon adsorption saturation term to the model will make the model more realistic. This is because China has a vast territory, diverse climates, wide latitudes, and a large distance from the sea. In addition, the terrain is different, and the terrain types and mountain ranges are diverse, which leads to various combinations of temperature and precipitation and different combinations of water temperatures form different types of forest vegetation. This is because the net carbon absorbed by each vegetation is the same under certain conditions every year. Furthermore, from the technical point of view, we know that the progress of carbon storage technology promotes the increase of carbon absorption, but with the relative backwardness of technology, the carbon storage technology will improve relatively slowly, resulting in the decrease of the change rate of carbon storage.
Considering that the dual wheels of optimizing industrial structure and energy structure proposed by Guo [
10] could make great contributions to national emission reduction, we can use the quadratic function simulated by previous articles to express the relationship between time and carbon emissions. We can assume that the distance between annual carbon emissions and peak carbon emissions represents the speed of carbon emissions, which is reasonable, because the smaller the distance between them, the larger the carbon emissions, and the smaller the slope of the curve. In practice, it shows that as the industrial structure is gradually transferred from the secondary industry to the tertiary industry, the energy structure is also changed from coal-based primary energy to natural gas-based clean energy, and the carbon emission changes slowly.
We know that there will be a series of processes from the transformation of industrial structure and technology research and development to the application of technology in time production, which will take a certain amount of time. If the relationship between carbon emission and carbon absorption in 2022 is simulated, the carbon emission reduction technology in 2022 will increase compared with the carbon emission when the technology is mature, because the carbon emission reduction technology is just successful but immature. Therefore, we should choose the distance from the peak to the carbon emissions before 2022 as the factor that will affect the carbon emissions in 2022. Increased investment from the government in carbon emission reduction technologies and the rapid development of carbon emission reduction technologies will accelerate the transformation of industrial structure and energy structure, as well as increasing the efficiency during the period of putting into use. For this reason, we establish the following model, the descriptions of these parameters are given in
Table 1, and we note that these parameters are all positive,
For convenience, we denote that
then, model (
2) becomes
According to the initial condition of the system (
3), we present a theorem about the nonnegtivity of solution of the system (
3).
Theorem 1. If , the solution of the system (3) with is positive. Proof. First, we prove
when
under the positive initial condition of the system (
3) with
.
We assume that
is not always positive for
and make
be the first time that
. According to the second equation of the system (
3), we can obtain
. The two conclusions we obtain are contradictory. Therefore, under the positive initial condition, the solution
of the system (
3) is positive for
. Then, we prove
when
under the positive initial condition of the system (
3) with
. We assume that
is not always positive for
and make
be the first time that
. According to the first equation of the system (
3), we can obtain
. The two conclusions we reach are contradictory. Therefore, under the positive initial condition, the solution
of the system (
3) with
is also positive for
. □
Remark 1. We prove if , the solution of the system (3) with is positive. It is also not easy for us to prove the solution of the system (3) is positive when . However, according to the numerical simulation of a group of real parameters, we find that the solution of the system (3) is always positive, which is not contradictory to the positivity of the solution of the system (3). Next, we will consider the dynamics phenomena of the system (
3).
4. Normal Form of Hopf Bifurcation
In
Section 3, we have shown that the equilibrium
is unstable when
, and the equilibrium
is locally asymptotically stable when
. To reflect the actual situation, we focus on the delay from technological innovation to practical production. Therefore, we consider the time-delay
as a bifurcation parameter and denote the critical value
, where
is given in (
23). When
, characteristic Equation (
21) has a pair of pure imaginary roots
. Therefore, system (
3) undergoes a Hopf bifurcation at equilibrium
. In this section, we derive the normal form of Hopf bifurcation for the system (
3) by using the multiple time scales method given in [
23,
24].
In order to normalize the delay, we first re-scale the time
t by using
, then translate the equilibrium
to the origin, that is,
For convenience, we still use
x and
y to represent
and
respectively, so Equation (
3) is transformed into:
Equation (
25) can also be written as:
where
and
We let
h be eigenvector corresponding to eigenvalue
of Equation (
26), and
be the eigenvector corresponding to eigenvalue
of adjoint matrix of Equation (
26), satisfying
We treat the delay
as the bifurcation parameter, let
, where
is the Hopf bifurcation critical value,
is perturbation parameter,
is dimensionless scale parameter. Suppose system (
26) undergoes a Hopf bifurcation from the trivial equilibrium at the critical point
, and then, by the MTS method, the solution of (
26) is assumed as follows:
where
and the derivative with regard to
t is transformed into
where
is differential operator, and
We expand
at
by the Taylor expansion, that is,
where
We substitute Formulas (
29)–(
32) into Equation (
26), then comparing the coefficients of
and
on both sides of the equation, respectively. Then, we obtain the following expressions, respectively,
Equation (
33) has the solution with following form,
where
h is given by (
28). Equation (
34) is a linear non-homogeneous equation, and the non-homogeneous equation has a solution if and only if a solvability condition is satisfied. That is, the right-hand side of (
34) should be orthogonal to every solution of the adjoint homogeneous problem. Thus, we solve (
36) into the right part of equation (
34), and the coefficient vector of
is noted as
, by
so we can solve
, namely,
where
with
We solve Equation (
34), as
is a small disturbance coefficient, and we only consider the influence of
on low-order terms, thus, we obtain its solutions with following form:
where
stands for the complex conjugate of the preceding terms, then, we substitute solutions (
38) into (
34), and we obtain
with
Next, substituting solution (
36) and (
38) into (
35), and with the coefficient vector of
noted as
, by solvability condition, we have
Note that
is a disturbance parameter, and
has little influence for small unfolding parameter, thus, we can ignore the
term, then
, can be solved to yield
where
with
Let
, we obtain the normal form of Hopf bifurcation of system (
26) truncated at the cubic order terms:
where
M is given in (
37), and
H given in (
41).
With the polar coordinate
, substituting that expression into Equation (
42), we obtain the amplitude equation of Equation (
42) on the center manifold as
Theorem 4. When , system (26) has periodic solutions. - (1)
If , then the periodic solution reduced on the center manifold is unstable;
- (2)
If , then the periodic solution reduced on the center manifold is stable.
6. Conclusions
In this paper, considering the competitive relationship between carbon emission and carbon absorption, we set up a new two-dimensional differential equation model with time delay to make some predictions and analyze whether China can achieve it under the current policy, which depends on China’s technology research and development level and government policy investment. In practice, the parameters of the model are variable. In order to simplify the problem, the parameters of the model (
3) are constant coefficients. At the same time, the model in this paper does not consider too many factors. The simulation process of the model may be different from the real process. For example, the peak value of the simulation will be higher than the future peak value because we have not fully considered the carbon emission reduction measures, but on the whole, the stability of our model is consistent with the actual one. In addition, we theoretically analyzed the existence and stability of the equilibrium and the existence of Hopf bifurcation, and we also derive the normal form of Hopf bifurcation for the system (
3) by using the multiple time scales method. After that, we selected a set of data for numerical analysis to verify our theoretical analysis results, we find that equilibrium
of system (
3) is locally asymptotically stable when
. When
, the system (
3) has a stable periodic solution near the equilibrium
and we find from
Figure 4 that the optimization and adjustment of industrial structure and energy structure has an important impetus to China’s realization of peak carbon dioxide emissions and carbon neutrality. When the industrial structure is optimized and the energy structure is improved, the time for China to reach peak carbon dioxide emissions will be shortened (see
Figure 4).
Next, our numerical analysis also shows that when the natural growth rate of carbon absorption increases, the time for China to achieve carbon peak carbon dioxide emissions will be shortened and the peak value will also decrease (see from
Figure 5). From
Figure 5, we predict that when the natural growth rate of carbon absorption is 0.3, China will achieve carbon neutrality before 2060. As the natural growth rate of carbon absorption is actually too high, we also need to deepen the industrial reform and optimize the energy structure to reduce the natural growth rate of China’s carbon emissions. Therefore, based on the above research, this paper emphasizes planting trees and improving the level of carbon storage technology to improve the natural growth rate of carbon absorption and carbon emission reduction technology, and improving the development and application technology of new energy to achieve in-depth industrial structure adjustment and energy structure optimization. The above measures are of great significance to China’s realization of peak carbon dioxide emissions by 2030 and carbon neutrality by 2060.