Skip to main content

Relative effects of climate factors and malaria control interventions on changes of parasitaemia risk in Burkina Faso from 2014 to 2017/2018

Abstract

Background

In Burkina Faso, the prevalence of malaria has decreased over the past two decades, following the scale-up of control interventions. The successful development of malaria parasites depends on several climatic factors. Intervention gains may be reversed by changes in climatic factors. In this study, we investigated the role of malaria control interventions and climatic factors in influencing changes in the risk of malaria parasitaemia.

Methods

Bayesian logistic geostatistical models were fitted on Malaria Indicator Survey data from Burkina Faso obtained in 2014 and 2017/2018 to estimate the effects of malaria control interventions and climatic factors on the temporal changes of malaria parasite prevalence. Additionally, intervention effects were assessed at regional level, using a spatially varying coefficients model.

Results

Temperature showed a statistically important negative association with the geographic distribution of parasitaemia prevalence in both surveys; however, the effects of insecticide-treated nets (ITNs) use was negative and statistically important only in 2017/2018. Overall, the estimated number of infected children under the age of 5 years decreased from 704,202 in 2014 to 290,189 in 2017/2018. The use of ITNs was related to the decline at national and regional level, but coverage with artemisinin-based combination therapy only at regional level.

Conclusion

Interventions contributed more than climatic factors to the observed change of parasitaemia risk in Burkina Faso during the period of 2014 to 2017/2018. Intervention effects varied in space. Longer time series analyses are warranted to determine the differential effect of a changing climate on malaria parasitaemia risk.

Peer Review reports

Background

Malaria is a major global health problem affecting hundreds of millions of people every year and accounted for 619,000 deaths in 2021 [1]. Sub-Saharan Africa bears the largest burden contributing more than 80% of the total cases, mainly children under the age of 5 years [2]. Over the past two decades, an increasing number of partners and financial resources have enhanced malaria control efforts. The scale-up of interventions has saved millions of lives globally and cut malaria mortality by an estimated 36% from 2010 to 2020 [3]. Yet, malaria continues to cause severe morbidity and mortality, particularly in six sub-Saharan African countries, including Burkina Faso [2, 4]. Indeed, Burkina Faso contributed to 3.3% of global malaria cases in 2021 [1]. In 2020, 11,311,560 people were affected with a diagnostic confirmation rate of 92%, including 508,097 severe cases and 3,966 deaths [5]. As in many other malaria endemic countries, Burkina Faso mitigates the malaria burden through interventions such as insecticide-treated nets (ITNs), artemisinin-based combination therapies (ACTs), indoor residual spraying (IRS) and chemoprophylaxis. National ITN distribution campaigns and introduction of seasonal malaria chemoprevention (SMC) in children aged 36–59 months were carried out by the Ministry of Health since 2010 [6, 7]. However, malaria remained the main cause of outpatient attendance at health centres in all age groups in 2020 [6].

In Burkina Faso, malaria transmission is seasonal with a peak of incidence observed during and just after the rainy season [8]. This seasonal peak varies across the three major geographic zones linked with the duration of the rainy season; for instance, up to 3 months in the North, 6 months in the Centre, and 9 months in the South of the country [9]. Studies in different parts of the world have shown a relative impact of climatic factors on malaria, but the effect of climate variables varies according to regions. The development of malaria parasites in the mosquito depends on climatic and environmental factors [10]. As regards climatic factors, temperature and humidity are particularly relevant [10]. Environmental changes may result in conditions that favour the spread of malaria, offsetting the impact of preventive measures, such as ITNs, IRS, ACT and SMC, thus putting more people at risk of malaria [11]. Altitude and rainfall have been significantly associated with malaria cases in south Sumatra, Indonesia [12]. Beloconi et al. (2023), using transmission models, have reported that temperature had a significant negative effect on malaria incidence in the lowland of Kenya [13]. Another study in West Africa, using a Bayesian modelling approach, has shown that the average annual temperature had negative association with malaria incidence [14]. It has also been demonstrated that, besides the important effects of interventions, living in urban areas was associated with a lower odds of malaria, while rainfall, land surface temperature (LST) day and night were associated with a higher odds in Uganda [15]. Studies conducted in Kenya and Uganda found that LST and rainfall showed a negative and positive association with malaria, respectively [16, 17]. Using Bayesian geostatistical modelling, previous studies have identified several environmental factors that were significantly associated with malaria transmission in Burkina Faso. For example, Samadoulougou et al. (2014), have reported that rainfall and maximum temperature were respectively positively and negatively associated with malaria after accounting for spatial dependence [18]. Likewise, Diboulo et al. (2016) and Rouamba et al. (2019) found that night LST and rainfall, respectively, showed a negative and positive association with malaria [19, 20].

While there is a growing body of literature studying the impact of climatic factors on malaria using spatial and temporal models at the individual level, the relationship between changes in climate and interventions on malaria is not well established. In this study, we estimated the differential effect of malaria interventions and climatic factors on the temporal changes of malaria parasite prevalence of children under the age of 5 years. The analysis is based on data from two Malaria Indicator Surveys (MIS) conducted in 2014 and 2017/2018, employing Bayesian geostatistical modelling. The effects of malaria interventions are estimated at country and regional levels, readily adjusted for climatic confounders. Furthermore, we predicted the geographical distribution of parasitaemia risk at high spatial resolution in 2014 and 2017/2018 and the posterior probability of reduction in parasitaemia risk between the two surveys and of the number of infected children.

Methods

Country profile

Burkina Faso is a landlocked country in West Africa, which belongs to the Sahel region. It is mainly located on a savannah plateau. On average, Burkina Faso receives between 500 mm and 1,000 mm of rainfall mainly between June and September and has an average annual temperature of around 28 °C [21]. The country had a population slightly more than 20 million in 2020, all of whom are at risk of malaria [9]. Malaria was responsible for 60% of the hospitalisations and was the largest contributor to mortality for children under the age of 5 years in 2021 [9, 22]. Plasmodium falciparum is the parasite responsible for 90% of malaria cases and Anopheles gambiae is the most efficient malaria vector species in Burkina Faso [23]. The LST day and LST night have increased in Burkina Faso over the past decade (see Fig. 1B in additional file 1). In recent years, precipitation increased from May to the end of August, while the temperature (day and night) decreased. However, from September to November, precipitation decreased and temperature increased in the recent past (see Fig. 1B in additional file 1).

Survey data

The National Institute for Statistics and Demography (INSD) of Burkina Faso conducted a first MIS from October to November 2014 and a second MIS from November 2017 to March 2018. Hence, the MIS 2014 was carried out shortly after the end of the rainy season, while the MIS 2017/2018 in the dry season. For each survey, a random sample of 252 clusters (200 in rural and 52 in urban settings) with 6,552 households were selected through a stratified two-stage sampling procedure. The clusters were the census units established by INSD in 2006 [24]. Each region is divided into urban and rural areas to form the two sampling clusters within the region and the samples were drawn independently in each cluster. At the first stage, 252 clusters were drawn with a probability proportional to the size of households in the cluster during the mapping operation for the 2006 census. At the second sampling stage, a count of households in each of these 252 clusters provided a list of households from which 26 households were selected with equal probability systematic sampling. As part of the final sampling, one in two households were randomly selected and all children under the age of 5 years tested for parasitaemia. Two malaria diagnostic tests were performed; (i) a rapid diagnostic test (RDT) and (ii) a blood film examined under a microscope [25, 26]. Our data analysis included 224 clusters for 2014 and 194 clusters for 2017/2018 that had geo-referenced information available. The number of households used was 6,449 in 2014 and 6,318 in 2017/2018. The surveys included information about interventions and ACT at the cluster locations. We created intervention coverage indicators, namely the proportion of children and adults using ITNs, the proportion of household owning ITNs and the proportion of children under the age of 5 years with fever 2 weeks prior to the survey treated by ACT.

Climatic data

Climatic data were downloaded from remote sensing sources for the full year prior to the end of the surveys, i.e. from November 2013 to November 2014 for MIS 2014 and from March 2017 to March 2018 for MIS 2017/2018. Day and night LST were extracted from Moderate Resolution Imaging Spectroradiometer (MODIS) at a spatial resolution of 1 × 1 km2 and a temporal resolution of 8 days. Night light, distance to permanent water and crops coverage were obtained from MODIS at 1 × 1 km2 spatial resolution and a temporal resolution of 12 months. Dekadal rainfall data were obtained from the US Early Warning and Environmental Monitoring System (EWES) at 8 × 8 km2 spatial resolution and a temporal resolution of 5 days. The climatic data were extracted at their original spatial and temporal resolution and summarised by their average for the year prior to the end of the surveys. They were matched spatially to the MIS locations and to the locations of the predictions. Altitude was obtained from the Shuttle Radar Topographic Mission (SRTM) using the digital elevation model. The high spatial resolution population data were downloaded from WorldPop at 1 × 1 km2 spatial resolution. Table 1 describes the sources and the spatio-temporal resolution of the data used.

Table 1 Description of the sources, spatio-temporal resolution of the climatic data for the period of 2014 to 2017 of Burkina Faso

Statistical analysis

Descriptive statistics were employed to examine the characteristics of the data and for summarising the variables. Prior to fitting the Bayesian model, collinearity between all pairs of predictors was assessed using non-spatial regression methods based on values of variance inflation factor (VIF); predictors with VIF > 5 were excluded from further analysis.

Bayesian geostatistical bivariate logistic regression was used first to assess the association between all predictors and parasitaemia risk. We assumed that the number of children under the age of 5 years, tested positive by blood microscopy, follows a binomial distribution. A multivariate model was developed to predict the prevalence at the two survey years using only the important climatic factors and interventions, as predictors that were identified by the bivariate analysis. We standardised the predictors by subtracting their mean and dividing the difference by their standard deviation (SD) to facilitate comparison of their effects at the country level. The spatial correlation was taken into account by introducing location random effects. They were modelled by a multivariate normal distribution with an exponential function of distance between any pair of locations. Furthermore, the model allowed us to estimate the minimum distance at which the spatial correlation is lower than 5% (range parameter). Bayesian kriging was carried out to predict parasitaemia risk over a 2 × 2 km2 resolution grid. The posterior probability of a decline in parasitaemia risk between the two survey years was determined using samples from the corresponding posterior predictive distributions (PPD). We draw a sample of size 100 from each PPD at grid cell level and computed their differences (i.e. predicted parasitaemia risk in 2017/2018 minus parasitaemia risk in 2014). The proportion of samples with a positive difference at each grid cell estimates the posterior probability of a decline in parasitaemia risk.

We assessed the effects of the changes in the coverage of interventions and climatic factors on changes in parasitaemia prevalence between the two surveys by fitting a Bayesian geostatistical model to the malaria survey data of 2017/2018. The changes in interventions (i.e. the differences in 2017/2018 coverage compared to the 2014 coverage) and the changes in climatic factors (i.e. differences in 2017/2018 values from those in 2014) were included as covariates in the logit of parasitaemia risk. The logit of parasitaemia in 2014 was considered as an offset term to enable changes in the risk between the survey years to be modelled on the logit scale. The geographical misalignment of the locations between the two surveys was taken into account by predicting the prevalence of the first survey at the locations of the second survey.

The number of children infected with Plasmodium in each survey year was estimated by multiplying population data obtained from WorldPop [27] at 1 km2 spatial resolution with the predicted pixel-level malaria parasitaemia prevalence. The number of children under the age of 5 years in a specific pixel was estimated by multiplying the population counts by 0.183 [28] and 0.180 [29], which are the proportions of population under the age of 5 years in 2014 and 2017/2018, respectively. Regional and national estimates of the number of infected children under 5 years old were obtained by aggregating the pixel-level estimates within each region and across the country. The relative reduction of the number of infected children was calculated by dividing the difference of the numbers of infected of that age group between the two surveys by the number of infected in the first survey. The estimates of the total number of infected children under 5 years old, both by region and for the entire country, were divided by the total population of that age group. This allowed us to obtain population-adjusted prevalence estimates at both regional and national levels, respectively.

A spatially varying coefficients model was fitted to estimate the effects of malaria interventions at a regional level. A conditional autoregressive (CAR) prior distribution was considered to introduce a neighbour-based spatial structure on the regression coefficients related to each intervention effect in the study. Neighbours were defined as the adjacent areas regions and used to create a matrix of spatial weights taking the value of one for a direct neighbour and zero otherwise [19, 30]. We produced maps of the random effects of the interventions. The random effect was estimated as the difference between the regional effect and the national effect. Following this approach, a negative value indicates that the effect of intervention at regional level is statistically lower than the national average.

Geostatistical models were fitted in Integrated Nested Laplace Approximation (INLA) [31] package in R (version 4.1.2) [32]. The spatially varying coefficients model was implemented in OpenBugs (version 3.2.3, Imperial College and Medical Research Council; London, UK) [33]. We summarised the parameter estimates by computing the posterior medians of odds ratio (OR) and their corresponding 95% Bayesian credible intervals (BCIs). The effect of a covariate was considered statistically important when the 95% BCI did not include one. A detailed description of the analysis is provided in supplementary file (see Additional file 2).

Results

Descriptive analysis

A description of the survey data and climatic factors is given in Table 2. Maps of the observed prevalence for both surveys are shown in Fig. 1. The number of screened children under the age of 5 years were 5,082 and 4,417 in the 2014 and 2017/2018 MIS, respectively. The prevalence of malaria was reduced from 44.7% in 2014 to 16.2% in 2017/2018. The lowest prevalence was observed in the Centre (9.3% in 2014 and 6.6% in 2017/2018) and the highest in the region of South-West (60.2% in 2014 and 39.0% in 2017/2018) (see Fig. 1A in additional file 1, and Table 3A in additional file 3). Among the regions, the highest decrease in the prevalence of malaria from 2014 to 2017/2018 occurred in the region of North with a relative reduction of 75.7%. The proportion of the population that slept under an ITN the night before the survey decreased by 33.7% comparing 2014 with 2017/2018. However, the proportion of fever episodes treated with ACT increased by 65.2%.

Fig. 1
figure 1

Observed malaria prevalence and survey locations of MIS 2014 (A) and 2017/2018 (B); predicted parasitaemia risk in 2014 (C) and 2017/2018 (D); distribution of estimated number of infected children per km2 pixel in 2014 (E) and 2017/2018 (F); posterior probability of reduction in parasitaemia risk (G); location diagram of Burkina Faso (H)

Table 2 Survey information, climatic factors and malaria intervention coverage in 2014 and 2017/2018 in Burkina Faso

Effect of interventions and climatic factors on the geographical distribution of parasitaemia risk from 2014 to 2017/2018

Results from the geostatistical bivariate model (Table 3) showed that the relationship between the intervention coverage indicators and the malaria prevalence was not statistically important in 2014, except of the proportion of households with at least one ITN for every two people. However, there is evidence that in 2017/2018, ITN ownership, ITN use and ACTs were associated with a decrease of parasitaemia risk. Night LST, night light and distance to permanent water were negatively associated with parasitaemia risk in both periods. Altitude was statistically important in 2014 but not in 2017/2018.

The final multivariate geostatistical models showed that more climatic factors were associated with the geographical distribution of parasitaemia risk in 2014 than in 2017/2018 (Table 4). In 2014, none of the intervention indicators included in the analysis was statistically important. In 2017/2018, only the proportion of children who slept under an ITN was important with an OR of 0.82 (95% BCI: 0.70–0.96). The spatial range parameter was smaller, while the spatial variance was higher in 2014 compared to 2017/2018. The LST night was statistically important for both time points with an OR of 0.72 (95% BCI: 0.64–0.93) in 2014 and 0.89 (95% BCI: 0.81–0.96) in 2017/2018. There were no statistically important interactions between malaria interventions and LST night (see Table 3B in additional file 3). Night light was also negatively associated to parasitaemia risk in 2014 (OR = 0.47, 95% BCI: 0.38–0.59) and 2017/2018 (OR = 0.49, 95% BCI: 0.37–0.64). An altitude of above 307 m was associated with a 28% lower odds of parasitaemia in 2014. Figure 1C and D depict the predicted prevalence in 2014 and 2017/2018 over a 2 × 2 km2 resolution grid. In 2014, we observed a high-predicted parasitaemia risk (above or equal to 50%) in some areas, such as South-West, Centre-West and Hauts-Bassins. The posterior probability of reduction in parasitaemia risk was more than 0.50 in most part of the country. The highest decline in malaria prevalence was reported in the regions of Centre-East, Plateau-Central and Centre, while the least decline occurred in the South-West region (Fig. 1G).The number of infected children decreased from 704,202 in 2014 to 290,189 in 2017/2018 (Table 5), corresponding to a relative reduction in the country of 58.8%. The number of infected children was reduced in all regions (Fig. 1E and F). Reduction rates above the national average were observed in the regions of Centre-East, Hauts-Bassins, Boucle du Mouhoun, North and Centre-Southern, whereas the lowest was in the South-West (38.3%) (Table 5). In both surveys, the highest and lowest estimated numbers of children infected were observed in the regions of Hauts-Bassins and Centre, respectively. In addition, a reduction of population-adjusted prevalence of 63.1% was observed overall.

Effects of changes in intervention coverage and climatic factors on the decline in parasitaemia prevalence from 2014 to 2017/2018

The effects of changes in intervention coverage indicators and climatic factors between the two surveys on the parasitaemia risk decline are shown in Table 6. The results suggest that a 1% increase in the proportion of children who slept under an ITN was associated with a decrease of parasitaemia odds by 24%. The change in ACT coverage and in the proportion of households with at least one ITN for every two persons showed no statistically important effects in the parasitaemia risk decline. A one-unit increase in the change in night light was associated with a 15% lower odds of parasitaemia between the two surveys. Changes in LST night and distance to permanent water between 2014 and 2017/2018 were not associated with a decrease of parasitaemia risk. The residual spatial variation in the changes of parasitaemia risk accounted for 71% of total variation. The geographical distribution of the effects of interventions at regional level (Fig. 2) indicates that ITN ownership had a statistically important effect lower than the national effect in the regions of South-West, Boucle du Mouhoun, Centre and East. In North, Sahel, Centre-North, Centre-West, Centre-South and Hauts-Bassins, this effect was stronger than the national average. The effect of the proportion of children who slept under an ITN was statistically different from the national average in all regions, except for the regions of North and East. ACT was statistically important and more effective in the regions of Hauts-Bassins, Boucle du Mouhoun and Centre-South. ACT was less effective than nationally in the regions of South-West and Centre-East.

Fig. 2
figure 2

Effect of proportion of households with at least one ITN for every two people at regional level (A); effect of proportion of children who slept under an ITN at regional level (B) and effect of ACT coverage at regional level (C)

*Statistically important effect less than national effect

+Statistically important effect higher than national effect

Table 3 Posterior odds ratios (ORs) and 95% Bayesian credible intervals (BCIs) estimated by bivariate geostatistical logistic regression models, fitted to MIS 2014 and 2017/2018 data from Burkina Faso
Table 4 Posterior odds ratios (ORs) and 95% Bayesian credible intervals (BCIs) estimated by multivariate geostatistical logistic regression models fitted to MIS 2014 and 2017/2018 data from Burkina Faso
Table 5 Estimated number of children infected with Plasmodium and population adjusted prevalence in 2014 and 2017/2018 in Burkina Faso
Table 6 Posterior estimates of the effect of changes in interventions and climatic factors on parasitaemia risk decline from 2014 to 2017/2018 in Burkina Faso

Discussion

This study assessed the effects of malaria control interventions and climatic factors simultaneously on the geographical distribution and on the change of malaria parasitaemia risk in Burkina Faso between 2014 and 2017/2018. ITN use had a statistically important effect on the geographical distribution of parasitaemia risk in 2017/2018 but not in 2014. ITN use was also associated with parasitaemia risk decline from 2014 to 2017/2018. Similar results have been reported before [34, 35]. A study conducted in Africa reported also that bednets were effective, mainly because they prevent Anopheles mosquitoes from feeding blood, even in case nets are damaged and also because of their capacity to eliminate the disease-carrying vectors [36]. Bednets act as a physical barrier between mosquitoes and the people sleeping under the nets. The nets are treated to repel mosquitoes or kill them when they land on the nets [37]. Bednets have been found to have a significant effect on malaria transmission [38, 39] in children under 5 years of age in additional research conducted in Africa and Papua New Guinea.

The statistical important effect of a decreased coverage of ITNs use between the two time periods could be explained by the fact that the good use of ITNs protect not only individual protection but also contribute to community level protection by reducing the overall population of mosquitoes in an area. This observation is supported by previous studies in Malawi [40] and Mozambique [41].

There was a variation in the effect of ITN across the country. Previous studies in Burkina Faso [19], in Haiti [42] and Uganda [30] reported varying effect of ITN use. A study has showed that even with bednets usage by the entire population, the reduction of malaria transmission would be constrained to a certain minimum level. This is due to the fact that bednets are typically used during the night [43], and studies have documented instances of early mosquito bites occurring when people are outdoors [42, 44]. Interventions’ effects might be affected by many factors, such as human behaviour, mosquito behaviour, pyrethroid vector resistance or the physical integrity of bednets. It should also be noted that bednets have a variable lifespan in terms of physical condition, residual insecticidal effect and perceived usefulness, depending on the geographical and cultural context [34]. For instance, bednets can be effective for up to 20 washings, but their efficacy can be significantly reduced if certain measures are not taken, such as avoiding the use of bleach or hot water when washing them, and avoiding drying them in direct sunlight [45]. Hence, low compliance of bednets use may reduce its potential impact on malaria transmission [19]. A previous study also found that bednets may not always be effective because children are vulnerable during the early evening and early morning hours when they are not protected by the nets or inside the house [46].

ACT coverage was not associated with the geographical distribution and decline of parasitaemia risk from 2014 to 2017/2018. Similar findings have been reported in Angola, Liberia, Senegal, Mozambique, Rwanda, Tanzania [34], Burkina Faso [19] and Uganda [30]. This non-protective effect of ACT might be explained by the fact that this intervention was measured by the proportion of fever episodes rather than clinical malaria cases treated by ACT. As the aetiology of fever in children under 5 years old is multifactorial, causes other than malaria must be considered, such as respiratory viral infections or bloodstream infection, and hence, the effect of ACT is not clear [47, 48]. Nevertheless, at a regional level, this intervention had a stronger effect on malaria reduction in some regions comparing to the national effect. Other studies in Tanzania [49] and Uganda [17] also revealed an important and negative association of ACTs and malaria prevalence. Similarly, Diboulo et al. (2016), in a study conducted in Burkina Faso, have reported that ACT coverage was not important at national level but was statistically important at regional level. This finding could be explained by the fact that the effectiveness of a given intervention is related to both its coverage and transmission levels [19]. The effect of IRS could not be estimated because the survey of 2017/2018 did not have coverage related information.

Rainfall was not found to be an important predictor for the observed change of parasitaemia risk. This result might be related to the data collection period, which was at the beginning of dry season in 2017/2018 and shortly after the end of the rainy season in 2014. Although Roll Back Malaria (RBM) has recommended that surveys take place during the highest transmission season (rainy season), this is not always feasible because of logistical issues and a scarcity of human resources to reach remote areas [34]. LST night was negatively associated with parasitaemia risk in 2014 and 2017/2018. Similar results have been reported before from Burkina Faso [18, 19] and Kenya [13, 16]. This finding might be explained through the optimum temperature at which the mosquito and the parasite develop. Both Anopheles and Plasmodium are sensitive to temperature. Because mosquitoes are ectotherms, each life stage is dependent on temperature in the developmental and mortality rates [50]. Agyekum et al. (2021) reported that the longevity and survival rate of An. gambiae mosquitoes were higher in the rainy season, which is associated with cooler temperatures (17.5 ± 2.9 days) and higher humidity (84.5% ± 10.5%), than in the dry season which is associated with hotter temperatures and lower humidity (7.3 ± 2.8 days and 57.5% ± 14.9%, respectively) [51]. In effect, Anopheles infection rates depend on the duration of the reproduction cycle of the parasites in the mosquito, which in turn is related to the temperature [52]. Increasing temperature decreases the longevity of mosquitoes and increases mosquito mortalities, and hence, the transmission decreases. The relationship between temperature and longevity could be explained in two ways. First, higher temperatures may decrease the longevity by speeding the reaction rate of various metabolic processes that affect development and life history. Second, higher temperatures might heighten the damage caused by the by-products of metabolism, such as reactive oxygen species [51]. It has also been reported that P. falciparum only develop between 15 and 35 °C, the optimum being at 25 °C. The more one deviates from this optimum, the longer the duration of the parasite’s reproduction cycle, and hence, the infection rate of Anopheles decreases [52]. Additional studies have demonstrated that malaria parasites have a temperature-dependent development curve that does not match up with the mosquito temperature curve. This suggest, for example, that if the optimal temperature for parasite development is higher than the range preferred by mosquitoes, this could result in slower development of the parasite in the mosquito, and hence, reducing the malaria transmission. In contrary, if the mosquito’s preferred temperature is higher than the optimal range for parasite development, the mosquito may survive longer, leading to an increase in the transmission rate [50]. However, other studies have reported that an increase of LST night was associated with an increase of parasitaemia risk in Uganda and Cameroon [15, 53].

An increase of night light, which is a proxy for urbanization [54], was associated with a reduction of malaria parasite prevalence. This result is consistent with previous studies in Africa [15, 19, 55] that have shown that children living in urban areas are at a lower risk of malaria compared to children residing in rural areas where An. gambiae is the predominate vector. Urbanization can change malaria transmission patterns over longer time scales by creating more or less favourable microclimatic conditions and habitats for malaria vectors, depending on the species [56, 57]. For example, urban development can eliminate mosquito habitats and reduce malaria incidence. In some cases, increases have been observed in peri-urban malaria as a result of economic migration and the creation of novel breeding habitats [58]. Urbanization affects anopheline species diversity, mosquito density, survival rates and infection rates with P. falciparum and the frequency with which mosquitoes bite people are all affected. Hence, fewer people acquire malaria infection, become ill or die of its consequences in urban areas [59]. Furthermore, urban areas have stronger networks of better equipped higher level health facilities both public and private, while rural area often lack medical facilities in term of quantity and quality.

At higher altitudes, malaria transmission decreased. Similarly, a study in Tanzania has shown that the primary effect of increasing altitude was a log-linear reduction in vector abundance and, to a lesser extent, a reduction in the proportion of infective mosquitoes [60]. In the lowland regions of sub-Saharan countries, malaria is holoendemic with perennial transmission. On the contrary, highland areas are known to be malaria hypoendemic, due to climate (low temperature and relative humidity), which is not suitable for anopheline development and their reproductive fitness [61]. At higher elevations, temperatures have also historically been too cold for Plasmodium to be transmitted [62]. Therefore, the risk of malaria parasite decreases with altitude due to an impact of temperature on the cycle of the Anopheles mosquito of the parasite. This risk is not linear due to the many other determinants that influence the epidemiology of malaria at a given altitude [63].

Our results showed a strong reduction in the number of under-5-year-old children infected with Plasmodium in Burkina Faso when comparing MIS 2014 with 2017/2018, which may underline the effect of ITN use and changing climatic conditions. The highest posterior probability of reduction in parasitaemia risk occurred in the region of Centre and the lowest in the region of South-West. In effect, the Centre is the capital and the most urbanized setting in Burkina Faso where the effect of ITN use was important. In contrast, the region of South-West is less urbanized with high endemicity and lower average temperature [64]. South-West has also a longer rainy season and is characterised by many swamps, which act as reservoirs for malaria parasites [65]. Finally, the decrease of spatial variance associated to an increase of the spatial range from 2014 to 2017/2018 can be explained by the fact that more climatic factors were linked to malaria parasite prevalence in 2014 and therefore the residual spatial structure in the model was lower compared to 2017/2018.

A main limitation in this study is that the MIS of 2014 and of 2017/2018 were carried out during the early dry and dry season, respectively. Therefore, it is likely that the difference in transmission seasons, overestimates the reduction in malaria parasitaemia prevalence.

Conclusion

Parasitaemia prevalence was considerably lower in 2014 compared to 2017/2018. During the same time the coverage of bednets interventions decreased and that of ACTs increased. Annual mean temperature remained the same but rainfall was lower mainly as the second survey was carried out in the dry season. Temperature had a statistically important negative association with the geographical distribution of parasitaemia prevalence in both surveys, however the effects of malaria interventions (i.e. proportion of children who slept under an ITN) were negative and statistically important only in 2017/2018. The latter was also related with the geographical distribution of the changes in malaria prevalence after adjusting for changes in temperature. The decline of the parasite prevalence distribution was related to the ACT coverage at regional but not national level.

The malaria control programme should pay attention to the potential hotspot regions of South-West where the posterior probability of reduction in parasitaemia risk was lower than 20%. The effectiveness of all interventions were lower than that of national average or not statistically important in South-West and Boucle du Mouhoun. To accurately evaluate the relative effects of climate and malaria interventions to malaria parasite declines, MIS surveys should be carried out at the same time of the year. In addition, longer time series analysis can help to more fully understand the differential effects of a changing climate on malaria parasite risk, to develop forecasting models, and to improve the ability to anticipate and also mitigate the future effects of climate change on malaria outbreaks.

Data availability

The study data are available upon request from the Demographic and Health Surveys programmes (https://dhsprogram.com/).

References

  1. World malaria report. 2022 [https://cdn.who.int/media/docs/default-source/malaria/world-malaria-reports/world-malaria-report-2022.pdf?sfvrsn=40bfc53a_4].

  2. World Malaria Report. 2021 [https://cdn.who.int/media/docs/default-source/malaria/world-malaria-reports/world-malaria-report-2021-global-briefing-kit-eng.pdf?sfvrsn=8e5e915_23&download=true]

  3. CDC-Malaria.: Worldwide Impact of Malaria [https://www.cdc.gov/malaria/malaria_worldwide/impact.html].

  4. Musoke D, Miiro G, Ndejjo R, Karani G, Morris K, Kasasa S, Nakiyingi-Miiro J, Guwatudde D, Musoke MB. Malaria prevention practices and associated environmental risk factors in a rural community in Wakiso district, Uganda. PLoS ONE. 2018;13(10):e0205210.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Lutte contre le paludisme.: le Burkina Faso élabore son plan national de plaidoyer Speak Up Africa. [https://www.speakupafrica.org/fr/lutte-contre-le-paludisme-le-burkina-faso-elabore-son-plan-national-de-plaidoyer/ ].

  6. Yaro JB, Tiono AB, Ouedraogo A, Lambert B, Ouedraogo ZA, Diarra A, Traore A, Lankouande M, Soulama I, Sanou A, et al. Risk of Plasmodium falciparum infection in south-west Burkina Faso: potential impact of expanding eligibility for seasonal malaria chemoprevention. Sci Rep. 2022;12(1):1402.

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  7. Millogo AA, Yaméogo L, Kassié D, de Charles Ouédraogo F, Guissou C, Diabaté A. Spatial modelling of malaria prevalence associated with geographical factors in Houet province of Burkina Faso, West Africa. GeoJournal 2022.

  8. Sangaré I, Ouattara CA, Soma DD, Soma D, Assogba BS, Namountougou M, Tougri G, Savadogo LB. Spatial-temporal pattern of malaria in Burkina Faso from 2013 to 2020. Parasite Epidemiol Control. 2022;18:e00261.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Malaria in Burkina Faso.: Statistics & Facts| Severe Malaria Observatory [https://www.severemalaria.org/countries/burkina-faso].

  10. CDC-Malaria.: About Malaria Biology [https://www.cdc.gov/malaria/about/biology/index.html].

  11. Moore S. Effect of Climate Change on Malaria. 2022. Available at: https://www.news-medical.net/health/The-Effect-of-Climate-Change-on-Malaria.aspx.

  12. Hasyim H, Nursafingi A, Haque U, Montag D, Groneberg DA, Dhimal M, Kuch U, Müller R. Spatial modelling of malaria cases associated with environmental factors in South Sumatra, Indonesia. Malar J. 2018;17(1):87.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Beloconi A, Nyawanda BO, Bigogo G, Khagayi S, Obor D, Danquah I, Kariuki S, Munga S, Vounatsou P. Malaria, climate variability, and interventions: modelling transmission dynamics. Sci Rep. 2023;13(1):7367.

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  14. Arab A, Jackson MC, Kongoli C. Modelling the effects of weather and climate on malaria distributions in West Africa. Malar J. 2014;13:1–9.

    Article  Google Scholar 

  15. Ssempiira J, Nambuusi B, Kissa J, Agaba B, Makumbi F, Kasasa S, Vounatsou P. The contribution of malaria control interventions on spatio-temporal changes of parasitaemia risk in Uganda during 2009–2014. Parasites & Vectors. 2017;10(1):450.

    Article  Google Scholar 

  16. Nyawanda B, Beloconi A, Khagayi S, Bigogo G, Obor D, Otieno N, Lange S, Franke J, Sauerborn R, Utzinger J, et al. The relative effect of climate variability on malaria incidence after scale-up of interventions in western Kenya: a time-series analysis of monthly incidence data from 2008 to 2019. Parasite Epidemiol Control. 2023;21:e00297.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Ssempiira J, Kissa J, Nambuusi B, Kyozira C, Rutazaana D, Mukooyo E, Opigo J, Makumbi F, Kasasa S, Vounatsou P. The effect of case management and vector-control interventions on space–time patterns of malaria incidence in Uganda. Malar J. 2018;17(1):1–11.

    Article  Google Scholar 

  18. Samadoulougou S, Maheu-Giroux M, Kirakoya-Samadoulougou F, De Keukeleire M, Castro MC, Robert A. Multilevel and geo-statistical modeling of malaria risk in children of Burkina Faso. Parasites & Vectors. 2014;7(1):350.

    Article  Google Scholar 

  19. Diboulo E, Sié A, Vounatsou P. Assessing the effects of malaria interventions on the geographical distribution of parasitaemia risk in Burkina Faso. Malar J. 2016;15(1):228.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Rouamba T, Nakanabo-Diallo S, Derra K, Rouamba E, Kazienga A, Inoue Y, Ouédraogo EK, Waongo M, Dieng S, Guindo A, et al. Socioeconomic and environmental factors associated with malaria hotspots in the Nanoro demographic surveillance area, Burkina Faso. BMC Public Health. 2019;19(1):249.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Burkina Faso - Climat. [https://www.pik-potsdam.de/en/institute/departments/climate-resilience/projects/project-pages/agrica/giz_climate-risk-profile_burkina-faso_en].

  22. Burkina, Faso. Target Malaria [https://targetmalaria.org/fr/where-we-operate/burkina-faso/].

  23. Démographie-INSD/Burkina INdlSedl.: Burkina Faso Enquête Démographique et de Santé et à Indicateurs Multiples (EDSBF-MICS IV) 2010. 2012.

  24. Institut National de la Statistique et de la Démographie-INSD/Burkina.: Recensement generale de la population et de l’habitation du Burkina de 2006 [www.cns.bf/IMG/pdf/RGPH_2006.pdf].

  25. Demographic and Health Surveys Program: Burkina Faso MIS 2014-MIS Final Report. 2015. [https://dhsprogram.com/pubs/pdf/MIS19/MIS19.pdf].

  26. Demographic and Health Surveys Program.: Burkina Faso MIS, 2017-18-MIS Final Report 2018. [https://dhsprogram.com/pubs/pdf/MIS32/MIS32.pdf].

  27. WorldPop. Open Spatial Demographic Data and Research. [https://www.worldpop.org/].

  28. Institut National de la Statistique et de la. Démographie:Tableau de Bord Demographique 2015. [https://fr.readkong.com/page/tableau-de-bord-demographique-u-1243866].

  29. Institut National de la Statistique et. de la Démographie: Tableau de bord social. 2018. [https://www.insd.bf/sites/default/files/2021-12/TBS_2018.pdf].

  30. Ssempiira J, Nambuusi B, Kissa J, Agaba B, Makumbi F, Kasasa S, Vounatsou P. Geostatistical modelling of malaria indicator survey data to assess the effects of interventions on the geographical distribution of malaria prevalence in children less than 5 years in Uganda. PLoS ONE. 2017;12(4):e0174948.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Blangiardo M, Cameletti M. Spatial and spatio-temporal bayesian models with. R-INLA: John Wiley & Sons; 2015.

    Book  Google Scholar 

  32. Brunsdon C, Comber L. An introduction to R for spatial analysis and mapping. Sage; 2014.

  33. Haneuse S. An Introduction to WinBUGS - PIMS Collaborative Research Group Summer School on Bayesian Modeling and Computation. Graphical models:111.

  34. Giardina F, Kasasa S, Sié A, Utzinger J, Tanner M, Vounatsou P. Effects of vector-control interventions on changes in risk of malaria parasitaemia in sub-saharan Africa: a spatial and temporal analysis. The Lancet Global Health. 2014;2(10):e601–15.

    Article  PubMed  Google Scholar 

  35. Pryce J, Richardson M, Lengeler C. Insecticide-treated nets for preventing malaria. Cochrane Database of Systematic Reviews 2018.

  36. Okumu FO, Moore SJ. Combining indoor residual spraying and insecticide-treated nets for malaria control in Africa: a review of possible outcomes and an outline of suggestions for the future. Malar J. 2011;10(1):1–13.

    Article  Google Scholar 

  37. Malaria prevention through insecticides treated nets. [https://www.malariaconsortium.org/media-downloads/802/Malaria%20prevention%20through%20insecticide%20treated%20nets].

  38. Lines J, Myamba J, Curtis C. Experimental hut trials of permethrin-impregnated mosquito nets and eave curtains against malaria vectors in Tanzania. Med Vet Entomol. 1987;1(1):37–51.

    Article  CAS  PubMed  Google Scholar 

  39. Hii J, Smith T, Vounatsou P, Alexander N, Mai A, Ibam E, Alpers M. Area effects of bednet use in a malaria-endemic area in Papua New Guinea. Trans R Soc Trop Med Hyg. 2001;95(1):7–13.

    Article  CAS  PubMed  Google Scholar 

  40. Escamilla V, Alker A, Dandalo L, Juliano JJ, Miller WC, Kamthuza P, Tembo T, Tegha G, Martinson F, Emch M. Effects of community-level bed net coverage on malaria morbidity in Lilongwe, Malawi. Malar J. 2017;16:1–9.

    Article  Google Scholar 

  41. Plucinski MM, Chicuecue S, Macete E, Colborn J, Yoon SS, Kachur SP, Aide P, Alonso P, Guinovart C, Morgan J. Evaluation of a universal coverage bed net distribution campaign in four districts in Sofala Province, Mozambique. Malar J. 2014;13(1):1–8.

    Article  Google Scholar 

  42. Steinhardt LC, St Jean Y, Impoinvil D, Mace KE, Wiegand R, Huber CS, Alexandre JSF, Frederick J, Nkurunziza E, Jean S. Effectiveness of insecticide-treated bednets in malaria prevention in Haiti: a case-control study. The Lancet Global Health. 2017;5(1):e96–e103.

    Article  PubMed  Google Scholar 

  43. Bala S, Gimba B. Global sensitivity analysis to study the impacts of bed-nets, drug treatment, and their efficacies on a two-strain malaria model. Math Comput Appl. 2019;24(1):32.

    MathSciNet  Google Scholar 

  44. Smith NR, Trauer JM, Gambhir M, Richards JS, Maude RJ, Keith JM, Flegg JA. Agent-based models of malaria transmission: a systematic review. Malar J. 2018;17(1):1–16.

    Article  CAS  Google Scholar 

  45. Boucle du Mouhoun, Zone BFC. Monthly Weather Averages and Historical Data [https://tcktcktck.org/burkina-faso/boucle-du-mouhoun].

  46. Exposure to Light May. Help Prevent Malaria Mosquitoes to Bite at Night [https://food.ndtv.com/health/short-light-treatment-may-reduce-mosquito-bites-at-night-1713970].

  47. D’Acremont V, Kilowoko M, Kyungu E, Philipina S, Sangu W, Kahama-Maro J, Lengeler C, Cherpillod P, Kaiser L, Genton B. Beyond Malaria — Causes of Fever in Outpatient Tanzanian Children. N Engl J Med. 2014;370(9):809–17.

    Article  PubMed  Google Scholar 

  48. Kiemde F, Tahita MC, Lompo P, Rouamba T, Some AM, Tinto H, Mens PF, Schallig HDFH, van Hensbroek MB. Treatable causes of fever among children under five years in a seasonal malaria transmission area in Burkina Faso. Infect Dis Poverty. 2018;7(1):60.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Froeschl G, Saathoff E, Kroidl I, Berens-Riha N, Clowes P, Maboko L, Assisya W, Mwalongo W, Gerhardt M, Ntinginya EN. Reduction of malaria prevalence after introduction of artemisinin-combination-therapy in Mbeya Region, Tanzania: results from a cohort study with 6773 participants. Malar J. 2018;17(1):1–11.

    Article  Google Scholar 

  50. Beck-Johnson LM, Nelson WA, Paaijmans KP, Read AF, Thomas MB, Bjørnstad ON. The Effect of temperature on Anopheles Mosquito Population Dynamics and the potential for Malaria Transmission. PLoS ONE. 2013;8(11):e79276.

    Article  PubMed  PubMed Central  ADS  Google Scholar 

  51. Agyekum TP, Botwe PK, Arko-Mensah J, Issah I, Acquah AA, Hogarh JN, Dwomoh D, Robins TG, Fobil JN. A systematic review of the effects of temperature on Anopheles Mosquito Development and Survival: implications for Malaria Control in a future warmer climate. Int J Environ Res Public Health. 2021;18(14):7255.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Développement. et Santé| Le moustique [https://devsante.org/articles/le-moustique].

  53. Nyasa RB, Awatboh F, Kwenti TE, Titanji VPK, Ayamba NLM. The effect of climatic factors on the number of malaria cases in an inland and a coastal setting from 2011 to 2017 in the equatorial rain forest of Cameroon. BMC Infect Dis. 2022;22(1):461.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Proville J, Zavala-Araiza D, Wagner G. Night-time lights: a global, long term look at links to socio-economic trends. PLoS ONE. 2017;12(3):e0174610.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Kabaria CW, Gilbert M, Noor AM, Snow RW, Linard C. The impact of urbanization and population density on childhood Plasmodium Falciparum parasite prevalence rates in Africa. Malar J. 2017;16(1):49.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Tucker Lima JM, Vittor A, Rifai S, Valle D. Does deforestation promote or inhibit malaria transmission in the Amazon? A systematic literature review and critical appraisal of current evidence. Philosophical Trans Royal Soc B: Biol Sci. 2017;372(1722):20160125.

    Article  Google Scholar 

  57. Wayant NM, Maldonado D, Arias ARd, Cousiño B, Goodin DG. Correlation between normalized difference vegetation index and malaria in a subtropical rain forest undergoing rapid anthropogenic alteration. Geospat Health. 2010;4(2):179–90.

    Article  PubMed  Google Scholar 

  58. Padilla JC, Chaparro PE, Molina K, Arevalo-Herrera M, Herrera S. Is there malaria transmission in urban settings in Colombia? Malar J. 2015;14(1):453.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW. Urbanization, malaria transmission and disease burden in Africa. Nat Rev Microbiol. 2005;3(1):81–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Bødker R, Akida J, Shayo D, Kisinza W, Msangeni HA, Pedersen EM, Lindsay SW. Relationship between Altitude and Intensity of Malaria Transmission in the Usambara Mountains, Tanzania. J Med Entomol. 2003;40(5):706–17.

    Article  PubMed  Google Scholar 

  61. Tchuinkam T, Nyih-Kong B, Fopa F, Simard F, Antonio-Nkondjio C, Awono-Ambene H-P, Guidone L, Mpoame M. Distribution of Plasmodium Falciparum gametocytes and malaria-attributable fraction of fever episodes along an altitudinal transect in Western Cameroon. Malar J. 2015;14(1):96.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Malaria risk zones. expand to higher elevations in Ethiopian Highlands| NOAA Climate.gov [ http://www.climate.gov/malaria_ethiopia].

  63. Chappuis F, Loutan L. Malaria titude. Médecine et hygiène 2003:1029–32.

  64. Théodore KS. Disparités régionales de la morbidité palustre des enfants de moins de cinq ans au Burkina Faso.20.

  65. La région du Sud-. ouest du Burkina Faso [https://www.burkinatourism.com/La-region-du-Sud-ouest-du-Burkina-Faso.html].

Download references

Acknowledgements

The authors are grateful to the Demographic and Health Surveys (DHS) for making the data available and all collaborating institutions whose contributions have enabled this study.

Funding

This research was funded by the Swiss National Science Foundation (project no. 310030E_186574).

Open access funding provided by University of Basel

Author information

Authors and Affiliations

Authors

Contributions

PV, AS and NT conceived the study. NT and PV analysed and interpreted the data. PV prepared the datasets. NT and TS drafted the manuscript. All authors contributed to the study design and implementation, critically reviewed the manuscript for its intellectual content and approved the final version prior to submission.

Corresponding author

Correspondence to Penelope Vounatsou.

Ethics declarations

Ethics approval and consent to participate

In this study we analysed secondary data made available by the Demographic and Health Survey (DHS) MEASURE. According to survey protocols and related documents of the two surveys, ethical approval was obtained from the Institutional Review Board of International Consulting Firm (ICF) of Calverton, Maryland, and the “Comité National d’Éthique” of Burkina Faso.

Consent for publication

Not applicable.

Disclaimer

The findings and conclusions in this article are those of the authors and do not necessarily represent the official position of the funders or institutions.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Supplementary Material 3

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Traoré, N., Singhal, T., Millogo, O. et al. Relative effects of climate factors and malaria control interventions on changes of parasitaemia risk in Burkina Faso from 2014 to 2017/2018. BMC Infect Dis 24, 166 (2024). https://doi.org/10.1186/s12879-024-08981-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-024-08981-2

Keywords