This paper presents knowledge-based planning of the most appropriate treatment schemes for drinking water supply (DWS) systems, jointly applying statistical tools and experimental tests. Milan City was chosen as the case study, its DWS system being composed of more than 20 DWS units, widespread in the urban area. First, multivariate statistical techniques (factor analysis and cluster analysis) were applied to groundwater monitoring data, to identify specific contaminations of the captured groundwater and their spatial distribution. In detail, eight typical compositions of captured groundwater and their distribution in DWS units were identified, leading to the selection of different treatment processes to be adopted in DWS units. Then, process performances were experimentally evaluated on representative groundwater samples. Heterotrophic denitrification was tested at pilot-scale on groundwater contaminated by nitrate, identifying optimal operating and design parameters for complying with regulation limits. Granular activated carbon (GAC) adsorption was tested at laboratory-scale for the removal of volatile organic compounds (VOCs) and pesticides, present in various ratios in groundwater samples: removal efficiencies were obtained and the occurrence of competition phenomena was highlighted. Air stripping was tested at pilot-scale for removing the most volatile VOCs (trichloroethylene, tetrachloroethylene), highlighting its usefulness as pre-treatment before adsorption to increase GAC lifetime.
Drinking water supply (DWS) in industrialized and densely populated metropolitan areas is often assured through multiple DWS units that intercept various punctual and diffuse pollution plumes.
In such a complex context, planning the appropriate treatment trains can be a hard task. The knowledge of raw groundwater composition and its variability, in order to individuate the water parameters that need to be removed, is critical. The water quality parameters of interest being many and characterized by a complex spatial distribution, data-mining techniques may be extremely supportive (Olsen et al. 2012; Selle et al. 2013) in this perspective. In the literature, different statistical tools have been applied for the analysis of water quality data; among these, data driven techniques such as principal components analysis (PCA), factor analysis (FA) and cluster analysis (CA) are reported as effective and unbiased methods to extract meaningful information from a large dataset (Papaioannou et al. 2010). Despite the fact that the use of multivariate statistics for water quality assessment is widely reported in the literature (Papaioannou et al. 2010; Olsen et al. 2012; Page et al. 2012), its application as a supporting tool for the planning and management of DWS systems is not conventional.
The planning of a DWS system requires the selection of the appropriate treatment schemes, based on the comparison of the available technologies, considering technical and non-technical factors (Hamouda et al. 2009). Since knowledge acquisition is often case-specific, experimental and pilot studies are needed in order to obtain representative information about the technical and economic efficiency of the alternatives (Joksimovic et al. 2006).
The present work considers a case study of the Milan City DWS system. The system provides about 230×106 m3 y−1 of drinking water, collecting groundwater from three aquifers by 522 wells spread across the city area and organized in clusters of about 20 wells, each one conveying water to a DWS unit. This study applies a knowledge-based concept to the planning of the proper treatment train to have in the DWS units. As the first step of the knowledge-based concept, a multivariate statistical analysis was applied to the available groundwater monitoring data in order to identify specific contaminations of the captured groundwater and their spatial distribution. Then the second step of the knowledge framework, the experimental testing of some treatment processes, was performed at laboratory- and pilot-scale to assess the optimal operating conditions. Statistical analyses were used as a supporting tool to guide the selection of (1) the most appropriate treatment technologies and (2) the most representative groundwater samples to be tested, in order to obtain results that could be extended to a high number of cases. Specifically, as treatments are concerned, activated carbon (AC) adsorption, air stripping and heterotrophic biodenitrification were evaluated respectively to address the removal of organics (i.e. volatile organic compounds (VOCs) and pesticides) and nitrate.
HCA was run based on the FA extracted varifactors and Ward's method was used as the clustering method. All the data analysis was performed through the IBM SPSS Statistics 21 statistical package.
As far as the experimental setup is concerned, a series of laboratory-scale batch experiments was performed to assess organics removal in groundwater samples by adsorption on AC. A commercial AC (Cecacarbon GAC 1240) was tested through bottle-point isotherm measurements (ASTM D5919-96): different doses (from 1 to 35 mg L−1) of AC were dosed in 1 L pyrex bottles and kept mixed on a magnetic stirrer (400 RPM) at ambient temperature (20–25 °C) for 48 h before measuring the residual concentration of VOCs and pesticides in the liquid phase.
Air stripping was tested through a pilot-scale steel column (0.15 m diameter, 1.8 m height, 3.2 min HRT, 15.8 m3 m−2 h−1 air flow rate), analyzing inlet and outlet concentrations of VOCs and pesticides.
Finally, heterotrophic denitrification was tested in a pilot-scale biofilter consisting of a PVC column (0.33 m inner diameter, 2.40 m height, 21 min EBCT, 0.38 porosity) filled with expanded clay media (Biolite®). After a start-up period of about 50 days, aimed at allowing biomass colonization, the biofilter was fed up-flow with groundwater extracted by a supply well at a flow rate of 0.6 m3 h−1. Average influent concentration of NO3- during the experimentation was 8.9 ± 0.65 mg N L−1. The process was evaluated for 10 weeks by monitoring NO3-, NO2-, volatile suspended solids (VSS) and total organic carbon (TOC) in water at the biofilter inlet, outlet and at four intermediate biofilter heights. During this period, different operating conditions were tested: (i) nutrient availability was varied between 40% and 135% of the stoichiometric request, evaluated according to Henze et al. (2008), by supplying different doses of sodium acetate and potassium phosphate; (ii) cycle durations were varied (24, 48 and 72 h), by periodically backwashing the reactor in order to remove excess biomass.
Regarding analytical methods, the concentrations of VOCs, pesticides, nitrite and nitrate ions and VSS were measured according to EPA 524.3-09, APAT-IRSA 5060/03, EPA 300.1 and APHA 2540 methods, respectively; TOC was measured photometrically by a kit method (LCK 380, Hach Lange).
RESULTS AND DISCUSSION
FA was applied to the groundwater quality dataset in order to reduce the number of variables and to identify sets of parameters with the same pattern of variability. Starting from the 88 analyzed parameters, 16 relevant factors were extracted, explaining 57.6% of data variance and identifying different kinds of contamination in the studied area. Among factors summarizing parameters usually present in natural groundwater, such as physical–chemical parameters or alkaline metals, 13 of the extracted factors could be attributed to pollution of anthropic origin, such as fertilizer constituents (nitrate, sulphate and chloride), PAHs, pesticides, heavy metals (Sb, Cd, V, Se, As and Cr) and VOCs. So FA allowed the identification of the main pollution components affecting the groundwater captured in the studied DWS system.
Therefore FA and CA results provided useful indications for identifying the most proper technologies to be adopted locally for treating groundwater in the studied DWS system.
Nitrate, pesticides and VOCs were the species higher or very close to the national regulation limits so different types of treatments were experimentally tested to address the removal of these pollutants, particularly: (1) heterotrophic denitrification, for nitrate removal, (2) AC adsorption, for both pesticide and VOC removal, and (3) air stripping for the removal of the most volatile VOCs, i.e. trichloroethylene (TCE), tetrachloroethylene (PCE). Experiments were carried out on a limited number of groundwater samples, extracted from representative DWS units which were selected based on the results of the statistical analysis.
Groundwater belonging to cluster E showed a nitrate concentration above the Italian quality standard for drinking waters (i.e. 50 mg L−1). The suitability of a heterotrophic denitrification treatment was then experimentally assessed using a pilot-scale biofilter, evaluating the main aspects of the process: (1) nitrate removal efficiency, (2) residual nitrite concentration in the effluent, and (3) risk of bacterial regrowth in the distribution network triggered by the release of biomass and the availability of readily biodegradable organic carbon. All these aspects were evaluated with respect to the main process parameters, namely organic carbon availability and nitrogen volumetric load.
In order to assess the risk of bacterial regrowth in distribution networks, the release of organic carbon and biomass during pilot tests was also evaluated: residual TOC and VSS in the effluent were equal to respectively 4.6 ± 1.34 mg C L−1 and 3.2 ± 1.57 mg VSS L−1. No significant differences (Student's t-test, p-value > 0.05) among the tested carbon dosage conditions were found in TOC and VSS effluent concentrations.
The effect of the volumetric load on the denitrification process was evaluated considering the process efficiency along the biofilter height. Under C-limiting conditions, none of the tested BV provided nitrate removal efficiencies higher than 60%, resulting in the highest nitrite concentrations in the effluent (up to 6.2 mg N L−1). Under carbon oversupply conditions, removal efficiencies were strongly affected by the treatment cycle duration (tTR), which was modified by varying the biofilter backwashing frequency: for tTR of 24 and 72 h and volumetric loads between 0.64 and 1.07 kg N m-3 d-1, complete nitrate removal and nitrite concentrations below the regulation limit (0.5 mg L−1) were obtained; on the other hand, at tTR of 48 h, nitrate removal efficiencies >90% were obtained only for BV up to 0.80 kg N m-3 d-1. Moreover, nitrite concentrations in the effluent below the regulation limit were observed at tTR of 48 h only when BV was set to 0.64 kg N m-3 d-1. This behavior may probably be attributed to transitory or less stable conditions achieved at the intermediate tTR of 48 h; however, further investigations are needed to draw conclusions about this behavior.
Finally, no significant removal was observed for the three VOCs (PCE, TCE, vinyl chloride) and nine pesticides (2,6-dichlorbenzamide, 3,6-dichloropiridazine, atrazine, atrazine desethyl, atrazine desisopropyl, bromacile, hexazinone, terbuthylazine desethyl, tris-2(chloroethyl)-phosphate) which were initially present in the inlet water. Therefore, additional treatments, such as AC adsorption or air stripping, should be designed for groundwater requiring the removal of these organics.
These results suggested that adsorption onto AC is a suitable process for the removal of organic micropollutants. Nonetheless, when poorly adsorbable organics are present, an upgrade of the process is required, in order to increase the bed lifetime, for example, adopting higher EBCTs or a two-stage in-series configuration.
In some cases, especially for DWS units treating groundwater with high levels of VOCs (such as clusters C, D and F), air stripping could be a suitable solution both as a stand-alone treatment or as a pre-treatment before adsorption onto AC. Stripping efficiency was evaluated using a pilot-scale aeration column, fed on groundwater containing various VOCs, mainly PCE and TCE at average initial concentrations of 16.2 μg L-1 and 2.4 μg L−1, respectively. Significant removal was observed for PCE and TCE, on average 36% ± 11% and 18% ± 13%, respectively. PCE removal was higher than TCE removal, since a lower air–water ratio is required as expected by their values of Henry's constant. Moreover, results for TCE showed a significant correlation between removal efficiency and initial concentration (ρ = 0.713, p-value = 0.001, N = 10), while the same behavior was not observed for PCE. Results suggested that the adoption of an air-stripping step before GAC could be a valid solution, especially in the presence of PCE at high concentrations: as a matter of fact, the preliminary removal of PCE could reduce competition for adsorption sites, increasing the adsorption of the other organic micropollutants and reducing the frequency of GAC regeneration. However, the choice and the proper design of such a layout must also be based on economic evaluations, considering the costs related to both the treatment phases, depending on the specific water composition.
Data mining applied to groundwater quality data permitted the identification of the main kinds of contamination affecting the DWS system of Milan City, defining typical water composition contributing to each DWS unit. According to these results, different treatment processes were selected to be adopted in DWS units and their performances were experimentally evaluated. Experimental tests were performed on the most representative groundwater samples to obtain information generalizable to a high number of DWS units.
Heterotrophic denitrification resulted in being effective for nitrate abatement despite the significant organic contamination due to pesticides and VOCs. Carbon supply over 90% of the stoichiometric requirements and volumetric loads ≤1.07 kg N m-3 d-1 were found to assure effective and complete denitrification.
Adsorption onto AC was effective for the removal of several VOCs and pesticides. However, variable efficiencies were highlighted towards contaminants due to (1) their physical–chemical properties and (2) the influence of water composition, which magnifies competition phenomena. An air-stripping step before GAC resulted in a successful solution to removing the most volatile compounds (e.g. PCE), favoring the following adsorption of low-adsorbable compounds (e.g. TCE).
Reported results showed that data-mining techniques can effectively support the optimal design of a drinking water treatment train, allowing (1) identification of the main contaminants to be controlled and (2) reduction of efforts in testing the performances of the most suitable treatment processes, supporting the selection of the most representative water samples to be used in the experiments, in order to obtain results of general validity. The applied knowledge-based approach resulted in a suitable supporting tool for the planning of a DWS system in an urban area, such as Milan City, coping with a composite source water quality, a complex groundwater capturing system and many fingerprint constraints.
This research has been funded by Metropolitana Milanese S.p.A. (MM), which manages the integrated water service for the city of Milan (I).
The authors to thank Carlo Carrettini and Fabio Marelli for having made this study possible, and Angela Manenti, Chiara Pagano and the MM laboratory staff for their helpful cooperation during the experimental work.