Skip to main content

Spatio-temporal distribution of malaria and its association with climatic factors and vector-control interventions in two high-risk districts of Nepal



Over the last decade, the incidence of confirmed malaria has declined significantly in Nepal. The aim of this paper is to assess the spatio-temporal distribution of malaria and its association with climatic factors and vector control interventions in two high-risk districts of Nepal.


Hotspot analysis was used to visualize the spatio-temporal variation of malaria incidence over the years at village level and generalized additive mixed models were fitted to assess the association of malaria incidence with climatic variables and vector control interventions.


Opposing trends of malaria incidence were observed in two high-risk malaria districts of eastern and far-western Nepal after the introduction of long-lasting insecticidal nets (LLINs). The confirmed malaria incidence was reduced from 2.24 per 10,000 in 2007 to 0.31 per 10,000 population in 2011 in Morang district but increased from 3.38 to 8.29 per 10,000 population in Kailali district. Malaria hotspots persisted mostly in the same villages of Kailali district, whereas in Morang district malaria hotspots shifted to new villages after the introduction of LLINs. A 1° C increase in minimum and mean temperatures increased malaria incidence by 27% (RR =1.27, 95% CI =1.12-1.45) and 25% (RR =1.25, 95% CI =1.11-1.43), respectively. The reduction in malaria incidence was 25% per one unit increase of LLINs (RR =0.75, 95% CI =0.62-0.92). The incidence of malaria was 82% lower in Morang than in Kailali district (RR =0.18, 95% CI =0.11-0.33).


The study findings suggest that LLIN coverage should be scaled up to entire districts rather than high-incidence foci only. Climatic factors should be considered for malaria micro-stratification, mosquito repellents should be prescribed for those living in forests, forest fringe and foothills and have regular visits to forests, and imported cases should be controlled by establishing fever check posts at border crossings.


Nepal has made significant progress in reducing the incidence of malaria by more than 84% over the last decade [1]. It has already achieved and exceeded the target of the Millennium Development Goals (MDGs) and universal coverage of malaria control interventions, and the Roll Back Malaria (RBM) targets of 2010 [15]. However, an increasing proportion of Plasmodium falciparum and imported malaria has been observed in the last decade [1]. Although the analysis of malaria data from 31 malaria-risk districts has shown a significant decline of confirmed malaria and Plasmodium vivax malaria incidences, no decline was found for P. falciparum and clinically suspected malaria incidences [1]. Malaria cases in Nepal are highly clustered and vector control interventions are not uniformly implemented. Accurate identification of malaria foci or clusters at the local level can greatly increase the effectiveness of vector control interventions whereas not identifying foci can cause effective control measures to fail [6]. Hence, an accurate mapping of malaria-endemic foci at village development committee (VDC) or household levels improves the effectiveness and cost-effectiveness of vector control interventions.

Although about 84% (23 million) of the people in Nepal were estimated to be at risk of malaria in 2012, with 4% at high-risk [7], the recent microstratification of malaria-risk areas in 2012 at the VDC level shows only approximately 13.02 million people (48%) living in malaria-endemic VDCs [2]. Out of the estimated total population living in endemic areas, 0.98 million (3.6%) live in high-risk VDCs, 2.7 million (9.8%) live in moderate-risk VDCs, and 9.4 million (34.5%) in low-risk VDCs. Similarly, there are fewer districts classified as high- and moderate-risk (25 instead of 31) and the overall population living in VDCs at risk (estimated at 1,254 VDCs out of 3,972) is declining [2], which indicates the progress made to eliminate malaria-endemic VDCs (foci) and to achieve the malaria elimination goal by 2026.

The high-risk areas consist of forest, forest fringe, foothills, river belts, hills, and river valleys. Malaria incidence has declined in districts regardless of whether vector control interventions (i.e., implementation of indoor residual spraying (IRS) and distribution of long-lasting insecticidal nets (LLINs)) have been implemented [1]. Moreover, reports of malaria have increased in high-risk districts (Kailali, Nawalparasi, Dhanusha, and Mahottari) where vector control interventions were in place [5]. The Global Fund to Fight AIDS, Tuberculosis and Malaria (GFATM) began supporting a malaria control programme in 13 high-risk districts in 2004 with support for rapid diagnostic test (RDT) kits, artemisinin combination therapy (ACT), LLINs and information, education and communication/behaviour change communication (IEC/BCC) for LLIN use [5]. The distribution of LLINs started in 2006 in both districts.

Climatic factors such as temperature, humidity and rainfall play important roles in malaria transmission [8]. Temperature is the major determinant of malaria risk [9], but rainfall is important where temperature is not a limiting factor (e.g., in tropical areas of Africa and Bangladesh) [10, 11]. Humidity is a direct product of temperature and rainfall and affects the distribution and longevity of vector mosquitoes, which in turn affects malaria transmission [8]. The effect of climate change on malaria transmission in temperate regions and tropical highlands has been reported in many recent studies [1214].

The main aim of this study was to assess the spatio-temporal distribution of malaria and its association with climatic factors and vector control interventions in two high-risk districts of Nepal.


Study area

