Next Article in Journal
Post-Traumatic Stress Disorder in the Military Police of Rio de Janeiro: Can a Risk Profile Be Identified?
Previous Article in Journal
Six-Minute Walk Distance in Breast Cancer Survivors—A Systematic Review with Meta-Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mapping Risk of Malaria as a Function of Anthropic and Environmental Conditions in Sussundenga Village, Mozambique

1
Instituto Superior de Ciências e Educação a Distância, Beira 2102, Mozambique
2
School of Public Health, University of Minnesota, Minneapolis, MN 55455, USA
3
Direcção Distrital de Saúde de Sussundenga, Sussundenga 2207, Mozambique
4
Centro de Informação Geográfica-Faculdade de Economia da UCM, Beira 2102, Mozambique
5
Faculdade de Ciência de Saúde da UCM, Beira 2102, Mozambique
6
Faculdade de Engenharia da UCM, Chimoio 2203, Mozambique
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(5), 2568; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18052568
Submission received: 5 January 2021 / Revised: 18 February 2021 / Accepted: 20 February 2021 / Published: 5 March 2021

Abstract

:
Mozambique is a country in Southern Africa with around 30 million inhabitants. Malaria is the leading cause of mortality in the country. According to the WHO, Mozambique has the third highest number of malaria cases in the world, representing approximately 5% of the world total cases. Sussundenga District has the highest incidence in the Manica province and environmental conditions are the major contributor to malaria transmission. There is a lack of malaria risk maps to inform transmission dynamics in Sussundenga village. This study develops a malaria risk map for Sussundenga Village in Mozambique and identifies high risk areas to inform on appropriate malaria control and eradication efforts. One hundred houses were randomly sampled and tested for malaria in Sussundenga Rural Municipality. To construct the map, a spatial conceptual model was used to estimate risk areas using ten environmental and anthropic factors. Data from Worldclim, 30 × 30 Landsat images were used, and layers were produced in a raster data set. Layers between class values were compared by assigning numerical values to the classes within each layer of the map with equal rank. Data set input was classified, using diverse weights depending on their appropriateness. The reclassified data outputs were combined after reclassification. The map indicated a high risk for malaria in the northeast and southeast, that is, the neighborhoods of Nhamazara, Nhamarenza, and Unidade. The central eastern areas, that is, 25 de Junho, 1 and 2, 7 de Abril, and Chicueu presented a moderate risk. In Sussundenga village there was 92% moderate and 8% high risk. High malaria risk areas are most often located in densely populated areas and areas close to water bodies. The relevant findings of this study can inform on effective malaria interventions.

1. Background

Mozambique is a country in Southern Africa with around 30 million inhabitants. Malaria is the leading cause of mortality in the country. According to the WHO, Mozambique has the third highest cases worldwide, approximately 5% of the world cases [1,2]. In Mozambique malaria is endemic across the country, ranging from hyper-endemic in areas along the coast, meso-endemic in inland flatlands, and hypo-endemic areas in the highlands. Malaria transmission occurs all year round, peaking during the rainy season. Anopheles funestus is the most frequent, with 90% of the cases, followed by Plasmodim ovale with 9%, and Plamodium malariae with 1% of the cases [3].
Manica Province, located in the central region of Mozambique, has the fourth largest number of cases among all provinces in Mozambique. In 2019 the province recorded 821,775 malaria cases. With Manica Province, Sussundenga District had an increased number of malaria cases and in the first six months of 2020, recorded 63,526 cases accounting for around 19% of the province cases (out of 12 districts) [4].
Several factors contribute to malaria endemicity, including climate and environmental conditions such as increased temperatures and rainfall patterns, favorable places for vector development, and sociodemographic and economic conditions. Prior study in Chimoio, Mozambique [5] indicated that environmental conditions, account for 73% of malaria cases while, Global Fund [6] indicated a figure of 90%. A study in Sussundenga reported that social determinants contribute to 15% of the positive malaria cases.
The spatial distribution of transmission intensity and associated variables have become an urgent requirement, especially in endemic areas. GIS is an innovative technique to understand malaria spatiality and identify high-risk zones to develop appropriate interventions [7]. A malaria risk map does not exist for the Sussundenga village, as well as local scale maps to localize risk areas based on the distance between the vector breeding sites and host households, the vector dispersion, and case clusters [8]. The maps could also be useful long term for supporting elimination efforts. The objective of this study was to develop a malaria map of for Sussundenga Village in Mozambique, identifying zones at risk to inform malaria eradication efforts, and collect environmental, anthropic, and clinic data.

2. Methods

2.1. Study Area

Sussundenga village (Figure 1) is located in Sussundenga District, 42 km from Chimoio, the capital of Manica Province, around 45 km from the Zimbabwe border.
The rural and agrarian region is about 156.9 square kilometers and has 41,354 inhabitants. The village is administratively divided into 17 residential areas called “Bairros” (Figure 2) [9].
There are two seasons, the rainy season starts in November and ends in March, the dry season is from April to October. The average rainfall is 1067.6 mm, which varies significantly in quantity and distribution within and between years. The average highest temperatures are in October and November (30 °C), the lowest in July (12 °C) [10].
The predominant vegetation consists of deciduous forests, evergreen forests, prairie, shrubs, and savanna. The hydrographic network is comprised of 6 main rivers, Revué, Munhinga, Mussapa, Lucite, Chicueu, and Muzória with permanent flow, and several waterbodies built for irrigation. The village has one public health center and the majority of residents have low access to health services. Bairros 25 Junho 1, 2, and Chicueu have a health center in close proximity [9].

2.2. Data Collection and Analysis

Spatial, anthropic, and environmental variables were used to model the risk map. A four-step strategy was applied to develop the map as described by Ferrao et al. [11].
The factors of risk used can be seen in Figure 3 as well as the included population density (inhabitants/km2), point prevalence per Bairro (%), average temperature (°C), precipitation (mm), altitude (meters), slope (degrees), NDVI (Normalized digital index), distance of the road (km), distance to waterbodies (km), and Land usage and cover (LULC).
Ten risk factors were used and Table 1 demonstrates the risk factors of malaria, with the weights and classes of anthropic and environmental conditions.

2.2.1. Point Prevalence

One Hundred houses were randomly sampled and tested for malaria in Sussundenga Village using a malaria rapid diagnostic test (RDT), RightSign Biotest® (Biotest, Hangzhou Biotech Co, China). The field team visited communities (Bairros) in different locations in the village, based on samples of the individual at malaria risk in each Bairro. They recorded the individual test result and, each individual positive case was used to calculate the point-prevalence, that is, the proportion of individuals in a population that has the disease at a single point in time, regardless of the duration of time that the individual might have that disease. This was calculated by dividing the number of positive cases per residential area by the total population of the residential area and multiplying by 100 [12] using the following formula:
P o i n t   P r e v a l e n c e   % = P e r s o n s   m a l a r i a   p o s i t i v e B a i r r o   P o p u l a t i o n   d u r i n g   t h e   p e r i o d 100

2.2.2. Average Temperature (°C)

Long term minimum, maximum temperature as presented in Figure 4 was extracted from WorldClim-Global Climate Data (European Environmental Agency), from 1970 to 2000 [13]. The average temperature was calculated and thematic map was produced using ArcGis 10.7.1 (Esri, Redlands, CA, USA). [14].

2.2.3. Precipitation (mm)

Precipitation data as presented in Figure 4 were extracted from the Bioclim (WorldClim) from 1970 to 2000 [13], processed at Diva GIS 7.5 (Centro International da Papa & Food Agricultural Program, Lima-Peru, Rome-Italy) [15] and the thematic map produced.

2.2.4. Altitude (Meters)

A digital elevation model (DEM) [16] with a resolution of 30 × 30 m was used to estimate the altitude and thematic map.

2.2.5. Slope (Degrees)

The slope was derived from the 30 m×30 m digital elevation model [16], which was obtained from the ArcGIS and thematic map spatial analysis tool.

2.2.6. Land Cover and Land Use (LULC)

Land cover and land use data were retrieved from the most recent Landsat 8 satellite image (September 2019), Appendix A [17] and thematic map. The image was reclassified into different LULC classes using the manual training sampling technique and the maximum likelihood algorithm. The LULC classes were as follows:
  • Agricultural crop area, grass. and water body.
  • Shrubland and mosaic cover vegetation
  • Forest, bare and urban settlement areas [18].

2.2.7. Distance to the Road (km)

The Euclidean distance to the nearest road was calculated using ArcGIS, classifying a 2019 Landsat image, 30 m × 30 m and thematic map. Distances from road locations were calculated using the distance measurement function in the ArcGIS 10.7.2 software.

2.2.8. Distance to Waterbodies (km)

The distance to the nearest body of water was calculated with ArcGIS 10.7.2, classifying a 2019 Landsat image of 30 m 30 m for the water and undefined zones [19] and the thematic map.

2.2.9. Vegetation Index by Normalized Differences (NDVI)

Vegetation vigor is indicated by NDVI, and greater amounts of green vegetation in the soil indicate higher NDVI. Non-vegetation classes are generally lower than the vegetation classes in NDVI and the thematic map. The resulting expression of NDVI is as follows:
NDVI = NIR RED NIR + RED
where:
NIR is the reflectance in the near infrared band.
RED is the reflectance in the red band [20].
NDVI was taken from the Landsat image.

2.2.10. Determining Risk Factor Weights (Analytical Hierarchical Process)

Analytical Hierarchical Process (AHP), a simple way to assign relative weights for different factors through applying pairwise comparison was applied [21,22]. In pairwise comparisons each factor is evaluated according to a scale ranging from 1 to 9 (Table 2). As a result of the pairwise comparisons process, a reciprocal matrix was produced, where each factor in the matrix represented the dominance of a certain factor over another in terms of their contribution to malaria risk. Thereafter, each factor in the reciprocal matrix was divided by the sum of its column. Finally, the weight of each factor was calculated as an average across the rows [23]. The true consistency ratio was used after the pairwise matrix derivation, which was calculated by dividing the consistency index for the set judgments by the index for the corresponding random matrix. Saaty suggests that, if that ratio exceeds 0.1, the set of judgments may be too inconsistent to be reliable [24].

2.2.11. Mapping Risk of Malaria

Figure 5 represents the schematic representation of the data flow and analysis to generate a malaria risk map in the Sussundenga Village.
In this step weighted overlay analysis, a modeling method for suitability, that analyses multiclass maps based on the relative importance of each thematic layer and layers class was used [25]. Raster layers were overlayed and multiplied by their assigned weight. Raster layers were overlayed and multiplied by each raster cell’s suitability value and its layer weight to derive a suitability value (Table 2), as presented in the following formula:
S =   WiSij Wi
where:
Wi is the weight of the i-th factor map,
Sij is the i-th spatial class weight of j-th,
S is the spatial unit value in the output map [26]. As a result, a new raster surface was generated representing different levels of malaria risk based on the anthropic and environmental conditions.

2.2.12. Accuracy Assessment of the Produced Map

The intention in this step was to evaluate the accuracy of the produced risk map. Data analysis of the 2019 record books was also carried out to access location of the malaria patients as a tool to check for accuracy. The risk map was also discussed with health practitioners and malaria program officers in Sussundenga to assess accuracy.

3. Results

The average temperature ranges from 21.07 °C to 21.74 °C. Precipitation ranges from 1092.01 mm to 1124 mm. The altitude ranges from 534 m to 711 m. The slope varies from 0% to 24.32% (Figure 6).
Figure 7 shows the maps of water bodies (DTWB), distance to the road (DTR), population density (PD), malaria point prevalence (MPP), and normalized difference index (NDVI). The DTWB varied from 0 to 4116.06 m, the DTR varied from 0 to 8153.21 m, the PD varied from 1639 to 41,770 inhabitants/km2. MPP varied from 0 to 47.29% and NDVI from −0.02 to 0.49.
Table 3 represents the comparison matrix of malaria risk factors used for weighting. A value of 1 means that the comparison factors have the same weight and that they equally affect malaria occurrence. A value of five means that the factor in the column has five times the risk of malaria occurrence than the comparison in the row.
The weights of each factor used for the spatial model to produce the malaria risk map were: average temperature (22.6%). precipitation (20.0%). altitude (12.7%). distance to the water body (9.36%). slope (8.28%). land use and cover (7.69%). population (6.46%). distance from the road (4.86%). malaria prevalence (4.73%). NDVI (3.30%). The consistency ratio for the pairwise matrix was 0.08.
Figure 8 shows the maps of average temperature. precipitation. altitude. slope and LULC. In terms of temperature. 100% of the Sussundenga village area was at low risk. For precipitation. the entire village (100%) had a high risk. For altitude. 100% of the area of the village of Sussundenga was at high risk. For slope. 0.21% had low risk. 14.51% had moderate risk. and 85.28% high risk. For LULC. 30.74% of the area had a low risk. 7.73% a moderate risk. and 61.53% a high risk.
Figure 9 shows the maps of DTWB. DTR. population density. prevalence of malaria. and use and normalized difference vegetation index (NDVI). For DTWB. 28.52% of the area had low risk. 48.68% moderate risk. and 22.80% high risk. For DTR. 58.19% had low risk. 21.96% had moderate risk. and 19.85% high risk. For population density. 70.57% had a low risk. 24.13% had a moderate risk. and 5.30% a high risk. For malaria prevalence. 92% had low risk. 3% had moderate risk. and 5% high risk. For NDVI. 0.5% had low risk. 90% had moderate risk. and 9.5% high risk.
Figure 10 represents the malaria risk map. Neighborhoods of Nhamaraza. Nhamarenza. Chicueu. and Unidade have higher incidence of malaria in the northeast and southeast areas. In contrast. the central eastern part. had a lower risk of malaria. in the neighborhoods 25 de Junho 1. 25 de Junho 2. and 7 de Abril.
Table 4 represents the area under different malaria risk in Sussundenga village. Around 8% of the area is high malaria risk and 92% is moderate malaria risk in Sussundenga village.

Accuracy Check

Comparing the map results and the data extracted from the book records (Appendix B) it can be noted that only one neighborhood (14%) failed the prediction. The consistency ratio for the pairwise matrix was 0.08.

4. Discussion

There are currently no studies available on the mapping risk of malaria in Sussundenga village. In this study precipitation. slope. and LULC were related to malaria risk. Similar results have been reported in Chimoio [11]. in Zimbabwe [22]. Burundi [23]. Zambia [27]. Ghana [28]. Colombia [28,29]. and China [30]. Rainfall increases the presence of stagnant water for malaria vectors for breeding and a weekly precipitation of over 10 mm will propitiate mosquito development [31]. Sussundenga has an average weekly precipitation of 22 m. from June to September the weekly precipitation is below 10 mm [10]. which may explain the high malaria risk in Sussundenga. Flat landscape areas are more prone to accumulate water from rain which increases the malaria risk [32] and may explain the increased risk in Sussundenga.
In this study. LULC seems also to play an important role in malaria occurrence. Modified man-made landscape is conducive to malaria vectors. which could be attributed to the elevated population. [8,28].
The malaria risk map in this study. suggests that places near bodies of water and higher population density have increased malaria risk compared to the other areas. Large and small dams can increase malaria cases in villages that are located near them. since dams create additional breeding habitats for vectors and higher malaria transmission [33,34]. In Maputo. Mozambique. the risk of malaria was 6.2 times greater for individuals living less than 200 m from the breeding sites than that of individuals living 500 m or more from the breeding sites [34]. Rapid urbanization of areas within villages or on the outskirts of existing urban centers is common and often developed in without proper planning. Conditions are crowded and developers tend to dig pits to extract stones and soil for house construction. creating numerous breeding sites for mosquitoes. [8,28]. In densely populated areas. mosquitoes do not have to travel far to feed [35].
In this study. Sussundenga village had malaria risk ranging from high to moderate. Similar results were reported in Chimoio. Mozambique [11]. The consistency ratio for the pairwise matrix was 0.08. which is well below Saaty’s suggested 0.1 ratio to delineate a reliable result [24]. The map accuracy in this study was 86% and similar results were reported in Swaziland [36]. In Madagascar [37] an accuracy of 74% was reported [38].

5. Limitations of the Study

This paper was an attempt to validate a previous study carried out in Chimoio in 2018 and which would have been expected to generate a predictive model that could predict malaria prevalence without the necessity of prevalence as an input.

6. Conclusions

Sussundenga village had 8% with high risk and 92% moderate risk. High malaria risk areas are most often located in densely populated areas and areas close to water bodies. The relevant findings of this study can inform on effective malaria interventions.

Author Contributions

J.L.F. contributed to data analysis and manuscript writing. D.E. contributed to English and manuscript revisions. A.N. contributed to data provision. R.M. contributed to mapping. A.T. contributed to data collection. K.M.S. contributed to English and manuscript revision. All authors have read and agreed to the published version of the manuscript. M.F.B. contributed to writing—original draft.

Funding

This study was funded through a SEED grant from the University of Minnesota Center for Global Health and Social Responsibility (CGHSR).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki. Ethical review and approval for the study was completed by the Institutional Review Board (IRB) at the University of Minnesota [STUDY00007184] and from A Comissão Nacional de Bioética em Saúde (CNBS) at the Ministry of Health of Mozambique [IRB00002657].

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Acknowledgments

We would like to thank the Provincial and district directorates of health to grant permission to carry out this study specially Firmino Jaqueta, Serafina Benesse, Filipe Murgorgo and Elsa Trabuco. To Gabriel Viegas for images formatting.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The table presents the spectral characteristics of imaging device instruments OLI & TIRS.
BandsWavelength (µm)Resolution (m)
Band 1—Coastal aerossol0.43–0.4530
Band 2—Blue0.45–0.5130
Band 3—Green0.53–0.5930
Band 4—Red0.64–0.6730
Band 5—Near Infrared (NIR)0.85–0.8830
Band 6—SWIR 11.57–1.6530
Band 7—SWIR 22.11–2.2930
Band 8—Panchromatic0.50–0.6815
Band 9—Cirrus1.36–1.3830
Band 10—Thermal Infrared (TIRS) 110.60–11.19100
Band 11—Thermal Infrared (TIRS) 211.50–12.51100
OLI = Operational Land Imager. TIR = Thermal Infra-red Sensor; Source: https://www.usgs.gov/land-resources/nli/landsat/landsat-8?qt (accessed on 10 September 2020).

Appendix B. Malaria Origin in Sussundenga Village

NeighborhoodNumber Malaria Patients%Accuracy
Nhamarenza388314.0Ok
Nhamizara761927.4Ok
25 de Junho 2280710.1Ok
25 de Junho 121337.7Ok
Chicueu657523.6Ok
Unidade20777.5Failed
7 de Abril27359.8Ok
Total27.829100.0
Source: Sussundenga Rural Centre Malaria Data.

References

  1. Word Health Organization. World Malaria Report 2019. Available online: https://www.who.int/publications/i/item/9789241565721 (accessed on 20 September 2020).
  2. Making Strides against Malaria in Mozambique. Available online: https://www.who.int/malaria/making-strides-against-malaria-in-mozambique (accessed on 19 April 2020).
  3. Instituto Nacional de Saúde (INS). Inquérito Nacional sobre Indicadores de Malária (IIM) 2018. Available online: https://microdata.worldbank.org/index.php/catalog/3488 (accessed on 20 September 2020).
  4. DPS. Relatório do Primeiro Semestre de 2020; Direccao Provincial de Saude de Manica: Chimoio, Mozambique, 2020. [Google Scholar]
  5. Ferrão, J.L.; Mendes, J.M.; Painho, M. Modelling the influence of climate on malaria occurrence in Chimoio Municipality. Mozambique. Parasites Vectors 2017, 10, 260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Global Health Group. Malaria Elimination Country Briefings Published by the UCSF Global Health Group. Available online: https://malariaworld.org/blog/malaria-elimination-country-briefings-published-ucsf-global-health-group (accessed on 22 September 2020).
  7. Musa, G.J.; Chiang, P.H.; Sylk, T.; Bavley, R.; Keating, W.; Lakew, B.; Tsou, H.C.; Hoven, C.W. Use of GIS Mapping as a Public Health Tool-From Cholera to Cancer. Health Serv. Insights 2013. [Google Scholar] [CrossRef] [PubMed]
  8. de Oliveira, E.C.; dos Santos, E.S.; Zeilhofer, P.; Souza-Santos, R.; Atanaka-Santos, M. Geographic information systems and logistic regression for high-resolution Mapping risk of malaria in a rural settlement of the southern Brazilian Amazon. Malar. J. 2013. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Plano Estratégico do Município de Sussundenga. Public Presentation Workshop. 2018. Sussundenga. Manica. Mocambique. Unpublished work 2018.
  10. Mozambique Sussundega: Climate. Weather. Temperatures. Available online: http://www.levoyageur.net/weather-city-SUSSUNDENGA.html (accessed on 5 October 2020).
  11. Ferrao, J.L.; Niquisse, S.A.; Mendes, J.M.; Painho, M. Mapping and Modelling Malaria Risk Areas Using Climate. Socio-Demographic and Clinical Variables in Chimoio, Mozambique. Int. J. Environ. Res. Public Health 2018, 15, 795. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. CDC. Principles of Epidemiology in Public Health Practice, 3rd ed.; CDC: Atlanta, GA, USA, 2006. Available online: https://www.cdc.gov/csels/dsepd/ss1978/SS1978.pdf (accessed on 5 October 2020).
  13. BIOCLIM: WorldClim Bioclimatic Variables (2000. 20. 50. 80). Available online: https://www.edenextdata.com/?q=content/bioclim-worldclim-bioclimatic-variables-2000-20-50-80 (accessed on 10 October 2020).
  14. ArcGIS Desktop 10.7.1 Quick Start Guide. Available online: https://desktop.arcgis.com/en/arcmap/10.7/get-started/setup/arcgis-desktop-quick-start-guide.htm (accessed on 10 October 2020).
  15. DIVA-GIS. Available online: https://www.gbif.org/tool/81420/diva-gis (accessed on 15 October 2020).
  16. USCGS. Digital Elevation Models (DEMs). Available online: https://www.usgs.gov/faqs/what-are-digital-elevation-models-dems?qt-news_science_products=0#qt-news_science_products (accessed on 15 October 2020).
  17. Lansat 8. Landsat-8 Imagery in the Period 2019. Available online: https://eoportal.org/web/eoportal/satellite-missions/l/landsat-8-2019 (accessed on 15 October 2020).
  18. NOAA. National Centers for Environmental Information. 2020. Available online: https://www.ncdc.noaa.gov/ (accessed on 15 October 2020).
  19. Acharya, T.D.; Lee, D.H.; Yang, I.T.; Lee, J.K. Identification of Water Bodies in a Landsat 8 OLI Image Using a J48 Decision Tree. Sensors 2016, 16, 1075. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Rouse, J.W. Monitoring vegetation systems in the great plains with ERTS. Third Earth Resour. Technol. Satell. Symp. 1973, 1, 309–317. Available online: http://0-doi-org.brum.beds.ac.uk/citeulike-article-id:12009708 (accessed on 2 November 2020).
  21. Tatjana, A.P.; Pachemska Martin, L.; Riste, T. Analytical Hierarchical Process (AHP) Method Application in the Process of Selection and Evaluation. 2014. Available online: http://rad.ihu.edu.gr/fileadmin/labsfiles/decision_support_systems/lessons/ahp/AHP_Lesson_1.pdf (accessed on 9 September 2020).
  22. Chikodzi, D. Spatial Modelling of Malaria Risk Zones Using Environmental. Anthropogenic Variables and GeograPhical Information Systems Techniques. J. Geosci. Geomat. 2013, 1, 8–14. [Google Scholar] [CrossRef]
  23. Hassaan, M.A.; Abdrabo, M.A.; Masabarakiza, P. GIS-Based Model for Mapping Malaria Risk under Climate Change Case Study: Burundi. J. Geosci. Environ. Prot. 2017, 5, 102–117. [Google Scholar] [CrossRef] [Green Version]
  24. Coile, G. The Analitycal Hierachical Process (AHP). Available online: https://training.fws.gov/courses/references/tutorials/geospatial/CSP7306/Readings/AHP-Technique.pdf (accessed on 5 January 2021).
  25. ESRI. Available online: https://www.esri.com/about/newsroom/wp-content/uploads/2019/05/understanding-weighted-overlay.pdf (accessed on 5 October 2020).
  26. Carter, B. Geographic Information System for Geoscientists: Modelling with GIS; Pergamon: Ottawa, ON, Canada, 1991; pp. 319–470. [Google Scholar]
  27. Bridges, D.J.; Chishimba, S.; Mwenda, M.; Winters, A.M.; Slawsky, E.; Mambwe, B.; Mulube, C.; Searle, K.M.; Hakalima, A.; Mwenechanya, R.; et al. The use of spatial and genetic tools to assess Plasmodium falciparum transmission in Lusaka. Zambia between 2011 and 2015. Malar. J. 2020, 19, 20. [Google Scholar] [CrossRef] [PubMed]
  28. Amratia, P.; Psychas, P.; Abuaku, B.; Ahorlu, C.; Millar, J.; Oppong, S.; Koram, K.; Valle, D. Characterizing local-scale heterogeneity of malaria risk: A case study in Bunkpurugu-Yunyoo district in northern Ghana. Malar. J. 2019, 18, 81. [Google Scholar] [CrossRef] [PubMed]
  29. Schubert, H.; Caballero Calvo, A.; Rauchecker, M.; Rojas-Zamora, O.; Brokamp, G.; Schütt, B. Assessment of Land Cover Changes in the Hinterland of Barranquilla (Colombia) Using Landsat Imagery and Logistic Regression. Land 2018, 7, 152. [Google Scholar] [CrossRef] [Green Version]
  30. Piedrahita, S.; Altamiranda-Saavedra, M.; Correa, M.M. Spatial fine-resolution model of malaria risk for the Colombian Pacific region. Trop. Med. Int. Health 2020, 25, 1024–1031. [Google Scholar] [CrossRef] [PubMed]
  31. Zhao, X.; Thanapongtharm, W.; Lawawirojwong, S.; Wei, C.; Tang, Y.; Zhou, Y.; Sun, X.; Cui, L.; Sattabongkot, J.; Kaewkungwal, J. Malaria Risk Map Using Spatial Multi-Criteria Decision Analysis along Yunnan Border During the Pre-elimination Period. Am. J. Trop. Med. Hyg. 2020, 103, 793–809. [Google Scholar] [CrossRef] [PubMed]
  32. Ferrao, J.L. Characterising, Modelling and Mapping Malaria Occurrence and Its Mortality Trend for Precision Public Health. Repositorio Universidade Nova. Available online: https://run.unl.pt/handle/10362/64068?locale=en (accessed on 2 January 2021).
  33. Sanchez-Ribas, J.; Parra-Henao, G.; Guimaraes, A.E. Impact of dams and irrigation schemes in Anopheline (Diptera: Culicidae) bionomics and malaria epidemiology. Rev. Inst. Med. Trop. São Paulo. 2012, 54, 179–191. Available online: http://www.scielo.br/scielo.php?script=sci_arttext&pid=s0036-46652012000400001&nrm=iso (accessed on 2 September 2020). [CrossRef] [PubMed] [Green Version]
  34. Kibret, S.; Lautze, J.; McCartney, M.; Nhamo, L.; Yan, G. Malaria around large dams in Africa: Effect of environmental and transmission endemicity factors. Malar. J. 2019, 18, 303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Thompson, R.; Begtrup, K.; Cuamba, N.; Dgege, M.; Mendes, C.; Gamage-Mendes, A.; Enosse, S.M.; Barreto, J.; Sinden, R.E.; Hogh, B. The Matola malaria Project: A temporal and spatial transmission of malaria disease in a suburban area of Maputo. Mozambique. Am. J. Trop. Hyg. 1997, 57, 550–559. [Google Scholar] [CrossRef] [PubMed]
  36. Krefis, A.C.; Schwarz, N.G.; Nkrumah, B.; Acquah, S.; Loag, W.; Oldeland, J.; Sarpong, N.; Adu-Sarkodie, Y.; Ranft, U.; May, J. Spatial Analysis of Land Cover Determinants of Malaria Incidence in the Ashanti Region, Ghana. PLoS ONE 2011, 6, e17905. [Google Scholar] [CrossRef] [PubMed]
  37. Rakotoarison, H.A.; Rasamimalala, M.; Rakotondramanga, J.M.; Ramiranirina, B.; Franchard, T.; Kapesa, L.; Razafindrakoto, J.; Guis, H.; Tantely, L.M.; Girod, R.; et al. Remote Sensing and Multi-Criteria Evaluation for Malaria Risk Mapping to Support Indoor Residual Spraying Prioritization in the Central Highlands of Madagascar. Remote Sens. 2020, 12, 1585. [Google Scholar] [CrossRef]
  38. Sturrock, H.J.; Cohen, J.M.; Keil, P.; Tatem, A.J.; Le Menach, A.; Ntshalintshali, N.E.; Hsiang, M.S.; Gosling, R.D. Fine-scale Mapping risk of malaria from routine aggregated case data. Malar. J. 2014, 13, 421. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Map of Mozambique, Sussundenga Village and the study area.
Figure 1. Map of Mozambique, Sussundenga Village and the study area.
Ijerph 18 02568 g001
Figure 2. Sussundenga village map and Bairro’s. 1. Nhamatiquite, 2. Samora Machel, 3. Bapua, 4. Congrosso, 5. Chizizira, 6. 3 de Fevereiro, 7. Mussacumbira, 8. Nhamarenza, 9. Dhoa, 10. Chissamba, 11. Nhaguzue, 12. Chissanba, 13. Chipendeque, 14. Tave, 15. Muzoria, 16. Unidade, 17. 7 de Abril, 18. 25 de Junho 2, 19. 25 de Junho 1, 20. Nhamawaia, e 21. Nhamezara.
Figure 2. Sussundenga village map and Bairro’s. 1. Nhamatiquite, 2. Samora Machel, 3. Bapua, 4. Congrosso, 5. Chizizira, 6. 3 de Fevereiro, 7. Mussacumbira, 8. Nhamarenza, 9. Dhoa, 10. Chissamba, 11. Nhaguzue, 12. Chissanba, 13. Chipendeque, 14. Tave, 15. Muzoria, 16. Unidade, 17. 7 de Abril, 18. 25 de Junho 2, 19. 25 de Junho 1, 20. Nhamawaia, e 21. Nhamezara.
Ijerph 18 02568 g002
Figure 3. Representation of malaria risk factors and their flow for malaria risk map generation for Sussundenga Village. DEM/SRTM = Digital elevation model/Shuttle Radar Topographic Mission, NDVI = Normalized difference vegetation index, LULC = Land use and land cover.
Figure 3. Representation of malaria risk factors and their flow for malaria risk map generation for Sussundenga Village. DEM/SRTM = Digital elevation model/Shuttle Radar Topographic Mission, NDVI = Normalized difference vegetation index, LULC = Land use and land cover.
Ijerph 18 02568 g003
Figure 4. Long term minimum and maximum temperature and, precipitation for Sussundenga 1970/2000. Source: Worldclim (European Environmental Agency).
Figure 4. Long term minimum and maximum temperature and, precipitation for Sussundenga 1970/2000. Source: Worldclim (European Environmental Agency).
Ijerph 18 02568 g004
Figure 5. Schematic representation of the data flow and analysis to generate a malaria risk map in the Sussundenga Village.
Figure 5. Schematic representation of the data flow and analysis to generate a malaria risk map in the Sussundenga Village.
Ijerph 18 02568 g005
Figure 6. Presents the map of average temperature, precipitation, altitude, slope, and land in Sussundenga Village.
Figure 6. Presents the map of average temperature, precipitation, altitude, slope, and land in Sussundenga Village.
Ijerph 18 02568 g006
Figure 7. Maps of water bodies (DTWB). distance to the road (DTR). population density (PD). malaria point prevalence (MPP) and normalized difference index (NDVI) in Sussundenga Village.
Figure 7. Maps of water bodies (DTWB). distance to the road (DTR). population density (PD). malaria point prevalence (MPP) and normalized difference index (NDVI) in Sussundenga Village.
Ijerph 18 02568 g007
Figure 8. Risk for malaria. (A). Temperature. B) Precipitation. (C). Altitude. (D). Slope. (E). LULC for Sussundenga Village.
Figure 8. Risk for malaria. (A). Temperature. B) Precipitation. (C). Altitude. (D). Slope. (E). LULC for Sussundenga Village.
Ijerph 18 02568 g008
Figure 9. Risk for malaria. (A). Distance to waterbodies. (B). Distance to roads. (C). Pop. Density. (D). Malaria prevalence. (E). NDVI for Sussundenga Village.
Figure 9. Risk for malaria. (A). Distance to waterbodies. (B). Distance to roads. (C). Pop. Density. (D). Malaria prevalence. (E). NDVI for Sussundenga Village.
Ijerph 18 02568 g009
Figure 10. Malaria risk map for Sussundenga Village.
Figure 10. Malaria risk map for Sussundenga Village.
Ijerph 18 02568 g010
Table 1. Malaria risk factor, weights, and classes of anthropic and environmental conditions.
Table 1. Malaria risk factor, weights, and classes of anthropic and environmental conditions.
Risk FactorWeight %ClassRankRisk Degree
Average Temperature °C22.422–323High
>322Moderate
<221Low
Precipitation (mm)20.8>7003High
450–7002Moderate
<4501Low
Altitude (mm)10.4<2003High
201–5002Moderate
>5001Low
Slope (degrees)7.30–53High
5–152Moderate
>151Low
LULC8.2Agric. crop area, grass and water body.3High
Shrubland & mosaic cover vegetation2Moderate
Forest, bare and urban settlement1Low
DTWB (km) <5003High
500–15002Moderate
>15001Low
DTR (km)3.8>53High
2.5–52Moderate
<2.51Low
Pop. Density5.1>90003High
6001–90002Moderate
<60001Low
Prevalence (%)5.1>213High
14–212Moderate
<141Low
NDVI4.70.255–0.9863High
0–0.252Moderate
−0.288–01Low
LULC = Land use and land cover, DTWB = Distance do Waterbody, DTR = Distance do Road, Pop. = Population, NDVI = Normalized difference vegetation index. Source: dated from: Ferrao et al. Mapping and Modelling Malaria Risk Areas Using Climate, Socio-Demographic and Clinical Variables in Chimoio, Mozambique.
Table 2. Scale to develop the pairwise comparison matrix.
Table 2. Scale to develop the pairwise comparison matrix.
1Equal importanceTwo factors also contribute equally to the objective.
3Moderate importanceExperience and judgment slightly favor one factor in relation to the other.
5Much more importantExperience and judgment strongly favor one factor over the other
7Very importantExperience and judgment very strongly favor one over the other factor.
9Absolutely importantThe evidence favoring one over the other is the highest possible validity
2,4,6,8Intermediate valuesWhen compromise is needed
Table 3. 10 × 10 Comparison matrix of Risk Factors used in the study.
Table 3. 10 × 10 Comparison matrix of Risk Factors used in the study.
Risk FactorT AveragePrepAltSLPLULCDTWBDTRPopPrevNDVI
T average1.001.003.004.004.002.007.004.004.005.00
PP1.001.003.004.003.001.007.004.004.003.00
Alt0.330.331.003.003.001.004.002.002.003.00
SLP0.250.250.331.001.002.001.003.002.002.00
LULC0.250.330.331.001.002.002.003.001.001.00
DTWB0.501.001.000.500.501.002.003.003.002.00
DTR0.140.140.251.000.500.501.001.001.002.00
Pop den0.250.250.500.330.330.331.001.002.004.00
Prevalence0.250.250.500.501.000.331.000.501.002.00
NDVI0.200.330.330.501.000.500.500.250.501.00
T average = mean temperature. Prep = precipitation. Alt = Altitude. SLP = Slope. LULC = Land use and land cover. DTWB = Distance to waterbodies. DTR = Distance to Road. Prev = Prevalence.
Table 4. Areas under different malaria risk in Sussundenga Village.
Table 4. Areas under different malaria risk in Sussundenga Village.
ValueRisk ClassificationArea (Hectare)Percentage
1Hight244.447.59
2Moderate2972.3492.41
Total3216.78100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ferrão, J.L.; Earland, D.; Novela, A.; Mendes, R.; Ballat, M.F.; Tungaza, A.; Searle, K.M. Mapping Risk of Malaria as a Function of Anthropic and Environmental Conditions in Sussundenga Village, Mozambique. Int. J. Environ. Res. Public Health 2021, 18, 2568. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18052568

AMA Style

Ferrão JL, Earland D, Novela A, Mendes R, Ballat MF, Tungaza A, Searle KM. Mapping Risk of Malaria as a Function of Anthropic and Environmental Conditions in Sussundenga Village, Mozambique. International Journal of Environmental Research and Public Health. 2021; 18(5):2568. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18052568

Chicago/Turabian Style

Ferrão, João L., Dominique Earland, Anísio Novela, Roberto Mendes, Marcos F. Ballat, Alberto Tungaza, and Kelly M. Searle. 2021. "Mapping Risk of Malaria as a Function of Anthropic and Environmental Conditions in Sussundenga Village, Mozambique" International Journal of Environmental Research and Public Health 18, no. 5: 2568. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18052568

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop