Next Article in Journal
Using Supra-Covered Bonds to Enhance Liquidity in the Euro Area: Assessment of Advantages for the Banking Sector
Next Article in Special Issue
Nonparametric Estimation of the Ruin Probability in the Classical Compound Poisson Risk Model
Previous Article in Journal
Global FDI Inflow and Its Implication across Economic Income Groups
Previous Article in Special Issue
Propensity Score Weighting with Mismeasured Covariates: An Application to Two Financial Literacy Interventions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Environmental Kuznets Curve: A Semiparametric Approach with Cross-Sectional Dependence

Department of Economics, University of Cantabria, Avda. de los Castros 56, 39005 Santander, Spain
*
Authors to whom correspondence should be addressed.
J. Risk Financial Manag. 2020, 13(11), 292; https://doi.org/10.3390/jrfm13110292
Submission received: 21 October 2020 / Revised: 12 November 2020 / Accepted: 16 November 2020 / Published: 23 November 2020
(This article belongs to the Special Issue Nonparametric Econometric Methods and Application II)

Abstract

:
This paper proposes a new approach to examine the relationship between CO 2 emissions and economic developing. In particular, we propose to test the Environmental Kuznets Curve (EKC) hypothesis for a panel of 24 OECD countries and 32 non-OECD countries by developing a more flexible estimation technique which enables to account for functional form misspecification, cross-sectional dependence, and heterogeneous relationships among variables, simultaneously. We propose a new nonparametric estimator that extends the well-known Common Correlated Effect (CCE) approach from a fully parametric framework to a semiparametric panel data model. Our results corroborates that the nature and validity of the income–pollution relationship based on the EKC hypothesis depends on the model assumptions about the functional form specification. For all the countries analyzed, the proposed semiparametric estimator leads to non-monotonically increasing or decreasing relationships for CO 2 emissions, depending on the level of economic development of the country.

1. Introduction

The growing interest of citizens and governments in environmental degradation, generated mainly by economic growth and the excessive use of natural resources, has led economies around the world to take measures to mitigate the effects of global warming and climate change, according to what was established in the Kyoto Protocol and the 2015 Paris Agreement on Climate Change Mitigation. However, despite the efforts made, the burning fossil fuels (i.e., carbon, oil and gases) used for the production of energy necessary for economic development continues to significantly contribute to CO 2 emissions.
This situation has led the main policymakers to focus their efforts on promoting policies that pretend to achieve sustainable economic growth, that is, to achieve an increase in economic growth that is compatible with the environment. Since the relationship between economy and environment is highly complex and controversial, much of the empirical studies carried out in the literature are based on the so-called Environmental Kuznets Curves (EKC), which are still subject to an intense debate nowadays.
The EKC is based on the concept of the Kuznets curve, initially proposed by Kuznets (1955), which describes an inverted U-shaped relationship between income inequality and income per capita. This hypothesis argues that income inequality initially increases as per capita income rises and then begins to decrease from a certain threshold point. In the 1990s, the concept of the Kuznets curve was applied for the first time to the environmental quality trying to corroborate whether the relationship between income per capita and environmental degradation follows a similar relationship in terms of inverted U-shaped pattern (see Grossman and Krueger 1993, 1995 and Stern 1998, for example). In this framework, the EKC postulates that low income levels are correlated with a growing deterioration of the environment, but after a certain threshold point of the income per capita, the relationship between the two variables becomes negative again. In this way, this hypothesis is supported by the argument that higher levels of development imply a change in the economic structure in favor of industry and services, where production processes turn to more efficient and environmentally friendly technologies that help to preserve natural resources and reduce significantly the environmental deterioration. Similarly, the EKC shows the development of an economy over time. In a first stage, the economy is based on the agricultural sector with a strong impact on the quality of the environment. In a second stage, the industry is developed and, although it generates a higher level of wealth, it causes a great damage in the environmental’s quality. After a threshold point, the economy sustains its growth in efficient and cleaner technologies, mainly in the services sector. Thus, the EKC hypothesis emphasizes that the economic growth is a precondition for reducing environmental pollution (see Beckerman 1992).
Since the pioneering work of Grossman and Krueger (1991), several empirical studies have emerged in the literature with the aim of corroborating the EKC hypothesis in the income–pollution relationship. However, so far there is no consensus on this relationship (see Dinda (2004); Galeotti (2007) and Kaika and Zervas (2013a, 2013b) for intensive reviews about this issue). On the one hand, Selden and Song (1994), Grossman and Krueger (1995), List and Gallet (1999) and Stern and Common (2001), among others, point out that the intensity of pollutant emissions initially increases with the income per capita in the early stages of the economic development, but eventually falls as the income per capita rises above certain threshold point, at least in the case of the developed countries. On the other hand, Harbaugh et al. (2002) and Effiong and Oriabije (2018) highlight that there is no evidence that this relationship is robust for a number of emission pollutants.
A possible reason for this lack of consensus may be the presence of several specification errors in the models used to corroborate the EKC hypothesis. Firstly, it is common in this literature to assume ex-ante specific functional forms (i.e., quadratic or cubic polynomials) that may not fully explain the complexity of the EKC relationship. In the last two decades, nonparametric and semiparametric models are gaining popularity in order to deal with this problem (see Harbaugh et al. 2002; Bertinelli and Strobl 2005, and Azomahou et al. 2006 among others). Secondly, most of these studies completely ignore the unobservable heterogeneity of countries or regions due to economic, social, political, structural, and biophysical differences that can have variable effects on environmental quality (see Dinda 2004 or Purcel 2020 and the references therein, among others). Hence, assuming homogeneity of parameters (i.e., suggesting that the income–pollution relationship is the same for all countries) can lead to misleading inferences. Thirdly, most studies with longitudinal data use standard panel data techniques that crucially depend on the independence assumption between the individuals. However, the growing process of globalization and the importance of economic and social interconnections between economic agents make necessary to control for the unobserved cross-sectional dependence (see Apergis et al. 2017 and Martinez-Zarzoso and Bengochea-Morancho 2004, for example).
Given the lack of consensus about the income–pollution relationship based on the results obtained with fully parametric estimation techniques, this article reviews the EKC hypothesis trying to overcome the three potential misspecification problems discussed above: (i) functional misspecification; (ii) cross-sectional dependence from common factors; (iii) heterogeneous relations among variables. More precisely, in this article a new nonparametric estimation procedure is developed that allows to detect the true shape of the income–pollution relationship, allowing heterogeneous relationships among countries and controlling for the possible cross-section dependence related to the presence of a finite number of unobservable common factors (see Bai 2009 and Pesaran 2006, among others). To the best of our knowledge, this is the first paper that try to solve the three misspecification problems simultaneously with the aim of modelling the income–pollution relationship in a more efficient way and leading to more suitable political prescriptions. Finally, with the aim of being able to fully characterize the Kuznets curve, we will carry out an empirical study with both OECD and non-OECD countries. In this way, we will be able to draw conclusions for both developing countries and those that are developing.
The paper is organized as follows. In Section 2, we present a review of the literature on EKCs, focusing mainly on issues related to econometric specifications. In Section 3, the model specification and the data description are presented. In Section 4, we present the main results and in Section 5 we discuss the main empirical findings. Finally, in Section 6 we offer concluding remarks. In the Appendices, we present details about the econometric methodology of the semiparametric model that we propose and the necessary assumptions to obtain the main asymptotic properties of the proposed estimators.

2. Literature Overview on the EKC

After the seminal work of Grossman and Krueger (1991, 1995); Shafik and Bandyopadhyay (1992) and Panayotou (1993) on the EKC, the relationship between pollution and economic growth has been intensively analyzer over the past two decades from both theoretical and empirical point of view.
Despite the enormous amount of studies devoted in this topic, the agreement of the EKC hypothesis are far. Evidence of an EKC has been found for several indicators, but these findings are not unanimously accepted in the literature as it can be seen in the case of the carbon dioxide emissions. For this particular environmental measure, an EKC was corroborated in the studies of Holtz-Eakin and Selden (1995), Galeotti and Lanza (1999), Heil and Selden (2001), Cole (2004), Galeotti et al. (2006), and Parajuli et al. (2019), among others, in contradiction with the results obtained in Shafik and Bandyopadhyay (1992), Galeotti and Lanza (2001), York et al. (2003), and Nutakor et al. (2020), for example. Other researches such as Sengupta (1996) and Martinez-Zarzoso and Bengochea-Morancho (2004) find a cubic N-shaped relationship, whereas Barra and Zotti (2018) points out that preliminary evidence validates the Kuznets’s hypothesis, although this relationship turned out to be misleading once the issue of (non)-stationarity has been taking into account.
Most of these empirical studies have relied on parametric specifications (i.e., quadratic or cubic specifications). However, although most of them resort to specific panel data technique to deal with the presence of unobservable individual heterogeneity. In this situation, the risk of misspecification is considerably high for two reasons. Firstly, most of the above studies crucially depend on the independence assumption between the individuals. Hence, the presence of cross-sectional independence may invalidate most of their conclusions. Secondly, fully parametric models are subject to restrictive conditions about the functional form. When the assumptions of the empirical model and parameters are inconsistent with the real data generating process, there exists model specification problem and it can lead us to misleading inference and conclusions. In order to avoid this problem, nonparametric and semiparametric regression models have become very popular in the recent literature since they do not require the specification of a specific functional form in order to investigate the existence of EKC in the relationship between pollution and economic growth.
Among the papers that use nonparametric estimation techniques is that of Harbaugh et al. (2002), where they use a nonparametric pooled regression to examine the robustness of the evidence for the existence of an inverted U-shaped relationship between national income and pollution for a panel of countries. They find that the results are highly sensitive both to slight variations in the data and to reasonable permutations of the econometric specifications. Hence, they conclude that there is little empirical support for an inverted U-shaped relationship between pollution and national income in these data. The work in Bertinelli and Strobl (2005) proposes a more flexible semi-parametric specification to overcome the above problems and obtain that they are unable to reject a linear relationship between economic growth and environmental degradation for a panel of 122 countries over the period 1950–1990.
Azomahou et al. (2006) examine the empirical relationship between CO 2 emissions per capita and GDP per capita during the period 1960–1996 using a panel of 100 countries. In order to avoid any ad hoc choice of a parametric functional form, they consider a fully nonparametric model. Further, in order to avoid any bias related with the presence of unobserved individual heterogeneity, the propose to estimate this regression using marginal integration technique, since they propose to use the first differencing transformation in order to avoid any bias problem related to the presence of unobserved individual heterogeneity. In this framework they obtain a monotonically-increasing relationship between CO 2 emissions and economic development of the existence of an EKC for CO 2 emissions has been clearly contradicted in this more flexible specification.
Bertinelli et al. (2012) investigated the relationship between the CO 2 emissions per capita and the GDP per capita using a kernel regression estimation in a set of countries individually. They obtain that the relationship between output and pollution between countries is different. For some developed countries the relationship has been heterogeneous after 1960, whereas for almost all the developing countries it was always upward sloping. Chen and Chen (2015) analyze the nonlinear relationship between industrial pollution and economic development level based on nonparametric method for testing and verifying the EKC of carbon dioxide in China. Using standard nonparametric techniques they verify the pattern of EKC in the relationship between the industrial carbon dioxide emissions and the level of economic development. Kalaitsidakis et al. (2018) examine the relationship between total factor productivity growth and emissions using a semiparametric smooth coefficient model that allow to directly estimate the output elasticity of emissions using a panel of 17 countries for the period of 1981–1998. They are able to show that there exists a monotonically-increasing relationship between emissions and total factor productivity growth. Further, they find that the CO 2 emissions effect varies depending on a country’s emissions level.
Although the literature on the relationship between pollution and economic growth is extensive, there is still no consensus on the EKC hypothesis even when we resort to more flexible specifications such as nonparametric or semiparametric specifications. Further, more of the above studies ignores the presence of cross-sectional dependence. In the following section, we present an alternative estimation technique that enables us to overcome all these difficulties simultaneously.

