Consumption of Free Chlorine in an Aqueduct Scheme with Low Protection : Case Study of the New Aqueduct Simbrivio-Castelli ( NASC ) , Italy

The safety of high quality drinking water supply relies on the quantities to be delivered, on the complexity of the water supply systems, and on the widespread phenomena of the contamination of water bodies. These parameters indicate the need for the development of an application that will allow the quick acquisition of data on strategic management. This is requires both the analysis of factors related to the hydraulic operation of the plants and the characteristics of water quality. The present paper aims to evaluate the use of models that predict data for water quality in a distribution system. The assessment is made in order to consider the use of the model as a support tool for the management system of a supply network and to optimize the quality of the provided service. The improvement of the control system related to the operations of disinfection, in particular, in the case of long pipelines, is absolutely mandatory in order to ensure the safety of public health and respect for the environment at high levels.


Introduction
In water supply network systems, when water is supplied from poorly protected catchment areas or where a proper treatment phase is not provided, periodic or accidental degradation of drinking water quality is often observed.The monitoring of water quality is therefore intended to play a proactive role.The spread of systems based on online analytical techniques for the measurement of the quality of different piping system parameters, combined with efficient computational models, can enable the prediction of the characteristics of the water in a fast and accurate way and thus provide data not only to circumvent the supply of water of low quality but, more importantly, improve the

Issues Related to the Use of Chlorine in Water Supply Networks
The quality of the water distributed through the supply network is influenced by different factors and the transformation processes that occur during the treatment process and the distribution in the piping [8,11].
The following issues have great interest: (a) the consumption of residual disinfectant in the network; (b) the formation of toxic substances from the reaction of the oxidants used in the disinfection; (c) the establishment of microbial growth processes in the network.

Consumption of Free Chlorine in the Network
The presence of a residual concentration of free chlorine (hypochlorous acid (HOCl) and ion hypochlorite (OCl − )) in treated water is indicative of an efficient disinfection (i.e., inefficient destruction of pathogenic microorganisms).The technical literature indicates that concentration values of residual free chlorine in the range 0.2-0.5 ppm ensure a good level of protection [12].Following conditions of quality degradation in the water transported downstream of the treatment, a reduction in the activity of the disinfectant (chlorine) decay may occur.This is due to different reactions that lead to the reduction of free forms of chlorine from both chemical and biological transformations of material present in suspended form in the liquid phase or immobilized on the solid-liquid interface in contact with the walls of the pipes and tanks, interfering with the fouling and corrosion processes.
The reaction of chlorine in water follows the reaction shown below: Hydrochloric acid dissociates in turn to form hydrogen and chloride ions Water 2018, 10, 127 3 of 13 Hypochlorous acid, however, dissociates only partially Several studies have been conducted to evaluate the factors that influence the process of consumption of free chlorine [13][14][15].A non-negligible influence of corrosion phenomena, in particular in pipes of cast iron, was also detected [16].

Final Products of Chlorine Reactions
The oxidation of the soluble natural organic fraction (NOM) present in the water favors the formation of organohalogen compounds.Chlorine can effectively convert humic substances to trihalomethanes (THMs), and other oxidation by-products of organohalogen compounds under the reaction conditions which occur in the liquid phase within the treatment systems [17].
In order to identify the organochlorine compounds produced in the disinfection phase, certain parameters, with which different fractions of halogenated compounds are detected, have been introduced.The main ones are the total organic halides (TOX), the adsorbable organic halides (AOX), and the extractable organic halides (EOXs).The choice is made on the basis of the possibility of evaluating the parameter, its correspondence to the physical phenomenon that is to be analyzed, and the ease of repeating the analysis in other contexts [18].

