Depletion of Forest Resources and Wildlife Population with Habitat Complexity : A Mathematical Model

Nowadays, management and regulation of natural resources like agriculture, fisheries, forestry and wildlife is one of the popular topics in research. The evolution of humankind is largely dependent on the quality of the environment and the resources it provides; but numerous human-induced factors, and climate change may drastically alter the conditions of human sustainability. This paper deals with effect of numerous human-induced activities on the depletion of forestry resources and wildlife population with habitat complexity. A nonlinear mathematical model is proposed and analyzed. In modeling process, we assume that the growth rate of wildlife population wholly depends on forestry biomass. It is depleted by human-induced activities. Local stability analysis of the mathematical model along with the persistence of the system is checked by using theory of nonlinear ordinary differential equations and Butler-McGhee lemma. Analytical results obtained are justified numerically through numerical simulation. Important parameters are investigated and variation of variables with change in these parameters is determined.


Introduction
India is known for its rich heritage of biological diversity, having already documented over 91,000 species of animals and 45,500 species of plants in its ten bio-geographic regions.Nearly 6500 native plants are still used prominently in indigenous healthcare systems.Thousands of locally-adapted crop varieties, grown traditionally since ancient times, and nearly 140 native breeds of farm livestock, continue to thrive in its diversified farming systems.The country is recognized as one of the eight Vavilovian Centres of Origin and Diversity of Crop Plants, having more than 300 wild ancestors and close relatives of cultivated plants still growing and evolving under natural conditions.But Forests are being destroyed and degraded at alarming rates.Deforestation comes in many forms, including fires, clear-cutting for agriculture, ranching and development, unsustainable logging for timber, and degradation due to climate change.This impacts people's livelihoods and threatens a wide range of plant and animal species.Forests are disappearing at an alarming rate-18.7 million acres of forests are lost annually, equivalent to 27 soccer fields every minute.
Mathematical modeling is an essential tool in studying a diverse range of such ecological situations.Proposed work deals with the study of mathematical model on the sustainable management of renewable resources in India.The depletion of forest biomass by human population and industrialization has been investigated both theoretically as well as experimentally by many researchers ([1]- [7]).In particular, [8] [9] have proposed a mathematical model for the depletion of forest resources where they have shown that the forest resources may be doomed to extinction under growing population living in the forest or by industrialization which wholly depends on forest resources.[8] [9] results have also been found by other researchers showing their concern about biodiversity caused by deforestation.[10] has proposed and analyzed a mathematical model for conservation of forestry biomass and wildlife population.In above paper, we assumed that the growth rate of wildlife conservation is proportional to the depletion of forestry biomass due to wildlife population.We have also studied the affect of illegal trade in forestry biomass and wildlife population.[11] has proposed and analysed the mathematical model and shown the effect of population and industrialization on the forestry biomass and wildlife population.This model is motivated by our food chain model given by [12].We proposed a nonlinear mathematical model and analyzed to see the effect of harvesting on food chain model with habitat complexity.We again study this model in the view of effect of human-induced activities on the forestry biomass and wildlife population with habitat complexity.
The rest of this paper is organized as follows: in Section 2, we introduce our mathematical model.In Section 3 and 4, our model is analyzed with regard to equilibria and their stabilities, respectively.In Section 5, we find the conditions for persistence of the system.In Section 6, we present a numerical example to illustrate the applicability of results obtained and also investigate the occurrence of Hopf bifurcation for certain set of parameters.We conclude with a short discussion in Section 7.

Mathematical Model
The mathematical model for effecting human-induced activities on the forestry biomass and wildlife population with habitat complexity can be described by the ( ) ( ) ( ) ( ) The flowchart of the model (2.1) as shown in Figure 1.
Here, B be the forestry biomass, x W be the density of wild grazer population and y W be the density of wild predator population.It is assumed that the dynamics of forestry biomass follow logrithimic growth rate.It is further assumed that the rate of depletion of B due to human-induced activities is proportional to the product of forestry biomass and human-induced activities coefficient 1 q .As wild grazer population is wholly dependent on forestry biomass, its growth rate is proportional to the interaction term ( ) ( ) The depletion rate of wild grazer population due to human activities is proportional to the product of concentration of x W and human activities coefficient 2 q .The growth rate of wild predator population is proportional to the interaction term (i.e.

( ) ( )
).The natural depletion rate of wild predator population is considered proportional to y W .The depletion rate of wild predator population due to human activities is proportional to the product of concentration of y W and human activities coefficient 3 q .0 r is initial growth rate of the forestry biomass, K is the carrying capacity of forestry biomass, 1 a is the attack rate of grazers, 1 b is the half-saturation constant, ξ is the forestry bio- mass-grazer conversion rate, η is the zero population growth grazer intake, 2 a is the saturation killing rate, θ is the prey (grazer)-predator conversion rate, µ is the zero population growth predator intake, h is the handling time and 0 1 c < < is a dimension less parameter that measures the degree or strength of habitat complexity which is effected by human-induced activities like industrialization, development, framing etc.It is to be noted, when when there is no complexity, we get back the original Holling Type II response function.It may be pointed out that for feasibility of the model (2.1), the growth rates of wild grazer population should be positive.
To analyze the model (2.1), we need the bounds of dependent variables involved.For this we find the region of attraction in the following lemma.
Lemma 2.1: The set ( ) 0 , , : 0 is the region of attraction for all solutions initiating in the interior of the positive octant.Proof: Let x y B t W t W t be any solution with positive initial conditions ( ) Define a function Computing the time derivative of ( ) V t along the solutions of a system (2.1), we get Applying a theorem of differential inequalities, we obtain Therefore all solutions of system (2.1) enter into the region ( ) 0 , , : 0 .
This completes the proof of lemma.

Equilibrium Analysis
There exist following four equilibria of the system (2.1) 1) , where Existence of ( ) , x B W are the positive solutions of the system of algebraic equations given below Here * * * , , x y B W W are the positive solutions of the system of algebraic equations given below After solving (3.2), we have conditions of existence for equilibrium 3 E as follows and * B is the unique positive root of the following equation ( )

Local Stability Analysis
To discuss the local stability of the system (2.1), we compute the variational matrix of the system (2.1).The entries of the general variational matrix are given by differentiating the right side of the system (2.1) with respect to , x B W and y W , i.e.
( ) Accordingly, the linear stability analysis about the equilibrium points , = 0,1, 2, 3 i E i gives the following results: 1) The equilibrium point 0 E is unstable manifold in B-direction.
2) The equilibrium point 1 E is unstable manifold in W x -direction.
3) The equilibrium point 2 E is unstable manifold in W y -direction.
4) The characteristic equation corresponding to variational matrix ( ) where
Proof: Suppose u is a point in the positive octant and ( ) O u is the orbit through u and ω is the omega limit set of the orbit through u Note that ( ) We claim that ( ) ( ) ( ) Next, we show that ( ) Finally, since only the closed orbits and the equilibria form the omega limit set of the solutions on the boundary of 3 R + and system (2.1) is deceptive, by Butler-McGhee lemma , implies system is uniformly persist.

Numerical Simulation for Bifurcation
In this section, we present numerical simulation to illustrate results obtained in previous sections.The system (2.1) is solved using fourth order Runge-Kutta Method with the help of MATLAB software package.
We shall vary c in system (2.1) so as to obtain a Hopf bifurcation.Now, we write autonomous system (2.1) in the form ( ) x y n B W W k r K a b q H a c h q q ξ η θ µ
We consider the system ( ) ( ) ( ) ( ) The system (6.2) has always non-negative equilibrium point ( ) All eigenvalues have negative real parts, so equilibrium point is locally asymptotically stable at 0.8 c = .
From Figures 2-6, we analyse that habitat complexity play important role in stability of ecosystem.Human-induced activities effect the habitat complexity.
We also see that from Figure 3 and Figure 4, increasing the value of coefficient of human-induced activities, forestry biomass become extinct and wild life wholly depended on the forestry biomass also become extinct.It shows numerically that the Hopf point is found when 0.7 c = . 3 E is unstable when 0.7 c < and stable when 0.7 c > . From this result, we can say that if we conserve habitat complexity for wildlife population, we can reduce the rate of depletion of forestry biomass.Human being depend on the forestry biomass but if they use forestry biomass keeping in the range of control of habitat complexity of wildlife population, we conserve the forestry biomass as well as wildlife population.The numerical study presented here shows that, using the parameter c as control, it is possible to break the stable behavior of the system (6.2) and drive it in an unstable state.Also, it is possible to keep the population levels in a required state using the above control.

Conclusion
In this paper, we have proposed a nonlinear mathematical model to study the     effect of human activities on forestry biomass and its effect on wildlife species living in the forest.The proposed model has been analyzed by the stability theory of differential equations.The conditions of existence of equilibrium points and their stability in have been obtained.We have also obtained conditions under which system persists by using differential inequality.Our numerical study shows that habitat complexity behaves as a control; it is possible that it breaks the stable behavior of the system (2.1) and drives it to an unstable state.
Also, it is possible to keep the population levels in a required state using the above control.We see that effect of human-induced activities on forestry biomass and wild grazer population and the wild predator population.The effect of depletion of forestry biomass on the survival of wildlife species has also been shown graphically.
following system of nonlinear differential equations.
positive octant and the system (2.1) are persistent.