Periodic dynamic systems for infected hosts and mosquitoes

Sistemas dinâmicos periódicos para hospedeiros e mosquitos infectados

Abstracts

A mathematical model for the purpose of analysing the dynamic of the populations of infected hosts anf infected mosquitoes when the populations of mosquitoes are periodic in time is here presented. By the computation of a parameter lambda (the spectral radius of a certain monodromy matrix) one can state that either the infection peters out naturally) (lambda <= 1) or if lambda > 1 the infection becomes endemic. The model generalizes previous models for malaria by considering the case of periodic coefficients; it is also a variation of that for gonorrhea. The main motivation for the consideration of this present model was the recent studies on mosquitoes at an experimental rice irrigation system, in the South-Eastern region of Brazil.

Malaria; Culicidae; Population dinamics; Ecology, vectors

Desenvolveu-se um modelo matemático para analisar a dinâmica das populações de indivíduos e mosquitos infectados quando as populações de mosquitos são periódicas no tempo. Pela determinação de um parâmetro lambda (o raio espectral de uma matriz de monodromia) pode-se estabelecer que a infecção termina naturalmente (lambda <= 1) ou se lambda > 1 que a infecção torna-se endêmica. O modelo generaliza, para o caso de coeficientes periódicos, modelos anteriores para malária; como também é uma variação de modelo para a gonorréia. A principal motivação para a consideração do modelo proposto foram os recentes estudos sobre mosquitos numa estação experimental de arroz irrigado, na região Sudeste do Brasil.

Malária; Culicidae; Dinâmica populacional; Ecologia de vetores

Periodic dynamic systems for infected hostsand mosquitoes

Sistemas dinâmicos periódicos para hospedeiros e mosquitos infectados

W. M. Oliva and E. M. Sallum

Departamento de Matemática Aplicada do Instituto Superior Técnico da Universidade Técnica de Lisboa. Portugal (W.M.O.), Instituto de Matemática e Estatística da Universidade de São Paulo. São Paulo, SP - Brasil (E.M.S.)

ABSTRACT

A mathematical model for the purpose of analysing the dynamic of the populations of infected hosts anf infected mosquitoes when the populations of mosquitoes are periodic in time is here presented. By the computation of a parameter l (the spectral radius of a certain monodromy matrix) one can state that either the infection peters out naturally) (l £ 1) or if l > 1 the infection becomes endemic. The model generalizes previous models for malaria by considering the case of periodic coefficients; it is also a variation of that for gonorrhea. The main motivation for the consideration of this present model was the recent studies on mosquitoes at an experimental rice irrigation system, in the South-Eastern region of Brazil.

Malaria, epidemiology. Culicidae. Population dinamics. Ecology, vectors.

RESUMO

Desenvolveu-se um modelo matemático para analisar a dinâmica das populações de indivíduos e mosquitos infectados quando as populações de mosquitos são periódicas no tempo. Pela determinação de um parâmetro l (o raio espectral de uma matriz de monodromia) pode-se estabelecer que a infecção termina naturalmente (l £ 1) ou se l > 1 que a infecção torna-se endêmica. O modelo generaliza, para o caso de coeficientes periódicos, modelos anteriores para malária; como também é uma variação de modelo para a gonorréia. A principal motivação para a consideração do modelo proposto foram os recentes estudos sobre mosquitos numa estação experimental de arroz irrigado, na região Sudeste do Brasil.

Malária, epidemiologia. Culicidae. Dinâmica populacional. Ecologia de vetores.

INTRODUCTION

A model was constructed for the purpose of analysing the dynamic of the populations of infected hosts and infected mosquitoes, when the populations of mosquitoes are periodic in time. A deterministic estimator l was established to predict whether the infection will peter out naturally (when l £ 1) or, if l > 1, the infection will become endemic (see Appendix). In the model, incubation and immunity are neglected.

The special case of one patch of hosts and two patches of mosquitoes are considered and the main results concerning the dynamics are given (see Theorems A and B). The estimator l is the spectral radius of the monodromy matrix C, corresponding to the periodic matrix A(t) of the linear system associated with the complete system (1). It is a non trivial matter to compute the matrix C as also to discover its spectral radius l.

The case of one patch of hosts and one (resp. two) patch(es) of mosquitoes is studied and it is assumed that the populations of mosquitoes behave according to periodic step functions. The monodromy matrix C is the product of exponentials of computable matrices; moreover, using suitable changes of coordinates, it was possible to obtain explicitly all the exponentials that appear in the expression of C, so that its spectral radius can be easily determined.

The general case (system (*)) that can be used for any number of patches of hosts and mosquitoes is presented in the Appendix.

The main motivation for the construction of the model here presented as well as the applications made was the recent studies by Forattini et al.⁴ on mosquitoes at an experimental rice irrigation system^{5, 6, 7, 8, 9}.

The constant and periodic coefficients that appear in the equations of the mathematical model should be determined experimentally. If they are available, reasonable predictions about possible mosquito-borne diseases that can appear at irrigation systems can be made.

ONE PATCH OF HOSTS AND TWO PATCHES OF MOSQUITOES

A deterministic model is here presented with a view to describing the dynamic of the population of hosts and mosquitoes infected by malaria when there is a homogeneous group of hosts with a constant population H and two groups of mosquitoes of different types i = 1, 2, with population V_i= V_i(t) periodic in the time t with periodic T > 0.

Let us considere the following system of ordinary differential equations:

where:

H : population of hosts;

V_i= V_i(t) = V_i(t + T) : population of mosquitoes at instant t, periodic of periodic T > 0;

S = S(t) : population of infected hosts at instant t;

I_i= I_i(t) : population of infected mosquitoes of type i at instant t;

d_i : death rate of mosquitoes of type i;

x : cure rate of sick hosts;

b'_i : bites by one mosquito of type i on hosts per unit of time;

b_i = b_i(t) : bites by mosquitoes of type i taken on one person, per unit of time, which is a periodic function with period T > 0;

f'_iI_i : population of infected mosquitoes of type i which are infective;

f_iS : population of infected hosts which are infective.

Since b_i(t)H = b'_iV_i(t), the system above can be written as:

where H, x, b'_i, d_i , f_i are positive constants and V_i = V_i(t) are continuous periodic functions of period T > 0.

That system (1) corresponds to a generalization of the Ross model (Lotka¹²) and of that of Dye-Hasibeder^{2, 3} for malaria. Moreover it is also a variation of the Aronson-Mellander¹ model that describes the dynamics of gonorrhea.

The main results A and B, stated below, for system (1), are special cases of general results proved in the Appendix.

Theorem A

Let (S(t), I₁(t), I₂(t)) be a non zero solution of (1) such that 0 £ S(t₀) £ H, 0 £ I_i(t₀) £ V_i(t₀), i = 1, 2, and some t₀³ 0. Then 0 < S(t) < H, 0 < I_i(t) < V_i(t), i = 1, 2, for all t > t₀.

We may write (1) in a matricial form

Let ø(t) be the matricial solution of = A(t)y such that ø(0) = I_d. The monodromy matrix C = ø(T) is positive (Aronson, Mellander¹, Lemma 2) and, by Perron's theorem (Gantmacher¹⁰), it has a simple positive eigenvalue l such that

Theorem B

There are two possibilities for the non zero solutions (S(t), I(t)) = (S(t), I₁(t), I₂(t)) of (1) such that for some t₀³ 0, 0 £ S(t₀) £ H and 0 £ I_i(t₀) £ V_i(t₀), i = 1, 2:

a) If l £ 1 then (S(t), I(t)) tends to the zero solution as t ® ¥ ;

b) If l > 1, then there exists a unique T -periodic solution (S*, I*) such that for any t > t₀ we have 0 < S*(t) < H and 0 < I*_i(t) < V_i(t), i = 1, 2. In this case (S(t) - S*(t), I(t) - I* (t)) tends to zero as t ® ¥ .

In other words, Theorem B says that the infection peters out naturally when l £ 1; or, if l > 1, the infection becomes endemic if the initial number of infected cases in at least one group is positive.

As usual, it is a non trivial matter to obtain the matrix C = ø(T) and an expression for l.

THE PERIODIC POPULATIONS OF MOSQUITOES AS STEP FUNCTIONS

One patch of hosts and one patch of mosquitoes

Let us consider the following system:

This describes the dynamic of malaria when we have one group of mosquitoes with a periodic population V = V(t) interacting with a homogeneous group of individuals of a fixed population H. It will be shown that for V = V(t) periodic with period T = 12 (months) with V(t) = V_i positive and constant, i < t < (i + 1), i = 1, ..., 12, the eigenvalues of the fundamental solution of the associated linear system can be computed directly showing their dependence on the data of system (2) that is, on the parameters x, d, H, b', f and on V =V(t).

Let us write (2) in the following form

For each i = 1, ..., 12 one considers ø_i(t) = , the fundamental solution of

If we make U = - dS - b'fI, (2)' becomes

one obtains

Let ø(t) be the matricial solution of such that ø(0) = I_d . Then one can write

So from (2)' we have

and so, from (6)

where the are the real eigenvalues of A(i) = A_i, j = 1, 2, i = 1, ..., 12. From (3),

so we have

and then

Finally we are able to obtain the eigenvalues of ø(T) that are those of the following matrix:

One Patch of Hosts and Two Patches of Mosquitoes

In this section we consider system (1) with d₁ = d₂ = d and V_i =V_i (t) as a periodic step function, which is positive and constant during the month i, i = 1, 2, ..., 12, and periodic with period equal to T = 12 (months). We write

The three invariants (coefficients of the characteristic polynomium) of A(t) are:

det A(t) = - xd ² +dj(t)

TrA(t) = - x - 2d (TrA(t) does not depend on t)

D₂A(t) = 2xd +d ²- j(t)

where

It is easy to check that (- d ) is a eigenvalue of A(t), say l₃ = - d . From the expressions of the invariants, the other two eigenvalues l₁ and l₂ satisfy.

These two last equations imply that

Note that the three eigenvalues of A(t) are real and distinct; moreover, we have l₂ <- d < l₁.

For each i = 1, 2, ..., 12, one considers ø_i(t) = e^tA(i) which is the fundamental solution of

If we make

system (9) becomes

or

where

From (10) we see that A(i) = M(i)M^{- 1} with M^{- 1} =

Let F(t) be the matricial solution of

such that I(0) = I_d. Then, as above, one can write

Let us write the matrix (i) of system (11) in the form

Then

and also

The eigenvalues of F(T) are e^-Td together with the eigenvalues of e^B(12)e^{B (11)}...e^B(1). So we only need to compute e^B(i) for i = 1, ..., 12. For that one has to solve the system

or, equivalently, to solve the second order equation

Since the characteristic roots of (13) are the eigenvalues of B(i), that is,

j = 1, 2, and i = 1, ..., 12, then, as we did in section 3.1, e^B(i) is given by the second member of equality (8).

APPENDIX

Let us consider the following system

which can be written as

y' = A(t)y + N(t, y) where

Assume a_i to be positive constants; c_i= c_i(t) to be positive continuous periodic functions of period T > 0; b_ji = b_ji(t) to be non negative continuous periodic functions of period T > 0 such that A(t) is an irreducible matrix (Gantmacher¹⁰), for all t. Moreover, assume there exists e > 0 such that b_ji(t) ³ e for all t provided that b_ji(t) are not identically zero.

For x = (x₁, ..., x_n) and y = (y₁, ..., y_n) we denote x £ y (x < y) if, for all i, x_i £ y_i ( x_i < y_i ). With the same arguments used in Aronson, Mellander¹ one can state:

I) If y(t) and z(t) are non-zero solutions of (*) such that for some t₀³ 0 we have 0 £ y(t₀) £ z(t₀) £ c(t₀) with y(t₀) ¹ z(t₀), then 0 < y(t) < z(t) < c(t) for all t > t₀.

For fixed t₀ ³ 0 one considers the map f(y₀) = y(t₀ + T, t₀, y₀) for 0 £ y₀£ c(t₀). When t₀ = 0 we take the sequence of positive numbers c(0) = c₀ > c₁ > c₂ > ... > c_n > ... > 0, where c_n = f0 (c_{n- 1}) and Q = c_n ³ 0. When Q = 0 we have y(t, t₀, y₀) = 0 for 0 £ y₀£ c(t₀). If Q 0 we have y(t) = y(t, 0, Q) positive, periodic with period T > 0 and, by I), Q > 0.

Let F(t) the matrix solution of = A(t)y, F(0) = I_d; since A(t) = (a_ij(t)) is irreducible with a_ij(t) ³ 0 for i ¹ j, then the matrix C = F(T) is positive and so, by Perron's theorem, it has a simple positive eigenvalue l = max{Rel_i; det(C - l_iI_d) = 0}.

As in theorem 1 of Aronson, Mellander¹ one has analogously:

II)If l < 1, there exists K > 0 such that | y(t)| £ |y(t₀)| for all t ³ t₀³ 0 and any solution y(t) = y (t, t₀,y₀) of (*) such that 0 £ y₀£ c(t₀).

Consider now the case l > 1. Let E(t₀) = {y Î IRⁿ: 0 £ y £ c(t₀)}, w > 0 eigenvector of C^tcorresponding to l, v(t₀)^t= w^tF(t₀)^{- 1} and E_q(t₀) = {y Î E(t₀) : v(t₀)^t y ³ q}. We claim that for q > 0 sufficiently small we have (E_q(t₀ )) Ì E_q(t₀) where 0 £ t₀£ T. In fact

Since

is continuous on the compact set

{(t₀, y) | t₀Î [0, T], y Î E(t₀)},

then there exist d > 0 and N > 0 such that ³ v(t₀)^ty + N | y | for t₀ Î [0, T] and | y | £ d. Using the Brower fixed point theorem (Hönig¹¹) one concludes that has a fixed point in E_q(t₀).

For any two solutions y(t) = y(t, t₀, y₀) and z(t) = y(t, t₀, z₀) of (*) such that 0 < y₀, z₀ < c(t₀), one has

D⁺u(t) £ - (u - 1) min(y_z, y_y) < 0 for t ³ t₀ ³ 0

(see Aronson and Mellander¹; Lemma 5, [1]).

So, one can conclude that for q > 0 sufficiently small, f : E_q(t₀) ® E_q(t₀) has only one fixed point Q, and then f : E(t₀) ®E(t₀) has only one fixed point Q > 0, besides the origin, that corresponds to a periodic orbit of period T > 0 for system (*).

For y(t) = y(t, t₀, y₀) with 0 £ y₀ £c(t₀) and z(t) = y(t, t₀, Q) we have

where and p is continuous and T-

periodic, T > 0. So, for each x₀ > 0 sufficiently close to the origin, there exists a constant M(x₀) > 0 such that

for all y₀, x₀ £ y₀ £c(t₀) and all t ³ t₀ ³ 0 because the c_k(t) are bounded functions (see Aronson, Mellander¹, Th. 2).

Given a compact set K Ì E(t₀), one considers a point x₀, 0 < x₀ <c(t₀), such that x₀ < y(t₀ + T, t₀, y₀) <c(t₀) for all y₀Î K. Since there are constants M(x₀) > 0 and N > 0 such that for all y₀ Î K we have

then one has the following result.

III) For l > 1, system (*) admits a unique non zero periodic solution y(t, t₀, Q), which has period T, and a constant a > 0 such that for each compact set K Ì E(t₀) there corresponds a constant M_K> 0 and we have

Moreover, as in theorem 3 of Aronson, Mellander¹ we state:

IV) For l = 1 there is a constant L > 0 such that for any solution y(t, t₀, y₀) = y(t) of (*) with 0 £ y₀ £c(t₀), t₀ ³ 0, we have

Correspondence to:

Elvia Mureb Sallum

Departamento de Matemática Aplicada do Instituto de Matemática e Estatística da Universidade de São Paulo.

R. do Matão, 1010. Cidade Universitária - 05508-900 São Paulo, SP - Brasil.

Fax: [phone]

E-mail: [e-mail]

The publication of this article was sponsored by FAPESP (Process 95/2290-6).

Received in 5.17.1995. Approved in 1.29.1996.

• 1. ARONSON, G. & MELLANDER, I. A deterministic model in biomathematics. Asymptotic behavior and threshold conditions. Math. Biosci., 49:207-22, 1980.
• 2. DYE, C & HASIBEDER, G. Population dynamics of mosquito-borne deseases: effects of flies which bite some peoples more frequently than others. Trans. R. Soc. Trop. Med. Hyg, 80-:69-77, 1986.
• 3. DYE, C. & HASIBEDER, G. Population dynamics of mosquito-borne diseases: persistence in a completely heterogeneous environment. Theoret. Pop. Biol., 33:31-53, 1988.
• 4. FORATTINI, O.P.; KAKITANI, I.; MASSAD, E.; GOMES, A. de C. Studies on mosquitoes (Diptera: Culicidae) and anthropic environment. 1. Parity of blood seeking Anopheles (Kerteszia) in South-Eastern Brazil. Rev. Saúde Pública, 27:1-18, 1993.
• 5. FORATTINI, O.P.; KAKITANI, I.; MASSAD, E.; MARUCCI, D. Studies on mosquitoes (Diptera: Culicidae) and anthropic environment. 2. Immature stages research at a rice irrigation system location in South-Eastern Brazil. Rev. Saúde Pública, 27:227-36, 1993.
• 6. FORATTINI, O.P.; KAKITANI, I.; MASSAD, E.; MARUCCI, D. Studies on mosquitoes (Diptera: Culicidae) and anthropic environment. 3. Survey of adult stages at the rice irrigation system and the emergence of Anopheles albitarsis in South-Eastern Brazil. Rev. Saúde Pública, 27:318-411, 1993.
• 7. FORATTINI, O.P.; KAKITANI, I.; MASSAD, E.; MARUCCI, D. Studies on mosquitoes (Diptera: Culicidae) and anthropic environment. 4. Survey of resting adults and synanthropic behavior in South-Eastern Brazil. Rev. Saúde Pública, 27:398-411, 1993.
• 8. FORATTINI, O.P.; KAKITANI, I.; MASSAD, E.; MARUCCI, D. Studies on mosquitoes (Diptera: Culicidae) and anthropic environment. 5. Breeding of Anopheles albitarsis in flooded rice fields in South-Eastern Brazil. Rev. Saúde Pública, 28:329-31, 1994.
• 9. FORATTINI, O.P.; KAKITANI, I.; MASSAD, E.; MARUCCI, D. Studies on mosquitoes (Diptera: Culicidae) and anthropic environment. 6. Breeding in empty conditions of rice fields in South-Eastern Brazil. Rev. Saúde Pública, 28:395-9, 1994.
• 10. GANTMACHER, F. R. The theory of matrices Vol. 2. New York, Chelsea, 1959.
• 11. HÖNIG, C. S. Aplicações da topologia à análise IMPA, Rio de Janeiro, 1976.
• 12. LOTKA, A. J. Contrib. to the analysis of malaria epidemiology. Am. J. Hyg., 3, 1923.

Publication Dates

• Publication in this collection
  03 Jan 2003
• Date of issue
  June 1996

History

• Accepted
  29 Jan 1996
• Received
  17 May 1995