vii southern-summer school on mathematical biology200.145.112.249/webcast/files/lecture_5_18.pdf ·...
TRANSCRIPT
![Page 1: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/1.jpg)
VII Southern-Summer School on Mathematical Biology
Roberto André Kraenkel, IFT
http://www.ift.unesp.br/users/kraenkel
Lecture V
São Paulo, January 2018
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 1 / 29
![Page 2: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/2.jpg)
Outline
1 Density and Diffusion
2 Reaction-Diffusion
3 Bibliography
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 2 / 29
![Page 3: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/3.jpg)
Outline
1 Density and Diffusion
2 Reaction-Diffusion
3 Bibliography
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 2 / 29
![Page 4: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/4.jpg)
Outline
1 Density and Diffusion
2 Reaction-Diffusion
3 Bibliography
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 2 / 29
![Page 5: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/5.jpg)
Space
Up to this point, all models that we have studied assume implicitly that all individuals are
in certain region of space.
This region has been supposed not to be very important .
We think of homogeneous regions.
Well-mixed populations.
HOWEVER...Individuals move, possibly generating the spatial redistribution of the population.
And space may be heterogeneous due to several factors :
Iclimate
Isoil
Ivegetation
Icomposition
Isalinity....
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 3 / 29
![Page 6: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/6.jpg)
Density
Let us consider a population in space.Let space be homogeneous. How do populations spread over space?.First point: we will not speak of number of individuals.Instead we will speak of density of individuals.The number of individuals per unit space.The usual notation is ⇢(~x , t) for density. It is a function of time and space.In some contexts, we use the term concentration.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 4 / 29
![Page 7: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/7.jpg)
Diffusion
Our main hypothesis is that individuals move randomly.
In some sense, they behave as molecules in a gas.If we look at such population from a space scale much larger than thetypical scale of the movement of the individuals we will see themacroscopic phenomenon called diffusion.Particles in a gas obey Fick’s law.We will assume the same for a population.So, what’s Fick’s law?
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 5 / 29
![Page 8: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/8.jpg)
Fick
The Fickian diffusion law states that:I
The flux
~J of "material"( animals, cells,..) is proportional to to the
gradient of the density of the material:
~J = �D
~r⇢ ⌘ �D(@⇢
@x,@⇢
@y)
Iwhere we took a two-dimensional space.
IBut to simplify the calculations let us consider the one-dimensional
case:
J ⇠ �@⇢
@x
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 6 / 29
![Page 9: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/9.jpg)
Mass/number of individuals conservation
Let us impose a conservation law:I
The rate of change in time of the quantity of individuals in a region of
space is equal to the flux through the borders.
that is, (in one dimension, (x0
� x
1
) being the size of the region):
@
@t
Zx1
x0
⇢(x, t)dx = J(x0, t) � J(x1, t)
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 7 / 29
![Page 10: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/10.jpg)
The diffusion equation
@@t
Rx
1
x
0
⇢(x , t)dx = J(x0
, t)� J(x1
, t)
We can write the previous equation in a differential form:I Take x1 = x0 + �x .I So that for �x ! 0:
FR
x1x0
⇢(x, t)dx ! ⇢(x0, t)�x
FJ(x1, t) ! J(x0, t) + �x
⇣@J(x,t)
@x
⌘
x=x0I Which implies::
@⇢@t
�x = ��x
✓@J(x, t)
@x
◆
I and using Fick’s law
@⇢@t
= �@J(x, t)@x
= D
@2⇢@x2
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 8 / 29
![Page 11: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/11.jpg)
The diffusion equation
@⇢@t
= D
@2⇢@x2
The above equation is known as the diffusion equation.In two dimensions we would have:
@⇢
@t= Dr2⇢
where r2⇢ ⌘ @2⇢@x2
+ @2⇢@y2
It is the same equation that describes heat diffusion if ⇢ is taken astemperature.Let us recall some facts about it.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 9 / 29
![Page 12: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/12.jpg)
Diffusion Equation
The diffusion equation is a partial differential equation, a PDE.
It is linear, and the coefficients are constants.
It can be solved analytically.
Mathematical comment
In order to speak of a solution of a differential equation, we need to specify
supplementary conditions.
In the case of the diffusion equation we should give an initial condition
⇢(x , 0) and the values of either ⇢(x , t) or
@⇢(x,t)@x at the borders or for
x ! ±1.
To solve it analytically, means that we can find a formula connecting ⇢(x , t)to ⇢(x , 0).
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 10 / 29
![Page 13: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/13.jpg)
Gauss
There is a distinctive solution: a Gaussian function.In one dimension we have, for t > 0 and x 2 [�1,+1]:
⇢(x, t) =Q
2(⇡Dt)1/2 e�x
2/(4Dt)
where Q os a constant.It is a Gaussian that "widens"with time.Corresponds to an initial condition concentrated in x = 0.Here is a plot.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 11 / 29
![Page 14: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/14.jpg)
Gauss: plots
Solution to the 1D diffusion equation
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 12 / 29
![Page 15: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/15.jpg)
Gauss: 2D plot
Solution to the 2D diffusion equation
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 13 / 29
![Page 16: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/16.jpg)
Diffusion, biology
Let us put some biology in this lecture!
Let us give a biological sense to all that.Suppose that at t = 0 a population of N individuals is released atx = 0.After a certain amount of time we want to know the the extensionoccupied by the population.Let’s be more specific: we want the extension of the region containing95% of the population.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 14 / 29
![Page 17: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/17.jpg)
Diffusion, biology
Knowing the density of a population allows us to calculate the totalpopulation in a given area. In the 1D case, we have:
Population between �L and L = N
L
=
Z +L
�L
⇢(x , t)dx .
If we use the Gaussian for ⇢(x , t),perform the integral, we obtain that95% of the population is a region of size 2
p2Dt.
Which grows in time proportional to t
1/2.Or, at a speed which goes like t
�1/2. Decreasing.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 15 / 29
![Page 18: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/18.jpg)
Diffusion + Growth
The previous case corresponds to a non-growing population.
Let us incorporate growth:
@⇢@t
= D
@2⇢@x2 + a⇢(x, t)
Still linear.
But, as we already learning, some saturation mechanism should become relevantfor large enough populations Say:
@⇢@t
= D
@2⇢@x2 + a⇢(x, t) � b⇢2(x, t)
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 16 / 29
![Page 19: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/19.jpg)
Fisher-Kolmogorov
@⇢@t
= D
@2⇢@x2 + a⇢(x, t) � b⇢2(x, t)
Figura : Robert. A. Fisher
Figura : Andrei N. Kolmogorov
The above equation is calledFisher-Kolmogorov equation.It is the simplest equation with diffusion,growth and self-regulation of a species.It is nonlinear.
It is a representative of the class of“reaction-diffusion” equations.
I This name comes from chemistry.The 2D version is obvious:
@⇢@t
= Dr2⇢ + a⇢ � b⇢2
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 17 / 29
![Page 20: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/20.jpg)
Fisher-Kolmogorov
@⇢@t
= D
@2⇢@x2 + a⇢(x, t) � b⇢2(x, t)
Let us again look at the problem of a population released at a point (x = 0).Suppose it obeys the Fisher-Kolmogorov equation (and not anymore the simplediffusion equation).No explicit formula.But look at the plot::
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 18 / 29
![Page 21: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/21.jpg)
Fisher-Kolmogorov
@⇢@t
= D
@2⇢@x2 + a⇢(x, t) � b⇢2(x, t)
We can see that there is a wave-front. And it moves with speed v = 2
paD.
Constant.
In the case of simple diffusion the speed decreased with time.This pattern can be made the basis of experimental verification.Our observations should concentrate on the front’s speed. .
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 19 / 29
![Page 22: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/22.jpg)
Skellam
The speed does not depend on b.Therefore, the constant wavefront speed is not related to densitydependence. The nonlinear term is there to avoid infinities.A equation
@⇢
@t= D
@2⇢
@x2 + a⇢(x, t)
is called the Skellam equation.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 20 / 29
![Page 23: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/23.jpg)
The classic example
Muskrat
The muskrat, a species native of North-america, was introduced inEurope.In 1905, five individuals were introduced in Prague.Today, there are millions in EuropeIn what follow, we see the expansion of the muskrat’s range aroundPrague over 17 years..
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 21 / 29
![Page 24: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/24.jpg)
Muskrat
1905
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 22 / 29
![Page 25: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/25.jpg)
Muskrat
1909
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 22 / 29
![Page 26: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/26.jpg)
Muskrat
1913
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 22 / 29
![Page 27: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/27.jpg)
Muskrat
1917
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 22 / 29
![Page 28: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/28.jpg)
Muskrat
1921
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 22 / 29
![Page 29: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/29.jpg)
Skellam !
From these observations we can estimate the speed of invasion as afunction of time.Here it is:
A straight line. Constant speed. Skellam dixit! REJOICE!.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 23 / 29
![Page 30: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/30.jpg)
Micro X macro
From the theory of the Brownian motion we can see D as the meansquare displacement per unit of time.We could try to track individuals and calculate it .Beware!, it is likely that you get a wrong value for D. Too large.Why?
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 24 / 29
![Page 31: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/31.jpg)
Home range effects
Many species have home ranges.This comes from several factors: the need to find food, the need tofind shelter .This slows down the diffusion process.In general, a mechanistic study of D is difficult. In most studies it istaken as a phenomenological constant.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 25 / 29
![Page 32: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/32.jpg)
Example: Hantavirus
In 2000, a new species of Hantavirus was discovered, being theetiological agent of f a respiratory syndrome. It is fatal in up to 60%of casesThe host is Oligoryzomys fulvescens. Take a look at him:
Where you find the rat, you find the HantavirusThe disease "follows"the spread of the rat.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 26 / 29
![Page 33: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/33.jpg)
Hantavirus II
The diffusion of the hosts is well modeled by the usual models,.But D is small.Oligoryzomys fulvescens has a limited home-range.The population spreads through juvenile migrants.A statistically rare event.But determinant for the spatial redistribution of the population.The diffusion coefficient appearing in the equations is a proxy of allthese processes.
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 27 / 29
![Page 34: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/34.jpg)
Bibliography
J.D. Murray: Mathematical Biology I e II (Springer, 2002)N.F. Britton: Essential Mathematical Biology ( Springer, 2003).R.S. Cantrell e C. Cosner: Spatial Ecology via Reaction-DiffusionEquations (Wiley, 2003).A. Okubo e S.A. Levin: Diffusion and Ecological Problems (Springer,2001).P. Turchin: Quantitative Analysis of Movement (Sinauer, 1998)
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 28 / 29
![Page 35: VII Southern-Summer School on Mathematical Biology200.145.112.249/webcast/files/lecture_5_18.pdf · Diffusion, biology Let us put some biology in this lecture! Let us give a biological](https://reader034.vdocument.in/reader034/viewer/2022042321/5f0a9d827e708231d42c80c4/html5/thumbnails/35.jpg)
Online Resources
http://ecologia.ib.usp.br/ssmb/
Thank you for your attention
Roberto A. Kraenkel (IFT-UNESP) VII SSSMB São Paulo, Jan 2018 29 / 29