Climate Change-Greenhouse Gas Emissions Analysis and Forecast in Romania

: Greenhouse gases (GHG), such as carbon dioxide, methane, nitrous oxide, and other gases, are considered to be the main cause of global climate change, and this problem has received signiﬁcant global attention. Carbon dioxide has been considered the most signiﬁcant gas contributing to global climate change. Our paper presents an analysis of the greenhouse gas emissions in Romania along with a forecast for the years to come. For the study, data from the National Institute of Statistics and Eurostat were gathered and used for the analysis in order to present the results. To obtain the results, the data gathered were analyzed using forecasting methods that can be of help in solving some uncertainties that surround the future. The greenhouse gas (GHG) emissions trends in Romania were analyzed both for linear and exponential function methods. The obtained results showed that the linear function analysis of total GHG emissions in Romania had a forecast accuracy higher than the exponential function method. From the analytical methods used we can draw the conclusion that the emissions are on a descending scale and choosing a proper method is important in analyzing data.


Introduction
The biggest threat to the environment and also to a country's economic situation comes from the climate changes that occur. An analysis over the past 150 years revealed an increase in median temperature by one degree Celsius [1]. The highest temperatures ever recorded since 1850 until now happened after 1995. It is estimated that the annual median temperature by 2100 will increase by 4 to 5 degrees Celsius [2,3]. The Fourth Assessment Report by Working Group I of the Intergovernmental Panel on Climate Change (IPCC) estimated that global climate change is forecasted to result from six scenarios of GHG aerosol emissions, taken from a larger set of scenarios from the IPCC Special Report on Emissions Scenarios. Working Group I assessed climate change under these scenarios, including ocean-atmosphere general circulation models (GCMs) and simple climate models (SCMs), with some models that include carbon-cycle feedback for the climate-change projections. These scenarios are identified as A1FI, A1B, A1T, A2, B1, and B2, with the lowest cumulative emissions being projected by the B1 scenario and the highest by A1FI. They estimated that the A1FI emissions scenario would lead to a warming of 4 • C in the 2070s relative to pre-Industrial temperatures [4].
According to Liu and Raftery [5], based on a fully statistical probabilistic model for forecasting future fossil fuel and industry carbon emissions, the probability of staying below a warming of 2 • C is only 5%.
The impact of climate change can be seen mostly through extreme climate events like tornados, heatwave, floods, etc., all happening more and more often both in Europe as well as in other parts of the world [6]. Natural ecosystems, water reserves, and human health are highly affected in each part of the globe. For example, in Europe, rains have increased in the north/northwest side. The economic sector, once developed in accordance with the region and the climate, now will be affected and produce suffering for the natural ecosystem, with devastating consequences, some implying the loss of biodiversity.
Sustainable climate change systems in the context of decreasing environmental impact is a key subject included in the broader topic of environmental protection in Europe. The Sustainable Development Report 2019 mentions the need to take urgent action to combat climate change and its impact [7]. In this respect, the European Commission created the legal framework and policies for climate and energy from 2020 to 2030-COM (2014), which set one important target, reducing greenhouse emissions by 2050 to under 80-85%, which is actually a level below the level recorded in 1990 [8]. In addition, the new targets for Europe set by the EU mentions the point of reaching a climate-neutral continent by 2050, a goal set in the European Green Deal. The related legislation is binding and involves a commitment to climate neutrality [9].
The climatic sensitivity of global wildfire activity was analyzed by researchers [10]. Using satellite-derived wildlife products and processes, they were able to present the climate factors that influence wildfires. Process-oriented fire models helped to simulate past and future wildfires. Thus, their results present the spatiotemporal dynamics of burned areas along with their possible response to climatic factors [10].
Different areas around the world were analyzed in order to highlight the features of regional climate change. King et al. [11] chose to analyze the glaciers draining from the Geladandong ice caps since the 1960s. Fu et al. [12] presented the results related to climate change in the Tibetan Plateau using global climate change model simulation.
Surface air temperature and precipitation are the climatic factors that are usually analyzed [13]. Global temperature trend analysis was the main topic in several scientific articles [14][15][16]. Mudelsee M. [16] used a system to generate a time series of global surface temperature.
The earth system is influenced highly by living organisms because they impact the formation of clouds and atmospheric dust and alter the concentration of GHG in the atmosphere [17][18][19]. The World Meteorological Organization specified that the period between the years 2011 and 2020 was the warmest decade on record [20]. According to scientists and environmental experts, to be able to stop this level of change, to stop biodiversity change, it is necessary to reduce GHG emissions along with applying at a worldwide level a proper legal framework and directive that contains strategies and actions for the future. The difference between being able to fulfil the legal framework set in place and moving ahead as the situation is now will be determined by the financial capacity of each country.
Regarding the studies conducted that analyzed GHG emissions, many concluded that the main responsible factor for climate change disaster is GHGs of anthropogenic origin [21][22][23][24].
According to the Kyoto Protocol, GHGs include six gases, namely, carbon dioxide (CO 2 ), methane (CH 4 ), nitrous oxide (N 2 O), hydrofluorocarbons (HFCs), perfluorocarbons (PFCs), and sulfur hexafluoride (SF6) [25]. The most common GHG emitted by human activities is carbon dioxide. Anthropogenic GHGs include carbon dioxide from burning fossil fuels, deforestation, land clearing for agriculture, degradation of soils, methane from agricultural activities, waste management, energy use and biomass burning, nitrous oxide from agricultural fertilizers, and fluorinated gases from industrial processes and refrigeration [26,27].
In State of the Global Climate 2020 report it is mentioned that globally averaged mole fractions of carbon dioxide (CO 2 ) have already exceeded 410 parts per million (ppm), and if the CO 2 concentration follows the same pattern as in previous years, it could reach or exceed 414 ppm in 2021 [20].
Alhindawi et al. [28] investigated the GHG emissions from the transportation sector and developed a multivariate linear regression and double exponential smoothing model for the projection of GHG emissions from this sector.
Andrejiova et al. [29] presented in their research the current situation regarding GHG emissions in the air in 27 European Union (EU-27) member states. The authors compared the European Union member states on the basis of 15 examined parameters (population, area, the number of commercial airports with more than 15,000 passengers per year, the total GHG emissions, the amount of aviation GHGs, the amounts of individual aviation GHG emissions (CO 2 , N 2 O, CH 4 , HFC), the total amount of PM2.5 particulate matter in the air, etc.) by applying the principal component analysis method and a cluster analysis. They reported that in the years 2012 to 2018 in almost all EU-27 member states, there was an increase in aviation GHG emissions.
In another study, Saleh et al. [30] investigated CO 2 emissions using a prediction support vector machine (SVM) model, which considered variable energy consumption that has an impact on the emergence of CO 2 emissions. The authors considered electrical energy and burning coal as energy consumption that directly affects increasing CO 2 emissions. The trial-and-error approach was applied in order to obtain a better prediction model with lower error.
Li et al. [31] proposed a new forecasting model for carbon emissions related to energy consumption in the Beijing-Tianjin-Hebei region from 2017 to 2030 based on the extreme learning machine (ELM) algorithm optimized by grey prediction theory and a support vector machine (SVM). They proved that the SVM-ELM model has higher prediction accuracy than the SVM model and ELM model through the analysis of empirical research. In the same way, Zhou et al. [32] proposed a combined forecasting model based on the rough set and grey prediction models as well as the support vector machine (SVM) model to forecast carbon emissions from 2012 to 2017. They used the data on Chinese carbon emissions from 1992 to 2011.
Another study that used the support vector machine method for the prediction of CO 2 was conducted by Sun et al. [33]. The authors presented a method based on principal component analysis (PCA) and an improved least squares support vector machine (LSSVM) prediction.
Lee et al. [34] analyzed the uncertainty of GHG emissions using the parametric Monte Carlo simulation method and the non-parametric bootstrap method, whereas Akyol and Ucar estimated the greenhouse gas emissions of Turkey in the year 2030 using a time series forecasting algorithm in the WEKA (Waikato Environment for Knowledge Analysis) data mining software [35]. Kijewska and Bluszcz [36] analyzed the level of differentiation of European Union member states in terms of the emissions of four greenhouse gases and identified groups of similar countries based on these criteria. The authors used a taxonomic method-cluster analysis, namely, the agglomerative algorithm, which enables the extraction of objects that are similar to each other from the data and then merges them into groups. Their conclusion was that among the largest emitters of greenhouse gases are Germany, the United Kingdom, France, Turkey, Poland, Italy, and Spain.
In their research, Ding et al. [37] explored the dynamic relationship between business growth and carbon emissions performance by constructing and using a time series model to predict the trend of carbon emissions. They used the time series method (ADF unit root test, cointegration test, and VAR model) to investigate 805 companies listed on the Taiwan Stock Exchange from 2012 to 2017. Their results showed a long-term dynamic relationship between business growth and carbon emissions performance. Tubiello et al. [38] assessed the GHG emissions from drained organic sols. To better highlight the emissions, they were able to organize the data in two important domains-land use and agriculture-and then just the N 2 O and CO 2 emissions per unit area were used for the analysis. Their conclusion showed a significant presence of CO 2 in the emissions from drained organic soils. GHG emissions were also analyzed from a positive point of view by seeing whether there is a connection between agricultural GHG emissions and productivity growth [39]. The analysis continued and analyzed GHG emissions in correlation to renewable energy [40] and non-renewable energy. It was observed that the CO 2 emissions decreased in the case of renewable energy.
To summarize, the EU efforts and thus, Romanian efforts, focus on reducing GHG emissions, keeping in mind along the way that this process should not interfere with a country's economic state.
The present manuscript's objective was to conduct an analysis on GHG emissions in Romania in the EU context starting with the legislative framework in both the EU and Romania, as well as aspects of GHG emissions in the EU and Romania, and then conduct a forecast analysis of GHGs using data from Eurostat and the Romanian National Statistic Institute.
According to the 4th Biennial Report of Romania, the general trend of GHG emissions in Romania shows a strong decrease compared to the base year. Thus, between 1989 and 2018, total GHG emissions (excluding Land Use, Land-Use Change, and Forestry-LULUCF) decreased by 62.10% and net GHG emissions (including LULUCF) decreased by 68.32%. CO 2 had the largest share of total GHG emissions, followed by CH 4 and N 2 O. The emissions trends by sector in 2018 was as follows: Total GHG emissions from the energy sector accounted for the largest share (66.32%), followed by those from the agriculture sector, with a share of 17.1%, and those from industrial processes and the product-use sector, with a share of 11.58% [41].

Materials and Methods
In order to reach the research objective, we used data from the official statistics platforms Eurostat and INS (Romanian National Statistic Institute, Bucharest, Romania), focusing on the total GHGs emissions in Romania and the EU, GHG emissions by sector (energy, industrial process and product usage, waste and agriculture), GHG emissions per capita, and GHG emissions intensity of energy consumption. Our intention was to deepen the correlation between the two databases. The forecast analysis was done using MSOffice Excel statistics using the forecasting tools available. Forecasting theory was chosen based on the fact that the databases had past information that was of help for future prediction. Before applying the forecasting tool, the data were pre-processed because we observed missing data from some states in the EU.
Since we had a set timeline, year by year, a forecast analysis using chronological series was used. The chronological time series (SCr) is a set of values that a certain value can take (quantitative and qualitative) at different moment or during successive timeframes. In our case, the timeframes were successive and had the same length, so we used SCr of moments with equal timeframes between moments. Knowing this, the median indicators of the chronological time series were calculated.

Median level for SCr intervals
This can be calculated as a simple arithmetic average by the ratio between the sum of the series terms and the number of terms: where we have the following terms: n-number of terms; y t -the series terms; • SCr moments with equal intervals between moments-the median level is calculated as a chronological simple average, a particular case of the chronological weighted average where h 1 = h 2 = . . . = h n − 1 = h n = h and h represents the interval length (in time units).
Then we have: where: hi represents the length of the intervals between moments t i and t i+1 , i = 1, . . . , n − 1, expressed in time units. -Ai represents the simple moving chronological average (for example: A 1 = ((y 1 + y 2 ) × h)/2).
• Absolute average change (absolute average gain) This indicator can be determined as a simple average based on absolute modifications (changes) with a mobile base using the equation: y n -the last term of the SCr; y 1 -the first term of the SCr; n-number of terms.
This indicator/parameter shows with how many concrete unit measures the level of process is modified, on average, during the timeframe analyzed.

•
Dynamic average parameter (increase or decrease) This parameter can be determined as a geometrical average of the mobile base coefficients (parameters) using the following equation: This parameter shows how many times, on average, the level of the analyzed process has been modified considering the timeframe considered for calculating.

•
Dynamic rhythm (relative average change or relative median gain or median increase/ decrease ratio) This shows the percentage of increase/decrease, on average, of the process level from one period to another over the entire time horizon.
The analytical method-linear trend analysis and exponential trend analysis were used because they offer a more exact adjustment of the chronological series because they consider all the SCr terms.
The equation that describes the linear trend analysis is: Based on this, the parameters a and b of the function were calculated as follows: Regarding the exponential trend, the equation that best describes the analysis is: For the exponential trend after the logarithm of the relationship and applying the least squares method, we can determine parameters a and b as follows:

Brief Overview of GHG Emissions in the European Union
In order to achieve a sustainable future for the next generations, the EU is constantly developing methods and tools to monitor the progress of European countries in reducing the impact of GHG emissions. A reduction in the negative impact of GHG emissions can be achieved through an increase in public and industry awareness. The strategic EU documents that are considered to have an impact on GHG emissions are included in Table 1. Long-term low GHG emission development strategy of the EU and its member states "The endorsement of the climate neutrality objective was reached following an inclusive institutional and societal debate based on the strategic long-term vision proposed by the European Commission which includes a detailed analysis of solutions that could be pursued for the transition to a net zero GHG emissions economy and insights regarding the corresponding strategic priorities and an enabling framework that would allow reaching climate neutrality by 2050." https://unfccc.int/documents/210328, accessed on 14 October 2021 A Clean Planet for All-A European strategic long-term vision for a prosperous, modern, competitive, and climate-neutral economy "Research, development and demonstration will significantly reduce costs of breakthrough technologies. This will lead to genuinely new products replacing today's industrial products, such as carbon fibre or stronger cements reducing the volume of production while increasing product value. A net-zero GHG emissions economy will see new business concepts develop with re-use and additional services at its core." According to Eurostat and the scientific community, the GHG emissions resulting from human activity are the main cause of environmental decline, mainly Earth's average temperature increase over the last 250 years. The GHG emissions generated in 2019 in the EU were 8.4%, which was less than the GHG emissions in 2000, according to Eurostat. There was a period when the emissions slightly increased from 2002 to 2008, when the trend started to descend. The GHG emissions recorded by Eurostat presented the highest value for Luxemburg, with 30.8 metric tons per capita, in 2005, which in 2019 decreased to 20.3, whereas the smallest value recorded was in Latvia.
Romania is part of the new GHG emissions system that needs to be properly implemented in all of Europe. Efforts and significant investments have been made to align the local policies and practices to the European Community regulations. The situation is evolving and changes are made constantly.
In Figure 1, the total national emissions of the so-called "Kyoto basket" of GHGs are presented for the EU states. This includes carbon dioxide (CO 2 ), methane (CH 4 ), nitrous oxide (N 2 O), and the so-called F-gases (hydrofluorocarbons, perfluorocarbons, nitrogen trifluoride (NF3), and sulfur hexafluoride (SF6)). The EU member states report the emissions annually according to the United Nations Framework Convention on Climate Change (UNFCCC). As can be seen in Figure 1, there are still some differences among EU states. In Figure 1, the cumulative data are presented for each country, whereas for the European Union-27 countries, only the cumulative average is presented. Based on data reported by the member states, the EU climate action and the European Green Deal set a plan to cut emissions by 2030 by at least 55%. The ultimate goal is for Europe to become by 2050 the world's first climate-neutral continent. The initiatives and the measures put in place can have an impact if the GHG emissions are cut as the package of measures states and also by investing in research and innovation. In addition, if we compare the EU with other countries/continents, we can say that in 2020 the US with- Based on data reported by the member states, the EU climate action and the European Green Deal set a plan to cut emissions by 2030 by at least 55%. The ultimate goal is for Europe to become by 2050 the world's first climate-neutral continent. The initiatives and the measures put in place can have an impact if the GHG emissions are cut as the package of measures states and also by investing in research and innovation. In addition, if we compare the EU with other countries/continents, we can say that in 2020 the US withdrawal from the Paris Agreement and the projections were different from the EU. At the same time, China has adopted renewable energy much faster, thus leading to a decrease in GHG emissions [42].
GHG emissions from all operators covered by the European Union Emissions Trading Scheme (ETS) declined by 13.3% in 2020 compared to 2019. Additionally, in the energy sector the emissions decreased by 14.9% due to reduced electricity consumption along with using renewable energy that replaced the energy produced using coal. The industrial emissions also decreased by 7%. Considering the fact that in 2020 the COVID-19 pandemic started, the European Commission is not sure whether the emissions decrease can be attributed just to the regulation's efficiency.
Using the data available on Eurostat and applying Excel forecast tools, Figures 2 and 3 were created.
There is a significant decrease in GHG emissions by 2050, as can be seen in Figure 2. If the current plan for the European Union is implemented, then the emissions will decrease by 35% in 2030 and by 61% in 2050 compared to 2000. Additionally, if the analysis is done by taking into consideration 1990, the base year for many researchers, we could see a decrease of 60.25% by 2050, so just a 0.75% difference between the years 1990 and 2000. A great difference can be seen in Luxembourg, with less than 37% GHG emissions in 2030 and 76% in 2050 compared to 2000.
As can be seen from Figure 3, the countries that will change the most are Spain and Portugal from dark blue (values that are very high), and the country that actually will not move forward if the legislation, regulations, and implementation of GHG emissions reduction strategies do not change is Turkey [43]. We chose to add Turkey due to the fact that it demonstrates one of the worst scenarios, where regulations and the implementation of GHG reduction are not properly addressed. There is a significant decrease in GHG emissions by 2050, as can be seen in Figure 2. If the current plan for the European Union is implemented, then the emissions will decrease by 35% in 2030 and by 61% in 2050 compared to 2000. Additionally, if the analysis is done by taking into consideration 1990, the base year for many researchers, we could see a decrease of 60.25% by 2050, so just a 0.75% difference between the years 1990 and 2000. A great difference can be seen in Luxembourg, with less than 37% GHG emissions in 2030 and 76% in 2050 compared to 2000.
As can be seen from Figure 3, the countries that will change the most are Spain and Portugal from dark blue (values that are very high), and the country that actually will not   It is necessary to keep in mind that we are now in the fourth phase of EU ETS, which applies to 2021-2030 period, in correlation with the Commission Implementing Regulation (EU) 2018/2066 of 19 December 2018, aimed at monitoring and reporting GHG emissions. In addition, during the fourth phase the principles of Stage 3 of the EU ETS remain the same. The main potential value of the EU ETS is to make the price of GHG emissions visible throughout the production and consumption chain, thus paving the way for an efficient low-carbon economy. Only full bidding for emissions can determine their real price, as it is the necessary economic incentive to contribute to this transition. Since 2013, companies must bid in order to obtain an emissions share, shares that are fewer each year than the previous year. For example, in 2020 there were 70.9 million certified emissions reductions (carbon credits) for developing countries and those that are conducting clean development mechanism activities. There was a high number of requests for these credits.
In the same context, Dritsaki and Dritsaki [44] investigated the optimum model to forecast CO 2 emissions in the EU-28 based on annual data (from 1960 until 2014). The authors forecasted CO 2 emissions for six years (2015-2020) and used an autoregressive integrated moving average (ARIMA) (1,1,1)-autoregressive conditional heteroscedasticity (ARCH) (1) model combined with the linear ARIMA model and the conditional variance of the ARCH model. Their findings support the fact that the year 2020 presented a considerable decrease in CO 2 emissions, reaching 33.8% less than in the year 1990 (Kyoto Protocol) [44].

GHG Emissions Analysis and Forecasting in Romania
A total of 28 states submitted to the EU in 2020 their strategies for lowering GHG emissions in the long term. The ambitious goals for climate change were set while working with a global coalition of countries and international institutions. Despite all this, it seems that the EU needs better knowledge of the GHG emissions reduction plan and regulations.
Romania carries out an annual inventory of GHGs in order to comply with its European obligations and international regulations. Both the annual reporting obligations to the UNFCCC and the EU annual reporting requirements applied to all member states must be met. Specialized GHG services are being developed in Romania as a result for adoption of the Kyoto Protocol [45]. The aspects that are most often considered and can be used by industries of all size operations are: -Conformity verification, which is done under the international standard ISO 14065. This gives credit for a series of local, regional, and global emissions schemes. The companies that are developing in order to provide this kind of service usually have experience in EU ETS. -GHG emissions inventory and product verification against internationally accepted standards. -Carbon footprint-the first step for any business, which means making a verifiable carbon footprint for the organization or product and identifying opportunities to reduce costs by developing low-carbon business strategies and a roadmap for reducing GHG emissions.
Considering the data available on Eurostat, in Figure 4 we have the Romanian GHG emissions that present a descending slope from the first available data in 2000. It can be seen that since 2007, when the emissions had a value of 131,532.46 metric tons, the highest value reported by Romania, the GHG emissions decreased to 91,656.49 metric tons in 2008, which represents an improvement of 30%. Using the data available on Eurostat and applying Excel forecast tools, Figures 5 and  6 were developed. Since forecast accuracy is a key indicator, it is necessary to mention that at the end of the analysis scientists usually use relative approximation error, root mean squared error, mean absolute percentage error, or mean absolute deviation. Figure 5 shows the forecast and the boundary values for the confidence interval. To apply this method, a 95% confidence interval was selected in order to see the future values within this range (calculated as a normal distribution). The smoothing factor α in our forecast was 0.9, which indicates that the amount of smoothing was reduced, thus more weight could be placed on the more recent observations. The same forecast tool was also used for the analysis of the emissions intensity of the energy consumption, as shown in Figure 6. The symbol used was IGES and the data were measured in a thousand-metric-tons CO2 equivalent. As can be seen, the highest value  Using the data available on Eurostat and applying Excel forecast tools, Figures 5 and 6 were developed. Since forecast accuracy is a key indicator, it is necessary to mention that at the end of the analysis scientists usually use relative approximation error, root mean squared error, mean absolute percentage error, or mean absolute deviation. Using the data available on Eurostat and applying Excel forecast tools, Figures 5 and  6 were developed. Since forecast accuracy is a key indicator, it is necessary to mention that at the end of the analysis scientists usually use relative approximation error, root mean squared error, mean absolute percentage error, or mean absolute deviation. Figure 5 shows the forecast and the boundary values for the confidence interval. To apply this method, a 95% confidence interval was selected in order to see the future values within this range (calculated as a normal distribution). The smoothing factor α in our forecast was 0.9, which indicates that the amount of smoothing was reduced, thus more weight could be placed on the more recent observations. The same forecast tool was also used for the analysis of the emissions intensity of the energy consumption, as shown in Figure 6. The symbol used was IGES and the data were measured in a thousand-metric-tons CO2 equivalent. As can be seen, the highest value recorded happened in 2007, with a 90,370.00 metric ton CO2 equivalent, whereas the lowest value of 66,259.00 metric ton CO2 equivalent was registered in 2016. A 36% improvement can be seen in 2018 compared to 2007, but a 1% decrease compared to 2016. Although from 2000 to 2018 the values slightly decreased, the improvement in GHG emissions registered was 18%.   Figure 5 shows the forecast and the boundary values for the confidence interval. To apply this method, a 95% confidence interval was selected in order to see the future values within this range (calculated as a normal distribution). The smoothing factor α in our forecast was 0.9, which indicates that the amount of smoothing was reduced, thus more weight could be placed on the more recent observations. The same forecast tool was also used for the analysis of the emissions intensity of the energy consumption, as shown in Figure 6. The symbol used was I GES and the data were measured in a thousand-metric-tons CO 2 equivalent. As can be seen, the highest value this regulation and others that were issued since then contributed to decreasing G emissions due to their implementation among member states.  The data reported to Eurostat from 2000 to 2018 present that the energy sector the highest pollutant of GHGs. Their reported value never decreased under a 76,404 t sand metric ton CO2 equivalent, the highest value being recorded in 2003. The emis quantity decreases were visible mainly between 2007 and 2018, an improvement of The difference between the sectors is significant. The next sectors that were close in reporting numbers (the CO2 equivalent) were industrial processes and product usage agriculture.
As can be seen in Figure 7, their reported values were in the range of an app mately 19,000 metric ton CO2 equivalent. The highest value recorded in industrial cesses and product usage was 22,457.98 in the year 2004, whereas the lowest value 11,769 in the year 2013. Unfortunately, the decreasing trend was not consistently li with the correlation coefficient in this case being R 2 = 0.64. In 2018, the values repo were 14% higher than in 2013 and 30% lower than in 2000. In agriculture, 2018 reco higher emissions than 2010, when the lowest emissions were recorded. The increase by 13% compared to 2010, still 7% higher than 2000 but a 4% decrease when compar 2007. Since 2010, the values of GHG emissions continuously increased year by year a we do not change the regulations in this sector the emissions registered now will n much different in 2035 or 2050.
If we analyze the last sector in Figure 7, the waste sector, we can see that altho the values reported were lower than in previous sector, the trend was actually increa The values registered in 2018 were 5809 metric tons of CO2 equivalent and in 2000, metric tons of CO2 equivalent. As can be seen, a 9.5% increase in GHG emissions recorded. The highest value reported during this time was in 2017, whereas the lo was in 2011. The GHG emissions from waste in the EU fell by 42% since 1995, but a be seen from Romania's data, it is necessary to apply better waste management and ulations. By focusing on recycling and reducing disposal of waste as landfill, many member states were able to register lower values for emissions from waste. If Rom changes its waste treatment policies or applies the European regulations as it should, creating a more circular economy, then emissions will be reduced and climate change present a lower impact. The waste sector accounts for about 3% of total GHG emis [46].  The difference between the sectors is significant. The next sectors that were close in the reporting numbers (the CO 2 equivalent) were industrial processes and product usage and agriculture.
As can be seen in Figure 7, their reported values were in the range of an approximately 19,000 metric ton CO 2 equivalent. The highest value recorded in industrial processes and product usage was 22,457.98 in the year 2004, whereas the lowest value was 11,769 in the year 2013. Unfortunately, the decreasing trend was not consistently linear, with the correlation coefficient in this case being R 2 = 0.64. In 2018, the values reported were 14% higher than in 2013 and 30% lower than in 2000. In agriculture, 2018 recorded higher emissions than 2010, when the lowest emissions were recorded. The increase was by 13% compared to 2010, still 7% higher than 2000 but a 4% decrease when compared to 2007. Since 2010, the values of GHG emissions continuously increased year by year and if we do not change the regulations in this sector the emissions registered now will not be much different in 2035 or 2050.
If we analyze the last sector in Figure 7, the waste sector, we can see that although the values reported were lower than in previous sector, the trend was actually increasing. The values registered in 2018 were 5809 metric tons of CO 2 equivalent and in 2000, 5303 metric tons of CO 2 equivalent. As can be seen, a 9.5% increase in GHG emissions was recorded. The highest value reported during this time was in 2017, whereas the lowest was in 2011. The GHG emissions from waste in the EU fell by 42% since 1995, but as can be seen from Romania's data, it is necessary to apply better waste management and regulations. By focusing on recycling and reducing disposal of waste as landfill, many EU member states were able to register lower values for emissions from waste. If Romania changes its waste treatment policies or applies the European regulations as it should, thus creating a more circular economy, then emissions will be reduced and climate change will present a lower impact. The waste sector accounts for about 3% of total GHG emissions [46].

Chronological Series Analysis
As previously mentioned, the first analysis that was conducted was done by using time series techniques, specifically chronological series (SCr) analysis. The time series techniques were used to look back on patterns of the history of GHG emissions reporting; thus, if we can identify and project the pattern then we have our forecast. Chronological series component analysis can explain the past and present evolution of processes and phenomena but also offers the possibility to predict their behavior in the near or distant future through future projection of evolution models and past behavior. This implies a fundamental supposition that past models and trends that happened in the past will also be valid in the future. This assumption is closer to the truth for the immediate period, whereas for a longer period of time some factors might intervene that can modify the preestablished behavior. Despite the apparently random nature of variable manifestation within a timeframe, usually there is a basic pattern, either theoretical or fundamental, that can be identified, and that contains different elements and reflects the influence of some factors in such a manner that the original series can be reproduced as accurately as possible.
As said, the timeframes are successive and have the same length, so we used SCr of moments with equal timeframes between moments. Knowing this, the median indicators of the chronological time series were calculated.
Medium level for SCr intervals: Then the equation for SCr moments with equal intervals between moments is: The interval length is 1 year. The absolute average change (absolute average gain) result was: Regarding the dynamic average parameter (increase or decrease), the result was:

Chronological Series Analysis
As previously mentioned, the first analysis that was conducted was done by using time series techniques, specifically chronological series (SCr) analysis. The time series techniques were used to look back on patterns of the history of GHG emissions reporting; thus, if we can identify and project the pattern then we have our forecast. Chronological series component analysis can explain the past and present evolution of processes and phenomena but also offers the possibility to predict their behavior in the near or distant future through future projection of evolution models and past behavior. This implies a fundamental supposition that past models and trends that happened in the past will also be valid in the future. This assumption is closer to the truth for the immediate period, whereas for a longer period of time some factors might intervene that can modify the preestablished behavior. Despite the apparently random nature of variable manifestation within a timeframe, usually there is a basic pattern, either theoretical or fundamental, that can be identified, and that contains different elements and reflects the influence of some factors in such a manner that the original series can be reproduced as accurately as possible.
As said, the timeframes are successive and have the same length, so we used SCr of moments with equal timeframes between moments. Knowing this, the median indicators of the chronological time series were calculated.
Medium level for SCr intervals: y = ∑ n t=1 y t n = 110, 746.87 thousand tonnes CO 2 equivalent Then the equation for SCr moments with equal intervals between moments is: The interval length is 1 year. The absolute average change (absolute average gain) result was: Regarding the dynamic average parameter (increase or decrease), the result was: The last indicator of the chronologic series was the dynamic rhythm (relative average change or relative median gain or median increase/decrease ratio): Thus, based on the data used we can say that the average GHG emissions is the equivalent of 111,079.16 metric tons of CO 2 . The emissions decreased in the analyzed period by an equivalent of 1456.18 metric tons CO 2 on average in a year, which means 0.986 times. The dynamic rhythm result shows us that the GHG emissions only decreased by 1.4%. It is necessary to further analyze this dynamic rhythm in a few more years to see whether things have started to change for the better, by not just 1.4% but, if possible, by a higher percentage.

GHG Emissions Analytical Method Analysis
Analytical methods offer a more exact adjustment of the chronological series because they consider all the chronological series terms. These methods assume using the appropriate mathematical/analytical functions to determine the trend values. The fundamental problem is choosing the right type of function, also called the trend adjustment function. Most used functions are linear, polynomial, hyperbolic, exponential, logistics, and parabolic. Choosing one function can be done by graphic representation. The data gathered from Eurostat is a chronological series formed from an odd number of terms, so in this case the central point of the series was t = 0 and the rest of the terms were symmetrical to that one.
Using graphical representation, in Figures 8 and 9 we can see the GHG emissions trend in Romania both for linear and exponential analysis. As can be observed from the correlation coefficient, the linear trend analysis was a better fit because R 2 = 0.7557, unlike R 2 = 0.7376 for the exponential trendline analysis. The difference between the correlation coefficients was not significant, which is why this both methods are acceptable but still will need further analysis in the future. The last indicator of the chronologic series was the dynamic rhythm (relative average change or relative median gain or median increase/decrease ratio): Thus, based on the data used we can say that the average GHG emissions is the equivalent of 111,079.16 metric tons of CO2. The emissions decreased in the analyzed period by an equivalent of 1456.18 metric tons CO2 on average in a year, which means 0.986 times. The dynamic rhythm result shows us that the GHG emissions only decreased by 1.4%. It is necessary to further analyze this dynamic rhythm in a few more years to see whether things have started to change for the better, by not just 1.4% but, if possible, by a higher percentage.

GHG Emissions Analytical Method Analysis
Analytical methods offer a more exact adjustment of the chronological series because they consider all the chronological series terms. These methods assume using the appropriate mathematical/analytical functions to determine the trend values. The fundamental problem is choosing the right type of function, also called the trend adjustment function. Most used functions are linear, polynomial, hyperbolic, exponential, logistics, and parabolic. Choosing one function can be done by graphic representation. The data gathered from Eurostat is a chronological series formed from an odd number of terms, so in this case the central point of the series was t = 0 and the rest of the terms were symmetrical to that one.
Using graphical representation, in Figures 8 and 9 we can see the GHG emissions trend in Romania both for linear and exponential analysis. As can be observed from the correlation coefficient, the linear trend analysis was a better fit because R 2 = 0.7557, unlike R 2 = 0.7376 for the exponential trendline analysis. The difference between the correlation coefficients was not significant, which is why this both methods are acceptable but still will need further analysis in the future.
Thus, the equation that describes the linear trend analysis is: The values were calculated and are presented in Table 2.
The values were calculated and are presented in Table 2. As can be seen, the resulting equation was the linear trend equation for which the correlation coefficient is equal to 1. It can be observed that the analysis was correct because  Table 2. Linear trend analysis of total GHG emissions in Romania.

Years
y t t t 2 T*y tŷ t y t −ŷ t (y t −ŷ t ) 2 As can be seen, the resulting equation was the linear trend equation for which the correlation coefficient is equal to 1. It can be observed that the analysis was correct because the sum of the GHG emissions reported and the sum of the GHG emissions adjusted are equal. After calculating the relative approximation error for the linear trend analysis, we can say that the forecast accuracy is high because it is below 10%.
For the exponential method we used the equation above to determine the parameters, so they were: log a = ∑ log y t n = 95.76 19 = 5.04 (21) log b = ∑ t × log y t ∑ t 2 = −5.46 570 = −0.00957 (22) Thus, a = 109, 711.5 In this case the exponential equation will be: The values were calculated based on Table 3. Table 3. Exponential trend analysis of total GHG emissions in Romania.
Years y t t t 2 T*y tŷ t y t −ŷ t (y t −ŷ t ) 2 In addition, for the exponential trend analysis the relative approximation error was calculated, but the result was above 10%, which shows that this method is not accurate.
Considering the results, we can make a choice regarding the optimum method using the following criteria: - We select the graphic representation of the adjusted chronological series that best fits the real values; -We select the chronological series trend function for which the sum of the adjusted values is closest to the real values; -We select the adjusted function for which the sum of the squares of the differences between the adjusted and real values is minimum; - We select the function that has the minimum mean square deviation.
Considering all the criteria between the two analytical methods used for analyzing the total GHG emissions, we can say that the most exact method is the linear trend method. The criteria selection is also sustained by the relative approximation error that was calculated.
ε linear < ε exponential 5.44% < 12.32% Although the results were achieved using the basics tools for a statistic analysis, they revealed that if we do not apply the proper method for analysis then we will not be able to predict the emissions as close to what actually happens as we would like. For example, as seen from the relative approximation error, the exponential forecast leads to different results. The methods used in the paper were selected to be expressed in a simple way in order to be easily understood by scientists and interested people in the changes in GHG emissions.

Study Limitations
For a great quality analysis it is necessary to mention the limitations. Although the forecasting method used in the article is conventional and expressed good results, future research should consider combining both conventional and newly developed methods of forecast analysis. The recommendation to combine them comes from the fact that each method has its own advantages and limitations and only together can they present accurate results.
Significant results can be obtained using conventional forecast methods that use statistical modeling/analysis based on mathematical models that provide statistical diagnostic significance tests and using machine learning models that can accelerate data processing, automate forecasts, and increase adaptability to changes.
Due to the fact that the data collected had a clear period of time, the chronological analysis was a good fit, but to sustain the results, more combined methods should be considered because every analysis has its own limitations, including this method.
In this respect, in our future analysis we intend to analyze future data related to GHGs using IT tools, but for the start of our research we limited the analysis to the basics (forecasting tool from Excel).

Conclusions
Due to technological development in the past 30 years, the environment suffered changes that have had an impact on the environment in the long term, thus affecting the life of future generations. Because people's needs increased in this time, the GHG emissions that resulted from these processes also increased.
Thus, this paper highlights the GHG emissions and the importance of good emissions management. Based on the results of our research, it can be seen that statistics are a tool that can be used in predicting GHG emissions for the conditions that the current factors do not change. The regulations and policies that are adapted year by year influence any analysis. Environmental management and analysis tools are needed and should constitute an important part of the sustainable development plan. Assessing the situation based on known data will help stakeholders make better and more informed decisions.
A better system that can be used for calculating and monitoring GHG emissions is at the moment a requirement not just for Romania but for other countries as well. This conclusion can be drawn from the dynamic rhythm calculated that showed only a 1.4% improvement in the reduction of GHG emissions.
In addition, from the analytical methods used we can draw the conclusion that emissions are on a descending scale and choosing a proper method is important in analyzing data.
Regarding the results presented in the paper, it must be mentioned that one weakness that the forecast method should have addressed is using more analysis methods. Using more methods could have ensured a larger comparison between the results, thus resulting in a more accurate forecast. At the same time, the use of chronological series and the analytical methods analysis both showed a decrease in GHG emissions.
The forecasting tool from Excel showed that in the time analyzed, some improvements were made, and if we move forward in the future using the current regulations and policies, it is possible to achieve the Green Deal goal. Of course, the stakeholders' attitude towards changing and adapting to the EU directives will highly influence the future outcome.
Although in Europe there are states that are making considerable progress, there are still other countries like Romania that need support to implement the proper tools. Therefore, in order to see improvement in the numbers reported to the EU, in Romania it is necessary to: -Ensure transparency in the emissions quantity for each stakeholder, disregarding the sector; -Include better rights and obligations for stakeholders; -Provide the tools/means to change the current processes to sustainable, zero emissions ones; -Offer state grants for companies adapting to green energy and zero emissions; -Implement proper legislation for selecting, managing, and monitoring waste, using green energy for recycling; -Maintain a decrease of at least 1.4% in GHG emissions, as shown in our analysis; -Maintain a descending trend in the values GHG emissions generated and gathered, which are then reported to the EU; -Use the proper statistical methods to analyze the trends in GHGs emissions (from our research, linear trend analysis is one option); -Implement, on a national scale, carbon recovery machines or high-carbon-using trees in every part of the country.
In addition, the improvements mentioned before can be achieved only by conducting future scientific research on GHG emissions. One important indicator that should be analyzed is the carbon footprint, which is easy to understand since there is global interest in it and it is measured in quantitative units, but still, it is necessary to keep in mind that one weakness is the insufficient accuracy of the data. Future research should analyze each industry by itself and its GHG emissions, and from the results develop the proper tools for successful GHG emissions reduction.