Mitigation potential of sanitation infrastructure on groundwater contamination by nitrate in maputo

: In Maputo, the capital of Mozambique, nitrate concentrations above 250 mg L − 1 in groundwater have been reported. This happens due to the widespread use of latrines and septic tanks that allow for constant inﬁltration of its content into the soil and eventually to groundwater sources, a situation that is widespread in the Global South and represents a serious threat for human health and for the environment. This is a reﬂection of limited access to safe and adequate sanitation services, which the local authorities have set to improve in the forthcoming decades with a recently commissioned city-wide sanitation masterplan serving as a basis for the works. In this article, we aimed at understanding whether the infrastructure projected in the masterplan would lead to a reduction of nitrogen reaching groundwater. Currently, according to our calculations, almost 500 onnes of nitrogen reach the city’s groundwater sources each year, with the masterplan potentially resulting in a 14% reduction, a small reduction due to its reliance on maintaining and expanding fecal sludge services, without considering investments to improve domestic systems (e.g., construction of contained systems). An alternative, not presented in the Masterplan and put forward by the authors, could be the construction of simpliﬁed sewers in two of the city’s most densely populated neighborhoods, with a potential 29% reduction in nitrogen reaching groundwater.


Introduction
Nitrate is a compound of nitrogen that occurs naturally in moderate concentrations, in many natural environments [1]. In higher concentrations, nitrate becomes a water contaminant that can affect humans and ecosystems alike: high levels of nitrate in drinking water have been linked to health problems, such as blue-baby syndrome, blood disorders, hypertension, and cancer, with WHO recommending 50 mg L −1 as the absolute maximum nitrate concentration in water for human consumption [2]. Nitrate concentrations, above WHO's limit, in surface water can also lead to eutrophication and aquatic biodiversity loss [3].
For the particular case of groundwater contamination, nitrate sources include agriculture fertilizers [4][5][6];, animal waste from farming [7][8][9], and human fecal sludge, urine, and feces, deposited into on-site sanitation systems [1], henceforth defined as OSS systems. The particular impact of up, it can be assumed that there is no driving force for the household to request an emptying service which means that the fecal sludge is lost and that nitrate will continuously reach groundwater. Of the facilities reported to fill up, most (54%) are unsafely emptied by informal operators, with fecal sludge being assumed to be discharged, untreated, into the environment. According to estimates, only 3% of the systems that fill up are safely emptied by formal operators while 43% are considered to be "safely abandoned after use" [47]. For contextualization, it is important to mention that the existing wastewater treatment plant is composed of two anaerobic ponds and two facultative ponds operating in parallel. However, the system is currently overloaded. On the one hand, both anaerobic ponds are full of sludge as these are the discharge points for formal fecal sludge operators. On the other hand, the facultative ponds are 80-90% covered in water hyacinths. Unpublished work made by the authors has shown that the treatment capacity of the plant, as BOD and COD removal rates, is limited. Thus, it is clear that the largest percentage of wastewater and fecal sludge produced in the city is not adequately treated, which leads to serious environmental and health issues [44,46].
Regarding water supply, while half of Maputo's population gets its drinking water from the large, public water provider (Águas da Região de Maputo, ARM), the rest, mainly in the peri-urban areas which are the most densely populated [49] and where the dwellers mostly rely on OSS [47], are supplied by small-scale, independent water providers (SSIPs). These providers depend entirely [50] on a superficial, semiconfined aquifer that has a maximum width of 15 m, and depths ranging from 20 to 40 m [49,51,52]. This aquifer has been reported to have high nitrate concentrations, with many boreholes throughout the city displaying concentrations above 50 mg L −1 , sometimes reaching 250 mg L −1 [50], especially in the neighborhoods where the population mostly relies on OSS [53]. An up-to-date and extensive description of the geo-hydrology of the area, including nitrate concentrations, can be found elsewhere [51].
This situation clearly demonstrates the link between inadequate sanitations infrastructure and services, environmental contamination, and urban health and wellbeing. Therefore, having in mind the needed improvements to the existing sanitation infrastructure in Maputo, the National Water and Sanitation Infrastructure Board (AIAS, in Portuguese), the primary agency responsible for ensuring both water supply in small towns and sanitation in all urban areas in Mozambique [54], commissioned a city-wide sanitation masterplan. This document proposes improvements to the existing sewer network, its expansion, and the construction of new wastewater treatment plants, together with a network of fecal sludge transfer stations and fecal sludge treatment plants. Based on the data collected for Maputo city, the potential effect of the implementation of the City-Wide Sanitation Masterplan (henceforth defined as Masterplan) was evaluated with MFA. The implementation of the Masterplan was evaluated as the baseline scenario for the Maputo case. Additionally, two other options were studied: whereas one assumed the further expansion of sewer infrastructure, the other assumed improvements of domestic sanitation systems and of fecal sludge services in the city. These two optional solutions were not budgeted or projected by the authorities, and were devised by the authors for the purpose of exploring the possible alternatives to improve the effectiveness of the Masterplan. To the knowledge of the authors, it is the first time that MFA is used to assess the environmental impact of planned sanitation infrastructure in the Global South.

Material Flow Analysis Model Development
The system boundary was defined to cover Maputo, for which the flows of water and nitrogen were modelled for the year 2015. In MFA, flows can be represented as substances, such as nitrogen, which, in this model, is classified as a pollution indicator, or goods, which can have a positive (drinking water) or negative (wastewater) economic value for society. The conceptual MFA model is depicted in Figure 1. The model includes a network of 11 processes, representing the most relevant activities for the transformation, storage, and/or transport of water and nitrogen within the system boundary and Sustainability 2018, 10, 858 4 of 18 the 33 flows that connect them. Within the modelled processes, three are natural sinks (Atmosphere (9), Surface water (10), and Groundwater (11)), representing where the water and associated nitrogen will finally settle; nitrogen ending up in (10) and (11) indicates pollution. These three processes were considered sinks and were not balanced, meaning that the total inputs and outputs of these processes were not determined. The other eight processes are the main social-economic processes in the study area that mobilize and transform water and nitrogen. For instance, Agriculture (8) transfers nitrogen, from fertilizers that are imported to the system, to fresh food produce. Water supply (1) and Markets (7) provide drinking water and fresh produce to the city dwellers. These, represented as Households (2), then transform drinking water and food into wastewater and fecal sludge. While wastewater reaches Sewer and drainage (5) and is partially conveyed to the wastewater treatment plant (6), fecal sludge is discharged into OSS (3) and then either infiltrates the groundwater or is transported via truck to (6). Finally, the Landfill (4) receives all solid waste produced in the system.
To characterize the flows of nitrogen in the city, it is important to first calculate the flows of wastewater (WW WWTP ) and fecal sludge (FS WWTP ) that are treated at the wastewater treatment plant (WWTP). Both flows will be directly influenced by infrastructure projected in the Masterplan and are the main sources of nitrogen. This is done using Equations (1) and (2).
where f WWTP WW is the fraction [-] of wastewater that is conveyed to the WWTP; WW wet and WW dry are, respectively, the volumes of wastewater flowing in the sewer network during the wet and dry periods of the year m 3 year −1 .
where FS OSS is the volume of fecal sludge that reaches OSS m 3 year −1 ; and f Fill OSS , f Emptied OSS , f Formal OSS are, respectively, the fraction of OSS that fill up, are emptied, and are emptied by formal operators [-]. Finally, f WWTP FS is the fraction of fecal sludge that is transported to the WWTP [-]. Then, to quantify the amount of nitrogen reaching groundwater (N GW ), Equation (3) was used.
where f N−GW is the fraction of N sourced from the different sources that reaches groundwater [-]; N GrW is the amount of nitrogen sourced from greywater (Equation (4)); N LF is the amount of nitrogen in the city's landfill leachate (Equation (5)); N OAD is sourced from open-air defecation (Equation (6)), N OSS is the amount of nitrogen originating in OSS (Equation (7)); N Leak is the amount of nitrogen in the drinking water lost due to leakage (Equation (8)); and N WWTP is the amount of nitrogen originating in wastewater exfiltrated from the sewer infrastructure (Equation (9)). All these variables have units of [kg year −1 ].
where N conc GrW is the nitrogen concentration in greywater mg N L −1 ; Pop is the city's population [capita]; GR OSS GrW is the greywater generation rate for the population with OSS, and GR Sewer GrW is the greywater generation rate for the population covered by the sewer network, both in L capita −1 day −1 ; f OSS Pop is the fraction of the population with OSS; f Sewer Pop is the fraction of the population served by the sewer network [-]; and f seep is the fraction of wastewater that seeps for the sewer network [-].
where N conc LF is the nitrogen concentration in the landfill's leachate mg N L −1 , A LF is the area of the landfill m 2 , PP Maputo is the average annual precipitation in Maputo mm year −1 ≡ L·m −2 year −1 , and f rain−leachate is the fraction of the rain that becomes leachate [-].
where N cont Excreta is the nitrogen concentration in human excreta mg N kg Excreta −1 , Excreta is the daily excreta production kg Excreta capita −1 day −1 .
where N cont FS is the nitrogen concentration in fecal sludge mg N kg FS −1 ; and FS I,A,NF is the mass of fecal sludge in OSS that is informally emptied, abandoned, or that does not fill up kg FS day −1 . Please note that the nitrogen sourced from OSS, as defined in Equation (7), takes only into consideration nitrogen content in fecal sludge; nitrogen reaching groundwater and sourced from greywater discharged into OSS is calculated using Equation (4).
where N cont DW is the nitrogen concentration in drinking water mg N L −1 ; and Leak is the leakage in the water supply network m 3 year −1 .
where N cont Eff is the nitrogen concentration in the WWTP's effluent mg N L −1 ; Eff is the effluent flow at the wastewater treatment plant m 3 year −1 . This component was modelled as a best-case scenario: the nitrogen concentration from seepage is the same as the nitrogen concentration at the WWTP outlet.
For all calculations, SIMBOX, a program developed at Eawag, was used. SIMBOX can partially overcome the problems of limited data availability and uncertainty by defining the value of input parameters as probability distributions instead of point values [55]. With SIMBOX, each input parameter has a set of attributes, average, standard deviation, and distribution. These parameters are used to calculate output variables using the aforementioned equations.
In the program there are five options for parameter distribution: normal, truncated normal, lognormal, truncated lognormal, or uniform [36,56]. Uniform distribution was adopted for parameters for which little information was available; this included authors' assumptions, validated through discussions with local experts, and secondary data obtained from other geographical regions (e.g., dietary data obtained for other African countries). Truncated normal distribution was used for all fractions used in the calculations, with the exception of the ones that were assigned a uniform distribution; the use of truncated normal distributions avoids obtaining negative values for calculated output variables that cannot be negative (e.g., flows). Lognormal distribution was used for parameter concentrations, as these have the tendency to be non-negative and positively skewed [57]. All other parameters were assumed to follow normal distribution. A list with the most relevant input parameters is given in Table 1.
used for all fractions used in the calculations, with the exception of the ones that were assigned a uniform distribution; the use of truncated normal distributions avoids obtaining negative values for calculated output variables that cannot be negative (e.g., flows). Lognormal distribution was used for parameter concentrations, as these have the tendency to be non-negative and positively skewed [57]. All other parameters were assumed to follow normal distribution. A list with the most relevant input parameters is given in Table 1.   Once the uncertainties of input parameter values were defined, the uncertainties in output variables were calculated using Monte Carlo simulation [70]. Monte Carlo simulation uses random parameter values that are within the error margin previously defined for the input parameter values to calculate a distribution of probable output values [52]. For this work, each scenario was run 1000 times.
Sensitivity analysis, which consists in quantifying the effect of a defined input parameter change on an output variable [57], was also evaluated. This analysis enables identifying which input parameters exert a significant influence on a given output variable. To determine the sensitivity of a certain variable, this variable's value was calculated based on an initial set of input parameter values; afterwards, the value of one of the parameters was altered by 100%, with all other parameters being kept unaltered. The difference between the initial variable mean value and the value obtained after changing the input parameter was then analyzed. The procedure was repeated for all parameters influencing the variable of interest [57]. SIMBOX also displays the ranking of parameter knowledge, which consists of a list of ranked input parameters with respect to their impact on the uncertainty of the variable calculated. Based on this model, the flows of water and nitrogen in Maputo for the year of 2015 were calculated. Most input data was obtained from interviews with local water and sanitation experts and extracted from city-wide surveys that covered around 1200 households each [47,48]. It is important to underline that the survey [48] was specifically designed to provide input data for this model, in particular, fecal sludge flows in the city. While the first survey aimed at characterizing the spatial distribution of OSS systems in Maputo, the second studied attitudes and practices regarding the management of such systems. Data that could not be collected was calculated from surrogate parameters; e.g., the amount of fecal sludge that is not transported to the WWTP (6) was calculated considering the Mozambican diet, the average excreta generation in developing countries, and the rate of filled-up, replaced, or OSS systems that are (informally) emptied in Maputo. The sources for the most relevant parameters are given in Table 1.
In this manuscript, data are presented as box plots depicting median (middle of the data set), upper and lower quartiles (25% of data greater/lower), and minimum and maximum values; outliers are presented as circles.

Current Situation
According to recent surveys [47], the sewer system covers around 5% of Maputo's population, with most of the collected wastewater being discharged untreated into the sea. In addition, about 1% of the population practice open-air defecation, which happens throughout the city [47,48]. In the model, the overarching assumptions were:

1.
Only fecal sludge that is emptied by formal operators is transported to the wastewater treatment plant, 2.
Fecal sludge is homogeneously produced throughout the city, and 3.
The fraction of nitrogen that infiltrates the ground and reaches groundwater sources is the same, irrespective of the source.

Development Scenarios
Three development scenarios were formulated for this study. All scenarios represent improvements to the status quo and the impact of each scenario was quantified through the variation of input parameters that are related to the access to sanitation services and to its usage ( Table 2). The underlying assumptions regarding behavioral changes and management practices made for each scenario were discussed with local water and sanitation experts.
Scenario MP (from Masterplan) was formulated based on the Masterplan. Firstly, it aims at expanding the existing sewer system to reach around 14% of the population. This scenario includes the installation of more than 250 km of new sewers and the construction of five pumping stations and six WWTPs. Secondly, there is a focus on improving the quality of the actual fecal sludge management services (projected 86% of the population), with the construction of four fecal sludge treatment plants and 24 fecal sludge transfer stations [71]. Under this scenario, the Masterplan does not include plans for improvements to the domestic systems, therefore we assumed that its breakdown among the city's population will not change. Furthermore, we assumed that the construction and improvement of fecal sludge infrastructures and the increased access to efficient fecal sludge management services along the supply chain will allow more affordable fecal sludge management services and consequently raise the fraction of OSS systems that are emptied by formal operators by 25% (Table 2).
Scenario MPS (from Masterplan plus sewer) consists of the (unplanned) expansion of the sewer network to about 40% of the population, which would allow also covering two of the poorest urban districts of Maputo, Nhlamankulo, and KaMaxaquene, which are the districts with the highest population density in the city (almost 200 hab km −2 ). These are also the urban districts with the highest demand for emptying services where, currently, more than 40% is done unsafely (by hand with buckets) by informal operators [47]. For this scenario we have also assumed a 25% increase in the number of OSS systems that will be emptied and emptied by formal operators (Table 2).
Scenario MPO (from Masterplan plus OSS) consists of the (not yet planned) improvement of OSS. Over 72% of the OSS systems in Maputo are reported not to fill up [47], which means that currently only a small fraction of these systems work properly and are (formally) emptied (Table 2). A possible path for development could be to partially replace the existing systems (e.g., various types of latrines that represent more than 50% of the total [47]) with contained systems that require emptying. For this scenario we assumed that this replacement will occur homogenously throughout the city and that there will be an increase of 25% in the amount of systems that fill up, which will lead to an increase in the fractions of systems (formally) emptied ( Table 2). 1 % of the total population covered by sewer network; 2 % of the total population covered by onsite sanitation systems; 3 As % of the total; 4 As % of the systems that fill up; 5 As % of the systems that are emptied; 6 In the current situation, 1% of the population practices open-air defecation.

Fecal Sludge and Wastewater Flows
All three development scenarios will lead to an increased volume of wastewater ( Figure 2) and fecal sludge (Figure 3) conveyed to the WWTP (Table 3), when compared to the current situation.

Mass Flows
Currently, almost 500 tonnes of nitrogen reach groundwater each year (Figure 4). According to our results, Scenario MP could potentially result in a slight reduction (14%) of the nitrogen flow into the groundwater, Scenario MPS allows a 29% reduction, and Scenario MPO, a 19% reduction ( Table  2).
The sources of nitrogen reaching groundwater, as calculated by the model, are depicted in Figure 5. These include greywater discharged unsafely onto streets, a common practice among households relying on OSS [47,48], and directly into OSS systems, common among septic tank owners; the leachate from the city's landfill; leakage from the drinking water supply systems; open-air defecation; fecal sludge disposed of into OSS systems; and exfiltration from the sewer infrastructure, including network and WWTP. OSS infiltration is for all scenarios the largest source of nitrogen reaching groundwater, followed by nitrogen by leakage from the city's drinking water networks. The impact of nitrogen in leakage is related to large volumes of water lost each day in the

Mass Flows
Currently, almost 500 tonnes of nitrogen reach groundwater each year (Figure 4). According to our results, Scenario MP could potentially result in a slight reduction (14%) of the nitrogen flow into the groundwater, Scenario MPS allows a 29% reduction, and Scenario MPO, a 19% reduction ( Table 2).
The sources of nitrogen reaching groundwater, as calculated by the model, are depicted in Figure 5. These include greywater discharged unsafely onto streets, a common practice among households relying on OSS [47,48], and directly into OSS systems, common among septic tank owners; the leachate from the city's landfill; leakage from the drinking water supply systems; open-air defecation; fecal sludge disposed of into OSS systems; and exfiltration from the sewer infrastructure, including network and WWTP. OSS infiltration is for all scenarios the largest source of nitrogen reaching groundwater, followed by nitrogen by leakage from the city's drinking water networks. The impact of nitrogen in leakage is related to large volumes of water lost each day in the city's network (around 50% of 200,000 m 3 ).

Sensitivity Analysis and Parameter Knowledge Ranking
The sensitivity analysis for nitrogen reaching groundwater, fecal sludge transported, and wastewater conveyed to the WWTP is depicted in Figure 6. The corresponding change for each variable's output value is given for a 100% increase in several input parameters. Only parameters that lead to a variation of at least 30% are presented.

Sensitivity Analysis and Parameter Knowledge Ranking
The sensitivity analysis for nitrogen reaching groundwater, fecal sludge transported, and wastewater conveyed to the WWTP is depicted in Figure 6. The corresponding change for each variable's output value is given for a 100% increase in several input parameters. Only parameters that lead to a variation of at least 30% are presented.

Sensitivity Analysis and Parameter Knowledge Ranking
The sensitivity analysis for nitrogen reaching groundwater, fecal sludge transported, and wastewater conveyed to the WWTP is depicted in Figure 6. The corresponding change for each variable's output value is given for a 100% increase in several input parameters. Only parameters that lead to a variation of at least 30% are presented. Regarding fecal sludge, this variable is most sensitive to the household practices regarding septic tank management, in particular the selection of formal or informal operators. The volume of wastewater conveyed to the treatment plant is most sensitive to the number of city dwellers and to the amount of rain that reaches the network, given by the fraction of the year that is dry or wet. Finally, the amount of nitrogen that reaches groundwater is most sensitive to the fraction of nitrogen that reaches groundwater from all nitrogen discharged into the soil, to the size of the population and the food consumed by the population which will influence the content of nitrogen in fecal sludge. Population size is the only parameter that has a significant impact on the output of the three variables.
Regarding parameter knowledge and consequent uncertainty of the output results, when the uncertainty of a certain variable is too large, the analysis shows which parameter uncertainties would have to be reduced (knowledge gathering) in order to achieve a lower uncertainty [72]. For both fecal sludge and wastewater, most uncertainty is explained by population size, emptying practices of septic tanks, the amount of water used for flushing, and greywater generation rate. Finally, for nitrogen reaching groundwater, most uncertainty comes from the percentage of nitrogen that reaches the soil and finally groundwater sources and from the type of food used by the city's population. These parameters will thus have to be studied in detail in order to refine the model's output-produce less output uncertainty. It should also be mentioned that despite "nitrogen reaching groundwater" being very sensitive to Population, this parameter explains very little uncertainty on the output (approximately 8%; Table 4). In fact, Population is one of the best known input parameters to the model as this comes from a recent house-to-house census. This means that this parameter has little uncertainty and has little impact on the total uncertainty of "nitrogen reaching groundwater".  Regarding fecal sludge, this variable is most sensitive to the household practices regarding septic tank management, in particular the selection of formal or informal operators. The volume of wastewater conveyed to the treatment plant is most sensitive to the number of city dwellers and to the amount of rain that reaches the network, given by the fraction of the year that is dry or wet. Finally, the amount of nitrogen that reaches groundwater is most sensitive to the fraction of nitrogen that reaches groundwater from all nitrogen discharged into the soil, to the size of the population and the food consumed by the population which will influence the content of nitrogen in fecal sludge. Population size is the only parameter that has a significant impact on the output of the three variables.
Regarding parameter knowledge and consequent uncertainty of the output results, when the uncertainty of a certain variable is too large, the analysis shows which parameter uncertainties would have to be reduced (knowledge gathering) in order to achieve a lower uncertainty [72]. For both fecal sludge and wastewater, most uncertainty is explained by population size, emptying practices of septic tanks, the amount of water used for flushing, and greywater generation rate. Finally, for nitrogen reaching groundwater, most uncertainty comes from the percentage of nitrogen that reaches the soil and finally groundwater sources and from the type of food used by the city's population. These parameters will thus have to be studied in detail in order to refine the model's output-produce less output uncertainty. It should also be mentioned that despite "nitrogen reaching groundwater" being very sensitive to Population, this parameter explains very little uncertainty on the output (approximately 8%; Table 4). In fact, Population is one of the best known input parameters to the model as this comes from a recent house-to-house census. This means that this parameter has little uncertainty and has little impact on the total uncertainty of "nitrogen reaching groundwater".

Discussion
The flow of fecal sludge that could reach the WWTP, calculated with our model, has been overestimated (2000 m 3 day −1 , Figure 3). In fact, truck-counting campaigns at the WWTP set the actual transported fecal sludge to approximately four times less, about 500 m 3 day −1 [47,48]. We hypothesize one possible reason for this discrepancy: in the surveys [47,48], the households were asked "Did you select a formal or informal operator to empty your OSS?", with the subsequent assumption that formal operators always legally discharge at the WWTP, which might not be the case. On the other hand, estimations of wastewater conveyed to the WWTP (2300 m 3 day −1 ) were in the same order of magnitude of field measurements done by the authors 4000-5000 m 3 day −1 .
All three development scenarios could result in an improvement when compared to the current situation, with almost all wastewater produced in the city being conveyed to the WWTP, adequately treated, consequently reducing the levels of pollution in Maputo Bay. Also, the volumes of fecal sludge transported to the treatment plants would greatly increase. We expect that the increase in wastewater and fecal sludge, adequately treated, would translate into improved urban health and wellbeing for the city dwellers [46].
According to our calculations, OSS is the largest source of nitrogen reaching groundwater, a conclusion similar to that of other authors [35,36]. In the current situation, throughout the city, almost 500 tonnes of nitrogen reach the groundwater sources each year, a situation that leads to groundwater contamination by nitrate. We assume that this flow will have larger impacts in neighborhoods with higher population densities and where the population mostly relies on OSS [53]. Nevertheless, a comment is due to be made regarding a few of the underlying assumptions made in this manuscript and, in particular, the fraction of nitrogen that reaches groundwater from OSS. Nevertheless, it has been shown that groundwater recharge rate and hydraulic head have a significant influence on nitrate concentrations in aquifers [22]. The aquifer around Maputo shows a recharge rate of around 0.0002 to 0.0003 m 3 m −2 day −1 and a hydraulic head gradient in the range of 0.0015 to 0.002 [51]. Both parameters are an order of magnitude lower than the values assumed by [22] as "best case condition"-lower potential for nitrogen reaching groundwater sources. Finally, according to [51], no information concerning hydraulic conductivity of the aquitard unit is currently available for the area. Despite the lack of data, we argue that from a hydrology perspective, the potential for groundwater contamination by nitrate in the city of Maputo is high.
According to our results, Scenario MP would have the smallest impact on nitrogen reaching groundwater (14% reduction when compared to the current situation). Scenario MPS, on the other hand, which encompasses the expansion of the existing sewer network to a restricted and densely populated area, could potentially yield better results (29% reduction). Classic sewer network construction/expansion is often considered prohibitively expensive in the context of countries in the Global South [73]. However, an option for expanding the sewer network in peri-urban Maputo are simplified sewers. [74] argues that "values for sanitation technologies can vary widely and should be treated with caution and costs are always site specific depending on topography, housing layout, customer choice, materials used and population density". However, the author continues, "on-site sanitation and simplified sewerage are invariably cheaper than conventional sewerage". Finally, the author goes on stating that "for peri-urban areas, economies of scale can make simplified sewerage cheaper than even on-site sanitation technologies, in terms of total annual costs per household, at population densities greater than 160 hab km −2 ", giving the example of a case study in Brazil. In Maputo, the urban districts, Nhlamankulo and KaMaxaquene, both with almost 200 hab km −2 , are the only districts in the city above that threshold and also the ones with the highest demand of emptying services [47].
A second alternative development scenario (Scenario MPO) includes improvements on domestic OSS, an approach that is mostly based on fecal sludge management services and with potential for a 19% reduction in nitrogen reaching groundwater. Scenario MPO assumes the replacement of part of the existing OSS systems that are reported not to fill-up with systems that fill up and allow emptying.
The true impact of the baseline improvements in the three Scenarios, which are connected to the expansion of fecal sludge infrastructure and services, is not only difficult to make, but might have been overestimated. One of the underlying assumptions in all three scenarios is that there will be an increased request of formal operators due to the construction of decentralized fecal sludge infrastructure, which would translate into more fecal sludge being transported to the treatment plant and adequately treated. However, the impact of fecal sludge infrastructure depends greatly on behavioral changes and subsidization schemes [75] that were not planned or budgeted in the existing version of the Masterplan.

Conclusions
Cities throughout the Global South struggle to offer its dwellers safe and inclusive sanitation services. This is also the case in Maputo, capital of Mozambique, where most of the fecal sludge generated is not treated [47,48]. In part due to this situation, more than half of the city's population is supplied with groundwater where nitrate concentrations as high as 250 mg L −1 have been reported [50]. To change this situation, a city-wide sanitation Masterplan, proposing improvements to sanitation and storm water drainage infrastructure, was commissioned.
In this work, we quantified the amount of nitrogen reaching groundwater each year, and up to which extent the Masterplan would be an improvement. Using Material Flow Analysis (MFA), we concluded that each year more than 500 tonnes of nitrogen, mostly from OSS systems, reach the city's groundwater sources.
We demonstrated that the current version of the Masterplan falls short on improving the ongoing groundwater contamination by nitrate and proposed simplified sewer networks instead of only relying on OSS. Such approach not only demonstrates the added value of MFA as a planning and decision support tool, proving that such tools should be exploited more often to characterize the impact of planned urban sanitation in the Global South, but also adds to the argument of [76], according to whom "more attention is needed to how the planning process of sanitation infrastructure in cities across the Global South is designed and conducted". In preliminary studies for such city-wide plans, when sanitation infrastructure is absent, external factors, such as environmental and health indicators, should also be taken into consideration.
These results will be the starting point for a cost-benefit analysis that will take into consideration not only the costs of initial investments and the costs of maintenance, but also direct health benefits through improved sanitation, and external factors such as environment protection and indirect health risk, through groundwater consumption, to further support local decision-making. Furthermore, a better characterization of the area, including more geo-hydrological data, and bringing together the results of our simulations and climatological, soil, and land-use information [77] could further lead to insightful conclusions about the situation in Maputo.