Acessibilidade / Reportar erro

Evidências sobre curva ambiental de Kuznets e convergência das emissões Os autores agradecem os comentários de um parecerista anônimo.

Resumo

Recentemente, diversos pesquisadores se dedicaram a estudar a relação existente entre crescimento econômico e poluição. Desses estudos surgiram dois enfoques distintos: a Curva Ambiental de Kuznets (CAK) e a convergência das emissões. Neste contexto, Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010. desenvolveram um modelo que faz uma ligação entre estas duas metodologias, à medida que a convergência de emissões implica a CAK. O objetivo deste artigo é estimar a CAK (de um modo tradicional) e o modelo de Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010. com o intuito de verificar se os resultados se reforçam mutuamente. As estimações da curva foram sensíveis ao modelo escolhido. A estimação do modelo apontou evidências de convergência condicional das emissões per capita e da CAK.

Palavras-Chave
Curva Ambiental de Kuznets; Convergência das emissões; Crescimento econômico

Abstract

Recently many researchers started to study the relationship between economic growth and pollution. Two distinct approaches followed from these studies: the Environmental Kuznets Curve (EKC) and the convergence of emissions. In this context, Brock and Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010. built a model that makes a bridge between both methodologies, as convergence of emissions implies the EKC. The purpose of this paper is to estimate the EKC (in a traditional way) and the model proposed by Brock and Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010. in order to verify if the results reinforce each other. The estimates of the curve were sensitive to the model chosen. The estimates of the model pointed to evidence of conditional convergence of per capita emissions and of the EKC.

Keywords
Environmental Kuznets Curve; Convergence of emissions; Economic growth

1. Introdução

A energia ocupa um papel de destaque no desenvolvimento econômico. A partir da Revolução Industrial, as inovações tecnológicas favoreceram o uso crescente de combustíveis fósseis. A consequência direta desse uso é o aumento da concentração de CO2 na atmosfera (ultrapassando 400 ppm no Observatório de Mauna Loa, em 2013). Assim, atribui-se este resultado ao crescimento econômico, visto como o principal responsável pela poluição gerada que provocou danos ao meio ambiente. Discussões a respeito dos conflitos entre crescimento econômico e meio ambiente se iniciaram na década de 1970, quando se defendeu que o meio ambiente, por causa da finitude de recursos naturais, cria limites ao crescimento econômico (Meadows et al., 1972MEADOWS, D. H. et al. The limits to growth. New York: Universe Books, 1972.).

Porém, com o avanço dos estudos referentes ao crescimento econômico e meio ambiente, uma nova visão surgiu. Por essa ótica, seria possível conciliar a riqueza de um país com um ambiente sustentável. Neste contexto, duas metodologias podem ser destacadas: a Curva Ambiental de Kuznets, desenvolvida por Grossman e Krueger (1991)GROSSMAN, G. M.; KRUEGER, A. B. Environmental impacts of a North American free trade agreement. National Bureau of Economic Research, 1991., Shafik e Bandyopadhyay (1992)SHAFIK, N.; BANDYOPADHYAY, S. Economic growth and environmental quality: a time series and cross-country evidence. Journal of Environmental Economics and Management, v. 4, p. 1-24, 1992. e Panayotou (1993)PANAYOTOU, T. Empirical tests and policy analysis of environmental degradation at different stages of economic development. Technology and Employment Programme, 1993.; e a convergência das emissões, proposta por List (1999)LIST, J. A. Have air pollutant emissions converged among U.S. regions? Evidence from unit root tests. Southern Economic Journal, v. 66, n. 1, p. 144-155, 1999..

A Curva Ambiental de Kuznets (CAK deste ponto em diante) postula uma relação entre emissões e renda no formato de U invertido. Em níveis iniciais do crescimento econômico, renda e emissões guardam uma relação direta. Porém, à medida que a renda cresce e chega a um determinado nível, as emissões tendem a apresentar trajetória descendente. Isto ocorre porque, a níveis maiores de renda, as pessoas estão dispostas a buscar crescimento associado a um meio ambiente menos degradado. De acordo com Grossman e Krueger (1995)______. Economic growth and the environment. Quarterly Journal of Economics, v. 110, p. 353-377, 1995., é possível explicar o formato da CAK com base em três efeitos distintos: efeito escala, quando um aumento da renda está associado a um aumento das emissões de poluentes para uma dada composição do produto e uma dada tecnologia; efeito composição, quando uma mudança da composição do produto associado a um aumento da renda leva a uma alteração do peso dos produtos menos poluentes no produto total, conduzindo a uma redução das emissões quando o peso dos produtos menos poluentes aumenta; efeito técnica, quando a adoção de tecnologias menos poluentes associada a um aumento da renda leva a uma redução das emissões. Olhando para a CAK, a relação direta entre renda e emissões é explicada pela predominância do efeito escala em relação aos demais. Esta relação direta ocorre nas fases iniciais do crescimento de um país. A relação inversa entre renda e emissões é explicada pela predominância dos efeitos composição e técnica sobre o efeito escala. Neste último caso, provavelmente o efeito técnica é o mais importante. Esta relação inversa ocorre no caso de países mais avançados. Segundo Grossman e Krueger (1995)______. Economic growth and the environment. Quarterly Journal of Economics, v. 110, p. 353-377, 1995., o nível de renda em que as emissões atingem o máximo é conhecido na literatura como turning point (que pode ser traduzido de forma livre como pico ou ponto de inflexão). O ponto de inflexão pode ocorrer em um nível de renda per capita factível de ser obtido no futuro próximo. Nesse caso, a diminuição das emissões ocorreria de uma forma natural, sem a necessidade de uma intervenção da autoridade ambiental e a introdução de ineficiências dela decorrentes. Caso contrário, pode-se forçar a ocorrência do ponto de inflexão por meio de politicas ambientais apropriadas, conforme sugerido em Dasgupta et al. (2002)DASGUPTA, S.; LAPLANTE, B.; WANG, H.; WHEELER, D. Confronting the Environmental Kuznets Curve. Journal of Economic Perspectives, v. 16, p. 147-168, 2002.. Por exemplo, a adoção de tecnologias menos poluentes pode ser o resultado de políticas ambientais que introduzam os incentivos corretos para a adoção de tecnologias mais limpas.

Os impactos sobre o meio ambiente associados ao crescimento econômico produziram uma série de estudos empíricos. Graças a uma maior disponibilidade de bancos de dados, a CAK passou a ser estimada para diversos tipos de poluentes. Evidências empíricas no formato de U invertido foram obtidas por Grossman e Krueger (1991)GROSSMAN, G. M.; KRUEGER, A. B. Environmental impacts of a North American free trade agreement. National Bureau of Economic Research, 1991., Shafik e Bandyopadhyay (1992)SHAFIK, N.; BANDYOPADHYAY, S. Economic growth and environmental quality: a time series and cross-country evidence. Journal of Environmental Economics and Management, v. 4, p. 1-24, 1992. e Panayotou (1993)PANAYOTOU, T. Empirical tests and policy analysis of environmental degradation at different stages of economic development. Technology and Employment Programme, 1993.. O produto comum desses estudos é a noção de que a degradação ambiental, durante os primeiros estágios do desenvolvimento, intensifica-se e, posteriormente, reduz-se, seguindo uma curva no formato de U invertido. Autores como Selden e Song (1994)SELDEN, T. M.; SONG, D. Environmental quality and development: is there a Kuznets Curve for air pollution emissions? Journal of Environmental Economics and Management, v. 27, n. 2, p. 147-162, 1994., Auci e Becchetti (2006)AUCI, S.; BECCHETTI, L. The instability of the adjusted and unadjusted Environmental Kuznets Curves. Ecological Economics, v. 60, p. 282-298, 2006., Fonseca e Ribeiro (2005)FONSECA, L. N.; RIBEIRO, E. P. Preservação ambiental e crescimento econômico no brasil. In: Encontro de economia da Região Sul, VII. Rio Grande do Sul: [s.n.], 2005., Carvalho e Almeida (2010)CARVALHO, T. S.; ALMEIDA, E. A hipótese da Curva de Kuznets Ambiental Global: uma perspectiva econométrico-espacial. Estudos Econômicos, v. 40, n. 3, p. 587-615, p. 1149-1159, 2010. encontraram evidências favoráveis à existência da CAK.

Porém, como ressaltam Hilton e Levinson (1998)HILTON, F. G. H.; LEVINSON, A. Factoring the Environmental Kuznets Curve: evidence from automotive lead emissions. Journal of Environmental Economics and Management, v. 35, p. 126-141, 1998., através de um exame das evidências de uma CAK para as emissões automotivas de chumbo, temos que a melhoria ambiental que acompanha o crescimento da renda depende dos tipos de política ambiental que induzem o efeito técnica, ligados aqui ao que os autores chamam de intensidade da poluição. É necessário cautela para interpretar a CAK. Segundo De Bruyn (2000)DE BRUYN, S. M. Economic Growth and the Environment. Dordrecht: Kluwer Academic Publishers, 2000., houve uma interpretação errônea de que a implementação de políticas de controle de emissões só seria bem-sucedida quando a renda ultrapassasse certo nível. Além do mais, conforme ressaltado em Grossman e Krueger (1995)______. Economic growth and the environment. Quarterly Journal of Economics, v. 110, p. 353-377, 1995., mesmo que as emissões apresentem uma trajetória na forma de U invertido, isso não garante que este processo de ascensão e queda das emissões se sustente no longo prazo. Logo, o crescimento econômico sozinho pode não levar ao que é postulado nesta teoria. Assim, são necessárias políticas ambientais que, juntamente com o crescimento econômico, venham a induzir a redução da pressão ambiental.

Apesar dos resultados favoráveis, muitas críticas foram feitas à hipótese da CAK. A principal reclamação reside no uso incorreto da técnica econométrica, que pode levar a evidências enganosas. Stern (2004)STERN, D. I. The rise and fall of the environmental Kuznets curve. World Development, v. 32, p. 1419-1439, 2004. argumenta que, sob o rigor de uma econometria correta, a CAK não existe. Além disso, a possível relação obtida entre as variáveis pode ser espúria (Wagner, 2008WAGNER, M. The carbon Kuznets Curve: a cloudy picture emitted by bad econometrics? Resource and Energy Economics, v. 30, p. 388-408, 2008.). Os resultados das estimações são sensíveis ao tipo de poluente considerado. O ponto de inflexão e os próprios coeficientes também se mostram sensíveis à forma funcional e ao tipo de modelo utilizado (De Bruyn, 2000DE BRUYN, S. M. Economic Growth and the Environment. Dordrecht: Kluwer Academic Publishers, 2000.).

O estudo acerca da convergência das emissões de poluentes está inspirado no modelo de Solow (1956)SOLOW, R. A contribution to the theory of economic growth. Quarterly Journal of Economics, v. 70, p. 65-94, 1956., em que a taxa de crescimento da renda per capita é inversamente proporcional ao nível inicial da renda per capita. Uma implicação do resultado de Solow é que os países mais pobres devem alcançar os países mais ricos ao longo do tempo (desde que os parâmetros dos países sejam semelhantes). Na nova vertente que examina a convergência das emissões de poluentes, List (1999)LIST, J. A. Have air pollutant emissions converged among U.S. regions? Evidence from unit root tests. Southern Economic Journal, v. 66, n. 1, p. 144-155, 1999. avalia se a convergência da renda pode ocorrer juntamente com a convergência das emissões de poluentes para regiões dos EUA no período 1929-1994. List (1999)LIST, J. A. Have air pollutant emissions converged among U.S. regions? Evidence from unit root tests. Southern Economic Journal, v. 66, n. 1, p. 144-155, 1999. obtém o resultado de convergência dos indicadores de qualidade ambiental das regiões dos EUA. A explicação desse resultado é que as regiões mais pobres tendem a se especializar em produtos intensivos em poluição, enquanto as regiões mais ricas tendem a se especializar em bens cuja produção utiliza uma tecnologia menos poluente.

Cientes das limitações da CAK tradicional, Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010., desenvolveram um modelo teórico baseado no modelo de Solow que mostra implicitamente a existência da CAK quando ocorre a convergência das emissões. Segundo os autores, pelo fato de a equação a ser estimada ter sido obtida diretamente da teoria, a análise é mais robusta e o modelo não é tão sensível aos dados quanto a análise tradicional da CAK.

Este trabalho propõe-se a fazer a estimação de dois modelos diferentes: a análise tradicional da CAK e a convergência das emissões de Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010..1 1 Embora a estimação da equação de convergência tenha sido feita por Brock e Taylor (2010), neste trabalho será adotada uma técnica mais moderna, com o uso da metodologia de painel dinâmico. Além disso, os dados que serão utilizados neste trabalho são mais atuais do que aqueles utilizados pelos referidos autores. Espera-se, com isso, obter resultados sobre a existência da CAK que se complementem. Essa análise mais completa é útil porque, no modelo de Brock e Taylor, a existência da CAK é implícita, não sendo possível calcular o nível de renda associado ao ponto de inflexão. Logo, é razoável supor que os resultados da convergência das emissões e da CAK sejam consistentes entre si, utilizando-se o mesmo conjunto de países.

2. Modelo teórico

2.1 Curva Ambiental de Kuznets Tradicional

O conceito de Curva Ambiental de Kuznets surgiu no começo da década de 1990 para descrever os impactos que o desenvolvimento econômico exerce sobre os níveis de emissões de poluentes, uma vez que estes estão intimamente relacionados com a produção de um país. O crescimento econômico em um país pobre vem acompanhado por emissões em níveis crescentes associadas ao crescimento da atividade produtiva. Devido à condição de pobreza do país, é dada uma prioridade menor sobre a qualidade ambiental, o que implica menos controle e poucas políticas ambientais sobre os níveis de emissões.

A partir do momento em que a renda cresce e alcança o ponto de inflexão, a prioridade dada pelos indivíduos muda em direção a um ambiente menos degradado. Se esse efeito é forte o suficiente, ocorrerá um declínio nas emissões de poluentes, de acordo com a Figura 1.

Figura 1
Curva Ambiental de Kuznets

Conforme vimos na seção anterior, existem três forças que agem por trás da CAK: efeito escala, efeito composição e efeito técnica. A hipótese teórica em que se baseia a CAK é a de que o impacto negativo do efeito escala tende a prevalecer nos estágios iniciais do crescimento econômico, sendo compensado ao longo do processo de desenvolvimento do país pelos efeitos composição e técnica, os quais tendem a reduzir os níveis das emissões.

Porém, conforme ressalta Diniz (2001)DINIZ, E. M. Crescimento, poluição e o Protocolo de Quioto: uma avaliação do caso brasileiro. São Paulo: Banco Santos / Universidade de Oxford, 2001., se o efeito escala predomina sobre os demais não ocorrerá a segunda fase da CAK, uma vez que um nível maior de produto sempre produz um nível maior de poluição. Segundo Diniz (2007)______. Lessons from the Kyoto protocol. Ambiente e Sociedade, v. 10, n. 1, p. 27-38, 2007., o efeito escala é predominante nas nações menos desenvolvidas, uma vez que o crescimento econômico é o objetivo mais desejável a priori. Já nos países desenvolvidos que estão na parte negativa da curva, os efeitos composição e técnica predominam sobre o efeito escala.

A discussão a respeito da existência da CAK se desenvolveu principalmente no plano econométrico. Esta relação, por ser empírica, fundamenta-se na seguinte expressão:

em que o subscrito t se refere ao tempo, α é uma constante, Eit representa o logaritmo das emissões e Yit o logaritmo da renda per capita. A curva em formato de U invertido surge quando β1>0, β2<0, e β3= 0. Neste caso o ponto de inflexão é dado por .

Questionamentos foram feitos de que a CAK seria um artefato empírico sem o embasamento de uma teoria econômica. A ocorrência de uma relação inversa entre poluição e renda foi apontada em diversos modelos de crescimento econômico, como, por exemplo, os relatados em Stokey (1998)STOKEY, N. L. Are there limits to growth? International Economic Review, v. 39, n. 1, p. 1-31, 1998. e Dinda (2005)DINDA, S. A theoretical basis for the Environmental Kuznets Curve. Ecological Economics, v. 53, p. 403 - 413, 2005.. Stokey (1998)STOKEY, N. L. Are there limits to growth? International Economic Review, v. 39, n. 1, p. 1-31, 1998. examinou os principais modelos de crescimento econômico e obteve uma condição única para o aparecimento de uma CAK e que é válida para todos os modelos. Logo, a CAK pode ocorrer em modelos de crescimento econômico microfundamentados. Segundo Stokey (1998)STOKEY, N. L. Are there limits to growth? International Economic Review, v. 39, n. 1, p. 1-31, 1998., um método mais poluente de produção é utilizado se a capacidade produtiva estiver abaixo do ponto de inflexão. À medida que a economia cresce, o ponto de inflexão é ultrapassado e as tecnologias menos poluentes podem ser usadas. Nesse contexto, abaixo do ponto de inflexão a poluição cresce linearmente com a renda, porém, acima dele, a relação entre emissões e renda depende da elasticidade de substituição intertemporal do consumo. Se esta for menor do que um, encontra-se a relação no formato de U invertido. Através do exame de dados empíricos em Campbell e Mankiw (1989)CAMPBELL, J. Y.; MANKIW, N. G. Consumption, income, and interest rates: reinterpreting the time series evidence. NBER Macroeconomics Annual, v. 4, p. 185-216, 1989., temos que a condição de Stokey (1998)STOKEY, N. L. Are there limits to growth? International Economic Review, v. 39, n. 1, p. 1-31, 1998. é satisfeita, de modo que pode ocorrer uma CAK e a análise empírica precisa ser feita caso a caso. Logo, a CAK pode surgir pelo próprio processo do crescimento econômico e, assim, sua existência torna-se apenas uma questão empírica.

Por outro lado, Dinda (2005)DINDA, S. A theoretical basis for the Environmental Kuznets Curve. Ecological Economics, v. 53, p. 403 - 413, 2005. mostra a importância de políticas de controle e redução das emissões no contexto da CAK. O autor utilizou um modelo de crescimento endógeno em que uma parte do capital é utilizada para a produção de mercadorias, as quais geram poluição. Porém, outra parte do capital é usada para práticas mitigadoras no custo da poluição. Desta forma, a mudança de uma alocação insuficiente de capital para uma alocação suficiente do capital destinado à mitigação é a base para a relação em que se sustenta a CAK no modelo apresentado.

2.2 Modelo de convergência das emissões

Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010. desenvolveram um modelo para emissões de poluentes, intimamente ligado ao modelo de Solow. No modelo proposto, o progresso tecnológico sobre a produção de bens e redução das emissões é exógeno, o que leva a um crescimento contínuo da economia e também a uma melhoria na qualidade ambiental. Além disso, tanto a poupança quanto as práticas de mitigação escolhidas são, também, determinadas exogenamente.2 2 Uma discussão sobre o problema de decisão acerca da adoção de práticas de mitigação pode ser encontrada, por exemplo, em Baumol e Oates (1988) e Perman et al. (2011). As medidas para mitigar a poluição englobam impostos de Pigou e permissões comercializáveis de emissões. O imposto de Pigou consiste em um imposto sobre as emissões de poluentes que procura internalizar os custos e fazer com que a firma venha a considerar as externalidades em seu processo interno de decisão de produção. A alíquota do imposto de Pigou é igual ao dano social marginal. As permissões comercializáveis de emissões são uma tentativa de obter o mesmo resultado do imposto de Pigou com o auxílio do mercado. A autoridade ambiental determina a quantidade agregada de emissões no montante em que o dano social marginal é igual ao custo marginal de mitigação. A seguir, emite permissões comercializáveis de emissões nesse montante e deixa a alocação das emissões entre fontes para o mercado. Considerando as curvas de dano social marginal e custo marginal de mitigação em um ambiente de certeza, a alíquota do imposto de Pigou é aquela que produz o montante permitido de emissões, de modo que o imposto de Pigou e as permissões comercializáveis de emissões produzem exatamente o mesmo resultado em termos práticos. Em um contexto de incerteza, a equivalência não permanece, e as duas práticas de mitigação podem produzir diferenças substanciais nos valores esperados do bem-estar social. A autoridade ambiental irá preferir, em um contexto de incerteza, uma política ou outra, dependendo da incerteza afetar o dano social marginal, o custo marginal de mitigação ou ambos. Para maiores detalhes, consultar Baumol e Oates (1988), cap. 5 e Perman et al. (2011), cap. 6. Desta forma têm-se as seguintes relações:

O produto Y possui uma função de produção estritamente côncava, com retornos constantes de escala, e utiliza os insumos trabalho efetivo e capital. O capital se acumula a uma taxa de poupança fixa s e deprecia a uma taxa δ. A taxa do progresso tecnológico é dada por gB. A força de trabalho cresce à taxa n.

A introdução das emissões no modelo é dada por Ω, em que cada unidade F da atividade econômica gera Ω unidades de poluição. Porém, a quantidade de poluição produzida vai depender da existência de práticas mitigadoras nas emissões de poluentes. Supõe-se que a prática mitigadora seja uma função com retornos constantes de escala e que a quantidade de poluição mitigada seja uma função estritamente côncava da atividade econômica F e dos esforços de abatimento F A. Se a prática de mitigação no nível ΩA remove unidades de poluição do total criado, então a poluição emitida é igual a poluição gerada menos a poluição mitigada. Assim,

Desta forma, combinando as hipóteses sobre a poluição e redução das emissões com o modelo de Solow, o produto Y disponível para consumo ou investimento torna-se Y=[1 – θ]F. Neste caso θ representa a fração da atividade econômica destinada à mitigação. Supõe-se, ainda, que a taxa de progresso tecnológico seja exógena em relação à mitigação das emissões, o que faz com que Ω se reduza à taxa gA>0.

As taxas de crescimento devem ser iguais ao longo da trajetória de crescimento equilibrado. À medida que a economia se encontra na trajetória de crescimento equilibrado, o produto agregado, consumo e capital crescem à taxa gB + n, enquanto suas versões por trabalhador crescem à taxa gB, isto é, gy = gk = gc = gB > 0.

E quando a economia se encontra em seu estado estacionário, a taxa de crescimento das emissões agregadas é dada por gE = gB + n - gA. A taxa de crescimento das emissões depende de um efeito escala (gB + n), e de um efeito técnica criado pelo impacto do progresso tecnológico no abatimento das emissões (gA). O crescimento sustentável ocorre quando a taxa de crescimento da emissões agregadas no equilíbrio é negativa (gE< 0).

O modelo de Solow verde relaciona o nível de renda per capita e a qualidade ambiental produzindo uma trajetória conhecida como Curva Ambiental de Kuznets. Isso pode ser verificado por meio das taxas de crescimento das emissões agregadas e de variação do capital por trabalhador efetivo, dadas por:

As duas expressões acima valem para qualquer ponto do tempo. Por meio delas, pode-se ligar a dinâmica da acumulação do capital à evolução dos níveis de poluição. A curva gerada desta relação é negativamente inclinada, e pode ser considerada um locus da poupança, uma vez que se desloca quando ocorre uma variação na taxa de poupança s. À medida que o capital por trabalhador efetivo aumenta, o locus da poupança se aproxima de zero.

A análise gráfica (Figura 2) permite mostrar que a taxa de variação do capital por trabalhador efetivo é positiva em pontos à esquerda de B, e negativa em pontos à direita. A intersecção no ponto B gera o estado estacionário do capital por trabalhador efetivo k*. Níveis de mais baixos e mais distantes do equilíbrio implicam um crescimento mais acelerado. Quando a economia entra em uma trajetória de crescimento equilibrado, a taxa de variação do capital por trabalhador efetivo é zero. Nesse ponto, o produto agregado da economia e o capital crescem à taxa gB + n. Com relação às emissões, a taxa de crescimento é zero em T, positiva à esquerda e negativa à direita. Dessa forma, o ponto T pode ser associado ao ponto de inflexão das emissões. Sob a hipótese de crescimento sustentável, o modelo gera a Curva Ambiental de Kuznets retratada no segundo gráfico da Figura 2.

Figura 2
Trajetória das emissões de poluentes

Se o estoque inicial de capital em uma economia for baixo, então as emissões, em um primeiro momento, crescerão e, posteriormente, diminuirão. Se o estoque inicial de capital for alto, pode ocorrer uma queda monotônica das emissões à medida que a economia se move em direção ao equilíbrio.

Na especificação longa da equação, supõe-se que as taxas de progresso tecnológico são comuns a todos os países, e o conjunto de fatores que determinam o estado estacionário e a taxa de crescimento ao longo da trajetória de crescimento equilibrado varia entre os países. Assim, permitindo que si, θi e (δ + n + gB)i e sejam uma média temporal específica para cada país, a seguinte equação pode ser obtida:

em que o termo à esquerda é a taxa de crescimento das emissões em cada período, sit é a poupança, θi é a taxa de mitigação e (δ + n + gB)it é a taxa de depreciação efetiva.

De acordo com Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010., poupança e mitigação afetam a trajetória de crescimento das emissões e do produto e mudam as magnitudes do estado estacionário. Porém, não exercem impacto algum sobre a taxa de crescimento de longo prazo. Por sua vez, o crescimento tecnológico e o crescimento populacional, além de exercer impactos sobre a trajetória de crescimento das emissões e do produto e sobre os valores de equilíbrio, também alteram as taxas de crescimento de longo prazo.

3. Metodologia e dados

Neste trabalho serão utilizadas duas metodologias diferentes para estimação de dados em painel, uma vez que os modelos usam dimensões temporais diferentes. De acordo com Levin et al. (2002)LEVIN, A.; LIN, C. F.; CHU, C. S. J. Unit roots tests in panel data: asymptotic finite sample properties. Journal of Econometrics, v. 108, p. 1-24, 2002., para painéis de tamanho médio (25 < T < 250, 10 < N < 250), em que T se refere a dimensão temporal e N corresponde ao número de indivíduos, é recomendado o uso de testes de raiz unitária e demais procedimentos de séries temporais específicos para metodologia de painel. No entanto, para painéis em que T é pequeno (considerando T < 25) e N é grande (10 < N < 250) pode-se aplicar a metodologia padrão de painel e não é necessário utilizar os procedimentos de séries temporais. Para o caso da CAK padrão, em que N = 99, e T = 37, será necessário o uso da metodologia de séries temporais. Já para modelo de Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010., como N = 99 e T = 6, não será necessário o uso das técnicas de séries temporais. Segue uma breve descrição dos procedimentos que serão aplicados neste trabalho.

3.1 Testes de raiz unitária 3 3 Maiores detalhes em Breitung e Pesaram (2005).

Os testes de raiz unitária desenvolvidos para dados em painel, podem ser divididos em duas gerações de testes. A primeira geração assume a hipótese de independência em corte transversal, já a segunda geração permite que ocorra correlação entre estas unidades.

A ideia básica dos testes de raiz unitária em painel, assim como na análise convencional de séries temporais, é encontrar o valor do coeficiente autorregressivo de um modelo AR(1) como o descrito abaixo:

em que i = 1,...,N, t = 1,...,T, yit é a variável a ser testada e εit é um termo de erro estacionário. O termo xit representa os componentes estocásticos e não estocásticos. O modelo (8) pode ser reescrito como uma regressão do tipo Dickey-Fuller Aumentado (ADF), da seguinte forma:

em que Δyit = yit - yit-1, φi = 1 - ρi e as p defasagens variam para cada setor.

A hipótese nula é a mesma para todos os testes:4 4 Com exceção do teste Hadri (2000).

Ou seja, a hipótese nula é a de que todas as unidades do painel possuem raiz unitária. A peculiaridade dos testes de raiz unitária em painel é a presença de duas hipóteses alternativas diferentes:

Sob a hipótese alternativa H1a, o componente autorregressivo é o mesmo para todas as unidades do painel, por isso é chamada de hipótese alternativa homogênea. Os testes de Levin, Lin e Chu (2002)LEVIN, A.; LIN, C. F.; CHU, C. S. J. Unit roots tests in panel data: asymptotic finite sample properties. Journal of Econometrics, v. 108, p. 1-24, 2002. e Breitung (2000) adotam este tipo de hipótese. A alternativa H1b, por sua vez, é chamada de alternativa heterogênea, já que não impõe um coeficiente autorregressivo específico para cada unidade de corte transversal. Os testes IPS, de Im, Pesaram e Shin (2003), Fisher ADF e Fisher PP, de Madalla e Wu (1999), Choi ADF e Choi PP, de Choi (2001), e CIPS, de Pesaram (2003), adotam este tipo de hipótese alternativa.

3.2 Cointegração

Uma vez verificado se as séries são estacionárias ou não, procede-se aos testes de cointegração em painel, a fim de verificar se a relação entre as séries não é espúria. Pedroni (1999)PEDRONI, P. Critical values for cointegration tests in heterogeneous panels with multiple regressors. Oxford Bulletin of Economics and Statistics, v. 61, p. 653-670, 1999. desenvolveu testes estatísticos seguindo a metodologia proposta por Engle e Granger (1987). O autor considerou a seguinte regressão:

para t = 1,...,T; i = 1,...,N; m = 1,...,M.

Neste caso, T se refere ao número de observações no tempo, N se refere ao número de indivíduos no painel e M se refere ao número de regressores incluídos na equação. Supõe-se que as variáveis yit e xi,t sejam integradas de ordem um (I(1)). Assim, como na análise de cointegração em séries temporais, o estabelecimento de relações de cointegração entre as séries é equivalente a mostrar que o termo de erro na Equação (13) é estacionário.

A hipótese nula associada ao procedimento do teste de Pedroni é que o resíduo possui raiz unitária, que é equivalente a supor que não existem relações de cointegração. Aqui também há duas hipóteses alternativas diferentes. A hipótese alternativa homogênea considera um mesmo coeficiente autorregressivo para o coeficiente ρ. Já a hipótese alternativa heterogênea não impõe que o coeficiente tem que ser o mesmo para todas as unidades de corte transversal. Uma limitação do teste de Pedroni é a hipótese de independência das unidades do painel. Assim como nos testes de raiz unitária em painel, o teste de Pedroni pode apresentar distorções de tamanho e poder na presença deste tipo de correlação.

A fim de evitar este tipo de problema, Westerlund (2007)WESTERLUND, J. Testing for error correction in panel data. Oxford Bulletin of Economics and Statistics, v. 69, p. 709-748, 2007. desenvolveu um teste de cointegração baseado no método de correção dos erros, no qual a presença de dependência não afeta o tamanho e o poder do teste. O autor considerou o seguinte processo:

em que t = 1,...,T e i = 1,...,N. O escalar yit contém tanto termos deterministas (Φ1i + Φ2it) quanto um termo estocástico Zit.

A construção das estatísticas é feita sob a seguinte equação:

em que dt = (1,t)',αi, é o parâmetro de correção dos erros e λi = -αiβi. O parâmetro αi não é afetado pela imposição de βi, o que sugere a validade do teste sobre αi para verificar a cointegração das séries. O teste faz uso de quatro estatísticas: duas desenvolvidas na versão pooled e duas construídas na versão group mean. A hipótese alternativa difere para as duas versões de estatísticas.

Neste trabalho, optou-se por fazer dois testes de cointegração. Um baseado na primeira geração de testes e outro baseado na segunda geração. Espera-se, com isso, evitar que distorções no teste de Pedroni afetem a robustez das estimações.

3.3 Estimação do vetor de cointegração

Engle e Granger (1987)5 5 Engle, R.F.; Granger, C.W.J. Cointegration and error correction: representation, estimation, and testing. Econometrica, v.55, n.2, pp.251-76, 1987. argumentam que se uma combinação linear de duas ou mais séries I(1) for estacionária, pode-se dizer que as séries cointegram, ou seja, apresentam uma relação de longo prazo. Neste trabalho serão utilizados estimadores de equações simples, a fim de se estimar a CAK.

Para uma melhor definição dos estimadores que serão utilizados, considere o seguinte modelo simples de efeitos fixos em painel:

em que {yit} é uma sequência de escalares, β é um vetor (K X 1) de parâmetros e {αi} são os interceptos individuais de cada unidade de corte transversal. Se o termo de erro {uit} for estacionário, então a cointegração de yit e xit prevalece para todos os indivíduos do painel. Neste caso, a Expressão (17) descreve um sistema de regressões cointegradas, em que {xit} é um processo integrado de ordem 1 para todo i:

em que εit é i.i.d com média zero e variância constante.

Breitung e Pesaram (2005)BREITUNG, J.; PESARAM, M. H. Unit roots and cointegration in panels. Frankfurt: Deutsche Bundesbank, 2005. argumentam que, mesmo sob as fortes hipóteses que garantem a consistência, o estimador MQO se torna ineficiente, devido à endogeneidade dos regressores e à correlação serial dos resíduos.

A fim de corrigir o problema de endogeneidade dos regressores e correlação serial dos resíduos6 6 Decorrente da estimação da Equação (25) via Mínimos Quadrados Ordinários (MQO). e fornecer estimativas consistentes e eficientes do vetor de cointegração, foram desenvolvidos estimadores alternativos. Entre eles, pode-se citar os estimadores Fully-Modified OLS (FM-OLS) e Dynamic OLS (DOLS), que serão utilizados neste trabalho. O estimador FM-OLS, proposto por Pedroni (2000)______. Fully modified ols for heterogeneous cointegrated panels. Advances in Econometrics, v. 15, p. 93-130, 2000., ajusta os efeitos da endogeneidade dos regressores e da dinâmica de curto prazo dos resíduos. Para corrigir o efeito da endogeneidade de longo prazo dos regressores, a variável dependente é ajustada pela parte do erro que é correlacionada com o regressor.

Outra forma de correção da endogeneidade e correlação serial dos resíduos é feita por meio do estimador DOLS. A versão group-means do estimador foi proposta por Pedroni (2000)______. Fully modified ols for heterogeneous cointegrated panels. Advances in Econometrics, v. 15, p. 93-130, 2000.. A correção paramétrica feita pelo estimador ocorre pela inclusão de p valores defasados e p valores precedentes de Δxit, como regressores adicionais.

3.4 Painel dinâmico

Considere a seguinte especificação de um modelo de painel dinâmico:

em que yit representa a variável dependente, yi,t-1 é a variável dependente defasada em um período e xit são os regressores adicionais. Finalmente, αi e ηt são termos de efeitos individual e temporal, e νit é um termo de erro.

Em um painel dinâmico, a estimação, via OLS, dos coeficientes γ e β é enviesada e inconsistente. Isto ocorre pelo fato de que o regressor yi,t-1 é correlacionado com αi e possivelmente com o termo de erro composto (αi + νit). Mesmo em um painel com efeitos fixos, a estimação é enviesada, uma vez que a hipótese de exogeneidade não pode ser verificada (Cameron e Trivedi, 2005CAMERON, A. C.; TRIVEDI, P. K. Microeconometrics: Methods and Applications. New York: Cambridge University Press, 2005.).

Para contornar esse problema, diversos autores sugerem o uso de variáveis instrumentais, a fim de lidar com o problema da endogeneidade causada pelo regressor defasado. Os estimadores de painel dinâmico utilizam um modelo transformado em primeiras diferenças dado por:

Anderson e Hsiao (1981)7 7 ApudIslam (1995). sugerem o uso de yi,t-2 como instrumento para (yi,t-1 - yi,t-2). Segundo os autores, este é um instrumento válido, já que não é correlacionado com (νit - νit-1), desde que os resíduos não exibam correlação serial. Segundo Islam (2003)______. What have we learnt from the convergence debate? Journal of Economic Surveys, v. 17, n. 3, p. 309-362, 2003. a desvantagem deste método é a perda de duas unidades de corte transversal. Outra forma alternativa e mais eficiente, segundo Anderson e Hsiao (1981), é utilizar Δyi,t-2 como instrumento para Δyi,t-1. Já Arellano e Bond (1991)8 8 ApudIslam (2003). propõem o uso de todas as defasagens da variável dependente como instrumento no modelo de painel dinâmico, por meio de um GMM em primeira diferença.

Porém, Arellano e Bover (1995)ARELLANO, M.; BOVER, O. Another look at the instrumental variable estimation of error components models. Journal of Econometrics, v. 68, n. 1, p. 29-51, 1995. e Blundell e Bond (1998)BLUNDELL, R.; BOND, S. Initial conditions and moment restrictions in dynamic panel data models. Journal of Econometrics, v. 87, p. 115-143, 1998. argumentam que, apesar de eficientes, os instrumentos utilizados pelos modelos anteriores são fracos, principalmente se a variável apresentar forte persistência. Uma alternativa é utilizar um sistema que combina equações em nível e equações em primeira diferença (System-GMM). Os instrumentos para a equação em primeira diferença são os mesmos descritos anteriormente e os instrumentos para a equação em nível são as diferenças defasadas das respectivas variáveis.

O estimador obtido depende da ausência de correlação serial de segunda ordem nos resíduos. Então se aplica um teste de autocorrelação de primeira ordem e outro de segunda ordem nos resíduos em primeira diferença. E para testar a validade dos instrumentos aplica-se o teste de Hansen e o teste de Sargan.

3.5 Dados

Para se fazer a análise pretendida neste trabalho, os dados anuais obtidos correspondem ao período 1972-2008. A amostra contém 99 países que foram escolhidos em função da disponibilidade dos dados. Para o cálculo do modelo de convergência foram considerados seis9 9 Conforme Islam (2003), a frequência estimada em vários estudos de painel é definida por um intervalo de cinco anos, porém, em função da quantidade de dados, optou-se por definir um intervalo de seis anos a fim de incorporar na análise o maior intervalo temporal possível entre a primeira e a última observação. pontos de corte transversal, referentes aos anos 1978, 1984, 1990, 1996, 2002 e 2008.

  • Emissões per capita de CO2 (eit): foram obtidas por meio da International Energy Agency (IEA), Medidas em toneladas métricas. As estimações de CO2 calculadas pela IEA representam o total das emissões de combustíveis fósseis.10 10 O banco de dados CAIT (Climate Analysis Indicators Tool), elaborado pelo World Resources Institute (WRI), poderia ser usado alternativamente para essa análise. Ele possui a vantagem de estar disponível publicamente pela internet no sítio http://cait2.wri.org/wri, ao contrário da base da IEA. Ambas as bases possuem alta confiabilidade. Neste artigo optamos por usar o banco de dados do IEA.

  • Renda per capita (yit): é medida pelo PIB real per capita, corrigido pela paridade do poder de compra. A variável foi coletada na Penn World Table, versão 7.0, divulgada em 2011.

  • Taxa média de crescimento das emissões per capitaeit): foi construída seguindo a especificação proposta por Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010.. Nesta especificação, obtém-se uma aproximação da taxa média de crescimento exata.

  • Poupança per capita (sit): foi medida pela taxa de investimento disponível na Penn World Table, versão 7.0. O uso desta taxa pode ser justificado pela hipótese utilizada no modelo de Solow de que a economia é fechada. A mesma taxa é utilizada por Barro e Sala-i-Martin (2004, p. 519)BARRO, R. J.; SALA-i-MARTIN. Economic Growth - second edition. Cambridge, MA: MIT Press, 2004..

  • Taxa de depreciação efetiva (δ + n + gB)it: também foi construída seguindo Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010.. A taxa de crescimento populacional foi obtida junto ao Banco Mundial. Segundo Romer (2001)ROMER, D. Advanced Macroeconometrics. 2. ed. Nova Iorque: McGraw-Hill, 2001., a taxa de depreciação (δ) varia entre 1% e 2%. Já a taxa de crescimento tecnológico (gB) varia entre 3% e 4 %. Assim, Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010. consideraram (δ + gB) equivalente a 5%. Desta forma, a taxa de depreciação efetiva será composta da média do crescimento populacional de cada período considerado, acrescido de 0,05 para cada média calculada.

A seguir serão reportadas informações a respeito da base de dados. A Tabela 1 mostra a estatística descritiva das variáveis e o quadro 1 mostra quais foram os países utilizados na amostra.

Tabela 1
Estatística Descritiva das Variáveis
Quadro 1
Países constantes da amostra

4. Resultados

A presente seção pretende mostrar os resultados do trabalho. Uma vez que foram utilizados dois modelos diferentes, CAK e modelo de Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010., foi necessário o uso de duas metodologias distintas: séries temporais em painel e modelos de painel curto. Por isso, optou-se por apresentar os resultados divididos em duas seções. A primeira seção traz todos os resultados e estimações para a CAK. Já a segunda seção traz os resultados do modelo de Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010..

4.1 Curva Ambiental de Kuznets

Quando se trata da análise de cointegração, o primeiro passo é verificar se as séries são não estacionárias por meio de testes de raiz unitária. Os testes aqui reportados adotam a hipótese alternativa heterogênea com a especificação mais geral (tendência e intercepto).11 11 Foram feitos testes nas demais especificações (intercepto e passeio aleatório), que apontaram para a presença de raiz unitária. Sob esta hipótese, a rejeição de H0, implica que parte das séries é estacionária.12 12 Pois H0: as séries possuem raiz unitária. Desta forma, a Tabela 2 reporta os resultados dos seguintes testes: IPS, Fisher ADF, Fisher PP, Choi ADF, Choi PP, representando a primeira geração de testes. Já o teste CIPS é o teste representativo da segunda geração.

Tabela 2
Testes de Raiz Unitária: intercepto e tendência

De modo geral, os testes apontam para a presença de raiz unitária nas séries. A exceção fica por conta dos testes Fisher ADF e Fisher PP, que apontaram para a estacionariedade da variável emissões. No entanto, esses testes pressupõem independência entre as unidades do painel e, conforme ressalta Baltagi (2005)BALTAGI, B. H. Econometric Analysis of Panel Data. 3. ed. Inglaterra: John Wiley & Sons, 2005., apresentam distorções de tamanho quando N = 100. E é exatamente o que se verifica para amostra em questão, uma vez que N = 99. Além do mais, o teste CIPS, que relaxa a suposição de independência, não rejeitou a hipótese de raiz unitária. Desta forma, para a especificação geral (intercepto e tendência), os testes apontaram para a presença de raiz unitária nas séries.

Uma vez verificado que as séries apresentam raiz unitária e são integradas de mesma ordem,13 13 Foram feitos os mesmos testes sobre a primeira diferença das variáveis, porém optou-se por não reportá-los neste trabalho. Os testes apontaram que todas as variáveis são estacionárias na primeira diferença. parte-se agora para os testes de cointegração. O teste de Pedroni apresenta três estatísticas de painel e três estatísticas de grupo.14 14 Optou-se por não reportar a estatística v, uma vez que os resultados foram próximos aos da estatística rho, e também porque podem existir distorções de tamanho e poder devido à dimensão das séries utilizadas neste trabalho. A hipótese nula do teste é ausência de cointegração. Logo, verifica-se cointegração entre as variáveis quando se rejeita a hipótese nula. A equação de cointegração foi especificada admitindo-se apenas uma constante.

Daqui em diante serão considerados dois modelos: um modelo com os termos renda e renda ao quadrado; e um modelo com os termos renda, renda ao quadrado e renda ao cubo.

Tabela 3
Teste de cointegração de Pedroni

