Modelling entomological-climatic interactions of Plasmodium falciparum malaria transmission in two Colombian endemic-regions: contributions to a National Malaria Early Warning System
- Daniel Ruiz^{1, 2}Email author,
- Germán Poveda^{1},
- Iván D Vélez^{3},
- Martha L Quiñones^{3, 4},
- Guillermo L Rúa^{3},
- Luz E Velásquez^{3} and
- Juan S Zuluaga^{5}
https://doi.org/10.1186/1475-2875-5-66
© Ruiz et al;licensee BioMed Central Ltd. 2006
Received: 27 March 2006
Accepted: 01 August 2006
Published: 01 August 2006
Abstract
Background
Malaria has recently re-emerged as a public health burden in Colombia. Although the problem seems to be climate-driven, there remain significant gaps of knowledge in the understanding of the complexity of malaria transmission, which have motivated attempts to develop a comprehensive model.
Methods
The mathematical tool was applied to represent Plasmodium falciparum malaria transmission in two endemic-areas. Entomological exogenous variables were estimated through field campaigns and laboratory experiments. Availability of breeding places was included towards representing fluctuations in vector densities. Diverse scenarios, sensitivity analyses and instabilities cases were considered during experimentation-validation process.
Results
Correlation coefficients and mean square errors between observed and modelled incidences reached 0.897–0.668 (P > 0.95) and 0.0002–0.0005, respectively. Temperature became the most relevant climatic parameter driving the final incidence. Accordingly, malaria outbreaks are possible during the favourable epochs following the onset of El Niño warm events. Sporogonic and gonotrophic cycles showed to be the entomological key-variables controlling the transmission potential of mosquitoes' population. Simulation results also showed that seasonality of vector density becomes an important factor towards understanding disease transmission.
Conclusion
The model constitutes a promising tool to deepen the understanding of the multiple interactions related to malaria transmission conducive to outbreaks. In the foreseeable future it could be implemented as a tool to diagnose possible dynamical patterns of malaria incidence under several scenarios, as well as a decision-making tool for the early detection and control of outbreaks. The model will be also able to be merged with forecasts of El Niño events to provide a National Malaria Early Warning System.
Background
The World Health Organization estimates that malaria parasites infect from 200 to 300 million persons and kill more than 1 million people each year, primarily children under the age of five [1]. The efficacy of control measures has decreased over the past decades because mosquitoes and parasites are becoming more resistant to the commonly used insecticides and anti-malarial drugs [2]. As a result, malaria kills more people today than three decades ago. To further complicate matters, the international community suggests that diseases relayed by mosquitoes, such as malaria, are among those infectious diseases most likely to spread dramatically as global temperatures head upward [3, 4]. Forecasts from different models suggest that by the end of the 21st century ongoing warming will have enlarged the zone of potential malaria transmission from an area containing 45 percent of the world's population to an area containing about 60 percent [4]. Although this expansion certainly fits the predictions, the cause of that growth may not be attributed convincingly to global warming. Other factors (maybe the dominant contributors) may have been involved as well, for instance: disruption of the environment in ways that favour the mosquitoes' proliferation, declines in vector-control and in other public health programs, and rises in drug and insecticide resistance [4]. To make matters worse, there are a number of 'macro-factors' such as economic inequalities, continuous human migratory patterns and vegetation patterns that may increase or decrease vulnerability and exposure to vector-borne infections and play an important role in regional morbidity and mortality profiles. Clearly, it is impossible to treat climate isolated from the other biological, socioeconomic and demographic contributors if the determination of the potential impact of climate change on malaria incidence is urgently need to be made. However, the hypothesis for a climatic contribution becomes stronger: even though malaria is a highly complex multi-factorial disease, previous studies have identified environmental factors and climate variability as going a considerable way in helping to explain the fluctuations of disease incidence [5–7].
In Colombia, currently, the population exceeds 44 million people; more than five million live in endemic-prone regions. Malaria has recently re-emerged as a significant public health burden in the country: incidence during epidemic years (total positive cases for both Plasmodium falciparum and Plasmodium vivax malaria scaled by the total population at risk per 1,000 inhabitants) increased from less than 2.5 in 1963 to 6 in 1983 and almost 10 in 1998. During 1996, disease transmission reached 42 cases per 1,000 inhabitants in high-risk areas. In the Chocó Department, on the Colombian Pacific Coast, more than 80,000 cases were reported during 1998, when the population at risk was 380,000 people [8]. During the year 2000 malaria cases were reported in 69% of the 125 towns of the Department of Antioquia, in the Andean region. Besides the aforementioned strongly increasing trend, the Annual Parasite Index also exhibited a significant association between the increase in the number of malaria cases and the occurrence of the El Niño warm event, which is considered the main forcing mechanism of Colombia's hydroclimatology at inter-annual timescales [8]. There is strong evidence that the El Niño event intensifies the annual cycle of malaria cases in endemic rural areas as a consequence of concomitant anomalies in the normal annual cycle of temperature and precipitation. During 'normal' years, endemic malaria exhibits a clear-cut 'normal' annual cycle, which is associated with prevalent climatic conditions. During El Niño events, characterized by increasing temperatures and decreasing monthly rainfalls, river flows and soil moistures, malaria outbreaks were found to be enhanced by climate anomalies. The annual cycle of malaria seems to be affected in its amplitude (increase in the number of cases), although the phase (or timing) remains unaltered in most malaria prone-regions of the country [8].
As recent scientific and technological advances have permitted much better predictions of the El Niño event, there is an evident opportunity to incorporate this climate forecasting capability into programs, campaigns, control measures and mitigation plans to reduce the human health impact of malaria outbreaks in Colombian prone-regions. Nevertheless, there remain significant gaps of knowledge in the understanding of the interactions between climatic factors and the dynamics of malaria transmission, which have motivated attempts to explore and develop mathematical comprehensive models.
Detailed dynamic models of malaria epidemiology have been developed extensively and have made important contributions to understanding the transmission of malaria and other diseases [3, 9–21]. There are limitations, however. Although these models can provide adequate approximations to some biological and epidemiological characteristics, most of them have not been able to describe the overall transmission dynamics [21]. Even though mosquitoes and vertebrate host dynamics are generally included, most models do not include changes in environmental and climatic patterns that could affect disease incidence in both spatial and temporal scales, nor do they include the social, economical and demographical conditions prevailing in the communities that were studied. In this article experiences in developing and implementing a vector-borne disease model, which is based on previously proposed differential-equation 'compartment' models, are reported. The aim was to validate the mathematical tool by comparing its results under conditions of intense transmission in two specific malaria prone-regions of Colombia. To do so, malaria prevalence changes were monitored in these areas and model predictions were tested to assess their consistency with field observations.
In developing the mathematical model, three major steps were followed: (i) First, the parasite transmission cycle was studied to define those endogenous variables strongly affected by climatic conditions. (ii) The vector ecology, behavioural patterns of mosquitoes and entomological parameters were analysed in order to determine those exogenous variables (and their relationships with climatic anomalies) that might be relevant for representing malaria transmission. Analyses were focused on the necessity to represent the observed fluctuations in vector density. (iii) The vector population dynamics during pre-imago stages was represented, including the availability of adequate breeding sites and the predator-prey interactions during larval stage.
In applying and validating the mathematical tool, five major steps were followed: (a) homogeneity analysis of hydrological time series, to detect possible changes in mean and variance, as well as significant trends in historical series of both temperature and precipitation. (b) Analysis of total positive malaria cases, to determine the behaviour of epidemiological time series, to establish the initial values of state variables, and for fine-tuning the parameters. (c) Correlation analysis, to determine the contribution of climate factors on the temporal variance of malaria transmission. (d) Simulation-base scenario, to analyse the behaviour of the proposed differential-equation model and its preliminary results. And (e) simulation-alternate scenarios, to analyse several simulation results of the set of non-linear differential equations when perturbations of several exogenous variables, changes in control parameters and changing climate scenarios were considered.
The mathematical tool allowed improving the understanding of the linkages between climate patterns and malaria outbreaks, mainly during the onset of ENSO warm events. Experiences suggest that forecasts of future El Niño events will be able to be merged with mathematical models to provide a Malaria Early Warning System (MEWS) to enhance disease surveillance, control and response to epidemics, and facilitate early, coupled and environmentally sound public health interventions [8].
Methods
Main module
By ignoring the secondary exo-erythrocytic schizogony, that takes place in the parenchymal cells of the liver, simulation was limited to representing P. falciparum malaria transmission. Based upon the widely known biological and epidemiological aspects of malaria transmission, a susceptible person may become a member of the infected class if the individual is inoculated with sporozoites during an adult female Anopheles mosquito infective bite. The exo-erythrocytic schizogony and the erythrocytic schizogony take place in the infected vertebrate host. During this stage, sexually differentiated forms (gametocytes) are present in the host's blood circulation and can be taken by a susceptible female mosquito during a new blood meal. Parasites can then be transmitted from this mosquito to a susceptible person in some latter human-vector interaction. The progressive increase of parasitaemia in the infectious host as a result of multiple erythrocytic cycles may be slowed down by the immune response of the human. Finally, the immune individual loses immunity and returns to the reservoir of susceptible hosts.
In the vector ecology, a susceptible adult female Anopheles mosquito becomes a member of the infected class when it takes a blood meal from an infectious host and ingests the gametocytes present in the individual's blood stream. The vector finally enters the reservoir of infectious mosquitoes after the completion of the sporogonic cycle.
Since vector density is commonly assumed to be constant by most researchers, this variable has played a minor role in the construction of their mathematical infectious disease models. It is generally accepted that it is very difficult to arrive at a reliable estimate of the changes in vector density over large areas as a result of changes in temperature, precipitation and humidity. The effort of this research is aimed at fully representing the vector density in human-mosquito interactions, by considering three additional stages (virtual compartments) in the mosquito ecology: the pre-imago stages eggs, larvae, and pupae. Interests included the analysis of average development periods from eggs to adult mosquitoes, expressed as time delays and the related survival probabilities, expressed as successes in hatching, larval development, and adult emergence. All these entomological variables seem to be strongly dependent on the ambient temperature or water temperature of the breeding site.
Breeding places availability model
The availability of larval habitats (Figure 1) has been based on level fluctuations of the quantity of water in breeding places of varying capacities. Table functions affecting the egg laying (oviposition) of a vector population -into the main module- have been estimated through Fuzzy Logic. This methodology was used to generalize the discrete events 'larval habitat completely available-completely full' and 'larval habitat not available-completely dry' to a continuous (fuzzy) form or partial truth. We considered the collections of water as 'reservoirs' ranging from 10 to 60 mm in depth opened to observed rainfall, estimated actual evaporation, and assumed constant infiltration. Daily evapotranspiration was estimated using the models introduced by: Turc (1945), assuming a homogeneous distribution of evapotranspiration throughout the year; Coutagne (1974) using Budyko's equation (1974); Cenicafé (1997), a specific evaporation model suited for Colombia; and Thornthwaite (1948), assuming a relation between actual and potential evapotranspiration of about 0.70 and a homogeneous distribution throughout each month [22].
As ovipositing anopheline females often prefer standing fresh water with low organic contents [23], this model determine the number of successive days that a breeding site remained between 90 and 100% of its capacity. This period of time is needed to complete all stages including egg batching, growth and metamorphosis in aquatic existence, and adult (imago) emergence. A 'favourable' condition is given when this period exceeds 12–15 successive days [2], and is normally affected by the water temperature of the breeding site. Finally, as anopheline eggs cannot survive prolonged desiccation, we assumed that a breeding site is 'not favourable' for the development of the species during pre-imago stages if it remains dry for more than two successive days. For the sake of simplicity, other relevant factors such as exposure to sunlight, emergent vegetation, temperature, salinity, and organic content [23], were not included in the model.
Predator-prey interaction models
This module allowed us to represent the interaction between the pre-imago stages of the mosquito population (mainly during larvae stage) and natural predators (Figure 1). As a first effort, the Lotka-Volterra predator-prey type model [24–26] has been introduced into the differential-equation system to analyse the dynamics of both the prey and predator populations. Parameters for other types of models that were included in the model, but are not currently used, including Henize predator-prey model [27, 9], Kaibab III prey-food model [10], Ecological System with External Forcing [9], and Autonomous Ecological System [9], are being calibrated through laboratory experiments.
System of coupled non-linear differential equations
The mathematical model represents the feedback mechanisms and interactions depicted in Figure 1. The comprehensive tool is based on a system of coupled non-linear differential equations, which is solved through a fourth order Runge-Kutta numerical algorithm. Simulations are being run using the computer software Powersim Constructor Version 2.51. The system of differential equations for the human population component and hence, the dynamics of malaria transmission in human hosts, is the following:
$\begin{array}{l}\frac{d\left(HUS\right)}{dt}={R}_{na}\cdot {P}_{hu}+\gamma \cdot HUM-{R}_{mo}\cdot HUS-VC\cdot \frac{HUF}{{P}_{hu}}\cdot HUS,\\ \frac{d\left(HUI\right)}{dt}=VC\cdot \frac{HUF}{{P}_{hu}}\cdot HUS-{\left[HUI\right]}_{t-\left({k}_{in}+{k}_{er}\right)}-{R}_{mo}\cdot HUI,\\ \frac{d\left(HUF\right)}{dt}={\left[HUI\right]}_{t-\left({k}_{in}+{k}_{er}\right)}-V\cdot HUF-{R}_{mo}\cdot HUF,\text{and}\\ \frac{d\left(HUM\right)}{dt}=V\cdot HUF-\gamma \cdot HUM-{R}_{mo}\cdot HUM.\end{array}$
The variables HUS, HUI, HUF, and HUM denote the total number of susceptible, malaria-infected, infectious, and immune human hosts at a given time t, respectively. In this component, Martens' equations [3] were followed, although considering the infective stage of the human hosts.
For the vector ecology component, the model uses the following system of coupled non-linear differential equations to describe the dynamics of the pre-imago stages of the vector population:
$\begin{array}{l}\frac{d\left(E\right)}{dt}{T}_{1}\left({f}_{L}\right)\cdot {f}_{L}\cdot {f}_{L\_S}\cdot {\u3008\left(VS+VI+VF\right)\cdot {R}_{po}\u3009}_{FI}-\frac{E}{{k}_{E}}-E\cdot {\mu}_{E}\cdot {T}_{2}\left(E\right),\\ \frac{d\left(L\right)}{dt}=\frac{E}{{k}_{E}}-\frac{L}{{k}_{L}}-L\cdot {\mu}_{L}\cdot {T}_{3}\left(L\right)-{a}_{S}\cdot L\cdot PD,\text{and}\\ \frac{d\left(PU\right)}{dt}=\frac{L}{{k}_{L}}-PU\cdot {\mu}_{PU}\cdot {T}_{4}\left(PU\right)-\frac{PU}{{k}_{em}}.\end{array}$
The dynamic variables E, L, and PU denote the total number of eggs, larvae, and pupae at a given time t, respectively. The table functions T_{1}, T_{2}, T_{3}, and T_{4} represent the multiplier factors affecting oviposition (competition during egg-laying), and eggs, larvae and pupae mortalities, respectively.
To describe the dynamics of the imago stages of the vector population and hence, the dynamics of malaria transmission in adult female mosquitoes, the model uses a second set of differential equations:
$\begin{array}{l}\frac{d\left(VS\right)}{dt}=\frac{PU}{{k}_{em}}-\left({\mu}_{m}+{\alpha}_{m}\right)\cdot VS-VS\cdot f\cdot \frac{HUF}{{P}_{hu}}\cdot S\_V,\\ \frac{d\left(VI\right)}{dt}=VS\cdot f\cdot \frac{HUF}{{P}_{hu}}\cdot S\_V-\left({\mu}_{m}+{\alpha}_{m}\right)\cdot VI-\frac{VI}{n},\text{and}\\ \frac{d\left(VF\right)}{dt}=\frac{VI}{n}-\left({\mu}_{m}+{\alpha}_{m}\right)\cdot VF.\end{array}$
The dynamic variables VS, VI, and VF represent the susceptible or non-infectious, malaria-infected, and malaria infectious mosquito (adult females) population sizes at a given time t, respectively.
To represent the dynamics of natural predators, the following equation is initially considered in the overall malaria model:
$\frac{d\left(PD\right)}{dt}={b}_{p}\cdot L\cdot PD-{m}_{p}\cdot PD.$
Exogenous variables considered for mosquito population
Exogenous variable | Used variable | Depending on | Default value | |
---|---|---|---|---|
Hatching time delay | k_{E} or R_E | [days] | Water T | (A) |
Eggs becoming non-viable | φ_{E} or E_NV | [days] | Water T | (B) |
Larvae developing time delay | k_{L} or R_L | [days] | Water T | (C) |
Larvae becoming non-viable | φ_{L} or L_NV | [days] | Water T | (D) |
Adult emergence time delay | k_{em} or R_EM | [days] | Water T | (E) |
Pupae becoming non-viable | φ_{PU}, PU_NV | [days] | Water T | (F) |
Induced mortality of mosquitoes | IM_M | [days] | SEC | 98–191.8 [21] |
Plasmodium species | P. falciparum | |||
Human Blood Index | HBI | [dec] | Species | 0.38–0.46 [3] |
Rate of oviposition (eggs per batch) | R_{po} or R_O | [eggs/vector] | T, SEC | 75–150 |
Degree days required for digestion of blood | D_{bd} | [°C-day] | RH | 36.5 at RH 70–80% [3] |
Minimum T required for digestion of blood meal | T_{min, bd} | [°] | 9.9 [3] | |
Degree days required for parasite development | D_{m} | [°-day] | Species | 100–120 [3] |
Min T required for parasite development | T_{min, p} | [°] | Species | 16 [3] |
Endogenous variables considered for mosquito population
Endogenous variable | Used variable | Depending on | Function | |
---|---|---|---|---|
Rate of eggs becoming non-viable | μ_{E} or Mu_E | [1/day] | φ_{E} | $\frac{1}{{\phi}_{E}}$ |
Rate of larvae becoming non-viable | μ_{L} or Mu_L | [1/day] | φ_{L} | $\frac{1}{{\phi}_{L}}$ |
Rate of pupae becoming non-viable | μ_{PU} or Mu_PU | [1/day] | φ_{PU} | $\frac{1}{{\phi}_{PU}}$ |
Daily survival probability | p or P_S | [dec] | T | ${e}^{-\frac{1}{\left(-4.4+1.31*T-0.03*{T}^{2}\right)}}$ [3] or longevity under CLC. |
Natural mortality rate of mosquitoes | μ_{m} or Mu_M | [1/day] | p | 1 - p |
Induced mortality rate of mosquitoes | α_{m} or Alfa_M | [1/day] | IM_M | $\frac{1}{IM\_M}$ |
Feeding interval | FI | [days] | T, D_{bd}, T_{min, bd} | $\frac{{D}_{bd}}{T-{T}_{min,bd}}$ [35,31] or length of the gonotrophic cycle under CLC. |
Sporogonic cycle | n or EIP | [days] | T, D_{m}, T_{min, p} | $\frac{{D}_{m}}{T-{T}_{min,p}}$ or length of the sporogonic period under CLC. |
Frequency with which human blood meals are taken | a | [1/day] | FI, HBI, Risk | $\frac{HBI}{FI}$ |
Exogenous variables considered for the human population
Exogenous variable | Used variable | Depending on | Default value | |
---|---|---|---|---|
Natural per-capita natality rate | R_{na} | [/day] | SEC | For simplicities' sake, R_{na} equals the natural death rate |
Natural per-capita death rate | R_{mo} | [/day] | SEC | 1/(52.5*365) for intermediate conditions |
Average infected period | ν | [years] | Parasite species | 0.9–1.5 for Plasmodium falciparum; 3 for Plasmodium vivax |
Mean duration of immunity | τ | [years] | 1.5 |
Endogenous variables considered for human population
Endogenous variable | Used variable | Depending on | Function | |
---|---|---|---|---|
Total human population at risk | P_{hu} | [individual] | HUS, HUI, HUF, HUM | HUS + HUI + HUF + HUM |
Primary exo-erythrocytic schizogony | k_{in} | [days] | Parasite species, SEC | 5–13 for P. falciparum; 2 for P. vivax |
Erythrocytic schizogony | k_{er} | [days] | Parasite species, SEC | 2 for P. falciparum; 2 for P. vivax |
Average infectious period | 1/V or P_F | [years] | HUF, ν, VC, P_{hu} | $\frac{EXP\left[VC*\frac{HUF}{{P}_{hu}}*\nu \right]-1}{VC*\frac{HUF}{{P}_{hu}}}$ [3] |
Average immune period | 1/γ or P_M | [years] | HUF, τ, VC, P_{hu} | $\frac{EXP\left[VC*\frac{HUF}{{P}_{hu}}*\tau \right]-1}{VC*\frac{HUF}{{P}_{hu}}}$ [3] |
Prevalence** | G or Prev | [dec] | HUF, P_{hu} | $\frac{HUF}{{P}_{hu}}$ |
Exogenous variables considered in the infectious disease models
Exogenous variable | Used variable | Depending on | Default value | |
---|---|---|---|---|
Human Biting Density | ma or HBD | [blood-meals/person/night] | Total mosquitoes captured (indoor and outdoor landing captures 06:00–09:00 pm) | |
Sporozoite Rate | s or SR | [dec] | 0.07–0.12 | |
Human susceptibility (efficiency with which an infective mosquito infects a susceptible human) | S_H | [dec] | 1.00 | |
Mosquito susceptibility (efficiency with which an infective human infects a susceptible mosquito) | S_V | [dec] | 1.00 | |
Proportion of infective female vectors (proportion of those anophelines with sporozoites in their salivary glands which are actually infective) | b | [dec] | 0.01 | |
Recovery rate in man | C1 | [dec] | 1.00 | |
Transmission rate | f | [1/days] | f (Risk) | (G) |
Endogenous variables considered in the infectious disease models
Endogenous variable | Used variable | Depending on | Function | |
---|---|---|---|---|
Critical density | DCR | [mosquitoes/human] | S_H, S_V, a, p, n, C1 | $C1*\frac{-LN\left(p\right)}{S\_H*S\_V*{a}^{2}*{p}^{n}}$ [3] |
Simulated density | DC | [mosquitoes/human] | VS, VI, VF, P_{hu} | $\frac{VS+VI+VF}{{P}_{hu}}$ |
Calibration parameter (VC) | K1 | HBD, DCR, DC, D_0, D_1 | ||
Entomological Inoculation Rate | EIR | [infected bites/human/night] | m, a, s, b | m*a*s*b |
Basic Reproduction Rate | Z | VC, 1/V | VC*(1/V) | |
Interaction Susceptibility | IN_F | [dec] | R_IN | Decay function: 1.2 if R_IN = 0; 0.853 if R_IN = 2; 0.812 if R_IN = 5; and 0.802 if R_IN = 15 |
Vectorial Capacity | VC | [infected bites/human/night] | IN_F, K1, a, p, n | $IN\_F*K1*\frac{{a}^{2}*{p}^{n}}{-LN\left(p\right)}$ [3] |
The function ${\u3008\left(VS+VI+VF\right)\cdot {R}_{po}\u3009}_{FI}$ represents the net oviposition that occur every FI pulse (see feeding interval in Table 2). That is, it is assumed that the interval between blood meals (i.e. the length of the gonotrophic cycle) is equivalent to the time interval between successive ovipositions. In the predator-prey interactions module, PD represents the predator population, a_{S} the predation rate coefficient, b_{p} the reproduction rate of predators per 1 prey eaten, and m_{p} the predator mortality rate. In the breeding places availability model, f_{L} and f_{L_S} represent, respectively, the water availability and the level of desiccation (table functions FL_P and FL_P_S) affecting the egg laying in vector population into main module. To consider the preferences of ovipositing females for specific larval habitats, weighted factors of 0.40, 0.30, 0.10, 0.10, 0.05, and 0.05 were assumed for the breeding sites of 10, 20, 30, 40, 50, and 60 mm of capacity, respectively. Finally, the carrying capacity of the mosquito population was defined assuming a maximum vector-human ratio of 30:1.
Study sites and demographic data
Climate data, missing periods and homogeneity analysis
Description of the climate data
Endemic region | Time period | Historical time series | Daily timescale | Monthly timescale | ||
---|---|---|---|---|---|---|
Total | Missing | Total | Missing | |||
Nuqui | Nov/01/1997 – Dec/31/2004 | Mean daily temperatures | 2,618 days | 459 | 86 months | 14 |
Total daily rainfall | 672 | 22 | ||||
Mean daily relative humidity values | 462 | 14 | ||||
El Bagre | Jan/01/1990 – Dec/31/2004 | Mean daily temperatures | 5,479 days | 1,234 | 180 months | 23 |
Total daily rainfall | 1,018 | 30 | ||||
Mean daily relative humidity values | 1,561 | 27 |
As mentioned above, homogeneity analyses were conducted to detect changes in the mean, variance, as well as significant trends in the climatic data sets. To detect the point of change in the mean and/or the variance in each single hydrological time series the Bayesian Analysis (BA) and the Abbe' Criterion for Homogeneity Test were used, both at a 0.05 significance level. The BA allowed determining the mean and the mode of the strongest point of change (NC point) and the mean of the total amount of change. The Abbe test allowed determining whether the hypothesis of homogeneity of each historical time series could be accepted or rejected.
To detect changes in the variance of each single hydrological time series (using the point of change estimated by the BA), the Simple F Test, the Simple F Test with corrections of dependence, the Modified F Test using both Chi-square and F distributions, the Ansari-Bradley Test, the Bartlett Test, and the Levene Test were conducted, all at a 0.05 significance level.
To detect changes in the mean of each single hydrological time series (using the point of change estimated by the BA), the Mann-Whitney/Wilcoxon Rank Sum Test, the Simple T Test assuming change and no change in the variance, the Modified T Test assuming change and no change in the variance, the Simple T Test with corrections of dependence assuming change and no change in the variance, and the Kruskal-Wallis Test were conducted, all at a 0.05 significance level.
To detect trends in the historical time series (free of seasonality) the T Test for the detection of linear trends, the Hotelling-Pabst Test, the Man-Kendall Test, and the Sen Test were conducted, all at a 0.05 significance level.
Finally, to detect trends in the seasonal time series the T Test for the detection of linear trends, the Hotelling-Pabst Test, the Seasonal Kendall Test, and the Seasonal Homogeneity of Trends Test were conducted, all at a 0.05 significance level.
All these hypotheses tests were run using the ASH (Analysis of Historical Time Series) Software developed by J.D. Salas and R.A. Smith at the Hydrology and Water Resources Program, Colorado State University, and recently modified at the Water Resources Graduate Program, National University of Colombia at Medellin.
The hypotheses tests for detecting non-homogeneities in time series recorded in the Nuqui Region showed that the historical records are homogeneous for the entire available period (results not shown). Both time series of mean monthly temperatures and mean monthly relative humidity values registered in the El Bagre region exhibit, on the other hand, strong non-stationarities in the mean and variance in May, 1992 (Figure 3B). As the shift showed not to be climate-induced, it was considered that only historical records can be assumed homogenous from July 1^{st}, 1992 onwards. Noteworthy, mean monthly temperatures observed in both regions showed strong increases in mean annual temperatures (approximately +1.0 – +1.5°C) during the onset of the El Niño warm event, which comprised the periods 1994–1995, 1997–1998 and 2002–2003 (Figure 3).
Finally, all daily missing records for both temperature and relative humidity were substituted by an estimated value based on observed monthly records. Missing daily precipitation records were reconstructed assuming no rain on each day.
Epidemiological data, statistical and correlation analyses
Even though historical time series seem to be homogeneous in the Nuqui region, the continuous changes in the responsibilities of the local authorities have lead, on the contrary, to several non-homogeneities in epidemiological time series of the El Bagre region (Figure 4B). In this case, strong changes in the mean were detected in the epidemiological records. Although the errors or omissions appear to be important in such historical time series, datasets were only checked for accuracy, logic and range of values.
In the Nuquí region, monthly epidemiological data comprises the continuous period of January, 1994 through June, 2005 (138 months or 4,199 days). The Antioquia Health Service and the National Institute of Health provided regional malaria data for the El Bagre region, for the period spanning from January 1^{st}, 1990 through December 31^{st}, 2004 (1st EP 1990-13th EP 2004; 180 months or 5,479 days).
These data sets were collected according to the type of infection (P. falciparum, P. vivax, and mixed malaria) and age groups (<1, 1–4, 5–14, 15–44, 45–59, and >60 years), and were used to develop different types of statistical analyses for understanding the linear and no-linear degree of correlation between climatic variables and malaria transmission indexes. Statistical and correlation analyses were conducted using both the Microsoft Excel Data Analysis Tool and the SPSS version 11.5. The central tendency and variability of all these data sets were examined using descriptive statistics. Finally, the relationship between climatic variables and malaria cases (or incidences) was explored using linear regression and cross correlation with various time lags.
Correlation analysis
Nuqui region | El Bagre region | ||||||||
---|---|---|---|---|---|---|---|---|---|
Lag (month) | IVM (T) | IVM (P) | IFM (T) | IFM (P) | Lag (EP) | IVEP (T) | IVEP (P) | IFEP (T) | IFEP (P) |
0 | 0.772 | -0.250 | 0.810 | -0.332 | 0 | 0.315 | -0.155 | 0.385 | -0.067 |
1 | 0.833 | -0.293 | 0.827 | -0.384 | 1 | 0.323 | -0.317 | 0.487 | -0.315 |
2 | 0.818 | -0.301 | 0.788 | -0.409 | 2 | 0.312 | -0.349 | 0.466 | -0.436 |
3 | 0.832 | -0.400 | 0.733 | -0.447 | 3 | 0.342 | -0.416 | 0.472 | -0.491 |
4 | 0.790 | -0.441 | 0.608 | -0.413 | 4 | 0.280 | -0.435 | 0.398 | -0.496 |
Finally, correlation analyses show that epidemiological time series seem to be strongly associated with mean temperatures, particularly during the onset of the El Niño warm event, as shown in Figure 4. Noteworthy, similar increments in the number of malaria cases during the El Niño event have been found throughout Colombia [5, 6, 33, 34].
Base scenarios
Roughly speaking, it is generally accepted that there are two scenarios during which malaria transmission could be possible in Colombian endemic-regions: (1) High levels of disease transmission could be expected during periods of 'excellent' conditions following long rainy seasons, when mosquito densities are generally high and temperatures are not low enough to inhibit a successful parasite development within the mosquito host. This scenario is primarily controlled by vector density and suggests that temperature and precipitation variables might have synergistic effects on malaria transmission. Temperature affects those entomological variables relevant to the biology of the mosquito host mainly during its aquatic stages; precipitation controls the availability (and, in turn productivity) of adequate breeding sites. This scenario requires linking and simulating the dynamics of vector ecology and malaria transmission in adult mosquitoes and human hosts. (2) High levels of disease transmission could be also expected during periods of 'good' conditions following warm dry seasons, periods when vector densities tend to be extraordinarily low. Although this scenario is 'good' in terms of temperature for the development of parasites in the mosquito host, it seems not to be particularly favourable for disease transmission in terms of low vector densities. Under this scenario, transmission is primarily controlled by those entomological variables (sporogonic cycle and feeding interval) that are strongly affected by temperature. High temperatures could lead to a shortening of the duration of both the extrinsic incubation period and the gonotrophic cycle. The simulation of the dynamics of malaria transmission could be then simplified by estimating a constant vector density (or a smooth seasonal fluctuation of vector density) and separately simulating the disease transmission in human hosts.
As the analysis of annual cycles of rainfall, temperature and malaria incidence patterns shows peaks in malaria transmission following the periods of high temperatures (and when rainy seasons starts to occur), as shown in Figures 5 and 6, it was proposed the second epidemiological scenario discussed above for the Nuqui and El Bagre regions. Hence, malaria transmission was simulated assuming the Vectorial Capacity infectious disease model for specific constant mosquito densities, which were estimated through the analysis of the availability of larval habitats in the selected localities.
Based on the continuous and homogeneous periods of climatic and epidemiological data sets available in the El Bagre region, the observed P. falciparum malaria incidence in the study site was modelled for the period from January 3^{rd}, 1994 to December 28^{th}, 2002 (first epidemiological period 1994 – 13th EP 2002). This period corresponds to a simulation exercise of about 3,282 days or 108 months, or 9 years (Figures 7B and 7C). The initial total human population at risk was assumed to be equal to 15,862 people, the population of rural areas at the end of 1994. The initial number of total infectious, infected and immune individuals was assumed to be equal to HUF(0) = 326, HUI(0) = 261, and HUM(0) = 453 people, equivalent to the total number of positive cases reported during the 13^{th} EP 1993, 1^{st} EP 1994, and 12^{th} EP 1993, respectively. Accordingly, the total number of susceptible individuals HUS(0) was assumed to be equal to 14,822 individuals. In the endemic area, malaria incidence was simulated assuming temperature values equal to either mean daily indoor temperatures or a constant value of 27.7°C, equivalent to the mean annual temperature of this malaria prone-region.
The system of coupled differential equations was solved by using a fourth order Runge-Kutta numerical algorithm for a 1-day time step at both study sites. Values of exogenous variables are shown in Tables 1, 3 and 5.
Preliminary analysis of instability cases
Sensitivity to initial conditions
Climate scenarios
Four changing climate scenarios were simulated for the El Bagre region assuming increases in mean daily temperatures for the entire simulation period from January 3^{rd}, 1994 to December 28^{th}, 2002. Based on the Special Report on Emission Scenarios and the predicted Colombian Climate Change Scenarios for year 2050 (particularly for the study site), increases of 1.0°C, 1.6°C, 1.9°C, and 2.8°C in mean daily ambient temperatures were assumed. Although it does not simulate a future scenario per se, it represents the possible dynamics of malaria transmission of an endemic area in which a population with similar characteristics of the El Bagre region lives. Hence, simulation of a changing climate scenario will only allow for the realization of a strong non-linearity involved in the mathematical model.
Results
Vector densities
During a 1,800-day period assumed for the El Bagre region (January 3rd 1994 – December 7th 1998) mean daily temperatures ranged from 24.0 to 32.5°C, and at least 83% of the observed values were in the interval (27-30°C]. Within these temperature ranges, the critical density necessary to maintain parasite transmission ranged from 16.9 to 6.6, and from 8.9 to 6.7 mosquitoes per human host, respectively. The assumed simulation period also showed five wet seasons, which are generally associated with peaks in vector densities, and four dry seasons, which commonly lead to low levels in this entomological variable. Simulation results of the dynamics of mosquito ecology showed high vector densities of about 15, 17, 12, 17, and 12 mosquitoes per human host during these wet seasons, and low values of about 3, 5, 4, and 3 vectors during dry periods in this study site.
Similar analyses were conducted for the Nuqui region although results are not discussed here. However, it is worth mentioning that during the simulation period of November 1^{st}, 1997 to December 31^{st}, 2003 (2,252 days), mean daily temperatures in this area ranged from 22.6 to 29.1°C, and the critical density necessary to maintain parasite transmission reached values ranging from 7.5 to 28.6 mosquitoes per human host.
The combined analyses of availability of breeding sites, predator-prey interactions, critical densities for malaria transmission, and dynamics of vector ecology, showed that constant vector densities, ranging from 7 to 12 mosquitoes per human host, might be considered for simulating the dynamics of malaria transmission and conducting the analysis of instability cases in the El Bagre region. Similar simulation exercises carried out for the Nuqui region showed constant vector densities ranging from 10 to 16 mosquitoes per individual.
Finally, simulation results also showed average vector densities of about 11, 10, and 9 mosquitoes per human host for the periods January 3^{rd}, 1994 to March 12^{th}, 1996 (time steps 1–800), March 13^{th}, 1996 to November 2^{nd}, 1997 (801–1,400), and November 3^{rd}, 1997 to December 7^{th}, 1998 (1,401–1,800), if a crude seasonality of vector density is desired for simulating malaria transmission in the El Bagre region.
Base scenarios
Under the base scenario, the Vectorial Capacity (VC) and the Basic Reproduction Rate (R_{o}) estimated for El Bagre showed values ranging from 0.6 to 1.8 potential infective bites per day, and from 0.5 to 1.8 arising secondary cases, respectively, for mean daily outdoor temperatures and the assumed vector density. For the interval 27–30°C, VC ranged from 1.0 up to 1.8 daily potential infective bites. In Nuqui, VC and R_{o} reached values in the ranges 0.5–1.9 and 0.5–2.1, respectively.
The estimated R_{o} for the El Bagre region showed a broad range of values although all of them remained at lower levels. For temperatures below 27°C, R_{o} showed values below the transmission threshold (R_{o} = 1) meaning that the disease cannot develop in this community. For temperatures above 27°C, most of R_{o} values ranged between 1.0 and 2.0, suggesting that the disease can develop in the population at risk but at low endemic levels. In Nuqui, R_{o} reached values below the transmission threshold for temperatures below 25°C. For temperatures above this threshold, R_{o} values ranged linearly between 1.0 and 2.0.
Instability cases
In El Bagre, simulation results for changes in vector density show low correlation coefficients (less than 0.351) and high mean square errors (above 0.0030) for densities greater than 11 mosquitoes per human host (Figure 8). For instance, if vector density is set to a constant value of 12 mosquitoes per host, the model significantly overestimates observed malaria incidence (R = 0.293 and MSE = 0.0055; see solid line 1). For values of this entomological variable lower than 7–8 (below the critical density threshold necessary to maintain parasite transmission), the model shows, as expected, that malaria disease recedes and dies out due to the low values of the Basic Reproduction Rate. If vector density is set to 9 mosquitoes per human host, the model represents (R = 0.662 and MSE = 0.0001, see solid line 3) the observed malaria incidence but mainly for the period from January 1994 to July 1999. If vector density is set to 10 mosquitoes per human host (base scenario, see solid line 2), the correlation coefficient drops to 0.455 and the mean square error reaches 0.0010. In this case, the model adequately represents the malaria incidence for the period from January 1994 to November 2000. Modelling results suggest that a density of about 9–10 mosquitoes per human host are required for the development of malaria in the community at risk of the El Bagre region.
Simulation results (not included) for changes in Human Blood Index (HBI) show correlation coefficients and mean square errors reaching 0.455–0.395 and 0.0010–0.0019, respectively, for values of HBI in the range 0.41–0.42. For values above 0.44, the model considerably overestimates observed malaria incidences. In this case, simulation results show low correlation coefficients (less than 0.318) and high mean square errors (greater than 0.0042). If HBI is set to 0.38, the model is not able to reproduce either the observed malaria incidence or the recorded outbreaks. Modelling results suggest that a Human Blood Index of about 0.40–0.41 is required for the development of malaria in the community at risk. If this exogenous variable is set to these values, the model adequately represents malaria incidence but mainly for the period from January 1994 to November 2000.
Simulation results (not included) for changes in the number of Degree Days for Parasite Development (DD) in the mosquito host show high correlation coefficients (0.455–0.529) and low mean square errors (0.0010–0.0005) for values of DD in the range 111–115°C-days. If DD is set to 111°C-days, the Extrinsic Incubation Period (EIP) or sporogonic cycle decreases from 13.9 days at 24°C to 7.0 days at 32°C. If this exogenous variable is set to 115°C-days, the EIP is shortened from 14.4 days at 24°C to 7.2 days at 32°C. If DD is set to 120°C-days (EIP ranging between 15 and 7.5 days) the correlation coefficient increases to 0.675 and the mean square error decreases to 0.0001. If DD is set to 105°C-days (EIP ranging between13.1 and 6.6 days), R drops to 0.383 and MSE increases to 0.0021. In this case, the mathematical model considerably overestimates observed malaria incidences.
Simulation results (not included) for changes in the Human Average Infected Period (ν) show correlation coefficients and mean square errors reaching 0.455–0.315 and 0.0010–0.0059, respectively, for values of ν in the range 0.95–1.10 years. If ν is set to 1.50 years, the correlation coefficient drops to 0.199 and the mean square error significantly increases to 0.0403. In this case, the model overestimates the recorded malaria incidence. If this exogenous variable is set to 0.80 years, the model shows that malaria disease recedes and dies out due to low values of R_{o}.
Finally, simulation results (not included) for changes in the time lag regarding the primary Exo-erythrocytic Schizogony (k_{in}) show high correlation coefficients and low mean square errors for values of k_{in} in the range 8–10 days. For values above or below this threshold, correlation coefficients and mean square errors do not change significantly.
Sensitivity to initial conditions
Under the assumed first scenario and for a daily survival probability (p) of 0.899, simulation results show that the mathematical model is able, although not incredibly so (R = 0.396 and MSE = 0.0049), to reproduce increasing and decreasing trends observed in malaria incidences. Nevertheless, the model considerably overestimates recorded disease outbreaks. If p is set to 0.890, when temperatures are in the interval (27–30°C] (as shown in Figure 9), R increases to 0.566 and MSE drops to 0.0009. In this case, the model adequately represents the observed malaria incidence for the period from January 1994 to January 2002 and overestimates malaria outbreaks during February-November 2002. For daily survival probabilities below 0.880, when temperatures are above 30°C, the model shows that the malaria disease dies out. For values of p ranging from 0.895 to 0.901, when temperatures are in the favourable range 29–24°C for vector survivorship, R ranges from 0.448 to 0.374, and MSE ranges from 0.0027 to 0.0062.
Under the assumed second scenario and for a feeding interval (FI) of 2.0 days, simulation results show that the mathematical model is not able (R = 0.196 and MSE = 0.0011) to represent seasonality, prevalence and time of occurrence of recorded outbreaks. The model briefly simulates the observed values for the period from January 1994 to June 1999 but significantly overestimates malaria incidence from July 1999 onward. If FI is set to 2.1 days, when temperatures are in the favourable range for vector survivorship (Figure 9), the model adequately represents the observed incidence (R increases to 0.472 and MSE decreases to 0.0001). For values of FI below 2.0 days, i.e. 1.7 days, R dramatically drops to 0.008 and MSE increases to 0.0115. As expected, malaria incidence and outbreaks are significantly overestimated.
Finally, simulation results for the assumed third scenario show that the mathematical model is not able (R = 0.151 and MSE = 0.0042) to represent seasonality, prevalence and time of occurrence of recorded outbreaks. For a constant sporogonic cycle (n) of 8.0 days, the model significantly overestimates malaria incidence throughout the entire simulation period. If n is set to 10 days, when temperatures are in the favourable range for vector survivorship (Figure 9), the model adequately represents the observed disease incidence (R increases to 0.576 and MSE decreases to 0.0001). Finally, for values of n below 8 days, i.e. n = 7 days, R decreases to 0.113 and MSE increases to 0.0082. As expected, malaria incidence and outbreaks are significantly overestimated.
Final simulation scenarios
Figure 12 also depicts the time series of model results for a first final simulation scenario considered for El Bagre region (see solid line 1 in Figure 12C), in which (a) a crude seasonality of mosquito density (values between 9–11 mosquitoes per human host) is assumed, according to the availability of adequate breeding sites; (b) the HBI exogenous variable is set to a constant value of 0.41; (c) the DD exogenous variable is set to 114°C-days; (d) the ν exogenous variable is set to 0.95 years; (e) the k_{in} endogenous variable is set to 8 days; (f) the daily survival probability of the mosquito host is set to a constant value of 0.890; and (g) the endogenous variables FI and n are assumed to be function of mean daily temperatures.
Under this final scenario, the proposed mathematical model adequately reproduced (R = 0.628 and MSE = 0.0005) both seasonality and prevalence of P. falciparum malaria recorded in the area during the period from January 3^{rd}, 1994 (first EP 1994) to December 28^{th}, 2002 (13th EP 2002). Noteworthy, is that the system of differential equations was able to replicate malaria outbreaks that were recorded in this endemic prone-region as a result of climatic and environmental anomalies associated with the occurrence of the El Niño in the Tropical Pacific (Figure 12C). The period from January 1994 to December 1996, which comprised the event of the El Niño 1994–1995, was better simulated. Nevertheless, modelling results still exhibit peaks in transmission two epidemiological periods before recorded outbreaks (see solid line 1 in Figure 12C).
Finally, Figure 12C depicts the time series of model results for a second final simulation scenario considered for the El Bagre region, in which the Vectorial Capacity was delayed one epidemiological period (see solid line 2 in Figure 12C). In this case, the correlation coefficient between observed and simulated malaria incidence rates increased to 0.668, whereas the mean square error remained constant.
Climate scenarios
Discussion
Application of the model during the validation periods from November 1^{st}, 1997 to December 31^{st}, 2003 and from January 3^{rd}, 1994 to December 28^{th}, 2002 for the Nuqui and El Bagre regions, produced adequate results, with correlation coefficients and mean square errors between observed and modelled malaria incidences reaching R = 0.897–0.668 (P > 0.95) and MSE = 0.0002–0.0005, respectively. Under the described epidemiological scenarios, modelling results suggest that malaria can develop in the populations at risk at low endemic levels and minor efforts can lead to the eradication of the disease.
For the El Bagre region it was found that a constant density of about 9–10 mosquitoes per human host are required for the development of malaria in the community. The results of the analysis of instability cases seem to suggest that, in order to counteract malaria in the selected region, it is not necessary to eliminate Anopheles mosquitoes from there entirely: field campaigns aimed at the control and eradication of the disease must reduce the number of vectors to values below 9 mosquitoes per human host. The model also suggests that if environmental conditions favour the development of the mosquito population (letting vector densities reach values above 11 mosquitoes per human host), malaria disease will propagate in the community and reach endemic levels of 9% of the total population at risk. In this case, malaria outbreaks could affect almost 15% of the total number of people living in the area.
The results of the analysis of changes in Human Blood Index seem to suggest that An. darlingi, the common primary vector found in El Bagre, is basically zoophilic in its feeding behaviour. It is of interest to note that, if this characteristic of the vector population drastically increases to 60% of total bites on humans, malaria disease could propagate in the community at high endemic levels and reach 20% of the total population at risk.
Based on modelling results it can therefore be argued that if an infection of P. falciparum lasts for a period of time greater than 1.1–1.2 years in humans (when left untreated), malaria could propagate rapidly in the population at risk and reach high-endemic levels. If the average infected period is shortened to periods below 0.9 years, the model suggests that the disease could be controlled through the ecosystems themselves and minor efforts can lead to its eradication. The Human Average Infected Period might be considered one of the critical variables for understanding the malaria cycle in the vertebrate host. Under our conditions in the field (research is still ongoing), most patients are treated immediately or within one month after the detection of malaria parasites in their blood streams.
Simulation results also suggest, as expected, that changes in the time lag representing the Primary Exo-erythrocytic Schizogony do not appreciably affect malaria transmission. In fact, the duration of the parasite cycle in human hosts is more dependent on host characteristics than on external ecological factors.
So far, the mathematical model suggests those entomological variables more strongly affected by climatic conditions and that in turn finally impact transmission potential of mosquito population. These parameters include the sporogonic cycle of the malaria parasite, the daily rate of vector biting, the daily rate of vector natality and the daily survival probability of the mosquito host. Field research must be undertaken in order to estimate how changes in microclimatic and environmental conditions of breeding and resting sites affect these relevant entomological variables. Other variables such as hatching, larval developing, and adult emerging time delays, which seem also to be affected by ambient temperatures, must be included as well.
However, simulation results show that the degree days required for digestion of a portion of ingested blood and for parasite development, which in turn influence the feeding interval and the sporogonic cycle inside vectors, respectively, might be the key variables that cause sensitivity to initial conditions and affect or control malaria transmission. According to these results, vector survivorship is NOT the most important element in the Basic Reproduction Rate of malaria. This conclusion is contrary to Macdonald's statement: 'the formula for the Basic Reproduction Rate holds that the influence of vector survivorship is greater than the influence of the average number of men bitten by one mosquito in one day or the sporogonic cycle, which are in turn greater than the influence of the proportion of anophelines with sporozoites in their salivary glands which are actually infective, or the proportion of affected people, who have received one infective inoculum only, who revert to the unaffected state in one day' [19]. Nevertheless, it has to be said that research is still ongoing.
According to the model and based on the simulated scenarios, malaria outbreaks in the selected regions are possible during the favourable periods following the onset of the El Niño warm event. This epidemiological scenario is primarily controlled by the sporogonic cycle and the feeding interval, which are strongly affected by ambient temperatures. Hence, under the conditions in the field, it appears that temperature becomes the most relevant parameter driving the final malaria incidence and must be considered an essential variable for an adequate representation of malaria transmission in the selected endemic areas. It has to be argued that precipitation, which controls the availability of adequate breeding sites, is responsible for representing observed fluctuations in vector density. Even though mosquito's density is not the key variable for understanding the patterns of malaria outbreaks in the selected region, the model results show that a crude seasonality of vector density is also needed for understanding disease transmission. The model suggests that between periods of high temperatures, malaria disease is maintained at low endemic levels due to the role of vector density. Field research is currently under way in order to consider possible synergistic effects of temperature and precipitation on malaria transmission.
Simulation results suggest that assessing the impact of climate on malaria transmission not only requires consideration of the changes in annual mean temperatures, but also, and more importantly, the extent of inter-annual variability in both temperature and rainfall.
Most of the mathematical models used to derive the standard formulae in malaria epidemiology assume a homogeneous mosquito population, random and uniform vector bites on humans, and a constant adult mosquito population size. Although the first two useful approximations were followed, simple differential equations and basic entomological exogenous variables were included for representing seasonal fluctuations of mosquito density. Thus, the intrinsic differences in mortality rates, the mosquito biting preferences, and the proximity to larval habitats were deliberately ignored. But the dynamics of vectors during pre-imago stages, the predator-prey-food interactions during larval stage, and the relationships between environmental factors and survivorship-behaviour of vectors were purposely included. In fact, the analysis of human biting rates available from field campaigns conducted up to this point shows that, as expected, population densities might be significantly associated with observed rainfall patterns. These vector dynamics, following the preliminary analyses, are being clearly captured by the mathematical model.
New analysis will focus on more detailed examination of key variables that control the dynamic response of malaria transmission in the study sites. Particular interests include testing the model for parsimony and simplicity. A thorough dissection of all the components is currently being conducted in order to capture the most basic interactions between climatic and entomological variables. The mathematical tool will be improved by the inclusion of new subroutines for predator-prey interactions, as well as new field and laboratory data for the entomological parameters.
New analysis will focus also on additional stage variables (following Yang's equations [20, 21]) that will be used to describe different levels of protection or anti-malarial immunity, which are commonly observed among people exposed to continuous and intense disease transmission. Migratory patterns and exogenous introduction of infected individuals at certain times are to be included in the model.
Conclusion
The observed monthly values of malaria incidence are characterized by seasonal oscillations alternating irregularly with high incidence periods. Modelling results must show at least such fluctuations of disease transmission. Simulations using only mean annual temperatures do not represent this characteristic. In contrast, simulations using mean monthly/daily temperatures adequately represent the seasonality of malaria prevalence. Apparently, results seem to suggest that the intra-annual cycle of mean temperatures is good enough to represent the seasonal oscillations of monthly malaria incidence rates in the selected study sites. However, simulation results show that seasonality of vector density also becomes an important factor towards understanding disease transmission, thus suggesting that seasonal and inter-annual variability of both temperature and rainfall are required to assess the impact of climate on malaria dynamics.
In the applied model described here, attempts were made to include the seasonality of vector density. The module representing the availability of breeding sites was intentionally used to estimate the seasonal fluctuations of such an important entomological variable. Modelling the pre-imago cohorts for different larval habitats resulted in an estimation of the availability of reservoirs. The BPAM model simulations suggested that (a) oviposition and the development of mosquitoes during pre-imago stages are significantly slower during dry periods, and (b) egg-laying by the mosquito population is appreciably enhanced during periods of notable availability of adequate larval habitats. The system of differential equations also included the simulation of vector dynamics in the field on a daily time scale. In addition to climate data, biological information was required to incorporate species-specific values for each of the parameters of the mathematical model. As a result, seasonal fluctuations in vector density were predicted, with values in the range of the estimated critical densities necessary to maintain parasite transmission.
Malaria is a complex disease that cannot be studied in a uni-disciplinary and exclusively qualitative-descriptive manner. The research project has used an interdisciplinary approach for understanding malaria epidemiology by studying, concurrently, the mosquito vector ecology, the malaria parasite, the human population dynamics, and the characteristics of our heterogeneous environment. Thus, this effort summarizes the joint work of diverse disciplines including Mathematics, Hydrology, Climatology, Entomology, Epidemiology, Field Study and Public Health Medicine.
The mathematical model and the proposed simulation process constitute promising tools to deepen the understanding of the entomological, epidemiological, and climatic interactions related to malaria transmission conducive to disease outbreaks. Mathematical models could help in comprehending how climatic (and non-climatic) factors affect the dynamics of malaria transmission in order to evaluate its spatial and temporal risks in each endemic area. A profound understanding of all these linkages is required for basic research, and for important practical mitigation and human health control interventions. In this context, biological models are significantly necessary to have a quantitative understanding of malaria transmission. As this model is still in its early stages of development, so far it has not been used for operational purposes. However, in the foreseeable future, biological/eco-epidemiological models could be implemented as powerful tools to diagnose possible dynamic patterns of malaria incidence under several entomological, demographic and climatic scenarios, as well as decision-making tools, for the early detection and control of disease outbreaks.
Transmission models could permit (a) to estimate the time of occurrence of unexpected malaria outbreaks; (b) to evaluate the possible magnitude of the concomitant sharp rises in the incidence of the disease; and (c) to pose and answer "what if" questions for remote areas potentially vulnerable to changing scenarios, for which data is commonly sparse. Mathematical models could also identify those elements or appropriate components of the complex eco-epidemiological systems for which there exist proven and effective interventions, particularly in the areas of mosquito control, education in environmental health, and treatment. That is, comprehensive models could provide quantitative goals for effective interventions adapted to the specific ecological circumstances of each endemic area. Then, mathematical models will be able to determine the most appropriate preventative actions that have to be taken in order to reduce the vulnerability of populations to climate-induced epidemics and to prevent outbreaks before they begin. Finally, they are also able to identify 'when' these interventions need to be implemented in each area.
Although malaria is a highly complex multi-factorial disease related to diverse socio-economic and demographic factors, it has been demonstrated here that environmental factors and climate variability go a considerable way in explaining fluctuations of disease incidence in these particular endemic areas. Forecasts of relevant climatic variables (future ENSO events) can be merged with selected mathematical tools (statistical and biological/eco-epidemiological models), creating a Malaria Early Warning System (MEWS) to facilitate early, coupled and environmentally sound public health interventions [8].
A complete MEWS framework will permit a continuous evaluation of the local risk of malaria transmission in the face of multiple changing scenarios. The MEWS could then be combined with several public health strategies in order to develop Integrated Malaria Surveillance and Control Systems. Other components should be included, such as: (a) epidemiological surveillance and control activities, (b) early diagnosis and treatment of primary cases, and (c) entomological surveillance and control activities.
Declarations
Acknowledgements
D. Ruiz thanks G. Yan from the Department of Biological Sciences (State University of New York at Buffalo, USA), C. Berberich from the Kumasi Centre for Collaborative Research in Tropical Medicine (Kumasi, Ghana), and S. Masaaki from the Research Centre for Tropical Infectious Diseases (Institute of Tropical Medicine, Nagasaki University, Japan) for their valuable comments. D. Ruiz also thanks Vivian Rojas Quintero for her patience and support. The research group thank 'Mineros de Antioquia' for their assistance in collecting and storing climate data, the Instituto de Hidrologia, Meteorologia y Estudios Ambientales-Colombia for providing climate data, Direccion Seccional de Salud de Antioquia and Instituto Nacional de Salud for providing epidemiological data and support.
The group also thank The Inter-American Institute for Global Change Research (Collaborative Research Network "Diagnostics and Prediction of Climate Variability and Human Health Impacts in the Tropical Americas") and COLCIENCIAS-Colombia for funding and supporting activities during the period 2000–2002. G. Poveda's research was partially supported by CIRES (University of Colorado, Boulder) through a Visiting Fellowship during 2001–2002. D. Ruiz was partially supported by the project 'PDF-B Integrated National Adaptation Pilot: High Mountain Ecosystems, Colombia's Caribbean Insular Areas, and Human Health/sub-activity Biological Modelling, activity Risk Assessment, sub-component Early Warning System, component Human Health', funded by The Global Environment Facility (GEF)/World Bank and Conservation International during Nov-Dec 2005. Current modelling activities are being supported by Programa Ingeniería Ambiental, Escuela de Ingeniería de Antioquia.
Authors’ Affiliations
References
- World Health Organization: World malaria situation in Part I. Wkly Epidemiol Rec. 1994, 72: 269-274.Google Scholar
- Gilles HM, Warrell DA: Bruce Chwatt's Essential Malariology. 1993, London: William Heinemann Medical Books LtdGoogle Scholar
- Martens WJM: Health impacts of climate change and ozone depletion. An eco-epidemiological modelling approach. 1997, The Netherlands: Maastricht University PressGoogle Scholar
- Epstein PR: Is global warming harmful to health?. Scientific American. 2000, AugustGoogle Scholar
- Poveda G, Graham NE, Epstein PR, Rojas W, Quinones ML, Velez ID, Martens WJM: Climate and ENSO Variability Associated with Vector-Borne Diseases in Colombia. El Niño and the Southern Oscillation, Multiscale Variability and Global and Regional Impacts. Edited by: Diaz HF and Markgraf V. 2000, Cambridge University Press, 183-204. 1Google Scholar
- Bouma M, Poveda G, Rojas W, Quinones ML, Cox J, Patz J: Predicting high-risk years for malaria in Colombia using parameters of El Niño-Southern Oscillation. Trop Med Int Health. 1997, 2: 1122-1127. 10.1046/j.1365-3156.1997.d01-210.x.View ArticlePubMedGoogle Scholar
- Poveda G, Rojas W: Evidencias de la asociación entre brotes epidémicos de malaria en Colombia y el fenómeno El Niño-Oscilación del Sur. Revista de la Academia Colombiana de Ciencia. 1997, 21: 421-429.Google Scholar
- Poveda G, Rojas W, Quinones ML, Velez ID, Mantilla RI, Ruiz D, Zuluaga JS, Rua GL: Coupling between annual and ENSO timescales in the malaria-climate association in Colombia. Environ Health Perspect. 2001, 109: 489-493.PubMed CentralPubMedGoogle Scholar
- Toro M, Aracil J: Qualitative analysis of system dynamics ecological models. System Dynamics Review. 1988, 4: 56-80.View ArticleGoogle Scholar
- Aracil J, Toro M: Métodos cualitativos en Dinámica de Sistemas. 1993, Sevilla: Secretariado de Publicaciones de la Universidad de SevillaGoogle Scholar
- Aron J, May RM: The population dynamics of malaria. The population dynamics of infectious diseases: theory and applications. Edited by: Anderson RM. 1982, London: Chapman and Hall, 139-179.View ArticleGoogle Scholar
- Blower SM, Small PM, Hopewell PC: Control strategies for tuberculosis epidemics: new models for old problems. Science. 1996, 273: 497-500.View ArticlePubMedGoogle Scholar
- Blower SM, Porco TC, Darby G: Predicting and preventing the emergence of antiviral drug resistance in HSV-2. Nat Med. 1998, 4: 673-678. 10.1038/nm0698-673.View ArticlePubMedGoogle Scholar
- Blower SM, Gerberding JL: Understanding, predicting and controlling the emergence of drug-resistant tuberculosis: a theoretical framework. J Mol Med. 1998, 76: 624-636. 10.1007/s001090050260.View ArticlePubMedGoogle Scholar
- Dietz K: Mathematical models for transmission and control of malaria. Malaria: principles and practice of malariology. Edited by: Wernsdorfer WH, McGregor I. 1988, Churchill Livingstone, New Cork, 2: 1091-1133.Google Scholar
- Dyner I: Dinámica de sistemas y simulación continua en el proceso de planificación. Primera edición. 1993, Medellín: Publicación Universidad Nacional de ColombiaGoogle Scholar
- Haile DG: Computer simulation of the effects of changes in weather patterns on vector-borne disease transmission. The potential effects of global climate change in the United States. Edited by: Smith JB, Tirpak DA. 1989, Document No. 230-05-89-057, Appendix G. Washington, D.C. U.S. Environmental Protection AgencyGoogle Scholar
- Mckenzie FE, Wong RC, Bossert WH: Discrete-event simulation models of Plasmodium falciparum malaria. Simulation. 1998, 71: 213-217.View ArticleGoogle Scholar
- Macdonald G: The epidemiology and control of malaria. 1957, London UK: Oxford University PressGoogle Scholar
- Yang H: Malaria transmission model for different levels of acquired immunity and temperature-dependent parameters (vector). Rev. Saúde Pública. 2000, 34: 223-231.View ArticlePubMedGoogle Scholar
- Yang H, Ferreira M: Assessing the effects of global warming and local social and economic conditions on the malaria transmission. Rev. Saúde Pública. 2000, 34: 214-22.View ArticlePubMedGoogle Scholar
- De Bruin HAR: Evapotranspiration in humid tropical regions. Hydrology of Humid Tropical Regions with Particular Reference to the Hydrological Effects of Agriculture and Forestry Practice. Proceedings of the Hamburg Symposium. 1983, IAHS Publ, 140:Google Scholar
- Gabinaud A: Ecologie de deux Aedes halophiles du Littoral Mediterraneen Français: Aedes caspius (Pallas, 1771), Aedes detritus (Haliday, 1833). Utilisation de la végétation comme indicateur biotique pour l'établissement d'une carte écologique. Aplication en dynamique des populations. PhD thesis. 1975, Universitte des Sciences et Techniques du Languedoc, Centre Universitaire de PerpignanGoogle Scholar
- Lotka A: Elements of physical biology. 1925, Baltimore: Williams&Wilkins CoGoogle Scholar
- Volterra V: Variazioni e fluttuazioni del numero d'individui in specie animali conviventi. Mem R Accad Naz dei Lincei Ser. 1926, 2: 31-113.Google Scholar
- Sharov A: Quantitative population ecology. 1996, On-line lectures, Department of Entomology, Virginia Tech, Blacksburg, VA, [http://www.ento.vt.edu/~sharov/PopEcol/popecol.html]Google Scholar
- Henize J: Course materials for principles of Dynamic Systems. 1971, MITGoogle Scholar
- Garrett-Jones C, Shidrawi GR: Malaria vectorial capacity of a population of Anopheles gambiae . Bull World Health Organ. 1969, 40: 531-545.PubMed CentralPubMedGoogle Scholar
- Service MW: A guide to medical entomology. 1980, London: The Macmillan Press LtdView ArticleGoogle Scholar
- Ruiz D, Poveda G, Quiñones ML, Velez ID, Rua G, Rojas W, Zuluaga JS: Modelación sistémica para el diagnóstico de la interacción clima-malaria en Colombia. Aplicación durante El Niño 1997–1998 y La Niña 1998–2000. Meteorologia Colombiana. 2002, 5: 41-48.Google Scholar
- Ruiz D, Poveda G, Quiñones ML, Vélez ID, Rúa G, Rojas W, Zuluaga JS: Modeling entomological-climatic interaction of malaria transmission at Nuquí (Colombian Pacific Coast). Case of study: El Niño 1997–1998 and La Niña 1998–2000. Conference Climate Variability and Change and their Health Effects in the Caribbean: Information for Climate Variability and Change Adaptation Planning in the Health Sector. 2002, PAHO/WHO. Bridgetown, BarbadosGoogle Scholar
- Ruiz D, Poveda G, Mantilla RI, Quiñones ML, Vélez ID, Rúa G, Rojas W, Zuluaga JS: Modelación de la interacción entomológica-climática de la transmisión de la malaria mediante Dinámica de Sistemas. Revista Colombiana de Entomología. 2003, 29: 191-201.Google Scholar
- Poveda G, Rojas W: Evidencias de la asociación entre epidemias de malaria en Colombia y el fenómeno El Niño/Oscilación del Sur. Revista Academia Colombiana de Ciencias. 1997, 21: 421-429.Google Scholar
- Poveda G, Rojas W, Quiñones ML: Relación entre brotes epidémicos de malaria en Colombia y el Niño-Oscilación del Sur. 1999, Informe COLCIENCIAS. Medellín, ColombiaGoogle Scholar
- Detinova TS, Beklemishev WN, Bertram DS: Age-grouping methods in diptera of medical importance. 1962, Geneva: WHO Monograph, 47:Google Scholar
- Boyd MF: Epidemiology: factors related to the definitive host. Malariology. Edited by: Boyd MF. 1949, Philadelphia and London: WB Saunders Company, 1: 608-697.Google Scholar
- Gillies MT, De Meillon B: The anophelinae of Africa south of the Sahara (Ethiopian zoogeographical region). 1968, South African Institute of Medical Research, Publication No. 54, 2Google Scholar
- White GB: Malaria vector ecology and genetics. Brit Med Bull. 1982, 38: 207-212.PubMedGoogle Scholar
- Molineaux L: The epidemiology of human malaria as an explanation of its distribution, including some implications for its control. Malaria: principles and practice of malariology. Edited by: Wernsdorfer WH, McGregor I. 1988, New York USA: Churchill Livingstone, 2: 913-998.Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.