Next Article in Journal
Controllability and Hyers–Ulam Stability of Differential Systems with Pure Delay
Next Article in Special Issue
Pursuit Differential Game with Slow Pursuers on the 1-Skeleton Graph of the Icosahedron
Previous Article in Journal
The Research on the Interactions between the Emerging and Developed Markets: From Region and Structural Break Perspectives
Previous Article in Special Issue
A Note on the Concept of Time in Extensive Games
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Qualitative Analysis of a Model of Renewable Resources and Population with Distributed Delays

by
Tiziana Ciano
1,2,
Massimiliano Ferrara
1,2,3,* and
Luca Guerrini
2,4
1
Department of Law, Economics and Human Sciences, University “Mediterranea” of Reggio Calabria, Via dell’Università, 25, 89124 Reggio Calabria, Italy
2
Decision Lab, University “Mediterranea” of Reggio Calabria, Via dell’Università, 25, 89124 Reggio Calabria, Italy
3
Department of Management and Technology, ICRIOS-The Invernizzi Centre for Research in Innovation, Organization, Strategy and Entrepreneurship Bocconi University, Via Sarfatti, 25, 20136 Milano, Italy
4
Department of Management, Polytechnic University of Marche, Piazzale Martelli 8, 60121 Ancona, Italy
*
Author to whom correspondence should be addressed.
Submission received: 18 March 2022 / Revised: 5 April 2022 / Accepted: 7 April 2022 / Published: 10 April 2022
(This article belongs to the Special Issue Differential Games and Its Applications)

Abstract

:
This work generalizes Matsumoto et al.’s dynamic model of population and renewable resources by substituting a distributed delay for the time delay. It is proved that the equilibrium point may lose or gain local stability, allowing for the observation of alternating stability/instability areas if some conditions hold.
MSC:
37H20; 37N40

1. Introduction

Following the findings of archeological studies, many Pacific islands follow similar evolutionary patterns in terms of natural resources and population dynamics. The research by Brander and Taylor [1] (BT henceforth) examines the archeological and anthropological findings from Easter Island as a case in point from an economic standpoint. They present a general equilibrium model of renewable resources and population dynamics to explain how and why Easter Island grew and fell throughout the 1400 years between the 4th century and the middle of the 18th century. Their findings suggest that an economic model linking resources and population dynamics may be able to explain not only the sources of past historical evolution discovered on these small islands, but also the possibility of sustainable growth for our global economy in an era in which a rapidly increasing population and a rapidly degrading environment are becoming serious problems. The analysis included inside the BT model has been expanded in several ways. Coats and Dalton [2] investigated the influence of market institutions and various property-rights frameworks on the distribution of wealth. According to Reuveny and Decker [3], the long-run dynamics of Easter Island were altered by technical advancement and population management reform in the past. However, there has not been much disclosed about the “history” of Easter Island. It is thought that a small group of Polynesians came on the island about the year 400, that deforestation began around the year 1000, that the majority of the sculptures were carved between the years 1000 and 1400, and so on. BT and other academics are attempting to replicate a dynamic pattern of natural resources and population based on this “common knowledge” in their study. However, “wisdom” is still just one of the numerous plausible possibilities, and it has not yet been proven beyond reasonable doubt. For instance, new reconsiderations of archeological material on the island, according to Intoh [4], indicate that the arrival time of the Polynesians on the island cannot be determined with certainty. It is merely a guess that a small number of people arrived between 410 and 1270 AD. Compared to traditional knowledge, this discovery is a jarring departure. In other words, we may have a different history and an evolutionary pattern of Easter Island than what is now known, although the historical data is consistent. It is thus essential to develop a model for small islands capable of generating a variety of dynamic patterns to cope with the confusing qualities of archeological data. Agricultural output may be critical to economic activity in a preindustrial economic civilization. It is well-known that a vital aspect of this kind of agricultural production is the substantial time lag between when producers decide to sow seeds in the fields and when they harvest the crops. Thus, it is logical to ask how delayed output has influenced the evolutionary pattern of a small island economy. Following this consideration, Matsumoto et al. [5] incorporate a delay in manufacturing into the BT model and discover the destabilizing effect caused by the delay in production on the evolution of a small island economy; namely, there is a critical value of the delay for which a loss of stability occurs. In the current literature, time delays are treated as fixed or continuously distributed (distributed delay henceforth). The former relates to economic situations when there is a specified time gap for the actors involved. This latter is suited for economic scenarios where actors’ delays vary. A fundamental issue is the unknown nature of time delays. Distributed delays, on the other hand, are the weighted average of all prior data from time zero to the present. Thus, distributed delays better describe time-delayed economic systems (see MacDonald [6]). Additionally, Caperon [7] shows there is some experimental evidence that they are more precise than those with instantaneous time delays. Cushing [8] pioneered and popularized distributed delays in mathematical biology, while Invernizzi and Medio [9] brought distributed delays to mathematical economics. In this paper, the authors reconsider Matsumoto et al.’s model [5] replacing time delay with distributed delay and investigate the dynamic effects of such lags on the adjustment process in population and stock of natural resources. Analytically, our system of functional delay differential equations has a single equilibrium point. The qualitative study demonstrates that this critical point loses or gains local stability when the lag duration increases or decreases. If certain circumstances are met, a series of intervals in which zones of stability and instability alternate may occur. This indicates that it may have significant challenges stabilizing the economic system. Another significant possibility that emerges from the qualitative investigation is the occurrence of Hopf bifurcation-induced limit cycles. The paper is organised as follows. Section 2 considers the distributed delay version of the time-delayed dynamic model of Matsumoto et al. [5]. Section 3 discusses the dynamics of weak delays in which the growth rate’s greatest weighted response is to present population density, whereas previous densities have an exponentially declining effect. Section 4 deals with strong delays in the sense that the largest effect on growth rate response at any time t is related to population density at the prior time tT. Section 5 and Section 6 analyze the cases of a weak delay together with a fixed delay given by the average length of the continuous delay. Finally, Section 7 concludes.

2. The Model

Matsumoto et al. [5] introduce a discrete-time delay into a dynamic economic model for a small island based on the Easter Island study by BT [1]. The model is an economy with two items and three economic agents (two producers and one consumer). The agricultural good is the harvest of renewable resources, while the manufactured good is something else. The natural resource stock S ( t ) and population L ( t ) are supplied. Producers set labor and supply needs to optimize profitability. A manufacturing producer offers the manufactured item created by labor alone. BT presume that per capita use of the resource good increases fertility and/or reduces death in accordance with Malthusian population dynamics. Including a delay τ 0 in production to account for the reality that agricultural items require time to mature from sowing to harvesting, they obtain the following two-dimensional system of delay differential equations:
S . ( t ) = r S ( t ) 1 S ( t ) K α β L ( t ) S ( t τ ) , L ˙ ( t ) = L ( t ) c + α β ϕ S ( t τ ) ,
where α is a positive constant indicating the harvesting efficiency; β ( 0 , 1 ) indicates the agricultural good’s preference, ϕ is a positive constant. K is the maximum feasible size of the resource stock, while r denotes the natural resource’s intrinsic growth rate, and both are positive constants. The difference between the underlying birth and death rates determines population change. The net rate, represented as c, is expected to be negative. When modeling the delay distribution across a large population, using a discrete delay might be considered a crude approximation at times. Of course, if the delay is continuously distributed by a continuous distribution function with a mean delay equal to the discrete delay and a positive variance to account for the delay difference across people, it may seem to be much more realistic. As a result, we suggest a generalizing model (1) with distributed delays. The dynamics of the economy are then governed by the following integro-differential equations system
S . ( t ) = r S ( t ) 1 S ( t ) K α β L ( t ) t g 1 ( t r ) S ( r ) d r , L ˙ ( t ) = L ( t ) c + α β ϕ t g 2 ( t r ) S ( r ) d r ,
where the delay kernel g is given by a gamma distribution function, i.e.,
g j ( u ) = l T l u l 1 e l T u ( l 1 ) ! ,
with ( j , l ) = ( 1 , m ) or ( j , l ) = ( 2 , n ) , m , n positive integers, and T 0 is a parameter associated with the mean time delay of the distribution. The weighting function’s form is determined by the parameter l. The two cases l = 0 and l = 1 in (3) are the weak delay kernel and the strong delay kernel, respectively. Notice that the distribution function approaches the Dirac distribution as T 0 , while it converges to a fixed delay T as the shape parameters m and n approach infinity. It is clear that the delayed system has the same equilibrium point as the basic system (1). The coordinates of an interior equilibrium are obtained by solving
r 1 S K = α β L , c + α β ϕ S = 0 ,
and so they are
S e = c α β ϕ , L e = r α β ϕ K + c α 2 β 2 ϕ K , with α β ϕ K + c > 0 .
To analyze the stability of system (2), we should first convert its coordinates to create a new system centered on the equilibrium point ( S e , L e ) , and then linearize the resulting system at the origin to obtain its characteristic equation. In order to render the eigenvalue analysis analytically tractable, we shall focus on some special cases and apply the linear chain trick technique (see MacDonald [6]), which allows an equation with gamma distributed kernels to be replaced by an equivalent system of ordinary differential equations.

3. Case m = n = 1

System (2) with weak kernels, i.e.,
g 1 ( t r ) = g 2 ( t r ) = 1 T e 1 T ( t r ) ,
takes the form
S . ( t ) = r S ( t ) 1 S ( t ) K α β L ( t ) t 1 T e 1 T ( t r ) S ( r ) d r , L ˙ ( t ) = L ( t ) c + α β ϕ t 1 T e 1 T ( t r ) S ( r ) d r .
Introducing the new variable
x ( t ) = t 1 T e 1 T ( t r ) S ( r ) d r ,
then, according to the law of solving the derivative for an integral with parameterized variables, system (4) can be rewritten as the following three-dimensional system of ODEs:
S . ( t ) = r S ( t ) 1 S ( t ) K α β L ( t ) x ( t ) , x . ( t ) = 1 T S ( t ) x ( t ) , L ˙ ( t ) = L ( t ) c + α β ϕ x ( t ) .
The equilibrium ( S e , L e ) of system (2) is transformed into the equilibrium ( S e , x e , L e ) of (5), where x e = S e . Thus, the stability study of equilibrium ( S e , L e ) of (2) is equivalent to the stability study of equilibrium ( S e , x e , L e ) of (5). Linearizing the system at the equilibrium point, one finds that the associated characteristic equation is
r 1 2 S e K λ α β L e α β x e 1 T 1 T λ 0 0 α β ϕ L e λ = 0 ,
i.e.,
λ 3 + a 1 λ 2 + a 2 λ + a 3 = 0 ,
where
a 1 = a 1 ( T ) = 1 T ( α β ϕ K + 2 c ) r α β ϕ K , a 2 = a 2 ( T ) = c r α β ϕ K T > 0
and
a 3 = a 3 ( T ) = ( α β ϕ K + c ) c r α β ϕ K T > 0 .
To study the stability of the equilibrium, we need to investigate the distribution of roots in the complex plane of the characteristic Equation (6). The stability will be determined by the real parts of the roots of Equation (6). If all roots of Equation (6) are located in the left-half complex plane, then the equilibrium is locally asymptotically stable. If Equation (6) has a root with positive real part, the equilibrium is unstable. According to the Routh-Hurwitz stability criterion, the necessary and sufficient conditions of asymptotic stability are a 1 > 0 , a 3 > 0 and a 1 a 2 > a 3 . From (7) and (8), these conditions reduce to a 1 a 2 > a 3 . Hence, we must have
1 T > r ( α β ϕ K + 2 c ) + α β ϕ K + c α β ϕ K α β ϕ K M .
Notice that the inequality (9) is verified for all T if M 0 and for T < 1 / M if M > 0 . Using the Hopf bifurcation theorem, we check the possibility of the emergence of a limit cycle at T = T * , where
T * = α β ϕ K r ( α β ϕ K + 2 c ) + α β ϕ K + c α β ϕ K .
At this value of T , one has a 1 * a 2 * = a 3 * , where a j * = a j ( T * ) ( j = 1 , 2 , 3 ) . As a result, Equation (6) becomes ( λ + a 1 * ) ( λ 2 + a 2 * ) = 0 , yielding the existence of a pair of purely imaginary roots λ 1 , 2 = ± i ω * , with ω * = a 2 * , and a real root λ 3 = a 1 * . Differentiating Equation (6) with respect to T , we obtain
( 3 λ 2 + 2 a 1 λ + a 2 ) d λ d T = a 1 λ 2 + a 2 λ + a 3 ,
where
a 1 = a 1 ( T ) = 1 T 2 , a 2 = a 2 ( T ) = c r α β ϕ K T 2 , a 3 = a 3 ( T ) = ( α β ϕ K + c ) c r α β ϕ K T 2 .
Recalling that ω * 2 = a 2 * , from (11) we derive
R e d λ d T λ = i ω * = a 1 * a 2 * + a 1 * a 2 * a 3 * 2 a 2 * + a 1 * 2 .
Then,
s i g n d R e λ d T λ = i ω * = s i g n a 1 * a 2 * a 1 * a 2 * + a 3 * = s i g n c r α β ϕ K T * 2 2 T * α β ϕ K + 2 c α β ϕ K r α β ϕ K + c > 0 .
Consequently, one pair of complex roots of (6) crosses through the imaginary axis transversally at T = T * from the left half-plane to the right half. Finally, notice from (11) that λ = i ω * is a simple root of (6). Otherwise, (11) would imply a 1 * i ω * 2 + a 2 * i ω * + a 3 * = 0 , and so the absurd ω * = 0 . By virtue of the previous analysis, we have the following conclusion:
Theorem 1.
Let T * be defined as in (10).
(1)
If α β ϕ K + 2 c < 0 and r α β ϕ K + c α β ϕ K / α β ϕ K + 2 c , the equilibrium point ( S e , x e , L e ) of (5) is locally asymptotically stable regardless of the values of T .
(2)
If α β ϕ K + 2 c 0 or α β ϕ K + 2 c < 0 and r < α β ϕ K + c α β ϕ K / α β ϕ K + 2 c , the equilibrium point ( S e , x e , L e ) of (5) is locally asymptotically stable for T < T * while it becomes locally unstable for T > T * and bifurcates to a limit cycle for T = T * .

4. Case m = n = 2

System (2) with strong kernels, i.e.,
g 1 ( t r ) = g 2 ( t r ) = 2 T 2 ( t r ) e 2 T ( t r ) ,
becomes
S . ( t ) = r S ( t ) 1 S ( t ) K α β L ( t ) t 2 T 2 ( t r ) e 2 T ( t r ) S ( r ) d r , L ˙ ( t ) = L ( t ) c + α β ϕ t 2 T 2 ( t r ) e 2 T ( t r ) S ( r ) d r .
Introducing the new variables
x ( t ) = t 2 T 2 ( t r ) e 2 T ( t r ) S ( r ) d r , y ( t ) = t 2 T 2 e 2 T ( t r ) S ( r ) d r ,
Equation (12) is rewritten as the following four-dimensional system of ODEs:
S . ( t ) = r S ( t ) 1 S ( t ) K α β L ( t ) x ( t ) , x . ( t ) = 2 T y ( t ) x ( t ) ) , y ˙ ( t ) = 2 T S ( t ) y ( t ) , L ˙ ( t ) = L ( t ) c + α β ϕ x ( t ) .
System (13) has a unique interior equilibrium point ( S e , x e , y e , L e ) , where x e = y e = S e . In this case, the characteristic equation associated to the linearized system is of the form
r 1 2 S e K λ α β L e 0 α β x e 0 2 T λ 2 T 0 2 T 0 2 T λ 0 0 α β ϕ L e 0 λ = 0 ,
i.e.,
λ 4 + a 1 λ 3 + a 2 λ 2 + a 3 λ + a 4 = 0 ,
where
a 1 = a 1 ( T ) = 4 T ( α β ϕ K + 2 c ) r α β ϕ K , a 2 = a 2 ( T ) = 4 T 1 T ( α β ϕ K + 2 c ) r α β ϕ K
and
a 3 = a 3 ( T ) = 4 c r α β ϕ K T 2 > 0 , a 4 = a 4 ( T ) = 4 ( α β ϕ K + c ) c r α β ϕ K T 2 > 0 .
For Equation (14), according to the Routh–Hurwitz criterion, the equilibrium is locally asymptotically stable if and only if a 1 > 0 , a 3 > 0 , a 4 > 0 and a 1 a 2 a 3 > a 3 2 + a 1 2 a 4 . Looking at the signs of the coefficients of (14), it follows these conditions reduce to a 1 > 0 and a 1 a 2 a 3 > a 3 2 + a 1 2 a 4 . The first of these two inequalities is equivalent to 4 α β ϕ K ( α β ϕ K + 2 c ) r T > 0 , which is verified for α β ϕ K + 2 c 0 or α β ϕ K + 2 c > 0 and T < 4 α β ϕ K / [ ( α β ϕ K + 2 c ) r ] . The second inequality instead means
ψ ( T ) = a 1 a 2 a 3 a 3 2 a 1 2 a 4 > 0 ,
where
ψ ( T ) = 4 α β ϕ K 2 ( α β ϕ K + 2 c ) 2 r 2 T 5 + 20 α β ϕ K 3 ( α β ϕ K + 2 c ) r T 4 + 4 α β ϕ K 2 4 α β ϕ K 2 + c r T 3 + ( α β ϕ K + c ) α β ϕ K + 2 c 2 r 2 T 2 + 8 α β ϕ K ( α β ϕ K + c ) ( α β ϕ K + 2 c ) r T 16 α β ϕ K 2 ( α β ϕ K + c ) .
Lemma 1.
(1)
If α β ϕ K + 2 c = 0 and 16 c + r < 0 , then ψ ( T ) > 0 for T > 4 / [ ( 16 c + r ) ] 3 = T 1 * , and ψ ( T 1 * ) = 0 .
(2)
If α β ϕ K + 2 c < 0 , then ψ ( T ) = 0 has exactly one positive root, say T 2 * .
(3)
If α β ϕ K + 2 c > 0 and T < 4 α β ϕ K / [ ( α β ϕ K + 2 c ) r ] , then ψ ( T ) = 0 may have one, three or five positive roots.
Proof. 
Let α β ϕ K + 2 c = 0 . Then ψ ( T ) = 16 c 3 [ 16 c + r T 3 + 4 ] . The conclusion is immediate. Let α β ϕ K + 2 c 0 . Then, one has the cases α β ϕ K + 2 c < 0 or α β ϕ K + 2 c > 0 and T < 4 α β ϕ K / [ ( α β ϕ K + 2 c ) r ] . Since ψ ( + ) = + and ψ ( 0 ) < 0 , there exists at least T = T * > 0 such that ψ ( T * ) = 0 . In the latter case, we must also have T * < 4 α β ϕ K / [ ( α β ϕ K + 2 c ) r ] . One can actually count the roots of ψ ( T ) = 0 using Descartes rule of signs, that says the number of positive roots of the polynomial is either equal to the number of sign differences between consecutive nonzero coefficients, or is less than it by an even number. □
Henceforth, let us assume α β ϕ K + 2 c = 0 or ψ ( T ) > 0 together with α β ϕ K + 2 c < 0 or α β ϕ K + 2 c > 0 and T < 4 α β ϕ K / [ ( α β ϕ K + 2 c ) r ] . Since there exists T = T * > 0 such that ψ ( T * ) = 0 , i.e., a 1 * a 2 * a 3 * a 3 * 2 = a 1 * 2 a 4 * , the characteristic Equation (14) can be factored as
λ 2 + a 3 * a 1 * λ 2 + a 1 * λ + a 1 * a 2 * a 3 * a 3 * 2 a 1 * = 0 .
We have three characteristic roots, two purely imaginary
λ 1 , 2 = ± i a 3 * a 1 *
and two roots,
λ 3 , 4 = a 1 * ± a 1 * 2 4 a 1 * a 2 * a 3 * a 3 * 2 a 1 * 2 ,
which have real parts different from zero since λ 3 + λ 4 = a 1 * < 0 and λ 3 λ 4 = a 1 * a 2 * a 3 * a 3 * 2 / a 1 * > 0 . Choosing T as a bifurcation parameter, we apply the Hopf bifurcation theorem to establish the existence of a cyclical movement. According to this theorem, one can establish the existence of a cyclic solution at T = T * . It remains to check that the root λ = i ω * , with ω * = ω ( T * ) > 0 , is simple as well as the transversality condition holds. A differentiation of (14) with respect to T yields
4 λ 3 + 3 a 1 λ 2 + 2 a 2 λ + a 3 d λ d T = a 1 λ 3 + a 2 λ 2 + a 3 λ + a 4 ,
where
a 1 = a 1 ( T ) = 4 T 2 < 0 , a 2 = a 2 ( T ) = 4 T 2 2 T ( α β ϕ K + 2 c ) r α β ϕ K < 0
and
a 3 = a 3 ( T ) = 8 c r α β ϕ K T 3 < 0 , a 4 = a 4 ( T ) = 8 ( α β ϕ K + c ) c r α β ϕ K T 3 < 0 .
If λ = i ω * is a repeated root of (14), then we get from (15) that a 1 λ 3 + a 2 λ 2 + a 3 λ + a 4 evaluated at T = T * must be zero. Thus, separating the real and imaginary parts, we find a 1 * ω * 2 = a 3 * ,   a 2 * ω * 2 = a 4 * , where a j * = a j ( T * )   ( j = 1 , 2 , 3 , 4 ) . However, λ = i ω * is a root of (14), and so a 1 * ω * 2 = a 3 * . In other words, one obtains a 3 * a 1 * = a 1 * a 3 * , and so the contradiction 0 = ( 2 / T * ) [ ( α β ϕ K + 2 c ) r / ( α β ϕ K ) ] = a 2 * T * 2 / 4 > 0 . Finally, from (15) we derive
d λ d T = a 1 λ 3 + a 2 λ 2 + a 3 λ + a 4 4 λ 3 + 3 a 1 λ 2 + 2 a 2 λ + a 3 ,
leading to
R e d λ d T λ = i ω * = a 1 * ψ ( T * ) 2 a 1 * 3 a 3 * + ω * 2 a 1 * a 2 * 2 a 3 * 2 ,
where
ψ ( T * ) = a 1 * a 2 * a 3 * + a 1 * a 2 * a 3 * + a 1 * a 2 * a 3 * 2 a 3 * a 3 * 2 a 1 * a 1 * a 4 * a 1 * 2 a 4 * .
A negative sign of ψ ( T * ) , and so a positive sign of (16), implies that one pair of complex roots of (14) crosses through the imaginary axis transversally at T = T * from the left half plane to the right half plane. On the other hand, a positive sign of ψ ( T * ) yields that the roots can cross the imaginary axis only from right to left as T increases. Summarizing the above analysis, we have the following results.
Theorem 2.
Let T 1 * and T 2 * be defined as in the previous Lemma.
(1)
Let α β ϕ K + 2 c = 0 and 16 c + r < 0 . The equilibrium point ( S e , x e , y e , L e ) of (13) is locally asymptotically stable for T < T 1 * , unstable for T > T 1 * , and it bifurcates to chaos at T = T 1 * .
(2)
Let α β ϕ K + 2 c < 0 . The equilibrium pointof ( S e , x e , y e , L e ) of (13) is locally asymptotically stable if ψ ( T ) > 0 and T < T 2 * . As T > T 2 * , it remains stable ψ ( T 2 * ) < 0 and becomes unstable if ψ ( T 2 * ) < 0 . In this latter case, a Hopf bifurcation appears at T = T 2 * .
(3)
Let α β ϕ K + 2 c > 0 and T < 4 α β ϕ K / [ ( α β ϕ K + 2 c ) r ] . The equilibrium pointof ( S e , x e , y e , L e ) of (13) is locally asymptotically stable if ψ ( T ) > 0 . and T < T * , where T = T * is the smallest value such that ψ ( T * ) = 0 . and ψ ( T * ) 0 , then a Hopf bifurcation occurs at the equilibrium point as T passes through T * . Moreover, stability switches may take place at values of T where ψ ( T ) = 0 .

5. Case m = 1 and n

System (2) becomes
S . ( t ) = r S ( t ) 1 S ( t ) K α β L ( t ) t 1 T e 1 T ( t r ) S ( r ) d r , L ˙ ( t ) = L ( t ) c + α β ϕ S ( t T ) .
Letting
x ( t ) = t 1 T e 1 T ( t r ) S ( r ) d r ,
the dynamical system (17) assumes the form
S . ( t ) = r S ( t ) 1 S ( t ) K α β L ( t ) x ( t ) , x . ( t ) = 1 T S ( t ) x ( t ) , L ˙ ( t ) = L ( t ) c + α β ϕ S ( t T ) .
The characteristic equation of the linearized system around the interior equilibrium point ( S e , x e , L e ) , where x e = S e is
r 1 2 S e K λ α β L e α β x e 1 T 1 T λ 0 α β ϕ L e e λ T 0 λ = 0 ,
i.e.,
λ 3 + A λ 2 + B λ + C + D λ e λ T = 0 ,
where
A = A ( T ) = 1 T α β ϕ K + 2 c α β ϕ K r , B = B ( T ) = c r α β ϕ K T > 0
and
C = C ( T ) = α β ϕ K + c c r α β ϕ K T > 0 , D = α β ϕ K + c c r α β ϕ K > 0 .
We now turn our attention to the issue of stability switching. The stability analysis concerns with whether all roots are located in the left half of the complex plane. When the delay changes, the study of stability switching is concerned with whether the roots cross the imaginary axis. Obviously, λ = i ω ( ω > 0 ) is a root of Equation (19) if and only if ω satisfies
i ω 3 A ω 2 + B i ω + C + D i ω cos ω T i sin ω T e i ω T = 0 .
Since e i ω T = cos ω T i sin ω T , separating the real and imaginary parts, we become
A ω 2 = D ω sin ω T + C cos ω T ,
ω 3 B ω = D ω cos ω T C sin ω T .
Squaring and adding these two equations yields
ω 6 + p ω 4 + q ω 2 + s = 0 ,
where
p = A 2 2 B , q = B 2 D 2 , s = C 2 < 0 .
Let z = ω 2 , then Equation (22) rewrites as follows:
f ( z ) = z 3 + p z 2 + q z + s = 0 .
Since f ( 0 ) = s < 0 and f ( + ) = + , we conclude that Equation (23) has at least one positive real root.
Lemma 2.
Let Δ = p 2 q 2 + 18 p q s 4 q 3 4 p 3 s 27 s 2 be the discriminant of the cubic equation. The following cases can be discerned.
(1)
If p 0 or q 0 , Equation (23) has only one real positive root.
(2)
If p < 0 , q > 0 and Δ > 0 , Equation (23) has only one real positive root.
(3)
If p < 0 , q > 0 and Δ 0 , Equation (23) has three real positive roots that are distinct if Δ 0 .
Proof. 
Since s < 0 , in virtue of Descartes rule, the three enunciated cases follow immediately. □
Let us denote by z 0 a solution of (23). Then, Equation (22) has a positive root ω 0 = z 0 . To determine the values of T at which Equation (19) has purely imaginary roots λ = ± i ω 0 , we derive from (20) and (21)
sin ω 0 T = D ω 0 2 + A C B D ω 0 3 D 2 ω 0 2 + C 2 , cos ω 0 T = A D + C ω 0 3 B C D 2 ω 0 2 + C 2 .
Therefore, from (24) we obtain
T 0 ( j ) = 1 ω 0 cos 1 D ω 0 2 + A C B D D 2 ω 0 2 + C 2 + 2 j π ω 0 , if A D + C ω 0 3 B C 0 , 2 ( j + 1 ) π ω 0 1 ω 0 cos 1 D ω 0 2 + A C B D D 2 ω 0 2 + C 2 , otherwise ,
where j = 0 , 1 , 2 , According to the Hopf bifurcation theorem, to verify stability switching, we need to determine the sign of the derivative of R e ( λ ) at T = T 0 ( j ) .
Proposition 1.
Let λ T be the root of (19) such that R e ( T 0 ( j ) ) = 0 , I m ( T 0 ( j ) ) = ω 0 . Then λ = ± i ω 0 are simple roots of (19) at T = T 0 ( j ) and
s i g n d R e λ d T T = T 0 ( j ) = s i g n f ( z 0 ) .
Proof. 
Substituting λ = λ T into (19) and by differentiating with respect to T yields
3 λ 2 + 2 A λ + B + D e λ T C + D λ T e λ T d λ d T = C + D λ λ e λ T .
For convenience, we check the sign of d λ / d T 1 , that is written as
d λ d T 1 = 3 λ 2 + 2 A λ + B λ C + D λ e λ T + D λ C + D λ T λ ,
Using (19), we obtain
d λ d T 1 = 3 λ 2 + 2 A λ + B λ λ 3 + A λ 2 + B λ + D λ C + D λ T λ .
Plugging λ = i ω 0 and T = T 0 ( j ) into (27), together with (23), we have
s i g n d R e λ d T T = T 0 ( j ) = s i g n R e d λ d T T = T 0 ( j ) 1 = s i g n 3 ω 0 4 + 2 A 2 2 B ω 0 2 + B 2 D 2 ω 0 6 + A 2 2 B ω 0 4 + B 2 ω 0 2 = s i g n f ( z 0 ) .
It remains to prove the simplicity of the root λ = i ω 0 . Assuming that this root is a repeated root of (19), then (26) implies ( C + D i ω 0 ) ω 0 e i ω 0 T 0 ( j ) = 0 , i.e., D ω 0 sin ω 0 T 0 ( j ) + C cos ω 0 T 0 ( j ) = 0 and D ω 0 cos ω 0 T 0 ( j ) C sin ω 0 T 0 ( j ) = 0 . Hence, from (20) and (21), we derive that A = 0 and B = ω 0 2 . Consequently, (22) yields the absurd B D 2 + C 2 = 0 . This concludes the proof. □
As T grows, all roots cross the imaginary axis from left to right if the sign is positive, indicating that stability is lost or instability is maintained. On the other hand, the negative sign indicates that the axis is being crossed in the opposite direction, implying that stability may be achieved. Bearing all the above discussion in mind, the following theorem holds.
Theorem 3.
Let f ( z ) and T 0 ( j ) ( j = 0 , 1 , 2 , ) be defined as in (23) and (25), respectively. If f ( z ) = 0 has one or two solutions, there is only one stability switch. If f ( z ) = 0 has three solutions, then at least a stability switch exists. System (18) may undergo a Hopf bifurcation at ( S e , x e , L e ) for those values of T at which a stability switch occurs.
Proof. 
Suppose that Equation (23) has a unique solution and let ω 0 be the corresponding unique positive root of Equation (22). Recalling f ( 0 ) < 0 and f ( + ) = + , the polynomial f ( z ) is an increasing function in a neighborhood of ω 0 , and so its the derivative at ω 0 is positive. Consequently, λ = ± i ω 0 crosses the imaginary axis from left to right. Hence, a stability switch occurs. Suppose now that Equation (23) has two solutions, with ω 1 < ω 2 the resulting positive roots of (22). Then, the polynomial f ( z ) is a decreasing function in a neighborhood of ω 1 and an increasing function in a neighborhood of ω 2 . Thus, one has a crossing of imaginary axis from left to right in correspondence of ω 2 and a crossing from right to left in correspondence of ω 1 . As a result, there is just one stability switch. Finally, suppose that Equation (23) has three solutions and ω 1 < ω 2 < ω 3 are the three positive roots of (22). In this case, the polynomial f ( z ) is an increasing function in a neighborhood of ω 1 , ω 3 , and a decreasing function in a neighborhood of ω 2 . Crossing of imaginary axis is from left to right in correspondence of ω 1 and ω 3 , from right to left in correspondence of ω 2 . There is at least one stability switch present. □

6. Case m and n = 1

System (2) takes the form
S . ( t ) = r S ( t ) 1 S ( t ) K α β L ( t ) S ( t T ) , L ˙ ( t ) = L ( t ) c + α β ϕ t 1 T e 1 T ( t r ) S ( r ) d r .
Setting
x ( t ) = t 1 T e 1 T ( t r ) S ( r ) d r ,
system (28) rewrites as
S . = r S ( t ) 1 S ( t ) K α β L ( t ) S ( t T ) , x . = 1 T S ( t ) x ( t ) , L ˙ = L ( t ) c + α β ϕ x ( t ) .
To examine local dynamics of the above system in a neighborhood of the equilibrium point ( S e , x e , L e ) , where x e = S e , we consider the linearized system and derive the following characteristic equation
r 1 2 S e K λ α β L e e λ T 0 α β S e 1 T 1 T λ 0 0 α β ϕ L e λ = 0 ,
i.e.,
λ 3 + M λ 2 + N λ + P + Q λ + R λ 2 e λ T = 0 ,
where
M = M ( T ) = 1 T α β ϕ K + 2 c α β ϕ K r , N = N ( T ) = α β ϕ K + 2 c r α β ϕ K T
and
P = P ( T ) = α β ϕ K + c c r α β ϕ K T > 0 , Q = Q ( T ) = α β ϕ K + c r α β ϕ K T > 0 , R = R ( T ) = α β ϕ K + c r α β ϕ K > 0 .
Let λ = i ω with ω > 0 be a solution of (30). Substituting it into (30), and separating the real and imaginary parts, we have
M ω 2 P = Q ω sin ω T + R ω 2 cos ω T ,
ω 3 N ω = Q ω cos ω T + R ω 2 sin ω T .
The sum of the squares of these two equations yields the sextic equation in ω
ω 6 + p ω 4 + q ω 2 + s = 0 ,
where
p = M 2 R 2 2 N , q = N 2 Q 2 2 M P , s = P 2 > 0 .
Letting z = ω 2 , Equation (33) reduces to a cubic equation
g ( z ) = z 3 + p z 2 + q z + s = 0 .
Lemma 3.
Let Δ = p 2 q 2 + 18 p q s 4 q 3 4 p 3 s 27 s 2 denote the discriminant of the cubic equation. The following scenarios may be identified.
(1)
If p 0 and q 0 , Equation (34) has no real positive roots.
(2)
If p < 0 or q < 0 and Δ > 0 , Equation (34) has no real positive roots.
(3)
If p < 0 , q < 0 and Δ 0 , Equation (34) has two real positive roots, which are distinct if Δ 0 .
Proof. 
The statement follows from Descartes rule, noticing that s > 0 .
From (31) and (32), we gain the equation
sin ω T = R ω 4 + M Q + N R ω 3 P Q R 2 ω 3 + Q 2 ω , cos ω T = Q M R ω 2 + P R N Q R 2 ω 2 + Q 2 ,
Let ω 0 be the positive root of (33), corresponding to the positive solution z 0 of (34). Then, we conclude that Equation (30) has purely imaginary roots λ = ± i ω 0 for
T 0 ( j ) = 1 ω 0 cos 1 Q M R ω 0 2 + P R N Q R 2 ω 0 2 + Q 2 + 2 j π ω 0 , if R ω 0 4 + M Q + N R ω 0 3 P Q 0 , 2 ( j + 1 ) π ω 0 1 ω 0 cos 1 Q M R ω 0 2 + P R N Q R 2 ω 0 2 + Q 2 , otherwise ,
where j = 0 , 1 , 2 , . . . To verify whether a Hopf bifurcation occurs, we study the transversality condition as follows.
Proposition 2.
The pair of pure imaginary roots λ = ± i ω 0 at T 0 ( j ) , j = 1 , 2 , , crosses the imaginary axis from left to right if g ( z 0 ) > 0 , and crosses the imaginary axis from right to left if g ( z 0 ) < 0 .
Proof. 
Differentiating Equation (30) for T. we obtain
3 λ 2 + 2 M λ + N + Q + 2 R λ e λ T Q λ + R λ 2 T e λ T d λ d T = Q λ + R λ 2 λ e λ T ,
Therefore,
d λ d T 1 = 3 λ 2 + 2 M λ + N λ λ 3 + M λ 2 + N λ + P + Q + 2 R λ λ Q λ + R λ 2 T λ .
by using (30). Thus,
s i g n d R e λ d T T = T 0 ( j ) = s i g n R e d λ d T T = T 0 ( j ) 1 = s i g n 3 ω k 4 + 2 M 2 R 2 2 N ω k 2 + N 2 Q 2 2 M P ω k 6 + M 2 2 N ω k 4 + N 2 2 M P ω k 2 + P 2 = s i g n g ( z 0 ) .
The conclusion holds. □
Then, we have the following results.
Theorem 4.
Let g ( z ) and T 0 ( j ) ( j = 0 , 1 , 2 , . . . ) be defined as in (34) and (35) respectively. If g ( z ) = 0 has no solutions, then no stability switches exist. If g ( z ) = 0 has one or two solutions, then there is at least a stability switch. For those values of T at which a stability transition occurs, system (29) may exhibit a Hopf bifurcation.
Proof. 
If Equation (34) has no solutions, the proof is trivial. The second part of the statement follows the same justifications for the theorem presented in the previous section. □

7. Conclusions

This paper generalizes the BT’s delayed continuous-time model for small islands presented by Matsumoto et al. [5]. Time lags modeled by way of distributed delays allow the reduction of the dynamics to a set of ordinary differential equations. The gamma distribution function has been considered to have only weak and strong kernels to simplify our analysis. Sufficient conditions of existence and stability of equilibria in four different cases are derived. In contrast with the delayed model, it is found that the positive equilibrium becomes unstable for all large delay values, or the stability of equilibrium switches back, leading to the occurrence of multiple stability switches.

Author Contributions

Conceptualization, M.F. and L.G.; methodology, T.C.; software, T.C.; validation, M.F., L.G. and T.C.; formal analysis, T.C.; investigation, L.G.; resources, M.F.; data curation, T.C.; writing—original draft preparation, L.G.; writing—review and editing, M.F.; visualization, L.G.; supervision, M.F.; project administration, M.F.; funding acquisition, M.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the anonymous reviewers for their careful reading of our manuscript and their many insightful comments and suggestions. The usual disclaimers apply.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brander, J.A.; Taylor, N.S. The simple economics of Easter Island: A Ricardo–Malthus model of renewable resource use. Am. Econ. Rev. 1998, 88, 119–138. [Google Scholar]
  2. Coats, R.M.; Dalton, T.R. Could institutional reform have saved Easter Island? J. Evol. Econ. 2000, 10, 489–505. [Google Scholar]
  3. Reuveny, R.; Decker, C.S. Easter Island: Historical anecdote or warning for the future. Ecol. Econ. 2000, 35, 271–287. [Google Scholar] [CrossRef]
  4. Intoh, M. Prehistoric oceania. In Oceania History; Yamamoto, M., Ed.; Yamakawa Publishing Co.: Tokyo, Japan, 2000; pp. 17–45. [Google Scholar]
  5. Matsumoto, A.; Suzuki, M.; Saito, Y. Time-delayed dynamic model of renewable resource and population. In Time and Space in Economics; Asada, T., Ishikawa, T., Eds.; Springer: Tokyo, Japan, 2007; pp. 145–159. [Google Scholar]
  6. MacDonald, N. Time Lags in Biological Systems; Springer: New York, NY, USA, 1978. [Google Scholar]
  7. Caperon, J. Time lag in population growth response of Isochrysis galbana to a variable nitrate environment. Ecology 1969, 50, 188–192. [Google Scholar] [CrossRef]
  8. Cushing, J.M. Integro-Differential Equations and Delay Models in Population Dynamics; Springer: Berlin, Germany, 1977. [Google Scholar]
  9. Invernizzi, S.; Medio, A. On lags and chaos in economic dynamic models. J. Math. Econ. 1991, 20, 521–550. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ciano, T.; Ferrara, M.; Guerrini, L. Qualitative Analysis of a Model of Renewable Resources and Population with Distributed Delays. Mathematics 2022, 10, 1247. https://0-doi-org.brum.beds.ac.uk/10.3390/math10081247

AMA Style

Ciano T, Ferrara M, Guerrini L. Qualitative Analysis of a Model of Renewable Resources and Population with Distributed Delays. Mathematics. 2022; 10(8):1247. https://0-doi-org.brum.beds.ac.uk/10.3390/math10081247

Chicago/Turabian Style

Ciano, Tiziana, Massimiliano Ferrara, and Luca Guerrini. 2022. "Qualitative Analysis of a Model of Renewable Resources and Population with Distributed Delays" Mathematics 10, no. 8: 1247. https://0-doi-org.brum.beds.ac.uk/10.3390/math10081247

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop