Hostname: page-component-69cd664f8f-gvmbd Total loading time: 0 Render date: 2025-03-13T05:58:52.146Z Has data issue: false hasContentIssue false

Pollution, public debt, and growth: the question of sustainability

Published online by Cambridge University Press:  12 March 2025

Marion Davin
Affiliation:
CEE-M, Univ Montpellier, CNRS, INRAE, Institut Agro, Montpellier, France
Mouez Fodha
Affiliation:
Paris School of Economics, University Paris 1 Panthéon-Sorbonne and Aix Marseille Univ, CNRS, AMSE, Marseille, France
Thomas Seegmuller*
Affiliation:
Aix Marseille Univ, CNRS, AMSE, Marseille, France
*
Corresponding author: Thomas Seegmuller; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

This paper examines an endogenous growth model that allows us to consider the dynamics and sustainability of debt, pollution, and growth. Debt evolves according to the financing adaptation and mitigation efforts and to the damages caused by pollution. Three types of features are important for our analysis: the technology through the negative effect of pollution on TFP; the fiscal policy; the initial level of pollution and debt with respect to capital. Indeed, if the initial level of pollution is too high, the economy is relegated to an endogenous tipping zone where pollution perpetually increases relatively to capital. If the effect of pollution on TFP is too strong, the economy cannot converge to a stable and sustainable long-run balanced growth path. If the income tax rates are high enough, we can converge to a stable balanced growth path with low pollution and high debt relative to capital. This sustainable equilibrium can even be characterized by higher growth and welfare. This last result underlines the role that tax policy can play in reconciling debt and environmental sustainability.

Type
Articles
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2025. Published by Cambridge University Press

1. Introduction

In the context of growing public debt and rising costs related to global pollution, it is crucial to have a clear understanding of the interplay between debt and global pollution dynamics. Today, several countries that are particularly vulnerable to the impacts of climate change also find themselves burdened with high levels of debt. The economic repercussions of the COVID-19 pandemic further exacerbated this situation (see e.g Dibley et al. Reference Dibley, Wetzer and Hepburn2021). This issue is widely brought up for developing countries,Footnote 1 but is also a major concern for developed countries. Indeed, since 2008, we have observed in all groups of countries an increasing trend in the share of debt in GDP (Figure 1).

Figure 1. Evolution of the central government debt to GDP ratio by groups of countries (World bank classification in 2023). IMF global debt database.

Concerning expenses related to global pollution, in addition to the substantial investments linked to the transition to a less polluting economy (i.e. mitigation), the expenses associated with adaptation are expected to grow in all countries (IPCC, Reference Birkmann, Liwenga, Pandey, Boyd, Djalante, Gemenne, Filho, Pinho, Stringer, Wrathall, Pörtner, Roberts, Tignor, Poloczanska, Mintenbeck, Alegría, Craig, Langsdorf, Löschke, Möller, Okem and Rama2022). Meanwhile, the connections between public debt and global environmental challenges can be illustrated by the positive and substantial impacts of climate vulnerability on debt (Buhr et al. Reference Buhr, Donovan, Kling, Lo, Murinde, Pullin and Volz2018), suggesting the existence of a detrimental cycle wherein vulnerability perpetuates itself through public debt management. Moreover, the rise in disaster-related losses will result in reduced tax revenues. Zenios (Reference Zenios2024) gives an overview of the direct and indirect channels and suggests a possible doom loop between climate change and sovereign debt. We consider in this paper these different dimensions to examine the interplay between debt and pollution dynamics.

As reported by IMF (2023), policymakers face a fundamental tradeoff. On the one hand, relying on spending-based measures to achieve emission goals and to adapt could lead to a substantial increase in public debt. On the other hand, limited environmental action exposes the world to adverse consequences from global pollution, increasing the cost of adaptation. The recent proposals from the European Commission regarding new economic governance rules highlight the interconnected nature of questions surrounding fiscal sustainability and environmental concerns. In particular, the Commission considers climate change as a structural trend representing a challenge to the financial stability of Member States’ public finances.Footnote 2 The Pisani-Ferry and Mahfouz (Reference Pisani-Ferry and Mahfouz2023) report is in line with this argument. In the case of France, it recommends using debt to finance the investments needed for the ecological transition. The use of debt should be limited to “green” investments that have a positive impact on the climate and generate long-term economic returns. This policy must be accompanied by more progressive taxation.

This paper contributes to this debate. Within an endogenous growth framework, we study the dynamic path of pollution, debt, and economic growth when public authorities finance mitigation and adaptation to tackle the damages caused by the pollution stock. From a normative perspective, we look at how fiscal and environmental policy instruments can be used to guarantee sustainability and improve welfare.

We develop an overlapping generations (OLG) model where debt, pollution, and growth are endogenous. Households live for two periods and save through two assets, capital, the source of growth, and public debt. The government issues debt securities because taxes on capital and labor incomes do not cover public expenses for pollution mitigation and adaptation, and the service of past debt. The pollution stock evolves with productive capital and mitigation measures, and is a source of damage by reducing aggregate productivity (TFP). Technology plays an important role: the final good is produced using an AK-type production function to have a simple engine of growth, but considering that the adaptation policy of the government can dampen the negative effect of pollution on TFP. In other words, the government can affect pollution using two instruments: government spending of mitigation reduces the pollution flow, whereas the government spending of adaptation offsets the negative effect of pollution on TFP. The first type of expenditure could be on waste management, the installation of filters to reduce harmful emissions or carbon removal solutions, while the second one includes building infrastructure to protect against rising sea levels or creating drought-resistant agriculture.

The long-run equilibria analysis shows that two balanced growth paths (BGPs) may exist, defined by constant ratios of debt and pollution over capital. One is characterized by a low pollution-to-capital ratio and a high debt-to-capital ratio while the reverse is observed for the other. When TFP vulnerability to pollution is not too high, the growth rate is higher at the BGP with low pollution and high debt, which also means higher welfare at this BGP. In that case, the higher long-term growth rate is at the expense of a high level of debt per unit of capital. This suggests a crowding-in effect of debt on growth.

The analysis of dynamics allows us to show that depending on fiscal policy, the TFP vulnerability to pollution, and initial conditions of capital, debt, and pollution stocks, the economy either converges to the BGP with low pollution over capital, collapses or experiences a perpetual increase in pollution over capital. More precisely, we show that the BGP with a high pollution-to-capital ratio is not sustainable, as the economy cannot converge to this state, while the convergence to the BGP with low pollution-to-capital may be possible for a sufficiently high labor taxation and a reasonable TFP vulnerability to pollution. We identify an extreme case in which sustainability is completely excluded. None of the two BGPs is stable and the economy either collapses or is characterized by a perpetual increase in pollution over capital. This unfavorable situation occurs when the tax rates on labor and capital income are low and TFP vulnerability to pollution is high.

If the BGP with high pollution over capital is the only unstable one, it is a saddle toward which the stock variables cannot converge. Therefore, the stable manifold of this equilibrium defines an endogenous tipping zone (ETZ). If the initial conditions on debt, pollution, and capital are such that the economy is in this zone, the dynamics of pollution relative to capital explode.Footnote 3 Interestingly, the higher the debt relative to pollution, the easier the economy can be relegated to the ETZ. It corresponds to a form of debt vulnerability: higher debt favors unsustainable dynamic paths for pollution.

Finally, in the case of a long-run sustainability, we investigate if the policy can still improve welfare. Therefore, we analyze the effect of the fiscal policy, mitigation and adaptation on the level of the stable BGP. When TFP vulnerability to pollution is not excessive, increasing taxation enhances welfare along the stable BGP, while the effects of environmental policy instruments depend greatly on their efficiency. Insufficiently efficient mitigation in reducing pollution may worsen the pollution-to-capital ratio, whereas adaptation improves welfare only if TFP responds strongly to such expenditure. In general, this analysis shows the crucial role played by adaptation.

Previous contributions studied the macroeconomic implications of the interplay between public debt and environmental factors without (Heijdra et al. Reference Heijdra, Kooiman and Ligthart2006) and with public mitigation Fodha and Seegmuller (Reference Fodha and Seegmuller2012, Reference Fodha and Seegmuller2014).Footnote 4 Nonetheless, debt is often considered as an exogenous instrument, and its dynamic evolution of financing adaptation and mitigation together with the dynamic path of pollution is left aside, meaning that the question of sustainability is not properly addressed. In Andersen et al. (Reference Andersen, Bhattacharya and Liu2020), a dynamic debt scheme is designed to finance public mitigation at each period and to improve welfare. However, in contrast to our framework, they cannot discuss the implications on economic growth. In Baret and Menuet (Reference Baret and Menuet2024), debt allows financing mitigation expenditure. However, this paper cannot address the question of sustainability by assuming a constant long-term debt-to-output ratio and a stabilizing rule ensuring convergence towards this objective. Moreover, they leave aside adaptation expenditures while they must be significantly integrated into environmental spending. We go one step further by developing a model that incorporates endogenous public debt dynamics and acknowledges the impact of environmental issues on fiscal sustainability. This dimension seems essential to address the economic consequences of global pollution and highlights how it exerts its influence on sovereign debt. Our paper is related to some recent papers that consider endogenous public debt and its dynamics. For example, Boly et al. (Reference Boly, Combes, Menuet, Minea, Motel and Villieu2022) examine the relationship between public and environmental debt in a framework with public mitigation but no adaptation. We depart from this paper focusing on fiscal sustainability and considering the economic damage entailed by pollution stock and its impact on debt. Catalano et al. (Reference Catalano, Forni and Pezzolla2020) examine the role of fiscal policy in climate change adaptation. They use a calibrated macroeconomic model of an open economy, that does not allow them to explicitly identify the interplay between debt dynamics, growth, and the environment. More generally, and in contrast to the literature that assumes limits on long-term public debt (Baret and Menuet, Reference Baret and Menuet2024; Boly et al. Reference Boly, Combes, Menuet, Minea, Motel and Villieu2022; Seghini and Dees, Reference Seghini and Dees2024), we do not impose any restrictions on sovereign debt. This is crucial for studying environment-debt interactions and sustainability.

The rest of this paper is organized as follows. Section 2 presents an OLG model in which pollution is proportional to capital stock, and the government issues debt and imposes taxes on capital and labor incomes for financing adaptation and mitigation expenditures. Section 3 defines the equilibrium. Section 4 studies the balanced growth paths and examines the existence and multiplicity of BGPs. Section 5 considers the dynamics and the possibility of an endogenous tipping zone. Section 6 presents some policy implications. The final section provides the conclusion.

2. The model

We consider a dynamic model with pollution and three types of agents, firms, consumers, and a government. Time is discrete, $t=0, 1,\ldots, +\infty$ , and there is no uncertainty.

2.1 Production

We consider an AK model of economic growth in which TFP decreases with pollution stock $P_{t}$ .Footnote 5 Considering that pollution or climate change is detrimental to production is particularly relevant in addressing debt and environmental issues, as it allows to focus on funding adaptation efforts, extending beyond mere mitigation. The need for adaptation strategies will increase with the intensification of climate change impacts. These adaptations come with associated costs, such as building infrastructure to protect against rising sea levels, creating drought-resistant agriculture or developing air conditioning systems to cope with heat waves. Public action can provide the necessary financial resources to implement these adaptation measures, reducing the vulnerability of countries to environmental shocks. We thus assume that the capacity for adaptation reduces the incremental damage caused by pollution stock. This ability to adapt is ensured by the public authorities, who devote an amount $G_{1t}$ , specifically for this purpose.Footnote 6

Therefore, the production $Y_t$ is given by:

(1) \begin{equation} Y_t=A\left (\frac {G_{1t}}{P_t}\right )K_t^\alpha (\bar {K}_tL_t)^{1- \alpha } \end{equation}

with $K_t$ the capital, $\bar {K}_t$ the aggregate level of capital, $L_t$ the labor, $\alpha \in (0,1)$ and:

Assumption 1. $A(X)$ is a strictly increasing function, with $ A(0)=A_0\geqslant 0$ , $A(\!+\infty )=A_1\lt +\infty$ , and $A^{\prime }(X)X/A(X)\in (0,1)$ .

This last assumption implies an elasticity of production to the adaptation to pollution ratio lower than one. When the adaptation to pollution ratio goes up it causes a less-than-proportional increase in productivity.

Note that $A(X)$ may capture the fact that climate change destroys a part of aggregate output at each period (Golosov et al. Reference Golosov, Hassler, Krusell and Tsyvinski2014; Dietz and Stern, Reference Dietz and Stern2015). It can also represent the health effects of global pollution stock or the impacts of a change in temperature, which results in reduced aggregate productivity (Dasgupta et al. Reference Dasgupta, van Maanen, Gosling, Piontek, Otto and Schleussner2021; Burke et al. Reference Burke, Hsiang and Miguel2015).

Example: we can consider the following specifications for $A(X)$ :

(2) \begin{equation} A(X)=\frac {A_1 X}{1+X} \end{equation}

with $A(0)=0$ . This function is increasing and concave, with:

(3) \begin{equation} \frac {A^{\prime }(X)X}{A(X)}=\frac {1}{1+X}\in (0, 1) \end{equation}

Let $r_t$ be the interest rate and $w_t$ the wage. At equilibrium, we have $ \bar {K}_t=K_t$ and assuming that the population in this economy is equal to one, labor input is $L_t=1$ . Therefore, profit maximization gives:

(4) \begin{align} r_t&= \alpha A\left (\frac {G_{1t}}{P_t}\right ) \end{align}
(5) \begin{align} w_t &= (1-\alpha ) A\left (\frac {G_{1t}}{P_t}\right )K_t \\[6pt]\nonumber\end{align}

Returns of factor being a positive function of productivity, they decrease with pollution.

2.2 Government

We consider public actions to tackle environmental issues. Public spending $G_t$ linearly increases with GDP:

(6) \begin{equation} G_t=g Y_t \end{equation}

and are divided into public spending which attenuates the effect of pollution on production $G_{1t}$ , i.e. adaptation to climate effect, and mitigation $G_{2t}$ :

(7) \begin{eqnarray} G_{it}=g_i Y_t \end{eqnarray}

for $i=1,2$ , with $g_1+g_2=g$ .

Environmental policy instruments consist of public spending on both mitigation and adaptation. While most of the literature has focused on their potential substitutability, these two strategies are now seen as simultaneously needed in the face of climate emergencies. This need is reflected in international ambitions to balance climate finance spending between the two strategies (Sadler et al. Reference Sadler, Ranger, Fankhauser, Marotta and O’Callaghan2024).

Since we are in an endogenous growth framework, we assume that the government determines its public spending for adaptation and mitigation by fixing their amount per GDP unit. This means that the policy of adaptation will be determined by $g_1$ and the policy of mitigation by $g_2$ .

To finance these spending, the government collects taxes on labor and capital incomes, $\tau _L$ and $\tau _K$ , and issues debt $B_t$ . Therefore, its expenditures include repayment of debt and interest payments. The government faces the following budget constraint at each period:

(8) \begin{equation} R^b_tB_t+G_t=B_{t+1}+\tau _L w_t+\tau _K r_tK_t \end{equation}

with $R^b_t$ the interest factor of debt and $B_0\gt 0$ the initial stock of debt. We are in an economy with a positive initial stock of public debt. The different policy parameters as well as the interest factor of debt and the income will determine how public debt evolves through time. We will precisely study the interplay between debt accumulation, dynamics of pollution stock, and growth.

2.3 Pollution

The stock of pollution increases with the emission flow and partly leaves the atmosphere through a natural process in a share $0\lt m\lt 1$ . Emission flow is assumed to be proportional to the stock of capital. Mitigation measures $G_{2t}$ are implemented by public authorities to further reduce pollution flows. Mitigation expenditure may include investment in carbon capture and sequestration or carbon dioxide removal solutions. The stock of pollution evolves according to:

(9) \begin{equation} P_{t+1}=(1-m)P_{t}-\psi G_{2t}+\mu K_{t} \end{equation}

The parameter $\psi \gt 0$ captures the efficiency of public mitigation and $\mu \gt 0$ the pollution flow resulting from capital stock, as in Heijdra et al. (Reference Heijdra, Kooiman and Ligthart2006) or Chiroleu-Assouline and Fodha (Reference Chiroleu-Assouline and Fodha2014). Capital is the source of pollution. Indeed, capital accumulation favors the production of pollution-intensive goods. This is consistent with the evidences (e.g. Cole et al. Reference Cole, Elliott and Shimamoto2005; Andersen, Reference Andersen2017). For example, Cole et al. (Reference Cole, Elliott and Shimamoto2005) find that industrial processes that tend to be physical capital intensive generate more emissions than less capital-intensive processes. Andersen (Reference Andersen2017) obtains that pollution emissions are higher for industries that use more intensively tangible assets, such as physical capital, than intangible assets, such as labor.

Note that since $G_{2t}$ will increase with income, public adaptation will also have a direct negative effect on the pollution stock by increasing TFP, income and therefore mitigation. In addition, polluting capital which will be derived from a portion of savings will increase with labor income and therefore with TFP and public adaptation.

This stock of pollution only affects the real side of the economy through its negative effect on the TFP. We will not consider a direct negative effect of pollution on households welfare.

2.4 Consumers

Consumers are in overlapping generations. The population size of each generation is constant and normalized to one. Each consumer lives for two periods, consumes in both periods, and saves through two assets, public debt and capital. Capital depreciates at rate $\delta \in (0,1)$ , meaning that return on capital is given by $1-\delta +(1-\tau _K) r_{t}$ .

The utility function of the generation born in $t$ is given by:

(10) \begin{equation} U(c_t, d_{t+1})= \ln c_t + \beta \ln d_{t+1} \end{equation}

with $c_t$ the consumption when young, $d_{t+1}$ the consumption when old, and $\beta \in (0,1)$ . The household maximizes her utility under the two budget constraints:

(11) \begin{align} c_t+K_{t+1}+B_{t+1}&= (1-\tau _L)w_t \end{align}
(12) \begin{align} d_{t+1}&= [1-\delta +(1-\tau _K) r_{t+1}]K_{t+1}+R^b_{t+1}B_{t+1} \\[9pt]\nonumber\end{align}

The household can solve her optimal behavior in two steps. Given her saving, she determines her portfolio choice between capital and debt holding. Given this choice, she chooses her optimal saving. Maximizing the utility with respect to $K_{t+1}$ and $B_{t+1}$ , we obtain the following equation:

(13) \begin{equation} 1-\delta +(1-\tau _K) r_{t+1}=R^b_{t+1}\equiv R_{t+1} \end{equation}

Bonds and capital assets provide the same return, which means that they are perfect substitutes.Footnote 7 Given this result, the utility maximization with respect to $K_{t+1}+B_{t+1}$ gives the saving function:Footnote 8

(14) \begin{equation} K_{t+1}+B_{t+1}=\frac {\beta }{1+\beta } (1-\tau _L)w_t \end{equation}

3. Equilibrium

We define an equilibrium as a function of capital, debt, and pollution stocks. Market clearing is obtained substituting (4), (5) and (13) in (8), (9) and (14). We get the following functions:

(15) \begin{align} K_{t+1}+B_{t+1}&=\frac {\beta }{1+\beta } (1-\tau _L)(1-\alpha ) A\left (\frac { G_{1t}}{P_t}\right )K_t \end{align}
(16) \begin{align} B_{t+1}&= R_tB_t+g A\left (\frac {G_{1t}}{P_t}\right )K_t-(\tau _L (1-\alpha )+\tau _K \alpha ) A\left (\frac {G_{1t}}{P_t}\right )K_t \end{align}
(17) \begin{align} P_{t+1}&=(1-m)P_t-\psi g_2 A\left (\frac {G_{1t}}{P_t}\right )K_t +\mu K_t \\[6pt]\nonumber\end{align}

We introduce the following variables to conduct our analysis: the growth factor $\gamma _{t+1}\equiv K_{t+1}/K_t$ , debt per unit of capital $b_t \equiv B_t/K_t$ , and pollution per unit of capital $\pi _t \equiv P_t/K_t$ . Using the government budget (7), public adaptation per unit of capital can thus be written as a function of $\pi _t$ :

(18) \begin{equation} \frac {G_{1t}}{P_t}= \frac {g_1 }{\pi _t}A\left (\frac {G_{1t}}{P_t}\right ) \end{equation}

Since $A^{\prime }(x)x/A(x)\in (0,1)$ , this equation implicitly defines a decreasing function $G_{1t}/P_t= \varepsilon (\pi _t)$ if $\lim _{x \rightarrow 0}A(x)/x\gt \pi _t/g_1$ (see Appendix A for details). Using Assumption1 and equation (18), we further have $ \varepsilon (0)=+\infty$ and $\varepsilon (\!+\infty )=0$ . As pollution per unit of capital increases, the amount of public spending dedicated to addressing pollution-related issues for each unit of pollution decreases. Because of the productivity cost entailed by pollution, the higher the pollution per unit of capital, the lower the resources that can be allocated to adaptation. This may be counteracted if the government decides to allocate a larger proportion of its budget to climate adaptation efforts, $i.e$ if $g_1$ increases.

The total factor productivity can thus be expressed as a function of pollution per capital. We have $A(G_{1t}/P_t)=A[\varepsilon (\pi _t)]\equiv a(\pi _t)$ , with $a^{\prime }(\pi _t)\lt 0$ , $a(0)=A[\varepsilon (0)]=A(\!+\infty )=A_1$ , and $a(\!+\infty )=A[\varepsilon (\!+\infty )]=A(0)=A_0$ . Similarly, the interest factor can be written as $R_{t+1}=1-\delta +(1-\tau _K)\alpha a(\pi _{t+1})\equiv R(\pi _{t+1})$ , with $R^{\prime }(\pi _{t+1})\lt 0$ .

Example (continued): considering our example given by equation (2), we have $\varepsilon (\pi _t)=g_1 A_1/\pi _t-1$ , which implies that:

(19) \begin{equation} a(\pi _t)=A_1-\frac {\pi _t}{g_1} \end{equation}

which requires that $\pi _t\lt g_1A_1$ .

Now, we can rewrite the dynamic system (15)–(17) as follows:

(20) \begin{align} \gamma _{t+1}+b_{t+1}\gamma _{t+1}&= \Sigma (1-\tau _L)(1-\alpha ) a(\pi _t) \end{align}
(21) \begin{align} b_{t+1}\gamma _{t+1} &= R(\pi _t) b_t+g a(\pi _t)-(\tau _L (1-\alpha )+\tau _K \alpha ) a(\pi _t) \end{align}
(22) \begin{align} \pi _{t+1}\gamma _{t+1} &=(1-m)\pi _t-\psi g_2 a(\pi _t) +\mu\\[12pt]\nonumber \end{align}

with $\Sigma \equiv \frac {\beta }{1+\beta }$ the saving rate.

Rearranging equations (20)–(22), we finally obtain:

(23) \begin{align} \gamma _{t+1}&= a(\pi _t)\left [\Sigma (1-\tau _L)(1-\alpha ) +(\tau _L (1-\alpha )+\tau _K \alpha )-g-\frac {R(\pi _t)}{a(\pi _t)} b_t\right ] \end{align}
(24) \begin{align} b_{t+1} &= \frac {\frac {R(\pi _t)}{a(\pi _t)} b_t+g -(\tau _L (1-\alpha )+\tau _K \alpha ) }{\Sigma (1-\tau _L)(1-\alpha ) +(\tau _L (1-\alpha )+\tau _K \alpha )-g-\frac {R(\pi _t)}{a(\pi _t)} b_t} \end{align}
(25) \begin{align} \pi _{t+1} &=\frac {(1-m)\frac {\pi _t}{a(\pi _t)}-\psi g_2 +\frac {\mu }{a(\pi _t)} }{\Sigma (1-\tau _L)(1-\alpha ) +(\tau _L (1-\alpha )+\tau _K \alpha )-g-\frac {R(\pi _t)}{a(\pi _t)} b_t}\\[12pt]\nonumber \end{align}

Equations (24) and (25) give the dynamics of $(b_t, \pi _t)$ for $t\geq 0$ , taking into account that both $b_t$ and $\pi _t$ are predetermined variables. The dynamics of debt per capital is coupled with the dynamics of pollution because the depreciation of capital is not complete. This implies that the cost of capital falls with productivity loss but not in the same proportion as growth $\left (\frac {R(\pi _t)}{a(\pi _t)}\text { still depends on }\pi _t\right )$ . The value of the growth factor $\gamma _t$ is deduced from these two variables using equation (23). To conduct our analysis, we focus on relevant situations in which the growth factor, debt per unit of capital, and pollution per unit of capital are all positive. To ensure $\gamma _{t+1}\gt 0$ , $b_{t+1}\gt 0$ and $\pi _{t+1}\gt 0$ , we assume a primary deficit and the following restrictions:

\begin{eqnarray*} \Sigma (1-\tau _L)(1-\alpha ) &\gt &\frac {R(\pi _t)}{a(\pi _t)} b_t+g - (\tau _L (1-\alpha )+\tau _K \alpha ) \\ \mu &\gt &a(\pi _t) \psi g_2 \end{eqnarray*}

They are satisfied under the next assumption:

Assumption 2.

\begin{eqnarray*} b_t &\lt &\frac {a(\!+\infty )}{R(\!+\infty )}[\Sigma (1-\tau _L)(1-\alpha ) -g + \tau _L (1-\alpha )+\tau _K \alpha ] \\ g&\gt & \tau _L (1-\alpha )+\tau _K \alpha \\ \mu &\gt & a(0) \psi g_2 \end{eqnarray*}

The first inequality characterizes an upper bound for public debt, which increases with the amount of savings and decreases with the primary deficit, the second a primary deficit, and the third emission intensity per unit of capital higher than the efficiency of pollution abatement.

4. Balanced growth paths: multiplicity and main features

We focus here on long-run equilibria. We first show the existence and multiplicity of BGPs. Then, we investigate the main features of these equilibria. We will in particular understand how they are ranked according to the levels of debt per capital, pollution per capital, and growth.

4.1 Existence and multiplicity of BGPs

Along a balanced growth path, capital, debt, and pollution grow at a constant rate $\gamma -1$ . A balanced growth path is thus characterized by $b_t=b_{t+1}=b$ and $\pi _t=\pi _{t+1}=\pi$ solving (24) and (25). Hence, it is a stationary solution $(b, \pi )$ to:

(26) \begin{align} b &= \frac {\frac {R(\pi )}{a(\pi )} b+g -(\tau _L (1-\alpha )+\tau _K \alpha ) }{\Sigma (1-\tau _L)(1-\alpha ) +(\tau _L (1-\alpha )+\tau _K \alpha )-g-\frac {R(\pi )}{a(\pi )} b} \end{align}
(27) \begin{align} \pi &=\frac {(1-m)\frac {\pi }{a(\pi )}-\psi g_2 +\frac {\mu }{a(\pi )} }{\Sigma (1-\tau _L)(1-\alpha ) +(\tau _L (1-\alpha )+\tau _K \alpha )-g-\frac {R(\pi )}{a(\pi )} b} \\[12pt]\nonumber \end{align}

Given such a solution, the growth factor corresponds to:

(28) \begin{equation} \gamma = a(\pi )\left [\Sigma (1-\tau _L)(1-\alpha ) +(\tau _L (1-\alpha )+\tau _K \alpha )-g-\frac {R(\pi )}{a(\pi )} b\right ] \end{equation}

The ratio of (26) and (27) gives:

(29) \begin{equation} b=\frac {g -(\tau _L (1-\alpha )+\tau _K \alpha ) }{\frac {1-m-R(\pi )}{a(\pi )}- \frac {\psi g_2}{\pi } +\frac {\mu }{a(\pi )\pi } } \end{equation}

Using (13), it is equivalent to:

(30) \begin{equation} b =\frac {a(\pi )\pi [g -(\tau _L (1-\alpha )+\tau _K \alpha )] }{\pi (\delta -m )-(1-\tau _K)\alpha a(\pi )\pi - \psi g_2 a(\pi ) +\mu }= \frac {g -\tau _L (1-\alpha )-\tau _K \alpha }{X(\pi )-\frac {1-\delta }{a(\pi )}-(1-\tau _K)\alpha }\equiv B_1(\pi ) \end{equation}

with $X(\pi )\equiv \frac {1-m}{a(\pi )}- \frac {\psi g_2}{\pi } +\frac {\mu }{a(\pi )\pi }$ .

Moreover, (27) can be rewritten as:

(31) \begin{eqnarray} b&=&\frac {a(\pi )}{R(\pi )} \left [ \Sigma (1-\tau _L)(1-\alpha ) +\tau _L (1-\alpha )+\tau _K \alpha -g - X(\pi )\right ] \notag \\ &\equiv & B_2(\pi ) \end{eqnarray}

In the following, we assume:

Assumption 3. $\delta \geq m$ and $\epsilon _a(\pi )\lt -1$ .

The first part of the assumption implies that the rate of pollution absorption is lower than the depreciation rate of capital. This is consistent in our context, as pollution stock can refer to greenhouse gases whose some will remain in the atmosphere for thousands of years. The second part of the assumption implies that total factor productivity is elastic to pollution over capital ratio, illustrating an important vulnerability to climate change (see IPCC 2022). When pollution per unit of capital increases total factor productivity falls more than proportionally. It implies that $a(\pi )\pi$ is decreasing in $\pi$ . As a result, $X(\pi )$ is an increasing function of $\pi$ .

Example (continued): note that with the example defined in equation (19), $\epsilon _a(\pi )\lt -1$ implies that $\pi \gt g_1A_1/2$ .

This example illustrates the fact that $\epsilon _a(\pi )\lt -1$ could introduce a lower bound $\underline {\pi }\gt 0$ defined by $\epsilon _a(\underline {\pi })=-1$ such that $\epsilon _a(\pi )\lt -1$ for all $\pi \gt \underline {\pi }$ .

Under Assumptions2 and 3, the numerator of (30) is positive and the denominator is increasing in $\pi$ . We thus have $B_1^{\prime }(\pi )\lt 0$ . In addition, to ensure a positive debt along the balanced growth path, we restrict our attention to cases where:

(32) \begin{equation} X(\pi )-\frac {1-\delta }{a(\pi )}\gt (1-\tau _K)\alpha \end{equation}

New debt emissions should be higher than the cost of existing debt. Indeed, using (27) and (28), $X(\pi )=\gamma /a(\pi )$ . It implies that inequality (32) is equivalent to $\gamma \gt R(\pi )$ . Using (30), we have $b[\gamma -R(\pi )]=a(\pi )[g -\tau _L (1-\alpha )-\tau _K \alpha ]$ . Since we assume that the government budget is characterized by a primary deficit, a BGP should be characterized by a growth factor larger than the interest factor.

Since the left-hand side of inequality (32) is increasing in $\pi$ , there exists $\pi _1\gt 0$ such that $X(\pi _1)=\frac {1-\delta }{a(\pi _1)}+(1-\tau _K)\alpha$ if there is a value $\widetilde {\pi }\gt 0$ such that $X(\widetilde {\pi })\lt \frac {1-\delta }{a(\widetilde {\pi })}+(1-\tau _K)\alpha$ . Then, inequality (32) is satisfied for all $\pi \gt \pi _1$ , with $\displaystyle \lim _{\pi \to \pi _1^+}B_1(\pi )=+\infty$ . Note that the existence of $\pi _1\gt 0$ can be compatible with Assumption2.

Lemma 1. Under Assumptions 1 3 , assume that there exists a value $\widetilde {\pi }\gt 0$ such that $X(\widetilde {\pi })\lt \frac {1-\delta }{a(\widetilde {\pi })}+(1-\tau _K)\alpha$ . Then, there exists $\pi _1\gt 0$ such that $X(\pi )\gt \frac {1-\delta }{a(\pi ) }+(1-\tau _K)\alpha$ for $\pi \gt \pi _1$ .

Our example illustrates that this lemma is satisfied for a non-empty set of parameters.

Example (continued): using the example defined in equation (19), we illustrate the existence of the bound $\pi _1$ . Inequality (32) writes $F(\pi )\gt (1-\tau _K)\alpha$ , with:

(33) \begin{equation} F(\pi )\equiv \frac {\pi [g_1A_1(\delta -m)+\psi g_2]+g_1A_1(\mu -\psi g_2)}{\pi (g_1A_1-\pi )} \end{equation}

Since $F(\pi )$ is an increasing function and $F(g_1A_1)=+\infty$ , there exists $\pi _1\in (g_1A_1/2, g_1 A_1)$ if $F(g_1A_1/2)\lt (1-\tau _K)\alpha$ . This happens if $\delta -m \lt A_1(1-\tau _K)\alpha /2$ and

(34) \begin{eqnarray} g_1\gt \frac {4 \mu -2\psi g_2 A_1}{A_1 [(1-\tau _K)\alpha A_1-2 (\delta -m)]}\equiv \underline {g_{1}} \left ( g_{2}\right ) \end{eqnarray}

Using (31), a steady state with positive debt per unit of capital should satisfy:

(35) \begin{equation} \Sigma (1-\tau _L)(1-\alpha ) \gt g-\tau _L (1-\alpha )-\tau _K \alpha +X(\pi ) \end{equation}

Since the left-hand side of this inequality is constant and the right-hand side is increasing in $\pi$ , there exists $\pi _2\gt 0$ such that inequality (35) is satisfied for all $\pi \lt \pi _2$ , with $B_2(\pi _2)=0$ . In this case, we also deduce that under Assumption3, we have $B_2^{\prime }(\pi )\lt 0$ . For this, we need to have a value $\widehat {\pi }\gt 0$ such that $\Sigma (1-\tau _L)(1-\alpha ) \lt g-\tau _L (1-\alpha )-\tau _K \alpha +X(\widehat { \pi })$ .

Lemma 2. Under Assumptions 1 3 , assume that there exists a value $\widehat {\pi }\gt 0$ such that $\Sigma (1-\tau _L)(1-\alpha ) \lt g-\tau _L (1-\alpha )-\tau _K \alpha +X(\widehat {\pi })$ . The interval $(\pi _1, \pi _2)$ is non-empty if the following inequality is satisfied:

(36) \begin{equation} \Sigma (1-\tau _L)(1-\alpha ) \gt g-\tau _L (1-\alpha )-\tau _K \alpha +\frac { 1-\delta }{a(\pi _1)}+(1-\tau _K)\alpha \end{equation}

where $\pi _2$ is defined by $\Sigma (1-\tau _L)(1-\alpha ) =g-\tau _L (1-\alpha )-\tau _K \alpha +X(\pi _2)$ .

Example (continued) : in our example, let $\widehat {\pi }=g_1A_1$ be such that $a (\widehat {\pi })=0$ . In this case, we have $X(\widehat {\pi })=+\infty$ , which ensures the first inequality in the lemma.

Note that in our example, we have $\pi _1$ higher but arbitrarily close to $g_1A_1/2$ if $g_1$ tends to $\frac {4 \mu -2\psi g_2 A_1}{A_1 [(1-\tau _K)\alpha A_1-2 (\delta -m)]}$ . Therefore, inequality (36) is satisfied if $\Sigma (1-\tau _L)(1-\alpha ) \gt g-\tau _L (1-\alpha )-\tau _K \alpha +\frac { 2(1-\delta )}{A_1}+(1-\tau _K)\alpha$ . This last inequality is satisfied if $\tau _K$ and $A_1$ are high enough and the primary deficit is not too important. It proves the existence of $\pi _2$ and of a non-empty interval $(\pi _1, \pi _2)$ .

Since $B_1(\pi _1)=+\infty \gt B_2(\pi _1)$ and $B_1(\pi _2)\gt B_2(\pi _2)=0$ , the economy may be characterized by an even number (two) of steady states.

Proposition 1. Under Assumptions 1 3 , and inequality ( 36 ), there exists $\overline {g}\gt \tau _L (1-\alpha )+\tau _K \alpha$ such that for $g$ $\in$ ( $\tau _L (1-\alpha )+\tau _K \alpha, \overline {g}$ ), there are (at least) two BGPs, $(\pi _I, b_I)$ and $(\pi _{II}, b_{II})$ , with $\pi _I\lt \pi _{II}$ and $b_I \gt b_{II}$ .

Proof. See Appendix B.

A primary deficit ( $g\gt \tau _L (1-\alpha )+\tau _K \alpha$ ) ensures the possibility of having a positive stationary debt per unit of capital in our context where growth is higher than the interest factor. At the same time, if the environmental expenditure is too high ( $g\gt \overline {g}$ ), the primary deficit is too significant to observe stationarity in debt per unit of capital. The share of GDP devoted to environmental issues has to take an intermediate value to observe stationary solutions. In addition, condition (36) has to be satisfied.

Example (continued): using the example defined in equations (2) and (19), we illustrate what the conditions for existence of BGPs imply in terms of policy.

Using Lemmas1 and 2, we have conditions ensuring that the interval $ (\pi _{1},\pi _{2})$ , necessary for the existence of a steady-state $\pi _{i}$ , is non-empty. Adaptation expenditures must satisfy the constraint (34), $g_{1}\gt \underline {g_{1}}\left ( g_{2}\right )$ , where adaptation $g_1$ must not deviate significantly from the minimum value $\underline {g_{1}}(g_2)$ .

Thus, for any tax rates $\tau _{L}$ et $\tau _{K}$ , existence of a primary deficit and (34) determine the precise level of environmental policy expenditures $g_{1}$ and $g_{2}$ . Adaptation $g_{1}$ has to be close to $\underline {g_{1}}\left ( g_{2}\right )$ and mitigation $g_{2}$ adjusts so as to have a weak primary deficit. Using Assumption3, mitigation is also constrained by the condition imposing a maximum threshold $\mu /A_{1}\gt \psi g_{2}$ . Considering that $\mu$ and $A_1$ are close or equal, this last inequality is never binding.

We also need to ensure that savings are high enough to finance the current deficit, investment, and the cost of debt, i.e. inequality (36) holds:

(37) \begin{equation} \frac {\beta }{1+\beta }\left ( 1-\tau _{L}\right ) \left ( 1-\alpha \right ) \gt g_{1}+g_{2}-(\tau _{L}(1-\alpha )+\tau _{K}\alpha )+2\frac {\left ( 1-\delta \right ) }{A_{1}}+\alpha \left ( 1-\tau _{K}\right ) \end{equation}

To understand how considering these different conditions create constraints on the policy instruments, let us define $\kappa \gt 0$ small enough as being the primary deficit: $\kappa =g_{1}+g_{2}-\left [ \tau _{L}\left ( 1-\alpha \right ) +\tau _{K}\alpha \right ]$ . Substituting this expression in inequality (37), we obtain:

\begin{equation*} \tau _{L}\lt \frac {\alpha }{1-\alpha }\frac {1+\beta }{\beta }\tau _{K}+1-\frac {1}{1-\alpha }\frac {1+\beta }{\beta }\left ( \kappa +\alpha +2\frac {1-\delta }{A_{1}}\right ) \end{equation*}

We show that these conditions impose interdependencies between the tax rates $\tau _L$ and $\tau _K$ , and the size of the primary deficit $\kappa$ . If the tax rate on labor income is positive, it requires a sufficiently high tax rate on capital income. Moreover, the higher the primary deficit is, the higher the capital tax rate. This positive relationship allows to keep the debt burden low enough such as the debt over capital ratio can be stationnary.

The BGP characterized by the lowest pollution to capital ratio ( $\pi _{I}$ ) has the highest level of debt over capital ( $b_{I}$ ), while the one with the highest pollution-to-capital ratio ( $\pi _{II}$ ) is also defined by the lowest debt per unit of capital ( $b_{II}$ ). We can note that the decreasing relationship we observe between $\pi$ and $b$ is ensured by the sufficient TFP vulnerability ( $\epsilon _a(\pi )\lt -1$ ). This property is specific to our analytical framework and is explained in detail in the following section.

4.2 Balanced growth and the role of TFP sensitivity to pollution

We want to clearly understand the links between pollution to capital, debt to capital, and growth that come from the comparison of the two BGPs.

First, we turn our attention to the growth factor $\gamma$ . From (28), we see that it is a declining function of the debt per capital ratio $b$ , through a priori a usual crowding-out effect of debt on investment. Moreover, as pollution generates negative external effects on production and therefore also on savings, the growth factor also depends negatively on the pollution per capital ratio $\pi$ . In our context of TFP vulnerability ( $\epsilon _a(\pi )\lt -1$ ), a BGP with a higher level of $\pi$ is characterized by a lower debt per capital ratio $b$ . Therefore, at the BGPs, the relationship between $\gamma$ and $\pi$ (or $b$ ) seems ambiguous.

To avoid this ambiguity, we exploit the fact that the growth of capital is equal to the growth of the pollution stock. Then, using (27) and (28), the growth factor can be expressed as a function of $\pi$ :

(38) \begin{equation} \gamma = 1-m+\frac {\mu -\psi g_2 a(\pi )}{\pi } \end{equation}

Hence, we note that the growth factor is higher than 1, i.e. growth is positive, as soon as $m$ is low enough. The growth factor increases with the pollution flow but decreases with the current pollution stock. Indeed, since mitigation decreases with pollution through its effect on TFP, $\pi$ has two opposite effects on growth, a positive one through pollution flows and a negative one through the pollution stock.

Therefore, growth is a decreasing function of pollution over capital ratio (and hence an increasing function of $b$ ) if and only if the elasticity of TFP with respect to $\pi$ satisfies:

(39) \begin{equation} -\epsilon _a(\pi )\lt \frac {\mu -\psi g_2 a(\pi )}{\psi g_2 a(\pi )} \end{equation}

which may be in accordance with Assumptions2 and 3. We thus deduce that:

Corollary 1. Under Assumptions 1 3 , inequality ( 36 ), and $\overline {g}\gt g\gt \tau _L (1-\alpha )+\tau _K \alpha$ , we have:

  1. 1. $\gamma _I\gt \gamma _{II}$ iif $-\epsilon _a(\pi )\lt \frac {\mu -\psi g_2 a(\pi ) }{\psi g_2 a(\pi )}$ for all $\pi \in (\pi _1, \pi _2)$ ;

  2. 2. $\gamma _I\lt \gamma _{II}$ iif $-\epsilon _a(\pi )\gt \frac {\mu -\psi g_2 a(\pi ) }{\psi g_2 a(\pi )}$ for all $\pi \in (\pi _1, \pi _2)$ .

In case 1, the growth rate decreases with $\pi$ , so it is lower in the state with low debt and high pollution ( $\pi _{II}$ , $b_{II}$ ). This configuration is characterized by a not excessive TFP vulnerability to pollution. When the production is not too sensitive to pollution through the TFP, a higher level of $\pi$ means a lower pollution flow over pollution stock, which implies lower growth.

In case 2, the high TFP vulnerability explains that growth is higher in the state ( $\pi _{II}$ , $b_{II}$ ). Indeed, a higher level of pollution over capital implies a strong increase in the pollution flow because of the decrease in public mitigation. Then, the pollution flow over the pollution stock increases, which implies higher growth.

This result is important as it reveals that as long as we consider a not excessive TFP vulnerability to pollution (case 1 of Corollary1), all other things being equal, a BGP characterized by a lower pollution level per unit of capital is associated with higher capital growth. Recall that this BGP also has a higher level of debt over capital. A direct implication of Corollary1 is that when TFP vulnerability is not excessive, a BGP with higher growth means a BGP with higher debt over capital. In contrast, when TFP vulnerability is very high, a BGP with lower growth means a BGP with higher debt over capital. Therefore, there is a crowding-in effect of debt on growth in the first case and a crowding-out effect in the second one.

Using (26) and (28), the intertemporal budget constraint evaluated at a BGP can be written:

(40) \begin{equation} b [\gamma -R(\pi )]=a(\pi )[g-\tau _L (1-\alpha )-\tau _K \alpha ] \end{equation}

When pollution over capital is low, the TFP and, therefore, the interest factor are high. This means that both the primary deficit and debt services are high. Debt over capital is high, even if growth is higher than at the steady state with a higher ratio of pollution over capital (case 1 of Corollary1). This explains that higher debt can be compatible with higher growth. This is an interesting result regarding the macroeconomic literature that mainly finds that public debt usually has a crowding-out effect on growth (see the seminal contribution by Diamond, Reference Diamond1965), except in the presence of some financial imperfections (see Woodford, Reference Woodford1990). Using (40), if the TFP is constant and, therefore, the interest factor too, we immediately deduce that public debt over capital and growth are inversely related. Public debt always has a crowding-out effect on growth.

Now, we want to understand precisely why a high TFP vulnerability to pollution is the source of a negative link between $\pi$ and $b$ when we compare both BGPs. We start by examining the extreme case in which TFP is not sensitive to pollution damage (i.e. $\epsilon _a(\pi )=0$ ). It implies that $a(\pi )=a$ and $R(\pi )=R$ are constant. Equations (30) and (31) thus become:

(41) \begin{align} B_1(\pi ) &=\frac {a[g -(\tau _L (1-\alpha )+\tau _K \alpha )] }{\delta -m -(1-\tau _K)\alpha a +\frac {\mu - \psi g_2 a}{\pi } } \end{align}
(42) \begin{align} B_2(\pi )&=\frac {a}{R} \left [ \Sigma (1-\tau _L)(1-\alpha ) +\tau _L (1-\alpha )+\tau _K \alpha -g - \frac {1-m}{a}- \frac {1}{\pi }\left (\frac {\mu }{a}-\psi g_2\right )\right ] \\[8pt]\nonumber\end{align}

with $B_1^{\prime }(\pi )\gt 0$ and $B_2^{\prime }(\pi )\gt 0$ . In that case, if there still exist several steady states, the one with the highest level of pollution over capital will be also characterized by the highest level of debt over capital. Comparing equations (41) and (42) with equations (30) and (31) provides insights into the differences that arise.

Equation (41) represents debt over capital as the ratio of the primary deficit over the new debt emission, which increases with growth, minus the cost of debt services measured by the interest factor. As previously mentioned, the growth of capital is equal to the growth of pollution at a BGP, which explains that it is decreasing in the pollution stock over capital and hence that there is a positive relationship between $\pi$ and $b$ . When productivity is negatively affected by pollution over capital, several adding effects may imply a reversal of this link. These effects can be perceived by using equation (30). Production being affected negatively by $\pi$ , a higher $\pi$ implies a lower primary deficit $a(\pi )[g-\tau _L(1-\alpha )-\tau _K \alpha ]$ . In addition, the cost of debt reimbursement goes down with productivity loss. These adding effects, due to endogenous productivity, are important when vulnerability to pollution is high enough. In that case, we have a negative relationship between debt over capital and pollution over capital. Equation (42) comes from the equilibrium on the asset market taking into account that capital growth is equal to pollution growth and considering the government budget constraint. Debt per unit of capital is equal to the difference, discounted by the interest factor, between savings and the sum of the primary deficit and the increase in capital (which is here equal to pollution growth). When productivity is constant, the only effect of $\pi$ on $b$ is positive and is due to its negative effect on pollution growth, as previously mentioned. When the TFP decreases with pollution over capital, several adding effects may overturn the link between debt over capital and pollution over capital. These effects can be perceived by using equation (31). First, a lower TFP decreases income discounted by the interest factor. Second, growth per unit of TFP increases with $\pi$ , which has a crowding-out effet on debt.

To summarize, TFP vulnerability to pollution implies that higher pollution to capital reduces product so that the debt to capital ratio reduces too because less deficit has to be financed and less saving is available to buy public debt. Turning to the analysis of the dynamics is now essential to determine toward which equilibria the economy will converge.

5. Dynamics, endogenous tipping zone and sustainability

The first main question we ask is whether the economy might converge to a BGP. Only in such a case, the economy will be sustainable in the long run. Otherwise, either the economy will collapse, or pollution over capital will follow an explosive dynamic path. We will especially identify a zone in terms of initial conditions, that we call the endogenous tipping zone, such that the economy will not be sustainable.

Proposition1 states that two BGPs with positive debt, capital, and pollution may coexist. We analyze the dynamics in this interesting case. The question is to know toward which BGP the economy will converge. We aim to identify the conditions for a sustainable or explosive and unsustainable dynamic path. The objective is to highlight the respective roles of fiscal instruments, TFP vulnerability to pollution, and initial conditions on pollution debt and capital.

5.1 Stability and sustainability

The dynamics are driven by equations (24) and (25). They can alternatively be driven by a combination of (24) and (25), and equation (25).

Let $\Phi _t\equiv \frac {b_t}{\pi _t}=\frac {B_t}{P_t}$ . Using (24) and (25), we obtain:

(43) \begin{equation} \Phi _{t+1} = \frac {\frac {R(\pi _t)}{a(\pi _t)} \Phi _t+\frac {g -\tau _L (1-\alpha )-\tau _K \alpha }{\pi _t}}{X(\pi _t)} \end{equation}

Then, $\Phi _{t+1}\geqslant \Phi _t$ is equivalent to:

(44) \begin{equation} \Phi _t \leqslant \frac {g -\tau _L (1-\alpha )-\tau _K \alpha }{\pi _t [X(\pi _t)- \frac {1-\delta }{a(\pi _t)}-(1-\tau _K)\alpha ]} \equiv B_3(\pi _t) \end{equation}

with $B_3(\pi _t)=B_1(\pi _t)/\pi _t$ , $B_3(\pi _1)=\infty$ , $B_3(\pi _2)\gt 0$ and $B_3^{\prime }(\pi _t)\lt 0$ .

Using (25) and (31), $\pi _{t+1}\geqslant \pi _t$ rewrites $b_t\geqslant B_2 (\pi _t)$ . This is equivalent to $\Phi _t\geqslant B_2 (\pi _t)/\pi _t \equiv B_4 (\pi _t)$ , given by:

(45) \begin{eqnarray} B_4(\pi _t) &=&\frac {a(\pi _t)}{R(\pi _t)\pi _t} \left [ \Sigma (\pi _t)(1-\tau _L)(1-\alpha ) +\tau _L (1-\alpha )+\tau _K \alpha -g - X(\pi _t)\right ] \end{eqnarray}

with $B_4^{\prime }(\pi _t)\lt 0$ , $B_4 (\pi _1)\gt 0$ and $B_4 (\pi _2)=0$ .

We can draw a phase diagram using these different ingredients and the results of the previous section. The stationary values of debt over pollution are $\Phi _I= \frac {b_I}{\pi _I}$ and $\Phi _{II}= \frac {b_{II}}{\pi _{II}}$ . Since $B_3(\pi _t)$ and $B_4 (\pi _t)$ are both decreasing and $\pi _I\lt \pi _{II}$ , we deduce that $\Phi _I\gt \Phi _{II}$ . We further note that $B_3^{\prime }(\pi _I)\lt B_4^{\prime }(\pi _I)$ , while $B_3^{\prime }(\pi _{II})\gt B_4^{\prime }(\pi _{II})$ .

A qualitative picture of the dynamics is represented in Figure 2. We conjecture that the steady state $(\pi _I, \Phi _I)$ is stable, whereas the steady state $(\pi _{II}, \Phi _{II})$ is a saddle. Since the two dynamic variables $\pi _t$ and $\Phi _t$ are predetermined, a saddle is generically unstable. We will now confirm this conjecture by the analysis of local dynamics.

Figure 2. Dynamics with sustainability.

Proposition 2. Under Assumptions 1 3 , inequality ( 36 ), and $g$ $\in$ ( $\tau _L (1-\alpha )+\tau _K \alpha, \overline {g}$ ), we have the following:

  1. 1. The steady state $(\pi _{II}, b_{II}, \Phi _{II})$ is a saddle;

  2. 2. The steady state $(\pi _{I}, b_{I}, \Phi _{I})$ is locally stable if $\tau _L$ is high enough and $\epsilon _a (\pi )$ is not too negative.

Proof. See Appendix C.

This proposition shows that the economy never converges to the steady state with a high level of pollution over capital. As pollution, debt, and capital are all predetermined, the saddle ( $\pi _{II}, b_{II}$ ) is never achieved and therefore delimits a pollution trap, as we will discuss later. Indeed, when initial conditions are characterized by too high levels of pollution and debt with respect to capital, the economy cannot converge to a long-run BGP. Both $b_t$ and $\pi _t$ increase across time. The too high level of pollution over capital implies a too low TFP and GDP growth to be compatible with the convergence to a stable and sustainable BGP. When pollution over capital is not too high, the economy would converge to the stable BGP characterized by the low level of pollution over capital and the highest level of debt over capital or pollution $(\pi _I, b_I)$ , as represented in Figure 2. Such a dynamic path could experience oscillations converging to the steady state.Footnote 11 Depending on the level of pollution over capital, the convergence to this steady state does not require a so low debt relative to capital. The stability of this BGP requires a not too strong TFP vulnerability to pollution and a high tax rate on labor income. If the first condition is not fulfilled, any increase in pollution over capital implies a strong decrease in GDP and, therefore, investment in future capital due to the high TFP vulnerability to pollution. Future pollution over capital increases even more, generating an explosive dynamic path (see for instance equation (25)). Under the last condition, savings are low enough to prevent an unsustainable buildup of debt and capital.

5.2 Endogenous tipping zone (ETZ)

Since the steady state $(\pi _{II}, b_{II}, \Phi _{II})$ is a saddle, it has one stable and one unstable manifold. By inspection of Figure 2, we observe that the stable manifold of this steady state delineates a zone such that when the economy is on the right side of this manifold, pollution will be explosive. We call this zone an Endogenous Tipping Zone (ETZ), because it is endogenously defined by the dynamic behavior of the economy. In the following proposition, we show that:

Proposition 3. Under Assumptions 1 3 , inequality ( 36 ), and $g\in (\tau _L (1-\alpha )+\tau _K \alpha, \overline {g}$ ), the stable manifold of the steady state $(\pi _{II}, b_{II}, \Phi _{II})$ has a negative slope at least in the neighborhood of the steady state, while the unstable manifold has a negative slope but higher than the stable one.

Proof. See Appendix D.

Figure 3. Endogenous tipping zone (ETZ).

The stable manifold $(SM)$ of the steady state $(\pi _{II}, b_{II}, \Phi _{II})$ is clearly represented in Figure 3. Since the two dynamic variables are predetermined, the economy is generically on the left or the right side of this decreasing curve. If the initial conditions are such that the economy is on the left side of $(SM)$ , the dynamics could be characterized by convergence to the stable steady state with an increase in the long run of debt over pollution and capital and a decrease in pollution per unit of capital.

If the initial conditions are such that we are on the right side of $(SM)$ , pollution per unit of capital will increase a priori indefinitely. Since the unstable manifold is negatively sloped, pollution will also increase with respect to debt. The economy is in the ETZ if the stock of pollution is sufficiently high with respect to capital. Interestingly, at least around the BGP, the ETZ is delimited by a negatively slopped relationship between $\Phi _t$ and $\pi _t$ . This means that the economy will experience explosive paths for lower pollution levels over capital when the debt relative to pollution is higher. When an economy is already burdened with high debt levels, it may struggle to allocate additional resources toward pollution reduction efforts and adaptation. A vicious cycle is triggered, where a high level of debt does not translate into significant spending on pollution control. Instead, it illustrates an inability to adapt and mitigate adequately, thereby increasing damages and making it more difficult to stabilize debt and pollution per unit of capital. In contrast, the lower the debt over pollution, the higher the level of pollution over capital to have an explosive path. This means that lower public debt and higher capital allow to reach a sustained growing economy more easily.

Considering the damage of pollution stock on production, we provide a theoretical mechanism explaining why, for a given pollution-to-capital ratio, the likelihood of observing explosive paths for pollution is higher in countries with higher debt. Such a detrimental situation is pointed out in the literature, but as we can remark in Zenios (Reference Zenios2024), the papers that attempt to formalize it focus mainly on the impact of climate risks on public spending; he highlights that these risks increase the cost of public debt, making public finances even more vulnerable. We obtain a similar conclusion of debt vulnerability considering a purely deterministic context.

The importance of our result is all the greater as public debt-to-GDP ratios have risen for decades and have reached record levels in a significant number of both developed and developing countries(see Figure 1 and WorldBank, 2023). This trend, combined with the heightened vulnerability to environmental issues, emphasizes the need to propose adapted policy tools to avoid a vicious circle of pollution over debt.

5.3 Unsustainability

The system can also be completely unsustainable. This happens if no steady state is stable, i.e. if the equilibrium $(\pi _{I}, b_{I}, \Phi _{I})$ is unstable. In such a situation, pollution and debt over capital will follow an explosive dynamic paths. As we will see, it can occur if the TFP vulnerability to pollution is high and the tax rates are low. The next proposition provides sufficient conditions for such an undesirable configuration:

Figure 4. Dynamics with unsustainability.

Proposition 4. Under Assumptions 1 3 , inequality ( 36 ), and $g$ $\in$ ( $\tau _L (1-\alpha )+\tau _K \alpha, \overline {g}$ ), we have the following:

  1. 1. The steady state $(\pi _{II}, b_{II}, \Phi _{II})$ is a saddle;

  2. 2. The steady state $(\pi _{I}, b_{I}, \Phi _{I})$ is unstable if $\tau _L$ is low enough, $\epsilon _a(\pi )\lt -2$ , $1+m\gt 2\delta$ and $1-\tau _K\gt \mu /[2\alpha a(\pi _2)\pi _2]$ .

Proof. See Appendix E.

This proposition gives sufficient conditions to have all steady states saddle or unstable. It particularly requires sufficiently low tax rates on capital and labor incomes and a strong TFP vulnerability to pollution. Of course, the BGPs are affected by the level of the tax rates, which means that the stationary values of $\pi _i$ and $b_i$ are not similar in Proposition4 and in Proposition2. This explains that we do not conduct an analysis of bifurcations, but in both configurations, two BGPs coexist. When the tax rates are low, the primary deficit is important, as is the cost of debt. Both effects directly deteriorate public finance. The positive impact of low tax rates on aggregate savings is not sufficient to compensate. Low taxation creates conditions promoting unsustainable debt levels and hence hinders the government’s ability to tackle pollution issues. Moreover, when TFP is highly vulnerable to pollution damage, the economy can never converge to a long-run BGP. As we have already seen, a small increase in pollution over capital implies a strong decrease in production and capital, which implies a larger future increase in pollution over capital. Either pollution and debt will go to infinity, characterizing a vicious circle of debt and pollution or the economy will collapse. There may exist a dynamic trajectory that diverges from the BGP with low pollution over capital to converge to the BGP with high pollution over capital. However, since the variables are predetermined, the economy will not generically experience such a dynamic path.

To summarize this section, high levels of tax rates are key ingredients to rule out an explosive accumulation of pollution and debt. However, the economy may be unsustainable for technological reasons, i.e. a high TFP vulnerability to pollution. Finally, high initial levels of debt and pollution with respect to capital promote instability of the dynamic path. Note that a high initial debt over capital is not a priori a source for unsustainability since it may reinforce the possibility of converging to the stable steady state with low pollution over capital. It will depend on the level of pollution over capital.

6. Policy implications

We are now interested in a configuration of possible sustainability. This means that we consider the case in which the BGP $(\pi _{I}, b_{I}, \Phi _{I})$ is locally stable. In this encouraging situation, we examine which public policy can improve welfare in the long run. Instead of focusing on first best policies,Footnote 14 we examine policies which improve welfare in the long run or after a given number of periods. Such analysis can be performed by analyzing comparative statics at the stable BGP. Starting at or close to this BGP, the economy will converge to it after the policy change. Therefore, analyzing the change of welfare at such a BGP is appropriate to evaluate the change of welfare after a certain number of periods. For this aim, we precisely study the effect of policy variables that allow the management of environmental adaptation and mitigation, $g_1$ and $g_2$ , and the fiscal revenue, $\tau _L$ and $\tau _K$ , on the welfare at the BGP $i=I$ .

We start by evaluating the welfare at a BGP $i=\{I, II\}$ . The consumptions are given by:

(46) \begin{align} c_{it}&=(1-\Sigma )(1-\tau _L)w_{it} \end{align}
(47) \begin{align} d_{it+1}&=R(\pi _i)\Sigma (1-\tau _L)w_{it} \\[9pt]\nonumber\end{align}

with $w_{it}=(1-\alpha )a(\pi _i)\gamma _i^t K_0$ . Substituting these two consumptions into the utility function (10), we get:

\begin{equation*} U(c_{it}, d_{it+1})=\ln (1-\Sigma )+\beta \ln \Sigma R(\pi _i) +(1+\beta )\ln (1-\tau _L)(1-\alpha )a(\pi _i)\gamma _i^t K_0 \end{equation*}

Therefore, the main factor determining the utility evaluated at a BGP is the growth factor $\gamma _i$ . A higher growth factor has a dominant positive effect on the welfare. Since $R(\pi )$ and $a(\pi )$ are decreasing in $\pi$ , a lower pollution over capital also positively affects welfare. Based on (38), we see that policy variables affect $\gamma _i$ through their impact on $\pi _i$ . In addition, we note that the policy variable $g_2$ has also a direct negative effect on the growth rate.

One necessary step to determine how policy variables modify welfare is to examine their effects on ( $\pi _i, b_i$ ). For this aim, we redefine all relevant equations as functions of ( $\pi _i, b_i$ ) and policy instruments ( $ \tau _L, \tau _K, g_1, g_2$ ). Productivity can thus be defined as a function of $g_1$ and $\pi _i$ :

\begin{equation*} A(G_{1}/P)=A[\theta (\pi _i, g_1)]\equiv \Theta (\pi _i, g_1) \end{equation*}

where $\theta (\pi _i, g_1)$ is defined by (18), and is increasing in $g_1$ and decreasing in $\pi _i$ . We deduce the sign of the two following derivatives: $\Theta _{g_1}(\pi _i, g_1)\gt 0$ , $\Theta _{\pi _i}(\pi _i, g_1)\lt 0$ . The function $\Theta _{g_1}(\pi _i, g_1)$ captures how productivity responds to a variation in mitigation.

We present the influence of fiscal policy, via $\tau _L$ and $\tau _K$ , and environmental policy, through $g_1$ and $g_2$ , focusing on the potentially stable and sustainable steady state $(b_I, \pi _{I})$ (see Proposition2).

Lemma 3. Under Assumptions 1 3 and inequality ( 36 ), we have the following properties:

  • An increase in $\tau _L$ and/or $\tau _K$ increases $b_I$ and reduces $\pi _{I}$ ;

  • An increase in $g_1$ reduces $\pi _{I}$ if $\Theta _{g_1}(\pi, g_1)$ is high enough, $\tau _L$ is high and $g-\tau _L (1-\alpha )-\tau _K\alpha$ is not too small;

  • There exists $\bar {\psi }\gt 0$ , such that if $\psi \lt \bar {\psi }$ , an increase in $g_2$ reduces $b_I$ and increases $\pi _{I}$ .

Proof. See Appendix F.

Lemma3 emphasizes that at the low pollution BGP ( $b_I,\pi _{I}$ ), more stringent taxation reduces pollution over capital. This is because more fiscal revenue reduces the primary deficit and hence promotes capital investment. However, even though the government earns more revenue from taxation, the reduction in pollution to capital ratio goes with an increase in the debt ratio. Higher fiscal pressure leading to a higher debt over capital can appear counter-intuitive. It is driven by two key assumptions: productivity reacts to the fall in pollution per unit of capital and there is a primary deficit ( $g\gt \tau _L(1-\alpha )+\tau _K\alpha$ ). In this context, the increase in the debt ratio is a positive side effect of productivity gains, which is facilitated by the reduction in pollution damage. Note that even if the qualitative effect of tax on capital and labor is similar, it masks different mechanisms. Labor tax reduces labor income while capital taxation reduces the cost of capital and debt.

Concerning environmental policy variables, an increase in the share of the public budget allocated to adaptation leads to the following consequences: the rise in $g_1$ generates competing effects on productivity and hence on pollution to capital ratio. First, it puts pressure on the primary deficit, leading to a reduction in capital accumulation due to a crowding-out effect. This increases pollution per unit of capital ( $\pi$ ) and reduces productivity. Second, the increase in $g_1$ has a direct positive effect on productivity, as a larger share of public spending is allocated to adaptation. When productivity is sufficiently sensitive to $g_1$ ( $\Theta _{g_1}(\pi, g_1)$ high enough), the negative effect driven by the crowding-out is surpassed by the direct positive effect of $g_1$ . Pollution per unit of capital ( $\pi$ ) goes down while productivity goes up.

Concerning the effect of mitigation, the efficiency of public spending to reduce pollution flow (captured by $\psi$ ) is crucial, because it determines how pollution and debt respond to an increase in $g_2$ at the BGP. More precisely, we can observe a backfire effect of mitigation policy in the sense that it can increase pollution per unit of capital and reduce debt. Indeed, as for $g_1$ , the rise in $g_2$ puts pressure on the primary deficit, leading to a reduction in capital accumulation due to a crowding-out effect. This, in turn, leads to an increase in pollution per unit of capital ( $\pi$ ), resulting in decreased productivity $a(\pi )$ . At the same time, the increase in $g_2$ has a direct negative effect on pollution, whose importance depends on $\psi$ , as a larger share of public spending is allocated to mitigation. However, when mitigation is not sufficiently efficient to reduce pollution flow, this direct effect of an increase in $g_2$ on pollution is low. The negative feedback effects of an increase in mitigation on capital accumulation, and hence on productivity, surpasses the initial effect of the policy.

Now, we can identify welfare-improving policy scenarios along the BGP $i=I$ , i.e. policy scenarios that increase the growth factor $\gamma _I$ and decrease pollution over capital. In addition, using Corollary1, we recall that as $-\epsilon _a(\pi )$ is lower than $\frac {\mu -\psi g_2 a(\pi )}{\psi g_2 a(\pi )}$ , the growth rate evaluated at a BGP decreases with pollution over capital. Using this relationship, the one between $\pi _I$ and policy instruments presented in Lemma3, and equation (38), we have the following result:

Proposition 5. Under Assumptions 1 3 , inequality ( 36 ) and $-\epsilon _a(\pi )\lt \frac {\mu -\psi g_2 a(\pi )}{\psi g_2 a(\pi )}$ , the welfare along the BGP $i=I$

  • increases with taxation $\tau _L$ and $\tau _K$

  • increases with the share of budget allocated to adaptation $g_1$ if $\Theta _{g_1}(\pi, g_1)$ is high enough, $\tau _L$ is high and $g-\tau _L (1-\alpha )-\tau _K\alpha$ is not too small

  • decreases with the share of budget allocated to mitigation $g_2$ if $\psi$ is sufficiently low

Proof. See Appendix G.

When the TFP vulnerability to pollution is still not too high, a more stringent fiscal policy is a means to increase welfare as it allows for a reduction in the pollution-to-capital ratio. This is not so evident for environmental instruments because their financing costs can heavily burden the primary budget. This creates a crowding-out effect that outweighs the benefits when the effectiveness of these policies in achieving their main objectives (reducing pollution or increasing adaptation) is limited. However, when TFP reacts sufficiently to adaptation expenditures, this spending can improve growth and welfare despite the fact the primary deficit enlarges.

To illustrate in particular that the condition $\Theta _{g_1}(\pi, g_1)$ high can be in accordance with an elasticity of productivity with respect to pollution ratio satisfying $-\epsilon _a(\pi )\lt \frac {\mu -\psi g_2 a(\pi )}{\psi g_2 a(\pi )}$ , let us come back on our example.

Example (continued): in our example, $-\epsilon _a(\pi )=\frac {\pi /g_1}{A_1-\pi /g_1}=\frac {A_1}{a(\pi )}-1$ . Thus, condition (39) is equivalent to $\mu \gt A_1\psi g_2$ , which holds under Assumption2. This means that $-\epsilon _a(\pi )\lt \frac {\mu -\psi g_2 a(\pi )}{\psi g_2 a(\pi )}$ . Moreover, in our example, we have $\Theta _{g_1}(\pi, g_1)=\pi /g_1^2$ which is high if $g_1$ is low. In such a case, an increase in $g_1$ can improve welfare, as shown in Proposition5 .

Finally, according to Assumption2, mitigation is not highly efficient. Consequently, relying solely on this instrument is not a viable option for increasing welfare, as it does not directly dampen the fall in productivity caused by global pollution.

Showing that an increase in a tax (on labor or capital) improves welfare on a stable BGP means that, starting near this BGP before the policy change, the economy converges to the modified BGP. After a certain number of periods, the generations will benefit from the welfare increase because the economy is not too far from the new BGP. Of course, when the increase in tax (on labor or capital) is implemented, this is not beneficial to the current and first generations after the policy change takes place. This is a well-known result, shown in a different context by Fodha et al. (Reference Fodha, Seegmuller and Yamagami2018). Indeed, by inspection of equations (46) and (47), we easily see that an increase in the tax rates has a direct negative effect on at least one consumption, while the increase in growth and the reduction of pollution over capital will take some times to have a significant effect. However, after a certain number of periods, the economy will converge and become closer to the BGP with a higher level of growth and a lower pollution over capital. Hence, after a certain time, the welfare of all the successive generations will be higher. This does not imply that intergenerational equity issues are absent during the convergence process because of potential oscillations in the dynamics. Nonetheless, the cost of these fluctuations can be viewed as relatively small compared with the positive effect on welfare associated with the new BGP.

Two last remarks are at stake. First, we have in mind to consider policies that do not change the stability conditions of the BGP with the lowest level of pollution over capital. This means that we do not consider too drastic changes, although we can also note that increasing the tax rates promotes stability. Second, when the policy reduces pollution over capital at the lowest BGP $(b_I, \pi _{I})$ , it increases pollution over capital at the saddle BGP $(b_{II}, \pi _{II})$ . This suggests that the ETZ will occur for higher levels of pollution over capital. Therefore, this zone is reducing and the convergence to the stable BGP $(b_I, \pi _{I})$ could take place for a larger range of initial conditions.

7. Conclusion

This paper examines the complex interplay between public debt, environmental quality, and economic growth, particularly pertinent in the context of growing public debt and global pollution concerns. We address these issues within an endogenous growth framework that incorporates public debt dynamics, adaptation and mitigation spending, and feedback effects of pollution on productivity. We identify two balanced growth paths characterized by varying levels of pollution and debt relative to capital. Depending on fiscal policy, initial conditions, and the responsiveness of productivity to pollution and adaptation, the economy either converges to a sustainable BGP, collapses, or experiences perpetual increases in debt and pollution. Unsustainable debt and pollution arise particularly when pollution-induced damage severely impacts productivity, underscoring the importance of policy interventions for environmental and fiscal stability.

Our results suggest that technological efficiency is fundamental: improving technology to adapt to the impacts of climate change is important for welfare at a sustainable path. Government should start tackling global emissions (mitigation and adaptation) before a pollution threshold is reached; crossing this threshold, on the other hand, pushes the economy into a tipping zone of unsustainability. Finally, since environmental policies can be financed through both public debt and income tax revenues, we find conditions under which fiscal policy is a key element in reconciling debt sustainability and environmental sustainability.

Our findings suggest that increasing taxation, particularly in contexts where productivity is highly sensitive to pollution damage, can enhance welfare along stable growth paths. This emphasizes the crucial role of efficient environmental policy instruments in promoting adaptation.

Acknowledgements

We thank two anonymous referees for their helpful suggestions. We also thank participants in the conferences LORDE 2024 (Grenoble, March 2024), 28th ICMAIF (Crete, May 2024), Macroeconomic Green Dynamics (Urbino, June 2024), ASSET 2024 (Venice, November 2024) for their remarks and comments. All remaining errors are own.

Funding statement

Thomas Seegmuller acknowledges financial support from the French government under the France 2030 investment plan managed by the French National Research Agency Grant ANR-17-EURE-0020, and by the Excellence Initiative of Aix-Marseille University - A*MIDEX.

Appendix.

A. The relationship between $G_{1t}/P_t$ and $\pi _t$

Let $x_t=G_{1t}/P_t$ . Then, equation (18) writes:

\begin{equation*} x_t=\frac {g_1}{\pi _t}A(x_t) \end{equation*}

Differentiating this equation with respect to $x_t$ and $\pi _t$ , we obtain:

\begin{equation*} \frac {dx_t}{d\pi _t}=\frac {g_1 A(x_t)/\pi _t^2}{\frac {g_1}{\pi _t}A^{\prime }(x_t)-1} \end{equation*}

Since $\frac {g_1}{\pi _t}=\frac {x_t}{A(x_t)}$ , we finally have:

\begin{equation*} \frac {dx_t}{d\pi _t}=\frac {x_t/\pi _t}{\frac {x_t A^{\prime }(x_t)}{A(x_t)}-1}\lt 0 \end{equation*}

because $x_t A^{\prime }(x_t)\lt A(x_t)$ under Assumption1. This also means that $A(x_t)/x_t$ is decreasing in $x_t$ . Since $A(\!+\infty )\lt +\infty$ under Assumption1, $\lim _{x_t \rightarrow + \infty }\frac {A(x_t)}{x_t}\lt \frac {\pi _t}{g_1}$ . There is a unique solution $x_t=\varepsilon (\pi _t)$ to the equation $\frac {\pi _t}{g_1}=\frac {A(x_t)}{x_t}$ if $\lim _{x_t \rightarrow 0}\frac {A(x_t)}{x_t}\gt \frac {\pi _t}{g_1}$ .

B. Proof of Proposition1

There exist (at least) two BGPs if there are two solutions $\pi \in (\pi _1, \pi _2)$ to the equation $B_1(\pi )=B_2(\pi )$ . Since $B_1(\pi _1)\gt B_2(\pi _1)$ and $B_1(\pi _2)\gt B_2(\pi _2)$ , it requires the existence of at least one value of $\pi \in (\pi _1, \pi _2)$ such that $B_1(\pi )\lt B_2(\pi )$ .

Using equations (30) and (31), the inequality $B_1(\pi )\lt B_2(\pi )$ is equivalent to:

(B.1) \begin{equation} \Omega (\pi )\gt g-\tau _L (1-\alpha )-\tau _K \alpha \end{equation}

with

(B.2) \begin{eqnarray} \Omega (\pi ) &\equiv &\frac {a(\pi )}{R(\pi )}\left [X(\pi )-\frac {1-\delta }{a(\pi )}-(1-\tau _K)\alpha \right ] \notag \\ &&\left [ \Sigma (1-\tau _L)(1-\alpha ) +\tau _L (1-\alpha )+\tau _K \alpha -g - X(\pi )\right ] \end{eqnarray}

By construction, we have $\Omega (\pi _1)=\Omega (\pi _2)=0$ and $\Omega (\pi )\gt 0$ for all $\pi \in (\pi _1, \pi _2)$ . Let us consider $\widetilde {\pi }=\lambda \pi _1$ , with $\lambda \in (1, \pi _2/\pi _1)$ a constant independent of $g$ . Then, $\widetilde {\pi }$ does not depend on $g$ and $\Omega (\widetilde {\pi })\gt 0$ . When $g$ tends to $\tau _L (1-\alpha )+\tau _K \alpha$ , inequality (B.1) evaluated at $\pi =\widetilde {\pi }$ is satisfied. By continuity, there exists $\overline {g}\gt \tau _L (1-\alpha )+\tau _K \alpha$ such that $\Omega (\widetilde {\pi })\gt g-\tau _L (1-\alpha )-\tau _K \alpha$ for all $g\in (\tau _L (1-\alpha )+\tau _K \alpha, \overline {g})$ . This proves the existence of two solutions $\pi _I$ and $\pi _{II}$ , with $\pi _1\lt \pi _I\lt \pi _{II}\lt \pi _2$ . Since $B_1(\pi )$ and $B_2(\pi )$ are decreasing functions, the associated stationary values $b_I$ and $b_{II}$ are ranked in the following way: $b_I\gt b_{II}$ .

C. Proof of Proposition2

The dynamic system we consider is given by equation (43) and equation (25) which rewrites:

(C.1) \begin{eqnarray} \pi _{t+1} &=&\frac {X(\pi _t)\pi _t}{\Sigma (1-\tau _L)(1-\alpha ) +(\tau _L (1-\alpha )+\tau _K \alpha )-g-\frac {R(\pi _t)}{a(\pi _t)}\pi _t \Phi _t} \end{eqnarray}

Let us note $R(\pi )/a(\pi )=(1-\delta )/a(\pi ) +(1-\tau _K)\alpha \equiv \widetilde {R}(\pi )$ . We have $X^{\prime }(\pi )\gt \widetilde {R}^{\prime }(\pi )$ and, using (32), $X(\pi )\gt \widetilde {R}(\pi )$ .

Differentiating equation (43), we get:

(C.2) \begin{eqnarray} \frac {d \Phi _{t+1}}{\Phi }&=&\frac {\widetilde {R}(\pi )}{X(\pi )}\frac {d \Phi _{t}}{\Phi } +\left [\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}\frac {\widetilde {R}(\pi )}{X(\pi )}-\frac {g-\tau _L (1-\alpha )-\tau _K\alpha }{\pi \Phi X(\pi )}-\frac {X^{\prime }(\pi )\pi }{X(\pi )}\right ]\frac {d \pi _t}{\pi } \end{eqnarray}

Differentiating (C.1), we get:

(C.3) \begin{eqnarray} && \frac {d \pi _{t+1}}{\pi }=\frac {R(\pi )\pi \Phi }{\gamma }\frac {d\Phi _t}{\Phi } + \left [1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {R(\pi )\pi \Phi }{\gamma } \left (\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}+1\right )\right ]\frac {d \pi _t}{\pi } \end{eqnarray}

The trace $T$ and the determinant $D$ of the associated Jacobian matrix are given by:

(C.4) \begin{align} T & = \frac {\widetilde {R}(\pi )}{X(\pi )} +1+\frac {X^{\prime }(\pi )\pi }{X(\pi ) }+\frac {R(\pi )\pi \Phi }{\gamma } \left (\frac {\widetilde {R}^{\prime }(\pi )\pi }{ \widetilde {R}(\pi )}+1\right ) \gt 1 \end{align}
(C.5) \begin{align} D & = \frac {\widetilde {R}(\pi )}{X(\pi )}\left (1+\frac {X^{\prime }(\pi )\pi }{ X(\pi )}+\frac {R(\pi )\pi \Phi }{\gamma }\right )+\frac {R(\pi )\pi \Phi }{\gamma } \left ( \frac {g-\tau _L (1-\alpha )-\tau _K\alpha }{\pi \Phi X(\pi )}+\frac {X^{\prime }(\pi )\pi }{X(\pi )}\right ) \\[9pt]\nonumber\end{align}

and define the characteristic polynomial $P(\lambda )\equiv \lambda ^2-T\lambda +D =0$ . We have $P(0)=D\gt 0$ . Since $P(\!-\infty )=+\infty$ , $P(\!+\infty )=+\infty$ , and $T\gt 1$ , the two roots are positive or complex conjugates.

Using $\widetilde {R}(\pi )=R(\pi )/a(\pi )$ , $\gamma =a(\pi )X(\pi )$ and equation (43) at a steady state, $T$ and $D$ rewrite:

(C.6) \begin{align} T&=\frac {\widetilde {R}(\pi )}{X(\pi )} + 1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {\widetilde {R}(\pi )\pi \Phi }{X(\pi )} \left (\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}+1\right ) \gt 1 \end{align}
(C.7) \begin{align} D&=\frac {\widetilde {R}(\pi )}{X(\pi )}\left (1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {\widetilde {R}(\pi )\pi \Phi }{X(\pi )}\right )+\frac {\widetilde {R}(\pi )\pi \Phi }{X(\pi )}\left ( 1-\frac {\widetilde {R}(\pi )}{ X(\pi )}+\frac {X^{\prime }(\pi )\pi }{X(\pi )}\right ) \notag \\ & =\frac {\widetilde {R}(\pi )}{X(\pi )}\left (1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}\right )(1+\pi \Phi ) \\[9pt]\nonumber\end{align}

We deduce that:

(C.8) \begin{eqnarray} P(1)&=&1-T+D \notag \\ &=& -\left (1-\frac {\widetilde {R}(\pi )}{X(\pi )}\right )\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {\widetilde {R}(\pi )\pi \Phi }{X(\pi )} \left (\frac {X^{\prime }(\pi )\pi }{X(\pi )}-\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}\right ) \end{eqnarray}

Using (44) and (45), we rewrite $B_3(\pi )$ and $B_4(\pi )$ as follows:

(C.9) \begin{align} B_3(\pi )&=\frac {g -\tau _L (1-\alpha )-\tau _K \alpha }{\pi [X(\pi )-\widetilde {R}(\pi ) ]} \end{align}
(C.10) \begin{align} B_4(\pi ) &=\frac {1}{\widetilde {R}(\pi )\pi } \left [ \Sigma (1-\tau _L)(1-\alpha ) +\tau _L (1-\alpha )+\tau _K \alpha -g - X(\pi )\right ] \\[9pt]\nonumber\end{align}

We deduce that:

(C.11) \begin{align} \frac {B_3^{\prime }(\pi )\pi }{B_3(\pi )}&=-1-\frac {X^{\prime }(\pi )-\widetilde {R}^{\prime }(\pi )}{X(\pi )-\widetilde {R}(\pi )} \pi \end{align}
(C.12) \begin{align} \frac {B_4^{\prime }(\pi )\pi }{B_4(\pi )} &=-1-\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}-\frac {X(\pi )}{\Phi \widetilde {R}(\pi )\pi }\frac {X^{\prime }(\pi )\pi }{X(\pi )} \\[9pt]\nonumber\end{align}

After some computations, we can show that $\frac {B_3^{\prime }(\pi )\pi }{B_3(\pi )}\lt \frac {B_4^{\prime }(\pi )\pi }{B_4(\pi )}$ is equivalent to:

(C.13) \begin{eqnarray} &&\left (1-\frac {\widetilde {R}(\pi )}{X(\pi )}\right )\left (-\frac {X^{\prime }(\pi )\pi }{X(\pi )}\right )+\frac {\widetilde {R}(\pi )\pi \Phi }{X(\pi )} \left (\frac {X^{\prime }(\pi )\pi }{X(\pi )}-\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}\right )\gt 0 \end{eqnarray}

We recall that $B_3^{\prime }(\pi _I)\lt B_4^{\prime }(\pi _I)$ and $B_3^{\prime }(\pi _{II})\gt B_4^{\prime }(\pi _{II})$ . By inspection of equations (C.8) and (C.13), we deduce that $P(1)\gt 0$ at the steady state $(\pi _I, \Phi _I)$ and $P(1)\lt 0$ at the steady state $(\pi _{II}, \Phi _{II})$ . Therefore, the steady-state $(\pi _{II}, \Phi _{II})$ is a saddle, with one eigenvalue between 0 and 1 and one higher than 1. The steady-state $(\pi _I, \Phi _I)$ is stable if $D\lt 1$ and unstable if $D\gt 1$ .

Using (C.7) and $\pi \Phi =b$ , $D\lt 1$ is equivalent to:

(C.14) \begin{equation} \widetilde {R}(\pi )(1+b)\lt \frac {X(\pi )}{1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}} \end{equation}

Using the expression of $\widetilde {R}(\pi )$ and equation (31), the left-hand side of inequality (C.14) is given by:

(C.15) \begin{align} \widetilde {R}(\pi )(1+b)&=\Sigma (\pi )(1-\tau _L)(1-\alpha )-[g-\tau _L (1-\alpha )-\tau _K \alpha ]-[X(\pi )-\widetilde {R}(\pi )]\nonumber\\&\lt \Sigma (\pi )(1-\tau _L)(1-\alpha ) \end{align}

which is arbitrarily low if the tax rate $\tau _L$ is high. Using $X(\pi )= \frac {1-m}{a(\pi )}- \frac {\psi g_2}{\pi } +\frac {\mu }{a(\pi )\pi }$ , we also have:

(C.16) \begin{eqnarray} \frac {X^{\prime }(\pi )\pi }{X(\pi )}=\frac {-\frac {1-m}{a(\pi )}\epsilon _a (\pi )+ \frac {\psi g_2}{\pi } -\frac {\mu }{a(\pi )\pi }(1+\epsilon _a (\pi ))}{\frac {1-m}{a(\pi )}- \frac {\psi g_2}{\pi } +\frac {\mu }{a(\pi )\pi }} \end{eqnarray}

Therefore, $X^{\prime }(\pi )\pi /X(\pi )$ is not too high if $\epsilon _a (\pi )$ is not too negative.

We deduce that we have $D\lt 1$ , which ensures the stability of the steady state $(\pi _I, \Phi _I)$ , if $\tau _L$ is sufficiently high and close to 1 and $\epsilon _a (\pi )$ is not too negative. Note that at least when $D$ is close to 1, the eigenvalues are complex conjugates, meaning that the dynamic path converges with oscillations around the steady state.

D. Proof of Proposition3

The Jacobian matrix of the linearized dynamic system is given by:

\begin{equation*} J= \left (\begin{array}{c@{\quad}c} \frac {\widetilde {R}(\pi )}{X(\pi )} & \frac {\widetilde {R}^{\prime }(\pi )\pi }{\vphantom{\frac{1}{2}}\widetilde {R}(\pi )}\frac {\widetilde {R}(\pi )}{X(\pi )}-\frac {g-\tau _L (1-\alpha )-\tau _K\alpha }{\pi \Phi X(\pi )}-\frac {X^{\prime }(\pi )\pi }{X(\pi )} \\ \frac {R(\pi )\pi \Phi }{\gamma } & 1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {R(\pi )\pi \Phi }{\gamma } \left (\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}+1\right )\end{array}\right ) \end{equation*}

Let $E_s=(1,e_s)$ be the eigenvector associated to the stable eigenvalue $\lambda _s\in (0,1)$ . We have $J E_s=\lambda _s E_s$ . Using the second equation of this system, we deduce that:

\begin{equation*} \frac {R(\pi )\pi \Phi }{\gamma }=e_s\left [\lambda _s-\left [1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {R(\pi )\pi \Phi }{\gamma } \left (\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}+1\right )\right ]\right ] \end{equation*}

Since we have:

\begin{equation*} \lambda _s \lt 1\lt 1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {R(\pi )\pi \Phi }{\gamma } \left (\frac {\widetilde {R}^{\prime }(\pi )\pi }{ \widetilde {R}(\pi )}+1\right ) \end{equation*}

we deduce that $e_s\lt 0$ . This means that the eigenvector associated with the stable eigenvalue has a negative slope. On the stable manifold, we have $\frac {d\Phi _t}{\Phi }/\frac {d \pi _t}{\pi }=1/e_s\lt 0$ in the neighborhood of the steady state.

Let $E_u=(1,e_u)$ be the eigenvector associated to the unstable eigenvalue $\lambda _u \gt 1$ . $det (J-\lambda _uI)=0$ is equivalent to:

\begin{eqnarray*} &&\left ( \frac {\widetilde {R}(\pi )}{X(\pi )}-\lambda _u \right )\left [1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {R(\pi )\pi \Phi }{\gamma } \left (\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}+1\right )- \lambda _u \right ] \\ &&=-\frac {R(\pi )\pi \Phi }{\gamma } \left [\frac {g-\tau _L (1-\alpha )-\tau _K\alpha }{\pi \Phi X(\pi )}+\frac {X^{\prime }(\pi )\pi }{X(\pi )}- \frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}\frac {\widetilde {R}(\pi )}{X(\pi )}\right ] \end{eqnarray*}

Since $X^{\prime }(\pi )\gt \widetilde {R}^{\prime }(\pi )$ , $X(\pi )\gt \widetilde {R}(\pi )$ and $\lambda _u \gt 1$ , the right-hand side of this equation is strictly negative, which implies that:

\begin{equation*} \lambda _u \lt 1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {R(\pi )\pi \Phi }{\gamma } \left (\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}+1\right ) \end{equation*}

Using

\begin{equation*} \frac {R(\pi )\pi \Phi }{\gamma }=e_u\left [\lambda _u-\left [1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}+\frac {R(\pi )\pi \Phi }{\gamma } \left (\frac {\widetilde {R}^{\prime }(\pi )\pi }{\widetilde {R}(\pi )}+1\right )\right ]\right ] \end{equation*}

we deduce that $e_u\lt 0$ , which means that the eigenvector associated to the unstable eigenvalue has a negative slope. On the unstable manifold, we have $\frac {d\Phi _t}{\Phi }/\frac {d \pi _t}{\pi }=1/e_u\lt 0$ in the neighborhood of the steady state. Since $\lambda _u\gt \lambda _s$ , we even have $1/e_u\gt 1/e_s$ , which means that around the steady state, the negative slope of the unstable manifold is greater than the one of the stable manifold.

E. Proof of Proposition4

Using the proof of Proposition2, the steady state $(\pi _{II}, \Phi _{II})$ is saddle. We know that the steady state $(\pi _I, \Phi _I)$ is unstable if $D\gt 1$ , i.e.

(E.1) \begin{equation} \widetilde {R}(\pi )(1+b)\gt \frac {X(\pi )}{1+\frac {X^{\prime }(\pi )\pi }{X(\pi )}} \end{equation}

which requires $\tau _L$ sufficiently low and $\epsilon _a (\pi )$ sufficiently negative.

Using $X(\pi )= \frac {1-m}{a(\pi )}- \frac {\psi g_2}{\pi } +\frac {\mu }{a(\pi )\pi }$ and (C.16), the inequality $X^{\prime }(\pi )\pi /X(\pi )\gt 1$ can be written:

\begin{eqnarray*} -\frac {1-m}{a(\pi )}\epsilon _a (\pi )+ \frac {\psi g_2}{\pi } -\frac {\mu }{a(\pi )\pi }(1+\epsilon _a (\pi ))\gt \frac {1-m}{a(\pi )}- \frac {\psi g_2}{\pi } +\frac {\mu }{a(\pi )\pi } \end{eqnarray*}

which is equivalent to:

\begin{eqnarray*} -\frac {1-m}{a(\pi )}(1+\epsilon _a (\pi ))+ 2\frac {\psi g_2}{\pi } -\frac {\mu }{a(\pi )\pi }(2+\epsilon _a (\pi ))\gt 0 \end{eqnarray*}

This inequality is satisfied for $\epsilon _a (\pi )\lt -2$ . Using the fact that $X^{\prime }(\pi )\pi /X(\pi )\gt 1$ , inequality (E.1) is satisfied if:

(E.2) \begin{equation} \widetilde {R}(\pi )b\gt \frac {X(\pi )}{2}-\widetilde {R}(\pi ) \end{equation}

The left-hand side of this inequality is positive, whereas the right-hand side is strictly negative if $X(\pi )\lt 2\widetilde {R}(\pi )$ , i.e.

(E.3) \begin{equation} -\frac {\psi g_2}{\pi }\lt \frac {1-2\delta +m}{a(\pi )}+2(1-\tau _K)\alpha -\frac {\mu }{a(\pi )\pi } \end{equation}

Since $a(\pi )\pi$ is decreasing, this last inequality is satisfied if $1+m\gt 2\delta$ and $1-\tau _K\gt \mu /[2\alpha a(\pi _2)\pi _2]$ . It gives a sufficient condition to have $D\gt 1$ .

F. Proof of Lemma3

Using (26)–(28), the system that defines the stationary solutions can be rewritten as:

(F.1) \begin{align} b_i &=\frac {g_1+g_2 -(\tau _L (1-\alpha )+\tau _K \alpha ) }{z(\pi _i, g_1, g_2)-y(\pi _i, g_1, \tau _K)}\equiv \mathcal {B}_1(\pi _i, g_1, g_2, \tau _K, \tau _L) \end{align}
(F.2) \begin{align} b_i&=\frac {\left [\Sigma (1-\tau _L)(1-\alpha ) +\tau _L (1-\alpha )+\tau _K \alpha -(g_1+g_2) -z(\pi _i, g_1, g_2)\right ]}{y(\pi _i, g_1, \tau _K)} \notag \\ &\equiv \mathcal {B}_2(\pi _i, g_1, g_2, \tau _K, \tau _L)\\[8pt]\nonumber \end{align}

with

\begin{equation*} y(\pi _i, g_1, \tau _K)=(1-\delta )/\Theta (\pi _i, g_1)+\alpha (1-\tau _K). \end{equation*}

The sign of partial derivatives are: $y_{\pi _i}\gt 0$ ; $y_{g_1}\lt 0$ ; $y_{\tau _K}\lt 0$ ; and

\begin{equation*} z(\pi _i, g_1, g_2)=\frac {1-m}{\Theta (\pi _i, g_1)}- \frac {\psi g_2}{\pi _i} +\frac {\mu }{\Theta (\pi _i, g_1)\pi _i}. \end{equation*}

The sign of partial derivatives are: $z_{\pi _i}\gt 0$ ; $z_{g_1}\lt 0$ ; $z_{g_2}\lt 0$ .

Total differentiation of Equations (F.1) and (F.2) gives:

\begin{equation*} \mathcal {C} \times \begin{bmatrix} \mathrm {d}b \\ \mathrm {d} \pi \end{bmatrix}= \mathcal {D} \times \begin{bmatrix} \mathrm {d}g_1 \\ \mathrm {d} g_2 \\ \mathrm {d}\tau _L \\ \mathrm {d}\tau _K\end{bmatrix}\end{equation*}

where

\begin{equation*} \mathcal {C} =\begin{bmatrix} 1 & -B^{\prime }_1(\pi ) \\ 1 & -B^{\prime }_2(\pi )\end{bmatrix}\end{equation*}

with $B^{\prime }_1(\pi )=\frac {\partial \mathcal {B}_1(\pi, g_1, g_2, \tau _K, \tau _L)}{\partial \pi }$ and $B^{\prime }_2(\pi )=\frac {\partial \mathcal {B}_2(\pi, g_1, g_2, \tau _K, \tau _L)}{\partial \pi }$ and where

\begin{equation*} \mathcal {D}=\begin{bmatrix} \mathcal {D}_1 & \mathcal {D}_2 & \mathcal {D}_3 & \mathcal {D}_4 \\ \mathcal {D}_5 & \mathcal {D}_6 & \mathcal {D}_7 & \mathcal {D}_8\end{bmatrix}\end{equation*}
(F.3) \begin{equation} = \begin{bmatrix} \frac {1+\frac {b\Theta _{g_1}}{(\Theta )^2}\left (\delta -m+\mu /\pi \right )}{z-y} & \frac {\pi +\psi b}{\pi (z-y)} & \frac {\alpha -1}{z-y} & \frac {-\alpha (1+b)}{z-y} \\ \frac {-1+\frac {\Theta _{g_1}}{(\Theta )^2}\left (b(1-\delta )+\mu /\pi +1-m\right )}{y} & \frac {\psi -\pi }{\pi y} & \frac {(1-\alpha )(1-\Sigma )}{y} & \alpha \left (\frac {1+b}{y}\right )\end{bmatrix} \end{equation}

The determinant of matrix $\mathcal {C}$ is $\det \mathcal {C}(\pi )= B^{\prime }_1(\pi )-B^{\prime }_2(\pi )$ . We recall that $B_1^{\prime }(\pi _I)\lt B_2^{\prime }(\pi _I)$ and $B_1^{\prime }(\pi _{II})\gt B_2^{\prime }(\pi _{II})$ . We thus have $\det \mathcal {C}(\pi _I)\lt 0$ and $\det \mathcal {C}(\pi _{II})\gt 0$ . Moreover, under Assumptions13, we have $\mathcal {D}_1, \mathcal {D}_2, \mathcal {D}_7, \mathcal {D}_8 \gt 0$ , $\mathcal {D}_3, \mathcal {D}_4\lt 0$ . The sign of $\mathcal {D}_5$ and $\mathcal {D}_6$ , referring to the impact of environmental policy variables, $g_1$ and $g_2$ , depends on policy and model parameters and is discussed later.

We obtain the effects of policy variables on stationary variables examining:

\begin{equation*} \begin{bmatrix} \mathrm {d}b \\ \mathrm {d}\pi \end{bmatrix}=\frac {1}{\det \mathcal {C}(\pi )}\times \begin{bmatrix} -B^{\prime }_2(\pi ) & B^{\prime }_1(\pi ) \\ -1 & 1\end{bmatrix} \times \mathcal {D} \times \begin{bmatrix} \mathrm {d}g_1 \\ \mathrm {d}g_2 \\ \mathrm {d}\tau _L \\ \mathrm {d}\tau _K\end{bmatrix}\end{equation*}
  • Effect of $g_1$

    \begin{equation*} \frac {\mathrm {d}b}{\mathrm {d}g_1}=\frac {-B^{\prime }_2(\pi )\mathcal {D}_1+B^{\prime }_1(\pi )\mathcal {D}_5}{\det \mathcal {C}(\pi )} \end{equation*}
    \begin{equation*} \frac {\mathrm {d}\pi }{\mathrm {d}g_1}=\frac {-\mathcal {D}_1+\mathcal {D}_5}{\det \mathcal {C}(\pi ) } \end{equation*}
    Using the expression for $\mathcal {D}_1$ and $\mathcal {D}_5$ given in (F.3), we have:
    \begin{eqnarray*} \frac {\mathrm {d}\pi }{\mathrm {d}g_1}&=&\left (\frac {-1+\frac {\Theta _{g_1}}{\Theta ^2}\left (b(1-\delta )+\frac {\mu }{\pi }+1-m\right )}{y}\right . \\ &&\left .-\frac {1+\frac {b\Theta _{g_1}}{\Theta ^2}\left (\delta -m+\frac {\mu }{\pi }\right )}{z-y} \right )\times \frac {1}{\det \mathcal {C}(\pi )} \end{eqnarray*}
    Rewriting the previous equation, we obtain:
    (F.4) \begin{eqnarray} \frac {\mathrm {d}\pi }{\mathrm {d}g_1}&=& \left (\frac {-z}{y(z-y)} +\frac {\Theta _{g_1}}{y\Theta ^2}\left [ b(1-\delta )+\frac {\mu }{\pi }+1-m \right .\right . \notag \\ && \left . \left . -b\frac {y}{z-y}\left (\frac {\mu }{\pi }+\delta -m\right )\right ]\right )\times \frac {1}{\det \mathcal {C}(\pi )} \end{eqnarray}
    Using (F.1) and (F.2), we have:
    \begin{equation*} \frac {y}{z-y}\lt \frac {x(1-\tau _L)(1-\alpha )}{g-\tau _L (1-\alpha )-\tau _K\alpha } \end{equation*}
    The right-hand side of this inequality is small if $\tau _L$ is high and $g-\tau _L (1-\alpha )-\tau _K\alpha$ is not too small. In this case, the term into brackets in equation (F.4) is positive such that for $\Theta _{g_1}$ high enough, we have $\textrm {sgn}\left (\frac {\mathrm {d}\pi _I}{\mathrm {d}g_1}\right )\lt 0$ and $\textrm {sgn}\left (\frac {\mathrm {d}\pi _{II}}{\mathrm {d}g_1}\right )\gt 0$ . Given the expression for $\frac {\mathrm {d} b}{\mathrm {d}g_1}$ , we do not conclude concerning the impact of $g_1$ on $b$ in this configuration. We can note that as long as $\Theta _{g_1}$ is low enough, $\mathcal {D}_5\lt 0$ . This implies $\frac {\mathrm {d}b}{\mathrm {d}g_1}\lt 0$ and $\frac {\mathrm {d}\pi }{\mathrm {d}g_1}\gt 0$ for the BGP $i=I$ and the reverse for the BGP $i=II$ .
  • Effect of $g_2$

    \begin{equation*} \frac {\mathrm {d}b}{\mathrm {d}g_2}=\frac {-B^{\prime }_2(\pi )\mathcal {D}_2+B^{\prime }_1(\pi )\mathcal {D}_6}{\det \mathcal {C}(\pi )} \end{equation*}
    \begin{equation*} \frac {\mathrm {d}\pi }{\mathrm {d}g_2}=\frac {-\mathcal {D}_2+\mathcal {D}_6}{\det \mathcal {C}(\pi )}=\frac {1}{\pi }\left (\frac {\psi -\pi }{y}-\frac {\pi +\psi b}{z-y}\right ) \end{equation*}
    Note that the two steady states exist with strictly positive values $\pi _I$ and $\pi _{II}$ even if $\psi$ tends to 0. This means that there exists $\bar {\psi }\gt 0$ such that $\psi \lt \pi _I\lt \pi _{II}$ for all $\psi \lt \bar {\psi }$ . In this case, we have $\mathcal {D}_6\lt 0$ . Then, an increase in $g_2$ means that $\pi _I$ and $b_{II}$ increases while $\pi _{II}$ and $b_I$ decrease.
  • Effect of $\tau _L$

    \begin{equation*} \frac {\mathrm {d}b}{\mathrm {d}\tau _L}=\frac {-B^{\prime }_2(\pi )\mathcal {D}_3+B^{\prime }_1(\pi )\mathcal {D}_7}{\det \mathcal {C}(\pi )} \end{equation*}
    \begin{equation*} \frac {\mathrm {d}\pi }{\mathrm {d}\tau _L}=\frac {-\mathcal {D}_3+\mathcal {D}_7}{\det \mathcal {C}(\pi )} \end{equation*}
    We have $\frac {\mathrm {d}\pi }{\mathrm {d}\tau _L}\lt 0$ and $\frac {\mathrm {d}b}{\mathrm {d}\tau _L}\gt 0$ for the BGP $i=I$ and the reverse, $\frac {\mathrm {d}\pi }{\mathrm {d}\tau _L}\gt 0$ and $\frac {\mathrm {d}b}{\mathrm {d} \tau _L}\lt 0$ , for the BGP $i=II$ .
  • Effect of $\tau _K$

    \begin{equation*} \frac {\mathrm {d}b}{\mathrm {d}\tau _K}=\frac {-B^{\prime }_2(\pi )\mathcal {D}_4+B^{\prime }_1(\pi )\mathcal {D}_8}{\det \mathcal {C}(\pi )} \end{equation*}
    \begin{equation*} \frac {\mathrm {d}\pi }{\mathrm {d}\tau _K}=\frac {-\mathcal {D}_4+\mathcal {D}_8}{\det \mathcal {C}(\pi )} \end{equation*}
    We have $\frac {\mathrm {d}\pi }{\mathrm {d}\tau _K}\lt 0$ and $\frac {\mathrm {d}b}{\mathrm {d}\tau _K}\gt 0$ for the BGP $i=I$ and the reverse, $\frac {\mathrm {d}\pi }{\mathrm {d}\tau _K}\gt 0$ and $\frac {\mathrm {d}b}{\mathrm {d} \tau _K}\lt 0$ , for the BGP $i=II$ .

G. Proof of Proposition5

We examine how policy variables affect $\gamma _I$ to deduce the welfare effect along the stable BGP $i=I$ . Using (38), we have:

(G.1) \begin{equation} \gamma _I = 1-m+\frac {\mu -\psi g_2 a(\pi _I)}{\pi _I} \end{equation}

The proof of Lemma3 reveals that $\frac {\mathrm {d}\pi _I}{\mathrm {d}\tau _L}\lt 0$ and $\frac {\mathrm {d}\pi _I}{\mathrm {d}\tau _K}\lt 0$ and gives the conditions to have $\frac {\mathrm {d}\pi _I}{\mathrm {d}g_1}\lt 0$ or $\gt 0$ . In the configuration where the fall in $\pi _I$ is good for growth $\left (-\epsilon _a(\pi )\lt \frac {\mu -\psi g_2 a(\pi )}{\psi g_2 a(\pi )}\right )$ , we directly have the policy scenarios for fiscal and adaptation policies that increase welfare.

As regards mitigation $g_2$ , it exerts a direct negative effect on growth, in addition to its effect through $\pi$ . Lemma3 gives a condition such that $\frac {\mathrm {d}\pi _I}{\mathrm {d}g_2}\gt 0$ , meaning that in the configuration where the fall in $\pi _I$ is good for growth, the increase in $g_2$ has a double negative effect on welfare.

Footnotes

1 In December 2023, among 67 low-income countries, the Debt Sustainability Analysis provided by the World Bank Group and the IMF identified 28 countries with a high risk of overall debt distress and 11 countries already in distress. This worrisome observation is coupled with the increase in the frequency and intensity of extreme weather events in these vulnerable countries.

3 Note that Bacchiocchi et al. (Reference Bacchiocchi, Bellocchi, Bischi and Travaglini2024) also exploit the saddle property of a steady state to exhibit a form of unsustainability of public debt. However, they are interested on the interplay between debt and inflation, whereas we are concerned with pollution, debt, and growth.

4 There is also a literature analyzing the link between environmental quality and production in the absence of public debt. For a recent reference, see Menuet et al. (Reference Menuet, Minea, Villieu and Xepapadeas2024) where public pollution abatement is introduced into the production function and positively affects the efficiency of pollution as an input.

5 Burke et al. (Reference Burke, Hsiang and Miguel2015) finds a nonlinear decline in macroeconomic productivity following a change in temperature, across sectors, in 166 rich and poor countries since the 1960s. This finding is confirmed by Kalkuhl and Wenz (Reference Kalkuhl and Wenz2020) which looks at the effects of rising global mean surface temperature on production.

6 Adaptation is formalized as a flow, that reduces the damage of pollution stock. This implies that adaptation expenses are recurrent during each period. In this time perspective, we suppose that changing climatic conditions necessitate ongoing adjustments, updates, or improvements.

7 As it is usually the case in models with pollution, production and public debt, we consider that there is no market imperfection on the asset market. Then, our results will not depend on a specific financial market imperfection.

8 In the working paper version of the paper Davin et al. (Reference Davin, Fodha and Seegmuller2024), the analysis is extended to the case of a CES utility function where consumption in both periods are substitutable. In such a case, the saving function is increasing in the real interest factor $R_{t+1}$ .

9 Note that it is always satisfied if $A(0)\gt 0$ . If $A(0)=0$ , the condition is equivalent to $\lim _{x \rightarrow 0}A^{\prime }(x)\gt \pi _t/g_1$ .

10 We can notice that the negative relationship between $b$ and $\pi$ may be observed because the depreciation of capital is not complete ( $\delta \ne 1$ ). Otherwise, when $\delta =1$ , $R(\pi )/a(\pi )$ is constant. From equation (26), the debt ratio $b$ does no longer evolve with $\pi$ . The reduction in the primary deficit and the cost of debt induced by the increase in $\pi$ are proportional to the decrease in savings so that the debt-to-capital ratio does not depend on productivity and thus not on $\pi$ .

11 Oscillatory convergence is one possible scenario we may observe, and discussion concerning such a convergence is not the central focus of our study.

12 Note that the assumption $1+m\gt 2\delta$ fit with $\delta \geq m$ . Indeed, since $m\lt 1$ , the interval $[m, (1+m)/2)$ is non-empty and we can choose $\delta$ in this interval.

13 Based on Proposition2, this implies a TFP vulnerability to pollution not too important. We also consider in this section the case in which growth is decreasing with pollution ratio $\pi$ , i.e. $-\epsilon _a(\pi )\lt \frac {\mu -\psi g_2 a(\pi )}{\psi g_2 a(\pi )}$ .

14 The analysis of the optimal allocation and its decentralization using fiscal instruments is out of the scope of this paper. It could create several problems linked to the existence of a solution because of the Ak production function already characterized by increasing returns and where the TFP is a function of both the adaptation policy and the stock of pollution.

15 The effects of comparative statics on the other steady state $(b_{II}, \pi _{II})$ are the opposite.

References

Andersen, D. C. (2017) Do credit constraints favor dirty production? Theory and plant-level evidence. Journal of Environmental Economics and Management 84, 189208.CrossRefGoogle Scholar
Andersen, T. M., Bhattacharya, J. and Liu, P.. (2020) Resolving intergenerational conflict over the environment under the pareto criterion. Journal of Environmental Economics and Management 100, 102290.CrossRefGoogle Scholar
Bacchiocchi, A., Bellocchi, A., Bischi, G. I. and Travaglini, G.. (2024) A nonlinear model of public debt with bonds and money finance. Economia Politica 41(2), 457498.CrossRefGoogle Scholar
Baret, M. and Menuet, M.. (2024) Fiscal and environmental sustainability: Is public debt environmentally friendly? Environmental and Resource Economics 87(6), 124.CrossRefGoogle Scholar
Boly, M., Combes, J.-L., Menuet, M., Minea, A., Motel, P. C. and Villieu, P.. (2022) Can public debt mitigate environmental debt? Theory and empirical evidence. Energy Economics 111, 105895.CrossRefGoogle Scholar
Buhr, B., Donovan, C., Kling, G., Lo, Y., Murinde, V., Pullin, N. and Volz, U.. (2018). Climate change and the cost of capital in developing countries. Imperial College London; SOAS University of London; UN Environment. https://unepinquiry.org/publication/climate-change-and-the-cost-of-capital-in-developing-countries/.Google Scholar
Burke, M., Hsiang, S. M. and Miguel, E.. (2015) Global non-linear effect of temperature on economic production. Nature 527(7577), 235239.CrossRefGoogle ScholarPubMed
Catalano, M., Forni, L. and Pezzolla, E.. (2020) Climate-change adaptation: the role of fiscal policy. Resource and Energy Economics 59, 101111.CrossRefGoogle Scholar
Chiroleu-Assouline, M. and Fodha, M.. (2014) From regressive pollution taxes to progressive environmental tax reforms. European Economic Review 69, 126142.CrossRefGoogle Scholar
Cole, M. A., Elliott, R. J. and Shimamoto, K.. (2005) Industrial characteristics, environmental regulations and air pollution: An analysis of the UK manufacturing sector. Journal of Environmental Economics and Management 50(1), 121143.CrossRefGoogle Scholar
Dasgupta, S., van Maanen, N., Gosling, S. N., Piontek, F., Otto, C. and Schleussner, C.-F.. (2021) Effects of climate change on combined labour productivity and supply: an empirical, multi-model study. The Lancet Planetary Health 5(7), e455e465.CrossRefGoogle ScholarPubMed
Davin, M., Fodha, M. and Seegmuller, T.. (2024). Pollution, public debt, and growth: the question of sustainability. AMSE Working Papers No. 2418, Aix-Marseille School of Economics, France.Google Scholar
Diamond, P. (1965) National debt in a neoclassical growth model. American economic review 55(5), 11261150.Google Scholar
Dibley, A., Wetzer, T. and Hepburn, C.. (2021) National COVID debts: climate change imperils countries’ ability to repay. Nature 592(7853), 184187.CrossRefGoogle ScholarPubMed
Dietz, S. and Stern, N.. (2015) Endogenous growth, convexity of damage and climate risk: how nordhaus’ framework supports deep cuts in carbon emissions. The Economic Journal 125(583), 574620.CrossRefGoogle Scholar
Fodha, M. and Seegmuller, T.. (2012) A note on environmental policy and public debt stabilization. Macroeconomic Dynamics 16(3), 477492.CrossRefGoogle Scholar
Fodha, M. and Seegmuller, T.. (2014) Environmental quality, public debt and economic development. Environmental and Resource Economics 57(4), 487504.CrossRefGoogle Scholar
Fodha, M., Seegmuller, T. and Yamagami, H.. (2018) Environmental tax reform under debt constraint. Annals of Economics and Statistics 129, 3352.CrossRefGoogle Scholar
Golosov, M., Hassler, J., Krusell, P. and Tsyvinski, A.. (2014) Optimal taxes on fossil fuel in general equilibrium. Econometrica 82(1), 4188.Google Scholar
Heijdra, B. J., Kooiman, J. P. and Ligthart, J. E.. (2006) Environmental quality, the macroeconomy, and intergenerational distribution. Resource and Energy Economics 28(1), 74104.CrossRefGoogle Scholar
IMF. (2023) International monetary fund. fiscal monitor: Climate crossroads: Fiscal policies in a warming world. Washington, DC: IMF.Google Scholar
Kalkuhl, M. and Wenz, L.. (2020) The impact of climate conditions on economic production. evidence from a global panel of regions. Journal of Environmental Economics and Management 103, 102360.CrossRefGoogle Scholar
Menuet, M., Minea, A., Villieu, P. and Xepapadeas, A.. (2024) Environmental quality along the process of economic growth: a theoretical reappraisal. Economic Theory 77(4), 12191258.CrossRefGoogle Scholar
Pisani-Ferry, J. and Mahfouz, S.. (2023). The economic implications of climate action. France Stratégie. https://www.strategie.gouv.fr/en/publications/economic-implications-climate-action.Google Scholar
Sadler, A., Ranger, N., Fankhauser, S., Marotta, F. and O’Callaghan, B.. (2024) The impact of COVID-19 fiscal spending on climate change adaptation and resilience. Nature Sustainability 7(3), 270281.CrossRefGoogle Scholar
Seghini, C. and Dees, S.. (2024) The green transition and public finances. Banque de France Working Paper No. 949.CrossRefGoogle Scholar
Birkmann, J., E., Liwenga, R., Pandey, E., Boyd, R., Djalante, F., Gemenne, W., Filho, Leal, Pinho, P.F., Stringer, L. and Wrathall, D.. (2022) Poverty, Livelihoods and Sustainable Development. In: Pörtner, H.-O., Roberts, D. C., Tignor, M., Poloczanska, E. S., Mintenbeck, K., Alegría, A., Craig, M., Langsdorf, S., Löschke, S., Möller, V., Okem, A., Rama, B. (eds.), Climate Change 2022: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK and New York, NY, USA: Cambridge University Press, pp. 11711274.Google Scholar
Woodford, M. (1990) Public debt as private liquidity. The American Economic Review 80(2), 382388.Google Scholar
WorldBank. (2023). International Debt Report 2023. Washington, DC: World Bank. http://hdl.handle.net/10986/40670.Google Scholar
Zenios, S. A. (2024) The climate-sovereign debt doom loop: what does the literature suggest? Current Opinion in Environmental Sustainability 67, 101414.CrossRefGoogle Scholar
Figure 0

Figure 1. Evolution of the central government debt to GDP ratio by groups of countries (World bank classification in 2023). IMF global debt database.

Figure 1

Figure 2. Dynamics with sustainability.

Figure 2

Figure 3. Endogenous tipping zone (ETZ).

Figure 3

Figure 4. Dynamics with unsustainability.