Abstract
It is essential to study CO2 emissions intensity as the most critical factor affecting temperature increase and climate change in a country like Iran, which ranked seven regarding CO2 emissions intensity. Investigating the convergence of CO2 emissions intensity is essential in recognizing its dynamics in identifying the effectiveness of government environmental policies. In this paper, using the Markov chain and spatial Markov chain methods, the convergence of CO2 emissions intensity from fossil-fuel consumption has been investigated in 28 provinces of Iran from 2002 to 2016. The empirical results showed that convergence clubs and neighbors significantly influenced the transition probability of regions to clubs with high and low CO2 emissions. Therefore, if a province had a neighbor with low (high) CO2 emissions intensity, the transition probability of this province to the club with low (high) CO2 intensity increased. Therefore, in provinces that have neighbors with low (high) CO2 emissions intensity, the transition probability to the club with low (high) CO2 intensity increases.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
CO2 emissions increased from 22 billion tons in 1990 to more than 36 billion tons in 2019 (ourworldindata.org). However, in 2020, with the COVID 19 pandemic, which reduced energy demand by 4%, the global CO2 emissions related to energy fell by 5.8% to about 2 billion tons per year, unprecedented in human history (www.iea.org). However, energy consumption and CO2 emissions are expected to return to normal with the prevention of the COVID 19 pandemic. Figure 1 shows the highest CO2 intensity (CO2-I) countries in the world. South Africa, Kazakhstan, Taiwan, Russia, and Ukraine have the highest CO2-I. Iran ranks seventh in the world regarding CO2-I. A comparison of the trend of CO2-I between Iran and the global average showed that although CO2-I in the world has been decreasing from 1990 to 2018, this trend is increasing in Iran (Fig. 2). This is in the context that Iran has the second-largest gas reserves after Russia and the fourth-largest oil reserves after Venezuela, Saudi Arabia, and Canada (Dudley, 2020).
One of the most critical effects of CO2 emissions is climate change and global warming. This is essential for survival on earth because “Every bit of warming matters” (IPCC, 2018). Figure 3 shows the effect of emissions on heating imbalance, in watts per square meter, caused by all major artificial greenhouse gases. It is known that the highest share and highest growth of greenhouse gases are related to CO2 emissions. As a result, the absorbed heat per watt per square meter had almost doubled from 1980 to 2019. Given the high effects of CO2 emissions on climate change and temperature, it is essential to identify the factors affecting CO2 emissions. Many studies, including Wang et al. (2018), Kang et al. (2019), Chen et al. (2019), and Shahnazi and Shabani (2021), have been done in this regard. One of the factors related to CO2 emissions is CO2 emissions convergence.
This research aims to investigate the convergence of CO2 emissions intensity in Iranian provinces. The convergence of CO2 emissions can be explained theoretically based on three foundations, i.e., environmental Kuznets theory, the orientation of global environmental treaties, and economic growth convergence.
The environmental Kuznets theory is one of the most important models used to explain CO2 emissions. Accordingly, the relationship between economic development and environmental quality is an inverse-U relationship. According to Grossman and Krueger (1991), economic development affects environment in three ways: (i) scale effect in the sense that increasing the scale of the economy leads to more destruction of natural resources and more emissions; (ii) composition effect: with the increasing development, production structures change and move towards cleaner activities with less pollution; (iii) technological effect resulting from increasing research and development costs in the development process which leads to new clean technologies. Therefore, considering these effects, economic growth has negative environmental effects in the short run and leads to increase pollution. However, in the long run, with the growth of clean technology development, the quality of environment improves, and different regions converge in CO2 emissions. The second explanation for the convergence of CO2 emissions is related to international environmental treaties such as the Kyoto Protocol or the Paris Agreement (Haider and Akram 2019). International environmental treaties can be the basis for CO2 emissions convergence to provide unified policies for improving environmental conditions in different countries. Regarding the treaties of universal laws and regulations from 1850 up to now, current content includes 1368 Multilateral environment agreements (MEAs), 2264 bilateral environment agreements (BEAs), and 250 other environmental agreements (https://iea.uoregon.edu). It should be noted that Iran is present in 320 of the agreements described in Table 1. The third basis for explaining the convergence of CO2 emissions is related to the principles of convergence of economic growth. The concept of convergence is stemming from economic growth literature (Barro, 1991). Due to the fact that part of CO2 emissions is related to energy consumption for economic growth, as convergence in economic growth, CO2 emissions associated with economic growth also converge.
This study has two significant contributions: first, a non-parametric method that considers the spatial effects was used to assess CO2 emission convergence. Although parametric methods have been used in different studies, which may provide misleading results as multimodality and distribution dynamics are ignored (Quah, 1993, 1997). On the other hand, if research deals with spatial data with locational information, spatial effects must be considered in data analysis as ignoring the spatial effects may lead to biased results. Second, as far as researchers examined, it is the first study examining CO2-I convergence in Iran. Previous studies investigated CO2 emissions in Iran (e.g., Mirzaei and Bekri 2017, Hajilary et al. 2018, Hosseini et al. 2019), but CO2-I convergence was not investigated in Iran. It is important to examine CO2-I convergence in Iran regarding international environmental goals to fulfill Iran’s commitments in international environmental treaties and agreements such as the Kyoto Protocol and the Paris Agreement.
This study has the following sections: literature review is provided in “Literature review.” “Methodology and data” presents methodology. The empirical results are reported in “Empirical results.” Conclusion and policy implications are drawn in “Conclusion and policy implications.”
Literature review
Many studies explored the convergence of CO2 emissions is shown in Table 2. These studies differ in the areas, duration, method, and type of convergence. Some studies such as Yu et al. (2018) and Boussemart et al. (2015) focused on studying the convergence of CO2 emissions in the interior of a country, and some like Romero-Avila (2008), Yan et al. (2017), and Li et al. (2020) examined the convergence between countries. The single country study focuses mainly on the study of the USA and China. In multi-country studies, some have examined countries located in a region, mainly including EU and OECD countries. And some articles have reviewed a large number of countries in more than 100 countries. The study period ranges from 11 years (Li et al., 2017) to more than 100 years (Presno et al., 2018).
Previous studies differed in terms of the method used. For example, Romero-Ávila (2008) used the bootstrap methods, Van (2005) used the Kernel density estimates, Li et al. (2017) applied the spatial panel data models, Lee et al. (2008) used the panel unit root stationarity analysis, Huang and Meng (2013) employed the spatial panel data, and Morales-Lage et al. (2019) applied the nonlinear dynamic factor model to investigate the convergence of CO2-I. Another difference between the studies was related to the type of convergence. Accordingly, some studies such as Boussemart et al. (2015) showed the existence of non-club convergences, and others such as Kounetas (2018), Morales-Lage et al. (2019), and Tiwari et al. (2021) showed the existence of club convergences. Although most studies conclude that there is convergence in emission intensity, in some studies, such as Yan et al. (2017), convergence is not confirmed.
A review of previous studies indicates the existence of two research gaps that the purpose of this article is to fill two. First, previous studies have mostly been performed with parametric methods or have not seen spatial effects. And second, CO2-I convergence has not been studied in Iran.
Some previous studies have used parametric methods that may provide misleading results as multimodality and distribution dynamics are ignored (Quah, 1993, 1997). On the other hand, Van (2005) and Kounetas (2018) investigated the convergence of CO2-I using a non-parametric method while not considering the spatial or neighbor effects on CO2-I convergence; ignoring the spatial effects may lead to biased results. This paper uses a non-parametric method of spatial effects to fill this research gap. As Table 2 shows, there is no study on CO2-I convergence in Iran. Although Iran has seventh place globally, CO2 emissions intensity and the trend of CO2-I are increasing in Iran, contrary to the global trend from 1990 to 2018 (https://yearbook.enerdata.net); identifying CO2-I convergence is very important for Iran in terms of policies to reduce emissions.
Methodology and data
Theoretical models
The concept of convergence was first introduced in the economic growth literature in 1991 and is now being used in other fields such as carbon dioxide emissions. Three concepts of beta, sigma, and random convergence can be used to check convergence. Sigma convergence examines the decreased difference in carbon dioxide emissions between regions over time. The increased difference in carbon dioxide emissions between regions means divergence. The sigma convergence is investigated using the standard deviation of carbon dioxide between regions or using the kernel density function. Beta convergence is divided into unconditional and conditional convergence. The convergence of CO2 emissions of all regions into a common and stable steady state is called unconditional beta convergence. Their convergence to the steady-state conditions in proportion to the initial level of factors in each region is called conditional beta convergence. Based on the conditional beta convergence, each region converges to its own specific steady state, and regions farther from their equilibrium point converge more rapidly to their steady state. This type of convergence led to forming a new concept called club convergence or multiple equilibrium models. In this case, regions with similar carbon dioxide emissions (e.g., high carbon dioxide emissions and low carbon dioxide emissions) converge to a steady state. Thus, in the long run, not only the difference between the regions does not disappear, but the carbon dioxide emissions of the regions converge towards similar peaks, and carbon dioxide emissions distribution of the regions moves towards two or more peaks. In this case, homogeneity between regions or countries increases, and differences between groups increase (Le Gallo, 2004). An examination of economic literature showed that the concept of club convergence was first proposed by Baumol (1986) and then was completed by Durlauf and Johnson (1992, 1995). Next, Quah (1997) developed it in other ways. Kova argued that the econometric methods, a common method for examining beta convergence, often have misleading results, produce Galton’s problem, and do not describe distribution behavior and movements. Hence, he introduces the Markov chain as a powerful tool for studying and identifying convergence clubs (Quah, 1993). Stochastic convergence indicates that shocks to carbon dioxide emissions are temporary (Fallahi, 2017). A randomness unit root test is used to investigate random convergence.
In this paper, due to limited information, it was not possible to investigate Stochastic convergence. Moreover, since the provinces of Iran are not homogeneous and similar in terms of economic and geographical factors, conditional beta convergence (club convergence) should be considered. According to Quah (1993), the use of measurement methods to investigate the conditional beta convergence (club convergence) leads to misleading results. Hence, in this paper, the Markov chain method was used to investigate club convergence. In addition, the kernel density function has also been used to investigate Sigma convergence.
Empirical models
Kernel density
The kernel density estimator is a weighted moving average of the sample frequency distribution. The function K is the kernel if it is as follows:
In most cases, K is a symmetric probability density function. The kernel density estimator at point x with K is as follows:
where n is the number of observations, xi is the observation of i, and h is the smoothing parameter or bandwidth.
Markov chain
A random process with a set of finite values xr, t + 1, xr, t, xr, t − 1, …, xr, 0 is called a Markov chain, if this process is in state i at time t (xr, t = i) and the probability Pij is in state j at time t + 1(xr, t + 1 = j). So, for all values t ≥ 0 and all states,j, i−1, i−2, …, i0, we can write:
P ij is the probability of transition of the process from state i now to state j in the next time (Bickenbach and Bode 2003). If the Markov process is time-independent, the Markov chain is completely defined by the Markov transfer matrix (M). M is called the transition probability matrix; each of its values in row i and column j is Pij, which is the probability of moving from state i to j. Assuming the number of states is 2, the Markov transition probability matrix is a 2×2 matrix as follow:
The number of states in the Markov chain must be chosen so that the initial distribution of data in the different states is approximately the same (Quah 1993). The initial distribution of the data is explained by the initial vector (Monfort 2008). If the Markov chain is irreducible and periodical, the process reaches a steady-state or stationary-state vector. The steady-state vector is called the ergodic distribution vector and calculated using an eigenvector associated with the eigenvalues of 1 in the transpose of transition probability matrix (Ibe 2009). The ergodic vector is used for considering convergence in the literature (Liddle 2012, Li et al. 2019). If the data concentration is in one state in the ergodic vector, there is absolute convergence. If the data is divided between different states with equal frequency, there is a divergence. Convergence clubs exist when the data are concentrated in some states (Le Gallo, 2004). According to transition probability matrix, we can calculate stability and half-life index. The stability index is as follows:
where N is the number of states, and Tr(M) is the trace of the matrix M (Monfort 2008). The half-life index shows the speed of convergence toward the ergodic vector and calculates the years to cover the half distance to the ergodic vector. This index is as follows:
where α2 is the largest eigenvalue after 1 of transpose M, T is year interval (Shorrocks 1978). In this study, the year interval is 1. Matrix M presents the transition probability matrix but does not take into account spatial autocorrelation between observations of regions.
The Moran’s I statistic (1950) is the common method to investigate spatial autocorrelation of regions as follow:
where Ci is the variable of region i,\(\overline{C}\) presents the average of variable, n is the number of regions, and Wij is ith and jth components of the spatial weights matrix.
For considering the effect of neighbors on the transition probability matrix, Rey (2001) proposed spatial Markov transition probability. In this matrix, the transition probabilities depended on the state of neighbors (spatial lag). If the number of states is 2, the spatial Markov matrix is a 2×2×2 matrix (Rey 2001).
Table 3 presents the spatial Markov transition probability with two states. Spatial lags are determined based on the average CO2 in the neighboring provinces in the starting year (e.g., 2002). For example, provinces are in state one according to spatial lag when the average CO2 of neighbors is in a state of 1 in 2002. In Table 3, pij/k is the probability shifting from state i to j when the neighboring provinces are in k state. The spatial effects can be understood by comparing the elements of matrix M with the corresponding elements in Table 2. If the transition probabilities are as follows:
There is no spatial effect. Spatial effects exist when, for example, P12/1 < P12, meaning that if the neighbors of a province are in state 1, then the province moves from state 1 to state 2 with a lower probability.
Data
Data of CO2-I for 28 Iranian provinces in the period of 2002 to 2016 was used in this study. To calculate the CO2-I, first, energy consumption was extracted from different energy carriers. Second, the CO2 emissions was calculated based on the carbon dioxide coefficient for each energy carrier documented in the Iranian energy balance sheet. Third, CO2 emissions for each province is divided into the real GDP to calculate CO2-I. The data of GDP was collected from the statistical yearbook of provinces (www.amar.org.ir). The real GDP was calculated using the provincial consumer price index.
In the Markov chain method, relative variables are used and the relative CO2-I is applied. To do this, CO2-I in each province was divided into the average CO2-I in Iran. Table 4 shows the descriptive statistics of CO2-I, GDP, and the relative CO2-I in Iranian provinces.
Empirical results
The kernel density in Fig. 4 shows how CO2-I in 2002 evolves into 2016. The probability density function was not symmetric and was left-skewed and had a positive kurtosis. The positive kurtosis indicates that the distribution is peaked and has thick tails relative to normal distribution. As shown in Fig. 4, the distribution of the relative CO2-I of Iranian provinces in 2002 was not normal and was in 105% of the country’s mean. Therefore, evidence suggested that the CO2-I of most provinces in Iran was higher than the country’s average. In addition, kurtosis of distribution showed more data at this point than the normal state. Figure 4 shows that the concentration of provincial CO2-I was 120% of the country’s average in 2016. Moreover, the negative skewness of distribution has increased, which shows that the CO2-I of some provinces has been decreased compared to 2002.
A comparison of the four density distributions in 2002, 2007, 2012, and 2016 showed that the distribution kurtosis increased in 2016, and the skewness shifted more to the left than in 2002, 2007, and 2012. Thus, over the 15 year, two peaks became apparent, one small peak was formed at around 30% of the country’s average, and a major peak was formed above the average. So, there are individual convergence clubs in CO2 emissions. The reason for this can be inequality in GDP per capita (Padilla and Duro 2013), and different energy mix and efficiency (Barassi et al. 2008; Wang et al. 2005), different industrial structures (Talukdar and Meisner 2001), and demographic variation (Knapp and Mookerjee 1996) among the provinces of Iran.
Although the distribution density function showed the shape of the distribution and formation of the convergence clubs, it cannot explain the intra-distribution movements and the evolution of distribution. Therefore, to investigate the intra-distribution dynamics, steady state, club convergence, and convergence speed of CO2-I, the first-order static Markov chain model, can be used.
In this study, the number of observations was 420, classified into four states. State 1 includes observations whose CO2-I was 0 to 75% of Iran’s average CO2-I, which is called low CO2-I. In state 2, observations whose relative CO2-I were between 75 and 100% of average are called medium-low CO2-I. The relative CO2-I is between 100 and 120% of Iran’s CO2-I average in state 3, a medium CO2-I state. Finally, observations with CO2-I above 120% of the average are in state 4. This state is called high CO2-I. The initial distribution of observations between different states using the initial distribution vector showed that 24% of the observations were in state 1, 25% of the data were in state 2, 25% of the observations were in state 3, and 26% were in state 4 (Table 5). In this study, four states are considered, so the transition probability matrix was a 4×4. The transition probabilities matrix is presented in Table 5, the elements of which indicate the probability of transfer from one state to another. As the results in Table 5 show the element on the main diagonal of the transition matrix were high and, on average, 80%. Therefore, the stability index was 80%. Moreover, the probability of movements outside the main diagonal was 7% on average. Therefore, according to the stability index and mobility outside the main diagonal, as well as the high probabilities staying in the states 1 (95%) and 4 (80%), it can be concluded that the provinces with high and low relative CO2-I do not see changes in their relative position over time.
Table 5 showed that 87 observations were in state 1 during 15 years which remained in state 1 with a 95% probability the next year and will shift up to states 2 and 3 with a probability of 4% and 1%, respectively. Therefore, provinces with low CO2-I are very unlikely to leave the low CO2 emissions club. According to the data (Appendix Table 10 in appendix), Ilam, Tehran, Khuzestan, Kohkiluyeh, and Boyer-Ahmad provinces have been in the state 1 in 15 years. Golestan has been in state 1 for 8 years and then has moved to state 2; Lorestan has been in state 1 for 7 years. Moreover, Ardabil province has been in state 1 for 1 year and has shifted to state 2. Furthermore, 95 observations were in state 2, of which 78% remained in state 2 in the next year, and the highest probability of movement was 17% into state 3. Chaharmahal Bakhtiari province for 15 years; Ardabil and Zanjan for 12 years; Azerbaijan Sharghei, Kerman, and Lorestan for 8 years; Semnan, Fars, and Golestan for 6 years; Khorasan for 5 years; Kurdistan for 3 years; Gilan for 2 years; Bushehr, Mazandaran, and Markazi for 1 year were in state 2. Most observations (120 observations) were in state 3, which had a CO2-I of 100 to 120% of Iran’s average. According to Table 5, the observations in state 3 will remain in this state with a probability of 66% in the next year. The probability of shifting down to state 2 is 12%, and the probability of shifting up to state 4 is 22%. Azerbaijan Gharbei, Khorasan, Gilan, and Mazandaran have been in state 2 for more than 10 years.
As shown in Table 5, 118 observations were in state 4, with a probability of 80% remaining in state 4 in the next year and with a probability of 18% and 2% moving down to states 2 and 3 in the next year, respectively. Isfahan, Sistan and Baluchistan, and Qazvin provinces have been in state 4 for 15 years; Hormozgan for 14 years; Hamedan for 10 years, Markazi for 9 years; Kurdistan for 8 years; Yazd, Kermanshah, and Qom for 6 years; Azerbaijan Gharbei for 4 years; Bushehr, Semnan, and Gilan provinces for 2 years; and Mazandaran, Kerman, Zanjan, and Azerbaijan Sharghei provinces for 1 year. The ergodic distribution vector is presented in Table 5. The results of this vector were used to investigate the convergence. The results of ergodic vectors showed that in the long run, 87% of the observations are in states 1 and 2. The number of regions in state 1 or low CO2-I club will increase from 24% in the initial vector to 27% in the ergodic vector. Moreover, the observations in state 2, i.e., in the medium-low CO2-I club, will increase from 25 to 59%. According to the results of the ergodic vector, the observations in state 3, i.e., the medium CO2-I club, will decrease from 25 to 12%, and the number of observations in state 4 will decrease from 26 to 2%. However, the results indicate that there is club convergence. This finding was contradicted Boussemart et al. (2015) and Huang and Meng (2013) that support convergence in CO2 emissions across provinces in a country, and consistent with Rios and Gianmoena (2018), Kounetas (2018), Apergis and Payne (2017), and Tiwari et al. (2021) that predict the club convergence of CO2 emissions. A possible explanation can be attributed to the fact that the provinces of Iran have different climatic conditions and structures of economies. Therefore, provinces with similar climates, production composition, fuel consumption, and energy efficiency will converge towards the same provinces.
The half-life index was 14 years, indicating that 14 years are needed to cover the half distance between the initial vector to the ergodic vector. The results of half-life and stability indices indicated that the move towards long-run equilibrium was very slow and stability of the process was high. The results of Table 5 examined the evolution and dynamics of CO2-I distribution without considering the neighborhood effects.
However, the data used in this study were spatial data containing the locational information, so the spatial effects should be considered in the data analysis. Figure 5 shows the continuous state in which each province is located from 2002 to 2016. As Fig. 5 shows, most of the provinces that had a similar continuous state were located next to each other. A spatial correlation test must confirm this visual impression of spatial correlation.
The common test for the spatial autocorrelation was Moran’s I statistic. According to the results of Table 6, the Moran value was positive and was statistically significant in all years except 2009, showing that there was a positive spatial autocorrelation between CO2-I in Iranian provinces. As a result, provinces with high (low) CO2-I were adjacent to each other. This finding is consistence with Shahnazi and Shabani 2019, 2021, Liu et al. 2019, You and Lv 2018.
Table 7 shows two transition probability matrices. The first matrix was constructed using the neighbor-relative CO2-I distribution, in which CO2-I of each province is normalized to the average CO2-I of the neighboring provinces. The second matrix included the probability of transition between Iran-relative CO2-I distribution and the neighbor-relative CO2-I distribution. This matrix showed probability transition between two different distributions at a given time. If there are no spatial effects, the first matrix in Table 7 is similar to the matrix in Table 5, and the second matrix in Table 7 is a unity matrix. The results of Table 6 showed that the first matrix is not similar to the matrix in Table 5, and the second matrix is not a unity matrix.
N, number of observations, upper endpoint is a ratio of the province to the national average CO2-I
Table 8 shows the spatial transition probability matrix, which is conditioned on the spatial lag in 2002. The spatial lag indicates the state of the average CO2-I of neighbors. It seems that the spatial lag of a province affects the transition probability of a province between different states. For example, provinces with low CO2-I are negatively affected by the provinces with high CO2-I. In fact, by increasing the CO2-I of neighbors, the transition probability of moving up to the high states increase. According to Table 5, provinces with low CO2-I (state 1) with an average probability of 5% of this state shift up to states 2 and 3. However, according to Table 8, if the neighboring provinces with low CO2-I (state 1 in spatial lag) surrounded these provinces, the probability of transition upwards is only 4%. But, if the neighboring provinces had a medium-low CO2-I (state 2 in spatial lag), the probability increases to 6%, and if the neighboring provinces were in state 3, the probability increases to 13%.
If the province is in state 2, as shown in Table 5, the probability of moving to states 3 and 4 is 18%. But if this province has neighbors with low CO2-I, then the probability of moving upwards will be 0%. Moreover, if the neighbors are in low-medium CO2-I, the probability of moving upwards increases to 12%, and if it has neighbors with medium CO2-I, the probability of moving upwards is 24%, and if the neighbors are in high CO2-I, then, the probability of moving upwards will be 100%. Furthermore, provinces in state 3 have been affected by the neighbors. As seen in Table 5, if a province is in state 3, the probability of transition to state 4 is 22%. But, according to the results of Table 8, this province will not have neighbors with low CO2-I. If this province is surrounded by the provinces that are in state 2, the probability of moving upwards is 17%, and if it has neighbors with the medium CO2-I (state 3), the probability of moving up is 20%, and in the presence of neighbors that they are in state 4, the probability of moving upwards is 38%. The results of Table-8 indicated that the provinces with high CO2-I had not been affected by their neighbors. The results of ergodic distribution vectors of spatial lag are presented in Table 8. The results showed that if neighbors are in low or medium-low CO2-I states, the frequencies of observations in state 1 were 53% and 46%, respectively. When neighbors were in state 3, the concentration of observations in state 3 was 42%, and when neighbors were in state 4, the data concentration in state 4 was 100%. As a result, the results of the ergodic vectors indicated that their neighbors influenced the provinces.
Finally, we can summarize the information of the Markov spatial chain by considering the relationship between the probability of a specific transition (down, none, or up) and the CO2-I of neighbors (spatial lag). Table 9 showed that neighbors had an important effect on the probability of moving down or up. For example, when a province was surrounded by provinces with low CO2-I, the probability of moving downward was higher than moving upward. Moreover, if the provinces had neighbors with the high CO2-I, then the probability of shifting upwards was higher than downwards. From the perspective of interregional convergence, all these results highlighted that the geographical environment affected interregional convergence. In this way, the neighbors affected the location of each province in the convergence clubs. This finding can be considered a result of spatial or social interactions in the provinces. These interaction processes can be considered to result from economic linkages, knowledge spillovers, policy copying among provinces, and trade (Huang and Meng 2013).
This result contradicts Huang and Meng 2013, who state that spatial autocorrelation will lead to convergence. Also, our finding is consistent with Li et al. (2017) and Rios, V., and Gianmoena, L. (2018). In these studies, spatial autocorrelation has caused divergence and convergence clubs.
Conclusion and policy implications
Iran is one of the countries with the highest CO2-I globally, making it essential to study CO2-I convergence to encourage the government’s CO2 policies to decrease the external effects of CO2 emissions. On the other hand, assessing CO2-I convergence is essential as the CO2-I gap between provinces (or regions) can assist policymakers design and implementing effective policies to decrease CO2-I and the loss of the environment. As a result, the Markov chain and spatial Markov chain methods examined the CO2-I convergence in 28 Iranian provinces from 2002 to 2016. The results indicate convergence clubs in CO2 emissions in Iranian provinces. Therefore, the results emphasized that the geographical environment affects interregional convergence. In this way, each province’s location in the convergence clubs is influenced by its neighbors.
Figure 6 presented a conditional map of CO2-I. These figures demonstrated the CO2-I of provinces in 2016 compared to 2002. Accordingly, the dynamics of CO2-I are different in different provinces. While the provinces of Ardabil, Zanjan, West Azerbaijan, Kurdistan, and Yazd have improved, in Fars, Chaharmahal and Bakhtiari, Khorasan, Semnan, and Qom provinces, the situation has worsened. The intensity of CO2 emissions in Iran is increasing. Also, the results showed that there was a fundamental difference in the overall CO2-I and dynamics of CO2-I in Iranian provinces.
These findings have policy implications for emissions reduction targets, as recognizing differences in convergence patterns in Iranian provinces helps design and implement emission reduction policies.
The amount and intensity of carbon dioxide emissions in Iran are high for two main reasons: macro and micro. In both categories, based on the results obtained in this paper, it is important to pay attention to regional spillover effects. At the macro level, the high CO2-I in Iran is related to national and regional reasons. Factors such as the low price of fossil fuels are important at the national level. A different and unequal technology pattern between regions can lead to the high CO2-I and CO2-I divergence at the regional level. So, the spillover effects of CO2 emissions should be considered by the policymakers, and paragraph 3 of Article 6 of the Sixth Development Plan of Iran (2017–2021), in which the distribution of funds obtained from pollution tax is limited to population proportion should be changed by considering the spillover effects of CO2 emissions for neighboring cities or provinces. At the micro-level, the main factors of high energy intensity and subsequent increase in CO2 emissions in Iran are the following: low and out-of-date power plants (Iranian parliament research center, 2015a). About 1.5 times higher electricity transmission and distribution losses in Iran compared to the global average (data.worldbank.org), vehicles with old technologies and high energy consumption, more in low-income areas, weakness of rail transportation system, disadvantages of intercity and intercity road networks and heavy traffic loads, a high share of old technology industries and low share of high-tech industries in GDP (Iranian parliament research center, 2015b; Shahnazi 2021), and imbalance in the type of industry in different provinces of Iran.
In order to promote regional convergence, the adoption of new technologies and greater use of renewable energy, making older equipment more efficient, and changing consumer behavior are effective to achieve long-term convergence. Therefore, the Iranian government must provide policy tools to introduce and disseminate technologies throughout the provinces. At the same time, the provinces, especially those that use inefficient technologies and have the high CO2 emission density (including Hormozgan, Qazvin, Hamadan, Isfahan, Kermanshah, Sistan and Baluchestan, Qom, Gilan, and South Khorasan), must be active in promoting energy-intensive and low-emission technologies. Of course, adopting new technologies must consider local characteristics and regional compatibility. Providing incentives for technological development and innovation, efficient energy production, and distribution is also effective in convergence. Hence, as a complementary policy, supportive fiscal policy can be effective. The government can reduce taxes on purchasing new, energy-intensive, and low-emission technologies or provide financial subsidies to adopt efficient technologies.
Data availability
The datasets generated and/or analyzed during the current study are available in the https://pep.moe.gov.ir/.
References
Ahmed M, Khan AM, Bibi S, Zakaria M (2017) Convergence of per capita CO2 emissions across the globe: insights via wavelet analysis. Renew Sustain Energy Rev 75:86–97
Apergis N, Payne JE (2017) Per capita carbon dioxide emissions across US states by sector and fossil fuel source: evidence from club convergence tests. Energy Econ 63:365–372
Apergis N, Payne JE (2020) NAFTA and the convergence of CO2 emissions intensity and its determinants. Intl Econ 161:1–9
Barro RJ (1991) Economic growth in a cross section of countries. Quart J Econ 106(2):407–443
Barassi MR, Cole MA, Elliott RJ (2008) Stochastic divergence or convergence of per capita carbon dioxide emissions: re-examining the evidence. Environ Resour Econ 40(1):121–137
Baumol WJ (1986) Productivity growth, convergence and welfare: what the long-run data show. Am Econ Rev 76(5):1702–1785
Bickenbach F, Bode E (2003) Evaluating the Markov property in studies of economic convergence. Intl Reg Sci Rev 26(3):363–392
Boussemart JP, Leleu H, Shen Z (2015) Environmental growth convergence among Chinese regions. China Econ Rev 34:1–18
Chen S, Jin H, Lu Y (2019) Impact of urbanization on CO2 emissions and energy consumption structure: a panel data analysis for Chinese prefecture-level cities. Struct Change Econ Dynam 49:107–119
Churchill SA, Inekwe J, Ivanovski K (2020) Stochastic convergence in per capita CO2 emissions: evidence from emerging economies, 1921–2014. Energy Econ 86:104659
Dudley B (2020) BP statistical review of world energy. BP Statistical Review, London, UK
Durlauf SN, Johnson PA (1992) Local versus global convergence across national economies (No. w3996). National Bureau of Economic Research
Durlauf SN, Johnson PA (1995) Multiple regimes and cross-country growth behaviour. J Appl Econ 10(4):365–384
Fallahi F (2017) Stochastic convergence in per capita energy use in world. Energy Economics 65:228–239
Friedman M (1992) Do old fallacies ever die? Journal of EconomicLiterature 30:2032–2129
Grossman GM, Krueger AB (1991) Environmental impacts of a North American free trade agreement (No. w3914). National Bureau of economic research
Haider S, Akram V (2019) Club convergence of per capita carbon emission: global insight from disaggregated level data. Environ Sci Pollut Res 26(11):11074–11086
Hajilary N, Shahi A, Rezakazemi M (2018) Evaluation of socio-economic factors on CO2 emissions in Iran: factorial design and multivariable methods. J Clean Prod 189:108–115
Herrerias MJ (2012) World energy intensity convergence revisited: a weighted distribution dynamics approach. Energy Policy 49:383–399
Herrerias MJ (2013) The environmental convergence hypothesis: carbon dioxide emissions according to the source of energy. Energy Policy 61:1140–1150
Hosseini SM, Saifoddin A, Shirmohammadi R, Aslani A (2019) Forecasting of CO2 emissions in Iran based on time series and regression analysis. Energy Reports 5:619–631
Huang B, Meng L (2013) Convergence of per capita carbon dioxide emissions in urban China: a spatio-temporal perspective. Appl Geog 40:21–29
Ibe, O. (2009). Markov processes for stochastic modeling (2nd ed). Elsevier Academic Press.
IPCC (2018). Global warming of 1.5°C. An IPCC Special Report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty [Masson-Delmotte, V., P. Zhai, H.-O. Pörtner, D. Roberts, J. Skea, P.R. Shukla, A. Pirani, W. Moufouma-Okia, C. Péan, R. Pidcock, S. Connors, J.B.R. Matthews, Y. Chen, X. Zhou, M.I. Gomis, E. Lonnoy, T. Maycock, M. Tignor, and T. Waterfield (eds.)]. In Press.
Jobert T, Karanfil F, Tykhonenko A (2010) Convergence of per capita carbon dioxide emissions in the EU: legend or reality? Energy Econ 32(6):1364–1373
Kang SH, Islam F, Tiwari AK (2019) The dynamic relationships among CO2 emissions, renewable and non-renewable energy sources, and economic growth in India: evidence from time-varying Bayesian VAR model. Struct Change Econ Dynam 50:90–101
Kang YQ, Zhao T, Yang YY (2016) Environmental Kuznets curve for CO2 emissions in China: a spatial panel data approach. Ecol Indicators 63:231–239
Knapp T, Mookerjee R (1996) Population growth and global CO2 emissions: a secular perspective. Energy Policy 24(1):31–37
Kounetas KE (2018) Energy consumption and CO2 emissions convergence in European Union member countries. A tonneau des Danaides? Energy Econ 69:111–127
Lee CC, Chang CP, Chen PF (2008) Do CO2 emission levels converge among 21 OECD countries? New evidence from unit root structural break tests. Appl Econ Lett 15(7):551–556
Lee CC, Chang CP (2008) New evidence on the convergence of per capita carbon dioxide emissions from panel seemingly unrelated regressions augmented Dickey–Fuller tests. Energy 33(9):1468–1475
Le Gallo J (2004) Space-time analysis of GDP disparities among European regions: a Markov chains approach. Intl Reg Sci Rev 27(2):138–163
Li J, Huang X, Yang H, Chuai X, Wu C (2017) Convergence of carbon intensity in the Yangtze River Delta, China. Habitat Intl 60:58–68
Li W, Zhao T, Wang Y, Zheng X, Yang J (2019) How does foreign direct investment influence energy intensity convergence in China? Evidence from prefecture-level data. J Clean Prod 219:57–65
Li C, Zuo J, Wang Z, Zhang X (2020) Production-and consumption-based convergence analyses of global CO2 emissions. J Clean Prod 264:121723
Liddle B (2012) OECD energy intensity. Energy Efficiency 5(4):583–597
Liu H, Nie J, Cai B, Cao L, Wu P, Pang L, Wang X (2019) CO2 emissions patterns of 26 cities in the Yangtze River Delta in 2015: Evidence and implications. Environ Pollut 252:1678–1686
Mirzaei M, Bekri M (2017) Energy consumption and CO2 emissions in Iran, 2025. Environ Res 154:345–351
Monfort P (2008) Convergence of EU regions: Measures and evolution. European Commission, Regional Policy, Brussels
Morales-Lage R, Bengochea-Morancho A, Camarero M, Martínez-Zarzoso I (2019) Club convergence of sectoral CO2 emissions in the European Union. Energy Policy 135:111019
Moran PA (1950) Notes on continuous stochastic phenomena. Biometrika 37(1/2):17–23
Padilla E, Duro JA (2013) Explanatory factors of CO2 per capita emission inequality in the European Union. Energy Policy 62:1320–1328
Panopoulou E, Pantelidis T (2009) Club convergence in carbon dioxide emissions. Environ Resour Econ 44(1):47–70
Presno MJ, Landajo M, González PF (2018) Stochastic convergence in per capita CO2 emissions. An approach from nonlinear stationarity analysis. Energy Econ 70:563–581
Quah D (1993) Empirical cross-section dynamics in economic growth. Eur Econ Rev 37:426–434
Quah DT (1997) Empirics for growth and distribution: stratification, polarization, and convergence clubs. J Econ Growth 2(1):27–59
Rey SJ (2001) Spatial empirics for economic growth and convergence. Geog Anal 33(3):195–214
Rios V, Gianmoena L (2018) Convergence in CO2 emissions: a spatial economic analysis with cross-country interactions. Energy Econ 75:222–238
Romero-Ávila D (2008) Convergence in carbon dioxide emissions among industrialised countries revisited. Energy Econ 30(5):2265–2282
Shabani ZD, Shahnazi R (2021) Energy intensity convergence in Iranian provinces: evidence from energy carriers’ consumption intensity. Environ Sci Pollut Res, 1-20
Shahnazi R, Shabani ZD (2019) The effects of spatial spillover information and communications technology on carbon dioxide emissions in Iran. Environmental Science and Pollution Research 26(23):24198–24212
Shahnazi R, Shabani ZD (2021) The effects of renewable energy, spatial spillover of CO2 emissions and economic freedom on CO2 emissions in the EU. Renewable Energy 169:293–307
Shahnazi R (2021) Iranian energy and environmental policy: investigation of governance frameworks. Global Encyclopedia of Public Administration, Public Policy, and Governance. https://link.springer.com/referenceworkentry/10.1007/978-3-319-31816-5_4223-1
Shi X, Yu J, Cheong TS (2020) Convergence and distribution dynamics of energy consumption among China’s households. Energy Policy 142:111496
Shorrocks, A. F. (1978). The measurement of mobility. Econometrica: Journal of the Econometric Society, 1013-1024.
Solarin SA (2014) Convergence of CO2 emission levels: evidence from African countries. J Econ Res 19(1):65–92
Solarin SA (2019) Convergence in CO2 emissions, carbon footprint and ecological footprint: evidence from OECD countries. Environ Sci Pollut Res 26(6):6167–6181
Solarin SA, Tiwari AK, Bello MO (2019) A multi-country convergence analysis of ecological footprint and its components. Sustain Cities Soc 46:101422
Sun J, Su CW, Shao GL (2016) Is carbon dioxide emission convergence in the ten largest economies? Intl J Green Energy 13(5):454–461
Talukdar D, Meisner CM (2001) Does the private sector help or hurt the environment? Evidence from carbon dioxide pollution in developing countries. World Dev 29(5):827–840
Tiwari AK, Nasir MA, Shahbaz M, Raheem ID (2021) Convergence and club convergence of CO2 emissions at state levels: a nonlinear analysis of the USA. J Clean Prod 288:125093
Ulucak R, Apergis N (2018) Does convergence really matter for the environment? An application based on club convergence and on the ecological footprint concept for the EU countries. Environ Sci Policy 80:21–27
Van PN (2005) Distribution dynamics of CO 2 emissions. Environ Resour Econ 32(4):495–508
Wang C, Chen J, Zou J (2005) Decomposition of energy-related CO2 emission in China: 1957–2000. Energy 30(1):73–83
Wang S, Wang J, Zhou Y (2018) Estimating the effects of socioeconomic structure on CO2 emissions in China using an econometric analysis framework. Struct Change Econ Dynam 47:18–27
Yan Z, Du K, Yang Z, Deng M (2017) Convergence or divergence? Understanding the global development trend of low-carbon technologies. Energy Policy 109:499–509
You W, Lv Z (2018) Spillover effects of economic globalization on CO2 emissions: a spatial panel approach. Energy Econ 73:248–257
Yu S, Hu X, Fan JL, Cheng J (2018) Convergence of carbon emissions intensity across Chinese industrial sectors. J Clean Prod 194:179–192
Author information
Authors and Affiliations
Contributions
Conceptualization (Rouhollah Shahnazi), data curation (Zahra Dehghan Shabani), formal analysis (Zahra Dehghan Shabani and Rouhollah Shahnazi), investigation (Rouhollah Shahnazi), methodology (Zahra Dehghan Shabani), project administration (Zahra Dehghan Shabani), resources (Muhammad Shahbaz, Rouhollah Shahnazi, Xuan Vinh Vo, and Zahra Dehghan Shabani), software (Zahra Dehghan Shabani), supervision (Zahra Dehghan Shabani), validation (Muhammad Shahbaz, Rouhollah Shahnazi, Xuan Vinh Vo, and Zahra Dehghan Shabani), visualization (Rouhollah Shahnazi and Zahra Dehghan Shabani), writing—original draft (Muhammad Shahbaz, Rouhollah Shahnazi, Xuan Vinh Vo, and Zahra Dehghan Shabani), writing—review and editing (Muhammad Shahbaz, Rouhollah Shahnazi, Xuan Vinh Vo, and Zahra Dehghan Shabani).
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Responsible Editor: Ilhan Ozturk
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Shahbaz, M., Shabani, Z.D., Shahnazi, R. et al. The spatial distribution dynamic and convergence of CO2 emissions in Iran’s provinces. Environ Sci Pollut Res 29, 69573–69587 (2022). https://doi.org/10.1007/s11356-022-20552-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11356-022-20552-z