Next Article in Journal
Global Solar Radiation and Its Interactions with Atmospheric Substances and Their Effects on Air Temperature Change in Ankara Province
Previous Article in Journal
A Systematic Review of Agroecology Strategies for Adapting to Climate Change Impacts on Smallholder Crop Farmers’ Livelihoods in South Africa
Previous Article in Special Issue
Emission and Reduction of Air Pollutants from Charcoal-Making Process in the Vietnamese Mekong Delta
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Analysis and Modeling of the CO2 Series Emitted by Thirty European Countries

by
Alina Bărbulescu
Department of Civil Engineering, Transilvania University of Brașov, 5 Turnului Str., 500152 Brasov, Romania
Climate 2024, 12(3), 34; https://doi.org/10.3390/cli12030034
Submission received: 20 December 2023 / Revised: 25 February 2024 / Accepted: 27 February 2024 / Published: 29 February 2024

Abstract

:
In recent decades, an increase in the earth’s atmospheric temperature has been noticed due to the augmentation of the volume of gases with the greenhouse effect (GHG) released into the atmosphere. To reduce this effect, the European Union’s directives indicate the action directions for reducing these emissions, among which carbon dioxide (CO2) recorded the highest amount. In this context, the article analyzes the CO2 series reported in 1990–2021 by 30 European countries. The Kruskal-Wallis test rejected the hypothesis that the series comes from the same underlying distribution. The Anderson-Darling test rejected the normality hypothesis for seven series out of thirty, and Sen’s procedure found a decreasing trend slope only for 17 series. ARIMA models have been built for all individual series. Grouping the series (by the k-means and hierarchical clustering) provided the base for building the Regional series (RegS), which describes the CO2 pollution evolution over Europe. The advantage of this approach is to provide the synthetic image of the regional evolution of the CO2 emission volume (mt), incorporating information from 30 series (one for each country) in only one—RegS. It is also shown that selecting the number of clusters involved in building RegS and assessing their stability is essential for the model’s goodness of fit.

1. Introduction

The Intergovernmental Panel on Climate Change (IPCC) [1] shows that the GHG emissions from human activities reached the highest levels compared to the previous 800,000 years. The US Inventory [2] indicates that in 2021, in the USA, the GHG volume exceeded 6340.2 mil mt of CO2 equivalents, higher by 6% than in 2020 but lower by 17% than in 2005. The GHG emissions from the EU economy rose to 941 mt CO2 equivalents in the first quarter of 2023 [3].
Carbon dioxide, the GHG with the highest weight in the total GHG (followed by CO2, NH4, and N2H), results from natural and anthropic sources. The first category includes respiration, ocean release, and decomposition. The second is represented by burning fossil fuels, deforestation, industrial production (like the cement industry), and agricultural activities [4,5,6]. Transportation emitted about a quarter of the EU’s CO2 volume in 2019, of which more than two-thirds were produced by road transportation [4,5]. CO2 emissions from transportation have increased by 33.5% from 1990 to 2019 in UE [4].
According to the Intergovernmental Panel on Climate Change (IPCC) document [1], the main cause of temperature augmentation from the middle of the twentieth century has been the enormous amount of GHG emitted into the atmosphere. The Sixth Assessment Report of IPCC shows that to limit this increase, immediate measures must be implemented [7,8,9,10]. According to [11,12], natural sinks contribute to removing some GHG from the atmosphere. However, when the GHG volume is high, these gases remain and accumulate for decades or even hundreds of years (like in the case of CO2). Therefore, conditions to stimulate the natural sinks’ activity must be created, and artificial sinks should be designed to balance the GHG’s effect [12,13,14].
The experimental results show that people’s constant exposure to CO2 may harm their health [15,16,17,18,19]. The review of Jacobson et al. [20] and the articles on which it is based also highlight some effects of chronic exposure to CO2, for example, oxidative stress, diminishing cognitive abilities, kidney calcification, bone demineralization, etc. [21]. Moreover, the adverse effects of the CO2 excess on the ecosystems, plants, and animals are presented in [22,23]. In this context, Hadipoor et al. [24] indicate some measures to control and reduce CO2 emissions.
Despite the necessity of emissions’ reduction being understood [25,26,27,28,29,30], there exist differences among countries related to carbon emissions’ responsibility [31,32]. So, most studies in the field discuss the responsibility for embodied CO2 emissions in international trade [25,31]. Some scientists analyzed the results of the international conventions’ implementation [33], while others focused on the source attribution of GHG emissions [34].
Techniques like ARIMA for modeling the CO2 series from 1972 to 2013 in Bangladesh [35], exponential smoothing and MLP for forecasting the CO2 series in Pakistan [36], and the use of the Cobb-Douglas function to estimate the CO2 sectorial amount in the total CO2 emission in Indonesia [37] proved to be efficient for this purpose. A review of the methods utilized for emphasizing the correlations between the emissions of CO2, the consumption of energy, and economic growth is presented in [38]. A survey on approaches utilized for modeling the CO2 emissions from stationary sources, GIS, and economic assessment can be found in [39]. An overview of the CO2 globally averaged concentrations series since 1830 is performed in [40].
In the context of the climate change negotiations, which started with the adoption of the Convention on Climate Change in 1992, followed by the Kyoto Protocol in 1997, and the launching of the EU’s Emissions Trading System in 2005, the European Union has been a key player in the fight against climate change. Investigating the CO2 series evolution is essential for assessing the results of implementing the measures for reducing climate change and diminishing its environmental effects [40,41]. Therefore, different articles studied the CO2 series’ temporal [35,40,42] or spatial [43,44] evolution in different countries, but fewer performed spatiotemporal [45] analyses, mostly for countries in Asia. In this study, we propose such an approach for investigating the trend of the CO2 series recorded in thirty European countries in Europe. First, we perform a deep statistical analysis, then build models (utilizing the Box–Jenkins methodology) for the series recorded in each country, and test the hypothesis that all series have the same underlying distribution. The rest of this study is dedicated to concisely presenting the regional and temporal evolution of the CO2 emission series in Europe and incorporating the information collected into one representative series—the Regional series (RegS) — significant from a spatial viewpoint. RegS is constructed using a selection from the 30 series recorded in the studied countries based on two clustering algorithms—k-means and hierarchical clustering [46,47]. Based on our knowledge, this approach is new in the study field.
It is also shown that finding the optimal number of clusters using different selection algorithms and testing the clusters’ stability is essential for obtaining the best RegS. So, even if it appears at a superficial analysis that clustering is the research goal, it is only one of the most critical steps in the proposed algorithm.

