key: cord-0754355-4n1a1c8z
authors: Wang, Huiping; Wang, Yi
title: Estimating CO(2) emissions using a fractional grey Bernoulli model with time power term
date: 2022-02-17
journal: Environ Sci Pollut Res Int
DOI: 10.1007/s11356-022-18803-0
sha: ac91d64cb6f0e402fc029dead8123236e66f2085
doc_id: 754355
cord_uid: 4n1a1c8z

Global warming caused by CO(2) emissions will directly harm the health and quality of life of people. Accurate prediction of CO(2) emissions is highly important for policy-makers to formulate scientific and reasonable low-carbon environmental protection policies. To accurately predict the CO(2) emissions of the world’s major economies, this paper proposes a new fractional grey Bernoulli model (FGBM(1,1,[Formula: see text] )). First, this paper introduces the modeling mechanism and characteristics of the FGBM(1,1,[Formula: see text] ) model. The new model can be transformed into other grey prediction models through parameter adjustment, so the new model exhibits high adaptability. Second, this paper employs four carbon emission datasets to establish a grey prediction model, calculates model parameters with three optimization algorithms, adopts two evaluation criteria to evaluate the accuracy of the model results, and selects the optimization algorithm and model results that yield the highest model accuracy, which verifies that the FGBM(1,1,[Formula: see text] ) model is more feasible and effective than the other six grey models. Finally, this paper applies the FGBM(1,1,[Formula: see text] ) model to predict the CO(2) emissions of the USA, India, Asia Pacific, and the world over the next 5 years. The forecast results reveal that from 2020 to 2024, the CO(2) emissions of India, the Asia Pacific region, and the world will gradually rise, but that in USA will slowly decline over the next 5 years.

Global warming caused by the greenhouse effect is one of the factors that seriously threatens human survival and development, and the increase in CO 2 emissions is considered the main cause of the greenhouse effect. According to the Global Climate 2015-2019 report released by the World Meteorological Organization, the growth rate of CO 2 in the atmosphere from 2015 to 2019 was 18% higher than that during the previous 5 years, and the average temperature was 0.2 °C higher than that during the previous 5 years. Hence, this period encompasses the hottest 5 years on record. British Petroleum (BP) World Energy Statistics (2020) indicated that in 2019, the global CO 2 emissions reached 34,169 million tons. Among the various countries, the USA and India rank second and third worldwide, respectively, in terms of their CO 2 emissions, and the total CO 2 emissions of these two countries account for 21.8% of the world emissions. The Asia Pacific region accounts for 50.5% of the global CO 2 emissions. Upon entering the twenty-first century, the growth of global CO 2 emissions is soaring. Thus, from 2009 to 2019, the CO 2 emissions in the USA basically remained at approximately 5000 MT, revealing a trend of continuous fluctuation. The CO 2 emissions in India grew the fastest, by 55.4%. The growth rate of the CO 2 emissions in the Asia Pacific region reached 30.4%, which is also one of the fastest growing economies in the world. From the perspective of global CO 2 emissions, although the growth rate is not high, with an average annual growth rate of 1.1%, CO 2 emissions remain on the rise. To achieve the target of the Paris Agreement, more than 20 countries put forward carbon neutrality goals at the 2020 Climate Ambition Summit, and more than 40 countries established new commitments to independently improve national contributions. Therefore, accurate prediction of future CO 2 emission data of the USA, India, Asia Pacific region, and the world can help policy-makers formulate more scientific and reasonable environmental policies to truly achieve the goal of carbon emission reduction.

The existing models for CO 2 emission forecast can be classified into three categories. The first category includes nonlinear intelligent models, such as the least squares support vector machine (Sun and Liu 2016) , extreme learning machine (Sun and Sun 2017) , generalized regression neural network (Heydari et al. 2019) ; and improved chicken swarm optimization (ICSO-SVM) model using the ICSO algorithm to optimize support vector machine parameters (Wen and Cao 2020) . The second category involves the statistical analysis model, which studies the quantitative relationship between CO 2 emissions and influencing factors and applies the relationship equation to predict CO 2 emissions. These approaches include trend analysis (Köne and Büke 2010) , logistic equations (Meng and Niu 2011) , improved Gaussian process regression (Fang et al. 2018) , panel quantile regression (Zhu et al. 2018) , log-average decomposition index (Xu et al. 2019) , and comprehensive methods combining multiple regression analysis, input-output techniques, and structural decomposition analysis (Xia et al. 2019) . The last category is the grey prediction model, which was first proposed by Deng (1982) . The grey prediction model was established founded on a small amount of incomplete information to describe the development trend of objects more accurately. Compared to machine learning-based and statistical prediction methods founded on large data samples, the grey prediction model can realize the simulation and prediction of small data samples. The classical grey model GM(1,1) and its extended grey prediction model are widely used in the fields of energy, environment, and social management (Tsai 2016; Liu et al. 2020; Liu et al. 2021a; Wang et al. 2020b; Liu et al. 2021b ).

