Applied Mathematics
Vol.05 No.10(2014), Article ID:46590,10 pages
10.4236/am.2014.510147

Formulation of a Vector SIS Malaria Model in a Patchy Environment with Two Age Classes

Josephine Wairimu1,2*, Sallet Gauthier2, Wandera Ogana1

1School of mathematics, University of Nairobi, Nairobi, Kenya

2INRIA, Metz and University of Lorraine, Metz, France

Email: *jwndirangu@uonbi.ac.ke

Copyright © 2014 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/   Received 6 March 2014; revised 6 April 2014; accepted 13 April 2014

ABSTRACT

We formulate an SIS model describing transmission of highland malaria in Western Kenya. The host population is classified as children, age 1- 5 years and adults, above 5 years. The susceptibility and infectivity of an individual depend on age class and residence. The large scale system with 6n equations is reduced into a compact form of 3n equations by a change of variables. Then 3n equations are vectorialized using the matrix theory to get a one dimension, compact form of the system, equation in . Using Vidyasagar theorem  , the graph of the reduced system is shown to be strongly connected and the system is a monotone dynamical system. This means that circulation of malaria parasites among the species and among the patches is strongly connected, hence transmission is sustained. We show that for the n-dimensional age structured system the positive orthant is positively invariant for all positive values of the variables.

Keywords:

Highland Malaria, Differentiated Susceptibility and Infectivity, Monotone Dynamical Systems, Age Structure 1. Introduction

In Kenya, malaria is the leading cause of morbidity and mortality. It accounts for 30% of all outpatient attendances and 19% of all admissions to health facilities. The infection ranges from intense in the lowland to endemic in the highlands causing havoc to the public health system. About 20% of all deaths in children under five result from malaria (Ministry of Health-Kenya, 2006). People living around the lake, the coast, and the Western highlands epidemic-prone districts have 20% risk to be infected.

In the Western Kenya highlands, the risk of infection is 70%. The climate and topography influence the epidemic magnitudes. Drainage quality and rainfall determine vector breeding. The U-shaped valleys are broad and with slow moving rivers with poor drainage. This favor mosquitoes and high malaria infection and incidence.

The V-shaped valleys with narrow bottoms and fast flowing rivers with good drainage are less favorable to mosquitoes. The plateaus are flat, but have good drainage. Their ecosystem resembles the V-shaped valleys without dams. The terrain can modify the transmission of malaria. The V-shaped and the U-shaped valleys are separated from each other. In Kenya malaria is a “traveling disease”. 80% percent of the people treated for malaria in Kibera had travelled out of Nairobi. The neighboring estate is inhabited by people who originate from Lake Victoria region where the disease is widespread. By migrating to Nairobi, they are less exposed to malaria, thus loose the semi-immunity they used to have. This is why they contract it easily when traveling upcountry.

Children between 1 and 5 years of age are easily infected  . They are not bitten in same way as adults,   . Most deaths occur in infants and parasitemia levels of infected individuals decrease with age  . The large areas affected by malaria make the spraying of every house impossible,  . Hyman  , formulated a general differential susceptibility and differential infectivity model to prove that the disease free equilibrium is globally stable when and unstable otherwise. When there exists an endemic equilibrium, which solutions approach asymptotically. The model can only be used for direct transmission. Pongsumpun  modelled the influence of age structure in an SIS model for dengue hemorrhagic fever (DHF). They showed that age structure reduces the periods of oscillations on the susceptible and infected human population and on the infected mosquito population. The difference with us is the metapopulation setting, the differentiated patch, age susceptibility and infectivity. Gao and Ruan  examined how population dispersal affects malaria spread between patches. The residents migrate to other patches instead of making short visits as we shall assume. There is no age structure. They showed that travel can lead the disease to become endemic in both patches, even though the disease dies out in each isolated patch. Auger  modified Ross  model to patches without vectors migration. They model assumed that susceptibility and infectivity are similar in all patches.

Motivated by the work of Auger  and Pongsumpun  , we formulate an age structured model of malaria with susceptibility and infectivity depending on residence patch.

The region is subdivided into homogenous patches, and we distinguish children, aged 1-5, from adults, over 5 years of age. The two age groups are allowed to visit other patches other than their residences. Children move with the same rate as adults. Mosquitoes fly between patches less than 2 km apart, the approximate distance a mosquito can travel, Lutambi  and Kelly  .