O teste de cointegração de Pedroni apresentou alguns resultados aparentemente contraditórios. Enquanto a estatística rho não foi estatisticamente significante, as demais estatísticas apontam para rejeição da hipótese nula, pelo menos a 10% de significância para as estatísticas de painel. Já utilizando as estatísticas de Grupo, rejeitou-se a hipótese nula a 1%, com exceção novamente da estatística rho.

Pedroni (2004)______. Panel Cointegration: Asymptotic and Finite Sample Properties of Pooled Time Series Tests with an Application to the PPP Hypothesis. Econometric Theory, v. 20, p. 597-625, 2004. argumenta que a estatística rho tende a ter bom poder para painéis com N pequeno e T relativamente grande. Desta forma, o autor adverte que este é um caso extremo, e que em situações diferentes desta, as demais estatísticas tendem a se adequar melhor. Além do mais, os testes assumem independência entre as unidades do painel, caso isso não ocorra os testes podem ter problemas de tamanho e poder.

Uma vez que podem existir problemas de distorção e tamanho no teste de Pedroni em virtude da suposição de independência das unidades do painel, optou-se por realizar outro teste de cointegração. O teste de Westerlund (2007)WESTERLUND, J. Testing for error correction in panel data. Oxford Bulletin of Economics and Statistics, v. 69, p. 709-748, 2007. permite que se faça a correção deste problema usando estatísticas robustas, obtidas por meio do método de bootstrap.15 15 Foram feitas 1.000 replicações.

Tabela 4
Testes de cointegração de Westerlund

Os resultados indicam a existência de cointegração a 1% de significância. Uma vez que o teste permite a correlação entre os indivíduos do painel, os resultados são mais robustos do que aqueles reportados pelo teste de Pedroni. De qualquer forma, ambos os testes apresentaram indícios de que as variáveis são cointegradas.

Para avaliar a ocorrência da relação no formato de U invertido, foram utilizados os estimadores FMOLS e DOLS, ambos na versão group-mean. De acordo com Tabela 5, pode-se verificar que todos os coeficientes estimados foram significantes e apresentaram os sinais esperados pela teoria.

Tabela 5
Estimação do vetor de cointegração

O procedimento de análise da estimação de uma CAK compreende a análise dos sinais das variáveis explicativas e da magnitude do coeficiente da renda ao cubo. Segundo Grossman e Krueger (1995)______. Economic growth and the environment. Quarterly Journal of Economics, v. 110, p. 353-377, 1995. nenhuma inferência pode ser feita a respeito da relação entre os valores dos coeficientes das variáveis explicativas com a variável dependente.

Em ambos modelos, os sinais dos coeficientes são os esperados pela teoria. No modelo 2, que inclui a renda ao cubo como variável explicativa adicional, temos que o coeficiente dessa variável está bem próximo de zero, mas significativo, o que aponta para uma curva do tipo N com a segunda inflexão bem suave, quando a equação é estimada pelo método FMOLS. Porém, quando a equação é estimada pelo método DOLS, o valor é ainda mais próximo de zero.

Já os pontos de inflexão apresentaram uma grande diferença entre os dois modelos. No modelo 1, a renda per capita necessária para que as emissões passem a apresentar uma trajetória declinante é de US$ 792.805,6 quando estimado pelo método FMOLS. Já a estimação por meio do método DOLS apresentou um ponto de inflexão ainda mais elevado: US$ 2.513.755,00. Ambos os valores são extremamente elevados e muito distantes da realidade. Porém, o ponto de inflexão do modelo 2 é mais razoável do que o anterior. O nível de renda per capita de U$S 6.168,88 não é tão elevado quando se confronta com os resultados de Cole et al. (1997) que foi de U$S 25.100,00. Porém, foi superior aos resultados de Shafik e Bandyopadhyay (1992)SHAFIK, N.; BANDYOPADHYAY, S. Economic growth and environmental quality: a time series and cross-country evidence. Journal of Environmental Economics and Management, v. 4, p. 1-24, 1992., que oscilou entre US$ 3.000,00 e U$S 4.000,00. Mesmo quando se considera o ponto de inflexão obtido por meio do método DOLS, um valor US$ 57.837,04 é razoável quando se compara com resultados obtidos anteriormente. Por exemplo, Wagner (2008)WAGNER, M. The carbon Kuznets Curve: a cloudy picture emitted by bad econometrics? Resource and Energy Economics, v. 30, p. 388-408, 2008. encontrou diversos pontos de inflexão de acordo com o método utilizado. Os pontos de inflexão obtidos variaram entre US$ 13.090,00 e US$ 394.449,00.

Desta forma, a ocorrência do movimento descrito pela CAK está condicionada ao nível de renda. Ou seja, de acordo com a teoria, as emissões cairão a partir do momento em que a renda alcançar o ponto de inflexão.

O objetivo desta estimação é apenas confrontá-la com os resultados obtidos do modelo de Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010.. Por isso, foram utilizados apenas dois tipos de estimadores e um modelo tradicional, cuja escolha dos estimadores estava condicionada às equações a serem testadas. Como o foco do trabalho é estimar um modelo alternativo à CAK, a discussão a respeito desta metodologia não será levada adiante.

4.2 Modelo de Convergência das emissões (Brock e Taylor, 2010BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010.)

Preliminarmente, foram feitos os testes de poolabilidade, o teste de Hausman e correlação serial. Quanto ao teste de poolabilidade,16 16 Poolabilidade pode ser entendida como a reunião ou empilhamento dos dados. Sendo assim, consiste em impor um mesmo parâmetro para todas as unidades do painel ao longo do tempo. segundo Baltagi (2005)BALTAGI, B. H. Econometric Analysis of Panel Data. 3. ed. Inglaterra: John Wiley & Sons, 2005., o procedimento do teste consiste em comparar um modelo de MQO com um modelo que utiliza a metodologia de painel.17 17 No teste feito neste trabalho, a hipótese alternativa foi um modelo de Efeitos Fixos. A rejeição da hipótese nula implica a necessidade de se aplicar métodos de painel na estimação do modelo a ser estudado. Já o teste de Hausman consiste na escolha entre dois estimadores: Efeitos Fixos e Efeitos Aleatórios. Sob H0 ambos estimadores são consistentes e eficientes. Já sob a hipótese alternativa, apenas o estimador de Efeitos Fixos é consistente.

Conforme a Tabela 6, os resultados dos testes de poolabilidade apontaram para o uso da metodologia de painel por meio do estimador de efeitos fixos, já que o teste de Hausman apontou para a rejeição da hipótese nula. Quanto à correlação serial, pode-se verificar que ela ocorre no modelo de convergência absoluta. Por isso, foram estimadas regressões robustas à heterocedasticidade e correlação serial pelo procedimento de Windmeijer (2005). Na estimação do modelo de convergência com efeitos fixos foram usadas as Equações (23) para convergência absoluta e (24) para convergência condicional:18 18 Todas as variáveis estão expressas em logaritmos, exceto Δeit.

Tabela 6
Testes Preliminares

Para calcular o modelo por meio do método System-GMM, as Equações (22) e (23) foram reescritas da seguinte forma:19 19 Conforme Islam (2003).

em que γ = 1 + β1 e a convergência ocorre quando γ < 1, pois para que ocorra a convergência é necessário que β1 seja negativo.20 20 Reescrevendo de outra forma: β1 =γ – 1. Neste trabalho foram reportados os valores do coeficiente β1.

De acordo com a Tabela 7, as estimações apontaram para a existência de convergência. Embora, a estimação da convergência absoluta tenha sido estatisticamente significante, prevalece a convergência condicional, uma vez que a adição de variáveis explicativas com significância estatística é um indício a favor da convergência condicional.

Tabela 7
Estimação Modelo de Convergência

Nas estimações, constatou-se a relação positiva que a poupança possui com o nível das emissões, assim como prevê o modelo de Solow verde. Além disso, a relação negativa da taxa de depreciação efetiva apontada no modelo teórico se manteve. E a convergência, apontada pelo sinal negativo do termo defasado das emissões, também se verificou. A taxa de convergência implícita λ pode ser considerada uma taxa relativamente baixa, enquanto a fração da renda que remunera o capital α é mais elevada no caso de estimadores fixos e convergência absoluta. Nos demais casos, a estimativa de α é consistente com aquela normalmente encontrada na literatura (ao redor de 1/3) ou com a estimativa referente a um conceito mais amplo de capital (ao redor de 2/3).

Com relação à estimação utilizando o método System-GMM, a validade das estimações está sujeita aos resultados dos testes de Sargan e do teste de correlação serial sobre o resíduo transformado. O teste de Sargan verifica se os instrumentos utilizados nas regressões são válidos. A rejeição da hipótese nula implica que as restrições de sobreidentificação do modelo não são válidas. Todos os testes apontaram para a não rejeição de H0, o que implica que os instrumentos são válidos. O outro conjunto de testes busca verificar se os resíduos transformados de 1ª e 2ª ordem apresentam correlação serial. A hipótese nula, neste caso é ausência de correlação serial. Com exceção do modelo de convergência absoluta, todos os demais testes apontaram para a não rejeição de H0.

De uma forma geral, as estimações através do System-GMM acompanharam os resultados do painel de Efeitos Fixos, ainda que de forma implícita, pois as magnitudes dos coeficientes foram diferentes. Conforme ressaltam Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010., o modelo de Solow verde é mais robusto do que a estimação de uma CAK tradicional. Tal fato, de certa forma, foi comprovado neste trabalho. Utilizando dois modelos diferentes na estimação da CAK, os resultados também foram diferentes, pelo menos pelo método FMOLS. Já na estimação do modelo de convergência, utilizando-se dois estimadores diferentes, os resultados foram muito próximos.

Comparando o modelo de convergência de Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010. com os resultados obtidos por meio da estimação da CAK, pode-se concluir que, apesar dos resultados apontarem para a convergência das emissões, não se pode afirmar que ocorre uma curva no formato de U invertido, uma vez que os pontos de inflexão próximos à realidade foram obtidos por meio de uma curva no formato de N. Este resultado não implica que o modelo de Solow Verde não é capaz de descrever a CAK, conforme apontam Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010., mas sim que a estimação econométrica da CAK pode realmente ser passível de falhas, como ressaltam Stern (2004)STERN, D. I. The rise and fall of the environmental Kuznets curve. World Development, v. 32, p. 1419-1439, 2004. e Wagner (2008)WAGNER, M. The carbon Kuznets Curve: a cloudy picture emitted by bad econometrics? Resource and Energy Economics, v. 30, p. 388-408, 2008..

Sendo assim, constata-se que no período considerado as emissões de dióxido de carbono apresentaram uma taxa de crescimento equilibrado, o que representa uma taxa constante de crescimento rumo ao estado estacionário. Mas isso não implica que o nível de emissões de dióxido de carbono tenha sido baixo, mas apenas que o crescimento das emissões é dado por uma taxa constante.

5. Conclusão

Os resultados obtidos neste trabalho estão em conformidade com aqueles descritos em Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010., apontando para a convergência das emissões e, implicitamente, a existência de uma CAK. Porém, quando se utilizou a metodologia de estimação da CAK, foi apontada uma curva no formato N,21 21 Considerando o modelo que gerou estimativas para os pontos de inflexão cujos valores da renda per capita implicados estão mais próximos dos valores que a renda per capita assume na realidade para os países considerados. o que difere da curva no formato de U invertido, implícita no modelo proposto por Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010..

O modelo de convergência de Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010. possui duas limitações importantes. A primeira é considerar nulos os custos de abatimento do carbono. A outra limitação é que o modelo por si só não permite o cálculo empírico do ponto de inflexão. Embora os autores tenham obtido em sua análise do modelo teórico o instante em que ocorre o ponto de inflexão, bem como a renda per capita e as emissões associadas a esse ponto, sua estimação não foi implementada no artigo original. Enfim, mesmo que as emissões de dióxido de carbono tenham apresentado uma taxa de crescimento constante, é essencial o uso de políticas ambientais que criem os incentivos para redução do nível de emissão de poluentes, em especial pela adoção de tecnologias menos poluentes. Ratificando o que autores como Grossman e Krueger (1995)______. Economic growth and the environment. Quarterly Journal of Economics, v. 110, p. 353-377, 1995. e Hilton e Levinson (1998)HILTON, F. G. H.; LEVINSON, A. Factoring the Environmental Kuznets Curve: evidence from automotive lead emissions. Journal of Environmental Economics and Management, v. 35, p. 126-141, 1998. salientam, temos que o crescimento econômico por si só não é capaz de provocar necessariamente a redução das emissões a um nível de renda per capita factível. É preciso que haja uma ação conjunta de governos, países, indústrias e pessoas para que possamos ter um crescimento econômico associado a um ambiente mais limpo, obtendo dessa forma a sustentabilidade.

  • 1
    Embora a estimação da equação de convergência tenha sido feita por Brock e Taylor (2010)BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010., neste trabalho será adotada uma técnica mais moderna, com o uso da metodologia de painel dinâmico. Além disso, os dados que serão utilizados neste trabalho são mais atuais do que aqueles utilizados pelos referidos autores.
  • 2
    Uma discussão sobre o problema de decisão acerca da adoção de práticas de mitigação pode ser encontrada, por exemplo, em Baumol e Oates (1988)BAUMOL, W. J.; OATES, W. E. The theory of environmental policy - second edition. Cambridge: Cambridge University Press, 1988. e Perman et al. (2011)PERMAN, R.; MA, Y.; COMMON, M.; MADDISON, D.; McGILVRAY, J. Natural Resource and Environmental Economics - 4th edition. Essex: Pearson, 2011.. As medidas para mitigar a poluição englobam impostos de Pigou e permissões comercializáveis de emissões. O imposto de Pigou consiste em um imposto sobre as emissões de poluentes que procura internalizar os custos e fazer com que a firma venha a considerar as externalidades em seu processo interno de decisão de produção. A alíquota do imposto de Pigou é igual ao dano social marginal. As permissões comercializáveis de emissões são uma tentativa de obter o mesmo resultado do imposto de Pigou com o auxílio do mercado. A autoridade ambiental determina a quantidade agregada de emissões no montante em que o dano social marginal é igual ao custo marginal de mitigação. A seguir, emite permissões comercializáveis de emissões nesse montante e deixa a alocação das emissões entre fontes para o mercado. Considerando as curvas de dano social marginal e custo marginal de mitigação em um ambiente de certeza, a alíquota do imposto de Pigou é aquela que produz o montante permitido de emissões, de modo que o imposto de Pigou e as permissões comercializáveis de emissões produzem exatamente o mesmo resultado em termos práticos. Em um contexto de incerteza, a equivalência não permanece, e as duas práticas de mitigação podem produzir diferenças substanciais nos valores esperados do bem-estar social. A autoridade ambiental irá preferir, em um contexto de incerteza, uma política ou outra, dependendo da incerteza afetar o dano social marginal, o custo marginal de mitigação ou ambos. Para maiores detalhes, consultar Baumol e Oates (1988)BAUMOL, W. J.; OATES, W. E. The theory of environmental policy - second edition. Cambridge: Cambridge University Press, 1988., cap. 5 e Perman et al. (2011)PERMAN, R.; MA, Y.; COMMON, M.; MADDISON, D.; McGILVRAY, J. Natural Resource and Environmental Economics - 4th edition. Essex: Pearson, 2011., cap. 6.
  • 3
    Maiores detalhes em Breitung e Pesaram (2005)BREITUNG, J.; PESARAM, M. H. Unit roots and cointegration in panels. Frankfurt: Deutsche Bundesbank, 2005..
  • 4
    Com exceção do teste Hadri (2000).
  • 5
    Engle, R.F.; Granger, C.W.J. Cointegration and error correction: representation, estimation, and testing. Econometrica, v.55, n.2, pp.251-76, 1987.
  • 6
    Decorrente da estimação da Equação (25) via Mínimos Quadrados Ordinários (MQO).
  • 7
    ApudIslam (1995)ISLAM, N. Growth empirics: A panel data approach. The Quarterly Journal of Economics, v. 110, n. 4, p. 1127-1170, 1995..
  • 8
    ApudIslam (2003)______. What have we learnt from the convergence debate? Journal of Economic Surveys, v. 17, n. 3, p. 309-362, 2003..
  • 9
    Conforme Islam (2003)______. What have we learnt from the convergence debate? Journal of Economic Surveys, v. 17, n. 3, p. 309-362, 2003., a frequência estimada em vários estudos de painel é definida por um intervalo de cinco anos, porém, em função da quantidade de dados, optou-se por definir um intervalo de seis anos a fim de incorporar na análise o maior intervalo temporal possível entre a primeira e a última observação.
  • 10
    O banco de dados CAIT (Climate Analysis Indicators Tool), elaborado pelo World Resources Institute (WRI), poderia ser usado alternativamente para essa análise. Ele possui a vantagem de estar disponível publicamente pela internet no sítio http://cait2.wri.org/wri, ao contrário da base da IEA. Ambas as bases possuem alta confiabilidade. Neste artigo optamos por usar o banco de dados do IEA.
  • 11
    Foram feitos testes nas demais especificações (intercepto e passeio aleatório), que apontaram para a presença de raiz unitária.
  • 12
    Pois H0: as séries possuem raiz unitária.
  • 13
    Foram feitos os mesmos testes sobre a primeira diferença das variáveis, porém optou-se por não reportá-los neste trabalho. Os testes apontaram que todas as variáveis são estacionárias na primeira diferença.
  • 14
    Optou-se por não reportar a estatística v, uma vez que os resultados foram próximos aos da estatística rho, e também porque podem existir distorções de tamanho e poder devido à dimensão das séries utilizadas neste trabalho.
  • 15
    Foram feitas 1.000 replicações.
  • 16
    Poolabilidade pode ser entendida como a reunião ou empilhamento dos dados. Sendo assim, consiste em impor um mesmo parâmetro para todas as unidades do painel ao longo do tempo.
  • 17
    No teste feito neste trabalho, a hipótese alternativa foi um modelo de Efeitos Fixos.
  • 18
    Todas as variáveis estão expressas em logaritmos, exceto Δeit.
  • 19
    Conforme Islam (2003)______. What have we learnt from the convergence debate? Journal of Economic Surveys, v. 17, n. 3, p. 309-362, 2003..
  • 20
    Reescrevendo de outra forma: β1 =γ – 1. Neste trabalho foram reportados os valores do coeficiente β1.
  • 21
    Considerando o modelo que gerou estimativas para os pontos de inflexão cujos valores da renda per capita implicados estão mais próximos dos valores que a renda per capita assume na realidade para os países considerados.
  • Os autores agradecem os comentários de um parecerista anônimo.

