Acessibilidade / Reportar erro

A Model for Aedes aegypti Infestation According to Meteorological Variables: Case of Caratinga (Minas Gerais - Brazil)

ABSTRACT

Aedes aegypti mosquitoes are the vector of diseases such as dengue, zika, chikungunya, and yellow fever among others. All the stages of development, eggs, larvae, pupa, and the adult of the species have its population modulated by meteorological variables, such as precipitation and temperature through affecting the productivity of breeding sites, metabolic processes, and others. Since adult females are responsible for transmitting the virus, the population of females becomes a direct indicator of the risk of infection. For this reason, some ongoing vector surveillance programs are based on adult female capture. In turn, all the stages of development have its population modulated by meteorological variables, such as precipitation and temperature, through productivity of breeding sites, metabolic processes and others. In this work, field data of capture of females was used to evaluate if a population dynamics model of Aedes aegypti under the effect of weather would be able to forecast field population. The nonlinear dynamic system model comprises: (1) four equations for the populations of the stages of development of the mosquito, designed for the ongoing surveillance program; (2) parametric dependencies of the rates of development on mean temperature and weekly accumulated precipitation. The dependencies on temperature and precipitation are modelled with aim of simplicity with the fewer number of parameters as possible. Temperature dependence is modelled based on values of the related literature under the assumption of existence of a optimum temperature for the rates, getting worse for extreme temperatures. The dependence on precipitation which is barely treated in experiments is modelled under the assumption of a monotonic dependence described by a power law with values estimated in orders of magnitude from data in the literature. By comparison with field data of an entomological indicator based on the number of Ae. aegypti females captured by a public health program in the city of Caratinga (Minas Gerais, Brazil), the model showed a significant correlation (R = 0.75). The result shows that the approach, if refined, can provide forecasting for of the population size.

Keywords:
Aedes aegypti; meteorological variables; dynamical systems

RESUMO

Os mosquitos Aedes aegypti são vetores de doenças como dengue, zika, chikungunya, febre amarela, entre outras. Como as fêmeas adultas são responsáveis pela transmissão do vírus, sua população torna-se um indicador direto do risco de infecção. Por esse motivo, alguns programas de vigilância de vetores baseiam-se na captura de fêmeas adultas. Por sua vez, todas as fases de desenvolvimento (ovos, larvas, pupas e o vetor adulto) têm sua população modulada por variáveis meteorológicas, como precipitação e temperatura, através da produtividade dos criadouros, processos metabólicos e outros. Neste trabalho, dados de captura de fêmeas em campo foram utilizados para avaliar se um modelo de dinâmica populacional de Aedes aegypti , sob o efeito de variáveis meteorológicas, seria capaz de prever a população de campo. O modelo de sistema dinâmico não linear compreende: quatro equações para as populações dos estágios de desenvolvimento do mosquito, projetadas para um programa de vigilância em andamento; (2) dependências paramétricas das taxas de desenvolvimento em relação a` temperatura média e precipitação acumulada semanalmente. As dependências de temperatura e precipitação são modeladas, por simplicidade, com o menor número de parâmetros possível. A dependência da temperatura é modelada com base em valores da literatura relacionada sob o pressuposto da existência de uma temperatura ótima para as taxas, piorando para temperaturas extremas. A dependência da precipitação, que mal é tratada em experimentos, é modelada sob o pressuposto de uma dependência monotônica descrita por uma lei de potência com valores estimados em ordens de magnitude a partir de dados da literatura. Por comparação com dados de campo de um indicador entomológico baseado no número de fêmeas Ae. aegypti capturadas por um programa de saúde pública na cidade de Caratinga (Minas Gerais, Brasil), o modelo apresentou correlação significativa (R = 0,75). O resultado mostra que a abordagem, se refinada, pode fornecer previsões para o tamanho da população.

Palavras-chave:
Aedes aegypti; variáveis meteorológicas; sistemas dinâmicos

1 INTRODUCTION

The Aedes (Stegomyia) aegypti (hereafter referred to as Ae. aegypti ), a species of the Culicidae family of mosquitoes, is the main vector of arboviruses like dengue fever, yellow fever, zika, and chikungunya 2727 World Health Organization. Dengue and severe dengue. http://www.who.int/mediacentre/factsheets/fs117/en/ (2017). Acessado em 20/03/2018.
http://www.who.int/mediacentre/factsheet...
. The species is found mainly in tropical and subtropical regions of the world and its distribution is closely related to climate patterns determined by variables such as temperature 11 E.B. Beserra, C.R. Fernandes, S.A.D.O. Silva, L.A.D. Silva & J.W.D. Santos. Efeitos da temperatura no ciclo de vida, exigências térmicas e estimativas do número de gerações anuais de Aedes aegypti (Diptera, Culicidae). Iheringia. Série Zoologia, 99 (2009), 142-148.),(1111 M.U.G. Kraemer, M.E. Sinka, K.A. Duda, A.Q. Mylne, F.M. Shearer, C.M. Barker, C.G. Moore, R.G. Carvalho, G.E. Coelho, W. Van Bortel et al. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. elife, 4 (2015), e08347.),(2121 G. Smith, D. Eliason, C. Moore & E. Ihenacho. Use of elevated temperatures to kill Aedes albopictus and Ae. aegypti. Journal of the American Mosquito Control Association, 4(4) (1988), 557-558.),(3030 H.M. Yang, M.L.G. Macoris, K.C. Galvani, M.T.M. Andrighetti & D.M.V. Wanderley. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and infection, 137(08) (2009), 1188-1202. , precipitation 1414 V.S.G. Neto & J.M.M. Rebêlo. Aspectos epidemiológicos do dengue no Município de São Luís, Maranhão, Brasil, 1997-2002. Cad. Saúde Pública, 20(5) (2004), 1424-1431. , humidity 55 E.A.P.D.A. Costa, E.M.D.M. Santos, J.C. Correia & C.M.R.D. Albuquerque. Impact of small variations in temperature and humidity on the reproductive activity and survival of Aedes aegypti (Diptera, Culicidae). Revista Brasileira de Entomologia, 54(3) (2010), 488-493. , and winds 99 A.D.S. Gomes, C.J. Sciavico & Á .E. Eiras. Periodicity of oviposition of females of Aedes aegypti (Linnaeus, 1762)(Diptera: Culicidae) in laboratory and field. Revista da Sociedade Brasileira de Medicina Tropical, 39(4) (2006), 327-332.),(2222 R. Souza-Santos & M.S. Carvalho. Análise da distribuição espacial de larvas de Aedes aegypti na Ilha do Governador, Rio de Janeiro, Brasil. Cadernos de saúde Pública, 16(1) (2000), 31-42. . Those factors act in the productivity of breeding sites, metabolic processes, behaviour, and other characteristics 1818 A.F. Ribeiro, G.R.A.M. Marques, J.C. Voltolini & M.L.F. Condino. Associação entre incidência de dengue e variáveis climáticas. Rev Saúde Pública, 40(4) (2006), 671-6. .

The development of Ae. aegypti takes place by complete metamorphosis 88 FUNASA. Programa Nacional de Controle da Dengue. Ministério da Saúde Vigilância Epidemiológica, (2002). . The stages of development occur in two distinct environments: (a) aquatic, composed by the egg, larva and pupa stages, and (b) aerial with the winged stage of the vector 1515 D.P. Neves. “Parasitologia humana”. Atheneu (2005). . To obtain nutrients for maturation and development of eggs, and energy reserves, adult gravid females perform blood feeding almost exclusively on humans 1919 T.W. Scott, P.H. Amerasinghe, A.C. Morrison, L.H. Lorenz, G.G. Clark, D. Strickman, P. Kittayapong & J.D. Edman. Longitudinal studies of Aedes aegypti (Diptera: Culicidae) in Thailand and Puerto Rico: blood feeding frequency. Journal of medical entomology, 37(1) (2000), 89-101. . Transmission of arboviruses to humans takes place when a female infected by the virus bites a susceptible human. Among the arboviruses, there are dengue fever, zika, chikungunya and yellow fever and others. According to the World Health Organization (WHO), just considering dengue fever, there occur 50 to 100 million new infections per year in more than 100 countries 2828 World Health Organization. Neglected tropical diseases. http://www.searo.who.int/entity/vector_borne_tropical_diseases/data/data_factsheet/en/ (2018). Acessado em 10/09/2018.
http://www.searo.who.int/entity/vector_b...
.

To avoid social impacts with epidemics of diseases transmitted by Ae. aegypti , the approach recommended by the WHO are the actions of entomological surveillance and control 2727 World Health Organization. Dengue and severe dengue. http://www.who.int/mediacentre/factsheets/fs117/en/ (2017). Acessado em 20/03/2018.
http://www.who.int/mediacentre/factsheet...
. The MIAEDES (Intelligent Aedes Monitoring), a system that make use of an adult mosquito trap called MosquiTRAP®. The MosquiTRAP® captures gravid females for monitoring the population of Ae. aegypti females.

In the surveilled municipalities, the traps are installed in the study area, spaced from each other by 200-250 meters, and monitored weekly. The surveillance data is summarised in terms of an entomological indicator MFAI (Mean Females Aedes Index), which is defined as the ratio of the total number of females caught during an epidemiological week to the number of positive traps inspected.

There is a long and well-known connection between meteorological and climatic factors and the rates of development stages of Ae. aegypti77 D.A. Focks, D.G. Haile, E. Daniels & G.A. Mount. Dynamic life table model for Aedes aegypti (Diptera: Culicidae): analysis of the literature and model development. Journal of medical entomology, 30(6) (1993), 1003-1017.),(1010 D.J. Gubler. The global emergence/resurgence of arboviral diseases as public health problems. Archives of medical research, 33(4) (2002), 330-342. . Such a connection produces a seasonal population pattern for Ae. aegypti , and consequently, for the number of cases of the infections. High precipitations and high temperatures increase the reproduction and survival of Ae. aegypti , markedly increasing the mosquito population size, and, consequently, the risk of infections 33 D.D. Chadee, F.L. Williams & U.D. Kitron. Impact of vector control on a dengue fever outbreak in Trinidad, West Indies, in 1998. Tropical Medicine & International Health, 10(8) (2005), 748-754.), (1010 D.J. Gubler. The global emergence/resurgence of arboviral diseases as public health problems. Archives of medical research, 33(4) (2002), 330-342. . The knowledge of the influence of the meteorological factors is important for the effectiveness of the environmental control. However, the adaptability of Ae. aegypti to urban environments and to human way of life, together with social conditions, allow the mosquito to survive throughout the year 33 D.D. Chadee, F.L. Williams & U.D. Kitron. Impact of vector control on a dengue fever outbreak in Trinidad, West Indies, in 1998. Tropical Medicine & International Health, 10(8) (2005), 748-754.), (1212 G. Kuno. Review of the factors modulating dengue transmission. Epidemiologic reviews, 17(2) (1995), 321-335. .

Mathematical and/or statistical models that consider the associations between meteorological variables and Ae. aegypti life cycle are developed with the purpose of describing the population dynamics of the vector, as well as the epidemics caused by that 11 E.B. Beserra, C.R. Fernandes, S.A.D.O. Silva, L.A.D. Silva & J.W.D. Santos. Efeitos da temperatura no ciclo de vida, exigências térmicas e estimativas do número de gerações anuais de Aedes aegypti (Diptera, Culicidae). Iheringia. Série Zoologia, 99 (2009), 142-148.), (22 M. Burattini, M. Chen, A. Chow, F. Coutinho, K. Goh, L. Lopez, S. Ma & E. Massad. Modelling the control strategies against dengue in Singapore. Epidemiology & Infection, 136(3) (2008), 309-319.), (77 D.A. Focks, D.G. Haile, E. Daniels & G.A. Mount. Dynamic life table model for Aedes aegypti (Diptera: Culicidae): analysis of the literature and model development. Journal of medical entomology, 30(6) (1993), 1003-1017.), (1717 M. Otero, H.G. Solari & N. Schweigmann. A stochastic population dynamics model for Aedes aegypti: formulation and application to a city with temperate climate. Bulletin of mathematical biology, 68(8) (2006), 1945-1974.), (2424 A. Tran, G. L’Ambert, G. Lacour, R. Benôıt, M. Demarchi, M. Cros, P. Cailly, M. Aubry-Kientz, T. Balenghien & P. Ezanno. A rainfall-and temperature-driven abundance model for Aedes albopictus populations. International journal of environmental research and public health, 10(5) (2013), 1698-1719.), (3030 H.M. Yang, M.L.G. Macoris, K.C. Galvani, M.T.M. Andrighetti & D.M.V. Wanderley. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and infection, 137(08) (2009), 1188-1202. .

In this work, a mathematical model based on differential equations is used to describe the Ae. aegypti life cycle as a function of meteorological data (temperature and precipitation). The main goal is evaluate if a simple model under the effect of weather would be able to forecast field population. Hence it was presented a nonlinear dynamic system model that comprises: (1) four equations for the populations of the stages of development of the mosquito, designed for the ongoing surveillance program; (2) parametric dependencies of the rates of development on mean temperature and weekly accumulated precipitation.

The dependencies of the entomological rates on temperature and precipitation are modelled with aim of simplicity with the fewer number of parameters as possible. Temperature dependence is modelled based on values of the related literature with polynomials of degree 2, under the assumption of existence of a optimum temperature for the rates, and getting worse for extreme temperatures 11 E.B. Beserra, C.R. Fernandes, S.A.D.O. Silva, L.A.D. Silva & J.W.D. Santos. Efeitos da temperatura no ciclo de vida, exigências térmicas e estimativas do número de gerações anuais de Aedes aegypti (Diptera, Culicidae). Iheringia. Série Zoologia, 99 (2009), 142-148.), (2121 G. Smith, D. Eliason, C. Moore & E. Ihenacho. Use of elevated temperatures to kill Aedes albopictus and Ae. aegypti. Journal of the American Mosquito Control Association, 4(4) (1988), 557-558.), (3030 H.M. Yang, M.L.G. Macoris, K.C. Galvani, M.T.M. Andrighetti & D.M.V. Wanderley. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and infection, 137(08) (2009), 1188-1202. . The dependence on precipitation, which is barely found in experiments, is modelled under the assumption of a monotonic dependence through a power law 2525 J. Waldock, N.L. Chandra, J. Lelieveld, Y. Proestos, E. Michael, G. Christophides & P.E. Parham. The role of environmental variables on Aedes albopictus biology and chikungunya epidemiology. Pathogens and global health, 107(5) (2013), 224-241. . The power law constants are estimated by orders of magnitude from the literature that refers to entomological rates, regardless if it refers to precipitation itself 11 E.B. Beserra, C.R. Fernandes, S.A.D.O. Silva, L.A.D. Silva & J.W.D. Santos. Efeitos da temperatura no ciclo de vida, exigências térmicas e estimativas do número de gerações anuais de Aedes aegypti (Diptera, Culicidae). Iheringia. Série Zoologia, 99 (2009), 142-148.), (77 D.A. Focks, D.G. Haile, E. Daniels & G.A. Mount. Dynamic life table model for Aedes aegypti (Diptera: Culicidae): analysis of the literature and model development. Journal of medical entomology, 30(6) (1993), 1003-1017.), (1717 M. Otero, H.G. Solari & N. Schweigmann. A stochastic population dynamics model for Aedes aegypti: formulation and application to a city with temperate climate. Bulletin of mathematical biology, 68(8) (2006), 1945-1974.), (2424 A. Tran, G. L’Ambert, G. Lacour, R. Benôıt, M. Demarchi, M. Cros, P. Cailly, M. Aubry-Kientz, T. Balenghien & P. Ezanno. A rainfall-and temperature-driven abundance model for Aedes albopictus populations. International journal of environmental research and public health, 10(5) (2013), 1698-1719.), (2929 H.M. Yang. The transovarial transmission in the dynamics of dengue infection: Epidemiological implications and thresholds. Mathematical biosciences, 286 (2017), 1-15.), (3030 H.M. Yang, M.L.G. Macoris, K.C. Galvani, M.T.M. Andrighetti & D.M.V. Wanderley. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and infection, 137(08) (2009), 1188-1202. . Although the determination of those parameters from experiments is of limited realism, it is intend to verify the forecasting power of the model as a evaluation of that approach for eventual future refinement. For that purpose, the results of the computational model implementation will be compared to experimental data from Aedes females (MFAI) from the municipality of Caratinga, Minas Gerais, Brazil.

In Section 2, the entomological data of captures and meteorological data are characterised. In Section 3, the mathematical model for the life cycle of Aedes aegypti is detailed, considering the dependence with the meteorological variables. A description of the simulation and validation of the results is shown in Section 4. The final considerations are made in Section 5.

2 DATA BASE

The data used in this work refers to the epidemiological weeks 23-52 of 2009 and epidemiological weeks 1-51 of 2010 in the city of Caratinga (Minas Gerais, Brazil). The climate of the city is classified as tropical semi-humid. Average temperatures are around 22C, with a mean average of 16 . 7C and average maximum of 27 . 7C (averages considering the period from 1961 to 2014). The average annual rainfall is 1 , 122 . 6mm 1313 A.P.M. Moreira, V.A.L. C. & S.J. C. Tendências climáticas e anomalias de precipitação em CaratingaMG. In “Os Desafios da Geografia Física na Fronteira do Conhecimento” (2018), pp. 2000-2009. .

The weekly meteorological data, accumulated precipitation p (mm) and average temperature T (C) for the study period were obtained through the MDBTR (Meteorological Data Bank for Teaching and Research) on the website of the National Meteorological Institute INMET (www.inmet.gov.br).

Experimental data from the Ae. aegypti females population were obtained by using MosquiTRAP® in a public heath surveillance from the study area. During the period of study of 96 epidemiological weeks, the MFAI was generated in the MI-Aedes program developed in the municipality by the biotechnology company Ecovec S. A. (Belo Horizonte, MG, Brazil).

3 MATHEMATICAL MODEL

To evaluate the viability of a predictive population model with dependencies on mean temperature and precipitation, it is presented a nonlinear dynamic system model that comprises: (1) four equations for the populations of the stages of development of the mosquito, designed for the ongoing surveillance program; (2) parametric dependencies of the rates of development on mean temperature and weekly accumulated precipitation. The population model (3.1) written as a system of differential equations is proposed to describe the dynamics of the population size of the four stages of development Ae. aegypti , considering the weekly variation of precipitation and average temperature. The population of eggs, aquatic forms (larvae and pupae), pre-bloodmeal and post-bloodmeal females are represented by E ( t ), A ( t ), F 1( t ) e F 2( t ), respectively (see Figure 1 ).

Figure 1:
Diagram of the model populations related to the life cycle of Aedes aegypti with rates of development as parametric dependent on temperature T and precipitation p .

d E d t = ϕ p , T 1 - E t C p , T F 2 t - α 1 p , T E t - μ E p , T E t , d A d t = α 1 p , T E t - α 2 p , T A t - μ A p , T A t , d F 1 d t = α 2 p , T A t - α 3 p , T F 1 t - μ F 1 p , T F 1 t , d F 2 d t = α 3 p , T F 1 t - - μ F 2 p , T F 2 t , C , ϕ , α 1 , α 2 , α 3 , μ E , μ A , μ F 1 , μ F 2 0 , p , T , t + . (3.1)

The coefficients of the model are parametric functions of the weekly precipitation index p and the mean temperature T . The oviposition rate carried out by post-bloodmeal females is symbolised by ϕ. The rate α 1 corresponds to the development of eggs in an aquatic population. In turn, the rate of individuals of the aquatic population that develop to blood pre-bloodmeal females is represented by α 2. Rate α 3 represents the maturation rate of pre-bloodmeal females to blood post-bloodmeal females. The natural mortality rates of the populations involved are given by µ E , µ A , μF1 e μF2 , respectively, linear term, ϕ1-ECF2 includes itself the term -ϕEF2C , which stands to mitigate rate ϕ if the population of stage E becomes large enough compared to the value of C66 L. Esteva & H.M. Yang. Mathematical model to assess the control of Aedes aegypti mosquitoes by the sterile insect technique. Mathematical biosciences, 198(2) (2005), 132-147. .

The relationship between the life cycle of Ae. aegypti and the temperature is strongly consolidated in the literature 11 E.B. Beserra, C.R. Fernandes, S.A.D.O. Silva, L.A.D. Silva & J.W.D. Santos. Efeitos da temperatura no ciclo de vida, exigências térmicas e estimativas do número de gerações anuais de Aedes aegypti (Diptera, Culicidae). Iheringia. Série Zoologia, 99 (2009), 142-148.), (3030 H.M. Yang, M.L.G. Macoris, K.C. Galvani, M.T.M. Andrighetti & D.M.V. Wanderley. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and infection, 137(08) (2009), 1188-1202. . In turn, the relationship between the population of Ae. aegypti and precipitation has not been sufficiently studied in laboratories or field experiments 2525 J. Waldock, N.L. Chandra, J. Lelieveld, Y. Proestos, E. Michael, G. Christophides & P.E. Parham. The role of environmental variables on Aedes albopictus biology and chikungunya epidemiology. Pathogens and global health, 107(5) (2013), 224-241. . The association between those dependencies is therefore ignored so far. Let Π=ϕ,α1,α2,α3,μE,μA,μF1,μF2 be a general representation of the model (3.1) rates, and let the parametric dependencies of the rates on precipitation p and temperature T be, respectively, Π=Ψp and Γ=ΓT . Then, the design for the dependence of the rate Π on p and T could be constructed by supposing it a well behaved enough function (analytic or class C n ), such that ΠΨ,Γ=a0+a1Ψ+b1Γ+OΨ2,Γ2 . In this case, the association of the dependencies of the model rates on p and T would be considered in a general form. By simplicity, as a working hypothesis, and because the model (daily) rates Π are typically smaller than one, we retain the linear part to the purpose of the association:

Π p , T = Ψ p + Γ T . (3.2)

Hence, by such simplifying consideration, the dependencies on p and T turn to be additive, as if they were rates itself in parallel association. Furthermore, such approach helps in minimisation of the number of free parameters.

To provide a insight on the parametric dependency of rates on precipitation Ψ ( p ), let us focus in the Figure 2a and notice that the peaks pattern resembles the precipitation peaks pattern. This fact is also manifested in Figure 2b , where there are significant correlations ( p < 0 . 05) between precipitation and MFAI for lags ranging from 0 to 12. Because of that similarity, the dependency of the rates on precipitation is supposed to be a monotonously increasing function. Other premises adopted here are simplicity and the fewest possible parameters. For those reasons, the parametric function Ψ=Ψp was written in terms of a power law, as

Ψ p = Ψ 0 + p - p 0 p 1 - p 0 r Ψ 1 - Ψ 0 ; with p p 0 . (3.3)

Figure 2:
(a) Time series of precipitation, temperature and MFAI data from Caratinga-MG, Brazil, and (b) cross correlation between precipitation and MFAI time-series. The crossing curve stands for the significance level. Lags between 1 to 12 are associated to significant correlations ( p < 0 . 05)

The values Ψ 0 and Ψ 0 are associated to p 0 = 0 mm and p 1 = 34 . 62 mm, that correspond in weekly base to the threshold of rainfall volume between subtropical and tropical rainforest (1800mm/year) 1616 S. Noguchi, A.R. Nik & M. Tani. Rainfall characteristics of tropical rainforest at Pasoh forest reserve, Negeri Sembilan, Peninsular Malaysia. In “Pasoh”. Springer (2003), pp. 51-58. . This power law parametric equation can provide increasing, decreasing, constant functions which can be also concave up and down. In the present work, we consider r > 0, focusing on increasing dependencies.

More exact relations for the dependencies between rates and precipitation should be obtained by optimisation or more complex experiments. The values adopted in this work ( Table 1 ) belong to the spectrum of values from other studies, regardless if they are related to precipitation. Due to this lack of exact values the constants of Equation (3.3) are set here in terms of orders of magnitude.

Table 1:
Free parameters of the dependencies of the model rates on temperature (3.4) and precipitation (3.3) with references which based their choice.

Although there are many studies for the dependence of the rates on temperature 11 E.B. Beserra, C.R. Fernandes, S.A.D.O. Silva, L.A.D. Silva & J.W.D. Santos. Efeitos da temperatura no ciclo de vida, exigências térmicas e estimativas do número de gerações anuais de Aedes aegypti (Diptera, Culicidae). Iheringia. Série Zoologia, 99 (2009), 142-148.), (2121 G. Smith, D. Eliason, C. Moore & E. Ihenacho. Use of elevated temperatures to kill Aedes albopictus and Ae. aegypti. Journal of the American Mosquito Control Association, 4(4) (1988), 557-558.), (3030 H.M. Yang, M.L.G. Macoris, K.C. Galvani, M.T.M. Andrighetti & D.M.V. Wanderley. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and infection, 137(08) (2009), 1188-1202. , in a first approach and as a working hypothesis, we adopt simple functions according to the following assumptions: (1) simple relation with a small number of free parameters; (2) the existence of an optimal temperature for the metabolic processes. As a consequence, there arises a maximum value for the rates of development ( Figures 3a , 3b and 3f ), or a minimum, for the mortality rates ( Figures 3c , 3d and 3e ). Once the temperature gets far from the optimum value, it is then expected that there decreases the fitness in the life cycle. Hence, for the case of rates of development ( α 1, α 2 and ϕ), as far as the temperature is from the optimum value, their values decrease until vanishing for very high or very low temperatures. In the case of mortalities ( µ E , µ A and µ F ), their values increase as the temperature becomes far from the optimum value. To fulfill these conditions, the rates were modelled by simple polynomials of degree 2, concave up for mortalities, and down for development rates with vanishing values in the outer intervals from the roots to avoid negative values.

Γ T = a T - T o p 2 + b . (3.4)

Figure 3:
Parameter setting α1, α2, µE , µA and µF , considering experimental data available in 11 E.B. Beserra, C.R. Fernandes, S.A.D.O. Silva, L.A.D. Silva & J.W.D. Santos. Efeitos da temperatura no ciclo de vida, exigências térmicas e estimativas do número de gerações anuais de Aedes aegypti (Diptera, Culicidae). Iheringia. Série Zoologia, 99 (2009), 142-148.), (2121 G. Smith, D. Eliason, C. Moore & E. Ihenacho. Use of elevated temperatures to kill Aedes albopictus and Ae. aegypti. Journal of the American Mosquito Control Association, 4(4) (1988), 557-558.), (3030 H.M. Yang, M.L.G. Macoris, K.C. Galvani, M.T.M. Andrighetti & D.M.V. Wanderley. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and infection, 137(08) (2009), 1188-1202. .

The values of a , b , and T op were obtained by fitting the degree 2 polynomials to experimental data present in some works 11 E.B. Beserra, C.R. Fernandes, S.A.D.O. Silva, L.A.D. Silva & J.W.D. Santos. Efeitos da temperatura no ciclo de vida, exigências térmicas e estimativas do número de gerações anuais de Aedes aegypti (Diptera, Culicidae). Iheringia. Série Zoologia, 99 (2009), 142-148.), (2121 G. Smith, D. Eliason, C. Moore & E. Ihenacho. Use of elevated temperatures to kill Aedes albopictus and Ae. aegypti. Journal of the American Mosquito Control Association, 4(4) (1988), 557-558.), (3030 H.M. Yang, M.L.G. Macoris, K.C. Galvani, M.T.M. Andrighetti & D.M.V. Wanderley. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and infection, 137(08) (2009), 1188-1202. . The idea is more to find the best polynomials of degree 2 according to data, than to provide a high goodness of fit. Figure 3 depicts the modelled parameters Γ( T ) together with the set of experimental data in related literature. The parameters of rate α 3 were determined by imposing three conditions: the inverse of the mean lifetime of F 1 of 5 days was set as the mean value of the function (3.4), and the two roots were chosen to be the temperatures 15 C and 40 C. Because outside of such interval the females do not bite, the values of T < 15 C and T > 40 C are set to be zero 44 S.R. Christophers. “Aedes aegypti: the yellow fever mosquito”. CUP Archive (1960).), (2121 G. Smith, D. Eliason, C. Moore & E. Ihenacho. Use of elevated temperatures to kill Aedes albopictus and Ae. aegypti. Journal of the American Mosquito Control Association, 4(4) (1988), 557-558. . In the cases of rates of development, outside the roots of the polynomials, the values in the model are fixed to zero, although it never happens in the present range of data.

The system (3.1) is non-autonomous because it has coefficients dependent on rainfall and temperature which, in turn, depend o time 2020 L.D.S.B. Silva, R.T.N. Cardoso, J.L.A. Fernandes, C.A. Silva & A.E. Eiras. Modelo Entomológico Determinístico sob Efeito da Pluviosidade para o Aedes aegypti e o Aedes albopictus. TEMA (São Carlos), 19(2) (2018), 289-303. . However the temperature and precipitation data is updated only weekly and the time evolution is hourly. So the evolution of the system along the yer can be viewed as a succession of autonomous systems. Hence, in each week, the system parameters become constant in time. Under such assumption, by applying the condition of equilibrium dXxt=0 with X=E,A,F1,F2 in model (3.1), we identify a trivial equilibrium point P 0 and a non-trivial equilibrium point P 1, described by

P 0 = E * , A * , F 1 * , F 2 * = 0 , 0 , 0 , 0 (3.5)

P 1 = E * * , A * * , F 1 * * , F 2 * * = 1 - 1 Q 0 C , α 1 α 2 + μ A E * * , α 2 α 3 + μ F 1 A * * , α 3 μ F 2 F 1 * * (3.6)

in which Q 0 is the vector reproduction number, given by

Q 0 = α 1 α 1 + μ E α 2 α 2 + μ A α 3 α 3 + μ F 1 ϕ μ F 2 . (3.7)

The equilibrium points (see Appendix A Appendix A STABILITY OF THE EQUILIBRIUM POINTS Let the system of equations (3.1) rewritten as d X d t = f X , where X = E , A , F 1 , F 2 . (A.1) The condition fX=0 revealed two critical points, one trivial, P 0(Equation 3.5), and one nontrivial, P 1(Equation 3.6). In the following it is performed the analysis of stability of the system. After the expansion near the equilibrium points P k ; k = 0;1 and computation of the Jacobian matrices BPk=∂f∂XPk , k = 0 , 1, their characteristic polynomials are given, respectively, by Φ P k λ = a k 4 λ 4 + a k 3 λ 3 + a k 2 λ 2 + a k 1 λ + a k 0 , (A.2) with coefficients for P 0, a 04 = 1 , a 03 = α 1 + μ E + α 2 + μ A + α 3 + μ F 1 + μ F 2 , a 02 = α 1 + μ E μ F 2 + α 2 + α 3 + μ A + μ F 1 α 1 + μ E + μ F 2 + α 2 + μ A α 3 + μ F 1 , a 01 = α 1 + μ E α 2 + μ A α 3 + μ F 1 + μ F 2 + α 1 + α 2 + μ E + μ A α 3 + μ F 1 μ F 2 , a 00 = α 1 + μ E α 2 + μ A α 3 + μ F 1 μ F 2 1 - Q 0 ; and for P 1, a 14 = 1 , a 13 = ( α 1 + μ E ) Q 0 + α 2 + α 3 + μ A + μ F 1 + μ F 2 , a 12 = ( α 1 + μ E ) ( α 2 + α 3 + μ A + μ F 1 + μ F 2 ) Q 0 + ( α 2 + μ A ) ( α 3 + μ F 1 + μ F 2 ) + ( α 3 + μ F 1 ) μ F 2 , a 11 = ( α 1 + μ E ) ( α 2 + μ A ) ( α 3 + μ F 1 + μ F 2 ) Q 0 + ( ( α 1 + μ E ) Q 0 + α 2 + μ A ) ( α 3 + μ F 1 ) μ F 2 , a 10 = ϕ α 1 α 2 α 3 ( 1 − 1 Q 0 ) . The stability of each equilibrium points P k can be assessed from the sign of the real part of their respective eigenvalues. The polynomial ΦPk (A.2) with have roots with real part strictly less than zero if and only if all elements in the first column of the Routh Table (A.3) are nonzero and have the same sign 26 . R k = 1 a k 2 a k 0 a k 3 a k 1 0 a k 2 a k 3 - a k 1 a k 3 a k 0 0 a k 1 - a k 3 2 a k 0 a k 2 a k 3 - a k 1 0 0 a k 0 0 0 (A.3) The elements R11k=1>0 , k = 0 , 1 are positive. Since all model rates remain positive by definition, R21k>0 , k = 0 , 1 is positive. As consequence, R31k , k = 0 , 1, are positive if ak2ak3-ak1>0 . For k = 0, it reads a 02 a 03 − a 01 = 2 ( α 1 + μ E ) ( α 2 + α 3 + μ A + μ F 1 ) μ F 2 + 2 ( α 2 + μ A ) ( α 3 + μ F 1 ) ( α 1 + μ E + μ F 2 ) + ( α 1 + μ E ) 2 ( α 2 + α 3 + μ A + μ F 1 + μ F 2 ) + ( α 2 + μ A ) 2 ( α 1 + α 3 + μ E + μ F 1 + μ F 2 ) + ( α 3 + μ F 1 ) 2 ( α 1 + α 2 + μ E + μ A + μ F 2 ) + ( μ F 2 ) 2 ( α 1 + α 2 + α 3 + μ E + μ A + μ F 1 ) > 0 . In the case k = 1, it is expressed as a 12 a 13 − a 11 = ( α 1 + μ E ) 2 ( α 2 + α 3 + μ A + μ F 1 + μ F 2 ) Q 0 2 + 2 ( α 2 + μ A ) ( α 3 + μ F 1 ) μ F 2 + ( α 1 + μ E ) ( α 2 + α 3 + μ A + μ F 1 + μ F 2 ) ( α 2 + α 3 + μ A + μ F 1 ) Q 0 + ( α 2 + μ A ) 2 ( α 3 + μ F 1 + μ F 2 ) + μ F 2 2 ( ( α 1 + μ E ) Q 0 + α 2 + α 3 + μ A + μ F 1 ) + ( α 3 + μ F 1 ) 2 ( α 2 + μ A + μ F 2 ) + ( α 1 + μ E ) ( α 2 + α 3 + μ A + μ F 1 ) μ F 2 Q 0 > 0 . Hence, R13k>0 , k = 0 , 1. Since ak2ak3−ak1>0 , k = 1 , 2, the analysis of R14k restricts in the analysis of whether ak1ak2ak3−ak012-ak32ak0>0 , k = 0 , 1. For k = 0 the expression reads a 01 a 02 a 03 - a 01 2 - a 03 2 a 00 = α 1 + μ E 3 α 2 + μ A 2 α 3 + μ F 1 + μ F 2 + α 1 + μ E 3 α 3 + μ F 1 2 α 2 + μ A + μ F 2 + α 1 + μ E 3 μ F 2 2 α 2 + α 3 + μ A + μ F 1 + 2 α 1 + μ E 3 a 2 + μ A α 3 + μ F 1 μ F 2 + α 1 + μ E 2 a 2 + μ A 3 α 3 + μ F 1 + μ F 2 + α 2 + μ A 3 a 3 + μ F 1 2 α 1 + μ E + μ F 2 + α 2 + μ A 3 μ F 2 2 α 1 + α 3 + μ E + μ F 1 + 2 α 1 + μ E α 2 + μ A 3 α 3 + μ F 1 μ F 2 + α 1 + μ E 2 α 3 + μ F 1 3 α 2 + μ A + μ F 2 + α 2 + μ A 2 α 3 + μ F 1 3 α 1 + μ E + μ F 2 + α 3 + μ F 1 3 μ F 2 2 ( α 1 + α 2 + μ E + μ A + 2 α 1 + μ E α 2 + μ A α 3 + μ F 1 3 μ F 2 + α 1 + μ E 2 μ F 2 3 α 2 + α 3 + μ A + μ F 1 + α 2 + μ A 2 μ F 2 3 α 1 + α 3 + μ E + μ F 1 + α 3 + μ F 1 2 μ F 2 3 α 1 + α 2 + μ E + μ A + 2 α 1 + μ E α 2 + μ A α 3 + μ F 1 μ F 2 3 + 4 α 1 + μ E 2 α 2 + μ A α 3 + μ F 1 μ F 2 α 2 + μ A + μ F 2 + 4 α 1 + μ E α 2 + μ A 2 α 3 + μ F 1 μ F 2 α 3 + μ F 1 + μ F 2 + 4 α 1 + μ E α 2 + μ A α 3 + μ F 1 2 μ F 2 α 3 + μ E + μ F 2 + ϕ α 1 α 2 α 3 α 1 + α 2 + α 3 + μ A + μ E + μ F 1 + μ F 2 2 > 0 . The expression is positive because all parameters of the model are positive. In the case k = 1, it reads a 11 a 12 a 13 - a 11 2 - a 10 a 13 2 = α 1 + μ E 3 α 2 + μ A 2 α 3 + μ F 1 + μ F 2 Q 0 3 + α 1 + μ E 3 α 3 + μ F 1 2 α 2 + μ A + μ F 2 Q 0 3 + α 1 + μ E 3 μ F 2 2 α 2 + α 3 + μ A + μ F 1 Q 0 3 + α 1 + μ E 3 a 2 + μ A α 3 + μ F 1 μ F 2 Q 0 2 2 Q 0 + 1 + 2 α 1 + μ E 2 a 2 + μ A 2 α 3 + μ F 1 2 Q 0 2 + 2 α 1 + μ E 2 a 2 + μ A 2 μ F 2 2 Q 0 2 + 2 α 1 + μ E 2 a 3 + μ F 1 2 μ F 2 2 Q 0 2 + 4 α 1 + μ E 2 α 2 + μ A α 3 + μ F 1 μ F 2 α 2 + α 3 + μ A + μ F 1 + μ F 2 Q 0 2 + α 1 + μ E 2 α 2 + μ A 3 α 3 + μ F 1 + μ F 2 Q 0 2 + α 1 + μ E 2 α 3 + μ F 1 3 α 2 + μ A + μ F 2 Q 0 2 + α 1 + μ E 2 μ F 2 3 α 2 + α 3 + μ A + μ F 1 Q 0 2 + α 1 + μ E α 2 + μ A 3 α 3 + μ F 1 + μ F 2 2 Q 0 + α 1 + μ E α 3 + μ F 1 3 α 2 + μ A + μ F 2 2 Q 0 + α 1 + μ E μ F 2 3 α 2 + α 3 + μ A + μ F 1 2 Q 0 + 2 α 1 + μ E α 2 + μ A 2 α 3 + μ F 1 2 μ F 2 2 Q 0 + 1 + 2 α 1 + μ E α 2 + μ A 2 α 3 + μ F 1 μ F 2 2 2 Q 0 + 1 + 2 α 1 + μ E α 2 + μ A α 3 + μ F 1 2 μ F 2 2 2 Q 0 + 1 + α 2 + μ A 3 α 3 + μ F 1 μ F 2 α 1 + α 3 + μ E + μ F 1 + μ F 2 + α 2 + μ A α 3 + μ F 1 3 μ F 2 α 1 + α 2 + μ E + μ A + μ F 2 + α 2 + μ A α 3 + μ F 1 μ F 2 3 α 1 + α 2 + α 3 + μ E + μ A + μ F 1 + 2 α 2 + μ A 2 α 3 + μ F 1 2 μ F 2 2 + 2 ϕ α 1 α 2 α 3 α 1 + μ E α 2 + α 3 + μ A + μ F 1 + μ F 2 > 0 . It is also positive because all parameters of the model are positive. Finally, for k = 0, R150=a00=α1+μEα2+μAα3+μF1μF21-Q0 and for k = 1, R151=a10=ϕα1α2α31-1Q0 . In both cases, the signal of the matrix element will depend on Q 0to be 0<Q0<1 or Q0>1 . Moreover, the factors 1-Q0 and 1-1Q0 have opposite signals in those cases. Therefore, in the case 0<Q0<1,R150>0 and the first column of the Routh Table for k = 0 will be strictly positive, causing P 0to be stable. Yet, for k = 1, in the same case, R151<0 , and there will have a change of sign in the last element of the first column of the Routh Table, causing P 1to be unstable. Now, turning to the case Q0>1 , we now have R150<0 for k = 0 and R151>0 for k > 1. Hence the change of sign in the first column of the Routh Tables switches and now P 0becomes unstable and P 1becomes stable. In the case Q 0= 1, R15k=0 for k = 0; 1, P 0and P 1are coincident and trivial,the independent terms a k0 of the characteristic polynomials ΦPkλ , k = 0;1 vanish and there take place a null eigenvalue, λ = 0. This fact together with the change of stability in ρ = 1 is sufficient to characterise the point as an transcritical bifurcation point. ) are degenerated and trivial, are P1=P0 , if Q0=1 . The point P 0 is asymptotically stable and P 1, unstable, if 0<Q0<1 , while P 1 turns to be asymptotically stable and P 0 unstable, if Q0>1 . The system shows transcritical bifurcation at Q0=1 2323 S.H. Strogatz. “Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering”. CRC Press (2018). .

As the meteorological parameters vary weekly along the ranges TminTTmax×0ppmax , the model parameters (3.2) vary. Therefore, Q 0 has different values each week but, according the the parameters defined, it remains positive Q03.9187 along the range of ( p, T ) data. Consequently, the non trivial equilibrium point P 1 varies, but remains asymptotically stable. In particular, the equilibrium point coordinate F2** values may be represented by a surface according to p and T (see Fig. 4 ).

Figure 4:
Surface of the females coordinate F2**=F2**p,T of the non trivial asymptotically stable equilibrium point, according to meteorological variables. The dependence on temperature reaches a maximum at a optimum temperature and increases monotonically with increasing precipitation.

4 NUMERICAL SOLUTION AND MODEL VERIFICATION

Model (3.1) was solved numerically with a 4th order Runge-Kutta Method using Matlab® software with time step 0.001, corresponding to 7000 points per week of data. In the evaluation of the results, we are interested only in the evolution of the population F 2( t ), since the data refer to the capture of pregnant females who seek the trap to make oviposition. The model is evaluated, minimising the mean square error given by:

S 2 = 1 2 N k = 1 N I k - λ F 2 k - l 2 , (4.1)

where N is the number of epidemiological weeks of the data and I ( k ) is the value of MFAI data in the week k . Since the weekly monitoring of traps takes time, the lag between model and data is given by l . The scale of the model is arbitrarily related to the environmental capacity, so the λ is a scale factor to adjust the scale to the sample capture data (MFAI). Subject to condition S2λ=0 , one finds

λ l = k = 1 N I k F 2 k - l k = 1 N F 2 k - l 2 . (4.2)

By setting l max as the value for which the cross-correlation between MFAI series I ( k ) and F 2 reaches its maximum, the scale factor between data and model can be fixed so that λ=λlmax . The value of parameter r in equation (3.3) was determined empirically by a simple one-dimensional optimisation as the process is repeated for values or r such that 0<r2 , with step of 0 . 1 searching for the lowest value of S 2. In the present case, it was obtained for r = 1 . 2, value used in the simulation. For simplicity, the carrying capacity C ( p, T ) of the medium was set to be constant, C = 1, defined as relative to the best environmental conditions, when all existing water containers are available, filled of water and larval nutrients.

To evaluate the result, it was used cross-correlation between the MFAI series I ( k ) and model female population corrected by the scale factor f2=λF2k (see Figure 5a ). The maximum of the correlation is 0.75 at lag l max = 6. The cross correlation between data and model remains significant ( p < 0 . 05) for lags ranging from 0 to 14. Figure 5b depicts a comparison between MFAI and f2=k-6 . The simple model reproduces with certain quality the field data of captures. More complex details like peaks are not fitted in general.

Figure 5:
(a) Cross-correlation between model population f 2 and infestation index MFAI. The crossing line stands for significant level. Maximum correlation occurs for lag of l = 6 weeks (p¡0.05) (b) Plots of the infestation index MFAI together with the lagged, l = 6, and scaled f 2 model population.

5 CONCLUSIONS

After the comparison between model and data of female captures, one can conclude the curves are similar, indicating that the proposed model is able to reproduce the data in an acceptable manner. y the population of F 2( t ). In some intervals, the model population becomes larger than the MFAI indicator, and, in others, the opposite occurs. Although the modelled peaks are not as pronounced as the infestation data, it is possible to notice their existence. The provided infestation data is obtained from a a public health program and by this reason, it is subjected to entomological control with insecticides. However, with few simple hypotheses on the dependence of the entomological parameters with temperature and precipitation, the model was able to reproduce the baseline of the infestation. Although the determination of dependencies of model rates on temperature and precipitation were subjected to simplifying assumptions and, therefore, of limited realism, the capability of the model to reproduce and forecast the data according to meteorological parameters is good. The task to evaluate the approach of those dependencies on meteorological data was successful. Eventual future refinement on those parametric dependencies with more realistic temperature parameters could provide more accuracy. Since the modelled precipitation parameters are set only as order of magnitude, an approach via optimisation could search for candidates for the values of the dependence of the entomological rates of development with precipitation. The present parameters can be used as initial set of parameters for optimisation. The result of the optimised parameters could guide biological experiments. The model is able to predict infestation according to meteorological parameters in public heath control programs.

Acknowledgments

The authors express their gratitude to the financial support of the Minas Gerais State Research Support Foundation (FAPEMIG), to the Federal Centre for Technological Education of Minas Gerais (CEFET), and Ecovec S.A.

REFERENCES

  • 1
    E.B. Beserra, C.R. Fernandes, S.A.D.O. Silva, L.A.D. Silva & J.W.D. Santos. Efeitos da temperatura no ciclo de vida, exigências térmicas e estimativas do número de gerações anuais de Aedes aegypti (Diptera, Culicidae). Iheringia. Série Zoologia, 99 (2009), 142-148.
  • 2
    M. Burattini, M. Chen, A. Chow, F. Coutinho, K. Goh, L. Lopez, S. Ma & E. Massad. Modelling the control strategies against dengue in Singapore. Epidemiology & Infection, 136(3) (2008), 309-319.
  • 3
    D.D. Chadee, F.L. Williams & U.D. Kitron. Impact of vector control on a dengue fever outbreak in Trinidad, West Indies, in 1998. Tropical Medicine & International Health, 10(8) (2005), 748-754.
  • 4
    S.R. Christophers. “Aedes aegypti: the yellow fever mosquito”. CUP Archive (1960).
  • 5
    E.A.P.D.A. Costa, E.M.D.M. Santos, J.C. Correia & C.M.R.D. Albuquerque. Impact of small variations in temperature and humidity on the reproductive activity and survival of Aedes aegypti (Diptera, Culicidae). Revista Brasileira de Entomologia, 54(3) (2010), 488-493.
  • 6
    L. Esteva & H.M. Yang. Mathematical model to assess the control of Aedes aegypti mosquitoes by the sterile insect technique. Mathematical biosciences, 198(2) (2005), 132-147.
  • 7
    D.A. Focks, D.G. Haile, E. Daniels & G.A. Mount. Dynamic life table model for Aedes aegypti (Diptera: Culicidae): analysis of the literature and model development. Journal of medical entomology, 30(6) (1993), 1003-1017.
  • 8
    FUNASA. Programa Nacional de Controle da Dengue. Ministério da Saúde Vigilância Epidemiológica, (2002).
  • 9
    A.D.S. Gomes, C.J. Sciavico & Á .E. Eiras. Periodicity of oviposition of females of Aedes aegypti (Linnaeus, 1762)(Diptera: Culicidae) in laboratory and field. Revista da Sociedade Brasileira de Medicina Tropical, 39(4) (2006), 327-332.
  • 10
    D.J. Gubler. The global emergence/resurgence of arboviral diseases as public health problems. Archives of medical research, 33(4) (2002), 330-342.
  • 11
    M.U.G. Kraemer, M.E. Sinka, K.A. Duda, A.Q. Mylne, F.M. Shearer, C.M. Barker, C.G. Moore, R.G. Carvalho, G.E. Coelho, W. Van Bortel et al. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. elife, 4 (2015), e08347.
  • 12
    G. Kuno. Review of the factors modulating dengue transmission. Epidemiologic reviews, 17(2) (1995), 321-335.
  • 13
    A.P.M. Moreira, V.A.L. C. & S.J. C. Tendências climáticas e anomalias de precipitação em CaratingaMG. In “Os Desafios da Geografia Física na Fronteira do Conhecimento” (2018), pp. 2000-2009.
  • 14
    V.S.G. Neto & J.M.M. Rebêlo. Aspectos epidemiológicos do dengue no Município de São Luís, Maranhão, Brasil, 1997-2002. Cad. Saúde Pública, 20(5) (2004), 1424-1431.
  • 15
    D.P. Neves. “Parasitologia humana”. Atheneu (2005).
  • 16
    S. Noguchi, A.R. Nik & M. Tani. Rainfall characteristics of tropical rainforest at Pasoh forest reserve, Negeri Sembilan, Peninsular Malaysia. In “Pasoh”. Springer (2003), pp. 51-58.
  • 17
    M. Otero, H.G. Solari & N. Schweigmann. A stochastic population dynamics model for Aedes aegypti: formulation and application to a city with temperate climate. Bulletin of mathematical biology, 68(8) (2006), 1945-1974.
  • 18
    A.F. Ribeiro, G.R.A.M. Marques, J.C. Voltolini & M.L.F. Condino. Associação entre incidência de dengue e variáveis climáticas. Rev Saúde Pública, 40(4) (2006), 671-6.
  • 19
    T.W. Scott, P.H. Amerasinghe, A.C. Morrison, L.H. Lorenz, G.G. Clark, D. Strickman, P. Kittayapong & J.D. Edman. Longitudinal studies of Aedes aegypti (Diptera: Culicidae) in Thailand and Puerto Rico: blood feeding frequency. Journal of medical entomology, 37(1) (2000), 89-101.
  • 20
    L.D.S.B. Silva, R.T.N. Cardoso, J.L.A. Fernandes, C.A. Silva & A.E. Eiras. Modelo Entomológico Determinístico sob Efeito da Pluviosidade para o Aedes aegypti e o Aedes albopictus. TEMA (São Carlos), 19(2) (2018), 289-303.
  • 21
    G. Smith, D. Eliason, C. Moore & E. Ihenacho. Use of elevated temperatures to kill Aedes albopictus and Ae. aegypti. Journal of the American Mosquito Control Association, 4(4) (1988), 557-558.
  • 22
    R. Souza-Santos & M.S. Carvalho. Análise da distribuição espacial de larvas de Aedes aegypti na Ilha do Governador, Rio de Janeiro, Brasil. Cadernos de saúde Pública, 16(1) (2000), 31-42.
  • 23
    S.H. Strogatz. “Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering”. CRC Press (2018).
  • 24
    A. Tran, G. L’Ambert, G. Lacour, R. Benôıt, M. Demarchi, M. Cros, P. Cailly, M. Aubry-Kientz, T. Balenghien & P. Ezanno. A rainfall-and temperature-driven abundance model for Aedes albopictus populations. International journal of environmental research and public health, 10(5) (2013), 1698-1719.
  • 25
    J. Waldock, N.L. Chandra, J. Lelieveld, Y. Proestos, E. Michael, G. Christophides & P.E. Parham. The role of environmental variables on Aedes albopictus biology and chikungunya epidemiology. Pathogens and global health, 107(5) (2013), 224-241.
  • 26
    S. Wiggins. “Introduction to applied nonlinear dynamical systems and chaos”, volume 2. Springer Science & Business Media (2003).
  • 27
    World Health Organization. Dengue and severe dengue. http://www.who.int/mediacentre/factsheets/fs117/en/ (2017). Acessado em 20/03/2018.
    » http://www.who.int/mediacentre/factsheets/fs117/en/
  • 28
    World Health Organization. Neglected tropical diseases. http://www.searo.who.int/entity/vector_borne_tropical_diseases/data/data_factsheet/en/ (2018). Acessado em 10/09/2018.
    » http://www.searo.who.int/entity/vector_borne_tropical_diseases/data/data_factsheet/en/
  • 29
    H.M. Yang. The transovarial transmission in the dynamics of dengue infection: Epidemiological implications and thresholds. Mathematical biosciences, 286 (2017), 1-15.
  • 30
    H.M. Yang, M.L.G. Macoris, K.C. Galvani, M.T.M. Andrighetti & D.M.V. Wanderley. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and infection, 137(08) (2009), 1188-1202.

Appendix A STABILITY OF THE EQUILIBRIUM POINTS

Let the system of equations (3.1) rewritten as

d X d t = f X , where X = E , A , F 1 , F 2 . (A.1)

The condition fX=0 revealed two critical points, one trivial, P 0(Equation 3.5), and one nontrivial, P 1(Equation 3.6). In the following it is performed the analysis of stability of the system.

After the expansion near the equilibrium points P k ; k = 0;1 and computation of the Jacobian matrices BPk=fXPk , k = 0 , 1, their characteristic polynomials are given, respectively, by

Φ P k λ = a k 4 λ 4 + a k 3 λ 3 + a k 2 λ 2 + a k 1 λ + a k 0 , (A.2)

with coefficients for P 0,

a 04 = 1 ,

a 03 = α 1 + μ E + α 2 + μ A + α 3 + μ F 1 + μ F 2 ,

a 02 = α 1 + μ E μ F 2 + α 2 + α 3 + μ A + μ F 1 α 1 + μ E + μ F 2 + α 2 + μ A α 3 + μ F 1 ,

a 01 = α 1 + μ E α 2 + μ A α 3 + μ F 1 + μ F 2 + α 1 + α 2 + μ E + μ A α 3 + μ F 1 μ F 2 ,

a 00 = α 1 + μ E α 2 + μ A α 3 + μ F 1 μ F 2 1 - Q 0 ;

and for P 1,

a 14 = 1 ,

a 13 = ( α 1 + μ E ) Q 0 + α 2 + α 3 + μ A + μ F 1 + μ F 2 ,

a 12 = ( α 1 + μ E ) ( α 2 + α 3 + μ A + μ F 1 + μ F 2 ) Q 0 + ( α 2 + μ A ) ( α 3 + μ F 1 + μ F 2 ) + ( α 3 + μ F 1 ) μ F 2 ,

a 11 = ( α 1 + μ E ) ( α 2 + μ A ) ( α 3 + μ F 1 + μ F 2 ) Q 0 + ( ( α 1 + μ E ) Q 0 + α 2 + μ A ) ( α 3 + μ F 1 ) μ F 2 ,

a 10 = ϕ α 1 α 2 α 3 ( 1 1 Q 0 ) .

The stability of each equilibrium points P k can be assessed from the sign of the real part of their respective eigenvalues. The polynomial ΦPk (A.2) with have roots with real part strictly less than zero if and only if all elements in the first column of the Routh Table (A.3) are nonzero and have the same sign 2626 S. Wiggins. “Introduction to applied nonlinear dynamical systems and chaos”, volume 2. Springer Science & Business Media (2003). .

R k = 1 a k 2 a k 0 a k 3 a k 1 0 a k 2 a k 3 - a k 1 a k 3 a k 0 0 a k 1 - a k 3 2 a k 0 a k 2 a k 3 - a k 1 0 0 a k 0 0 0 (A.3)

The elements R11k=1>0 , k = 0 , 1 are positive. Since all model rates remain positive by definition, R21k>0 , k = 0 , 1 is positive. As consequence, R31k , k = 0 , 1, are positive if ak2ak3-ak1>0 . For k = 0, it reads

a 02 a 03 a 01 = 2 ( α 1 + μ E ) ( α 2 + α 3 + μ A + μ F 1 ) μ F 2 + 2 ( α 2 + μ A ) ( α 3 + μ F 1 ) ( α 1 + μ E + μ F 2 ) + ( α 1 + μ E ) 2 ( α 2 + α 3 + μ A + μ F 1 + μ F 2 ) + ( α 2 + μ A ) 2 ( α 1 + α 3 + μ E + μ F 1 + μ F 2 ) + ( α 3 + μ F 1 ) 2 ( α 1 + α 2 + μ E + μ A + μ F 2 ) + ( μ F 2 ) 2 ( α 1 + α 2 + α 3 + μ E + μ A + μ F 1 ) > 0 .

In the case k = 1, it is expressed as

a 12 a 13 a 11 = ( α 1 + μ E ) 2 ( α 2 + α 3 + μ A + μ F 1 + μ F 2 ) Q 0 2 + 2 ( α 2 + μ A ) ( α 3 + μ F 1 ) μ F 2 + ( α 1 + μ E ) ( α 2 + α 3 + μ A + μ F 1 + μ F 2 ) ( α 2 + α 3 + μ A + μ F 1 ) Q 0 + ( α 2 + μ A ) 2 ( α 3 + μ F 1 + μ F 2 ) + μ F 2 2 ( ( α 1 + μ E ) Q 0 + α 2 + α 3 + μ A + μ F 1 ) + ( α 3 + μ F 1 ) 2 ( α 2 + μ A + μ F 2 ) + ( α 1 + μ E ) ( α 2 + α 3 + μ A + μ F 1 ) μ F 2 Q 0 > 0 .

Hence, R13k>0 , k = 0 , 1. Since ak2ak3ak1>0 , k = 1 , 2, the analysis of R14k restricts in the analysis of whether ak1ak2ak3ak012-ak32ak0>0 , k = 0 , 1. For k = 0 the expression reads

a 01 a 02 a 03 - a 01 2 - a 03 2 a 00 = α 1 + μ E 3 α 2 + μ A 2 α 3 + μ F 1 + μ F 2 + α 1 + μ E 3 α 3 + μ F 1 2 α 2 + μ A + μ F 2 + α 1 + μ E 3 μ F 2 2 α 2 + α 3 + μ A + μ F 1 + 2 α 1 + μ E 3 a 2 + μ A α 3 + μ F 1 μ F 2 + α 1 + μ E 2 a 2 + μ A 3 α 3 + μ F 1 + μ F 2 + α 2 + μ A 3 a 3 + μ F 1 2 α 1 + μ E + μ F 2 + α 2 + μ A 3 μ F 2 2 α 1 + α 3 + μ E + μ F 1 + 2 α 1 + μ E α 2 + μ A 3 α 3 + μ F 1 μ F 2 + α 1 + μ E 2 α 3 + μ F 1 3 α 2 + μ A + μ F 2 + α 2 + μ A 2 α 3 + μ F 1 3 α 1 + μ E + μ F 2 + α 3 + μ F 1 3 μ F 2 2 ( α 1 + α 2 + μ E + μ A + 2 α 1 + μ E α 2 + μ A α 3 + μ F 1 3 μ F 2 + α 1 + μ E 2 μ F 2 3 α 2 + α 3 + μ A + μ F 1 + α 2 + μ A 2 μ F 2 3 α 1 + α 3 + μ E + μ F 1 + α 3 + μ F 1 2 μ F 2 3 α 1 + α 2 + μ E + μ A + 2 α 1 + μ E α 2 + μ A α 3 + μ F 1 μ F 2 3 + 4 α 1 + μ E 2 α 2 + μ A α 3 + μ F 1 μ F 2 α 2 + μ A + μ F 2 + 4 α 1 + μ E α 2 + μ A 2 α 3 + μ F 1 μ F 2 α 3 + μ F 1 + μ F 2 + 4 α 1 + μ E α 2 + μ A α 3 + μ F 1 2 μ F 2 α 3 + μ E + μ F 2 + ϕ α 1 α 2 α 3 α 1 + α 2 + α 3 + μ A + μ E + μ F 1 + μ F 2 2 > 0 .

The expression is positive because all parameters of the model are positive. In the case k = 1, it reads

a 11 a 12 a 13 - a 11 2 - a 10 a 13 2 = α 1 + μ E 3 α 2 + μ A 2 α 3 + μ F 1 + μ F 2 Q 0 3 + α 1 + μ E 3 α 3 + μ F 1 2 α 2 + μ A + μ F 2 Q 0 3 + α 1 + μ E 3 μ F 2 2 α 2 + α 3 + μ A + μ F 1 Q 0 3 + α 1 + μ E 3 a 2 + μ A α 3 + μ F 1 μ F 2 Q 0 2 2 Q 0 + 1 + 2 α 1 + μ E 2 a 2 + μ A 2 α 3 + μ F 1 2 Q 0 2 + 2 α 1 + μ E 2 a 2 + μ A 2 μ F 2 2 Q 0 2 + 2 α 1 + μ E 2 a 3 + μ F 1 2 μ F 2 2 Q 0 2 + 4 α 1 + μ E 2 α 2 + μ A α 3 + μ F 1 μ F 2 α 2 + α 3 + μ A + μ F 1 + μ F 2 Q 0 2 + α 1 + μ E 2 α 2 + μ A 3 α 3 + μ F 1 + μ F 2 Q 0 2 + α 1 + μ E 2 α 3 + μ F 1 3 α 2 + μ A + μ F 2 Q 0 2 + α 1 + μ E 2 μ F 2 3 α 2 + α 3 + μ A + μ F 1 Q 0 2 + α 1 + μ E α 2 + μ A 3 α 3 + μ F 1 + μ F 2 2 Q 0 + α 1 + μ E α 3 + μ F 1 3 α 2 + μ A + μ F 2 2 Q 0 + α 1 + μ E μ F 2 3 α 2 + α 3 + μ A + μ F 1 2 Q 0 + 2 α 1 + μ E α 2 + μ A 2 α 3 + μ F 1 2 μ F 2 2 Q 0 + 1 + 2 α 1 + μ E α 2 + μ A 2 α 3 + μ F 1 μ F 2 2 2 Q 0 + 1 + 2 α 1 + μ E α 2 + μ A α 3 + μ F 1 2 μ F 2 2 2 Q 0 + 1 + α 2 + μ A 3 α 3 + μ F 1 μ F 2 α 1 + α 3 + μ E + μ F 1 + μ F 2 + α 2 + μ A α 3 + μ F 1 3 μ F 2 α 1 + α 2 + μ E + μ A + μ F 2 + α 2 + μ A α 3 + μ F 1 μ F 2 3 α 1 + α 2 + α 3 + μ E + μ A + μ F 1 + 2 α 2 + μ A 2 α 3 + μ F 1 2 μ F 2 2 + 2 ϕ α 1 α 2 α 3 α 1 + μ E α 2 + α 3 + μ A + μ F 1 + μ F 2 > 0 .

It is also positive because all parameters of the model are positive.

Finally, for k = 0, R150=a00=α1+μEα2+μAα3+μF1μF21-Q0 and for k = 1, R151=a10=ϕα1α2α31-1Q0 . In both cases, the signal of the matrix element will depend on Q 0to be 0<Q0<1 or Q0>1 . Moreover, the factors 1-Q0 and 1-1Q0 have opposite signals in those cases. Therefore, in the case 0<Q0<1,R150>0 and the first column of the Routh Table for k = 0 will be strictly positive, causing P 0to be stable. Yet, for k = 1, in the same case, R151<0 , and there will have a change of sign in the last element of the first column of the Routh Table, causing P 1to be unstable. Now, turning to the case Q0>1 , we now have R150<0 for k = 0 and R151>0 for k > 1. Hence the change of sign in the first column of the Routh Tables switches and now P 0becomes unstable and P 1becomes stable. In the case Q 0= 1, R15k=0 for k = 0; 1, P 0and P 1are coincident and trivial,the independent terms a k0 of the characteristic polynomials ΦPkλ , k = 0;1 vanish and there take place a null eigenvalue, λ = 0. This fact together with the change of stability in ρ = 1 is sufficient to characterise the point as an transcritical bifurcation point.

Publication Dates

  • Publication in this collection
    05 Apr 2021
  • Date of issue
    Jan-Mar 2021

History

  • Received
    12 Sept 2018
  • Accepted
    15 Oct 2020
Sociedade Brasileira de Matemática Aplicada e Computacional - SBMAC Rua Maestro João Seppe, nº. 900, 16º. andar - Sala 163, Cep: 13561-120 - SP / São Carlos - Brasil, +55 (16) 3412-9752 - São Carlos - SP - Brazil
E-mail: sbmac@sbmac.org.br