A detailed description of Nepal and its administrative and geographic divisions has been provided in a previous publication [1]. Out of 13 high-risk malaria districts in Nepal, one from the eastern development region (Morang) and one from the far-western development region (Kailai) of Nepal were selected for this study. The study areas are shown in Figure 1. These two districts contain a substantial proportion of all malaria cases in the respective regions. According to the 2011 Census of Nepal, the total population of Morang district is 964,709 and that of Kailali district, 775,709 [15]. Morang district has 66 VDCs including one submetropolitan city, Biratnagar. Similarly, Kailali district has 44 VDCs including three municipalities. Although a new classification of VDCs and municipalities was established in 2014, the old classification was used for this study because the available data conform to the old classification. The area of Morang district covers 1,855 sq km and that of Kailali district 3,235 sq km. Both comprise lowlands of the terai, forest and foothills, and extend to the hills. The principal malaria vectors of study areas are Anopheles fluviatilis, Anopheles annularis and Anopheles maculatus complex members [1, 16, 17].

Figure 1
figure 1

Classification of malaria risk districts in Nepal and study areas. This figure is updated from Dhimal et al. 2014 [1].

Research design

A retrospective study was designed to assess the association between malaria incidence, climatic variables and vector control interventions in two high-risk malaria districts of Nepal between 2004 and 2012.

Data collection

A data collection tool was used to enter the number of houses and population covered by LLINs and IRS in each VDC or municipality. These figures were then submitted to District (Public) Health Office (DPHO) who in turn produced VDC level reports that were latter consolidated at both district and national level for each year. The VDCs were selected for LLIN distribution considering the following criteria [3]:

  1. 1)

    LLINs were provided to population of all wards of high-to-moderate risk VDCs except the population of VDCs protected by IRS.

  2. 2)

    Forecasting of LLINs was done as population projected to be covered/1.8.

Two rounds of IRS were carried out routinely in high-risk VDCs. The first round of IRS (IRS1) was undertaken during pre-monsoon (April-May) and the second round of IRS (IRS2) in monsoon (July-August) each year [1, 3]. In the last five years, synthetic pyrethroid insecticides such as lambda-cyhalothrin (25 mg/m2), deltamethrin (20 mg/m2) and alpha-cypermethrin (30 mg/m2) have been used in rotation for IRS, and susceptibility tests using the WHO acute test procedure showed that all anopheline mosquitoes tested were susceptible to these insecticides [3, 5]. Although LLINs and IRS are complementary interventions, the present policy is to increase the coverage of LLINs and to reduce IRS. During 2007–2008, LLINs were distributed with a policy of one LLIN per household but this changed to one LLIN per two persons in a house in 2009 [1, 3]. In addition, LLINs were also provided free of cost by DPHO to all pregnant women attending antenatal care clinics at the public health institutions. The LLIN distribution strategy is to target one third of VDCs per year in each district and cover all high-risk VDCs by the end of the third year [18]. In GFAMT program districts, LLINs were distributed through mass campaigns. The LLINs distribution started in both districts from 2006 to 2008 by Population Service International (PSI)/Nepal in partnership with local NGO through social marketing. After 2008, social marketing stopped and LLINS were distributed freely by PSI in the GFAMT program areas. The IRS was done by DPHO. The DPHO collected and stored IRS and LLINs data. Finally, these data were reported to Epidemiology and Disease Control Division (EDCD) through health management information system (HMIS). Monthly aggregated district-level malaria data between 2004 and 2012 were obtained from the EDCD, Department of Health Services (DoHS), Ministry of Health and Population (MoHP), Government of Nepal. Similarly, annually aggregated VDC or municipality level (i.e., the smallest administrative unit in each district) malaria indicators and vector control intervention data were also obtained from the EDCD. These data were collected during a malaria microstratification study in 2012–2013. However, VDC level data were compiled only between 2007 and 2011 because the complete data for 2012 were not available at the time of data collection and vector control interventions (mainly distributing LLINs only) started in 2006 in these two districts. The details of the malaria surveillance system in Nepal have been described in a previous publication [1]. The data were verified by personnel of EDCD and cleaned independently by two members of this study team.

All of the malaria cases used in this study were confirmed cases, either by microscopy or RDT kits. Data on the population at risk of malaria at both VDC and district levels were obtained from the EDCD and were projected using national population census data of 2001 and 2011 [15, 19]. Records of both imported and indigenous cases were based on travel history as recorded in the Nepalese HMIS system. Indigenous malaria is defined as “Malaria acquired by mosquito transmission in an area where malaria is a regular occurrence” [20]. Geo-environmental data were collected based on classification of VDCs using topographic maps of the study districts. The monthly accumulated rainfall (mm), air temperature (minimum and maximum) (°C) and average relative humidity (RH) (%) of respective districts were obtained from the Department of Hydrology and Meteorology (DHM), Ministry of Science, Technology and Environment, Government of Nepal.

Hotspot analysis

Annual malaria incidence at the VDC level was analysed independently for spatial clustering (or hotspots) using the Getis-Ord Gi* statistic [2123] in the Arc GIS software version 10 (ESRI, Redland, CA, USA). This statistic reflects whether differences between the local mean (i.e., the incidence for a VDC and its nearest neighbouring VDC within a district) was significantly different from the global mean (i.e., the overall incidence of all VDCs for that particular district) [22, 24, 25]. A statistically significant positive z-score value shows a hotspot for high incidence rates while a statistically significant negative z-score value for a VDC specifies local spatial clustering of low incidence rates [22, 2426]. Two separate cluster analysis was performed and presented with one using all malaria cases and one with just imported malaria cases.

Statistical analysis

The data were entered in Microsoft Excel and analysed in R computing software [27]. Assuming an average net life of three years and one net per two household residents, VDC level coverage of LLINs were calculated per person [1, 3]. Similarly, VDC level IRS coverage was calculated per person using number of households covered by IRS and population residing in those households. The incidence of malaria was expressed per 10,000 population at risk of malaria. Generalized additive mixed models (GAMM) were used to assess the effects of climatic variables on malaria incidence using the district-level monthly aggregated data, and the effects of vector control interventions on malaria incidence using annually aggregated VDC level data. A GAMM is an extension of generalized linear models that allows random effects to be included in the predictor, and non-parametric smoothing terms in the place of the constant parameters [28, 29].

The use of the generalized additive model (GAM) approach is useful for this study because it provides a flexible method to identify nonlinear covariate effects in exponential family and other likelihood-based regression methods [30]. Furthermore, instead of estimating a single parameter, GAM provides a general unspecified (non-parametric) function that relates the predicted (transformed) response values to the predictor values. The Poisson distribution with a log link function for the effects of climatic factors and vector-control interventions on the malaria incidence was used. As the collected data of malaria, vector control interventions and climate were of different spatial (district and village level) and temporal (annual and monthly) scales, the following two separate models were fitted to the data and can be summarised as:

District level malaria incidence ~ 1 | Year + District + s Rainfall + s Temperature + s Relative Humidity + s Month Model I and V D C level malaria incidence ~ 1 | V D C + s I R S 1 + s I R S 2 + s LLINs + Year + District + G e o environmental region Model I I

The random effects are denoted by (1|X). In the first model (Model I), year was used to model variation between districts and months. In the second model (Model II), VDC was used to model variation between year, districts and geo-environmental regions. Similarly, a spline (denoted s(x)) was fitted to assess the association between malaria incidence and climatic factors (monthly data) and vector control interventions (IRS and LLINs) using the annually aggregated data. Climatic factors and vector control interventions seem to have linear effects as their effective degrees of freedom (‘edf’) were estimated as 1, which indicates that the spline is not distinguishable from a straight line. In order to assess changes in malaria incidence, risk ratios (RR) were calculated to observe differences in means between discrete time periods or places [31]. The log of the malaria incidence rate is expected to be linearly associated with the vector control interventions or climatic variables of different time periods, and the model parameters after exponentiation can be interpreted as RR, which is similar to relative risk or relative incidence [1, 32]. Models were fitted in R using its ‘mgcv’ package [28].

Ethical considerations

The Ethical Review Board of the Nepal Health Research Council (NHRC) approved the conduct of this study. Only data that had been approved and documented by the EDCD, DoHS, Government of Nepal, were used for this study.


Changes in the burden of malaria in two malaria high-risk districts

Opposing trends in malaria incidences were observed in the districts after the introduction of LLINs (Table 1). Confirmed malaria incidence was reduced from 2.24 per 10,000 in 2007 to 0.31 per 10,000 population in the Morang district of eastern Nepal. In contrast, confirmed malaria incidence increased from 3.38 per 10,000 to 8.29 per 10,000 population in Kailali district of far-western Nepal. More than 50% of the cases in Kailali district were attributed to imported malaria cases. The proportion of both indigenous malaria cases in general and that of indigenous P. falciparum cases were higher in Morang than Kailali district. The LLINs and first IRS coverage were higher in Kailali district compared to Morang district. The incidence of malaria was 82% lower in Morang than in Kailali district (RR = 0.18, 95% CI = 0.11-0.33).

Table 1 Reported confirmed malaria cases and vector control interventions in two high-risk districts in Nepal (2007–2011)


Applying hotspot analysis to malaria incidence data showed that significant malaria hotspots (P <0.01, z > 2.58) were present at the VDC level in both districts. In Morang district, significant malaria hotspots (all malaria cases) shifted over the years from the eastern to the north-western part of the district as shown in Figure 2. The VDCs with new hotspots in 2011 were Belbari and Indrapur (terai plains), and Kerbari and Yangsila (forest, forest fringe and foothills) of Morang dstrict. Similarly, imported malaria hotspots also shifted over the years in new VDCs (Figure 3). In contrast, despite vector control interventions, significant malaria hotspots persisted in the same VDCs over the years, mainly in Malakheti, Godwari and Sahajpur (forest, forest fringe and foothill regions) of Kailali district which are located in the western part of the district neighbouring Kanchanpur district (Figures 4 and 5).

Figure 2
figure 2

Confirmed malaria hotspots in Morang district (2007–2011).

Figure 3
figure 3

Imported malaria hotspots in Morang district (2007–2011).

Figure 4
figure 4

Confirmed malaria hotspots in Kailali district (2007–2011).

Figure 5
figure 5

Imported malaria hotspots in Kailali district (2007–2011).

Associations between climatic factors and the burden of malaria

The minimum temperature, mean temperature and average RH were linearly associated with malaria incidence. The time series of monthly malaria cases of Morang and Kailaili district with temperature, humidity and rainfall are shown in Figures 6 and 7, respectively. Among the climatic variables only minimum temperature and RH were significant predictors of malaria incidence (Figure 8). Overall, a 1°C increase in minimum temperature increased malaria incidence by 27% (RR = 1.27, 95% CI = 1.12-1.45) and a 1% increase in mean RH decreased malaria incidence by 9% (RR = 0.91, 95% CI = 0.83-1.00). However, when mean temperature was used in model instead of minimum temperature, only mean temperature was a significant predictor. Overall, a 1°C increase in temperature increased malaria incidence by 25% (RR = 1.25, 95% CI = 1.11-1.43) (see Additional file 1). The effects of maximum temperature and rainfall on malaria incidence were not significant.

Figure 6
figure 6

Monthly malaria cases of Morang district with temperature, relative humidity and rainfall from 2004 to 2012.

Figure 7
figure 7

Monthly malaria cases of Kailali district with temperature, relative humidity and rainfall from 2004 to 2012.

Figure 8
figure 8

Effects of climatic factors and month on malaria incidence in Morang and Kailali districts (2004 to 2012). The solid line is the estimated effect, grey polygon is the 95% confidence region. The y-axis of each plot represents risk ratio of malaria incidence in log scale.

Associations between vector control interventions and the burden of malaria

The number of confirmed malaria cases and vector-control interventions (average IRS and LLINs coverage) in Morang and Kailai district is shown in Figures 9A and B, respectively. The effects of vector-control interventions and year on malaria incidence are presented in Figure 10. The decline in the burden of malaria was associated with LLIN coverage at the VDC level. Malaria incidence was reduced by 25% per one unit increase of LLINs (RR = 0.75, 95% CI = 0.62-0.92).The effect of both rounds of IRS was not significantly associated with malaria incidence. The combined effect of differences in changes in intervention over time shows significant effect of year on malaria incidence. The incidence of s malaria was significantly higher in 2010 (RR = 1.79, 95% CI = 1.38-2.32) and in 2011 (RR = 1.48, 95% CI = 1.16-1.9) compared to 2007. The incidence of malaria was significantly lower in hills and river valleys compared to plain terai (RR = 0.29, RR = 0.11-0.76). Although statistically not significant, the incidence of malaria was significantly higher in forest, forest fringe and foothills compared to plain terai (RR = 1.68, 95% CI = 0.90 – 3.11).

Figure 9
figure 9

Confirmed malaria cases and vector-control interventions in Morang and Kailali districts (2007–2011). Panel A shows confirmed malaria cases, LLINs and average IRS coverage of Morang district. Panel B shows confirmed malaria cases, LLINs and average IRS coverage of Kailali district.

Figure 10
figure 10

Effects of vector-control interventions and year on malaria incidence in Mornag and Kailali district (2007–2011). The y-axis of each plot represents risk ratio of malaria incidence in log scale.


Nepal has prepared for malaria pre-elimination since 2011 with the ambitious goal of malaria elimination by 2026. However, shifts of malaria hotspots to new VDCs in Morang district and stable malaria hotspots persist in specific VDCs of Kailali district, despite more than six years of continuous vector control interventions using IRS and LLINs, indicates that there are problems with the malaria elimination efforts in the country.

Despite declining trends of malaria in Morang district, malaria hotspots have shifted to new VDCs, which in the past had been regarded as low-risk VDCs so that vector control interventions were not in place. The hotspots were reported in eastern villages in 2009 since LLINs started to be distributed in high endemic VDCs of Morang district with one LLIN per two household persons policy in 2009. Malaria spread west in 2010 because eastern village hotspots disappeared and gradually shifted into forest and hill areas of the western part of the district. In contrast, despite vector control interventions since 2006, hotspots remained in the same VDCs in Kailali district. Both districts share a border with high-endemic districts in Nepal as well as with India. The scaling-up of LLINs started in Kailali and Morang districts in 2006. However, IRS activity was regular before and during the study period. Only LLINs is found to be associated with a decline of malaria incidence in the present study, which is consistent with findings from Bangladesh [24] Rwanda [33], and Zambia [34]. In contrast to findings of this study, malaria incidence was associated with IRS coverage only in Botswana [22]. These findings imply that vector control interventions with a low coverage, or those that are focussed on selected VDCs only, cannot break down malaria transmission in VDCs or districts bordering malaria-endemic areas because the presence of vectors, a suitable climate and the continuous import of malaria cases facilitate the transmission cycle. Furthermore, the higher incidence of malaria in the forest, forest fringe and foothills, especially in the high-risk VDCs of Kailali district where regular movement of people is high, indicates that forest-related malaria cannot be controlled using vector control interventions such as LLINs and IRS alone. For this, effective BCC strategies are needed to promote avoiding mosquito-man contact, e.g., by the use of mosquito repellents and/or protective clothing.

The proportion of imported malaria cases has almost constant over the years in both districts, consisting of more than 50% in Kailali district and 30% in Morang district (Table 1), which is consistent with previous findings [1]. Imported malaria has appeared as a major challenge for many countries embarking on malaria elimination [22, 3537]. Large-scale migration within and outside the district, across the southern border of Nepal with India, displacement of the population after natural disasters such as floods, and the continuous introduction of malaria cases among adult migrant workers could be seen as factors for the persistence of malaria, especially in Kailali district. Hence, a cross-border malaria strategy that is well coordinated with malaria-endemic neighbouring counties such as India is very urgently needed. In the absence of the implementation of such a strategy, regular border screening among travellers should be carried out by establishing check posts for acute febrile illness at the border crossings, which may contribute to a reduction in malaria transmission among other benefits.

Climatic factors can be important for increasing the risk of malaria transmission especially in sub-tropical and temperate regions where the minimum temperature is a limiting factor for malaria transmission. The positive association of minimum and mean temperature with malaria incidence in the present study is consistent with the findings of many previous studies [1214, 3845]. An increase in the minimum (and also mean) temperature increases mosquito abundance and biting rates and shortens the incubation period of Plasmodium parasites, thereby increasing the malaria transmission cycle. Minimum temperature is the most influential environmental variable for malaria transmission since it occurs at night [46]. This is because malaria vector mosquitoes remain active for biting during the night when people sleeping without bed nets are exposed to the bite of infected mosquitoes, resulting malaria infection. Furthermore, when exposed to high temperatures at night people usually do not cover themselves and some people (usually adult men) sleep outside the house under the trees to avoid the heat, which in turn increases the risk of malaria [4648]. These observations are consistent with the findings of this study that the incidence of malaria was higher among adults. In contrast, maximum temperature has a complex relationship to malaria because an increase in maximum temperature above a certain range interrupts mosquito and parasite development [8, 14]. Similarly, the effect of rainfall on malaria incidence is complex. In areas where their breeding sites are produced by rainfall, increasing rainfall increases mosquito populations. However, too little rain, or drought, affects the mosquito life cycle as well as too much rainfall which can flush away the breeding places and thus decrease mosquito abundance [11, 4951].

No significant effect of the month of the year on malaria incidence was observed, which indicates a perennial distribution of malaria in both districts. Thus, IRS spraying in the pre-monsoon (April-May) and monsoon (July-August) seasons only may be ineffective for controlling malaria vectors. This observation can be explained by the fact that the minimum temperature increased rapidly in the terai region in all seasons [52, 53], so that the transmission of malaria was possible even in winter. In contrast, an increase in mean temperature above 28°C [14] mainly influenced by maximum temperature drastically reduces malaria transmission, particularly in the summer, which is the warmest season in Nepal. As a result of these two contrasting effects, coupled with an influx of infected migrant workers in the post-monsoon and winter seasons, the effect of the month on malaria incidence may be insignificant in this setting. Alternatively, there might be no effect of the month because temperature alone is sufficient to explain seasonal variation. Based on the association found between malaria incidence and temperature, malaria transmission in the temperate regions of Nepal can be predicted as revealed elsewhere by many previous studies [1214, 45].

The latest microstratification of malaria at the VDC level may be very useful for an effective utilization of limited resources. However, malaria microstratification should be updated regularly according to the progress with an aim to target new hotspots. As the country progresses towards elimination, finer scale mapping, i.e., at the ward or household level, is needed to identify residual foci [5456]. For example, a study from Bangladesh identifies stable malaria hotspots and risk factors at the household level which guide for cost-effective targeting of malaria intervention that may finally contribute to potential elimination of malaria from the country [51]. The application of spatial decision support tools such as geographic position system (GPS), geographical information system (GIS) and mobile computing technology helps to identify spatial clusters of malaria transmission and provides effective monitoring, evaluation and surveillance tools to cope with the complexities that are associated with the spatial variability of malaria transmission and associated risk factors [57].

This study, like other studies that use secondary data, has several limitations so the findings should be interpreted with caution. First, the analyses were based on routinely collected passive surveillance data from public health institutions only. As reported elsewhere [2], malaria cases reported to the EDCD and malaria cases reported through the (HMIS) vary slightly. Monthly malaria data used in this study were aggregated at the district level and collected through the HMIS while yearly aggregated VDC-level malaria and vector control intervention data were collected from visits to each health institution of a district and from the EDCD, which resulted in slight differences in data in some years. Second, the observed association between malaria incidence with LLINs and IRS coverage and climatic factors is ecological and not at the level of individuals. Third, malaria incidence was calculated based on confirmed malaria cases at public health institutions only and therefore misses out on possible cases from the private healthcare sector. As a result, the results of this study may not represent the situation of malaria transmission at the population level. Furthermore, a single model could not be developed since the collected data of malaria, vector control interventions and climate factors were of different spatial (district and village level) and temporal (annual and monthly) scales. Despite these challenges, this study provides important information about the malaria situation at district and subdistrict (VDC) levels after the scaling-up of malaria control interventions from the GFAMT support, and this will be important for preparing the malaria pre-elimination phase in Nepal.


Despite a significant decline in malaria cases at the national level, an increasing trend of malaria incidence in Kailali district, with persistence of malaria cases in the same villages where vector control interventions had been in place, a shift of malaria hotspots to new villages in Morang district without vector control interventions, and a positive association of malaria incidence with temperature indicates worries about the elimination of malaria from the country. However, the malaria elimination goal can be achieved if hotspots of malaria can be identified accurately and vector control interventions such as LLIN coverage can be scaled up in the entire endemic districts rather than focusing only on selected VDCs. This is important because the movement of people within and between districts and across the border to India is high. Second, imported malaria cases should be controlled by establishing health check posts at the border crossings screening for people with acute febrile illness. Third, community-based prevalence surveys should be carried out to detect asymptomatic malaria cases, identify spatial clusters of malaria hotspots and determine the real malaria transmission situation at the population level. Fourth, data recording, reporting and surveillance systems should be strengthened and a case-based surveillance system should be started in all malaria-endemic districts. Fifth, microstratification of malaria transmission areas should be carried out integrating observed climatic data and high resolution remote sensing images, and differential diagnosis of any fever case coming from probable malaria transmission areas should be performed. Finally, continuous efforts are crucial to maintain and sustain the gains that have already been achieved.


  1. Dhimal M, Ahrens B, Kuch U: Malaria control in Nepal 1963–2012: challenges on the path towards elimination. Malar J. 2014, 13: 241-10.1186/1475-2875-13-241.

    Article  PubMed Central  PubMed  Google Scholar 

  2. DOHS: Annual Report 2069/70 (2012/2013). 2014, Kathmandu: Department of Health Services, Ministry of Health and Population, Goverment of Nepal

    Google Scholar 

  3. EDCD: Nepal Malaria Report 2010–2012. 2013, Kathmandu: Nepal Epidemiology and Disease Control Division (EDCD), Department of Health Services, Ministry of Health and Population, Goverment of Nepal

  4. WHO/SEARO: Nepal Malaria Programme Review (7–16 June, 2010). 2011, New Delhi: World Health Organization Regional Office for South - East Asia

    Google Scholar 

  5. EDCD: Nepal Malaria Strategic Plan 2011–2016 (Revised Version- December 2011). 2011, Kathmandu: Epidemiology and Diseases Control Division, Department of Health Services, Ministry of Health and Population, Government of Nepal

  6. Carter R, Mendis KN, Roberts D: Spatial targeting of interventions against malaria. Bull World Health Organ. 2000, 78: 1401-1411.

    PubMed Central  CAS  PubMed  Google Scholar 

  7. WHO: World Malaria Report 2013. 2013, Geneva: World Health Organization

    Google Scholar 

  8. WHO: Manual on practical entomology in malaria part 1: vector bionomicts and organization of anti-malarial activities. 1975, Geneva: World Health Organization

    Google Scholar 

  9. Paaijmans KP, Read AF, Thomas MB: Understanding the link between malaria risk and climate. Proc Natl Acad Sci U S A. 2009, 106: 13844-13849. 10.1073/pnas.0903423106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999, 15: 105-111. 10.1016/S0169-4758(99)01396-4.

    Article  CAS  PubMed  Google Scholar 

  11. Haque U, Hashizume M, Glass GE, Dewan AM, Overgaard HJ, Yamamoto T: The role of climate variability in the spread of malaria in Bangladeshi highlands. PLoS One. 2010, 5: e14341-10.1371/journal.pone.0014341.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Siraj AS, Santos-Vega M, Bouma MJ, Yadeta D, Ruiz Carrascal D, Pascual M: Altitudinal changes in malaria incidence in highlands of Ethiopia and Colombia. Science. 2014, 343: 1154-1158. 10.1126/science.1244325.

    Article  CAS  PubMed  Google Scholar 

  13. Caminade C, Kovats S, Rocklov J, Tompkins AM, Morse AP, Colon-Gonzalez FJ, Stenlund H, Martens P, Lloyd SJ: Impact of climate change on global malaria distribution. Proc Natl Acad Sci U S A. 2014, 111: 3286-3291. 10.1073/pnas.1302089111.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Mordecai EA, Paaijmans KP, Johnson LR, Balzer C, Ben-Horin T, de Moor E, McNally A, Pawar S, Ryan SJ, Smith TC, Lafferty KD: Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol Lett. 2013, 16: 22-30. 10.1111/ele.12015.

    Article  PubMed  Google Scholar 

  15. CBS: National Population and Housing Census 2011. 2012, Kathmandu: Central Bureau of Statistics, National Planning Commission Secretariat, Goverment of Nepal, 1:

    Google Scholar 

  16. Dhimal M, Ahrens B, Kuch U: Species composition, seasonal occurrence, habitat preference and altitudinal distribution of malaria and other disease vectors in eastern Nepal. Parasit Vectors. 2014, 7: 540-10.1186/s13071-014-0540-4.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Darsie RF, Pradhan S: The mosquitoes of Nepal: their identification, distribution and biology. Mosquito Syst. 1990, 22: 69-130.

    Google Scholar 

  18. Adhikari SR: Towards universal health coverage: an example of malaria intervention in Nepal. WHO South-East Asia J Public Health. 2014, 3: 103-112.

    Google Scholar 

  19. CBS: National Population Census 2001. 2002, Kathmandu: Central Bureau of Statistics, National Planning Commission Secretariat, Government of Nepal

    Google Scholar 

  20. Malaria (Plasmodium spp.) 2014 case definition.,

  21. Ord JK, Getis A: Local spatial autocorrelation statistics: distributional issues and an application. Geogr Anal. 1995, 27: 286-306.

    Article  Google Scholar 

  22. Simon C, Moakofhi K, Mosweunyane T, Jibril HB, Nkomo B, Motlaleng M, Ntebela DS, Chanda E, Haque U: Malaria control in Botswana, 2008–2012: the path towards elimination. Malar J. 2013, 12: 458-10.1186/1475-2875-12-458.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Haque U, Scott LM, Hashizume M, Fisher E, Haque R, Yamamoto T, Glass GE: Modelling malaria treatment practices in Bangladesh using spatial statistics. Malar J. 2012, 11: 63-10.1186/1475-2875-11-63.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Haque U, Overgaard HJ, Clements AC, Norris DE, Islam N, Karim J, Roy S, Haque W, Kabir M, Smith DL: Malaria burden and control in Bangladesh and prospects for elimination: an epidemiological and economic assessment. Lancet Glob Health. 2014, 2: e98-e105. 10.1016/S2214-109X(13)70176-1.

    Article  PubMed  Google Scholar 

  25. Haque U, Huda M, Hossain A, Ahmed SM, Moniruzzaman M, Haque R: Spatial malaria epidemiology in Bangladeshi highlands. Malar J. 2009, 8: 10.1186-

    Article  Google Scholar 

  26. Getis A, Ord JK: The analysis of spatial association by use of distance statistics. Geogr Anal. 1992, 24: 189-206.

    Article  Google Scholar 

  27. R Development Core Team: R: a language and environment for statistical computing 2.15.2. 2012, Vienna, Austria: R Foundation for Statistical Computing

    Google Scholar 

  28. Wood SN: Generalized Additive Models: An Introduction with R. 2006, Chapman and Hall/CRC Press

    Google Scholar 

  29. Guisan A, Edwards TC, Hastie T: Generalized linear and generalized additive models in studies of species distributions: setting the scene. Ecol Model. 2002, 157: 89-100. 10.1016/S0304-3800(02)00204-1.

    Article  Google Scholar 

  30. Hastie T, Tibshirani R: Generalized Additive Models. Stat Sci. 1986, 2: 297-318.

    Article  Google Scholar 

  31. Ceesay SJ, Casals-Pascual C, Erskine J, Anya SE, Duah NO, Fulford AJ, Sesay SS, Abubakar I, Dunyo S, Sey O, Palmer A, Fofana M, Corrah T, Bojang KA, Whittle HC, Greenwood BM, Conway DJ: Changes in malaria indices between 1999 and 2007 in The Gambia: a retrospective analysis. Lancet. 2008, 372: 1545-1554. 10.1016/S0140-6736(08)61654-2.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Nelson KE, Galai N, Safaeian M, Strathdee SA, Celentano DD, Vlahov D: Temporal trends in the incidence of human immunodeficiency virus infection and risk behavior among injection drug users in Baltimore, Maryland, 1988–1998. Am J Epidemiol. 2002, 156: 641-653. 10.1093/aje/kwf086.

    Article  PubMed  Google Scholar 

  33. Karema C, Aregawi MW, Rukundo A, Kabayiza A, Mulindahabi M, Fall IS, Gausi K, Williams RO, Lynch M, Cibulskis R, Fidele N, Nyemazi J, Ngamije D, Umulisa I, Newman R, Binagwaho A: 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.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Kamuliwo M, Chanda E, Haque U, Mwanza-Ingwe M, Sikaala C, Katebe-Sakala C, Mukonka VM, Norris DE, Smith DL, Glass GE, Moss WJ: The changing burden of malaria and association with vector control interventions in Zambia using district-level surveillance data, 2006–2011. Malar J. 2013, 12: 437-10.1186/1475-2875-12-437.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Abeyasinghe RR, Galappaththy GN, Smith Gueye C, Kahn JG, Feachem RG: Malaria control and elimination in Sri Lanka: documenting progress and success factors in a conflict setting. PLoS One. 2012, 7: e43162-10.1371/journal.pone.0043162.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Wickramage K, Galappaththy GN: Malaria burden in irregular migrants returning to Sri Lanka from human smuggling operations in West Africa and implications for a country reaching malaria elimination. Trans R Soc Trop Med Hyg. 2013, 107: 337-340. 10.1093/trstmh/trt009.

    Article  CAS  PubMed  Google Scholar 

  37. Yangzom T, Gueye CS, Namgay R, Galappaththy GN, Thimasarn K, Gosling R, Murugasampillay S, Dev V: Malaria control in Bhutan: case study of a country embarking on elimination. Malar J. 2012, 11: 9-10.1186/1475-2875-11-9.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Bhandari GP, Dhimal M, Gurung S, Bhusal CL: Climate change and malaria in Jhapa district of Nepal: emerging evidences from Nepal. J Health Manag. 2013, 15: 141-150. 10.1177/0972063413486026.

    Article  Google Scholar 

  39. Wardrop NA, Barnett AG, Atkinson JA, Clements AC: Plasmodium vivax malaria incidence over time and its association with temperature and rainfall in four counties of Yunnan Province. China Malar J. 2013, 12: 452-10.1186/1475-2875-12-452.

    Article  PubMed  Google Scholar 

  40. Bouma MJ, Dye C, van der Kaay HJ: Falciparum malaria and climate change in the northwest frontier province of Pakistan. Am J Trop Med Hyg. 1996, 55: 131-137.

    CAS  PubMed  Google Scholar 

  41. Lindsay SW, Birley MH: Climate change and malaria transmission. Ann Trop Med Parasitol. 1996, 90: 573-588.

    CAS  PubMed  Google Scholar 

  42. Patz JA, Hulme M, Rosenzweig C, Mitchell TD, Goldberg RA, Githeko AK, Lele S, McMichael AJ, Le Sueur D: Climate change: Regional warming and malaria resurgence. Nature. 2002, 420: 627-628. 10.1038/420627a.

    Article  CAS  PubMed  Google Scholar 

  43. Berrang-Ford L, Maclean JD, Gyorkos TW, Ford JD, Ogden NH: Climate change and malaria in Canada: a systems approach. Interdiscip Perspect Infect Dis. 2009, 2009: 385487-

    PubMed Central  CAS  PubMed  Google Scholar 

  44. Hui FM, Xu B, Chen ZW, Cheng X, Liang L, Huang HB, Fang LQ, Yang H, Zhou HN, Yang HL, Zhou XN, Cao WC, Gong P: Spatio-temporal distribution of malaria in Yunnan Province, China. Am J Trop Med Hyg. 2009, 81: 503-509.

    PubMed  Google Scholar 

  45. Jetten TH, Martens WJ, Takken W: Model stimulations to estimate malaria risk under climate change. J Med Entmol. 1996, 33: 361-371.

    Article  CAS  Google Scholar 

  46. Nkurunziza H, Gebhardt A, Pilz J: Geo-additive modelling of malaria in Burundi. Malar J. 2011, 10: 234-10.1186/1475-2875-10-234.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Banjara MR, Sirawaraporn W, Petmitr S, Imwong M, Joshi AB, Chavalitshewinkoon-Petmitr P: Characteristics and risk factors of Plasmodium falciparum malaria in Eastern and Central Nepal. Kathmandu Univ Med J. 2009, 7: 378-382.

    Article  CAS  Google Scholar 

  48. Parajuli K, Ghimire P: Epidemiology of malaria in two Eastern districts of Nepal. Kathmandu Univ Med J. 2010, 8: 45-50.

    Article  CAS  Google Scholar 

  49. Duncombe J, Clements A, Davis J, Hu W, Weinstein P, Ritchie S: Spatiotemporal patterns of Aedes aegypti populations in Cairns, Australia: assessing drivers of dengue transmission. Trop Med Int Health. 2013, 18: 839-849. 10.1111/tmi.12115.

    Article  PubMed  Google Scholar 

  50. Gubler DJ, Reiter P, Ebi KL, Yap W, Nasci R, Patz JA: Climate variability and change in the United States: potential impacts on vector- and rodent-borne diseases. Environ Health Perspect. 2001, 109 (Suppl 2): 223-233.

    Article  PubMed Central  PubMed  Google Scholar 

  51. Haque U, Glass GE, Bomblies A, Hashizume M, Mitra D, Noman N, Haque W, Kabir MM, Yamamoto T, Overgaard HJ: Risk factors associated with clinical malaria episodes in Bangladesh: a longitudinal study. Am J Trop Med Hyg. 2013, 88: 727-732. 10.4269/ajtmh.12-0456.

    Article  PubMed Central  PubMed  Google Scholar 

  52. PAN: Temporal and Spatial Variability of Climate Change over Nepal (1976–2005). 2009, Kathmandu: Practical Action Nepal

    Google Scholar 

  53. Baidya SK, Shrestha ML, Sheikh MM: Trends in daily climatic extremes of temperature and precipitation in Nepal. J Hydrol Meteorol. 2008, 5: 38-51.

    Google Scholar 

  54. Snow RW, Schellenberg JR, Peshu N, Forster D, Newton CR, Winstanley PA, Mwangi I, Waruiru C, Warn PA, Newbold C, Marsh K: Periodicity and space-time clustering of severe childhood malaria on the coast of Kenya. Trans R Soc Trop Med Hyg. 1993, 87: 386-390. 10.1016/0035-9203(93)90007-D.

    Article  CAS  PubMed  Google Scholar 

  55. Gaudart J, Poudiougou B, Dicko A, Ranque S, Toure O, Sagara I, Diallo M, Diawara S, Ouattara A, Diakite M, Doumbo OK: Space-time clustering of childhood malaria at the household level: a dynamic cohort in a Mali village. BMC Public Health. 2006, 6: 286-10.1186/1471-2458-6-286.

    Article  PubMed Central  PubMed  Google Scholar 

  56. malEra Consultative Group on Monitoring E, Surveillance: A research agenda for malaria eradication: monitoring, evaluation, and surveillance. PLoS Med. 2011, 8: e1000400-

    Article  Google Scholar 

  57. Kelly GC, Tanner M, Vallely A, Clements A: Malaria elimination: moving forward with spatial decision support systems. Trends Parasitol. 2012, 28: 297-304. 10.1016/

    Article  PubMed  Google Scholar 

Download references


We are grateful to the Nepal Health Research Council (NHRC) for providing ethical clearance to conduct this study. This research was financially supported by a German Academic Exchange Services (DAAD) PhD scholarship to MD and by the research funding programme LOEWE Landes-Offensive zur Entwicklung Wissenschaftlich-ökonomischer Exzellenz of the Ministry of Higher Education, Research and the Arts of the State of Hesse, Germany. We are grateful to EDCD, DoHS, MoHP, Government of Nepal for making data on malariometric indicators available for the study. We are also grateful to DHM, Government of Nepal, for making meteorological data available for this study. We are thankful to two anonymous reviewers for helpful comments.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Meghnath Dhimal.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MD designed the study, compiled and analysed the data, and wrote the manuscript; RBO’H wrote the code for statistical analysis and contributed in manuscript writing; RCK provided meteorological data required for this study; GDT provided malaria and vector control intervention data required for this study; UK contributed to the writing and critically reviewed the manuscript and BA conceived the study and critically reviewed the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dhimal, M., O’Hara, R.B., Karki, R. et al. Spatio-temporal distribution of malaria and its association with climatic factors and vector-control interventions in two high-risk districts of Nepal. Malar J 13, 457 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: