Next Article in Journal
Distribution Characteristics of Phosphorus in the Yarlung Zangbo River Basin
Next Article in Special Issue
Developing a Statistical Model to Improve Drinking Water Quality for Water Distribution System by Minimizing Heavy Metal Releases
Previous Article in Journal
N2O Fluxes and Rates of Nitrification and Denitrification at the Sediment–Water Interface in Taihu Lake, China
Previous Article in Special Issue
PATs Operating in Water Networks under Unsteady Flow Conditions: Control Valve Manoeuvre and Overspeed Effect
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vulnerability Assessment to Trihalomethane Exposure in Water Distribution Systems

1
Department of Civil and Mechanical Engineering, University of Cassino and Southern Lazio, 03043 Cassino, Italy
2
Department of Civil Engineering and Architecture, University of Pavia, 27100 Pavia, Italy
*
Author to whom correspondence should be addressed.
Water 2018, 10(7), 912; https://doi.org/10.3390/w10070912
Submission received: 28 May 2018 / Revised: 29 June 2018 / Accepted: 6 July 2018 / Published: 10 July 2018
(This article belongs to the Special Issue Advances in Water Distribution Networks)

Abstract

:
Chlorination is an effective and cheap disinfectant for preventing waterborne diseases-causing microorganisms, but its compounds tend to react with the natural organic matter (NOM), forming potentially harmful and unwanted disinfection by-products (DBPs) such as trihalomethanes (THMs), haloacetic acids (HAAs), and others. The present paper proposes a methodology for estimating the vulnerability with respect to users’ exposure to DPBs in water distribution systems (WDSs). The presented application considers total THMs (TTHMs) concentration, but the methodology can be used also for other types of DPBs. Five vulnerability indexes are adopted that furnish different kinds of information about the exposure. The methodology is applied to five case studies, and the results suggest that the introduced indexes identify different critical areas in respect to elevated concentrations of TTHMs. In this way, the use of the proposed methodology allows identifying the higher risk nodes with respect to the different kinds of exposure, whether it is a short period of exposure to high TTHMs values, or chronic exposure to low concentrations. The application of the methodology furnishes useful information for an optimal WDS management, for planning system modifications and district sectorization taking into account water quality.

1. Introduction

The disinfection of drinking water increases users’ protection by reducing their risk of exposure to pathogenic contamination and limiting the proliferation of microbial species [1,2]. Among the various technologies available nowadays, chlorination is still the major worldwide disinfection strategy, due mainly to its effectiveness [3,4], low cost, and simplicity (easy to produce, store, transport, and use).
Despite the proved effectiveness of chlorine in preventing waterborne disease-causing microorganisms, its compounds tend to react with the natural organic matter (NOM) naturally occurring in the resource through reactions in bulk water or close to the pipe wall, unintentionally forming potentially harmful and unwanted disinfection by-products (DBPs). It has been observed that some DBPs have possible adverse health effects. Indeed, they are carcinogenic potential substances [5,6,7], and some epidemiologic studies have raised the issue of potential adverse reproductive effects [8,9,10]. The exposure via non-ingestion routes, such as inhalation and dermal contact, may also pose risks to human health [11,12]. For this reason, a good water distribution management with respect to chlorination requires the solution of two conflicting objectives: to maintain sufficiently high chlorine residual and low enough DPBs concentrations in the whole system.
In the presented study, a methodology to assess the vulnerability respect to the DBPs exposure in a water distribution system (WDS) is proposed. The vulnerability assessment is useful for an optimal system operation for many purposes, such as for the optimization of chlorine booster functioning, to design public health strategies to reduce risks [13,14], for the optimal position of water quality monitoring stations [15,16], to obtain indications for epidemiological investigations [17], or for the district sectorization and skeletonization of a network based on water quality considerations [18,19].
Over the last 30 years, numerous DBPs have been studied and classified, including trihalomethanes (THMs), haloacetic acids (HAAs), haloacetonitriles (HANs), and haloketones (HKs). Even if it is unclear which specific DBPs are responsible for the adverse health effects [20,21], only two classes of DBPs—THMs and HAA—are regulated. The four THMs species (chloroform, bromodichloromethane, dibromochloromethane, and bromoform) are produced from sodium hypochlorite reacting with the NOM. Even if the total THMs (TTHMs) concentration does not give information on the HAA levels and it does not consider the relative role of the four THMs species in contributing to adverse health outcomes, it is often used as a surrogate for the measurement of DBPs [22].
Guidelines such as those from the World Health Organization [23] suggest maintaining in the system a free chlorine residual concentration in the range of 0.2–0.5 mg/L, while the TTHMs concentration is bounded by an upper regulation limit that is different in each country. For example, the United States (US) standard for the TTHMs concentration is 80 µg/L, the Canadian [24] and the Australian–New Zealand [25] guidelines set the maximum value to 100 μg/L and 250 μg/L, respectively, while in Italy, the regulation is very restrictive, with the limit value equal to 30 μg/L.
In the presented application, the methodology is performed considered the TTHMs concentration, furnishing indications to water companies for respecting the regulation limit. However, by adopting the adequate formation model, it can be used for studying other DPBs.
In the existing literature, water age (represented by residence time of the water in the system), contaminants, chlorine, and DBPs concentrations (such as bromide, TTHMs, and HAAs) are the main elements considered to evaluate indexes representative of a network vulnerability in terms of water quality.
Many studies estimated the exposure of the users to contaminants using different indexes. For example, Propato and Uber and Murrey et al. [26,27] quantified the contaminated water ingested by individuals during an event, evaluating the volume of polluted water that was delivered to the users. Nilson et al. [28] considered as a vulnerability index the total mass contaminated load in each node, and performed the sensitivity of the nodal mass distribution to various system characteristics and stochastic demands representation. Other studies have considered the impacts of contamination events without quantifying the ingestion by individuals, such as in Khanal et al. [29], which introduced an exposure index based on the percentage of the users exposed. Davis and Janke [30] accounted for the influence of ingestion timing in addition to the volume, and in a second work [31], the same authors also presented a probabilistic model for timing ingestion based on data collected by the American Time Use Survey in order to advise more realistic exposures.
Thompson et al. [32] estimated the vulnerability of a demand node in terms of exposure to low residual chlorine concentration, while Baoyu et al. [33] presented a vulnerability assessment model that also took into account the impact of water age, as well as the uncertainty related to bulk and wall reaction coefficients, demand multipliers, and pipe roughness.
Other works performed vulnerability analyses to identify the areas of the system and the population with a higher risk respect to DBPs exposure. More recently, Islam et al. [34] used a non-compliance potential index predicted using TTHMs, haloacetic acids (HAAs), and free residual chlorine as variables. Successively, Islam et al. [35] presented an index-based approach to locate chlorine booster stations in WDSs. The maximization of a Water Quality Index (WQI) based on the formation of TTHMs is used as objective of the optimization problem, and the required number of the stations is recommended based on a trade-off analysis of risk potentials, WQI, and the life cycle cost of a booster.
To perform the vulnerability analysis with respect to the consumers’ exposure to TTHMs, it is fundamental to use models that are able to predict their formation within WDSs. Many mathematical approaches have been proposed that have been obtained by means of laboratory and field scaled experiments on raw and pre-treated water [11,36]. They usually account for the correlations existing between the formation of TTHMs and the elements that influence them, such as precursors, operational parameters (pH, contact time), and environmental conditions (temperature and season variability). Many predicting models are represented by empirical relationships between TTHMs concentrations and the main parameters affecting the TTHMs formation [37,38], while others are based on kinetic formulation that are, involved during chlorine bulk and wall reactions [39,40].
The present paper proposes a methodology for estimating the vulnerability respect to the TTHMs exposure in a WDS. Five vulnerability indexes furnishing different kinds of information about the exposure are adopted. The individually or aggregated use of these indexes allow identifying critical points in a WDS with respect to different kinds of exposure: high concentration exposure for a short period, or chronic exposure to low TTHMs concentrations, which are both important. The methodology is easy to apply and is useful for the water companies to evaluate the more critical zone in terms of water quality and also for analyzing future scenarios.
In this paper, the preliminary analysis performed on a literature network in Di Cristo et al. [41] is completed with further investigation, and extended to two more literature schemes and two more real water distribution systems.
The paper is structured in the following way: the following paragraph describes the methodology, and then the vulnerability analysis is presented for the five considered cases of study. Finally, some conclusions are drawn.

2. The Methodology

The five introduced indexes, which were evaluated in each node of the system, differ regarding the kind of information they can provide, and the way that they are formulated. The first considered parameter is the maximum TTHMs concentration value reached in each node during the entire simulation:
  [ T T H M s ] i m a x = m a x { [ T T H M s ] i , t }
where [TTHMs]i,t is the actual TTHMs concentration at time t in the node i.
The second one is the TTHMs average concentration CiTTHMs weighted on nodes demand, which is defined as:
  C i T T H M s = 0 T O T q i ( t ) · [ T T H M s ] i · d t 0 T O T q i ( t ) · d t
where TOT is the total observation time, qi(t) is the actual demand node (L/s), and i represents the node index.
Two other useful indexes are the dimensionless exposure time TiE and the normalized contaminated volume Vic,n, which are expressed by:
  T i E = 0 T O T I ( t ) · d t T O T
  V i c , n = 0 T O T q i ( t ) · I ( t ) · d t 0 T O T q i ( t ) · d t
where the factor I(t) is defined as:
  I ( t ) = { 1 i f [ T T H M s ] i , t > [ T T H M s ] l i m 0 i f [ T T H M s ] i , t [ T T H M s ] l i m
[TTHMs]lim is a fixed attention threshold. TiE and Vic,n, represent the percentage of the time and the demand distributed in presence of contaminated water in a node during the observation period, respectively. They can assume values between 1, when all of the water that is delivered is contaminated, and 0, which means that the threshold is never exceeded in the node. Note that with fixed base demands in the nodes, there are no differences between the TiE and Vic,n, values. Moreover, in a node, both values are equal to 1 if the limit concentration is exceeded for the entire simulation period.
The last index is the contaminated mass load MiTTHMs, which corresponds to the mass of TTHMs delivered in each demand node when their concentration is above the threshold [TTHMs]lim, and is formulated as:
  M i T T H M s = 0 T O T ( [ T T H M s ] i , t [ T T H M s ] l i m ) · q i ( t ) · d t
As already mentioned, the five indexes furnish different and complementary information. In particular, both the maximum instantaneous and the average TTHMs concentrations are evidence of a likely dangerous exposure, without considering number of users involved and the duration of the exposition. The dimensionless exposure time and the normalized contaminated volume consider the duration of the exposition to dangerous values, while the contaminated mass is also related to the number of users exposed. The last three indexes indicate if there is a chronicle or sporadic TTHMs exposure above the alert bound.
For the evaluation of the indexes defined by Equations (1)–(6), the chlorine concentration is computed using the first-order kinetic equation, and assuming that it is consumed only through bulk reaction:
  [ C l ] i , t = [ C l 2 ] i , 0 exp ( k b · t )
where [Cl2]i,0 (mg/L) and [Cl2]i,t (mg/L) are the chlorine concentration initially and at time t for the node i, respectively; and kb (1/h) is the bulk chlorine decay coefficient.
The TTHMs concentrations are computed as a linear function of consumed chlorine [37,42], using the following equation:
  [ T T H M s ] i , t = [ T T H M s ] i , 0 + D · ( [ C l 2 ] i , 0 [ C l 2 ] i , t )
where [TTHMs]i,0 (μg/L) and [TTHMs]i,t (μg/L) are the TTHMs concentrations initially and at time t in the node i, respectively; and D (μg/mg) is the TTHMs yield coefficient, which is defined as the ratio of μg of TTHMs formed to mg of chlorine consumed.
Both the chlorine decay and the TTHMs formation coefficients have to be properly assigned. Practically, they can be fixed considering literature suggestions or calibrated using measured data.
Hydraulic and water quality simulations are performed using the EPANET [43] and EPANET Multi-Species Extension (EPANET-MSX) codes. EPANET-MSX [44,45] enables the evaluation of chlorine and TTHMs concentrations through just providing the formulations governing the reaction dynamics. Using the EPANET and EPANET-MSX toolkits working simultaneously, the whole methodology is implemented in a code written in the C++ programming language.

3. Results

The proposed methodology is applied to five case studies (Figure 1): the example network Net-3 introduced in the EPANET User’s Manual, the PW06 network by Prasad and Walter [46], the Anytown water distribution system [47], the main trunk of AVWSS (Aurunci-Valcanneto Water Supply System) located in Lazio region, Italy [36], and the real network of Cimitile, which is located in South Italy.
In all of the cases, water is supposed to be disinfected through sodium hypochlorite. For the reaction rate coefficients to model chlorine decay, kb (Equation (7)) and TTHMs formation, D (Equation (8)), unique global values are assumed, because not significant variations of the environmental conditions have been registered in the networks. In particular, for the cases of Net-3, PW06, Anytown, and Cimitile, the following values are chosen from the literature: kb = 0.071 × 1/h [48] and D = 45 μg/mg are included in the range identified by Bocelli et al. [49]. In the performed tests, sodium hypochlorite is supposed to be injected with a constant dosage equal to 0.68 mg/L, and the initial chlorine and TTHMs concentrations at any demand node are fixed equal to 0.68 mg/L and 10.1 μg/L, respectively. As regards AVWSS, the values of the rate coefficients for chlorine decay and TTHMs formation have been calibrated, which is explained in more detail in Section 3.4.
Considering that the Italian regulatory limit for TTHMs is of 30 μg/L, the TTHMs concentration threshold is fixed equal to 25 μg/L as a precautionary measure. A different threshold can be adopted in a different country or when using the methodology for studying another DPB, through modifying the identified vulnerable areas and number of users exposed. For the hydraulic and water quality simulations, the time step is 5 min. The vulnerability indexes are computed considering 24 h of simulation (TOT = 24 h).

3.1. Net-3 Network

The literature case study Net-3 is composed of 117 pipes, 92 junctions, two pumps, and three tanks, and it is supplied by two sources that are modeled as reservoirs (Figure 1a). Geometric data, base demands, and patterns are reported in the EPANET User’s Manual. The Hazen–Williams law is adopted as the resistance formula. Sodium hypochlorite is supposed to be injected in the two reservoirs.
Figure 2 presents the vulnerability contour plot maps obtained for the five considered indexes. The maps furnish complementary information about the different kind of exposure that they describe. Considering [TTHMs]imax and CiTTHMs (Figure 2a,b), the green regions represent the “safe” zones, where the [TTHMs]lim is never exceeded. The vulnerable nodes identified with CiTTHMs, from yellow to red areas, are less than the ones identified considering [TTHMs]imax, because the latter parameter is more restrictive and considers the node as critical even if the imposed limit has been exceeded only at one time step. In particular, 59 nodes reach during the day a TTHMs concentration above the limit, while CiTTHMs values above the threshold are registered in only 18 nodes. In particular, the nodes 131, 166, and 243 have the maximum CiTTHMs, which is approximately equal to 37 μg/L. Moreover, the highest values of both [TTHMs]imax and CiTTHMs are registered in node 243.
As expected, the maps obtained with Vic,n and TiE (Figure 2c,d) are very similar, indicating that there are eight nodes (the red ones) in which the limit concentration is exceeded for all of the simulation period (Vic,n = TiE = 1). All of these nodes are characterized by high water age values. In this group, the node 243 is also included.
Finally, the computed MiTTHMS values range from 277 μg to about 485 × 104 μg, as reported in Figure 2e. The vulnerability map is slightly different from the other ones, and node 15, which is not included in the ones selected through Vic,n and TiE, is individuated as the most vulnerable one. In fact, in this node, the water age is lower, but the base demand, representing the daily mean values of the outputs in the node, is higher, indicating a large number of users.

3.2. PW06 Network

The PW06 network (Figure 1b) has 47 pipes and 33 demand nodes with elevations that vary between 10–30 m, and it is supplied from a single source (reservoir), where sodium hypochlorite is supposed to be injected. The base demands are half of the ones reported in the original paper, and the demand pattern of Figure 3a has been assigned in all of the nodes. The Hazen–Williams law is adopted as the resistance formula.
Unlike the Net-3 system, in this case study, the [TTHMs]imax and CiTTHMs values are very similar, as shown from the heat maps reported in Figure 4a,b. For both parameters, the limit is exceeded in the same 13 nodes, and they identify node 15 as the most vulnerable one.
As in the Net-3 case, the heat maps of TiE and Vic,n, as reported in Figure 4c,d, are similar, with 10 nodes in which the limit concentration is exceeded for the entire simulation period (Vic,n = TiE = 1). They are located in the opposite side with respect to the source, where the water ages have larger values. Also, the mass parameters individuate the same vulnerable area (Figure 4e), and the maximum MiTTHMS value (MiTTHMS = 14,053 × 103 µg) is registered in node 7, which is included in the critical nodes selected through Vic,n and TiE.

3.3. Anytown Network

The Anytown water distribution system (Figure 1c) has 16 junctions, two tanks, 34 pipes, three pumps, and one reservoir, where sodium hypochlorite is supposed to be injected. The network has been partially modified for the purposes of the paper. In particular, in the nodes 20, 50, 60, 80, and 130, the base demands have been reduced by one third with respect to the original values, and the pattern in Figure 3b has been assigned, while for the rest of the nodes, the base demands have been reduced by one fifth with respect to the original values, and the demand pattern of Figure 3a has been used. The Hazen–Williams law is adopted as the resistance formula for the simulation.
Figure 5 depicts the heat maps obtained with the five considered parameters. In this case study, unlike PW06, [TTHMs]imax and CiTTHMs identify the same critical nodes, but the values that they assume are very different, as indicated from heat maps in Figure 5a,b. In 13 nodes, the TTHMs concentration exceeds the limit, while the CiTTHMs reaches the threshold in only six nodes. For both parameters, the most critical node is the number 170, which is located far from the source. Similarly to the Net-3 network, the comparison between Figure 5a and Figure 5b indicates that the maximum concentration is a more restrictive parameter.
The parameters TiE and Vic,n individuate the same nodes as critical, but their values are different, as indicated from the heat maps of Figure 5c,d. In particular, the values of the dimensionless exposure time are lower than the ones of the normalized contaminated volume. The limit concentration is exceeded for the entire simulation period (Vic,n = TiE = 1) in only nodes 140 and 170.
The heat map corresponding to the parameter MiTTHMS reported in Figure 5e is slightly different from the other ones, and the maximum value (MiTTHMS = 41,756 × 103 µg) is observed at node 80, which is not included in the ones individuated from the other parameters.

3.4. Aurunci Valcanneto Water Supply System (AVWSS)

The considered main trunk of the AVWSS is located in the Lazio region, and it serves about 110,000 inhabitants distributed among 25 towns. As shown in Figure 1d, the system is supplied by three sources (inlet nodes 1, 12, and 13), and water is delivered in nine outlet nodes, which represent the insertion into nine different WDSs.
In November 2008, a campaign of water quality measurements was realized, and the measured values of chlorine and TTHMs concentrations ([Cl2]m and [TTHMs]m) are reported in Table 1. A better description of the AVWSS can be found in Di Cristo et al. [36], where the system has been used as a case study to compare and validate some empirical and kinetics models for the prediction of the TTHMs concentration.
The demands in the delivering nodes are reported in Table 1, and no pattern is assigned, since the AVWSS is a main trunk. The simulation is performed in steady condition, and the obtained water age in the nodes are reported in Table 1. The limit considered on TTHMs concentration is 1.5 µg/L, which is fixed considering the extension of the WDSs connected to the outlet nodes. In other words, this threshold corresponds to about the 25 μg/L at the more distant outputs nodes of the connected WDSs. The initial chlorine concentration is 0.30 mg/L, while the TTHMs one is null. The reaction rate coefficients to model chlorine decay, kb (Equation (7)), and TTHMs formation, D (Equation (8)), are fixed through a specific calibration equal to 0.02 1/h and 43 µg/mL, respectively.
Without a pattern in the nodes, the parameter CiTTHMs does not make sense, and TiE and Vic,n are coincident. As shown in Table 1, the results of the vulnerability analysis identify different critical nodes according to the different indexes. In particular, [TTHMs]imax selects node 9 as the most vulnerable one (with a maximum value of 3.91 µg/L), which corresponds to the node with the maximum water age. TiE assume the value 0 in four nodes where the [TTHMs]lim has never been exceeded, and 1 in the nodes 5, 7, 8, 9, 10, corresponding to ages higher than 6 h. In terms of mass, the maximum value is delivered in node 10, with the maximum base demand (MiTTHMs = 16,460.44 × 103 µg), which is included in the group identified by TiE.

3.5. Cimitile Network

Figure 1e reports the scheme of the water distribution network of Cimitile, a small town located in the Campania region in Southern Italy. The system, which serves about 3700 residents located in 2.74 km2, consists of 433 pipes, with a total length of approximately 20 km, five regulation valves, and 404 junction nodes. Pipe diameters range from 20 mm to 260 mm. Predominant pipe materials are cast iron and steel, while there is only a small percentage of high-density polyethylene (HDPE) pipes. The water comes from a regional main trunk (not reported in the scheme) and it goes into the system through three inlet nodes (Gescal, Galluccio, and Pisacane).
Base demands values, as shown in the scheme of Figure 6 have been obtained working out utility consumption billing data, and the daily total flow rate supplied to the users is about 34.25 m3/s. The same water demand pattern, as characterized by two main peaks (Figure 3a), has been assigned to all of the nodes. The Darcy–Weisbach head loss formula is used for the simulation, and the hydraulic parameters have been previously calibrated using tank levels and pressure measurements. The chlorination of the water occurs in the three inlet nodes.
The vulnerability contour plot maps, as reported in Figure 7, indicate that also in this case, the vulnerable indexes tend to identify different risk areas considering different kinds of exposure. For a better understanding, Figure 6 reports the scheme of the network with the more vulnerable nodes selected using the different indexes. Moreover, Figure 8 illustrates the temporal variability of the TTHMs concentration in nodes 211, 290, and 328.
In terms of [TTHMs]imax and CiTTHMs (Figure 7a,b), the most critical nodes are located in the west part of the network. The number of nodes characterized by the excess of the attention limit in terms of maximum actual TTHMs concentration are more copious, because, as already mentioned, this parameter leads to also classifying as vulnerable those nodes where the [TTHMs]lim is exceeded only once during the simulation time. Both parameters individuate the node 211 as the most critical one (Figure 6), in which the maximum CiTTHMs is about 27 μg/L. Furthermore, the same node 211 is located in the area of the network that is characterized by the highest water ages.
The heat maps show that the exposure time and the normalized contaminated volume (Figure 7c,d) enable identifying the same critical areas in the system. In particular, the zones with high TiE and Vic,n values are located in the north, central, and south parts of the network, and the lower and higher residence times nodes are equally involved. In five nodes 7, 211, 212, 242, and 328 (Figure 6), the concentration exceeds the limit for the entire simulation period (Vic,n = TiE = 1).
The MiTTHMs values range from 600 × 103 μg at node 214 to about 23 × 1012 μg at node 290 (Figure 7e). The high values of this parameter coincide with nodes characterized by large base demands, namely with a high number of users served. Using this index, the most vulnerable is node number 290, which is located in northwest region of the network, and is not included in the group individuated by the other parameters (Figure 6).
Figure 8 reports the TTHMs concentration time variability during the simulating day for the nodes 211, 290, and 328, identified as the most vulnerable during the analysis.
In node 290, which is identified as the most critical one according to MiTTHMs, the TTHMs concentration exceeds the attention limit for about 8 h at the beginning of the simulation and in the last few minutes. In the other considered nodes, 211 and 328, which are vulnerable according to all of the parameters except MiTTHMs, the TTHMs concentration are always above the threshold fixed at 25 μg/L, reaching a maximum value of around 27 μg/L in node 211. These results indicate that in the area of node 290, many users were exposed to a high concentration for a short period, while in the areas of nodes 211 and 328, fewer consumers were always exposed to TTHMs concentrations above the fixed threshold. This kind of analysis is very useful for optimizing the network management for better water quality.

4. Conclusions

Water distribution vulnerability to elevated TTHMs concentrations is explored through the use of five different indexes that enable evaluating the population’s exposure. The methodology is applied to five case studies: three literature schemes—namely Net-3, PW06, and Anytown—and two real systems, namely AVWSS and Cimitile, which are both located in Italy.
The main outcome of the obtained results is that the introduced indexes identify different critical areas with respect to elevated TTHMs concentrations, depending on the kind of exposure to which they refer. They enable identifying the higher risk nodes with respect to the different kinds of exposure: high concentration exposure for a short period, or chronic exposure to low TTHMs concentrations. The comparison between the maps obtained from the different indexes suggests that their performances strongly depend on the behavior of the system.
The presented research is significant because vulnerability analyses are important for water companies and administrators to design public health and strategies to reduce risks. The proposed methodology provides valuable information for assisting water utilities in the WDSs management considering water quality. In particular, the results of the presented study and the information embedded are also effective in the implementation of methodologies for the formulation of optimization problems such as chlorine dosage or chlorine booster allocation in the system or for planning network modifications. The vulnerability indexes can also be adopted for performing the district sectorization and skeletonization of a network accounting for water quality. Moreover, even if it is presented considering the TTHMs, the proposed methodology can be used for evaluating the exposure to others DPBs, adopting the adequate modelling for predicting the concentration and selecting the relative threshold limit.

Author Contributions

Conceptualization, all authors; Methodology, C.Q. and A.L.; Software, C.Q.; Validation, A.L., Formal Analysis, C.Q. and C.D.C. Data Curation, C.Q. and A.L.; Writing-Original Draft Preparation, C.Q. and C.D.C.; Writing-Review & Editing, C.Q. and C.D.C.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vasconcelos, J.J.; Rossman, L.A.; Grayman, W.M.; Boulas, P.F.; Clark, R.M. Kinetics of chlorine decay. J. Am. Water Works Assoc. 1997, 89, 54–65. [Google Scholar] [CrossRef]
  2. Gang, D.C.; Clevenger, T.E.; Banerji, S.K. Modeling Chlorine Decay in Surface Water. J. Environ. Inform. 2003, 1, 21–27. [Google Scholar] [CrossRef] [Green Version]
  3. Greenberg, A.E. Public Health Aspects of Alternative Water Disinfectants. J. Am. Water Works Assoc. 1981, 73, 31–33. [Google Scholar] [CrossRef]
  4. White, G.C. Ozone. In Handbook of Chlorination Alternative Disinfectants; Van Nostrand Reinhold: New York, NY, USA, 1992; pp. 1046–1110. [Google Scholar]
  5. Morris, R.D.; Audet, A.M.; Angelillo, I.F. Chlorination, chlorination by-products and cancer: A metanalysis. Am. J. Public Health 1992, 2, 955–963. [Google Scholar] [CrossRef]
  6. Krewski, D.; Balbus, J.; Butler-Jones, D.; Haas, C.; Isaac-Renton, J.; Roberts, K.J.; Sinclair, M. Managing Health Risks from Drinking Water—A Report to the Walkerton Inquiry. J. Toxicol. Environ. Health Part A 2002, 65, 1635–1823. [Google Scholar] [CrossRef] [PubMed]
  7. Weinberg, H.S.; Krasner, S.W.; Richardson, S.W.; Thurston, A.D. The Occurrence of Disinfection By-Product (DBPs) of Health Concern in Drinking Water: Results of a Nationwide DBP Occurrence Study; EPA International Report 2002, EPA/600/R-02/068; EPA: Washington, DC, USA, 2002.
  8. Waller, K.; Swan, S.H.; DeLorenze, G.; Hopkins, B. Trihalomethanes in Drinking Water and Spontaneous Abortion. Epidemiology 1998, 9, 134–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Dodds, L.; King, W.D. Relation between Trihalomethane Compounds and Birth Defects. Occup. Environ. Med. 2001, 58, 443–446. [Google Scholar] [CrossRef] [PubMed]
  10. Richardson, S.D.; Simmons, J.E.; Rice, G. Disinfection Byproducts: The next Generation. Environ. Sci. Technol. 2002, 36, 198–205. [Google Scholar] [CrossRef]
  11. Chowdhury, S.; Champagne, P. Risk from exposure to trihalomethanes during shower: Probabilistic assessment and control. Sci. Total. Environ. 2009, 407, 1570–1578. [Google Scholar] [CrossRef] [PubMed]
  12. Villanueva, C.M.; Cantor, K.P.; Grimalt, J.O.; Castaño-Vynials Malats, N.; Silverman, D.; Tardon, A.; Garcia-Closas, R.; Serra, C.; Carrato, A.; Rothman, N.; et al. Bladder cancer and exposure to disinfection byproducts in water through ingestion, bathing, showering and swimming in pools: Findings from the Spanish Bladder Cancer Study. Am. J. Epidemiol. 2007, 165, 148–156. [Google Scholar] [CrossRef] [PubMed]
  13. Murray, R.; Janke, R.; Uber, J. The Threat Ensemble Vulnerability Assessment (TEVA) Program for Drinking Water Distribution System Security. In Proceedings of the World Water and Environmental Resources Congress, Salt Lake City, UT, USA, 27 June–1 July 2004. [Google Scholar]
  14. Gunnarsdottir, M.J.; Gardarsson, S.M.; Elliott, M.; Sigmundsdottir, G.; Bartram, J. Benefits of Water Safety Plans: Microbiology, Compliance, and Public Health. Environ. Sci. Technol. 2012, 46, 7782–7789. [Google Scholar] [CrossRef] [PubMed]
  15. Propato, M.; Piller, O.; Uber, J.G. A Sensor Location Model to Detect Contaminations in Water Distribution Networks. Impacts Glob. Clim. Chang. 2005, 1–12. [Google Scholar] [CrossRef]
  16. Janke, R.; Murray, R.; Uber, J.; Taxon, T. Comparison of Physical Sampling and Real-Time Monitoring Strategies for Designing a Contamination Warning System in a Drinking Water Distribution System. J. Water Resour. Plan. Manag. 2006, 132, 310–313. [Google Scholar] [CrossRef]
  17. Maslia, M.L.; Sautner, J.B.; Aral, M.M.; Reyes, J.J.; Abraham, J.E.; Williams, R.C. Using Water-Distribution System Modeling to Assist Epidemiologic Investigations. J. Water Resour. Plann. Manag. 2000, 126, 180–198. [Google Scholar] [CrossRef]
  18. Perelman, L.; Maslia, M.L.; Ostfeld, A.; Sautner, J.B. Using Aggregation/Skeletonization Network Models for Water Quality Simulations in Epidemiologic Studies. J. Am. Water Works Assoc. 2008, 100, 122–133. [Google Scholar] [CrossRef]
  19. Saldarriaga, J.G.; García, S.; Leön, N. A Methodology to Preserve Water Quality Modeling in Skeletonized Network Models. In Proceedings of the World Environmental and Water Resources Congress, Albuquerque, NM, USA, 20–24 May 2012; pp. 2932–2942. [Google Scholar]
  20. Weisel, C.P.; Kim, H.; Haltmeier, P.; Klotz, J.B. Exposure estimates to disinfection byproducts of chlorinated drinking water. Environ. Health Perspect. 1999, 107, 3–10. [Google Scholar] [CrossRef]
  21. Li, X.-F.; Mitch, W.A. Drinking water disinfection byproducts (DBPs) and human health effects: Multidisciplinary challenges and opportunities. Environ. Sci. Technol. 2018, 52, 1681–1689. [Google Scholar] [CrossRef] [PubMed]
  22. Krasner, S.W.; Cantor, K.P.; Weyer, P.J.; Hildesheim, M.; Amy, G. Case study approach to modeling historical disinfection by-product exposure in Iowa drinking waters. J. Environ. Sci. 2017, 58, 183–190. [Google Scholar] [CrossRef] [PubMed]
  23. World Health Organization. Guidelines for Drinking-Water Quality; World Health Organization: Geneva, Switzerland, 2004; Volume 1. [Google Scholar]
  24. Health Canada, Guidelines for Canadian Drinking Water Quality. 2007. Available online: http://www.hc-sc.gc.ca/ewhsemt/alt_formats/hecs-sesc/pdf/pubs/water-eau/2010-sum_guide-res_recom/sum_guide-res_recom-eng.pdf (accessed on 9 July 2018).
  25. National Health and Medical Research Council, Australian Drinking Water Guidelines; 2004; Volume 6. Available online: http://www.nhmrc.gov.au/guidelines/publications/eh34 (accessed on 9 July 2018).
  26. Propato, M.; Uber, J.G. Vulnerability of Water Distribution Systems to Pathogen Intrusion: How Effective Is a Disinfectant Residual? Environ. Sci. Technol. 2004, 38, 3713–3722. [Google Scholar] [CrossRef] [PubMed]
  27. Murray, R.; Uber, J.; Janke, R. Model for Estimating Acute Health Impacts from Consumption of Contaminated Drinking Water. J. Water Resour. Plan. Manag. 2006, 132, 293–299. [Google Scholar] [CrossRef]
  28. Nilsson, K.A.; Buchberger, S.G.; Clark, R.M. Simulating Exposures to Deliberate Intrusions into Water Distribution Systems. J. Water Resour. Plan. Manag. 2005, 131, 228–236. [Google Scholar] [CrossRef]
  29. Khanal, N.; Buchberger, S.G.; McKenna, S.A. Distribution system contamination events: Exposure, influence, and sensitivity. J. Water Resour. Plann. Manag. 2006, 132, 283–292. [Google Scholar] [CrossRef]
  30. Davis, M.J.; Janke, R. Importance of Exposure Model in Estimating Impacts When a Water Distribution System Is Contaminated. J. Water Resour. Plan. Manag. 2008, 134, 449–456. [Google Scholar] [CrossRef]
  31. Davis, M.J.; Janke, R. Development of a Probabilistic Timing Model for the Ingestion of Tap Water. J. Water Resour. Plan. Manag. 2009, 135, 397–405. [Google Scholar] [CrossRef]
  32. Thompson, S.L.; Casman, E.; Fischbeck, P.; Small, M.J.; VanBriesen, J.M. Vulnerability Assessment of a Drinking Water Distribution System: Implications for Public Water Utilities. In Proceedings of the World Environmental and Water Resources Congress, Tampa, FL, USA, 15–19 May 2007. [Google Scholar]
  33. Baoyu, Z.; Xinhua, Z.; Yuan, Z. Vulnerability Assessment of Regional Water Distribution Systems. In Proceedings of the 2009 International Conference on Environmental Science and Information Application Technology, Wuhan, China, 4–5 July 2009; pp. 473–477. [Google Scholar]
  34. Islam, N.; Sadiq, R.; Rodriguez, M.J.; Legay, C. Assessing Regulatory Violations of Disinfection By-Products in Water Distribution Networks Using a Non-Compliance Potential Index. Environ. Monit. Assess. 2016, 188, 304. [Google Scholar] [CrossRef] [PubMed]
  35. Islam, N.; Sadiq, R.; Rodriguez, M.J. Optimizing Locations for Chlorine Booster Stations in Small Water Distribution Networks. J. Water Resour. Plan. Manag. 2017, 143, 04017021. [Google Scholar] [CrossRef]
  36. Di Cristo, C.; Esposito, G.; Leopardi, A. Modelling trihalomethanes formation in water supply systems. Environ. Technol. 2013, 34, 61–70. [Google Scholar] [CrossRef] [PubMed]
  37. Kavanaugh, M.C.; Trussell, A.R.; Cromer, J.; Trussell, R.R. An Empirical Kinetic Model of Trihalomethane Formation: Applications to Meet the Proposed THM Standard. J. Am. Water Works Assoc. 1980, 72, 578–582. [Google Scholar] [CrossRef]
  38. Brown, D.; Bridgeman, J.; West, J.R. Understanding data requirement for trihalomethane formation in water supply systems. Urban Water J. 2011, 8, 41–56. [Google Scholar] [CrossRef]
  39. Adin, A.; Katzhendler, J.; Alkaslassy, D.; Rav-Acha, C. Trihalomethane Formation in Chlorinated Drinking Water: A Kinetic Model. Water Res. 1991, 25, 797–805. [Google Scholar] [CrossRef]
  40. Di Cristo, C.; Leopardi, A.; de Marinis, G. Assessing measurement uncertainty on trihalomethane prediction through kinetic models in water supply systems. J. Water Supply Res. Technol. AQUA 2015, 64, 516–528. [Google Scholar] [CrossRef]
  41. Di Cristo, C.; Leopardi, A.; Quintiliani, C.; de Marinis, G. Drinking Water Vulnerability Assessment after Disinfection through Chlorine. Procedia Eng. 2015, 119, 389–397. [Google Scholar] [CrossRef]
  42. Gang, D.C.; Segar, R.L.; Clevenger, T.E.; Banerji, S.K. Using Chlorine Demand to Predict TTHM and HAA9 Formation. J. Am. Water Works Assoc. 2002, 94, 76–86. [Google Scholar] [CrossRef]
  43. Rossman, L. EPANET User’s Manual; Risk Reduction Engineering Laboratory, U.S. Environmental Protection Agency: Cincinnati, OH, USA, 1994.
  44. Shang, F.; Uber, J.G.; Rossman, L.A. Modeling Reaction and Transport of Multiple Species in Water Distribution Systems. Environ. Sci. Technol. 2007, 42, 808–814. [Google Scholar] [CrossRef]
  45. Shang, F.; Uber, J.G.; Rossman, L.A. EPANET Multi-Species Extension User’s Manual; Risk Reduction Engineering Laboratory, U.S. Environmental Protection Agency: Cincinnati, OH, USA, 2008.
  46. Prasad, T.D.; Walters, G.A. Minimizing residence times by rerouting flows to improve water quality in distribution networks. Eng. Optim. 2006, 38, 923–939. [Google Scholar] [CrossRef]
  47. Walski, T.M.; Brill, E.D.; Gessler, J.; Goulter, I.C.; Jeppson, R.M.; Lansey, K.; Lee, H.; Liebman, J.C.; Mays, L.; Morgan, D.R.; et al. Battle of networks models: Epilogue. J. Water Resour. Plan. Manag. 1987, 113, 191–203. [Google Scholar] [CrossRef]
  48. Elshorbagy, W.A. Kinetics of THM Species in Finished Drinking Water. J. Water Resour. Plan. Manag. 2000, 126, 21–28. [Google Scholar] [CrossRef]
  49. Boccelli, D.; Tryby, M.E.; Uber, J.G.; Summers, R.S. A reactive species model for chlorine decay and THM formation under rechlorination conditions. Water Res. 2003, 37, 2654–2666. [Google Scholar] [CrossRef]
Figure 1. Scheme of case studies. (a) Net-3, (b) PW06, (c) Anytown, (d) Aurunci-Valcanneto Water Supply System (AVWSS), and (e) Cimitile.
Figure 1. Scheme of case studies. (a) Net-3, (b) PW06, (c) Anytown, (d) Aurunci-Valcanneto Water Supply System (AVWSS), and (e) Cimitile.
Water 10 00912 g001
Figure 2. Vulnerability contour plot maps for Net-3: (a) maximum actual total trihalomethanes (TTHMs) concentration, (b) demand weighted average TTHMs concentration, (c) dimensionless exposure time, (d) normalized contaminated volume, and (e) contaminated mass load.
Figure 2. Vulnerability contour plot maps for Net-3: (a) maximum actual total trihalomethanes (TTHMs) concentration, (b) demand weighted average TTHMs concentration, (c) dimensionless exposure time, (d) normalized contaminated volume, and (e) contaminated mass load.
Water 10 00912 g002
Figure 3. Demand patterns used in case studies. (a) PW06 and Cimitile, (a,b) Anytown.
Figure 3. Demand patterns used in case studies. (a) PW06 and Cimitile, (a,b) Anytown.
Water 10 00912 g003
Figure 4. Vulnerability contour plot maps for PW06: (a) maximum actual TTHMs concentration, (b) demand weighted average TTHMs concentration, (c) dimensionless exposure time, (d) normalized contaminated volume, and (e) contaminated mass load.
Figure 4. Vulnerability contour plot maps for PW06: (a) maximum actual TTHMs concentration, (b) demand weighted average TTHMs concentration, (c) dimensionless exposure time, (d) normalized contaminated volume, and (e) contaminated mass load.
Water 10 00912 g004aWater 10 00912 g004b
Figure 5. Vulnerability contour plot maps for Anytown: (a) maximum actual TTHMs concentration, (b) demand weighted average TTHMs concentration, (c) dimensionless exposure time, (d) normalized contaminated volume, and (e) contaminated mass load.
Figure 5. Vulnerability contour plot maps for Anytown: (a) maximum actual TTHMs concentration, (b) demand weighted average TTHMs concentration, (c) dimensionless exposure time, (d) normalized contaminated volume, and (e) contaminated mass load.
Water 10 00912 g005aWater 10 00912 g005b
Figure 6. Scheme of the Cimitile network with the more critical nodes.
Figure 6. Scheme of the Cimitile network with the more critical nodes.
Water 10 00912 g006
Figure 7. Vulnerability contour plot maps for Cimitile: (a) maximum actual TTHMs concentration, (b) demand weighted average TTHMs concentration, (c) dimensionless exposure time, (d) normalized contaminated volume, and (e), contaminated mass load.
Figure 7. Vulnerability contour plot maps for Cimitile: (a) maximum actual TTHMs concentration, (b) demand weighted average TTHMs concentration, (c) dimensionless exposure time, (d) normalized contaminated volume, and (e), contaminated mass load.
Water 10 00912 g007
Figure 8. Temporal variation of the TTHMs concentration in nodes 211, 290, and 328.
Figure 8. Temporal variation of the TTHMs concentration in nodes 211, 290, and 328.
Water 10 00912 g008
Table 1. Water quality measurements of November 2008 and the results of the vulnerability analysis.
Table 1. Water quality measurements of November 2008 and the results of the vulnerability analysis.
Node NumberAge
(h)
Q
(L/s)
[Cl2]m
(mg/L)
[TTHMs]m
(µg/L)
[TTHMs]imax
(µg/L)
TiEMiTTHMs
(103 µg)
20.1720.271.970.0400
31.8980.271.490.4600
43.922--0.9400
56.8520.261.321.65126.28
65.988--1.4500
76.372.5-2.001.5419.33
88.573-2.452.031138.36
918.06470.243.023.9119823.64
108.21120.25-3.20116,460.44

Share and Cite

MDPI and ACS Style

Quintiliani, C.; Di Cristo, C.; Leopardi, A. Vulnerability Assessment to Trihalomethane Exposure in Water Distribution Systems. Water 2018, 10, 912. https://doi.org/10.3390/w10070912

AMA Style

Quintiliani C, Di Cristo C, Leopardi A. Vulnerability Assessment to Trihalomethane Exposure in Water Distribution Systems. Water. 2018; 10(7):912. https://doi.org/10.3390/w10070912

Chicago/Turabian Style

Quintiliani, Claudia, Cristiana Di Cristo, and Angelo Leopardi. 2018. "Vulnerability Assessment to Trihalomethane Exposure in Water Distribution Systems" Water 10, no. 7: 912. https://doi.org/10.3390/w10070912

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop