Mathematical modeling and analysis of the depletion of dissolved oxygen in water bodies
|
|
- Violet Greer
- 5 years ago
- Views:
Transcription
1 Nonlinear Analysis: Real World Applications Mathematical modeling and analysis of the depletion of dissolved oxygen in water bodies A.K. Misra a, Peeyush Chandra a,, J.B. Shukla b a Department of Mathematics, Indian Institute of echnology, Kanpur 20806, India b Centre for Modeling Environment and Development, Kanpur , India Received 4 January 2005; accepted 5 September 2005 Abstract In this paper, a nonlinear mathematical model is proposed to study the depletion of dissolved oxygen in a water body caused by industrial and household discharges of organic matters pollutants. he problem is formulated as a food chain model by considering various interaction processes biodegradation and biochemical involving organic pollutants, bacteria, protozoa, an aquatic population and dissolved oxygen. Using stability theory, it is shown that as the rate of introduction of organic pollutants in a water body increases, the concentration of dissolved oxygen decreases due to various interaction processes. It is found that if the organic pollutants are continuously discharged into water body, the concentration of dissolved oxygen may become negligibly small, thus threatening the survival of aquatic populations. However, by using some effort to control the cumulative discharge of these pollutants into the water body, the concentration of dissolved oxygen can be maintained at a desired level Elsevier Ltd. All rights reserved. Keywords: Mathematical model; Water pollution; Dissolved oxygen; Monod interaction. Introduction It is well known that pollutants and toxicants discharged in water bodies are harmful to biological species such as a fish population 0,5,9,22,23. When household and industrial wastes are discharged into water, the organic matter present in them are uptaken by micro-organisms which convert them into inorganic matter using dissolved oxygen in the process. he density of these micro-organisms bacteria increases as the cumulative concentration of organic pollutants increases. As these are part of a food chain involving protozoa and other biological populations in a water body, the level of dissolved oxygen decreases further due to various interactive biochemical and biodegradation processes. It may be noted here that the input of the dissolved oxygen in a water body is mainly due to atmospheric diffusion through the water surface and to a certain extent due to its production by photosynthesis 3,0,8 20. herefore, while modeling the depletion of dissolved oxygen in an aquatic ecosystem, all these processes should be taken into account 9. he study of effect of the organic pollutants on dissolved oxygen in a water body started with the well-known model of Streeter and Phelps 2. Several workers including Dobbins 3, O Connor 4, Beck and Young have generalized this model. It may be pointed out here that these models are only linear. However, there exist various nonlinear Corresponding author. address: peeyush@iitk.ac.in P. Chandra /$ - see front matter 2005 Elsevier Ltd. All rights reserved. doi:0.06/j.nonrwa
2 A.K. Misra et al. / Nonlinear Analysis: Real World Applications processes in water bodies involving interactions such as organic pollutants with bacteria 6,,5, phytoplankton with zooplankton in presence of nutrients 2,8, organic chemical distribution in aquatic food chains 22, dissolved oxygen with aquatic populations 5, etc. For a case related to organic pollutants, Rinaldi et al. 5 proposed ecological type nonlinear models involving fast and slow biodegradable pollutants, which implicitly consider the depletion of dissolved oxygen by bacteria and protozoa in a water body. In this study, where simulation analysis has been carried out, several questions regarding the behavior of the nonzero equilibrium including nonlinear stability, etc. remained unanswered. Further some other aspects, such as natural depletion rates of pollutants, crowding of species, growth of oxygen from algae and macrophytes, survival of aquatic population dependent on dissolved oxygen, etc. have not been taken into account in various nonlinear modeling studies. It may be also pointed out here that, the level of dissolved oxygen can be increased by pumping air into the water body, being used for production of fish, and therefore, this aspect needs also to be considered in the modeling process. In this paper, therefore, a nonlinear ecological type of mathematical model for the depletion of dissolved oxygen in a water body is proposed and analyzed by considering some of the aspects mentioned above involving degradable organic pollutants, bacteria, protozoa and an aquatic population. his model takes into account monod interactions between organic pollutants and bacteria, bacteria and protozoa and a bilinear interaction between protozoa and the aquatic population. It is assumed that a constant rate of input of oxygen is applied into the water body, to increase the level of dissolved oxygen 3,8,9. A model to control the discharge of organic pollutants and its effect on dissolved oxygen is also proposed and analyzed. 2. Mathematical model We consider a water body where organic pollutants are discharged in the form of wastes. It is assumed that these pollutants are part of a food chain consisting of bacteria, protozoa and an aquatic population, for example a fish population using dissolved oxygen in this water body for various biochemical and biodegradation processes. Let, be the cumulative concentration of organic matter degradable pollutants, B be the density of the bacteria, P be the density of the protozoa, N be the density of an aquatic population, which depends wholly on protozoa and C be the concentration of the dissolved oxygen DO. It is assumed that the cumulative discharge of organic pollutants into the water body is at a constant rate Q and that the rate of decrease of concentration due to biochemical and other factors is proportional to. It is further assumed that the cumulative rate of depletion of due to bacteria is given by a monod type of interaction involving the density of bacteria B as well as the concentration i.e. B/K 2 + K. As bacteria wholly depend upon organic pollutants, the growth rate of density of its population is proportional to this term. It is considered that the natural depletion rate of density of bacteria is proportional to B and the decrease of its growth rate due to crowding is proportional to B 2. Further, the depletion rate of density of bacteria by its predator protozoa is assumed to be given by the monod interaction between bacteria and protozoa i.e. BP/K 2 + K 22 B. hus, the growth rate of protozoa density is also proportional to this interaction term. he natural depletion rate of protozoa density is assumed proportional to P while its decrease due to crowding is proportional top 2. he depletion rate of protozoa density by its predator such as a fish population is assumed proportional to PN and therefore the corresponding growth rate of density of this predator population is proportional to this product. he natural depletion rate of density of predator population is considered proportional to N and its decrease due to crowding is proportional to N 2. It is also assumed that the rate of growth of concentration of dissolved oxygen by various sources including diffusion, pumping of air into the water body, etc. is q a constant and its natural depletion rate is proportional to C.Itis assumed, further, that the rate of depletion of dissolved oxygen is proportional to various terms K B/K 2 + K, K 2 BP/K 2 + K 22 B, K 3 PN, α B,α 2 P,α 3 N representing growth and depletion rates 5. Keeping in mind the above considerations, the system is governed by the following differential equations: d db = Q α 0 K B K 2 + K, K B = λ K 2 + K α B λ 0 B 2 K 2 BP K 2 + K 22 B,
3 982 A.K. Misra et al. / Nonlinear Analysis: Real World Applications dp dn dc = λ 2 K 2 BP K 2 + K 22 B α 2P λ 20 P 2 K 3 PN, 2. = λ 3 K 3 PN α 3 N λ 30 N 2, K B = q α 4 C λ 34 K 3 PN λ α B λ 22 α 2 P λ 33 α 3 N λ 2 K 2 + K λ K 2 BP 23 K 2 + K 22 B, where 0 = 0 0, B0 = B 0 0, P0 = P 0 = P 0 0, N0 = N 0 0, C0 = C 0 0. In model 2. the constants λ 0, λ 20, λ 30 are crowding coefficients flaking off coefficients 5 of bacteria, protozoa and the aquatic population, respectively. he parameters K,K 2,K 3,K,K 2,K 2,K 22 involved in the monod type interactions are all constants. he coefficients λ, λ 2, λ 3, λ, λ 22, λ 33, λ 2, λ 23, λ 34 are positive constants and their magnitudes are less than equal to unity. It may be pointed out that for feasibility of model 2., the growth rates of bacteria and protozoa should be positive. Hence, from the second and third equation of the model, it follows that and λ K K α > λ 2 K 2 K 22 α 2 > Following, Freedman and So 5, the region of attraction Ω for all solutions initiating in the positive octant is given by } Ω = {0 Qα0, 0 B R B, 0 P R P, 0 N R N, 0 C qα4, 2.4 where R B = R P = λ K Q λ 0 K Q + K 2 α 0, 2.5 λ λ 2 K K 2 Q λ 20 K 22 λ K Q + λ 0 K 2 K Q + K 2 α 0, 2.6 λ λ 2 λ 3 K K 2 K 3 Q R N = λ 20 λ 30 K 22 λ K Q + λ 0 K 2 K Q + K 2 α It is assumed here that all the initial values of the variables belong to set Ω. In the following, we analyze this model 2. by using the qualitative theory of differential 4,7,2,6,7. 3. Equilibrium analysis System 2. has following four nonnegative equilibria in Ω. hey are listed below: i E Q/α 0, 0, 0, 0,q/α 4 which always exists. ii E 2 2,B 2, 0, 0,C 2 exists, provided the following conditions are satisfied: λ K K α Q K 2 α 0 α > 0, 3. q λ 2 K 2 B 2 K 2 + K 2 λ α B2 >
4 A.K. Misra et al. / Nonlinear Analysis: Real World Applications iii E 3 3,B 3,P 3, 0,C 3 exists, provided the following conditions are valid: λ 20 λ K K α Q K 2 α 0 α + K 2α 2 K Q + K 2 α 0 >0, K λ 0 + K 2 λ 2 K 2 K 2 λ2 K 2 R B λ 20 K 2 + K 22 R B 2 K 22 α 2 > 0, K 2 + K 22 R B K 2 + K 22 R B 3.4 λ 2 K 2 B 3 K 2 + K 22 B 3 α 2 > 0, 3.5 q λ 2 K 3 B 3 K 2 + K 3 λ 23 K 2 3 B 3 K 2 + K 22 B 3 λ α B3 λ 22α 2 P3 > iv E 4,B,P,N,C exists whenever the following conditions are satisfied: λ K K α Q K 2 α 0 α + K 2m K 2 l K Q + K 2 α 0 >0, 3.7 K 2 λ30 λ 2 K 2 K 2 λ 0 + K 2 + K 22 R B 2 l K 2 + K 22 R B K λ30 λ 2 K 2 R B 22 K 2 + K 22 R B m > 0, 3.8 λ 30 λ 2 K 2 B K 2 + K 22 B m>0, 3.9 λ 3 K 3 P α 3 > 0, 3.0 K B q λ 2 K 2 + K λ K 2 B P 23 K 2 + K 22 B λ 34K 3 P N λ α B λ 22 α 2 P λ 33 α 3 N > 0, 3. where l = λ 20 λ 30 + λ 3 K 2 3, m = λ 30α 2 K 3 α 3. Proof. he equilibria E i i =, 2, 3, 4 of system 2. are obtained by solving the following set of simultaneous algebraic equations: K B Q α 0 = 0, 3.2 K 2 + K λ K B K 2 + K α B λ 0 B 2 K 2 BP = 0, 3.3 K 2 + K 22 B λ 2 K 2 BP K 2 + K 22 B α 2P λ 20 P 2 = 0, 3.4 λ 3 K 3 PN α 3 N λ 30 N 2 = 0, 3.5 K B q α 4 C λ 2 K 2 + K λ K 2 BP 23 K 2 + K 22 B λ 34K 3 PN λ α B λ 22 α 2 P λ 33 α 3 N = he equilibrium E Q/α 0, 0, 0, 0,q/α 4 exists obviously. We show the existence of other equilibria as follows. Existence of E 2 : for the equilibrium E 2 2,B 2, 0, 0,C 2, the values of 2,B 2 and C 2 are obtained by solving the following algebraic equations: B = Q α 0 K + K 2, 3.7 K B = λ K λ 0 K 2 + K α. 3.8
5 984 A.K. Misra et al. / Nonlinear Analysis: Real World Applications B Fig.. Existence of 2,B 2. From the isocline given by 3.7, it is easy to note the following: i B = 0at = Q/α 0 ii B as 0, iii db d = K 2 Q α 0 2 K α 0 K K + K 2, which is negative for 0 Q/α 0. hus B decreases as increases. From the isocline given by 3.8, it is noted that B<0at = 0 and B = 0at = K 2 α /λ K K α, which is positive. Further db/d >0 and hence B increases for all >0. hus, the two isoclines 3.7 and 3.8 intersect at 2, B 2 see Fig. provided λ K K α Q K 2 α 0 α > 0. Using these values of 2 and B 2,wegetC 2 from 3.6 as follows: C 2 = α 4 q λ 2 K 2 B 2 K 2 + K 2 λ α B 2, 3.9 which is positive provided the right-hand side of 3.9 is positive. his gives condition 3.2. hus, E 2 2,B 2, 0, 0,C 2 exists, provided conditions 3. and 3.2 are satisfied. Existence of E 3 : for the equilibrium E 3 3,B 3,P 3, 0,C 3, the values of 3 and B 3 are given by the following algebraic equations: B = Q α 0 K + K K and λ K K 2 + K α = λ 0 B + K 2 K 2 + K 22 B λ 20 λ2 K 2 B K 2 + K 22 B α Isocline 3.20 is the same as 3.7, hence its behavior is also the same. From the isocline given by 3.2, it is noted that B increases with in 0 B R B R B is given by 2.5, provided inequality 3.4 is satisfied. Also B = 0at = K 2 α K / 2α 2 K 2 α 2 λ K K α + K, λ 20 K 2 λ 20 K 2
6 A.K. Misra et al. / Nonlinear Analysis: Real World Applications B Fig. 2. Existence of 3,B 3. which is positive or negative depending on the sign of ᾱ = α K 2 α 2 /λ 20 K 2. It is noted that the value of B at = 0 is given by the following cubic equation: λ 0 λ 20 K 2 22 B3 + 2λ 0 λ 20 K 2 K 22 + λ 20 K 2 22 α B 2 + λ 0 λ 20 K K 2λ 2 K 2 K 22 α 2 + λ 20 α K 2 K 22 B + λ 20 K 2 2ᾱ = It is clear that one root of Eq is positive or negative depending on whether ᾱ is negative or positive, respectively. hus, from the above analysis we conclude that the two isoclines 3.20 and 3.2 intersect at a unique point 3,B 3 in the interior of first quadrant see Fig. 2, which is drawn for ᾱ > 0 provided conditions 3.3 and 3.4 are satisfied. Using these values of 3,B 3 in Eqs. 3.4 and 3.6, we get P3 = λ2 K 2 B3 λ 20 K 2 + K 22 B3 C3 = K 3 q λ B 3 2 α 4 K 2 + K 3 α 2, 3.23 λ 23 K 2 B 3 P 3 K 2 + K 22 B 3 λ α B 3 λ 22α 2 P which are positive under conditions 3.5 and 3.6, respectively. hus E 3 exists provided conditions are satisfied. Keeping in mind model 2., it is noted on physical consideration that the growth rate of bacteria in the absence of protozoa is more than when protozoa is present as protozoa predate bacteria. hus B2 >B 3. his result can also be proved from model 2. by using a comparison theorem 7. Existence of E 4 : for the equilibrium E 4,B,P,N,C, and B are obtained by solving the following algebraic equations: B = Q α 0 K K + K and λ K K 2 + K α K 2 λ30 λ 2 K 2 B λ 0 B K 2 + K 22 B l K 2 + K 22 B m = 0, 3.26 where l = λ 20 λ 30 + λ 3 K3 2 and m = λ 30α 2 K 3 α 3.
7 986 A.K. Misra et al. / Nonlinear Analysis: Real World Applications B Fig. 3. Existence of,b. Isocline 3.25 is the same as 3.7; hence, its behavior is also the same. From isocline 3.26, it is noted that under condition 3.8 B increases with in 0 B R B. It is also noted that B = 0at = K 2 α K / { 2m λ K K α + K } K 2 m K 2 l K 2 l which is positive or negative depending upon the sign of α = α K 2 /K 2 m/l. Further, for this isocline, the value of B at = 0 is given by the following cubic equation: λ 0 K22 2 B3 + 2λ 0 K 2 K 22 + K22 2 α B 2 + λ 0 K2 2 + K 2 λ 30 λ 2 K 2 mk K 2 K 22 α B l + K2 2 α = One root of Eq is positive or negative depending on the negative or positive value of α, respectively. hus, from the above discussion, we conclude that the two isoclines 3.25 and 3.26 intersect at a unique point,b in the interior of first quadrant see Fig. 3, drawn for α > 0 provided 3.7 and 3.8 are satisfied. Using these values of and B we get, P = λ30 λ 2 K 2 B l K 2 + K 22 B m, 3.28 N = λ 3 K 3 P α 3, 3.29 λ 30 C = K B q λ 2 α 4 K 2 + K λ K 2 B P 23 K 2 + K 22 B λ 34K 3 P N λ α B λ 22 α 2 P λ 33 α 3 N 3.30 which are positive under conditions hus E 4 exists, provided conditions are satisfied. On physical consideration from model 2. it is noted here that growth rate of protozoa in the absence of lower order aquatic population is greater than that in its presence as protozoa is predated by this aquatic population. Hence P 3 >P. his result can also be proved from model 2. by using a comparison theorem 7.
8 A.K. Misra et al. / Nonlinear Analysis: Real World Applications Remark.. Using Eqs and 3.26 we can check that under condition 3.8, d /dq>0 and db /d > 0. his gives db /dq = db /d d /dq > 0. Also, from Eq. 3.28, we note that dp dq = l λ 30 λ 2 K 2 K 2 db K 2 + K 22 B 2 dq > 0. Similarly from Eq. 3.29, we get dn /dq = λ 3 K 3 /λ 30 dp /dq > 0. Hence,B,P and N increase as Q increases. Also from Eq it can be seen that dc /dq<0, dc /dλ ij < 0 and dc /dq>0. he first inequality implies that as the rate of discharge of organic pollutants Q increases, the equilibrium level of dissolved oxygen C decreases. he second inequality suggests that C further decreases due to presence of various degradation processes. Also, it can be noted from the third inequality that C increases as the rate of input of dissolved oxygen q increases. 4. Stability analysis In our analysis, we assume that all the above equilibria exist. he local stability behavior of these equilibria has been studied and the result is stated in the following theorem: heorem. he equilibrium E i i =, 2 or 3 is unstable whenever E i+ exists. he equilibrium E 4 is locally asymptotically stable provided the following condition is satisfied: λ 0 K 2K 22 K 2 2 he proof is given in Appendix A. P > heorem 2. he equilibrium E 4 is nonlinearly stable in Ω provided the following two conditions are satisfied: λ 0 K 2K 22 P > 0, 4.2 K 2 2 α 0 λ 0 K 2K 22 λ K2 2 P K K 2 Q 2 K 2 + K > K 2 α 0 + K Q he proof is given in Appendix B. We note from the above that if λ 0 = 0, the above conditions are never satisfied. his shows that crowding coefficient of bacteria population density i.e. flaking off term stabilizes the system. he above theorems imply that the concentration of dissolved oxygen decreases as the discharge rate of organic pollutants increases. hey also show that if the cumulative rate of discharge of pollutants is very high, the concentration of dissolved oxygen may become negligible caused by various degradation processes. 5. A model for control of organic pollutants Here, we propose and analyze a nonlinear mathematical model to study the effect of control of discharge rate of organic pollutants on dissolved oxygen by using some effort 5. Let e be the effort applied to control the discharge of organic pollutants. It is assumed that the growth rate of effort applied to control the organic pollutants is proportional to the cumulative concentration of organic pollutants i.e. r, and that the rate of natural depletion of effort is proportional to e i.e. r 0 e. It is assumed further that the rate of decrease of organic pollutants in the water body is proportional to effort e i.e. γe. In the following model the other notations and assumptions for growth rates and depletion rates
9 988 A.K. Misra et al. / Nonlinear Analysis: Real World Applications of organic pollutants, bacteria, protozoa, lower order aquatic population and dissolved oxygen are the same as in model 2.. Keeping the above considerations in mind, the control dynamics of the system is governed by the following system of differential equations: d db dp dn dc = Q α 0 K B K 2 + K γ e, = λ K B K 2 + K α B λ 0 B 2 K 2 BP K 2 + K 22 B, = λ 2 K 2 BP K 2 + K 22 B α 2P λ 20 P 2 K 3 PN, = λ 3 K 3 PN α 3 N λ 30 N 2, = q α 4 C λ 34 K 3 PN λ α B λ 22 α 2 P λ 33 α 3 N λ 2 K B K 2 + K λ 23 K 2 BP K 2 + K 22 B, de = r r 0 e, 5. where 0>0, B0>0, P0>0, N0>0, C0>0, e0>0. In 5. γ is the control rate of cumulative density of organic pollutants, r is the growth rate of effort and r 0 is its depletion rate. As before we can find that the set Ω c, defined below, is a region of attraction for model 5.. { Ω c = 0 Q, 0 B R B, 0 P R P, 0 N R N, 0 C q, 0 e r } Q, 5.2 α 0 α 4 r 0 α 0 where all the initial values of the variables belong to set Ω c. 5.. Equilibrium analysis Model 5. has four equilibria in Ω c as described below the proof of existence of each of these equilibria is similar as discussed for model 2.. i Ẽ r0 Q r 0 α 0 +γr, 0, 0, 0, q r α 4, Q r 0 α 0 +γr. ii Ẽ 2 2, B 2, 0, 0, C 2, ẽ 2 exists provided the following are satisfied: and λ K K α Q K 2 α α 0 + γr /r 0 >0, 5.3 q λ 2 K 2 B 2 K 2 + K 2 λ α B 2 > 0, 5.4 where 2, B 2, C 2 and ẽ 2 are given by the following equations: Q α 0 + γr 2 K 2 B 2 = 0, r 0 K 2 + K 2 λ K 2 B 2 K 2 + K 2 α λ 0 B 2 = 0,
10 A.K. Misra et al. / Nonlinear Analysis: Real World Applications C 2 = K 2 B 2 q λ 2 λ α B 2, α 4 K 2 + K 2 ẽ = r r 0 2. iii Ẽ 3 3, B 3, P 3, 0, C 3, ẽ 3 exists provided the following conditions are satisfied: λ 20 λ K K α Q K 2 α α 0 + γr + K 2α 2 K Q + K 2 α 0 + γr > 0, 5.5 r 0 K 2 r 0 λ 0 + K 2 λ 20 K 2 + K 22 R B 2 λ 2 K 2 K 2 K 2 + K 22 R B K 22 λ2 K 2 R B α 2 > 0, 5.6 K 2 + K 22 R B λ 2 K 2 B 3 K 2 + K 22 B 3 α 2 > 0, 5.7 K 3 B 3 K 2 B 3 P 3 q λ 2 λ 23 λ α B 3 λ 22 α 2 P 3 > 0, 5.8 K 2 + K 3 K 2 + K 22 B 3 where 3, B 3, P 3, C 3 and ẽ 3 are given by the following equations: Q α 0 + γr 3 K 3 B 3 = 0, r 0 K 2 + K 3 λ K 3 K 2 + K 3 α = λ 0 B + P 3 = λ 2 K 2 B 3 α 2, λ 20 K 2 + K 22 B 3 K 2 λ 20 K 2 + K 22 B 3 λ 2 K 2 B 3 K 2 + K 22 B 3 α 2 C 3 = K 3 B 3 K 2 B 3 P 3 q λ 2 λ 23 λ α B 3 λ 22 α 2 P 3, α 4 K 2 + K 3 K 2 + K 22 B 3 ẽ 3 = r r 0 3., iv Ẽ 4, B, P,Ñ, C,ẽ exists, provided the following conditions are satisfied: λ K K α Q K 2 α α 0 + γr + K 2m K Q + K 2 α 0 + γr > 0, 5.9 r 0 K 2 l r 0 λ 0 + K 2 λ30 λ 2 K 2 K 2 K 2 + K 22 R B 2 l K 2 + K 22 R B K λ30 λ 2 K 2 R B 22 K 2 + K 22 R B m > 0, 5.0 λ 30 λ 2 K 2 B m>0, 5. K 2 + K 22 B
11 990 A.K. Misra et al. / Nonlinear Analysis: Real World Applications λ 3 K 3 P α 3 > 0, 5.2 q λ 2K B K 2 + K λ 23K 2 B P K 2 + K 22 B λ 34K 3 P Ñ λ α B λ 22 α 2 P λ 33 α 3 Ñ>0, 5.3 where l = λ 20 λ 30 + λ 3 K 2 3, m = λ 30α 2 K 3 α 3 and, B, P,Ñ, C,ẽ are given by the following equations: Q α 0 + γr r 0 K B K 2 + K λ K K 2 + K α = λ 0 B + P = l = 0, 5.4 K 2 lk 2 + K 22 B λ 30 λ 2 K 2 B K 2 + K 22 B m, 5.5 λ 30 λ 2 K 2 B K 2 + K 22 B m, 5.6 Ñ = λ 30 λ 3 K 3 P α 3, 5.7 C = q α 4 λ 2K B K 2 + K λ 23K 2 B P K 2 + K 22 B λ 34K 3 P Ñ λ α B λ 22 α 2 P λ 33 α 3 Ñ, 5.8 ẽ = r r In this case also it can be checked that as the rate of control of organic pollutants γ increases, the equilibrium level of the dissolved oxygen concentration C increases. As in the previous case, for model 5. also we can easily show that Ẽ i i =, 2or3 is unstable whenever Ẽ i+ exists. he local stability behavior of the equilibrium Ẽ 4, B, P,Ñ, C,ẽ is stated in the following theorem. heorem 3. he equilibrium Ẽ 4 is locally asymptotically stable if λ 0 K 2K 22 K 2 2 he proof is given in Appendix C. P> heorem 4. he equilibrium Ẽ 4 is nonlinearly stable in Ω c if the following conditions are satisfied: λ 0 K 2K 22 K 2 2 P>0, 5.2 α 0 λ 2 λ 0 K 2K 22 K K 2 Q P K2 2 > K 2 + K K 2 α 0 + K Q In this case also, we may note the stabilizing effect of λ 0. Further, these theorems imply that if an appropriate effort is applied to control organic pollutants, the equilibrium level of dissolved oxygen can be maintained at a desired level. he proof is given in Appendix D.
12 A.K. Misra et al. / Nonlinear Analysis: Real World Applications *, B * B Fig. 4. Global stability in B plane. 6. Numerical example o check the feasibility of our analysis regarding the existence of E 4 and corresponding stability conditions, we conduct some numerical computation by choosing the following values of the parameters in model 2.: Q = 3.0, α 0 =.0, K =.0, K 2 = 0., K =.0, λ =.0, α = 0.5, λ 0 = 0., K 2 =.0, K 2 = 2.0, K 22 =.0, λ 2 =.0, α 2 = 0.02, λ 20 = 0.25, K 3 =.0, λ 3 =.0, α 3 = 0.2, λ 30 = 0.25, q = 0.0, α 4 = 2.0, λ 2 = λ 23 = λ 34 = λ = λ 22 = λ 33 = It is found that under the above set of parameters, conditions for the existence of interior equilibrium E 4,B,P,N,C are satisfied and E 4 is given by = , B = , P = , N = , C = he eigenvalues of the Jacobian matrix M corresponding to this equilibrium E 4 are 2, i, i, i, i, which are either negative or have negative real parts. Hence E 4 is locally stable. It is pointed out here that for the above set of parameters, the conditions for local stability 4. and nonlinear stability 4.2 and 4.3 are also satisfied. Further, for the above set of parameters model 2. computer generated graphs of verses B and verses C are shown in Figs. 4 and 5, which indicates the global stability of,b in the B-plane and,c in C-plane. 7. Conclusion In this paper, we have proposed and analyzed a nonlinear mathematical model for the depletion of dissolved oxygen due to discharge of organic pollutants in a water body by considering biodegradation and biochemical processes in
13 992 A.K. Misra et al. / Nonlinear Analysis: Real World Applications C 4 *, C * Fig. 5. Global stability in C plane. the food chain involving bacteria, protozoa and an aquatic population for example a fish population in a water body. It has been shown that if the cumulative rate of discharge of pollutants increases, the equilibrium level of dissolved oxygen decreases, the amount of which depends on the rate of growth of oxygen as well as on various biodegradation and biochemical processes in the water body. It has also been shown that under certain conditions, if the cumulative rate of introduction of water pollutants is too high, the equilibrium concentration of the dissolved oxygen may become negligibly small, threatening the survival of biological species in a water body. A model to control the cumulative discharge of organic pollutants in a water body is also proposed and analyzed. It has been shown that the equilibrium concentration of the dissolved oxygen can be maintained at a desired level by using an appropriate effort to control the cumulative discharge of organic pollutants. Appendix A. Local stability analysis he general variational matrix M for system 2. is given as follows: K f, B K 2 + K λ K K 2 B K 2 B K M = 2 + K 2 f 2,B,P 0 0 K 2 + K 22 B λ 2 K 2 K 2 P, 0 K 2 + K 22 B 2 f 3 B,P,N K 3 P λ 3 K 3 N f 4 P, N 0 g, B g 2,B,P g 3 B, N g 4 P α 4 where f, B = α 0 + f 2,B,P= f 3 B,P,N= K K 2 B K 2 + K 2, λ K K 2 + K α 2λ 0 B K 2 K 2 P K 2 + K 22 B 2, λ 2 K 2 B K 2 + K 22 B α 2 2λ 20 P K 3 N, f 4 P, N = λ 3 K 3 P α 3 2λ 30 N,
14 g, B = g 3 B, N = A.K. Misra et al. / Nonlinear Analysis: Real World Applications λ 2K K 2 B K 2 + K 2, g 2,B,P= λ 2K K 2 + K + λ 23K 2 K 2 P K 2 + K 22 B 2 + λ α, λ 23K 2 B K 2 + K 22 B + λ 22α 2 + λ 34 K 3 N, g 4 P = λ 34 K 3 P + λ 33 α 3. A. We use eigenvalue method to study the local stability behavior of equilibria E i i =, 2, 3. Let M i be the matrix obtained from M after substituting for E i. For the equilibrium E, we note that one of the eigenvalues of M is λ K K α Q K 2 α 0 α /K 2 α 0 +K Q, which is positive whenever E 2 exists condition 3.. Hence, E is unstable whenever E 2 exists. For the equilibrium E 2 2,B 2, 0, 0,C 2, we note that one of the eigenvalues of M 2 is λ 2 K 2 B 2 /K 2+K 22 B 2 α 2, which is positive whenever E 3 exists as B 3 <B 2. hus, E 2 is unstable whenever E 3 exists. For the equilibrium E 3 3,B 3,P 3, 0,C 3, we note that one of the eigenvalues of M 3 is λ 3 K 3 P 3 α 3, which is positive whenever E 4 exists as P <P 3. his shows that E 3 is unstable whenever E 4 exists. Since in our study E 4 is the most interesting equilibrium from ecological point of view and its behavior cannot be described in a simple manner from M 4, we discuss its behavior by using Lyapunov s method. Proof of heorem. We linearize system 2. by using the transformations = + τ, B= B + b, P = P + p, N = N + n, C = C + c. A.2 Now we consider the following positive definite function: V = 2 τ 2 + m B b2 + m 2 P p2 + m 3 N n2 + m 4 c 2, A.3 and use the linearized system of 2. to get dv = f,b τ 2 m f 5 B,P b 2, m 2 λ 20 p 2 m 3 λ 30 n 2 m 4 α 4 c 2, m 4 g,b τc m 4 g 2,B,P bc m 4 g 3 B,N pc m 4 g 4 P nc, A.4 where the functions f and g i i =, 2, 3, 4 are given by A., and f 5 B, P = λ 0 K 2 K 22 P /K 2 + K 22 B 2. Now choosing m = K 2 + K /λ K 2, m 2 = m K 2 + K 22 B /λ 2 K 2, m 3 = m 2 /λ 3 we note that the coefficient of b 2 is negative under condition 4.. hus dv/ can be made negative definite by appropriately choosing a positive value for m 4 following usual methods of nonlinear analysis 4,7,6,7. Appendix B. Proof of heorem 2 o prove this theorem we consider the following positive definite function: V = 2 /2 + m B B B lnb/b + m 2 P P P lnp /P + m 3 N N N lnn/n + m 4 C C 2 /2. B.
15 994 A.K. Misra et al. / Nonlinear Analysis: Real World Applications Now using system 2. we get, dv K K 2 B 2 = K 2 + K K 2 + K m K 2 K 22 P K2 2 K 2 + K 22 BK 2 + K 22 B B B 2 α 0 2 K 2 K 22 m λ 0 K 2 + K 22 B 2 P B B 2 m 2 λ 20 P P 2 m 3 λ 30 N N 2 m 4 α 4 C C 2 + B B K K 2 K 2 + K K 2 + K m 4 C C λ 2 K K 2 B K 2 + K K 2 + K m 4 B B C C λ2 K K 2 + K + λ 23 K 2 K 2 P K 2 + K 22 BK 2 + K 22 B + λ α m 4 P P C C λ23 K 2 B K 2 + K 22 B + λ 34K 3 N + λ 22 α 2 m 4 N N C C λ 34 K 3 P + λ 33 α 3. Note that K2 2 K 2 + K 22 BK 2 + K 22 B > 0 inside the region of attraction Ω. Now choosing m = /λ,m 2 = m K 2 + K 22 B /λ 2 K 2, m 3 = m 2 /λ 3, we note that dv/ can be made negative definite inside Ω if λ 0 K 2K 22 P > 0, K 2 2 α 0 λ 0 K 2K 22 λ K2 2 P K K 2 Q 2 K 2 + K > 0 B.4 K 2 α 0 + K Q by making an appropriate choice for m 4 as pointed out earlier. B.2 B.3 Appendix C. Proof of heorem 3 We linearize system 5. by using the following transformations: = + τ, B = B + b, P = P + p, N = Ñ + n, C = C + c, e =ẽ + e. Now using the following positive definite function: V = τ 2 + m 2 B b2 + m 2 P p2 + m 3 Ñ n2 + m 4 c 2 + m 5 e 2, C. and use the linearized system of 5. to get dv = f, Bτ 2 m f 5 B, Pb 2 m 2 λ 20 p 2 m 3 λ 30 n 2 m 4 α 4 c 2 γr 0 e 2 r m 4 g, Bτc m 4 g 2, B, Pbc m 4 g 3 B,Ñpc m 4 g 4 P nc and m s chosen as positive constants in the following manner: m = K 2 + K /λ K, m 2 = m K 2 + K 22 B/λ 2 K 2, m 3 = m 2 /λ 3,m 5 = γ/r. C.2
16 If the following condition is satisfied: λ 0 K 2K 22 K 2 2 P>0. A.K. Misra et al. / Nonlinear Analysis: Real World Applications then dv/ can be made negative definite by appropriate choice of m 4 > 0. Hence Ẽ 4 is locally stable under above condition. Appendix D. Proof of heorem 4 We prove this theorem by using the following positive definite function: V = m B B B ln B B + m 2 P P P ln P P + m 3 N Ñ Ñ ln Ñ + N 2 m 4C C m 5e ẽ 2, D. where m,m 2,m 3,m 4,m 5 are some positive constants to be chosen appropriately. Using system 5. with the constants m,m 2,m 3 and m 5 chosen as follows: m = K 2 + K /λ K, m 2 = m K 2 + K 22 B/λ 2 K 2, m 3 = m 2 /λ 3,m 5 = γ/r,weget dv K K 2 B 2 = K 2 + K K 2 + K m K 2 K 22 P K2 2 B B 2 m 5 r 0 e ẽ 2 K 2 + K 22 BK 2 + K 22 B α 0 2 m λ 0 K 2K 22 P K2 2 B B 2 m 2 λ 20 P P 2 m 3 λ 30 N Ñ 2 m 4 α 4 C C 2 K K 2 + B B K 2 + K K 2 + K λ 2 K K 2 B m 4 C C K 2 + K K 2 + K λ 2 K m 4 B BC C K 2 + K + λ 23 K 2 K 2 P K 2 + K 22 BK 2 + K 22 B + λ α m 4 P P λ 23 K 2 B C C K 2 + K 22 B + λ 34K 3 Ñ + λ 22 α 2 m 4 N ÑC Cλ 34 K 3 P + λ 33 α 3. Now dv/ can be made negative definite in Ω c provided the following conditions as given in the theorem are satisfied: λ 0 K 2K 2 K 2 2 P>0 and α 0 λ 0 K 2K 22 K K 2 P λ K2 2 K 2 + K for a suitable choice of m 4 > 0. 2 Q > 0 K 2 α 0 + K Q D.2
17 996 A.K. Misra et al. / Nonlinear Analysis: Real World Applications References M.B. Beck, P.C. Young, A dynamic model for DO BOD relationship in a non-tidal stream, Water Res S. Busenberg, S. Kishore Kumar, P. Austin, G. Wake, he dynamics of a plankton nutrient interaction, Bull. Math. Biol W.E. Dobbins, BOD and oxygen relationship in streams, J. San. Eng. Div., Proc. ASCE B. Dubey, P. Chandra, P. Sinha, A model for fishery resource with reserve area, Nonlinear Anal.: Real World Prob H.I. Freedman, J.W.H. So, Global stability and persistence of simple food chains, Math. Biosci W.E. Gates, J.. Marler, J.D. Westfield, he application of bacterial process kinetics in stream simulation and stream analysis, Water Res J.K. Hale, Ordinary Differential Equations, Wiley-Interscience, New York, G. Hallam, Structural sensitivity of grazing formulations in nutrient controlled plankton models, J. Math. Biol M. Jayweera,. Asaeda, Modeling of biomanipulation in shallow, eutrophic lakes: an application to Lake Bleiswijkse Zoom, the Netherlands, Ecol. Model S.E. Jorgenson, Fundamental of Ecological Modeling, Elsevier Science Publishers B.V., Amsterdam, Netherlands, 988. R.A. Kelly, Conceptual ecological model of the Delaware estuary, in: B.C. Patten Ed., Systems Analysis and Simulation in Ecology, vol. IV, Academic Press, New York, J. LaSalle, S. Lefschetz, Stability by Liopunov s Direct Method with Applications, Academic Press, New York, A.J. McDonel, Oxygen budgets in macrophyte impacted streams, Water Res D.J. O Connor, he temporal and spatial distribution of dissolved oxygen in streams, Water Resour. Res S. Rinaldi, R. Soncini-sessa, H. Stehfest, H. amura, Modeling and Control of River Quality, McGraw-Hill, UK, J.B. Shukla, A.K. Agrawal, B. Dubey, P. Sinha, Existence and survival of two competing species in a polluted environment, J. Biol. Syst J.B. Shukla, B. Dubey, H.I. Freedman, Effect of changing habitat on survival of species, Ecol. Model J.H. Steele, he Structure of Marine Ecosystems, Harvard University Press, Cambridge, J.H. Steele, B.W. Frost, he structure of plankton communities, Philos. rans. R. Soc. London B M. Straskraba, A. Gnauck, Freshwater ecosystems modeling and simulation. Developments in Environmental Modeling, vol. 8, Elsevier, Amsterdam, M.W. Streeter, E.E. Phelps, A study of the pollution and natural purification of the Ohio rivers, U.S. Public Health Service Bill, 925, pp R.V. homann, Bioaccumulation model of organic chemical distribution in aquatic food chains, Environ. Sci. echnol F.-L. Xu, S. ao, R.W. Dawson, P.-G. Li, J. Cao, Lake ecosystem health assessment indicators and methods, Water Res
Mathematical Modeling and Analysis of Eutrophication of Water Bodies Caused by Nutrients
Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 511 524 Mathematical Modeling and Analysis of Eutrophication of Water Bodies Caused by Nutrients A. K. Misra Department of Mathematics,
More informationStability Analysis of Plankton Ecosystem Model. Affected by Oxygen Deficit
Applied Mathematical Sciences Vol 9 2015 no 81 4043-4052 HIKARI Ltd wwwm-hikaricom http://dxdoiorg/1012988/ams201553255 Stability Analysis of Plankton Ecosystem Model Affected by Oxygen Deficit Yuriska
More informationHarvesting Model for Fishery Resource with Reserve Area and Modified Effort Function
Malaya J. Mat. 4(2)(2016) 255 262 Harvesting Model for Fishery Resource with Reserve Area and Modified Effort Function Bhanu Gupta and Amit Sharma P.G. Department of Mathematics, JC DAV College, Dasuya
More informationModelling the spread of bacterial infectious disease with environmental effect in a logistically growing human population
Nonlinear Analysis: Real World Applications 7 2006) 341 363 www.elsevier.com/locate/na Modelling the spread of bacterial infectious disease with environmental effect in a logistically growing human population
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 informationStability analysis of a prey-predator model with a reserved area
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 4, 5(3):93-3 ISSN: 976-86 CODEN (USA): AASRFC Stability analysis of a prey-predator model with a reserved area Neelima
More informationPrey Predator Model with Asymptotic. Non-Homogeneous Predation
Int. J. Contemp. Math. Sciences, Vol. 7, 0, no. 4, 975-987 Prey Predator Model with Asymptotic Non-Homogeneous Predation Harsha Mehta *, Neetu Trivedi **, Bijendra Singh * and B. K. Joshi *** *School of
More informationA mathematical model to study the simultaneous effect of pollutants emitted from wood and non-wood based industries on forest resources
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 12 2016 No. 1, pp. 34-47 A mathematical model to study the simultaneous effect of pollutants emitted from wood and non-wood based
More informationHARVESTING IN A TWO-PREY ONE-PREDATOR FISHERY: A BIOECONOMIC MODEL
ANZIAM J. 452004), 443 456 HARVESTING IN A TWO-PREY ONE-PREDATOR FISHERY: A BIOECONOMIC MODEL T. K. KAR 1 and K. S. CHAUDHURI 2 Received 22 June, 2001; revised 20 September, 2002) Abstract A multispecies
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 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 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 informationModeling the effect of deforestation caused by human population pressure on wildlife species
Nonlinear Analysis: Modelling and Control, Vol. 23, No. 3, 33 32 ISSN 1392-5113 https://doi.org/1.15388/na.218.3.2 Modeling the effect of deforestation caused by human population pressure on wildlife species
More informationECOLOGICAL MODELLING
ECOLOGICAL MODELLING II. CONCEPTUAL AND MODELING THEORY 1 1. Models - Mankind has always used models as tools to solve problems as they give a simplified picture of reality - The model never contain all
More informationAnalysis of Plankton Populations
Analysis of Plankton Populations Sean Li September 21, 25 Abstract I use numerical and analytical techniques to study the dynamics of various models of plankton food webs with and without resource fluctuations.
More informationGlobal Properties for Virus Dynamics Model with Beddington-DeAngelis Functional Response
Global Properties for Virus Dynamics Model with Beddington-DeAngelis Functional Response Gang Huang 1,2, Wanbiao Ma 2, Yasuhiro Takeuchi 1 1,Graduate School of Science and Technology, Shizuoka University,
More informationSENSITIVITY ANALYSIS OF STREETER-PHELPS MODELS
SENSITIVITY ANALYSIS OF STREETER-PHELPS MODELS S. RINALDI~AND R. SONCINI-SESSA* FEBRUARY 1977 * Centro Teoria dei Sistemi. C.N.R.. Via Ponzio. 3415, Milano, Italy. This work was supported jointly by IIASA
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 informationAnalysis of a Prey-Predator System with Modified Transmission Function
Research Paper American Journal of Engineering Research (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-3, Issue-9, pp-194-202 www.ajer.org Open Access Analysis of a Prey-Predator System with Modified
More informationStability Analysis of a SIS Epidemic Model with Standard Incidence
tability Analysis of a I Epidemic Model with tandard Incidence Cruz Vargas-De-León Received 19 April 2011; Accepted 19 Octuber 2011 leoncruz82@yahoo.com.mx Abstract In this paper, we study the global properties
More informationApplied Mathematics Letters. Stationary distribution, ergodicity and extinction of a stochastic generalized logistic system
Applied Mathematics Letters 5 (1) 198 1985 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Stationary distribution, ergodicity
More informationPERSISTENCE AND EXTINCTION OF SINGLE POPULATION IN A POLLUTED ENVIRONMENT
Electronic Journal o Dierential Equations, Vol 2004(2004), No 108, pp 1 5 ISSN: 1072-6691 URL: http://ejdemathtxstateedu or http://ejdemathuntedu tp ejdemathtxstateedu (login: tp) PERSISTENCE AND EXTINCTION
More informationExistence and uniqueness for optimal control of Oxygen absorption in aquatic system
Existence and uniqueness for optimal control of Oxygen absorption in aquatic system F.B. Agusto and O. M. Bamigbola Abstract. In this paper, we consider the mathematical formulation and analysis of an
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 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 informationPersistence and global stability in discrete models of Lotka Volterra type
J. Math. Anal. Appl. 330 2007 24 33 www.elsevier.com/locate/jmaa Persistence global stability in discrete models of Lotka Volterra type Yoshiaki Muroya 1 Department of Mathematical Sciences, Waseda University,
More informationBehaviour of simple population models under ecological processes
J. Biosci., Vol. 19, Number 2, June 1994, pp 247 254. Printed in India. Behaviour of simple population models under ecological processes SOMDATTA SINHA* and S PARTHASARATHY Centre for Cellular and Molecular
More informationROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 6, Number 1, Winter 1997 ROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM J. CHATTOPADHYAY, E. BERETTA AND F. SOLIMANO ABSTRACT. The present paper deals with
More informationGENERAL ECOLOGY STUDY NOTES
1.0 INTRODUCTION GENERAL ECOLOGY STUDY NOTES A community is made up of populations of different organisms living together in a unit environment. The manner in which these organisms relate together for
More informationGHS S.4 BIOLOGY TEST 2 APRIL Answer all the questions in Section A and B. in the spaces provided
GHS S.4 BIOLOGY TEST 2 APRIL 2016 TIME: 1 HOUR Instructions: Answer all the questions in Section A and B. in the spaces provided ANSERS TO SECTION A 1 6 11 16 21 26 2 7 12 17 22 27 3 8 13 18 23 28 4 9
More information4 Insect outbreak model
4 Insect outbreak model In this lecture I will put to a good use all the mathematical machinery we discussed so far. Consider an insect population, which is subject to predation by birds. It is a very
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 informationVariational iteration method for solving multispecies Lotka Volterra equations
Computers and Mathematics with Applications 54 27 93 99 www.elsevier.com/locate/camwa Variational iteration method for solving multispecies Lotka Volterra equations B. Batiha, M.S.M. Noorani, I. Hashim
More informationWater-vegetation interactions in a sea of water: Plankton in mesoscale turbulence
Water-vegetation interactions in a sea of water: Plankton in mesoscale turbulence A. Provenzale ISAC-CNR, Torino and CIMA, Savona, Italy Work done with: Annalisa Bracco, Claudia Pasquero Adrian Martin,
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationDroop models of nutrient-plankton interaction with intratrophic predation
Droop models of nutrient-plankton interaction with intratrophic predation S. R.-J. Jang 1, J. Baglama 2 1. Department of Mathematics, University of Louisiana at Lafayette, Lafayette, LA 754-11 2. Department
More informationDYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS. Wei Feng
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 7 pp. 36 37 DYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS Wei Feng Mathematics and Statistics Department
More informationAnalysis of bacterial population growth using extended logistic Growth model with distributed delay. Abstract INTRODUCTION
Analysis of bacterial population growth using extended logistic Growth model with distributed delay Tahani Ali Omer Department of Mathematics and Statistics University of Missouri-ansas City ansas City,
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 informationThe construction and use of a phase diagram to investigate the properties of a dynamic system
File: Phase.doc/pdf The construction and use of a phase diagram to investigate the properties of a dynamic system 1. Introduction Many economic or biological-economic models can be represented as dynamic
More informationThe effect of emigration and immigration on the dynamics of a discrete-generation population
J. Biosci., Vol. 20. Number 3, June 1995, pp 397 407. Printed in India. The effect of emigration and immigration on the dynamics of a discrete-generation population G D RUXTON Biomathematics and Statistics
More informationSuccess Criteria Life on Earth - National 5
Success Criteria Life on Earth - National 5 Colour the box at the side of each objective: RED I don t know much about this or am confused by it. AMBER I know a bit about this but do not feel I know it
More information9693 MARINE SCIENCE. Mark schemes should be read in conjunction with the question paper and the Principal Examiner Report for Teachers.
CAMBRIDGE INTERNATIONAL EXAMINATIONS GCE Advanced Subsidiary Level and GCE Advanced Level www.xtremepapers.com MARK SCHEME for the May/June 01 series 969 MARINE SCIENCE 969/01 Paper 1 (AS Structured Questions),
More informationHomework 2. Due Friday, July We studied the logistic equation in class as a model of population growth. It is given by dn dt = rn 1 N
Problem 1 (10 points) Homework Due Friday, July 7 017 We studied the logistic equation in class as a model of population growth. It is given by dn dt = rn 1 N, (1) K with N(0) = N 0. (a) Make the substitutions
More informationBIOGEOCHEMICAL CYCLES
BIOGEOCHEMICAL CYCLES BASICS Biogeochemical Cycle: The complete path a chemical takes through the four major components, or reservoirs, of Earth s system (atmosphere, lithosphere, hydrosphere and biosphere)
More informationTuring and Non-Turing patterns in diffusive plankton model
Computational Ecology and Software 05 5(): 6-7 Article Turing and Non-Turing patterns in diffusive plankton model N. K. Thakur R. Gupta R. K. Upadhyay National Institute of Technology Raipur India Indian
More informationUNIT 5. ECOSYSTEMS. Biocenosis Biotope Biotic factors Abiotic factors
UNIT 5. ECOSYSTEMS 1. Define: ecosystem, biocenosis, biotope, abiotic factor, biotic factor 2. Complete using this word: ecosphere, biosphere, ecology, ecosystem a) The is all of the living thing on Earth.
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 informationStability and Hopf bifurcation of a ratio-dependent predator prey model with time delay and stage structure
Electronic Journal of Qualitative Theory of Differential Equations 6, No 99, 3; doi: 43/ejqtde699 http://wwwmathu-szegedhu/ejqtde/ Stability and Hopf bifurcation of a ratio-dependent predator prey model
More informationSome New Results on Lyapunov-Type Diagonal Stability
Some New Results on Lyapunov-Type Diagonal Stability Mehmet Gumus (Joint work with Dr. Jianhong Xu) Department of Mathematics Southern Illinois University Carbondale 12/01/2016 mgumus@siu.edu (SIUC) Lyapunov-Type
More informationChapter 10. Marine Ecology
Chapter 10 Marine Ecology Copyright 2016 McGraw-Hill Education. All rights reserved. No reproduction or distribution without the prior written consent of McGraw-Hill Education. Marine Ecology Ecology is
More informationMath 2930 Worksheet Introduction to Differential Equations. What is a Differential Equation and what are Solutions?
Math 2930 Worksheet Introduction to Differential Equations Week 1 January 25, 2019 What is a Differential Equation and what are Solutions? A differential equation is an equation that relates an unknown
More informationApplied Mathematics Letters
Applied Mathematics Letters 24 (211) 219 223 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Laplace transform and fractional differential
More informationMathematical Modelling of Plankton-Oxygen Dynamics under the Climate Change. Yadigar Sekerci & Sergei Petrovskii 1
Mathematical Modelling of Plankton-Oxygen Dynamics under the Climate Change Yadigar Sekerci & Sergei Petrovskii Department of Mathematics, University of Leicester, University Road, Leicester LE 7RH, U.K.
More informationRosenzweig-MacArthur Model. Considering the Function that Protects a Fixed. Amount of Prey for Population Dynamics
Contemporary Engineering Sciences, Vol. 11, 18, no. 4, 1195-15 HIKAI Ltd, www.m-hikari.com https://doi.org/1.1988/ces.18.8395 osenzweig-macarthur Model Considering the Function that Protects a Fixed Amount
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 informationNonlinear Dynamic Analysis of an Entermorpha prolifera Population Model
The Fourth International Conference on Computational Systems Biology (ISB2010) Suzhou, China, September 9 11, 2010 Copyright 2010 ORSC & APORC, pp. 29 34 Nonlinear Dynamic Analysis of an Entermorpha prolifera
More informationAvailable online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 19 (2016 ) 11 18 IUTAM Symposium Analytical Methods in Nonlinear Dynamics A model of evolutionary dynamics with quasiperiodic forcing
More informationInternational Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: Vol.8, No.7, pp , 2015
International Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: 0974-4304 Vol.8, No.7, pp 99-, 05 Lotka-Volterra Two-Species Mutualistic Biology Models and Their Ecological Monitoring Sundarapandian
More informationAnalysis of a predator prey model with modified Leslie Gower and Holling-type II schemes with time delay
Nonlinear Analysis: Real World Applications 7 6 4 8 www.elsevier.com/locate/na Analysis of a predator prey model with modified Leslie Gower Holling-type II schemes with time delay A.F. Nindjin a, M.A.
More informationImpulsive Stabilization and Application to a Population Growth Model*
Nonlinear Dynamics and Systems Theory, 2(2) (2002) 173 184 Impulsive Stabilization and Application to a Population Growth Model* Xinzhi Liu 1 and Xuemin Shen 2 1 Department of Applied Mathematics, University
More information4. Ecology and Population Biology
4. Ecology and Population Biology 4.1 Ecology and The Energy Cycle 4.2 Ecological Cycles 4.3 Population Growth and Models 4.4 Population Growth and Limiting Factors 4.5 Community Structure and Biogeography
More informationGLOBAL DYNAMICS OF A MATHEMATICAL MODEL OF TUBERCULOSIS
CANADIAN APPIED MATHEMATICS QUARTERY Volume 13, Number 4, Winter 2005 GOBA DYNAMICS OF A MATHEMATICA MODE OF TUBERCUOSIS HONGBIN GUO ABSTRACT. Mathematical analysis is carried out for a mathematical model
More informationStability Analysis of Predator- Prey Models via the Liapunov Method
Stability Analysis of Predator- Prey Models via the Liapunov Method Gatto, M. and Rinaldi, S. IIASA Research Memorandum October 1975 Gatto, M. and Rinaldi, S. (1975) Stability Analysis of Predator-Prey
More informationPOINT SOURCES OF POLLUTION: LOCAL EFFECTS AND IT S CONTROL Vol. II - Contaminant Fate and Transport Process - Xi Yang and Gang Yu
CONTAMINANT FATE AND TRANSPORT PROCESS Xi Yang and Gang Yu Department of Environmental Sciences and Engineering, Tsinghua University, Beijing, P. R. China Keywords: chemical fate, physical transport, chemical
More informationCommunities Structure and Dynamics
Communities Structure and Dynamics (Outline) 1. Community & niche. 2. Inter-specific interactions with examples. 3. The trophic structure of a community 4. Food chain: primary, secondary, tertiary, and
More informationStability and optimal harvesting in a stage structure predator-prey switching model
Stability and optimal harvesting in a stage structure predator-prey switching model Q. J. A. Khan and Lakdere Benkherouf Department of Mathematics and Statistics, College of Science, Sultan Qaboos University,
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 informationNON-CONSTANT POSITIVE STEADY STATES FOR A STRONGLY COUPLED NONLINEAR REACTION-DIFFUSION SYSTEM ARISING IN POPULATION DYNAMICS
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 45, Number 4, 2015 NON-CONSTANT POSITIVE STEADY STATES FOR A STRONGLY COUPLED NONLINEAR REACTION-DIFFUSION SYSTEM ARISING IN POPULATION DYNAMICS ZIJUAN WEN
More information10/6/ th Grade Ecology and the Environment. Chapter 2: Ecosystems and Biomes
7 th Grade Ecology and the Environment Chapter 2: Ecosystems and Biomes Lesson 1 (Energy Flow in Ecosystems) Each organism in an ecosystem fills an energy role. Producer an organism that can make its own
More informationMODELLING AND ANALYSIS OF THE SPREAD OF MALARIA: ENVIRONMENTAL AND ECOLOGICAL EFFECTS
Journal of Biological Systems, Vol. 13, No. 1 (2005) 1 11 c World Scientific Publishing Company MODELLING AND ANALYSIS OF THE SPREAD OF MALARIA: ENVIRONMENTAL AND ECOLOGICAL EFFECTS SHIKHA SINGH,J.B.SHUKLA
More informationLiving Things and the Environment
Unit 21.1 Living Things and the Environment Section 21.1 Organisms obtain food, water, shelter, and other things it needs to live, grow, and reproduce from its environment. An environment that provides
More informationAbsolute value equations
Linear Algebra and its Applications 419 (2006) 359 367 www.elsevier.com/locate/laa Absolute value equations O.L. Mangasarian, R.R. Meyer Computer Sciences Department, University of Wisconsin, 1210 West
More informationOcean facts continued
Ocean Facts A dynamic system in which many chemical and physical changes take place Formed over millions of years as precipitation filled low areas on Earth called basins and now covers 70% of the Earth
More informationYour web browser (Safari 7) is out of date. For more security, comfort and the best experience on this site: Update your browser Ignore
Your web browser (Safari 7) is out of date. For more security, comfort and the best experience on this site: Update your browser Ignore PL ANKTO N REVEALED A critical component of life on Earth For the
More informationCommunities Structure and Dynamics
Communities Structure and Dynamics (Outline) 1. Community & niche. 2. Inter-specific interactions with examples. 3. The trophic structure of a community 4. Food chain: primary, secondary, tertiary, and
More informationOBSERVABILITY AND OBSERVERS IN A FOOD WEB
OBSERVABILITY AND OBSERVERS IN A FOOD WEB I. LÓPEZ, M. GÁMEZ, AND S. MOLNÁR Abstract. The problem of the possibility to recover the time-dependent state of a whole population system out of the observation
More informationA 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 informationHopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd
Hopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd College of William and Mary Williamsburg, Virginia 23187 Mathematical Applications in Ecology and Evolution Workshop
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 informationCommunities Structure and Dynamics
Communities Structure and Dynamics (Outline) 1. Community & niche. 2. Inter-specific interactions with examples. 3. The trophic structure of a community 4. Food chain: primary, secondary, tertiary, and
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 informationInterrelationships. 1. Temperature Wind Fire Rainfall Soil Type Floods Sunlight Altitude Earthquake
Interrelationships Abiotic Factors A. A Partial List 1. Temperature Wind Fire Rainfall Soil Type Floods Sunlight Altitude Earthquake B. Aquatic Adaptations 1. Pumping salt out a. Salt water fish 2. Pumping
More informationPersistence of a generalized prey-predator model with prey reserved and herd behaviour of prey in
International Journal of Engineering Science Invention (IJESI) ISSN (Online): 239 6734, ISSN (Print): 239 6726 Volume 7 Issue 3 Ver. I March 208 PP. 40-48 Persistence of a Generalized Prey-Predator Model
More informationLOTKA-VOLTERRA SYSTEMS WITH DELAY
870 1994 133-140 133 LOTKA-VOLTERRA SYSTEMS WITH DELAY Zhengyi LU and Yasuhiro TAKEUCHI Department of Applied Mathematics, Faculty of Engineering, Shizuoka University, Hamamatsu 432, JAPAN ABSTRACT Sufftcient
More informationClimate change, ocean acidification and individual-based models: Why the little things matter
Climate change, ocean acidification and individual-based models: Why the little things matter Author Richards, Russell, Meynecke, Olaf, Chaloupka, M, Tomlinson, Rodger Benson Published 2012 Conference
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 information/Q~ =/V~ (b, + ~. a~inj), j=l
Bulletin of Mathematical Biology, Vol. 44, No. 4, pp. 579-585, 1982. Printed in Great Britain. 0092-8240[82/0405794)7503.00/0 Pergamon Press Ltd. 9 1982 Society for Mathematical Biology NOTE GRAPH THEORETCAL
More informationPERMANENCE IN LOGISTIC AND LOTKA-VOLTERRA SYSTEMS WITH DISPERSAL AND TIME DELAYS
Electronic Journal of Differential Equations, Vol. 2005(2005), No. 60, pp. 1 11. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) PERMANENCE
More informationStudents will be able to identify vocabulary related to Florida water and wetlands
Pre Visit 1 What s in your Water 6-8 Vocabulary Purpose/Objective Students will be able to identify vocabulary related to Florida water and wetlands Materials Word search copies Pencils Dictionaries Procedure
More informationNONSTANDARD NUMERICAL METHODS FOR A CLASS OF PREDATOR-PREY MODELS WITH PREDATOR INTERFERENCE
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 67 75. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationBrief Intro on Sea Plants
1 Brief Intro on Sea Plants This website will help you learn about the major varieties of underwater sea plants beneath the vast, deep oceans and seas. You will learn their names, colors, types, how they
More informationModeling Microbial Populations in the Chemostat
Modeling Microbial Populations in the Chemostat Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Modeling Microbial Populations in the Chemostat MBI, June 3, 204 / 34 Outline Why
More informationDynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting
Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting August 2012 Overview Last Model ẋ = αx(1 x) a(1 m)xy 1+c(1 m)x H(x) ẏ = dy + b(1 m)xy 1+c(1 m)x (1) where
More informationRegulation of Combined Harvesting of an Inshore- Offshore Fishery by Taxation
ISSN (Online): 39-7064 Index Copernicus Value (03): 6.4 Impact Factor (03): 4.438 Regulation of Combined Harvesting of an Inshore- Offshore Fishery by Taxation Tarun Pradhan Department of Mathematics,
More informationPhysiography Ocean Provinces p. 1 Dimensions p. 1 Physiographic Provinces p. 2 Continental Margin Province p. 2 Deep-Ocean Basin Province p.
Physiography Ocean Provinces p. 1 Dimensions p. 1 Physiographic Provinces p. 2 Continental Margin Province p. 2 Deep-Ocean Basin Province p. 2 Mid-Ocean Ridge Province p. 3 Benthic and Pelagic Provinces
More informationMath 2930 Worksheet Introduction to Differential Equations
Math 2930 Worksheet Introduction to Differential Equations Week 1 August 24, 2017 Question 1. Is the function y = 1 + t a solution to the differential equation How about the function y = 1 + 2t? How about
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 information115.3 Assignment #9 Solutions
115. Assignment #9 Solutions-1 115. Assignment #9 Solutions 8.1-12 Solve the differential equation d dx = 2(1 ), where 0 = 2 for x 0 = 0. d 1 = 2dx d 1 = 2dx ln 1 =2x + C Find C b inserting the Initial
More informationFatty acid and stable isotope (δ13c, δ15n) signatures of particulate organic matter in the Lower Amazon River: Seasonal contrasts and connectivity between floodplain lakes and the mainstem Mortillaro J.M.,
More information