- Open Access
Assessing the social vulnerability to malaria in Rwanda
Malaria Journal volume 14, Article number: 2 (2015)
Since 2004, malaria interventions in Rwanda have resulted in substantial decline of malaria incidence. However, this achievement is fragile as potentials for local malaria transmissions remain. The risk of getting malaria infection is partially explained by social conditions of vulnerable populations. Since vulnerability to malaria is both influenced by social and environmental factors, its complexity cannot be measured by a single value. The aim of this paper is, therefore, to apply a composite indicator approach for assessing social vulnerability to malaria in Rwanda. This assessment informs the decision-makers in targeting malaria interventions and allocating limited resources to reduce malaria burden in Rwanda.
A literature review was used to conceptualize the social vulnerability to malaria and to select the appropriate vulnerability indicators. Indicators used in the index creation were classified into susceptibility and lack of resilience vulnerability domains. The main steps followed include selection of indicators and datasets, imputation of missing values, descriptive statistics, normalization and weighting of indicators, local sensitivity analysis and indicators aggregation. Correlation analysis helped to empirically evidence the association between the indicators and malaria incidence.
The high values of social vulnerability to malaria are found in Gicumbi, Rusizi, Nyaruguru and Gisagara, and low values in Muhanga, Nyarugenge, Kicukiro and Nyanza. The most influential susceptibility indicators to increase malaria are population change (r = 0.729), average number of persons per bedroom (r = 0.531), number of households affected by droughts and famines (r = 0.591), and area used for irrigation (r = 0.611). The bed net ownership (r = −0.398) and poor housing wall materials (0.378) are the lack of resilience indicators that significantly correlate with malaria incidence.
The developed composite index social vulnerability to malaria indicates which indicators need to be addressed and in which districts. The results from this study are salient for public health policy- and decision makers in malaria control in Rwanda and timely support the national integrated malaria initiative. Future research development should focus on spatial explicit vulnerability assessment by combining environmental and social drivers to achieve an integrated and complete assessment of vulnerability to malaria.
In Rwanda, malaria is mesoendemic in the lowlands and hypo-endemic in highlands. Malaria is a public health concern in Rwanda because the entire population of is at risk for malaria , including an estimated 2.2 million children under five years of age and 443,000 pregnant women . Nearly 63% of the country is epidemic-prone, while the remaining of the country is characterized by a stable and endemic malaria transmission. Malaria transmission occurs throughout the year with two peaks in rainy seasons. In addition to climate suitability, other factors that influence malaria transmission, including human settlements near the marshlands, internal population movement and migrations, cross-border movement of people and irrigation schemes .
In 2005, the Government of Rwanda benefited from US Presidential Malaria Initiative (PMI) to reduce malaria-related deaths in Africa. In Rwanda, PMI is expanding coverage of insecticide-treated mosquito nets (ITNs), indoor residual spraying (IRS) with insecticides, and prompt use of artemisinin-based combination therapy (ACT) for malaria treatment . These interventions resulted in substantial decline in malaria transmission. For example, after falling down between 2006 and 2008 owning to an increase in ITNs coverage, malaria increased again 2009 because of limited coverage with ITNs. The number of malaria cases declined again in 2010 following an ITN distribution campaign . This achievement is, however, fragile as potentials for local malaria transmission remain. The entire population is at risk particularly children aged under five years and pregnant women . Highland communities suffer from epidemic malaria attributed to low immunity . Despite the decreasing malaria through intensified interventions in Rwanda, the environmental changes induced by high population pressure might expose the highland populations to an increase risk of malaria and its epidemic particularly if the current interventions are not sustained .
A study in the southern Rwanda highlands revealed that asymptomatic malaria was associated with low social conditions and ineffective use of bed nets in local communities . Failure to sustain malaria control and reduction in bed nets use may result in malaria resurgence in Rwanda . For instance, an upsurge of malaria incidence in 2009 was linked to a short delay in bed net provision . This demonstrates the fragility of progress made and the need to maintain malaria control despite decreasing malaria in Rwanda .
The declining malaria provides new challenges for its elimination, including for example the resistance of malaria vectors to insecticides and livelihood activities that increase exposure to mosquito bites . With intensive use of bed nets, immunity to malaria would develop more slowly under the reduced transmission, leading to a longer period of susceptibility . A declining malaria transmission in Rwanda may therefore impair the development of immunity, which would increase the populations vulnerability to severe malaria if the control measures are not maintained . Because of population movements and migrations, malaria is imported from endemic areas to low-transmission areas. Land use changes are clustering malaria where populations share the same social and environmental factors , such as high population densities and pressure associated with land use changes and economic activities which increase mosquito breeding sites . Therefore, approaches to malaria reduction need to be aligned with these changes by adopting new strategies .
Recently, an integrated vector management strategy was adopted in Rwanda as a framework for interventions based on local ecology, malaria epidemiology and social factors . IVM strategy targets multiple vectors and different ecological and socio-economic settings. To be effective, IVM strategy should encompass environmental modifications through infrastructural development and sanitation services to regulate not only the vectors, but also the mosquito biting exposure. Additionally, IVM should improve public health and quality of life while minimizing the social disparities . The spatial assessment of social vulnerability to malaria is, therefore, a well-timed support to the IVM initiative related to public health improvement and social disparities among the populations at risk of malaria infection.
In any society, there are groups of people who have limited control over their ability to admit to illness, mobilize resources, access health facilities and services, and make decisions. This lack of personal control places them in a position in which they can be considered to be socially vulnerable to disease infection . Social vulnerability assessment in the arena of vector-borne diseases highlights the importance of social factors that make some groups or individuals more susceptible to infection and more limited in their ability to respond to illness than others . Until now, the public health approach to malaria in sub-Saharan Africa (including Rwanda) is based on Global Health Initiative strategy, which concentrates on reducing malaria burden by funding specific interventions to strengthening the health care system . Although these programmes have temporarily reduced the overall malaria infection rates countrywide, they are unlikely to be effective in producing sustainable reductions of malaria without addressing the proximate causes of malaria transmission and ultimate or efficient causes of malaria incidence in social structure, agro-ecological settings and demographic pressure of the country . In addition, if the efforts to reduce malaria are solely concentrated on health care sector and malaria disease control, this may fail to address other factors that shape the health of individuals or community, such as rural housing, food security and employment . The reduction of malaria burden requires, therefore, an integrated design of interventions that are placed within the broader social context of malaria transmission and incidence.
While the vulnerability to malaria is a multidimensional concept encompassing both environmental and social factors, its complexity cannot be measured by a single indicator value. The aim of this paper is to apply an integrated composite indicator approach for assessing the social vulnerability to malaria in Rwanda, independent of the spatial distribution of malaria based on a purely statistical approach and a composite indicator approach. Contextually, the social vulnerability to malaria encompasses a broader social context of malaria incidence in addition to environmental factors of malaria transmission. It reflects the predisposition of the populations or individuals to malaria infection and theirs ability to mitigate the risk of malaria. The results from this study indicate the most vulnerable districts in Rwanda as a combination of susceptibility to not withstand malaria infection and lack of resilience to anticipate, to cope with or to recover from malaria episodes in Rwanda.
The social vulnerability to malaria in Rwanda is timely supporting the national integrated malaria initiatives, which seek to improve the efficacy, effectiveness and sustainability of malaria control interventions through the inter-sectors collaboration, resources allocation and health infrastructure development . This assessment is also intending to provide information on vulnerability indicators that should be given priority in targeting intervention strategies and allocating limited resources in order to reduce the existing susceptibilities, strengthen community resilience, and thus reduce malaria burden in Rwanda.
Vulnerability is a well-documented concept in the fields of disaster risk reduction and climate change adaptation [22, 23]. Vulnerability to natural hazards refers to the conditions determined by physical, social, economic, and environmental factors that increase the susceptibility of a community to hazards . Füssel  stressed the need of an integrated approach for climate change vulnerability assessment which would also consider non-climatic factors. In line with an integrative approach towards vulnerability assessments, the European-funded research project MOVE (Methods for the Improvement of Vulnerability Assessment in Europe) identified key factors and dimensions that need to be addressed . The MOVE vulnerability framework considers society and environment as a coupled system between environment, hazard and society. The core elements in MOVE framework are exposure, susceptibility and lack of resilience. Vulnerability is further characterized by different scales (local, regional, national) and dimensions (environmental, ecological, social, economic, institutional, cultural) .
While the MOVE framework is oriented towards natural hazards and climate change adaptation, some common features can emerge when assessing the social vulnerability to vector-borne diseases like malaria. Recently, there has been growing recognition that social and cultural factors significantly influence the distribution of health and illness, and that issues of inequity affect how disease incidences are distributed and treated . From an international perspective, public health community involved in malaria reduction is also increasingly receptive to broader and more encompassing definitions of the malaria problem. The biomedical community is expanding its vision on malaria because existing tools of known efficacy to anticipate mosquito bites and cope with malaria infection may be of limited value because of social barriers to their effective implementation . Moreover, ignoring the social determinants of malaria is likely to allow public decision- and policy-makers to concentrate on malaria mosquitoes and not to be concerned with thorny problems of poverty and inequalities in health care facilities, access to health treatment, distribution of land and capital resources that play also a key role in malaria incidence [21, 29]. Packard presented the similar argument when highlighting that it is an unlikely possibility to eradicate malaria by only distributing mosquito nets, while sanitation is not ensured. The burden of malaria cannot be reduced by simply attacking anopheline mosquitoes and malaria parasites. Rather, removing the barriers that prevent people from achieving and maintaining health is also an imperative task that can help to combat malaria from multiple fronts . Consequently, there is a need to move away from a narrow biomedical approach which viewed malaria as a problem of mosquitoes, parasites and vector control activities toward a more integrated perspective which is tied to social conditions that increase the rates of malaria morbidity and mortality .
As pointed out by Ribera and Hausmann-Muela, malaria should not be seen as single illness episodes, but as successive sequences of interconnected events that, in a context of poverty and social inequalities, determine the course of the illness as well as the health-seeking process . Similarly, Stratton and colleagues advocated for an action paradigm that links the traditional proximal arenas of interventions with malaria’s fundamental causes by addressing the environmental, economic, and political dimensions of malaria risk. They concluded that strategies that exclusively focus on reducing exposure to mosquitoes or treating malaria may provide tangible health benefits, but their outcome will be less effective in the long run at reducing total health burden than approaches aimed at underlying causes of differential vulnerability given the high poverty rates in most malaria endemic regions . This was also supported by Jones and Williams, when advocating for a comprehensive approach for diseases control. In their discussion on social burden of malaria, they urge for taking a broader social perspective to malaria burden and advocate for a shift to integrated designs of interventions that are placed within the broader social, cultural, political and economic context .
In the assessment of vulnerability to diseases, there is a growing literature which defines the vulnerability to malaria based on biological and disease-related factors, such as age, pregnancy and immunity [32–34]. Other studies go beyond individual behaviour or characteristics, to larger-scale social and economic conditions, because some societies have been more successful in addressing malaria disease than others . In line with this integrative assessment of social vulnerability to malaria, this study is drawn on a holistic risk and vulnerability framework which was developed by Kienberger and Hagenlocher in the context of vulnerability to vector-borne diseases . Figure 1 shows the adopted and applied framework of the domains of social vulnerability to malaria and examples of relevant indicators.
Based on recent publications by Kienberger and Hagenlocher  and Hagenlocher et al. , the vulnerability to malaria is defined as the predisposition of the population to malaria burden. It is characterized by different interrelated domains, generic and biological susceptibility, and lack of resilience. The distinction between generic and biological factors pre-disposing a community to malaria burden versus those that determine its ability to anticipate responding to, to cope with or to recover from malaria burden should be relevant for decision makers [36, 38].
Vulnerability assessments often employ indicators to better reflect complexity and real-world information into a format that is relevant for decision-making. Despite their alleged weaknesses, composite indicators provide an essential tool at the science–policy interface . The use of composite indicators is an important way to communicate and monitor the vulnerability, and allows comparisons to be made across different geographical areas . By combining social and environmental indicators, geographic variations of vulnerability to malaria can be illustrated. The vulnerability map can, therefore, indicate the relative levels of vulnerability, where there is uneven capacity for preparedness and response, and where resources might be used effectively to reduce prevailing vulnerability. This provides information to decision-makers for targeting vulnerable populations while improving their resilience . While recognizing that vulnerability is a multidimensional and dynamic condition at different spatial and temporal scales, there is a need to assess how social drivers influence the vulnerability to malaria among the districts of Rwanda using a composite indicator approach.
Selection and justification of indicators
One of the key parts of this study is the selection of relevant indicators for assessing the social vulnerability to malaria. Guided by the vulnerability framework, 19 indicators have been identified. Table 1 shows vulnerability domains, indicators used and their corresponding weight from principal component analysis (PCA) and information on sources of datasets as well.
Generally, high density and population pressure in highlands results in limited land resources and increases human susceptibility to diseases . This is relevant for Rwanda highlands, where demographic pressure has significantly modified the local environment during the past decades . Swamps were reclaimed by agriculture to feed an ever-increasing population. Water requirements for irrigation have led to modifications of surface waters . Since malaria incidences are rooted in livelihood activities which interplay with the ineffective use and non-use of bed nets , irrigation-based agricultural practices may increase the susceptibility to malaria infection . These environmental changes are likely to influence malaria incidence if adequate protection measures are not implemented in nearby communities . In additionally, high population density and pressure have significantly influenced the environmental degradation and declining landholdings [47, 48], which have pushed people to settle near unsuitable sites with more exposure to mosquito bites .
Populations with little immunity may move into malaria-prone areas where they are more vulnerable . The Rwandan population increased in past decades from approximately 2.6 million in 1960 to 8.2 million in 2002 and 10.6 millions in 2012 . This was accompanied by environmental degradation and decreasing landholding size, pushing people to settle near unsuitable sites with more exposure to mosquito bites . The resettlement of non-immune people in endemic zones was accompanied by sporadic malaria epidemics .
The household size is often associated with malaria incidence . The more people are sleeping together in the same room, the higher the probability of spreading infection to a new person. Below a certain threshold number of persons sleeping together, infection rates drop below the replacement rate and disappear even without other control measures. The threshold is likely crossed when the average household size drops below four persons . Moreover, the increasing number of persons per room results in improper use of bed nets and increased probability of being infected by malaria parasites.
From biological perspective, the susceptibility to malaria reflects the efficiency with which an infective mosquito infects humans. It is a function of immunity which depends on age, pregnancy or co-infection with other diseases . Pregnant women and children under five years of age are more susceptible to severe malaria since women’s immunity is temporarily reduced during pregnancy, while the immune system of small children is not fully developed . When HIV and malaria co-infect, a severe malaria should be expected . In non-pregnant women, HIV was found to roughly double the risk of malaria by impairing the immune response and decreasing the ability to withstand malaria infection , and reducing the efficacy of drugs . A recent study in Rwanda confirmed higher malaria prevalence among HIV-positive, pregnant women . With regard to age, adults would be able to withstand malaria infection because of acquired immunity from previous exposure to mosquito bites . Famines and food shortage induced by drought also result in population movements and migration, increasing exposure to malaria in endemic lowlands, thus rising the malaria incidence in returning population to the highlands . Recurrent drought in Rwanda caused crop failures and food shortages, then threatening the most vulnerable populations with malnutrition and famine, especially in the Eastern Province of Rwanda .
The lack of resilience relates to the capacity to anticipate mosquito-biting exposure and to limited access to health infrastructure and means to recover from malaria episodes. This capacity may be influenced by the level of education, awareness about malaria transmission and prevention, and access to protection measures. The coping capacity relates to access to health care services and to adequate and effective treatment . In this study, this encompasses protection measures, housing conditions, education level and improved knowledge about malaria which lead to better use of malaria interventions . The use of bed nets is the most accepted protection method against malaria. For bed nets to be effective, coverage must be high, bed nets should be retreated promptly and individuals should properly deploy their bed nets each night. The more households own and use bed nets, the greater the benefit to neighbouring households without bed nets . Indoor residual spraying (IRS) application on houses’ walls and roofs and on domestic animal shelters kills the adult mosquitoes that rest on these surfaces. It reduces the longevity and density of mosquitoes so that they can no longer transmit malaria parasites . Information on IRS was not however factored in this study because the related data was not available for the entire study area. Once this information is available for the entire country, this could additionally reduce uncertainties in vulnerability assessment.
Poor housing quality provides less protection against mosquitoes . Houses with more malaria-infected children are more likely to have mud walls, open eaves and absent ceilings than those with uninfected children . Despite the significant reductions in malaria transmission by high coverage of bed net use in Rwanda, high numbers of host-seeking malaria vectors in rural areas may be found indoors due to poor housing quality. Education levels may influence malaria incidence since it affects the knowledge about malaria prevention and control. Over recent years there has been emphasis on that level of education and improved knowledge about malaria leads to better use of malaria interventions . Nevertheless, association between malaria and education may be due to its role as a proxy for poverty. Media communication plays a strong role in malaria control to ensure that bed nets are used appropriately . Mobile phone technology is also an efficient method for rapidly detecting malaria patients and reduce malaria deaths in remote rural areas .
In terms of coping and recovering capacity, the vulnerability to malaria may be influenced by access to treatment and prompt access to effective malaria treatment is central to the success of malaria control worldwide. The Roll Back Malaria partnership has set for 2010 a target of ensuring that 80 percent of those suffering from malaria have prompt access to, and are able to correctly use, affordable and appropriate treatment within 24 hours of symptoms onset . A strategy to provide such access should take into account poor populations in malaria-endemic zones who are particularly inadequately served by the health system . In such areas, convergence of malaria prevalence and poor health care infrastructure can result in high malaria incidence rates. Moreover, primary care facilities are increasingly becoming the focal point for distribution of intervention strategies, but physical access to these health facilities may limit the extent to which communities can be reached .
Construction of a composite index of social vulnerability to malaria
To construct a composite index of social vulnerability, a methodology developed by Nardo et al.  was adopted. The steps followed were: the selection of relevant indicators, identification of appropriate datasets, imputation of missing data, and descriptive statistics for data exploration (distribution and multicollinearity among indicators, outliers’ detection), data transformation and normalization. The PCA was performed to weight indicators for further aggregation. Prior to aggregation and final vulnerability index creation, the local sensitivity analysis was applied to evaluate the influence of indicators on final index. Then Spearman correlation the regression analysis was applied to provide evidence of association between the used indicator and malaria at district level by highlighting the most sensitive indicators. A final composite index of vulnerability was then computed and visualized for each district. Figure 2 illustrates the main steps followed in constructing a composite index of social vulnerability to malaria.
Descriptive statistics were used for each indicator to evaluate the degree of missing data and potential outliers. Skewness with absolute values above 0.2 were taken as highest threshold for multicollinearity to avoid double counting among indicators . Indicators with kurtosis values greater than 3.5 were treated for asymmetric distribution . Outliers were removed by changing outliers’ values to next highest/lowest non-outlier number . Afterwards, outlier analysis was re-run to check if all outliers were eliminated .
Since indicators were at different measurement scales, standardization was required before their aggregation . Sometimes high ‘raw’ indicator data values mean low vulnerability, so therefore theses indicators have been transformed to have high values for high vulnerability. The minimum-maximum transformation method was used by accounting for the direction of indicators  and using the flowing formula:
Where Vi is equal to the standardized indicator i; Xi represents the indicator value before its transformation; Xi, min is the minimum score of indicator i before its transformation; and Xi, max as maximum score of indicator i before its transformation. All indicator values were transformed into a relative score ranging from 0 to 1, where higher values imply high vulnerability . A positive sign implies that high indicator values increase the vulnerability (+), while low values decrease the vulnerability and vice versa.
Prior to PCA, the correlation coefficient matrix of two vulnerability domain indicators was scanned to check for the values greater than 0.92 as an indication of collinearity as an indication of collinearity . After scanning the correlation coefficient matrix, no value greater than this threshold was found. Since many correlation coefficient values were above 0.300, the PCA would yield the acceptable results . The Kaiser-Meyer-Olkin (KMO) criteria and Bartlett’s test were also performed to examine the data suitability for PCA. The general rule is that the overall KMO should be 0.60 or higher to proceed with PCA. The high values for overall KMO for both susceptibility and lack of resilience indicators were 0.655 and 0.615, respectively, meaning that they are acceptable for PCA.
PCA for weighting indicators
A PCA was used to assign the weights to individual indicators based on common variance explained. Using a varimax orthogonal rotation, components with eigenvalues larger than one; which contribute individually to overall variance by more than 10%; and cumulatively to more than 60% were chosen . Table 2 illustrates the two extracted components explain 64.54% of total variance for lack of resilience indicators.
As earlier mentioned, components with eigenvalues larger than one, individually contributing to overall variance by more than 10%, and cumulatively by more than 60%, were chosen to weight the individual of indicators. Table 3 shows the highest indicator scores in each component and their weights for each single indicator are highlighted in bold.
The lack of resilience indicators that have the highest squared loading factors in the first principal component are poor housing roof materials, poor housing wall materials, low literacy rate and nurse ratio to populations. The second principal component is heavily loaded by the households without mobile phone, the households without radio, households without bed nets and number of health facilities.
With regard to the susceptibility indicators, three extracted components were able to explain 73.04% of the total variance as shown in Table 4.
The susceptibility indicators with the highest indicator scores in each principal component are highlighted in bold as illustrated in Table 5.
The susceptibility indicators that have the highest squared loadings in the first principal component are the number of poor populations, number of arriving populations, HIV prevalence, and population density. Apart from HIV prevalence, most of the highest scored indicators in first principal component are related to generic susceptibility sub-domain. With the exception of the average number of persons per bedroom, the majority of the biological susceptibility sub-domain (children under five years of age, women of child-bearing age, households affected by droughts and famines and elder population above 65 years) dominate the second principal component. The land area used for irrigation and the population change between 2002 and 2012 have an excessive influence in the last principal component.
Approaches for indicator weighting and aggregation have been subject to debate, but so far no approach is without limitations . Some approaches employ equal weights which ensure transparency and are straightforward, but have been criticized for assignment of implicit equal weights . A study by Hagenlocher et al.  for modeling the vulnerability to disease based on administrative boundaries concluded expert-based and purely statistical-based modeling approaches reveal similar outputs, indicating that in the absence of local expertise, statistical approaches could be used. The PCA was therefore used since it has the advantage of determining weights which explains the largest variation in original indicators . The weighted sum algorithm was chosen as the most widespread aggregation method. The normalized indicators were first aggregated according to their respective domains. The normalized weighted indicators for each domain were summed using the following formula:
SUS and LoR represent the value of vulnerability for susceptibility and lack of resilience domains respectively; W refers to the weight of a single indicator in each domain; and I is the normalized value of the indicator. The final composite index was calculated by aggregating the two domains and taking into account the number of indicators in each domain so that the domains grouping the larger number of indicators will have higher weight as follow:
Where n represents the number of indicators for a given domain; d refers to the value of each vulnerability domain while N is equal to the total number of indicators. From nineteen indicators that have been identified from literature, eleven indicators were assigned to susceptibility domain, and eight indicators to the lack of resilience domain. For easy visualization of the results, the final index values were normalized within a new range from zero to one, where zero reflects a very low and one a very high social vulnerability to malaria. The higher the values of the vulnerability index, the more the district is vulnerable.
Sensitivity analysis evaluates the contribution of individual sources of the uncertainty to the output. For this study, spearman correlation analysis was first used to validate the appropriateness of using indicators using International Business Machines Corporation SPSS statistics 22.0 for windows software. The general aim was to provide empirical evidence of association between indicators and malaria incidence by highlighting the most appropriate to influence malaria incidence at district level. In addition to correlation analysis, local sensitivity analysis helped to assess the influence of the input vulnerability indicators on the vulnerability index. This was achieved by targeting one construction stage at a time, while all other stages are held constant . Consequently, the use of box plots helped to assess the influence of the input vulnerability indicators by discarding one input at the time while keeping all other setting (normalization, weighting and aggregation) equal . This resulted in a series of alternative vulnerability indices. For each district, the alternative index was compared with the reference vulnerability index that takes into account the susceptibility and lack of resilience indicators. The results are displayed in the box plots showing the interquartile range, the minimum and maximum values . The larger the interquartile range, the higher is the influence of the respective input indicator .
Visualization of the results
Because public health decision-makers need the information on the most vulnerable districts and the social drivers of the vulnerability to malaria, a cartographic visualization method was adopted to translate the vulnerability index into a geographic map. ESRI ArcGIS10.2 software was used to display and map the final index of social vulnerability to malaria for each district. To ease interpretation of the assessment results, the final vulnerability index values were normalized within the zero to one range, where zero reflects a very low and one a very high social vulnerability to malaria. The values of the vulnerability index reflect the relative levels of social vulnerability among the district of Rwanda, which means that an index value of 0.00 does not imply the absence of vulnerability. The developed composite index shows the most vulnerable districts to integrated social indicators in terms of susceptibility to not withstand malaria and lack of resilience to anticipate, to cope with or to recover from malaria infection in Rwanda.
Vulnerability to malaria and underlying factors
As shown in Figure 3, high values of vulnerability index to malaria are found in Gicumbi (1.00), Rusizi (0.81), Nyaruguru (0.79), Gisagara (0.71), and Burera (0.67), and districts, and low values in Muhanga (0.00), Nyarugenge (0.10), Kicukiro (0.13), and Nyanza (0.17).
It is very important know that the vulnerability index values reflect the relative levels of vulnerability, which means that the value of 0.00 does not mean the absence of vulnerability to malaria within the district.
The bar charts for the most vulnerable districts in Figure 4 indicate which districts are less resilient, the relative share and contribution of the underlying factors of vulnerability. For lack of resilience domain, the poor housing wall and roof materials and low literacy rate are the main factors that are hampering the capacity of Gicumbi, Ruzizi, Nyaruguru, Gisagara and Nyamagabe districts to anticipate malaria mosquito-biting exposure (Figure 3). The vulnerability of Gicumbi, Nyaruguru and Nyamagabe districts is also exacerbated by the low rate of bed nets ownership low immune highland populations. Because the most vulnerable districts are mostly located in remote rural areas near the borders, access to communication through mobile technology is an impeding factor. With regard to capacity to cope with or to recover from malaria infection, the limited number of health facilities and insufficient medical personnel call for more improved interventions in Gicumbi, Rusizi and Gisagara districts. In these highland districts prone to epidemics malaria induced by climate change and variability, malaria incidence may combine with limited health infrastructure and poverty to result in high morbidity and mortality.
With regard to susceptibility domain, the district of Gicumbi is the most vulnerable because of poverty, high number children under five years of age, high number of child-bearing age women, and high number of elder populations (Figure 5).
The vulnerability of Gicumbi is also exacerbated by high number of arriving and mobile populations, particularly the refugees from North Kivu in Democratic Republic of Congo who are hosted in Gihembe camp. The vulnerability of Rusizi district is largely associated with a high number of child-bearing age women, high number of children under five years, high population density and high rate of HIV prevalence compared with other districts . The number of arriving populations was also noticed in Rusizi district. The susceptibility of Nyaruguru district is mostly explained by poverty, high number of persons per bedroom, high number households affected by famines, and high number of elder populations. The most influencing indicators in Gisagara district are poverty, high number of child-bearing age women, the irrigated livelihoods that increase the exposure of local population to mosquitoes’ bites, and high population density. The indicators that make Nyamagabe district more susceptible are mainly the poverty and demographic pressure.
Influence of input indicators on vulnerability index
The composition of susceptibility and lack of resilience was evaluated in terms of their influence on final composite indicator of social vulnerability to malaria. The influence of each of the single indicator was calculated by discarding one input at the time while keeping all other all settings remain constant. For susceptibility domain, the results show that the average number of persons per bedroom, elder populations above 65 years, population changes and population density have less impact on final vulnerability index. Conversely, the number of poor populations, number of arriving populations, number of children under five years and women in child-bearing age have an excessive influence on the final composite vulnerability index. In Figure 6, the box plots were used to show the influence of the single indicator on the composite vulnerability index.
For lack of resilience domain, the number of households without mobile phones and the number of households without radio are the less influential indicators. On the other hand, poor housing wall materials, poor housing roof material, number of households without bed nets have a marked impact of vulnerability index.
Association between the used indicators and malaria incidence
The regression analysis was used to validate the relevance of indicators by highlighting the most sensitive to influence malaria occurrence. Table 6 shows the spearman correlation coefficients between the used indicators and malaria positivity rate for the year 2010.
The susceptibility indicators that are scientifically correlated with malaria incidence are highlighted in bold. The population change 2002–2012, land area used for irrigation, households affected by drought and famine and average number of persons per bedroom have the high correlation coefficient values (r) raining from 0.729 to 0.531. The susceptibility indicators that have low values of correlation coefficient are the number of poor populations, number of arriving populations and, HIV prevalence and (r value range from 0.018 to 0.130).
This association was expected since the increase in populations augment their contact with an infective mosquito feeding on humans in high transmission areas . Owing to population increase in Rwanda, people moved to the lowlands and –malaria endemic districts such as Nyagatare, Gatsibo and Kayonza  with an increasing vulnerability of non-immune migrants. A moderate positive association between the number of persons per bedroom and malaria positivity rate was also found. When household size increases, bed nets are not properly used. Evidence has shown that people sharing a bed net with more than five people are more likely to have malaria than those sharing a bed net with up to two people . The districts with the highest average number of persons per bedroom are Nyagatare, Ngoma and Bugesera. This exacerbates the vulnerability of these most endemic-malaria districts of Rwanda.
A strong relationship between the number of households affected by drought and malaria is not surprising. While prolonged droughts reduce malaria transmission or turn rivers into strings of pools, preferred mosquito breeding sites after short rain, they may also reduce the food security, increase malnutrition and people susceptibility to malaria . The drought-prone districts in Rwanda are Bugesera, Kirehe, Ngoma, Rwamagana, Kayonza, Gatsibo and Nyagatare, Nyanza, Gisagara, Huye and Rusizi . This relationship between malaria and droughts is justified as the drought–prone districts are also malaria-prone areas.
The land area used for irrigation is also strongly correlated with malaria due to the creation of vector breeding sites by irrigation projects . Currently, most of the wetlands in Rwanda are being reclaimed for irrigated crops . A moderate correlation was also found between the number of children under five years and elder populations.
The lack of resilience indicators that have the high correlation coefficient values (r) are poor housing wall materials and the bed nets ownership with the correlation coefficient values of 0.378 and −0.398 respectively. Once again, the positive association between housing conditions and malaria incidence at district level in Rwanda highlights the need for effective malaria intervention that also target the housing improvement in rural areas in addition to bed net provision. Malaria decreases as the number of households with bed nets increases because intense malaria interventions with mass distribution of bed nets led to substantial malaria decline in Rwanda since 2006.
The present work was conducted to assess the social vulnerability to malaria in Rwanda independent of the spatial distribution of the disease based on purely statistical approach and a composite indicator framework. The results from this study are salient for public health policy- and decision makers in malaria control in Rwanda. This spatial assessment is a well-timed support to the national integrated malaria initiative which to seeks to improve the efficacy, effectiveness and sustainability of malaria control interventions through advocacy, social mobilization and inter-sectors collaboration to optimize the allocation of limited resources and health infrastructure at national scale level.
By displaying the lack of resilience versus susceptibility to malaria, the adopted approach provides useful information for decision-makers and a way of communicating the complex interactions between relevant factors of vulnerability to malaria. Squeezing the complex system of social conditions into a single vulnerability index, the developed approach yields a powerful comparative assessment tool capable of capturing societal conditions in a given district that drive people’s vulnerability to malaria infection. This has important policy implications as a successful reduction of malaria burden requires combining the best set of strategies that address the most important vulnerability factors in the most vulnerable districts.
By the fact that most of the vulnerable districts are located in the highlands where unstable malaria transmission would be limited by low temperatures, with the exception of Nyagatare and Bugesera this has an implication for malaria control program in Rwanda. In these low malaria transmission settings, malaria burden is difficult to assess, its impact is significant, and the cost-effectiveness of interventions to predict and respond is doubtful . As an example, malaria epidemics have been reported in Rwanda highlands since 1940s. . Later in the 1980s a steady malaria increase was attributable to low immune population movements from Nyarutovu highlands within Nemba Hospital catchment’s area to the lowlands in the East . Likewise, malaria increased during 1983 to 1987 at Gikonko health centre (Gisagara District), as a result of climate warming near altitude limits of malaria transmission . Unexpected rainfall also resulted in change of malaria patterns in northern Rwanda . In Byumba district hospital catchment’s area in Gicumbi district, a sharp increase of malaria among pregnant women and children was linked to rainfall anomalies in 1998 . It was therefore hypothesized that higher monthly rainfall increased malaria incidence in Byumba highlands, averaging at 2,300 m.
According to Rwanda Ministry of Health, ten border districts have the highest malaria prevalence . The remoteness of rural communities near the borders makes it difficult to provide adequate health services. When border control is inadequate, population’s vulnerability to malaria increases . Near the border areas, malaria infection may not be linked to location at which population came into contact with infective bite(s) of Anopheles vectors . Although long-term migrations are limited in Rwanda, temporary migrations with a significant heterogeneity were pointed out by Blumenstock , where cross-border movements are responsible for malaria transmission . Consequently, providing adequate health care facilities and appropriate and effective treatment to mobile communities near the borders would be an efficient malaria intervention. Malaria elimination needs to identify those migrant streams in the most vulnerable districts with potential to transport malaria and to target interventions accordingly to prevent potential epidemics malaria .
Likewise, conflicts are likely to precipitate the movement of refugees across borders and disruption of local infrastructure . They amplify the vulnerability to malaria owing to breakdown of health systems, mass population displacements, and resettlement of refugees in camps in malaria prone areas. Figure 7 shows the location of refugee camps in relation to the levels of social vulnerability to malaria among the district of Rwanda.
In emergency situations, malaria flourishes as a result of breakdown health facilities, displacement of non-immune populations to malaria-prone areas, malnutrition, inappropriate access to treatments, and lack of appropriate shelter . Since 1996, conflicts have resulted in displacement of refugees seeking sanctuary in Rwanda.
Based on recent statistics of the Ministry of Disaster Management and Refugee Affairs in Rwanda, the refugee camps population was estimated to 73,786 until the end of April 2014. The largest camps are Kigeme (18,430), Kiziba (16,461), Gihembe (14,707), Nyabiheke (14,213) and Mugombwa (6,688). As livelihood opportunities are very limited for refugees who generally live in poor housing , refugees are more vulnerable to malaria infection. Particular attention should be paid to Mugombwa refugee camp, which is located in the most malaria-endemic area within the most vulnerable district of Gisagara.
Additionally, it is apparent that high population numbers for the most vulnerable age groups (children under five years, women in child-bearing age and elder population) largely explains the susceptibility to malaria in highlands districts. Higher population density and pressure in the highland districts result in stressed productive land as a result of land fragmentation and declining landholding size, pushing people to settle in unsuitable locations or to migrate to malaria-endemic areas with an increasing human exposure to malaria. In past five years, malaria-endemic districts attracted internal migrants are Kayonza, Gatsibo, Rwamagana, Kirehe, and Bugesera in Eastern Province, but also Gasabo and Kicukiro districts within Kigali . The population increase and high average household size may act as impeding factor for effective use of bed nets in Nyagatare and in Nyaruguru districts . This research finding implies that in regions with large households or large populations sharing sleeping room such as refugee camps or worker migration camps near the rice farming areas, the elimination of malaria will require segmenting sleeping quarters into smaller units, such as with mosquito nets.
The low rate of bed nets ownership exacerbates the exiting vulnerability of Gicumbi and Nyaruguru districts. These districts are characterized by low malaria endemicity justifying the low level of malaria interventions with IRS and LLITNs. However, preventive measures targeted to potential hotspots of malaria transmission in highland districts could effectively help in for preventing epidemics and highland malaria. In Nyaruguru, Gisagara and Rusizi District, poverty and poor housing conditions suggest that housing characteristics should be considered in malaria interventions to prevent mosquitoes entering houses and reduce mosquito-biting exposures. This result is in strong agreement with the fourth population and housing census report in Rwanda  which revealed that the high percentage of households with poor housing walls in Southern Province is a result of poverty in rural populations. The population increase coupled with land use changes through the extensive marshlands reclamation for irrigation in Gisagara and Nyagatare districts should also attract the attention of public health planners .
As the area allocated to irrigation is strongly associated with malaria incidence, this has an implication for malaria control program, to take into account the livelihood activities which interplay with the ineffective use and non-use of bed nets to increase the vulnerability to malaria. This is more important in the country like Rwanda where marshlands irrigation as adaptation to climate variability and food shortage has favoured the development of mosquitoes and malaria expansion in highlands . In responding to food shortage, fish ponds in valley bottoms shaped the creation of vector breeding sites . Since 2008, marshland irrigation is being promoted for sustainable food security . Currently, most wetlands in Rwanda are reclaimed for growing rice  and rice production increased from 11,949 tons in 2000 to 72,000 tons in 2009 . Rice cropping is a promising solution to food insecurity but also increases malaria in local communities  as reported by farmers in Cyabayaga and Rugeramigozi wetlands . Association between malaria and marshlands cleaning is also confirmed by the coincidence of spatial distribution of malaria parasite prevalence rate and area affected by irrigation per district as illustrated in Figure 8.
Ijumba and Lindsay  argued that irrigation increases an unstable malaria transmission, where people have little immunity. The vulnerability of irrigation communities may however, be low in some areas owing to better social conditions and farmers’ awareness of bed net use . Planners should consider this opportunity to improve health care facilities when planning marshland irrigation in unstable malaria transmission zones. Integrated malaria control and health impact assessments are recommended near irrigation schemes.
The proposed indicator framework is however without limitations and comes along with some challenges. First, malaria vulnerability assessment based on administrative boundaries may be incomplete for making decisions in complex situations. The composite index based on district administrative boundaries in Rwanda may not be always appropriate for explicitly displaying spatial vulnerability to malaria because boundaries rarely correspond to variations in prevailing social vulnerability at country level . Moreover, the scale is often determined by what is conveniently collected rather than by what is appropriate . Beside, the transmission of malaria does not respect the artificial boundaries which are defined for administrative and political purposes. The PCA weighting scheme does not consider how indicator influence might change across the district. Since PCA analysis is a global statistical analysis, factors and components loadings will not vary at the local scale, and therefore extracted components using the standard PCA do not depend on location. Thus, a strategy for integrating indicators for assessing the social vulnerability to malaria in explicit manner  and to examine how the influence of indicators changes across the districts through statistical methods that correct for spatial processes may be required. Secondly, additional indicators such as extent and coverage of IRS campaigns were not available and thus not integrated into the vulnerability analysis because the related data was not available for the entire study area. Once this information is available for the entire country, this could additionally reduce uncertainties in vulnerability. Additionally, some indicators such as social networks and behavior change are difficult to quantitatively measure in malaria vulnerability assessment and therefore not considered in the in this paper. Finally, future research development should focus on an integrated and spatial explicit vulnerability assessment by combining environmental and social drivers is important to achieve an integrated and complete assessment of vulnerability to malaria in order to target malaria interventions that are responsive to the needs of the most vulnerable people in Rwanda.
This paper applied a composite indicator approach for assessing social vulnerability to malaria among the districts of Rwanda. It is drawn on published works and simplifies the complex information from multisource indicators of vulnerability to malaria into a format that is relevant for decision-making. It shows the most vulnerable districts to integrated social indicators in terms of susceptibility to not withstand malaria and lack of resilience to anticipate, to cope with or to recover from malaria infection. By decomposing vulnerability into its underlying factors, it indicates which factors need to be addressed in each district. The developed composite indicator framework supports the prioritization of appropriate interventions in Gicumbi, Burera, Nyaruguru, Nyagatare, Gisagara, Bugesera, Rusizi, and Burera districts. Being located in the highlands, the prevailing vulnerability in these districts may be exacerbated by cross-border migrations where malaria can be imported from outside.
The health of people in most vulnerable district will not improve unless poverty and expanding inequality are reduced and this includes the effort to control malaria on a large scale. But nothing can be accomplished without positioning the problem in social and cultural contexts of Rwanda. Increasing the community resilience in terms of bed net provision, housing improvement, poverty reduction and access to health care facilities and treatment can be seen as a promising approach for policy makers to be proactive towards malaria in the most vulnerable districts. Besides, improved land-use planning and environmental management can reduce community susceptibility to malaria in Rwandan highlands. Policies enabling activities towards border-crossing populations and enhancement of vector control is needed. As the indicators were aggregated at district level, a study using disaggregated data at household level may be needed. Finally, the approach used provides a comparative assessment and generalizes the relative levels of vulnerability at among the district of Rwanda, but does not provide useful information about what areas within the districts are most vulnerable to malaria infection. A lack of spatial details on prevailing vulnerability within the same district can result in the implementation of uniform interventions that do not necessarily translate to uniform vulnerability reduction inside the district.
President’s Malaria Initiative: Malaria Operational Plan (MOP). Kigali, Rwanda: Rwanda FY 2014; 2014.
National Institute of Statistics of Rwanda: Fourth Population and Housing Census. Thematic Report: Characteristics of households and housing. Ministry of Finance and Economic Planning; 2012.
President’s Malaria Initiative: Malaria Operational Plan (MOP). Kigali, Rwanda: Rwanda FY 2013; 2013.
Karema C, Aregawi MW, Rukundo A, Kabayiza A, Mulindahabi M, Fall IS, et al.: Trends in malaria cases, hospital admissions and deaths following scale-up of anti-malarial interventions, 2000–2010, Rwanda. Malar J 2012, 11:236. 10.1186/1475-2875-11-236
Government of Rwanda, President’s Malaria Initiative: Malaria Operational Plan (MOP). Kigali, Rwanda: Rwanda FY 2013; 2012.
Government of Rwanda: Rwanda Service Provision Assessment Survey 2007. Kigali, Rwanda: National Institute of Statistics Ministry of Finance and Economic Planning; 2008.
Himeidan YE, Kweka EJ: Malaria in East African highlands during the past 30 years: impact of environmental changes. Front Physiol 2012, 3:315.
Gahutu JB, Steininger C, Shyirambere C, Zeile I, Cwinya-Ay N, Danquah I, et al.: Prevalence and risk factors of malaria among children in southern highland Rwanda. Malar J 2011, 10:134. 10.1186/1475-2875-10-134
WHO: World Malaria Report. Geneva, Switzerland: World Health Organization; 2010.
WHO: World Malaria Report. Geneva, Switzerland: Global Malaria Programme; World Health Organization; 2011.
Mboera LEG, Mazigo HD, Rumisha SF, Kramer RA: Towards malaria elimination and its implication for vector control, disease management and livelihoods in Tanzania. Malaria World Journal 2013, 4:19.
Lengeler C: Insecticide-treated nets for malaria control: real gains. Bull World Health Organ 2004, 82:2.
Githeko AK, Ototo EN, Guiyun Y: Progress towards understanding the ecology and epidemiology of malaria in the western Kenya highlands: opportunities and challenges for control under climate change risk. Acta Trop 2012, 121:19–25. 10.1016/j.actatropica.2011.10.002
Rulisa S, Kateera F, Bizimana JP, Agaba S, Dukuzumuremyi J, Baas L, et al.: Malaria prevalence, spatial clustering and risk factors in a low endemic area of eastern Rwanda: a cross sectional study. PLoS One 2013, 8:e69443. 10.1371/journal.pone.0069443
Thaxton M: Integrating Population, Health, and Environment in Rwanda. Rwanda February: Policy Brief of the Population Reference Bureau; 2009.
Cotter C, Sturrock HJW, Hsiang MS, Liu J, Phillips AA, Hwang J, et al.: The changing epidemiology of malaria elimination: new strategies for new challenges. Lancet 2013, 382:900–11. 10.1016/S0140-6736(13)60310-4
Lizzi KM, Qualls WA, Brown SC, Beier JC: Expanding integrated vector management to promote healthy environments. Trends Parasitol 2014, 30:394–400. 10.1016/j.pt.2014.06.001
Jones CO, Williams HA: The social burden of malaria: what are we measuring? Am J Trop Med Hyg 2004, 71:156–61.
Parker RG: Empowerment, community mobilization and social change in the face of HIV/AIDS. AIDS 1996, 10:S27–31.
Packard RM: The Making of a Tropical Disease: A Short History of Malaria. Baltimore: John Hopkins University Press; 2007.
Morgan L: Community participation in health: perpetual allures, persistent challenge. Health Policy Plan 2001, 16:221–30. 10.1093/heapol/16.3.221
Birkmann J: Measuring Vulnerability to Promote Disaster-Resilient Societies: Conceptual Frameworks and Definitions. In Measuring Vulnerability to Natural Hazards. Edited by: Birkmann J. Tokyo-New York-Paris: United Nations University Press; 2006.
Birkmann J, Cardona OA, Carreno L, Barbat A, Pelling M, Schneiderbauer S, et al.: Framing vulnerability, risk and societal responses: the MOVE framework. Nat Hazards 2013, 67:193–211. 10.1007/s11069-013-0558-5
UNISDR: Terminology on Disaster Risk Reduction. Geneva, Switzerland: United Nations International Strategy for Disaster Reduction; 2009.
Füssel HM: Adaptation planning for climate change: concepts, assessment approaches, and key lessons. Sustain Sci 2007, 2:265–75. 10.1007/s11625-007-0032-y
Field CB, Barros V, Stocker TF, Qin D, Dokken DJ, Ebi KL, et al.: Managing the risks of extreme events and disasters to advance climate change adaptation. Cambridge Melbourne, Madrid, Cape Town, Singapore, São Paulo, Delhi, Tokyo, Mexico City, New York: Special Report of the Intergovernmental Panel on Climate Change; Cambridge University Press; 2012.
Heggenhougen HK: More than just “interesting!” Anthropology, health and human rights. Soc Sci Med 2000, 50:1171–5. 10.1016/S0277-9536(99)00416-5
Packard RM, Brown PJ: Rethinking health, development, and malaria: historicizing a cultural model in international health. Med Anthropol 2010, 17:181–94.
Brown PJ: Culture and the Global Resurgence of Malaria. In The Anthropology of Infectious Disease: International Health Perspectives. Edited by: Inhorn MC, Brown PJ. Amsterdam: Gordon and Breach Science Publishers; 1997:119–41.
Ribera JM, Hausmann-Muela S: The straw that breaks the camel’s back: redirecting health-seeking behavior studies on malaria and vulnerability. Med Anthropol Q 2011, 25:103–21. 10.1111/j.1548-1387.2010.01139.x
Stratton L, O’Neill MS, Kruk ME, Bell ML: The persistent problem of malaria: addressing the fundamental causes of a global killer. Soc Sci Med 2008, 67:854–62. 10.1016/j.socscimed.2008.05.013
Protopopoff N, Bortel WV, Speybroeck N, Van Geertruyden J-P, Baza D, D’Alessandro U, et al.: Ranking malaria risk factors to guide malaria control efforts in African Highlands. PLoS One 2009, 4:e8022. 10.1371/journal.pone.0008022
Bates I, Fenton C, Gruber J, Lalloo D, Medina LA, Squire SB, et al.: Vulnerability to malaria, tuberculosis, and HIV/AIDS infection and disease. Part I: determinants operating at individual and household level. Lancet Infect Dis 2004, 4:267–77. 10.1016/S1473-3099(04)01002-3
Sutherst RW: Global change and human vulnerability to vector-borne diseases. Clin Microbiol Rev 2004, 17:136–73. 10.1128/CMR.17.1.136-173.2004
Bates I, Fenton C, Gruber J, Lalloo D, Medina LA, Squire SB, et al.: Vulnerability to malaria, tuberculosis, and HIV/AIDS infection and disease. Part II: determinants operating at environmental and institutional level. Lancet Infect Dis 2004, 4:368–75. 10.1016/S1473-3099(04)01047-3
Kienberger S, Hagenlocher M: Spatial-explicit modeling of social vulnerability to malaria in East Africa. Int J Health Geogr 2014, 13:29. 10.1186/1476-072X-13-29
Hagenlocher M, Delmelle E, Casas I, Kienberger S: Assessing socioeconomic vulnerability to dengue fever in Cali, Colombia: statistical vs expert-based modeling. Int J Health Geogr 2013, 12:36. 10.1186/1476-072X-12-36
Hagenlocher M, Kienberger S, Lang S, Blaschke T: Implications of spatial scales and reporting units for the spatial modelling of vulnerability to vector-borne diseases. In Jekel T, Car A, Strobl J, Griesebner G, editors. GI_Forum 2014: Geospatial Innovation for Society. Salzburg, Austria: Heidelberg: Wichmann; 2014:197–206.
Ravallion M: Troubling Tradeoffs in the Human Development Index; Policy Research Working Paper 5484. Washington DC: World Bank; 2010.
Brooks N, Adger WN, Kelly PM: The determinants of vulnerability and adaptive capacity at the national level and the implications for adaptation. Glob Environ Chang 2005, 15:151–63. 10.1016/j.gloenvcha.2004.12.006
Borderon M: Why Here and not There? Developing a Spatial Risk Model for Malaria in Dakar, Senegal. In From Social Vulnerability to Resilience: Measuring Progress Toward Disaster Risk Reduction; Publication Series of UNU-EHS No 17/2012, vol. 108–120. Edited by: Cutter S, Corendea C. Bonn, Germany: United Nations University Institute for Environment and Human Security (UNU-EHS); 2013.
Lindsay SW, Martens WJM: Malaria in the African highlands: past, present and future. Bull World Health Organ 1998, 76:33–45.
Malesu MM, Oduor AR, Chrogony K, Nyolei D, Gachene CKK, Biamah EK, et al.: Rwanda Irrigation Master Plan. Nairobi, Kenya: The Government of Rwanda, Ministry of agriculture and animal resources, Ebony Company Limited and World Agro-forestry Centre (ICRAF); 2010.
Dunn CE, Le Mare A, Makungu C: Malaria risk behaviours, socio-cultural practices and rural livelihoods in southern Tanzania: implications for bed net usage. Soc Sci Med 2011, 72:408–17. 10.1016/j.socscimed.2010.11.009
Mboera LEG, Mukabana WR, Njunwa KJ, Kabbale FG: Integrated research partnerships for malaria control through an ecohealth approach in East Africa: Kenya. Tanzania and Uganda Projects. Final Report: Rwanda; 2014.
Ijumba JN, Lindsay SW: Impact of irrigation on malaria in Africa: paddies paradox. Med Vet Entomol 2001, 15:1–11. 10.1046/j.1365-2915.2001.00279.x
Clay DC: Fighting an Uphill Battle: Demographic Pressure, the Structure of Landholding, and Land Degradation in Rwanda; ASPAP/DAI Rapport No 192U.S. Kigali, Rwanda: Agency for International Development and the Rwanda Ministry of Agriculture and Livestock; ASPAP/DAI Rapport No 192; 1992.
Clay DC, Johnson NE: Size of farm or size of family: which comes first? Popul Stud 1992, 46:491–505. 10.1080/0032472031000146476
Deressa W, Ali A, Berhane Y: Review of the interplay between population dynamics and malaria transmission in Ethiopia. Eth J Health and Dev 2006, 20:1–8.
Government of Rwanda: The 4th General Census of Population and Housing in Rwanda: Provisional Results. Kigali, Rwanda: Ministry of Finance and Economic Planning; 2012.
WHO: Prevention and Control of Malaria Epidemics: Tutor’s Guide. Geneva: World Health Organization; 2003.
Huldén L: Household size explains successful malaria eradication. Malar J 2010, 9:2. 10.1186/1475-2875-9-2
Huldén L, McKitrick R, Huldén L: Average household size and the eradication of malaria. J R Stat Soc 2013, 177:725–42.
WHO: Integrating Poverty and Gender into Health Programmes: A Sourcebook for Health Professionals (Module on Malaria). Geneva, Switzerland: World Health Organization; 2006.
Saracino A, Nacarapa EA, da Costa Massinga ÉA, Martinelli D, Scacchetti M, de Oliveira C, et al.: Prevalence and clinical features of HIV and malaria co-infection in hospitalized adults in Beira. Mozambique Malar J 2012, 11:241. 10.1186/1475-2875-11-241
Cohen C, Karstaedt A, Frean J, Thomas J, Govender N, Prentice E, et al.: Increased prevalence of severe malaria in HIV-infected adults in South Africa. Clin Infect Dis 2005, 41:1631–37. 10.1086/498023
González R, Ataíde R, Naniche D, Menéndez C, Mayor A: HIV and malaria interactions. Where do we stand? Expert Rev Anti-Infect Ther 2012, 10:153–65. 10.1586/eri.11.167
Ivan E, Crowther NJ, Rucogoza AT, Osuwat LO, Munyazesa E, Mutimura E, et al.: Malaria and helminthic co-infection among HIV-positive pregnant women: prevalence and effects of antiretroviral therapy. Acta Trop 2012, 124:179–84. 10.1016/j.actatropica.2012.08.004
Langhorne J, Ndungu FM, Sponaas AM, Marsh K: Immunity to malaria: more questions than answers. Nat Immunol 2008, 9:725–32.
Bangs MJ, Subianto DB: El Niño and associated outbreaks of severe malaria in highland populations in Irian Jaya, Indonesia: a review and epidemiological perspective. Southeast Asian J Trop Med Public Health 1999, 30:608–19.
MIDIMAR: National Disaster Risk Management Plan; Ministry of Disaster Management and Refugee Affairs. September: Kigali; 2013.
Appiah-Darkwah I, Badu-Nyarko SK: Knowledge of malaria prevention and control in a sub-urban community in Accra, Ghana. Int J Trop Med 2011, 6:61–9.
Janssen W: Insecticide-treated nets: efficacy, impact and operational constraints. Humanitarian Exchange 2005, 31:11–3.
WHO: Indoor Residual Spraying: Use of Indoor Residual Spraying for Scaling up Global Malaria Control and Elimination. Geneva: Switzerland: World Health Organization; Global Malaria Programme; 2006.
Mwangangi JM, Midega J, Kahindi S, Njoroge L, Nzovu J, Githure J, et al.: Mosquito species abundance and diversity in Malindi, Kenya and their potential implication in pathogen transmission. Parasitol Res 2012, 110:61–71. 10.1007/s00436-011-2449-6
Yé Y, Hoshen M, Louis V, Séraphin S, Traoré I, Sauerborn R: Housing conditions and Plasmodium falciparum infection: protective effect of iron-sheet roofed houses. Malar J 2006, 5:8. 10.1186/1475-2875-5-8
Bowen HL: Impact of a mass media campaign on bed net use in Cameroon. Malar J 2013, 12:36. 10.1186/1475-2875-12-36
Prue CS, Shannon KL, Khyang J, Edwards LJ, Ahmed S, Ram M, et al.: Mobile phones improve case detection and management of malaria in rural Bangladesh. Malar J 2013, 12:48. 10.1186/1475-2875-12-48
Roll Back Malaria: Global Strategic Plan 2005–2015. Geneva: Switzerland: Roll Back Malaria Partnership, World Health Organization; 2005.
WHO: The African Summit on Roll Back Malaria, Abuja. Geneva: World Health Organization; 2000. WHO/CDS/RBM/2000.17
O’Meara WP, Noor A, Gatakaa H, Tsofa B, McKenzie FE, Marsh K: The impact of primary health care on malaria morbidity - defining access by disease burden. Trop Med Int Health 2009, 14:29–35. 10.1111/j.1365-3156.2008.02194.x
Nardo M, Saisana M, Saltelli A, Tarantola S: Tools for Composite Indicators Building. Ispra, Italy: European Commission; Joint Research Centre; Institute for the Protection and Security of the Citizen Econometrics and Statistical Support to Antifraud and the Yale Center for Environmental Law and Policy; 2005.
Hildebrand DK: Statistical Thinking for Behavioral Scientists. Boston: Duxbury; 1986.
Saisana M: A do-it-Yourself Guide in Excel for Composite Indicator Development. Joint Research Centre, Italy: European Commission; 2012.
Osborne JW, Overbay A: The power of outliers (and why researchers should always check for them). Practical Assessment, Research and Evaluation 2004, 9:6.
Landau S, Everitt BS: A Handbook of Statistical Analyses Using SPSS. Boca, Raton, London, New York, Washington, D.C: Chapman & Hall/CRC Press LLC; 2004.
Malczewski J: GIS and Multicriteria Decision Analysis. New York, NY: John Wiley and Sons; 1999.
Munda G, Saisana M: Methodological considerations on regional sustainability assessment based on multicriteria and sensitivity analysis. Reg Stud 2011, 42:261–76.
Field AP: Discovering Statistics Using SPSS for Windows. London: SAGE Publication; 2005.
Tabachnick BG, Fidell LS: Using Multivariate Statistics. MA: Allyn & Bacon: Needham Heights; 2001.
OECD: Handbook on Constructing Composite Indicators: Methodology and User Guide. Paris: Organization for Economic Co-operation and Development; 2008.
Giné R, Pérez-Foguet A: Improved method to calculate the water poverty index at local scale. J Environ Eng 2010, 136:1287–98. 10.1061/(ASCE)EE.1943-7870.0000255
Slottje DJ: Measuring the quality of life across countries. Review of Economics and Statistics 1991, 73:684–93. 10.2307/2109407
Tate E: Social vulnerability indices: a comparative assessment using uncertainty and sensitivity analysis. Nat Hazards 2012, 63:325–47. 10.1007/s11069-012-0152-2
Lung T, Lavalle C, Hiederer R, Dosio A, Bouwer LM: A multi-hazard regional level impact assessment for Europe combining indicators of climatic and non-climatic change. Glob Environ Chang 2013, 23:522–36. 10.1016/j.gloenvcha.2012.11.009
National Insitute of Statistics of Rwanda: Rwanda Demographic and Health Survey. Kigali, Rwanda: National Institute of Statistics of Rwanda; Preliminary Report; 2011.
Martens P, Hall L: Malaria on the move: human population movement and malaria transmission. Emerg Infect Dis 2000, 6:103–9. 10.3201/eid0602.000202
Government of Rwanda: The Third Integrated Household Living Conditions Survey -EICV3; Main Indicators Report. Kigali, Rwanda: Ministry of Finance and Economic Planning; 2011.
Nonaka D, Laimanivong S, Kobayashi J, Chindavonsa K, Kano S, Vanisaveth V, et al.: Is staying overnight in a farming hut a risk factor for malaria infection in a setting with insecticide-treated bed nets in rural Laos? Malar J 2010, 9:372. 10.1186/1475-2875-9-372
WHO: Malaria Early Warning Systems: Concepts, Indication and Partners. “A Framework for Field Research in Africa”. Geneva: World Health Organization; 2001.
National Institute of Statistics of Rwanda: Rwanda: Comprehensive Food Security and Vulnerability Analysis and Nutrition Survey. Kigali, Rwanda: National Institute of Statistics of Rwanda; 2009.
Nabahungu NL, Visser SM: Farmers’ knowledge and perception of agricultural wetland management in Rwanda. Land Degrad Dev 2013, 24:363–74. 10.1002/ldr.1133
Abeku TA: Response to malaria epidemics in Africa. Emerg Infect Dis 2007, 13:681–86. 10.3201/eid1305.061333
Vincke IH, Jadin JB: Contribution à l’étude de l’anophélisme en pays d’altitude. Ann Soc Belge de Méd Trop 1946, 26:483–500.
Gascon J, Merlos A, Madrenys N: Coll: Epidemiology of malaria in Nyarutovu (Rwanda): a clinical, parasitological and serological study. Trans R Soc Trop Med Hyg 1988, 82:421–2.
Loevinsohn ME: Climatic warming and increased malaria incidence in Rwanda. Lancet 1994, 343:714–18. 10.1016/S0140-6736(94)91586-5
Hammerich A, Campbell OMR, Chandramohan D: Unstable malaria transmission and maternal mortality – experiences from Rwanda. Trop Med Int Health 2002, 7:573–76. 10.1046/j.1365-3156.2002.00898.x
Rwanda National Malaria Control Programme: Towards Malaria Pre-Elimination: How to Sustain Achievements and get to Zero Malaria Deaths in Rwanda. Kigali, Rwanda: Concept note for the first Rwanda Malaria Forum; 26th to 28th September 2012 Serena Hotel; 2012.
Mamel WW: Health dilemma at the borders-global challenge. World Health Forum 1997, 16:9–16.
Bhumiratana A, Intarapuk A, Sorosjinda-Nunthawarasilp P, Maneekan P, Koyadun S: Border malaria associated with multidrug resistance on Thailand-Myanmar and Thailand-Cambodia borders: transmission dynamic, vulnerability, and surveillance. Biomed Res Int 2013, 2013:363417.
Blumenstock JE: Inferring patterns of internal migration from mobile phone call records: evidence from Rwanda. Inf Technol Dev 2012, 18:107–25. 10.1080/02681102.2011.643209
Lynch C, Roper C: The transit phase of migration: circulation of malaria and its multidrug-resistant forms in Africa. PLoS Med 2011, 8:e1001040. 10.1371/journal.pmed.1001040
Rowland M: Malaria epidemiology and control in refugee camps and complex emergencies. Ann Trop Med Parasitol 2001, 95:741–54. 10.1080/00034980120103405
UNEP: Rwanda from Post-Conflict to Environmentally Sustainable Development. Nairobi, Kenya: United Nations Environment Programme; 2011.
National Institute of Statistics of Rwanda: Fourth Population and Housing Census. Thematic Report: Migration and spatial mobility. Kigali Rwanda: National Institute of Statistics of Rwanda; 2012.
National Institute of Statistics of Rwanda: Fourth Population and Housing Census. Kigali Rwanda: Thematic Report: Characteristics of Households and Housing; National Institute of Statistics of Rwanda; 2012.
MINAGRI: Étude technique d’exécution des travaux d’aménagements hydro agricoles du périmètre 8 de la Vallée Muvumba – Kagitumba, District de Nyagatare, Province de l’Est. Kigali, Rwanda: Rapport final - Evaluation de l’impact sur l’environnement; 2008.
Coosemans M, Mouchet J: Consequences of rural development on vectors and their control. Ann Soc Belge de Méd Trop 1990, 70:5–23.
Deutsch WG: Environmental assessment of ten aquaculture sites in Rwanda, Africa. Alabama: Research Development Series No 36; International Centre for Aquaculture and Aquatic Environment; Department of Fisheries and Allied Aquaculture; Auburn University; 1992.
MINAGRI: Strategic Plan for the Transformation of Agriculture in Rwanda (SPAT II). Kigali, Rwanda: Ministry of Agriculture and Animal Livestock, Kigali, Final report; 2008.
Government of Rwanda: Enabling Self Sufficiency and Competitiveness of Rwanda Rice: Issues and Policy Options. Kigali Rwanda: Ministry of Agriculture and Animal Resources Republic of Rwanda; 2010.
Larmarange J, Vallo R, Yaro S, Msellati P, Meda N: Methods for mapping regional trends of HIV prevalence from Demographic and Health Surveys (DHS). European Journal of Geography 2011., 558: doi:10.4000/cybergeo.24606 http://cybergeo.revues.org/24606;
Meade MS: Geographic analysis of disease and care. Annu Rev Public Health 1986, 7:313–35. 10.1146/annurev.pu.07.050186.001525
Hongoh V, Gatewood Hoen A, Aenishaenslin C, Waaub JP, Belanger D, Michel P: Spatially explicit multi-criteria decision analysis for managing vector-borne diseases. Int J Health Geogr 2011, 10:70. 10.1186/1476-072X-10-70
The research leading to these results has received funding from the European Union’s Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 266327 (HEALTHY FUTURES, http://www.healthyfutures.eu/). The DHS Program of the US Agency for International Development (USAID) for proving the DHS data for Rwanda. The authors thank all staff from Malaria and Other Parasitic Diseases Division of Rwanda Biomedical Centre for their helpful cooperation especially Dr Corine Karema, Emmanuel Hakizimana, Dr Noella Umulisa, Dr Monique Mulindahabi, Alphonse Mutabazi, Alphonse Rukundo, and Alain Kabayiza.
The authors declare they have no competing interests.
JPB and SK have made substantial contributions to conception and design. JPB was responsible for data collection, processing, analysis, interpretation of the results and drafting the manuscript. ET and SK provided substantial input to the manuscript by revising it critically for important intellectual content. All authors read and approved the final manuscript.