{'fraction_non_alphanumeric': 0.03582972285262362, 'fraction_numerical': 0.02957155628911354, 'mean_word_length': 4.518785578747628, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Trend Analysis of NOX and SO2 Emissions in Indonesia from the Period of 1990 -2015 using Data Analysis Tool

A R T I C L E I N F O A B S T R A C T Article history: Received: 09 October, 2020 Accepted: 09 January, 2021 Online: 15 January, 2021 NOX and SO2 gas pollution have a direct impact on health problems and environmental damage. Therefore, to map the emission patterns and predict the resulting impacts, complete data and information on emissions of the two pollutants are needed. In Indonesia, data on NOX and SO2 emissions that are recorded over a long period of time (for example 5 decades) are very difficult to obtain. Meanwhile, REASv3.1 is a global emission inventory that provides complete data on air emissions in Asia during 1950 2015. Therefore, this study aimed to analyze NOX and SO2 emission trends, forecast data for 2016 2020, and compare the accuracy of calculations from the method used. The processing of both emission data used two methods, namely trend analysis based on exponential and polynomial approaches, and smoothing methods based on Double Moving Average (DMA) and Double Exponential Smoothing (DES). Furthermore, validation of the accuracy from both methods used the value of Mean Absolute Deviation (MAD), Root Mean Square Error (RMSE), and Mean Absolute Percentage Error (MAPE). The results showed that for the smoothing method, DMA was more accurate than DES. Meanwhile, the indicators are MAD, RMSE, and MAPE values, which are smaller and at a very good category. For forecasting results for 2016 – 2020, it was shown that the total emissions of both NOX and SO2 showed an increase, but with different gains. Furthermore, the total NOX emission gain was two times greater than the total of SO2. The road transportation and power plant sectors in NOX emissions showed an increasing trend, with an emission gain ratio of 3:20. Meanwhile, for SO2, the power plant sector experienced a significant increase, while the industrial sector actually showed a downward trend.


Introduction
Air quality is a measure of how much air is free from pollution and not harmful when inhaled by humans [1]. Air is stated to be polluted when it contains substances or energy, and other components that exceeds the quality standard [2]. Meanwhile, the sources of air pollution can come from natural processes or human activities (anthropogenic), and the air that is soiled by pollutants causes its quality to be poor. In fact, air pollution can have a serious impact on environmental damage, human health, and ecological balance [3]. This can be formed from chemical reactions with other pollutants and physical elements in the atmosphere.
According to The Environmental Protection Agency (EPA), 6 types of air pollutants can cause serious impacts on health and environmental damage, namely CO, NO2, Pb, PM, O3, and SO2 [4]. Therefore, monitoring air pollution is very necessary, because the data obtained can provide a lot of information about air quality in an area and within a certain time. REAS (Regional Emissions inventory in Asia) is an emission inventory owned by the Frontier ASTESJ ISSN: 2415-6698 Research Center for Global Change (FRCGC), Japan Agency for Marine-Earth Science and Technology (JAMSTEC), which provide data sets and various information about air emissions in the Asian region [5]. Furthermore, the REASv3.1 (latest version) provides complete anthropogenic emission data for the period of 1950 -2015.
The three anthropogenic emissions that causes respiratory problems, such as airway irritation, bronchitis, asthma, and pneumonia, are NOX, SO2, and PM [6][7][8]. The spread of these emissions is very broad, fast, and has a direct impact on health and the environment [9]. In Indonesia, complete and up to date NOX and SO2 emission data is rather difficult to obtain, because many ISPU stations (Air Pollutant Standard Index) that are tasked with monitoring air pollution are not operating properly. Therefore complete data on air emissions in Indonesia are mostly obtained on websites of global emission data providers, even until 2015. This data can be used to analyze air quality trends, and as a basis for forecasting data for 2016 -2020. Trend analysis is an empirical approach used to determine changes in the values of random variables, whether increasing or decreasing over a period of time in statistical terms. In addition to knowing future air quality conditions, forecasting air emission data is often used to anticipate risks in the event of exposure to poor air quality, as well as to formulate environmental pollution control strategies. [10,11].
Therefore, this research aimed to (1) analyze trends and forecasts of total NOX and SO2 emissions in Indonesia from 1950 -2015, and to compare the methods used to determine the bestperforming methods, (2) perform smoothing and forecasting of the dominant sector data from NOX and SO2, as well as compare the means used to find out which method has the best performance.

Materials and Methods
This study used data from REAS (Regional Emission inventory in Asia) ver.3.1. The data include 10 types of air pollutants (BC, CO, CO2, NH3, NOX, NMNVOC, OC, PM2.5, PM10, and SO2), as well as emission sources, both from the producing sector and the type of fuel used [12]. The data used are 2 types of pollutants (NOX and SO2), and the 2 largest emission-producing sectors of each pollutants types.
In this research, air emission data processing used two methods, namely the smoothing method with Double Moving Average (DMA) and Double Exponential Smoothing (DES) for emission data per sector. Furthermore, trend analysis was conducted with the exponential approach and polynomial orders of 2 and 3 for total emissions data. Meanwhile, to determine the accuracy of both methods, the results were validated using the calculation of MAD (Mean Absolute Deviation), RMSE (Root Mean Squared Error), and MAPE (Mean Absolute Percentage Error) [13].
where is the real data of each emission ̅ is forecasting data n is the amount of data used.
The focus of this study is in Indonesia, and the research methodology can be seen in Figure 1. To determine the MAPE accuracy criteria, it can be referred to Table 1   Table 1: The accuracy criteria of MAPE [14] Criteria The limit of MAPE percentage Very Good <10% Good 10% -20% Enough 20% -50% Not Accurate >50% All data were processed using the Analysis Toolpak, which is a set of "data analysis" tools in the processing group in Microsoft Excel.
The last stage was forecasting data for the next 5 years, starting from 2016 -2020. For DMA and DES, data forecasting was performed only for the best performance based on the values of MAD, RMSE, and MAPE that met the criteria. For trend analysis, data forecasting was carried out based on the line equation formed from the approach used.

Trend analysis of air pollutants (NOX, and SO2)
Estimation or forecasting of total emission data from NOX and SO2 pollutants in the coming year was carried out using trend analysis with three approaches, namely, exponential, polynomial order 2, and order 3. Furthermore, the selection of this approach was based on the suitability of the graphical patterns formed between real data and the 3 approaches.    This study used a lot of data obtained in the period of 1970-2015 or about 45 years. This relatively large data can provide accurate results to minimize errors. Furthermore, this study validated the results using MAD, RMSE, and MAPE to determine the best measurement accuracy of the approach used. Figure 2 showed a graph of trend analysis results on both pollutants, whereas Table 2. showed the performance evaluation results of both graphs Based on the data validity evaluation, the trend analysis that showed the best performance is the smallest MAD, RMSE, and MAPE values. For NOX emissions, the use of polynomial order 2 produced a good performance, while for SO2 emissions, it used a polynomial order 3. Based on the MAPE value, the accuracy of NOX emission measurements was in the very good category because the value was <10% (5.77%). Meanwhile, for SO2 emission, it was in a good category because the MAPE value was >10% (10.21%).
For the next stage, the research forecasted data for the next 5 years, between 2016 -2020, based on graphical equations that have been obtained from trend analysis that showed its best performance. NOX emission used a polynomial order 2 graph, while SO2 emission used order 3. The graph equation used for data forecasting was   Figure 3 showed that in Indonesia, there is an upward trend in emissions for both types of pollutants. This is different from what happened in developed countries, where NOX and SO2 emissions showed a downward trend, even though there was growth in the economic and industrial sectors [15]. The quantization graph of forecasting emissions of both pollutants in numerical data can be seen in Table 3, which showed the forecasting results of both pollutants for the period of 2016-2020. Meanwhile, SO2 emission is only 227.8 tons/year. This means that the NOX emission gain is more than two times the gain of SO2. Furthermore, the SO2 / NOX ratio has decreased from 0.878 in 2016 to 0.812 in 2020. This showed that the contribution of SO2 emissions to air pollution in Indonesia is not too significant compared to NOX

Smoothing and forecasting data of NOX and SO2 emissions
NOX and SO2 emission data in this article were collected from 1970 -2015, and were obtained from REAS version 3.1 [16]. Data smoothing is a method used to reduce randomness from time-series data in order to obtain relatively regular data [17,18].      To determine the data smoothing performance, the study validated the results obtained. Table 4 showed the performance validation results based on the MAD, RMSE, and MAPE measurements. Figure 5 showed forecasting NOX emission data for the Road and PP sectors in the 2016-2020 period. This data forecasting was based on the best performance obtained, namely at n = 2 (Road) and n = 3 (IND). Table 4 showed that data smoothing using the DMA was better than the DES methods. This can be seen from the smaller MAD, RMSE, and MAPE values for the DMA than the DES methods. For the road transportation sector, the best performance was in the 2 (n = 2) interval, while for the power plant, it was in the 3 (n = 3) interval. For the MAPE value, NOX emission measurement in the road transportation sector was 2.72%, while the power plant sector was 6.68%. Therefore, the accuracy was in the very good category because the value was <10%. Figure 5 showed that the forecasting graph for the transportation sector in 2016-2020 is relatively constant or the changes are not too significant, with a gain of 17.46 kton / year. Meanwhile, the power plant sector showed a significant increase with a gain of 120.02 kton / year Figure 6 showed the results of smoothing SO2 emission data using DMA and DES techniques for the power plant (PP) and the industrial sectors (IND). The results of smoothing the data from the two methods were then validated using MAD, RMSE, and MAPE. Subsequently, the results were compared to determine the best performance. Table 5 showed the validation results of both methods used Figure 6: Graph of SO2 emission data smoothing for (a). PP sector using DMA (b). PP sector using DES, (c). IND sector using DMA and (d). IND sector using DES  Table 5 showed that the best performing data smoothing for the PP and IND sectors using DMA was at n = 2 and n = 4, while for DES, it lies at weights =0.4 and =0.2. When the results are compared, the DMA method was better than DES, because the values of MAD, RMSE, and MAPE were smaller. For MAPE, the measurement of SO2 emissions in the power plant sector had a value of 8.41%, while the industrial sector was 8.45%. Therefore, the accuracy was in the very good category because the value was <10%.

SO2 emissions
Forecasting data for SO2 emissions in the Road and PP sectors in the 2016 -2020 period can be seen in Figure 7. This forecasting is based on the best performance obtained, namely at n = 2 (PP) and n = 4 (IND). Figure 7 showed that the data forecasting graph for the power plant sector in 2016-2020 has increased significantly with a gain of 798.5 kton / year, while for the industrial sector showed a decline with a gain of -73.54 kton/year.
In 2015, the installed power capacity for Steam Power Plants (PLTU) was around 53.1% of the total power generation in Indonesia [19]. This showed that the need for coal fuel to supply PLTU needs is very large. Furthermore, the increase in the amount of coal is proportional to the increasing demand for electricity. Therefore, SO2 gas pollution in the period of 2016 -2020, appears to have increased significantly. In contrast to the power plant, SO2 emissions in the industrial sector have decreased, because many industries have implemented the clean production principle by reducing the use of coal as their industrial fuel [20].

Conclusions
The use of Double Moving Average (DMA) method in smoothing NOX and SO2 emission data is better than the Double Exponential Smoothing (DES). This is based on the validation results of both methods, where the MAD, RMSE, and MAPE values for the DMA method are smaller than DES. For NOx and SO2 emissions, the resulting MAPE value in the DMA method was in the very good category, because the value was <10%. Furthermore, the best variation of the moving average interval for NOX emissions lies at n = 2 (road transportation sector) and n = 3 (power plant sector), while for SO2 emissions, it lies at n = 2 (power plant sector) and n = 4 (Industry sector).
The forecasting stage for the period of 2016 -2020 generally showed that the total emissions of NOX and SO2 have increased, with the gains for each emission being 533.9 kton / year and 227.8 kton / year. This means that the total emission gain of NOX is two times greater than SO2. Furthermore, in the road transportation sector that produces NOX emissions, changes are not too big, with an emission gain of 17.46 kton / year. Meanwhile, in the power plant sector, it has a quite significant increase, with an emission gain of 120.02 kton / year. SO2 emissions showed that there is an increase in emissions of power plant sector with a gain of 798.5 kton / year, while for the industrial sector, there is a decrease with a gain of -73.54 kton / year. In addition, the reduction in SO2 emissions in the industrial sector showed that there is an industrial policy to reduce the use of coal as its energy source.