An urban-to-rural continuum of malaria risk: new analytic approaches characterize patterns in Malawi
Malaria Journal volume 20, Article number: 418 (2021)
The urban–rural designation has been an important risk factor in infectious disease epidemiology. Many studies rely on a politically determined dichotomization of rural versus urban spaces, which fails to capture the complex mosaic of infrastructural, social and environmental factors driving risk. Such evaluation is especially important for Plasmodium transmission and malaria disease. To improve targeting of anti-malarial interventions, a continuous composite measure of urbanicity using spatially-referenced data was developed to evaluate household-level malaria risk from a house-to-house survey of children in Malawi.
Children from 7564 households from eight districts throughout Malawi were tested for presence of Plasmodium parasites through finger-prick blood sampling and slide microscopy. A survey questionnaire was administered and latitude and longitude coordinates were recorded for each household. Distances from households to features associated with high and low levels of development (health facilities, roads, rivers, lakes) and population density were used to produce a principal component analysis (PCA)-based composite measure for all centroid locations of a fine geo-spatial grid covering Malawi. Regression methods were used to test associations of the urbanicity measure against Plasmodium infection status and to predict parasitaemia risk for all locations in Malawi.
Infection probability declined with increasing urbanicity. The new urbanicity metric was more predictive than either a governmentally defined rural/urban dichotomous variable or a population density variable. One reason for this was that 23% of cells within politically defined rural areas exhibited lower risk, more like those normally associated with “urban” locations.
In addition to increasing predictive power, the new continuous urbanicity metric provided a clearer mechanistic understanding than the dichotomous urban/rural designations. Such designations often ignore urban-like, low-risk pockets within traditionally rural areas, as were found in Malawi, along with rural-like, potentially high-risk environments within urban areas. This method of characterizing urbanicity can be applied to other infectious disease processes in rapidly urbanizing contexts.
Considerable and diverse evidence demonstrates that health profiles differ between urban and rural areas across the globe, but especially in developing countries [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. Underlying these differences are diverse, context-specific factors involving the social and built environment that characterize and differentiate urban and rural spaces, and thereby affect health profiles. Studies often use dichotomous classifications of urban and rural to explain differences in disease prevalence or incidence, but such designations have been usually created for governmental administrative purposes , and may be unrelated to topographical, infrastructural and economic patterns that might serve as important determinants of disease risk. Some studies have evaluated specific disease risk factors in relation to a derived urbanicity scale , while others have created multi-country urbanicity metrics based on composite geographic, economic and political variables . Studies of health and urbanicity gradients have primarily considered chronic diseases [18,19,20], but recent research has started to examine infectious outcomes as well [21, 22]. As developing countries urbanize and potentially create new opportunities for infectious disease transmission and health , developing and testing new measures of urbanicity becomes critical.
The need for a multifactorial, continuous "urbanicity" metric that can be generalized across different settings has been recognized [22, 24], as it provides a more nuanced framework that captures the complexities of urban and rural environments and their impacts on health . Vlahov and Galea  developed a theoretical framework for studying connections between urbanicity and health that included the social environment (spaces where people interact with one another through social and market exchanges), the physical environment (nature of the physical spaces where people live and work), and health and social services (access to medical care and interventions). Lourenco  argued for recognizing a mosaic of rural and urban spaces, and for the need to consider a rural-urban continuum. Although the framework of a rural-urban continuum is not widely used in disease pattern studies, the features that contribute to a more health-relevant definition of urbanicity have been noted. Cuba, for example, understood the need to provide health care services in isolated and remote areas that were comparable in quality to those found in more affluent urban centres . This has resulted in improved health profiles in isolated areas. Similarly, rural areas in Kenya with access to telecommunications services (e.g. cell phones and mobile banking) have significantly better health profiles than other, nearby, equally rural areas that lack these services . Indeed, features associated with urbanicity may have more complex, indirect relationships to health. The building of improved transportation infrastructure in rural areas, for example, might increase the economic profile of rural communities, but at the same time lead to increased road accidents [29,30,31,32]. Thus, a comprehensive, continuous measure of urbanicity should comprise relevant factors that differentiate rural and urban spaces in a manner that is associated with disease risks or health outcomes of interest. Some generic characteristics of such a metric are likely to be similar; however, details of what is relevant will vary with the ecology and epidemiology of each disease. Ultimately, such an urbanicity measure might help improve forecasting capacity and specificity to prevention or control.
Any robust measure of urbanicity would ideally allow for both within- and among-country comparisons. Studies with a scaled measure of urbanicity typically have relied on local survey data to create this construct [33, 34]. While comprehensive, these studies lack generalizability and thus require local resources to collect the necessary data. A metric that uses readily available, spatially-referenced, Geographic Information System (GIS)-based data layers (e.g. census, environmental, economic) could be generalizable, while providing increased specificity for targeted interventions to high-risk regions.
Malaria is an example of a disease whose risks differ between rural and urban settings, and whose control strategies would benefit from a more nuanced understanding of which “rural” features might actually be protective, and which “urban” features could represent risk. Incidence of malaria is generally lower in urban compared to rural areas [25, 35, 36]. Differences in malaria incidence between urban and rural spaces are partly explained by fewer opportunities for mosquito vector reproduction in urban areas and reduced blood feeding on humans [37,38,39]. Urban habitats are considered generally less favourable for most competent Anopheles species due to the paucity of suitable breeding sites [40, 41]. In addition, urban house structures tend to restrict access of adult mosquitoes to humans at night, thus reducing transmission . However, despite historical reports indicating that malaria is a “rural” problem, this disease is considered to be an emerging threat in rapidly urbanizing areas of sub-Saharan Africa [43, 44].
These complex relationships between malaria risk and urbanicity suggest that, in addition to malaria control programs that have promoted insecticide treated nets (ITNs), improved diagnostics and vector control, recent reductions in malaria incidence and mortality might be partly attributable to rapid urbanization and development . The lines between what were formerly considered rural and urban have become increasingly blurred; stronger connections between the city and country through human movement, increased “citylike” infrastructure in rural areas, and increasing population density in peri-urban areas have rendered a dichotomous classification obsolete . Factors that characterize very large urban settings are found within smaller urban and partly rural settings (e.g., improved housing, accessible and available health care, greater population density, as well as access to transportation and markets). Conversely, characteristics associated with rural settings, such as small-scale cropping, standing water, poor housing, and inadequate health services, can also be found in some urban areas. A more thorough understanding of how specific features of urbanization could influence Plasmodium transmission across the urban-to-rural landscape gradient should improve the effectiveness and efficiency of targeted interventions.
To develop and apply a new urbanicity measure relevant to countrywide malaria risk in Malawi, the Lourenco  and the Vlahov/Galea  frameworks were integrated by creating a continuous, composite, scaled measure of urbanicity that comprises aspects of the social environment, the physical environment, as well as health and social services access. Using this composite urbanicity measure, the manner in which rural-like environments exist within, and proximate to, urban areas was analysed, as well as how urban-like contexts have emerged within otherwise rural settings. Malaria, a disease known to vary between urban and rural contexts, was evaluated for whether there are consistent and predictable patterns of risk along an urbanicity gradient.
Household survey and child infection data
The survey involved 7564 households in eight Districts throughout Malawi, conducted during April/May 2007 . Within each District, 16 to 30 Enumeration Areas (EAs) were chosen at random. Within each EA, 10 to 81 households were randomly selected. The geographic location of each sampled house was recorded using a GPS unit. After obtaining informed consent from the head of household, demographic information on all household members was collected using a standard questionnaire. In addition, one child between 6 and 59 months of age was identified for inclusion in the survey. Although household and individual malaria risk data were gathered from all households, only 4684 households where eligible children were present were selected for testing for Plasmodium infection. In households with more than one eligible child, only one was selected for malaria testing. A finger-prick blood sample was taken and later examined by a trained microscopist for presence of any species of Plasmodium. There was no attempt to identify the Plasmodium species, but a majority of infections in Malawi are known to be P. falciparum . See Fig. 2 for locations of included households and parasitaemia status of selected children.
Information on household-level material assets collected by project enumerators included: type of house construction, water and light sources, type of toilet facilities, and presence of livestock, electronic goods and vehicles (bicycles, motorcycles and cars). All information was entered into a Microsoft Access database. Principal Components Analysis (PCA) was applied to these household-level material items. PCA  is a data dimension reduction technique that, given a set of p, possibly correlated, variables yields p new, uncorrelated, and mutually orthogonal variables, called Principal Components (PCs). Each component is defined as a linear combination of the original p variables, with the coefficients of the original p variables in the linear combination called factor loadings. The PCs are ordered so that the first PC accounts for the largest amount of variation in the data, and the subsequent components explain a lesser amount. Dimension reduction via PCA is achieved by replacing the original p variables with an appropriate number r, where r ≤ p, of PCs that explain a significant percentage of the total variation in the data. The first PC was selected as a composite measure of material wealth and socio-economic status (SES) and divided into quintiles. The SES quintiles, representing SES "classes" for each household were then added to the database .
Data sources and urbanicity measure creation
To create a measure of urbanicity that does not require detailed, on-the-ground surveys, easily obtained, publicly available data were analysed, based on sources that could allow for cross-country comparisons. Six environmental, social, or infrastructural features commonly considered relevant to malaria risk, and that often differentiate urban and rural areas, were considered: population density, transportation infrastructure, and location of health services along with elevation and proximity to surface water (rivers and lakes). These components relate to the social environment, to the social/health services associated with the urban-to-rural continuum of interest and to environmental factors that characterize rural regions in this context . Increased population, proximity to roads and health services all characterize development, access to areas of economic activity and the ability to receive government services. Low elevation and proximity to water characterize extremely rural and undeveloped regions in this context such as areas along Lake Malawi or in the low lying regions along the border with Zambia .
Spatially-referenced data for elevation, water bodies, and transportation networks in Malawi were downloaded from the website of DIVA-GIS . Mean population within a 1 km buffer around each household location was extracted from a 30 m resolution raster-based composite of census and remotely sensed data on locations of human settlements obtained from the WorldPop Project . These environmental and population data were integrated into a GIS database using ArcGIS  (Fig. 1). In addition, both public and private health facility locations in Malawi, that had been comprehensively surveyed in 2003 by the Japan International Cooperation Agency (JICA), were obtained. The geo-coordinates for each health facility, as well as the health facility type, ownership and funding source, were included in the database .
For point features (e.g. health facilities), point-to-point distances were calculated. For linear features (e.g. roads), distance to the nearest point on the closest line was calculated. Computing distances using the Euclidean distance formula provides measures that some research has shown to be reasonable proxies for distances generated by other, more complex methods that measure actual travel routes [55, 56]. However, other studies have suggested that measures of access to health services based on Euclidean distances can overestimate access for a significant percentage of a given population compared with more complex metrics (e.g. road networks) due to underestimation of true travel distances . Though Euclidean distances tend to underestimate true travel distances, they are correlated with true travel distances [58, 59]. Given this finding, and lacking information on actual travel routes from households to health facilities, Euclidean distances were used as a proxy for travel times required to receive health services.
From the human population data layer, 1 km buffers around each point of interest were created, and the mean population within each buffer was extracted and recorded. To construct a measure of urbanicity for the entire country of Malawi, a 1 km grid of points covering Malawi was constructed. Distances from a household, or any other point of interest, to any other point, line, or polygon feature were calculated, and raster values were extracted based on GPS location.
To derive the composite measure of urbanicity, PCA was applied to the six gridded variables that were considered relevant to the multi-dimensional concept of urbanicity, namely, distance to health facility, road, lake and river, as well as population and elevation, all appropriately normalized. Each household was assigned the values of the PCs relative to the grid centroid closest to its location.
To investigate specific patterns of association between features of urbanicity and Plasmodium infection, controlling for potential confounders, multivariate logistic regression models were developed. To assess whether the proposed composite measure of urbanicity provided a better explanation of Plasmodium infection than other readily available indicators of urbanicity, such as population density or a binary indicator of urban vs. rural areas obtained from government statistics, three logistic regression models were fitted. The three models regressed the logit of observed malaria occurrence on: (i) the composite urbanicity measure defined by the first and second Principal Components; (ii) population density; and (iii) a binary urban/rural indicator. The three models were compared based on their Akaike’s Information Criteria (AIC) value , recognizing that the model with the lowest AIC is the model that best fits the data and better explains the variability in the observed occurrence of infection. In addition, receiver operating characteristic (ROC) curves  were developed to compare model performance between the continuous measure produced in this research and the dichotomous measure. All analyses were performed using R (ver. 3.5.1) statistical software .
Plasmodium infection survey
Of 4684 children tested for Plasmodium infection, 966 (20.6%) were found to be slide-positive (Fig. 2). More than half (57%) of these children reported sleeping under an ITN the previous night (Table 1). The mean age of children tested for infection was 17.1 months (Table 1). Children who were tested were younger and slightly poorer compared to those who were not tested (Additional file 1: Table S1).
Infection and geographic features
Children with Plasmodium infection lived further from the nearest health facility and closer to a lake body compared to those who tested negative (5.1 km vs. 3.8 km and 13.1 km vs. 27.7 km, p ≤ 0.001, respectively). Infected children, however, did not live closer to a road (2.6 km vs. 2.7 km, p = 0.6) nor to a river (2.3 km vs. 2.3 km, p = 0.9). Infected children resided in significantly less densely populated areas (z = 1.12 vs. z = 1.64, p ≤ 0.001), and at a lower elevation (z = -0.57 vs. z = − 0.28, p ≤ 0.001) than children who were not infected (Table 1).
Urbanicity measure and government classification
Results of the PCA showed that all six urbanicity variables were strongly represented in the first and second PCs (Table 2). In particular, the six variables tended to fall into two distinct groups. The first PC weighted more heavily those aspects of urbanicity that are related to human and infrastructural factors, such as population density, roads and health facilities. On the other hand, the second PC is a combination of variables involving biogeographic aspects of urbanicity, e.g. elevation, proximity to lakes, and proximity to rivers and streams, which should be important to vector mosquito abundance. The first two PCs together account for 49% of the total variation in the data and represent different aspects of urbanicity. Because these two PCs are both dominant and more interpretable than the remaining four, a composite measure of urbanicity was defined as the first two PCs. Figure 3 provides maps of how this urbanicity measure, representing the sum of the first two principal components, is geographically distributed throughout Malawi.
This PCA-derived composite urbanicity measure was compared with official government urban and rural designations to assess whether the new measure identified “urban-like” spaces in areas designated as “rural” by the official classification, and conversely, whether there were “rural-like” spaces in areas officially designated as “urban.” Each household location was re-classified as “rural” if the urbanicity measure was below the median, or as "urban" if above the median. This new, PCA-derived, urban–rural re-classification was then compared with the official governmental designation of urban or rural to evaluate which household locations were classified differently. Nearly half (45.6%) of the household locations defined as rural by the Malawian government were classified as urban by the PCA-derived urbanicity measure. Likewise, there was heterogeneity within the survey areas that were officially defined as urban, with 4.4% of study households (all within the two largest cities) being classified as rural according to the urbanicity-measure driven classification.
Relationships of urbanicity to Plasmodium infection
To assess the predictive utility of the PCA-derived composite measure, it was compared with the dichotomous urban–rural classification in the Malawi census data using receiver operating characteristic (ROC) curves (Fig. 4). For the government dichotomous classification, the curve mostly rests on the diagonal, while for the PCA-derived composite, the curve extends closer to the upper left corner. This indicates that the PCA composite classification is a superior predictor of Plasmodium infection compared to the official census-based measure.
Logistic regression models indicated that ITN use was highly predictive of reduced infection risk, and there was no difference by gender (Table 3). Using the poorest SES quintile as a reference, SES was inversely related to Plasmodium parasitaemia, with the wealthier groups experiencing the lowest likelihood of testing positive. In other univariate logistic regression models, age, distance to nearest health facility, distance to lake, and population density were all significantly associated with parasitaemia, but distance to river and road were not (Table 3).
To evaluate these relationships further, all covariates were included in a single multivariate logistic regression model (Table 3). Again, gender and distance to river were not significantly associated with parasitaemia risk when controlling for all the other covariates. ITN use was still protective against infection, although at a reduced magnitude. After accounting for other spatial indicators, the association between SES and infection disappeared. Older age, increased distance to health facilities and roads, decreased distance to lake, and increased population density were all associated with an elevated risk. A model including only the composite urban–rural measure indicated that increased levels of urbanicity were associated with decreased risks of being positive for Plasmodium infection (Table 3). Based on AIC values for logistic regression models with binary outcomes of parasitaemia status, the composite urbanicity measure performed better (AIC 4347) than one with a single predictor for population density (AIC 4758) and also better than one with only the urban/rural indicator variable (AIC 4760.)
Malaria risks in Malawi are heterogeneous within politically defined urban and rural areas. For example, 23% of rural survey locations exhibited low risk characteristics of an urban designation (based on the lower 25th percentile probability of parasitaemia). On the other hand, 4% of the urban survey locations exhibited high risk characteristics of a rural designation (based on the upper 25th percentile of the probability of parasitaemia.).
Predicted malaria risk based on the composite measure
Using the composite measure of urbanicity and logistic regression, infection risk was calculated at each grid point over the entire map of Malawi (Fig. 5). Although areas along the lake were “high risk,” low lying areas of low population density such as those in and around the wildlife parks also had elevated levels of malaria risk. However, the prediction errors in these areas, where data are lacking, were very high. Finally, the estimated infection risk for households based on the PCA-derived composite measure was compared with the official census-based urban/rural designation of household location (Fig. 6). While the mean composite measure risk for rural areas was higher than that for urban areas, there was considerable overlap between the two, particularly at the second quartile. Many outlier household location points designated as “urban” had risk equivalent to points at the high end of “rural” risk. Indeed, 45 household locations designated as “urban” had Plasmodium infection risk equivalent to that of the highest quantile of risk among “rural” areas. Likewise, 295 households designated as “rural” had infection risk equivalent to that of the lowest risk among households designated as “urban.”
The analyses of this study applied different methods that created a more informative and relevant geosocial classification of malaria risk along an "urban-to-rural" continuum, indicate that there is a graded relationship with urbanicity. Indeed, the results suggest the value of creating a more complex, nuanced characterization of urbanicity when assessing Plasmodium transmission risk than what a classical binary classification offers. There are a number of special features in the results that are worth noting.
The analyses found that while Plasmodium parasitaemia status was strongly associated with distance to lakes, this was not so with proximity to rivers. Lakefronts may provide breeding opportunities for Anopheles, particularly Anopheles funestus [63, 64]. In addition, lake water is generally more seasonally stable in volume than that of rivers. However, the water layers used for this study did not include information on river width, depth or degree of water flow. A more comprehensive map of waterways for Malawi could yield more detailed results. Different river widths create environments that are more or less amenable to Anopheles species-specific reproduction, and might include conditions that are similar to lake and other still water bodies.
Other community-level factors such as greater distances to health services and roads were also found to be important risk factors (Table 3). Appreciating these community contexts, and thereby identifying specific areas at higher potential risk, will be important to developing more focused prevention. This process will lead to a better understanding of the causal mechanisms underlying risk and may unveil the manner in which closer proximity to certain features of more urbanized environments mitigates against Plasmodium infection in children.
In addition to the community-based urbanicity metric, various household-level factors, such as greater material wealth (a measure of SES) and the use of ITNs, were also found to be associated with reduced malaria risk in children. The analysis suggests that higher SES, greater ITN use, and living in more urbanized contexts each acted independently to decrease malaria risk. Thereby, the important household-level causal pathways appear to be more than house construction and condition, as is often posited. Considerably more work is needed to untangle the relationships between household- and community-level factors, and their relationships with malaria. The issue of urban–rural classification, a community-level construct, will become ever more salient as sub-Saharan African (SSA) countries continue to develop, and urbanization continues to expand.
This study and others have demonstrated the limitations of using a politically defined urban–rural construct. Within a politically defined urban or rural environment, there exists heterogeneity with respect to geographic, economic and socio-political risk factors. In general, the urban designations in use are not consistently defined across the countries of SSA. As this study in Malawi has shown, where fewer areas are classified as “urban”, people might be considered residents of rural areas even though their communities share indicators of urban settings. The reverse is also true. As other researchers have noted, independent of population density or remoteness, some “rural” communities experience better standards of living, increased economic opportunities, and city-like services such as health care and schools . Inversely, small-area environmental conditions (family gardens, lowland springs, river edges) in urban settings might create highly suitable breeding sites for Anopheles vectors of malaria, much like what is typically thought of as rural.
The limitations of using a dichotomous definition of urbanicity have practical implications. When providing anti-malaria interventions, some “rural” communities may have areas with low disease risk, yet receive extensive prevention support, leading to a waste of resources that could be better targeted elsewhere. Similarly, rural-like areas near large and dense human settlements may be overlooked during anti-malaria interventions. Uneven geographic distribution of malaria risk and inadequate knowledge of the locations of transmission foci have been recognized as challenges to the spatial targeting of malaria control  in settings of heterogeneous transmission across geographic locations . Results from the present study confirm that the transmission profile throughout a holoendemic country like Malawi varies considerably over space. These same analyses offer measurable and easily obtained markers that could potentially be used to aid in efficient targeting of anti-malaria interventions.
Another important aspect of the urban–rural continuum involves the level of connectivity between more rural and more developed contexts. In Malawi, however, ease of transport between areas along the developmental gradient is challenging. One study showed that people in extremely rural areas must rely on animal-drawn carts to receive basic medical services and must spend money to reach better secondary and tertiary care facilities located in larger cities and towns . Transportation expenses have been shown to be a major barrier to receiving prompt care for serious health conditions in other African contexts . Travel distance has been shown to be an obstacle to obtaining HIV care in Malawi . While the present research used a simple measure of access to roads and health facilities as a proxy for spatial access to areas that provide economic opportunities and government services, future efforts might employ more detailed measures that take into account different types of transport infrastructure or travel times.
To further refine the use of a continuous urbanicity metric, at least three major areas need further consideration. First, although this study created a malaria-relevant urbanicity measure based on data that are readily available for any developing country, the lack of more specific markers of social and ecological contexts and economic development hampered the analysis. Data on the locations of schools or markets, for example, might have helped to make the results more focally accurate. Given the experiences and knowledge of the authors in Malawi, however, this analysis suggests that the composite measure is a reasonable representation of the gradient of urbanicity in Malawi. Second, associations with river and lake locations may not to be generalizable for cross-country comparisons. With some exceptions, colonial powers tended to establish SSA cities far from swamps and fresh water bodies. Residents of rural areas, on the other hand, often live near surface water. Even within cities, wealthier areas tend to be located at higher elevations and away from areas where water usually collects. These SSA generalities stand in contrast to European contexts which tended to favor locations amenable to trade by water routes. Third, the cross-sectional malaria data that were used may ignore important temporal effects of seasonal malaria. Seasonal transmission of Plasmodium has been demonstrated in Malawi in other studies [69, 70], although transmission does occur throughout the year. Data used in the present study were collected to assess infection prevalence toward the end of peak transmission, but seasonal patterns of transmission may also differ by urbanicity and developmental context. Regardless, seasonal variation cannot confound the relationship between urbanicity and parasitaemia, because the urbanicity measure does not vary in time in this analysis. It is recognized, however, that the data were collected in a single year and may not be representative of other years.
Policy makers need more accurate classification of urban and rural spaces to make better use of limited intervention resources . Control and prevention strategies that inappropriately target all “rural” areas may unnecessarily expend costly resources where they are not needed, while underserving truly remote and disconnected areas that face crushing malaria incidence and mortality [36, 43, 71,72,73,74,75,76,77]. The complement to this involves “urban” areas that are ignored because of a governmental designation. Malaria control policies are likely to be more cost-effective if governments encouraged a more scientific, context and disease-specific analysis of the complex urban–rural continuum. The analysis of measures associated with urbanicity in this Malawi study demonstrated a rural–urban gradient of associated malaria risk, but with urban-like pockets in areas traditionally classified as rural, and vice versa. Additionally, malaria risk showed a graded association with levels of the new urbanicity metric that were developed. Governments might encourage use of such an approach in developing policies to more effectively target anti-malaria interventions to populations with scarce resources.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. All code used for analysis is available in the kambanane/Urbanicity repository (URL: https://github.com/kambanane/Urbanicity).
Akaike’s Information Criterion
Geographic Information System
Insecticide treated net
Principal Components Analysis
Receiver Operating Characteristics
Abdesslam B. Evolution of rural-urban health gaps in Morocco: 1992–2011. BMC Res Notes. 2012;5:381.
Fotso J. Urban-rural differentials in child malnutrition: trends and socioeconomic correlates in sub-Saharan Africa. Health Place. 2007;13:205–23.
Fotso JC. Child health inequities in developing countries: differences across urban and rural areas. Int J Equity Health. 2006;5:9.
Fox K, Heaton TB. Child nutritional status by rural/urban residence: a cross-national analysis. J Rural Health. 2012;28:380–91.
Govoetchan R, Gnanguenon V, Azond´ekon R, Agossa RF, Sovi A, Ok´e-Agbo F et al. Evidence for perennial malaria in rural and urban areas under the Sudanian climate of Kandi, Northeastern Benin. Parasit Vectors. 2014;7:79.
Harris B, Goudge J, Ataguba JE, McIntyre D, Nxumalo N, Jikwana S, et al. Inequities in access to health care in South Africa. J Public Health Policy. 2011;32(Suppl 1):S102–23.
Johnson GZ. Health conditions in rural and urban areas of developing countries. Popul Stud. 1964;17:293.
Liu H, Fang H, Zhao Z. Urban-rural disparities of child health and nutritional status in China from 1989 to 2006. Econ Hum Biol. 2013;11:294.
Matthews Z, Channon A, Neal S, Osrin D, Madise N, Stones W. Examining the ”urban advantage” in maternal health care in developing countries. PLoS Med. 2010;7:e1000327.
Molyneux C, Harpham T. Urban health in developing countries: a review. Progr Development Studies. 2001;1:113.
Nandita S, Abhishek S, Domantas J, Faujdar R. Explaining the rural-urban gap in infant mortality in India. Demographic Res. 2013;29:473.
Okeke TA, Okeibunor JC. Rural–urban differences in health-seeking for the treatment of childhood malaria in south-east Nigeria. Health Policy. 2010;95:62.
Menashe-Oren A, Stecklov G. Urban-rural disparities in adult mortality in sub-Saharan Africa. Demographic Res. 2018;39:136–76.
Beatriz ED, Molnar BE, Griffith JL, Salhi C. Urban-rural disparity and urban population growth: a multilevel analysis of under-5 mortality in 30 sub-Saharan African countries. Health Place. 2018;52:196–204.
Vlahov D, Galea S. Urbanization, urbanicity, and health. J Urban Health. 2002;79(4 Suppl 1):S1–12.
Dahly DL, Adair LS. Quantifying the urban environment: a scale measure of urbanicity outperforms the urban–rural dichotomy. Soc Sci Med. 2007;64:1407–19.
Novak NL, Allender S, Scarborough P, West D. The development and validation of an urbanicity scale in a multi-country study. BMC Public Health. 2012;12:530.
Allender S, Wickramasinghe K, Goldacre M, Matthews D, Katulanda P. Quantifying urbanization as a risk factor for noncommunicable disease. J Urban Health. 2011;88:906–18.
Popkin BM, Gordon-Larsen P. The nutrition transition: worldwide obesity dynamics and their determinants. Int J Obes Relat Metab Disord. 2004;28:S2–9.
Hermosillo-Gallardo ME, Jago R, Sebire SJ. Association between urbanicity and physical activity in Mexican adolescents: the use of a composite urbanicity measure. PLoS One. 2018;13:e0204739.
Dal’Asta AP, Lana RM, Amaral S, Codeco CT, Monteiro AMV. The urban gradient in malaria-endemic municipalities in Acre: revisiting the role of locality. Int J Environ Res Public Health. 2018;15:1254.
Dear NF, Kadangwe C, Mzilahowa T, Bauleni A, Mathanga DP, Duster C, et al. Household-level and surrounding peri-domestic environmental characteristics associated with malaria vectors Anopheles arabiensis and Anopheles funestus along an urban–rural continuum in Blantyre. Malawi Malar J. 2018;17:229.
Alirol E, Getaz L, Stoll B, Chappuis F, Loutan L. Urbanisation and infectious diseases in a globalised world. Lancet Infect Dis. 2011;11:131–41.
Wilson ML, Krogstad DJ, Arinaitwe E, Arevalo-Herrera M, Chery L, Ferreira MU, et al. Urban malaria: understanding its epidemiology, ecology, and transmission across seven diverse ICEMR network sites. Am J Trop Med Hyg. 2015;93(3 Suppl):110–23.
Robert V, Macintyre K, Keating J, Trape JF, Duchemin JB, Warren M, et al. Malaria transmission in urban sub-Saharan Africa. Am J Trop Med Hyg. 2003;68:169–76.
Lourenco AE. The meaning of ’rural’ in rural health: a review and case study from Brazil. Glob Public Health. 2012;7:1–13.
Spiegel JM, Yassi A. Lessons from the margins of globalization: appreciating the Cuban health paradox. J Public Health Policy. 2004;25:85.
Jack W, Suri T. Risk sharing and transactions costs: evidence from Kenya’s mobile money revolution. Am Econ Rev. 2014;104:183.
Khanbhai M, Lutomia MBL. Motorcycle accident injuries seen at Kakamega Provincial Hospital in Kenya. East Central Afr J Surg. 2012;17:43–6.
Lamont M. Speed governors: road safety and infrastructural overload in post-colonial Kenya, c. 1963–2013. Africa. 2013;83:367–84.
Mogaka EO, Ng’ang’a Z, Oundo J, Omolo J, Luman E. Factors associated with severity of road traffic injuries, Thika, Kenya. Pan Afr Med J. 2011;8:20.
Otieno T, Woodfield JC, Bird P, Hill AG. Trauma in rural Kenya. Injury. 2004;35:1228.
Jones-Smith J, Popkin B. Understanding community context and adult health changes in China: development of an urbanicity scale. Soc Sci Med. 2010;71:1436–46.
Cyril S, Oldroyd J, Renzaho A. Urbanisation, urbanicity, and health: a systematic review of the reliability and validity of urbanicity scales. BMC Public Health. 2013;13:513.
Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW. Urbanization, malaria transmission and disease burden in Africa. Nat Rev Microbiol. 2005;3:81–90.
Keiser J, Utzinger J, Caldas de Castro M, Smith T, Tanner M, Singer B. Urbanization in sub-Saharan Africa and implication for malaria control. Am J Trop Med Hyg. 2004;71:118–27.
Kelly-Hope LA, McKenzie FE. The multiplicity of malaria transmission: a review of entomological inoculation rate measurements and methods across sub-Saharan Africa. Malar J. 2009;8:19.
Killeen GF, McKenzie FE, Foy BD, Schieffelin C, Billingsley PF, Beier JC. A simplified model for predicting malaria entomologic inoculation rates based on entomologic and parasitologic parameters relevant to control. Am J Trop Med Hyg. 2000;62:535–44.
Smith DL, McKenzie FE. Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar J. 2004;3:13.
Trape JF. The public health impact of chloroquine resistance in Africa. Am J Trop Med Hyg. 2001;64(1–2 Suppl):12.
Lindsay SW, Campbell H, Adiamah JH, Greenwood AM, Bangali JE, Greenwood BM. Malaria in a peri-urban area of The Gambia. Ann Trop Med Parasitol. 1990;84:553–62.
Lindsay SW, Emerson PM, Charlwood JD. Reducing malaria by mosquito-proofing houses. Trends Parasitol. 2002;18:510–4.
Alemu A, Tsegaye W, Golassa L, Abebe G. Urban malaria and associated risk factors in Jimma town, south-west Ethiopia. Malar J. 2011;10:173.
Mzilahowa T, Luka-Banda M, Uzalili V, Mathanga DP, Campbell CH, Mukaka M, et al. Risk factors for Anopheles mosquitoes in rural and urban areas of Blantyre District. Southern Malawi Malawi Med J. 2016;28:151–8.
Tatem AJ, Gething PW, Smith DL, Hay SI. Urbanization and the global malaria recession. Malar J. 2013;12:133.
Larson PS, Mathanga DP, Campbell CH, Wilson ML. Distance to health services influences insecticide-treated net possession and use among six to 59 month-old children in Malawi. Malar J. 2012;11:18.
WHO. World Malaria Report 2008. Geneva: World Health Organization; 2008.
Dommenget D, Latif M. A cautionary note on the interpretation of EOFs. J Climate. 2002;15:216–25.
Filmer D, Pritchett LH. Estimating wealth effects without expenditure data—or tears: an application to educational enrollments in states of India. Demography. 2001;38:115–32.
Sumanik-Leary J, Org M, Little M, Persson J, Kalonga C, Eales A, et al. Market assessment for locally manufactured PV-wind hybrid systems in Malawi. University of Strathclyde; 2016.
Hijmans R, Guarino L, Bussink C, Mathur P, Cruz M, Barrentes. A geographic information system for the analysis of species distribution data. Bioinformatics. 2004;19.
Tatem AJ, Gething P, Weiss D, Bhatt S, Pezzulo C. Pilot high resolution poverty maps. University of Southampton, 2013.
ArcGIS Desktop: Release 10.
Carter R, Mendis K, Roberts D. Spatial targeting of interventions against malaria. Bull World Health Organ. 2000;78:1401–11.
Nesbitt RC, Gabrysch S, Laub A, Soremekun S, Manu A, Kirkwood BR, et al. Methods to measure potential spatial access to delivery care in low- and middle-income countries: a case study in rural Ghana. Int J Health Geogr. 2014;13:25.
Jones SG, Ashby AJ, Momin SR, Naidoo A. Spatial Implications associated with using euclidean distance measurements and geographic centroid imputation in health care research. Health Serv Res. 2010;45:316–27.
Noor AM, Amin AA, Gething PW, Atkinson PM, Hay SI, Snow RW. Modelling distances travelled to government health services in Kenya. Trop Med Int Health. 2006;11:188–96.
Shahid R, Bertazzon S, Knudtson M, Ghali W. Comparison of distance measures in spatial analytical modeling for health service planning. BMC Health Serv Res. 2009;9:200.
Boscoe FP, Henry KA, Zdeb MS. A nationwide comparison of driving distance versus straight-line distance to hospitals. Prof Geogr. 2012. https://doi.org/10.1080/00330124.2011.583586.
Akaike H. Information theory and an extension of the maximum likelihood principle. New York: Springer; 1998. p. 199–213.
Fawcett T. An introduction to ROC analysis. Pattern Recogn Lett. 2006;27:861–74.
R Core Team; R Foundation for statistical computing. R: A Language and environment for statistical computing. 2020.
Mendis C, Jacobsen JL, Gamage-Mendis A, Bule E, Dgedge M, Thompson R, et al. Anopheles arabiensis and An. funestus are equally important vectors of malaria in Matola coastal suburb of Maputo, southern Mozambique. Med Vet Entomol. 2000;14:171–80.
Vezenegho SB, Chiphwanya J, Hunt RH, Coetzee M, Bass C, Koekemoer LL. Characterization of the Anopheles funestus group, including Anopheles funestus-like, from Northern Malawi. Trans R Soc Trop Med Hyg. 2013;107:753–62.
Kreuels B, Kobbe R, Adjei S, Kreuzberg C, von Reden C, B¨ater K et al. Spatial variation of malaria incidence in young children from a geographically homogeneous area with high endemicity. J Infect Dis. 2008;197:85–93.
Varela C, Young S, Mkandawire N, Groen RS, Banza L, Viste A. Transportation barriers to access health care for surgical conditions in Malawi: a cross sectional nationwide household survey. BMC Public Health. 2019;19:264.
Eide AH, Mannan H, Khogali M, Van Rooy G, Swartz L, Munthali A, et al. Perceived barriers for accessing health services among individuals with disability in four African countries. PLoS One. 2015;10:e0125915.
Palk L, Okano JT, Dullie L, Blower S. Travel time to health-care facilities, mode of transportation, and HIV elimination in Malawi: a geospatial modelling analysis. Lancet Glob Health. 2020;8:e1555–64.
Lowe R, Chirombo J, Tompkins AM. Relative importance of climatic, geographic and socio-economic determinants of malaria in Malawi. Malar J. 2013;12:416.
Roca-Feltrer A, Kwizombe CJ, Sanjoaquin MA, Sesay SS, Faragher B, Harrison J, et al. Lack of decline in childhood malaria, Malawi, 2001–2010. Emerg Infect Dis. 2012;18:272–8.
De Silva PM, Marshall JM. Factors contributing to urban malaria transmission in sub-Saharan Africa: a systematic review. J Trop Med. 2012;2012:819563.
Nduka F, Kalu K. A comparative study of malaria parasitemia in urban and rural areas of South East Nigeria. Trop Med Int Health. 2011;16:129–30.
Rakotomanana F, Ratovonjato J, Randremanana R, Randrianasolo L, Raherinjafy R, Rudant JP, et al. Geographical and environmental approaches to urban malaria in Antananarivo (Madagascar). BMC Infect Dis. 2010;10:173.
Ngom R, Siegmund A. Urban malaria in Africa: an environmental and socio-economic modelling approach for Yaound´e. Cameroon Natural Hazards. 2010;55:599–619.
Thwing JI, Mihigo J, Fernandes AP, Saute F, Ferreira C, Fortes F, et al. How much malaria occurs in urban Luanda, Angola? A health facility-based assessment. Am J Trop Med Hyg. 2009;80:487–91.
Siri JG, Lindblade KA, Rosen DH, Onyango B, Vulule JM, Slutsker L, et al. A census-weighted, spatially-stratified household sampling strategy for urban malaria epidemiology. Malar J. 2008;7:39.
Wang SJ, Lengeler C, Mtasiwa D, Mshana T, Manane L, Maro G, et al. Rapid Urban Malaria Appraisal (RUMA) II: epidemiology of urban malaria in Dar es Salaam (Tanzania). Malar J. 2006;5:28.
The authors would like to thank the parents, guardians and all children involved in this study for their cooperation and participation. The authors also acknowledge the following people for their support of this project: Dr. Carl H. Campbell Jr. of the Center for Tropical and Emerging Global Diseases at the University of Georgia, Storn Kabuluzi of the Malawi Ministry of Health, Prof. Cameron Bowie of the University of Malawi College of Medicine, Dr. Leo Zulu of the Department of Geography, Michigan State University, and Prof. Edward Ionides of the Department of Statistics, University of Michigan.
Funding for initial data gathering was provided by the US Centers for Disease Control and the President Malaria Initiative, and supported through the Cooperative Agreement number 5 U01 CI000189 from the CDC. Additional support came from the Department of African and Afro-American Studies and the Global Public Health Program of the University of Michigan. Analysis supported by Grant 1U19AI089683-01 from NIAID for the International Centers of Excellence for Malaria Research (ICEMR) in Malawi.
Ethics approval and consent to participate
The Institutional Review Boards at the College of Medicine, University of Malawi, University of Michigan and the US Centers for Disease Control and Prevention (CDC) approved the study protocol. After explaining the goals, risks and benefits of the study, verbal and written consent for participation was received from household heads.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Univariate and multivariate models of parasitaemia risk. Multivariate model 1 includes the individual components of the composite urbanicity measure except for the last two variables. Here, component 1 and 2 were included in one model. Multivariate model 2 includes the continuous urbanicity composite measure
About this article
Cite this article
Larson, P.S., Eisenberg, J.N.S., Berrocal, V.J. et al. An urban-to-rural continuum of malaria risk: new analytic approaches characterize patterns in Malawi. Malar J 20, 418 (2021). https://doi.org/10.1186/s12936-021-03950-5