 Research
 Open Access
 Published:
Spatial panorama of malaria prevalence in Africa under climate change and interventions scenarios
International Journal of Health Geographics volume 17, Article number: 2 (2018)
Abstract
Background
Malaria is highly sensitive to climatic variables and is strongly influenced by the presence of vectors in a region that further contribute to parasite development and sustained disease transmission. Mathematical analysis of malaria transmission through the use and application of the value of the basic reproduction number (R_{0}) threshold is an important and useful tool for the understanding of disease patterns.
Methods
Temperature dependence aspect of R_{0} obtained from dynamical mathematical network model was used to derive the spatial distribution maps for malaria transmission under different climatic and intervention scenarios. Model validation was conducted using MARA map and the Annual Plasmodium falciparum Entomological Inoculation Rates for Africa.
Results
The inclusion of the coupling between patches in dynamical model seems to have no effects on the estimate of the optimal temperature (about 25 °C) for malaria transmission. In patches environment, we were able to establish a threshold value (about α = 5) representing the ratio between the migration rates from one patch to another that has no effect on the magnitude of R_{0}. Such findings allow us to limit the production of the spatial distribution map of R_{0} to a single patch model. Future projections using temperature changes indicated a shift in malaria transmission areas towards the southern and northern areas of Africa and the application of the interventions scenario yielded a considerable reduction in transmission within malaria endemic areas of the continent.
Conclusions
The approach employed here is a sole study that defined the limits of contemporary malaria transmission, using R_{0} derived from a dynamical mathematical model. It has offered a unique prospect for measuring the impacts of interventions through simple manipulation of model parameters. Projections at scale provide options to visualize and query the results, when linked to the human population could potentially deliver adequate highlight on the number of individuals at risk of malaria infection across Africa. The findings provide a reasonable basis for understanding the fundamental effects of malaria control and could contribute towards disease elimination, which is considered as a challenge especially in the context of climate change.
Background
Malaria is one of the oldest and deadliest human vectorborne diseases. It is well known that malaria is transmitted among humans by the female Anopheles mosquito species, and the transmission cycle is essentially driven by the humanbiting habit of the mosquito [1]. Malaria vectors are found in tropical and subtropical areas of the world and in subSahara Africa, this disease is one of the major problems of public health [2, 3]. In 2015 an estimated 212 million (range 148–304 million) cases of malaria occurred worldwide and 429,000 people died, mostly children in Africa [2, 3]. In 2016, there were an estimated 216 million cases of malaria in 91 countries, an increase of nearly 5 million cases over 2015 [3]. The disease is endemic and/or epidemic, depending on the climatic parameters and ecological characteristics of the geographic regions [2, 3]. In general, occurrence of malaria vectors is organized into various patterns under the influence of temperature and the presence of breeding sites, which facilitates their reproduction and the increase of vector populations [2, 3]. These parameters therefore influence the epidemiology of malaria in a particular region.
The use of dynamic mathematical models to describe the time evolution of epidemiological systems dates back to centuries [4, 5]. The time evolution of the interactions between human population and malaria vectors are modeled using several methods, which include differential equations, discrete time maps, metapopulation networks and geospatial approaches [6, 7]. Differential equations where space is ignored and the total number of individuals in the population of vectors or hosts is constant are mostly applied. In reality, the dynamics of populations are usually effected by a variety of interactions (both inter and intra); and as a consequence, it is appropriate to consider a grouping of small local population of vectors occupying small habitat zone “network” distributed in uniform matrix [8, 9]. In certain recent models of spreading epidemic, the location of the patches in space is treated explicitly by taking into account the number of connections k (degree) that any given patch within the network may have [10]. A common feature in the majority of such studies is the inclusion of temperature as a climatic variable [11]. Temperaturedependent models have been abundantly developed and the effects of temperature on population dynamics of malaria vectors have been well analyzed [12,13,14]. The proliferation of such modeling approach was motivated by the fact that temperature is a key element in the development of malaria vectors and parasites. A simulation model that includes the four life stages of mosquito life cycle (egg, larva, pupa, and adult) using delayed differential equations is presented in Depinay et al. [12]. Through, this model, the authors showed that the ambient temperature conditions play important roles both on the mosquito lifehistory processes and the parasites. In the study carried out by Mordecai et al. [13] temperature was introduced on the basic reproduction number (R_{0}) and the results showed that the optimal temperature for malaria transmission is about 25 °C, whereas the work done by Parham and colleagues [14] estimated identical variable with a magnitude of 31 °C that is 6 °C higher than predicted in Mordecai et al. [13].
Mathematical analysis of the threshold values of R_{0} is an important and useful tool for the understanding of disease patterns. The basic reproduction number R_{0} depends on the following variables: mosquito density and biting rate, vector competence and survival rate as well as parasite development time within the mosquito (i.e. extrinsic incubation time) and human recovery rate, which are related to mosquito abundance, biology or physiology, and are linked to environmental conditions [13]. Mathematically, R_{0} is employed to characterize the possibility of a disease outbreak to occur if it exceeds 1 (R_{0} > 1), and the possibility of the disease to die out when R_{0} < 1. The quantity R_{0} can also help in determining the initial exponential increase in the number of infection during the disease outbreak [15]. In regions with endemic vectorborne diseases, it possible to determine potential control actions by varying the magnitude of the threshold value of R_{0} (usually 1). Such analysis provides appropriate guidance in the context of public health initiatives directed towards the reduction of the disease burden [14].
The control and eradication of infectious diseases are among the most important goals for improving public health. Although global eradication of communicable diseases (e.g. smallpox) was achieved [16], it has been difficult to eradicate malaria and other vectorborne diseases. Several countries in Africa, through diverse intervention measures [usage of artemisinincombination therapy (ACT), long lasting insecticidetreated nets (LLINs), indoor residual spraying (IRS), mass screening and treatment (MSAT), etc.] have strengthened malaria control programs with perceptible triumph in decreasing both the disease incidence and sustained transmission. However, climate variability and changes associated to other factors are potentially redefining the conditions and areas of suitability and competence of malaria vectors and the risk of the disease as well as the sustenance of residual transmission [17, 18]. The most direct way in which climate change is projected to effect malaria is associated with the variation of ambient temperature. This climate variable has an impact on both the mosquitoes’ lifehistory processes and the parasites development; which both contribute to malaria transmission [13].
Geospatial science is currently offering new opportunities to infer spatial and temporal knowledge linking models that incorporate important disease characteristic patterns to modeling and mapping [19]. The outcomes of such models often yield potential zones of disease distribution and vector occurrence [19]. In this context, different modeling techniques have been used in mapping the spatial distribution of malaria vectors and transmission areas within local, regional and global scales [20,21,22,23,24,25]. In [26] advanced statistical techniques (nonlinear discriminant analysis, random forest, and generalized linear model) were employed to investigate the environmental suitability in the Netherlands for three indigenous mosquito species. In [27] the annual entomological inoculation rate for the Plasmodium falciparum transmitted by Anopheles gambiae was utilized as a proxy to map the intensity of malaria transmission in Uganda. Climate parameters coupled with areas of seasonal abundance were used for spatial simulations of key malaria vectors (A. gambiae and Anopheles arabiensis) in Africa [28, 29]. Justifications in undertaking these studies were to reduce the disease burden through knowledge and information dissemination to empower public health workers and policymakers in better approaches to management of malaria. Still in this direction, the present study uses temperaturedependence aspect of the reproductive number obtained from dynamical mathematical network model to derive spatial distribution maps under different climatic and intervention scenarios to understand and forecast malaria transmission dynamics. The roles played by different interventions to control the disease and changes in temperature on the basic reproduction number are well analyzed.
Methods
Temperature data
The spatial simulation of R_{0} is made with temperature data obtained from WorldClim (http://www.worldclim.org/) and CCAFS (Climate Change, Agriculture and Food Security; http://www.ccafsclimate.org) databases. Data for the year 2000 were used to analyze the present situation whereas simulated data for the year 2050 were used to represent future scenario. The data are layers (grids) with a spatial resolution of 2.5 min containing average minimum and maximum temperatures.
Mathematical model
Among existing mathematical models developed for malaria, the model of Ross and Macdonald [30], which has laid the foundation of majority of current epidemiological models has attracted a lot of considerations. In conducting the spatial panorama of malaria prevalence in Africa under climate change and interventions scenarios, the Ross and Macdonald model [30] was selected. This model is entrenched on the assumption that at any given moment, an entire population of either humans or vectors can be divided into distinct compartments made of the susceptible (those who are vulnerable to infection) and infectious (those that have acquired infection and are able to infect others). The infection spreads by random contact between susceptible and infectious “compartments” of the human and mosquito populations [31]. The Ross–Macdonald model has two equations for describing the changes in the number of infectious hosts and vectors. The changes in the number of susceptible hosts and vectors are implicitly modeled since the host and vector population sizes are kept constant [32]. The patchy model of the system was proposed in [33] by considering vectorial transmission. Such transmission is considered to be strongly correlated to climate variables, which heavily influence the distribution of vectors across geographical scales [11]. However, the extent to which current and future projections of temperature contribute to the patterns of disease prevalence is yet to be well established. This understanding is of great importance in evaluating and better preparing for future malaria transmission/risk, and would inform strategies of addressing current residual malaria transmission in Africa. A multimalaria model intercomparison that considered the impact of climate change on the disease transmission at global scale is presented in [18]. Herein, we use a patchy processbased model that accounts for both the vector and parasite influences. From the model, the basis basic reproduction number R_{0} was derived and used to project the disease transmission at scale. Then the result was reproduced at country level to both analyze climate change and intervention impacts on the population. The latter factor, which is often neglected in most studies was considered here as the element that brings originality into the present research. In this framework, the Ross–MacDonald model in patch environment is described by the following equations [33].
The variables u_{ i } and v_{ i } are the numbers of infected host and vector respectively, in each patch. Without the loss of generality, the following variables were kept constant in all patches and are defined as: a is the biting rate of the vector, p is the probability that the bite of infectious vector will lead to the successful infection of a susceptible host, q is the probability that a susceptible vector that bites an infectious host will become infected, ξ is the rate at which infectious hosts recover, δ is the mortality rate of infected vectors and, M_{ Hi } represents the total host population in patch i, and M_{ Vi } is the vector population in patch I; m_{ ij } is the migration rate from patch i to patch j with i ≠ j.
Basic reproduction number in patches environment
The basic reproduction number generally represented by R_{0} is a fundamental quantity in the study and analysis of mathematical models in epidemiology. R_{0} is commonly defined as the number of secondary cases that can be produced by one case of an infected individual in a completely susceptible population. The general framework for the determination of the basic reproduction number can be found in [15]. R_{0} is presented for different disease transmission models based on a system of differential equations. The expression of \( R_{0} \) for the patchy Ross–Macdonald model presented above is given in [33]. To observe the effect of the coupling model components, we sequentially analyzed the model with single and double patches. The expression of the basic reproduction number for a single patch is given by [32]:
With double patches, the expression of the basic reproduction number is more complex and is given as follows [33]:
The m_{12} and m_{21} are the migration rates from patch 1 to patch 2 and from patch 2 to patch 1 respectively, \( M_{{V_{1} }} \) and \( M_{{V_{2} }} \) are the vector population densities in patch 1 and 2 respectively, while M_{ H } is the total host population per patch.
Two cases were considered in the analysis:

1.
If the population is distributed uniformly, then the migration parameters are set to m_{12} = m_{21} = k. In this instance, the identical proportion of vector population migrates from one patch to another, while keeping the overall populations at individual patches constant. Under this assumption the formulation of the basic reproduction number is given by:
$$ R_{0}^{2} = 2\hat{R}_{0}^{2} $$(4)$$ \begin{aligned} R_{0}^{2} & = a^{2} \frac{pq}{{\delta \xi M_{H} }}\frac{2k}{{2\left( {\xi + 2k} \right)}}\left[ {\left( {1 + \frac{\xi }{k}} \right)M_{{V_{1} }} + \left( {1 + \frac{\xi }{k}} \right)M_{{V_{2} }} } \right. \\ & \quad \left. { +\,\sqrt {\left( {\left( {1 + \frac{\xi }{k}} \right) (M_{{V_{1} }} + M_{{V_{2} }} )} \right)^{2} + 4M_{{V_{1} }} M_{{V_{2} }} } } \right] \\ \end{aligned} $$(5) 
2.
Here, the population is distributed uniformly, however, the migration parameters are different (m_{12} ≠ m_{21}) and the population of infected vectors in one of the patches is set to 0. This case depicts a situation whereby the population of infectious mosquitoes in one patch is maximum and null in the other patch. With these considerations, the expression of the basic reproduction number, which depends on coupling and other parameters of the system, is given by:
$$ R_{0}^{2} = a^{2} \frac{pq}{\delta \xi }*\frac{{M_{{V_{1} }} }}{{M_{H} }}\left( {1 + \frac{{\xi m_{21} }}{{m_{12} \left( {m_{12} + m_{21} + \xi } \right)}}} \right) $$(6)
The existence of an endemic equilibrium in the model is guaranteed when \( R_{0} > 1 \) [34]. The control of vectors contribute to the decrease of vectorhost ratio, which reduces the basic reproductive ratio of the pathogen; then decreases the equilibrium number of infectious hosts and vectors. Therefore, for successful disease control, the vector population density has to be reduced at least to a level below the entomological threshold that is mathematically expressed as:
Parameterization of the reproductive number
To parameterize the reproductive number, we opted to undertake a literature search. In Eq. (3), there are 10 variables, among which five are found to be dependent on temperature. They have been well described by authors in [34]. Therefore, under the assumptions that vectors and parasites are dependent on temperature, the function to represent the reproductive number as temperature dependent for a single patch is given in Eq. (8) below:
Similarly R_{0} for two patches is obtained by:
In this Eq. (9), other parameters ξ, m_{12}, m_{21}, and M_{ H } are kept constant. In [13], the optimal temperature for malaria transmission was estimated at around 25 °C. This estimate is close to the values obtained from a study on entomological inoculation rates [35]. Hence, for our parameterization we opted for the approach proposed by [13] in which a constant number of vector population density is assumed in each patch. The quantity M_{ Vi }(T) (i = 1,2) is estimated as:
EFD_{ i }(T) is the number of eggs laid by a vector per day, P_{ EAi }(T) is the probability that an egg survives to become an adult mosquito vector, and τ_{ EA }i(T) is the duration of egg to adult development in each patch i. The relation \( p\left( T \right)q(T) = \phi (T)e^{  \delta \left( T \right)EIP(T)} \) expresses the product of the probability that the bite of an infectious vector will lead to the successful infection of a susceptible host p, and the probability that a susceptible vector that bites an infectious host will become infected q with the vector competence ϕ(T). The mortality rate is δ(T) and the extrinsic inoculation rate is EIP. Using relation (9) and replacing the parameter time by rate, the complete temperature dependent R_{0} for malaria can be described as per the following cases:
When the system has only one patch Eq. (2) is given by:
When the system possesses two patches, then two subcases are generated and their equations are:

1.
\( R_{0}^{2} (T) = 2\hat{R}_{0}^{2} (T) \) with \( \hat{R}_{0}^{2} (T) \) defined by Eq. (10),

2.
\( R_{0}^{2} = a^{2} (T)\frac{{\phi (T)e^{  \delta (T)/PDR(T)} }}{{\delta (T)^{3} \xi }}\frac{{EFD\left( T \right)P_{EA} \left( T \right)MDR(T)}}{{M_{H} }}\left( {1 + \frac{{\xi m_{21} }}{{m_{12} \left( {m_{12} + m_{21} + \xi } \right)}}} \right) \)
Table 1 is the summary of the temperature dependent expressions of R_{0} parameters. They were obtained from the analysis described in Mordecai et al. [13].
Spatial projection of the reproductive number
Many countries in Africa experience a spatial variation of temperatures and the predictions of humidity and precipitation are highly uncertain, here temperature was considered as a key factor. The selection of temperature alone is further justified by the fact that a small increase in the value of this variable leads to a significant increase in malaria vector development time and in the frequency of blood feeding in adults [36]. It has been well established that the parasites responsible for malaria occurrence are transmitted to their hosts during blood meals by female vectors [36]. Therefore, to regionally map R_{0}, the monthly max and min temperatures are simultaneously extracted from the databases and their average estimated. The obtained information is organized into matrices with longitude as the column and latitude as the row. A point object, representing the temperature dependent mathematical expression of R_{0} is generated and applied in each geographical coordinate. A new matrix is formed with the values of the R_{0} in the respective geographical coordinates. The results are converted to American Standard Code for Information Interchange (ASCII) files and transferred to an open source software QGIS [37] for visualization of the potential risk areas of malaria transmission under selected climate scenarios. A similar approach was used for the intervention scenarios. Overall, the method exploits the spatial variation of temperature to provide geographical predictions.
Climate change and intervention scenarios
The choice of scenarios was guided by the availability of datasets, which were obtained from WorldClim (http://www.worldclim.org/) and CCAFS (Climate Change, Agriculture and Food Security) (http://www.ccafsclimate.org) databases. Downscale datasets for the year 2050 (future climate) of the SRESA1B [38] was chosen for the analysis.
To take into account the contribution of different interventions on the burden of malaria in Africa for our projections, we identified studies that recorded trends in the disease indicator over a period of time [39]. The study demonstrated that vector control has immensely contributed to the decline of the burden caused by malaria in subSaharan Africa. With the selected modeling framework, four parameters (mosquito biting rate a, vector competence bc, adult mosquito mortality rate µ, and the probability that mosquito eggs survive to become adult P_{ EA }), which are directly linked to malaria vectors population were adjusted based on intervention effects. The application of the interventions criteria on these parameters has a direct impact on the amplitude of the basic reproduction number R_{0}. The following two intervention scenarios were adopted: (1) uniform application of 40% reduction in the amplitude of a, bc,µ and P_{ EA } respectively; and (2) selection of two African countries and reducing the magnitude of a, bc,µ and P_{ EA } by up to 40% (Cameroon) and 80% (Kenya) [40, 41]. The baseline for the beginning of interventions was set to year 2000. Year 2050 represents the future scenario to measure the impacts of the interventions. In the analysis, the temperature values of the years 2000 and 2050 were successively replaced by the variable T in the mathematical expression of the reproductive number R_{0}.
Model validation
The process of determining the degree to which the simulated regional maps are accurate was conducted through comparison of observed risk areas for malaria transmission as described in [17, 42] with the maps of the basics reproduction rate R_{0} obtained from the present analysis. The P. falciparum entomological inoculation rate (PfEIR) is a measure of the proportion of exposure to infectious mosquitoes [27]. It is often translated as the amount of P. falciparum infective bites received by each human within a season or annually [27]. The synthesis of Annual P. falciparum Entomological Inoculation Rates (APfEIR) data found in the literature of 22 spatially distinct records in Africa, post1980, across the continent was used to compare the values of R_{0} obtained from our model [17, 42].
Results
In the present context, the basic reproduction number R_{0} is used as a tool in evaluating and predicting risk zones for malaria transmission. High values of R_{0} (> 1) correspond to malaria risk zones. Our results indicate that for a population density of M_{ H } = 100 and an immunity coefficient of ξ = 0.1, the values of R_{0} range between 0 and 9. Such values provide an overall view of the intensity of malaria transmission in the region. The areas in Africa where R_{0} was close to the value 9 correspond to locations with high malaria prevalence. The value of R_{0} tends to zero in areas with low risk of malaria transmission. The maximum values of R_{0} are obtained at the locations where the optimal temperature for malaria transmission is about 25 °C.
In the case of one patch, the graph of R_{0} is represented in Fig. 1a. For two patches, if m_{12} = m_{21}, the value of the basic reproduction number doubles. When m_{12} ≠ m_{21}, we set m_{12} = α m_{21} and substituted it in Eq. (3) and thenplotted R_{0} in three dimensions as presented in Fig. 1b. The result of the simulation was quite interesting as highlighted herein. In Fig. 1b, it was observed that as the ratio of the migration between the two patches increases, the basic reproduction number decreases up to a threshold value α = 5. Furthermore, an increase of the value of α does not change the value of R_{0}. In this circumstance, the numerical value of R_{0} obtained when α > 5 (i.e. where the basic reproduction number converge) is identical to the value estimated by Ross–Macdonald model with a single patch. Whenα = 1, the basic reproduction number is the double of the critical value of R_{0} obtained at α > 5.
It was observed that the introduction of the coupling between the patches only has effects on the amplitude of R_{0} and not the optimal temperature for the transmission of malaria. In order words, a model with one or two patches will produce almost similar outcomes when projected spatially. Based on these results, we decided to limit the production of the spatial distribution map of the basic reproduction number R_{0} to a single patch model given by Eq. (11). The distribution maps predicting the values of R_{0} in Africa are shown in Figs. 2 and 3. The predictions comprehensively match the known distribution map for malaria transmission and endemic zones in Africa.
The map representing the present distribution (Fig. 2a) inferred by the degree of magnitude of R_{0} indicates greater values in most malaria endemic regions such as west and central Africa. The magnitude of R_{0} considerably reduced as we move northward and southward and became almost close to zero in the temperature regions of northern and southern parts of Africa, which are currently characterized by a very low level of malaria transmission. A similar trend is observed in the highlands of eastern Africa, South Africa, western Cameroon, central Angola and the plateau of Madagascar.
For future projection (2050), we observe a shift in areas of transmission of the disease towards the northern and southern regions of Africa (Fig. 2b). A projection of the developed model indicates a high risk of transmission and establishment of the disease in these regions. Malaria transmission will remain present in regions close to the equator. The future scenario also indicates that the level of malaria transmission in the Sahara Africa will highly reduce.
The changes in the values of R_{0} due to the temperature change between the years 2010 and 2050 are presented in Fig. 2c. The map displays the difference ΔR_{0} = R_{01} − R_{02}, where R_{01}, is the basic reproduction of the baseline scenario obtained by replacing the variable temperature with the values of the year 2000 and, R_{02} the malaria transmission map for future scenario corresponding to substituting temperature with the values of the year 2050. Two classes of singularities are observed on the map: (1) areas where there will be an increase in the value of R_{0}, and (2) areas in which R_{0} values decrease considerably. The increases in the values are found in regions mores close to the equator and at the Horn of Africa. The southern and northern regions of Africa are also of great concern given the observed increase in value of R_{0} (Fig. 2c), which demonstrates a possibility of disease endemicity in these areas.
The effects of interventions were computed from the subtraction of the baseline (year 2000) map to future scenario (2050) map, in which changes in the amplitude of certain models parameters as described above were made. By including intervention measures into the model, most of the regions are likely to experience a decrease in the values of R_{0} (Fig. 2d). The central and eastern African highlands and the southern regions of Africa are with the least decline in R_{0} values. Current steady areas where malaria is endemic are characterized by values of R_{0} that tend to zero; justifying the efficiency of the interventions in reducing the disease burden.
The maps of two selected countries with different intervention levels are presented in Fig. 3. This analysis was carried out in order to observe the effect of the intervention at the local scale. For Cameroon, a large variation in the value of R_{0} is observed in western highlands and Adamawa Region (Fig. 3a, b). The magnitude of the difference ΔR_{0} varies from − 4.5 to 4. We defined, ΔR_{0I} = R_{01} − R_{02I} as the difference between the values of the reproduction number obtained for the baseline (the year 2000) and future scenario (the year 2050) with and without interventions; respectively. The application of the intervention measures considerably decreased the magnitude of ΔR_{0I} as many points on the map have negative values. These values range between − 6.4736 and − 0.0029. The smaller the value of ΔR_{0I}, the higher the impact of the intervention(s).
The maps of Kenya (Fig. 3c, d) indicate the values of ΔR_{0} ranging between − 4.00 and 4.00. By applying the intervention measures, the changes in the areas of prevalence between 2000 and 2050 (ΔR_{0I}) ranges from − 8.4188 to 0.0052. The lowest magnitude in the change of the values of the basic reproduction number was observed in highland areas of the countries. Overall, the changes in the values of R_{0} are more pronounced in Kenya than in Cameroon as Kenya may have benefited from high level of interventions than Cameroon.
To further confirm the performance of our approach, we overlaid the entomological APfEIR survey data with the magnitudes of the basic reproduction number. By considering a region with the APfEIR value greater than 0 as endemic, and if APfEIR < 1, the region is moderately endemic. APfEIR > 10 were predicted for areas with very high risk of the disease transmission. In Table 2 the summary of the basic reproduction numbers computed for sites specific to the values of APfEIR are presented. The values of the basic reproduction number are large enough in all sites corresponding to high values of APfEIR. Overall, the maps of R_{0} for malaria endemicity revealed minimal bias with predicted APfEIR map. Nevertheless, the present map with R_{0} has the general tendency to overestimate prevalence by R_{0}. The matching between predicted map of R_{0} and the estimates of APfEIR map indicate consistent agreement in the 2 parameters.
Discussion and conclusion
Climate variables play an important role in the dynamics, distribution, and transmission of vectorborne diseases such as malaria. Although rainfall is critical in providing suitable habitats for mosquitoes to breed, it explains little additional variance in malaria transmission; and therefore temperature is a key driver that affects the essential processes of mosquito biology and parasite life cycle [13]. To be more specific, temperature determines transmission intensity, including mosquito development rate, biting rate, and survival of the parasite within the mosquito [11, 43]. Because of the compatibility and relevance of including temperature data into dynamical mathematical models, the study here only accounted for temperature as the main climatic variable.
Mathematical models are generally an abstraction of reality. They have the ability to address the problem in form of equations that capture important linkages of complex transmission dynamics of the disease, which cannot be unraveled using laboratory and ecological experiments [44]. The basic reproduction number R_{0} determines the threshold values for which, disease models exhibit changes in their stabilities [45]. The analyses presented in [45] highlighted that, R_{0} can provide a reasonable estimate of the reduction level of malaria transmission intensity and therefore it offers a good measure to use as proxy for understanding and analyzing possible options to eliminate the disease [45]. However, the transition from malariafree to an endemic area is generally associated to different types of bifurcations. Bifurcation arises when a slight increase/decrease in the magnitude of model parameter values triggers a sudden topological change in the disease trend. Traditional mathematical models often consider one endemic equilibrium when R_{0} > 1, which is translated to a stable diseasefree equilibrium for R_{0} < 1 and unstable when R_{0} > 1 [46]. In such context, the bifurcation leading from a diseasefree equilibrium to an endemic equilibrium is considered forward. Backward bifurcation occurs when the condition R_{0} < 1 cannot be used alone to explain all the required disease elimination efforts. In such case, a backward bifurcation point should be identified and the thresholds of the values of R_{0} defined for the effective control of the disease. A backward bifurcation point coexists with the stable diseasefree equilibrium [46,47,48]. Such phenomenon has important public health implications because the condition of having R_{0} below unity will no longer be sufficient to guarantee the disease elimination; hence a new threshold for the basic reproduction number should be estimated and considered as the critical point, in which the elimination of the disease can be possible [47, 48]. The hypothesis of making m_{12} = α m_{21} helped to estimate the magnitude of the threshold (α) whose slight increase in the magnitude of its estimated value produces a decrease of R_{0} and hence create the possibility to obtain a globally stable endemic area of malaria. In contrast, the disease persists endemically in all cases when R_{0} > 1. In the process of using R_{0} to estimate the spatial representation of malaria transmission areas, an increase in the difference of R_{0}(ΔR) does not necessarily indicate that the areas in the map are an endemic zone of the disease. In fact, changes of the state of malaria transmission could be noticed for values of R_{0} less than unity.
Temperature was used to parameterize the expression of R_{0} in patches environment for understanding the dynamics of malaria transmission and later generate maps for the disease risk areas in Africa. Biting rate, vector competence, adult mortality rate, parasite development rate, eggs laid per adult female per day, eggtoadult survival probability, and mosquito development rate were expressed as dependent on temperature, hence the parameters help to include a climatic variable into the expression of R_{0} [13]. Combining these thermal responses constrained from high to low limits into a spatial distribution of temperature values yielded maps showing areas of potential transmission of malaria. The obtained maps of R_{0} agree to a certain level with existing malaria risk maps of Africa obtained from an experimental study such as the MARA maps [42]. The present results further agree with the spatial distribution map of P. falciparum malaria endemic areas in Africa [49]. The model predicted numerous regions with R_{0} greater than unity because of the clumping effects of spatial models [50]. The accuracy of our model is enshrined in the capacity of reproducing almost all known regions at risk of malaria, which were the same as that identified by various other models [42]. The estimates of the basic reproduction number magnitudes for randomly selected points in Africa compared to APfEIR further confirmed the accuracy of the approach used. But the variations between APfEIR and the values of R_{0} in endemic areas were substantial due to the shortrange heterogeneity used in the computation of R_{0} compared to the patchy distribution of the field datasets. Usually, APfEIR are computed with data obtained from disparate health units within the region to characterize the point location, whereas the values of R_{0} is estimated based on the exact values of the temperature at point location represented by its geographical coordinates.
In literature, several studies have used different approaches to tackle problems similar to the current study. For instance, Craig et al. [42] have proposed a “fuzzy logic” model for the distribution of stable malaria transmission in subSaharan Africa. However, the model development required detailed information, not always available in all localities of Africa. The study in [49] interpolated the probabilities of P. falciparum entomological inoculation rate (PfEIR), P. falciparum basic reproductive number (PfR), and P. falciparum parasite rate (PfPR) to generate the map for each quantity, which were then combined to yield the malaria transmission map. Although the outcome was realistic, the method did not include the dynamics of malaria vectors and parasites. Other studies [28, 29] used CLIMEX model, a platform designed to infer species responses to environmental conditions and climatic parameters to produce potential distributional maps of disease occurrence. The authors only focused on A. gambiae and A. arabiensis distributions to infer malaria transmission, with no consideration of the parasites. Overall, malaria parasite transmission intensity is spatially heterogeneous and this heterogeneity has important implications for risks and age patterns of progression from malaria infection to disease, infirmity, morbidity, and death [51]. The present study is among the very first investigations that define the limits of contemporary malaria transmission, using the basic reproduction number derived from a dynamical mathematical model. The authors in [14] used a similar approach but limited the study to expected changes in transmission at the Republic of Tanzania. In addition, the models used by Parham and Michael [14] failed to predict accurately the optimal value of temperature for malaria transmission, thus producing maps with false high magnitude areas of malaria transmission. Studies using approaches similar to what is presented here were applied to map bluetongue virus in the Netherlands [52]; to develop a temperaturedriven map of the emerging tick vector of Lyme disease Ixodes scapularis in Canada [53], and to assess the effect of climate change in the risk of Chagas disease transmission in Colombia [54].
The link between climate change and vectorborne diseases such as malaria is well established [55]. Although other studies highlighted the importance of climate change on the burden of malaria in Africa [56, 57], few researches have considered temperature as a paramount climatic factor and used it to predict future changes in malaria transmission. Our findings showed that changes in the magnitude of temperatures are likely to create a shift in the distribution of malariaendemic areas than to expand its geographic ranges. This finding concurs with the results reported in [57]. A similar observation was also reported in studies by [28, 29]. It is important to note that, in cases where temperature extremes set boundaries on the vectors distributions, climate change might alter the range (in altitude or latitude) of favorable environmental conditions for the malaria vectors and the parasites. We noted that the greatest effect of climate change on malaria is likely to be observed at a temperature equal to 25 °C, corresponding to the favorable conditions for disease transmission. The optimum temperature value for malaria transmission has been estimated based on the vectors and parasites biological characteristics. With species evolution and under certain circumstances, the vectors and parasites could undergo some genetic modifications in order to adapt to current environments with different optimum temperatures. With this hypothesis, the range for the shift of malaria transmission/endemic areas will be less than predicted. In addition to the direct influence of temperature on the biology of vectors and parasites, changing precipitation patterns could also have short and longterm effects on vector development [17]. A recent analysis of global mean surface precipitation over the period 1901–1995 indicates that precipitation trends vary across Africa. Precipitation appears to be increasing in east Africa but decreasing in western and northern part of the continent [17]. These are common observations that if actualized, may reduce the level of certainty on our findings.
The study reported in [20] assessed the contribution of different malaria interventions and revealed a considerable fall in malaria burden in subSaharan Africa beginning to the year 1980. The decline of malaria burden is partly due to a number of vector control tools [insecticidetreated nets (ITNs), indoor/aerial spraying and other], which have been developed and extensively used. Improvement in diagnosis and treatment has also contributed to the decline of the malaria burden. However, the hypothesis supporting the introduction of interventions here emphasized more on the vector control aspect, which is justified by selecting model parameters (mosquito biting rate, vector competence, adult mosquito mortality rate, and the probability that mosquito eggs survive to become adult) that directly affect malaria vectors for manipulations. Depending on the disease ecology and the policies of the government in each country, the control of malaria often uses multiple intervention measures packaged in an integrated vector management strategy. With the developed model, it was found that if the rate of interventions continues as present, a considerable reduction of malaria transmission is likely to happen by the year 2050 in Africa. However, focusing at an individual country such as Cameroon, in which the latest report is dated 2004 [41]; our findings only revealed a little decline in the prevalence of malaria. Such an outcome may be attributed to the lack of adoption and failure to implement the Abuja plan of action [40, 41] and the low number of projects and research activities directed towards reducing malaria burden as compared to other regions of Africa. On the contrary the application of intervention measures has substantially contributed to the reduction of the burden of malaria in several areas of Kenya, with up to 70% decline in malaria morbidity. An identical projection was obtained by the current study. Moreover, it is important to acknowledge that the existence of diverse malaria environments in Africa, each requiring focal intervention packages to achieve success in disease control. Herein, interventions were applied identically for each country creating a possibility of bias in the overall results.
In summary, by applying the basic reproduction number R_{0} derived from a dynamical model in patches environment, the map of malaria transmission intensity was obtained. The findings in this research could constitute a realistic basis for understanding the interactions and complexities between the disease (malaria), its vectors and the parasites. Including interventions in the analysis allowed to measure the level at which, continued efforts made by different governments in Africa have so far contributed to the reduction of the malaria burden. Considerable, complementary and concurrent efforts are still needed in the drive towards malaria eradication, especially in the context of climate change.
References
 1.
Steketee RW. Measuring malaria transmission reduction en route to elimination. s. 2014;13:O26.
 2.
CDC. CDC (Centers for Disease Control and Prevention)—Malaria—Travelers—Malaria Map Application. 2015.
 3.
WHO. WHO (World Health Organization) Global Malaria Programme: World malaria Report 2016. (http://www.who.int/malaria/media/worldmalariareport2016/en/).
 4.
Ruan S, Wolkowicz GSK, Wu J, et al. Dynamical systems and their applications in biology. Providence: American Mathematical Society; 2003.
 5.
Gillies MT, Coetzee M. A supplement to the Anophelinae of Africa South of the Sahara. Publ S Afr Inst Med Res. 1987;55:1–143.
 6.
Bailey NTJ. The mathematical theory of infectious diseases and its applications. Bucks: Charles Griffin & Company Ltd; 1975.
 7.
Reluga TC, Medlock J, Galvani AP. A model of spatial epidemic spread when individuals move within overlapping home ranges. Bull Math Biol. 2006;68:401–16.
 8.
Keeling MJ, Eames KTD. Networks and epidemic models. J R Soc Interface. 2005;2:295–307.
 9.
Eames KTD, Read JM, Edmunds WJ. Epidemic prediction and control in weighted networks. Epidemics. 2009;1:70–6.
 10.
Juher D, Ripoll J, Saldaña J. Analysis and Monte Carlo simulations of a model for the spread of infectious diseases in heterogeneous metapopulations. Phys Rev E. 2009;80:041920.
 11.
Mabaso ML, Vounatsou P, Midzi S, Da Silva J, Smith T. Spatiotemporal analysis of the role of climate in interannual variation of malaria incidence in Zimbabwe. Int J Health Geogr. 2006;5:20.
 12.
Depinay JMO, Mbogo CM, Killeen G, Knols B, Beier J, Carlson J, et al. A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malar J. 2004;3:29.
 13.
Mordecai EA, Paaijmans KP, Johnson LR, Balzer C, BenHorin T, Moor E, et al. Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol Lett. 2013;16:22–30.
 14.
Parham EP, Edwin M. Modelling climate change and malaria transmission. Adv Exp Med Biol 2010;673:184–99.
 15.
Hartemink NA, Randolph SE, Davis SA, Heesterbeek JAP. The basic reproduction number for complex disease systems: defining R0 for tickborne infections. Am Nat. 2008;171:743–54.
 16.
Breman JG, Henderson DA. Diagnosis and management of smallpox. N Engl J Med. 2002;346:1300–8.
 17.
Githeko AK, Lindsay SW, Confalonieri UE, Patz JA. Climate change and vectorborne diseases: a regional analysis. Bull World Health Organ. 2000;78:1136–47.
 18.
Martens WJM, Jetten TH, Rotmans J, Niessen LW. Climate change and vectorborne diseases: a global modelling perspective. Glob Environ Change. 1995;5:195–209.
 19.
Tanser FC, Le Sueur D. The application of geographical information systems to important public health problems in Africa. Int J Health Geogr. 2002;1:4.
 20.
Caminade C, Kovats S, Rocklov J, Tompkins AM, Morse AP, ColónGonzález FJ, et al. Impact of climate change on global malaria distribution. Proc Natl Acad Sci. 2014;111:3286–91.
 21.
Yeshiwondim AK, Gopal S, Hailemariam AT, Dengela DO, Patel HP. Spatial analysis of malaria incidence at the village level in areas with unstable transmission in Ethiopia. Int. J. Health Geogr. 2009;8:5.
 22.
SouzaSantos R, de Oliveira MV, Escobar AL, Santos RV, Coimbra CE. Spatial heterogeneity of malaria in Indian reserves of Southwestern Amazonia, Brazil. Int J Health Geogr. 2008;7:55.
 23.
Kienberger S, Hagenlocher M. Spatialexplicit modeling of social vulnerability to malaria in East Africa. Int J Health Geogr. 2014;13:29.
 24.
Stefani A, Roux E, Fotsing JM, Carme B. Studying relationships between environment and malaria incidence in Camopi (French Guiana) through the objective selection of bufferbased landscape characterisations. Int J Health Geogr. 2011;10:65.
 25.
Kabaria CW, Molteni F, Mandike R, Chacky F, Noor AM, Snow RW, et al. Mapping intraurban malaria risk using high resolution satellite imagery: a case study of Dar es Salaam. Int J Health Geogr. 2016;15:26.
 26.
Cianci D, Hartemink N, IbáñezJusticia A. Modelling the potential spatial distribution of mosquito species using three different techniques. Int J Health Geogr. 2015;14:10.
 27.
Kilama M, Smith DL, Hutchinson R, Kigozi R, Yeka A, Lavoy G, et al. Estimating the annual entomological inoculation rate for Plasmodium falciparum transmitted by Anopheles gambiae using three sampling methods in three sites in Uganda. Malar J. 2014;13:111.
 28.
Tonnang HE, Kangalawe RY, Yanda PZ. Review—predicting and mapping malaria under climate change scenarios: the potential redistribution of malaria vectors in Africa. Malar J. 2010;9:1–10.
 29.
Tonnang HE, Tchouassi DP, Juarez HS, Igweta LK, Djouaka RF. Zoom in at African country level: potential climate induced changes in areas of suitability for survival of malaria vectors. Int J Health Geogr. 2014;13:12.
 30.
Macdonald G, et al. The epidemiology and control of malaria. London: Oxford University Press; 1957.
 31.
McKenzie FE, Samba EM. The role of mathematical modeling in evidencebased malaria control. Am J Trop Med Hyg. 2004;71:94–6.
 32.
Magori K, Drake JM. The population dynamics of vectorborne diseases. Nat Educ Knowl. 2013;4:14.
 33.
Auger P, Kouokam E, Sallet G, Tchuente M, Tsanou B. The Ross–Macdonald model in a patchy environment. Math Biosci. 2008;216:123–31.
 34.
Van Den Driessche P, Watmough J. Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002;180:29–48.
 35.
Uppala SM, Kållberg PW, Simmons AJ, Andrae U, Bechtold V, Fiorino M, et al. The ERA40 reanalysis. Q J R Meteorol Soc. 2005;131:2961–3012.
 36.
Blanford JI, Blanford S, Crane RG, Mann ME, Paaijmans KP, Schreiber KV, Thomas MB. Implications of temperature variation for malaria parasite development across Africa. Sci Rep. 2013;3:1300. https://doi.org/10.1038/srep01300.
 37.
QGIS Development Team. QGIS geographic information system. Open Source Geospatial Foundation. 2009. http://qgis.osgeo.org.
 38.
Solomon S. Climate change 2007—the physical science basis: working group I contribution to the fourth assessment report of the IPCC. Cambridge: Cambridge University Press; 2007.
 39.
O’Meara WP, Mangeni JN, Steketee R, Greenwood B. Changes in the burden of malaria in subSaharan Africa. Lancet Infect. Dis. 2010;10:545–55.
 40.
Malaria RB, Organization WH, et al. The Abuja Declaration and the plan of action. An extract from the African summit on roll back malaria. Abuja, 25 April 2000.
 41.
Omumbo JA, Noor AM, Fall IS, Snow RW. How well are malaria maps used to design and finance malaria control in Africa. PLoS ONE. 2013;8:e53198.
 42.
Craig MH, Snow RW, Le Sueur D. A climatebased distribution model of malaria transmission in subSaharan Africa. Parasitol Today. 1999;15:105–11.
 43.
Kazembe LN, Kleinschmidt I, Holtz TH, Sharp BL. Spatial analysis and mapping of malaria risk in Malawi using pointreferenced prevalence of infection data. Int J Health Geogr. 2006;5:41.
 44.
Mandal S, Sarkar RR, Sinha S. Mathematical models of malaria—a review. Malar J. 2011;10:10–1186.
 45.
Smith DL, McKenzie FE, Snow RW, Hay SI. Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007;5:e42.
 46.
Gumel AB. Causes of backward bifurcations in some epidemiological models. J Math Anal Appl. 2012;395:355–65.
 47.
Buonomo B. A note on the direction of the transcritical bifurcation in epidemic models. Nonlinear Anal Model Control. 2015;20:38–55.
 48.
Buonomo B, Cerasuolo M. Stability and bifurcation in plant–pathogens interactions. Appl Math Comput. 2014;232:858–71.
 49.
Gething PW, Elyazar IR, Moyes CL, Smith DL, Battle KE, Guerra CA, et al. A long neglected world malaria map: Plasmodium vivax endemicity in 2010. PLoS Negl Trop Dis. 2012;6(9):e1814.
 50.
Riley S, Eames K, Isham V, Mollison D, Trapman P. Five challenges for spatial epidemic models. Epidemics. 2015;10:68–71.
 51.
Hay SI, Guerra CA, Gething PW, Patil AP, Tatem AJ, Noor AM, et al. A world malaria map: Plasmodium falciparum endemicity in 2007. PLoS Med. 2009;6:286.
 52.
Hartemink NA, Purse BV, Meiswinkel R, Brown HE, De Koeijer A, Elbers ARW, et al. Mapping the basic reproduction number (R 0) for vectorborne diseases: a case study on bluetongue virus. Epidemics. 2009;1:153–61.
 53.
Wu X, Duvvuri VR, Lou Y, Ogden NH, Pelcat Y, Wu J. Developing a temperaturedriven map of the basic reproductive number of the emerging tick vector of Lyme disease Ixodes scapularis in Canada. J Theor Biol. 2013;319:50–61.
 54.
Cordovez JM, Rendon LM, Gonzalez C, Guhl F. Using the basic reproduction number to assess the effects of climate change in the risk of Chagas disease transmission in Colombia. Acta Trop. 2014;129:74–82.
 55.
Anyamba A, Chretien JP, Small J, Tucker CJ, Linthicum KJ. Developing global climate anomalies suggest potential disease risks for 2006–2007. Int J Health Geogr. 2006;5:60.
 56.
Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IR, Johnston GL, et al. A new world malaria map: Plasmodium falciparum endemicity in 2010. Malar J. 2011;10:1475–2875.
 57.
Lafferty KD. The ecology of climate change and infectious diseases. Ecology. 2009;90:888–900.
Authors’ contributions
FMMK, HEZT designed the research, FMMK, RYAG, FTN, JM, HEZT conducted the analysis, FMMK, RYAG, FTN, HEZT, SAP, JM wrote the manuscript. All authors contributed to the final manuscript and approved its submission. All authors read and approved the final manuscript.
Acknowledgements
Francois M. Moukam Kakmeni would like to express his gratitude to the Deutscher Akademischer Austausch Dienst (DAAD) and the International Centre of Insect Physiology and Ecology (ICIPE), for financial and material support in this research.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The paper used already published information.
Consent for publication
All authors approved its submission.
Ethics approval and consent to participate
Not applicable.
Funding
No real research funding was obtained to carry out this study; only the first author benefited from a fellowship from DAAD for his upkeep.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Vectorborne disease
 Transmission
 Network model
 Basic reproduction number
 Geographical information system (GIS)