Microbial Growth
All types of waters, even those that come from deep aquifers practically well protected, contain some microorganisms, which continuously grow, even at relatively low temperatures [15].Low levels of dissolved oxygen (DO) and significant amounts of soluble iron and sulfur compounds were found in some groundwaters [19].These conditions are favorable to the growth of unwanted bacterial species, such as iron bacteria and sulfur bacteria [20,21].
The microorganisms already present in the water introduced into the supply network, even in small amounts, can promote the phenomena of microbial regrowth in pipes when the water has a sufficient content of organic substrate and nutrients (e.g., phosphorus), and the microorganisms are in the presence of DO, and when there are the thermal conditions stimulating bacterial replication [22,23].These conditions generally are present in all types of pipes [24].
In order to assess the presence of organic matter in drinking water treatment units and water supply systems that can be metabolized by microorganisms, Biodegradable Dissolved Organic Carbon (BDOC) and the Assimilable Organic Carbon (AOC) were used as indicators.Even if other possible alternatives might have been available, these represented the simple options that could be used in order to give an operative approach to the research.The AOC is the fraction of the organic carbon quickly usable by microorganisms and is generally a fraction of the total organic carbon (TOC) in variable percentage values between 0.1% and 9% [25].The BDOC is the fraction of the Dissolved Organic Carbon (DOC) which can be metabolized by bacteria within a period of a few days to few months, which is of particular importance for drinking water, as it can be the starting point of bacteria regrowth leading to the development of undesirable higher organisms [26].The results of a survey conducted in the United States have shown that the ratio between BDOC and DOC in water distribution networks can have values greater than 0.1 mg/L.In river water, that value may increase up to 0.5 mg/L [27,28].

Data Analysis Section-Modeling of the Chlorine Consumption with Semi-Empirical Methods
The calculation models that simulate chlorine decay are valuable tools for the management of chlorine residuals and for reproducing the evolution of quality parameters in the network.They can assist water network managers in identifying the correct action to be pursued [29][30][31][32][33][34].The models typically operate with two distinct modules: the first provides the calculation of characteristic quantities Water 2018, 10, 127 4 of 13 of the hydraulic motion processes and the second one, on the basis of the reconstruction of the velocity field, analyzes the changes in the concentrations of contaminants, taking into account diffusive phenomena, mixing, and decay.In particular, the process of propagation of contaminants is represented on the basis of a scheme consisting of elements or volumetric finite segments in which the pipeline was divided or analyzing calculation blocks of variable sizes that move along the pipeline.
EPANET is one of the most widely used computing models [35,36] and this is the reason why this software was chosen.It is a software that enables the modelling of multiple constituents throughout the water distribution system (WDS) [37].The chemical reactions between the constituents are described by differential-algebraic equations and are solved by numerical methods.In the present study, the consumption of the simulated substance, i.e., chlorine, in the network is described through the development of two processes that develop in the liquid phase and at the liquid-solid contact with the inner walls of the pipeline.The reaction rate is set with a first-order model, which assumes that the chlorine concentration decreases exponentially: where: C, chlorine concentration; t, time; k b constant reaction rate, in h −1 or day −1 .
The k b values were estimated for distribution networks [22,38]: the values vary from 0.01 h −1 to 0.3 h −1 .For the determination of k b as a function of temperature and concentration of the organic carbon, Kiene et al. (1998) [38] proposed the following experimental procedure: Temperature between 5 • C and 25 • C and TOC concentration between 1 and 3 mgC L −1 : where: T, temperature in Kelvin degrees; C TOC , concentration of TOC in mgC L −1 ; a and b, coefficients: a = 1.8 × 10 6 and b = 6050.The kinetics of the reactions that occur at the solid-liquid interface of the walls of the pipes are influenced by the specific surface and the process of mass transfer between the liquid medium and the solid phase [39].The consumption of chlorine on the walls of the duct is determined by applying first-order kinetics.The rate of reaction is expressed by the following relationship: where: R, radius of the pipe; k w , constant reaction rate to the wall The k w value was determined according to pipe roughness by using the following equation: where: C, is the roughness coefficient used in the Hazen-Williams formula, F, is a calculation coefficient.

of 13
Applying F values between 1.5 and 15, k w varied between 0.001 and 0.01 m day −1 .The determination of k f is carried out in accordance to the characteristics of the analyzed species, by means of the reactivity conditions, the molecular diffusivity, and the hydrodynamic conditions of the flow, using the Reynolds number [40].

Characterization of the Distribution Network for Water Supply
The study aims to analyze the process of consumption of free chlorine in the New Aqueduct Simbrivio-Castelli (NASC), which supplies drinking water to about sixty (60) municipalities and three industrial areas in an area located in Central Italy, serving a population of about 350,000 inhabitants.Figure 1 shows the quality index of those areas, by comparing the water quality data.
Water 2018, 10, x FOR PEER REVIEW 5 of 13

Characterization of the Distribution Network for Water Supply
The study aims to analyze the process of consumption of free chlorine in the New Aqueduct Simbrivio-Castelli (NASC), which supplies drinking water to about sixty (60) municipalities and three industrial areas in an area located in Central Italy, serving a population of about 350,000 inhabitants.Figure 1 shows the quality index of those areas, by comparing the water quality data.The flow conveyed to the distribution centers is about 0.7 m 3 s −1 .The pipelines extend for about 350 km; they are built mostly from steel.In some places there are pipelines made of gray and ductile iron.The diameters range from ND50 to ND800 (nominal diameter equal to 50 and 800 mm, respectively).The supply of the aqueduct, originally secured by groups of springs, has been integrated with the waters from other springs and groundwater sources.
The operative scheme is briefly described in the following points: 1. the water from the different sources is mixed in a longitudinal tank, called a "hydraulic gallery"; in this basin, which is the source of supply pipelines, sodium hypochlorite is added as a disinfection agent; 2. the main duct branches into numerous secondary pipelines.
The water scheme has been characterized according to its different elements: tanks, dividers, pipelines, lifting equipment.Pipes were characterized based upon their construction material, year of commissioning, length, and diameter.The data relating to the management of the water in the years 2008-2010 were collected with the following information considered: 1. flows from the different sources; 2. days of interruption of the supply service; 3. repair of hydraulic losses occurred on pipes; 4. interruptions in the disinfection procedure.The flow conveyed to the distribution centers is about 0.7 m 3 s −1 .The pipelines extend for about 350 km; they are built mostly from steel.In some places there are pipelines made of gray and ductile iron.The diameters range from ND50 to ND800 (nominal diameter equal to 50 and 800 mm, respectively).The supply of the aqueduct, originally secured by groups of springs, has been integrated with the waters from other springs and groundwater sources.
The operative scheme is briefly described in the following points: 1. the water from the different sources is mixed in a longitudinal tank, called a "hydraulic gallery"; in this basin, which is the source of supply pipelines, sodium hypochlorite is added as a disinfection agent; 2.
the main duct branches into numerous secondary pipelines.
The water scheme has been characterized according to its different elements: tanks, dividers, pipelines, lifting equipment.Pipes were characterized based upon their construction material, year of commissioning, length, and diameter.The data relating to the management of the water in the years 2008-2010 were collected with the following information considered: 1.
flows from the different sources; 2.
days of interruption of the supply service; 3.
repair of hydraulic losses occurred on pipes; 4.
interruptions in the disinfection procedure.

Collection of Data on Water Quality
The monitoring of water quality was conducted in 34 sampling stations.
The analyses relating to the water quality refer to a time period of two years, in which 1117 samples were analyzed.The samples were divided into the following classes: The methods of analysis comply with the Unichim Manual, which describes the water analysis procedures [41].The technical standards adopted are those used as reference by State where the case study was conducted.

Method of Processing Data
In the first phase, the water quality time series data (temperature, pH, nitrates, turbidity, Kubel oxidability, free chlorine, total and fecal coliforms, fecal streptococci colonies on agar at 22 • C and 36 • C) were reconstructed and the main statistical parameters were calculated (Figures 2 and 3).
In the second phase, a monitoring network, constituted by a reduced number of sampling stations, was selected, in accordance with appropriate access to the sampling points.On the basis of availability and significance of the data, and taking into account the location along one of the main pipelines and the distance between the stations, in the order of 15-25 km, seven sampling stations were selected.
Water 2018, 10, x FOR PEER REVIEW 6 of 13

Collection of Data on Water Quality
The monitoring of water quality was conducted in 34 sampling stations.
The analyses relating to the water quality refer to a time period of two years, in which 1117 samples were analyzed.The samples were divided into the following classes: The methods of analysis comply with the Unichim Manual, which describes the water analysis procedures [41].The technical standards adopted are those used as reference by the State where the case study was conducted.

Method of Processing Data
In the first phase, the water quality time series data (temperature, pH, nitrates, turbidity, Kubel oxidability, free chlorine, total and fecal coliforms, fecal streptococci colonies on agar at 22 °C and 36 °C) were reconstructed and the main statistical parameters were calculated (Figures 2 and 3).
In the second phase, a monitoring network, constituted by a reduced number of sampling stations, was selected, in accordance with appropriate access to the sampling points.On the basis of availability and significance of the data, and taking into account the location along one of the main pipelines and the distance between the stations, in the order of 15-25 km, seven sampling stations were selected.

Processing of Free Chlorine Data
The data were measured over two years (June 2010-June 2012).Table 1 shows the values of the main statistical parameters, and Figure 4 shows the trends of chlorine concentrations and ccumulative trends observed in the survey period in the sampling stations.

Processing of Free Chlorine Data
The data were measured over two years (June 2010-June 2012).Table 1 shows the values of the main statistical parameters, and Figure 4 shows the trends of chlorine concentrations and ccumulative trends observed in the survey period in the sampling stations.

Simulations with EPANET Model
The last phase of the survey consisted of the reconstruction of the field of free chlorine concentration using the 2.0 EPANET model [35].The simulations were carried out considering an ordinary hydraulic regime, in stationary conditions, inside the supply network [42][43][44].The hydraulic layout is constituted by: two tanks working with constant hydraulic load (tanks receiving water from springs and wells); four tanks in variable hydraulic load (hydraulic gallery and hydraulic dividers); 107 nodes (dividers in pressure and flow delivery points); 105 pipes with an internal diameter ranging between 50 mm to 800 mm; 11 water pumping systems.The delivered flow rate was set equal to the annual average value.
The pressure losses were calculated with the formula of Hazen-Williams by varying the coefficients of roughness between 130 and 135 for cast-iron pipes and between 140 and 150 for steel pipelines [35].The flow rates delivered to the nodes were set equal to the annual average values of flow at the delivery points of municipal and cooperatives utilities.The trend of the concentration of chlorine was studied by applying a first-order kinetic model for the reactions in the liquid medium and at the walls of the pipe.The values of kb were made to vary between 0.1 and 0.9 day −1 .The values of the coefficient kw were determined according to the roughness of the pipeline and kw ranged between 0.001 and 0.01 m day −1 .Table 2 shows the hydraulic residence times in the network, calculated in the sampling sections.Figure 5

Simulations with EPANET Model
The last phase of the survey consisted of the reconstruction of the field of free chlorine concentration using the 2.0 EPANET model [35].The simulations were carried out considering an ordinary hydraulic regime, in stationary conditions, inside the supply network [42][43][44].The hydraulic layout is constituted by: two tanks working with constant hydraulic load (tanks receiving water from springs and wells); four tanks in variable hydraulic load (hydraulic gallery and hydraulic dividers); 107 nodes (dividers in pressure and flow delivery points); 105 pipes with an internal diameter ranging between 50 mm to 800 mm; 11 water pumping systems.The delivered flow rate was set equal to the annual average value.
The pressure losses were calculated with the formula of Hazen-Williams by varying the coefficients of roughness between 130 and 135 for cast-iron pipes and between 140 and 150 for steel pipelines [35].The flow rates delivered to the nodes were set equal to the annual average values of flow at the delivery Water 2018, 10, 127 9 of 13 points of municipal and cooperatives utilities.The trend of the concentration of chlorine was studied by applying a first-order kinetic model for the reactions in the liquid medium and at the walls of the pipe.The values of k b were made to vary between 0.1 and 0.9 day −1 .The values of the coefficient k w were determined according to the roughness of the pipeline and k w ranged between 0.001 and 0.01 m day −1 .Table 2 shows the hydraulic residence times in the network, calculated in the sampling sections.Figure 5 illustrates the comparison between the chlorine concentration values calculated in the sampling sections with the model compared with the statistical variables related to the measured data.The model provides a good fit and a good capacity to describe the degradative processes active in the aqueduct, in particular in hydraulic steady conditions.The diagram is constructed by putting the values of the hydraulic residence time in the network, and the chlorine concentration on the abscissa axis and on the ordinate axis, respectively.
Further evidence for the good fit in terms of agreement between the observed and calculated data can also be observed in Figure 6.The diagram is constructed by putting the values of the hydraulic residence time in the network, and the chlorine concentration on the abscissa axis and on the ordinate axis, respectively.
Further evidence for the good fit in terms of agreement between the observed and calculated data can also be observed in Figure 6.The diagram is constructed by putting the values of the hydraulic residence time in the network, and the chlorine concentration on the abscissa axis and on the ordinate axis, respectively.
Further evidence for the good fit in terms of agreement between the observed and calculated data can also be observed in Figure 6.

