Discrete Dynamics in Nature and Society
Volume 2020 (2020), Article ID 1053283, 10 pages
https://doi.org/10.1155/2020/1053283

Research Article

Complex Behavior Analysis of a Fractional-Order Land Dynamical Model with Holling-II Type Land Reclamation Rate on Time Delay

Li Wu,1,2 Zhouhong Li,³ Yuan Zhang,³ and Binggeng Xie¹

¹College of Resources and Environmental Sciences, Hunan Normal University, Changsha, Hunan 410028, China
²Department of Geography and Land Engineering, Yuxi Normal University, Yuxi, Yunnan 653100, China
³Department of Mathematics, Yuxi Normal University, Yuxi, Yunnan 653100, China

Correspondence should be addressed to Binggeng Xie; xbgyb1961@163.com

Received 14 May 2020; Accepted 25 June 2020; Published 25 July 2020

Guest Editor: Jinliang Wang

Copyright © 2020 Li Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In this paper, a fractional-order land model with Holling-II type transformation rate and time delay is investigated. First of all, the variable-order fractional derivative is defined in the Caputo type. Second, by applying time delay as the bifurcation parameter, some criteria to determine the stability and Hopf bifurcation of the model are presented. It turns out that the time delay can drive the model to be oscillatory, even when its steady state is stable. Finally, one numerical example is proposed to justify the validity of theoretical analysis. These results may provide insights to the development of a reasonable strategy to control land-use change.

1. Introduction

Since the 20th century, the global environment has changed at an unprecedented speed, and a series of major global environmental problems pose a serious threat to the survival and development of mankind. From the perspective of earth system science, the generation of global environmental problems can be recognized as the result of the interaction between the earth’s atmosphere, hydrosphere, biosphere, lithosphere, and human activities [1, 2]. Land use/land cover change (LUCC) is as a key link connecting the four circles of the earth system, being an ideal entry point for the study of natural and human processes, and has rapidly become the focus of global change research [3, 4]. LUCC dynamic change is to describe, evaluate, interpret, and predict the quantity, quality, spatial distribution, classification, rate, and process of land use/land cover change by accumulating spatiotemporal continuous land use/land cover data and using the mathematical model [5]. The driving forces of land use change are mainly divided into natural and social systems. Notably, hydrology, topography, and geomorphology [6, 7] are the main driving forces in the natural system, while economic development, population growth, and policies are the main driving forces in the social and economic systems [8]. Generally speaking, it is the interwoven factors in these two systems that cause land use changes [9], but compared with the natural factors (such as climate, soil, and topography), human factors (such as regional policy, economic development, and population growth) play a dominant role [10].

Since the key of LUCC research is to understand the driving force and driving mechanism [11], a variety of system analysis and mathematical statistic methods have been well applied in the research of driving force. For example, the internal and human driving force of land use change is studied through the analysis of the multiple linear regression model [12]. The system dynamic model was used to quantitatively diagnose the contribution of each driving factor to land use change [13]. The CLUE comprehensive model was used to establish the interaction between land-use change and its influencing factors, and land change and its spatial distribution were simulated [14–16]. Therefore, the comprehensive mathematical model will be a development trend of LUCC research at present and in the future.

In 1997, Dobson et al. considered the land model with four coupled differential equations [17]:

where stands for an original area of pristine forest habitat, denotes the agricultural land, unused land is , and is the population density with time , for more detail, one can read [17] and the references cited therein. In [18], Chen et al. studied the dynamical analysis of the land model with Holling-II type land reclamation rate, and they found the positive equilibrium point, the basic reproductive rate , and, furthermore, some sufficient conditions for the global stability of the positive equilibrium and one of the boundary equilibrium.

On the one hand, it is well known that time delays are unavoidable for population modelling. It is therefore important to consider the dynamics for the population models with time delays, and dynamic complex analysis is obviously one of the most important problems. For example, delay differential equations have been applied in the stability of population dynamic systems [19–25], impulsive effects, and control systems [26–28].

On the other hand, there are two types of differential equations, namely, integer-order and fractional-order differential equations. Traditionally, the fractional-order differential equations have enjoyed a preference over integer-order differential equations because of the mathematical tractability of fractional systems. Over the past 30 decades, fractional-order differential equations have been applied in many fields, such as design and control of various ecological systems [29–32], secure communication [33–35], and system control [36–42].

Recently, the existence of Hopf bifurcation of differential equations has been studied as an important qualitative behavior of integer-order differential equations [43]. Moreover, fractional calculus is merged into complicated, dynamical systems which extremely renovate the theory of the design and control performance for complex systems. It has been discovered that physical phenomena in nature can be depicted more accurately by fractional-order models in comparison with classical integer-order ones [44, 45]. Some scholars introduced fractional calculus into population models and constructed fractional population or epidemic models [46–50], fractional dynamic systems, and fractional neural networks [51–53]. However, to the best of our knowledge, there are few studies to investigate the existence of Hopf bifurcation to the fractional-order land model with time delay.

Motived by the above ideas, in this paper, we will consider the following fractional-order land model with Holling-II type response and time delay:

where is the fractional order, represent the land of wood and grass, survival land, and unused land, respectively, and is the density of the predator population at time ; (land area of survival) of a period of time becomes unused land (area ), which in turn recovers through natural succession or ecological restoration to become a forest after a time interval . The unused land may also be restored to survival viable land after a time interval . Average clearing ability is described by constant . is the growth rate into the population. is the time delay required for the gestation of the mature population.

The main contributions can be summed up in three aspects:

(1)A novel fractional-order land model with Holling-II type land reclamation rate and time delay is formulated(2)Two primary dynamical properties—stability and oscillation—of the delayed fractional-order land model are investigated(3)The influences of the order on the Hopf bifurcation are obtained

Throughout this paper, we address the following assumptions.

Assumption 1. .

Assumption 2. .
Suppose Assumptions 1 and 2 hold, then system (2) has an unique positive equilibrium point , which is described by

The remainder of the current paper is organized as follows: In Section 2, some definitions and lemmas of fractional calculus are recalled. In Section 3, the fractional land model with Holling-II type land reclamation rate and time delay is investigated by using time delay as the bifurcation parameter, and the conditions of Hopf bifurcation are presented. One numerical example is given to illustrate the effectiveness of our main results in Section 4. Finally, conclusions are drawn in the last section.

2. Basic Tools of Fractional Calculus

In this section, we recall some definitions and lemmas of fractional calculus, which can be used in the proofs of the main results in Section 3.

There are several definitions of fractional derivatives. The Grnwald–Letnikov definition, the Riemann–Liouville definition, and the Caputo definition are usually used to deal with fractional-order systems. Since the Caputo derivative only requires the initial conditions which are based on the integer-order derivative and represents well-understood features of physical situation, it is more applicable to real-world problems. Hence, the Caputo fractional-order derivative is adopted in this paper.

Definition 1. (see [44]). The fractional integral of order for a function is defined as

where , , is the gamma function, and .

Definition 2. (see [44]). The Caputo fractional derivative of order for a function is defined by

where and is a positive integer such that .
Moreover, when ,

Lemma 1. (see [41]). The following autonomous system

where , , and , is asymptotically stable if and only if . In this case, each component of the states decays towards 0 like . Also, this system is stable if and only if , and those critical eigenvalues that satisfy have geometric multiplicity one.

3. Main Results

In this section, by applying the previous analytic technique, we shall investigate the stability and bifurcation of system (2) with time delay by taking time delay as the bifurcation parameter. The conditions of delay-induced bifurcation will be derived.

From reference [18], the dynamics of system (2) crucially depend on the basic reproduction number , which is given by

Let , and then we have

Taking advantage of the Taylor expansion formula, the linearized system of system (9) at the zero equilibrium is

where .

The associated characteristic equation of system (10) can be obtained as

which leads to

where

Let and , and from equation (12), we have

Assume that is a root of equation (12), . Substituting into equation (12) and separating the real and imaginary parts, one can have

Applying equation (15), direct calculation yields

where

From (16), we obtain

Define

where is defined by equation (18).