Referências

  • ARELLANO, M.; BOVER, O. Another look at the instrumental variable estimation of error components models. Journal of Econometrics, v. 68, n. 1, p. 29-51, 1995.
  • AUCI, S.; BECCHETTI, L. The instability of the adjusted and unadjusted Environmental Kuznets Curves. Ecological Economics, v. 60, p. 282-298, 2006.
  • BALTAGI, B. H. Econometric Analysis of Panel Data 3. ed. Inglaterra: John Wiley & Sons, 2005.
  • BANCO MUNDIAL. World Development Indicators Estados Unidos, 2011.
  • BARRO, R. J.; SALA-i-MARTIN. Economic Growth - second edition Cambridge, MA: MIT Press, 2004.
  • BAUMOL, W. J.; OATES, W. E. The theory of environmental policy - second edition. Cambridge: Cambridge University Press, 1988.
  • BLUNDELL, R.; BOND, S. Initial conditions and moment restrictions in dynamic panel data models. Journal of Econometrics, v. 87, p. 115-143, 1998.
  • BREITUNG, J.; PESARAM, M. H. Unit roots and cointegration in panels Frankfurt: Deutsche Bundesbank, 2005.
  • BROCK, W. A.; TAYLOR, M. S. The green Solow model. Journal of Economic Growth, v. 15, p. 127-153, 2010.
  • CAMERON, A. C.; TRIVEDI, P. K. Microeconometrics: Methods and Applications New York: Cambridge University Press, 2005.
  • CAMPBELL, J. Y.; MANKIW, N. G. Consumption, income, and interest rates: reinterpreting the time series evidence. NBER Macroeconomics Annual, v. 4, p. 185-216, 1989.
  • CARVALHO, T. S.; ALMEIDA, E. A hipótese da Curva de Kuznets Ambiental Global: uma perspectiva econométrico-espacial. Estudos Econômicos, v. 40, n. 3, p. 587-615, p. 1149-1159, 2010.
  • CHENG, S. Y.; WEI, J. S.; HOU, H. A cointegration analysis of purchasing power parity and country risk. International Journal of Business and Economics, v. 7, p. 199-211, 2008.
  • COLE, M. A.; ELLIOT, R. J. R.; AZAR, A. K. The determinants of trade in pollution intensive industries: North South evidence University of Birmingham, 2001.
  • DASGUPTA, S.; LAPLANTE, B.; WANG, H.; WHEELER, D. Confronting the Environmental Kuznets Curve. Journal of Economic Perspectives, v. 16, p. 147-168, 2002.
  • DE BRUYN, S. M. Economic Growth and the Environment Dordrecht: Kluwer Academic Publishers, 2000.
  • DINDA, S. A theoretical basis for the Environmental Kuznets Curve. Ecological Economics, v. 53, p. 403 - 413, 2005.
  • DINIZ, E. M. Crescimento, poluição e o Protocolo de Quioto: uma avaliação do caso brasileiro São Paulo: Banco Santos / Universidade de Oxford, 2001.
  • ______. Lessons from the Kyoto protocol. Ambiente e Sociedade, v. 10, n. 1, p. 27-38, 2007.
  • FONSECA, L. N.; RIBEIRO, E. P. Preservação ambiental e crescimento econômico no brasil. In: Encontro de economia da Região Sul, VII. Rio Grande do Sul: [s.n.], 2005.
  • GROSSMAN, G. M.; KRUEGER, A. B. Environmental impacts of a North American free trade agreement. National Bureau of Economic Research, 1991.
  • ______. Economic growth and the environment. Quarterly Journal of Economics, v. 110, p. 353-377, 1995.
  • HESTON, A.; SUMMERS, R.; BETTINA, A. Penn World Table Version 7.0 Center for International Comparisons of Production, Income and Prices at the University of Pennsylvania, Maio, 2011.
  • HILTON, F. G. H.; LEVINSON, A. Factoring the Environmental Kuznets Curve: evidence from automotive lead emissions. Journal of Environmental Economics and Management, v. 35, p. 126-141, 1998.
  • INTERNACIONAL ENERGY AGENCY. CO2 Emissions from Fuel Combustion Paris: IEA, 2010.
  • ISLAM, N. Growth empirics: A panel data approach. The Quarterly Journal of Economics, v. 110, n. 4, p. 1127-1170, 1995.
  • ______. What have we learnt from the convergence debate? Journal of Economic Surveys, v. 17, n. 3, p. 309-362, 2003.
  • LEVIN, A.; LIN, C. F.; CHU, C. S. J. Unit roots tests in panel data: asymptotic finite sample properties. Journal of Econometrics, v. 108, p. 1-24, 2002.
  • LIST, J. A. Have air pollutant emissions converged among U.S. regions? Evidence from unit root tests. Southern Economic Journal, v. 66, n. 1, p. 144-155, 1999.
  • MEADOWS, D. H. et al. The limits to growth New York: Universe Books, 1972.
  • PANAYOTOU, T. Empirical tests and policy analysis of environmental degradation at different stages of economic development. Technology and Employment Programme, 1993.
  • PEDRONI, P. Critical values for cointegration tests in heterogeneous panels with multiple regressors. Oxford Bulletin of Economics and Statistics, v. 61, p. 653-670, 1999.
  • ______. Fully modified ols for heterogeneous cointegrated panels. Advances in Econometrics, v. 15, p. 93-130, 2000.
  • ______. Panel Cointegration: Asymptotic and Finite Sample Properties of Pooled Time Series Tests with an Application to the PPP Hypothesis. Econometric Theory, v. 20, p. 597-625, 2004.
  • PERMAN, R.; MA, Y.; COMMON, M.; MADDISON, D.; McGILVRAY, J. Natural Resource and Environmental Economics - 4th edition. Essex: Pearson, 2011.
  • ROMER, D. Advanced Macroeconometrics 2. ed. Nova Iorque: McGraw-Hill, 2001.
  • SELDEN, T. M.; SONG, D. Environmental quality and development: is there a Kuznets Curve for air pollution emissions? Journal of Environmental Economics and Management, v. 27, n. 2, p. 147-162, 1994.
  • SHAFIK, N.; BANDYOPADHYAY, S. Economic growth and environmental quality: a time series and cross-country evidence. Journal of Environmental Economics and Management, v. 4, p. 1-24, 1992.
  • SOLOW, R. A contribution to the theory of economic growth. Quarterly Journal of Economics, v. 70, p. 65-94, 1956.
  • STERN, D. I. The rise and fall of the environmental Kuznets curve. World Development, v. 32, p. 1419-1439, 2004.
  • STOKEY, N. L. Are there limits to growth? International Economic Review, v. 39, n. 1, p. 1-31, 1998.
  • WAGNER, M. The carbon Kuznets Curve: a cloudy picture emitted by bad econometrics? Resource and Energy Economics, v. 30, p. 388-408, 2008.
  • WESTERLUND, J. Testing for error correction in panel data. Oxford Bulletin of Economics and Statistics, v. 69, p. 709-748, 2007.

Datas de Publicação

  • Publicação nesta coleção
    Mar 2015

Histórico

  • Recebido
    16 Maio 2013
  • Aceito
    11 Set 2014
Departamento de Economia; Faculdade de Economia, Administração, Contabilidade e Atuária da Universidade de São Paulo (FEA-USP) Av. Prof. Luciano Gualberto, 908 - FEA 01 - Cid. Universitária, CEP: 05508-010, Tel.: (55 11) 3091-5803/5947 - São Paulo - SP - Brazil
E-mail: estudoseconomicos@usp.br