3. Methodology and Data Source

3.1. Model Specification

Following the standard approach in the literature, we use the Stochastic Impacts by Regression on Population, Affluence and Technology (STIRPAT) model proposed in Dietz and Rosa (1997) as the reference analytical framework for evaluating the anthropogenic forces behind environmental change. Hence, the environmental impacts (I) can be understood as a multiplicative function of population size (P), affluence described per capita of economic activity (A), and the level of technology per unit of consumption and production (T) whose model specification is given by
I i = a P i b A i c T i d ε i ,
where ( a , b , c , d ) are the parameters of the model, ε denotes the idiosyncratic error term, and the subscript i are observational units (i.e., countries, regions) in a cross-section data.
Taking the natural logarithm of (1) and adding a quadratic term of the affluence (A) variable in line with the EKC hypothesis to capture possible existence of an inverted U-shaped relationship, the resulting panel data model is specified as follows:
l n ( E i t ) = β 0 + β 1 l n ( g d p c i t ) + β 2 l n ( g d p c i t 2 ) + β 3 l n ( p o p i t ) + β 4 l n ( e n i t i t ) + α i d t + ε i t ,
where E i t is a measure of environmental quality of country i at time t; p o p denotes the population size; g d p c is GDP per capita; e n i t denotes technology which is proxied by energy intensity to capture technology damaging effect on the environment. α i represents country-specific effect that is constat with time, and a time-specific effect d t to account for time-varying stochastic shocks that are common to all countries (including deterministics such as intercepts or seasonal dummies). All the variables in Equation (2) are expressed in natural logarithms so the estimated coefficients are interpreted as elasticities.
Further, the above setup is sufficiently general and renders a variety of panel data models as special cases. For example, if d t = 1 , the familiar fixed or random effects model is obtained. As the reader can notice, a very interesting information is obtained depending on the sign and statistical significance of the slope parameters of the income variable ( g d p c ). On the one hand, if β 1 > 0 and β 2 = 0 , then the relationship income–pollution is monotonically increasing (or decreasing if β 1 < 0 and β 2 = 0 ). On the other hand, if β 1 > 0 and β 2 < 0 , then and inverted U-shaped curve is observed for that relationship with the turning point given as E * = β 1 2 β 2 .
In this framework, standard panel data techniques can be used in order to obtain consistent estimates for the slope parameters. Nevertheless, as it was discussed previously, the empirical model (2) exhibits several weaknesses. As noted in Yatchew (1998), it is quite common that Economic Theory does not provide enough information regarding the functional form existing between the dependent variable and the covariates of the model. Therefore, assuming an ex-ante specific functional form and ignoring the possible parameters heterogeneity among countries can lead to misleading inference. Furthermore, the regression model (2) completely ignores the potential cross-sectional dependence.
In order to avoid the possible functional form misspecification problem in the above parametric framework, we propose an alternative approach that enables to relax the functional form assumptions and allows the data generating process to determine the true shape of the income–pollution relationship. Furthermore, with the aim of controlling for the potential cross-sectional dependence, we assume that it is due to the presence of a finite number of unobservable common factors (see Pesaran 2006 and Bai 2009, for example). Therefore, the semiparametric heterogeneous panel data model with cross-sectional dependence that we specify is as follows:
l n ( E i t ) = β 0 i + m i ( l g d p c i t ) + β 3 i l n ( p o p i t ) + β 4 i l n ( e n i t i t ) + α i d t + γ i f t + ϵ i t ,
where m i ( · ) is an unknown smooth function to estimate, l g d p c = l n ( g d p c ) , f t ( f 1 t , f 2 t , , f r t ) is a r × 1 vector of unobserved common effects, γ i ( γ i 1 , γ i 2 , , γ i r ) is the associated vector of factor loadings, and ϵ i t are the individual-specific (idiosyncratic) errors assumed to be independently distributed of ( d t , x i t , z i t ) . As the reader can notice, this specification suggests that the income–pollution relationship and the slope parameters can be different among individuals.
In order to obtain consistent estimates for model (3), we propose a new estimation procedure which combines nonparametric techniques with the Common Correlated Effects (CCE) approach proposed in Pesaran (2006). In the Appendix A, we present the proposed estimation procedure in detail.

3.2. Data Source

To investigate the empirical relationship between wealth-pollution, we used panel data sets consisting of 24 OECD countries and 32 non-OECD countries for the period 1980 to 2016. Countries with insufficient data on CO 2 emissions are dropped from the database. The detailed list of countries used for the estimation is collected in the Appendix C.
The data used in this study come from two main sources. Environmental degradation captured using the CO 2 emissions obtained from the International Energy Statistics of the U.S. Energy Information Administration (EIA)1. The CO 2 emissions (metric tones per capita) include burning of fossil fuels and cement manufacturing, but excludes emissions from land use such as deforestation.
The data for all other variables (population, affluence, and technology) are obtained from the World Development Indicators (WDI) of the World Bank. Taking the study in Effiong and Oriabije (2018) as benchmark, population is measured as total population, while affluence which captures economic prosperity is measured as real GDP per capita (constant 2015 US dollars). For its part, technology is measured using energy intensity. In this literature, energy intensity is usually expressed as total primary energy consumption per dollar GDP (1000 Btu per year 2015 999 US dollars).

4. Results

In order to begin our empirical analysis about the EKC in the OECD and non-OECD countries, we first examine if our panel variables contain cross-sectional dependence (CSD) using the CSD test proposed by Pesaran (2004) which follows a N ( 0 , 1 ) distribution under the null hypothesis of cross-section independence. The results of this test, which are collected in Table 1, indicate the null hypothesis of CSD is rejected at the 1 % level of significance. Therefore, it is possible to point out that all the series of the database used in this study contains CSD. At the light of this result, we also implemented the Pesaran (2007) panel unit root test, which is robust to the presence of CSD. The results are reported in Table 2 and indicate that all variables have a unit root in their levels and are stationary in their first differences. Hence, the unite root test results indicate the presence of CSD and non-stationary of the variables for both OECD and non-OECD countries.
Finally, we carry out the Friedman (1937); Frees (1995) and Pesaran (2004) tests for cross-sectional dependence. The results are reported in Table 3 and in all the cases they reject the null hypothesis of cross-sectional independence. Therefore, the presence of CSD has been corroborated for this study and ignoring this fact can seriously damage our results.
In order to account for CSD, a fully parametric model with heterogeneous slope coefficients as the following one is going to be estimated,
l n ( c o 2 i t ) = β 0 i + β 1 i l n ( g d p c i t ) + β 2 i l n ( g d p c i t 2 ) + β 3 l n ( p o p i t ) + β 4 l n ( e n i t i t ) + α i d t + γ i f t + ε i t ,
where we are working with the natural logarithm of c o 2 , g d p c , p o p , and e n i t . Therefore, if we are interesting in calculating the income level at which the turning point occurs to corroborate the EKC hypothesis, it can be done using the following expression: g d p c * = e β ^ 1 2 β ^ 2 .
Table 4 and Table 5 illustrates the results of fully parametric models with heterogeneous slopes using the Common Correlated Estimator (CCE) proposed in Pesaran (2006) for OECD and non-OECD countries, respectively. As it has been discussed previously, each country has its own economic structure so it is expected that CO 2 emissions-income relationship may differ across countries. One of the advantages of the CCE estimator is that they facilitate testing the EKC hypothesis for individual countries so we are able to identify country-specific determinants of environmental degradation.
However, in this paper we are not willing to impose restrictive conditions about the functional form, so we propose to estimate an alternate semiparametric regression model with heterogeneous slope coefficients and CSD as the following
l n ( c o 2 i t ) = β 0 i + m i ( l g d p c i t ) + β 3 i l n ( p o p i t ) + β 4 i l n ( e n i t i t ) + α i d t + γ i f t + ϵ i t .
We estimate this regression model using the alternative nonparametric estimator proposed in Appendix A, which enables us to combine nonparametric techniques and controlling for CSD, simultaneously.
Unlike the above fully parametric models which yields a unique coefficient estimator for the population parameter, one of the advantage of the non- and semi-parametric techniques is that they provide a regression plot that describes the true shape of the relationship between the dependent variable and the explanatory variable of interest while holding other regressors at a fixed point such as their means. Table 6 and Table 7 illustrates the results for the slope coefficients of OECD and non-OECD countries, respectively, while Figure 1 and Figure 2 presents the nonparametric estimates of the income-CO 2 emissions for the OECD countries and Figure 3 and Figure 4 for the non-OECD countries.

5. Discussion

Analyzing the empirical estimates in Table 4 and Table 5, we can point out that the energy intensity variable is positive and statistically significant in most of the OECD countries. This implies that higher consumption of fossil fuels in the production process will increase CO 2 emissions which in turn, will put further pressure on environmental quality. On its part, we obtain that this variable is not statistically significant in most of the non-OECD countries. For the population variable we do not get a definite conclusion. On the one hand, there are some OECD countries such as Canada, Denmark or Finland, for example, which exhibit a negative and significative effect, which implies that higher population reduces the pressure on environmental quality. On the other hand, there are some non-OECD countries such as Brazil, Costa Rica or Guatemala, among others, that show a positive and significative effect, which means that higher population exacerbates pressure on environmental quality.
Focusing now on the EKC hypothesis, it is expected that the coefficient on GDP will be positive because of a scale effect, whereas the coefficient on GDP 2 will be negative since emissions fall via the composition effects in the latter stages of development. When we analyze the results of Table 4 and Table 5 we find evidence that supports the EKC hypothesis only in six of the countries analyzed. In particular, the EKC hypothesis is supported in Finland, Korea, New Zealand (in the OECD database) and in Argentina, Ecuador, and Guatemala (in the non-OECD database). In particular, we find a peak turning point for Finland occurring at USD 30,001. Korea possesses a higher peak of USD 51,328 while New Zealand has a turning point at USD 32,728. On its part, Argentina exhibits a turning point at USD 1139, Ecuador has a peak of USD 1964, and Guatemala has a turning point of USD 466.
These findings are generally comparable to other empirical studies who find an EKC in the OECD. In particular, Galeotti (2007), Apergis et al. (2017), Churchill et al. (2018), and Dijkgraaf and Volleberg (2005), among others, find evidence supporting an inverted EKC relationship for some OECD countries, but not others. Our failure to find an EKC for Italy is consistent with Apergis (2016) and Churchill et al. (2018), but differs from Shabaz et al. (2017). Furthermore, Shabaz et al. (2017) find evidence of an EKC for Germany, which we do not. Other authors such as Churchill et al. (2018) try to examine the possibility of a N-shaped EKC for which they suggest to include the cube term of the per capita income. In this case, they obtain that there is a second turning point (i.e., the emissions begin to rise again when rich countries reach a second income tipping point2) in Australia, Canada and Japan, whereas an inverted N-shaped EKC was found for Denmark.
Despite what was obtained in Table 4, where only Finland, Korea, and New Zealand corroborated the EKC hypothesis, completely different conclusions are obtained if the plots in Figure 1 and Figure 2 are analyzed. In the 24 OECD countries analyzed, the presence of a nonlinear relationship between the CO 2 emissions and the economic development is clearly corroborated. Further, there is evidence supporting the validity of the EKC hypothesis in the majority of these OECD countries. However, for Denmark, France, Ireland, Italy, Turkey, and the United Kingdom the EKC hypothesis is rejected. In these particular cases our findings are consistent with the scale, technological and composition effects of economic growth (see Dinda 2004). These countries are characterized by a large agricultural sector which only have a marginal, or no, impact on the environment. However, globalization has allowed them to carry out the transition process from the agriculture to the industrial sector very quickly, and probably that is where it comes from the exponential growth that can be seen in the CO 2 emissions as the income increases. Furthermore, for Japan, South Korea, Luxembourg, Norway and Spain a N-shaped EKC was found.
Regarding the non-OECD countries, the nonlinear relationship between pollution and income is also corroborated. Likewise, countries can be grouped into two categories. On the one hand, the countries that offer insight into the existence of the EKC hypothesis such as Argentina, Costa Rica, Ecuador, El Salvador, Honduras, Hong Kong, Jamaica, Malaysia, Nicaragua, Nigeria, Pakistan, and Peru. On the other hand, there is another group of countries which exhibit an exponential growth as the income per capita increases such as The Bahamas, Brazil, Cyprus, Republica Dominicana, Gabon, Guatemala, Israel, Nepal, and Paraguay. This result suggests that these countries are still in an intermediate stage of development where the agricultural sector remains dominant and the industrial sector is less sophisticated, so that economic growth will normally have a scale effect on the environment. Further, a N-shaped EKC was found for China, Colombia, Panama, South Africa, and Uruguay.
Analyzing the estimators obtained for the slope coefficients of the semiparametric model, we can point that the significance of the population and energy intensity variables has increased considerably for both OECD and non-OECD countries (see Table 6 and Table 7), compared to the results for the fully parametric models (see Table 5). Furthermore, in most of the countries the effect of the energy intensity is positive, so we can conclude that higher consumption of fossil fuels in the production process will increase CO 2 emissions impacting negatively on the environmental quality. On its part, the effect of the population variable is more ambiguous and it will be depend on the particularities of the country.
In summary, the empirical results show that the nature and validity of the income–pollution relationship based on the EKC hypothesis depends on the model assumptions on functional form specification, corroborating the results obtained in Effiong and Oriabije (2018). The fully parametric estimate, that avoid the possible bias related to the CSD, only offers insight into the existence of an inverted U-shaped EKC curve in six of the 50 analyzed countries. On the contrary, the semiparametric analysis with CSD provides a more definite shape of the income–pollution relationship with flexibility in functional form specification as a non-monotonically increasing or decreasing relationship for CO 2 emissions in all the countries, depending on the level of economic development of the country. This inconsistency reiterates the econometric caveats in the literature surrounding ex-ante restrictions on the functional form specification and robustness issues. Finally, allowing heterogeneous unknown smooth functions and slope parameters enables us to characterize the individual behaviour of each of the countries directly. This result is very useful in this literature given the heterogeneity that characterizes the economic structure of the countries.
The conclusions from this study can be very useful for the policy makers in order to decide where they have to focus their effort to achieve sustainable economic growth. Our results suggest that not only poor countries, but also richer countries, face environmental pollution. It implies that economic development is not a sufficient condition to reduce CO 2 emissions. In this situation, all countries, especially developed countries because of their important resources (financial, technological, etc.), should make an effort to reduce these emissions in order to reduce global warming.

6. Conclusions

CO 2 emissions play an important role in global warming and much effort is being pout into reducing such emissions all over the world. Consequently, studying the relationship between CO 2 emissions and economic development is important in order to be able to correctly advice policy makers on how best to make economic growth compatible with the environment. With this aim, a huge international literature has emerged trying to test of the validity of the EKC hypothesis although there is still a great lack of consensus about the income–pollution relationship.
A possible reason for this lack of consensus may be the presence of several specification errors in the models used to corroborate the EKC hypothesis. In this paper, we contribute to the literature revisiting the EKC hypothesis trough an alternative more flexible estimation technique which enables to account for functional form misspecification, cross-sectional dependence, and heterogeneous relationships among variables, simultaneously. The empirical results show that the nature and validity of the income–pollution relationship based on the EKC hypothesis crucially depends the model assumptions on functional form specification.
The semiparametric estimators proposed in this paper enables to provide a more definite shape of the income–pollution relationship non-monotonically increasing or decreasing for the OECD and non-OECD countries, depending on the level of economic development of the country. On the contrary, the fully parametric estimate, only offers insight into the existence of an inverted U-shaped EKC curve in six of the 50 analyzed countries. This inconsistency reiterates the econometric caveats in the literature surrounding ex-ante restrictions on the functional form specification and robustness issues.

Author Contributions

Conceptualization, data curation, formal analysis, funding acquisition, investigation, methodology, A.S.; writing—original draft preparation, writing—review and editing, supervision, I.D. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge financial support from the Programa Estatal de Generación de Conocimiento y Fortalecimiento Científico y Tecnolótico del Sistema de I+D+i y del Programa Estatal de I+D+i Orientada a los Retos de la Sociedad/Spanish Ministry of Science and Innovation (Ref. PID2019-105986GB-C22) for the partial support of this work.

Acknowledgments

The authors are thankful for the constructive suggestions provided by the members of the evaluation board of the final project of the Economics Degree of the University of Cantabria, which improved the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Econometric Estimation. Semiparametric Approach

A semi-parametric model specifies the conditional mean of the dependent variable as two separate components, one parametric and one non-parametric. These types of models are very attractive from an empirical point of view due to their flexibility to balance precision and robustness. On the one hand, it allows us to incorporate prior information from economic theory or past experience while maintaining more flexibility in the specification of the model. On the other hand, although there is a nonparametric part that shows a slower convergence rate, the estimators obtained for the parametric part exhibit the same statistical properties as if the whole model were totally parametric. This is the so-called N -consistency property (see Robinson (1988) and Speckman (1988), for example) for cross-sectional models. Finally, semiparametric models allow to provide a solution to the well-known “curse of dimensionality” of the fully nonparametric models. Several reviews on this topic exist and we suggest the interested reader consult Ai and Li (2008), Henderson and Parmeter (2015), Parmeter and Racine (2019), Rodriguez-Poo and Soberon (2017), and Su and Ullah (2011), among others.
A semi-parametric panel data model with heterogeneous slopes and unknown functions and cross-sectional dependence is given by
y i t = m i ( z i t ) + x i t β i + α i d t + γ i f t + ϵ i t , i = 1 , , N , t = 1 , , T ,
where y i t denotes the dependent variable (i.e., the environmental quality measure in Equation (5)), x i t and z i t are p × 1 and q × 1 vectors of the explanatory variables of interest (i.e., in the case of Equation (5), x i t = ( l n ( p o p i t ) , l n ( e n i t i t ) ) and z i t = l n ( g d p c i t ) ), respectively, m i ( · ) is an unknown smooth function and β i is a p × 1 vector of unknown population parameters.
The aim of the researchers is to obtain consistent estimators of m i ( · ) and β i knowing that d t ( d 1 t , d 2 t , , d N t ) is a N × 1 vector of observed common effects (including deterministic regressors such as intercepts or seasonal dummies), α i is a N × 1 vector of unknown parameters, f t ( f 1 t , f 2 t , , f r t ) is a r × 1 vector of unobserved common factors, γ i ( γ i 1 , γ i 2 , , γ i r ) is the corresponding vector of the factor loadings, and ϵ i t are the individual-specific (idiosyncratic) errors assumed to be independently distributed of ( d t , x i t , z i t ) . In general, however, the unobserved factors f t could be correlated with ( d t , x i t , z i t ) , and to allow for such a possibility, we adopt the fairly general models for the individual specific regressors,
x i t = A i d t + Γ 1 i f t + v 1 i t ,
z i t = B i d t + Γ 2 i f t + v 2 i t ,
where A i , B i , Γ 1 i , and Γ 2 i are N × p , N × q , r × p , and r × q , factor loading matrices with fixed components, and v 1 i t and v 2 i t are the specific components of x i t and z i t , respectively, distributed independently of the common effects and across i, but assumed to follow general covariance stationary processes.
With the aim of obtaining consistent estimators for β i and m i ( · ) , in the following we will show how, with modifications, the Common Correlated Effects (CCE) approach proposed in Pesaran (2006) can be applied to a semiparametric regression model.
Let 0 p and 0 q be matrices of zero of p × p and q × q dimension, respectively, I p and I q identity matrices of p × p and q × q dimension. If we combine (A1)–(A3) and rearrange terms, we can write
y i t x i t z i t = C i d t + D i f t + m i ( z i t ) + v 1 i t + ϵ i t v 1 i t v 2 i t ,
where C i and D i are matrices of N × ( 1 + p + q ) and r × ( 1 + p + q ) dimension, respectively, of the form
C i = α i , A i , B i 1 0 0 β i I p 0 p 0 q 0 q I q and D i = γ i , Γ 1 i , Γ 2 i 1 0 0 β i I p 0 p 0 q 0 q I q .
In order to show that using suitable proxies for the unobserved factors is enough to avoid having to use initial estimates of β i , we take the cross-sectional sample averages of (A4) obtaining
y ¯ A t x ¯ A t z ¯ A t = C ¯ A t d t + D ¯ A f t + m ¯ A t + v ¯ 1 A t + ϵ ¯ A t v ¯ 1 A t v ¯ 2 A t ,
where y ¯ A t = N 1 i = 1 N y i t , x ¯ A t = N 1 i = 1 N x i t , and z ¯ A t = N 1 i = 1 N z i t . Furthermore, C ¯ A t , D ¯ A t , v ¯ 1 A t , and v ¯ 2 A t are the cross-sectional averages of C i t , D i t , v 1 i t , and v 2 i t , respectively, and m ¯ A t = N 1 i = 1 N m ( z i t ) . Let Γ ˜ ( γ ¯ , Γ ¯ 1 , Γ ¯ 2 ) , where γ ¯ = E ( γ i ) , Γ ¯ 1 = E ( Γ 1 i ) , and Γ ¯ 2 = E ( Γ 2 i ) . Following Pesaran (2006), we can premultiply both sides of (A5) by Γ ˜ and solve for f t ,
f t = Γ ˜ Γ ˜ 1 Γ ˜ y ¯ A t x ¯ A t z ¯ A t C ¯ A d t m ¯ A t + v ¯ 1 A t + ϵ ¯ A t v ¯ 1 A t v ¯ 2 A t
provided that
r a n k ( Γ ¯ ) = r ( 1 + p + q ) for sufficiently large N .
As N , v ¯ 1 A t p 0 , v ¯ 2 A t p 0 , ϵ ¯ A t p 0 , and m ¯ A t p 0 for each t under weak conditions. It follows,
f t ( Γ ˜ Γ ˜ ) 1 Γ ˜ y ¯ A t x ¯ A t z ¯ A t C ¯ A d t p 0 as N .
The result of this last line suggests that we can use λ t ( y ¯ A t , x ¯ A t , z ¯ A t , d t ) as observable proxies of the unobservable factors, f t . Therefore, we can conclude that effectively the Common Correlated Effect (CCE) approach proposed in Pesaran (2006) for fully parametric models can be applied in a semi-parametric setting with slight changes.
In this situation, we can estimate β i and m i ( · ) by augmenting the semiparametric regression of y i t on x i t and z i t with h t obtaining the following regression model
y i t = m ( z i t ) + x i t β i + λ t δ i + ϵ i t + o p ( 1 ) ,
where o p ( 1 ) captures possible approximation errors of the proxies. In addition, λ t ( y ¯ A t , x ¯ A t , z ¯ A t , d t ) is a × 1 vector of proxies, where = ( 1 + p + q + N ) .
In order to get a T -consistent estimator of β i , we follow Robinson (1988) to eliminate the unknown functional m i ( · ) . Taking conditional expectations of (A9) yields
E ( y i t | z i t ) = m i ( z i t ) + E ( x i t | z i t ) β i + E ( λ t | z i t ) δ i + o p ( 1 )
and subtracting (A10) from (A9) yields
y i t E ( y i t | z i t ) = x i t E ( x i t | z i t ) β i + [ λ t E ( λ t | z i t ) ] δ i + ϵ i t + o p ( 1 ) .
In order to obtain feasible estimators for β i , it is well-known that these conditional expectations are unknown and need to be estimated. With this aim, Robinson (1988) proposes to use (higher-order) Nadaraya-Watson kernel estimators. Later, Linton (1995) and Hamilton and Truong (1997), among others, pointed out that partial regression methods can be improved further by using local linear smoothers (see Fan and Gijbels (1996) to a deeper discussion about the desirable properties of these estimators). At the light of these results, we propose to use a local linear smoothers to estimate these conditional expectations.
Let μ 1 = E ( A i t | z i t = z ) , where ( A i t = y i t , or x i t , or λ t ). For a given point z I R q and for z i t in a neighbourhood of z, we propose to minimize the following weighted local linear least-squares (LLLS) problem for μ 1 ,
min t = 1 T [ A i t μ 1 ( z i t z ) μ 2 ] 2 K a ( z i t z ) ,
where K a ( · ) = a q K ( ( z i t z ) / a ) is a product kernel function such as K ( u ) = l = 1 q k ( u l ) , u l is the lth component of u, and a is a positive bandwidth term. Of course, a general diagonal or non-diagonal bandwidth matrix could be employed, but for the sake of simplicity, a single scalar bandwidth is used.
Using the resulting estimators for these conditional estimators in (A11) and writing the resulting expression in vectorial form yields
( I T S i ) Y i · = ( I T S i ) X i · β i + ( I T S i ) Λ + ϵ i · *
where Y i · ( Y i 1 , , Y i T ) and ϵ i · ( ϵ i 1 * , , ϵ i T * ) are T-dimensional vectors, X i · ( X i 1 , , X i T ) and Λ ( λ 1 , , λ T ) are matrices of dimension T × d and T × , respectively, and I T is a T × T diagonal matrix. Further, assuming that Z z i K a ( z i t ) Z z i is invertible, S i is a T × T smoothing matrix associated to the individual i of the form
S i e 1 ( Z z i K a ( z i 1 ) Z z i ) 1 Z z i K a ( z i 1 ) e 1 ( Z z i K a ( z i T ) Z z i ) 1 Z z i K a ( z i T ) , for Z z i = 1 ( z i 1 z ) 1 ( z i T z ) ,
where Z z i is a T × ( 1 + q ) matrix, e 1 is a ( 1 + q ) × 1 vector having 1 in the first entry and all other entries 0, K a ( z ) = d i a g ( K a ( z i 1 z ) , , K a ( z i T z ) ) is a T × T diagonal matrix. Note that ϵ i t * is the new error term which consists of three elements: (i) original error term, (ii) approximation error of the proxies, (iii) approximation error of the Taylor expansion.
By the formula for partitioned regression, the estimator of β i in (A13) is given by
β ^ i = X ^ i · M Λ ^ X ^ i · 1 X ^ i · M Λ ^ Y ^ i ·
where M Λ ^ = I T Λ ^ ( Λ ^ Λ ^ ) 1 Λ ^ is a projection matrix, for Λ ^ = ( I T S i ) Λ , X ^ i · = ( I T S i ) X i · , and Y ^ i · = ( I T S i ) Y i · .
Following a similar reasoning, the estimator of δ i in (A13) is given by
δ ^ i = Λ ^ M X ^ i Λ ^ 1 Λ ^ M X ^ i Y ^ i · ,
where M X ^ i = I T X ^ i · ( X ^ i · X ^ i · ) 1 X ^ i · .
Focusing now on the nonparametric estimation of the smooth unknown function m i ( · ) , we use the above estimator so the corresponding weighted local linear least-squares problem to minimize is of the following form
t = 1 T ( y i t x i t β ^ i λ t δ ^ i ) m i ( z ) ( z i t z ) D m i ( z ) 2 K h ( z i t z ) ,
where K h ( · ) is a product kernel defined as in (A12) and h is the new bandwidth term.
Then, assuming that ( Z z i K h ( z i ) Z z i ) is invertible, the resulting CCE nonparametric estimator of m i ( · ) is given by
m ^ i ( z , h ) = e 1 ( Z z i K h ( z i ) Z z i ) 1 Z z i K h ( z i ) [ Y i · X i · β ^ i Λ δ ^ i ] ,
where K h ( z ) is a T × T diagonal matrix defined as K a ( z ) , with h instead of a.
Under the conditions in Appendix B one can show that the semiparametric CCE estimator, β ^ i is consistent and asymptotically normal as N and T tends to infinity. More precisely, following a similar proof scheme as in Pesaran (2006), the following result is obtained.
Theorem A1.
Consider the panel data model (A1), and suppose that β i < C , Γ 1 i < C , Γ 2 i < C , Assumptions A1–A3 and A4–A7 hold, ( N , T ) (in no particular order), and the rank condition (A7) is satisfied. Then, β ^ i and δ ^ i are consistent estimators of β i and δ i , respectively. If it is further assumed that T / N 0 as ( N , T ) , then
T ( β ^ i β i ) d N ( 0 , Σ β i ) , Σ β i = Σ v 1 i 1 Υ β i Σ v 1 i 1 , T ( δ ^ i δ i ) d N ( 0 , Σ δ i ) , Σ δ i = Σ Λ 1 Υ β i Σ Λ 1 ,
where Υ β i = plim T 1 T X ˜ i · M G ˜ Ω ϵ i M G ˜ X ˜ i · , Υ δ i = plim T 1 T Λ ˜ M G ˜ Ω ϵ i M G ˜ Λ ˜ , Σ Λ = V a r ( λ t ) and Ω ϵ i = E ( ϵ i · ϵ i · ) are covariance matrices. Furthermore, M G ˜ = I T G ˜ ( G ˜ G ˜ ) 1 G ˜ for G ˜ = ( D ˜ , F ˜ ) , where D ˜ and F ˜ are vectors whose tth element are such as d ˜ t = d t E ( d t | z i t = z ) f z i t ( z ) and f ˜ t = f t E ( f t | z i t = z ) f z i t ( z ) .
Similarly, under the conditions in Appendix B one can show that the nonparametric CCE estimator, m ^ i ( z ; h ) is consistent and asymptotically normal as N and T tends to infinity.
Theorem A2.
Consider the panel data model (A1), and suppose that Assumptions A1–A9 hold. Given the T -consistency of β ^ i and δ ^ i , as T ,
T h q m ^ i ( z ; h ) m i ( z ) h 2 2 μ 2 q ( K ) t r { H m i ( z ) } d N 0 , R q ( K ) σ ϵ i 2 f z i t ( z ) ,
where σ ϵ i 2 = E ( ϵ i t 2 ) and H m i ( · ) is the Hessian matrix of m i ( · ) .
That theorem is proved following a similar proof scheme as in Musolesi et al. (2020), so it is omitted. The detailed proof of the theorem can be provided upon request.
Finally note that the estimate of the variances of the above theorems can be used to construct standard errors for β ^ i or confidence bands for m ^ i ( · ) . We use a standard multivariate kernel density estimator with a Epanechnikov kernel and the Silverman’s rule-of-thumb to chose the bandwidth.

Appendix B. Assumptions

In order to derive the asymptotic distribution of β ^ i and m ^ i ( · ) obtained in Appendix A, the following notation is used. Denoting X ˜ i · = X i · B X ( z ) and Λ ˜ = Λ B Λ ( z ) , where B X ( z ) = E ( X i · | z i t = z ) f z i t ( z ) and B Λ ( z ) = E ( Λ | z i t = z ) f z i t ( z ) . Furthermore, the following conditions are required.
Assumption A1.
(Common Effects). The ( N + r ) × 1 vector for common effects ( d t , f t ) , is covariance stationary with absolute summable autocovariances, distributed independently of the individual-specific errors ϵ i t , v 1 i t , and v 2 i t for all i, t, and t .
Assumption A2.
(Individual-Specific Errors). The individual-specific errors ϵ i t , v 1 j t , and v 2 j t are distributed independently for all i, j, t, and t . Furthermore, for each i, ϵ i t , v 1 i t , and v 2 i t follow linear stationary processes with absolute summable autocovariances ϵ i t = φ = 0 b i φ ζ i , t φ , v 1 i t = φ = 0 S 1 i φ ν 1 i , t φ , and v 2 i t = φ = 0 S 2 i φ ν 2 i , t φ , where ( ζ i t , ν 1 i t , ν 2 i t ) are ( p + q + 1 ) × 1 vectors of identically, independently distributed ( i . i . d . ) random variables with mean zero, variance matrix I p + k + 1 , and finite fourth-order cumulants. In particular,
V a r ( ϵ i t ) = φ = 0 b i φ 2 = σ i 2 σ ¯ 2 < , V a r ( v 1 i t ) = φ = 0 S 1 i φ S 1 i φ = Σ 1 i Σ ¯ 1 < , V a r ( v 2 i t ) = φ = 0 S 2 i φ S 2 i φ = Σ 2 i Σ ¯ 2 < ,
for all i and some constants σ ¯ 2 , Σ ¯ 1 , and Σ ¯ 2 , where σ i 2 > 0 and Σ 1 i and Σ 2 i are positive definite matrices.
Assumption A3.
(Identification of β i ) For each i, Ψ ^ i A = T 1 X ˜ i · M Λ ˜ X ˜ i · and Ψ ^ i G = T 1 X ˜ i M G ˜ X ˜ i · are nonsingular p × p matrices and have finite second-order moments for all i. Furthermore, sup i E T 1 X ˜ i · M G ˜ X ˜ i · 2 < C < .
Assumption A4.
(Density function). The density of z i t satisfies 0 < f z i t ( · ) < and is twice continuously differentiable in all its arguments with bounded second-order derivatives at any point of its support.
Assumption A5.
(Smoothness condition). Let Z I R q be the support of z i t . The unknown functions E ( λ t | z i t ) , E ( x i t | z i t ) , and E ( y i t | z i t ) are bounded and twice continuously differentiable at z in the interior of Z with second-order derivatives bounded.
Assumption A6.
(Kernel function). K ( u ) = l = 1 q k ( u l ) is a product kernel, and the univariate kernel function k ( · ) is compactly supported and bounded such that k ( u ) d u = 1 , u u k ( u ) d u = μ 2 ( K ) I q , and k 2 ( u ) d u = R ( K ) , where μ 2 ( K ) 0 and R ( K ) 0 are scalars and I q is a q × q identity matrix. All odd-order moments of k vanish, that is, u 1 ι 1 , , u q ι q k ( u ) d u = 0 , for all non-negative integers ι 1 , , ι q such that their sum is odd.
Assumption A7.
(Bandwidth). a and h are positive bandwidths such that as T , a + ( T a q ) 1 0 , T a 4 0 and T a 2 q / l o g ( T ) . Furthermore, as T , h + ( T h q ) 1 0 .
Assumption A8.
The map m ( · ) : Z I R is twice continuously differentiable at z in the interior of Z with second-order derivatives bounded.
Assumption A9.
(Lyapounov). For some ς > 0 , E [ | ϵ i t | ( 2 + ς ) ] exists and is bounded.

Appendix C. List of Countries

Table A1. List of OCDE countries.
Table A1. List of OCDE countries.
Australia Denmark Japan New Zealand Switzerland
Austria Finland South Korea Norway Turkey
Belgium France Luxembourg Portugal United Kingdom
Canada Ireland Mexico Spain United States
Chile Italy The Netherlands Sweden
Table A2. List of non-OCDE countries.
Table A2. List of non-OCDE countries.
Argentina Cyprus Hong Kong of China Pakistan South Africa
Bahamas, The  Dominican Republic Israel Panama Sri Lanka
Belarus Ecuador Jamaica Paraguay Thailand
Brazil El Salvador Malaysia Peru Uruguay
China Gabon Nepal Philippines
Colombia Guatemala Nicaragua Senegal
Costa Rica Honduras Nigeria Singapore

References

  1. Ai, Chunrong, and Qi Li. 2008. The Econometrics of Panel Data: Fundamentals and Recent Developments in Theory and Practice. Chapter Semi-Parametric and Non-Parametric Methods in Panel Data Models. Berlin/Heidelberg: Springer-Verlag, pp. 451–78. [Google Scholar]
  2. Apergis, Nicholas. 2016. Environmental Kuznets curves: New evidence on both panel and country-level CO2 emissions. Energy Economics 54: 263–71. [Google Scholar] [CrossRef]
  3. Apergis, Nicholas, Christina Christou, and Rangan Gupta. 2017. Are there Environmental Kuznets Curves for US state-level CO2 emissions? Renewable and Sustainable Energy Reviews 69: 551–58. [Google Scholar] [CrossRef] [Green Version]
  4. Azomahou, Théophile, François Laisney, and Phu Nguyen Van. 2006. Economic development and CO2 emissions: A nonparametric panel approach. Journal of Public Economics 90: 1347–63. [Google Scholar] [CrossRef] [Green Version]
  5. Bai, Jushan. 2009. Panel data models with interactive fixed effects. Econometrica 77: 1229–79. [Google Scholar]
  6. Barra, Cristian, and Roberto Zotti. 2018. Investigating the non-linearity between national income and environmental pollution: International evidence of kuznets curve. Environmental Economics and Policy Studies 20: 179–210. [Google Scholar] [CrossRef]
  7. Beckerman, Wilfred. 1992. Economic growth and the environment: Whose growth? Whose environment? World Development 20: 448–96. [Google Scholar] [CrossRef]
  8. Bertinelli, Luisito, and Eric Strobl. 2005. The environmental Kuznets curve semi-parametrically revisited. Economics Letters 88: 350–57. [Google Scholar] [CrossRef]
  9. Bertinelli, Luisito, Eric Strobl, and Benteng Zou. 2012. Sustainable economic development and the environment: Theory and evidence. Energy Economics 34: 1105–14. [Google Scholar] [CrossRef]
  10. Chen, Linna, and Shiyi Chen. 2015. The estimation of environmental Kuznets curve in China: Nonparametric panel approach. Computational Economics 46: 405–20. [Google Scholar] [CrossRef]
  11. Churchill, Sefa Awaworyi, John Inekwe, Kris Ivanovski, and Russell Smyth. 2018. The environmental Kuznets Curve in the OECD: 1870–2014. Energy Economics 75: 389–99. [Google Scholar] [CrossRef]
  12. Cole, Matthew A. 2004. Us environmental load displacement: examining consumption, regulations and the role of nafta. Ecological Economics 48: 439–50. [Google Scholar] [CrossRef]
  13. Dietz, Thomas, and Eugene A. Rosa. 1997. Effects of population and affluence on CO2 emissions. Proceedings of the National Academy of Sciences USA 94: 175–79. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Dijkgraaf, Elbert, and Herman R.J. Volleberg. 2005. A test for parameter homogeneity in CO2 emissions. Environmental and Resource Economics 32: 229–39. [Google Scholar] [CrossRef]
  15. Dinda, Soumyananda. 2004. Environmental Kuznets curve hypothesis: A survey. Ecology Economics 49: 431–55. [Google Scholar] [CrossRef] [Green Version]
  16. Effiong, Ekpeno L., and Alex O. Oriabije. 2018. Let the data speak: Semiparametric evidence on the environmental Kuznets curve in africa. Quality and Quantity 52: 771–82. [Google Scholar] [CrossRef]
  17. Fan, Jinqing, and Irene Gijbels. 1996. Local Polynomial Modelling and Its Applications. London: Chapman and Hall. [Google Scholar]
  18. Frees, Edward W. 1995. Assessing cross-sectional correlation in panel data. Journal of Econometrics 69: 393–414. [Google Scholar] [CrossRef]
  19. Friedman, Milton. 1937. The use of ranks to avoid the assumption of normality implicit in the analysis of variance. Journal of the American Statistical Association 32: 675–701. [Google Scholar] [CrossRef]
  20. Galeotti, Marzio. 2007. Economic growth and the quality of the environment: Taking stock. Environment, Development and Sustainability 9: 427–54. [Google Scholar] [CrossRef]
  21. Galeotti, Marzio, and Alessandro Lanza. 1999. Richer and cleaner? A study on carbon dioxide emissions in developing countries. Energy Policy 27: 565–73. [Google Scholar] [CrossRef] [Green Version]
  22. Roca, Jordi, Emilio Padilla, Mariona Farré, and Vittorio Galletto. 2001. Economic growth and atmospheric pollution in Spain: Discussion the environmental kuznets hypothesis. Ecological Economics 39: 85–99. [Google Scholar] [CrossRef] [Green Version]
  23. Galeotti, Marzio, Alessandro Lanza, and Francesco Pauli. 2006. Reassessing the environmental kuznets curve for CO2 emission: A robustness exercise. Ecological Economics 57: 152–63. [Google Scholar] [CrossRef]
  24. Grossman, Gene M., and Alan B. Krueger. 1991. Environmental impacts of a North American free trade agreement. National Bureau of Economic Research. [Google Scholar] [CrossRef]
  25. Grossman, Gene M., and Alan B. Krueger. 1993. The U.S-Mexico Free Trade Agreement. Chapter Environmental Impacts of the North American Free Trade Agreement. Cambridge: MIT Press, pp. 13–56. [Google Scholar]
  26. Grossman, Gene M., and Alan B. Krueger. 1995. Economic growth and the environment. The Quarterly Journal of Economics 110: 353–77. [Google Scholar] [CrossRef] [Green Version]
  27. Hamilton, Scott A., and Young K. Truong. 1997. Local linear estimation in partly linear models. Journal of Multivariate Analysis 60: 1–19. [Google Scholar] [CrossRef] [Green Version]
  28. Harbaugh, William T., Arik Levingson, and David M. Wilson. 2002. Reexamining the empirical evidence for an environmental Kuznets curve. Review of Economics and Statistics 84: 541–51. [Google Scholar] [CrossRef]
  29. Heil, Mark T., and Thomas M. Selden. 2001. Carbon emissions and economic development: Future trajectories based on historical experience. Environment and Development Economics 6: 63–83. [Google Scholar] [CrossRef]
  30. Henderson, Daniel J., and Christopher F. Parmeter. 2015. Applied Nonparametric Econometrics. New York: Cambridge University Press. [Google Scholar]
  31. Holtz-Eakin, Douglas, and Thomas M. Selden. 1995. Stoking the fires? CO2 emissions and economic growth. Journal of Public Economics 57: 85–101. [Google Scholar] [CrossRef] [Green Version]
  32. Kaika, Dimitra, and Efthimios Zervas. 2013a. The environmental Kuznerts curve (EKC) theory—Part A: Concept, causes and the co2 emissions case. Energy Policy 62: 1392–402. [Google Scholar] [CrossRef]
  33. Kaika, Dimitra, and Efthimios Zervas. 2013b. The environmental Kuznets curve (EKC) theory—Part B: Critical issues. Energy Policy 62: 1403–11. [Google Scholar] [CrossRef]
  34. Kalaitsidakis, Pantelis, Theofanis P. Mamuneas, and Thanasis Stengos. 2018. Greenhouse emissions and productivity growth. Journal of Risk and Financial Management 11: 38. [Google Scholar] [CrossRef] [Green Version]
  35. Kuznets, Simon. 1955. Economic growth and income inequality. American Economic Review 45: 1–28. [Google Scholar]
  36. Linton, Oliver. B. 1995. Second order approximation in the partially linear regression model. Econometrica 63: 1079–112. [Google Scholar] [CrossRef] [Green Version]
  37. List, John A., and Craig A. Gallet. 1999. The environmental Kuznets curve: Does one size fit all? Ecological Economics 31: 409–23. [Google Scholar] [CrossRef] [Green Version]
  38. Martinez-Zarzoso, Inmaculada, and Aurelia Bengochea-Morancho. 2004. Pooled mean group estimation of an environmental Kuznets curve for CO2. Economics Letters 82: 121–26. [Google Scholar] [CrossRef]
  39. Musolesi, Antonio, Juan M. Rodriguez-Poo, and Alexandra Soberon. 2020. A Semiparametric Panel Data Model With Common Factors and Spatial Dependence: Theory and Empirical Application to the Knowledge Production Function. Working paper. Santander, Spain: Department of Economics, University of Cantabria. [Google Scholar]
  40. Nutakor, Felix, Sylvestre Bizumuremyi, Jinke Li, and Wei Liu. 2020. Kuznets curve for CO2 emissions exits for Rwanda? Evidence from bootstapped rolling-window Granger causality test. Sustainability 12: 8636. [Google Scholar] [CrossRef]
  41. Panayotou, Theodore. 1993. Empirical Tests and Policy Analysis of Environmental Degradation at Different Stages of Economic Development. Working Paper WP238 Technology and Employment Programme. Geneva, Switzerland: International Labor Office. [Google Scholar]
  42. Parajuli, Rajan, Omkar Joshi, and Tek Maraseni. 2019. Incorporating forest, agriculture and energy consumption in the framework of the Environmental Kuznets Curve: A dynamic panel data approach. Sustainability 11: 2688. [Google Scholar] [CrossRef] [Green Version]
  43. Parmeter, Christopher F., and Jeffrey S. Racine. 2019. Panel Data Econometrics: Theory. Chapter Nonparametric Estimation and Inference for Panel Data Models. Cambridge, MA: Academic Press, pp. 97–129. [Google Scholar]
  44. Pesaran, M. Hashem. 2004. General Diagnostic Tests for Cross Section Dependence in Panels. Cambridge Working Papers in Economics. Cambridge, UK: University of Cambridge. [Google Scholar]
  45. Pesaran, M. Hashem. 2006. Estimation and inference in large heterogeneous panels with a multifactor error structure. Econometrica 74: 967–1012. [Google Scholar] [CrossRef] [Green Version]
  46. Pesaran, M. Hashem. 2007. A simple panel unit root test in the presence of cross-section dependence. Journal of Applied Econometrics 22: 265–312. [Google Scholar] [CrossRef] [Green Version]
  47. Purcel, Alexandra-Anca. 2020. New insights into the environmental kuznets curve hypothesis in developming and trasition economies: A literature survey. Environmental Economics and Policy Studies 22: 585–631. [Google Scholar] [CrossRef]
  48. Robinson, Peter M. 1988. Root-n-consistent semiparametric regression. Econometrica 56: 931–54. [Google Scholar] [CrossRef] [Green Version]
  49. Rodriguez-Poo, Juan M., and Alexandra Soberon. 2017. Nonparametric and semiparametric panel data models: Recent developments. Journal of Economic Surveys 31: 923–60. [Google Scholar] [CrossRef]
  50. Selden, Thomas M., and Daqing Song. 1994. Environmental quality and development: Is there a Kuznets curve for air pollution emission? Journal of Environmental Economics and Management 27: 147–62. [Google Scholar] [CrossRef]
  51. Sengupta, Ramprasad. 1996. CO2 emission-income relationship: Policy approach for climate control. Pacific and Asian Journal of Energy 7: 207–29. [Google Scholar]
  52. Shabaz, Muhammad, Muhammad Shafiullah, Vassilios G. Papavassiliou, and Shawkat Hammoudh. 2017. The CO2-growth nexus revisited: A nonparametric analysis for the G7 economies over nearly two centuries. Energy Economies 65: 183–93. [Google Scholar] [CrossRef]
  53. Shafik, Nemat, and Sushenjit Bandyopadhyay. 1992. Econometric Growth and Environmental Quality: Time Series and Cross Coountry Evidence. World Development Report Working Paper WPS 904. Washington, DC, USA: The World Bank. [Google Scholar]
  54. Speckman, Paul. 1988. Kernel smoothing in partial linear models. Journal of the Royal Statistical Socierty: Series B 50: 413–36. [Google Scholar] [CrossRef]
  55. Stern, David I. 1998. Progress on the environmental Kuznets curve? Environment and Development Economics 3: 173–96. [Google Scholar] [CrossRef]
  56. Stern, David I., and Michael S. Common. 2001. Is there an environmental Kuznets curve for shulphur? Environmental Economics and Management 41: 162–78. [Google Scholar] [CrossRef] [Green Version]
  57. Su, Liangjun, and Aman Ullah. 2011. Nonparametric and semiparametric panel econometric models: Estimation and testing. In Handbook of Empirical Economics and Finance. Edited by Aman Ullah and David E. A. Giles. New York: Taylor and Francis Group. [Google Scholar]
  58. Yatchew, Adonis. 1998. Nonparametric regression techniques in economics. Journal of Economic Literature 36: 669–721. [Google Scholar]
  59. York, Richard, Eugene A. Rosa, and Thomas Dietz. 2003. STIRPAT, IPAT and ImPACT: Analytic tools for unpacking the driving forces of environmental impacts. Ecological Economics 46: 351–65. [Google Scholar] [CrossRef]
1.
2.
According to the EKC, there is an initial increase in environmental pollution when a country’s per capita income is low, but this falls over time as income per capita grows. However, it is possible that when income increases beyond a particular threshold (i.e., the second turning point), emissions begin to rise again, suggesting an N-shaped relationship between income and environmental quality.
Figure 1. Nonparametric estimates of income-CO 2 emissions for the OCDE countries.
Figure 1. Nonparametric estimates of income-CO 2 emissions for the OCDE countries.
Jrfm 13 00292 g001
Figure 2. Nonparametric estimates of income-CO 2 emissions for the OCDE countries (Continued).
Figure 2. Nonparametric estimates of income-CO 2 emissions for the OCDE countries (Continued).
Jrfm 13 00292 g002
Figure 3. Nonparametric estimates of income-CO 2 emissions for the non-OCDE countries.
Figure 3. Nonparametric estimates of income-CO 2 emissions for the non-OCDE countries.
Jrfm 13 00292 g003
Figure 4. Nonparametric estimates of income-CO 2 emissions for the non-OCDE countries (Continued).
Figure 4. Nonparametric estimates of income-CO 2 emissions for the non-OCDE countries (Continued).
Jrfm 13 00292 g004
Table 1. Pesaran (2004)’s test for cross-section correlation of the errors in the augmented Dickey–Fuller, ADF(p), regression.
Table 1. Pesaran (2004)’s test for cross-section correlation of the errors in the augmented Dickey–Fuller, ADF(p), regression.
OECDNon-OECD
VariableCD-Testp-ValueCorrabs (corr)CD-testp-ValueCorrabs (corr)
co 2 18.850 ***0.0000.1870.44943.329 ***0.0000.3200.520
enit64.327 ***0.0000.6430.80223.423 ***0.0000.1970.452
pop95.964 ***0.0000.9490.949133.84 ***0.0000.9880.988
gdpc95.096 ***0.0000.9410.941119.88 ***0.0000.8860.886
gdpc 2 92.918 ***0.0000.9190.919119.95 ***0.0000.8870.887
Note: *** indicates statistical significance at the 1 % level.
Table 2. Pesaran (2007)’s panel unit root tests accommodating cross-sectional dependence.
Table 2. Pesaran (2007)’s panel unit root tests accommodating cross-sectional dependence.
OECDNon-OECD
LevelFirst diff.LevelFirst diff.
VariableCIPSp-ValueCIPSp-ValueCIPSp-ValueCIPSp-Value
co 2 −2.720 **0.040−4.609 ***0.000−2.118 ***0.000−4.639 ***0.000
enit−2.806 **0.021−4.408 ***0.000−2.819 ***0.000−4.848 ***0.000
pop−2.777 **0.029−2.780 **0.028−6.005 ***0.000−5.112 ***0.000
gdpc−2.748 **0.037−3.822 ***0.000−1.929 ***0.000−3.589 ***0.000
gdpc 2 −3.377 ***0.000−3.571 ***0.000−2.189 ***0.000−3.086 ***0.000
Note: *** and ** statistical significance at the 1 % and 5 % level.
Table 3. Pesaran (2004); Frees (1995) and Friedman (1937) cross sectional dependence tests.
Table 3. Pesaran (2004); Frees (1995) and Friedman (1937) cross sectional dependence tests.
Pesaran TestFrees TestFriedman Test
CountriesModelCD Testp-ValueCD (Q)p-ValueCDp-Value
OECDRE model14.237 ***0.0004.356 ***0.000164.69 ***0.000
FE model14.132 ***0.0004.350 ***0.000164.22 ***0.000
non-OECDRE model2.897 ***0.0043.913 ***0.00160.497 ***0.001
FE model4.143 ***0.0003.965 ***0.00069.389 ***0.000
Notes: RE and FE denote random and fixed effects models. *** indicates statistical significance at the 1 % level.
Table 4. Parametric estimates for the OECD countries. Fully parametric model.
Table 4. Parametric estimates for the OECD countries. Fully parametric model.
cteln(enit)ln(pop)ln(gdpc)ln(gdpc 2 )EKC HoldsTurning Points
(cte 2015 USD)
Australia−12.5270.481−0.9020.468−0.019NO
(121.47)(0.462)(2.382)(0.862)(0.043)
Austria186.1340.840 ***1.2610.207−0.010NO
(130.299)(0.230)(1.487)(2.412)(0.125)
Belgium−172.4451.009 ***3.666 **−3.232 *0.160 *NO24,343
(107.917)(0.213)(1.785)(1.790)(0.092)
Canada−1.229−0.210−8.469 ***1.762−0.092NO
(140.520)(0.252)(2.006)(1.836)(0.089)
Chile589.628 ***0.723 **−0.269−0.7580.081NO
(180.748)(0.324)(3.064)(1.163)(0.070)
Denmark675.823 ***1.761 ***−9.773 ***−5.434 ***0.286 ***NO13,359
(119.383)(0.187)(2.423)(1.839)0.094)
Finland−327.6501.055 ***−14.819 ***11.980 ***−0.581 ***YES30,001
(241.997)(0.368)(4.757)(2.871)(0.147)
France−347.018−2.689 *−7.2406.640−0.328NO
(361.849)(1.576)(11.131)(9.308)(0.495)
Ireland141.127 **0.476 **0.0920.2220.010NO
(58.307)(0.207)(0.409)(0.855)(0.041)
Italy−544.884−4.263 **9.9113.300−0.163NO
(440.651)(1.899)(11.729)(7.809)(0.404)
Japan256.529 *0.165−1.233−2.366 *0.114 *NO32,112
(147.414)(0.287)(1.874)(1.361)(0.066)
Korea52.1781.086 ***0.2231.692 **−0.078 *YES51,328
(110.724)(0.302)(1.648)(0.783)(0.042)
Luxembourg−586.2580.251 **−1.5595.030−0.212NO
(406.912)(0.107)(1.529)(0.544)(0.033)
Mexico335858 **0.892 **−3.399 ***0.542−0.025NO
(127.078)(0.333)(1.105)(0.544)(0.033)
The Netherlands−224.7190.446 *−3.2030.401−0.019NO
(138.904)(0.228)(3.476)(2.886)(0.141)
New Zealand355.141 ***0.318 *0.2641.892 **−0.091 **YES32,728
(88.296)(0.160)(0.696)(0.809)(0.042)
Norway139.5250.1943.2344.013−0.195NO
(440.084)(0.632)(4.205)(6.519)(0.292)
Portugal398.5110.223−1.767−6.166 ***0.371 ***NO4064
(320.103)(0.302)(2.739)(1.360)(0.076)
Spain37.1081.201 ***−1.287−4.932 **0.285 ***NO5727
(168.705)(0.223)(1.825)(2.057)(0.104)
Sweden−144.8400.1404.6526.014−0.285NO
(199.627)(0.176)(3.066)(3.631)(0.180)
Switzerland7.5990.560 **−5.011 *−0.1460.007NO
(135.368)(0.269)(2.505)(2.046)(0.098)
Turkey−252.856 **0.598 ***6.866 ***−2.688 ***0.179 ***NO1822
(115.558)(0.077)(0.984)(0.751)(0.046)
United Kingdom265.025 ***0.319 **−0.075−1.1620.057NO
(46.884)(0.133)(0.907)(1.007)(0.048)
United States158.253 ***0.567 ***0.893−4.785 ***0.262 ***NO9246
(51.021)(0.152)(0.702)(0.828)(0.043)
Note: ***, **, and * indicate statistical significance at the 1 % , 5 % , and 10 % levels, respectively.
Table 5. Parametric estimates for the non-OECD countries. Fully parametric model.
Table 5. Parametric estimates for the non-OECD countries. Fully parametric model.
cteln(enit)ln(pop)ln(gdpc)ln(gdpc 2 )EKC HoldsTurning Points
(cte 2015 USD)
Argentina34.066−0.178−0.6521.675 **−0.119 **YES1139
(52.373)(0.337)(4.664)(0.790)(0.054)
Bahamas279.150 ***0.590 **4.187−0.3730.018NO
(81.493)(0.265)(6.680)(5.799)(0.298)
Belarus36.7970.3110.181−1.080 *0.084 **NO618
(23.936)(0.300)(4.381)(0.573)(0.036)
Brazil−161.371 **0.744 **22.657 ***0.513−0.009NO
(74.899)(0.336)(4.840)(0.439)(0.032)
China−27.3470.733 ***−5.0840.783−0.028NO
(17.252)(0.200)(3.809)(0.872)(0.058)
Colombia172.1020.36513.6912.880 **−0.186 *NO
(116.556)(0.257)(9.669)(1.355)(0.093)
Costa Rica−135.116 ***−0.112 ***2.854 ***0.0430.017NO
(20.728)(0.021)(0.946)(0.993)(0.052)
Cyprus−130.219 ***−0.129−10.908−1.9300.140NO
(31.464)(0.144)(13.741)(1.171)(0.085)
Dominican Republic28.4060.321−43.363 **0.0660.009NO
(119.033)(0.887)(18.838)(3.103)(0.225)
Ecuador−3.0440.380 *3.74911.223 ***−0.740 ***YES1964
(34.738)(0.196)(7.211)(2.603)(0.178)
El Salvador14.6120.0362.2292.025−0.125NO
(98.948)(0.191)(2.025)(3.298)(0.201)
Gabon−19.9091.267 ***14.730 **-0.5650.036NO
(36.395)(0.290)(5.443)(3.800)(0.276)
Guatemala41.916 *0.0327.452 **10.890 **−0.886 **YES466
(22.727)(0.186)(2.734)(4.497)(0.363)
Honduras15.2400.174−1.0552.510−0.120NO
(55.567)(0.382)(1.915)(2.087)(0.112)
Hong Kong of China−127.346 ***0.1230.161−0.2080.012NO
(20.151)(0.165)(0.764)(1.456)(0.082)
Israel54.8420.680 **−30.795 ***−3.395 *0.275 *NO479
(42.606)(0.260)(10.826)(1.880)(0.140)
Jamaica−93.768 ***−0.120−0.4983.595−0.188NO
(34.106)(0.228)(2.847)(2.236)(0.133)
Malaysia−510.806 ***0.280−15.069 ***1.129−0.043NO
(72.052)(0.367)(2.607)(2.558)(0.257)
Nepal−47.8380.282 *−18.493 ***−0.8840.107NO
(52.107)(0.141)(4.298)(1.971)(0.179)
Nicaragua1051.098 *0.313−39.380−0.1340.016NO
(583.419)(0.323)(28.689)(1.075)(0.080)
Nigeria−52.200 ***0.1874.516 *2.123−0.186NO
(20.390)(0.161)(2.391)(1.339)(0.128)
Pakistan122.608−0.015−16.1782.923−0.142NO
(180.81)(0.097)(28.561)(2.939)(0.165)
Panama−178.660 ***0.147−6.672−2.268 ***0.133 **NO5044
(51.818)(0.322)(5.616)(0.732)(0.052)
Paraguay40.4140.777 *0.3240.845−0.090NO
(106.309)(0.396)(11.100)(1.625)(0.122)
Peru61.5260.362−12.416 ***1.348−0.089NO
(46.032)(0.388)(11.100)(0.925)(0.070)
Philippines186.920 ***0.3393.7724.490−0.349NO
(49.095)(0.257))(2.323)(5.178)(0.389)
Senegal−88.332−0.8824.648 **−4.3000.229NO
(149.700)(1.477)(2.076)(6.670)(0.353)
Singapore49.5801.353 ***−6.855 ***1.163−0.077NO
(42.931)(0.443)(1.537)(0.886)(0.056)
South Africa44.6010.507 *−14.216 ***−0.279−0.028NO
(56.440)(0.290)(4.032)(1.492)(0.100)
Sri Lanka−39.911 ***1.001 ***−0.8310.948−0.031NO
(14.340)(0.150)(1.258)(0.603)(0.041)
Thailand205.845***−0.276−18.2790.705−0.024NO
(57.738)(0.276)(12.234)(0.915)(0.068)
Uruguay−2.2340.566 **2.641−0.219 **0.050 **NO8935
(32.848)(0.275)(1.902)(0.095)(0.018)
Note: ***, **, and * indicate statistical significance at the 1 % , 5 % , and 10 % levels, respectively.
Table 6. Parametric estimates for the OECD countries. Semiparametric model.
Table 6. Parametric estimates for the OECD countries. Semiparametric model.
ln(enit)ln(pop) ln(enit)ln(pop)
Australia1.303 ***4.605 ***Luxembourg0.172 **−0.320 ***
(0.322)(0.921) (0.063)(0.014)
Austria0.924 ***−0.092 ***Mexico0.043−0.373 ***
(0.205)(0.025) (0.203)(0.030)
Belgium1.132 ***0.073The Netherlands0.614 ***−2.198 ***
(0.228)(0.132) (0.179)(0.307)
Canada−0.0181.317 ***New Zealand0.1970.088 **
(0.339)(0.467) (0.173)(0.034)
Chile0.720 **5.277 ***Norway0.100−0.074
(0.353)(0.467) (0.535)(0.134)
Denmark1.445 ***−0.441 ***Portugal−0.452 *1.560 ***
(0.295)(0.065) (0.248)(0.185)
Finland0.734 **−0.033Spain1.492 ***−0.012
(0.292)(0.074) (0.325)(0.064)
France−3.371 **−0.381 ***Sweden0.2270.531***
(1.340)(0.037) (0.206)(0.100)
Ireland−0.453 *0.027Switzerland0.532 **0.178
(0.240)(0.044) (0.240)(0.106)
Italy−5.387 **−1.695 ***Turkey0.578 ***−0.532 ***
(1.992)(0.529) (0.157)(0.062)
Japan0.582 ***0.067 ***United Kingdom0.194 *0.090 ***
(0.121)(0.003) (0.102)(0.004)
Korea0.270−0.287 ***United States0.712 ***0.202 ***
(0.207)(0.089) (0.238)(0.035)
Note: ***, **, and * indicate statistical significance at the 1 % , 5 % , and 10 % levels, respectively.
Table 7. Parametric estimates for the non-OECD countries. Semiparametric model.
Table 7. Parametric estimates for the non-OECD countries. Semiparametric model.
ln(enit)ln(pop) ln(enit)ln(pop)
Argentina−0.0230.780 ***Jamaica−0.563 ***2.169 ***
(0.194)(0.079) (0.210)(0.147)
Bahamas0.212−0.409 ***Malaysia−0.221−5.886 ***
(0.271)(0.012) (0.328)(0.521)
Belarus0.196−0.023Nepal0.1081.014 ***
(0.268)(0.021) (0.192)(0.043)
Brazil0.609 ***−0.003Nicaragua−0.480−0.655 ***
(0.140)(0.042) (0.365)(0.115)
China1.300 ***0.267 ***Nigeria0.161−0.389 ***
(0.243)(0.066) (0.104)(0.008)
Colombia0.544 **0.166 ***Pakistan−0.075−0.058 ***
(0.253)(0.048) (0.091)(0.017)
Costa Rica−0.098 ***−0.449 ***Panama−0.2010.989 ***
(0.028)(0.005) (0.461)(0.234)
Cyprus−0.2020.025Paraguay0.676 ***−0.269 **
(0.157)(0.083) (0.312)(0.133)
Dominican Republic0.176−5.300 ***Peru1.133 ***−0.334 ***
(0.680)(1.710) (0.275)(0.025)
Ecuador0.0751.023 ***Philippines0.433 *4.315 ***
(0.188)(0.075) (0.249)(0.360)
El Salvador0.122−0.408 ***Senegal−1.185−2.545 **
(0.165)(0.039) (1.066)(0.937)
Gabon0.887 ***3.458 ***Singapore−0.0511.366 ***
(0.278)(1.208) (0.494)(0.440)
Guatemala−0.0041.253 ***South Africa0.037−3.551 ***
(0.323)(0.034) (0.259)(0.505)
Honduras0.257−1.465 ***Sri Lanka0.1280.107 **
(0.226)(0.086) (0.287)(0.045)
Hong Kong of China0.022−1.909 ***Thailand0.0930.047
(0.151)(0.024) (0.291)(0.030)
Israel1.003 ***−0.184 ***Uruguay0.4081.180 ***
(0.275)(0.031) (0.272)(0.239)
Note: ***, **, and * indicate statistical significance at the 1 % , 5 % , and 10 % levels, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Soberon, A.; D’Hers, I. The Environmental Kuznets Curve: A Semiparametric Approach with Cross-Sectional Dependence. J. Risk Financial Manag. 2020, 13, 292. https://doi.org/10.3390/jrfm13110292

AMA Style

Soberon A, D’Hers I. The Environmental Kuznets Curve: A Semiparametric Approach with Cross-Sectional Dependence. Journal of Risk and Financial Management. 2020; 13(11):292. https://doi.org/10.3390/jrfm13110292

Chicago/Turabian Style

Soberon, Alexandra, and Irene D’Hers. 2020. "The Environmental Kuznets Curve: A Semiparametric Approach with Cross-Sectional Dependence" Journal of Risk and Financial Management 13, no. 11: 292. https://doi.org/10.3390/jrfm13110292

APA Style

Soberon, A., & D’Hers, I. (2020). The Environmental Kuznets Curve: A Semiparametric Approach with Cross-Sectional Dependence. Journal of Risk and Financial Management, 13(11), 292. https://doi.org/10.3390/jrfm13110292

Article Metrics

Back to TopTop