Predator-Prey Model with Ratio-dependent Food
|
|
- Rolf Hampton
- 6 years ago
- Views:
Transcription
1 University of Minnesota Duluth Department of Mathematics and Statistics Predator-Prey Model with Ratio-dependent Food Processing Response Advisor: Harlan Stech Jana Hurkova June 2013
2 Table of Contents 1. Introduction Known models Functional response Rosenzweig-MacArthur model Model Development Rescaling Nullclines Prey nullcline Predator nullcline Saddle node bifurcation Transcritical bifurcation Equilibria Extinction Monoculture equilibrium Coexistence equilibrium Hopf bifurcation Solution Homoclinic bifurcation Comparison to Rosenzweig-MacArthur Model Suggestions for Future Research Table of figures References Appendix A - Saddle node bifurcation Appendix B - Phase plane plots Appendix C - Hopf bifurcation
3 Acknowledgments I would like to express the deepest appreciation to my committee chair, Professor Harlan Stech, who has been the best advisor and professor I have ever had. Thank you for your guidance, patience, valuable advices and grammar checks. You introduced me the interesting side of mathematical modeling. I would like to also thank my committee members, Professor Bruce Peckham and Professor John Pastor for serving on my committee. In Addition, a thank you to Professor Dalibor Froncek, whose friendly approach makes all students feel home. You are a big support and the way you care about all students is amazing.
4 Abstract Predator-Prey models have been studied for a long time. The most famous ones are the Lotka- Volterra model, the Lotka-Volterra model with logistic function, and the Rosenzweig-MacArthur model. The properties of these models are well known. The purpose of this paper is to introduce a new model based on the Rosenzweig-MacArthur model using a new functional response that depends on the ratio of the predator and the prey populations. The model is analyzed and simulated, and its properties compared to the Rosenzweig-MacArthur model. The new model shows a surprising variety of behaviors, including multiple coexistence steady states, simultaneously attracting periodic and steady states, and homoclinic loss of periodic solution.
5 1. Introduction 1.1. Known models The predator-prey models, also known as producer-consumer models, describe the dynamics of biological systems in which two species interact, one a predator and one its prey. The simplest predator-prey model is the Lotka-Volterra model. This model was proposed in 1925 by the American biophysicist Alfred Lotka and the Italian mathematician Vito Volterra. The assumptions of this model are: density-independent exponential growth of the prey population in the absence of predators, a constant predation rate, a constant conversion rate and a constant per capita mortality rate of predators. Unfortunately, this model does not describe actual behavior observed in nature. One of the biggest problems is that the prey population is not self-limiting and, therefore, it can be arbitrary large. In order to fix this problem, a new version of the Lotka-Volterra model was introduced in This model uses a logistic growth instead of exponential growth for the prey. The logistic growth guarantees that the prey population is self-limiting. The prey can grow only to a certain saturation level. See for example (Pastor 2008) for more details. Our paper introduces a model related to the Rosenzweig- MacArthur model, which was proposed in 1963 in American Naturalist. See (Rosenzweig & MacArthur, 1963). The difference between the Lotka-Volterra model with logistic growth and the Rosenzweig-MacArthur model is that the predation rate is no longer assumed to be proportional to prey density Functional response There are many examples in nature which demonstrate that predators control the numbers of their prey. When predators are faced with increasing local density of their prey, they often respond by changing their consumption rate. This relationship of an individual predator's rate of food consumption to prey density was termed the "functional response" by C. S. Holling in He analyzed the act of predation and broke it onto behavioral units termed the components of 1
6 predation. There are many factors that might influence predation but we will consider only four major components: search, capture, handling and digestion. Theory and laboratory testing suggest that the functional response of an animal may take one of three forms called the Type I, Type II and Type III functional responses. See (Real, 1977), (Kuang, 2007), (Liu & Lou, 2010) or (Cosner, DeAngelis, Ault, & Olson, 1999). Figure 1: Three types of functional response Type I is the simplest functional response in which prey consumption rises linearly with prey density to a threshold level. This is called a filter feeder's functional response. However, increasing prey density beyond a threshold does not result in an increased per capita predation rate. This is due to digestion and handling components of predation. Figure 2: Type I functional response 2
7 The Type II functional response is a cyrtoid functional response in which the attack rate increases at a decreasing rate with prey density until it becomes constant at satiation. Cyrtoid behavioral responses are typical of predators that specialize on one or a few prey. Figure 3: Type II functional response In the Type III response, prey consumption remains low until a threshold density is reached. The predation rate then increases exponentially until levels out. The shape of this curve is described as sigmoidal. The predator ignores the prey when densities are very low in order to make hunting energetically profitable. Figure 4: Type III functional response 3
8 1.3. Rosenzweig-MacArthur model The Rosenzweig-MacArthur model uses the Holling II type functional response. See (Real, 1977). It belongs to the category of models with Predator Food Handling Times. After killing a prey, predator eats its captured food. This process consists of eating and digesting the meal portion. This fact was not considered in previous models (Lotka-Volterra) in which the predation rate was linearly proportional to prey density. The Rosenzweig-MacArthur model is inspired by behavior that can be found in nature. We cannot assume that the predator eats its prey at a constant rate. Some complications might occur, for example not having enough prey to hunt, delay in digesting, etc. The functional response in general is where interval, and is time interval necessary for a predator to catch its prey, also called the encounter is a time interval for chewing and digesting the captured prey. The Holling II type functional response (predation rate) is After using a substitution and, the functional response has the form 4
9 where is called the predation half-saturation constant and it is the density of prey at which the predator's killing rate reaches half its maximum value,. The aggregated model is where is a birth/death rate, is the carrying capacity, is the maximum per capita killing rate, is the value of that gives half the per capita killing rate (half-saturation constant), is the biomass conversion efficiency term that specify how much of the biomass of a killed individual is transformed into predator biomass. We assume that predators die with random death rate. Nullclines Nullclines are curves where the rate of change of the equations are 0, that is where or. The equilibria lie on intersection of these curves. There are two sets of nullclines. Predator nullclines: Prey nullclines: 5
10 The x-axis and y-axis are trivial nullclines. Non-trivial nullclines are a parabola and a vertical line. Figure 5: Rosenzweig-MacArthur model - nullclines Equilibria Three equilibria points can be found. Case I : This equilibrium describes the extinction of both species. Neither of the two species is alive. Case II : This equilibrium is called the monoculture equilibrium because only the prey population is alive, while predator population is extinct. This equilibrium is stable if, so in the case where predator nullcline is to the left of this equilibrium. It is unstable if, so when predator nullcline intersects the vertical (prey) nullcline. This behavior is called a transcritical bifurcation. Case III : Finally, this point is a coexistence equilibrium because both species are alive. The stability of this point also depends on values of. If the nullcline then lies to the right 6
11 of the maximum of the parabolic nullcline and the coexistence equilibrium is stable. Trajectories spiral into the coexistence equilibrium. If, the equilibrium is neutrally stable. If, the nullcline lies on the left of the maximum of the parabolic nullcline. The coexistence equilibrium is now unstable and trajectories spiral away from it. It can be shown by the Kolmogorv's theorem that the model contains either a stable equilibrium or a stable limit cycle, see (Pastor, 2008). It can be also shown by the Poincaré-Bendixson Theorem that if, the coexistence equilibrium is surrounded by a stable limit cycle. See (Hirsch, Smale, & Devaney, 2004), ( Kot, 2001) or (Rosen 1970 ). The change from stable to unstable coexistence equilibrium is called a Hopf bifurcation. 2. Model Development As mentioned before, the Rosenzweig-MacArthur model has a form This can be written more generally as where is predation term per capita (functional response). This term is determined by time intervals to catch and to process the prey 7
12 In the Rosenzweig-MacArthur model this term was The encounter rate is assumed to be proportional to prey density, therefore in this model This is the fact that we now focus on.. Also the handling time is assumed to be constant, so. Is it really always correct to assume that the time necessary to process food (chew, digest,...) is always constant in all situations that can occur in nature? Let us consider the situation where there is a pack of predators who hunts in a group. This behavior is very typical for many species, for example wolves, lions, hyenas, dogs, etc.. If there are more predators sharing a captured prey, the time to process will depend on the amount of food that each predator will get. Picture a scenario of wolves and a deer. Imagine a pack of 5 wolves hunting and eating a deer. After they finish eating their portions and resting they are ready to hunt again. Now imagine that there were 10 wolves instead. In that case each of them would get a smaller portion and it would be possible them to be still hungry. In that case they would be ready to hunt right away. Hunting in a group of individuals is a crucial factor that influence how soon predators are ready to hunt again. A slightly different situation with the same idea are rabbits (or cows) eating grass. The number of rabbits eating the grass in the same place influence the time they spend there and how long will it take them to move to a new area. Inspired by these natural situations, we introduce a new predation term where time interval to catch the prey stays the same, so we assume that the amount of time spend foraging does not depend on the size of the group, so 8
13 Therefore, the time to process food is not assumed constant, but it is assumed to be a saturating function of meal portion size. It can be noticed that when increases than the time to process decreases, and vice versa. With using this predation rate, the model we will study is where : is the density of prey(producer) population is the density of predator(consumer) population is the initial birth rate/death rate of prey is the carrying capacity for prey is the encounter rate is proportional to prey density is the maximum amount of food that each predator can eat is the half-saturation constant for the food processing time interval is the conversion efficiency constant is the death rate of predators. Predators are assumed to die for natural causes only. Note that if, then the model reduces to the Lotka-Volterra model. Note also that if then the model reduces to the Rosenzweig-MacArthur model. 9
14 3. Rescaling The new model includes 7 parameters. It is more efficient to decrease the number of parameters in order to analyze the model. The following rescaling was used. Let and. Then, and,. After substitution the system reads It is possible to simplify the model as follows After multiplying the 1st equation by and the 2nd equation by we obtain 10
15 We let, and, then the system becomes The number of parameters has been decreased to 5 but by change of time scale we can take and the number of parameters is reduced to 4. The parameter occurs only in new parameter. Therefore if we want to observe behavior of the system by increasing enrichment (increasing ), we need to vary only this one parameter. If we want to increase half-saturation to examine the effects of meal size, we will vary only the new parameter Similarly, varying predator death rate,, effects only the new parameter. 4. Nullcline analysis As usual, nullclines have been analyzed. Recall that nullclines are curves where a state variable derivative is zero. Hence, these are curves where or. Since is a common factor in the first equation of, it does not affect the nullclines for the system. In our case 11
16 There are 2 prey nullclines and 2 predator nullclines :. The x-axis and y-axis are zero nullclines. The non-zero nullclines are different than before. The predator nullcline is no longer a vertical straight line and the prey nullcline is no longer a parabola. The shapes of these curves are more complicated Predator nullcline ( ) Notice that the predator nullcline is dependent only on the two parameters and. By solving the equation for we obtain This curve is not a straight line anymore. It is not obvious what shapes can be expected. 12
17 Asymptotes: Since we assumed being positive, the curve has a vertical asymptote at. As the curve is asymptotic to a straight line with slope. x and y-intercepts: So the curve always passes through the origin However, if, so we are not interested in this branch. So, the predator nullcline intercepts the x-axis at the point. We are interested only in positive intercepts. Since we assumed all parameters are positive, the only way how this intercept can be positive is when, both conditions are satisfied when. We draw predator nullclines for all three cases; for and. 13
18 1) Figure 6: Predator nullcline, H=0.8, D=0.5 We can see that the curve crosses the x-axis at positive value of. Figure 7: Predator Nullcline, H=0.8, D=1 In this case, there is no intercept with x- axis. By increasing asymptotic behavior as. we lost it. We can observe The prey nullcline has a horizontal asymptote at 14
19 3) Figure 8: Predator Nullcline, H=0.8, D=1.3 If we increase further, the nullcline bends up. This nullcline has absolutely different behavior then predator nullclline in the Rosenzweig- MacArthur model. For small values of it changes. it is almost a vertical straight line but with increasing 4.2 Prey nullcline ( ) By solving the equation 15
20 for was obtained This nullcline is too complicated to be analyzed. It is not monotonically increasing or decreasing and is not just concave down or concave up. The behavior is more complicated. This nullcline depends only on the two parameters and. x-intercept This curve intercepts the x-axis at. y-intercept This curve intercepts the y-axis at. The following behavior can be observed when and varies. 16
21 K=0.5 K=0.8 K=1 K=1.5 K=2.5 K=5 Figure 9: Behavior of the prey nullcline for different values of K This nullcline is obviously not parabolic any more Coexistence equilibria We plotted the nullclines for and fixed and being varied. We observed several situations and saw that there are two dividing cases that determine number of coexistence equilibria. Case I is when two non-zero nullclines have a first tangent contact at one point. This is the situation when there is a creation of coexistence equilibria that splits then into 2 equilibria. This situation will corresponds to a saddle node bifurcation of coexistence equilibria. Case II occurs when there are two coexistence equilibria and one of them will be lost at. There are two nullcline intercepts and one of them is getting closer to x-axis with increasing and then moves 17
22 to the 4th quadrant. This corresponds to a transcritical bifurcation. Before this coexistence equilibrium is at the x-axis, the monoculture equilibrium is stable, but once the coexistenceequilibrium crosses the x-axis it becomes unstable. This will lead to a transcritical bifurcation Saddle node bifurcation curve We considered first the situation when two nullclines get closer to each other until they have one common tangent point. We wanted to find a set of parameters and for this case. This situation occurs when the gradient of one nullcline is a non-zero multiple of the gradient of the second one. If we denote non-zero nullclines equations as functions and then the necessary condition is. This condition can be replaced by setting, where is the Jacobian matrix The set of conditions that will yield one contact point saddle node bifurcation is By using Mathematica we got the following results: 18
23 So the set of all conditions is The solution of this system of 3 equations can be found in Appendix A since it is too complicated to state it inside the paper. We fixed parameter and plotted values of parameter as a function of Here are pictures for different values of On the horizontal axis are values of parameter and on the vertical axis are values of parameter 19
24 H=0.1 H=0.5 H=0.8 H=1 H=1.5 H=2. Figure 10: Saddle node bifurcation curve On the left side of the curve there are no coexistence equilibrium and on the right side there are 2 of them. Among the line, there is one coexistence equilibrium Transcritical bifurcation A second option was when the nullclines intersect at the x-axis. This is a crutial situation for the monoculture equilibrium because it determines if it is stable or unstable. It can be a transition between 0 and 1 coexistence equilibria, if the nullclines meet for first time at x-axis or between 2 and 1 equilibrium, as is in the case if we lose one equilibrium because it moved into the 4th quadrant. 20
25 The set of conditions for this situation is Notice that this curve does not depend on values of This curve always passes through the origin. and it has horizontal asymptote at. Figure 11: Transcritical bifurcation (intercept at the x-axis) 21
26 We combined both pictures of both cases to see how the whole 1st quadrant is divided into regions with different number of equilibria. Figure 12: Saddle Node and Transcritical Bifurcation Curves We found a common point of these two curves, points where nullclines have a first (tangential) contact at the x-axis. This is the situation that divides 0 and 1 coexistence equilibrium. 22
27 Solution is A common point of these two curves always exists since we assumed. We added this point into our picture. Figure 13: Part of the curve without biological meaning for H=0.5 We observed that the dashed part of the curve corresponds to the situation where there is a tangent contact for negative values of and. Since we are interested in positive values, we looked only at the segment of the curve for. 23
28 In summary, we show these bifurcation curves for different values of. H=0.2 H=0.5 H=0.8 H=1 H=1.5 H=2.5 Figure 14: Parameter region division Observe that common point is moving to the right side for getting closer to and moving to the left side for getting larger. If, this point is in first quadrant. If, it is located at the origin and if it is in the third quadrant. The parameter space is divided into 3 regions with different numbers of coexistence equilibria. 24
29 5. Equilibria The number of equilibria are summarized as 0 saddle node 2 transcritical 1 Figure 15: Parameter space picture with number of equilibria for H=0.8 Solving the systems of two equations for and gives equilibria points Extinctiction This is the case when both species have died out. 25
30 5.2. Monoculture equilibrium This is where only the prey population is alive, while all predators have died out. In order to determined stability of this equilibrium, we computed the determinant of the Jacobian matrix, which relys on the partial derivatives of equations, when and. Let and. Then This value is not new for us because this value of is the same as value of for which an intercept of two nonzero nullclines is at the x-axis. If, so if, the equilibrium is stable and if the equilibrium is unstable. So as long as lies to the right side of the prey nullcline,the monoculture equilibrium is stable. When the prey nullcline 26
31 crosses, we either create one stable coexistence equilibrium or loose one of two coexistence equilibria and monoculture equilibrium looses it stability at a transcritical bifurcation 5.3. Coexistence equilibria There appeared to be at most two coexistence equilibria.they depend on parameters. and Since these equilibria are difficult to analyze, we used Mathematica to show what behaviors and which bifurcations is possible to observe for the system (7). Without loss of generality we let. For simplicity we set. We will fixed and vary only and. Let us look first at different situations for and. 27
32 Case I: From Figure 16 it is obvious that we can observe 4 different scenarios. When. 0 saddle node 2 transcritical 1 Figure 16: Set of parameters for H=0.8 1) We chose and consider the phase planes for system (7) for a sequence of increasing K values. 28
33 K=0.1 K=0.11 K=0.16 K=1 K=2 K=2.75 K=3.3 K=3.345 K=3.45 Figure 17: Behavior of the model for H=0.8, D=0.1 For small values of there is no coexistence equilibrium. There is only monoculture equilibrium that is stable. Then nullclines intercept at one contact point on the x-axis. After increasing the intersection point is moves deeper into the first quadrant and represents a coexistence equilibrium that is stable. The monoculture equilibrium becomes unstable and we can observe a transcritical bifurcation. With even larger we can observe creation of a periodic solution and a lost of stability of the coexistence equilibrium at a Hopf bifurcation. The trajectory spirals out 29
34 of the equilibrium into a limit cycle. This state appears to hold for all larger values of behavior is typical of the Rosenzweig-MacArthur system.. This 2) K=0.15 K= K=0.4 K=1.5 K=2.25 K=2.5 K=2.515 K=2.55 K=2.7 Figure 18: Behavior for H=0.8, D=0.2 We can observe the same behavior for this set of parameters as for. The only difference is that the transcritical and saddle node bifurcation occur at the same value of. 30
35 3) K=0.2 K=2.5 K=2.96 K=3.3 K=3.5 K=3.7 K=4 K=4.3 K=5 As before we start with no inner coexistence equilibria; there is only a stable monoculture equilibrium. After increasing, nullclines are getting closer and eventually have a contact point inside the first quadrant. This is a creation of a coexistence equilibrium. With larger this equilibrium splits into two points. The upper one is stable (attractor) and the lower one is unstable (saddle). We have observed the saddle node bifurcation. As before, monoculture equilibrium is stable for low values of. When there are two inner equilibria, the monoculture 31
36 equilibrium is still stable. With increasing, the lower equilibrium approaches the x-axis and then collides with the monoculture equilibrium It then moves to the 4th quadrant and become biologically irrelevant. After losing this second equilibrium, the monoculture equilibrium becomes unstable and we can observe again transcritical bifurcation. The upper equilibrium stays stable for low values of. However, with increasing, we can observe the creation of a periodic solution again by way of a Hopf bifurcation. The periodic solution appears to persist for all larger values of. 4) K=0.1 K=4 K=10 K=21 K=22 K=30 In this case, we can observe a slightly different situation. For small values of, there is again only stable monoculture equilibrium and no coexistence equilibrium. With larger, the prey nullcline approaches the predator nullcline until they have a contact point. Enriching the system even more, this inner equilibirum splits again into two coexistence equilibria but the 32
37 monoculture equilibria stays stable for all values of K. One equilibrium is a saddle and the other one is a unstable repellor. We do not observe periodic solutions, therefore no Hopf bifurcation occurs. This characteristic of having two unstable coexistence equilibria persists for all large values of. This behavior is not ecologically desirable because the predator will always die out. This means that if is too large (greater than 1), there is no way how to prevent the extinction of predators, and only the prey population will stay alive. Notice how large values of must be in order to observe two coexistence equilibria for, compare to previous cases. Case II: Note that for the model does not support the case where there is only zero or one equilibrium as we saw for smaller values. We looked at three different cases. and 0 saddle node 2 transcritical 1 Figure 19: Three different scenarios for H=1.5 33
38 1) K=0.4 K=0.4 2 K=0.415 K=0.43 K=1 K=2.97 K=3 K=2.7 K=4 Figure 20: H=1.2, D=0.3 In this enrichment sequence, we can observe the same behavior as in case. There is no coexistence equilibrium for small, only a stable monoculture equilibrium. After nullclines have a contact point inside the first quadrant there is one coexistence equilibrium that splits into two. One of them is a stable attractor and second one is a saddle node. The saddle node will collide with the monoculture equilibrium and switch stability with the monoculture equilibrium, which becomes unstable and at a ranscritical bifurcation. One 34
39 equilibrium moves into biologically irrelevant 4th quadrant. With increasing a Hopf bifurcation again., we can observe The only difference in behavior for is that we will always observe first the creation of a coexistence equilibrium inside the first quadrant (saddle node bifurcation) and after that a transcritical bifurcation. There is no set of parameters of and that would give us first the creation of a coexistence equilibrium at the x-axis (transcritical bifurcation). 2) K=0.1 K=2 K=2.33 K=2.5 K=2.8 K=3.14 K=3.4 K=4 K=5 Figure 21: H=1.2, D=0.8 35
40 3) K=0.5 K=5 K=8 K=9.08 K=10 K=15 Figure 22: H=1.2, D=1.2 In this enrichment sequnce a pair of unstable coexistence equilibria are created, but the monoculture equilibrium will always be attracting. Notice that large values of for obtaining a coexistence equilibria when. are necessary Stability region Simulations suggest that the behavior of model solutions for smaller values of depends on the initial conditions for the simulation. We looked at this situatuation more carefuly and the following behavior was observed. 36
41 Figure 23: Dependence relationship on initial values for H=0.8, D=0.1 and K=3.34 These two pictures are for the same set of parameters and. They differ only in choice of inital conditons. We can see that the behavior depends on the initial condition. For one choice the trajectory spirals towards the attractive equilibrium and for another choice trajectories are spiral out into a limit cycle. This demonstrates that for certain parameter values the model can support "bistability" ( a locally attracting coexistence equilibrium surrounded by an attracting periodic orbit.) The same behavior was observed for different values of parameters. Figure 24: Dependence of behavior on initial conditions for H=0.8, D=0.2, K=
42 Figure 25: H=1.2, D=0.3, K= Hopf Bifurcation We already know values where transcritical bifurcation and saddle node bifurcation occur. We also want to know for which set of parameters we can have a Hopf bifurcation. A Hopf bifurcation occurs when an equilibrium looses its stability at the creation of a periodic solution. In order to have Hopf bifurcation the eigenvalues of the system's Jacobian matrix must be purely imaginary. If is eigenvalue of the Jacobian matrix and its real part is 0 then. In order to find such eigenvalues we want, where is the identity matrix. Let and be equations of non-zero nullclines. The determinant condition becomes 38
43 This complex number is zero if the real part and the imaginary parts are both equal to zero. Therefore, we have the conditions Since we still work with equilibria we need to have simultaneous solution of our nullclines. The complete set of conditions that must be satisfied in order to have a Hopf bifurcation is Since and 39
44 we can write this as In our case Simplifications of this system and its solution can be found in Appendix C. If we fix this is a system of 4 equations for 5 unknowns and. We decided to fix also and solve this system for and. Even Mathematica had a hard time to solve this system exactly with the Solve command. It took it about 35 minutes to find a solution which was too complicated to be useful. Therefore, we tried to find a numerical solution and compared it to the exact solutions to check that they were the same. Since the NSolve commmand was the faster way how to solve the system, we choose this approach. A typical calculation took about 2 seconds to get a solution. Once we found a solution we tried different values of and obtained many combinations of values of associated and. Some of them were negative or complex. We focused on positive real ones since they were biologically relevant. After plotting these points it was obvious what shape of curve we might typically expect. We also chose several points and used the Spline commend in Mathematica to obtain a curve where the Hopf bifurcation occurs. This curve was added to Saddle node and transcritical bifurcation curves and the final region divisions for and were obtained. 40
45 H=0.8 H=5 Figure 26: Hopf bifurcation, single points Parameter space pictures that covers all types of bifucartion are observed. 0 saddle node 2 Homoclinic transcritical 1 Hopf Figure 27: Parameter space region division for H=5 Points coordinates : 41
46 0 saddle node 2 transcritical 1 Hopf Figure 28: Parameter space region division for H=0.8 Points coordinates: In summary, the behavior of the model is more complicated then behavior of Rosenzweig- MacArthur model. There are several cases that can occur when. Case I: The creation of a coexistence equilibrium at a transcritical bifurcation followed by a Hopf bifurcation. Case II: The creation of a coexistence equilibrium at a saddle node bifurcation and a transcritical bifurcation simultaneously followed by a Hopf bifurcation. Case III: The creation of two coexistence equilibria at a saddle node bifurcation followed by the loss of one coexistence equilibria at a transcritical bifurcation, and a Hopf bifurcation. Case IV: The creation of two coexistence equilibria at a saddle node bifurcation followed by the loss of one coexistence equilibria at a transcritical and a Hopf bifurcation simultaneously. Case V: The creation of two coexistence equilibria at a saddle node bifurcation followed by a Hopf bifurcation and the loss of one coexistence equlibrium at a transcritical bifurcation. Case VI: The creation of two coexistence equilibria at a saddle node bifurcation. 42
47 5.5. Solution We verified observed behavior by plotting solutions. On the horizontal axis is time and on the vertical axis are the prey density and the predator density. We chose. Figure 29: Solution for H=0.8, D=0.2, K=0.1 For small values of out., monoculture equilibrium is stable and the predator population has died Figure 30: Solution for H=0.8, D=0.2, K=1.5 For larger values, there is a stable coexistence equilibrium. Trajectories spiral in. 43
48 Figure 31: Solution for H=0.8, D=0.2, K=2.25 For even larger values, oscilliations are getting bigger. Figure 32: Solution for H=0.8, D=0.2, K=2.55 As we stated before, it is possible to observe bistability. For the same parameter values but different initial values, we observe either oscillations getting smaller or oscillations getting smaller which corresponds to spiraling in and spiraling out in the phrase plane. 44
49 Homoclinic bifurcation We could also observe the case of homoclinic bifurcation. This case was not fully studied in this paper. Here are pictures for values that shows that this situation can happen in our model. K=3.2 K=3.5 K=4.5 K=4.9 K=5.5 Figure 33: Homoclinic bifurcation, H=5, D= Comparison to Rosenzweig-MacArthur model Nullclines In the Rosenzweig-MacArthur model there are 4 nullclines, the x-axis, the y-axis,a straight line predator nullcline and a parabolic prey nullcline. In our model, there are also 4 nullcllines. The x 45
50 and y axis are the same, but the others differs significantly. Both of them are more complicated curves with behavior depending on model parameter values. Equilibria In both models there is an extinction equilibrium and monoculture equilibrium. We should not forget that used in our model is a rescaled parameter while in Rosenzweig- MacArthur it is the original not rescaled value. The models differ significantly in their coexistence equilibriua. While in the Rosenzweig-MacArthur model we can observe only one coexistence equilibrium, in our model under certain conditions can observe two coexistence equilibria. Bifurcation In the model developed here it is possible to observe a transcritical bifurcation as in the Rosenzweig-MacArthur model. We can also see a Hopf bifurcation as in the Rosenzweig- MacArthur model. However, the behavior of our model strongly depends on the parameter values for and. Our model shows richer behavior including bistability and Homoclinic bifurcations in contrast to the Rosenzweig-MacArthur model. 7. Further Research Suggestions 1) The influence on different values of values of and could be analyzed. Recall we assumed that and. 2) For which values of parameters and is possible to observe a Homoclinic bifurcation? 3) What is the dependence relationship between initial values and stability of a coexistence equilibrium? For which values of and is it possible to see such a behavior? 46
51 Table of Figures Figure 1: Three types of functional response... 2 Figure 2: Type I functional response... 2 Figure 3: Type II functional response... 3 Figure 4: Type III functional response... 3 Figure 5: Rosenzweig-MacArthur model - nullclines... 6 Figure 6: Predator nullcline, H=0.8, D= Figure 7: Predator Nullcline, H=0.8, D= Figure 8: Predator Nullcline, H=0.8, D= Figure 9: Behavior of the prey nullcline for different values of K Figure 10: Saddle node bifurcation curve Figure 11: Transcritical bifurcation (intercept at the x-axis) Figure 12: Saddle Node and Transcritical Bifurcation Curves Figure 13: Part of the curve without biological meaning for H= Figure 14: Parameter region division Figure 15: Parameter space picture with number of equilibria for H= Figure 16: Set of parameters for H= Figure 17: Behavior of the model for H=0.8, D= Figure 18: Behavior for H=0.8, D= Figure 19: Three different scenarios for H= Figure 20: H=1.2, D= Figure 21: H=1.2, D= Figure 22: H=1.2, D= Figure 23: Dependence relationship on initial values for H=0.8, D=0.1 and K= Figure 24: Dependence of behavior on initial conditions for H=0.8, D=0.2, K= Figure 25: H=1.2, D=0.3, K= Figure 26: Hopf bifurcation, single points Figure 27: Parameter space region division for H= Figure 28: Parameter space region division for H= Figure 29: Solution for H=0.8, D=0.2, K= Figure 30: Solution for H=0.8, D=0.2, K= Figure 31: Solution for H=0.8, D=0.2, K= Figure 32: Solution for H=0.8, D=0.2, K= Figure 33: Homoclinic bifurcation, H=5, D=
52 References Cosner, C., DeAngelis, D. L., Ault, J. S., & Olson, D. B. (1999). Effects of Spatial Grouping on the Functional Response of Predators. Theoretical Population Biology, Hirsch, M. W., Smale, S., & Devaney, R. L. (2004). Differential Equations, Dynamical Systems, an An Introduction to Chaos. London: Elsevies Academic Press. Kot, M. (2001). Elements of Mathematical Ecology. Cambridge, UK: Cambridge University Press. Kuang, Y. (2007). Some mechanistically derived functional responses. Mathematical Biosciences and Engineering. Liu, X., & Lou, Y. (2010). Global dynamics of a predator-prey model. Journal of Mathematical Analysis and Applications, Pastor, J. (2008). Mathematical Ecology of Populations and Ecosystems. John Wiley and Sons. Real, L. A. (1977, March-April). The kinetics of functional response. The American Naturalist, pp Rosen, R. (1970). Dynamical Systems Theory in Biology. New York: John Wiley and Sons. Rosenzweig, M., & MacArthur, R. H. (1963). Graphical representation and stability conditions of predator-prey interactions. The American Naturalist 97,
A Producer-Consumer Model With Stoichiometry
A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More informationEnrichment in a Producer-Consumer Model with varying rates of Stoichiometric Elimination
Enrichment in a Producer-Consumer Model with varying rates of Stoichiometric Elimination Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationModels Involving Interactions between Predator and Prey Populations
Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate
More informationMATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation
MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation 1 Bifurcations in Multiple Dimensions When we were considering one-dimensional systems, we saw that subtle changes in parameter
More informationEnrichment in a General Class of Stoichiometric Producer-Consumer Population Growth Models
Enrichment in a General Class of Stoichiometric Producer-Consumer Population Growth Models Harlan Stech a, c, Bruce Peckham a and John Pastor b a Department of Mathematics and Statistics, University of
More informationIs chaos possible in 1d? - yes - no - I don t know. What is the long term behavior for the following system if x(0) = π/2?
Is chaos possible in 1d? - yes - no - I don t know What is the long term behavior for the following system if x(0) = π/2? In the insect outbreak problem, what kind of bifurcation occurs at fixed value
More informationLab 5: Nonlinear Systems
Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting
More informationTHE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL
THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More information3.5 Competition Models: Principle of Competitive Exclusion
94 3. Models for Interacting Populations different dimensional parameter changes. For example, doubling the carrying capacity K is exactly equivalent to halving the predator response parameter D. The dimensionless
More informationAN EXTENDED ROSENZWEIG-MACARTHUR MODEL OF A TRITROPHIC FOOD CHAIN. Nicole Rocco
AN EXTENDED ROSENZWEIG-MACARTHUR MODEL OF A TRITROPHIC FOOD CHAIN Nicole Rocco A Thesis Submitted to the University of North Carolina Wilmington in Partial Fulfillment of the Requirements for the Degree
More informationx 2 F 1 = 0 K 2 v 2 E 1 E 2 F 2 = 0 v 1 K 1 x 1
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 20, Number 4, Fall 1990 ON THE STABILITY OF ONE-PREDATOR TWO-PREY SYSTEMS M. FARKAS 1. Introduction. The MacArthur-Rosenzweig graphical criterion" of stability
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationLecture 20/Lab 21: Systems of Nonlinear ODEs
Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More information= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
More informationDynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More informationDIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS
DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University
More informationM.5 Modeling the Effect of Functional Responses
M.5 Modeling the Effect of Functional Responses The functional response is referred to the predation rate as a function of the number of prey per predator. It is recognized that as the number of prey increases,
More informationEvolutionary suicide in a two-prey-one-predator model with Holling type II functional response. Caterina Vitale
Evolutionary suicide in a two-prey-one-predator model with Holling type II functional response Caterina Vitale Department of Mathematics and Statistics Faculty of Science University of Helsinki May 2016
More informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationGerardo Zavala. Math 388. Predator-Prey Models
Gerardo Zavala Math 388 Predator-Prey Models Spring 2013 1 History In the 1920s A. J. Lotka developed a mathematical model for the interaction between two species. The mathematician Vito Volterra worked
More informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationMath 1280 Notes 4 Last section revised, 1/31, 9:30 pm.
1 competing species Math 1280 Notes 4 Last section revised, 1/31, 9:30 pm. This section and the next deal with the subject of population biology. You will already have seen examples of this. Most calculus
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationThe Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor
, pp. 35-46 http://dx.doi.org/10.14257/ijbsbt.2017.9.3.04 The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor Alemu Geleta Wedajo
More informationDIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS
DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University
More informationLocal Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid)
Biology Medicine & Natural Product Chemistry ISSN: 089-6514 Volume 5 Number 1 016 Pages: 9-14 DOI: 10.1441/biomedich.016.51.9-14 Local Stability Analysis of a Mathematical Model of the Interaction of Two
More informationSTABILITY OF EIGENVALUES FOR PREDATOR- PREY RELATIONSHIPS
Research Article STABILITY OF EIGENVALUES FOR PREDATOR- PREY RELATIONSHIPS Tailor Ravi M., 2 Bhathawala P.H. Address for Correspondence Assistant professor of Laxmi Institute of Technology, Sarigam, Valsad
More informationBIOS 6150: Ecology Dr. Stephen Malcolm, Department of Biological Sciences
BIOS 6150: Ecology Dr. Stephen Malcolm, Department of Biological Sciences Week 7: Dynamics of Predation. Lecture summary: Categories of predation. Linked prey-predator cycles. Lotka-Volterra model. Density-dependence.
More information1. Population dynamics of rabbits and foxes
1. Population dynamics of rabbits and foxes (a) A simple Lotka Volterra Model We have discussed in detail the Lotka Volterra model for predator-prey relationships dn prey dt = +R prey,o N prey (t) γn prey
More informationDynamics Analysis of Anti-predator Model on Intermediate Predator With Ratio Dependent Functional Responses
Journal of Physics: Conference Series PAPER OPEN ACCESS Dynamics Analysis of Anti-predator Model on Intermediate Predator With Ratio Dependent Functional Responses To cite this article: D Savitri 2018
More informationv n+1 = v T + (v 0 - v T )exp(-[n +1]/ N )
Notes on Dynamical Systems (continued) 2. Maps The surprisingly complicated behavior of the physical pendulum, and many other physical systems as well, can be more readily understood by examining their
More informationQualitative Analysis of Tumor-Immune ODE System
of Tumor-Immune ODE System LG de Pillis and AE Radunskaya August 15, 2002 This work was supported in part by a grant from the WM Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1 Simplified System of
More informationTHETA-LOGISTIC PREDATOR PREY
THETA-LOGISTIC PREDATOR PREY What are the assumptions of this model? 1.) Functional responses are non-linear. Functional response refers to a change in the rate of exploitation of prey by an individual
More informationLogistic Model with Harvesting
Logistic Model with Harvesting Model Here, we will talk about a model for the growth and harvesting of a fish population. Suppose that, absent any fishermen, the population grows logistically. That is,
More informationStudent Exploration: Food Chain
Name: Date: Student Exploration: Food Chain Vocabulary: consumer, ecosystem, equilibrium, food chain,, predator, prey, producer Prior Knowledge Questions (Do these BEFORE using the Gizmo.) The Food Chain
More informationModeling Prey and Predator Populations
Modeling Prey and Predator Populations Alison Pool and Lydia Silva December 15, 2006 Abstract In this document, we will explore the modeling of two populations based on their relation to one another. Specifically
More informationPredator-prey interactions
Predator-prey interactions Key concepts ˆ Predator-prey cycles ˆ Phase portraits ˆ Stabilizing mechanisms ˆ Linear stability analysis ˆ Functional responses ˆ The paradox of enrichment Predator-prey cycles
More informationMath 232, Final Test, 20 March 2007
Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.
More informationNUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD CHAIN WITH LOTKA-VOLTERRA LINEAR FUNCTIONAL RESPONSE
Journal of Sustainability Science and Management Volume 6 Number 1, June 2011: 44-50 ISSN: 1823-8556 Universiti Malaysia Terengganu Publisher NUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD
More informationStable Coexistence of a Predator-Prey model with a Scavenger
ANTON DE KOM UNIVERSITEIT VAN SURINAME INSTITUTE FOR GRADUATE STUDIES AND RESEARCH Stable Coexistence of a Predator-Prey model with a Scavenger Thesis in partial fulfillment of the requirements for the
More informationPopulation Dynamics in a Microfounded Predator-Prey Model. Thomas Christiaans. University of Siegen. Discussion Paper No
Volkswirtschaftliche Diskussionsbeiträge V W L Population Dynamics in a Microfounded Predator-Prey Model Thomas Christiaans University of Siegen Discussion Paper No. 8-04 ISSN 433-058x UNIVERSITÄT SIEGEN
More information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationFunctional Response to Predators Holling type II, as a Function Refuge for Preys in Lotka-Volterra Model
Applied Mathematical Sciences, Vol. 9, 2015, no. 136, 6773-6781 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.53266 Functional Response to Predators Holling type II, as a Function Refuge
More informationBifurcation and Stability Analysis of a Prey-predator System with a Reserved Area
ISSN 746-733, England, UK World Journal of Modelling and Simulation Vol. 8 ( No. 4, pp. 85-9 Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area Debasis Mukherjee Department
More informationA Stability Analysis on Models of Cooperative and Competitive Species
Research Journal of Mathematical and Statistical Sciences ISSN 2320 6047 A Stability Analysis on Models of Cooperative and Competitive Species Abstract Gideon Kwadzo Gogovi 1, Justice Kwame Appati 1 and
More informationDynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs
Dynamics of a opulation Model Controlling the Spread of lague in rairie Dogs Catalin Georgescu The University of South Dakota Department of Mathematical Sciences 414 East Clark Street, Vermillion, SD USA
More information1 2 predators competing for 1 prey
1 2 predators competing for 1 prey I consider here the equations for two predator species competing for 1 prey species The equations of the system are H (t) = rh(1 H K ) a1hp1 1+a a 2HP 2 1T 1H 1 + a 2
More information2 One-dimensional models in discrete time
2 One-dimensional models in discrete time So far, we have assumed that demographic events happen continuously over time and can thus be written as rates. For many biological species with overlapping generations
More informationApplications in Biology
11 Applications in Biology In this chapter we make use of the techniques developed in the previous few chapters to examine some nonlinear systems that have been used as mathematical models for a variety
More informationNumerical bifurcation analysis of a tri-trophic food web with omnivory
Mathematical Biosciences 177&178 (2002) 201 228 www.elsevier.com/locate/mbs Numerical bifurcation analysis of a tri-trophic food web with omnivory B.W. Kooi *, L.D.J. Kuijper, M.P. Boer 1, S.A.L.M. Kooijman
More information8 Ecosystem stability
8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationPopulation Models Part I
Population Models Part I Marek Stastna Successoribus ad Successores Living things come in an incredible range of packages However from tiny cells to large mammals all living things live and die, making
More informationA NUMERICAL STUDY ON PREDATOR PREY MODEL
International Conference Mathematical and Computational Biology 2011 International Journal of Modern Physics: Conference Series Vol. 9 (2012) 347 353 World Scientific Publishing Company DOI: 10.1142/S2010194512005417
More informationThe Effect of a Smart Predator in a One Predator, Two Prey System
The Effect of a Smart Predator in a One Predator, Two Prey System Elizabeth Green October 6, 4 Abstract This paper analyzes a food web with a predator and two non-competing preys where the predation follows
More informationModeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool
Modeling the Immune System W9 Ordinary Differential Equations as Macroscopic Modeling Tool 1 Lecture Notes for ODE Models We use the lecture notes Theoretical Fysiology 2006 by Rob de Boer, U. Utrecht
More informationWorkshop on Theoretical Ecology and Global Change March 2009
2022-3 Workshop on Theoretical Ecology and Global Change 2-18 March 2009 Stability Analysis of Food Webs: An Introduction to Local Stability of Dynamical Systems S. Allesina National Center for Ecological
More informationInteractions between predators and prey
Interactions between predators and prey What is a predator? Predator An organism that consumes other organisms and inevitably kills them. Predators attack and kill many different prey individuals over
More informationEcology 302: Lecture VII. Species Interactions.
Ecology 302: Lecture VII. Species Interactions. (Gotelli, Chapters 6; Ricklefs, Chapter 14-15) MacArthur s warblers. Variation in feeding behavior allows morphologically similar species of the genus Dendroica
More informationGrowth models for cells in the chemostat
Growth models for cells in the chemostat V. Lemesle, J-L. Gouzé COMORE Project, INRIA Sophia Antipolis BP93, 06902 Sophia Antipolis, FRANCE Valerie.Lemesle, Jean-Luc.Gouze@sophia.inria.fr Abstract A chemostat
More informationSEMELPAROUS PERIODICAL INSECTS
SEMELPROUS PERIODICL INSECTS BRENDN FRY Honors Thesis Spring 2008 Chapter Introduction. History Periodical species are those whose life cycle has a fixed length of n years, and whose adults do not appear
More informationBiological control does not imply paradox
Mathematical Biosciences 208 (2007) 26 32 www.elsevier.com/locate/mbs Biological control does not imply paradox Bo Deng a, *, Shannon Jessie b, Glenn Ledder a, Alex Rand c, Sarah Srodulski d a Department
More informationESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp DOI: /m2an:
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp. 339 344 DOI: 10.1051/m2an:2003029 PERSISTENCE AND BIFURCATION ANALYSIS
More informationDynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part II: Biology Applications Lecture 6: Population dynamics Ilya Potapov Mathematics Department, TUT Room TD325 Living things are dynamical systems Dynamical systems theory
More informationdv dt Predator-Prey Models
Predator-Prey Models This is a diverse area that includes general models of consumption: Granivores eating seeds Parasitoids Parasite-host interactions Lotka-Voterra model prey and predator: V = victim
More informationPredation. Vine snake eating a young iguana, Panama. Vertebrate predators: lions and jaguars
Predation Vine snake eating a young iguana, Panama Vertebrate predators: lions and jaguars 1 Most predators are insects Parasitoids lay eggs in their hosts, and the larvae consume the host from the inside,
More informationMA 138: Calculus II for the Life Sciences
MA 138: Calculus II for the Life Sciences David Murrugarra Department of Mathematics, University of Kentucky. Spring 2016 David Murrugarra (University of Kentucky) MA 138: Section 11.4.2 Spring 2016 1
More informationAge (x) nx lx. Population dynamics Population size through time should be predictable N t+1 = N t + B + I - D - E
Population dynamics Population size through time should be predictable N t+1 = N t + B + I - D - E Time 1 N = 100 20 births 25 deaths 10 immigrants 15 emmigrants Time 2 100 + 20 +10 25 15 = 90 Life History
More informationDYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD
Journal of Sustainability Science Management Volume 10 Number 2, December 2015: 16-23 ISSN: 1823-8556 Penerbit UMT DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD SARKER MD SOHEL
More informationContinuous time population models
Continuous time population models Jaap van der Meer jaap.van.der.meer@nioz.nl Abstract Many simple theoretical population models in continuous time relate the rate of change of the size of two populations
More informationPredation. Predation & Herbivory. Lotka-Volterra. Predation rate. Total rate of predation. Predator population 10/23/2013. Review types of predation
Predation & Herbivory Chapter 14 Predation Review types of predation Carnivory Parasitism Parasitoidism Cannabalism Lotka-Volterra Predators control prey populations and prey control predator populations
More informationStability Analysis And Maximum Profit Of Logistic Population Model With Time Delay And Constant Effort Of Harvesting
Jurnal Vol 3, Matematika, No, 9-8, Juli Statistika 006 & Komputasi Vol No Juli 006 9-8, Juli 006 9 Stability Analysis And Maximum Profit Of Logistic Population Model With Time Delay And Constant Effort
More informationInteractions of Living Things
CHAPTER 14 LESSON 2 Interactions of Living Things Populations and Communities Key Concepts How do individuals and groups of organisms interact? What are some examples of symbiotic relationships? What do
More informationChapter 7. Nonlinear Systems. 7.1 Introduction
Nonlinear Systems Chapter 7 The scientist does not study nature because it is useful; he studies it because he delights in it, and he delights in it because it is beautiful. - Jules Henri Poincaré (1854-1912)
More informationNonlinear Dynamics. Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna.
Nonlinear Dynamics Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna http://www.moreno.marzolla.name/ 2 Introduction: Dynamics of Simple Maps 3 Dynamical systems A dynamical
More informationBIOS 3010: ECOLOGY. Dr Stephen Malcolm. Laboratory 6: Lotka-Volterra, the logistic. equation & Isle Royale
BIOS 3010: ECOLOGY Dr Stephen Malcolm Laboratory 6: Lotka-Volterra, the logistic equation & Isle Royale This is a computer-based activity using Populus software (P), followed by EcoBeaker analyses of moose
More information8 Autonomous first order ODE
8 Autonomous first order ODE There are different ways to approach differential equations. Prior to this lecture we mostly dealt with analytical methods, i.e., with methods that require a formula as a final
More informationLokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt
Periodic solutions A periodic solution is a solution (x(t), y(t)) of dx = f(x, y) dt dy = g(x, y) dt such that x(t + T ) = x(t) and y(t + T ) = y(t) for any t, where T is a fixed number which is a period
More informationLong Term Dynamics for Two Three-Species Food Webs
Long Term Dynamics for Two Three-Species Food Webs B. Bockleman, B. Deng, E. Green, G. Hines, L. Lippett, J. Sherman September 2, 23 Abstract In this paper, we analyze two possible scenarios for food webs
More informationGetting Started With The Predator - Prey Model: Nullclines
Getting Started With The Predator - Prey Model: Nullclines James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 28, 2013 Outline The Predator
More informationDynamical Systems: Ecological Modeling
Dynamical Systems: Ecological Modeling G Söderbacka Abstract Ecological modeling is becoming increasingly more important for modern engineers. The mathematical language of dynamical systems has been applied
More informationQualitative Analysis of Tumor-Immune ODE System
of Tumor-Immune ODE System L.G. de Pillis and A.E. Radunskaya August 15, 2002 This work was supported in part by a grant from the W.M. Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1. Simplified System
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationFinding numerically Newhouse sinks near a homoclinic tangency and investigation of their chaotic transients. Takayuki Yamaguchi
Hokkaido Mathematical Journal Vol. 44 (2015) p. 277 312 Finding numerically Newhouse sinks near a homoclinic tangency and investigation of their chaotic transients Takayuki Yamaguchi (Received March 13,
More informationPredator-Prey Population Models
21 Predator-Prey Population Models Tools Used in Lab 21 Hudson Bay Data (Hare- Lynx) Lotka-Volterra Lotka-Volterra with Harvest How can we model the interaction between a species of predators and their
More informationOutline. Learning Objectives. References. Lecture 2: Second-order Systems
Outline Lecture 2: Second-order Systems! Techniques based on linear systems analysis! Phase-plane analysis! Example: Neanderthal / Early man competition! Hartman-Grobman theorem -- validity of linearizations!
More informationModeling and Simulation Study of Mutuality Interactions with Type II functional Response and Harvesting
American Journal of Applied Mathematics 201; 6(3): 109-116 http://www.sciencepublishinggroup.com/j/ajam doi: 10.116/j.ajam.2010603.12 ISSN: 2330-003 (Print); ISSN: 2330-006X (Online) Modeling and Simulation
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More information2 Discrete growth models, logistic map (Murray, Chapter 2)
2 Discrete growth models, logistic map (Murray, Chapter 2) As argued in Lecture 1 the population of non-overlapping generations can be modelled as a discrete dynamical system. This is an example of an
More informationChapter 4. Transition towards chaos. 4.1 One-dimensional maps
Chapter 4 Transition towards chaos In this chapter we will study how successive bifurcations can lead to chaos when a parameter is tuned. It is not an extensive review : there exists a lot of different
More informationThe Eco Pyramid By Michael Stahl
The Eco Pyramid The Eco Pyramid By Michael Stahl An ecosystem is a community of living organisms interacting with one another as well as with nonliving things. One very important aspect of an ecosystem
More informationName Student ID. Good luck and impress us with your toolkit of ecological knowledge and concepts!
Page 1 BIOLOGY 150 Final Exam Winter Quarter 2000 Before starting be sure to put your name and student number on the top of each page. MINUS 3 POINTS IF YOU DO NOT WRITE YOUR NAME ON EACH PAGE! You have
More informationLecture 3. Dynamical Systems in Continuous Time
Lecture 3. Dynamical Systems in Continuous Time University of British Columbia, Vancouver Yue-Xian Li November 2, 2017 1 3.1 Exponential growth and decay A Population With Generation Overlap Consider a
More informationKey words and phrases. Bifurcation, Difference Equations, Fixed Points, Predator - Prey System, Stability.
ISO 9001:008 Certified Volume, Issue, March 013 Dynamical Behavior in a Discrete Prey- Predator Interactions M.ReniSagaya Raj 1, A.George Maria Selvam, R.Janagaraj 3.and D.Pushparajan 4 1,,3 Sacred Heart
More information(Refer Slide Time: 00:32)
Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 12 Scilab simulation of Lotka Volterra
More informationLong Term Dynamics for Two Three-Species Food Webs
Long Term Dynamics for Two Three-Species Food Webs B. Bockelman, E. Green, L. Lippitt, J. Sherman October 28, 23 Abstract In this paper, we analyze two possible scenarios for food webs with two prey and
More information