There are two main kinds of grey single-variable forecasting models: the first-order grey differential model GM(1,1) and the grey Bernoulli model GBM(1,1). Deng first proposed the GBM(1,1) model in 1985; namely, a power exponent was introduced into the differential Bernoulli equation. When the exponent is equal to 2, the model is also referred to as the grey Verhulst model. Chen et al. (2008) proposed the nonlinear grey Bernoulli model NGBM(1,1) for the first time. Compared to the general GM(1,1) model, the NGBM(1,1) model can better reflect the nonlinear growth trend of data series. Chen et al. (2010) proposed a new grey prediction model NGBM(1,1), by optimizing the background value and power index simultaneously. Subsequently, many scholars improved the NGBM(1,1) model from different perspectives; e.g., Pao et al. (2012) proposed an iterative method to optimize the parameters of the NGBM(1,1) model. Wang (2017) optimized the background coefficient and initial conditions and considered the weighted method in analysis. Guo et al. (2016) proposed a new model by combining the self-memory principle of a dynamic system with the NGBM(1,1) model. Ma et al. (2019) constructed the NGBM(1, n) model by combining the GMC(1, n) model and Bernoulli equation. Wu et al. (2019b) extended the first-order accumulation operation and established the FANGBM(1,1) model based on fractional accumulation. Liu and Xie (2019) proposed the establishment of the WBGM(1,1) model by combining the fitting performance of the NGBM(1,1) model with the Weibull cumulative distribution. Şahin (2020) proposed the OFANGBM(1,1) model based on the integral mean value theorem. Wu et al. (2020a) established the NGBM(1,1, k, c) model by combining the NGM(1,1, k, c) model with the NGBM(1,1) model. Jiang and Wu (2021) constructed a nonlinear grey Bernoulli model based on the fractional order reverse accumulation, namely, the FANGBM(1,1) model. At present, a simpler and more convenient conformable fractional cumulative grey model (CFGM) was proposed (Ma et al., 2020) . On this basis, Xie et al. (2020) proposed the conformable fractional grey model in opposite direction (CFGOM). Zheng et al. (2021) proposed the conformable fractional nonhomogeneous Bernoulli model (CFNHGBM (1,1, K)).

The grey prediction model is widely implemented in CO 2 emission prediction. Lin et al. (2011) applied the grey model to predict CO 2 emissions in Taiwan. Pao et al. (2012) used the NGBM model to predict the CO 2 emissions and real GDP growth of China. Lotfall et al. (2013) forecast CO 2 emissions based on a grey model and an autoregressive integrated moving average (ARIMA) model. The prediction accuracy of these two methods was compared according to the root mean square error (RMSE), mean absolute error (MAE), and mean absolute percentage error (MAPE). Gao et al. (2015) established a new discrete fractional order cumulation model, i.e., FAGM(1,1, D), to predict CO 2 emissions. Hamzacebi and Karakurt (2015) used the grey prediction model to predict the energy-related CO 2 emissions in Turkey. Yuan et al. (2017) established a linear programming model reflecting the relationship between the economic development and CO 2 emissions in China and employed the GM(1,1) model to predict the parameters of the planning model. Wang and Ye (2017) developed a nonlinear multivariable grey model to discuss the relationship between economic growth and CO 2 emissions. Xu et al. (2019) combined an adaptive grey model with the buffer rolling method to predict the greenhouse gas emissions in China from 2017 to 2025. Wang and Li (2019) adopted the nonequidistant grey Verhulst model to analyze the relationship between CO 2 emissions and economic growth. Wu et al. (2020b) implemented a conformable fractional nonhomogeneous grey model to predict the CO 2 emissions of BRICS countries. Chiu et al. (2020) proposed a multivariate grey prediction model using neural networks based on feature selection and residual correction to predict China's carbon emissions. Based on the grey Verhulst model, Duan and Luo (2020) introduced an extrapolation method to optimize the background value and predicted the CO 2 emissions for three coal resources in China. Wang et al. (2020a) established the metabolic nonlinear grey model (MNGM)-ARIMA method, established a new MNGM-BPNN combination model based on the MNGM model and back-propagation (BP) neural network and analyzed the CO 2 emissions in China, the USA, and India with these two methods. Zhou et al. (2021) proposed a method to process original sequence data by means of an average weakening buffer operator based on the grey rolling mechanism of the new information priority principle and predicted the trend of CO 2 emissions in China. Xie et al. (2021) established a new continuous conformable fractional nonlinear grey Bernoulli model to forecast CO 2 emissions from fuel combustion in China.