Discussion
In more than 70% of the observed data from the seven selected stations, the processing of the observed aforementioned data showed that, in the period of investigation, the concentration values were below 0.2 mgCl 2 /L, which is the recommended minimum value at the delivery point to the user, according to the current regulations.Furthermore, the distribution of the average values of the concentration of free chlorine showed a decreasing trend with increasing distance from the disinfection section.A similar trend is observed for other statistical parameters such as the median and the 70th percentile.The consumption of free chlorine in the network can be attributed to the degradation phenomena which occurred periodically under certain conditions for one of the water supply sources.This kind of condition cannot affect a judgment of potability and involves the need to implement a new chlorination phase in the basins of municipal and cooperative utilities in order to reach compliance with the value indicated in the current legislation.The reduction of free chlorine is evident in the pipes closest to the origin of the aqueduct and in the pipes with smaller diameters characterized by more extended running times.
In the first case, such a situation can be attributed to the presence of substances which are rapidly oxidized and therefore consume chlorine.This means that the chlorine demand has not been covered adequately.In the second case, the phenomenon is certainly due to local conditions.
The applied calculation model EPANET showed a good ability to represent the processes of quality degradation in the pipes working in hydraulic stationary conditions.In particular, the constant rate of consumption in the liquid medium, k b , and the constant rate of chlorine consumption at the solid-liquid contact with the pipeline walls, k w , were estimated.The simulation results showed very good agreement with the experimental data, with k b values varying between 0.1 and 0.3 day −1 combined with values of k w varying between 0.01 and 0.07 m day −1 .These values determined for the examined water networks are lower than those indicated in the technical literature referring to the distribution networks.This consideration demonstrates that in most cases the unfolding of the phenomena of the consumption of free chlorine in water results less intense.This final consideration makes it even more evident how important the issue is, especially with the aim of preserving public health.Obviously, the obtained results should be verified in other detection campaigns and the application should be repeated in other similar contexts.

Concluding Remarks
From the results of the present study, certain conclusions can be drawn:

•
In more than 70% of the observed data the concentration values were below 0.2 mgCl 2 /L.

•
The reduction of free chlorine is evident in the pipes closest to the origin of the aqueduct and in the pipes with smaller diameters characterized by more extended running times.

•
The applied calculation model EPANET showed a good ability to represent the processes of quality degradation in the pipes working in hydraulic stationary conditions.• k b values varied between 0.1 and 0.3 day −1 and k w varied between 0.01 and 0.07 m day −1 .

•
The experience should be repeated in other similar contexts, with future studies also trying to use some additional key parameters, in order to test the generalizability of this study's findings to other similar situations.

Figure 1 .
Figure 1.Quality index of the supply network located in Central Italy.

Figure 1 .
Figure 1.Quality index of the supply network located in Central Italy.
(a) aqueduct water supply (springs and wells); (b) water procured to the aqueduct network; (c) disinfected water supplied to the supply network; (d) water transported by the supply network to the consumers.
(a) aqueduct water supply (springs and wells); (b) water procured to the aqueduct network; (c) disinfected water supplied to the supply network; (d) water transported by the supply network to the consumers.

Figure 2 .
Figure 2. Water quality (temperature, pH, nitrates, turbidity, Kubel oxidability, free chlorine, total and fecal coliforms, fecal streptococci colonies on agar at 22 °C and 36 °C) observed in sampling stations (one after the other).

Figure 2 .
Figure 2. Water quality (temperature, pH, nitrates, turbidity, Kubel oxidability, free chlorine, total and fecal coliforms, fecal streptococci colonies on agar at 22 • C and 36 • C) observed in sampling stations (one after the other).

Figure 3 .
Figure 3. Cumulative frequency of the chlorine concentration (Mg Cl2/L) observed in sampling stations.

Figure 5 .
Figure 5. Values of chlorine calculated by the model.

Figure 5 .
Figure 5. Values of chlorine calculated by the model.

Figure 5 .
Figure 5. Values of chlorine calculated by the model.

Figure 6 .
Figure 6.Comparison between calculated values and statistical parameters of the observed data.

1 KbFigure 6 .
Figure 6.Comparison between calculated values and statistical parameters of the observed data.

Table 1 .
Number of samples and statistical parameters related to the observed concentration of free chlorine.
Figure 3. Cumulative frequency of the chlorine concentration (Mg Cl 2 /L) observed in sampling stations.

Table 1 .
Number of samples and statistical parameters related to the observed concentration of free chlorine.

Table 2 .
Simulation results: hydraulic residence times in the network calculated in the sampling stations.