* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download OBU Template
Survey
Document related concepts
Ecological fitting wikipedia , lookup
Habitat conservation wikipedia , lookup
Unified neutral theory of biodiversity wikipedia , lookup
Biodiversity action plan wikipedia , lookup
Introduced species wikipedia , lookup
Latitudinal gradients in species diversity wikipedia , lookup
Maximum sustainable yield wikipedia , lookup
Island restoration wikipedia , lookup
Occupancyβabundance relationship wikipedia , lookup
Overexploitation wikipedia , lookup
Perovskia atriplicifolia wikipedia , lookup
Storage effect wikipedia , lookup
Reconciliation ecology wikipedia , lookup
Transcript
MATHEMATICAL MODELLING OF POPULATION DYNAMICS IN THEORETICAL ECOLOGY Department of HARRY GREEN Supervisors: DR. CRISTIANA SEBU, PROF. KHALED HAYATLEH, DR. TIM SHREEVE The research project explores the use of ordinary differential equations to model population change and the effects of predation in multiple species ecosystems. A system of π coupled ODEs can be used to model a system consisting of π interacting species. These systems are nonlinear and have the potential for chaos. The models are essentially a development on the Lotka-Volterra equations. We will examine accuracy of the models by comparison with real data using parameter fitting algorithms. To model single isolated populations in the absence of predation or competition we use a logistic growth curve given by π π = ππ 1 β πΎ (1) where πΎ is the carrying capacity, an upper limit on the species growth. In systems involving more species, logistic growth will be used to simulate the growth of the species at the bottom of the food chain (i.e., those that do not need to feed on others to grow. For predation terms, we use the βHolling Type II functional responseβ given by π π₯ = ππ₯ (2) π+π₯ Where a is the maximum predation rate. π(π₯) represents the biomass consumed by each member of the predator species y, which is converted into predator biomass with efficiency c. Assuming predators die at some rate d, for a two species food chain, we have the dynamical system π ππ₯ Mechanical Engineering and Mathematical Sciences The behaviour of a tritrophic food chain, when a superpredator π§ is introduced which feeds on the predator π¦, can be modelled by adding another dimension to (3). The solutions become three dimensional, and we encounter chaotic motion through a period doubling cascade. The solutions demonstrate a sensitivity to initial conditions, eliminating the possibility of any sort of long term prediction of the ecosystem. As the system is also dissipative i.e., all motion is bounded and all trajectories end in an attracting set, we frequently find strange attractors. π = ππ 1 β πΎ - π+π₯ π¦ (3) ππ₯ π=π π¦ β ππ¦ π+π₯ Which can exhibit a range of behaviour including spiraling and limit cycles. Bifurcation theory can reveal the existence of a Hopf Bifurcation and show for which parameter values the species can co-exist. The stability of the fixed point at the point πΎ, 0 determines the long term survival of the predator. π Modelling Human population growth with the π = ππ 1β (1) πΎ logistic curve (1). Data from UN estimates A strange chaotic attractor plotted from a tritrophic food chain model. The attractor demonstrates a sensitivity to initial conditions, and a small change in parameter values changes the shape unpredictably Modelling global human population growth with the logistic curve (1). Data from UN estimates and predictions. Human growth is affected little by other species, and is a good example of an isolated population. Modelling Human population growth with the Solutions of equation (3) for various initial conditions where the curve (1). Data from estimates redlogistic line is a boundary on the attracting set,UN calculated from the eigenvectors at (πΎ, 0) (K= 1 for this plot). The models given here are only a small and simple subset of the ODE based models used in theoretical ecology. Real world species display a much broader range of behaviour than simple Holling Type II predation. Using non-autonomous systems of ODEs we can account for seasonal variation in the systemβs dynamics. It is also frequently seen to combine the predator-prey models with epidemiological models designed to model the spreading of disease to investigate the effect of an epidemic disease on a multiple-species ecosystem.