Environmental Kuznets Curve and the Pollution-Halo/Haven Hypotheses: An Investigation in Brazilian Municipalities

: There is much discussion on the non-linear relationship between economic growth and carbon dioxide (CO 2 ) emissions. Additionally, the effects of Foreign Direct Investment (FDI) on the environment are ambiguous, as both beneﬁcial (i.e., pollution-halo) and harmful (i.e., pollution-haven) effects were found. Therefore, the literature presents no consensus on either of these topics. This is especially problematic for developing regions, as these regions represent growing economies interested in receiving foreign investments, and their CO 2 -related research is limited. This study aims to understand the impacts of economic growth and FDI on the CO 2 emissions of S ã o Paulo state, Brazil. To perform this study, a unique dataset on regional FDI was built, and 592 municipalities were included. The analyses combine linear and non-linear estimations, and the results suggest a non-linear relationship between Gross Domestic Product (GDP) per capita and CO 2 emissions, along with a negative association between FDI and CO 2 . Finally, this study discusses possible policy implications and contributes to the international literature.


Introduction
Economic growth plays a vital role in countries' economic and human development. Over the past 40 years, the world economy has more than tripled [1]. Despite economic growth raising living standards in most countries, it was also responsible for an increase in CO 2 emissions and a depletion of natural resources. The results presented in the IPCC report confirm that global warming is largely attributed to human activities, especially those related to CO 2 emissions [2].
According to Mardani et al. [3], the majority of previous studies for the past two decades have been intensively focused on the study of the nexus between economic growth and energy use and their effects on CO 2 emissions.
Dinda [4] and Liu et al. [5] states that the investigation of the relationship between economic growth and CO 2 emissions mainly focuses on whether this relationship follows the Environmental Kuznets Curve (EKC) hypothesis. According to this hypothesis, the environmental quality deteriorates in the early stage of economic development/growth and improves in a later stage as the economy develops.
Most studies following this approach confirm the existence of an inverted U-shaped (and other curve shapes) for the relationship between growth and emissions (i.e., EKC) in implies an extensive manual check of approximately 22,000 export records to identify the source of capital of a company in a given municipality. Such data allows an indepth analysis of FDI in the region and may serve as a base for future research worldwide.
We choose Sao Paulo for three reasons. First, Brazil has received little attention from the international literature. The absence of empirical studies in Brazil is particularly worrisome considering the country's inward FDI stock. The Sao Paulo state is particularly relevant for a more in-depth analysis of FDI in Brazil due to foreign investments in the region. According to the Central Bank of Brazil [30], the state accounted for roughly 39% of foreign investments in the country, at least considering industrial enterprises. Besides, Sao Paulo presented the highest Gross Domestic Product (GDP) in 2016 and is the richest and most economically complex region in Brazil [31]. Second, Brazil is the leading air polluter in Latin America and is a country that experienced significant growth in recent years [32,33], with 441.8 million tons of CO 2 emissions in 2018 [34]. Sao Paulo state is responsi-ble for 71.50 × 106 tCO 2 /year [35], which is remarkable to control environmental damage in the country. In addition, Brazilian policymakers do not fully understand the environmental impacts of such investments. Third, there is an absence of studies using subnational data for developing countries.
Accordingly, our study presents three contributions to the literature. First, we present an in-depth CO 2 -related analysis in a developing region, which is important for local policymaking, especially for policymakers being pressured to promote economic growth with minimum environmental damage. Policymakers in other parts of Brazil and the world may access our policy recommendations. Second, due to the scarcity of data, we built a unique dataset with information about multinational companies located in Sao Paulo state using an export-oriented FDI measure similar to Moralles and Moreno [36] and Polloni-Silva et al. [37]. Third, our findings reveal the impact of FDI on CO 2 emissions at a municipal level, therefore contributing to the international debate on this issue.
The rest of this paper is organized as follows. The literature review on the impacts of FDI and the possible non-linear relationship between economic growth, FDI, and emissions is presented in Section 2. Section 3 focuses on the data gathering process and the estimation strategy. Section 4 presents the main findings. Section 5 presents the robustness checks. Finally, Section 6 presents the policy implications of our findings, along with a brief conclusion.

Pollution-Halo and Pollution-Haven
According to the literature, FDI may boost the host country's economic growth while putting the environment at risk [16,20,38]. Although previous research has found inconclusive results, scholars have attempted to explain the influence of FDI on the environ-ment, as it is often assumed to be connected to economic growth [39,40].
Nonetheless, the literature points to two alternative theories that associate FDI and the environment. Namely, the "pollution-haven" and "pollution-halo" theories. The pollution-haven hypothesis defends that inward FDI is attracted to weak regulatory and environmental laws presented by the host country [17,41]. Thus, FDI should negatively affect the environment, especially in developing economies with weaker environmental standards [42]. In fact, several countries are eager to attract FDI [21] and may use their lax environmental regulations to compete in the "FDI tournament" [38]. On the contrary, the pollution-halo hypothesis claims that foreign firms, especially from developed economies, arrive with newer and cleaner technologies and better management practices capable of improving the environmental quality of the host country [17,21]. The pollution-halo hypothesis is especially relevant when domestic firms are considered pollution-intensive when compared to international standards [17].
Concerning the pollution-haven and pollution-halo studies, the results are debatable, once again. Zafar et al. [21] found evidence of FDI bringing advanced technology and in-novative products in the United States, which decreased environmental pollution as the pollution-halo theory suggests. Additionally, Jiang et al. [43] claim that FDI improved air quality in China, especially in populated cities.
By employing cointegration analysis, Nepal et al. [44] found that an increase in FDI inflows may reduce energy use, evidencing that energy-efficient techniques through FDI can reduce carbon emissions, which reinforces the importance of policies that encourage FDI inflow in renewable energy sectors. Evidence of pollution-halo in developing economies is also found by Tiba and Belaid [15] and Yi et al. [23].
However, FDI may present a significant positive impact on environmental degradation, thus supporting (especially in developing economies) the pollution-haven hypothesis. In this sense, He and Yao [16] found (at least to some extent) evidence to pollution-havens in China, while similar results are also presented by Jiang [8], Shahbaz et al. [20], and Opoku and Boachie [19]. Furthermore, Kastratović [42] studied 63 different de-veloping countries, including Brazil, and found evidence of pollution-havens. Zakarya et al. [45] studied the BRICS and also found evidence of pollution-havens, as FDI increased CO 2 emissions in the selected countries.
Yet, the results presented by Shao et al. [46] refuted the pollution-haven hypothesis for the BRICS. Thus, in essence, the literature still lacks conclusive results. As Nasir et al. [18] declared, FDI may have either a positive or negative impact on the host's environ-ment. It depends on the research method, timeframe, and region [16,18,25]. De facto, in-conclusive or mixed results for developing countries were also found by Nadeem et al. [47], Chen and Yang [48], and Xu et al. [49]. Therefore, further research is essential to better comprehend the link between FDI and the environment.
Thus, considering the current literature on FDI, environment, and location choice, we formulate the following hypothesis: Hypothesis 1 (H1). The effect of FDI on CO 2 emissions is statistically significant.

Non-Linear Relationship between Economic Growth, FDI, and CO 2 Emissions
Global warming is one of the main concerns of world globalization in terms of industrial development. There has been a steady growth in energy demands to sustain the development of countries. At the same time, there are increasing concentrations of greenhouse gases (GHGs) in the atmosphere, particularly carbon dioxide (CO 2 ), intensifying the natural greenhouse effect and the heat it generates. According to IPCC [2], human activities are estimated to have caused approximately 1.0 • C of global warming above preindustrial levels.
Since the 1990s, some authors have come to defend that the interaction between economic growth and environmental problems has a certain regularity. This empirical evidence first emerged by Grossman and Krueger [50,51], in which the relationship between GDP per capita and emission of pollutants takes the form of a "U-inverted" or Environmental Kuznets Curve (EKC). Therefore, their work linked economic growth to environmental degradation in a non-linear way.
The EKC theory asserts that environmental degradation increases in the early stages of economic growth, and decreases when the income surpasses a certain threshold level [16,50]. In other words, the EKC theory suggests that at an early stage of economic development, people will demand more energy and natural resources, while neglecting the environment. However, at a later stage of development, people may present a different mindset and demand a cleaner environment, energy-efficient products, and the protection of natural resources, as well as the environment, in general [18,21]. In summary, EKC represents the notion that despite the environmental degradation at first, economic growth can improve the environmental quality in the long-run [18]. As FDI is often assumed to be connected to economic growth [39,40], multiple scholars are now analyzing FDI and its effects on the environment from the EKC perspective.
However, results are still inconclusive [18]. Both country-specific and diversified studies have been published. Fosten et al. [52] found significant evidence of EKC in the UK. Esteve and Tamarit [53] studied Spain and found that, in the long-run, economic growth reduces CO 2 emissions. However, Ghosh [54] found no conclusive evidence of this longrun relationship between growth and emissions in India, and similar results were found for Turkey [55], where pollution was not found to be associated with FDI or economic growth. Instead, Haug and Ucal [56] examines the effects of foreign trade and FDI on CO 2 emissions in Turkey and the results show that FDI has no statistically significant long-run effects and CO 2 intensity is not influenced by FDI, but an environmental Kuznets curve is present for CO 2 .
More recently, Nasir et al. [18] found evidence of EKC in Asia, although their results were not statistically significant. Contrarily, the authors of [57] validated the EKC phenomenon in Malaysia. Furthermore, Onafowora and Owoye [58] analyzed the EKC phenomenon in selected economies such as Brazil, Mexico, China, Japan, and South Korea (among others) and only found evidence of EKC in developed economies. Malik et al. [59] studied the case of Pakistan and confirms the EKC hypothesis, and the symmetric results show that economic growth and FDI intensify carbon emission in both the long and short-run and Pakistan is a pollution haven for FDI.
The results of Liu et al. [5] revealed that there is a reverse U-shape relationship be-tween economic growth and carbon dioxide emissions in China, and an inverted Nshaped relationship between FDI and carbon dioxide emissions. However, he suggested reexamining the results by grouping different regions.
Considering the literature, it is possible to infer that the results on the non-linear relationship between economic growth and environmental degradation are still inconclusive [16].
Considering the current literature on economic growth, EKC, and FDI, the following hypotheses were formulated: Hypothesis 2 (H2). Economic growth presents a non-linear relationship with CO 2 emissions; Hypothesis 3 (H3). FDI presents a non-linear relationship with CO 2 emissions.

Data and Method
Our sample includes 592 out of the 645 municipalities in São Paulo state. We excluded 1 municipality for the lack of industry-related data and 52 small-sized municipalities (i.e., averaging 3000 inhabitants each) for the lack of energy-related data. Consequently, the employed sample comprises roughly 99.3% of the state's population. Also, the data ranges from 2010 to 2016, resulting in a seven-year period and 4144 observations.

CO 2 Emissions
This study employs CO 2 emissions per capita as the dependent variable. However, a detailed (i.e., municipal-level) CO 2 database is not available at this moment. Previous scholars also faced this issue when studying CO 2 at a more micro-level. For this reason, we manually calculate the emissions from energy consumption for the municipalities of São Paulo, which is a common procedure in CO 2 -related studies.
Therefore, we use data on electricity consumption and fossil fuels. Namely, the calculation includes gasoline, diesel, liquefied petroleum gas, fuel oil, aviation gasoline, and aviation kerosene. By following this calculation, this study is aligned with previous research, employing data from the Brazilian Energy Research Office (EPE) and the National Agency for Petroleum, Natural Gas and Biofuels (ANP). Yet, the Brazilian fuel legislature makes it inappropriate to calculate CO 2 emission directly. Thus, this study followed the IPCC [60] guidelines and subtracted the volume of ethanol added to gasoline C, which is commonly commercialized in the country, to use gasoline A in our calculations. The volume of biodiesel added to diesel oil was also subtracted.
Finally, the CO 2 emissions are calculated using the conversion factors presented in Table 1.

Explanatory and Control Variables
This study is interested in better understanding the general effects of the development of São Paulo state on the environment. Accordingly, the regions' GDP per capita is employed as a measure for economic growth, as done by previous CO 2 -studies [7,61,62]. Although alternative variables (e.g., income levels) could provide valuable information on the overall richness of a region, there is a general consensus on the importance of GDP as the primary measure of how an economy or region is performing.
Additionally, in recent years, FDI received considerable recognition in many CO 2related studies [5,15,62]. However, similar to the CO 2 -issue presented earlier, São Paulo state does not present a regional FDI database. De facto, there are no governmental statistics on FDI at the state or city-levels [36], and this lack of data virtually explains the current absence of Brazilian studies on FDI and its regional effects, except for Moralles and Moreno [36] and Polloni-Silva et al. [37].
Therefore, we build a proxy for the regional FDI presence in each municipality using export-data from the Brazilian Integrated System of Foreign Trade (SISCOMEX). This sys-tem provides information on all exporting companies operating in São Paulo state, including their names, addresses, and export band value (e.g., up to 1 million BRL, between 1 and 5 million BRL, between 5 and 10 million BRL). With this information, a research group manually verified if these companies presented international origins (i.e., if the headquarters is located in a foreign country). For the whole São Paulo state during the 2010-2016 period, the research group performed approximately 22,000 individual checks. Hence, we follow Moralles and Moreno [36] in building a proxy variable for the intensity of FDI in a region, following Equations (1) and (2) below. And, where FDI is the intensity of FDI for every region j at the year t, W is the adjustment weight regarding the exports from foreign companies, MNC is the number of multinational companies (i.e., foreign companies), TC is the total number of exporting companies, k is the number of foreign companies, SFE is the sum of exports from foreign companies, and STE is the sum of total exports (including national and foreign companies). The export band values (i.e., 1 million; 5 million; 10 million; 50 million, 100 million, and 320 million) are used to calculate SFE and STE, similar to Moralles and Moreno [36]. This FDI measure varies from 0 (no foreign companies) to 1 (all exports are dependent upon foreign companies). Thus, this measure represents the intensity of FDI for every municipality in every year of the sample and shows the relative importance of FDI for each region. With this information, our study will be the first to evaluate the municipal-level effects of FDI on the environment in a Brazilian region.
It should be noted that this FDI measure is not perfect as it excludes foreign companies commercializing their products exclusively within the country. However, this measure is arguably a valuable proxy for FDI considering that many foreign companies use Brazil as an export station to other Mercosul (Mercado Comum do Sul) countries [63]. As an example, roughly 22% of Paraguay's imports and 24% of Argentina's imports come from Brazil, and these products involve a multitude of industrial sectors (e.g., machinery, chemicals, clothing, cars) [64]. As the SISCOMEX system presents data on all exporting companies from all sectors, our FDI measure is arguably capable of including a considerable share of all foreign companies in the country. Additionally, a similar export-oriented proxy for FDI was already employed by Moralles and Moreno [36] and Polloni-Silva et al. [37].
Furthermore, regarding the control variables, previous scholars have employed a multitude of demographic, economic, and technology-related variables in previous studies. For this study, the models include the industrial and service sectors' share of total GDP, the share of electricity consumption in total consumption, and the population density. Therefore, our model is similar to previous employed models [9,61,62,65]. Table 2 presents the descriptive statistics and the sources used.

Explanatory and Control Variables
We start the analysis by verifying the linear effects of the explanatory variables on the regions' CO 2 emissions. Thus, the primary model can be written as: where CO 2 is the dependent variable for every municipality j at time t, GDPPC is the real GDP per capita, FDI is the aforementioned proxy for regional FDI intensity, X is the set of control variables, a j is the regional fixed effect, and e is the error term. All variables are employed in their natural-log form, except for FDI as the variable is already in percentage form. In addition to the linear effects of the other explanatory variables, this model formulation allows us to verify if FDI presents a pollution-haven or pollutionhalo behavior.
However, this study also aims to test for possible non-linear relationships between economic growth, FDI, and CO 2 emissions. Thus, by following many previous scholars, we add quadratic terms to our model: And, with these models, it is possible to test for U or inverted-U curves between GDP, FDI, and the CO 2 emissions, thus verifying the existence of EKC curves and better understanding how future development may affect the environment if São Paulo continues to grow and attract foreign investments. Furthermore, it is necessary to consider possible non-spherical disturbances that are commonly found in panel datasets. Initially, we verified if our baseline models present collinearity issues (see Appendix A), which is a common issue in environmental studies [66]. Moreover, the Hausman test was employed to examine which specification (i.e., random or fixed-effects) was better suited for our models. Additionally, the modified Wald test for group-wise heteroscedasticity [67], the Wooldridge test for autocorrelation [68], and the Pesaran test for cross-sectional dependence [69] were used. In sum, these tests' results show that the fixed-effects model was the preferable specification and that the sample shows heteroscedasticity, autocorrelation, and cross-sectional dependence problems.
To deal with these issues, CO 2 -studies have been using the Driskoll-Kraay [70] estimator [9,61,71], as this technique deals with heteroscedasticity and autocorrelation problems and is robust to general forms of cross-sectional (spatial) and temporal dependence with panels of all sizes, including an N > T dataset [72]. Therefore, the fixed-effects Driscoll-Kraay (DK) estimator is the main technique employed in the analyses. Yet, the sam-ple is limited as a higher time frame is not available at the moment. Thus, to take ad-vantage of the higher degrees of freedom as the random-effects specification has, the Feasible Generalized Least Squares (FGLS) technique was employed to double-check our findings. FGLS is usually a random-effects estimator capable of dealing with heterosce-dasticity, autocorrelation, and cross-sectional correlation [68,73].

Main Results
Our main findings are presented in Table 3, which shows the linear and quadratic relationships between our explanatory variables and the regions' CO 2 emissions. As shown, economic growth increases CO 2 emissions, as expected. Moreover, Table 3 shows that a non-linear relationship between GDP per capita and CO 2 emissions exists, as the quadratic term for GDP per capita is significant and negative. Thus, Table 3 suggests that after a certain threshold, the overall growth/richness of a region, here represented by the GDP per capita, helps to decrease the emission level.
Additionally, FDI has a linear and negative relationship with CO 2 emissions, which indicates that FDI presents a pollution-halo behavior. Thus, these results are aligned with some previous studies defending the pollution-halo hypothesis for FDI such as Jiang et al. [74], Zhang and Zhou [25], and Shao et al. [46]. Indeed, FDI seems to be environmentally beneficial to the region. Furthermore, a non-linear relationship between FDI and CO 2 is not found, as the quadratic term is not significant. Yet, before declaring support for any hypothesis, further tests will be made. Regarding the control variables, Table 3 shows that both the industrial and the service sectors contribute to higher emission levels. Again, this demonstrates that economic growth is environmentally costly. The results demonstrate that the industrial sector is positive and significant, while FDI is negative and significant. Thus, the overall industrial sector increases the emission levels, while a higher FDI intensity is beneficial.
Likewise, the service sector, which is ignored by many scholars, also increases pollution. This result is expected considering that the CO 2 measure employed by our study is calculated using energy consumption. De facto, the service sector is expanding in the region, and many services, particularly in highly urbanized regions, involve road transportation, which generates higher levels of pollution [75]. In addition, a higher share of residential electricity consumption in the total consumption presents negative coefficients, indicating that higher shares of consumption by other sectors (e.g., industry) result in higher CO 2 than the residential sector. In other words, the impact from the residential sector is lower than other sectors. Also, regions with higher population density levels present higher emissions, as these regions are arguably larger and more populated.
To further evaluate if the results from Table 3 are reliable, we split the sample into subsamples. First, we estimated the models exclusively with the municipalities with foreign companies operating during the 2010-2016 period, therefore excluding smaller (i.e., less populated/developed) regions. Moreover, we excluded the city of São Paulo and the São Paulo metropolitan region, as these municipalities represent significantly welldeveloped and populated regions. São Paulo city and its neighborhood are responsible for a significant share of the state's GDP and CO 2 emissions and therefore represent dissident regions compared to the rest of the state. In sum, results from Table 4 confirm the effects of GDP and FDI on the environment and are aligned with Table 3. The results from Table 4 confirm the findings in Table 3. In sum, GDP seems to present a non-linear relationship with CO 2 . An inverted U-shaped relationship between GDP and CO 2 seems to represent the Brazilian reality similar to the original Simon Kuznets [76] curve. Moreover, FDI is, again, linear and negatively associated with emissions, thus being environmentally beneficial. Once more, before declaring support for the hypotheses, we perform robustness checks.

Robustness Checks
The robustness checking section is twofold. First, we verify if endogeneity is an issue. Second, we evaluate if a non-linear relationship is significant using the threshold regression technique, as some scholars presented criticism towards the use of quadratic varia-bles to test for non-linearity.
The main explanatory variables (GDP per capita and FDI) were tested for endogeneity. Sustainability studies showed that economic growth and environmental degradation may suffer from simultaneity problems. Additionally, the international business literature argues that foreign companies do not choose their subsidiaries' locations randomly. Foreign firms may be attracted to some of the host's characteristics [26]. If foreign firms are choosing to invest in more developed regions, and these regions generate more pollution, then the model may suffer from reverse causality.
Therefore, the C-statistic test (also known as the difference-in-Sargan test) was employed to verify endogeneity issues [77] and we employed the lagged values as instruments, following Yue et al. [78] and Song [79]. Previous scholars argue that finding instruments to employ in CO 2 -related studies is difficult, but defend that past outputs should not impact present emissions. Therefore, one-year lagged values of GDP per capita and FDI were used. Results (see Appendix B) show that FDI is exogenous, but endogeneity is significant for GDP per capita.
To mitigate this issue, we employed the Limited Information Maximum Likelihood (LIML) estimation technique. Previous studies commented that the LIML technique is better than the traditional Two-Stages Least Squares (2SLS) technique, and more capable of ensuring consistency and asymptotic normality [80]. LIML presents good results even in smaller samples [81] and reduces the bias of the estimations [82].
The results presented in Table 5 demonstrate that our linear results are consistent, as  they are aligned with the previous estimations in Table 3. GDP per capita is positively associated with CO 2 emissions. Conversely, FDI is negatively associated with emissions.
Furthermore, there is the issue of possible non-linear relationships between our explanatory variables and CO 2 emissions. Previously, we employed the quadratic forms of GDP and FDI in the models, as it is done by the majority of EKC studies [7,15,62]. Recently, however, scholars have commented on the limitations of this method. Girma [83] argues that, although the quadratic specification tends to be more informative than the linear model, it establishes a very specific functional form on the non-linearity that may not be a good representation of the pattern in the data. He and Yao [16] pointed out that, although EKC is traditionally measured in quadratic or even cubic forms, this specification is arbitrary and subjective. For these reasons, scholars are now using more novel methods such as nonparametric regressions and the Non-linear Autoregressive Distributed Lag (NARDL) in non-linear studies [84]. Particularly, the threshold regression is being employed by CO 2 -related studies [16,85,86].
Thus, this study also investigates the non-linear behavior of the variables using the fixedeffect panel threshold regression proposed by Hansen [87] as presented in Equations (6) and (7). ln CO 2 jt = β 0 + β 1 ln GDPPC jt I GDPPC jt ≤ λ + β 2 ln GDPPC jt I(GDPPC jt > λ) + β 3 FDI jt + β 4 ln X jt + a j + e jt (6) And, ln CO 2 jt = β 0 + β 1 ln FDI jt I FDI jt ≤ λ + β 2 ln FDI jt I(FDI jt > λ) + β 3 GDPPC jt + β 4 ln X jt + a j + e jt (7) where I(·) is an indicator function. We follow Girma [83] and employ a set of quantiles (1%, 1.25%, 1.50%, . . . , 98.75%; 99%) to estimate the threshold parameter (λ). Then, we test the null hypothesis of no non-linearity (H 0 : β 1 = β 2 ) using the likelihood ratio test statistics and their bootstrapped p-values. Table 6 shows that a non-linear relationship is only significant for the GDP-CO 2 relationship, with one significant threshold point.  Therefore, the threshold tests confirm the linear relationship between FDI and CO 2 , and FDI seems to be negatively associated with emissions. Therefore, the results support H1a and reject H1b. Table 7 shows the estimated coefficients for GDP following two distinct regimes: lower and higher levels of GDP. These results suggest that, indeed, there is a non-linear relationship between GDP and CO 2 . After the threshold point, the impact of GDP on the environment is lower, meaning that richer regions (i.e., more developed) tend to have a smaller environmental impact, which is beneficial to the region. Nonetheless, the threshold regression does not support the EKC hypothesis, as GDP still presents positive effects in both regimes. Hence, Table 7 shows support for H2 and rejects H3. São Paulo state does not present the inverted-U curve that European [88] and Asian [5,7,62] regions present. In sum, the GDP-CO 2 relationship is not linear or quadratic. The policy implications of these findings will be discussed in the next section.

Conclusions and Policy Implications
The findings show that the relationship between economic growth and emissions is environmentally costly and non-linear. To some extent, GDP increases CO 2 emissions, which corroborates with previous studies' understanding [21,89]. Other studies proposed an inverted U-shape curve between economic growth and emissions [5,8,15,19]. Our study contributes by showing that GDP is non-linear with CO 2 emissions, although this relationship did not present the U-shape form previously proposed. In this sense, we found that GDP becomes less environmentally damaging after a certain threshold point. This result supports hypothesis 2 (H2: Economic growth presents a non-linear relationship with CO 2 emissions).
Additionally, the results reveal a negative impact of FDI on CO 2 emissions, therefore supporting hypothesis 1 (H1a: FDI is negatively associated with CO 2 emissions). This result corroborates with previous studies that advocate in favor of FDI to reduce pollutant emissions [6,21,23,78]. In contrast, our main finding presents evidence against studies arguing that foreign investments increase environmental damage [7,8,16,19,62,90]. The dataset we employed shows that the majority of foreign companies operating in São Paulo state have origins in developed economies (e.g., United States, Germany, and Japan), and are operating in medium and high-tech sectors (e.g., cars, energy, mobile technologies, and electronics). Therefore, the negative impact of FDI on CO 2 emissions might result from multinational companies with green practices and higher levels of productivity. This is relevant because public authorities in developing regions must create specific policies to attract international enterprises to generate green growth conditions. We argue that FDI must be one of the major drivers to promote economic growth along with sustainable development.
There are two main useful contributions for policymakers. First, the non-linear effect of GDP on emissions might encourage public authorities to find alternative economic growth strategies. Arguably, local policymakers have an essential role in boosting economic growth while reducing power consumption, as energy policies should be managed at the local level [91]. Second, FDI is an essential policy tool to create green development paths in Sao Paulo state. In this sense, local policies must focus on strategies to attract green technologies [92], clean energy sources [93], green jobs and products [94], and higher productivity [95]. Note that FDI is crucial to promote economic growth because it guarantees access to foreign firms' knowledge and technologies to transform the economic structure. For example, the industry is a relevant sector to the Sao Paulo region. Public authorities could reduce taxes to national and multinational companies that emit less CO 2 . A tax policy could also lead green companies to regions that are already presenting higher levels of CO 2 , such as the Sao Paulo city and densely populated regions [96]. New firms might also create local new jobs, increase the local income, and generate green economic growth paths. This new cycle would generate development with lower CO 2 emissions. Besides, alternative transport policies are required in metropolitan areas, such as efficient public transport for passengers, vehicle utilization improvements, and fuel efficiency to reduce pollutant emissions [75,97]. The public sector might invest in environmental projects to develop national green products, such as photovoltaic cells, smart grids, bioplastics, among others [98]. Also, the vast regions in Sao Paulo countryside could focus on alternative energy resources [99], such as the electric energy from sugarcane mills or by building photovoltaic power plants.
Noteworthy, Sao Paulo is the most developed state in Brazil. It means that Sao Paulo has adequate infrastructure (roads, railways, ports, and airports) to connect the region with the entire world. Also, Sao Paulo presents one of the highest human development in Brazil, which is crucial to create human capital (education, health, income, among others). These conditions favor the attraction of FDI to Sao Paulo state. Other emerging regions may use the policy recommendations presented here to facilitate their development.
Finally, some limitations of this study are expected to open multiple avenues for future research. First, this study measures CO 2 emissions through consumption. Future studies are encouraged to examine the impact of FDI and GDP on other types of emissions (e.g., sulfur dioxide emissions (SO 2 ), air quality, or ecological footprint). Second, this study did not verify how renewable and non-renewable energy sources can individually affect São Paulo state's emissions. Third, we did not control the effect of economic structure and the type of FDI on pollution.
Despite these limitations, this study contributes to the literature on growth, foreign investments, and CO 2 emissions. As previously described, there is no consensus on the effects of growth and FDI on the environment, particularly in emerging economies. To the best of our knowledge, this is the first study to apply Brazilian regional FDI data to better comprehend its environmental consequences. To perform such an analysis, a handmade and unique dataset was built, and more novel econometric techniques were employed. Our findings have important policy implications and are arguably capable of guiding future research in developing regions. Funding: This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior-Brasil (CAPES). This study was also supported by the São Paulo Research Foundation (FAPESP) (Process no. 2019/19905-0). The first author acknowledges the financial support from the Brazilian National Council for Scientific and Technological Development (CNPq, 140263/2020-9).
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Data will be made available upon reasonable request, following the embargo period on the FDI-related dataset, which will make the dataset publicly available on 30/11/2023.

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