Skip to main content

Spatiotemporal epidemiology, environmental correlates, and demography of malaria in Tak Province, Thailand (2012–2015)

Abstract

Background

Tak Province, at the Thai–Myanmar border, is one of three high malaria incidence areas in Thailand. This study aimed to describe and identify possible factors driving the spatiotemporal trends of disease incidence from 2012 to 2015.

Methods

Climate variables and forest cover were correlated with malaria incidence using Pearson’s r. Statistically significant clusters of high (hot spots) and low (cold spots) annual parasite incidence per 1000 population (API) were identified using Getis-Ord Gi* statistic.

Results

The total number of confirmed cases declined by 76% from 2012 to 2015 (Plasmodium falciparum by 81%, Plasmodium vivax by 73%). Incidence was highly seasonal with two main annual peaks. Most cases were male (62.75%), ≥ 15 years (56.07%), and of Myanmar (56.64%) or Thai (39.25%) nationality. Median temperature (1- and 2-month lags), average temperature (1- and 2-month lags) and average relative humidity (2- and 3-month lags) correlated positively with monthly total, P. falciparum and P. vivax API. Total rainfall in the same month correlated with API for total cases and P. vivax but not P. falciparum. At sub-district level, percentage forest cover had a low positive correlation with P. falciparum, P. vivax, and total API in most years. There was a decrease in API in most sub-districts for both P. falciparum and P. vivax. Sub-districts with the highest API were in the Tha Song Yang and Umphang Districts along the Thai–Myanmar border. Annual hot spots were mostly in the extreme north and south of the province.

Conclusions

There has been a large decline in reported clinical malaria from 2012 to 2015 in Tak Province. API was correlated with monthly climate and annual forest cover but these did not account for the trends over time. Ongoing elimination interventions on one or both sides of the border are more likely to have been the cause but it was not possible to assess this due to a lack of suitable data. Two main hot spot areas were identified that could be targeted for intensified elimination activities.

Background

There is a goal to eliminate malaria across countries in the Greater Mekong Subregion (Cambodia, Yunnan Province in China, Lao People’s Democratic Republic (PDR), Myanmar, Thailand and Viet Nam) by 2030. In Thailand, malaria remains as a major public health priority with over 13 million people (19% of the total population) currently at risk [1], most of whom are situated along border areas with its neighbouring countries [2]. There is significant geographical heterogeneity in spatial distribution of disease incidence, which exemplifies the “border malaria” type [3], characterized by high transmission along international borders [3,4,5].

Thailand has made significant progress towards its goal of eliminating malaria by 2024 [2]. Between 2012 and 2017, Thailand reduced its malaria cases by 67%, with an accelerated decrease of 39% during the period 2016 to 2017 [6] with scale-up of malaria elimination interventions. According to the 2018 World Malaria Report, the country had 11,440 reported confirmed cases in health facilities (which include 1075 cases at community level and 3023 cases from the private sector) and 11 reported deaths due to malaria in 2017. With evidence of emerging artemisinin drug resistance in Plasmodium falciparum malaria in the Thailand–Myanmar and Thailand–Cambodia–Lao PDR border areas [7,8,9,10,11,12], including insecticide resistance [13], adding complexity to the elimination agenda in Thailand, it is increasingly important to understand the local epidemiology of malaria to guide targeted elimination efforts. Strengthening malaria surveillance has been highlighted by the World Health Organization as a core intervention to accelerate progress towards elimination [14, 15].

Tak Province, at the Thai–Myanmar border, is one of three high incidence areas in Thailand. It has historically recorded one of the highest incidences of malaria cases across the country [16]. The province reported the highest number of malaria cases nationwide in 2012 and 2013 [17], and ranked second only to Ubon Ratchathani Province the following 2 years [18]. In 2015, Tak accounted for 21% of total reported cases in Thailand [19]. Plasmodium vivax is the predominant species, reflective of the national situation [1]. The availability of detailed climate and forest cover data presents an opportunity to better understand the epidemiology of malaria in the area.

This study analysed routine malaria surveillance data from 2012 to 2015 to describe and identify possible factors driving the spatiotemporal trends of disease incidence.

Methods

Study setting

The study location was Tak Province in western Thailand. It is subdivided into nine administrative districts (Fig. 1), 60 sub-districts [20], and 493 villages, with a total population of 526,381 as of the 2010 national census [21]. The climate in Tak is characterized as tropical with rainy season from May to October [22]. The province has a total land area of 16,406.6 km2 (fourth largest among 76 Thai provinces). At its western edge, Tak Province shares international borders with Kayin State of Myanmar.

Fig. 1
figure 1

Reference map. Place names are provided for the nine districts (amphoe) subdividing the province

Malaria diagnosis

In this area, malaria is diagnosed mostly through passive case detection by light microscopy by trained microscopists at malaria posts, health promoting hospitals and district hospitals. This is supplemented by active case detection (specifically by reactive case detection) [23]. Rapid diagnostic test (RDT) (SD BIOLINE Malaria Ag Pf/Pv®) are also used. There is also a system of village volunteer health workers who facilitate access to malaria diagnosis and treatment. Data on reported malaria cases is thought to be of high completeness and quality [24].

Data collection

Information about numbers of confirmed malaria cases at sub-district level, between 1 January 2012 and 31 December 2015 were extracted from the consolidated malaria database of the Department of Disease Control, Ministry of Public Health (MOPH). This included individual data gathered from both the Bureau of vector-borne diseases (BVBD) and Bureau of epidemiology (BOE) national malaria surveillance systems. These two databases have been combined since 2012 with removal of duplicate cases from these two sources through a collaborative project supported by the Global Fund [19, 24, 25]. A confirmed case was any person with a positive malaria blood smear or RDT result reported by government health workers [2]. Data included age, sex, and nationality (nationals of Thailand = “Thai”, nationals of Myanmar = “Myanmar”, nationals of Lao PDR = “Lao” and “Other”), as well as address to sub-district level, date of diagnosis and parasite species.

Aside from malaria incidence data, other information used in this study was (i) annual population estimates from the Ministry of Interior (http://www.moi.go.th), (ii) administrative boundary layers for Thailand from the Ministry of Interior (http://www.moi.go.th), with the most recent available version being from 2015; (iii) weather station data, collected as daily measurements of rainfall (in millimetres), temperature (in degrees Celsius,  °C) and relative humidity (in percentage), from the Thai Meteorological Department (https://www.tmd.go.th); (iv) forest cover layers from combining publicly available rasters characterizing global forest extent and change from 2000 through 2016 derived from Landsat data (http://earthenginepartners.appspot.com/science-2013-global-forest) [26]. These environmental variables were chosen as they are what was available from the climate monitoring stations in Thailand and previous studies had found correlation between each of them and malaria [27,28,29]. Forest cover was used as much of the malaria transmission in Thailand and the broader Greater Mekong Subregion is known to be in forest and forest fringe areas [30]. Spatial data management and processing were done in ArcCatalog 10.3.1 (Esri, Redlands, California, USA) and Excel 2013 (Microsoft, Redmond, Washington, USA).

Data analysis

Spatial and temporal trends in malaria incidence were examined at sub-district level. Thematic maps were produced using ArcMap 10.3.1 (Esri, Redlands, California, USA). Mixed infections were added to numbers of both P. falciparum-only and P. vivax-only cases prior to mapping of numbers of confirmed cases and annual parasite incidence (API) by species. API was calculated as the number of confirmed cases per 1000 total population as estimates of population at risk in each sub-district were not available. Plasmodium falciparum-to-P. vivax ratio was the number of confirmed P. falciparum divided by the number of confirmed P. vivax monoinfections. Hotspot maps [31] of API were produced using the Hot Spot Analysis (Getis-Ord Gi*) tool in ArcMap. Spatial relationships among features (i.e. sub-district boundaries) were defined by fixed distance band, and calculated with Euclidean (straight-line) distance. Correlations between malaria incidence rate, climate measurements and percent forest cover were examined using Pearson’s r in GraphPad Prism version 7.04 for Windows (GraphPad Software, La Jolla, California, USA). Age-sex distributions of malaria cases were plotted using R version 3.4.4 (The R Foundation for Statistical Computing, Vienna, Austria). Time series analysis of monthly malaria cases was done using stl() function (“seasonal and trend decomposition using LOESS) [32], also in R.

Forest cover was calculated by combining annual rasters of forest change and baseline forest cover in 2012 to derive annual forest cover layers at 30 m resolution. These were then combined with sub-district boundary files to calculate percentage forest cover by sub-district. Change in forest cover from 2012 to 2015 was calculated directly from the annual rasters. This was correlated with change in API using Pearson’s correlation on the logged values. Daily climate data from the five climate stations in Tak Province were averaged and then summarized by month as mean, median, maximum and minimum values. These were then correlated with monthly API at province level using Pearson’s correlation.

Ethical considerations

Permission to access consolidated and de-identified malaria surveillance data was provided by the MOPH. The protocol for this study was reviewed and approved by the Faculty of Tropical Medicine Ethics Committee (FTM-EC), Mahidol University.

Results

Total cases

From 2012 to 2015, there were 36,536 confirmed malaria cases in Tak Province (Fig. 2 and Table 1). In each surveillance year, the majority of cases were infected with P. vivax. On average, these two species combined comprised at least 97% of total cases annually with only 878 cases reported as unknown species (an overall percentage of 2.40% over the four years).

Fig. 2
figure 2

Monthly malaria cases and climate station measurements in Tak Province, 2012–2015. Each bar represents the total numbers of confirmed malaria cases in the province, including those with unknown species, from January 2012 to December 2015. Coloured lines represent corresponding monthly climate measurements (relative humidity; total rainfall; minimum, maximum, and average temperature)

Table 1 Confirmed malaria cases by type of species in Tak Province, 2012–2015

The general population of Tak Province represented almost 1% of the national total during the study period. Total API in Tak has steadily declined from 2012 to 2015 by 73.93%, with a marked decrease of at least 50% from 2013 to 2014 and again from 2014 to 2015. Similar downward trends were seen for P. falciparum and P. vivax infections. Plasmodium falciparum-to-P. vivax ratio decreased from 0.64 in 2012 to 0.44 in 2015. A summary of annual API by species and P. falciparum-to-P. vivax ratio is shown in Table 2.

Table 2 Population, annual parasite incidence and Pf:Pv ratio in Tak Province, 2012–2015

Exploration of seasonal trend

The time series decomposition of malaria cases in Tak Province shows a seasonal pattern and a declining trend from January 2012 to December 2015 (Fig. 3) with the same pattern being seen for each of P. falciparum (Additional file 1: Fig. S1) and P. vivax (Additional file 1: Fig. S2). The peak number of total cases occurred in the month of June, a seasonal trend similarly seen in both major malaria species (Fig. 4). P. falciparum-to-P. vivax ratio by month is shown in Fig. 5.

Fig. 3
figure 3

Time series decomposition of total confirmed malaria cases in Tak Province, 2012–2015. Top plot is decomposed into seasonal, trend and irregular components. The grey bar on the right indicates the relative magnitude of each of the decomposed components (length of bar = relative magnitude)

Fig. 4
figure 4

Confirmed malaria cases by month in Tak Province, 2012–2015. Plasmodium falciparum and P. vivax cases shown are numbers of mono-infections

Fig. 5
figure 5

Plasmodium falciparum-to-P. vivax ratio (mono-infections) by month in Tak Province, 2012–2015

Patient demographics

Demographic characteristics of people with malaria for the province are summarized by year in Table 3. Most cases were male (n = 22,926, 62.75%), ≥ 15 years (56.07%), and Myanmar (56.64%) or Thai (39.25%) nationality. Only 86.12% (n = 31,464) of total cases had their age recorded [all from BVBD, with four outliers (age > 150)].

Table 3 Demographic profile of malaria cases in Tak Province, 2012–2015

The age-sex distributions of malaria cases are presented as annual population pyramids in Fig. 6. These are compared to the provincial population age-sex structure as of the 2010 national census. Across all years and species, there were more males than females with malaria in Tak Province (62.75% males vs 37.25% females) whereas in the general population the proportion of each sex was roughly equal (49.53% males vs 50.47% females). Additionally, 47.41% of cases belonged to the “5–19” age groups in both sexes, whereas in the census, this group comprised only 26.24%, therefore they were overrepresented in the malaria group. There were no major differences in age-sex distribution between the years or between those with P. falciparum and P. vivax.

Fig. 6
figure 6

Age-sex distributions (pyramid) of malaria cases in Tak Province, 2012–2015. First column pyramids represent the population in Tak Province as of the most recent census in 2010 (N = 526,376). The second to fifth columns show the corresponding annual patterns of confirmed malaria cases, where the rows represent type of infection (top: Total cases, middle: P. falciparum mono-infections, bottom: P. vivax monoinfections)

Spatial distribution

There was a decrease in total API in 49 sub-districts and an increase in one from 2012 to 2015. For P. falciparum, API decreased in 45 sub-districts. For P. vivax, API decreased in 46 sub-districts and an increase in one (Additional file 1: Fig. S3 for confirmed cases). For each measure, the remaining sub-districts (out of the possible 60) were considered malaria-free when there had been no malaria cases in both 2012 and 2015, and therefore no change in API. Sub-districts with the highest API were in the Tha Song Yang and Umphang Districts along the Thai–Myanmar border (Fig. 7). Annual hot spots were mostly in the extreme north and south of the province (Fig. 8). In 2015, the five sub-districts with the highest API (for each of P. falciparum, P. vivax) were Tha Song Yang (23.09, 31.76), Mokro (10.62, 28.98), Nong Luang (6.48, 19.43), Mae Chan (21.41, 16.38), and Umphang (7.13, 16.37). At least 99% of malaria cases in Tak Province were reported in the five districts that border eastern Myanmar, with most cases reported in Tha Song Yang and Umphang (Table 4).

Fig. 7
figure 7

API at sub-district (tambon) level in Tak Province, 2012–2015. Mixed infections were added to both P. falciparum (Pf) and P. vivax (Pv) cases

Fig. 8
figure 8

API hot spots by species in Tak Province, 2012–2015. Spatial relationships among features (i.e. sub-district boundaries) were defined by fixed distance band, and calculated with Euclidean (straight-line) distance

Table 4 Confirmed malaria cases by district in Tak Province, 2012–2015

Figure 9 shows numbers of malaria cases by nationality for Myanmar and Thai (Additional file 1: Fig. S4 for P. falciparum plus mixed, Additional file 1: Fig. S5 for P. vivax plus mixed). The spatial distribution was similar for both nationalities with a significant correlation in all four years (Fig. 10). When mapped as proportion of cases by nationality, non-Thai cases (93% Myanmar nationals) were more common in the centre of the province in sub-districts nearer to the Myanmar border (Fig. 11).

Fig. 9
figure 9

Number of confirmed malaria cases by nationality at sub-district (tambon) level in Tak Province, 2012–2015. Top panel maps show total Myanmar cases, while bottom panel maps depict total Thai cases

Fig. 10
figure 10

Confirmed Myanmar and Thai malaria cases at sub-district (tambon) level in Tak Province, 2012–2015. Correlation test results and axes of scatterplots are in logarithmic scale

Fig. 11
figure 11

Percent non-Thai malaria (total cases) at sub-district (tambon) level in Tak Province, 2012–2015

While many non-Thai cases are seen along the international border with Myanmar, most of them were recorded consistently in specific sub-districts, namely Tha Song Yang (of the same district name), Mae Chan and Mokro (Umphang District) (Fig. 9).

Forest cover

Satellite images of global forest change available online allowed the derivation of sub-district forest cover expressed as percentage coverage. This was correlated with API. For the whole province, the total percentage forest cover decreased by 1.35%, from 63.42% in 2012 to 62.56% in 2015. The highest percent forest loss was in Tha Song Yang District (4.97%), from 81.14 to 77.10%. The most loss by sub-district being 5.96% in Mae Song and 5.92% in Mae Wa Luang, both in Tha Song Yang District. Percent forest cover decrease by sub-district from 2012 to 2015 did not correlate strongly with change in API over the same period (Fig. 12). At sub-district level, percentage forest cover had a low positive correlation with P. falciparum, P. vivax, and total cases, within most years (Additional file 1: Fig. S6).

Fig. 12
figure 12

Percent decrease in forest cover and decrease in API at sub-district (tambon) level. Correlation test results and axes of scatterplots are in logarithmic scale. Each plot (from left: Total cases, P. falciparum plus mixed, P. vivax plus mixed) shows the results of Pearson’s correlation test

Climate

Results of correlation of monthly cases per 1000 population and climate variables are shown in Table 5. There were five climate stations in Tak Province. Median temperature (1- and 2-month lags), average temperature (1- and 2-month lags) and average relative humidity (2- and 3-month lags) correlated positively with monthly total, P. falciparum and P. vivax API. Total rainfall in the same month correlated with API for total cases and P. vivax, but not P. falciparum. There was no consistent trend in climate variables for each of the five climate stations and no significant correlation between change in climate variables and change in malaria from 2012 to 2015.

Table 5 Relationship between climate and malaria incidence in Tak Province, 2012–2015

Discussion

From 2012 to 2015, Tak Province has experienced a substantial and consistent decline of both P. falciparum and P. vivax malaria with a decline in P. falciparum-to-P. vivax ratio over the same period. Around three quarters of the cases and areas with highest risk (indicated by API) were consistently in clusters of a few sub-districts in two out of the nine districts (Tha Song Yang and Umphang) at the northern and southern extremities of the province by the Thai–Myanmar border. There was some correlation with the percent decrease in forest cover and decline in API at sub-district level. This geographical pattern was similar for both species.

Although climate correlated reasonably well with API, there was no clear change in seasonal climate pattern over the period of the study. This suggests that change in climate was not responsible for the decrease in malaria. The relationship between trends in climate and malaria is not straightforward in that there is an optimal ecological niche of climate within which mosquito breeding and malaria transmission are more efficient [33]. In Thailand, there are seven Anopheles species known to transmit malaria, which have different environmental optima and geographic distributions and it is not clearly known to what degree they contribute to transmission in Tak Province [34]. Because the correlation was weak, the change in forest cover is unlikely to have been the main driver overall for the reduction in malaria. It is possible that small changes in forest cover may have a greater contribution in small foci where forest transmission is occurring but it was not possible to study this, as suitable detailed data was not available. More likely, it is a result of the Thai MOPH’s ongoing elimination efforts in the area with high coverage of long-lasting insecticidal nets, good access to diagnosis and treatment via malaria posts and health promoting hospitals, active case detection and village health volunteers existing with vector control and facilitating referral of symptomatic individuals [2]. However, it was not possible to assess this as suitable data on interventions were not available. This decline in cases mirrors that of what is seen elsewhere in Thailand, and in subsequent years, with overall strong progress towards the national goal of elimination [2]. Also unique to Tak Province is the presence of Shoklo Malaria Research Unit (SMRU) in Mae Sot District, which has been providing early diagnosis and treatment for malaria since 1986 [35]. SMRU operates at both sides of the Thai–Myanmar border, and reported 17,446 malaria cases between 2012 and 2015 (Additional file 1: Table S1 for annual reported cases by species) at sites in Tak Province with an 85.64% decline in reported malaria cases from 2012 to 2015. There was a similar decline in cases reported by the Thai MOPH at 75.78% in the same period. SMRU sites in Kayin State, Myanmar also saw a decline in malaria incidence, as brought about by recent gains in elimination efforts in Myanmar [36]. Due to migration of people between the two areas, the decline in reported cases on the Thai side may also be partly due to the decline on the other side of the border. Lack of functioning malaria posts in Tha Song Yang (2014) and Umphang (2015) Districts (source: SMRU unpublished reports) may have contributed to increased transmission in the areas identified as hotspots in this analysis.

A strength of this study is that it combined detailed data from all sources reporting to the MOPH in Thailand together [37, 38], with population data from the same area. This allowed examination of which groups of people had the highest rates of malaria in each major species. Similar to other studies in the region [39,40,41] and elsewhere [42], around two-thirds were male, particularly in young adults, possibly due to increased occupational risk and forest exposure in this group. More than half the cases were of Myanmar nationality. Data were not available to indicate whether they were short-term or long-term migrants, but this study mapped cases by place of residence where this information was available (i.e. which sub-district a patient confirmed to have malaria was reported to be from). When cases were mapped by nationality, those of Myanmar nationality were clustered in a few sub-districts in the centre of the province along the country border.

The study had several important limitations. Although the data was of high quality, there was some missing information, e.g. 2.40% of cases had unknown malaria species, 13.88% of patients had no age recorded, 4.08% had no nationality recorded, and 1.30% had missing sub-district codes (Additional file 1: Table S2 for yearly percentages, and 4-year average). The completeness of data for each field did not consistently increase over the period of the study. The study utilized data obtained mostly from passive case detection with a few from active case detection (Additional file 1: Table S3 for an annual breakdown of confirmed malaria cases by detection method), but it was not clear which applied to each case. There are other sources of surveillance data that are not fully captured in the MOPH database, including those from the military and non-governmental organizations [43]. PCR is not used for routine surveillance in Thailand and there was no available data from the government on asymptomatic infections, which could render asymptomatic cases to be undetected, as some research studies have revealed in both high transmission [44, 45] and low transmission [46, 47] settings at the Thai–Myanmar border. Data were not included from SMRU, and these additional case data could be substantial enough to affect the spatial and temporal patterns of disease found in this study. However, it is not clear to what degree cases recorded by SMRU overlap with the dataset used in this study due to, e.g. presentation of cases to more than one organization within the same malaria episode, and there is no system to effectively integrate the two databases accounting for this potential overlap. Perhaps more importantly, many of these cases are short-term migrants who are likely to have been infected in Myanmar, there being a porous border between the two countries in this area. Including cases from SMRU should thus have had little impact on the results as the analysis was based on place of residence being in Tak Province. Other data that were not available for the analysis were numbers of people tested, or detailed data on interventions including bed net distribution or treatment, which could have been used to confirm the downward trend in cases being due to successful elimination efforts. Finally, the data used to map forest cover had limitations in that there was no information on forest gain from 2013 to 2015 (stable forest gain was assumed during this period), only forest loss [26].

It would be informative to explore further about possible causes of the spatial pattern of disease. For example, delays in seeking treatment [48], importation or spread through human movement [49, 50], travel into the forest [51], locations of particular occupational risk groups, or the role of undiagnosed cases [52] and parasite genetics [53]. This would require additional surveys to collect this information. Another possible future analysis would be to use satellite remote sensing data as a higher resolution source of information on climate to examine the variation across space. More detailed mapping of the distribution of cases to assess the degree of clustering would also be helpful but this would require village level data on place of residence, which was not available for this study. This is particularly important as further progress is made towards elimination requiring targeting of these clusters to reach zero. It would also be informative to include data from neighbouring areas of Myanmar [54] to assess risk at either side of the border. Given the cross-border nature of the study setting [55], and how this could affect malaria transmission, there is an opportunity to advocate for a regional data system that would allow for a combined analysis of trends for both Tak Province in Thailand and Kayin State in Myanmar.

Conclusions

There has been a large decline in reported clinical malaria from 2012 to 2015 in Tak Province. Climate and forest correlated with malaria incidence rates but did not account for this decrease. Ongoing elimination interventions on one or both sides of the border are more likely to have been the cause but it was not possible to assess this due to a lack of suitable data. Two main hot spot areas were identified along the Thai–Myanmar border that should be studied in more detail so they can be targeted for elimination activities.

Availability of data and materials

The malaria surveillance dataset is generated and owned by the Department of Disease Control, Ministry of Public Health (MOPH), Thailand, but are available on reasonable request.

Abbreviations

API:

annual parasite incidence

BOE:

Bureau of Epidemiology

BVBD:

Bureau of Vector-borne Diseases

MOPH:

Ministry of Public Health, Thailand

PDR:

Lao People’s Democratic Republic

RDT:

rapid diagnostic test

SMRU:

Shoklo Malaria Research Unit

References

  1. WHO. World malaria report 2018. Geneva: World Health Organization; 2018.

    Google Scholar 

  2. Ministry of Public Health. National malaria elimination strategy 2017–2026. Nonthaburi: Ministry of Public Health; 2016.

    Google Scholar 

  3. Cui L, Yan G, Sattabongkot J, Cao Y, Chen B, Chen X, et al. Malaria in the Greater Mekong Subregion: heterogeneity and complexity. Acta Trop. 2012;121:227–39.

    Article  Google Scholar 

  4. 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.

    Article  Google Scholar 

  5. Corbel V, Nosten F, Thanispong K, Luxemburger C, Kongmee M, Chareonviriyaphap T. Challenges and prospects for dengue and malaria control in Thailand, Southeast Asia. Trends Parasitol. 2013;29:623–33.

    Article  Google Scholar 

  6. WHO. Mekong Malaria Elimination (MME) programme. Geneva: World Health Organization; 2018.

    Google Scholar 

  7. Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, et al. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009;361:455–67.

    Article  CAS  Google Scholar 

  8. WHO. Emergency response to artemisinin resistance in the Greater Mekong Subregion: regional framework for action 2013–2015. Geneva: World Health Organization; 2013.

    Google Scholar 

  9. Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014;371:411–23.

    Article  Google Scholar 

  10. Phyo AP, Nkhoma S, Stepniewska K, Ashley EA, Nair S, McGready R, et al. Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study. Lancet. 2012;379:1960–6.

    Article  Google Scholar 

  11. Wongsrichanalai C, Sirichaisinthop J, Karwacki JJ, Congpuong K, Miller RS, Pang L, et al. Drug resistant malaria on the Thai–Myanmar and Thai–Cambodian borders, Southeast Asian. J Trop Med Public Health. 2001;32:41–9.

    CAS  Google Scholar 

  12. Phyo AP, Ashley EA, Anderson TJC, Bozdech Z, Carrara VI, Sriprawat K, et al. Declining efficacy of artemisinin combination therapy against Plasmodium falciparum malaria on the Thai–Myanmar border (2003–2013): the role of parasite genetic factors. Clin Infect Dis. 2016;63:784–91.

    Article  CAS  Google Scholar 

  13. Chaumeau V, Cerqueira D, Zadrozny J, Kittiphanakun P, Andolina C, Chareonviriyaphap T, et al. Insecticide resistance in malaria vectors along the Thailand–Myanmar border. Parasit Vectors. 2017;10:165.

    Article  Google Scholar 

  14. WHO. Global technical strategy for malaria 2016–2030. Geneva: World Health Organization; 2015.

    Google Scholar 

  15. WHO. A framework for malaria elimination. Geneva: World Health Organization; 2017.

  16. Thimasarn K, Jatapadma S, Vijaykadga S, Sirichaisinthop J, Wongsrichanalai C. Epidemiology of malaria in Thailand. J Travel Med. 1995;2:59–65.

    Article  CAS  Google Scholar 

  17. Ministry of Public Health, Thailand. Bureau of vector-borne diseases: annual report 2013. Nonthaburi: Ministry of Public Health; 2013.

    Google Scholar 

  18. Ministry of Public Health, Thailand. Bureau of vector-borne diseases: annual report 2015. Nonthaburi: Ministry of Public Health; 2015.

    Google Scholar 

  19. Malaria Map. Department of Disease Control, Ministry of Public Health, Thailand, Nonthaburi. http://mgis.ddc.moph.go.th/. Accessed 15 Aug 2016.

  20. Ministry of Interior. Directory of agencies and positions. Nonthaburi: Department of Provincial Administration; 2015.

    Google Scholar 

  21. National Statistical Office, Thailand. The 2010 population and housing census: Tak Province. In: Population Statistics Group: Social Statistics Bureau. 2012.

  22. Carrara VI, Sirilak S, Thonglairuam J, Rojanawatsirivet C, Proux S, Gilbos V, et al. Deployment of early diagnosis and mefloquine–artesunate treatment of falciparum malaria in Thailand: the Tak malaria initiative. PLoS Med. 2006;3:e183.

    Article  Google Scholar 

  23. Smith Gueye C, Sanders KC, Galappaththy GN, Rundi C, Tobgay T, Sovannaroth S, et al. Active case detection for malaria elimination: a survey among Asia Pacific countries. Malar J. 2013;12:358.

    Article  Google Scholar 

  24. Khamsiriwatchara A, Sudathip P, Sawang S, Vijakadge S, Potithavoranan T, Sangvichean A, et al. Artemisinin resistance containment project in Thailand. (I): implementation of electronic-based malaria information system for early case detection and individual case management in provinces along the Thai–Cambodian border. Malar J. 2012;11:247.

    Article  Google Scholar 

  25. Khamsiriwatchara A, Kaewkungwal J, Lawpoolsri S, Akarasewi P, Premsri N, Chinanonwait N, et al. Estimating “real” malaria burden in Thailand: Applying automated fuzzy logic algorithm for cleaning redundancy data in combining two national malaria surveillance systems. American Society of Tropical Medicine and Hygiene 64th Annual Meeting; 2015; Philadelphia, USA.

  26. Hansen MC, Potapov PV, Moore R, Hancher M, Turubanova SA, Tyukavina A, et al. High-resolution global maps of 21st-century forest cover change. Science. 2013;342:850–3.

    Article  CAS  Google Scholar 

  27. Kotepui M, Kotepui KU. Impact of weekly climatic variables on weekly malaria incidence throughout Thailand: a country-based six-year retrospective study. J Environ Public Health. 2018;2018:8397815.

    Article  Google Scholar 

  28. Wiwanitkit V. Correlation between rainfall and the prevalence of malaria in Thailand. J Infect. 2006;52:227–30.

    Article  Google Scholar 

  29. Kiang R, Adimi F, Soika V, Nigro J, Singhasivanon P, Sirichaisinthop J, et al. Meteorological, environmental remote sensing and neural network analysis of the epidemiology of malaria transmission in Thailand. Geospatial Health. 2006;1:71–84.

    Article  Google Scholar 

  30. Lyttleton C. Deviance and resistance: malaria elimination in the Greater Mekong Subregion. Soc Sci Med. 2016;150:144–52.

    Article  Google Scholar 

  31. Noe A, Zaman SI, Rahman M, Saha AK, Aktaruzzaman MM, Maude RJ. Mapping the stability of malaria hotspots in Bangladesh from 2013 to 2016. Malar J. 2018;17:259.

    Article  Google Scholar 

  32. Cleveland RB, Cleveland WS, McRae JE, Terpenning I. STL: a seasonal-trend decomposition procedure based on Loess. J Off Stat. 1990;6:3–73.

    Google Scholar 

  33. Mordecai EA, Paaijmans KP, Johnson LR, Balzer C, Ben-Horin T, de Moor E, et al. Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol Lett. 2013;16:22–30.

    Article  Google Scholar 

  34. Tainchum K, Kongmee M, Manguin S, Bangs MJ, Chareonviriyaphap T. Anopheles species diversity and distribution of the malaria vectors of Thailand. Trends Parasitol. 2015;31:109–19.

    Article  Google Scholar 

  35. Carrara VI, Lwin KM, Phyo AP, Ashley E, Wiladphaingern J, Sriprawat K, et al. Malaria burden and artemisinin resistance in the mobile and migrant population on the Thai–Myanmar border, 1999-2011: an observational study. PLoS Med. 2013;10:e1001398.

    Article  Google Scholar 

  36. Landier J, Parker DM, Thu AM, Lwin KM, Delmas G, Nosten FH, et al. Effect of generalised access to early diagnosis and treatment and targeted mass drug administration on Plasmodium falciparum malaria in Eastern Myanmar: an observational study of a regional elimination programme. Lancet. 2018;391:1916–26.

    Article  Google Scholar 

  37. Soe TN, Yeesoonsang S, Jiraphongsa C, Nguyen HAP, Thammapalo S, Tao S, et al. Malaria surveillance at Thai–Myanmar border, Mae Sot District, Tak Province, Thailand, July 2013. OSIR J. 2016;9:12–9.

    Google Scholar 

  38. Ma S, Lawpoolsri S, Soonthornworasiri N, Khamsiriwatchara A, Jandee K, Taweeseneepitch K, et al. Effectiveness of implementation of electronic malaria information system as the national malaria surveillance system in Thailand. JMIR Public Health and Surveillance. 2016;2:e20.

    Article  Google Scholar 

  39. Wangroongsarb P, Satimai W, Khamsiriwatchara A, Thwing J, Eliades JM, Kaewkungwal J, et al. Respondent-driven sampling on the Thailand–Cambodia border. II. Knowledge, perception, practice and treatment-seeking behaviour of migrants in malaria endemic zones. Malar J. 2011;10:117.

    Article  Google Scholar 

  40. Siv S, Roca-Feltrer A, Vinjamuri SB, Bouth DM, Lek D, Rashid MA, et al. Plasmodium vivax malaria in Cambodia. Am J Trop Med Hyg. 2016;95:97–107.

    Article  Google Scholar 

  41. Wang RB, Zhang J, Zhang QF. Malaria baseline survey in four special regions of northern Myanmar near China: a cross-sectional study. Malar J. 2014;13:302.

    Article  Google Scholar 

  42. Cotter C, Sturrock HJ, 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.

    Article  Google Scholar 

  43. Mercado CE, Ekapirat N, Dondorp AM, Maude RJ. An assessment of national surveillance systems for malaria elimination in the Asia Pacific. Malar J. 2017;16:127.

    Article  Google Scholar 

  44. Imwong M, Nguyen TN, Tripura R, Peto TJ, Lee SJ, Lwin KM, et al. The epidemiology of subclinical malaria infections in South-East Asia: findings from cross-sectional surveys in Thailand–Myanmar border areas, Cambodia, and Vietnam. Malar J. 2015;14:381.

    Article  Google Scholar 

  45. Sattabongkot J, Suansomjit C, Nguitragool W, Sirichaisinthop J, Warit S, Tiensuwan M, et al. Prevalence of asymptomatic Plasmodium infections with sub-microscopic parasite densities in the northwestern border of Thailand: a potential threat to malaria elimination. Malar J. 2018;17:329.

    Article  Google Scholar 

  46. Imwong M, Stepniewska K, Tripura R, Peto TJ, Lwin KM, Vihokhern B, et al. Numerical distributions of parasite densities during asymptomatic malaria. J Infect Dis. 2016;213:1322–9.

    Article  CAS  Google Scholar 

  47. Baum E, Sattabongkot J, Sirichaisinthop J, Kiattibutr K, Jain A, Taghavian O, et al. Common asymptomatic and submicroscopic malaria infections in Western Thailand revealed in longitudinal molecular and serological studies: a challenge to malaria elimination. Malar J. 2016;15:333.

    Article  Google Scholar 

  48. Sonkong K, Chaiklieng S, Neave P, Suggaravetsiri P. Factors affecting delay in seeking treatment among malaria patients along Thailand–Myanmar border in Tak Province, Thailand. Malar J. 2015;14:3.

    Article  Google Scholar 

  49. Chang HH, Wesolowski A, Sinha I, Jacob CG, Mahmud A, Uddin D, et al. Mapping imported malaria in Bangladesh using parasite genetic and human mobility data. eLife. 2019. https://doi.org/10.7554/eLife.43481.001.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Sriwichai P, Karl S, Samung Y, Kiattibutr K, Sirichaisinthop J, Mueller I, et al. Imported Plasmodium falciparum and locally transmitted Plasmodium vivax: cross-border malaria transmission scenario in northwestern Thailand. Malar J. 2017;16:258.

    Article  Google Scholar 

  51. Bhumiratana A, Sorosjinda-Nunthawarasilp P, Kaewwaen W, Maneekan P, Pimnon S. Malaria-associated rubber plantations in Thailand. Travel Med Infect Dis. 2013;11:37–50.

    Article  Google Scholar 

  52. Parker DM, Matthews SA, Yan G, Zhou G, Lee MC, Sirichaisinthop J, et al. Microgeography and molecular epidemiology of malaria at the Thailand–Myanmar border in the malaria pre-elimination phase. Malar J. 2015;14:198.

    Article  Google Scholar 

  53. Wesolowski A, Taylor AR, Chang HH, Verity R, Tessema S, Bailey JA, et al. Mapping malaria by combining parasite genomic and epidemiologic data. BMC Med. 2018;16:190.

    Article  Google Scholar 

  54. Landier J, Parker DM, Thu AM, Carrara VI, Lwin KM, Bonnington CA, et al. The role of early detection and treatment in malaria elimination. Malar J. 2016;15:363.

    Article  Google Scholar 

  55. Parker DM, Carrara VI, Pukrittayakamee S, McGready R, Nosten FH. Malaria ecology along the Thailand–Myanmar border. Malar J. 2015;14:388.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the staff at the Center of Excellence for Biomedical and Public Health Informatics BIOPHICS for assistance with extracting malaria surveillance data used in this study. We show appreciation to Ipsita Sinha, Win Zaw, Vilasinee Yuwaree, Jin Young Cha, Nattawan Aunnarongrat, Chawarat Rotejanaprasert, Mavuto Mukaka and James Watson (Mahidol Oxford Tropical Medicine Research Unit MORU) for their support. We are also grateful to Daniel Parker (University of California, Irvine) and François Nosten (Shoklo Malaria Research Unit SMRU) for enlightening discussions and Steeve Ebener (Health GeoLab Collaborative HGLC) for helpful advice on spatial data management and thematic mapping.

Funding

Mahidol Oxford Tropical Medicine Research Unit (MORU) is funded by the Wellcome Trust of Great Britain (101148/Z/13/Z). Individual graduate support for CM was provided by the Rockefeller Foundation and the Faculty of Tropical Medicine, Mahidol University (Dean-MORU scholarship) and salary support by the Bill and Melinda Gates Foundation (OPP1110500). RJM is funded by the Bill and Melinda Gates Foundation (OPP1129596, OPP1110500 and CPT000390) and Asian Development Bank (TA-8656REG 47197-001). These funding bodies did not play a role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

CM SaL and RJM conceptualized and designed the study. JK AK PS and JW provided technical advice and data on malaria surveillance in Thailand. CM NE and RJM coordinated with Thai agencies. NE KH and RJM assisted in data processing. CM analysed the data and wrote the manuscript. SaL WP SY SiL and RJM supervised the work. SaL RRM VC ND AD and RJM critically reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Chris Erwin G. Mercado.

Ethics declarations

Ethics approval and consent to participate

Ethics approval for this study was obtained from the Faculty of Tropical Medicine Ethics Committee (FTM-EC), Mahidol University (MUTM 2014-061-01).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Additional file

Additional file 1.

Additional figures and tables.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mercado, C.E.G., Lawpoolsri, S., Sudathip, P. et al. Spatiotemporal epidemiology, environmental correlates, and demography of malaria in Tak Province, Thailand (2012–2015). Malar J 18, 240 (2019). https://doi.org/10.1186/s12936-019-2871-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-019-2871-2

Keywords