Policy Research Working Paper 10469 Cholera Risk in Lusaka A Geospatial Analysis to Inform Improved Water and Sanitation Provision Peter W. Gething Sophie Ayling Josses Mugabi Odete Duarte Muximpua George Joseph Solomon Sitinadziwe Kagulura Water Global Practice June 2023 Policy Research Working Paper 10469 Abstract Urbanization combined with climate change are exacerbat- to networked and non-networked Water and sanitation ing water scarcity for an increasing number of the world’s infrastructure, groundwater vulnerability, and drainage, to emerging cities. Water and sanitation infrastructure, which generate a high-resolution map of cholera risk across the city. in the first place was largely built to cater only to a small The analysis presents scenarios of standalone or combined subsector of developing city populations during colonial investments across sewerage coverage and maintenance, times, are increasingly coming under excessive strain. In on-site sanitation improvements, piped water network cov- the rapidly growing cities of the developing world, expan- erage and quality, and ensuring the safety of point source sion does not always keep pace with population demand, water. It identifies the investment most strongly correlated leading to waterborne diseases, such as cholera (Vibrio with the largest reduction in cholera risk as the provision cholerae) and typhoid (Salmonella serotype Typhi). Fund- of flush to sewer infrastructure citywide. However, it also ing gaps therefore make targeting for efficient spending on considers the trade-offs in terms of financial cost versus infrastructure upgrades essential for reducing the burden health benefits and takes note of where the next highest of disease. This paper applies geospatial analysis in Lusaka, health benefits could be achieved for a much lower cost. Zambia, in the context of the cholera outbreak of Octo- Finally, the analysis was done in the context of a considered ber 2017 to May 2018, to identify different water and restructuring of an existing World Bank investment, the sanitation infrastructure investment scenarios and their Lusaka Sanitation Program. It identifies what appears to be relative impact on reducing the risk of cholera in the city. the most efficient combined initiative as partial sanitation The analysis presented uses cholera case location data and investment scale-up and investment in piped water in 10 geospatial covariates, including the location of and access priority wards where the cholera risk was highest. This paper is a product of the Water Global Practice. It is part of a larger effort by the World Bank to provide open access to its research and make a contribution to development policy discussions around the world. Policy Research Working Papers are also posted on the Web at http://www.worldbank.org/prwp. The authors may be contacted at gjoseph@worldbank.org. The Policy Research Working Paper Series disseminates the findings of work in progress to encourage the exchange of ideas about development issues. An objective of the series is to get the findings out quickly, even if the presentations are less than fully polished. The papers carry the names of the authors and should be cited accordingly. The findings, interpretations, and conclusions expressed in this paper are entirely those of the authors. They do not necessarily represent the views of the International Bank for Reconstruction and Development/World Bank and its affiliated organizations, or those of the Executive Directors of the World Bank or the governments they represent. Produced by the Research Support Team Cholera Risk in Lusaka: A Geospatial Analysis to Inform Improved Water and Sanitation Provision A1 AUTHORS & AFFILIATIONS Peter W. Gethingi,ii § Sophie Aylingiii Josses Mugabi iii Odete Duarte Muximpua iii George Josephiii Solomon Sitinadziwe Kaguluraiv 1 JEL: L95;C21; C19: I18 Key Words: Water, Sanitation, Utilities, Urban, Spatial Modelling, Targeting, Infrastructure, Cholera, 1 i. John Curtin Distinguished Professor, Curtin School of Population Health, Faculty of Health Sciences, Curtin University, Kent St, Bentley, WA 6102, Australia ii. Telethon Kids Institute, Perth Childrens Hospital, Nedlands, WA6009, Australia. iii. Water Global Practice, World Bank Group, 1818 H Street NW, Washington DC, 20433, USA. § Corresponding author The authors gratefully acknowledge Lusaka Water and Sanitation Company for their collaboration namely Jonathan Kampata (Managing Director, LWSC), Kennedy Mayumbelo (Project Manager – LSP), Mwansa Nachula Mukuka (Sanitation Specialist, LSP) Kalikeka Malate (Senior Engineer, Project Planning, LSP), Lusungu Nyirenda (Wastewater Specialist, LSP). The cholera case data was made available with thanks to the Zambian Ministry of Health and Mazyanga Mazaba liwewe, Prof Victor Mukonka, Prof Nathan Kapata and Dr Muzala Kapina from the Zambian National Institute of Public Health (ZNPHI). The household water and sanitation survey data collection was made possible with then Field Coordinators Sensio Banda and Gertrude Namitala; survey firms Fibonacci Engineering and Palm Associates. From inside the World Bank, the work was made possible with the support of operational Task Team Leaders (TTLs) Josses Mugabi, Odete Muximpua, Ruth Kennedy Walker and Ai-ju Huang from the Water Global Practice of the World Bank as well as Collins Chansa, Senior Economist, Health Global Practice. Comments were gratefully received from Solomon Sitinadziwe Kagulura, Senior Health Specialist, Health Global Practice. This work was financed by the World Bank Global Water Security and Sanitation Partnership (GWSSP). The findings, interpretations and conclusions expressed in this paper do not necessarily reflect the views of the World Bank, the Executive Directors of the World Bank or the governments whom they represent. The World Bank does not guarantee the accuracy of the data included in this work. Official delimitation of areas and borders might not reflect the official position of the World Bank Group. Country borders or names do not necessarily reflect the World Bank Group's official position. These maps are for illustrative purposes and does not imply the expression of any opinion on the part of the World Bank, concerning the legal status of any country or territory or concerning the delimitation of frontiers or boundaries. A1 INTRODUCTION Due to a myriad of push and pull factors, such as promising economic opportunities in urban centers, as well as climate or declining wages in rural areas (Hoffmann et al. 2019; Patel et al. 2020), urban populations in developing countries continue to expand at a rapid rate. Such growth in developing cities in Africa and India is due to account for almost two-thirds of the projected growth in the world´s urban population by 2050 (Bryan, Glaeser, and Tsivanidis 2019). Furthermore, the number of large cities exposed to water scarcity is projected to increase from 193 to 284 by 2050 (He et al. 2021). Alongside such rapid growth come public health challenges stemming from the inability of older infrastructure to keep pace with that growth. This is particularly the case for water and sanitation infrastructure, as is commonly seen in informal settlements that grow up in the peri-urban areas of developing cities (Allen, Dávila, and Hofmann 2006; Patel, Shah, and Beauregard 2020; Patel and Shah 2021). This can have important repercussions on public health. The acute need for handwashing facilities as one of the first preventative measures against the spread of COVID-19 brought this issue to the fore globally once again. However, Water and Sanitation Services (WSS) infrastructure has been a long-standing challenge for much older waterborne diseases such as cholera and typhoid (Salmonella serotype Typhi). Meanwhile, the funding to make infrastructural improvements can rarely meet the full demand, with consequences for population health. The water and sanitation utility in Lusaka has taken US$131 million in loans (World Bank 2015; Lusaka Water and Sanitation Company 2022) to pay for upgrades and extensions to the sanitation systems available there, across three international funders – the African Development Bank (ADB), the World Bank, and the European Investment Bank (EIB). There are also other initiatives from the Millennium Challenge Corporation (MCZ) and the Gates Foundation. The first three investments are intended to reach 1.7 million beneficiaries. Despite the seemingly large investment across a number of funders, the budget needs to be stretched to meet required investment, especially for an aging sewer network. This makes targeting essential for reducing the burden of disease. Geospatial approaches can enable spatial targeting of infrastructure to take place. This paper applied such methods to Lusaka, Zambia, in the aftermath of the cholera outbreak of 2017-18, in order to identify different long-term investment scenarios and their relative impact on reducing the risk of cholera in the city. In Lusaka, the capital city of Zambia, recurrent cholera outbreaks have caused significant morbidity and mortality in recent decades (Olu et al. 2013; Kapata et al. 2018). The cholera outbreak in Zambia in 2017- 18 was declared on October 6, 2017, in Lusaka. The number of cases increased from several hundred in early December 2017 to peak at approximately 2,000 by early January 2018 and cumulatively resulted in at least 98 deaths (CFR = 1.8%) by May 2018 (Sinyange et al. 2018). Of the suspected cases, 91.7% occurred among Lusaka residents. In the short term, emergency response measures (Kapata et al. 2018), such as provision of tanked or bottled water, hyper chlorination of drinking water, or oral vaccine campaigns (Bi et al. 2017),(Ferreras et al. 2018; Azman et al. 2016), can reduce disease burden and shorten the duration of the outbreak. In Lusaka, multiple local and international partners engaged in the humanitarian response to include the Ministry of Health, Ministry of Water, the water utility Lusaka Water and Sewerage Company (LWSC), Lusaka City Council (LCC), the Center for Disease Control in Zambia (CDC) and the World Health Organization (WHO). It included provision of safe water supplies to compounds through bowsers, emergency works to water networks in affected areas as well as provision of water through additional boreholes, kiosks and elevated tanks. CDC conducted extensive water quality monitoring, LCC enforced pit latrine emptying and the burying of shallow wells and there was an intensified effort by LWSC to attend to sewer blockage complaints. Approximately 2 million doses of a cholera vaccine were also administered to residents over 1 year of age from January 10 to February 14, 2018. However, these vaccines only provide a short-term immunity and are not intended as a substitute for addressing underlying poor water and sanitation conditions. In fact, due to heavy flooding in March 2018, and associated widespread water shortages, there was a subsequent resurgence despite the vaccine (Sinyange et al. 2018). Longer-term prevention measures necessitate a robust underlying infrastructure so that contaminant sources including human waste are disposed of safely and local populations are able to access routinely safe and clean drinking water. Without such longer-term investments in infrastructural renovation and maintenance, similar outbreaks often recur during the rainy season. Indeed, Zambia has seen frequent cholera outbreaks, with its first reported in 1977 and its most recent in just January 2023, with a case fatality rate has ranged between 2% and 8.3% during the past three decades. The country has registered cases every year except for the 1984-1988, 1994-1995 and 2012- 2015 periods. Outbreaks usually occur between October and June during the rainy season, in rural fishing camps (particularly Lake Kariba, Lake Tanganyika, and Lukanga swamps) and peri-urban areas of Lusaka and the Copperbelt provinces. Zambia is not alone in the region, and cholera continues to be a global issue. Cholera is a diarrheal disease caused by the bacterium Vibrio Cholera when the host ingests contaminated food or water. According to the Global Task Force on Cholera Control (GTFCC), there are up to 143,000 deaths worldwide and up to 4.0 million cases of cholera each year (GTFCC 2022) due largely to inadequate access to safe water and sanitation. Without proper treatment, death can occur within hours due to severe dehydration. However, with proper treatment, the recovery from cholera is almost as dramatic as the disease onset. Patients get well rapidly and recover fully. According to WHO, the cholera Case Fatality Rate should be below 1%. Nonetheless, in 2020, in Sub-Saharan Africa the annual Case Fatality Rate was 1.6% - the highest in the world (WHO 2020) and in 2009, for instance, 98% of the 221,226 notified cases worldwide were from Africa. Further, in recent years, remarkable epidemics struck various African regions located far from the coast. For instance, in 2008–2009, Zimbabwe experienced the largest cholera outbreak ever recorded in Africa, with more than 100,000 cases and more than 4,300 deaths (World Bank 2016). These examples stress the need to better characterize cholera outbreaks in non-coastal regions of Africa and how to use various measures to reduce the cholera risk in susceptible environments. The analysis presented in this paper uses cholera case location data and geospatial covariates to include the location of networked and non-networked WSS infrastructure as well as groundwater data to generate a high-resolution map of cholera risk across the city. It has sought to explain patterns of risk in relation to known risk factors for the disease from the literature, identify priority areas and compare alternative strategies for improved infrastructure. The analysis presents scenarios of standalone or combined investments across sewerage coverage and maintenance; on-site sanitation improvements; piped water network coverage and quality; and ensuring the safety of point source water and the relative costing estimates for each, and number of population it would affect (i.e. cost efficiency). It identifies the investment most strongly correlated with the largest reduction in cholera risk would be the provision of flush to sewer infrastructure citywide. However, such an investment would also be disproportionately the most expensive. It identifies the next most impactful stand-alone investment would be the expansion of piped water city-wide, followed by addressing water quality in the existing network. Finally, the analysis was carried out in the context of a considered restructuring of an existing World Bank investment called the Lusaka Sanitation Program (LSP), and identifies the most efficient combined initiative was partial sanitation investment scale-up and investment in piped water in 10 priority wards, costing US$134 million, and reducing the risk of cholera city-wide by some 48%. In the following sections, this paper outlines the relevant data for Lusaka that was assembled to complete the analysis, such as case location data on cholera cases, water supply and sanitation infrastructure and other potential predictors of cholera to conduct a geospatial analysis for addressing three important objectives: (i) use cholera case location data and geospatial covariates to generate a high-resolution map of cholera risk across the city; (ii) seek to explain patterns of risk in relation to putative causal factors; (iii) compare current water and sanitation infrastructure and access to the pattern of cholera risk, identify priority areas and compare alternative strategies for improved infrastructure. In the results section the paper presents a high-resolution map of cholera risk followed by the results of scenario analysis exploring the difference in level of association between the covariates and risk reduction. It also presents the results of combined scenarios which mirror a series of alternatives proposed under a World Bank investment which was being restructured at the time of the analysis. It ends with discussion and limitations. This paper contributed to informing the targeting of existing investments, and the methods employed here can continue to be used to inform future investments by the Government of the Republic of Zambia once these projects are closed. The granularity of the analysis which has brought together multiple data sources enhances our understanding of the pattern of underlying risk across the city, the role of different putative risk factors, and the potential impact of different infrastructure provision or other mitigation strategies. High quality, spatially referenced, data have been generated on the location of the cholera cases themselves, the configuration and quality of existing infrastructure, population access to safe water and sanitation, and other relevant environmental variables. The modeling techniques used here contribute to a growing field of geospatial modeling work in disease modeling that aims to identify and estimate relative contributions of factors that contribute to outbreaks (Kulldorff and Nagarwalla 1995; Besag and Newell 1991), including cholera. For example in Bangladesh, Ali et al defined areas of high cholera risk based on environmental risk factors such as proximity to surface water, high population density and low educational status (Ali et al. 2002; Bi et al. 2016); in Ghana Anamzui-Ya measured proximity to refuse dumps and water reservoirs within communities of Kumasi and used a spatial conditional autoregressive (CAR) model to determine the spatial dependency of cholera prevalence on both digitized imagery and RapidEye image. They found an inverse spatial relationship between cholera prevalence and proximity to both refuse dumps and classified reservoirs (Anamzui-Ya 2012). In Zimbabwe (Luque Fernandez et al. 2012) identified a spatial pattern in the distribution of cholera cases in an epidemic with Harare, characterized by a lower cholera risk in the suburbs with highest elevation. A parallel study to this one was also conducted in Harare, Zimbabwe taking into account a similar suite of geospatial covariates (Joseph et al. 2023). There have also been other efforts to map cholera risk in Zambia using spatial methods. For example Tambatamba et al analyzed the distribution of cholera cases, the mode of cholera transmission, and the risk factors affecting cholera infection in a peri-urban area of Lusaka by using a Geographic Information System (GIS) and a matched case-control method (Tambatamba et al. 2008). Other authors looked explicitly at precipitation patterns, and the association between drainage network availability and cholera cases with a regression analysis (Sasaki et al. 2009). Mwamba et al took a broader look using a cholera risk in Zambia between 2008 and 2017, identifying 16 districts at higher risk using a discreate Poisson-based space-time scan statistic to allow for variation over the ten-year study period (Mwaba et al. 2020). However, this study provides a unique contribution in terms of the combination of data sources that it brings together, considering multiple known environmental risk factors, and its granular focus on Lusaka. A2 MATERIALS AND METHODS Studies of previous epidemics in Lusaka and elsewhere have implicated a range of different factors as contributing to cholera risk, including inadequate provision of safe water and sanitation (Mengel et al. 2014; Adedire et al. 2015), climatic conditions (Luque Fernández et al. 2009), inadequate drainage (Sasaki et al. 2009) and food-borne contamination (Dubois et al. 2006; Sinkala et al. 2004). We also see associated environmental factors where cholera outbreaks could be more prevalent such as areas with higher levels of poverty (Cowman et al. 2017) and densely populated peri-urban settlements (Chimusoro et al. 2018). In this work multiple data sources that reflect these covariates as far as possible were brought together to identify where the strongest correlates with risk were in the case of the 2017-18 Lusaka outbreak. Table 1 provides an overview of all the data sources that went into the analysis. They include the cholera case data itself, data on poverty rates by ward, population estimates, household size information, groundwater vulnerability data, distance from water and sewer networks and burial grounds and household reported access to treated water, improved sanitation, and handwashing facilities. Finally, data on the risk of flooding; frequency of complaints regarding sewer network; water supply/quality and prevalence of E. coli in drinking water sources were all input into the model. Table 1- Data sources used to generate the covariates in the analysis Data Description Data Type, Resolution, Dates, Obs, Source sampling technique Cholera case data De-identified GPS coordinates for Zambia National Public October 2017-May 2018 (5,444 Health Institute (ZNPHI) cases). Sampling must have followed http://znphi.co.zm/ the ZNPHI and GTFCC protocols for identifying a cholera outbreak. Pre-existing geospatial grids Poverty rates World Bank Group Small Area de la Fuente, A., Murr, A. Estimations – ward level and Rascon, E. (2015) Mapping Subnational Poverty in Zambia. Lusaka: World Bank Group and Zambia Central Statistical Office. Population estimates Household census, global dataset, 2015 Zambia census, Open global raster Street Map , Nightime light brightness Household size Household surveys Night-time light brightness and Global dataset Resolution NASA (2015) ‘ASTER Level 1 Digital Elevation Model Precision Terrain Corrected Registered At- Sensor Radiance Version 3’. Sioux Falls, South Dakota: NASA EOSDIS Land Processes DAAC, USGS Earth Resources Observation and Science (EROS) Center. Groundwater vulnerability to Raster KfW contaminants Covariate grids were generated based on GIS data describing relevant WASH infrastructure or natural features Distance from water network Shapefile LWSC Distance from sewer network Shapefile LWSC Distance from stream/river Raster file (based on hydrological NASA (2015) – as per Digital analysis of Digital Elevation Model Elevation Model data source data) Distance from cemetery/burial Global Raster file Open Street Map ground WASH covariate grids were generated by implementing a Bayesian geostatistical model applied to household or other survey data describing water and sanitation facilities, water quality, and flood risk % HHs with soap in toilet GPS point data from household World Bank Lusaka % HHs not treating water surveys Sanitation Assessment % households with survey (2015) and the unimproved sanitation Lusaka Sanitation Program Density of population with (LSP) baseline survey unimproved sanitation undertaken by Vision RI (2016). Composite risk index for water source Composite risk index for toilet facility Prevalence of E. coli in drinking 290 randomly selected water source NWASCO and CDC water sources locations (taps, boreholes, tanks etc) (CDC) Risk of flooding (poor Raster file Millennium Challenge drainage) Corporation (MCC) Frequency of complaints Complaints were aggregated by LWSC regarding sewer network, Township (the administrative level water supply/quality below the District Metering Area (DMA) used to administer water across the city, with 78 units across Lusaka District) and categorized as relating to water quality, water supply, or sewerage. Figure 1 shows a schematic summary of the main components of the analysis and how they interrelate. The analysis proceeded in three main stages: (i) assembly and standardization of a suite of geospatial covariates, gridded at a resolution of 100m × 100m and representing a range of environmental, sociodemographic and infrastructure covariates potentially related to Cholera risk; (ii) use of these covariates along with data on cholera case locations in a Log-Gaussian Cox Process model to generate a predicted map of underlying cholera risk across the city; and (iii) the use of a counterfactual modeling approach to evaluate and compare the possible impact of different mitigation strategies. The data and analytical steps are now discussed in more detail. Study area and Cholera case data Figure 2 shows the study area, consisting of the 34 Wards that make up the metropolitan District of Lusaka. Data on cholera cases were obtained from the Zambian National Institute of Public Health (ZNPHI) with the study dataset consisting of the 5,444 cases reported between October 2017 and May 2018, with each geopositioned by household location. Constructing covariates and assessing their linear association with cholera risk A suite of geospatial covariates was constructed from the data sources listed in Table 1 to include a range of factors potentially related to cholera risk and explaining some fraction of the observed spatial variation in case incidence. All covariates were defined on a spatial raster grid at 100m × 100m resolution across the study area, and fell into one of three categories, as follows. (i) Hydrology. A digital elevation model(NASA 2015) (DEM) was used to derive the path of main river and stream channels and a raster grid created denoting distance of each grid cell from those channels. Data on reported flooding at 128 point locations across the city were obtained from the Millennium Challenge Account – Zambia (MCA-Z), and these were converted into a raster grid denoting distance to flood prone areas. Raster data were also obtained on the vulnerability to pollution of the underlying groundwater aquifers across Lusaka, based on geological characteristics. This includes taking into account infiltration characteristics including consideration of rock type, composition, tectonic lineaments, catchments or drainage basins, groundwater flow, and surface vegetation (Nick, Mweene, and Baumle 2012). (ii) Water and sanitation infrastructure. Digitized data were obtained from the Lusaka Water and Sewerage Company (LWSC) on the piped water and sewerage infrastructure that underly the city, and these were converted into raster grids denoting distance from the water and sewerage networks. Two further types of data were assembled that reflected aspects of the quality of water and sewerage provision across these networks. First, data were obtained on recorded complaints made by Lusaka residents to LWSC about aspects of the sewer and piped water network. Complaints were aggregated by Township (the administrative level below the District Metering Area (DMA) used to administer water across the city, with 78 units across Lusaka District) and categorized as relating to water quality, water supply, or sewerage. Second, data were obtained from the Zambian office of the Centres for Disease Control ((CDC), 2018) who undertook water sampling at 290 randomly selected water source locations (taps, boreholes, tanks etc.) across the city and recorded whether each had evidence of contamination with Escherichia coli bacteria, which served as a proxy for unsafe drinking water. (iii) Household water and sanitation access and other characteristics. Data were compiled from two household surveys conducted in Lusaka: the World Bank Lusaka Sanitation Assessment survey (2015) and the Lusaka Sanitation Program (LSP) baseline survey undertaken by Vision RI (2016). Both were representative sample, questionnaire-based surveys in which selected households were asked broad ranging questions, including on their access to and use of improved water and sanitation facilities, and latitude and longitude coordinates of households were recorded via GPS. In order to use these detailed survey data as putative covariates of cholera risk across the city, two further steps were undertaken. First, questionnaire responses from each household were condensed into an index of risk for both water and sanitation, as detailed in Table 2. Second, the set of resulting risk level data at household locations was used in a Bayesian geostatistical model (P W Gething, Dasgupta, and Andres 2017; Peter W. Gething and Joseph 2017) to yield an interpolated raster grid (one each for sanitation and water risk). Three other geospatial covariates were also derived from the household survey data using the same approach: the percentage of households with soap in their toilet area (as a proxy for handwashing); the mean household size; and the percentage of households with unimproved sanitation. The latter variable was combined with an underlying population map to derive a population density of people without adequate sanitation across the city. Finally, geospatial data was obtained on the percentage of households classified as being in extreme poverty by municipal ward (de la Fuente, Murr, and Rascon 2015) and the location of graveyards across the city which is a potential risk factor for spread of cholera. Both were converted into raster grids. To assess the association between each covariate and cholera risk, the Pearson's product moment correlation coefficient was calculated. Table 2 Risk indices for water and sanitation derived from household survey responses. Risk levels are presented in ascending order of risk. Risk level Description Risk associated with household source of water 0 Piped to premises 1 Public standpipe/tap or borehole within 100m 2 Public standpipe /tap more than >100m from house 3 Well more than 100mm from the house 4 well within 1-10m of a pit latrine 5 well within 11-40m of a pit latrine 6 Surface water source Risk associated with household access to sanitation 0 Flush to sewer no blockages reported 1 Flush to sewer (reported blockages) 2 onsite improved + unshared + not broken/leaking 3 onsite improved + shared + not broken/leaking 4 onsite improved + unshared + broken/leaking 5 onsite improved + shared + broken/leaking 6 unimproved/ambiguous + shared + not broken 7 unimproved/ambiguous + unshared + not broken 8 unimproved/ambiguous + unshared + broken 9 open defecation Development of a cholera risk map To generate a high resolution map showing spatial variation in underlying cholera risk across Lusaka, a log-Gaussian Cox Process (LGCP) model was developed (Diggle et al. 2013). LGCP models are appropriate when the data in question can be considered the result of a spatial point process, denoting the location of events (in this instance, cholera cases), and the objective is to infer an underlying spatial surface (a map) denoting how the intensity of occurrence varies across the study area. Formally, the model was structured as follows: ( , ) = ( , )exp{Υ(, )} where is the intensity (in our application, cholera incidence rate) at a spatial location ( , ); is the fixed spatial component; and Υ is a spatial Gaussian Process (a multivariate normal Gaussian distribution characterized by a mean and spatial covariance function). The fixed spatial component is then defined as: = exp{} where is a vector of coefficients and is a vector of geospatial covariates selected from those described above. Covariate selection was performed via an exhaustive screening process (using the glmulti package in R, (Calcagno 2020)) that formulated every possible model based on combinations of the candidate covariates and compute the Akaike Information Criteria for each. Because this related only to the fixed spatial component, candidate models were constructed as standard generalized linear models using the log link function (i.e. Poisson regression), rather than each requiring full inference of the spatial Gaussian process. This procedure also yielded a measure of the relative ‘importance’ of each covariate as a predictor by averaging the performance of models (expressed as an Akaike weight) in which each covariate appeared across the set of all possible models. Having identified the best performing set of covariates, inference of the full LGCP model was performed via Markov Chain Monte Carlo (MCMC) simulation, and all analysis was conducted in R using the lgcp package (Taylor et al. 2015). The model output consisted of posterior realizations of the intensity surface defined on the same 100m × 100m resolution spatial grid as the input covariates. From these realizations, summary maps were generated including the posterior mean, which was used as the final point estimate map. Model fit was directly assessed by comparing observed incidence rates (calculated as the number of cases reported in each pixel, annualized, and divided by the population of the pixel) against predicted values from the point estimate surface. Correspondence between the observed and predicted value across the data set was summarized by calculating the correlation coefficient, mean error (as a measure of overall bias), and mean absolute error (as a measure of overall precision). Counterfactual analysis to explore impact and targeting of improved water and sanitation infrastructure One feature of the LGCP model described above is that it enabled exploration of the multivariate relationship between cholera risk and specific risk factors. Some of those risk factors relate to inadequate access to safe water and improved sanitation, or exposure to certain environmental risks, that can potentially be addressed by deliberate intervention strategies. To explore and compare the potential impact of different mitigation actions, a counterfactual analysis was conducted using the following steps. First, a set of six intervention scenarios were defined: (A) provision of piped water to house premises; (B) ensuring that households have, as a minimum, access to a public tap within 100m; (C) provision of flush to sewer by connecting households to sewer network; (D) ensuring that households have, as a minimum, access to a shared improved onsite sanitation facility; (E) reducing risk of flooding (eliminating risk, or reducing to ‘medium’ or ‘low’ risk levels); (F) eliminating E.coli contamination in LWSC water sources. Scenarios A, B, and C relate to enhancements to the existing networks of piped water and sewerage and would be progressively costlier to implement as their extent expanded further from the current networks. As such, three variants were defined for each of those scenarios whereby the enhancements would apply to areas within 500m of the existing networks, within 1000m, and finally extending to the entire city. Second, for each scenario the relevant geospatial covariate grid was modified to represent the counterfactual scenario. For example, provision of piped water to all premises in the city was represented by replacing the real-world covariate grid for water risk with one in which in which the risk level was set to the minimum level of 1 (see Table 2) throughout the city; elimination of E. coli risk was represented by replacing the real-world covariate grid with one in which risk was set to zero throughout, and so on. Third, for each scenario the LGCP model predicted using the real-world covariate was compared to a scenario model that used instead the counterfactual covariate to predict risk across the city. Fourth, the resulting level of risk (number of predicted cases across the city) predicted under each counterfactual scenario was compared to the real-world number, thus providing an estimate of the impact of the proposed intervention scenario. Finally, the impact for each scenario was computed within each individual Ward, thereby allowing Wards to be ranked for the purposes of potential targeting of interventions. A3 RESULTS Exploration of factors driving spatial risk patterns Table 3 shows the correlation coefficients between each putative risk factor and observed cholera cases per grid cell. These correlations are simply an indicator of the bivariate association between each factor and cholera risk. Covariates most strongly associated with increased risk included the density of population with unimproved sanitation, the sanitation risk index (see Table 2), and the prevalence of E. coli contamination in water sources. Decreasing risk was most strongly associated with increasing distance from flood-prone areas. Table 3 also shows results from the multivariate model selection procedure that compared models with every possible combination of these risk factors as covariates. Covariates consistently included in better performing models have a higher ‘relative importance’ score (standardized to 1 for the most important). The final subset of covariates identified as forming the best predicting model are also indicated. Table 3 Associations between putative risk factors and observed cholera risk, and their inclusion in the LGCP model. For each candidate geospatial covariate , the Pearson's product moment correlation coefficient is given, along with its associated P-value and upper and lower bounds of its 95% confidence interval. Also shows are the relative importance values indicating the predictive utility of each covariate and whether it was included in the identified best model. Covariates are listed in descending order of relative importance. DD = decimal degrees (1 DD is approximately 111km). Geospatial covariate Corr. p-value L95 U95 Importance Include d - Distance to flood-prone area (DD) -0.115 <0.001 -0.106 1 Yes 0.123 - Groundwater vulnerability -0.098 <0.001 -0.089 1 Yes 0.107 - Number of complaints - sewerage -0.074 <0.001 -0.062 1 Yes 0.086 Number of complaints - water supply 0.059 <0.001 0.047 0.071 1 Yes Prevalence of E. coli contamination (LWSC 0.147 <0.001 0.138 0.156 1 Yes sources) Population density with unimproved 0.249 <0.001 0.241 0.258 1 Yes sanitation (per 100m2) Water risk index 0.087 <0.001 0.079 0.096 0.999 Yes Sanitation risk index 0.139 <0.001 0.130 0.147 0.999 Yes - Distance to sewer (DD) -0.06 <0.001 -0.051 0.990 Yes 0.069 - Distance to stream (DD) -0.042 <0.001 -0.034 0.990 Yes 0.051 % HHs with unimproved sanitation 0.142 <0.001 0.134 0.151 0.853 Yes - Distance to graveyard (DD) -0.044 <0.001 -0.035 0.683 Yes 0.052 Combined water-sanitation risk index 0.15 <0.001 0.141 0.159 0.29 No Prevalence of E. coli contamination (non- - -0.068 <0.001 -0.059 0.245 No LWSC sources) 0.076 - Distance to piped network (DD) -0.088 <0.001 -0.08 0.231 No 0.097 - Mean size of HH -0.001 0.773 0.007 0.102 No 0.010 % HHs with soap in toilet 0.090 <0.001 0.081 0.099 0.043 No Number of complaints - water quality 0.012 0.052 0.000 0.023 0.032 No % HHs in extreme poverty 0.076 <0.001 0.067 0.085 0.002 No A high resolution map of cholera risk in Lusaka Figure 3 shows the 100m X 100m predicted map of cholera risk across the study area, along with the prediction uncertainty per pixel. The geographical pattern of risk is heterogeneous and predicted annual incidence rates vary from zero to in excess of five cases per 1,000 people per annum. The greatest concentration of elevated risk lies in the western Wards of Kanyama and Harry Nkumbula, which both house densely populated informal settlements, and Mwembeshi to the North West. Other notable foci of risk are predicted in Nkoloma, Chawama and Kamwala Wards to the south and, to the east, in Chainda and Mtendere. Areas of higher risk to the western and eastern extents of the city are bisected by a central strip of much lower risk running North-West to South-East. This central area broadly corresponds with the more affluent regions of the city which are also more well covered by the existing sewer network. Table 4 shows results of the model validation where predicted incidence was compared to observed values within each pixel. While linear association between observed and predictive values was relatively weak (correlation of 0.48), the mean absolute error indicates reasonably high precision (the average magnitude of absolute difference in predicted versus actual case incidence per 1,000 people per annum was 1.88) and the mean error indicates minimal overall bias (when over- and under-estimation are considered together, the overall tendency was to under-estimate case incidence by just 0.4 cases per 1,000 per annum). Table 3 Model performance statistics. Values summarize the comparison of predicted versus observed values. Mean absolute error and mean error are in the same units as the modelled metric i.e. cholera incidence per 1,000 people per annum. Validation statistic Value Correlation 0.48 Mean absolute error 1.88 Mean error -0.4 Exploration of potential impact and targeting of improved water and sanitation infrastructure Table 5 shows the results of the counterfactual analysis estimating the potential impact on cholera risk of six alternative mitigation scenarios, expressed as the estimated percentage reduction in cases city- wide if each scenario was enacted. The most impactful scenario was the provision of flush to sewer to all households, which yielded an estimated 90% reduction in cholera cases if implemented across the entire city. Implementing this intervention within only 1,000m and 500m of the existing sewer network was estimated to yield a 29% and 13% reduction in cases, respectively. The second most impactful scenario was the provision of piped water to all households which, if implemented universally, led to an estimated 61% drop in cases (31% and 25% if restricted to 1,000m and 500m of the existing water network, respectively). Ensuring all households have access to at least an improved shared sanitation facility yielded an estimated 56% reduction, and elimination of E.coli risk a 52% reduction. Ensuring all households had access to a public tap within 100m of their household yielded much smaller reductions in cases of just 6%, even if implemented universally. This could be because coverage of piped water is already quite high, but in addition is an interesting finding as it may reflect additional contamination that can occur when carrying water from the standpipe to homes, or reflect unsafe storage. Table 4 Estimated impact and efficiency of alternative risk mitigation scenarios. % reduction in Population Efficiency Action Location cases city-wide targeted ratio Provide piped water to Everywhere -60.89 1,761,855 0.35 premises <1000m from network -31.04 770,622 0.40 <500m from network -25.14 619,937 0.41 Ensure at least public tap Everywhere -5.95 532,065 0.11 within 100m <1000m from network -3.24 244,082 0.13 <500m from network -2.82 210,473 0.13 Provide universal flush to Everywhere -89.57 1,771,116 0.51 sewer <1000m from network -28.78 689,562 0.42 <500m from network -12.71 412,765 0.31 Ensure at least improved Everywhere -56.35 1,457,934 0.39 shared Eliminate E.coli risk Everywhere -52.00 1,771,116 0.29 Existing piped water network -21.03 738,937 0.28 Reduce flood risk Eliminate risk everywhere -29.55 1,771,116 0.17 Reduce to ‘low’ risk everywhere -10.57 748,480 0.14 Reduce to ‘medium’ risk everywhere -2.18 213,726 0.10 Alongside overall impact, a second consideration when evaluating different mitigation strategies is the relative cost and level of effort required to implement them. As a simple indicator of relative cost, the total population that would require provision under each mitigation scenario was identified. An ‘efficiency ratio’ was then computed for each scenario as simply the relative ratio of the reduction in cases versus the size of population requiring the intervention. Thus, scenarios with large efficiency ratios yielded relatively larger impact in case reductions given the number of people requiring the intervention. Under this metric, the provision to households of flush to sewer facilities city-wide (efficiency ratio 0.51) or within 1000m of the sewer network (0.42) were the two most efficient scenarios; with provision to households of piped water within 500m (0.41) and 1000m (0.40) of the water network next most efficient. Figure 4 shows how each of the six intervention scenarios could be geographically prioritized across the city, as measured in terms of the magnitude of overall case reduction that was estimated to be achieved within each municipal ward. As can be seen, the Kanyama ward to the far west of the city is identified as the Ward in which most impact could be achieved across all six scenarios, and the neighboring Harry Nkumbula Ward to the immediate east is second in terms of impact for five of the six scenarios. Other Wards commonly featuring in the top five ranks include Raphael Chota, Lima, and Mtendere. Figures 7, 8 and 9 also show how different wards would be targeted by different interventions based on their potential to reduce cholera risk, depending on whether they are piped water, onsite or networked sanitation options. Combined targeting scenarios A final round of counterfactual analyses worked in conjunction with the aforementioned investment, Lusaka Sanitation Program. The US$65 million investment in upgrading Lusaka’s sewer infrastructure also had a component to provide discounted onsite sanitation options to households in peri-urban areas. For this part of the exercise, three categories of scenarios were considered that included either maintaining the project as it was planned in 2015, much prior to the cholera outbreak (a “status quo” option), a partial scale-up of the existing investment operation, or a full scale-up. Each of the three had one or two sub- variants. For a summarized list of the scenarios, see Table 6. Table 5 Table showing investment scenarios considered for World Bank intervention expansions. Each of the scenarios pertain to scenarios considered for restructuring of the World Bank financing that was being discussed in April 2019. Scenario 1 pertained to a reduced scope of the investments to cover Onsite Sanitation (OSS) in WB targeted areas and the sewer network planned for the Year 1 catchment of the project. Scenario 2 referred to Scenario 1 plus sewer network planned for Years 2-5 of the project (Kanyama and a secondary catchment area called CSE- 20). Scenario 3A covered a slightly higher proportion of the Years 2-5 project areas. Scenario 3B was 3A plus improving water quality by addressing water and sanitation related complaints on the existing network. Scenario 3C was 3A plus improving water access via extension of the existing network with no upgrades to cover the whole population in the top 10 ranked high wards for piped on premises access. Scenario 4A involved covering the whole of the years 2-5 areas for the sewer network. Scenario 4B was 4A plus improving water via addressing e-coli in existing network + eliminating water and sanitation related complaints. Scenario 4C was 4A plus improving water access via extension of the water network in its current state in the top 10 ranked high impact target wards for pipon premises access. US$90 million in additional financing. Population targeted Ensure at % Provide least reduction Provide Eliminate Eliminate Intervention piped improved+ Eliminate in flush to sewer water Type water to shared E.coli risk cases city- sewer complaints complaints premises onsite wide facility Scenario 1 - Sanitation – WB OSS+Y1 -23.7 59,094 303,500 Scenario 2 - Status quo Sanitation – WB OSS+ YI+ Very Limited Y2 -28.3 89,233 298,152 Partial Scenario 3A - scale Sanitation – -28.3 89,432 298,152 WB OSS+ YI + Limited Y2 Scenario 3B - . Sanitation – Limited Y2 + Water Quality -35.2 89,432 298,152 738,937 1,021,299 1,021,299 Scenario 3C - Sanitation – limited Y2 + Increased Water Access -47.7 947,477 89,432 298,152 Scenario 4A - Sanitation – complete Y2 -31.2 172,723 298,152 Scenario 4B - Sanitation – complete Y2 + Water Quality -37.6 172,723 298,152 738,937 1,021,299 1,021,299 Scenario 4C - Full scale up Sanitation – complete Y2 + Increased Water Access -50.6 947,477 172,723 298,152 In terms of the sewer network, the status quo options included year 1 sewers planned in an area of the city, which was already broadly covered by an older network, was largely middle income or was made up of commercial buildings (Emmasdale, Chaisa and Kafue Road), plus very limited works in one peri-urban area (Kanyama) which had been affected by the outbreak. The partial scale-up included extending the sewer network also to an industrial area, with very little population impact but with revenue generation potential for LusakaF Water Supply and Sanitation Company (LWSC). The full scale-up from a sewer perspective involved extending the sewer network to cover several other wards in the city to include Chawama/Kuomboka and Garden, which did have the potential to impact on cholera risk areas. These were the full works originally planned as part of the intervention but were cost constrained for completion. Within the partial and full scale-up options, two sub-scenarios were considered to include either re- habilitating the existing water infrastructure to ensure improved water quality or expanding the network to cover the ten wards identified as highest priority for targeting piped water on premises interventions that would have the greatest impact on eliminating cholera. For all three categories of scenario – status quo, partial and full scale-up-, onsite sanitation was included to enable access to 10,000 toilets at a highly subsidized price in peri-urban communities. For each of these scenarios, the combined level of risk reduction was calculated using a multivariate model that enabled consideration of the complementary effects of combined interventions in provision of piped water/elimination of e-coli in the water supply and providing onsite sanitation/sewer network expansion as appropriate. The findings (Table 6) demonstrated that the existing World Bank investment intervention as planned at the time, valuing US$65 million, was associated with a reduction to the risk of cholera in the city of between 24% and 28%. Partial and full scale-up of the sewer network was associated with a reduction of the risk of cholera by 28% and 31% respectively. However, greater impact was seen in adding additional components of addressing water quality through providing maintenance to eliminate the risk of e-coli in the existing network which was associated with a 35%-38% reduction in cholera risk when combined with the partial and full scale-up of the water network. Finally, the strongest association with reduced cholera risk was observed on increasing access to the water network by targeting the top ten wards, prioritized as part of the mapping exercise for investment in piped water supply. The combination of the partial sanitation investment scale-up and investment in piped water in 10 priority wards would cost US$134 million, but was associated with a reduced risk of cholera city-wide by some 48%. Meanwhile investing in the full scale-up of sewer network investments could cost up to US$155 million, but was associated with a 51% reduction in cholera risk city-wide. On deciding which investment scenario would be most cost efficient, the study also assessed the cost per case reduced for every 10,000 people, so that relative cost efficiency and health efficiency of investments could be compared. Full results can be found in the annex, but each scenario was costed as between US$1.89 and US$3.76 per case reduced for every 10,000 people. A4 DISCUSSION The main contribution of this paper is the demonstration of how geospatial methods can be employed to target water supply and sanitation investments to attain desired outcomes – reduction of cholera risk in the city of Lusaka in Zambia. This study has shown how a wide variety of very detailed geospatial data can be brought together in a formal spatial statistical analysis to explore the geographic patterns and potential drivers of cholera risk during an urban epidemic. The high-resolution risk map provides a highly granular information source for decision makers considering how to reduce risk of future outbreaks, clearly delineating neighborhoods with elevated underlying risk from those where, even in the midst of the epidemic, risk was low. The analysis has also demonstrated how cholera risk arises from a complex set of interacting factors. Unsurprisingly, water and sanitation access play a key role, but this is mediated by environmental factors such as poor drainage and variations in the vulnerability to contamination of the underlying groundwater. The approach presented here has also enabled estimation of the possible impact of different cholera risk mitigation strategies to improve water and sanitation infrastructure, quality, and access, including their relative efficiency and potential to target the highest impact parts of the city. The fact that two Wards (Kanyama and Harry Nkumbula) were ranked as first and second highest impact for nearly all mitigation scenarios reflects the clear need to prioritize this region of the city. The potential for greatest impact here likely stems from the confluence of a number of factors: the very high concentration of cholera risk (meaning predicted incidence rates are high), the dense population (meaning larger numbers can potentially be averted), and the poor levels of access to water and sanitation at the time (meaning there was large scope to improve on current provision). All three factors, in turn, were driven by the nature of these neighborhoods which were characterized by high rates of poverty and informal settlements with inadequate infrastructure. The approach developed here has a number of limitations. All analysis stemmed from a model that sought to explain the spatial variation in reported cases across the city. These case data were geolocated according to the home residence of the patient, but of more interest is the location at which the infection was acquired. In many cases, this may be in or close to the home, but in others the exposure may have occurred somewhere else entirely – for example at a workplace, when visiting another household in a different parts of the city, or in a public venue. In future analyses, it may be possible to obtain data representing patterns of human movement (for example as collected via mobile phone records) and include these in analysis of spatial risk patterns (Bengtsson et al. 2015). A second limitation concerns the counterfactual analysis. This relied on the necessary assumption that the empirical relationships captured in the modeling between cholera risk and each putative risk factor were causal rather than mere associative. For example, it was assumed that the observed relationship between poor sanitation and cholera meant that improving sanitation would reduce cholera risk. This seems largely defensible since there are established causal pathways linking the two, and since other factors were controlled for in the analysis, but results nevertheless remain contingent on the assumption of causation. For some risk factors, the causal pathway may be more indirect. The elimination of E. coli contamination in piped water, for example, was estimated to yield reductions in cholera risk. While E. coli is, of course, not the pathogen responsible for cholera, it is plausible that actions taken to eliminate it by improving water infrastructure and treatment would also be effective against Vibrio cholerae itself. A third limitation pertains to the distinction between access to and use of water and sanitation infrastructure. This study deals solely with the question of access and does not claim to know the extent to which city residents are using the available WSS infrastructure. There is a lack of available data on the extent of household use of WSS infrastructure, but if it were later to become available, it could be incorporated in the model. To compare the possible effort required to implement different scenarios, we used the beneficiary population in each geographic area that would have to receive the intervention. This is an inevitably crude proxy that does not account for the different unit or per-capita costs of, for example, providing a household with flush to sewer versus providing piped water. Further analyses could attempt to draw up actual costings to enable more refined comparisons. A5 CONCLUSION Ongoing epidemics in cities across the developing world demonstrate the persistent threat of Cholera and its potentially devastating impact. International efforts led by the Global Task Force on Cholera Control (GTFCC) have been focused on rapid response through stepping up vaccine availability, while a long-term solution rests in the provision of clean water and safe sanitation infrastructure facilities. Basic improvements in infrastructure, alongside more specific time-limited interventions, can dramatically reduce risk but, in heavily resource constrained settings, must be effectively chosen and targeted in order to maximize their impact. This study has showcased how the growing availability of geospatial data on water and sanitation infrastructure, environmental characteristics and other risk factors can be coupled with spatial statistical analysis to provide granular and robust information to support more precise and impactful infrastructure interventions to improve public health. A6 REFERENCES Adedire, E. B., A. B. Usman, G. A. Abbass, I. O. Ajayi, and O. I. Fawole. 2015. “Descriptive Characterization of Cholera Epidemic Caused by Break Down of Public Pipe Borne Water Supply-Egbeda, Oyo State Nigeria September 2013.” International Journal of Epidemiology 44 (suppl_1): i153–i153. https://doi.org/10.1093/ije/dyv096.189. Ali, Mohammad, Michael Emch, J. P. Donnay, Mohammad Yunus, and R. B. Sack. 2002. “Identifying Environmental Risk Factors for Endemic Cholera: A Raster GIS Approach.” Health and Place. https://doi.org/10.1016/S1353-8292(01)00043-0. Allen, Adriana, Julio D Dávila, and Pascale Hofmann. 2006. “The Peri-Urban Water Poor: Citizens or Consumers?” Environment and Urbanization 18 (2): 333–51. https://doi.org/10.1177/0956247806069608. Anamzui-Ya, J.A. 2012. “Spatial Analysis and Mapping of Cholera Causing Factors in Kumasi, Ghana.” University of Twente. Azman, Andrew S., Lucy A. Parker, John Rumunu, Fisseha Tadesse, Francesco Grandesso, Lul L. Deng, Richard Laku Lino, et al. 2016. “Effectiveness of One Dose of Oral Cholera Vaccine in Response to an Outbreak: A Case-Cohort Study.” The Lancet Global Health. https://doi.org/10.1016/S2214-109X(16)30211-X. Bengtsson, Linus, Jean Gaudart, Xin Lu, Sandra Moore, Erik Wetter, Kankoe Sallah, Stanislas Rebaudet, and Renaud Piarroux. 2015. “Using Mobile Phone Data to Predict the Spatial Spread of Cholera.” Scientific Reports 5 (1): 8923. https://doi.org/10.1038/srep08923. Besag, Julian, and James Newell. 1991. “The Detection of Clusters in Rare Diseases.” Journal of the Royal Statistical Society. Series A (Statistics in Society) 154 (1): 143. https://doi.org/10.2307/2982708. Bi, Qifang, Andrew S. Azman, Syed Moinuddin Satter, Azharul Islam Khan, Dilruba Ahmed, Altaf Ahmed Riaj, Emily S. Gurley, and Justin Lessler. 2016. “Micro-Scale Spatial Clustering of Cholera Risk Factors in Urban Bangladesh.” PLoS Neglected Tropical Diseases. https://doi.org/10.1371/journal.pntd.0004400. Bi, Qifang, Eva Ferreras, Lorenzo Pezzoli, Dominique Legros, Louise C. Ivers, Kashmira Date, Firdausi Qadri, et al. 2017. “Protection against Cholera from Killed Whole-Cell Oral Cholera Vaccines: A Systematic Review and Meta-Analysis.” The Lancet Infectious Diseases. https://doi.org/10.1016/S1473-3099(17)30359-6. Bryan, Gharad, Edward Glaeser, and Nick Tsivanidis. 2019. “Cities in the Developing World.” w26390. Cambridge, MA: National Bureau of Economic Research. https://doi.org/10.3386/w26390. Calcagno, V. 2020. “Glmulti: Model Selection and Multimodel Inference Made Easy. R Package Version 1.0.8.” 2020. https://cran.r-project.org/web/packages/glmulti/index.html. Chimusoro, Anderson, Stephen Maphosa, Portia Manangazira, Isaac Phiri, Tonderai Nhende, Sydney Danda, Ottias Tapfumanei, Stanley Munyaradzi Midzi, and Juliet Nabyonga-Orem. 2018. “Responding to Cholera Outbreaks in Zimbabwe: Building Resilience over Time.” In Current Issues in Global Health, edited by David Claborn. IntechOpen. https://doi.org/10.5772/intechopen.79794. Cowman, Gretchen, Shikanga Otipo, Ian Njeru, Thomas Achia, Harsha Thirumurthy, Jamie Bartram, and Jackson Kioko. 2017. “Factors Associated with Cholera in Kenya, 2008-2013.” Pan African Medical Journal 28. https://doi.org/10.11604/pamj.2017.28.101.12806. Diggle, Peter J., Paula Moraga, Barry Rowlingson, and Benjamin M. Taylor. 2013. “Spatial and Spatio- Temporal Log-Gaussian Cox Processes: Extending the Geostatistical Paradigm.” Statistical Science. https://doi.org/10.1214/13-sts441. Dubois, A. E., M. Sinkala, P. Kalluri, M. Makasa-Chikoya, and R. E. Quick. 2006. “Epidemic Cholera in Urban Zambia: Hand Soap and Dried Fish as Protective Factors.” Epidemiology and Infection. https://doi.org/10.1017/S0950268806006273. Ferreras, Eva, Elizabeth Chizema-Kawesha, Alexandre Blake, Orbrie Chewe, John Mwaba, Gideon Zulu, Marc Poncin, et al. 2018. “Single-Dose Cholera Vaccine in Response to an Outbreak in Zambia.” New England Journal of Medicine. https://doi.org/10.1056/NEJMc1711583. Gething, P W, B Dasgupta, and Luis Andres. 2017. “Geospatial Modeling of Water and Sanitation Indicators in Nigeria in 2015. Report Prepared for the World Bank.” Washington DC: The World Bank Group. Gething, Peter W., and G Joseph. 2017. “Geospatial Analysis of Access to Safe Water and Sanitation in Tanzania and Their Association with Poverty and Health Outcomes. Report Prepared for the World Bank.” Washington DC: The World Bank Group. GTFCC. 2022. “About Cholera -A Disease of Inequity That Strikes the World’s Poorest and Most Vulnerable People.” 2022. https://www.gtfcc.org/about-cholera/. He, Chunyang, Zhifeng Liu, Jianguo Wu, Xinhao Pan, Zihang Fang, Jingwei Li, and Brett A. Bryan. 2021. “Future Global Urban Water Scarcity and Potential Solutions.” Nature Communications 12 (1): 4667. https://doi.org/10.1038/s41467-021-25026-3. Hoffmann, Ellen M., Verena Konerding, Sunil Nautiyal, and Andreas Buerkert. 2019. “Is the Push-Pull Paradigm Useful to Explain Rural-Urban Migration? A Case Study in Uttarakhand, India.” Edited by W. David Walter. PLOS ONE 14 (4): e0214511. https://doi.org/10.1371/journal.pone.0214511. Joseph, George, Sveta Milusheva, Hugh Sturrock, Tonderai Mapoko, Sophie Charlotte Amy, and Yi Rong Hoo. 2023. “The Importance of Maintenance: Geospatial Analysis of Cholera Risk and Water and Sanitation Infrastructure in Harare, Zimbabwe.” Kapata, Nathan, Nyambe Sinyange, Mazyanga Lucy Mazaba, Kunda Musonda, Raymond Hamoonga, Muzala Kapina, Khozya Zyambo, et al. 2018. “A Multisectoral Emergency Response Approach to a Cholera Outbreak in Zambia: October 2017-February 2018.” Journal of Infectious Diseases. https://doi.org/10.1093/infdis/jiy490. Kithiia, Justus, and Gamoyo Majambo. 2020. “Motion but No Speed: Colonial to Post-Colonial Status of Water and Sanitation Service Provision in Mombasa City.” Cities 107 (December): 102867. https://doi.org/10.1016/j.cities.2020.102867. Kulldorff, Martin, and Neville Nagarwalla. 1995. “Spatial Disease Clusters: Detection and Inference.” Statistics in Medicine 14 (8): 799–810. https://doi.org/10.1002/sim.4780140809. la Fuente, A de, A Murr, and E Rascon. 2015. “Mapping Subnational Poverty in Zambia.” Lusaka: World Bank Group and Zambia Central Statistical Office. Luque Fernandez, Miguel A., Michael Schomaker, Peter R. Mason, Jean F. Fesselet, Yves Baudot, Andrew Boulle, and Peter Maes. 2012. “Elevation and Cholera: An Epidemiological Spatial Analysis of the Cholera Epidemic in Harare, Zimbabwe, 2008-2009.” BMC Public Health. https://doi.org/10.1186/1471-2458-12-442. Luque Fernández, Miguel Ángel, Ariane Bauernfeind, Julio Díaz Jiménez, Cristina Linares Gil, Nathalie El Omeiri, and Dionisio Herrera Guibert. 2009. “Influence of Temperature and Rainfall on the Evolution of Cholera Epidemics in Lusaka, Zambia, 2003-2006: Analysis of a Time Series.” Transactions of the Royal Society of Tropical Medicine and Hygiene. https://doi.org/10.1016/j.trstmh.2008.07.017. Lusaka Water and Sanitation Company. 2022. “Lusaka Sanitation Programme Factsheet.” Mengel, Martin A., Isabelle Delrieu, Leonard Heyerdahl, and Bradford D. Gessner. 2014. “Cholera Outbreaks in Africa.” In Cholera Outbreaks, edited by G. Balakrish Nair and Yoshifumi Takeda, 379:117–44. Current Topics in Microbiology and Immunology. Berlin, Heidelberg: Springer Berlin Heidelberg. https://doi.org/10.1007/82_2014_369. Mwaba, John, Amanda K. Debes, Patrick Shea, Victor Mukonka, Orbrie Chewe, Caroline Chisenga, Michelo Simuyandi, et al. 2020. “Identification of Cholera Hotspots in Zambia: A Spatiotemporal Analysis of Cholera Data from 2008 to 2017.” Edited by Adam Akullian. PLOS Neglected Tropical Diseases 14 (4): e0008227. https://doi.org/10.1371/journal.pntd.0008227. NASA. 2015. “ASTER Level 1 Precision Terrain Corrected Registered At-Sensor Radiance Version 3.” Sioux Falls, South Dakota: NASA EOSDIS Land Processes DAAC, USGS Earth Resources Observation and Science (EROS) Center. Nick, A, R Mweene, and R Baumle. 2012. “Groundwater Resources of the Mwembeshi and Chongwe Catchments, Including the Lusaka Region. A Manual with Explanations for the Use of Hydrogeological Maps and Vulnerability Map.” Lusaka: Ministry of Lands, Energy and Water Development, Department of Water Affairs and Federal Institute for Geosciences and Natural Resources. Olu, O., O. Babaniyi, P. Songolo, B. Matapo, E. Chizema, M. Kapin’a-Kanyanga, E. Musenga, and O. Walker. 2013. “Cholera Epidemiology in Zambia from 2000 to 2010: Implications for Improving Cholera Prevention and Control Strategies in the Country.” East African Medical Journal. Patel, Amit, George Joseph, Namesh Killemsetty, and Sokha Eng. 2020. “Effects of Residential Mobility and Migration on Standards of Living in Dar Es Salaam, Tanzania: A Life-Course Approach.” Edited by Wenjia Zhang. PLOS ONE 15 (9): e0239735. https://doi.org/10.1371/journal.pone.0239735. Patel, Amit, and Phoram Shah. 2021. “Rethinking Slums, Cities, and Urban Planning: Lessons from the COVID-19 Pandemic.” Cities & Health 5 (sup1): S145–47. https://doi.org/10.1080/23748834.2020.1790252. Patel, Amit, Phoram Shah, and Brian E. Beauregard. 2020. “Measuring Multiple Housing Deprivations in Urban India Using Slum Severity Index.” Habitat International 101 (July): 102190. https://doi.org/10.1016/j.habitatint.2020.102190. Sasaki, Satoshi, Hiroshi Suzuki, Yasuyuki Fujino, Yoshinari Kimura, and Meetwell Cheelo. 2009. “Impact of Drainage Networks on Cholera Outbreaks in Lusaka, Zambia.” American Journal of Public Health. https://doi.org/10.2105/AJPH.2008.151076. Sinkala, M., M. Makasa, F. Mwanza, P. Mulenga, P. Kalluri, R. Quick, E. Mintz, R. M. Hoekstra, and A. DuBois. 2004. “Cholera Epidemic Associated with Raw Vegetables - Lusaka, Zambia, 2003- 2004.” Journal of the American Medical Association. https://doi.org/10.1001/jama.292.17.2077. Sinyange, Nyambe, Joan M Brunkard, Nathan Kapata, Mazyanga Lucy Mazaba, Kunda G Musonda, Raymond Hamoonga, Muzala Kapina, et al. 2018. “Cholera Epidemic - Lusaka, Zambia, October 2017-May 2018.” MMWR. Morbidity and Mortality Weekly Report. https://doi.org/10.15585/mmwr.mm6719a5. Tambatamba, Bushimbwa, Philip Mulenga, Satoshi Sasaki, Hiroshi Suzuki, and Kumiko Igarashi. 2008. “Spatial Analysis of Risk Factor of Cholera Outbreak for 2003–2004 in a Peri-Urban Area of Lusaka, Zambia.” The American Journal of Tropical Medicine and Hygiene 79 (3): 414–21. https://doi.org/10.4269/ajtmh.2008.79.414. Taylor, Benjamin M., Tilman M. Davies, Barry S. Rowlingson, and Peter J. Diggle. 2015. “Lgcp : An R Package for Inference with Spatial and Spatio-Temporal Log-Gaussian Cox Processes .” Journal of Statistical Software. https://doi.org/10.18637/jss.v052.i04. WHO. 2020. “World Health Organization, ‘Weekly Epidemiology Record’, No. 37.” https://apps.who.int/iris/bitstream/handle/10665/334241/WER9537-eng- fre.pdf?sequence=1&isAllowed=y. World Bank. 2015. “INTERNATIONAL DEVELOPMENT ASSOCIATION PROJECT APPRAISAL DOCUMENT Lusaka Sanitation Project (P149091).Pdf.” ———. 2016. “Project Appraisal Document (PAD) National Water Project - Zimbabwe P154861.” A7 FIGURES Figure 1 Schematic Diagram summarizing main components of analysis Figure 2-Study area and cholera case locations. Inset map show location of Lusaka in Zambia. Main map shows Lusaka study area with cholera case location data used in study, along with administrative Ward boundaries. Underlay taken from OpenStreetMap. Figure 3 Predicted map of cholera risk. (A) The predicted incidence rate during the 2017/18 epidemic, expressed as cases per year (although note that all cases fell in the eight months from October 2017 to May 2018) (B) Prediction uncertatiny expressed as Standard Error of the predicted incidence rate per pixel. (C) Wards within the study area. Non-residential areas (including industrial and commercial land, parks) are masked in Panels A and B. A B C D E F Figure 4 - Geographical targeting by Ward of potential interventions. Intervention scenarios shown are (A) elimination of flood risk; (B) elimination of E. coli risk in piped water; (C) provision of piped water to households; (D) ensuring access to a public tap within 100m of households; (E) provision to households of flush to sewer; (F) ensuring households have access to at least an improved shared on-site sanitation facility. For each scenario, Wards are ranked according to the hypothetical reduction in cholera cases resulting from the intervention if deployed in that Ward. Figure 5 - Graph displaying cost efficiency of interventions by Type and Spatial Scope vs. Number of Cholera cases reduced (Highest cost=least efficient) Figure 6 - Graph demonstrating the cost efficiency vs. number of cholera cases reduced by investment scenario explored (highest cost=least efficient) Figure 7 Map and graph showing where in the city interventions can be prioritized for piped water coverage based on the highest reduction in cholera risk Figure 8 Map and graph showing where in the city interventions can be prioritized for at least improved and shared onsite sanitation based on the highest reduction in cholera risk Figure 9 Map and graph showing where in the city should be targeted for sewer connected facilities based on the highest reduction in cholera risk 1 Manuscript for re-submittion to PLOS Water Call for Papers on ‘Safe and Sustainable Water in Cities’ edited 6th April 2023