Scholars have greatly promoted the optimization and application of grey models. However, although the existing studies have optimized the grey model considering the structure or parameters, each optimization method only improves the model performance to a certain extent, and the accuracy remains insufficient. In addition, most models only apply one optimization algorithm to determine the optimal parameters, and research on the application of multiple optimization algorithms is rare. Therefore, to better predict carbon emissions, based on the optimization of existing models, namely, NGBM(1,1) and FAGM(1,1,t ) (Wu et al. 2019a) , this paper proposes a new grey prediction model FGBM(1,1,t ) and considers a variety of optimization algorithms to determine the optimal structural parameters of the model. The main contributions of this paper are as follows:

(1). Based on the advantages of the grey Bernoulli prediction model NGBM(1,1) and FAGM(1,1, t ), a new model, i.e., FGBM(1,1,t ), is proposed. The new model can be transformed into other grey prediction models by changing its parameters. (2). In the FAGM(1,1,t ) model, application of the trapezoidal integral method to obtain an approximate solution yields errors, but this paper provides an analytical solution of the time response function in the FGBM(1,1,t ) model, and the result is more accurate. (3). Grey wolf optimization (GWO), particle swarm optimization (PSO), and quantum genetic algorithm (QGA) are applied to solve the parameters of the FGBM(1,1,t ) model. (4). The FGBM(1,1,t ) model is used to predict the CO 2 emissions of four countries and regions over the next 5 years.

This paper is organized as follows: the second section first introduces the basic theory of the NGBM(1,1) and FAGM(1,1,t ) models. Second, a new grey prediction model FGBM(1,1,t ) is constructed, which mainly includes the modeling mechanism, model characteristics, and solution method. Then, three common optimization algorithms (GWO, PSO, and QGA) are employed to solve the parameters of the FGBM(1,1,t ) model, and two widely adopted error measures are introduced. The third section considers the CO 2 emission data of four economies to verify the feasibility and effectiveness of the FGBM(1,1,t ) model. The fourth section applies the FGBM(1,1,t ) model to predict the CO 2 emissions of the USA, India, Asia Pacific region and the world over the next five years. The fifth section contains the conclusion of this paper.

This section mainly introduces the grey Bernoulli model NGBM(1,1) and the FAGM(1,1,t ) model as follows:

Definition 1: Given a nonnegative sequence

} is referred to as the firsto r d e r ge n e r a t i n g s e qu e n c e o f X (0) , w h e r e Chen et al. (2008) , it can be found that the expression of NGBM(1,1) is:

The above is a nonlinear equation, and the exponent exponent ≠ 1.

The parameters a, b of the NGBM(1,1) model are calculated as follows:

and n is the sample number of the modeling sequence.

The time response function of the NGBM(1,1) model is as follows:

The predicted values of the model are as follows:

Definition 2: Given a nonnegative sequence

, and its r -th order accumulation sequence is (Wu et al, (2013) . Denoted as matrix A r , the r-AGO matrix that satisfies X (r) = A r X (0) is:

Denoted as matrix D r , the r-IAGO matrix that satisfies

Matrices A r and D r satisfy A r D r = I n .

The whitening differential equation of the FAGM(1,1,t ) model is:

Applying integration to Eq. (7) in the time interval of [k − 1, k], we obtain: B a s e d o n t h e t r a p e z o i d e q u a t i o n a n d z (r) 

Equation (8) is written as:

Based on Eq. (9), unknown parameters a, b, c of FAGM(1,1,t ) satisfy:

The time response function of the FAGM(1,1,t ) model is:

Moreover, the restored value of x (0) (k), k = 2, 3, … , n is given as:

Description of the FGBM(1,1,t˛)

The expression of the FGBM(1,1,t ) model is as follows:

where 0 ≤ r ≤ 1, 0 ≤ ≤ 4, 0 ≤ ≤ 3, ≠ 1. Theorem 1. The time response function of FGBM (1,1,t ) is derived as:

The general solution formula of the equation is as follows:

To perform the definite integral on interval [1, t] , we know that:

Although is not an integer, ∫ k 1 e a � ( −1) d can be integral by numerical integration method to obtain the real number.

The predicted values x (0) (k) can be obtained by:

By integrating on [k − 1, k] both side of Eq. (16) at the same time, the following conclusion can be obtained:

According to the integral median theorem, we can get:

The results are as follows:

According to the commonly used method to solve the parameters of grey prediction model, the least square criterion of FGBM (1,1,t ) model can be described as the following unconstrained optimization problems:

The solution of this optimization problem is:

where In order to facilitate readers to understand the solution process of the FGBM (1,1,t ) model, Fig. 1 shows the algorithm implementation process of FGBM (1,1,t ) model.

The expression of the FGBM(1,1,t ) model indicates that the FGBM(1,1,t ) model can be transformed into other grey prediction models when different values are assigned to parameters and r, , such as GM(1,1) (Deng 1982 (Wu et al. 2020a ). The relationships between FGBM (1,1,t ) with other grey models are shown in the Fig. 2 . To minimize the model error, we must determine the optimal values of parameters r, , , and adopt the MAPE as the main criterion. The optimization problems when solving the optimal parameters are as follows: The above optimization problems are essentially nonlinear programming problems with equality constraints, which can be solved with intelligent optimization algorithms or heuristic algorithms. In this paper, GWO, PSO, and QGA are adopted to solve the parameters. The reason for choosing these three algorithms is: GWO algorithm has the characteristics of simple structure, has few parameters to be adjusted, and is easy to implement. Among them, there are adaptive convergence factors and information feedback mechanism, which can achieve a balance between local optimization and global search. Therefore, GWO algorithm has good performance in solving the problem accuracy and convergence speed. The advantage of PSO algorithm is that it has memory ability. In the implementation of intelligent search, it can combine the individual and global optimal location to realize location. It has a very fast speed of approaching the optimal solution and can effectively optimize the parameters of the system. The population coding method used by QGA greatly

, t = 1, 2, 3, ⋯ , n enriches the diversity of the population, at the same time, the quantum revolving gate is used to update the population and evolve based on the information of the current optimal individual. It has the characteristics of strong adaptability, fast convergence and suitable for global search. The parameters of the algorithm are shown in Table 1 .

This section considers two error metrics widely used for prediction models to test the effectiveness and applicability of the grey prediction model, as indicated in Table 2 .

In this section, four cases are presented to check the accuracy of the FGBM(1,1,t ), and the forecasting results are compared to those obtained with other grey prediction models, such as NGM(1,1), SIGM, GMP, FGM, NGBM(1,1,k,c), and FAGM(1,1,t ). In Sections 3.1 to 3.4, the raw data in Table 1 and the above seven grey models are used to simulate and 

According to Statistics Review of World Energy 2020, the USA produces the highest CO 2 emissions worldwide. In 2019, the CO 2 emissions of the USA reached 4964.69 million tons, accounting for 14.5% of the total global CO 2 emissions, ranking second in the world. The task of emission reduction in the USA is far from realization. In this regard, the USA has put forward the Zero Carbon Action Plan (ZCAP), which deeply considers stabilization of domestic employment, economic innovation, and coordinated development of the environment in the postepidemic era. Therefore, this section adopts the USA as an example to test the accuracy of the FGBM(1,1,t ) model. The parameters and MAPE values of FGBM(1,1,t ) calculated with the three optimization algorithms are listed in Table 3 . Figure 3 shows the number of iterations of the three optimization algorithms and the relationship between MAPE and the parameters. The MAPE and test_MAPE values for QGA are 1.1320% and 1.5805%, respectively, which are the smallest among those for the three optimization algorithms, and the performance is the best. Therefore, the QGA results are selected here. The fitting and prediction results for each model are shown in Fig. 4 and Table 4 , and the error measurement results for each model are shown in Fig. 5 and Table 5 . The prediction value of the FGBM(1,1,t ) model is the closest to the actual value. Based on the fitting values, the two error measurement indexes of the FGBM(1,1,t ) model are the best. This also indicates that the FGBM(1,1,t ) model is more accurate than the other grey models in predicting the CO 2 emissions in the USA.

According to Statistics Review of World Energy 2020, in 2019, the CO 2 emissions in India reached 2480.35 million tons, accounting for 7.3% of the global CO 2 emissions, ranking third in the world. From 2010 to 2019, the CO 2 emissions in India increased, at an annual growth rate of 5.48%. With the rapid growth of the Indian economy and population, its CO 2 emissions also increased notably, which is estimated to account for 11% of the global CO 2 emissions by 2030. In the face of the call of the international community to promote carbon emission reduction, the Indian government has promised to reduce greenhouse gas emissions by 33-35% from 2015 to 2030. Therefore, this section chooses India as an example to test the accuracy of the FGBM(1,1,t ) model. The parameters and MAPE values of the FGBM(1,1,t ) model calculated with the three optimization algorithms are listed in Table 6 . The MAPE and test_MAPE values for GWO are 0.5090% and 0.7130%, respectively, which are the smallest among the three optimization algorithms, and the performance is the best. Therefore, GWO is selected here. Figure 6 shows the number of iterations and the relationship between MAPE and the parameters. The fitting and prediction results for each model are shown in Fig. 7 and Table 7 , and the error measurement results for each model are shown in Fig. 8 and Table 8 . The prediction value of the FGBM(1,1,t ) model is the closest to the actual value. Based on the prediction and fitting values, the two error measurement indexes of the FGBM(1,1,t ) model are the best. This further demonstrates that the FGBM(1,1,t ) model is more accurate than the other grey models in predicting the CO 2 emissions of India. 

Mean absolute error respectively, which are the smallest and the best among the three optimization algorithms, so the PSO algorithm is selected here. Figure 9 shows the relationship between the number of iterations, MAPE, and parameters. The fitting and prediction results for each model are shown in Fig. 10 and Table 10 , and the error measurement results for each model are shown in Fig. 11 and Table 11 . The predicted value of the FGBM(1,1,t ) model is the closest to the actual value. According to the prediction and fitting values, the two error metrics of the FGBM(1,1,t ) model are the best. This also indicates that the FGBM(1,1,t ) model is more accurate than the other grey prediction models in predicting the CO 2 emissions in the Asia Pacific region.

According to Statistics Review of World Energy 2020, in 2019, the global CO 2 emissions reached 34,169.00 million tons, which has increased greatly over the past few decades. From 2010 to 2019, the global CO 2 emissions increased by 1.102% annually, exhibiting a slow rising trend. However, economic recovery after the COVID-19 pandemic will facilitate a global rebound in CO 2 emissions. Therefore, this section adopts the world as an example to test the accuracy of the FGBM(1,1,t ) model. The parameters and MAPE values of the FGBM(1,1,t ) model calculated with the three optimization algorithms are summarized in Table 12 . The MAPE and test_MAPE values obtained with QGA are 0.168% and 2.352%, respectively, which are the smallest and the best among the three optimization algorithms. Therefore, QGA is selected here. Figure 12 shows the relationship between the number of iterations, MAPE, and parameters. The fitting and prediction results for each model are shown in Fig. 13 and Table 13 , and the error measurement results for each model are shown in Fig. 14 and Table 14 . The predicted value of the FGBM(1,1,t ) model is the closest to the actual value. Based on the prediction and fitting values, the two error metrics of the FGBM(1,1,t ) model are the best. This further verifies that the FGBM(1,1,t ) model is more accurate than the other grey prediction models in predicting the global CO 2 emissions.

The above results reveal that although FGBM(1,1,t ) achieves the best effect in predicting the global CO 2 emissions among the many models, due to the sudden increase in the original data of CO 2 emissions in 2016, the fitting effect of the model is good, but the prediction error is large. The predicted global CO 2 emissions over the next 5 years may be far from the real value. To resolve this problem, global carbon emission data can be divided into two parts, namely, Organisation for Economic Co-operation and Development 16 2011 1735.15 1749.11 1749.12 1749.45 1751.62 1737.49 1751.26 1745.24 2012 1848.13 1849.31 1849.32 1849.41 1848.95 1846.63 1847.28 1845.79 2013 1929.35 1949.63 1949.64 1949.65 1948.07 1956.97 1944.88 1949.24 2014 2083.54 2050.06 2050.07 2050.17 2047.3 2060.50 2045.31 2051.17 2015 2149.38 2150.62 2150.63 2150.95 2145.88 (OECD) and non-OECD data. The trend of CO 2 emissions based on OECD and non-OECD data is more stable than that based on global data. Therefore, this paper applies the FGBM(1,1,t ) model to fit and predict the CO 2 emissions based on OECD and non-OECD data, as listed in Tables 15,   16 , and 17. Tables 15 and 16 reveal that between the OECD and non-OECD data, the prediction error obtained with QGA is the smallest, so QGA is selected to calculate the model parameters. The global CO 2 emissions are determined by adding the OECD and non-OECD-based CO 2 emissions, as summarized in Table 17 . The fitting and prediction errors of the global CO 2 emissions are listed in Table 18 and are lower than those listed in Table 13 . Table 19 provides the OECD and non-OECD based CO 2 emissions over the next 5 years. The global CO 2 emissions over the next 5 years can be obtained by adding the two types of emissions. The fitting and prediction results are more accurate than the direct prediction results of the global CO 2 emissions.

In this section, we employ the FGBM(1,1,t ) model to predict the CO 2 emissions of the USA, India, Asia Pacific region, and the world over the next 5 years (2020-2024). The prediction results are summarized in Table 20 and Fig. 15 . The CO 2 emissions in India, the Asia Pacific region, and the world will gradually increase over the 5 years. In addition, the CO 2 emissions in the USA will slowly decline.

To better describe the future CO 2 emissions of the USA, India, Asia Pacific region, and the world, a new grey prediction model, i.e., the FGBM(1,1,t ) model, is proposed based on the NGBM(1,1) and FAGM(1,1,t ) models, and a precise solution of the new model is obtained via the numerical integration method. Moreover, this paper applies three common optimization algorithms to calculate the model parameters. By changing the model parameters, the FGBM(1,1,t ) model can be transformed into other models, so the model achieves a strong adaptability. The CO 2 emission fitting and forecasting results in the above four economies indicate that the FGBM(1,1,t ) model is more effective and accurate than the existing NGM(1,1), SIGM, GMP, FGM, NGBM(1,1,k,c), and FAGM(1,1,t ). Moreover, we employ the FGBM(1,1,t ) model to predict the CO 2 emissions of the USA, India, Asia Pacific region, and the world over the next 5 years. The forecast results reveal that from 2020 to 2024, the CO 2 emissions of India, the Asia Pacific region, and the world will gradually rise. The CO 2 emissions of the USA will slowly decline over the next 5 years. Notably, the grey prediction FGBM(1,1,t ) model can be applied not only in the prediction of CO 2 emissions but also in the prediction of other data, with a high adaptability. In order to meet the requirements for carbon emission reduction in the sustainable development goals, several suggestions are put forward based on the above findings. First, promoting the green and low-carbon transformation of the energy system is the key. We should focus on energy conservation and improve the dual control system of total energy consumption and intensity. We should increase the utilization ratio of renewable energy such as wind power and photovoltaic power generation, and develop hydropower, geothermal energy, marine energy, hydrogen energy, biomass energy, and photothermal power generation according to local conditions. Second, the fiscal and tax support should be strengthened to encourage the research and development of green and low-carbon technologies. Financial funds should be used to support the development of green environmental protection industry and energy efficient utilization, and a number of scientific and technological projects should be arranged in the fields of energy conservation and environmental protection, cleaner production, and clean energy. Third, the recycling of renewable resources should be strengthened to realize a green and low-carbon lifestyle. We should promote the classification, reduction, and resource utilization of domestic waste according to local conditions, accelerate the construction of waste material recovery system, and strengthen the recycling of waste paper, waste plastics, waste tires and other resources.

This study still has some limitations. First, the model has many hyper-parameters, which may lead to over-fitting problems. Second, the model is a univariate model, which may ignore some important influencing factors (e.g., economic development, industrial policy, etc.), so the multivariate grey model can also be applied to predict carbon emissions, which is the optimization direction in the future.

Author contribution All authors contributed to the study conception and design. Material preparation, data collection, and analysis were performed by Huiping Wang and Yi Wang. The first draft of the manuscript was written by Yi Wang. All authors read and approved the final manuscript. Data availability The datasets used or analyzed during the current study are available from the corresponding author on reasonable request. 

Ethics approval Not applicable.

Consent to publish Not applicable.

The authors declare no competing interests.

Forecasting of foreign exchange rates of Taiwan's major trading partners by novel nonlinear grey Bernoulli model NGBM(1, 1)

Forecasting Taiwan's major stock indices by the Nash nonlinear grey Bernoulli model

A multivariate grey prediction model using neural networks with application to carbon dioxide emissions forecasting

Novel gray forecasting model and its modeling mechanism

Control problems of grey systems

Grey optimization Verhulst model and its application in forecasting coal-related CO 2 emissions

A novel method for carbon dioxide emission forecasting based on improved Gaussian processes regression

Estimation of Chinese CO 2 emission based on a discrete fractional accumulation grey model

Forecasting China's SO 2 emissions by the nonlinear grey Bernoulli self-memory model

Forecasting the energy-related CO 2 emissions of Turkey using a grey prediction model

Renewable energies generation and carbon dioxide emission forecasting in microgrids and national grids using GRNN-GWO methodology

Forecasting of CO 2 emissions from fuel combustion using trend analysis

Nonlinear grey Bernoulli Model with fractional-order opposite-direction accumulation and its application

Grey forecasting model for CO 2 emissions: a Taiwan study

Application of optimized fractional grey model-based variable background value to predict electricity consumption

Application of a novel fractional grey prediction model with time power term to predict the electricity consumption of India and China

Predicting Chinese total retail sales of consumer goods by employing an extended discrete grey polynomial model

A nonlinear grey forecasting model with double shape parameters and its application

Prediction of CO 2 emissions in iran using grey and arima models

Application of a novel nonlinear multivariate grey Bernoulli model to predict the tourist income of China

The conformable fractional grey system model

Modeling CO 2 emissions from fossil fuel combustion using the logistic equation

Forecasting of CO 2 emissions, energy consumption and economic growth in China using an improved grey model

Grey GM(1,1,) model with time power and its application

Prediction and analysis of the three major industries and residential consumption CO 2 emissions based on least squares support vector machine in China

Prediction of carbon dioxide emissions based on principal component analysis with regularized extreme learning machine

Projections of Turkey's electricity generation and installed capacity from total renewable and hydro energy using fractional nonlinear grey Bernoulli model and its reduced forms

Using grey models for forecasting China's growth trends in renewable energy consumption

Modeling carbon emission trajectory of China, US and India

A weighted nonlinear grey Bernoulli model for forecasting nonlinear economic time series with small data sets

Modelling the nonlinear relationship between CO 2 emissions and economic growth using a PSO algorithmbased grey Verhulst model

Forecasting the industrial solar energy consumption using a novel seasonal GM(1,1) model with dynamic seasonal adjustment factors

Forecasting Chinese carbon emissions from fossil energy consumption using non-linear grey multivariable models

Influencing factors analysis and forecasting of residential energy-related CO 2 emissions utilizing optimized support vector machine

Grey system model with the fractional-order accumulation

A novel Grey Bernoulli model for short-term natural gas consumption forecasting

Forecasting shortterm renewable energy consumption of China using a novel fractional nonlinear grey Bernoulli model

A novel conformable fractional non-homogeneous grey model for forecasting carbon dioxide emissiocns of BRICS countries

Analysis of novel FAGM(1,1,) model to forecast health expenditure of China

The indirect carbon emission from household consumption in China between 1995-2009 and 2010-2030: A decomposition and prediction analysis

Forecasting annual electricity consumption in China by employing a conformable fractional grey model in opposite direction

Forecasting fuel combustion-related CO 2 emissions by a novel continuous fractional nonlinear grey Bernoulli model with grey wolf optimizer

Forecasting Chinese greenhouse gas emissions from energy consumption using a novel grey rolling model

An investigation into the relationship between China's economic development and carbon dioxide emissions

A MFO-based conformable fractional nonhomogeneous grey Bernoulli model for natural gas production and consumption forecasting

Forecasting Chinese carbon emissions using a novel grey rolling prediction model

The heterogeneous effects of urbanization and income inequality on CO 2 emissions in BRICS ecomomies: evidence from panel quantile regression