To derive the condition of the occurrence for Hopf bifurcation, we further give the following assumption.

Assumption 3.

where

Lemma 2. Let be the root of equation (18) near satisfying , , then the following transversality condition holds

Proof. By using implicit function theorem and differentiating (12) with respect to , we have

where are the derivatives of .
Thus, we obtain

where

Let and be the real and imaginary parts of , respectively. We further suppose that and are the real and imaginary parts of , respectively, and , are the real and imaginary parts of , respectively.
Hence,

Therefore, based on Assumption 3, the transversality condition is satisfied. This completes the proof of Lemma 2.
Denote

where

Lemma 3. If Assumptions 1 and 2 hold, then the positive equilibrium point of the delayed fractional-order land model (2) is asymptotically stable when .

Proof. When , system (2) changes into

If the conditions of hold and, furthermore, if Assumptions 1 and 2 are satisfied, it is easy to check from the Routh–Hurwitz criterion that the four eigenvalues of the characteristic equation (29) have negative real parts. Hence, the positive equilibrium point of the delayed fractional-order land model (2) is asymptotically stable when .
From the above analysis, together with Lemma 3, we will give the following theorem.

Theorem 1. For the model (2):

(i)Under Assumptions 1–3, the equilibrium point is global asymptotically stable for (ii)Under Assumptions 1–3:

(a)The equilibrium point is locally asymptotically stable for (b)Model (2) undergoes a Hopf bifurcation at the when , i.e., it has a branch of periodic solutions bifurcating from the equilibrium near

4. Numerical Example and Application

In this section, we give one example to demonstrate the effectiveness of the proposed approach. The simulation results are based on theAdams–Bashforth–Moulton predictor-corrector scheme [54], and step length .

Consider the following fractional-order land dynamical model with time delay:

in this case, taking , and according to Yuxi city land and population data of Yunan Province, China, let . By computing, we obtain the positive equilibrium point as . Moreover, we have , , and . Furthermore, the initial values are selected as , and it is easy to verify that the conditions of (ii) in Theorem 1 are satisfied. Figures 1–8 depict the equilibrium point which is locally asymptotically stable when , while Figures 9–16 display the equilibrium point that is unstable, and Hopf bifurcation occurs when (see Supplementary Materials (available here)).

Figure 1: Time series solution of system (30) with , .

Figure 2: Time series solution of system (30) with , .

Figure 3: Time series solution of system (30) with , .

Figure 4: Time series solution of system (30) with , .

Figure 5: Portrait diagram of land model (30) with , .

Figure 6: Portrait diagram of land model (30) with , .

Figure 7: Portrait diagram of land model (30) with , .

Figure 8: Portrait diagram of land model (30) with , .

Figure 9: Time series solution of system (30) with , .

Figure 10: Time series solution of system (30) with , .

Figure 11: Time series solution of system (30) with , .

Figure 12: Time series solution of system (30) with , .

Figure 13: Portrait diagram of land model (30) with , .

Figure 14: Portrait diagram of land model (30) with , .

Figure 15: Portrait diagram of land model (30) with , .

Figure 16: Portrait diagram of land model (30) with , .

5. Conclusion

The problem of bifurcation for a fractional land model with time delay has been studied in this paper. Firstly, by comprehensively analyzing the characteristic equation, the stability criterion and the conditions of existence of Hopf bifurcation are established. It was shown that the presence of delay heavily affects the stability of the model. Finally, according to China, Yunan Province, Yuxi city land data, some simulation examples have been exploited to illustrate the applicability and usefulness of the developed theoretical result.

Data Availability

The labeled dataset used to support the findings of this study is available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant nos. 61903323 and 11361072, the Natural Scientific Research Fund Project of Yunnan Province under Grant nos. 2017FH001-045 and 2018FH001-012, and the National Key Research and Development Project under Grant no. 2016YFC0502406.

Supplementary Materials

In this paper, in “Section 4,” according to Yuxi city land and population data of Yunan Province, China, that is to say, from numerical simulation of data (see attachment of data.txt, basic data), we have a = 0.34, b = 0.012, r = 0.0048, s = 0.004, h = 0.05, and d = 0.08. By computing, we obtain the positive equilibrium point as (0.075617, 0.021169, 0.44984, 0.423379). Moreover, we have , ω₀ = 0.002358, and τ₀ = 728.403. Furthermore, by using software of mathematics, we obtain fitting results, that is, it is in the attachments data1.txt and data2.txt (data after numerical simulation). As a matter of fact, this shows that numerical simulation is consistent with the actual, see Figures 1–8 and Figures 9–16. (Supplementary Materials)

References

B. W. Huang, “On earth system science and sustainable development strategy (I),” Acta Geographica Sinica, vol. 51, pp. 350–354, 1996, in Chinese.

W. V. Reid, D. Chen, L. Goldfarb et al., “Earth system science for global sustainability: grand challenges,” Science, vol. 330, no. 6006, pp. 916-917, 2010.

B. L. Turner II, “Local faces, global flows: the role of land use and land cover in global environmental change,” Land Degradation and Development, vol. 5, no. 2, pp. 71–78, 1994.

X. B. Li, “A review of the international researches on land use/land cover change,” Acta Geographica Sinica, vol. 51, pp. 553–557, 1996, in Chinese.

A. M. Dewan and Y. Yamaguchi, “Land use and land cover change in greater Dhaka, Bangladesh: using remote sensing to promote sustainable urbanization,” Applied Geography, vol. 29, no. 3, pp. 390–401, 2009.

H. J. Geist and E. F. Lambin, “Proximate causes and underlying driving forces of tropical deforestation,” Bioscience, vol. 52, no. 2, pp. 143–150, 2002.

S. Rost, D. Gerten, A. Bondeau, W. Lucht, J. Rohwer, and S. Schaphoff, “Agricultural green and blue water consumption and its influence on the global water system,” Water Resources Research, vol. 44, pp. 137–148, 2008.

L. H. Hou and Y. L. Cai, “An essential analysis and review on land use/cover change research,” Progress in Geography, vol. 23, pp. 96–103, 2004, In Chinese.

E. O’Rourke, “Socio-natural interaction and landscape dynamics in the Burren, Ireland,” Landscape and Urban Planning, vol. 70, no. 1-2, pp. 69–83, 2005.

J. Ni and J. A. Shao, “The drivers of land use change in the migration area, three gorges project, China: advances and prospects,” Journal of Earth Science, vol. 24, no. 1, pp. 136–144, 2013.

N. K. Msofe, L. X. Sheng, and J. Lyimo, “Land use change trends and their driving forces in the Kilombero valley floodplain, Southeastern Tanzania,” Sustainability, vol. 11, no. 2, p. 505, 2019.

J. Kleemann, G. Baysal, H. N. N. Bulley, and C. Frst, “Assessing driving forces of land use and land cover change by a mixed-method approach in north-eastern Ghana, West Africa,” Journal of Environmental Management, vol. 196, pp. 411–442, 2017.

D. Y. Liu, X. Q. Zheng, and H. B. Wang, “Land-use simulation and decision-support system (landSDS): seamlessly integrating system dynamics, agent-based model, and cellular automata,” Ecological Modelling, vol. 417, pp. 1–10, 2020.

Y. Q. Chen, H. V. Peter, and B. Xu, “Spatial modeling of land use and its effects in China,” Progress in Geography, vol. 19, pp. 116–127, 2000, in Chinese.

J. Zhang, W. B. Zhu, S. Y. Wu, and S. C. Li, “Simulation of temporal and special land use changes in Jing-Jin-Ji urban agglomeration using CLUE-S model,” Acta Scientiarum Naturalium Universitatis Pekinensis, vol. 54, pp. 115–124, 2018, in Chinese.

J. J. Liu, M. Xia, Y. Z. Liu, K. L. Zhang, and Z. H. Zhang, “Driving mechanism of rural land use change based on multi-agent system and cellular automata,” Transactions of the Chinese Society of Agricultural Engineering, vol. 34, pp. 242–252, 2018, in Chinese.

A. P. Dobson, A. D. Bradshaw, and A. J. M. Baker, “Hopes for the future: restoration ecology and conservation biology,” Science, vol. 277, no. 5325, pp. 515–522, 1997.

S. S. Chen, J. Hua, F. X. Zhang, and Y. M. Li, “Dynamical analysis of land model with Holling II land reclamation rate,” Journal of Systems Science and Complexity, vol. 37, pp. 819–827, 2017.

I. Kubiaczyk and S. H. Saker, “Oscillation and stability in nonlinear delay differential equations of population dynamics,” Mathematical and Computer Modelling, vol. 35, no. 3-4, p. 295, 2002.

K. Gopalsamy and I. Leung, “Delay induced periodicity in a neural netlet of excitation and inhibition,” Physica D, vol. 89, no. 3-4, pp. 395–426, 1996.

Y. H. Xia and V. G. Romanovski, “Bifurcation analysis of a population dynamics in a critical state,” Bulletin of the Malaysian Mathematical Sciences Society, vol. 38, no. 2, pp. 499–527, 2015.

Z. Wei, Y. H. Xia, and T. H. Zhang, “Stability and bifurcation analysis of a amensalism model with weak Allee effect,” Qualitative Theory of Dynamical Systems, vol. 19, no. 1, pp. 19–23, 2020.

X. D. Li, J. H. Shen, and R. Rakkiyappan, “Persistent impulsive effects on stability of functional differential equations with finite or infinite delay,” Applied Mathematics and Computation, vol. 329, pp. 14–22, 2018.

X. Y. Yang, X. D. Li, X. Li, Q. Xi, and P. Y. Duan, “Review of stability and stabilization for impulsive delayed systems,” Mathematical Biosciences & Engineering, vol. 15, no. 6, pp. 1495–1515, 2018.

X. D. Li, X. Y. Yang, and T. W. Huang, “Persistence of delayed cooperative models: impulsive control method,” Applied Mathematics and Computation, vol. 342, pp. 130–146, 2019.

D. Yang, X. D. Li, J. Shen, and Z. J. Zhou, “State-dependent switching control of delayed switched systems with stable and unstable modes,” Mathematical Methods in the Applied Sciences, vol. 41, no. 16, pp. 6968–6983, 2018.

D. Yang, X. D. Li, and J. L. Qiu, “Output tracking control of delayed switched systems via state-dependent switching and dynamic output feedback,” Nonlinear Analysis: Hybrid Systems, vol. 32, pp. 294–305, 2019.

Y. Q. Wang, J. Q. Lu, and Y. J. Lou, “Halanay-type inequality with delayed impulses and its applications,” Science China: Information Sciences, vol. 62, no. 9, Article ID 192206, 2019.

E. Ahmed, A. M. A. El-Sayed, and H. A. A. El-Saka, “Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models,” Journal of Mathematical Analysis and Applications, vol. 325, no. 1, pp. 542–553, 2007.

C. D. Huang, J. D. Cao, M. Xiao, A. Alsaedi, and F. E. Alsaadi, “Controlling bifurcation in a delayed fractional predator-prey system with incommensurate orders,” Applied Mathematics and Computation, vol. 293, pp. 293–310, 2017.

S. Das, P. K. Gupta, and Rajeev, “A fractional predator-prey model and its solution,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 7, pp. 873–876, 2009.

H. A. El-Saka, E. Ahmed, M. I. Shehata, and A. M. A. El-Sayed, “On stability, persistence, and Hopf bifurcation in fractional order dynamical systems,” Nonlinear Dynamics, vol. 56, no. 1-2, pp. 121–126, 2009.

C. Luo and X. Y. Wang, “Chaos generated from the fractional-order complex chen system and its application to digital secure communication,” International Journal of Modern Physics C, vol. 24, no. 04, p. 1350025, 2013.

E. Kaslik and S. Sivasundaram, “Nonlinear dynamics and chaos in fractional-order neural networks,” Neural Networks, vol. 32, no. 1, pp. 245–256, 2012.

Y. Li, H. P. Wang, and Y. Tian, “Fractional-order adaptive controller for chaotic synchronization and application to a dual-channel secure communication system,” Modern Physics Letters B, vol. 33, no. 24, pp. 1097–1106, 2019.

L. X. Yang, W. S. He, and X. J. Liu, “Synchronization between a fractional-order system and an integer order system,” Computers & Mathematics with Applications, vol. 62, no. 12, pp. 4708–4716, 2011.

M. Odai and Y. Hori, “Controller design robust to nonlinear elements based on fractional order control system,” IEEE Transactions on Industry Applications, vol. 120, no. 1, pp. 11–18, 2010.

C. D. Huang, H. Liu, X. P. Chen, M. S. Zhang, J. D. Cao, and A. Alsaedi, “Dynamic optimal control of enhancing feedback treatment for a delayed fractional order predator-prey model,” Physica A: Statistical Mechanics and Its Applications, vol. 554, Article ID 124136, 2020.

Z. Wang, Y. K. Xie, J. W. Lu, and Y. X. Li, “Stability and bifurcation of a delayed generalized fractional-order prey-predator model with interspecific competition,” Applied Mathematics and Computation, vol. 347, pp. 360–369, 2019.

C. D. Huang, Z. H. Li, D. W. Ding, and J. D. Cao, “Bifurcation analysis in a delayed fractional neural network involving self-connection,” Neurocomputing, vol. 314, pp. 186–197, 2018.

D. Matignon, “Stability results for fractional differential equations with applications to control processing,” in Proceedings of the Computational Engineering in Systems and Application Multi-Conference, IMACS, vol. 2, pp. 963–968, Lille, France, 2016.

X. Li and R. C. Wu, “Hopf bifurcation analysis of a new commensurate fractional-order hyperchaotic system,” Nonlinear Dynamics, vol. 78, no. 1, pp. 279–288, 2014.

N. Laskin, “Fractional quantum mechanics,” Physical Review E, vol. 62, no. 3, pp. 3135–3145, 2000.

I. Podlubny, Fractional Differential Equations, Academic Press, Cambridge, MA, USA, 1999.

A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Application of Fractional Differential Equations, Elsevier, New York, NY, USA, 2006.

T. Zhao, Z. Z. Zhang, and R. K. Upadhyay, “Delay-induced Hopf bifurcation of an SVEIR computer virus model with nonlinear incidence rate,” Advances in Difference Equations, vol. 2018, no. 1, 2018.

E. Hanert, E. Schumacher, and E. Deleersnijder, “Front dynamics in fractional-order epidemic models,” Journal of Theoretical Biology, vol. 279, no. 1, pp. 9–16, 2011.

C. Vargas-De-León, “Volterra-type Lyapunov functions for fractional-order epidemic systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 24, no. 1–3, pp. 75–85, 2015.

H. AbdelhamedAsfour and M. Ibrahim, “On the differential fractional transformation method of MSEIR epidemic model,” International Journal of Computer Applications, vol. 113, no. 3, pp. 10–16, 2015.

W. Zhu, Y. H. Xia, B. Zhang, and Y. Bai, “Exact traveling wave solutions and bifurcations of the time fractional differential equations with applications,” International Journal of Bifurcation and Chaos, vol. 29, no. 3, Article ID 1950041, 2019.

A. S. Deshpande, V. Daftardar-Gejji, and Y. V. Sukale, “On Hopf bifurcation in fractional dynamical systems,” Chaos, Solitons & Fractals, vol. 98, pp. 189–198, 2017.

H. Wang, Y. Yu, G. Wen, and S. Zhang, “Stability analysis of fractional-order neural networks with time delay,” Neural Processing Letters, vol. 42, no. 2, pp. 479–500, 2015.

M. Xiao, W. X. Zheng, G. Jiang, and J. D. Cao, “Undamped oscillations generated by Hopf bifurcations in fractional-order recurrent neural networks with caputo derivative,” IEEE Transactions on Neural Networks and Learning Systems, vol. 26, no. 12, pp. 3201–3214, 2015.

S. Bhalekar and D. Varsha, “A predictor-corrector scheme for solving nonlinear delay differential equations of fractional order,” Journal of Fractional Calculus and Applications, vol. 1, pp. 1–9, 2011.