Geographic information systems and logistic regression for high-resolution malaria risk mapping in a rural settlement of the southern Brazilian Amazon
© de Oliveira et al.; licensee BioMed Central Ltd. 2013
Received: 21 May 2013
Accepted: 4 November 2013
Published: 15 November 2013
In Brazil, 99% of the cases of malaria are concentrated in the Amazon region, with high level of transmission. The objectives of the study were to use geographic information systems (GIS) analysis and logistic regression as a tool to identify and analyse the relative likelihood and its socio-environmental determinants of malaria infection in the Vale do Amanhecer rural settlement, Brazil.
A GIS database of georeferenced malaria cases, recorded in 2005, and multiple explanatory data layers was built, based on a multispectral Landsat 5 TM image, digital map of the settlement blocks and a SRTM digital elevation model. Satellite imagery was used to map the spatial patterns of land use and cover (LUC) and to derive spectral indices of vegetation density (NDVI) and soil/vegetation humidity (VSHI). An Euclidian distance operator was applied to measure proximity of domiciles to potential mosquito breeding habitats and gold mining areas. The malaria risk model was generated by multiple logistic regression, in which environmental factors were considered as independent variables and the number of cases, binarized by a threshold value was the dependent variable.
Out of a total of 336 cases of malaria, 133 positive slides were from inhabitants at Road 08, which corresponds to 37.60% of the notifications. The southern region of the settlement presented 276 cases and a greater number of domiciles in which more than ten cases/home were notified. From these, 102 (30.36%) cases were caused by Plasmodium falciparum and 174 (51.79%) cases by Plasmodium vivax. Malaria risk is the highest in the south of the settlement, associated with proximity to gold mining sites, intense land use, high levels of soil/vegetation humidity and low vegetation density.
Mid-resolution, remote sensing data and GIS-derived distance measures can be successfully combined with digital maps of the housing location of (non-) infected inhabitants to predict relative likelihood of disease infection through the analysis by logistic regression. Obtained findings on the relation between malaria cases and environmental factors should be applied in the future for land use planning in rural settlements in the Southern Amazon to minimize risks of disease transmission.
KeywordsMalaria Remote sensing Spatial analysis Epidemiology
Malaria is one a major public health problem and it affects more than three hundred million individuals per year. It severely impacts the African continent and affects more than one million people per year in the Amazon countries in South America. Brazil accounts for one-third of the reported malaria cases .
Malaria is one of the most serious and striking of the transmissible diseases, and it affects approximately 500 million people per year worldwide, causing more than one million deaths each year . In Brazil, 99% of the cases of malaria are concentrated in the Amazon region. The disease has been highly transmissible, maintaining high levels that are superior to those of the 1970s, when 3.9 cases per 1,000 inhabitants in the Amazon region were recorded. In 1999, 2001 and 2003, the annual parasite index (API) (number of positive slides for malaria/year per 1,000 inhabitants) of the region was of 31.9, 18.8 and 19.3 cases per 1,000 inhabitants, respectively . In 2005, 90% of the cases of malaria were registered from the recent occupation of rural areas, with activities such as manual gold mining, wood extraction and subsistence crop and cattle farming. On a regional scale, the distribution of malaria is commonly associated with environmental conditions and, mainly, with the tropical climate .
In 2005, Mato Grosso registered 9,774 cases of malaria, which corresponded to 2% of the total of cases in the Legal Amazon. When compared to 2004 (7,849 cases), there was a 38.7% increase of the cases in the State. The cities of Juruena and Rondolândia nearby in the northwest of the State deserve special attention for their high malaria incidences with APIs higher than 50 per 1.000 inhabitants . Climatic conditions guarantee the reproduction and longevity of the Anopheles mosquitoes that transmit the disease. Environmental factors, such as the periodicity of rainfall and the flood and recession of the Amazon tributaries, paired with social aspects such as the extension of the disorganized land occupation and installation of gold mining sites, are factors that have favoured and been responsible for the occurrences of malaria in the cities in the north of the State .
The city of Juruena presented 720 positive smears for malaria in 2004. This corresponded to an API of 116.8 per 1.000 inhabitants, which represented an increase of 284.9% in the incidence of positive slides when compared to the 2003 API of 41.0 per 1,000 inhabitants .
As it is part of the Amazon biome, this area originally had a predominance of semi-evergreen seasonal and evergreen forests with trees as tall as 50 m, and is located on a plateau with moderately hilly terrain in altitudes between about 205 and 295 m. The tropical lowland rainforest climate is characterized by an annual mean temperature of about 27°C and precipitations of 2,250 mm per year with two distinct periods; about 80% of the rainfall occurs between September and April, whereas monthly precipitation from May to August is mostly lower than 50 mm.
Remote sensing (RS) and geographic information systems (GIS) have proved to be an innovative and important component in studies of public health and epidemiology  and have been used for monitoring, surveillance and spatial modelling of diseases, provide examples of how Eearth observation satellites can be used in studies of ecology and prediction of malaria, and have contributed with examples for the mapping of malaria vectors, using mid-resolution RS-imagery such as Landsat ETM or SPOT [10–15].
In local-scale studies it has become common to use spatial analysis techniques, which allow the precise localization of risk areas, such as the distance between the vector’s breeding places and households, the dispersion of the vectors and case clusters [16–19].
Inasmuch as the risk of contracting malaria is related to diverse factors such as environmental alterations caused by human activities, this study aims to identify and analyse local-scale, spatial patterns of the disease in the Vale do Amanhecer settlement, an area with documented, elevated malaria incidences .
Data collection and analysis
The record of 585 malaria cases for 2005 were obtained at the Municipal Secretary of Health of Juruena. Out of the total of notifications in the Juruena municipality, 336 came from the Vale do Amanhecer land reform settlement. Notification forms recording positive slides for malaria were compared to those recorded in the national Epidemiological Surveillance Information System for malaria . Comparisons were made in order to verify data consistency as well as to complement the data when fields in the photocopied forms were illegible. The EpiInfo 3.3.2 (CDC) software was used to build the spatial database of malaria cases. Notification forms of positive slides were related to individuals domiciles, georeferenced in geographical coordinates by a GPS field survey (SAD 69 Datum). The database was then exported to the TerraView 3.1.4 (INPE) GIS software.
Spatial data processing
SPRING 4.2 (INPE) and ENVI (RSI) were applied for satellite imagery pre-processing such as band composition and calculation of spectral indices such as the normalized difference vegetation index (NDVI) and the soil and vegetation humidity index (VSHI) of the Tasseled Cap Transformation. Georeferencing and processing of vector layers, such as georeferenced field data and the settlement digital map (1:50.000 scale) obtained from INCRA (National Institute for Colonization and Agrarian Reform), visual interpretation for land-use mapping and spatial analysis procedures such as buffering were realized in the ArcGIS 9.2 (ESRI) and TerraView 3.1.4 (INPE).
Land Use and Cover (LUC) mapping was realized using a multispectral Landsat 5 TM-image provided by the Instituto Nacional de Pesquisas Espaciais (INPE) from 25/06/2005 (WRS 229/67).Visual interpretation of Landsat TM imagery by experts generally outperforms supervised classification techniques [22, 23]. Therefore the LUC map was elaborated by visual interpretation of a #3/#4/#5 colour composite, resulting in three classes coded ordinally to 1 - forest, 2 - secondary vegetation and 3 - agricultural area.
The NDVI was scaled into 256 grey levels and then three classes defined: low (0 to 100), intermediate (101 to 200) and high (> 200). The soil and vegetation humidity index (VSHI) was obtained as the third component of a Tasseled Cap Transformation , which represents the continuous distribution of the wetness of vegetation/soil during satellite passage (dry season) and then differentiated into three equally spaced classes of the dynamic range of the VSHI scene.
Slope and downslope direction maps, proxies for terrain shadowing and therefore potentially relevant for vector habitat sutiability [25, 26], were extracted from a 90 m resolution digital elevation model of the Shuttle Radar Topography Mission (SRTM). The slope map represents the continuous distribution of the terrain slope with values ranging from 0 to 16 degrees. The downslope direction map is a representation of nine classes of slope direction (eight main directions and “flat”), which were used in the logistic regression model as nominal coded.
Distance of domiciles to potential breeding sites [27–29] and mining areas [30, 31] has shown to have potential predictive power in household malaria risk mapping. Euclidian distance maps to water bodies with 100 m intervals, and to mining areas with 300 m intervals were generated using standard buffering techniques. These different intervals were chosen because of their maximum distances to the domiciles, ranging from 2,900 m at maximum in the case of most distant potential breeding site, and 6,900 m in the case of the most distant mining area.
This study is part of the malaria spatial analysis in rural land reform settlements project and was approved by the Julio Muller Hospital Ethics Committee, protocol number 326/CEP/HUJM/07 on the 9 May, 2007. This study is based on secondary data made available by the Health Sector of the Municipalities included in the investigation. The data banks do not include names or other information that allow for identifying the subjects.
Data analysis and spatial modelling
g(P i ) = link function
P i = likelihood of response for the –ith factor (or covariate)
β 0 = intercept
β i = coefficient
x i = independent variables
Exploratory variables for the analysis used in the logistic regression for malaria infection model in rural settlement in Brazilian Amazon
Number of classes
Vegetation index (NDVI)
Vegetation and soil humidity index (VSHI)
Mining area distance
In 2005, the Vale do Amanhecer settlement had 718 inhabitants, of which 394 (54.87%) were men and 324 (45.13%) were women. In the settlement, 359 cases of malaria were notified, distributed in 200 domiciles, the city of Juruena presented 720 positive smears for malaria in 2004. This corresponded to an API of 116.8 per 1.000 inhabitants, which represented an increase of 284.9% in the incidence of positive slides when compared to the 2003 API of 41.0 per 1,000 inhabitants .
Five cases were excluded because the notification forms of SIVEP-malaria did not inform their dwelling places, and another 18 because they belonged to domiciles outside the settlement limits. The 7,200 ha reserved for land owners and another 7,200 ha designated to permanent environmental reserve were considered the official limit area of the settlement.
Parameters of variables included in the logistic regression model for malaria infection in the rural settlement in Brazilian Amazon
Mining area distance
Mining area distance
Mining area distance
Mining area distance
Malaria risk and its predictive mapping in South America should be addressed under four interconnected spatial scales, which have to be considered even if a study focuses on a local scale. First, continental scale, high temperature, humid and semi-humid tropical, lowland climates favour elevated vector densities and malaria transmission in the Southern Amazon [32, 33]; second, (super) regional scale, which is little discussed in the literature yet, inside this climatic zone, extreme spatial heterogenities in malaria cases can be observed, ranging from less than 0.5 malaria incidences per 1,000 inhabitants in municipalities with consolidated land use, against up to more than 420 incidences in 2005 in municipalities with mining activities and/or land reform settlements , such as Juruena, where rural population densities are commonly higher than in areas with large crop farming activity. Inside high incidence municipalities, cases are further strongly clustered, frequently in areas with selective logging and mining activities [34, 35]. Mining areas in Mato Grosso have been found to have high rates of asymptomatic infection making disease control difficult . The assessment showed that malaria incidence in the Vale de Amanhecer settlement is more than ten times higher than in the rest of the municipality, probable favoured by an intense flux of populations with low immunity, poor housing conditions and precarious health care [36–39]. Due to the very high incidences in the settlement, logistic regression modeling results only in a reasonable overall performance, explained variance and sensitivity, if binarization is increased from one to two cases per household.
As reported by Rodrigues , the present study gives further evidence that GIS and logistic regression can be successfully applied for further local-scale zoning of high-risk areas. Greater likelihood of contracting malaria are significantly linked to areas close to gold mining sites , which present intense land use and occupation, high level in the VSHI and low NDVI values. That high incidences are related to the fine-grained spatial association between natural preserved environments of the vector habitat altering with heavily modified man-made landscape patches, characterized by the high exposition and abundant vector presence. In this context, emphasize that the positive association between malaria incidence in the Amazon with deforestation [42, 43] cannot be generalized, but may sharply differ or even be inverted as a function of observation scale and socio-environmental covariates. There is no evidence that malaria risk is elevated because of the presence of larger deforestation patches inside the settlement per se, but because these areas are characterized by elevated population fluxes and are frequently mining areas with man-made reproduction habitats favouring vector presence.
As deforestation tends to be higher near households, the output map must be interpreted as the potential risk, considering the spatial distribution of housings, which means that the geographical pattern of relative likelihood of disease infection would be altered if housing patterns were different (e g, allocation of new families). The relative importance of socio-environmental factors as determined by the logistic regression model however, is not biased as cases and non-cases are related to households and their surroundings.
Some local-scale malaria studies have reported a positive association between malaria incidences and proximity to watercourses , a variable not included as significant predictor in the logistic regression model in the present study. That these results have two reasons, mainly. First, most housing is located near watercourses, therefore, distance has a low amplitude of variation, and high distance, which hypothetically could have more non-cases, is under-represented. Second, mining areas are located along watercourses to facilitate gold extraction and include a mosaic of small man-made waterbodies. Therefore, both distance to watercourses and potential reproduction habitats have highly significant rank correlation with distance to mining areas and are therefore omitted as predictors in the final model.
Similarly, attribute lack of significant relationship of malaria incidences with DEM derived variables such as “slope” and “downslope” direction , with the spatial scale and the geographical configuration of study area. The terrain in the settlement is mostly flat or slightly undulated, so differences related to relief shadowing are low if compared to relief variations in the region on a coarser scale. Finally, the minimum mapping unit of TM imagery on the order of 0.5–1.0 ha, appropriate for thematic mapping not finer than 1:50,000 (Federal Geographic Data Committee, 1992, Wright et al. 2007, Collins and Stephens 2010), may explain the limited predictive power of land-use in the final logistic regression model.
Malaria is a focal disease, and even this small settlement area presented heterogeneity in the spatial distribution of incidences. These patterns are less related to the natural environment per se, than caused by land use, landscape modification due to human activities in the settlement and the proximity of individuals to places with elevated vector presence.
In a high malaria risk area, GIS and logistic regression could be successfully applied to predict relative likelihood of disease infection, which is positively related principally to proximity of gold mining areas and elevated nearby mining areas and, secondarily, to intense land use, lower vegetation density and higher soil humidity.
Findings on the relationship between malaria cases and environmental factors should be applied in the future for land use planning in rural settlements in the Southern Amazon to minimize risks of disease transmission.
- World Health Organization (WHO): Why tackle malaria?. 2001, Geneva: WHO, [serial on the Internet]. [cited 2001 Ago 10]. Available from: http://www.who.org/tdr
- Castro MC, Singer BH: Meio ambiente e saúde: metodologia para análise espacial da ocorrência de malária em projetos de assentamento. Revista Brasileira de Estudos Populacionais. 2007, 24: 247-262. 10.1590/S0102-30982007000200005.View Article
- Ministério da Saúde: Brasil, Secretaria de Vigilância em Saúde. 2004, Brasília: Programa Nacional de Controle da Malária - PNCM, 2
- Ministério da Saúde: Brasil, Secretaria de Vigilância em Saúde. 2005, Brasília: Programa Nacional de Controle da Malária - PNCM, 2
- Ministério da Saúde: Brasil, Secretaria de Vigilância em Saúde, Portal da Saúde. Relatório do sistema nacional de vigilância em saúde. 2006, Brasília, Acesso em novembro-2008; Disponível em http://www.saude.gov.br/relatório_snvs_mt
- Mato Grosso. Secretaria de Estado de Saúde: Boletim epidemiológico da malária. 2004, Cuiabá
- Mato Grosso. Secretaria de Estado de Saúde: Dados Básicos sobre Malária no Mato Grosso, período 1990 a 2005. 2005, Cuiabá
- INCRA – Instituto Nacional de Colonização e Reforma Agrária: Dados do Assentamento de Vale do Amanhecer – Juruena-MT. 2006, Cuiabá
- Clarke KC, Mclafferty SL, Tempalski BJ: On epidemiology and geographic information systems: a review and discussion of future directions. Emerg Infect Dis. 1996, 2: 85-92. 10.3201/eid0202.960202.PubMed CentralView ArticlePubMed
- Thomson MC, Connor SJ, Milligan PJM, Flasse SP: The ecology of malaria as seen from Earth observation satellites. Ann Trop Med Parasitol. 1996, 90: 243-264.PubMed
- Beck LR, Bradley ML, Byron LW: Remote sensing and human health: new sensors and new opportunities. Emerg Infect Dis. 2000, 6: 217-226. 10.3201/eid0603.000301.PubMed CentralView ArticlePubMed
- Rogers DJ, Randolph SE, Snow RW, Hay SI: Satellite imagery in the study and forecast of malaria. Nature. 2002, 415: 710-715. 10.1038/415710a.PubMed CentralView ArticlePubMed
- Pope KO, Rejmánková E, Savage HM, Arredondo Jimenez JI, Rodríguez MH, Roberts DR: Remote sensing of tropical wetlands for malaria control in Chiapas, Mexico. Ecol Appl. 1993, 4: 81-90.View Article
- Beck LR, Rodriguez MH, Dister SW, Rodriguez AD, Rejmankova E, Ulloa A, Meza RA, Roberts DR, Michael AS, Washino RK, Hacker C, Legters LJ: Remote sensing as a landscape epidemiologic tool to identify villages at high risk for malaria transmission. Am J Trop Med Hyg. 1994, 51: 271-280.PubMed
- Roberts DR, Paris JF, Manguin S, Harbach RE, Woodruff R, Rejmankova E, Polanco J, Wullschleger B, Legters LJ: Predictions of malaria vector distribution in Belize based on multispectral satellite data. Am J Trop Med Hyg. 1996, 54: 304-308.PubMed
- Booman M, Durrheim DN, La Grange K, Matrin C, Mabuza AM, Mbokazi FM: Using a geographical information system to plan a malaria control programme in South Africa. Bull World Health Organ. 2000, 78 (12): 1438-1444.PubMed CentralPubMed
- Carter R, Mendis KN, Roberts D: Spatial targeting of interventions against malaria. Bull World Health Organ. 2000, 78: 1401-1411.PubMed CentralPubMed
- Carvalho MS, Souza-Santos R: Análise de dados espaciais em saúde pública: métodos, problemas, perspectivas. Cad Saude Publica. 2005, 21 (2): 361-378. 10.1590/S0102-311X2005000200003.View ArticlePubMed
- Borror DJ, Delong DM: Introdução ao Estudo dos Insetos. 1988, São Paulo: Edgard Blücher
- De Oliveira EC, Santos ES, Zeilhofer P, Santos RS, Santos MA: Spatial patterns of malaria in a land reform colonization project, Juruena municipality, Mato Grosso, Brazil. Malar J. 2011, 10: 177-10.1186/1475-2875-10-177.PubMed CentralView ArticlePubMed
- Mato Grosso. Secretaria de Estado de Saúde de Mato Grosso: SIVEP-MALARIA. Acesso em novembro de 2007. Disponível em http://www.saude.gov.br/sivep_malaria
- Mas JS, Ramirez I: Comparison of land use classifications obtained by visual interpretation and digital processing. ITC Journal. 1996, 3 (4): 278-283.
- Lillesand T, Kiefer RW, Chipman J: Remote Sensing and Image Interpretation. 2007, New York: Wiley
- Crist EP, Cicone RC: Application of the Tasselled Cap Concept to Simulated Thematic Mapper Data, Photogramm. Eng. and Remote Sensing. 1994, 50: 343-352.
- Moss WJ, Hamapumbu H, Kobayashi T, Shields T, Kamanga A, Clennon J, Mharakurwa S, Thuma PE, Glass G: Use of remote sensing to identify spatial risk factors for malaria in a region of declining transmission: a cross-sectional and longitudinal community survey. Malar J. 2011, 10: 163-10.1186/1475-2875-10-163.PubMed CentralView ArticlePubMed
- Vanwambeke SO: Landscape and land cover factors influence the presence of Aedes and Anopheles larvae. J Med Entomol. 2007, 44 (1): 133-144. 10.1603/0022-2585(2007)44[133:LALCFI]2.0.CO;2.View ArticlePubMed
- Zhou S, Zhang S, Wang JJ, Zheng X, Huang F, Li WD, Zhang HW: Spatial correlation between malaria cases and water-bodies in Anopheles sinensis dominated areas of Huang-Huai plain, China. Parasites & Vectors. 2012, 5: 106-10.1186/1756-3305-5-106.View Article
- Peterson I, Borrell LN, El-Sadr W, Teklehaimanot A: Individual and household level factors associated with malaria incidence in a highland region of ethiopia: a multilevel analysis. Am J Trop Med Hyg. 2009, 80 (1): 103-111.PubMed
- Clarke SE, Bøgh C, Brown RC, Walraven GE, Thomas CJ, Lindsay SW: Risk of malaria attacks in Gambian children is greater away from malaria vector breeding sites. Trans R Soc Trop Med Hyg. 2002, 96: 499-506. 10.1016/S0035-9203(02)90419-0.View ArticlePubMed
- Moreno JE, Rubio-Palis Y, Páez E, Pérez E, Sánchez V: Abundance, biting behaviour and parous rate of anopheline mosquito species in relation to malaria incidence in gold mining areas of southern Venezuela. Med Vet Entomol. 2007, 21 (4): 339-349. 10.1111/j.1365-2915.2007.00704.x.View ArticlePubMed
- Patz JA, Daszak P, Tabor GM, Aguirre AA, Pearl M, Epstein J, Wolfe ND, Kilpatrick AM, Foufopoulos J, Molyneux D, Bradley DJ: Unhealthy landscapes: policy recommendations on land use change and infectious disease emergence. Environ Health Persp. 2004, 112 (10): 1092-10.1289/ehp.6877.View Article
- Rogers DJ, Randolph SE: The global spread of malaria in a future, warmer world. Science. 2000, 289: 1763-1766.View ArticlePubMed
- Olson SH, Gangon R, Elguero E, Durieux L, Guégan J-F, Foley JA: Links between climate and malaria in the Amazon Basin. Emerg Infect Dis. 2009,http://wwwnc.cdc.gov/eid/,
- Confalonieri UEC: Saúde na Amazônia: um modelo conceitual para a análise de paisagens e doenças. Estudos Avançados. 2005, 19 (53): 221-236.View Article
- Barbieri AF, Sawyer DO: Heterogeneity of malaria prevalence in alluvial gold mining areas in Northern Mato Grosso State, Brazil. Cad Saude Publica. 2007, 23 (12): 2878-2886. 10.1590/S0102-311X2007001200009.View ArticlePubMed
- Coura JR, Suárez-Mutis M, Ladeia-Andrade S: A new challenge for malaria control in Brazil: asymptomatic Plasmodium infection - a review. Mem Inst Oswaldo Cruz. 2006, 3: 229-237.
- Marques AC: Human migration and the spread of malaria in Brazil. Parasitol Today. 1987, 3: 166-170. 10.1016/0169-4758(87)90170-0.View Article
- Ladeia-Andrade S, Ferreira MU, Carvalho ME, Curado I, Coura JR: Age-dependent acquisition of protective immunity to malaria in riverine populations of the Amazon Basin of Brazil. Am J Trop Med Hyg. 2009, 80 (3): 452-459.PubMed
- Oliveira-Ferreira J, Lacerda MVG, Brasil P, Ladislau JLB, Tauil PL, Daniel-Ribeiro CT: Malaria in Brazil: an overview. Malar J. 2010, 9: 115-10.1186/1475-2875-9-115.PubMed CentralView ArticlePubMed
- Rodrigues AF, Escobar AL, Souza-Santos R: Análise espacial e determinação de áreas para o controle da malária no Estado de Rondônia. Rev Soc Bras Med Trop. 2008, 41 (1): 55-64. 10.1590/S0037-86822008000100011.View Article
- Duarte EC, Fontes CJF: Associação entre a produção anual de ouro em garimpo e incidência de malária em Mato Grosso-Brasil, 1995–1996. Revista Sociedade Brasileira de Medicina Tropical. 2002, 35 (6): 665-668. 10.1590/S0037-86822002000600020.View Article
- Vasconcelos CH, Novo EMLM, Donalisio MR: Uso do sensoriamento remoto para estudar a influência de alterações ambientais na distribuição da malária na Amazônia brasileira. Cad Saude Publica. 2006, 22 (3): 517-526. 10.1590/S0102-311X2006000300006.View ArticlePubMed
- Yasuoka J, Levins R: Impact of deforestation and agricultural development on anopheline ecology and malaria epidemiology. Am J Trop Med Hyg. 2007, 76 (3): 450-460.PubMed
- Kleinschmidt BL, Sharp GPY, Clarke B, Fraser CC: Use of generalized linear mixed models in the spatial analysis of small-area malaria incidence rates in KwaZulu Natal, South Africa. Am J Epidemiol. 2000, 153 (12): 13-21.
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.