STABILITY IN A MODEL FOR GENETICALLY ALTERED MOSQUITOS WITH PERIODIC PARAMETERS
|
|
- Oliver James
- 6 years ago
- Views:
Transcription
1 STABILITY IN A MODEL FOR GENETICALLY ALTERED MOSQUITOS WITH PERIODIC PARAMETERS Hubertus F. von Bremen and Robert J. Sacker Department of Mathematics and Statistics, California State Polytechnic University, Pomona, USA, hfvonbremen@csupomona.edu Department of Mathematics, University of Southern California, Los Angeles, USA, rsacker@math.usc.edu Abstract: In this paper we introduce seasonal variability into a discrete-time population dynamics model given by Jia Li for genetically altered and wild mosquitoes. This is done by allowing the birth and survival rates to vary periodically. Using the new model we show that under certain conditions the system will have a globally asymptotically stable solution. Keywords: Periodic difference equation, global stability, genetically altered mosquitoes.. INTRODUCTION Mosquitoes transmit many fatal and debilitating diseases. Malaria is the major killer in sub-saharan Africa with more than one million deaths attributed to it every year. Yellow fever, dengue fever, West Nile virus, encephalitis and filariasis are additional mosquito-borne diseases that have a significant impact on populations worldwide. Recently a significant amount of research has been conducted to genetically modify mosquitoes in order to block the transmission of disease generating parasites. Catteruccia et al [] state that success in genetically modifying the Mediterranean fruit fly, the yellow fever mosquito (Aedes aegypti) has been achieved. Researchers have proposed the release of genetically modified mosquitoes into the wild. The altered mosquitoes would interact with the wild ones and reproduce. This development has fueled the need to create and study mathematical models of the population dynamics of wild and genetically modified mosquitoes. Jia Li [] has recently proposed a promising discrete-time mathematical model for populations consisting of wild and genetically altered mosquitoes. Li [] developed a two species model having a hybrid Ricatti/Ricker type nonlinearity and provided sufficient conditions guaranteeing the existence of locally asymptotically stable fixed point for mosquitoes with a constant mating rate. We [] recently showed that under less restrictive conditions the fixed point is actually globally asymptotically stable with respect to initial populations in which both species are present. We also investigated the case when the mortality rates (fixed as well as density dependent) are different for the two species, and showed that the model becomes very sensitive to small changes in the reproductive parameters. More recently, we showed [] that for proportional mating rates, under a certain set of conditions, the system will have two positive fixed points. The positive fixed point with the largest value is locally asymptotically stable, provided certain conditions are satisfied, and the other positive fixed point is unstable. In the proportional mating rate case the origin is also a fixed point and it is locally asymptotically stable. This creates a "threshold" effect, or allee effect, so that very small initial populations tend to die out. In this paper we present an extension of the population dynamics model (see [], [], []) that takes into account the effect of seasonal variability by allowing the parameters to be periodic. We consider the cases of proportional and constant mating rates where the birth and mortality rates are periodic, but with equal mortality rates for the wild and the genetically altered mosquitoes. The population dynamics model for genetically altered and wild mosquitos with periodic parameters is given in section. In section we consider the ratio dynamics and we show that the under some conditions the ratio dynamics produces a globally asymptotically stable periodic solution. The periodic solutions of the original system are explored in sections and, for the constant mating rate case and the proportional mating rate, respectively.. MODEL FOR THE POPULATION OF MOSQUITOES An extension of the discrete-time mathematical model for populations consisting of wild and genetically altered mosquitoes proposed by Jia Li [] for systems with periodic parameters is considered in this paper. The following description closely follows []. Let x n be the number of wild mosquitoes present at generation n. The population dynamics of the wild mosquitoes is described by the difference equation x n+ = f(n, x n )s(n, x n )x n, () where f is the birth function (per-capita rate of offspring production) and s is the survival probability (fraction of the offspring that survive). The survival probability is assumed to have a Ricker-type form s(n, x n ) = e d(n) k(n)xn. In the model by Jia Li the survival parameters d and k were fixed and did not depend on the generation n. In the formulation here presented, it will be assumed that the parameters are generation dependent and periodic. Let y n be the number of genetically altered mosquitoes present at generation n, and assume that before the wild and altered mosquitoes interact, the dynamics of the altered mosquito population is similar to the one of the wild mosquitoes. Once the altered mosquitoes are released into the wild mosquito habitat, the mosquito populations are governed by the system of difference equations
2 x n+ = f (n, x n, y n )x n e d(n) k(n)(xn+yn), () y n+ = f (n, x n, y n )y n e d(n) k(n)(xn+yn). It is assumed that both wild and altered mosquitoes have the same survival probability e d(n) k(n)(xn+yn). For x n, y n and (x n, y n ) (, ) the birth rate functions f and f are given by f (n, x n, y n ) = c(n n ) α (n)x n + β (n)y n, () f (n, x n, y n ) = c(n n ) α (n)x n + β (n)y n, where c(n n ) is the number of matings per individual, per unit time with N n =. At generation n the number of matings, per individual, with wild mosquitoes is c(n n )x n /( ) and with altered mosquitoes, c(n n )y n /( ). Let α (n) be the number of wild offspring that a wild mosquito produces through mating with a wild mosquito, and β (n) be the number of wild mosquitoes produced through mating with an altered mosquito at generation n. Similarly, α (n) and β (n) are the number of altered mosquitoes produced by the mating of altered mosquitoes with wild and altered mosquitoes respectively at generation n. Combining () and () gives the following set of difference equations that govern the interacting populations of wild and altered mosquitoes x n+ = () c(n n ) α (n)x n + β (n)y n x n e d(n) k(n)(xn+yn), y n+ = c(n n ) α (n)x n + β (n)y n y n e d(n) k(n)(x n+y n ). The mating rate depends on the population density. When the population is relatively small the mating rate will be assumed to be proportional to the total population, N n =, that is, c(n n ) = c N n. Once the population size exceeds a certain level, we expect the number of matings to saturate, and we assume the mating rate is constant, that is,. In the case of small population sizes we assume that the mating rate is proportional to the total population, and the mating rate becomes c(n n ) = c ( ). Using () and a proportional mating rate, we get x n+ = (a (n)x n + b (n)y n )x n e d(n) k(n)(x n+y n ), () y n+ = (a (n)x n + b (n)y n )y n e d(n) k(n)(x n+y n ), where a i (n) = c α i (n) and b i (n) = c β i (n), for i =,. In the constant mating rate case letting c(n n ) = c, and letting a i (n) = cα i (n) and b i (n) = cβ i (n), for i =,, equation () becomes x n+ = a (n)x n + b (n)y n x n e d(n) k(n)(x n+y n ), () y n+ = a (n)x n + b (n)y n y n e d(n) k(n)(x n+y n ). In () and () we assume that x n, y n and (x n, y n ) (, ) for n.. THE RATIO DYNAMICS AND THE STABILITY OF PERIODIC SOLUTIONS In this section we will study the ratio z n = x n /y n and we will show that the under some conditions the ratio dynamics has a globally asymptotically stable periodic solution. Consider the equations in () or (), then the ratio z n = x n /y n for both systems is given by z n+ = a (n)z n + b (n) a (n)z n + b (n) z n. () Due to the decay survival probability term in () and in () the populations can not grow indefinitely ( is not possible). In the following Lemma we show that () has a unique globally asymptotically stable positive periodic solution. Lemma Suppose a(n), b(n), c(n), d(n) >, z, with a(n), b(n), c(n), d(n) being periodic with period k. Consider the difference equation z n+ = F (n, z n ), (8) F (n, z) = a(n)z + b(n) z c(n)z + d(n). If b(n)/d(n) > and c(n)/a(n) >, then (8) has a unique positive periodic s-cycle which is Globally Asymptotically Stable and is such that s divides k. Proof. The result of this Lemma follows directly from a result given by Elaydi and Sacker in []. Let f n (z) = F (n, z). In order to apply the result in [] we need to show that for each n the functions f n satisfy the following properties: )f n : R + R + is continuous, )f n is concave, and ) There exists z and z such that f n (z ) > z and f n (z ) < z. Property ) follows directly from the positivity of a(n), b(n), c(n), d(n), z. Direct computations show that f n(z) < for all z > (see []), showing that Property ) is satisfied. For each n, f n (z) has the unique positive fixed point z = (d(n) b(n))/(a(n) c(n)). Let z > z and z < z, then due to the fact that f n is concave, positive and continuous for z > we have that f n (z ) > z and f n (z ) < z. This shows that Property ) is satisfied. From Lemma() equation () has a unique globally asymptotically stable periodic orbit. Let z n zˆ n, so we have that y n /x n / zˆ n = r n, where r n is a unique globally asymptotically stable periodic orbit.
3 The periodic point r(n) for () gives rise to an invariant set of lines in the (x, y) plane, i.e., the union S = S n of the collection of lines {S n }, where S n = {(x, y) : y/x = r(n)}. From the global asymptotic stability of the periodic orbit r n we then have that y n /x n r(n), i.e. the ω-limit set of any point(x, y) with x >, y > lies in the union of the lines S n. Then to study the solutions in S we can set y n = r(n)x n in the first equation of () (or of (), and in the second equation we can set x n = (/r(n))y n. Using these substitutions we get two uncoupled equations (on S n ), as it will be shown in the next two sections.. CONSTANT MATING RATE CASE In this section we briefly explore the stability of the periodic solutions of systems with a constant mating rate. Using the substitutions y n = r(n)x n in the first equation of (), and x n = (/r(n))y n in the second equation,we get the following two uncoupled equations (on S). x n+ = a (n) + b (n)r(n) x n e d(n) k(n)(+r(n))xn, (9) + r(n) y n+ = a (n) + b (n)r(n) y n e d(n) k(n)(+/r(n))y n. + r(n) Using the transformations x n = u n /(k(n)(+r(n))) and y n = v n /(k(n)( + /r(n))) we can rewrite (9) in terms of u and v as u n+ = u n e pu(n) un, () v n+ = v n e pv(n) vn, Population 8 8 Altered Mosquitoes, y n Wild Mosquitoes, x n 8 9 Figure Population level of wild and genetically altered mosquitoes as a function of generation number Population ratio x n /y n where p u (n) and p v (n) are periodic with period k and given by. p u (n) = d(n) + () log[ k(n + )( + r(n + ))(a (n) + b (n)r(n)) k(n)( + r(n)) ], p v (n) = d(n) + log[ k(n + )r(n)( + r(n + ))(a (n) + b (n)r(n)) k(n)r(n + )( + r(n)) ]. The system in () is just two uncoupled Ricker s equations with periodic parameters. The stability of the solutions of the original coupled system () can be studied by considering the stability of the solutions of (). When periodic solutions of the ratio dynamics are stable, and the periodic solutions of the decoupled system are also stable (), then the periodic solutions of the original system are also stable. The following example illustrates the result presented in section (). Consider the system in () with the following period- parameters a = [,, ], a = [,,.], b = [8., 8., 8.], b = [8.,.,.], k = [.,.,.] and d = [.,.,.]. This set of parameters satisfy the condition in Lemma (), that is, a (n), a (n), b (n), b (n) >, with a (n)/a (n) > and b (n)/b (n) >. The population level of wild and genetically altered mosquitoes as a function of generation number Figure Ratio of population level of wild to genetically altered mosquitos as a function of generation number is shown in Figure. The figure was generated using the system in (), with the above-mentioned parameter values and initial conditions of x =. and y =.. The values of the population levels suggest that the system eventually settles to a periodic solution, with period. The ratio z n = x n /y n of the population levels of wild to genetically altered mosquitos for the above system is plotted on Figure as a function of generation number. As expected from Lemma (), the ratio z n = x n /y n eventually settles to a periodic solution with period. It should be noted that convergence of the ratio z n = x n /y n to a periodic solution is not sufficient to guarantee that x n and y n individually converge to a solution with the same period. Consider the above example with the same parameters, except for k = [.,.,.]. The ratio z n = x n /y n is identical to the previous one. However, x n and y n qualitatively behave significantly different, as shown in Figure. We now have a period- solution, as opposed to the earlier period- solution.
4 8 Altered Mosquitoes, y n. Population 8 Wild Mosquitoes, x n q.. q Figure Population level of wild and genetically altered mosquitoes as a function of generation number Figure Level curves of the positive fixed point of the composite map as a function of q and q Let q u (n) = e p u(n) in (), then the u equation can be written as u n+ = u n q u (n)e un. () Next we consider the special case of k =. Suppose q u (n) is period two, that is, q u (n + ) = q u (n) for all n. Let q = q u () and q = q u () with both q and q being positive. To study the period- solutions of () we can consider z n+ = g(z n ), and () g(z n ) = (g g )(z n ), with g i (z) = q i ze z, and i =,. The fixed points of () are the period- solutions of the u equation in (). Note that the v equation in () can be written in the same form as equation (). In the period- Ricker equation in the form given in (), the origin is a fixed point, and for positive values of q u (n) = q < there are no positive fixed points. In the range < q < e, the period- Ricker equation has unique globally asymptotically stable positive fixed points depending on the value of q. In the period- case the origin is also a fixed point of the composite map (), and for < q < e and < q < e the composite map will have at most one positive fixed point. Figure shows level curves of the positive fixed point of the composite map as a function of q and q. The blank region towards the lower left corner of the plot is the region in the q and q parameter space in which the system has no positive fixed point. Level curves of the derivative of the composite map at the positive fixed point as a function of q and q are shown in Figure. In the parameter range shown < q < e and < q < e the magnitude of the derivative is always less than unity. Thus the fixed point of the composite map is stable, and therefore () will have globally asymptotically stable period two solutions. For a q q Figure Level curves of the derivative of the composite map at the positive fixed point as a function of q and q given value of q and q one can use Figure to estimate the value of the fixed point of the composite map and Figure to estimate the value of the derivative of the composite map at that fixed point. We will present more detailed results on the stability of solutions of the periodic Ricker equation in []. These results give parameter regions guaranteeing stability that are larger than the ones provided by Z. Zhou and X. Zou []. See also [8].. PROPORTIONAL MATING RATE CASE In this section we will focus on the proportional mating rate case (for small population size). As before, using the substitutions y n = r(n)x n in the first equation of (), and x n = (/r(n))y n in the second equation,we get the following two uncoupled equations (on S)
5 x n+ = () (a (n) + b (n)r(n))x ne d(n) k(n)(+r)x n, y n+ = (a (n)(/r(n)) + b (n))yne d(n) k(n)(+/r(n))y n.. The x equation in () is of the form x n+ = α(n)x ne β(n)x n, () with α(n) = (a (n) + b (n)r(n))e d(n) and β(n) = k(n)( + r(n)). Note that the y equation in () has the same general form as (). For the special case of period- parameters (α(n) = α and β(n) = β) it was shown in [] that if α > βe, then () will have two positive fixed points. An unstable fixed point will occur at x < /β, and a stable fixed point will occur at x > /β. The fixed point, x, with x > /β will be locally asymptotically stable if < βx <. Since the origin is also a stable fixed point this gives rise to a threshold effect, i.e., very small initial populations tend to die out. Using the transformations x n = u n /β(n) we can rewrite () in terms of u as p p Figure Level curves of the first positive fixed point of the composite map as a function of p and p.8... u n+ = p u (n)u ne u n, () where p u (n) is periodic with period k and given by p u (n) = α(n)β(n+)/β(n). Again, the y equation in () can also be written in same form as (). The following numerical results deal with the special case of k =. Suppose p u (n) is period two, that is, p u (n + ) = p u (n) for all n. Let p = p u () and p = p u () with both p and p being positive. To study the -periodic solutions of () we can consider z n+ = g(z n ), and () g(z n ) = (g g )(z n ), with g i (z) = p i z e z, and i =,. The fixed points of () are the period- solutions of (). As for the case of period- parameters, the origin is also a stable fixed point of the composite map () for the period- case. For small enough parameter values of p and p, the composite map () has no positive fixed point. Thus in this parameter range, both populations of mosquitoes will become extinct. As the parameter values become sufficiently large, the composite map () will have two positive fixed points. Figure shows level curves of the first positive fixed point of the composite map as a function of p and p. The blank region on the plot corresponds to the set of parameter values of p and p for which the origin is the only fixed point. For a fixed combination of parameter values p and p, Figure can be used to obtain an estimate of the value of the first positive fixed point of the composite map corresponding to the chosen parameters. For a fixed combination of parameter values p and p, Figure can be used to estimate the value p p Figure Level curves of the derivative of the composite map at the first fixed point as a function of p and p..
6 . Wild Mosquitoes, x n. p. Population p. Altered Mosquitoes, y n Figure 8 Level curves of the the second fixed point of the composite map as a function of p and p Figure Population level of wild and genetically altered mosquitoes as a function of generation number p p Figure 9 Level curves of the derivative of the composite map at the second fixed point as a function of p and p of the derivative of the composite map at the first positive fixed point. Figure shows level curves of the derivative of the composite map at the first fixed point of as a function of p and p. The derivative of the composite map is always greater than unity in the range of p and p shown. Thus the first fixed point of the composite map is always unstable in the range of p and p shown. Figure 8 shows level curves of the second positive fixed point of the composite map as a function of p and p. As before, the blank region on the figure corresponds to the set of parameter values of p and p for which the origin is the only fixed point. Figure 9 shows level curves of the derivative of the composite map at the second fixed point of as a function of p and p. The derivative of the composite map is less than unity in the range of p and p shown. Thus the second fixed point of the composite map is stable in the range of p and p shown. The following example illustrates the results presented in this and in section (). Consider the system in () with the following period- parameters a = [.,.], a = [.,.], b = [.,.], b = [.,.], k = [.,.] and d = [.,.]. This set of parameters satisfy the condition in Lemma, that is, a (n), a (n), b (n), b (n) >, with a (n)/a (n) > and b (n)/b (n) >. And therefore we conclude that the ratio dynamics of the system has a unique asymptotically stable periodic solution. The population level of wild and genetically altered mosquitoes as a function of generation number is shown in Figure. The figure was generated using the system in (), with the above-mentioned parameter values and initial conditions of x =. and y =.. The values of the population levels suggest that the system eventually settles to a periodic solution, with period. Note that if the initial conditions are chosen sufficiently small, the system will converge to the zero solution. The ratio z n = x n /y n of the population levels of wild to genetically altered mosquitos is plotted on Figure as a function of generation number. As expected from Lemma, the ratio z n = x n /y n eventually settles to a periodic solution with period.. CONCLUSION An extension of the discrete-time population dynamics model (see [], [], []) for genetically altered and wild mosquitoes that takes into account the effect of seasonal variability is presented in this paper. The seasonal variability is introduced to the mathematical model by allowing the birth and survival rates vary periodically. In this paper we consider proportional and constant mating rates with periodic birth and survival (mortality) rates. The survival rates are allowed to be periodic but they are taken to be equal for the two types of mosquitoes. The stability of the solutions of the new model with periodic parameters is studied by considering the ratio dynamics. We provide conditions that guarantee that the ratio dynamics converges to a globally asymptotically stable periodic orbit. In the con-
7 Population ratio z n =x n /y n generation(n) Figure Ratio of population level of wild to genetically altered mosquitos as a function of generation number Control of Autonomous Decision Support Based Systems. Edited by E P Hofer and E Reithmeier, 8,. [] S. Elaydi and R. J. Sacker, "Global stability of periodic orbits of nonautonomous difference equations", J Differential Eq, 8():8-,. [] R.J. Sacker and H.F. von Bremen, "Stability in the periodic Ricker equation", in preparation. [] Z. Zhou and X. Zou, "Stable Periodic Solutions in a Discrete Periodic Logistic Equation", Applied Mathematical Letters, :-,. [8] R. Kon, "Attenuant cycles of population models with preriodic carrying capacity", Journal of Difference Equations and Applications, (-):-,. stant mating rate case we further show that under some additional conditions the system has globally asymptotically stable periodic orbits. In the proportional mating rate case with periodic parameters the origin is a locally asymptotically stable fixed point. This creates a "threshold" effect, or allee effect, so that very small initial populations tend to die out. In the period- parameter case we show that in a certain parameter region the system has one locally unstable and one locally stable positive periodic orbit. ACKNOWLEDGMENTS H. von Bremen is thankful for the support received from the California State Polytechnic University-Pomona, RSCA (Research, Scholarship, and Creative Activity) grant and the generous support from FUNDUNESP. Robert J. Sacker was supported by the University of Southern California, Letters Arts and Sciences Faculty Development Grant. REFERENCES [] Flaminia Catteruccia, Tony Nolan, Thanasis G. Loukeris, Claudia Blass, Charalambos Savakis, Fotis C. Kafatos, and Andrea Crisanti, "Stable germline transformation of the malaria mosquito Anopheles stephensi", Nature, :99-9, June,. [] J. Li, "Simple mathematical models for mosquito populations with genetically altered mosquitos", Math. Bioscience, 89:9 9,. [] R.J. Sacker and H.F. von Bremen, "Global Asymptotic Stability in the Jia Li Model for Gentically Altered Mosquitoes", Appeared in Difference Equations and Discrete Dynamical Systems, Edited by L J S Allen, B Aubach, S Elaydi and R Sacker, 8,. [] R.J. Sacker and H.F. von Bremen, "Some Stability Results in a Model for Genetically Altered Mosquitoes", Proceedings of the th International Workshop on Dynamics and Control. Appeared in Modeling and
Global Asymptotic Stability in the Jia Li Model for Genetically Altered mosquitoes
Global Asymptotic Stability in the Jia Li Model for Genetically Altered mosquitoes Robert J. Sacker Department of Mathematics, University of Southern California, Los Angeles, CA 90089-2532 USA, e-mail:
More informationDynamic Reduction with applications to Mathematical Biology and other areas
Journal of Biological Dynamics Vol. 00, No. 00, Month-Month 200x, 1 15 Dynamic Reduction with applications to Mathematical Biology and other areas Robert J. Sacker and Hubertus F. von Bremen Department
More informationGlobal stability of periodic orbits of nonautonomous difference equations in population biology and the Cushing-Henson conjectures
Global stability of periodic orbits of nonautonomous difference equations in population biology and the Cushing-Henson conjectures Saber Elaydi Department of Mathematics Trinity University San Antonio,
More informationPeriodic difference equations, population biology and the Cushing-Henson conjectures
Periodic difference equations, population biology and the Cushing-Henson conjectures Saber Elaydi Department of Mathematics Trinity University San Antonio, Texas 78212, USA E-mail: selaydi@trinity.edu
More informationpopulation size at time t, then in continuous time this assumption translates into the equation for exponential growth dn dt = rn N(0)
Appendix S1: Classic models of population dynamics in ecology and fisheries science Populations do not grow indefinitely. No concept is more fundamental to ecology and evolution. Malthus hypothesized that
More informationPeriodic difference equations, population biology and the Cushing-Henson conjectures
Periodic difference equations, population biology and the Cushing-Henson conjectures Saber Elaydi Department of Mathematics Trinity University San Antonio, Texas 78212, USA E-mail: selaydi@trinity.edu
More informationNonautonomous Beverton-Holt Equations and the Cushing-Henson Conjectures
Trinity University Digital Commons @ Trinity Mathematics Faculty Research Mathematics Department 4-2005 Nonautonomous Beverton-Holt Equations and the Cushing-Henson Conjectures Saber Elaydi Trinity University,
More informationMath 142-2, Homework 2
Math 142-2, Homework 2 Your name here April 7, 2014 Problem 35.3 Consider a species in which both no individuals live to three years old and only one-year olds reproduce. (a) Show that b 0 = 0, b 2 = 0,
More informationA New Mathematical Approach for. Rabies Endemy
Applied Mathematical Sciences, Vol. 8, 2014, no. 2, 59-67 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.39525 A New Mathematical Approach for Rabies Endemy Elif Demirci Ankara University
More informationExtinction and the Allee Effect in an Age Structured Population Model
AMS Special Session: Difference Equations and Applications Extinction and the Allee Effect in an Age Structured Population Model Nika Lazaryan and Hassan Sedaghat Department of Mathematics Virginia Commonwealth
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 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 information11. S. Jang, Dynamics of a discrete host-parasitoid system with stocking, Discrete Dynamics
Sophia Jang Department of Mathematics and Statistics Texas Tech University Office: MA202 Phone: (806) 834-7006 Fax: (806) 742-1112 E-mail: sophia.jang@ttu.edu Publications 1. M. De Silva, S. Jang, Period-doubling
More information7.2. Differential Equations. Phase Plots. Example 1 The logistic equation. Phase Plots. Phase Plots. Phase Plots, Equilibria, and Stability
Differential Equations 7 7.2, Equilibria, and Stability Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. Phase plots provide a way to visualize the dynamics
More informationToday. Introduction to Differential Equations. Linear DE ( y = ky ) Nonlinear DE (e.g. y = y (1-y) ) Qualitative analysis (phase line)
Today Introduction to Differential Equations Linear DE ( y = ky ) Nonlinear DE (e.g. y = y (1-y) ) Qualitative analysis (phase line) Differential equations (DE) Carbon dating: The amount of Carbon-14 in
More informationNonautonomous difference equations: Open problems and conjectures
Trinity University Digital Commons @ Trinity Mathematics Faculty Research Mathematics Department 7-2003 Nonautonomous difference equations: Open problems and conjectures Saber Elaydi Trinity University,
More informationAsynchronous and Synchronous Dispersals in Spatially Discrete Population Models
SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 7, No. 2, pp. 284 310 c 2008 Society for Industrial and Applied Mathematics Asynchronous and Synchronous Dispersals in Spatially Discrete Population Models Abdul-Aziz
More informationPopulation is often recorded in a form of data set. Population of Normal, Illinois
Population is often recorded in a form of data set Population of Normal, Illinois 1 Population of Venezuela 2 Population of world up to 1850 3 Population of world 4 Population of world (carton) 5 Population
More informationModeling the Existence of Basic Offspring Number on Basic Reproductive Ratio of Dengue without Vertical Transmission
International Journal on Recent and Innovation Trends in Computing and Communication ISSN: 232-869 Modeling the Existence of Basic Offspring Number on Basic Reproductive Ratio of Dengue without Vertical
More informationAccepted Manuscript. Backward Bifurcations in Dengue Transmission Dynamics. S.M. Garba, A.B. Gumel, M.R. Abu Bakar
Accepted Manuscript Backward Bifurcations in Dengue Transmission Dynamics S.M. Garba, A.B. Gumel, M.R. Abu Bakar PII: S0025-5564(08)00073-4 DOI: 10.1016/j.mbs.2008.05.002 Reference: MBS 6860 To appear
More informationEssential Ideas of Mathematical Modeling in Population Dynamics
Essential Ideas of Mathematical Modeling in Population Dynamics Toward the application for the disease transmission dynamics Hiromi SENO Research Center for Pure and Applied Mathematics Department of Computer
More informationHistory and meaning of the word Ecology A. Definition 1. Oikos, ology - the study of the house - the place we live
History and meaning of the word Ecology A. Definition 1. Oikos, ology - the study of the house - the place we live B. Etymology study of the origin and development of a word 1. Earliest - Haeckel (1869)
More informationPopulation modeling of marine mammal populations
Population modeling of marine mammal populations Lecture 1: simple models of population counts Eli Holmes National Marine Fisheries Service nmfs.noaa.gov Today s lecture topics: Density-independent growth
More informationAedes aegypti Population Model with Integrated Control
Applied Mathematical Sciences, Vol. 12, 218, no. 22, 175-183 HIKARI Ltd, www.m-hiari.com https://doi.org/1.12988/ams.218.71295 Aedes aegypti Population Model with Integrated Control Julián A. Hernández
More informationStability Analysis of a Population Dynamics Model with Allee Effect
Stability Analysis of a Population Dynamics Model with Allee Effect Canan Celik Abstract In this study, we focus on the stability analysis of equilibrium points of population dynamics with delay when the
More informationEvolution 1 Star. 6. The different tools used during the beaks of finches lab represented. A. feeding adaptations in finches
Name: Date: 1. ccording to modern evolutionary theory, genes responsible for new traits that help a species survive in a particular environment will usually. not change in frequency. decrease gradually
More informationExistence of Almost Periodic Solutions of Discrete Ricker Delay Models
International Journal of Difference Equations ISSN 0973-6069, Volume 9, Number 2, pp. 187 205 (2014) http://campus.mst.edu/ijde Existence of Almost Periodic Solutions of Discrete Ricker Delay Models Yoshihiro
More informationQualitative Analysis of a Discrete SIR Epidemic Model
ISSN (e): 2250 3005 Volume, 05 Issue, 03 March 2015 International Journal of Computational Engineering Research (IJCER) Qualitative Analysis of a Discrete SIR Epidemic Model A. George Maria Selvam 1, D.
More informationBIO S380T Page 1 Summer 2005: Exam 2
BIO S380T Page 1 Part I: Definitions. [5 points for each term] For each term, provide a brief definition that also indicates why the term is important in ecology or evolutionary biology. Where I ve provided
More informationPopulations Study Guide (KEY) All the members of a species living in the same place at the same time.
Populations Study Guide (KEY) 1. Define Population. All the members of a species living in the same place at the same time. 2. List and explain the three terms that describe population. a. Size. How large
More informationThursday. Threshold and Sensitivity Analysis
Thursday Threshold and Sensitivity Analysis SIR Model without Demography ds dt di dt dr dt = βsi (2.1) = βsi γi (2.2) = γi (2.3) With initial conditions S(0) > 0, I(0) > 0, and R(0) = 0. This model can
More informationSupplementary Information
Supplementary Information This document shows the supplementary figures referred to in the main article. The contents are as follows: a. Malaria maps b. Dengue maps c. Yellow fever maps d. Chikungunya
More informationPopulation Ecology Density dependence, regulation and the Allee effect
2/22/15 Population Ecology Density dependence, regulation and the Allee effect ESRM 450 Wildlife Ecology and Conservation Wildlife Populations Groups of animals, all of the same species, that live together
More informationA DISCRETE-TIME HOST-PARASITOID MODEL
A DISCRETE-TIME HOST-PARASITOID MODEL SOPHIA R.-J. JANG AND JUI-LING YU We study a discrete-time host-parasitoid model proposed by May et al. In this model, the parasitoid attacks the host first then followed
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 informationOn the global dynamics of periodic triangular maps
arxiv:1605.09631v1 [math.ds] 31 May 2016 On the global dynamics of periodic triangular maps Rafael Luís 1,2 2 University of Madeira, Funchal, Madeira, Portugal. 3 Center for Mathematical Analysis, Geometry,
More informationChapter 2: Growth & Decay
Chapter 2: Growth & Decay 107/226 Introduction In this chapter, model biological systems whose population is so large or where growth is so fine-grained that continuity can be assumed. These continuous
More informationLevels of Ecological Organization. Biotic and Abiotic Factors. Studying Ecology. Chapter 4 Population Ecology
Chapter 4 Population Ecology Lesson 4.1 Studying Ecology Levels of Ecological Organization Biotic and Abiotic Factors The study of how organisms interact with each other and with their environments Scientists
More informationChapter 4 Population Ecology
Chapter 4 Population Ecology Lesson 4.1 Studying Ecology Levels of Ecological Organization The study of how organisms interact with each other and with their environments Scientists study ecology at various
More informationMath 266: Autonomous equation and population dynamics
Math 266: Autonomous equation and population namics Long Jin Purdue, Spring 2018 Autonomous equation An autonomous equation is a differential equation which only involves the unknown function y and its
More informationSchool on Modelling Tools and Capacity Building in Climate and Public Health April Point Event Analysis
2453-12 School on Modelling Tools and Capacity Building in Climate and Public Health 15-26 April 2013 Point Event Analysis SA CARVALHO Marilia PROCC FIOCRUZ Avenida Brasil 4365 Rio De Janeiro 21040360
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 information1.Introduction: 2. The Model. Key words: Prey, Predator, Seasonality, Stability, Bifurcations, Chaos.
Dynamical behavior of a prey predator model with seasonally varying parameters Sunita Gakkhar, BrhamPal Singh, R K Naji Department of Mathematics I I T Roorkee,47667 INDIA Abstract : A dynamic model based
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 66, No. 5, pp. 1563 1587 c 2006 Society for Industrial and Applied Mathematics DISCRETE-TIME SIS EPIDEMIC MODEL IN A SEASONAL ENVIRONMENT JOHN E. FRANKE AND ABDUL-AZIZ YAKUBU Abstract.
More informationContinuous-Time Dynamical Systems: Sketching solution curves and bifurcations
Math 360 Winter 07 Section 0 HW#4 Continuous-Time Dynamical Systems: Sketching solution curves and bifurcations (Due in one week! Thursday Nov 9, 07) Materials contained in this assignment are of great
More informationChapter 2 Lecture. Density dependent growth and intraspecific competition ~ The Good, The Bad and The Ugly. Spring 2013
Chapter 2 Lecture Density dependent growth and intraspecific competition ~ The Good, The Bad and The Ugly Spring 2013 2.1 Density dependence, logistic equation and carrying capacity dn = rn K-N Dt K Where
More informationEcology! the study of the relationship between organisms and their environment
Ecology Notes Ecology! the study of the relationship between organisms and their environment An ecosystem is an organism s surroundings consisting of both living and nonliving things. Abiotic NONliving
More informationMathematical Model of Proportional Release on Population Suppression
Applied Mathematical Sciences, Vol. 12, 2018, no. 8, 407-415 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.711338 Mathematical Model of Proportional Release on Population Suppression Bo
More informationSYMBIOTIC MODELS WITH AN SIR DISEASE
SYMBIOTIC MODELS WITH AN SIR DISEASE A Thesis Presented to the Faculty of California State Polytechnic University, Pomona In Partial Fulfillment Of the Requirements for the Degree Master of Science In
More informationPopulation Models with Allee Effect: A New Model
Trinity University Digital Commons @ Trinity Mathematics Faculty Research Mathematics Department 7-2010 Population Models with Allee Effect: A New Model Saber Elaydi Trinity University, selaydi@trinity.edu
More informationMathematical Modeling and Analysis of Infectious Disease Dynamics
Mathematical Modeling and Analysis of Infectious Disease Dynamics V. A. Bokil Department of Mathematics Oregon State University Corvallis, OR MTH 323: Mathematical Modeling May 22, 2017 V. A. Bokil (OSU-Math)
More informationPopulation Ecology and the Distribution of Organisms. Essential Knowledge Objectives 2.D.1 (a-c), 4.A.5 (c), 4.A.6 (e)
Population Ecology and the Distribution of Organisms Essential Knowledge Objectives 2.D.1 (a-c), 4.A.5 (c), 4.A.6 (e) Ecology The scientific study of the interactions between organisms and the environment
More informationSlow and fast dynamics model of a Malaria with Sickle-Cell genetic disease with multi-stage infections of the mosquitoes population
Journal of Physics: Conference Series PAPER OPEN ACCESS Slow and fast dynamics model of a Malaria with Sickle-Cell genetic disease with multi-stage infections of the mosquitoes population To cite this
More informationMATH 415, WEEKS 14 & 15: 1 Recurrence Relations / Difference Equations
MATH 415, WEEKS 14 & 15: Recurrence Relations / Difference Equations 1 Recurrence Relations / Difference Equations In many applications, the systems are updated in discrete jumps rather than continuous
More informationGlobal Stability of a Computer Virus Model with Cure and Vertical Transmission
International Journal of Research Studies in Computer Science and Engineering (IJRSCSE) Volume 3, Issue 1, January 016, PP 16-4 ISSN 349-4840 (Print) & ISSN 349-4859 (Online) www.arcjournals.org Global
More informationIntroduction to SEIR Models
Department of Epidemiology and Public Health Health Systems Research and Dynamical Modelling Unit Introduction to SEIR Models Nakul Chitnis Workshop on Mathematical Models of Climate Variability, Environmental
More informationMA 137 Calculus 1 with Life Science Application A First Look at Differential Equations (Section 4.1.2)
MA 137 Calculus 1 with Life Science Application A First Look at Differential Equations (Section 4.1.2) Alberto Corso alberto.corso@uky.edu Department of Mathematics University of Kentucky October 12, 2015
More informationCh 5. Evolution, Biodiversity, and Population Ecology. Part 1: Foundations of Environmental Science
Ch 5 Evolution, Biodiversity, and Population Ecology Part 1: Foundations of Environmental Science PowerPoint Slides prepared by Jay Withgott and Heidi Marcum Copyright 2006 Pearson Education, Inc., publishing
More informationPredation is.. The eating of live organisms, regardless of their identity
Predation Predation Predation is.. The eating of live organisms, regardless of their identity Predation 1)Moves energy and nutrients through the ecosystem 2)Regulates populations 3)Weeds the unfit from
More informationPopulation Ecology & Biosystematics
Population Ecology & Biosystematics Population: a group of conspecific individuals occupying a particular place at a particular time This is an operational definition Compare with Deme: a population unevenly
More informationPopulation and Community Dynamics
Population and Community Dynamics Part 1. Genetic Diversity in Populations Pages 676 to 701 Part 2. Population Growth and Interactions Pages 702 to 745 I) Introduction I) Introduction to understand how
More informationEcology- Population Growth Rate Problems
AP Biology Name Ecology- Population Growth Rate Problems Birth rate b = B/N Death rate m = D/N Individual or Population Growth Rate (per capita) r=(b-d)/n or r = b-m Exponential Growth Rate Calculation
More informationA Mathematical Model for Transmission of Dengue
Applied Mathematical Sciences, Vol. 10, 2016, no. 7, 345-355 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.510662 A Mathematical Model for Transmission of Dengue Luis Eduardo López Departamento
More informationMODELING AND ANALYSIS OF THE SPREAD OF CARRIER DEPENDENT INFECTIOUS DISEASES WITH ENVIRONMENTAL EFFECTS
Journal of Biological Systems, Vol. 11, No. 3 2003 325 335 c World Scientific Publishing Company MODELING AND ANALYSIS OF THE SPREAD OF CARRIER DEPENDENT INFECTIOUS DISEASES WITH ENVIRONMENTAL EFFECTS
More informationDiscrete Time Markov Chain of a Dynamical System with a Rest Phase
Discrete Time Markov Chain of a Dynamical System with a Rest Phase Abstract A stochastic model, in the form of a discrete time Markov chain, is constructed to describe the dynamics of a population that
More informationMathematical Models of Biological Systems
Mathematical Models of Biological Systems CH924 November 25, 2012 Course outline Syllabus: Week 6: Stability of first-order autonomous differential equation, genetic switch/clock, introduction to enzyme
More informationMA 777: Topics in Mathematical Biology
MA 777: Topics in Mathematical Biology David Murrugarra Department of Mathematics, University of Kentucky http://www.math.uky.edu/~dmu228/ma777/ Spring 2018 David Murrugarra (University of Kentucky) Lecture
More informationAN ABSTRACT OF THE THESIS OF. Margaret-Rose W. Leung for the degree of Honors Baccalaureate of Science in Mathematics
AN ABSTRACT OF THE THESIS OF Margaret-Rose W. Leung for the degree of Honors Baccalaureate of Science in Mathematics presented on June 5, 2012. Title: A Vector Host Model for Coinfection by Barley/Cereal
More informationModeling with differential equations
Mathematical Modeling Lia Vas Modeling with differential equations When trying to predict the future value, one follows the following basic idea. Future value = present value + change. From this idea,
More informationDynamical Systems: Lecture 1 Naima Hammoud
Dynamical Systems: Lecture 1 Naima Hammoud Feb 21, 2017 What is dynamics? Dynamics is the study of systems that evolve in time What is dynamics? Dynamics is the study of systems that evolve in time a system
More information1 Mathematical models of the population growth
1 Mathematical models of the population growth 1.1 How many people can the Earth support? (Or a quick refresher on the least square method) I would like to start with a very simple and yet interesting
More informationEvolution - Unifying Theme of Biology Microevolution Chapters 13 &14
Evolution - Unifying Theme of Biology Microevolution Chapters 13 &14 New Synthesis Natural Selection Unequal Reproductive Success Examples and Selective Forces Types of Natural Selection Speciation http://www.biology-online.org/2/11_natural_selection.htm
More informationA branching process for the early spread of a transposable element in a diploid population
J. Math. Biol. (2008) 57:811 840 DOI 10.1007/s00285-008-0190-2 Mathematical Biology A branching process for the early spread of a transposable element in a diploid population John M. Marshall Received:
More informationThe fate of beneficial mutations in a range expansion
The fate of beneficial mutations in a range expansion R. Lehe (Paris) O. Hallatschek (Göttingen) L. Peliti (Naples) February 19, 2015 / Rutgers Outline Introduction The model Results Analysis Range expansion
More informationScholarship 2008 Biology
93101Q S 931012 Scholarship 2008 Biology 9.30 am Saturday 22 November 2008 Time allowed: Three hours Total marks: 24 QUESTION booklet There are THREE questions in this booklet. Answer ALL questions. Write
More informationLocal and Global Stability of Host-Vector Disease Models
Local and Global Stability of Host-Vector Disease Models Marc 4, 2008 Cruz Vargas-De-León 1, Jorge Armando Castro Hernández Unidad Académica de Matemáticas, Universidad Autónoma de Guerrero, México and
More informationWe have two possible solutions (intersections of null-clines. dt = bv + muv = g(u, v). du = au nuv = f (u, v),
Let us apply the approach presented above to the analysis of population dynamics models. 9. Lotka-Volterra predator-prey model: phase plane analysis. Earlier we introduced the system of equations for prey
More information3-LS1-1 From Molecules to Organisms: Structures and Processes
3-LS1-1 From Molecules to Organisms: Structures and Processes 3-LS1-1. Develop models to describe that organisms have unique and diverse life cycles but all have in common birth, growth, reproduction,
More informationMathematical Ecology
Mathematical Ecology Joachim Hermisson, Claus Rueffler & Meike Wittmann January 14, 2016 Literature and Software Sarah P. Otto, Troy Day: A Biologist s Guide to Mathematical Modeling in Ecology and Evolution,
More informationCHAPTER. Population Ecology
CHAPTER 4 Population Ecology Lesson 4.1 Studying Ecology Ernst Haeckel defined ecology in 1866 as the body of knowledge concerning the economy of nature the total relations of the animal to both its inorganic
More informationSimple Mathematical Model for Malaria Transmission
Journal of Advances in Mathematics and Computer Science 25(6): 1-24, 217; Article no.jamcs.37843 ISSN: 2456-9968 (Past name: British Journal of Mathematics & Computer Science, Past ISSN: 2231-851) Simple
More informationAnalysis of a Vaccination Model for Carrier Dependent Infectious Diseases with Environmental Effects
Nonlinear Analysis: Modelling and Control, 2008, Vol. 13, No. 3, 331 350 Analysis of a Vaccination Model for Carrier Dependent Infectious Diseases with Environmental Effects Ram Naresh 1, Surabhi Pandey
More informationOptimal Harvesting Models for Fishery Populations
Optimal Harvesting Models for Fishery Populations Corinne Wentworth St. Mary s College of Maryland Mentored by: Dr. Masami Fujiwara and Dr. Jay Walton July 28, 2011 Optimal Harvesting Models for Fishery
More informationCHAPTER 5. Interactions in the Ecosystem
CHAPTER 5 Interactions in the Ecosystem 1 SECTION 3.3 - THE ECOSYSTEM 2 SECTION 3.3 - THE ECOSYSTEM Levels of Organization Individual one organism from a species. Species a group of organisms so similar
More informationGeographical Information System (GIS)-based maps for monitoring of entomological risk factors affecting transmission of chikungunya in Sri Lanka
Geographical Information System (GIS)-based maps for monitoring of entomological risk factors affecting transmission of chikungunya in Sri Lanka M.D. Hapugoda 1, N.K. Gunewardena 1, P.H.D. Kusumawathie
More informationMultiple choice 2 pts each): x 2 = 18) Essay (pre-prepared) / 15 points. 19) Short Answer: / 2 points. 20) Short Answer / 5 points
P 1 Biology 217: Ecology Second Exam Fall 2004 There should be 7 ps in this exam - take a moment and count them now. Put your name on the first p of the exam, and on each of the ps with short answer questions.
More informationGlobal Analysis of an SEIRS Model with Saturating Contact Rate 1
Applied Mathematical Sciences, Vol. 6, 2012, no. 80, 3991-4003 Global Analysis of an SEIRS Model with Saturating Contact Rate 1 Shulin Sun a, Cuihua Guo b, and Chengmin Li a a School of Mathematics and
More informationMathematical Methods in Biology Part 3: Vectors and matrices
Mathematical Methods in Biology Part 3: Vectors and matrices Eva Kisdi Department of Mathematics and Statistics University of Helsinki c Eva Kisdi. Any part of this material may be copied or re-used only
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 informationMathematical modelling and controlling the dynamics of infectious diseases
Mathematical modelling and controlling the dynamics of infectious diseases Musa Mammadov Centre for Informatics and Applied Optimisation Federation University Australia 25 August 2017, School of Science,
More informationStudent Background Readings for Bug Dynamics: The Logistic Map
Student Background Readings for Bug Dynamics: The Logistic Map Figure 1 The population as a function of a growth rate parameter observed in the experiment of [Bugs]. Graphs such as this are called bifurcation
More informationHow to Use This Presentation
How to Use This Presentation To View the presentation as a slideshow with effects select View on the menu bar and click on Slide Show. To advance through the presentation, click the right-arrow key or
More informationA.P. Biology CH Population Ecology. Name
1 A.P. Biology CH. 53 - Population Ecology Name How many ants (shown below - 6 ants / cm 2 ) would there be in an ant colony that is flat and one meter long on each side? Dispersion Patterns Matching A
More informationFirst Order Differential Equations
First Order Differential Equations 1 Finding Solutions 1.1 Linear Equations + p(t)y = g(t), y() = y. First Step: Compute the Integrating Factor. We want the integrating factor to satisfy µ (t) = p(t)µ(t).
More informationMath 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations
Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations Junping Shi College of William and Mary, USA Equilibrium Model: x n+1 = f (x n ), here f is a nonlinear function
More informationOn some Concepts of (a, b, c)-trichotomy for Noninvertible Linear Discrete-Time Systems in Banach Spaces
DOI: 10.2478/awutm-2014-0016 Analele Universităţii de Vest, Timişoara Seria Matematică Informatică LII, 2, (2014), 145 155 On some Concepts of (a, b, c)-trichotomy for Noninvertible Linear Discrete-Time
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 informationA population is a group of individuals of the same species occupying a particular area at the same time
A population is a group of individuals of the same species occupying a particular area at the same time Population Growth As long as the birth rate exceeds the death rate a population will grow Immigration
More informationOptimal Translocation Strategies for Threatened Species
Optimal Translocation Strategies for Threatened Species Rout, T. M., C. E. Hauser and H. P. Possingham The Ecology Centre, University of Queensland, E-Mail: s428598@student.uq.edu.au Keywords: threatened
More informationSpotlight on Modeling: The Possum Plague
70 Spotlight on Modeling: The Possum Plague Reference: Sections 2.6, 7.2 and 7.3. The ecological balance in New Zealand has been disturbed by the introduction of the Australian possum, a marsupial the
More information