Orientation of global warming and its effects, especially in the Arctic, is quite obvious, and the main reason for this is the redundancy of CO₂ in the atmosphere [1] . To study the association between elevated atmospheric CO₂ and global warming, several studies have been conducted with different methods. Most of the approaches are the analysis of temperature data as well as the increase of atmospheric CO₂ extraction process. Due to the growing trend of global temperature and the rapid increase in atmospheric CO₂, the impression that there is a link between them is reinforced [2] . Assuming the establishment of such a link, Salahi et al. (2006) by simulating changes in temperature and precipitation in Tabriz city in terms of doubling of atmospheric carbon dioxide conclude that in this case, the temperature for all seasons increase significantly [3] . Carbon dioxide measurements in the atmosphere are very precise and show the annual growth of the accumulation in the atmosphere of the earth. Despite this certainty, the amount of carbon dioxide produced by the earth cannot be determined yet. Also the amount of gas absorbed by the land and oceans is not determined [4] , although a substantial sum is formed. For example, the most important way of reducing carbon dioxide in the atmosphere after reducing the consumption of fossil fuels is carbon storage in terrestrial ecosystems and land. Soil contains 75% of the land carbon source [5] . Determining the emission or absorption of carbon dioxide in the atmosphere requires a detailed knowledge of at least one of them. By subtracting the amount of gas in atmosphere we can obtain the other. Reaching this trust requires a wide range of global geophysical and economic data. Corinne Le Quéré (2009) prepared an overall rate of all production and absorption of carbon dioxide from 1959 to 2008 to develop and analyze the phenomena that affect them. The growth rate of carbon dioxide in the atmosphere was determined on the basis of direct measurements. The amount of carbon dioxide produced by fossil fuels was detected according to statistics from the consumption of fossil fuels in different countries. The amount of carbon dioxide that comes from land use changes was estimated by statistical deforestation, land use, spatial observations fires, and assumptions about the carbon density in the vegetation and soil. The result of the operation shows that the amounts of the estimated residuals were not compatible with the quantity of gas in the atmosphere. The inability to identify the source and extent of absorption and excretion of carbon dioxide prevents accurate understanding of the carbon cycle [4] . The purpose of this article is to present a method to identify locations where gas adsorption or desorption may have occurred. This method works by implementing spatial analysis in the global statistical information quantities of gas in the atmosphere.

Most studies on the precise knowledge of the factors that influence the uptake and release of carbon dioxide began around the 1990s. These studies can be divided into two groups: 1) Assessment of changes in carbon dioxide in the atmosphere, and understanding the effective factors that cause them. 2) Statistical analysis of climate data and maps to predict the spread and absorption of gases. If the amount of carbon dioxide in the atmosphere is the balancing result of known sources and sinks, according to the known factors that influence the absorption and release of carbon dioxide gas the budget can be shown as indicated by the following equation [6] :

The trend of carbon dioxide in the atmosphere is equivalent to the sum of the amount of gas emitted by the sources and absorbed by the sinks. According to measuring stations on the surface of the Earth, an average of 40% of the total annual CO₂ emissions between 1959 and 2008, remain in the atmosphere. But this number, as shown in Figure 1 [4] [7] , has also high annual variations too.

The amount of CO₂ emission values with the absorbed CO₂ is called residual values. The trend chart is shown in Figure 2 [4] [8] . The range in which the shadow appears, the calculation is approximate.

As the comparison between Figure 1 and Figure 2 shows, the trend of changes in atmospheric carbon dioxide does not justify the fluctuations of its calculated residual amounts. In other words, the growth rate of CO₂ in the atmosphere does not correspond with the calculated emissions and absorptions [4] .

Geostatistical analysis is another way that is used to learn more about the carbon cycle. Using two methods of neural networks and fuzzy network on meteorological data in Tehran, Khazaei et al. (1391) predicted the concentration of carbon monoxide emissions.

All arithmetic operations used in the studies mentioned in the previous section, are often limited to the statistical analysis of (non-spatial) descriptive data, while the spatial resolution with its good analytical equipment that is able to identify areas with certain characteristics in the environment is suitable to apply. The importance of the method presented in this study is its approach in using this kind of analysis.

The data used in this study were carbon dioxide values in the atmosphere between 2000 and 2011 that is measured by 63 sampling stations, through the database offered in GLOBAL VIEW. Information available on this site starts from 1999. Prior to 1988, the number of measurement stations was much less. Studies show that at least 10 stations for each study area are required to obtain good estimates [9] . The samples are offered both daily and monthly. By calculating the average of monthly data, annual values for each station were calculated and used. Using the geographical coordinates of stations and a view of world map, the working platform was prepared. Using polygon analysis of Thiessen, the functional areas of each station is shown in Figure 3.

Using annual data for each station and their geographical positions, a Triangular Irregular Network (TIN) model was produced. Figure 4 illustrates an example of this model for 2009.

Then by using the TIN models, the contour maps were generated from them. By implementing a method that will be explained in the next section, for each station a numerical value proportional to its role in the sink or source of the region obtained. These values are named in this paper by “Stations Volcanic Degree” (SVD). In conclusion, the amounts of these SVDs were compared during 2000-2011.

Figure 9 shows an example of the results of this paper in relation to 2009 data.

In this case, the strongest peak, with 107 SVD value is known at the BSC station in the “Black Sea” area. The strongest valley with SVD value of −41 is known at the BKT position near “Indonesia”. With these assumptions, the probability that in 2009 at the Black Sea region a strong CO₂ emission factor being founded is much more than any other place. In addition, the gas absorption factor in the position of Indonesia station will be more discoverable. By applying the proposed method, to CO₂ values of stations during 2000 to 2011, the 12 maps were produced. In each map the peaks and valleys were illustrated. The peaks and valleys were graded with their SVD values. Table 1 shows 10 stations that have the highest degree of positive and negative values of SVD.

Table 1 show that there is a strong domination in the BSC and BKT stations from all others. The station of the Black Sea (BSC), with its high SVD values, shows the highest degree of CO₂ emission zone, and can be a good option for discovering the CO₂ emission factors. This inference exists while the average values measured at the station is something close to the average of the others, and apparently no particular dominance over others, can be seen. And the Indonesia station (BKT) is recommended as a priority area for detecting CO₂ gas absorption factors. The SVD value of the PAL station has always been negative, and only in 2001 became largely positive. It can also be a sign of a temporary source of CO₂ emissions in 2001 in the region. Figure 10 shows the location of these stations.

In this paper, using the measuring stations annual data of carbon dioxide in the atmosphere and their locations, a TIN model was produced for each year. Then with a model that has been proposed based on the idea of the formation of volcanic mountains, the peaks and valleys of the TIN model were rated based on their similarities with forms of volcanic mountains. These grades were saved as Station Volcanic Degree (SVD). It is assumed that the special formation of volcanoes could be a suitable pattern for extracting those peaks and valleys of the TIN model which the related sources or sinks are located inside their own regions. By considering the computation results on data from 2000 to 2011, possible source locations of gas emission and adsorption were offered. As

mentioned in the previous sections, despite the extraordinary importance of preventing the growth of residual carbon dioxide in the atmosphere, the identification of factors that influence the emission and absorption of the gas is extremely low and rough. The method proposed in this paper is an innovative approach compared to previous studies, and expresses the results quickly and explicitly. It is good to use a method to evaluate the results of this study with cases of observation. In addition, because the residual quantity of atmospheric carbon dioxide is strongly dependent on the temperature, it is good to do this research, again with seasonal data, and the results reviewed.

I would like to deeply thank Dr. Ali Javidaneh who cooperated in the development of research ideas. And also Mr. Nima Ghasemloo and Dr. Hamid reza Ranjbar who helped me frankly in implementing the computer programs for this study.