Introduction

Carbon dioxide is the main greenhouse gas that causes global warming, and the burning of fossil fuels such as coal and petroleum will cause a large amount of carbon dioxide emissions. In the latest research, a large number of scholars have studied various energy issues, such as the industrial solar energy (Wang et al. 2020), clean energy (Wang 2015), coal (Shou et al. 2020), renewable and hydro energy (Utkucan Sahin 2020), and natural gas (Wang and Li 2020).

Energy consumption and carbon dioxide emissions are closely related, so scholars will pay attention to energy and carbon emissions at the same time. Pao et al. (2012) not only studied energy consumption, they also predicted carbon dioxide emissions and economic growth. There are also some scholars specializing in the study of carbon dioxide emissions in certain regions, such as the carbon dioxide emissions of the BRICS and OECD countries (Wu et al. 2020; Saidi and Omri 2020). The above scholars have done a lot of research, but there are some shortcomings. They both studied energy consumption and carbon dioxide emissions separately, and failed to establish a connection between them. The model in this paper makes up for this deficiency.

Carbon dioxide emissions are affected by many factors. Therefore, scholars have studied carbon dioxide emissions from different perspectives. For example, Lin and Agyeman (2020) used energy-related carbon dioxide emission dynamics to assess low-carbon development in sub-Saharan Africa. Tobelmann and Wendler (2020) study the impact of environmental innovation on carbon dioxide emissions. Li et al. (2019) used panel data to analyze the impact of modernization on carbon dioxide emissions. There are also studies on carbon dioxide emissions from the perspectives of urbanization (Zhou et al. 2019) and power generation (Eberle and Heath 2020). From the energy point of view, this paper establishes the connection between energy consumption and carbon dioxide emissions through the non-equigap GM(1,1) model, and analyzes the impact of energy changes on carbon dioxide emissions.

This paper is divided into 6 parts. The literature review in the “Literature review” section is studied from the perspective of grey accumulator operator and non-equigap grey model. The “Non-equigap GM(1,1) model with conformable fractional accumulation” section establishes a non-equigap GM(1,1) model with conformable fractional accumulation. The “Validity verification” section uses two cases to verify the rationality of the novel proposed model. The “Analyze the relationship between energy consumption and carbon dioxide emissions” section analyzes the relationship between energy consumption and carbon dioxide emissions in 53 countries and regions. The “Conclusion” section is a summary of the full text.

Literature review

The difference between the grey model and other models is that the grey model contains an accumulation operator. Since Professor Deng Julong proposed a grey model with first-order accumulation, many types of accumulation methods have emerged. For example, the fractional accumulation operator proposed by Wu et al. (2013), the conformable accumulation operator proposed by Ma et al. (2020), adjacent accumulation (Zhao and Wu 2020), and priority accumulation of new information (Wu and Zhang 2018). Many scholars have done a lot of research on this basis. For example, Liu and Yu (2017) established a non-equigap and non-homogeneous grey model from the fractional perspective. Continuous grey model with conformable fractional derivative is proposed by Xie et al. (2020). Xi et al. (2019) studied the NGM(1,1) model from the perspective of new information priority.

In 1993, the non-equigap GM(1,1) model (NGM(1,1)) was proposed by Shi (1993). Later, Professor Deng, the founder of the grey system theory, also conducted a series of studies on the NGM(1,1) model (Deng 1994; Deng 1997). The stability of non-equigap grey control model is studied by Xiao and Li (2009). On this basis, a series of new grey models was derived from non-equigap angles. For example, Wang et al. (2012) proposed the NGM(1,1) power model, which has a very flexible form. Both the NGM(1,1) model and the grey Verhulst model are special forms of the NGM(1,1) power model. Wang and Li (2019) used the non-equigap grey Verhulst model to study the relationship between CO2 and economic growth. Meanwhile, Li et al. (2020) also studied the non-equigap grey Bernoulli model to analyze the relationship between economic growth and pollutants. The above two latest studies have studied the relationship between the two factors, which inspired the inspiration for this paper. The model in this paper is different from previous time series models, replacing time series with energy consumption. Under other conditions unchanged, analyze the impact of energy consumption on carbon dioxide emissions.

This paper proposes a non-equigap GM(1,1) model with conformable fractional accumulation. The conformable fractional order accumulation was first proposed by Ma et al. (2020). This new type of fractional order is easier to calculate and has a better fitting effect. Some scholars have studied the non-equigap grey model from the perspective of cumulative generation. At present, no scholar has introduced this new type of accumulation operator into the NGM(1,1) model, so a non-equigap GM(1,1) model with conformable fractional accumulation has been proposed. This model considers both energy consumption and carbon dioxide emissions. The nonlinear relationship between energy consumption and carbon dioxide emissions can be reflected by this model.

Non-equigap GM(1,1) model with conformable fractional accumulation

There is an original sequence X(0)(tk) = (x(0)(t1), x(0)(t2), ⋯, x(0)(tn)). If the gap is Δtk = tk − tk − 1 ≠ c, k = 2, 3, ⋯, n, c is a constant, and X(0)(tk) is called a non-equigap sequence. X(r)(tk) is ther order cumulative sequence of X(0)(tk):

X(r)(tk) = (x(r)(t1), x(r)(t2), ⋯, x(r)(tn))

where:

$$ {x}^{(r)}\left({t}_k\right)=\Big\{{\displaystyle \begin{array}{c}{x}^{(0)}\left({t}_1\right)\begin{array}{cc}& \end{array}\kern1.00em \begin{array}{cc}& \end{array}k=1\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}\\ {}\begin{array}{cc}\sum \limits_{j=2}^k\frac{x^{(0)}\left({t}_j\right)\times \varDelta {t}_j}{{t_j}^{1-r}}& k=2,3,\cdots, n\end{array}\end{array}} $$

the range of r is (0, 1].

The contiguous mean generation sequence of the non-equigap sequence X(r)(tk) is Z(r) = (z(r)(2), z(r)(3), ⋯, z(r)(n)):

where:

$$ {z}^{(r)}\left({t}_k\right)=\frac{1}{2}\left({x}^{(r)}\left({t}_k\right)+{x}^{(r)}\left({t}_{k-1}\right)\right) $$
(1)
$$ {x}^{(0)}\left({t}_k\right)+a{z}^{(r)}\left({t}_k\right)=b $$
(2)

is called the mean value form of CFNGM(1,1). Its whitening differential equation is:

$$ \frac{\mathrm{d}{x}^{(r)}}{\mathrm{d}t}+a{x}^{(r)}=b $$
(3)

a is called the development coefficient and b is called grey action. The least square estimation of CFNGM(1,1) model satisfies:

$$ \left[\begin{array}{c}a\\ {}b\end{array}\right]={\left({B}^TB\right)}^{-1}{B}^TY $$
(4)

where:

$$ B=\left[\begin{array}{cc}-{z}^{(r)}\left({t}_2\right)& 1\\ {}-{z}^{(r)}\left({t}_3\right)& 1\\ {}\vdots & \vdots \\ {}-{z}^{(r)}\left({t}_n\right)& 1\end{array}\right],\kern0.5em Y=\left[\begin{array}{c}\frac{x^{(r)}\left({t}_2\right)-{x}^{(r)}\left({t}_1\right)}{t_2-{t}_1}\\ {}\frac{x^{(r)}\left({t}_3\right)-{x}^{(r)}\left({t}_2\right)}{t_3-{t}_2}\\ {}\vdots \\ {}\frac{x^{(r)}\left({t}_n\right)-{x}^{(r)}\left({t}_{n-1}\right)}{t_n-{t}_{n-1}}\end{array}\right] $$
(5)

The initial condition of the differential equation in Eq. (3) is:

$$ {x}^{(r)}(1)={x}^{(0)}(1) $$

The time response equation can be obtained as:

$$ {\hat{x}}^{(r)}\left({t}_k\right)=\left({x}^{(0)}\left({t}_1\right)-\frac{\hat{b}}{\hat{a}}\right){e}^{-\hat{a}\left({t}_k-{t}_1\right)}+\frac{\hat{b}}{\hat{a}} $$
(6)

So, the reduction sequence of \( {\hat{x}}^{(r)}(k) \) can be obtained as:

$$ {\hat{x}}^{(0)}\left({t}_k\right)=\Big\{{\displaystyle \begin{array}{cc}{\hat{x}}^{(1)}\left({t}_1\right)& k=1\\ {}\frac{t_k^{1-r}\left({\hat{x}}^{(r)}\left({t}_k\right)-{\hat{x}}^{(r)}\left({t}_{k-1}\right)\right)}{\varDelta {t}_k}& k=2,3,\cdots, n\end{array}} $$
(7)

This paper uses MAPE to measure the stability of the model, as shown in Eq. (8).

$$ \mathrm{MAPE}=\frac{1}{n}\sum \limits_{k=1}^n\frac{\mid {x}^{(0)}(k)-{\hat{x}}^{(0)}(k)\mid }{x^{(0)}(k)} $$
(8)

Validity verification

Verify the validity of the CFNGM(1,1) model

In this paper, two cases are used to verify the validity of the CFNGM(1,1) model. The model testing process is shown in the “Verify the validity of the CFNGM(1,1) model” section of the paper. In the “Applicability test of CFNGM(1,1) model in developed and developing countries” section, the USA and India are selected as representatives of developed and developing countries to verify the applicability of the model. The test results all show that the model has good predictive performance.

Case 1

The data in reference Wang et al. (2012) were selected to compare the fitting effects of the two models. In this set of data, the order is r = 0.99, parameters a = 0.001028, and b = 539.0238, and the final fitting results are shown in Table 1. It can be seen from Table 1 that the fitting error of the CFNGM(1,1) model is 0.32%, and the fitting error of the non-equigap GM(1,1) power model is 0.94%. Therefore, we say that the CFNGM(1,1) model is better than the non-equigap GM(1,1) power model.

Table 1 Compare the fitting results of different models in case 1

Case 2

The data in reference Xi et al. (2019) were selected for comparison between fitting and prediction effects. In case 2, the conformable fractional order of the CFNGM(1,1) model is r = 0.955122, the parameters a =  − 0.00058, b = 9.199959, and the fitting and prediction results are shown in Table 2. In Table 2, the fitting error of the CFNGM(1,1) model is 0.8448%, and the prediction error is 1.7219%. The fitting effect and prediction effect of CFNGM(1,1) model are better than reference Xi et al. (2019). It can be concluded from these two cases that the model proposed in this paper is reasonable and effective.

Table 2 Compare the fitting results of different models in case 2

Applicability test of CFNGM(1,1) model in developed and developing countries

The USA and India are selected to verify the applicability of the model. The verification results are shown in Table 3. In this verification process, data from 2011 to 2016 are selected for fitting, and data from 2017 to 2019 are used as predictions. Compare the 2017–2019 forecast results with the actual values to calculate MAPE. The MAPE in Table 3 is far less than 10%, which proves that the model is suitable for forecasting in developed and developing countries.

Table 3 Validation of CFNGM(1,1) model in developed and developing countries

Take the USA as an example, compare the nonlinear model with the linear model. Figure 1 is a linear model diagram of US carbon dioxide emissions. The horizontal axis is energy consumption, and the vertical axis is carbon dioxide emissions. Figure 1 is a scatter plot made by Excel 2010, which is generated by linear trend fitting. The prediction error of the linear model is 5.51%, which is higher than the prediction error of the nonlinear model.

Fig. 1
figure 1

Linear model of carbon dioxide emissions

Analyze the relationship between energy consumption and carbon dioxide emissions

The burning of a large number of fossil fuels such as coal, oil, and natural gas has caused the CO2 in the atmosphere to rise. Carbon emissions are mainly CO2 emissions. According to the information reviewed, CO2 contributes as much as 60% to the greenhouse effect. Therefore, it is necessary to study the relationship between energy consumption and carbon dioxide emissions.

The data in the “Analyze the relationship between energy consumption and carbon dioxide emissions” section comes from https://www.bp.com/en/global/corporate/energy-economics/statistical-review-of-world-energy/downloads.html. Although there are many data on energy and carbon dioxide emissions, we have used a grey model with the characteristics of “small samples and poor information” to solve the relationship between energy consumption and carbon dioxide emissions. The reason for this is that Professor Liu Sifeng, the developer of grey system theory, once said why we still use small samples to solve practical problems in the era of big data. There are four reasons:

  1. (1)

    Understanding that big data starts with small data,

  2. (2)

    Data with a certain shape in big data is often small data,

  3. (3)

    In the long history, data at a certain period of time is usually small data,

  4. (4)

    The information density of big data is low, and it has obvious characteristics of poor information.

Professor Liu believes that big data is like a sea of sand, and small data is the gold hidden in the sea of sand. Therefore, the grey model was selected to solve the relationship between energy consumption and carbon dioxide emissions.

North America

The “North America” section studies the relationship between energy consumption and carbon dioxide emissions in Canada, Mexico, and the USA. Firstly, according to the energy consumption and carbon dioxide emissions data from 2014 to 2019, the parameters a and b in the CFNGM(1,1) model are calculated to obtain the fitted sequence of carbon dioxide emissions. Secondly, according to the 2019 energy consumption growth rate announced by the BP World Energy Statistics Yearbook, assuming that energy consumption will continue to grow at this rate in the next three years, the energy consumption from 2020 to 2022 is calculated. Finally, on the basis of substituting parameters a and b into Eq. (6) to obtain the time response equation, the energy consumption data from 2020 to 2022 can be input into the time response equation to obtain the x(r) series. By using Eq. (7), carbon dioxide emissions for 2020-2022 can be calculated. Take Canada as an example to elaborate on the process.

Energy consumption is used as the independent variable, and carbon dioxide emissions as the dependent variable. There is an energy consumption sequence tk = (14.03, 13.99, 13.94, 14.11, 14.35, 14.21), and an initial sequence of carbon dioxide emissions X(0)(tk) = (553.46, 546.23, 537.78, 549.11, 565.64, 556.19).

  1. Step 1.

    The Δt is calculated.

$$ \varDelta t=\left(-0.04,-0.05,0.17,0.24,-0.14\right) $$
  1. Step 2.

    Particle swarm optimization algorithm is used to find the optimal fractional order r = 0.000001, so the cumulative sequence can be obtained.

$$ {x}^{(0.000001)}=\left(553.46,551.90,549.97,556.59,566.05,560.57\right) $$
  1. Step 3.

    The contiguous mean generation sequence of the non-equigap sequence X(0.000001)(tk) can be obtained as:

$$ {Z}^{(0.000001)}=\left(552.68,550.93,553.28,561.32,563.31\right). $$
  1. Step 4.

    The expression of B and Y is:

$$ B=\left[\begin{array}{cc}-552.68& 1\\ {}-550.93& 1\\ {}-553.28& 1\\ {}-561.32& 1\\ {}-563.31& 1\end{array}\right],Y=\left[\begin{array}{c}39.04\\ {}38.58\\ {}38.92\\ {}39.42\\ {}39.14\end{array}\right]. $$

The parameters \( \hat{a} \) and \( \hat{b} \) are calculated according to the least squares method.

$$ \hat{a}=-0.04339,\hat{b}=14.88381 $$
  1. Step 5.

    The parameters \( \hat{a}=-0.04339,\hat{b}=14.88381 \) are brought into Eq. (6):

$$ {\hat{x}}^{(0.000001)}=\left(553.46,551.91,549.97,556.58,565.99,560.49\right). $$
  1. Step 6.

    The restored sequence of the fitted value is:

$$ {\hat{x}}^{(0)}=\left(553.46,543.68,540.68,548.71,563.03,558.74\right). $$
  1. Step 7.

    The growth rate of energy consumption in 2018–2019 is − 0.9%. According to this growth rate, the energy consumption in 2020–2022 is calculated as t7 = 14.09, t8 = 13.96, t9 = 13.83.

  2. Step 8.

    The 2020–2020 energy consumption values 14.09, 13.96, and 13.83 are introduced into Eq. (6), and the time response sequence for the next 3 years can be obtained as:

$$ {\hat{x}}^{(0.000001)}\left({t}_7\right)=555.80,{\hat{x}}^{(0.000001)}\left({t}_8\right)=550.74,{\hat{x}}^{(0.000001)}\left({t}_9\right)=545.71 $$
  1. Step 9.

    The reduction sequence of the predicted value of carbon dioxide emissions from 2020 to 2022 is:

$$ {\hat{x}}^{(0)}\left({t}_7\right)=550.91,{\hat{x}}^{(0)}\left({t}_8\right)=542.87,{\hat{x}}^{(0)}\left({t}_9\right)=534.79 $$

The detailed calculation process of the CFNGM(1,1) model is given above by taking Canadian energy consumption and carbon dioxide emissions as examples. Table 4 shows the carbon dioxide emissions of Mexico and the USA. The fitting errors of Mexico and the USA are 0.61% and 0.75%, respectively. The fitting accuracy shows that the model is suitable for carbon dioxide emission prediction.

Table 4 Forecasting carbon dioxide emissions in Mexico and the USA. Energy consumption: Exajoules

Judging from the forecast results of the three countries in North America, Canada and Mexico have shown a downward trend in carbon dioxide emissions as energy consumption decreases. With the reduction of energy consumption in the USA, its carbon dioxide emissions have shown an upward trend. The USA shut down a large number of coal power plants, which should be one of the reasons for the reduction in energy consumption in the USA. Although the US coal use has fallen sharply, it is not enough to offset the rise in emissions from other economic sectors. Those neglected industrial plants and factories have also become a greater source of climate pollution. One of the reasons for the increase in carbon dioxide emissions in the USA may also be related to the weather. A relatively cold winter has led to a substantial increase in the demand for oil and natural gas heating in areas such as New England. Therefore, the USA is likely to see an increase in carbon dioxide emissions in the face of reduced energy consumption, so it should pay attention to this problem.

South America

The carbon dioxide emissions of four South American countries, including Argentina, Brazil, Colombia, and Venezuela, are predicted in the “South America” section. The prediction results are shown in Table 5. The energy consumption of each country from 2020 to 2022 is calculated based on the growth rate of energy consumption from 2018 to 2019. It can be seen from Table 5 that the carbon dioxide emissions of Argentina and Venezuela will be decreasing in the next 3 years. The carbon dioxide emissions of Brazil and Colombia will increase in the next 3 years. As energy consumption increases (decreases), carbon dioxide emissions show a corresponding increase (decrease) trend, but there is a nonlinear relationship between energy consumption and carbon dioxide emissions.

Table 5 Forecasting carbon dioxide emissions in Argentina, Brazil, Colombia, and Venezuela. Exajoules, million tonnes

Europe

In the “Europe” section, the carbon dioxide emissions of 14 European countries are predicted, and the predicted results are shown in Table 8. Table 6 and Table 7 are the original data of energy consumption and carbon dioxide emissions, respectively. Table 6 lists the growth rate of energy consumption in 2018–2019. Assuming that the next 3 years will continue to develop at this growth rate, the energy consumption for the next 3 years is calculated, as shown in Table 6.

Table 6 The original data of energy consumption. Exajoules
Table 7 The original data of carbon dioxide emissions. Million tonnes

The energy consumption in Table 6 is used as the independent variable, and the carbon dioxide data in Table 7 is used as the dependent variable to predict the carbon dioxide emissions of 14 European countries. In Table 8, Germany’s carbon dioxide emissions fitting error is the largest at 2.07%, which is far below 10%. Therefore, it shows that the model is suitable for carbon dioxide emission prediction in European countries.

Table 8 Forecast results of 14 European countries. Million tonnes

The carbon dioxide emission prediction results of these 14 European countries show that carbon dioxide emissions of Austria, Belgium, and Turkey will show an increasing as energy consumption increase. The carbon dioxide emissions of Finland, Germany, Italy, Poland, Spain, Ukraine, and the UK in the next 3 years will show a decreasing trend as energy consumption decreases. Greece’s carbon dioxide emissions have fluctuated up and down. Its carbon emissions are relatively stable.

France, Norway, and Sweden are quite special. The forecast of France and Norway shows that as energy consumption decreases, carbon dioxide emissions have increased. In April 2020, the French government launched the Energy Transition Action Schedule 2019–2028 “Multi-Year Energy Plan.” It plans to shut down 14 nuclear reactors by 2035 and reduce the proportion of nuclear power in France’s total electricity generation to 50%. The installed capacity of renewable energy power generation is expected to be significantly higher than the current level, with the new installations mainly coming from wind, electricity and solar energy. And the French government has also introduced a series of incentive policies to promote this plan. However, affected by the COVID-19, many of the previously announced renewable energy projects in France have not yet been released and tendered, and existing projects have also been delayed to varying degrees. Therefore, there may be an increase in carbon dioxide emissions in France. In order to reduce carbon emissions, starting from January 2020, Norwegian aviation fuel suppliers are required to add 0.5% biofuel to all aviation fuel. However, there are not many suppliers of this fuel. The demand for biofuels is strong, but the output is not enough. This way of reducing emissions has certain limitations. In Sweden, carbon dioxide decreases as energy consumption increases. This phenomenon in Sweden is in line with the facts. The energy consumption in this article refers to the consumption of primary energy. In addition to non-renewable energy such as coal and oil, primary energy also includes renewable energy such as solar energy, water power, and wind power. As early as 2012, hydropower in Sweden accounted for 48% of the total power generation, and wind power accounted for 4.4% of the total power generation. Hydropower and wind power have shown an upward trend, while fossil energy power generation has fallen by 8%. Therefore, although energy consumption is increasing in Sweden, its carbon dioxide may appear to decrease.

CIS

In the “CIS” section, the carbon dioxide emissions of six CIS countries are predicted. The original data of energy consumption and carbon dioxide emissions of CIS countries are shown in Table 9, and the growth rate data is also listed in Table 9.

Table 9 The original data of energy consumption and carbon dioxide emissions in CIS. Exajoules, million tonnes

The results of carbon dioxide emission forecasts for CIS countries are shown in Table 10. It can be seen in Table 10 that the fitting error of the six countries does not exceed 1%, indicating that the model is suitable for predicting the carbon dioxide emissions of the CIS. The forecast results show that the carbon dioxide emissions of Azerbaijan, Belarus, and Turkmenistan from 2020 to 2022 are on the rise with the increase in energy consumption. Kazakhstan, Russian Federation, and Uzbekistan show a downward trend in carbon dioxide emissions from 2020 to 2022 as energy consumption decreases.

Table 10 Forecast results of carbon dioxide emissions of the CIS. Million tonnes

Middle East

In 2019, the oil reserves in the Middle East accounted for 61.5% of the world’s oil reserves, and the oil production in the Middle East accounted for 30.7% of the world’s oil production. The Middle East is the region with the largest oil reserves and the most oil production and export in the world. In the “Middle East” section, the carbon dioxide emissions of six countries in the Middle East are predicted. The original data of energy consumption, carbon dioxide emissions, and the final prediction results are included in Table 11. From Table 11, it can be seen that the energy consumption of these six Middle Eastern countries has increased from 2018 to 2019. According to this growth rate, the energy consumption of 2020–2022 will also increase. The forecast results show that as energy consumption increases, carbon dioxide emissions in these six countries are also increasing year by year. Among them, Iran’s carbon dioxide emissions rank first among the six countries, followed by Saudi Arabia. Israel has the least carbon dioxide emissions.

Table 11 Forecast results of carbon dioxide emissions in the Middle East. Exajoules, million tonnes

Africa

Africa’s carbon dioxide emissions forecast results are shown in Table 12. The carbon dioxide emissions forecasts of South Africa, Eastern Africa, and Western Africa increase with the increase in energy consumption, and only Egypt’s carbon dioxide emissions show a downward trend. Although Egypt’s forecast shows a downward trend, it is still higher than the forecasts of Eastern Africa and Western Africa. South Africa’s carbon dioxide emission forecast results rank first among these four regions. The forecast value of South Africa in 2022 is 503.17, which is 3.77 times the forecast value of Eastern Africa in 2022. Therefore, South Africa should take measures to prevent carbon dioxide emissions from continuing to rise.

Table 12 Forecast results of carbon dioxide emissions in the Africa. Exajoules, million tonnes

Asia Pacific

The carbon dioxide emissions of 16 countries and regions in the Asia Pacific region are predicted in the “Asia Pacific” section. The original data of energy consumption in the Asia Pacific region is shown in Table 13. It can be seen from Table 13 that China’s energy consumption in the Asia Pacific region is the largest and is still on the rise. In China, energy consumption is inseparable from its population. As we all know, China is a country with a large population, and the electricity needed by these people needs to consume a lot of energy. However, in recent years, non-fossil energy power supply capacity in China has continued to increase.

Table 13 The original data of energy consumption in the Asia Pacific. Exajoules

The original data of carbon dioxide emissions from Asia Pacific regions from 2014 to 2019 are shown in Table 14. Energy consumption in the Asia Pacific regions ranks first is China, and it also has the highest carbon dioxide emissions. India and Japan’s carbon dioxide emissions are closely behind China. Carbon dioxide emissions of China are much higher than those of other countries in the Asia Pacific region. Carbon dioxide is the main greenhouse gas, so the forecast of carbon dioxide emissions can be warned in advance to remind the government whether to adjust the energy structure.

Table 14 The original data of carbon dioxide emissions in the Asia Pacific. Million tonnes

The forecast results of carbon dioxide emissions in the Asia Pacific region are shown in Table 15. The forecast results show that carbon dioxide emissions of 10 countries including Australia, Bangladesh, China, India, Indonesia, Pakistan, Philippines, Sri Lanka, Thailand, and Vietnam will show an upward trend from 2020 to 2022. The results of carbon dioxide emissions forecasts in China show an increase year by year, and will exceed 10 billion tonnes of carbon dioxide, which should be taken seriously. With the reduction of energy consumption in China Hong Kong SAR, Japan, Singapore, South Korea, and China Taiwan, carbon dioxide emissions have shown a downward trend. The forecast results show that as energy consumption increases, there is a trend of decreasing carbon dioxide emissions in Malaysia. Currently, the Ministry of Transport of Malaysia is promoting the development of green logistics to reduce the industry’s carbon footprint. The transportation industry is Malaysia’s second largest source of carbon emissions. Therefore, the goal of reducing carbon emissions is achieved through green logistics cooperation to improve efficiency.

Table 15 The forecast result of carbon dioxide emissions in the Asia Pacific. Million tonnes

Conclusion

The non-equigap GM(1,1) model with conformable fractional accumulation is proposed to predict the carbon dioxide emissions of 53 countries and regions around the world. In order to verify the effectiveness of the model in this paper, MAPE is used as a measurement standard, and two cases are compared to illustrate the validity of the model in this paper. The distribution intervals of all fractional orders are counted as shown in Fig. 2. In the context of the application of predicting carbon dioxide emissions, the proportion of conformable fractional orders less than 0.1 is 66.04%, and the conformable fractional orders in the range of 0.1–0.5 account for 15.09%, a total of 81.13%. This also shows that the conformable fractional order has a better prediction effect when the r is smaller.

Fig. 2
figure 2

Interval statistics of the conformable fractional distribution

In this article, energy consumption is used as input, and carbon dioxide emissions are used as output. The relationship between energy consumption and carbon dioxide emissions is analyzed. Energy consumption and carbon dioxide emissions show a nonlinear relationship. The carbon dioxide emission forecast is divided into seven parts, including North America, South America, Europe, CIS, Middle East, Africa, and Asia Pacific. In North America, carbon dioxide in the USA is on the rise. The forecast results of carbon dioxide emissions of South American countries Brazil and Colombia from 2020 to 2022 show an upward trend. With the growth of energy consumption, carbon dioxide emissions have also increased in three European countries, including Austria, Belgium, and Turkey. The forecast of France and Norway shows that as energy consumption decreases, carbon dioxide emissions have increased. In CIS, the countries with increasing carbon dioxide emissions are Azerbaijan, Belarus, and Turkmenistan. In the Middle East, countries with increased carbon dioxide emissions include Iran, Israel, Kuwait, Qatar, Saudi Arabia, and United Arab Emirates. In Africa, the forecast results of carbon dioxide emissions in the three regions have increased, including South Africa, Eastern Africa, and Western Africa. In Asia Pacific, countries with rising carbon dioxide emissions include Australia, Bangladesh, China, India, Indonesia, Pakistan, Philippines, Sri Lanka, Thailand, and Vietnam. The above are the countries with increased carbon dioxide emissions forecast results. There are also special circumstances in the prediction results. For example, the prediction results of the USA, France, and Norway show that as energy consumption decreases, carbon dioxide emissions have increased. In Sweden and Malaysia, energy consumption is increasing, but carbon dioxide is decreasing. The emergence of these situations is closely related to the country’s carbon dioxide emission policy. Enough attention should be paid to countries that predict increased carbon dioxide emissions.

There may be some uncertain factors in the forecasting process, such as technological changes and government policies. Assuming that certain uncertain factors appear in the historical trend of carbon dioxide that will affect future changes in carbon dioxide emissions, buffer operators can be used to eliminate this effect. And the energy consumption data used in the prediction of future carbon dioxide is calculated based on the latest data. To prevent the interference of uncertain factors, we can only forecast the data for the next year. After the actual data appears in the next year, the model data will be updated again to ensure the accuracy of the forecast. These suggestions are provided for reference when similar problems arise.