# Nonlinear Dynamics in Coevolution of Economic and Environmental Systems

^{1}Department of Economics and Law, University of Macerata, Macerata, Italy^{2}Department of Economics and Management, University of Pisa, Pisa, Italy

The goal of this paper is to review some of the works on the dynamics of economic models in which the environmental variable is introduced. In particular, we will focus on models that in addition to giving some analytical insights into the coevolution of economic and environmental systems, they can give rise to nonlinear dynamic effects as the emergence of chaotic dynamics, multistability and (local and global) indeterminacy. However, several questions on this issue remain open and we hope that our paper may attract researchers to contribute to this topic.

## 1. Introduction

The literature on economic growth has so far focused on highly stylized models^{1}. More recently, an increasing attention has been paid to the analysis of models which consider more realistic interplays among variables and the possible emergence of non-trivial dynamics.

In this paper, several models aimed to explore particular connections between economic activity and the environment will be examined. In particular, we will focus on the analysis of models described by dynamic systems able to highlights both the long-term trends of the variables and the emergence of short or medium-term phenomena arising from their interaction.

To this end, important models of the economic the ory, such as the Solow model, the Ramsey-Cass-Koopmans model, and the OLG structure are examined when environmental variables are introduced. In particular, we will analyse how the different approaches may affect the dynamic properties of the models. Indeed, the choice of considering the environment as productive input or as consumption good, the different hypotheses on the agents' rationality (and then the different allocative problems), and the choice between continuous and discrete time framework have relevant consequences on the models' dynamics. For example, depending on alternative assumptions, the productivity in the private sector could be effective or not in defining the complexity of the system. Moreover, another decisive axiom is related to the complementarity or substitutability between environmental and private good. Nevertheless, cyclical dynamics or multistability seem to be likely results whenever models take into account the interaction between economic and environmental systems.

The remainder of our paper is organized as follows: In section 2, we will analyse two seminal models on this issue. In section 3, we will review several models in a discrete time framework. In section 4, we will present the continuous time seminal model of Antoci et al. [2]. In section 5, we will discuss models in which a fixed time delay is introduced. Finally, section 6 will conclude.

## 2. Two General Frameworks

In this section, we review two milestone models in the literature on interactions between economic and environmental systems.

### 2.1. The Green Solow Model

One of the first models which highlights the possible existence of nonmonotonic relationships^{2} between economic activity and the environment has been proposed by Brock and Taylor [4].^{3} Specifically, the authors consider an extended version of the Solow model in which savings^{4} and environmental investment choices are fixed.

In the model, the output *Y* is generated by a strictly concave production function *F* with constant returns to scale and where the inputs are the effective labor *BL*, where *B* is the labor-augmenting technological progress, and the capital *K*. Capital accumulates via savings, *sY*, with *s* ∈ (0, 1) and depreciates at δ > 0. The rate of labor-augmenting technological progress, *g*_{B} > 0, determines the growth rate of labor efficiency:

The impact of pollution has been modeled assuming that each unit of economic activity, *F*, generates Ω units of pollution. Then, the amount of pollution differs from the amount generated by the economic activity only when abatement occurs. Moreover, abatement is considered as a constant returns to scale activity and the pollution abated is assumed as an increasing and strictly concave function of both the total economic activity *F* and the economy's ecological effort, *F*^{A}. Then, pollution is defined as:

where *A* is the abatement level. The specification can be rewritten also as:

where $a(\theta )\equiv [1-A(1,\frac{{F}^{A}}{F}\left)\right]$ and $\theta =\frac{{F}^{A}}{F}$. Combining environmental assumptions with the Solow model, the output *Y* becomes $\stackrel{~}{Y}=(1-\theta )F$ and, assuming also an exogenous technological process in abatement at the rate *g*_{A} > 0, it is derived:

where $k=\frac{K}{BL}$, $\stackrel{~}{y}=\frac{\stackrel{~}{Y}}{BL}$, $p=\frac{P}{BL}$ and *f*(*k*) = *F*(*k*, 1).

From a dynamic point of view, the model presents a differential equation in *k* from which, by assuming that Inada conditions^{5} hold for *F* and θ fixed, the unique interior fixed point *k** attracts every initial condition *k*(0) > 0.

Moreover, it can be noticed that, as *k* → *k**, the aggregate output, consumption and capital approach the same growth rate *g*_{B} + *n* and, consequently, their per capita levels grow at the rate *g*_{B} > 0.

At the steady state level *k**, the growth rate of emissions *g*_{P} is given by:

In order to describe the dynamic relationship between capital accumulation and emissions level, the authors introduce a definition of sustainable growth in these terms: a balanced growth path^{6} is sustainable if it is associated to both rising consumption per capita and an improving environment. In mathematical terms, it is characterized by the following conditions:

By differentiating the Equation (6) and assuming *f*(*k*) = *k*^{α}, the dynamics of the model are described by the following system:

Therefore, the system (9) allows to capture the evolution of emissions along time, as *k* → *k**. By assuming sustainable growth (that is, *g*_{P} < 0), the value *k*_{T} (called turning point), ensuring $\frac{\u1e56}{P}=0$, is lower than *k** and this implies that the time profile of emissions depends on the position of the initial value *k*_{0} related to the turning point. In particular, if the economy starts with an initial capital stock such that *k*_{0} < *k*_{T}, emissions first rise and then fall. Hence, an inverted U-shaped profile, recalling the environmental Kuznets curve, is obtained. Instead, if an initial capital stock *k*_{0} > *k*_{T} is assumed, emissions monotonically fall as *k* moves toward its steady state value *k**.

Finally, when the sustainable growth is not assumed, emissions grow for all *t* even as *k* approaches the steady state value.

### 2.2. The Ramsey-Cass-Koopmans Model With Environmental Pollution

The model in Xepapadeas [5] represents a natural evolution of the approach proposed in the previous subsection. In this context, consumption-investment choices are derived in a decentralized framework composed by intertemporal utility maximizing agents and perfectly competitive profit maximizing firms. The individual utility function depends on the per capita consumption flow *c*(*t*) and the pollution stock *P*(*t*).

The representative consumer considers the pollution level as fixed and solves the problem:

where ρ is the discount rate in the utility function, *k*(0) is the initial capital and $R(t)={\int}_{\tau =0}^{t}r(\tau )d\tau $, with *r*(τ) real interest rate at time τ and *e*^{−R(t)} an appropriate discount factor. Under standard assumptions of concavity on *U* and imposing that $\underset{c\to +\infty}{lim}}{U}_{c}(c,P)=0$, the solution of the optimization problem, obtained with the use of the maximum principle, is interior (that is, *c*(*t*) > 0 for every *t*) and the consumption path is defined by:

where $\eta =-\frac{{U}_{cc}}{{U}_{c}}c$.^{7} By assuming that (i) the production function *f*(*k*) satisfies the Inada conditions, (ii) markets are competitive and firms are profit maximizers (which imply *f*′(*k*) = *r* + δ), the economic dynamics is described by the following system:

By analysing the model, a unique steady state (*c**, *k**, *P**) exists and the dynamics, generated by an optimization process, evolve on the stable manifold of (*c**, *k**, *P**) and converge to it. As explained in Xepapadeas [5], by considering that pollution evolves, it can be noted that the coordinates (*c**, *k**) are generally not affected by the stationary state value of *P**, although the approach path to the steady state is affected [due to the presence of $\frac{{U}_{cP}}{{U}_{c}}\u1e56$ in the first equation of system (11)].

## 3. Discrete Time Models

If, on the one hand, the aforementioned models are suitable to describe long run trends, on the other hand they are not able to capture the occurrence of short or medium run nonlinear phenomena.

Regarding this point, the work of Day [6] provides an example of how complex phenomena may emerge from the dynamic analysis of a simple economic framework.

### 3.1. The Day's (1982) Model

In this work, the author considers a neoclassical growth model à la Solow (1956) in discrete time. Indeed, assuming no capital depreciation, it can be expressed by a first order difference equation in the capital-labor ratio ${k}_{t}=\frac{{K}_{t}}{{L}_{t}}$, that is

where σ is the saving ratio assumed as constant, *f*(·) is the production function and λ is the natural growth rate of population.

By assuming that productivity is reduced by a so called pollution effect, Day [6] introduces the following specification for *f*:

where *A* > 0 is a productivity parameter, β > 0 is the elasticity of capital, *m* > 0 is the state of environment if private production is not performed and γ > 0 weighs the effects of pollution.

Then, the map that describes the dynamics of the model becomes:

The map in (14) is unimodal, *C*^{1}, and admits a unique maximum point $\overline{k}=\left(\frac{\beta m}{\beta +\gamma}\right)$, from which it can be obtained the maximum capital-labor ratio ${k}^{m}=f(\overline{k})$, that is:

In particular, note that the slope of the production function indefinitely grows as *k* approaches zero and consequently, for sufficiently small initial conditions *k*_{0} > 0, growth is positive.

For values of *A* sufficiently small, the stationary state is monotonically achieved either from above or below. As *A* increases, the steady state capital stock increases until verify the equality ${k}^{m}=\overline{k}={k}^{s}$, where *k*^{s} is the steady state capital-labor ratio. This represents the bifurcation point from which oscillations in levels of *k* occur for even higher values of *A*.

The condition *k*^{m} ≤ *m* is assumed to prevent negative values of *k*. Then, the author derives the following sufficient condition for the existence of growth cycles, that is:

Indeed, for a parameterizations satisfying such condition, the capital-labor ratio exhibits bounded oscillations (perhaps after a period of growth).

In order to determine if such cycles can be chaotic, Day [6] shows the existence of a parameter set for the map such that the theorem in Li and Yorke [7] is satisfied. For example, when β = γ = *m* = 1, the map (14) becomes:

from which, following [8], the existence of a parameter set such that:

implies that irregular growth cycles appear (see Figure 1).

**Figure 1**. Parameter set: σ = 0.5, λ = 0.02, *A*_{l} = 5, *A*_{int} = 6.5, *A*_{h} = 7.5. **(A)** Changes in the shape of the map, as *A* varies. **(B)** Convergence to the fixed point for *A* = *A*_{l} (*k*_{0} = 0.3). **(C)** Convergence to a 2-cycle for *A* = *A*_{int} (*k*_{0} = 0.21). **(D)** Chaotic regime for *A* = *A*_{h} (*k*_{0} = 0.3).

### 3.2. The Zhang's (1999) Model

In this model, the author considers an economy à la [9], in which at every period *t* there are two overlapping generations of individuals. The representative agent of each generation gets utility from consumption, *c*_{t+1}, and environmental quality, *E*_{t+1}, when old. The utility function *U*(*c*_{t+1}, *E*_{t+1}) is assumed to be increasing with respect to each argument and strictly concave. Furthermore, some Inada like conditions are satisfied, in order to get an optimal bundle with strictly positive values for *c* and *E*.

The individual supplies inelastically a unit of labor only in the first period, earning a real wage *w*_{t}, and optimally allocate *w*_{t} between private savings *s*_{t}^{8} for the old age consumption and environmental improvement expenditures, *m*_{t}. In the old age, the agent earns (1 + *r*_{t+1} − δ)*s*_{t}, where *r* and δ are the real rate of return and the capital depreciation rate, respectively. Then, the individual life-cycle budget constraints are given by:

The environmental quality evolves according to the following rule, staken from John and Pecchenino [9]:

where *b* ∈ (0, 1) represents the degree of autonomous evolution of environmental quality, β*c*_{t} measures the consumption degradation of the environment, and γ*m*_{t} measures the environmental improvement (β and γ are assumed to be positive). The production function *f*(*k*_{t}):ℝ_{+} → ℝ_{+} is assumed to be *C*^{2} and strictly concave with respect to *k*_{t} (the capital-labor ratio). Firms maximize profits and then:

The representative agent maximizes his utility with respect to (19) and (20). Then, the first order necessary and sufficient condition for the maximization problem reads as:

Hence, a perfect foresight competitive equilibrium is characterized by the Equation (21), the condition (22) and the market clearing condition *k*_{t+1} = *s*_{t}.

To simplify the analysis, the author assumes a constant elasticity of substitution between consumption and environment, η_{E}, defined as:

and this allows to rewrite equation (22) as

From the market clearing condition, the author gets the following relationship between capital and environment:

and then the study of the dynamics can be reduced to the analysis of the following first-order difference equation:

where $\alpha (k)\equiv \frac{k{f}^{\prime}(k)}{f(k)}$ is the capital share of output. In order to simplify the analysis, α is assumed as constant. Therefore, by considering $f({k}_{t})=A{k}_{t}^{\alpha}$ (where *A* > 0 represents the productivity parameter), [10] gets the expression

where ${a}_{0}\equiv \frac{(1-b){\eta}_{E}\gamma -\beta (1-\delta )}{\gamma (1+{\eta}_{E})}$. It can be noticed that this parameter is less than 1 and it is independent with respect to α.

About the dynamics of the model, an interior fixed point exists if and only if γ(1 − α) − βα > 0. In particular, Zhang [10] shows that, when *a*_{0} ≥ 0, the unique stationary equilibrium is given by

and the following proposition holds:

Proposition 1. *Suppose* γ(1 − α) − βα > 0. *Then, if* *a*_{0} ≥ 0, *map* *G* *is monotonically increasing; for all* *E*_{0} ∈ (0, +∞), $\underset{t\to +\infty}{lim}{G}^{t}({E}_{0})={E}^{*}$. *That is, there exists a unique and asymptotically stable (attracting) positive steady state*.

When *a*_{0} < 0, dynamics may exhibit complexity. By assuming $-\frac{{\alpha}^{-\frac{\alpha}{1-\alpha}}}{1-\alpha}<{a}_{0}$, [10] proves that if ${a}_{0}\in [-\frac{1+\alpha}{1-\alpha},0)$, *E** is stable. Instead, if *a*_{0} decreases in the interval $(-\frac{{\alpha}^{-\frac{\alpha}{1-\alpha}}}{1-\alpha},-\frac{1+\alpha}{1-\alpha})$, the map *G* exhibits the classical period-doubling sequence, summarized by the bifurcation diagram in Figure 2. It can be noticed that *A* does not play any role in the dynamic properties of the model: it only influences the level of *E* at the stationary state.

### 3.3. The Naimzada-Sodini (2010) Model

In the case of discrete time unidimensional models with overlapping generations, a variation of the model proposed in Zhang [10], with a more general production function, has been provided by Naimzada and Sodini [11]. Remaining close to the paradigm introduced by John and Pecchenino [9], the authors propose a model in which a population of individuals is characterized by a utility function depending on the stock of an environmental good, *E*_{t}, and on the consumption of the private good, *c*_{t}. It is also assumed that *E*_{t} is negatively affected by the consumption but it is improved by specific environmental expenditures.

Differently from Zhang [10], Naimzada and Sodini [11] consider allocations in a decentralized economy. In this case, the choices of each agent between consumption and environmental expenditure generate externalities on the others. A consequence of this assumption is that the environmental maintenance, because of the environment is considered as a public good, is characterized by the classic free-riding problem. In particular, the functional form described by the authors is the following:

where *b* ∈ (0, 1) measures the autonomous evolution of environmental quality, β > 0 measures the consumption effect on environment, γ > 0 weighs the environmental expenditures efficiency, *N* is the number of the agents, and $\overline{E}>0$ represent the long run value of the environmental index in absence of anthropic activity. Individual preferences are described by the utility function *U*(*c*_{t+1}, *E*_{t+1}), with *U* assumed as twice continuous and differentiable. Also in this case, Inada like conditions are also assumed to avoid corner solutions in the optimization problem.

On the production side, compared with [10], the model considers a more general specification, that is a CES-technology:

where *k*_{t} is the physical capital level at *t*, *A* > 0 is a scaling parameter, α ∈ (0, 1) measures the degree of capital intensity of production, and $\theta =\frac{1}{1+\rho}$ represents the elasticity of substitution between labor and capital, with ρ > 0.

From the optimality conditions, equilibrium expressions of wage rate *w*_{t}, and interest rate *r*_{t} are derived as

By assuming that agents are identical, the problem faced by the individual born at *t* is:

in which

represents the expected environmental quality at time *t* + 1 (depending on ${m}_{t}^{e}$, that is, the expectation of agent *i* about strategies of others *N* − 1 identical agents) and the following constraints apply:

where *s*_{t} is the saving and δ represents the depreciation rate of capital.

By assuming that the young individual at *t* is able to perfectly foresee the environmental index *E*_{t+1} and recalling the Equations (31) - (32), the equilibrium conditions for all *t* become:

In order to analyse dynamics of the model, the authors, as in Zhang [10], introduce a constant elasticity of substitution between consumption and environment, defined as:

Hence, they characterize the intertemporal equilibrium conditions by means of a nonlinear difference equation in *E*_{t}:

Therefore, Naimzada and Sodini [11] describe the dynamics distinguishing between two possible cases, depending on the sign of ρ: (i) ρ < 0 and (ii) ρ > 0.

In the case (i), the authors show that the map admits a unique positive fixed point and highlights that qualitative results similar to the ones in Zhang [10] are achievable. Then, the unique fixed point *E** could be attracting or repelling and, in the latter case, limit cycles or a chaotic attractor arises. However, differently from Zhang [10], the scaling parameter *A* affects the stability of *E**. Indeed, by starting from the parameterizations $\alpha =0.831,\beta =0.55,\gamma =1.56,\delta =0.001,{\eta}_{E}=5,\rho =-0.4$, $A=3,b=0.52,\overline{E}=1,N=2,000$, an increase in *A* generates a loss of stability of *E** through a period-doubling bifurcation. In addition, for *A* ≃ 6, the map exhibits chaotic dynamics.

In the case (ii), the map admits an odd number of fixed points and the ones with even index are unstable. Then, differently from (i), multiple equilibria may exist. In particular, the authors prove that, for ρ > > 0, three steady states exist.

An other interesting phenomenon shown by the authors is the following. By considering a parameter set for which a unique fixed point ${E}_{1}^{*}$ exists, an increase of *A* is able to generate, through a fold bifurcation, two new steady states ${E}_{2}^{*}$ and ${E}_{3}^{*}$, where ${E}_{2}^{*}<{E}_{3}^{*}$ is repelling and separates the basins of attraction of the two attracting fixed points ${E}_{1}^{*}$ and ${E}_{3}^{*}$. Then, a change in *A* may be engine of a poverty trap (see [12]).

Considering another configuration of parameters (see the parameter set in Figure 3), the authors provide another example in which a different sequence of dynamic phenomena occurs, as *A* varies.

**Figure 3**. Parameter set: $\alpha =0.07,\beta =0.1,{\eta}_{E}=11,\rho =12,N=100,\delta =0.4$, $\gamma =0.04,b=0.58,\overline{E}=22$. The part of the bifurcation diagram depicted in black is generated starting from *E*_{0} = 0.6. The part depicted in red shows the existence of a second attractor for *A* ∈ [0.83, 0.88] (the initial condition is *E*_{0} = 0.29).

For *A* low, a unique repelling fixed point *E*^{1}, enclosed in an attracting limit 2-cycle, exists; by increasing *A*, the map undergoes the classic period-doubling sequence until a fold bifurcation generates two new fixed points, *E*^{2} (repelling) and *E*^{3} (attracting) with *E*^{3} < *E*^{2} < *E*^{1}. When *A* further increases, first ${E}_{3}^{*}$ loses its stability through a flip bifurcation and then a cascade of flip bifurcations appears. Moreover, there exists a region of parameters for which coexistence of attractors occurs (see Figure 3). For larger values of *A*, the lower attractor dies and all the feasible trajectories are attracted by the remaining attractor that lives at the right of the repelling steady state *E*^{2}.

Figure 4 shows an interesting global phenomenon analyzed by the authors, for even higher values of *A*: the remaining attractor enlarges and invades the space before occupied by the other attractor.

Finally, an uncommon phenomenon analyzed in the work is the following one. By fixing $\alpha =0.07,\beta =0.1,{\eta}_{E}=11,\rho =12$, $N=100,\delta =0.4,\gamma =0.04,b=0.58,\overline{E}=22$, the second iterate of the map is characterized by two humps or more precisely by a maximum and a minimum point. By investigating the second iterate *G*^{2}, it can be noticed that an increase of *A* may lead to a fold bifurcation of *G*^{2} inducing a stable 2-cycle and an unstable 2-cycle. In this case, the fold bifurcation arises far from the fixed point that maintains its instability for the whole process (see Figure 5).

**Figure 5**. Figure reproduced with the permission of the authors. Source: Naimzada and Sodini [11]. Evolution of the second iterate of *G* (drawn at the bifurcation value *A* = 0.132). For *A* = 0.10, *G*^{2} has no intersections with the main diagonal (dashed line); for *A* = 0.132, *G*^{2} creates 2 tangent points (solid line) and for *A* = 0.16, two fixed points for *G*^{2} are born (dotted line).

### 3.4. The Caravaggio-Sodini Model (2017)

Recently, by following some suggestions from empirical findings^{9}, Caravaggio and Sodini [14] have proposed a model in which, by assuming an overlapping generation framework and discrete time, the environmental resource is not an argument in the agent's utility but is taken as an input in the production of a private good. Moreover, differently from John and Pecchenino [9], the model on the one hand assumes that the economic activity damages the environment and that agents consider it only as an externality^{10}; on the other hand, it assumes that the government imposes a wage taxation with a fixed share devolved to the maintenance of a natural resource^{11}.

About the economic agent, at time *t* he supplies inelastically his time endowment and the earned wage *w*_{t} is allocated between the consumption *c*_{t} and the saving *s*_{t}. Hence, the individual preferences depend on consumption in the young age, *c*_{t}, consumption in the old age, *d*_{t+1} = *s*_{t}*R*_{t+1} and on the level of basic services provided by public sector. Thus, the utility function is given by:

where ρ, σ ∈ (0, 1) represent the intensity of consumption spillover in young and old age, respectively, ϕ > 0 determines the discount factor and ${\overline{c}}_{t}$, ${\overline{d}}_{t+1}$ are consumption externalities. Moreover, the government levies a tax on labor wage at the tax rate τ ∈ (0, 1). Saving, remunerated at the real interest factor *R*_{t+1}, is used to consume the final good in the old age. Then, the agent has to face the following intertemporal budget constraint:

Maximizing the utility function (42) with respect to (43) and considering the ex post conditions ${\overline{c}}_{t}={c}_{t}$, ${\overline{d}}_{t+1}={d}_{t+1}$, the saving function is obtained as:

On the production side, the private good is produced through the following Cobb-Douglas technology:

where ${\overline{k}}_{t}$ represents a production externality, *A* > 0 is the productivity scaling parameter and *E*_{t} is the stock of natural resource involved in the productive process. In that specification, authors assume both the possibility of decreasing social returns to scale and constant social returns to scale, that is α + β ≤ 1. From the known optimality conditions, the equilibrium expressions for wage and interest rate are derived as:

Regarding to the evolution of natural resource, a functional form similar to the one introduced in Smulders and Gradus [16] is adopted:

where $\overline{E}>0$ represents the natural level of environmental resource, that is the level when no productive activity exists, and λ, θ > 0 weigh the environmental impact of production *Y*_{t} and the impact of public environmental expenditures *G*_{t} = δ·(τ*w*_{t}).

By considering the aggregate consistency condition ${k}_{t}={\overline{k}}_{t}$ and the market clearing condition on the capital market *k*_{t+1} = *s*_{t}, the dynamics of the model are driven by the following dynamic system:

Introducing the change of variables:

the dynamics of the model can be studied through the analysis of the following one-dimensional map:

where *H*(*Y*_{t}) is defined as:

In order to analyse the dynamics of the system by studying the dynamic properties of the map *N*, some relationships between maps *M* and *N* are provided:

(a) 0 is a fixed point for *N* if and only if $(0,\overline{E})$ is a fixed point for *M*;

(b) *Y** > 0 is a fixed point for *N* if and only if there exist *k** and *E** such that (*k**, *E**) is a fixed point for *M*;

(c) an attracting fixed point for *N* corresponds to an attracting fixed point for *M* and vice versa.

For both mathematical and economic reasons, Caravaggio and Sodini [14] distinguish the two cases (i) α + β < 1 and (ii) α + β = 1, highlighting how different parameterizations and different initial conditions may lead the economy on a bounded growth path, on an unbounded growth path or on a poverty trap path.

In particular, in the case (i), a unique interior globally stable fixed point *Y** or three interior fixed points ${Y}_{1}^{*}<{Y}_{2}^{*}<{Y}_{3}^{*}$ with ${Y}_{1}^{*}$ and ${Y}_{3}^{*}$ locally stable and ${Y}_{2}^{*}$ unstable are admitted. Therefore, dynamics of *N*, for initial conditions on the left and the right of ${Y}_{2}^{*}$, are captured by ${Y}_{1}^{*}$ and ${Y}_{3}^{*}$ (see Figure 6A). In addition, investigating further properties of the relationship between *M* and *N*, it can be observed that starting from a positive stationary equilibrium *Y** for *N*, an invariant set in the plane (*k, E*), defined by $E=(\frac{{Y}^{*}}{A{k}^{\alpha +\beta}}{)}^{\frac{1}{\gamma}}$, is identified. Such expression, for the parameter set considered, describes the stable manifold of the saddle point $({k}_{2}^{*},{E}_{2}^{*})$ that separates the basins of attraction of $({k}_{1}^{*},{E}_{1}^{*})$ and $({k}_{3}^{*},{E}_{3}^{*})$, as depicted in Figure 6B.

**Figure 6**. Figure reproduced with the permission of the publisher. Source: Caravaggio and Sodini [14]. Parameter set: $\alpha =0.32,\beta =0.6,\gamma =0.4,\delta =0.7,\theta =3,\lambda =2.87$, $\rho =0.25,\sigma =0.7,\tau =0.2,\varphi =0.6,\overline{E}=0.0007,A=48.6$. **(A)** Time series of *Y* converging to ${Y}_{1}^{*}$ or to ${Y}_{3}^{*}$ for different initial conditions; **(B)** Dashed lines describe iso-product curves in the plane (*k, E*) associated with stable fixed points of *N*, ${Y}_{1}^{*}$ and ${Y}_{3}^{*}$. The iso-product curve related to the unstable ${Y}_{2}^{*}$ (solid line) is the stable manifold of the saddle $({k}_{2}^{*},{E}_{2}^{*})$ and separates the basin of attraction for $({k}_{1}^{*},{E}_{1}^{*})$ (green region) and $({k}_{3}^{*},{E}_{3}^{*})$ (yellow region), respectively.

Authors also highlight the particular role played by the scaling parameter *A*. Indeed, it can be noted that a decrease in *A* may cause the disappearance of the fixed points ${Y}_{2}^{*}$ and ${Y}_{3}^{*}$. Then ${Y}_{1}^{*}$ becomes the unique attractor of the system and a successive recovery of *A* at its original value may not allow the restoration of the previous path. Figures 7A,B provide a numerical example of this occurrence in terms of *k* and *E*.

**Figure 7**. Figure reproduced with the permission of the publisher. Source: Caravaggio and Sodini [14]. Parameter set: $\alpha =0.32,\beta =0.6,\gamma =0.4,\delta =0.7,\theta =3,\lambda =2.87,\rho =0.25$, $\sigma =0.7,\tau =0.2,\varphi =0.6,\overline{E}=0.0007,{A}_{h}=48.6,{A}_{l}=48.3$. **(A)** Evolution of capital accumulation over time when a shock in *A* occurs; **(B)** Evolution of the stock of environmental resources over time when a shock in *A* occurs. The red pointed vertical line represents the temporary (exactly 10 iterations) time series generated by *A*_{l}.

In the case (ii), Caravaggio and Sodini [14] underline two interesting results: (1) the noninvertibility and unimodality of the map may induce the birth of nonconnected basins of attraction (see Figure 8);

**Figure 8**. Figure reproduced with the permission of the publisher. Source: Caravaggio and Sodini [14]. Parameter set: $\alpha =0.6,\beta =0.4,\gamma =0.4,\delta =0.1,\theta =1.12,\lambda =8.92,\rho =0.1$, $\sigma =0.2,\tau =0.21,\varphi =0.19,\overline{E}=0.000001,{A}_{h}=5,000,{A}_{int}=250,{A}_{l}=150$. For *A* = *A*_{int}, the map admits three stationary equilibria: 0 (stable), ${Y}_{1}^{\star}=0.444$ (unstable) and ${Y}_{2}^{\star}=1.057$ (stable). The basin of 0 is nonconnected and consists of two parts: the immediate basin in $(0,{Y}_{1}^{*})$ and the other part (*Y*_{p}, +∞), where *Y*_{p} is the preimage of ${Y}_{1}^{*}$. On the contrary, the basin of ${Y}_{2}^{*}$ is connected being a single interval $({Y}_{1}^{*},{Y}_{p})$.

(2) For λ < θ, it is possible that trajectories, after a long transient, positively diverge. On the contrary, for λ > θ trajectories, after the transient phase, may converge to the unique attractor 0, that represents a nonproductive state for the economy (see Figures 9A,B, respectively).

**Figure 9**. Figure reproduced with the permission of the publisher. Source: Caravaggio and Sodini [14]. **(A)** Positive divergence when λ < θ; **(B)** convergence to a nonproductive state when λ > θ.

### 3.5. The Antoci-Gori-Sodini Model (2016)

Continuing with the literature on the discrete time OLG models, several works in which a different allocation problem for the agents is taken into account have been proposed, as in Antoci and Sodini [17] and in Antoci et al. [18]. In this context, differently from John and Pecchenino [9], where the investment in environmental defensive expenditures is internalized by agents in the allocation problem, economic agents are assumed to be unable of interven for the environmental maintenance. These models therefore consider an environmental quality that is deteriorated by economic activity and individuals who defend themselves from such worsening by simply increasing the consumption of the private good, perceived as a substitute of the environmental one. As shown in Antoci and Sodini [17], since the production has a negative impact on the environment, self-protection choices induce further environmental degradation. Then, a self-enforcing mechanism may be observed according to which environmental depletion leads to an increase in the consumption of private good, which in turn generates further environmental degradation and so on.

These works show also that such self-reinforcing substitution mechanism may fuel undesirable growth paths according to which an increase in consumption of the private good is associated with a reduction in the individual well-being. Moreover, this mechanism may lead to the occurrence of several global phenomena and the individuals' expectations about the future evolution of the environmental quality can induce local and global indeterminacy.^{12}

Before introducing the main features of the model in Antoci et al. [18], it can be useful to recall the definition of indeterminacy (local and global). Indeed, indeterminacy occurs if, for a given initial condition on the state variables, different values for the control variable can be chosen. In particular, if such values lead to the same Ω-limit set, then local indeterminacy occurs, otherwise if the different initial values of the control variable lead to different Ω-limit sets, then global indeterminacy arises. On the contrary, a steady state is said to be (a) locally determinate if, for an initial condition on the state variables close to the steady state, there exists, in a neighborhood of the steady state, only one feasible choice for the control variables and (b) globally determinate if, for a given initial condition on the state variables, there exists only one value of the control variables that can be chosen.

Moving on the model, the authors consider an OLG economy composed of identical and perfectly rational individuals who live for two periods, youth and old age. Each young agent allocates his time endowment between labor *l*_{t} ∈ (0, 2), remunerated at the wage rate *w*_{t}, and leisure 2−*l*_{t}. Then, the representative member of generation *t* earns a labor wage *w*_{t}*l*_{t} which is entirely saved at time *t*, that is *s*_{t} = *w*_{t}*l*_{t}. The saving *s*_{t} is invested in productive capital (*K*_{t+1} = *s*_{t}), that the individual rents to the representative firm at the interest factor *R*_{t+1}. The revenue *R*_{t+1}*s*_{t} is used at the second period of life to buy and consume the quantity *C*_{t+1}, so that *C*_{t+1} = *R*_{t+1}*s*_{t}. Therefore, the lifetime budget constraint of the individual of generation *t* is given by:

Individual preferences depend on leisure when individuals are young and consumption when individuals are old. In addition, the authors assume that the utility function is affected by the environmental quality, *E*_{t+1}. Then, the lifetime utility index is represented by a constant inter-temporal elasticity of substitution (CIES) function, that is:

where σ > 0 (σ ≠ 1) and γ > 0 (γ ≠ 1) determine the elasticity of utility with respect to consumption and leisure, respectively, and τ is a positive parameter reflecting the relative influence of environmental quality in defining the individual utility. Moreover, if σ ∈ (0, 1) (resp. σ > 1) then consumption and environment are complements (resp. substitutes). In particular, for σ ∈ (0, 1) (resp. σ > 1) the following derivative holds:

^{13}

The representative individual of generation *t* chooses labor supply *l*_{t} in order to maximize the value of utility function (54), subject to the intertemporal budget constraint (53) and the condition *l*_{t} ∈ (0, 2). Then, in order to get an interior solution, the necessary and sufficient first order condition for the maximization is:

On the production side, the authors assume that the firm produces output *Y*_{t} through a constant returns to scale Cobb-Douglas technology, that is:

where *A* > 0 is a scaling parameter and α ∈ (0, 1) represents the elasticity of the production function with respect to capital. By assuming that the capital stock *K*_{t} fully depreciates at the end of every period, the profit maximization implies that productive factors are remunerated at their marginal products, that is:

With regard to the environmental resource, the authors define the environmental quality index as a function that negatively depends on the economy-wide average production at time *t*, that is:

where $\overline{E}$ is a positive parameter symbolizing the natural value of the environmental quality, and ρ > 0 measures the environmental impact of production. It can be noticed that this specification eliminates the accumulation of environmental damages provoked by the production and, from a mathematical point of view, allows the authors to investigate the dynamical system in the plane (*K*_{t}, *l*_{t}).

By imposing the market clearing condition *K*_{t+1} = *s*_{t} = *w*_{t}*l*_{t}, assuming that individuals have perfect foresight and that the economy-wide average output coincides, ex post, with the output of the representative firm, equilibrium conditions can be described by the following map:

where *M* is defined in the set $D:=\left\{({K}_{0},{l}_{0})\in {\mathbb{R}}^{2}:{K}_{t}>0,0<{l}_{t}<2,\forall t>0\right\}$.

Regarding to existence, multiplicity and stability of fixed points, the map admits always at least one fixed point in the set *D*. In particular, if 0 < σ < 1, the map admits a unique fixed point, while if σ > 1, the map admits one or three fixed points, according to the natural value of the environmental quality, $\overline{E}$. Furthermore, the authors prove that for 0 < σ < 1 the fixed point is a saddle (globally determined), independently of the parameterizations; on the contrary, for σ > 1 two cases arise: (i) there exists a unique fixed point that can be a sink or a saddle independently of the parameterizations, or (ii) there exist three fixed points which have not any dynamic properties independent of the applied parameterizations. In the case (ii), given the three fixed points $({K}_{1}^{*},{l}_{1}^{*})$, $({K}_{2}^{*},{l}_{2}^{*})$ and $({K}_{3}^{*},{l}_{3}^{*})$, $({K}_{1}^{*},{l}_{1}^{*})$ and $({K}_{3}^{*},{l}_{3}^{*})$ can be sink or saddle whereas $({K}_{2}^{*},{l}_{2}^{*})$ can be source or saddle.

An important result shown by the authors is that, due to the absence of externalities on the production and the inclusion of the environmental variable only in the agents' utility, the map describing the dynamics of the model is invertible. This implies that, differently from Caravaggio and Sodini [14], basins of attraction are connected sets. Furthermore, the inverse map allows to numerically approximate the stable manifold of saddles, that is the trajectories converging to a determinate equilibrium.

The main issue investigated in the work is the possible rise of local and global indeterminacy. By recalling that in the model *l*_{t} is the control variable (that is, the initial value *l*_{0} is fixed after the optimization process), while *K*_{t} is the state variable, a stationary equilibrium (*K*_{t}, *l*_{t}) is locally indetermined if it is an attractor from a dynamic point of view (in other words, it holds that for every eigenvalue λ_{i}, |λ_{i}| < 1). On the contrary, (*K**, *l**) is locally determined if it is a saddle, that is one of the two eigenvalues (in modulus) is larger than 1 and the other (in modulus) is smaller than one. The economic dynamics are said to be globally indetermined if there is coexistence of attractors and/or saddles.

The main configurations of the dynamics can be summarized through four examples. In the example 1, the parameter set is $A=1.25,\overline{E}=12.2,\alpha =0.2,\gamma =1.12,\sigma =3,\tau =1$. In this case, for 0 < ρ < 1 there exists a unique fixed point $({K}_{1}^{*},{l}_{1}^{*})$ and it is a saddle (then locally determined). For ρ ≃ 1, $({K}_{1}^{*},{l}_{1}^{*})$ becomes stable and it implies local indeterminacy. This means that, given an initial value *K*_{0} sufficiently close to ${K}_{1}^{*}$, there exists a continuum of initial conditions *l*_{0} such that trajectories starting from (*K*_{0}, *l*_{0}) converge to $({K}_{1}^{*},{l}_{1}^{*})$ (see Figures 10A,B).

**Figure 10**. Figure reproduced with the permission of the publisher. Source: Antoci et al. [18]. **(A)** Stable manifold of the locally determined fixed point $({K}_{1}^{*},{l}_{1}^{*})$, obtained for ρ = 0.2. **(B)** Basin of attraction (in gray) of the locally indeterminate fixed point $({K}_{1}^{*},{l}_{1}^{*})$ for ρ = 1.

A further increase in ρ generates a saddle-node bifurcation from which $({K}_{1}^{*},{l}_{1}^{*})$ and $({K}_{3}^{*},{l}_{3}^{*})$ (attracting) and $({K}_{2}^{*},{l}_{2}^{*})$ (saddle) are obtained. Therefore, global indeterminacy arises and then there exists initial values *K*_{0} from which it is possible to reach all coexisting fixed points by opportunely choosing the initial value *l*_{0}, as shown in Figure 11. By further increasing the value of ρ, simultaneously $({K}_{1}^{*},{l}_{1}^{*})$ approaches the saddle $({K}_{2}^{*},{l}_{2}^{*})$ and $({K}_{3}^{*},{l}_{3}^{*})$ loses its stability in favor of 2-cycle via a supercritical flip bifurcation. After the collision between $({K}_{1}^{*},{l}_{1}^{*})$ and $({K}_{2}^{*},{l}_{2}^{*})$ almost all trajectories belonging to *D* are captured by the 2-cycle arisen from $({K}_{3}^{*},{l}_{3}^{*})$. Furthermore, the authors show that, in this case, (i) by starting from *K*_{0} sufficiently far from ${K}_{{3}^{*}}$, there exist several *l*_{0} such that trajectories converge to the locally determined $({K}_{3}^{*},{l}_{3}^{*})$ and (ii) given *K*_{0}, there exist other values of *l*_{0} for which trajectories converge to the 2-cycle. In the latter case, global indeterminacy occurs in a context in which a unique locally determined fixed point exists, as shown also in Coury and Wen [20].

**Figure 11**. Figure reproduced with the permission of the publisher. Source: [18]. Global indeterminacy (ρ = 3.2): two attracting fixed points, $({K}_{1}^{*},{l}_{1}^{*})$ and $({K}_{3}^{*},{l}_{3}^{*})$, coexist with the saddle $({K}_{2}^{*},{l}_{2}^{*})$.

In the example 2, considering the same parameter set of the example 1 but setting σ = 3.274, if ρ increases over the value 3.1, the system undergoes a sequence of flip bifurcations around $({K}_{3}^{*},{l}_{3}^{*})$ (see Figure 12A). Figure 12B shows, when ρ = 3.2465, the coexistence of $({K}_{1}^{*},{l}_{1}^{*})$ with a four-piece chaotic attractor. Because of the structure of the basins of attraction, predicting the long term behavior of the economy starting in the North-East portion of the state plane becomes rather difficult.^{14}

**Figure 12**. Figure reproduced with the permission of the publisher. Source: Antoci et al. [18]. **(A)** Bifurcation diagram for ρ. Due to the rise of attractors with rather tangled basins of attraction, the initial conditions were opportunely adapted to the change in ρ in order to follow the evolution of each attractor. **(B)** Basins of attraction when the attracting fixed point $({K}_{1}^{*},{l}_{1}^{*})$ coexists with a four-piece chaotic attractor (ρ = 3.2465).

Finally, In examples 3 and 4, the authors show the occurrence of quasi-periodic dynamics generated via Neimark-Sacker bifurcations, further phases of coexistence and the possible rise of the phenomenon of hysteresis.

## 4. Continuous Time Models

The models discussed in the previous section have been sketched the main strands of literature on discrete time dynamic modeling with environment. In parallel, the literature has been enriched by several contributions developed in a continuous time framework.

For example, the existence of growth dynamics affected by the substitutability between environmental depletion and consumption of a private good, highlighted in Antoci et al. [18], have been previously analyzed in Antoci [21] in a continuous time model. In this work, the author describes growth dynamics in an economy in which the individual well-being depends on leisure, an environmental resource and a private good that can be consumed as a substitute for the environmental one. Moreover, it is assumed that the production and consumption of the substitute good impose negative externalities through the environmental deterioration. The author shows that the ratio between the technical progress and the environmental quality plays a crucial role. In particular, different dynamic systems rule the equilibrium dynamics of the model depending on such ratio and this can provoke the birth of a bistable regime. Furthermore, the work anticipates the results provided in Antoci et al. [18], in which the substitutability mechanism is engine of undesirable growth paths.

In the same way, the dynamic effects induced by the inclusion of the environmental input in the production function, analyzed in Caravaggio and Sodini [14], have been discussed through a continuous time model in Antoci et al. [2]. The results of this last contribution are described in more detail in the next subsection.

### 4.1. The Antoci-Russu-Galeotti Model (2011)

In this work, the authors define a growth model in which negative environmental externalities are considered and where the environment is also assumed as a productive input.

The economy is assumed to be composed by a continuum of identical agents and the representative individual produces an output *Y*(*t*) through the following Cobb-Douglas technology:

where *K*(*t*) is the stock of physical capital accumulated, *L*(*t*) is the agent's labor input and *E*(*t*) is the stock of a natural resource employed as productive input.

It is assumed an additively nonseparable utility function, as in Bennett and Farmer [22], which depends on leisure 1 − *L*(*t*) and consumption of the produced good, *C*(*t*):

where η denotes the inverse of the intertemporal elasticity of substitution in consumption and ϵ, η > 0 with η ≠ 1. Moreover, the function is assumed to be concave with respect to *C* and 1 − *L*, and then the condition $\eta >\frac{\u03f5}{1+\u03f5}$ holds.

The evolution of *K*(*t*) is represented by the differential equation:

and the dynamics of the environmental resource are modeled with the well-known logistic equation augmented by the negative impact of production process:

in which $\overline{E}>0$ is the carrying capacity of the natural resource, $\overline{Y}$ represents the economy-wide average output and the parameter δ > 0 measures the negative impact of production on *E*.

The authors assume that the representative agent chooses *C*(*t*) and *L*(*t*) in order to solve the problem

where *K*(0) and *L*(0) are fixed, *K*(*t*), *L*(*t*), *C*(*t*)≥0, *L*(*t*) ∈ [0, 1] and θ > 0 is the discount factor.

Since agents are identical, the ex post condition $\overline{Y}=Y$ holds. This implies that trajectories generated by the model are not socially optimal but Nash equilibria.

Starting from the Hamiltonian function associated to the optimization problem:

in which Ω is the co-state variable associated to *K*, and solving the optimization problem, the authors derive the following dynamic system:

where $f(L)=\frac{\u03f5}{\beta}{(1-L)}^{\frac{\u03f5-\eta (1+\u03f5)}{\eta}}{L}^{1-\beta}={K}^{\alpha}{E}^{\gamma}{\Omega}^{\frac{1}{\eta}}$.

Regarding to the existence of stationary states, the authors show that the system admits (i) zero, one or two stationary states if α + γ < 1, (ii) one or zero stationary states if α + γ = 1, (iii) one stationary state if α + γ > 1.

In the case (i), the authors focus on the instance in which two stationary states ${P}_{1}^{*}$ and ${P}_{2}^{*}$ exist (with ${E}_{1}^{*}<{E}_{2}^{*}$ and ${K}_{1}^{*}<{K}_{2}^{*}$). This configuration is characterized by the assumption:

where $\stackrel{~}{E}$ is the unique positive value satisfying the condition ${g}^{\prime}(\stackrel{~}{E})=0$.

Analysing the Jacobian matrix *J* at ${P}_{1}^{*}$ and ${P}_{2}^{*}$ and using the Routh-Hurwitz criterion [23] to establish the number of eigenvalues with positive (or negative) real part, configurations of parameters for which ${P}_{1}^{*}$ is a sink, whereas ${P}_{2}^{*}$ is a saddle with a two-dimensional stable manifold are obtained.

By taking $\overline{E}$ as the bifurcation parameter, as $\overline{E}$ increases, ${P}_{2}^{*}$ does not change its nature, while ${P}_{1}^{*}$ can undergo one, two or no Hopf bifurcations.

These results show how there exist sufficient conditions to get local indeterminacy, and they are affected by the value of η,^{15} high enough values of δ and low enough values of θ (see the numerical example in Figure 13 where, given two initial values for the state variables (*K*_{0}, *E*_{0}) sufficiently close to their stationary values in ${P}_{1}^{*}$, there exists a continuum of initial values *l*_{0} such that the trajectories are attracted by the limit cycle around to ${P}_{1}^{*}$).

**Figure 13**. Parameter set: $\alpha =0.1,\beta =0.8,\gamma =0.58,\delta =0.05,\u03f5=1,\eta =1.5,\theta =0.001,\overline{E}=0.21$. Locally attracting limit cycle around the unstable stationary point ${P}_{1}^{*}$ depicted in red, arisen via Hopf bifurcation.

In the case (iii), where *P** is the unique stationary state, the trace of the Jacobian matrix *J** is negative for some values of $\overline{E}$, and *P** is a saddle with a two-dimensional stable manifold. By increasing the value of $\overline{E}$ (taken as bifurcation parameter) three chances can be obtained: (a) *P** remains a saddle with a two-dimensional stable manifold and no Hopf bifurcations occur; (b) *P** remains a saddle with a two-dimensional stable manifold but two Hopf bifurcations occur; (c) *P** becomes a source and only one Hopf bifurcation occurs.

Regarding to the global analysis, the main result shown by the authors is the emergence of global indeterminacy. As a matter of fact, in the case (i) with ${P}_{1}^{*}$ sink and ${P}_{2}^{*}$ saddle, there exists a neighborhood *N* of ${P}_{1}^{*}$ such that, for every (*E*_{0}, *K*_{0}, *L*_{0}) ∈ *N*, the half-line {*E* = *E*_{0}, *K* = *K*_{0}, *L* < *L*_{0}} intersects the stable manifold of ${P}_{2}^{*}$. This means that, for every initial couple (*E*_{0}, *K*_{0}) sufficiently close to $({E}_{1}^{*},{K}_{1}^{*})$, there exists a continuum of initial values ${L}_{0}^{1}$ for which the trajectory starting from $({E}_{0},{K}_{0},{L}_{0}^{1})$ converges to ${P}_{1}^{*}$ and a locally unique ${L}_{0}^{2}$ such that the trajectory starting from $({E}_{0},{K}_{0},{L}_{0}^{2})$ converges to ${P}_{2}^{*}$ (see Figure 14). Therefore, in this case global indeterminacy occurs.

**Figure 14**. Parameter set: $\alpha =0.1,\beta =0.8,\gamma =0.58,\delta =0.05,\u03f5=1,\eta =1.5,\theta =0.001,\overline{E}=0.17$. Global indeterminacy in the space (*E, K, L*). The trajectory in green approaches the sink ${P}_{1}^{*}$ (depicted in blue) starting from the point $({E}_{1}^{*},{K}_{1}^{*},{L}_{0}^{1}=0.3745910588598)$, while the trajectory in black, approaching the saddle ${P}_{2}^{*}$ (depicted in red), starts from the point $({E}_{1}^{*},{K}_{1}^{*},{L}_{0}^{2}=0.345910588598)$.

## 5. Delayed Dynamics

The interplay between economic activity and environmental quality can be finally analyzed through dynamic models where, although the time is assumed to be continuous, a time delay is introduced to underline the delayed effect of a variable on the others. Specifically, in a context in which the environmental quality is assumed, this approach turns out to be useful in showing both the noninstantaneous environmental effect of economic activity and the impact of the environmental evolution on the economy.

A seminal contribution on this issue is represented by the model described in Matsumoto and Szidarovszky [24], in which a fixed time delay is introduced in the productive sector (that is, a time to built technology is considered) and a negative environmental impact of production is assumed. The authors show that, in an extension of Day [6], chaotic dynamics may arise if a sufficiently large time delay is assumed. In particular, dynamics are described by the following equation:

where *x* is the capital per labor, α > 0 represents the depreciation ratio of capital and *h* defines the production function in which the fixed delay *T* > 0 is included.

Another significant contribution is given by Ferrara et al. [25]. In this work, the authors analyse a simplified version of the Green Solow model [4] in which the production negatively affects the environmental evolution but the environmental quality does not generate consequences on the production side. The model is defined by the following dynamic system

where *k*_{d} = *k*_{t−τ} is the individual's capital stock considering the introduction of a time lag τ > 0, *p* represents the stock of pollution and α ∈ (0, 1) is the capital's share. Moreover, in the first equation, *s* and δ > 0 denote the saving ratio and the capital depreciation rate, respectively, whereas in the second equation ϕ is the amount of emissions per unit of output and *m* > 0 describes the pollution decay. In the analysis, the authors provide conditions under which cyclical dynamics in the economic sector may occur and investigate how these ones may affect the environmental dynamics.

In Ferrara et al. [26], following the spirit of models provided in Wirl [27] and Antoci et al. [2], the authors consider an economy composed by agents that are not able to internalize the effect of their choices on the environmental dynamics (although they consider the environmental depletion in the utility function) but that perfectly foresee the evolution of economic variables. The novelty of the model is that the dynamics of the environmental good is assumed to be negatively affected by the economic activity that took place in an earlier date (*t* − τ).

More in depth, the model assumes that the representative agent produces an output *Y*(*t*) through the Cobb-Douglas technology:

where α ∈ (0, 1) and *K*(*t*) represents the stock of physical capital employed in the production process (for simplicity, *L*(*t*) = 1 is assumed). From the environmental point of view, the following linear specification is assumed:

where $\overline{Y}(t-\tau )$ is the economy-wide average output at time *t* − τ, $\overline{E}>0$ represents the carrying capacity of the natural resource, β > 0 measures the speed of convergence of the environmental quality toward its carrying capacity in absence of production and γ > 0 weighs the negative environmental impact of production.

The utility function *U*(*C*(*t*), *E*(*t*)) is dependent on the consumption of a private good and the environmental quality. Then, the following constant intertemporal elasticity of substitution (CIES) formulation is considered:

where σ > 0 (σ ≠ 1) represents the inverse of intertemporal elasticity of substitution in consumption. Such specification allows to capture both cases of complementarity (σ < 1) and substitutability (σ > 1) between consumption of the private good and environmental degradation.

By assuming no depreciation of capital, the dynamics of *K*(*t*) is given by the differential equation

Therefore, the representative agent faces the following intertemporal maximization problem:

where *r* > 0 is the intertemporal discount rate and ${\overline{Y}}_{d}=\overline{Y}(t-\tau )$.

By solving the maximization problem through the use of maximum principle and imposing, ex post, the equilibrium condition ${\overline{Y}}_{d}^{\alpha}={\overline{Y}}^{\alpha}$, the equilibrium dynamics are described by the following system of three nonadvanced differential equations with fixed delay, that is:

On the one hand, in absence of time delays (τ = 0), the model is saddlepath stable if private consumption and the environmental good are complements and cyclical behaviors appear only assuming substitutability, in line with [18]. On the other hand, with a positive time delay, in contrast with previous results, cyclical dynamics arise also in the case of complementarity between consumption of the private good and the environmental one.

## 6. Conclusions

In this work we have reviewed the literature highlighting the nonlinear interactions between economic activity and environmental dynamics. We have shown that, regardless of the time setting, the dynamic systems governing the evolution of both economic and environmental variables are able to exhibit the occurrence of cyclical dynamics, multistability or indeterminacy. In addition, we have analyzed the role of some decisive parameters in defining the dynamics of the model. This set of results points out that the analysis of the global properties of dynamic systems, i.e., the structure of both basins of attraction and Ω-limit sets, represents an essential instrument to deepen these issues.

In this field, the heterogeneity in the modeling of environmental resources has not been studied yet. These heterogeneities may concern (i) the introduction of environmental variables affecting the system both as a production input and a consumption good, or (ii) the inclusion of a multiplicity of environmental resources both in the utility and production functions. The effects of this kind of heterogeneities on the dynamic properties in models focusing on the coevolution of economic and environmental variables still represent open questions on which a scholar's agenda may concentrate upon.

Finally, a fruitful research question may concern the introduction of strategic interactions between agents. In particular, it could be relevant to analyse the dynamics driven by the recognition and the achievement of environmental targets, at the international level. Indeed, as described in Biancardi and Villani [28], the decision process related international environmental agreements may generate complex relationships (multistage games) which deserves further investigations.

## Author Contributions

AC, MS reviewed the literature on the coevolution of economic and environmental systems and performed the presentation (both theoretical and graphical) of all seminal contribution enclosed in the article.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

The authors gratefully acknowledge Angelo Antoci, Luca Gori, Federica Nieri, and Nicla Pasquini for their valuable suggestions and comments allowing to improve the quality of the work.

## Footnotes

2. ^See [3] for a complete review.

3. ^This work has been circulated as a working paper for many years and it was a sort of milestone for the literature of the 2000s on the relationship between growth and the environment.

4. ^This assumption is commonly used to simplify the analysis.

5. ^In economic literature, the Inada conditions are assumptions introduced in order to guarantee the stability of an interior equilibrium (see [1]).

6. ^In dynamic modeling, a balanced growth path is a trajectory such that all variables grow at constant (but potentially different) rates.

7. ^In order to simplify the notation, the subscripts denote partial derivatives. Therefore, the general equality ${V}_{x}(x,y)=\frac{\partial V(x,y)}{\partial x}={V}_{x}$ henceforth holds.

8. ^Savings are inelastically supplied to the firms.

10. ^Agents do not invest in defensive environmental expenditures.

11. ^About the role of government in public expenditure for the environmental maintenance, see [15].

12. ^With regard to local indeterminacy in OLG models where environmental quality is considered, see also [19]. In this work, the authors show the existence of a unique stationary state and a condition for which it is locally indeterminate. In addition, differently from Zhang [10], the work underlines that the possible emergence of endogenous fluctuations does not require high levels of pollution emission rate.

13. ^The increase in *U*_{t} due to an increase in *C*_{t+1} is positively related to the value of the index of environmental quality *E*_{t+1} if σ ∈ (0, 1), viceversa if σ > 1.

14. ^Small changes in *K*_{0} or *l*_{0} may generate a change in the Ω-limit set approached by the economy.

15. ^Indeed, a further condition is given by $\eta \ge \frac{\u03f5}{\u03f5+\alpha \beta}$.

## References

2. Antoci A, Galeotti M, Russu P. Poverty trap and global indeterminacy in a growth model with open-access natural resources. *J Econ Theory* (2011) **146**:569–91. doi: 10.1016/j.jet.2010.12.003

4. Brock W, Taylor MS. The green solow model. *J Econ Growth* (2010) **15**:127–53. doi: 10.1007/s10887-010-9051-0

5. Xepapadeas A. Economic growth and the environment. *Handb Environ Econ.* (2005) **3**:1219–71. doi: 10.1016/S1574-0099(05)03023-8

7. Li TY, Yorke JA. Period three implies chaos. *Am Math Monthly* (1975) **82**:985–92. doi: 10.1080/00029890.1975.11994008

8. Yorke JA, Yorke E. Chaotic behavior and fluid dynamics. In: Swinney HL, Gollub JR, editors. *Hydrodynamic Instabilities and the Transition to Turbulence*. (Berlin: Springer) (1981). p. 77–95.

9. John A, Pecchenino R. An overlapping generations model of growth and the environment. *Econ. J.* (1994) **104**:1393–410. doi: 10.2307/2235455

10. Zhang J. Environmental sustainability, nonlinear dynamics and chaos. *Econ Theory* (1999) **14**:489–500. doi: 10.1007/s001990050307

11. Naimzada A, Sodini M. Multiple attractors and nonlinear dynamics in an overlapping generations model with environment. *Discr Dyn Nat Soc.* (2010) **2010**:503695. doi: 10.1155/2010/503695

12. Azariadis C, Drazen A. Threshold externalities in economic development. *Q J Econ.* (1990) **105**:501–26. doi: 10.2307/2937797

13. Kozluk T, Zipperer V. Environmental policies and productivity growth. *OECD J.* (2015) **2014**:155–85. doi: 10.1787/eco_studies-2014-5jz2drqml75j

14. Caravaggio A, Sodini M. Multiple attractors and dynamics in an OLG model with productive environment. *Commun Nonlinear Sci Numer Simul*. (2017) **58**:167–84. doi: 10.1016/j.cnsns.2017.09.019

15. Raffin N, Seegmuller T. Longevity, pollution and growth. *Math Soc Sci.* (2014) **69**:22–33. doi: 10.1016/j.mathsocsci.2014.01.005

16. Smulders S, Gradus R. Pollution abatement and long-term growth. *Eur J Polit Econ.* (1996) **12**:505–32. doi: 10.1016/S0176-2680(96)00013-4

17. Antoci A, Sodini M. Indeterminacy, bifurcations and chaos in an overlapping generations model with negative environmental externalities. *Chaos Solitons Fract.* (2009) **42**:1439–50. doi: 10.1016/j.chaos.2009.03.055

18. Antoci A, Gori L, Sodini M. Nonlinear dynamics and global indeterminacy in an overlapping generations model with environmental resources. *Commun Nonlinear Sci Numer Simul.* (2016) **38**:59–71. doi: 10.1016/j.cnsns.2016.02.005

19. Seegmuller T, Verchère A. A note on indeterminacy in overlapping generations economies with environment and endogenous labor supply. *Macroecon Dyn.* (2007) **11**:423–9. doi: 10.1017/S1365100506060093

20. Coury T, Wen Y. Global indeterminacy in locally determinate real business cycle models. *Int J Econ Theory* (2009) **5**:49–60. doi: 10.1111/j.1742-7363.2008.00102.x

21. Antoci A. Environmental degradation as engine of undesirable economic growth via self-protection consumption choices. *Ecol Econ.* (2009) **68**:1385–97. doi: 10.1016/j.ecolecon.2008.09.009

22. Bennett RL, Farmer RE. Indeterminacy with non-separable utility. *J Econ Theory* (2000) **93**:118–43. doi: 10.1006/jeth.1999.2633

23. Hurwitz A. On the conditions under which an equation has only roots with negative real parts. *Math Ann.* (1895) **46**:273–84.

24. Matsumoto A, Szidarovszky F. Asymptotic behavior of a delay differential neoclassical growth model. *Sustainability* (2013) **5**:440–55. doi: 10.3390/su5020440

25. Ferrara M, Guerrini L, Sodini M. Stability and nonlinear dynamics in a Solow model with pollution. *Nonlinear Anal Model Control* (2014) **19**:565–77. doi: 10.15388/NA.2014.4.3

26. Ferrara M, Gori L, Sodini M. A continuous time economic growth model with time delays in environmental degradation. *J Inform Optim Sci.* (2018).

27. Wirl F. Stability and limit cycles in one-dimensional dynamic optimisations of competitive agents with a market externality. *J Evol Econ.* (1997) **7**:73–89. doi: 10.1007/s001910050035

Keywords: environmental quality, nonlinear dynamics, chaos, multistability, indeterminacy

Citation: Caravaggio A and Sodini M (2018) Nonlinear Dynamics in Coevolution of Economic and Environmental Systems. *Front. Appl. Math. Stat*. 4:26. doi: 10.3389/fams.2018.00026

Received: 23 December 2017; Accepted: 08 June 2018;

Published: 29 June 2018.

Edited by:

Gian Italo Bischi, Università degli Studi di Urbino Carlo Bo, ItalyReviewed by:

Villani Giovanni, Università degli Studi di Bari Aldo Moro, ItalyAnastasiia Panchuk, National Academy of Sciences of Ukraine, Ukraine

Copyright © 2018 Caravaggio and Sodini. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Andrea Caravaggio, a.caravaggio@studenti.unimc.it