First, five municipalities were selected to identify the potential for using treated WW. These study sites had the municipal public sewage collection and treatment services that generate reclaimed water for non-potable purposes. In addition, population growth was projected in the selected municipalities until the year 2035 in order to estimate the investment necessary to universalise access to the sewage service. Therefore, it was possible to obtain the period of return of this investment from the collection of the tariff paid by the users. In addition to the income collected by the service provider, in the situation where water reuse was assumed, an additional increase in the income of the service providers was estimated from the sale of recycled WW. Thus, the economic feasibility analysis aims to elucidate whether additional investment to promote a better quality of reclaimed water is feasible in view of the potential return on this investment from the commercialisation of this resource by service providers.
Therefore, variations in the investment required to expand access to sewage treatment services for the population by 2035 were assumed, mainly because of the proposed level of treatment, which can be secondary or tertiary (with disinfection systems). Consequently, along with the variation in the level of sewage treatment, the potential risk of microbiological contamination for users, operators, and the general population during the application and use of reclaimed WW was also important. For this reason, a quantitative microbiological risk assessment methodology was applied to identify the risk associated with contamination in each condition analysed. Therefore, methodologies of quantitative microbiological risk assessment and net present value were applied in a complementary manner to assess financial viability.
2.1 Selected municipalities as a case study and the study boundaries
The five most populous municipalities in MS State were used as study sites (Table 1).
Table 1
The most populous municipalities of Mato Grosso do Sul and their respective rankings in terms of gross value added at competing industry prices
Municipality | Population (inhabitants)a | Current ranking in terms of gross value added participation at current industry prices relative to the total state GDP (%)b | Sewage collection and treatment index (%)c |
---|
Campo Grande | 885,711 | 1 | 73.0 |
Dourados | 215,486 | 4 | 51.4 |
Três Lagoas | 119,465 | 2 | 37.3 |
Corumbá | 110,806 | 5 | 45.0 |
Ponta Porã | 91,082 | 10 | 27.6 |
Source: aIBGE (2017), bIBGE (2014), cANA (2017). |
The final uses for the treated WWs considered in the study were street cleaning, sewage network unblocking, crop irrigation, and industrial activities. These scenarios were compared with the scenario of discharging treated WW to water bodies. Figure 1 shows the analysis of water reuse from a financial perspective considering the risk associated with microbial contamination.
2.2. Economic feasibility study
The net present value (NPV) methodology was applied to revenues obtained under different investment scenarios from income generated by service user tariffs and the sale of reuse water by local concessionaires. Apart from assuming operation and maintenance costs, the NPV considered the investment required to universalise access to WW treatment services in the selected municipalities, as shown in Eq. 1:
\(NPV=\sum _{y=0}^{n}\frac{{CF}_{y}}{{(1+i)}^{y}}-{C}_{0}\) Eq. 1
where y is the year, CF is cash flow, i is the discount rate, and C 0 is the total initial investment cost. We assume a minimum attractiveness rate of 12%. The costs and benefits of the modelled interventions were calculated in the Brazilian currency Reais (R$), based on the January 2020 exchange rate, assuming an intervention period of 15 years, starting in 2020 and ending in 2035. All costs and benefits incurred after 2020 were converted to 2020-equivalent values using a discount rate of 2.5%.
Four panoramas of investment reflecting variations in WW treatment services were analysed in addition to the variations in the sensitivity analysis. In the first condition (baseline), receipts for the WW treatment were obtained solely from the collection of annual tariffs applied to users of the WW service. In the second, the annual tariff revenue was supplemented by the medium possible increment in annual revenue from the sale of treated WW, based on the demand for medium quality and medium cost water. Finally, in the third and fourth conditions, the annual tariff revenue was supplemented by the minimum and maximum possible revenue from low-to-high-quality water, respectively. The cost of the treated WW was based on the study of Sabesp (2017) and ranged from 1.50–5.08 R$/m3 depending on the required quality.
The estimated revenue of sewage services was calculated using data on billed sewage volumes from the Brazilian National Sanitation Information System (SNIS) from 2012–2017, considering projection trends up to 2035 including the revenue from the reuse of water, as shown in Eq. 2:
\(RT=\left(WTR\right)+\left(VR\times TR\right)\) Eq. 2
where RT is the total revenue from sanitary sewage, WTR is the WW treatment service revenue, VR is the volume of reused water billed, and TR is the average water reuse rate. Increases in the projected population served by the sanitation service over the duration of the project and readjustments to tariffs led to variations in the total annual receipts collected by the service provider. In the Online Resource 1, it is possible to verify all data considered in each municipality from 2013–2035, including estimated population, increase in WW treatment coverage by WWTPs (not including septic tanks), tariffs of WW treatment, operating revenue from WW treatment, operating expenditure from WW treatment, and water reuse volume according to different applications.
2.3. Recycled WW demand and supply
We estimated the water reuse demand in each municipality and conducted a sensitivity analysis considering a variation of ± 20% in the estimated average amount of water (Khan and Jain 2018; Hernandez-Sancho and Molinos-Senante 2015).
2.3.1. Street cleaning
We assumed that only 20% of the total streets paved per municipality would be cleaned using treated WW, as this practice is carried out in addition to simple sweeping and is only necessary on streets with high flows of people or economic activity, producing an excess of dirtiness. To estimate the volume of recycled water in each municipality necessary to meet this demand, SNIS data were used to determine the number of households supplied with the water distribution network. Data on the per-link extent of the water network from SNIS were also considered to obtain an estimate of the total length of the paved streets. It was assumed that the average width of the streets was 8 m, with a cleaning efficiency of approximately 100 m/m3, and a weekly street washing frequency.
Street cleaning was assumed to be performed by a team of four employees, including three workers and a driver, who wore all the necessary individual protection equipment (gloves, boots, proper clothing, and helmets). Materials, such as brooms and liquid detergents, were also used in the washing process. It was assumed that the treated WW was transported by tank trucks. Street cleanings should be scheduled at times when movements of people were minimal.
2.3.2. Sewage network cleaning
The need for unblocking the sewage network was assessed using information obtained from the United States Environmental Protection Agency (EPA 1999). This study adopted a method of simply discharging treated WW into the network at pressure under the control of a water truck (at a pump head of approximately 60 mca). The procedures adopted for clearing sewage networks are similar to those adopted for street cleaning, in which a tanker truck is used with employees wearing personal protective equipment.
Considering that the frequency of unblocking can vary from twice a month to once every six months, we assumed an average frequency of sewage network unblocking of once a month, over only 5% of the network length in each municipality. Based on the fact that for diameters of up to 500 mm, the average demand volume of treated WW is equivalent to three times the linear volume of the pipe, we assumed a required volume to remove obstructions of 2.4 m3 for each linear meter of sewage network, with a useful volume of approximately 0.8 m3. Note that these figures consider the growth of the sewage network of each municipality over time, and therefore show increases in the proportion of the demand for treated WW for cleaning the network.
2.3.3. Industrial activities
The system model used in this study considered a closed circuit for industrial water reuse. The model starts with the final effluent of an existing municipal WWTP, and the endpoint is the industry that receives this effluent from a water tanker. When the water tanker enters the facility, it is connected to a storage reservoir to supply the recycled water, which is later pumped to its final site of use. The pump is operated and maintained to ensure minimal risks for the workers.
We estimated the demand for treated WW in industry in the selected municipalities using data from SNIS on volumes of produced water, assuming that approximately 17% of the water consumed in Brazil is allotted to industry (FAO 2019), and that an average of 30% of this water supply could be replaced by treated WW. It is worth mentioning that SNIS data refer to water produced and distributed by public water supplies, and do not include other routes through which water can be supplied to the industry.
2.3.4. Crop irrigation
The agricultural sector is the largest water consumer globally, and consumes approximately 60% of the freshwater in Brazil (FAO 2019). We estimated the demand for treated WW in agriculture in the selected municipalities using data from SNIS on volumes of produced water, assuming that approximately 60% of the water consumed in Brazil is allotted to agriculture (FAO 2019), and that 30% of this amount could be used from treated WW.
2.4. Projection of sewage treatment demand
We analysed the demand for WW treatment in selected municipalities in MS State (Campo Grande, Dourados, Três Lagoas, Corumbá, and Ponta Porã) based on the population growth of each municipality for the year 2035. Therefore, we developed several investment conditions under the assumption that the entire urban population will be served by sewage treatment systems by 2035, with 90% of the population served by a centralised system and 10% served by individual solutions, such as septic tanks.
We applied mathematical and statistical methods to identify growth trends in the sanitation service industry. We estimated parameters such as tariffs, as well as operational and maintenance costs through 2035 by extrapolating historical data (covering 2012–2017) from the SNIS database. Appropriate parameters were obtained by applying the linear regression tool in Excel (Microsoft Corporation 2018) directly to data values (for linear growth) or to their logarithms (for exponential growth). Exponential growth curves were appropriate for the economic parameters, as tariffs and costs tend to keep pace with inflation, and therefore exhibit geometric growth. To calculate operational expenses for WW collection and treatment services, the total expenditure on both sewage and water supply services in the SNIS database was considered, as well as the percentage of participation of each service (water and sewage) in relation to the total operating revenue. Because the population growth rate in the state is currently decreasing, an exponential model of population growth would have been inappropriate. Instead, a linear model was applied and calibrated to match data from the last IBGE (2010) census with the 2035 population predictions produced by ANA (2017) as the starting and ending points, respectively.
The linear and exponential representations of the parameter growth are given in Equations 3 and 4, respectively:
\(Value=\left(Cst+rate\right) \times A \left(linear growth\right)\) Eq. 3
\(Value={e}^{\left(Cst+rate\right)A} \left(exponential growth\right)\) Eq. 4
where Value is the projected parameter, Cst is the historical value of the parameter, rate is the function-determined linear or exponential growth rate, and A is the year. These equations were used to model the costs of tariffs, operation, and maintenance, as well as the growth of population and sewage production for each year from 2020–2035.
2.5. Expansion of access to sanitation services
Costs for the expansion of access to WW collection and treatment were obtained from ANA (2017). These costs are related to conventional sewage systems, where raw effluent is collected separately from rainwater/drainage, and its treatment must achieve the standards for discharge into water bodies.
Tertiary treatment of WW for reuse is a common treatment level where close contact with the water is considered a possibility (Vojtěchovská Šrámková et al. 2018; Voulvoulis 2018). Thus, to guarantee the required quality, additional costs of 30% (based on the estimates by von Sperling 2016) were considered when introducing disinfection systems to decrease microbial contamination. It should be noted that some systems have lower investment and operating costs than others, with higher treatment levels generally corresponding to higher costs. In addition, per capita costs will decrease with increasing installed capacity, as a result of economies of scale. Finally, several local factors can influence system costs, including the quality of the materials used in construction, local personnel costs, energy costs, and potential land acquisition costs (depending on whether the land is owned by the government).
As the accuracy of the cost estimates in each locality was difficult to access, an additional variation of ± 20% in the cost of treatment for reuse was assumed in our sensitivity analysis (in addition to 30%). Table 2 lists the treatment systems used in each municipality and their respective characteristics (Sanesul 2017; Águas Guariroba 2017).
Table 2
Technologies used for the treatment systems in each selected municipality and their current respective characteristics
Wastewater treatment plant | Population project (inhab.) | Expected flowb (l/s) | Treatment process |
Três Lagoas WWTP 1 Planalto | 72,000 | 100 | UASB + physico-chemical |
Três Lagoas WWTP 2 Jupiá | 28,800 | 80 | UASB + Biodruma + Secondary decanter |
Corumbá WWTP 1 Olaria | 57,600 | 80 | UASB + Trickling filter + Secondary decanter |
Corumbá WWTP 2 Maria Leite | 57,600 | 80 | UASB + Trickling filter + Secondary decanter |
Dourados WWTP 1 Guaxini | 86,400 | 120 | UASB + Trickling filter + Secondary decanter |
Dourados WWTP 2 Água Boa | 79,200 | 110 | UASB + Trickling filter + Secondary decanter |
Dourados WWTP 3 Laranja Doce | 28,800 | 40 | UASB + Trickling filter + Secondary decanter |
Dourados WWTP 4 Harry Amorim | 28,800 | 40 | UASB + Facultative pond |
Ponta Porã WWTP 1 Estoril | 57,600 | 80 | UASB + Trickling filter + Secondary decanter |
Ponta Porã WWTP 2 São Thomaz | 28,800 | 40 | UASB + Trickling filter + Secondary decanter |
Campo Grande WWTP 1 Imbirussú | 67,000 | 120 | Compact aerobic system + Unitank |
Campo Grande WWTP 2 Los Angeles | 500,000 | 900 | UASB + contact tank for disinfection |
aConsisting of an anaerobic system, followed by an aerobic system with rotating microbial support type Tubular Reel; bDesign parameter: maximum expected flow when the system is in full operation. |
2.6. Microbial risk assessment
The risks associated with the use of treated WW for the aforementioned purposes were assessed using the quantitative microbial risk assessment (QMRA) approach. This involves the following four key steps.
2.6.1. Hazard identification
Escherichia coli 0157 was chosen as the pathogen for the QMRA. Infection with this pathogen from exposure to WW has been extensively reported (Forslund et al. 2012; Akiba et al. 2015; Tripathi et al. 2019). Data on E. coli concentrations in different scenarios can be found in the Online Resource 1. Data were collected monthly for a monitoring period of 18 months. Data from WWTPs in Campo Grande were obtained from Aristimunho (2019). The concentrations of E. coli in the final effluent from all the selected WWTPs ranged from 104 to 106. This was assumed to be the concentration before tertiary treatment, as observed in the treatment plants (see Online Resource 1).
2.6.2. Exposure assessment
Exposure to WW was considered based on the different uses of WW, currently or projected, within the various municipalities considered for this study. Exposure to treated WW was considered in six scenarios (Fig. 1). The main routes through which the exposed populations could ingest pathogens in treated WW are described in Table 3. For instance, during WW collection, risks were assessed for workers in the major exposure group. The routes of exposure were through inhalation of aerosols generated during collection, direct ingestion of aerosol droplets, and hand-to-mouth exposure (from contaminated hands). The duration of exposure was assumed to be 30 min during this stage. Similar to effluent collection, workers were the major exposed group during the transport of the treated WW. The exposure routes and durations were the same as those during the WW collection stage.
During the application of the treated WW for the different uses mentioned above (Fig. 1), two population groups were considered: the workers involved in the practice, and the public, who may be exposed during these practices. For the workers, the main routes of exposure considered were the same as mentioned above, with a duration of an hour. However, the exposure of the public was mainly through the route of inhalation of droplets for a maximum of 1 min. Application of WW for crop irrigation was also considered, with farmers as the main exposed population. Additionally, the risks of infection with the pathogen for consumers of farm produce were also considered. In this assessment, lettuce consumption was used to model the risk of infection for consumers. Table 3 presents the volume of water ingested used in the risk assessment.
Table 3
Volume of water ingested under different exposure scenarios
Exposure Scenario | Route of Exposure | Volume of effluent exposure |
Scenario I: Discharge into rivers | Direct ingestion during recreation (swimming) | 1–5 mLa |
Scenario II: Exposure to effluents during collection | Aerosol inhalation | Calculated using Equations 6–8 and input data from Table 6 |
| Ingestion of droplets |
| Hand-to-mouth ingestion |
Scenario III: Exposure during transport of effluents | Ingestion of droplets | Calculated using Equations 6–7 and input data from Table 6 |
| Hand-to-mouth ingestion | |
Scenario IV: Street and sewer network cleaning | Aerosol Inhalation for workers | Calculated using Equations 6–8 and input data from Table 6 |
| Ingestion of droplets for workers |
| Hand-to-mouth ingestion for workers |
| Aerosol inhalation by the general public (street cleaning) | Calculated using Equations 6–8 and input data from Table 6 |
Scenario V: Use of effluents for industrial activities | Aerosol inhalation | Calculated using Equations 6–8 and input data from Table 6 |
| Droplet ingestion |
| Hand-to-mouth ingestion |
| Aerosol inhalation (irrigation of green spaces) |
| Direct ingestion during irrigation (farmers) | 1–5 mLb |
| Consumption of vegetables (consumers) | Calculated using Equations 9 and input data from Table 7 |
aBased on assumptions; bWHO (2006). |
To determine the dose of pathogenic E. coli ingested during these different scenarios, Eq. 5 was used, as follows:
\(D=Craw\times V\) Eq. 5
Where “D” is the dose ingested, “Craw” the concentration of pathogenic E. coli per milliliter and “V” the volume (mL/day) ingested.
Equations 6–8 were used to determine the volume of WW ingested per person under each of the exposure scenarios.
Hand-to-mouth ingestion from wet hands:
\(QHM = h \times A \times fHM\) Eq. 6
where h is the thickness of the water film on the hands (mm), A is the skin surface area in contact with the mouth (mm2), and fHM is the frequency of the hand-to-mouth contact (n/min).
Ingestion of droplets of water:
\(QD = VD \times fD\) Eq. 7
where VD is the water droplet volume (µL) and fD is the frequency at which droplets splash into the mouth (n/min).
Inhaled volume of water per minute:
\(QI= IR \times VIWS\) Eq. 8
where IR is the air inhalation rate (m3/min) and VIWS is the fraction of inhalable water spray (µL water/m3 air). The values of the parameters used in Equations 7–9 are listed in Table 4.
Table 4
Additional data for QMRA calculations
Parameter | Value or distribution of values | Source |
Inhalation rate (m3/min) | | |
Children | Uniform (1.11E-02, 4.36E-02) | |
Adults | Uniform (1.03E-02, 7.77E-02) | USEPA 2011 |
VIWS, volume of inhalable water spray (µL/m3) | Average: 10.8 95% confidence interval: 1.76–36.3 | de Man et al. 2014 |
H, film thickness of water on hands (mm) | Uniform (1.97E-02, 2.34E-02) | USEPA 2011 |
A, surface area of hand that is mouthed (mm2) | Uniform (100, 2,000) | USEPA 2011 |
Fhm frequency of hand-to-mouth contact, (n/min) | Gamma (2, 0.5) | Freeman et al. 2001 |
Fd, frequency of water droplets landing in mouth (n/min) | Gamma (2.1, 0.17) | de Man et al. 2014 |
Vd, volume of a droplet (µL) | Uniform (0.5, 524) | de Man et al. 2014 |
The risk of infection from the indirect use of water, due to the consumption of vegetables irrigated with this water, was determined based on the dose of pathogenic E. coli ingested by consumers. This was modelled with lettuce as a surrogate vegetable, using the following formula:
\(DC=V\times l\times c\) Eq. 9
Where “V” is the volume of water caught on the lettuce in mL/g of lettuce, “l” the mean per capita intake of lettuce in grams per person per day, and “c” the concentration of pathogenic E. coli in the water used for irrigation. The different exposure scenarios and ingested volumes are listed in Table 5.
Table 5
Assumptions used in estimation of risks of consuming vegetables irrigated with recycled WW
Exposure scenario /Assumptions for dosage | Volume of water ingested (mL or g) | Frequency (days) | Reference |
Consumption of lettuce | | Uniform distribution (156,160) | |
Volume of water caught on lettuce | Normal distribution (0.108, 0.019) | | Hamilton et al. (2006) |
Per capita intake of lettuce | Pert distribution (25, 50, 75) | | Sant’Ana et al. (2014) |
To model the risk of infection from the discharge of the treated sewage into rivers, a baseline E. coli concentration of 1 colony forming unit/100 mL was assumed for the river water quality prior to discharge. Based on the assumption that workflow discharge was carried out at a fixed point, the reduction in the rate of auto-depurification was calculated as a function of time to determine the risk of infection in individuals swimming in the river between zero and five days after the release of the effluent (Eq. 10).
\(Nt= {N}_{0} \times {e}^{(-kt)}\) Eq. 10
where Nt is the coliform density at time t in running water (most probable number 100 mL− 1), N0 is the initial density of coliforms after dilution, and k is the first-order activation constant (d− 1).
2.6.3. Dose response assessment
Similarities in exposure risk arose from the fact that each type of activity assumed the same route of exposure and pathogen (pathogenic E. coli) concentration. The volumes of water ingested by direct contact were obtained from the 2006 WHO publication on WW reuse, and were assumed to be uniformly distributed over a range of 1–5 mL. Most of the E. coli load was disregarded, with only 8% assumed pathogenic. The dose of pathogenic E. coli ingested under each scenario considered in this study was calculated based on the volume of water consumed during these activities. We also assessed the potential reduction in the risk of infection when tertiary WW treatment was implemented. We modelled the risks with a 3 Log (99.9%) reduction in E. coli concentration after possible tertiary treatment (disinfection) with chlorination, ozonation, or UV radiation (von Sperling 2005). Thus, we assessed the risks of exposure to secondary-and tertiary-treated WW separately. The beta Poisson dose-response model was used for risk assessment (Haas et al. 2014), represented by the following equation:
\({p\left(d\right)=1-\left(1+\left(\frac{d}{{N}_{50}}\right)\left({2}^{\frac{1}{\alpha }}-1\right)\right)}^{-\alpha }\) Eq. 11
where “p(d)” is the probability of infection, “d” is the median infection dose representing the number of organisms that will infect 50% of the exposed population (N50), and α is the dimensionless infectivity constant. For this assessment, the risks were modelled for E. coli 0157 with an N50 value of 2.11 × 106, and α was set equal to 1.55 × 10− 1 (DuPont et al. 1971; Girardi et al. 2019).
2.6.4. Risk characterisation
Risks of infection from multiple exposures were determined using Eq. 12:
\(\text{P1(A)=1-}{\left(1-\text{P1}\left(d\right)\right)}^{n}\) Eq. 12
where ‘P1(d)’ is the risk of infection from a single exposure to a dose ‘d’ of the bacteria; and ‘n’ is the number of days of exposure to the single dose ‘d’(Sakaji and Funamizu 1998). The duration of exposure (n) for all scenarios was considered to be 30 min, except for the exposure to aerosols by the general public, which was modelled for 1 min.
Determination of harm by using the disability-adjusted life year (DALY) metric
The DALY metric can be used to widely compare different illnesses and other risks from daily life (Havelaar et al. 2000). The DALY aims to estimate the overall environmental burden of disease (Knol et al. 2009) and measures population health (Vocht et al. 2011) by estimating the loss of healthy life years of the population (Havelaar et al. 2000; Murray 1994). DALYs are calculated by considering the years of life lost due to premature death or mortality and years of life lived with a disability (Timm et al. 2016). In this study, the DALYs per year were calculated using Eq. 13:
\(DALY= \sum _{i=1}^{n}P\left(ill\backslash inf\right)\times P({outcome}_{i}\backslash ill)\times {Duration}_{i }\times {Severity}_{i}\) Eq. 13
where “n” refers to the total number of outcomes considered. Three outcomes were considered in this study, and all of them were considered to be consequences of diarrhoea. P(ill\inf) is the probability of illness given an infection, as we assessed the risks for E. coli 0157 to be 1. P(outcome\ill) provides the probability of an outcome given an illness. Durationi is the duration (years) of outcome i, and Severityi is the severity weight for outcome i. The severity, weight, and duration of the disability or disease used in this assessment were adjusted to developing country conditions (Katukiza et al. 2014). The severity weight scale ranged from 0 (healthy) to 1 (death) (Havelaar and Melse 2003; Katukiza et al. 2014). For each pathogen, the disease outcomes, duration, severity, and disease burden per case of infection were obtained from literature. An average life expectancy of 61 years was used for this study, and the years lost were based on death occurring at the age of 1 year (Howard and Pedley 2004). All input data taken from the literature for the DALY estimation are presented in Table 6. All the microbial risk (QMRA) and burden of disease (DALY) assessments were performed with Monte Carlo simulations with 10,000 iterations using the @Risk (Palisade, USA) add-on to Excel (Microsoft Corporation 2018).
Table 6
Health outcomes and the related probability inputs used in the DALY assessment
Relevant Health outcome | Probability of outcome (P(outcome\ill)) | Severity weight | Duration (Years) |
Watery diarrhoea | 0.53 | 0.067 | 0.009 |
Bloody diarrhoea | 0.47 | 0.39 | 0.015 |
Death | 0.007 | 1 | 60 |
*Input values taken from Havelaar and Melse (2003), Machdar et al. (2013), Katukiza et al. (2014), and Howard and Pedley (2004). |