Dynamical Models of Interactions between Herds Forage and Water Resources in Sahelian Region

and Applied Analysis 3 resources consumed by an animal in the herd per unit time, is ecologically plausible and provides a solid basis for theoretical development. As in [11–13], density dependence of the resources will be the starting point, which gives a functional response. In the simplest case, such a function is a linear function of the forage resources density (R), which is integrated into the classic Lotka-Volterra model. The linear functional response is a limited case and can only be seen over a short period. But one can use the Michaelis-Menten or Holling response function of type II defined by φ(R) = (Bω 0 R/(1 + Bω 1 R)), where R denotes the forage resources density, ω 0 and ω 1 denote, respectively, the time taken by a herd’s animal to search and consume forage resources, and B is the herd consumption rate per unit of time. With the Holling function response of type II, it is well known that the diminution of forage resources due to the herd consumption increases and the forage density decreases and becomes constant at the end. In the model formulated here, γ is the removal rate from the herd (death, off take. . .); parameter e denotes the conversion rate of the forage resources consumed by the herd’s animals into their biomass. When there is no herd, the dynamics of forage resources can be governed by the logistic equation, but the intrinsic production can be negative, such that (0,0) is a stable equilibrium. Setting H(t) as the tropical livestock unit of the herd at day t, the function g(R,H) = Bω 0 R(t)H(t)/(1 + Bω 1 R(t)) can also be written as g(R,H) = qR(t)H(t)/(1 + aR(t)), where q = Bω 0 denotes the utilization rate of forage resources and a = Bω 1 denotes the satiety rate of herd’s animals. The following differential equations model interactions between forage resources and animals of the herd:


Introduction
It is well known that pastoralism is the name given to the subsistence practice in which people care for and domesticate animals, usually ungulates such as camels, cattle, reindeer, sheep, and goats [1].
Pastoralism will continue for the near future in poor nations [2][3][4], especially in Africa, because it is generally an efficient, low energy requiring subsistence base for semiarid regions.During the 20th century, however, most national governments tried to force pastoralists to stop their migrations and to reduce the size of their herds in order to prevent overgrazing.These efforts at controlling them were consistently resisted by pastoralists.They usually saw large herds as symbols of wealth and as security against unpredictable climates and periodic epidemics among their animals [5].
Livestock in the Sahel is now basically a usurper such as land use regarding the vegetation that grows there naturally.The pastoral farming systems are those in which more than 90% of the dry matter consumed by livestock comes from grazing.The researchers focus in this paper on these pastoralfarming systems.The transition from one system to another is periodic and depends on climate and economic context [6].Part of pastoral population in the rural population varies greatly from country to country.Livestock numbers are not known with precision, because large part of statistics does not differentiate production systems.Based on national statistics, thus including agricultural areas, small ruminants are the most numerous, before cattle.Animal density can only be determined on well-defined regions, for which there are the number of both animals and areas used by livestock.The data always differ by authors.
The pastoral system is moving towards a systems approach mixed crop-livestock farming.This trend goes in North Africa with the grazing of steppe.In the Sahel, the environmental consequences are less obvious [2,7,8].The vast majority of feed consumed by ruminants in Sahelian region is still formed by natural pastures.Savannas, steppes, and training fallow provide the basis of livestock feed, even 2 Abstract and Applied Analysis in sedentary farming or being intensified.It is important to evaluate or to monitor their progress [1,6].
There have been many problems about pastoral livelihood vulnerability.The dry lands cannot support sustained and reliable agriculture because of low and variable rainfall and high temperatures.Pastoralism, however, is extremely well suited to this type of environment.Pastoralists make optimum use of the dry lands by practicing a mobile and extensive livestock-keeping system.They move according to where and when fodder becomes available and use different herd management strategies such as herd splitting, herd diversification, and herd maximization to ensure that they spread the risk of livestock loss from droughts, diseases, and theft.All the while, they make maximum use of the available vegetation without degrading the environment [1].
Different explanations have been advanced for the increasing vulnerability of pastoralists.Population growth in pastoral areas has created pressure on land.Climate change has increased the frequency of droughts, floods, and livestock diseases.However, these natural factors only exacerbate the effects of a harsh policy and legal environment that is focused on "modernizing" and settling pastoralists.The increased vulnerability of pastoral livelihoods to shocks and other drivers of change are in many ways a function of the cumulative effect of these policies [1].
Long-term sustainable pastoral development requires a good knowledge of the dynamics of multiple factors underlining pastoralism and here research has a crucial role to play.Alongside support to improve primary and secondary education is the need to strengthen institutions of higher learning and research in pastoral areas of Africa [5].Such institutions require support to conduct research on a range of political, social, and natural science subjects and at levels, which range from local adaptation to regional integration and global trends.The links between research and policy also need to be strengthened, so that policy responds to the dynamics of pastoral livelihoods in an efficient manner [5].This paper is our contribution to this research in pastoral areas of Africa.

Model Construction
In Africa, the Sahel has a wet season (June to October) and a dry season (November to May) which are very distinct, and livestock is greatly influenced by the amount of grass and shrubs.As we said previously, more than 90% of dry matter consumed by livestock comes from grazing.Some problems as overgrazing or scarcity of water are frequently observed.These problems force some pastoralists to migrate with their herd.According to these facts, the intrinsic production of forage resources can be a negative number, even the intrinsic disponibility of water resources.The intrinsic production  of forage resources can be written as  =  0 −  1 , where  0 =  01  +  02 as given in [8],  denotes the precipitation,  02 is the growth rate of the grasses and others without precipitation, and  1 is the disappearing rate due to many reasons (competition, respiration, and human harvesting) different from herd's nutrition.The particularity of this model, which is the difference with classic predator-prey models in the literature, is the fact that the intrinsic production rate can be negative ( < 0).This situation is frequently observed in Sub-Saharan Africa, particularly in Sahelian region when precipitations are scarce.Intensification of this situation can conduct pastoralists to migrate with their herds: this is the transhumance.We also know [9] that tropical livestock unit (TLU) is a basic criterion for a head of cattle weighing 250 kg.The daily volume of consumption of dry matter per TLU amounted to 6.5 kg.However, since this criterion varies, there is obviously the case with differences depending on the country, even within each country.The test described here is also currently used officially by the Ministry of Agriculture and Livestock in Niger [9].
Parameter  denotes the maximal number of livestock in this paper.When considering the question of the number of livestock that can be raised on a given surface, it is first important to know the amount of forage needed for livestock.To calculate the volume needed to feed livestock, measured by the volume of solids, the amount of assimilation of dry materials by livestock is generally estimated to be approximately 1.4% to 3.0% of the weight of livestock.In the Sahelian region, by experience, the following values are generally used on the basis of tropical livestock unit, in order to compute the daily volume of required dry matter per head [9].
One cattle: 5.2 kg [6.5 kg (TLU) × 0.8 (index of the considered space)].Two goats and sheep: 1.0 kg [6.5 kg (TLU) × 0.15 (index of the considered space)].We can calculate the number of heads that can be raised with the formula below [9].Rearing capacity = [(volume of feed supply of natural grassland × utilization rate) divided by (daily volume required per capita solids)] × number of days of reliance on natural grassland.
Let us give some hypotheses concerning the models in this paper.
(H1) The herds, forage, and water resources considered in this paper are in the same area.
(H2) When the intrinsic production is nonnegative, if feed resources are not consumed, they change their amount to the maximal production capacity and therefore storage.In this case the feed resources grow logistically.When the intrinsic production is negative, the feed resources will disappear at long time, and this situation conducts pastoralists to migrate with their herds: this corresponds to the transhumance situation.
(H3) If the herd is helpless, its number of animals or its tropical livestock unit (TLU) decreases and may disappear if nothing is done.
(H4) Interactions between the herd and forage resources are following the functional response of Michaelis-Menten or Holling function response of type II.
(H5) Consumption of resources has an instant effect on the reduction of forage resources and increased biomass of the herd, in proportion to their consumption.
In the study of interactions between herds, forage, and water resources, it is crucial to determine which specific form of functional response, describing the amount of resources consumed by an animal in the herd per unit time, is ecologically plausible and provides a solid basis for theoretical development.As in [11][12][13], density dependence of the resources will be the starting point, which gives a functional response.In the simplest case, such a function is a linear function of the forage resources density (), which is integrated into the classic Lotka-Volterra model.The linear functional response is a limited case and can only be seen over a short period.But one can use the Michaelis-Menten or Holling response function of type II defined by () = ( 0 /(1 +  1 )), where  denotes the forage resources density,  0 and  1 denote, respectively, the time taken by a herd's animal to search and consume forage resources, and  is the herd consumption rate per unit of time.With the Holling function response of type II, it is well known that the diminution of forage resources due to the herd consumption increases and the forage density decreases and becomes constant at the end.In the model formulated here,  is the removal rate from the herd (death, off take. ..); parameter  denotes the conversion rate of the forage resources consumed by the herd's animals into their biomass.When there is no herd, the dynamics of forage resources can be governed by the logistic equation, but the intrinsic production can be negative, such that (0,0) is a stable equilibrium.Setting () as the tropical livestock unit of the herd at day , the function (, ) =  0 ()()/(1 +  1 ()) can also be written as (, ) = ()()/(1 + ()), where  =  0 denotes the utilization rate of forage resources and  =  1 denotes the satiety rate of herd's animals.The following differential equations model interactions between forage resources and animals of the herd: According to a study conducted by the Japan Green Resources Corporation (JGRC) [9], in a region where the average annual precipitation is 500-600 mm, the volume of dry matter production of grass per hectare in natural grassland was in tons 1.54 in 1997, 1.6 in 1998, and 1.8 in 1999 (according to research conducted on mowing in October each year).Among wild herbs growing in this region are the grasses that cattle prefer, but they also appreciate the taste of some legumes.
Until now, in many cases (excluding reservoir dams) small artificial wells and natural ponds were used to supply livestock with drinking water.In all these cases, the storage capacity of water is low, and we find ourselves in a situation of chronic lack of water during the dry season.In addition, the storage capacity of water ponds in some places decreases by the reduction of vegetation in the water accumulation area and the influx of land in the ponds.There is sometimes also conflict when cattle in two villages must share the same drinking water source [9].It is well known that it is not easy to get water in Sahelian region.For sustainable use of water sources that are ponds, it is also necessary to take a measure of water conservation preventing the influx sand around ponds.For example, one study of the Japan Green Resources Corporation showed that establishment of lines of stones can retain the sand and promote the restoration of wild vegetation.The earth is accumulating inside line of stones (the side from which comes the influx); the storage capacity of water becomes higher and ability to grow fodder crops and planting fodder trees becomes larger.This allows obtaining fodder and shade for livestock, and it becomes possible to revive the natural clam that was impoverished around the pond, as a grazing area including a watering place for the cattle [9].
The pool of Eda (which is straddled between the villages of Magou and Eda in Niger) had long been used as a point of watering the cattle moving to the Burkina Faso, but due to the decrease in its subsequent storage capacity of water to the influx of sand in the pool, it could not provide the volume of water required to move livestock and livestock of the two villages.But thanks to the digging of the pond and the establishment of lines of stones (some with citizen participation), the revegetation of weeds is increasing, and livestock numbers even greater than in the past are not gathering [10].
The variable () is the tropical livestock unit of the herd at time , () is the forage resources density (kilograms of dry matter per hectare) at time , () denotes the water resources at time ,  is the removal rate from the herd,  1 is the conversing rate of forage resources consumed into herd biomass, and  2 is the conversing rate of water resources consumed into herd biomass; as previously, is the intrinsic production of forage resources and 1 the intrinsic disponibility of water;  is the precipitation rate in the area;  1 is the maximal capacity of the considered area to support forage resources;  2 is the maximal capacity of the area to contain large quantity of water;  1 and  2 are, respectively, the utilization rates of forage and water resources.
Following [14,15], we set   = ( 0 / 1 ) and   = (1/ 1 ).The dynamical model for interactions between herds, forage resources, and water is given by the following system of differential equations: We therefore have three trophic levels: one predator and two preys.The intrinsic production rates  1 and  2 defined, as previously, can be negative values.In the absence of forage and water resources ( 1 < 0 and  2 < 0), the herd is doomed if nothing is done.In this case, pastoralists will migrate with their herds (transhumance).Its workforce decreases and goes to extinction if nothing is done.But, in the absence of the herd, forage resources are stored and accumulate until the limit capacity when  1 > 0; forage resources can disappear when  1 < 0. The water resources can accumulate to exceed the carrying capacity, but it will very quickly go back to limit capacity when  2 > 0; water resources can disappear when  2 < 0.
Let us give now a mathematical analysis of System (1) and interpretations of the results in terms of pastoralism, using the data [10].
Remark 2. This theorem confirms the fact that forage resources quantity and the number of animals in the herd are always nonnegative numbers, since their trajectories are always in  2 + .Moreover, the forage resources quantity is bounded.This is ecologically plausible since temperatures are usually high and the stockage capacity is limited.

Lemma 3. The sign of the intrinsic production r has an incidence on the dynamics of forage resources and then on the dynamics of the herd's animals.
(1) If  > 0, then the forage resources converge to the maximal number of livestock K for a long time.
(2) If  < 0, then the forage resources will disappear for a long time if nothing is done.This situation conducts to the transhumance of pastoralists with their herds.
(1) Equilibrium  0 = (0, 0) and equilibrium  1 = (, 0) exist without any condition. ( ) is an ecologically acceptable equilibrium if  >  and the threshold Therefore, with these conditions satisfied, there are three equilibria for System (1). ( ) is an ecologically acceptable equilibrium if  >  and the threshold Therefore, with these conditions satisfied, there are three equilibria for System (1).
(1) The equilibrium  0 is a saddle-node when  > 0 and a stable node when  < 0.
(2) If  > 0,  1 is a globally asymptotically stable node when  1 < 1 and a saddle-node with stability for forage resources and instability for herds when  1 > 1.
If  < 0,  1 is a saddle-node when  1 < 1 and a locally unstable node when  1 > 1.  1 is a globally asymptotically stable node when  1 = 1.The centre manifold in this case is given by (3) The equilibrium  2 is not ecologically acceptable when , then  2 is a globally asymptotically stable focus and when  1 > 1 + (/(1 + )),  2 is an unstable focus and there exists a limit cycle for System (1).This phenomenon corresponds to supercritical Hopf bifurcation [16].
Remark 6.The stability of equilibrium  0 = (0, 0) means that there will be no more forage resources and no animal of the herd in the considered area for a long moment.There can be many explanations to this situation.Droughts can cause animals mortalities through starvation, emergency slaughtering, and sales, or definitive herds migrations (transhumance), which can create severe drops in the herd sizes.The Sahelian region is particularly affected by such climate shocks [17].Droughts can also cause forages' disappearance, and then the equilibrium  0 is stable.When this removal concerns only animals of the herd, depending on some climatic changes, the forage resources growth towards the maximal quantity needed for livestock and the equilibrium  1 = (, 0) is stable.
Concerning the coexisting equilibrium  2 = ( * ,  * ), there is a Hopf bifurcation.The pastoral interpretation of Hopf bifurcation is that animals in the herd will coexist with the forage resources, exhibiting oscillatory balance behavior.We have a peak for the herd trajectory, followed by a peak for forage trajectory.
(6) If  1 = 1, then  1 is a globally asymptotically stable node when  > 0 and just a locally asymptotically stable node when  < 0.
Remark 8 (Hopf bifurcation).There is a limit cycle when  1 passes through the value 1 + (/(1 + )).This phenomenon is known as Hopf bifurcation.Since the limit cycle is stable, it is a supercritical Hopf bifurcation.The pastoral interpretation of Hopf bifurcation as we said is that Herd's removal rate 1 per day [10] animals in the herd will coexist with the forage resources, exhibiting oscillatory balance behavior.This can be explained by the resilience of vegetation, which means the ability of the ecosystem to withstand unusual stress and recover spontaneously once they have disappeared.Thus, an investigation in [18] showed that cereal production in the Sahel went from a deficit of one million tons in 1987 to surplus 1 million tons in 1988.
Remark 10.Any trajectory with initial condition in the nonnegative orthant  3  + is trapped and will stay inside.The nonnegative orthant  3  + is then positively invariant and System (2) is mathematically well posed.We can then say that System (2) is well posed, since forage resources (), water resources (), and the tropical livestock unit of the herd () are always nonnegative quantities.

Proposition 12 (descartes criterion). The number of positive real roots of the polynomial equation is less than or equal to the number of changes in coefficient of ordered polynomial coefficients, and these two numbers have the same parity.
Proposition 13 (descartes rule of signs-I).The number of positive roots of the polynomial equation with real coefficients does not exceed the number of sign changes in its coefficients.A zero coefficient is not counted as a sign change.
(4) The equilibrium  13 = ( 1 , when it exists, can be locally asymptotically stable if the Routh-Hurwitz conditions are satisfied, where V 1 , V 2 , and V 3 are given as in the proof in Appendix B.Moreover, the system exhibits periodic oscillations with alternative peaks of forage, water, and the herd trajectories.
Remark 16.It is difficult for cattle to adequately live in Sahelian region when forage resources or water resources have dried up completely.When at least one of two completely lacks resources, pastoralists can be forced to migrate with livestock in order to avoid disaster.In addition, feed resources and water resources are often somewhat related since the fresh leaves have water content of about 80%.Transhumance can therefore take place when one of the two resources has failed to herd and in this case, pastoralists migrate while the second resource is available in the area.The stability of  10 ,  11 ,  12 , and  13 characterizes these situations.In usual predator-prey models, it is virtually impossible to have stability for  10 ,  11 , and  12 , becausestability of  10 , for example,means thatthe entire population will disappear, which is not the main objective when coupling preys and predators.However, this is easily explained in pastoralism as it corresponds to transhumance.

Numerical Simulations and Interpretations
The parameters values come from [9,10] and references therein, and we use the same method as in [9] to compute some parameters values.Let us recall these formulas, in order to have the livestock capacity in natural grassland.
(1) Assuming that the production volume of fresh grasses is 8000 kg/ha, the estimated water contained in fresh grasses is around 80% (this means 20% of dry matter).(2) The dry matter (DM) obtained for the livestock is 5.2 kg/day for cattle and 1 kg/day for goats and sheep.
We also have in Tables 1 and 2 the parameters values obtained in [9,10], which correspond to the situation where the intrinsic production  =  01  +  02 −  1 is negative.This means that, at this moment, the removal quantity of forage resources, which is due to competition or which is used to feed the livestock, is greater than the intrinsic production.This situation can force pastoralists to leave the area and migrate elsewhere in order to feed their herds.1) when  = 650 kgDM/ha,  = 0.05/day,  = 0.013/ha,  = 800 kgDM/ha,  = 0.4/day, and  = 1.0/day.This illustrates the global stability of the unique coexisting equilibrium  2 = ( * ,  * ).This means that the herd's animals and forage resources have periodic curves since there is a stable limit cycle.The two populations can exist together.

Conclusion
A mathematical model for pastoralism has been considered in order to explain and predict some situations in Sahelian region.This is one of the first deterministic models of pastoralism analyzed with interpretations.The results obtained in this paper have been simulated with the data in [9,10] and confirm many observations throughout pastoralism's  1 and 2. literature concerning interactions of forage resources and herds.At the end of this first part, we can say that since forage resources can be also considered sometimes as water resources, it will be better for us to add in the previous equations of System (1) another equation in order to take into account the water resources.
The statistics show the importance of livestock to the national economy in tropical arid countries.Two opposing trends develop: the continued contribution of livestock in national economy and the growing risk of deterioration of pastoral resource savings.
Sustainable pastoral resource management, equitable and secured access to pastoral resources, and peace and security which are always the aims of framework policy are also the guidelines of our contribution to pastoralism by mathematical models.Thus, it is necessary that these issues be addressed through a comprehensive pastoral policy, which confers full political, social, economic, and environmental benefits to the pastoral communities.The deterministic models proposed and analyzed in this paper are firstly to participate to these efforts.Although considerable efforts with some positive results are being made throughout Africa, a great deal more still needs to be done.This calls for both commitments from individual countries and cooperation at the regional and continental levels.We think that mathematical models can considerably contribute to increasing comprehension of pastoralism and therefore increase the quality of life of pastoralists.The models considered in this paper exhibit periodic oscillations when some conditions are satisfied.The comprehension of conditions to have periodic oscillations and then coexistence of herds, forage, and water resources can help pastoralists to better understand their environment.With this coexistence, there are alternative peaks of livestock tropical unit, forage resources, and water resources.This means that when these conditions are satisfied, pastoralists can live in some areas without transhumance, with disponibility of forage and water resources.

A. Expression of coefficients in Theorem 11
The points  10 ,  11 ,  12 , and  13 are obviously equilibria of System (2) without any condition.The coexisting equilibrium is given by the system The third equation implies The expressions of  * and  * are positives if  * >  1 ,  * >  2 , and  1  1 <  <  2  2 .In order to give the expression of  * , let us consider the second equation of the system
(5) At equilibrium ( * ,  * ,  * ), using the relations in System (2)  ) . (B.1) The characteristic polynomial is in the form where (B. 3) The Routh-Hurwitz conditions for stability of this equilibrium are Abstract and Applied Analysis 13 When these conditions are satisfied, a coexisting equilibrium when it exists is locally asymptotically stable and globally asymptotically stable if there is a unique coexisting equilibrium.This ends the proof.

Figure 1 :
Figure1: Phase portrait of the herd's animals and forage resources when r < 0. This illustrates the global stability of equilibrium  0 = (0, 0) for System (1).This is a new situation, which happens during migration of pastoralists with their herds, calling transhumance; the pastoralists leave the area because of scarcity of forage resources.The parameters in this case are given in Table1.