Analysis of Spatial Effects in the Relationship between CO 2 Emissions and Renewable Energy Consumption in the Context of Economic Growth

: The paper presents a spatial approach to the analysis of the relationship between air pollution, economic growth, and renewable energy consumption. The economic growth of every country is based on the energy consumption that leads to an increase in national productivity. Using renewable energy is very important for the environmental protection and security of the earth’s resources. Promoting environmentally friendly operations increases awareness of sustainable development, which is currently a major concern of state governments. In this study, we explored the inﬂuence of economic growth and the share of renewable energy out of total energy consumption on CO 2 emissions. The study was based on the classical environmental Kuznets curve (EKC) and enriched with the spatial dependencies. In particular, we determined the spatial spillovers in the form of the indirect effects of changes in renewable energy consumption of a speciﬁc country on the CO 2 emissions of neighboring countries. A neighborhood in this study was deﬁned by ecological development similarity. The neighborhood matrix was constructed based on the values of the ecological footprint measure. We used the spatio-temporal Durbin model, with which the indirect effects were determined in relation to the spatially lagged renewable energy consumption. The results of our study also show the strength of the effects caused by imitating actions from the states with high levels of environmental protection. The study was conducted using data for 75 selected countries from the period of 2013–2019. Cumulative spatial and spatio-temporal effects allowed us to determine (1) the countries with the greatest impact on others and (2) the countries that follow the leading ones.


Introduction
Caring for the natural environment should be an integral part of the economic development policy of each country. Unfortunately, state authorities have devoted too little attention to this issue so far, and the degradation of the environment caused by the overexploitation of natural resources and an excessive desire to become wealthy has been extreme. High levels of consumption of non-renewable energy sources and environmental pollution cause an increase in greenhouse gas emissions, mainly carbon dioxide (CO 2 ). Increasingly more emissions have a negative effect on the composition of the atmosphere and global climate [1].
To protect nature, the concept of sustainable development was created, the goal of which is economic development with consideration for the well-being of the present and future generations [2]. The most popular definition of sustainable development is the one formulated by the Brundtland Commission, which describes it as meeting the needs of the present generation without limiting the possibilities of meeting them for future generations.
In particular, sustainable development addresses the problem of reducing the consumption of limited resources of the Earth as well as reducing environmental pollution [3].
The relationship between economic growth and the amount of environmental degradation is usually described by the environmental Kuznets curve (EKC) [4][5][6][7]. In the basic version, the curve expresses the dependence between these processes in the form of an inverted "U" shape, that is, an increase in the level of income of states leads to ever greater environmental degradation, and then, when wealth reaches a certain level, the relationship is reversed. In addition to economic development, renewable energy consumption also has an impact on the natural environment. Increases in the levels of renewable energy consumption, as well as its share of total energy consumption, promote environmental protection [8]. The influence of other factors on the state of the natural environment has also been considered in the literature, for example, the level of trade openness, fossil fuel energy consumption, and the degree of urbanization or population density [9][10][11][12][13][14].
In many countries, an increase in the share of energy from renewable sources out of the total energy consumption has been observed. Moreover, the actions of some countries in this direction have influenced changes to the structure of energy consumption in others. The improvement of environmental conditions resulting from the increase in using renewable energy sources causes an imitation effect.
The aim of this study was to explore the influence of economic growth and the share of renewable energy out of the total energy consumption on the CO 2 emissions for 75 selected countries of the world in the period of 2013-2019. Our concern, in particular, was the impact of changes in renewable energy consumption in a specific country on the air pollution in neighboring countries (the so-called spatial spillovers). In the investigation, we used the spatio-temporal Durbin model (STDM) as a re-specification of the equation based on the concept of the classical environmental Kuznets curve (EKC) and determined the indirect effects in relation to spatially lagged renewable energy consumption. A neighborhood in this study was defined by ecological development similarity. The neighborhood matrix was constructed based on the values of the ecological footprint measure.
The results of our study also show the strength of the effect caused by imitating actions from states with high levels of environmental protection. In particular, the study allowed us to determine the countries with the greatest impact on others as well as the countries that follow the leading ones.
In this study, the following research hypotheses were formulated: (1) The neighborhood, in the sense of ecological similarity, is significant for the analysis of dependence between CO 2 emissions, economic growth, and the consumption of energy from renewable sources. (2) The countries characterized by a high share of energy from renewable sources out of the total energy consumption have less of an impact on the state of the natural environment in other countries than those wealthier but with a lower use of renewable energy.
The paper is organized as follows. In Section 2, we present a review of the literature related to the subject of our research. Section 3 presents a discussion on the tools and models that were used in the empirical analysis performed. The data are discussed in Section 4, as are the spatial distributions of the variables considered. Section 5 contains the details of the empirical results, and Section 6 summarizes the main results and presents the general conclusions. Finally, suggestions for further studies are presented.

Literature Review
The relationship between energy consumption and the emissions of pollutants has been analyzed by many researchers. Issues related to the effects of increasing total energy consumption as well as increasing the share of energy consumption from renewable sources have been discussed. These studies show that increases in energy consumption result in increases in the emissions of pollutants. Özokcu and Özdemir [15] consider this relationship on the basis of the cubic Kuznets curve, which was estimated for two groups of countries-26 highly developed OECD countries and 52 developing ones. Other authors, such as Aydin and Esen [16], Piłatowska and Włodarczyk [17], Presno et al. [18], and Yavuz and Yilanci [19] have also pointed out the negative impact of increased consumption energy on the environment. They used a nonlinear approach based on threshold analysis in their studies.
Studies that deal with the impact of renewable energy consumption on the environmental situation can be divided into two groups. The first group consists of studies in which the consumption of energy per capita was considered [20][21][22][23]; the second consists of those that considered the share of energy consumption from renewable sources [24]. In the work of Zoundi [25], 25 countries in the period of 1980-2012 were analyzed using the concept of co-integration. The same approach was presented by Zambrano-Monserrate et al. [26] with a discussion on the relationship in Brazil, by Jebli and Youssef [27], who considered the link between energy and the environment in Tunisia, as well as by Sahbi and Shahbaz [28], who focused on the countries of central-east and northern Africa. Similar analyses can be found in the works of Gill et al. [29], Sinha et al. [30], Dogan and Seker [31], and Bölük and Mert [32].
Despite the differences in the approaches to expressing energy consumption in the models used, the general results are the same. They show a positive effect of the increase in both the level and share of renewable energy consumption on the natural environment.
The research studies cited above were based on the environmental Kuznets curve, by which the role of the explanatory variable is played by an appropriate measure of economic growth. The models used were enriched with various additional explanatory variables. In a few works in this field, one can find a reference to the spatial connections between countries/regions. For example, Güçlü [33] incorporates spatial links into the Kuznets curve by analyzing the relationship between economic growth and environmental degradation for Turkish NUTS-3 regions in the years 2008-2013. The spatial environmental Kuznets curve was also used in the works by the following: Tan [34], Donfouet et al. [35], McPherson and Nieswiadomy [36], Burnett and Bergstrom [37], and Tevie et al. [38]. These researchers used simple spatial models, such as the spatial autoregressive model (SAR) and the spatial error model (SEM). In addition, Kang et al. [39], Wang et al. [40], Fong et al. [41], and Li et al. [42] used the spatial Durbin model (SDM). In their study, Li et al. [42] additionally determined the spatial direct and indirect effects resulting from changes in all explanatory variables included in the model.
In all of the above-mentioned studies, the significance of spatial connections for the relationship under investigation was indicated, and the authors formulated conclusions about the similarity of the environmental situation in the countries directly adjacent to each other. It should be emphasized that in these works, only the first-order neighborhood according to the common border criterion was considered.

Methodology
In the investigation, we used the models for pooled time series and cross-sectional data (TSCS), with particular reference to the spatial model. The basic space-time model was chosen, enriched only by spatial components, without any fixed or random effects that are characteristic of panel models. In this approach, we studied the heterogeneity of economies using the spatial trend, but for CO 2 emissions, it turned out to be statistically insignificant. We also considered the validity of using dynamic spatial models as well as dynamic spatial panel data models; however, given the insignificance of spatial effects and other diagnostics of these models, we decided to forgo them in further analysis. The justification for the use of the spatial models, that is, the models containing spatial lags of dependent or/and explanatory variables, comes from the specific interpretation of the parameters of these models, which measured the impact of changes in the variable values in neighboring observations/regions (i.e., y j , x kj ) on the dependent variable observation y i [43] (p. 34).
In classical terms, based on the concept of the environmental Kuznets curve in the variant of the quadratic function, the model describing the relationship between CO 2 Energies 2021, 14, 5829 4 of 16 emissions and GDP per capita as well as the share of energy consumption from renewable sources out of the total energy consumption takes the following form: where CO 2 denotes the carbon dioxide emissions per capita, GDP stands for the value of gross domestic product per capita, and RE is the share of renewable energy consumption. In turn, β 0 , β 1 , β 2 , and β 3 are the structural parameters of the model, and ε is its random component. All the variables have been expressed in logarithms to stabilize the variance.
A turning point can be determined for the last two of the above-mentioned relationships, indicating the level of GDP per capita at which CO 2 emissions reach the maximum value (iii) or the minimum value (iv). It is determined according to the following formula: In order to verify the validity of introducing spatial connections to our analysis, first for all the variables considered in every year the values of Moran's I have been calculated, using the following formula [44,45]: where y i denotes an observed value of the phenomenon in the region i, z means a column vector with elements z i = y i − y, S 0 = ∑ n i=1 ∑ n j=1 w ij is a sum of the corresponding elements of the weights' matrix W, and n stands for the number of regions. The matrix W of spatial connections in our study was defined based on the environmental development similarity of the countries.
The W matrix was derived from the 2017 ecological footprint value because this was the year for which the latest data was available. We decided to use the connectivity matrix constant in time due to the fact that in the period of 2013-2017, for the countries under consideration, there have been only minor changes in the ecological footprint values. Therefore, we concluded that this regularity was maintained in the following years. Thus, for the entire period of our study, the neighborhood structure remained unchanged.
The choice to use the ecological footprint as a criterion for determining the neighborhood of countries was dictated by its close relationship with the theory of sustainable development, in which special attention is paid to natural environmental protection. In addition, the level of CO 2 emissions, which was the subject of this study, is one of the main aspects of environmental pollution.
To construct the matrix W, we started by determining the distance between pairs of countries according to the following formula: where EF i and EF j are indicators of the ecological footprint for countries ith and jth, respectively. Then, the borderline level g of similarity between the countries was determined as the fifteenth percentile of all distances. This avoided the problem of excessive den- Subsequently, the non-zero elements of the distance matrix D were inverted as follows: and row-standardized to one. Finally, a block matrix of cross-sectional and temporal links between various countries in the field of environmental development was created.
In order to confirm the validity of introducing the spatial effects to model (1) the Lagrange multiplier tests (LM), in the basic and robust versions, were used. Thus, the following spatio-temporal Durbin model specification was considered: The models such as (6), thanks to the inclusion of spatial lags of the dependent variable and independent variables, allowed us to quantify the magnitude of the so-called direct and indirect effects in the short term [46] (p. 11). In this study, we were primarily interested in the indirect effects that were used to test the hypothesis whether in the area of the countries considered in terms of CO 2 emissions the spatial spillovers exist.
To see the way in which the mentioned effects were obtained, the general expression of the non-dynamic model was transformed into Equation (7) By transforming the equation to the form the following: and excluding from the matrix X t the vector regarding the variable X k , that is, X kt , the following equation was obtained: where . X t stands for the matrix from which the X kt has been removed. The expression (I − ρW) −1 (β k I N + θ k W) allowed us to determine the direct and indirect effects of the kth explanatory variable X k on the dependent variable Y. In our study, the indirect effects were determined in relation to the share of energy from renewable sources out of the total energy consumption in the neighboring regions.
The short-term effects were designated as the matrix of partial derivatives of Y with respect to the kth explanatory variable of X in spatial unit 1 up to unit N at a particular point in time, as shown in the following equation: which denotes the effect of a change of a particular explanatory variable in a particular spatial unit on the dependent variable of all other units in the short term. Similarly, the long-term effects could be determined from the dynamic model, which takes into account the time delays of the dependent and/or independent variables [46] (p. 11). The diagonal elements of the matrix (I − ρW) −1 (β k I N + θ k W) define the direct impacts of change in ith observation x k (denoted by x ik ) on y i , that is, on the values of the dependent variable in the same ith spatial unit. The average of the sum across the ith row of this matrix represents the average impact on the individual observation y i resulting from changing the kth explanatory variable by the amount across all observations-the Energies 2021, 14, 5829 6 of 16 average impact to an observation. In turn, the average of the sum in the jth column of the matrix yields the average impact over all y i observations from changing the kth explanatory variable by an amount in the jth observation-the average impact from an observation [43] (p. 37). To sum up, indirect effects as spatial spillovers were identified based on the non-diagonal elements of the matrix considered.

Data
The data used in the analysis came from three databases. First, the Our World in Data website (https://ourworldindata.org (accessed on 17 May 2021)) provided the data on CO 2 emissions per capita (CO 2 ) and the share of energy from renewable sources out of the total energy consumption (RE). Second, the World Bank (https://data.worldbank.org (accessed on 17 May 2021)) provided the GDP per capita (GDP). Third, the Global Footprint Network (https://data.footprintnetwork.org (accessed on 17 May 2021)) provided the ecological footprint by countries used to create a neighborhood matrix. All calculations and drawings were made in the program R-CRAN (version 4.0.2). Figure 1 presents the spatial distributions of carbon dioxide per capita in 2013 and 2019. In both years, the CO 2 values were distributed almost identically in the studied area. The lowest CO 2 emission values can be observed in the countries of South America, the southern part of Asia (on the Indian Peninsula and Indonesia), as well as in Southern Europe and the countries of Northern Africa. The highest values can be observed in North America (the US and Canada), northern and eastern parts of Asia, in Arab countries, as well as in Australia and New Zealand. Mostly, they are the relatively high development countries, which have a great impact on the world economy. average impact to an observation. In turn, the average of the sum in the jth column of the matrix yields the average impact over all observations from changing the kth explanatory variable by an amount in the jth observation-the average impact from an observation [43] (p. 37). To sum up, indirect effects as spatial spillovers were identified based on the non-diagonal elements of the matrix considered.

Data
The data used in the analysis came from three databases. First, the Our World in Data website (https://ourworldindata.org (accessed on 17 May 2021)) provided the data on CO2 emissions per capita ( ) and the share of energy from renewable sources out of the total energy consumption ( ) . Second, the World Bank (https://data.worldbank.org (accessed on 17 May 2021)) provided the GDP per capita ( ). Third, the Global Footprint Network (https://data.footprintnetwork.org (accessed on 17 May 2021)) provided the ecological footprint by countries used to create a neighborhood matrix. All calculations and drawings were made in the program R-CRAN (version 4.0.2). Figure 1 presents the spatial distributions of carbon dioxide per capita in 2013 and 2019. In both years, the CO2 values were distributed almost identically in the studied area. The lowest CO2 emission values can be observed in the countries of South America, the southern part of Asia (on the Indian Peninsula and Indonesia), as well as in Southern Europe and the countries of Northern Africa. The highest values can be observed in North America (the US and Canada), northern and eastern parts of Asia, in Arab countries, as well as in Australia and New Zealand. Mostly, they are the relatively high development countries, which have a great impact on the world economy.    Figures 1 and 2, it can be assumed that there is an inverse relationship between renewable energy consumption and CO 2 emissions in the areas of the surveyed countries. An exception may be highly developed countries, such as Canada, the United States, Australia, and China (in these countries, both variables have relatively high values), as well as less developed countries, such as Egypt, Morocco, and Algeria (in these countries, both variables are characterized by relatively low values).  Figure 2 shows the spatial distributions of the share of energy from renewable sources out of the total energy consumption in 2013 and 2019. The greatest share of renewable energy consumption characterized countries of both North and South America (excluding Mexico), most European countries (without Central and Eastern Europe), and China, Australia, and New Zealand. The lowest values were observed in Africa and in North and West Asia. By comparing the distributions of the variables under consideration in Figures 1 and 2, it can be assumed that there is an inverse relationship between renewable energy consumption and CO2 emissions in the areas of the surveyed countries. An exception may be highly developed countries, such as Canada, the United States, Australia, and China (in these countries, both variables have relatively high values), as well as less developed countries, such as Egypt, Morocco, and Algeria (in these countries, both variables are characterized by relatively low values).

Spatial Autocorrelation
The empirical analysis began with testing the spatial autocorrelation for the variables under consideration with established connections between countries based on the level of environmental development (ecological footprint). The level of significance was 0.05. Table 1 presents the values of Moran statistic (Moran's I) and the assessment of its statistical significance in the years 2013-2019.
Positive and statistically significant values of the spatial autocorrelation coefficient for the per capita carbon dioxide emission and the per capita GDP (expressed in natural logarithms and marked as l_CO2 and l_GDP, respectively) have been recorded for all the

Spatial Autocorrelation
The empirical analysis began with testing the spatial autocorrelation for the variables under consideration with established connections between countries based on the level of environmental development (ecological footprint). The level of significance was 0.05. Table 1 presents the values of Moran statistic (Moran's I) and the assessment of its statistical significance in the years 2013-2019. Positive and statistically significant values of the spatial autocorrelation coefficient for the per capita carbon dioxide emission and the per capita GDP (expressed in natural logarithms and marked as l_CO 2 and l_GDP, respectively) have been recorded for all the years. The positive spatial autocorrelation indicates similarity, in terms of CO 2 emissions as well as GDP, of countries with a similar level of environmental protection. The values of the Moran's I prove the strong links between countries with comparable levels of environmental development.
The situation is different in the case of the share of energy from renewable sources in total energy consumption (l_RE). The Moran statistics are statistically insignificant and indicate the lack of links, in this respect, between "neighboring" countries.
The results of spatial autocorrelation testing for the considered variables were the initial motivation for incorporating the spatial effects into the model of CO 2 emissions relative to GDP and renewable energy consumption using the Kuznets curve additionally.

Empirical Models
First, the space-time model (LM_pooled) in the form of Equation (1) was estimated and verified. The results obtained are presented in Table 2. The p-values for the parameters β 1 and β 2 indicate the significance of the impact of GDP per capita as well as its squares on CO 2 emissions. Moreover, the signs of the parameters (β 1 > 0 and β 2 < 0) allow us to conclude an inverse U-shaped relationship between GDP and CO 2 emissions. Thus, the considered relationship for selected countries of the world takes a classic shape of the Kuznets curve. The negative and statistically significant value of the β 3 parameter indicates an inverse relationship between renewable energy consumption and CO 2 emissions. Thus, an increase in the share of energy consumption from renewable sources in individual countries leads to Energies 2021, 14, 5829 9 of 16 improvement in their environmental situation. Based on the estimated Kuznets curve, its turning point was determined, amounting to $139,658.40 per capita. Taking into account the values of GDP per capita, it should be stated that none of the countries reached this ceiling during the period considered. Therefore, all the countries are on the path leading to the turning point, which may indicate a greater focus on economic development than on care for the natural environment. Figure 3 shows the shape of the Kuznets curve determined on the basis of model (1). increase in the share of energy consumption from renewable sources in individual countries leads to improvement in their environmental situation. Based on the estimated Kuznets curve, its turning point was determined, amounting to $139,658.40 per capita. Taking into account the values of GDP per capita, it should be stated that none of the countries reached this ceiling during the period considered. Therefore, all the countries are on the path leading to the turning point, which may indicate a greater focus on economic development than on care for the natural environment. Figure 3 shows the shape of the Kuznets curve determined on the basis of model (1). The Lagrange multiplier tests (the basic versions-, , and the robust ones-, ) indicate the legitimacy of supplementing model (1) with spatial connections between the countries. Therefore, the spatio-temporal Durbin model was specified (see Equation (6)). The results of estimation and verification of the model are reported in Table 3. The values of the and parameters, as in the case of the model without spatial effects, indicate an inverse U-shaped relationship between GDP and CO2 emissions. Importantly, these parameters are statistically significant. Moreover, the sign of the parameter has not changed, which, as in the previous model, indicates a positive impact of  (1)).
The Lagrange multiplier tests (the basic versions-LM SE , LM SAR , and the robust ones-RLM SE , RLM SAR ) indicate the legitimacy of supplementing model (1) with spatial connections between the countries. Therefore, the spatio-temporal Durbin model was specified (see Equation (6)). The results of estimation and verification of the model are reported in Table 3. The values of the β 1 and β 2 parameters, as in the case of the model without spatial effects, indicate an inverse U-shaped relationship between GDP and CO 2 emissions. Importantly, these parameters are statistically significant. Moreover, the sign of the parameter β 3 has not changed, which, as in the previous model, indicates a positive impact of renewable energy consumption on carbon dioxide emissions. Likewise, a negative and statistically significant parameter θ describing the effects of changes in renewable consumption in "neighboring" countries (with a similar level of environmental development) shows that its increase results in lower CO 2 emissions in a given country.
Compared to model (1), the GDP value at which CO 2 emissions started to decline decreased. In this case, the threshold value was estimated at $82,138.04. This is further evidence of a positive influence of pro-ecological neighbors' behavior on the environmental situation in a given country. It is worth emphasizing that only two countries have reached the threshold point, namely Luxembourg and Norway. The shape of the Kuznets curve, determined based on model (6), is presented in Figure 4.
Energies 2021, 14, x FOR PEER REVIEW 10 o renewable energy consumption on carbon dioxide emissions. Likewise, a negative a statistically significant parameter describing the effects of changes in renewable c sumption in "neighboring" countries (with a similar level of environmental developme shows that its increase results in lower CO2 emissions in a given country. Compared to model (1), the GDP value at which CO2 emissions started to decl decreased. In this case, the threshold value was estimated at $82,138.04. This is furt evidence of a positive influence of pro-ecological neighbors' behavior on the environm tal situation in a given country. It is worth emphasizing that only two countries ha reached the threshold point, namely Luxembourg and Norway. The shape of the Kuzn curve, determined based on model (6), is presented in Figure 4.  (6)).
The positive and statistically significant value of the autoregression paramete proves a similar level of CO2 emissions in countries with a similar level of environmen development.

Spatial Spillovers
In this subsection, we present the results of the empirical indirect effects analysis the years of the examined period, carried out on the basis of the following transformat of model (6), with respect to the spatially lagged renewable energy consumption, pressed in natural log ( ), that is, as the following formula: The indirect effects were determined in the form of the average values in the cro section of rows and, separately, in the cross-section of columns of the ( − ρ ) ( + ) matrix, excluding diagonal elements. In this way, measurements of the avera impacts (in terms of the analyzed variables) of individual countries on a given count and of a given country on other countries, respectively, were obtained. Due to the stabi of the spatial connectivity matrix over time, the spillover effects were the same in each the analyzed years. Figure 5 presents spatial distributions of indirect effects obtained. The first map Figure 5 shows the distribution of average inflows via the share of energy consumpt from renewable sources out of the total energy consumption in individual countries the CO2 emissions in a given country. The positive and statistically significant value of the autoregression parameter ρ proves a similar level of CO 2 emissions in countries with a similar level of environmental development.

Spatial Spillovers
In this subsection, we present the results of the empirical indirect effects analysis for the years of the examined period, carried out on the basis of the following transformation of model (6), with respect to the spatially lagged renewable energy consumption, expressed in natural log Wln(RE), that is, as the following formula: The indirect effects were determined in the form of the average values in the crosssection of rows and, separately, in the cross-section of columns of the (I − ρW) −1 (β 3 I N + θW) matrix, excluding diagonal elements. In this way, measurements of the average impacts (in terms of the analyzed variables) of individual countries on a given country, and of a given country on other countries, respectively, were obtained. Due to the stability of the spatial connectivity matrix over time, the spillover effects were the same in each of the analyzed years. Figure 5 presents spatial distributions of indirect effects obtained. The first map in Figure 5 shows the distribution of average inflows via the share of energy consumption from renewable sources out of the total energy consumption in individual countries on the CO 2 emissions in a given country. We can see that the countries of South America were among the ones that received transmission impulses from other countries with the highest strength. It should be noted that these countries were characterized by the lowest CO2 emissions and the highest share of renewable energy consumption. The countries that were least affected by all other countries through the transmission of renewable energy consumption included the United States, China, Russia, and Australia-the relatively highly developed economies.
The second map in Figure 5 shows the distribution of the average impacts of a given country's share of energy from renewable sources out of the total energy consumption on the CO2 emissions in all other economies. It is worth noting that countries that were the least influenced by others were the ones that most strongly affected other countries. Thus, renewable energy consumption in the United States, China, Russia, and Australia most strongly affected the CO2 emissions in other countries. Among the economies whose impact on other economies was the largest, there were also those of Italy and Norway. On the other hand, among the countries whose impact (through changes in the structure of energy consumption) on environmental pollution in other countries was the lowest, were Brazil, Algeria, Peru, and Venezuela. Figure 6 shows the impacts of two selected countries on other countries in the range of the variables considered. The maps in this figure present transmission impulse distributions resulting from changes in the structure of energy consumption in countries with We can see that the countries of South America were among the ones that received transmission impulses from other countries with the highest strength. It should be noted that these countries were characterized by the lowest CO 2 emissions and the highest share of renewable energy consumption. The countries that were least affected by all other countries through the transmission of renewable energy consumption included the United States, China, Russia, and Australia-the relatively highly developed economies.
The second map in Figure 5 shows the distribution of the average impacts of a given country's share of energy from renewable sources out of the total energy consumption on the CO 2 emissions in all other economies. It is worth noting that countries that were the least influenced by others were the ones that most strongly affected other countries. Thus, renewable energy consumption in the United States, China, Russia, and Australia most strongly affected the CO 2 emissions in other countries. Among the economies whose impact on other economies was the largest, there were also those of Italy and Norway. On the other hand, among the countries whose impact (through changes in the structure of energy consumption) on environmental pollution in other countries was the lowest, were Brazil, Algeria, Peru, and Venezuela. Figure 6 shows the impacts of two selected countries on other countries in the range of the variables considered. The maps in this figure present transmission impulse distributions resulting from changes in the structure of energy consumption in countries with the highest share of energy consumption from renewable sources, namely Norway and Brazil. the highest share of energy consumption from renewable sources, namely Norway and Brazil.
A change in the share of renewable energy in Norway had the strongest impact on environmental pollution in other Scandinavian countries, as well as in Central European countries and Russia. It can be assumed that this was due to the high degree of energy dependence on Norway of countries located close to each other in geographical space. In contrast, the countries of both North and South America, as well as South Asia, were least influenced by the changes in Norway. Changes in the share of renewable energy out of the total energy consumption in Brazil had the strongest impact on environmental pollution in most of other South American countries, China and Mexico, as well as in most Mediterranean countries. The reason for such dependencies may be the comparable, equatorial climate of the countries, where changes in the structure of energy consumption result in similar changes in terms of CO2 emissions. The least sensitive (from the environmental aspect) to changes in renewable energy consumption in the country were the United States, Canada, Norway, and Finland. Figure 7 presents the distributions of the average impacts of changes in the structure of energy consumption in countries with the strongest impact on others in terms of CO2 emissions. Based on the results obtained, it was established that such countries were Italy and China. A change in the share of renewable energy in Norway had the strongest impact on environmental pollution in other Scandinavian countries, as well as in Central European countries and Russia. It can be assumed that this was due to the high degree of energy dependence on Norway of countries located close to each other in geographical space. In contrast, the countries of both North and South America, as well as South Asia, were least influenced by the changes in Norway.
Changes in the share of renewable energy out of the total energy consumption in Brazil had the strongest impact on environmental pollution in most of other South American countries, China and Mexico, as well as in most Mediterranean countries. The reason for such dependencies may be the comparable, equatorial climate of the countries, where changes in the structure of energy consumption result in similar changes in terms of CO 2 emissions. The least sensitive (from the environmental aspect) to changes in renewable energy consumption in the country were the United States, Canada, Norway, and Finland. Figure 7 presents the distributions of the average impacts of changes in the structure of energy consumption in countries with the strongest impact on others in terms of CO 2 emissions. Based on the results obtained, it was established that such countries were Italy and China. In the cases analyzed, we found that a change in the share of renewable energy out of the total energy consumption in Italy had the strongest impact on the levels of CO2 emissions in most of the European economies taken into account, as well as in Chile. The North and East Asia countries, as well as Brazil and Argentina, were among the ones slightly less affected by Italy, whereas the smallest transmission impulses from Italy were received by the United States, Canada, Australia, and India.
China, in turn, had the strongest impact on Argentina, Brazil, Malaysia, as well as on most of the European Mediterranean countries. As with the impulses from Italy, the group of countries least sensitive to changes in the structure of energy consumption in China included the United States, Canada, and Australia.
Finally, it is worth emphasizing the similarity of the strength of influence of Italy and China on environmental pollution in South American countries as well as in most European countries. It is also worth noting the weak dependence of the level of environmental protection in the United States on changes in the structure of energy consumption in other considered world economies.

Conclusions
The results of this study underline the role of changes in the structure of energy consumption in the world economies for the improvement of the environmental situation. In the cases analyzed, we found that a change in the share of renewable energy out of the total energy consumption in Italy had the strongest impact on the levels of CO 2 emissions in most of the European economies taken into account, as well as in Chile. The North and East Asia countries, as well as Brazil and Argentina, were among the ones slightly less affected by Italy, whereas the smallest transmission impulses from Italy were received by the United States, Canada, Australia, and India.
China, in turn, had the strongest impact on Argentina, Brazil, Malaysia, as well as on most of the European Mediterranean countries. As with the impulses from Italy, the group of countries least sensitive to changes in the structure of energy consumption in China included the United States, Canada, and Australia.
Finally, it is worth emphasizing the similarity of the strength of influence of Italy and China on environmental pollution in South American countries as well as in most European countries. It is also worth noting the weak dependence of the level of environmental protection in the United States on changes in the structure of energy consumption in other considered world economies.

Conclusions
The results of this study underline the role of changes in the structure of energy consumption in the world economies for the improvement of the environmental situation.
Environmental protection has become an increasingly important element of the economic development of countries, which is reflected in the contemporary concept of sustainable development. Its purpose is to improve the state of the national economy while reducing the consumption of scarce resources.
The Kuznets curve determined for the selected countries pointed to the inverse Ushaped relationship between the per capita GDP and CO 2 emission. Including the share of energy from renewable sources in the total energy consumption as an additional explanatory variable in the models constructed confirmed the conclusions of other researchers that with the increase in this share, there was an improvement in the environmental situation, that is, the carbon dioxide emissions were reduced. Moreover, the inclusion of spatially lagged variables (i.e., the CO 2 emissions and energy from renewable sources consumption in "neighboring" countries) in the final model showed to what extent the pro-ecological actions of some economies affect others. Additionally, it can be seen that the impact of these variables on the dependent variable is smaller than their impact within a given territorial unit.
The spatial indirect effects determined based on the spatio-temporal Durbin model allowed us to identify, firstly, the countries that are most susceptible to the influence of other countries, and secondly, those with the strongest impact on others.
It is worth noting that relatively highly developed countries were among those in which the change in energy from renewable sources consumption had the greatest impact on the CO 2 emissions in other countries. This is mainly due to the fact that most of the economies failed to reach the turning point, that is, the level of GDP per capita at which the CO 2 emissions start to decline. The economies are still at a stage where the main focus is on economic development.
At the same time, the highly developed countries were minorly influenced by other countries in terms of the variables under consideration. The case of the United States should be distinguished as an economy independent of most others.
Undoubtedly, the positive impact of the changes taking place in the countries with a higher level of environmental development on the state of the environment in other countries was observed. This thesis was confirmed by the decline in GDP per capita at the turning point when relationships between neighbors were incorporated into the model.
The analysis of the spatial distributions of the impact of changes in the structure of energy consumption in Norway, Italy, Brazil, and China on air pollution in other countries leads to interesting conclusions. The mentioned European countries have a major influence on the others within the same continent, whereas impulses from economies such as Brazil and China, located on other continents, have a wider geographical scope. The countries influenced by them are not located in one cohesive area.
The results of the research show the importance of pro-ecological activities not only within a given country. The spatial spillovers in this regard are also significant.
The spatio-temporal Durbin model used in our study is only one of the possible specifications that turned out to be useful for the analysis of the phenomenon under consideration. Other model specifications should be used in further studies. Additionally, the use of other connectivity matrices should be verified. It is also worth determining the indirect effects in relation to other explanatory variables and establishing appropriate spatial regimes with regard to the wealth of the analyzed economies.

Conflicts of Interest:
The authors declare no conflict of interest.