2. The Model denotes the total number of susceptible children in patch , denotes the total number of susceptible adults in patch . denotes the total number of infected children in patch , denotes the total number of susceptible adults in patch. The vector population is likewise identified by and for the susceptible and infective vectors. The constant total human population comprises all hosts in all the patches: . So that denotes the total host population in patch. The total vector populations is given by the sum of all vectors in all the patches as and is the total vector population in patch. Migration of hosts is allowed in and out of the patch with an assumption that those who visit return to their resident patches (i.e. they make short visits sufficient for transmission to take place). denote the proportion of individuals travelling from from patch to patch. The flow of individuals into the different compartments is represented in Figure 1. We use to represent the U-shaped valleys and the V-shaped valleys. Due to the heterogenous nature of malaria in Western Kenya highlands, there are other patches characterized by other features not necessarily the ecosystem, hence the numbering.

The total host population on patch is defined as and its dynamics by, with

(1)

Figure 1. Transfer diagram for the different epidemiological groups in different patches.

where is the total number of residents of patch at time  .

We define the vectors, , and in. The matrix is

(2)

and the matrix analogously. Equation (1) becomes

The total population of mosquitoes satisfies,

If all the patches are sufficiently distant from each other, then mosquitoes do not migrate and

The dynamics of the total host and vector populations on patch is

(3)

where is the vector defined by

and the matrix is the diagonal block matrix

The migration equation becomes

(4)

We denote

is the recruitment of new born children in patch;

the susceptible host population and the infectious host population in patch;

the susceptible vector population and the infectious vector population in patch;

the man biting rate of vectors in patch;

the proportion of infectious bites on a host producing a patent infection;

the proportion of bites by susceptible vectors on an infected host producing a patent infection;

the head rate of children mortality in patch, the head rate of adult mortality in patch;

the mean duration of childhood;

the head rate of vector mortality;

the head rate of children recovery from infection, the head rate of adult recovery from infection and the mean duration of infection in patch.

The total host and vector populations are constant in all the patches. The parameter depends on the host age group in each patch. We also assume that the recovery rate does not depend on age  but on how quickly the host receives treatment and history of their malaria infection.

For patches we have the equations

(5)

The term represents the recruitment of newborns into the susceptible children’s class, while

in the first equation represents the differentiated infection of susceptible children by infec-

tions mosquitoes using frequency dependent transmission and a non constant host population in patch. The term defines the recovery of infected children in patch. The term is the per capita death rate of children in patch, is the rate children pass into adulthood, and the last term defines host migration. For the second and third equations the terms are defined similarly for Adults. In the sixth equation, for infectious vectors, the first term corresponds to the infection of susceptible mosquitoes by infected children and infected adults, while the last term caters for vector mortality.

We will assume in the sequel that for any index

The rationale for this assumption is that the rate of recovering, for an adult or a child, is considerably greater

than the mean sojourn time in the compartment of childhood. Actually and

The complete system is given by Equations (4) and (5).

3. Reduced System

It turns out that system given by Equation (5) can be rewritten in a triangular form so we need the following theorem to reduce such a system and thus study a smaller system.

Theorem 3.1 (Vidyasagar) Consider the following system:

(6)

If is globally asymptotically stable (GAS) in for the system, and if is GAS in, for the system, then is (locally) asymptotically stable for (6). Moreover if all the trajectories of (6) are forward bounded, then is a GAS for (6).

To apply the above Vidyasagar theorem we would need to prove the stability analysis of first equation, then we would only have to test the stability of the infection equation.

From Equation (5) if we add the first and third equation together we get

Which gives

We note that the matrix, is a Metzler matrix.

For a matrix the stability modulus, that we will denote, is the largest real part of the elements of the spectrum of

We have now

This proves that the stability modulus of satisfies

implying that this Metzler matrix is non singular, which in turn implies the opposite, that is, its inverse is nonnegative   . Therefore the equilibrium of this linear system is given by

and is globally asymptotically stable.

A similar result is obtained for the adult population with an equilibrium denoted by

and for the mosquito population with an equilibrium denoted by.

Reduction Process

We will now give different expressions for the equation of our system. Depending on the case at hand we will use the most convenient form to give the properties of this system and the corresponding proofs.

Using, and, we can rewrite system (5) in, with an immediate variable change

(7)

This system is clearly triangular if we consider the first variables. By application of Theorem (3.1), the stability analysis of (7) is now reduced to the stability analysis of the system

(8)

We set

Then Equation (8) can be written, in a vectorialized way, as

(9)

For another variable change, we set

Rewriting system (9) in terms of, and, we have

(10)

where we define the matrices

Finally, we will make a final “vectorization” of the system

We observe that where is the diagonal matrix

and are nonnegative matrices and we claim that is a stable Meztler matrix.

Indeed the stability modulus of is negative and

where the minimum is taken over the components of the 3 positive vectors.

We have

Since the matrices involved are Metzler matrices, this implies the following inequality for the corresponding stability modulus

The relation follows from and from which implies by Perron-Frobenius that.

Using the preceding matrices and the vector we rewrite Equation (10) in a compact form as

(11)

This system evolves on the unit cube of.

4. Basic Properties of the Model

For any index we shall use the classical notations: if and for some and we write if for all.

Proposition 4.1 (Positively Invariant Set)

The unit cube

is positively invariant for system (10).

Proof

To show the invariance of the unit cube, under the flow of the system (10), it suffices to show that each of the faces of the cube cannot be crossed.

On the patch we have

If

implying that cannot be crossed from positive to negative.

If we have

The relation is equivalent to. But is a demographic equilibrium satisfying

Taking this fact into account gives

or equivalently

which gives in turn, since

The equation for in patch is

If then

if a similar argument, as for the case gives

Finally the equation for in patch is given by

If, then for any,

and if a similar argument, as for the cases and gives

The proposition is proved.

Proposition 4.2

If the matrix is irreducible, then the system (11) is strongly monotone in the interior of the positively invariant set

Proof

We utilise the theory of monotone dynamical systems introduced by   , developed further in  and applied in  .

System (11) is monotone if its Jacobian is a Metzler matrix on the unit cube. The Jacobian of system (11) is given by

The Jacobian is clearly a Metzler matrix since, which implies that the system is cooperative in the unit cube.

Next, we show that the Jacobian is an irreducible matrix in the set. This will imply strong monotonicity of the system in the interior of the unit cube. In this set the diagonal terms of are positive, the same property is satisfied for, then the connectivity of the associated graph of is reduced to the connectivity of.

It is well known that a matrix is irreducible if its associated graph is strongly connected. Then only the off diagonal terms are concerned. Then it is sufficient to prove that the matrix

is irreducible. For the associated graph we distinguish three categories of vertices : the vertices corresponding to the (children vertices), (adults vertices), (mosquitoes vertices). Considering the matrix it is clear that a vertex is connected to the vertices of the same patch, i.e., and. In other words on the three vertices of patch are contained in a strongly connected component. The vertices from a patch is a subgraph and we can use the contraction of each in the general digraph to obtain a new digraph which represents movements between the patches. To prove the strong connectedness, it is sufficient to prove that this contracted digraph is strongly connected.

To the matrix, there is a corresponding digraph between the patches. This graph is equivalent to the associated digraph of. Since, and are Metzler matrices, we can consider the matrix and its corresponding multigraph on the graph of patches. This oriented multigraph represents actually the circulation of malaria parasites between the patches. Malaria parasites can be transported either by child, adults or mosquitoes.

Since we have assumed that the matrix is irreducible, our graph is strongly connected, consequently the Jacobian is irreducible. We note that we are in a context of multiple species (i.e., childs, adults and mosquitoes) in a metapopulation model as conceptualized in  .

Our assumption simply means that the circulation of parasites is strongly connected.

The establishment of the equilibrium solution and study of stability for the disease free and the endemic equilibrium will be tackled in a forth coming article.

5. Conclusions

In this study, we formulated an infinite model representing the spread of malaria in a heterogenous population classified. This population is classified as children (0 - 5 years) and adults (above 5 years). The infectivity and susceptibility of the population to highland malaria depend on: 1) the age class and; 2) the type of ecosystem the individual inhabits. The study captures the role played by the age of the individual and the ecosystem hetero- geneity in malaria epidemically spreading in the regions called patches. Hats are defined depending on their topography. Model properties are described to show that the solution set exists in the positive orthant which is positively invariant. The system is reduced to a single equation in, which is strongly connected and strongly monotone. The large scale system has been reduced using vectorialization into a compact form which is easier to study. The establishment of equilibrium solution and stability study will be done in a forthcoming article.

Acknowledgements

We wish to acknowledge of the Inria Metz, UMMISCO(IRD), the French Embassy in Nairobi and the Uni- versity of Nairobi, Kenya, for their financial, logistic and moral support during the writing of this article. We are very grateful to Dr. Githeko, KEMRI Kisumu for the great insight and literature he gave us during this study.

References

1. Vidyasagar, M. (1980) Decomposition Techniques for Large Scale Sytems with Nonadditive Interactions: Stability and Stabilization. IEEE Transactions on Automatic Control, 25, 773-779. http://dx.doi.org/10.1109/TAC.1980.1102422
2. Githeko, A.K., Ayisi, J.M., Odada, P.K., Atieli, F.K., Ndenga, B.A., Githure, J.I. and Yan, G. (2006) Topography and Malaria Transmission Heterogeneity in Western Kenya Highlands: Prospects for Focal Vector Control. Malaria Jour- nal, 5, 107. http://dx.doi.org/10.1186/1475-2875-5-107
3. Wolfgang, M.R., Willem, T., Richard, C. and Bart, K. (2002) Host-Specific Cues Cause Differential Attractiveness of Kenyan Men to the African Malaria Vector Anopheles Gambiae. Malaria Journal, 1, 17. http://dx.doi.org/10.1186/1475-2875-1-17
4. Smith, D.L., Guerra, C.A., Snow, R.W. and Simon, H.I. (2007) Standardizing Estimates of the Plasmodium Falciparum Parasite Rate. Malaria Journal, 6, 131. http://dx.doi.org/10.1186/1475-2875-6-131
5. Tumwiine, J., Mugisha, J.Y.T. and Luboobi, L.S. (2007) On Oscillatory Pattern of Malaria Dynamics in a Population with Temporary Immunity. Computational and Mathematical Methods in Medicine, 8, 191-203.
6. Wanjala, C.L., Waitumbi, J., Zhou, G. and Githeko, A.K. (2011) Identification of Malaria Transmission and Epidemic Hotspots in the Western Kenya Highlands: Its Application to Malaria Epidemic Prediction. Parasites and Vectors, 4, 81. http://dx.doi.org/10.1186/1756-3305-4-81
7. Hyman, J.M. and Li, J. (2006) Differential Susceptibility and Infectivity Epidemic Models. Mathematical Biosciences and Engineering, 3, 89-100. http://dx.doi.org/10.3934/mbe.2006.3.89
8. Pongsumpun, P. and Tang, I.M. (2003) Transmission of Dengue Haemorrhagic Fever in an Age Structured Population. Mathematical and Computer Modeling, 37, 949-961. http://dx.doi.org/10.1016/S0895-7177(03)00111-0
9. Gao, D. and Ruan, S. (2012) A Multipatch Malaria Model with Logistic Growth Populations. SIAM: SIAM Journal on Applied Mathematics, 72, 819-841. http://dx.doi.org/10.1137/110850761
10. Auger, P., Kouokam, E., Sallet, G., Tchuente, M. and Tsanou, B. (2008) The Ross-Macdonald Mode in a Patchy Envi- ronment. Mathematical Biosciences, 216, 123-131. http://dx.doi.org/10.1016/j.mbs.2008.08.010
11. Ross, R. (1911) The Prevention of Malaria. Springer-Verlag, Berlin.
12. Lutambi, A.M., Penny, M.A., Smith, N. and Chitnis N. (2013) Mathematical Modeling of Mosquito Dispersal in a He- terogenous Environment. Mathematical Biosciences, 241, 198-216.
13. Kelly, D.W. and Thompson, C.E. (2000) Epidemiology and Optimal Foraging: Modelling the Ideal Free Distribution of Insect Vectors. Parasitology, 120, 319-327.
14. Arino, J, Davis, J.R., Hartley, D., Jordan R., Miller, J. and van den Driessche, P. (2005) A Multispecies Epidemic Model with Spatial Dynamics. Mathematical Medical Biology, 22, 129-142.
15. Berman, A.P. and Robert, J. (1979) Nonnegative Matrices in the Mathematical Sciences, Volume 9 of Classics in Ap- plied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia. Revised Reprint of the 1979 Original.
16. Smith, H.L. (1995) Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems, Volume 41. Mathematical Surveys and Monographs.
17. Hirsch, M.W. (1982) Systems of Differential Equations That Are Competitive or Cooperative I: Limit Sets. SIAM: SIAM Journal on Applied Mathematics, 13, 167-179. http://dx.doi.org/10.1137/0513013
18. Hirsch, M.W. (1988) Systems of Differential Equations That Are Competitive or Cooperative III: Competing Species. SIAM: SIAM Journal on Applied Mathematics, 1, 51-71.
19. Hirsch, M.W. and Smith, H.L. (2005) Monotone Dynamical Systems. Handbook of Differential Equations: Ordinary Differential Equations. Vol. II, Elsevier B. V., Amsterdam, 239-357.
20. Arino, J. (2009) Diseases in Metapopulations. In: Ma, Z., Zhou, Y. and Wu, J., Eds., Modeling and Dynamics of Infec- tious Diseases in Series in Contemporary Applied Mathematics, World Scientific, Singapore City, 65-123. Also CDM Preprint Series Report 2008-04.

NOTES

*Corresponding author.