2. Materials and Methods

2.1. Data Series

The analyzed data set consists of the CO2 net emissions (in mt) series recorded from 1990 to 2021 in the EU—27 countries, Iceland, Norway, and Switzerland. The data series was downloaded in an .xlxs form from the official site of the European Union [48] and represents the series sent by the 30 countries to UNFCCC and the EU GHG Monitoring Mechanism. The series has no gaps. The series was processed and represented in Figure 1 on a logarithmic scale for clarity of representation. The international abbreviations of the countries’ names are utilized in this work: AT—Austria, BE—Belgium, BG—Bulgaria, CH—Switzerland, CY—Cyprus, CZ—Czechia, DE—Germany, DK—Denmark, EE—Estonia, EL—Greece, ES—Spain, FI—Finland, FR—France, HR—Croatia, HU—Hungary, IE—Ireland, IS—Island, IT—Italy, LT—Lithuania, LU—Luxembourg, LV—Latvia, MT—Malta, NL—Netherland, NO—Norway, PL—Poland, PT—Portugal, RO—Romania, SE—Sweden, SI—Slovenia, and SK—Slovakia.

2.2. Methodology

The study is split into two parts. The first focuses on evaluating the evolution of the time series recorded in each country and yearly CO2 series. The second one concerns building RegS.

2.2.1. Study the Time Series Recorded in Each Country

The steps in the individual series analysis are the following:
(1)
Test the hypothesis that the series is Gaussian against the hypothesis that it is not Gaussian by using the Anderson–Darling (AD) test [49].
(2)
Test the homoskedasticity hypothesis against the heteroskedasticity by the Fligner–Killeen (KF) test [50] (since it is less sensitive to the normality violation) [51].
(3)
Perform the Mann–Kendall (MK) [52] to test the randomness hypothesis against the existence of a monotonic trend. If the null is rejected, the slope of a linear trend will be computed by Sen’s [53]. The following series were subject to this analysis:
(a)
The CO2 series recorded in each country (30 series);
(b)
The total CO2 emissions during 1990–2021.
(4)
Perform the Augmented Dickey–Fuller (ADF) test [54] to assess the existence of a unit root vs. the time series stationarity for the series from (3).
(5)
Perform the Kruskal–Wallis (K-W) test [55] to test whether the series in a group originates from the same distribution against the alternative that at least one comes from a different distribution. When the null hypothesis was rejected, the post-hoc Dunn’s test [56], with the adjustment proposed by Hochberg [57], was run.
(6)
Modeling the time series from (3) using the ARIMA technique.
Considering a time process ( X t ), t∈Z, denote by d X t the difference of d-th order of X t . ( X t ) is an autoregressive integrated moving average process, ARIMA (p,d,q), if:
Φ ( B ) d X t = Θ ( B ) ε t ,
where Φ and Θ are polynomials of p and q orders, respectively, with roots higher than 1, and ( ε t ), t∈Z is white noise [58].
An autoregressive process of order p, AR (p), is a particular case of ARIMA, with q = d = 0. A moving average process of order q, MA (q), is an ARIMA (0,0,q).
If, after performing the ADF test, the null hypothesis was rejected, then d = 0. Otherwise, one may differentiate the series and perform the ADF test again until the null is rejected. Otherwise, one may differentiate the series and perform again the ADF test until the null is rejected. Thus, d will be equal to the number of differences taken on the series when getting the null rejected. The p and q orders were selected by analyzing the charts of the autocorrelation and partial autocorrelation functions (ACF and PACF, respectively). The model with the lowest value of the Akaike criterion [59] was kept. The Ljung–Box [60] test was used to check that the residuals form white noise.
All tests were performed at a significance level of 0.05. A p-value less than 0.05 computed in a test (except Dunn’s) leads to rejecting the corresponding null hypothesis. In Dunn’s test, the null hypothesis rejection is performed if the p-value < α / 2 .

2.2.2. Building RegS

The algorithm run for building RegS from the series recorded in the 30 European countries has the following steps [61,62]:
1.
Determine the optimal number of clusters, k, perform the k-means and hierarchical clustering, and choose the best clustering.
The k-means and hierarchical clustering [63,64,65] were utilized to group the series recorded in the 30 countries. The advantages and disadvantages of these techniques are presented in [66], shortly.
The silhouette [67], elbow [68], and gap statistics methods [69] were employed to determine the optimum k for running the k-means algorithm. Sometimes, these techniques do not provide the same k; thus, the study was performed for each possibility, and the best one was selected using some criteria that will be presented in the following.
The ratio WSS (the within-cluster sum of squares) and BSS/TSS (the between-clusters sum of squares by the total sum of squares) were computed to determine the best clustering in the k-means algorithm. The higher the BSS/TSS is, the better is the clustering. A smaller WSS indicates better groupings, as well [70].
Hierarchical clustering is a method whose graphical output is a dendrogram that indicates the series’ hierarchy. The distance between elements is computed in the first stage, and the distance matrix is built. Different distance functions can be used for this purpose, including Euclidean, Manhattan, Hamming, Jaccard, etc. In this study, the first one was utilized.
Hierarchical clustering can be agglomerative or divisive. In the first one utilized here, a cluster is initially created for each element, then, successively, the groups are merged until only one cluster is obtained. The selection of the clusters to be merged at an intermediate step is performed after checking the distances between the couples of clusters. The most similar clusters (couples with the lowest distances between them) are merged. The cluster similarity is defined using different linkage methods, among which are “average”, “complete”, “median”, “ward.D”, “ward.D2”, etc. [65,71]. In the “average” (“complete”) method, the average (maximum) distances between pairs of elements (one from a cluster and the other one from another cluster) are returned. These two methods performed the best on the CO2 dataset, taking into account the values of the cophenetic correlation coefficient [72,73]. Values greater than 0.9 (between 0.8 and 0.9) indicate a very good (good) clustering quality, whereas less than 0.8 show poor clustering results. After obtaining the clusters, the average Jaccard (AvgJaggard) measures and associated instabilities are computed to verify whether the algorithm (k-means and hierarchical clustering, in this case) provided a satisfactory representation in different groups of the studied dataset. AvgJaggard > 0.85 indicates a high stability of the cluster. AvgJaggard in [0.6, 0.85) (lower than 0.6) shows that the cluster is stable (unstable) [74].
2.
Select the cluster with the highest number of elements, Clmax. When at least two clusters have this property, Clmax is the one with the smallest WSS.
3.
Compute RegS, whose elements are the averages of the series from Clmax. More precisely, the value assigned to year j is the mean of values recorded in the same year in the countries from Clmax.
4.
Evaluate the modeling errors for each series by subtracting the values of RegS from the recorded values.
5.
Estimate the RegS’s goodness-of-fit of by computing the mean absolute percentage error (MAPE). MAPE was chosen because it is a non-dimensional index that can be utilized for comparing different models.
The R 4.3.2 software (https://www.r-project.org/) was used to perform the study.

3. Results and Discussion

3.1. Analysis of the CO2 Time Series

The Anderson–Darling test rejected the null hypothesis for the series recorded in BE, DK, EE, ES, FR, HU, IE, IS, IT, LT, MT, NL, PL, PT, SI, SK, and the total CO2 series. The homoscedasticity hypothesis was rejected for the IS, IT, LT, LU, NO, and SI series.
The MK trend test applied to the total CO2 series obtained by summing all series values recorded in each country during 1990–2021 rejected the null hypothesis. The nonparametric Sen’s procedure provided a negative slope of −2.9589 × 10 7 for the period 1990–2021. The polynomial trend determined for 1990–2002 has the equation:
y t = 766992 t 3 + 5 × 10 9 t 2 5 × 10 9 t + 6 × 10 11   ( R 2 = 0.8144 )
and after 2003, it has the equation:
y t = 5.3788 × 10 7 t + 1.1127 × 10 11   ( R 2 = 0.9012 ) ,
where t is the time and y t is the value of the series at the moment t.
Table 1 presents the p-values computed in the MK test for the CO2 series recorded in each country during the study period. If the p-value is less than 0.05, it is accompanied by the sign plus or minus between the brackets, meaning that the slope of the linear trend computed by Sen’s method is positive or negative, respectively. Out of 30 series, the null hypothesis was rejected for 20 s. A positive trend was found for the AT, CY, and IS series, whereas a negative trend was estimated for 17 countries.
The ADF test could not reject the hypothesis of a unit root existence for any series at the significance level of 5%. Therefore, there is not enough evidence for the series stationarity. Therefore, to reach stationarity, some transformations must be employed.
Table 2 contains the modeling results, as follows:
-
Column 2—the model type;
-
Column 3—the model’s coefficients (when a simple differentiation did not lead to the model) and the corresponding standard error (se) inside the brackets;
-
Column 4—the drift, if it exists, and the standard error (se) of its estimation inside the brackets;
-
The p-value computed in the Ljung–Box test applied to the model’s residual series;
-
The MAPE of the model.
The results from Table 2 indicate that all models except for ES were validated by the Ljung–Box test. Still, relatively high MAPEs were found in those of LV and SE. Therefore, we must find a better approach to model the ES, LV, and SE series. It will be the subject of another article. The AT, BG, CH, CZ, DE, HU, IS, LT, MT, PL, SE, and SK series are described by ARIMA (0,1,0), which is obtained by taking the first-order difference of the raw series. AR(1) models were determined for EE, HR, NO, and PT, AR(2) for LU, and MA(1) for FI and LV. The rest of the series, excluding ES, have models with d = 1 or 2. These findings are in concordance with those of the ADF test.
The model of the total CO2 series is of ARIMA (0,1,0) with the drift = −3,481,622 (se = 18,455,798) and MAPE = 2.5848. The p-value in the Ljung–Box test is 0.7458.

3.2. Building RegS

The critical step of the algorithm used for building RegS is to determine the optimal number of clusters, k. Different values were found for k—two, three, or one—by the silhouette, the elbow-knee method (Figure 2), and gap statistics, respectively. Since the result provided by the last method is influenced by the outlier existence (Germany, in this case—Figure 1) and the small distances between some clusters [61], it was neglected; thus, the study was performed for the remaining alternatives.
The clusters resulting after running the k-means algorithm are shown in Figure 3a for k = 2 and in Figure 3b for k = 3, where the names of the countries are replaced, for simplicity, by numbers from 1 to 30 assigned from AT to SK (in an alphabetical order).
WSSs are 86.7687 and 197.9709, respectively, and BSS/TSS = 69.3% for k = 2. WSS is high for the second cluster that contains DE(7), FR(13), IT(18), and PL(25). WSS has significantly lower values for k = 3 (0.000, 27.1215, and 34.1752), and BSS/TSS = 93.4%. Thus, the algorithm provides a significantly higher separation between the clusters when k = 3. In this case, Germany (7) forms a single cluster, ES(11), FR(13), IT(18), NL(23), and PL(25) belong to the third cluster, and the rest of the countries form the second one. This output can be explained by the series’ characteristics, like the average emissions levels or the existence of a certain type of trend of the time series in the same cluster. Indeed, the CO2 emissions in Germany are much higher than those of the countries from the third cluster: ES, FR, IT, NL, and PL. Moreover, the last four mentioned series have a decreasing trend. Still, the groups’ stability must be verified before selecting the best division of the countries in different subsets. The average Jaccard values and the instabilities for the two clusterings are shown in Table 3.
When k = 2, the first cluster is stable, and the second one is highly stable. When k = 3, the first is stable, and the other two are highly stable. In both cases, the highest instability is that of the first cluster (that contains Germany). Hierarchical clustering was also formed to cross-validate the above grouping and select the best one.
The cophenetic coefficient was computed to choose the linkage method. The highest value (0.9723) was obtained when using the “average”, followed by the “complete” (0.9473), compared to only 0.7474 in the case of “ward.D2”, or 0.7291 in the case of “ward” methods. Therefore, the “average” and “complete” linkage were utilized. The clustering obtained is provided in Figure 4a,b for k = 2 and Figure 4c,d for k = 3.
For k = 2, both methods provided the same clusters. The second cluster contained all countries but Germany. When k = 3, apart from Germany, which forms a cluster, there is a slight difference between the other two clusters when using “average” and “complete” methods (Figure 4c,d). In the first case, ES(11), FR (13), IT(18), and PL(25) belong to a cluster, while in the second one, NL(23) is added.
The clustering obtained by “complete” linkage, containing 1, 24, and 5 elements, coincides with that given by the k-means algorithm (Figure 3b). They cross-validate each other. The average Jaccard indicators computed after bootstrapping are given in Table 4.
Considering the ratio BSS/TSS in the k-means algorithm, the average Jaccard values and instabilities in the k-means, and hierarchical clustering, the best choice is k = 3. Indeed, when using “average”, two clusters are highly stable, and one is stable. If “complete” is employed, two clusters are stable, and one is highly stable. Both linkage methods provide good results in terms of AvgJaccard and instability.
That obtained by the “average” is better since the average instability is slightly lower than in the care of the “complete” method. Still, since the results are comparable, for building RegS, both alternatives were used, and the comparison is presented. Table 5 contains the MAPE (%). The MAPEs corresponding to the “complete” and “average” are denoted by MAPE_c and MAPE_av, respectively.
The values of MAPE_c (MAPE_av) are in the following intervals: 20 (20) are lower than 100, one (zero) is between 100 and 200, 3 (2) are between 200 and 300, 3 (3) are between 300 and 500, and the rest are higher than 500. The series values from the first and third clusters and most from the second cluster are well estimated. The highest values correspond to the countries with the lowest emissions, among which CY and IS have increasing trends. Considering the average MAPE for all countries (MAPE_c = 298.88, MAPE_av = 342.44), and that the lower the MAPE, the better the fitting is, we can conclude that the best result was obtained using the first linkage method.
To show that building RegS using the cluster with the highest number of elements is the best choice, we also built such a series using the series from the first and last clusters. MAPE was denoted by MAPE_DE (when using the first cluster—DE series), MAPE_c_III (when using the series from the third cluster with the “complete” method), and MAPE_av_III (when using the series from the third cluster with the “average” method). Based on the results from Table 6, the first choice is the worst, as it overestimates all the values from the other countries.
Building RegS with the series from the third cluster provides significantly (sometimes more than ten times) higher MAPE than when using the second cluster for this purpose (see MAPE_c_III from Table 6 vs. MAPE_c from Table 5, and MAPE_av_III from Table 6 vs. MAPE_av from Table 5).
The regional series obtained in all cases are depicted in Figure 5, where RegS_c (RegS_av) is the series obtained by using the elements from the cluster with the highest number of elements and the “complete” (“average”) method, DE is the series from the first cluster, and Reg_III_c (Reg_III_av) is the series obtained by using the elements from the cluster third cluster and the “complete” (“average”) method.
Apart from the overall decreasing trend captured by the models, one may note short augmentation periods followed by abrupt decreases, which could be explained by the continuous monitoring and regulations that appeared and were applied to control and restrict CO2 emissions. Comparing the values from Table 5 and Table 6 with the built regional series (Figure 5), it results a drastic overestimation of the regional evolution of CO2 series over Europe. So, the best approach is that proposed in the methodology.
Comparing the results from Table 1 with those obtained by the first procedure for the total GHG series recorded in Europe during the same period [47], the following groups of countries are determined as follows:
(a)
Countries (HR, EE, FI, EL, IE, PT, SI, ES) for which the MK test did not reject H0 for both CO2 and GHG series, so no significant monotonic trend can be emphasized;
(b)
Countries for which H0 was rejected and both CO2 and GHG series have the same type of trend: negative (BE, BG, CH, CZ, DE, DK, FR, HU, IT, LT, LU, NL, PL, RO, SE, and SK) or positive (AT, CY, and IS);
(c)
Countries for which the CO2 series has a negative slope of the trend, but H0 was not rejected for the total GHGs series (MT).
(d)
Countries for which H0 was not rejected, but the GHGs series have a monotonic increasing (LV) or decreasing (NO) trend.
So, most CO2 and total GHG series trends have a similar pattern. All series evolution but one can be described by ARIMA models. For most series, the correlation between the CO2 and GHG series (reported by the same country) is above 0.973, except for the Netherlands (0.847) and Norway (0.853).

4. Conclusions

This article studied the CO2 series reported by 30 countries from 1990 to 2021 and built the RegS using an original algorithm. RegS shows a decreasing trend of CO2 emissions, whereas, at the individual level, measures should be applied in about half of the countries to achieve the goal of pollution decrement.
Statistical analysis of big data series from trustworthy sources provides the background for making scientifically based decisions on implementing measures for reducing the pollution from anthropic sources and mitigating environmental disasters.
An analogous study will be developed on the CO2 volume per capita or GDP to assess the relationship between the economy, society, and environment in an interdisciplinary framework. A probabilistic approach to correlations between such variables will provide an in-depth analysis of the causality relationships. Other problems to be addressed are (a) the uncertainty, delays, and inertia in the proposed models; (b) integrating the socio-economic factors into environmental models; and (c) the risk evaluation.
Future research should answer the questions concerning the most appropriate measures and best implantation methods for sustainable development, and how to motivate society members to support actions towards a cleaner production and environment.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be available on request from the author.

Conflicts of Interest

The present research work carries no conflicts of interest.

References

  1. IPCC 2014. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Available online: https://www.ipcc.ch/report/ar5/syr/ (accessed on 14 September 2023).
  2. Inventory of U.S. Greenhouse Gas Emissions and Sinks. Available online: https://www.epa.gov/ghgemissions/inventory-us-greenhouse-gas-emissions-and-sinks (accessed on 16 November 2023).
  3. EU Economy Greenhouse Gas Emissions: −3% in Q1 2023. Available online: https://ec.europa.eu/eurostat/web/products-eurostat-news/w/ddn-20230816-1#:~:text=In%20the%20first%20quarter%20of,of%20CO2%2Deq (accessed on 14 September 2023).
  4. CO2 Emissions from Cars: Facts and Figures (Infographics). Available online: https://www.europarl.europa.eu/news/en/headlines/society/20190313STO31218/co2-emissions-from-cars-facts-and-figures-infographics (accessed on 7 August 2023).
  5. IEA. CO2 Emissions from Fuel Combustion; IEA: Paris, France, 2019. [Google Scholar] [CrossRef]
  6. Greenhouse Gas Emissions. Sources of Greenhouse Gas Emissions. Available online: https://www.epa.gov/ghgemissions/sources-greenhouse-gas-emissions (accessed on 7 July 2023).
  7. IPCC. Climate Change 2022: Mitigation of Climate Change; Cambridge University Press: Cambridge, UK, 2022. [Google Scholar]
  8. Paris Agreement. In Report of the Conference of the Parties to the United Nations Framework Convention on Climate Change (21st Session, 2015: Paris); HeinOnline: Paris, France, 2015.
  9. Calotă, R.; Antonescu, N.N.; Stănescu, D.-P.; Năstase, I. The Direct Effect of Enriching the Gaseous Combustible with 23% Hydrogen in Condensing Boilers’ Operation. Energies 2022, 15, 9373. [Google Scholar] [CrossRef]
  10. Antonescu, N.N.; Stănescu, D.-P.; Calotă, R. CO2 Emissions Reduction through Increasing H2 Participation in Gaseous Combustible—Condensing Boilers Functional Response. Appl. Sci. 2022, 12, 3831. [Google Scholar] [CrossRef]
  11. Denman, K.L.; Brasseur, G.; Chidthaisong, A.; Ciais, P.; Cox, P.M.; Dickinson, R.E.; Hauglustaine, D.; Heinze, C.; Holland, E.; Jacob, D.; et al. Couplings Between Changes in the Climate System and Biogeochemistry. In Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2007. [Google Scholar]
  12. Saini, M.; Saini, S.K. Identification of Major Sinks of Greenhouse Gases. In Greenhouse Gases: Sources, Sinks and Mitigation; Sonwani, S., Saxena, P., Eds.; Springer: Singapore, 2022; pp. 39–62. [Google Scholar]
  13. Greenhouse Gas Removal Methods and Their Potential UK Deployment. 2021. Available online: https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/1026988/ggr-methods-potential-deployment.pdf (accessed on 14 September 2023).
  14. Ming, T.; de Richter, R.; Oeste, F.D.; Tulip, R.; Caillol, S. A nature-based negative emissions technology able to remove atmospheric methane and other greenhouse gases. Atmos. Pollut. Res. 2021, 12, 101035. [Google Scholar] [CrossRef]
  15. Naiyer, S.; Abbas, S.S. Effect of Greenhouse Gases on Human Health. In Greenhouse Gases: Sources, Sinks and Mitigation; Sonwani, S., Saxena, P., Eds.; Springer: Singapore, 2022; pp. 85–106. [Google Scholar]
  16. Hajtar, L.; Herczeg, L. Influence of carbon-dioxide concentration on human well-being and intensity of mental work. Időjárás 2012, 116, 145–169. [Google Scholar]
  17. Satish, U.; Mendell, M.J.; Shekhar, K.; Hotchi, T.; Sullivan, D.; Strufert, S.; Fisk, W.J. Is CO2 an indoor pollutant? Direct effects of low-to-moderate CO2 concentrations on human decision-making performance. Environ. Health Perspect. 2012, 120, 1671–1677. [Google Scholar] [CrossRef] [PubMed]
  18. Allen, J.G.; MacNaughton, P.; Satish, U.; Santanam, S.; Vallarino, J.; Spengler, J.D. Associations of cognitive function scores with carbon dioxide, ventilation, and volatile organic compound exposures in office workers: A controlled exposure study of green and conventional office environment. Environ. Health Perspect. 2016, 124, 805–812. [Google Scholar] [CrossRef] [PubMed]
  19. Snow, S.; Boyson, A.S.; Paas, K.H.W.; Gough, H.; King, M.-F.; Barlow, J.; Noakes, C.J.; Schraefel, M.C. Exploring the physiological, neurophysiological and cognitive performance effects of elevated carbon dioxide concentrations indoors. Build. Environ. 2019, 156, 243–252. [Google Scholar] [CrossRef]
  20. Jacobson, T.A.; Kler, J.S.; Hernke, M.T.; Braun, R.K.; Meyer, K.C.; Funk, W.E. Direct human health risks of increased atmospheric carbon dioxide. Nat. Sustain. 2019, 2, 691–701. [Google Scholar] [CrossRef]
  21. Smith, D. Is Carbon Dioxide Harmful to People? Available online: https://learn.kaiterra.com/en/air-academy/is-carbon-dioxide-harmful-to-people (accessed on 15 September 2023).
  22. Hunter, P. The impact of CO2. The global rise in the levels of CO2 is good for trees, bad for grasses and terrible for corals. EMBO Rep. 2007, 8, 1104–1106. [Google Scholar] [CrossRef] [PubMed]
  23. Ehleringer, J.R.; Cerling, T.; Dearing, M.D. A History of Atmospheric CO2 and Its Effects on Plants, Animals, and Ecosystems; Ecological Studies 177; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  24. Hadipoor, M.; Keivanimehr, F.; Baghban, A.; Reza Ganjali, M.; Habibzadeh, S. Chapter 24—Carbon dioxide as a main source of air pollution: Prospective and current trends to control. In Sorbents Materials for Controlling Environmental Pollution; Núñez-Delgado, A., Ed.; Elsevier: Amsterdam, The Netherlands, 2021; pp. 623–688. [Google Scholar]
  25. Wang, T.; Chen, Y.; Zeng, L. Spatial-Temporal Evolution Analysis of Carbon Emissions Embodied in Inter-Provincial Trade in China. Int. J. Environ. Res. Public Health 2022, 19, 6794. [Google Scholar] [CrossRef]
  26. Li, Y.L.; Chen, B.; Chen, G.Q. Carbon network embodied in international trade: Global structural evolution and its policy implications. Energy Policy 2020, 139, 11131. [Google Scholar] [CrossRef]
  27. Chiritescu, R.-V.; Luca, E.; Iorga, G. Observational study of major air pollutants over urban Romania in 2020 in comparison with 2019. Rom. Rep. Phys. 2024, 76, 702. [Google Scholar]
  28. Dumitru, A.; Olaru, E.-A.; Dumitru, M.; Iorga, G. Assessment of air pollution by aerosols over a coal open-mine influenced region in southwestern Romania. Rom. J. Phys. 2024, 69, 801. [Google Scholar]
  29. Bărbulescu, A.; Dumitriu, C.S.; Popescu-Bodorin, N. On the aerosol optical depth series in the Arabian Gulf region. Rom. J. Phys. 2022, 67, 814. [Google Scholar]
  30. Bărbulescu, A.; Dumitriu, C.S.; Popescu-Bodorin, N. Assessing atmospheric pollution and its impact on the human health. Atmosphere 2022, 13, 938. [Google Scholar] [CrossRef]
  31. Ullah, I.; Ali, S.; Shah, M.H.; Yasim, F.; Rehman, A.; Al-Ghazali, B.M. Linkages between Trade, CO2 Emissions and Healthcare Spending in China. Int. J. Environ. Res. Public Health 2019, 16, 4298. [Google Scholar] [CrossRef] [PubMed]
  32. Jakob, M.; Marschinski, R. Interpreting trade-related CO2 emission transfers. Nat. Clim. Chang. 2013, 3, 19–23. [Google Scholar] [CrossRef]
  33. Slechten, A.; Verardi, V. Assessing the Effectiveness of Global Air-Pollution Treaties on CO2 Emissions; Working Papers 64981625; Lancaster University Management School, Economics Department: Lancaster, UK, 2014. [Google Scholar]
  34. Lindenmaier, R.; Dubey, M.K.; Henderson, B.G.; Butterfield, Z.T.; Herman, J.R.; Rahn, T.; Lee, S. Multiscale observations of CO2, and pollutants at Four Corners for emission verification and attribution. Proc. Nat. Acad. Sci. USA 2014, 111, 8386–8391. [Google Scholar] [CrossRef] [PubMed]
  35. Budiono, R.; Juahir, H.; Mamat, M.; Sukono; Supian, S.; Nurzaman, M. Modeling and analysis of CO2 emissions in million tons of sectoral greenhouse gases in Indonesia. IOP Conf. Ser. Mat. Sci. Eng. 2019, 621, 012020. [Google Scholar] [CrossRef]
  36. Debone, D.; Pazini Leite, V.; Georges El Khouri Miraglia, S. Modelling approach for carbon emissions, energy consumption and economic growth: A systematic review. Urban Clim. 2021, 37, 100849. [Google Scholar] [CrossRef]
  37. Yousefi-Sahzabi, A.; Sasaki, K.; Djamaluddin, I.; Yousefi, H.; Sugai, Y. GIS modeling of CO2 emission sources and storage possibilities. Energy Proceed. 2011, 4, 2831–2838. [Google Scholar] [CrossRef]
  38. Bâki iz, H. The evolution of large-scale variations in globally averaged atmospheric CO2 concentrations since 1830. All Earth 2022, 34, 16–26. [Google Scholar]
  39. Hansen, J.; Sato, M.; Kharecha, P.; Beerling, D.; Berner, R.; Masson-Delmotte, V.; Pagani, M.; Raymo, M.; Royer, D.L.; Zachos, J.C. Target Atmospheric CO2: Where Should Humanity Aim? Available online: https://arxiv.org/ftp/arxiv/papers/0804/0804.1126.pdf (accessed on 15 September 2023).
  40. Rahman, A.; Mahmudul Hasan, M.M. Modeling and Forecasting of Carbon Dioxide Emissions in Bangladesh Using Autoregressive Integrated Moving Average (ARIMA) Models. Open J. Stat. 2017, 7, 560–566. [Google Scholar] [CrossRef]
  41. Tawiah, K.; Daniyal, M.; Qureshi, M. Pakistan CO2 Emission Modelling and Forecasting: A Linear and Nonlinear Time Series Approach. J. Environ. Public Health 2023, 2023, 5903362. [Google Scholar] [CrossRef] [PubMed]
  42. Chang, S.C.; Lee, C.T. Evaluation of the temporal variations of air quality in Taipei City, Taiwan, from 1994 to 2003. J. Environ. Manag. 2008, 86, 627–635. [Google Scholar] [CrossRef] [PubMed]
  43. Kuttippurath, J.; Peter, R.; Singh, A.; Raj, S. The increasing atmospheric CO2 over India: Comparison to global trends. iScience 2022, 25, 104863. [Google Scholar] [CrossRef] [PubMed]
  44. Zhou, Y.; Zhang, J.; Hu, S. Regression analysis and driving force model building of CO2 emissions in China. Sci. Rep. 2021, 11, 6715. [Google Scholar] [CrossRef]
  45. Jiang, P.; Gong, X.; Yang, Y.; Tang, K.; Zhao, Y.; Liu, S.; Liu, L. Research on spatial and temporal differences of carbon emissions and influencing factors in eight economic regions of China based on LMDI model. Sci. Rep. 2023, 13, 7965. [Google Scholar] [CrossRef]
  46. Bărbulescu, A.; Nazzal, Y.; Howari, F. Statistical analysis and estimation of the regional trend of aerosol size over the Arabian Gulf Region during 2002–2016. Sci. Rep. 2018, 8, 9571. [Google Scholar] [CrossRef]
  47. Bărbulescu, A.; Postolache, F. New approaches for modeling the regional pollution in Europe. Sci. Total Environ. 2021, 753, 141993. [Google Scholar] [CrossRef]
  48. EEA Greenhouse Gases—Data Viewer. Available online: https://www.eea.europa.eu/data-and-maps/data/data-viewers/greenhouse-gases-viewer (accessed on 15 June 2023).
  49. Anderson, T.W.; Darling, D.A. Asymptotic theory of certain “goodness of fit” criteria based on stochastic processes. Ann. Math. Stat. 1952, 23, 193–212. [Google Scholar] [CrossRef]
  50. Fligner, M.A.; Killeen, T.J. Distribution-free two-sample tests for scale. J. Am. Stat. Assoc. 1976, 71, 210–213. [Google Scholar] [CrossRef]
  51. Conover, W.J.; Johnson, M.E.; Johnson, M.M. A comparative study of tests for homogeneity of variances, with applications to the outer continental shelf biding data. Technometrics 1981, 23, 351–361. [Google Scholar] [CrossRef]
  52. Kendall, M.G. Rank Correlation Methods, 4th ed.; Charles Griffin: London, UK, 1975. [Google Scholar]
  53. Sen, P.K. Estimates of the regression coefficient based on Kendall’s tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  54. Fuller, W.A. Introduction to Statistical Time Series; John Wiley and Sons: New York, NY, USA, 1976. [Google Scholar]
  55. Kruskal, W.H.; Wallis, W.A. Use of ranks in one-criterion variance analysis. J. Am. Stat. Assoc. 1952, 47, 583–621. [Google Scholar] [CrossRef]
  56. Dunn, O.J. Multiple comparisons using rank sums. Technometrics 1964, 6, 241–252. [Google Scholar] [CrossRef]
  57. Hochberg, Y. A sharper Bonferroni procedure for multiple tests of significance. Biometrika 1988, 75, 800–802. [Google Scholar] [CrossRef]
  58. Brockwell, P.J.; Davis, R.A. Introduction to Time Series and Forecasting; Springer International Publishing: Cham, Switzerland, 2016. [Google Scholar]
  59. Akaike, H. Time series analysis and control through parametric models. In Applied Time Series Analysis; Findley, D.F., Ed.; Academic Press: New York, NY, USA, 1978. [Google Scholar]
  60. Ljung, G.M.; Box, G.E.P. On a Measure of a Lack of Fit in Time Series Models. Biometrika 1978, 65, 297–303. [Google Scholar] [CrossRef]
  61. Bărbulescu, A.; Postolache, F.; Dumitriu, C.Ș. Estimating the precipitation amount at regional scale using a new tool, Climate Analyzer. Hidrology 2021, 8, 125. [Google Scholar] [CrossRef]
  62. Bărbulescu, A. On the Regional Temperature Series Evolution in the South-Eastern Part of Romania. Appl. Sci. 2023, 13, 3904. [Google Scholar] [CrossRef]
  63. K-Mean: Getting the Optimal Number of Clusters. Available online: https://www.analyticsvidhya.com/blog/2021/05/k-mean-getting-the-optimal-number-of-clusters/ (accessed on 16 June 2023).
  64. Hierarchical Clustering in, R. Available online: https://www.datacamp.com/tutorial/hierarchical-clustering-R (accessed on 18 June 2023).
  65. Kassambara, A. Practical Guide to Cluster Analysis in R. Unsupervised Machine Learning. Available online: https://www.datanovia.com/en/product/practical-guide-to-cluster-analysis-in-r/ (accessed on 10 May 2023).
  66. Daburra, I. K-Means Clustering: Algorithm, Applications, Evaluation Methods, and Drawbacks. Available online: https://towardsdatascience.com/k-means-clustering-algorithm-applications-evaluation-methods-and-drawbacks-aa03e644b48a (accessed on 18 June 2023).
  67. Rousseeuw, P. Silhouettes: A Graphical Aid to the Interpretation and Validation of Cluster Analysis. J. Comput. Appl. Math. 1987, 20, 53–65. [Google Scholar] [CrossRef]
  68. Thorndike, R.L. Who belongs in the family? Psychometrika 1953, 18, 267–276. [Google Scholar] [CrossRef]
  69. Löhr, T. K-Means Clustering and the Gap-Statistics. Available online: https://towardsdatascience.com/k-means-clustering-and-the-gap-statistics-4c5d414acd29 (accessed on 16 June 2023).
  70. Madsen, B.E.; Browning, S.R. A Groupwise Association Test for Rare Mutations Using a Weighted Sum Statistic. PLoS Genet. 2009, 5, e1000384. [Google Scholar] [CrossRef]
  71. Sokal, R.R.; Rohlf, F.J. The comparison of dendrograms by objective methods. Taxon 1962, 11, 33–40. [Google Scholar] [CrossRef]
  72. Farris, J.S. On the cophenetic correlation coefficient. System. Zool. 1969, 18, 279–285. [Google Scholar] [CrossRef]
  73. Myrphy, P. Clustering Data in R. Available online: https://rstudio-pubs-static.s3.amazonaws.com/599072_93cf94954aa64fc7a4b99ca524e5371c.html (accessed on 30 June 2023).
  74. Soetewey, A. Stats and R. The Complete Guide to Clustering Analysis: K-Means and Hierarchical Clustering by Hand and in R. Available online: https://statsandr.com/blog/clustering-analysis-k-means-and-hierarchical-clustering-by-hand-and-in-r/ (accessed on 20 June 2023).
Figure 1. The annual CO2 series [mt] during 1990–2021, represented in a log10 scale.
Figure 1. The annual CO2 series [mt] during 1990–2021, represented in a log10 scale.
Climate 12 00034 g001
Figure 2. The optimal number of clusters determined by (left) the silhouette and (right) the elbow method for the CO2 series recorded in the European countries.
Figure 2. The optimal number of clusters determined by (left) the silhouette and (right) the elbow method for the CO2 series recorded in the European countries.
Climate 12 00034 g002
Figure 3. Results of the k-means algorithm for grouping the series recorded in the European countries when (a) k = 2 and (b) k = 3. Numbers from 1 to 30 are assigned to the countries in the alphabetical order of names’ abbreviations.
Figure 3. Results of the k-means algorithm for grouping the series recorded in the European countries when (a) k = 2 and (b) k = 3. Numbers from 1 to 30 are assigned to the countries in the alphabetical order of names’ abbreviations.
Climate 12 00034 g003
Figure 4. Hierarchical clustering (a) k = 2 and “average” method; (b) k = 2 and “complete” method; (c) k = 3 and “average” method, (d) k = 3 and “complete” method. The countries are numbered from 1—AT to 30—SK, in the alphabetical order of the countries’ names’ abbreviations.
Figure 4. Hierarchical clustering (a) k = 2 and “average” method; (b) k = 2 and “complete” method; (c) k = 3 and “average” method, (d) k = 3 and “complete” method. The countries are numbered from 1—AT to 30—SK, in the alphabetical order of the countries’ names’ abbreviations.
Climate 12 00034 g004
Figure 5. RegS for CO2 series.
Figure 5. RegS for CO2 series.
Climate 12 00034 g005
Table 1. The p-values in the MK test on the CO2 series reported by the European countries.
Table 1. The p-values in the MK test on the CO2 series reported by the European countries.
CountryATBEBGCHCYCZDEDKEEEL
MK p−val0.000 (+)0.000 (−)0.007 (−)0.002 (−)0.001 (+)0.000 (−)0.000 (−)0.000 (−)0.57030.116
CountryESFIFRHRHUIEISITLTLU
p−val0.7830.2120.000 (−)0.2500.000 (−)0.4460.000 (+)0.001 (−)0.003 (−)0.034 (−)
CountryLVMTNLNOPLPTROSESISK
p−val0.5700.043 (−)0.001 (−)0.8840.017 (−)0.9100.000 (−)0.000 (−)0.0620.000 (−)
Table 2. ARIMA-type models for the series recorded in each country.
Table 2. ARIMA-type models for the series recorded in each country.
ModelCoefficient (se)Drift (se)Ljung–BoxMAPE
ATARIMA (0,1,0) 0.0689.540
BEARIMA (1,1,0)ar1 = −0.5379 (0.1541)−823,323.2 (434,892.5)0.5852.528
BGARIMA (0,1,0) 0.69610.242
CHARIMA (0,1,0) 0.1633.565
CYARIMA (1,2,0)ar1 = −0.4796 (0.1660) 0.8723.118
CZARIMA (0,1,0) −1,672,886.6 (912,138.5)0.7213.382
DEARIMA (0,1,0) −13,175,314 (5,529,230)0.6632.801
DKARIMA (0,1,1)ma1 = −0.3653 (0.1724)−1,071,914.1 (530,997.1)0.8635.936
EEARIMA (1,0,0),mean = 18,424,007 (5,355,966)ar1 = 0.9194 (0.0809) 0.97815.576
ELARIMA (0,2,1)ma1 = −0.7087 (0.1233) 0.5783.755
ESARIMA (0,0,0), mean = 240,898,858 (7,123,938) 0.00013.679
FIARIMA (0,0,1) with mean = 30,004,148 (1,645,259)ma1 = 0.4819 (0.2230) 0.66219.416
FRARIMA (0,1,1)ma1 = −0.4832 (0.1950)−3,180,106 (1,285,215)0.5752.903
HRARIMA (1,0,0), mean = 12,922,689 (1,229,364)ar1 = 0.7948 (0.0809) 0.8898.736
HUARIMA (0,1,0) −923,762.9 (441,331.6)0.4293.588
IEARIMA (1,0,2) with mean = 44,879,052 (2,889,765)ar1 = 0.8256 (0.0983)
ma1 = 0.2435 (0.1239)
ma2 = 0.6980 (0.1712)
0.9972.913
ISARIMA (0,1,0) 37689.71 (24201.84)0.5011.227
ITARIMA (0,1,0) 0.9563.974
LTARIMA (0,1,0) 0.25627.837
LUARIMA (2,0,0) with mean = 9,754,522.1 (622,592.4)ar1 = 1.1611 (0.1629)
ma1 = −0.3443 (0.1769)
0.6265.342
LVARIMA (0,0,1) with mean = 4,189,224.4 (525,627.8)ma1 = 0.5705 (0.1333) 0.586258.177
MTARIMA (0,1,0) 0.9346.524
NLARIMA (1,2,1)ar1 = −0.2340 (0.1926)
ma1 = 0.8333 (0.1073)
0.5672.651
NOARIMA (1,0,0), mean = 24,099,466 (1,405,598)ar1 = 0.6871 (0.1205) 9.482
PLARIMA (0,1,0) 0.7953.040
PTARIMA (1,0,0), mean = 49,480,302 (3,843,606)ar1 = 0.6729 (0.1330) 0.94510.135
ROARIMA (0,1,1)ma1 = 0.5272 (0.1465)−4,199,933 (2,188,561)0.56910.766
SEARIMA (0,1,0) 0.71455.183
SIARIMA (0,1,0) 0.3567.363
SKARIMA (0,1,0) −793,429.1 (483,098)0.7836.182
Table 3. The average Jaccard values and instabilities in the k-means clustering.
Table 3. The average Jaccard values and instabilities in the k-means clustering.
k = 2k = 3
ClustersAvgJaccardInstabilityClustersAvgJaccardInstability
10.75070.34310.64450.361
20.95680.00020.90610.007
30.85110.087
Table 4. The average Jaccard values and instabilities in the hierarchical clustering.
Table 4. The average Jaccard values and instabilities in the hierarchical clustering.
k = 2 k = 3
ClustersAvgJaccardInstabilityClustersAvgJaccardInstability
10.94530.00010.96850.003
Average20.63200.36820.63300.367
30.86070.100
10.94790.00010.95670.002
Complete20.65000.35020.65200.348
30.79940.161
Table 5. MAPE (%) in the RegS modeling built using the series form the second cluster. The cluster to which the country belongs is marked inside the bracket, with Roman numerals: I = 1, II = 2, and III = 3.
Table 5. MAPE (%) in the RegS modeling built using the series form the second cluster. The cluster to which the country belongs is marked inside the bracket, with Roman numerals: I = 1, II = 2, and III = 3.
CountryAT(II)BE(II)BG(II)CH(II)CY(II)CZ(II)DE(I)DK(II)EE(II)EL(II)
MAPE_c31.3767.3210.488.98461.9868.0995.6129.29162.0556.28
MAPE_av22.4562.4617.157.32544.3963.3394.9620.99201.1549.78
CountryES(III)FI(II)FR(III)HR(II)HU(II)IE(II)IS(II)IT(III)LT(II)LU(II)
MAPE_c84.3035.6789.35206.6827.7020.23322.7190.55460.03286.51
MAPE_av81.9853.3287.76252.0416.8713.55384.9889.14545.67343.94
CountryLV(II)MT(II)NL *NO(II)PL(III)PT(II)RO(II)SE(II)SI(II)SK(II)
MAPE_c3479.831515.6178.7359.1287.7024.2637.91790.94267.0313.25
MAPE_av3989.871756.4375.5882.6385.8816.6936.49926.02321.1129.31
* Note: NL belongs to the cluster II (III) when using complete (average) method.
Table 6. MAPE (%) when building regional series using the cluster I (DE) and III.
Table 6. MAPE (%) when building regional series using the cluster I (DE) and III.
CountryAT(II)BE(II)BG(II)CH(II)CY(II)CZ(II)DE(I)DK(II)EE(II)EL(II)
MAPE_DE1470.27645.882185.871999.8212,765.59628.3801514.205862.92899.55
MAPE_c_III495.42183.73773.61698.864755.16177.6961.84515.072191.58279.10
MAPE_av_III440.63157.82693.56625.744308.59152.3465.33459.301981.18244.53
CountryES(III)FI(II)FR(III)HR(II)HU(II)IE(II)IS(II)IT(III)LT(II)LU(II)
MAPE_DE259.422948.87143.006927.441550.891771.969549.77116.0412,757.8981,008.29
MAPE_c_III35.951054.967.642560.80529.14608.883564.4917.934820.683261.18
MAPE_av_III23.53948.9315.942316.76471.88543.933227.4725.404376.742953.39
CountryLV(II)MT(II)NL *NO(II)PL(III)PT(II)RO(II)SE(II)SI(II)SK(II)
MAPE_DE36,791.50385.233523.23180.291632.331504.6619,918.228266.202465.8636,791.50
MAPE_c_III31,259.2813,915.3884.611283.5610.63557.26511.897706.143090.75878.89
MAPE_av_III28,253.6012,639.5567.681156.2310.35497.13457.156998.472796.41789.56
* Note: NL belongs to the cluster II (III) when using complete (average) method.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bărbulescu, A. Statistical Analysis and Modeling of the CO2 Series Emitted by Thirty European Countries. Climate 2024, 12, 34. https://doi.org/10.3390/cli12030034

AMA Style

Bărbulescu A. Statistical Analysis and Modeling of the CO2 Series Emitted by Thirty European Countries. Climate. 2024; 12(3):34. https://doi.org/10.3390/cli12030034

Chicago/Turabian Style

Bărbulescu, Alina. 2024. "Statistical Analysis and Modeling of the CO2 Series Emitted by Thirty European Countries" Climate 12, no. 3: 34. https://doi.org/10.3390/cli12030034

APA Style

Bărbulescu, A. (2024). Statistical Analysis and Modeling of the CO2 Series Emitted by Thirty European Countries. Climate, 12(3), 34. https://doi.org/10.3390/cli12030034

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop