A Drinking Water Relevant Water Chemistry Model for the Free Chlorine and Cyanuric Acid System from 5°C to 35°C

Abstract

In the United States, approved methods to measure free chlorine concentrations in drinking water systems adding sodium dichloroisocyanurate (dichlor) or trichloroisocyanuric acid (trichlor) as chlorine sources exhibit measurement bias from the presence of chlorinated cyanurates, leading to overestimated free chlorine concentrations for regulatory compliance. One option to overcome this limitation is to estimate free chlorine concentrations using an established water chemistry model (full model), but the full model has only been determined for 25°C. The current research used a simplified version of the full model (simple model) and estimated the unknown temperature dependence (5°C to 35°C) of the two remaining equilibrium constants (K_7a and K_9a) required for the simple model. At 0 M ionic strength (μ), \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\ln { { \rm { K } } _ { 7 { \rm { a } } } } = - { \frac { 4 , 671 } { { { \rm { T } } _ { \rm { K } } } } } + 4.95$$ \end{document} or \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$ { \rm { p } } { { \rm { K } } _ { 7 { \rm { a } } } } = { \frac { 2 , 028 } { { { \rm { T } } _ { \rm { K } } } } } - 2.15$$ \end{document} , \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\Delta { \rm{H}}_{7{ \rm{a}}}^0$$ \end{document} = 38.8 ± 6.0 kJ/mol (95% confidence interval [CI]), \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\ln { { \rm { K } } _ { 9 { \rm { a } } } } = - { \frac { 5 , 133 } { { { \rm { T } } _ { \rm { K } } } } } + 3.79$$ \end{document} or \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$ { \rm { p } } { { \rm { K } } _ { 9 { \rm { a } } } } = { \frac { 2 , 229 } { { { \rm { T } } _ { \rm { K } } } } } - 1.65$$ \end{document} , and \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\Delta { \rm{H}}_{9{ \rm{a}}}^0$$ \end{document} = 42.7 ± 3.0 kJ/mol (95% CI). At 25°C and μ of 0 M, the simple model estimated pK_7a and pK_9a are 4.65 ± 0.059 (95% CI) and 5.83 ± 0.020 (95% CI), respectively. As an example of the impact of temperature, the free chlorine concentration for a 2 mg Cl₂/L dichlor addition (pH 7.0) decreases from 0.90 mg Cl₂/L at 25°C to 0.60 mg Cl₂/L at 5°C. If temperature was not considered, a system operating at 5°C would overestimate their free chlorine concentration by 50%, which could have significant implications for understanding disinfection efficacy, illustrating the developed model's significance.

Introduction

In July 2001, the use of dichlor (anhydrous sodium dichloroisocyanurate, C₃Cl₂N₃NaO₃, or sodium dichloroisocyanurate dihydrate, C₃Cl₂N₃NaO₃•2H₂O) and trichlor (trichloroisocyanuric acid, C₃Cl₃N₃O₃) for the routine treatment of drinking water in the United States was first approved under the Federal Insecticide, Fungicide, and Rodenticide Act (FIFRA) (Kuechler, 2009). Subsequently, several manufacturers have achieved NSF 60 certification for various dichlor or trichlor products for addition to drinking water (NSF International, 2018). Ultimately, in the United States, dichlor or trichlor may be approved by Safe Drinking Water Act (SDWA) primacy agencies for use in drinking water. In addition to use in U.S. drinking water systems, the World Health Organization (WHO) provides guidelines on the use of dichlor in drinking water (World Health Organization, 2007, 2011), and dichlor has been suggested as a form of chlorine for household point-of-use water treatment (Clasen and Edmondson, 2006). Kuechler (2009) summarizes the stated advantages for using dichlor or trichlor as free chlorine sources in drinking water, including that they are high percentage chlorine chemicals (55–92%) and solid chlorine forms that (1) require simpler handling than liquid or gas chlorine forms, (2) have storage stabilities on the order of years, and (3) do not add calcium to the water.

Addition of dichlor or trichlor to drinking water is equivalent to adding free chlorine (hypochlorous acid [HOCl] plus hypochlorite ion [OCl⁻]) and cyanuric acid (C₃H₃N₃O₃) in a fixed proportion. The C₃N₃O₃ moiety (126.05 molecular weight) is typically denoted as Cy; therefore, cyanuric acid is represented as H₃Cy, Cl₃Cy is trichlor, and NaCl₂Cy is dichlor. For context and assuming 100% chemical purity, each mg Cl₂/L added from dichlor or trichlor adds 0.89 or 0.59 mg Cy/L respectively. Regardless of how cyanuric acid becomes present in water with free chlorine (i.e., through addition of dichlor, trichlor, or free chlorine and cyanuric acid individually), a rapid equilibrium (Matte et al., 1990) is obtained between H₃Cy and its three dissociation products, HOCl, and chlorinated cyanurates (Fig. 1) that has been established at 25°C by O'Brien (1972) and O'Brien et al. (1974), referred to herein as the full model (Table 1). The resulting free chlorine concentration present in the water (assuming no ammonia is present) is dependent on (1) temperature, (2) pH, (3) chemical added (i.e., dichlor or trichlor), (4) chemical dose that determines the total chlorine (sum of chlorine in free chlorine and chlorinated cyanurates, TOTCl) and the total cyanurate added (sum of the cyanurate in H₃Cy and its three dissociation products and chlorinated cyanurates, TOTCy), and (5) extent of chlorine demand that has occurred, which determines the remaining TOTCl concentration (Wahman, 2018a).

FIG. 1.

Hypochlorous acid, cyanuric acid, cyanuric acid ionization products, and chlorinated cyanurate equilibria for the full model. Gray shaded regions indicate chemicals and reactions that are relevant under drinking water conditions and chemical dosages and are those included in the simple model. Hydrogen ion for dissociation reactions and water and hypochlorous acid for hydrolysis reactions are not shown for clarity, and required full reactions for the full and simple models are shown in Table 1. Adapted with permission from Brady et al. (1963) Copyright 1963 American Chemical Society. Dotted boxes denote those chemicals which are called out by name in the figure by the text next to them (e.g., trichlor, dichlor, or cyanuric acid).

Table 1.

Summary of Required Equilibria to Simulate Chlorinated Cyanurate Water Chemistry for the Full and Simple Models and Associated Equilibrium Constants as Noted

Stoichiometry	Constant	pK ^a	Equilibrium included
Stoichiometry	Constant	pK ^a	Full model	Simple model
Cl₃Cy + H₂O ⇋ HCl₂Cy + HOCl	K_1a	1.8	✓
HCl₂Cy ⇋ Cl₂Cy⁻ + H⁺	K₂	3.75	✓
H₂ClCy ⇋ HClCy⁻ + H⁺	K₄	5.33	✓
H₃Cy ⇋ H₂Cy⁻ + H⁺	K₆	6.88	✓	✓
Cl₂Cy⁻ + H₂O ⇋ HClCy⁻ + HOCl	K_7a	4.51	✓	✓
HClCy⁻ ⇋ ClCy²⁻ + H⁺	K₈	10.12	✓
HClCy⁻ + H₂O ⇋ H₂Cy⁻ + HOCl	K_9a	5.62	^*	✓
H₂Cy⁻ ⇋ HCy²⁻ + H⁺	K₁₀	11.40	✓
ClCy²⁻ + H₂O ⇋ HCy²⁻ + HOCl	K_11a	6.90	✓
HCy²⁻ ⇋ Cy³⁻ + H⁺	K₁₂	13.5	✓
HOCl ⇋ OCl⁻ + H⁺	K	7.54^b	✓	✓

Determined by O'Brien et al. (1974) but not required in the full model implementation.

O'Brien et al. (1974) at 0.020 M ionic strength and 25°C unless otherwise noted.

Morris (1966) at 0 M ionic strength and 25°C.

Because of the combination of a fast equilibrium (Matte et al., 1990) between H₃Cy and its three dissociation products, HOCl, and chlorinated cyanurates (Fig. 1 and Table 1) and the possibility of direct reactions with method reagents (Jensen and Johnson, 1990), typical wet chemistry methods (e.g., N,N–diethyl–p–phenylenediamine [DPD]-based methods and amperometric titration) used to measure the free chlorine concentration in drinking water systems adding dichlor or trichlor exhibit measurement bias from the presence of chlorinated cyanurates (Whittle, 1970; Canelli, 1980; Wajon and Morris, 1980; Jensen and Johnson, 1990; Wahman et al., 2017). Therefore, when TOTCy is greater than zero and no ammonia is present to form chloramines, the free chlorine concentration is overestimated, which is important because the effective disinfectant in systems wherein chlorinated cyanurates and free chlorine are present is considered to be free chlorine, specifically HOCl (Canelli, 1974; Engel et al., 1983; Yamashita et al., 1988; Wojtowicz, 1996; Murphy et al., 2015).

Currently, the inability to accurately determine free chlorine concentrations in systems adding dichlor or trichlor prohibits determination of disinfection efficacy, including compliance with U.S. primary disinfection requirements or numeric free chlorine distribution system residual (i.e., greater than detectable) requirements required by several states (Wahman and Pressman, 2015). The current issue with free chlorine measurement is highlighted in the State of Utah's plan review policy for dichlor and trichlor installations (Utah Division of Drinking Water, 2017) and a recent memorandum from the U.S. Environmental Protection Agency (USEPA, 2018).

To overcome the inability to directly measure free chlorine in systems using dichlor or trichlor, one option is to simulate free chlorine concentrations using the established equilibrium water chemistry of the free chlorine and cyanuric acid system. For example, a web-based application has recently been developed to estimate free chlorine concentrations using this concept (Wahman, 2017). A limitation with the water chemistry estimate is that it is currently only applicable to 25°C, limiting its usefulness in actual drinking water systems where temperature may vary over a wide range (e.g., 5°C to 35°C). In addition, O'Brien (1972) and O'Brien et al. (1974) used chemical concentrations one to two orders of magnitude greater than would be seen in drinking water (e.g., 149 mg Cl₂/L TOTCl and 378 mg Cy/L TOTCy) to estimate the equilibrium constants evaluated herein. Therefore, extending the water chemistry model to include temperatures from 5°C to 35°C and evaluating it under relevant drinking water conditions are an immediate need.

If pH, TOTCy, and TOTCl are known, the full model for the free chlorine and cyanuric acid system consists of 10 independent equilibrium expressions that are temperature dependent, 2 mass balance equations for TOTCl and TOTCy, and 12 unknown chemical concentrations (Table 1 and Fig. 1) (O'Brien et al., 1974). As described in Wahman (2018a) and for drinking water samples, pH and TOTCl can be directly measured, and TOTCy can be conservatively calculated based on the dichlor or trichlor mass originally dosed. Fortunately, under relevant drinking water conditions (i.e., 6.5 ≤ pH ≤9.5, TOTCl ≤10 mg Cl₂/L, and TOTCy ≤10 mg Cy/L), the system is reduced to six chemicals present at meaningful concentrations (Fig. 1) because of the pH range and TOTCl and TOTCy concentrations. Therefore, only four (K, K₆, K_7a, and K_9a) temperature-dependent equilibrium expressions (Table 1 and Fig. 1) are required along with the TOTCl and TOTCy mass balances. The resulting model is herein termed the simple model (Table 1).

Figure 2 displays the impact of using the simple model versus the full model to simulate free chlorine concentrations at 25°C and was generated using relevant drinking water combinations (13,671 total simulations) of pH (6.5 to 9.5 every 0.1 units), TOTCl (0 to 10 mg Cl₂/L every 0.5 mg Cl₂/L), and TOTCy (0 to 10 mg Cy/L every 0.5 mg Cy/L). For most conditions, the simulated free chlorine concentrations by the simple model were slightly greater than those from the full model, and the differences were generally <0.05 mg Cl₂/L, reaching a maximum of 0.20 mg Cl₂/L at pH 9.5 (Fig. 2A). In the worst-case scenario (pH 9.5), the deviation of the simple model was only 6% of the full model-simulated free chlorine concentration (Fig. 2B). Furthermore, Supplementary Fig. S1 shows that for the conditions evaluated in the current research (Table 2), the full and simple models simulate essentially the same free chlorine concentrations, assuming 25°C.

FIG. 2.

Comparison of the difference in simulated free chlorine concentrations at 25°C between the full model of O'Brien et al. (1974) and the proposed simple model over a range of relevant drinking water conditions (6.5 ≤ pH ≤9.5, TOTCl ≤10 mg Cl₂/L, and TOTCy ≤10 mg Cy/L) on a concentration (A) and percentage (B) basis. Δ free chlorine (mg Cl₂/L) = simple model − full model and Δ free chlorine (%) = (simple model − full model)/full model. Equilibrium constants used for all simulations are those shown in Table 1.

Table 2.

Summary of Experimental Conditions (208 Total Samples)

Value ± standard deviation (No. of samples with absorbance >0.020 with a 10-cm cell)
Temperature (°C)	pH	Total chlorine (mg Cl₂/L)	Dosed cyanuric acid (mg/L H₃Cy)
4.3 ± 0.13 (18)	7.50 ± 0.003 (54)	1.87 ± 0.054 (90)	1 (60)
5.3 ± 0.28 (18)	8.30 ± 0.004 (74)	3.70 ± 0.10 (118)	2 (55)
9.0 ± 0.19 (19)	9.00 ± 0.003 (80)		3 (48)
15.3 ± 0.24 (21)			4 (45)
19.8 ± 0.076 (21)
20.6 ± 0.10 (22)
24.4 ± 0.075 (22)
30.2 ± 0.14 (22)
34.2 ± 0.14 (22)
34.5 ± 0.26 (23)

In the simple model, two of the four required equilibrium expressions (Table 1 and Fig. 1) represent the acid–base chemistry of free chlorine (HOCl and OCl⁻, K) and cyanuric acid (H₃Cy and H₂Cy⁻, K₆) whose temperature dependence in Kelvin (T_K) have been established [Eqs. (1) and (2), respectively, 0 M ionic strength (μ)] (Morris, 1966; Wahman, 2018b). \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} { \rm { pK } } = { \frac { 3 , 000 } { { { \rm { T } } _ { \rm { K } } } } } - 10.0686 + 0.0253 { { \rm { T } } _ { \rm { K } } } \tag { 1 } \end{align*} \end{document} \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} { \rm { p } } { { \rm { K } } _6 } = { \frac { 1 , 743 } { { { \rm { T } } _ { \rm { K } } } } } + 1.12 \tag { 2 } \end{align*} \end{document}

Using the reported K and K₆ relationships, only the temperature dependence of two equilibrium constants (K_7a and K_9a) involving chlorinated cyanurates is required to expand the simple model to a range of relevant drinking water temperatures. Therefore, the goal of the current research was to determine the temperature dependence of K_7a and K_9a from 5°C to 35°C. With the temperature dependence of K_7a and K_9a determined, the simple model's temperature dependence is complete, providing the necessary information to simulate free chlorine concentrations in drinking water systems using dichlor or trichlor over a range of practical temperatures (e.g., updating the existing web-based application to incorporate temperature).

In addition to the previously described drinking water application, the simple model developed in this research may have broader applications. For instance, cyanuric acid, dichlor, and trichlor have been used as disinfectants and algicides in cooling towers (Hammond et al., 1986) and added to outdoor swimming pools since 1958 to stabilize chlorine residual by minimizing sunlight-related degradation (Canelli, 1974).

Materials and Methods

Reagent and glassware preparation

Solutions were prepared in ultra-pure water (18.2 MΩ-cm, Barnstead NANOpure Diamond). Stock cyanuric acid solutions (1,000 mg/L of H₃Cy, pH >9, 129.07 molecular weight) were prepared from reagent grade cyanuric acid (Acros Organics), accounting for reagent purity of 98%. Free chlorine and OCl⁻ solutions were prepared from reagent grade sodium hypochlorite (Fisher Chemical). Carbonate buffers were made from sodium bicarbonate (Fisher Chemical) and adjusted to the target pH with 0.05 and 0.5 N hydrochloric acid (HCl) or sodium hydroxide (NaOH) as required. All glassware was treated to make it chlorine-demand free before use (Summers et al., 1996).

Analytical methods

Spectrophotometric measurements (292 nm) were conducted with a 10-cm path length quartz cuvette using a Hach DR/4000U Spectrophotometer. The spectrophotometer was blanked with ultra-pure water before use. Triplicate absorbance measurements were taken and averaged for any measured sample.

The sample pH and temperature were measured immediately before each spectrophotometric measurement with an Orion pH electrode with automatic temperature compensation connected to an Orion Star A211 pH meter that had been calibrated with nominal pH standards of 4, 7, and 10 per manufacturer's instructions. Free and total chlorine colorimetric DPD samples were analyzed in parallel with each measurement using a Hach SL1000 Portable Parallel Analyzer (Hach Company, 2015, 2017).

Absorbance verification experiments

Morris (1966) estimated molar absorptivity values for HOCl and OCl⁻ from 230 to 300 nm between 5°C and 35°C, demonstrating that a linear relationship between absorbance and HOCl and OCl⁻ concentrations existed at each temperature and wavelength studied, and the relationship did not change substantially with temperature. In this research, absorbance was measured at 292 nm because (1) at the TOTCl and TOTCy concentrations used in the experiments, HOCl and OCl⁻ are the only two chemicals present in the experimental water that meaningfully absorb at 292 nm (Morris, 1966; O'Brien, 1972; O'Brien et al., 1974) (see Supplementary Data for confirmation, Absorbance Interference Check at 292 nm section), (2) the maximum molar absorptivity for OCl⁻ at 292 nm is ∼3.5 times that of HOCl at 230 nm (Morris, 1966), and (3) at the pH values evaluated in this research, OCl⁻ is present at approximately the same or greater concentration than HOCl. Therefore, 292 nm represented the wavelength providing the maximum absorbance that could be attributed to free chlorine and not subject to interference from other chemical species present.

As subsequently discussed (see Equilibrium Constant Experiments section), the molar absorptivity values for HOCl and OCl⁻ at 292 nm are not required for data analysis. Rather, only verification on the current spectroscopy equipment that a linear relationship exists between concentration and absorbance at 292 nm from 5°C to 35°C is required. Therefore, 0, 0.25, 0.50, 1, 2, and 4 mg Cl₂/L OCl⁻ (pH 10) standards were prepared in duplicate at 5, 20, and 35°C and measured at 292 nm to verify that the current equipment could produce results comparable with Morris (1966). OCl⁻ was chosen to demonstrate spectroscopy equipment linearity at the various temperatures because it was simpler to create solutions of OCl⁻ versus HOCl, and OCl⁻ has a greater molar absorptivity, allowing TOTCl concentrations (0 to 4 mg Cl₂/L) to be used that were relevant to this research.

Equilibrium constant determination experiments

Equilibrium constant experiments

All experiments and associated measurements were conducted in a temperature controlled chamber maintained nominally at either 5, 10, 15, 20, 25, 30, or 35°C. Duplicate experiments were conducted at 5, 20, and 35°C. To acquire the required data to allow estimation of K_7a and K_9a temperature dependence, solutions of free chlorine only and free chlorine with cyanuric acid additions of either 1, 2, 3, or 4 mg/L cyanuric acid were created in 0.004 M carbonate buffer at pHs 7.5, 8.3, and 9.0 at each experimental temperature and for two TOTCl concentrations (1.9 and 3.7 mg Cl₂/L).

A sample's free chlorine concentration was determined using a combination of spectrophotometry and colorimetric methods. In a solution containing only free chlorine (i.e., no cyanuric acid addition and designated as Solution #1) at a fixed pH and temperature, the 292-nm absorbance is proportional to HOCl and OCl⁻ concentrations [Eq. (3)]. \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} Ab{s_1} & = { \varepsilon _{HOCl}}l{\; \propto _0}\;TOTFre{e_1} + { \varepsilon _{OCl}}l{\; \propto _1}\;TOTFre{e_1} \\ & = { \varepsilon _{HOCl}}l{\; \propto _0}\;TOTCl + { \varepsilon _{OCl}}l{\; \propto _1}\;TOTCl \tag{3} \end{align*} \end{document}

In Equation (3), Abs₁ is the measured total absorbance at 292 nm in Solution #1; α₀ and α₁ are the fractions of free chlorine that are HOCl and OCl⁻, respectively; ɛ_HOCl and ɛ_OCl are the molar absorptivities of HOCl and OCl⁻ at 292 nm, respectively; l is the path length of the absorbance cell; and TOTFree₁ is the free chlorine concentration in Solution #1 that equals the TOTCl concentration in Solution #1.

Adding some mass of cyanuric acid to Solution #1 creates Solution #2 wherein the absorbance at 292 nm is now related by Equation (4). Note that the only difference between Solution #1 and Solution #2 is the presence of cyanuric acid in Solution #2 (i.e., Solution #2 has the same pH, temperature, and TOTCl concentration as Solution #1). Therefore, α₀ and α₁ are also the same between Solution #1 and Solution #2. In Equation (4), Abs₂ is the measured total absorbance at 292 nm in Solution #2 and TOTFree₂ is the free chlorine concentration in Solution #2. Because some TOTCl is now bound as chlorinated cyanurates in Solution #2, TOTFree₂ is less than TOTCl. \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} Ab{s_2} = { \varepsilon _{HOCl}}l{\; \propto _0}\;TOTFre{e_2} + { \varepsilon _{OCl}}l{\; \propto _1}\;TOTFre{e_2} \tag{4} \end{align*} \end{document}

Taking the absorbance ratio for Solution #1 and Solution #2 results in Equation (5), which relates the absorbance ratio to the fraction of free chlorine in Solution #2 and becomes independent of HOCl and OCl⁻ molar absorptivities. \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} {{Ab{s_2}} \over {Ab{s_1}}} & = {{{{ \varepsilon _{HOCl}}l{\; \propto _0}\;TOTFre{e_2} + { \varepsilon _{OCl}}}l{ \propto _1}\;TOTFre{e_2}} \over {{ \varepsilon _{HOCl}}l{\; \propto _0}\;TOTCl + { \varepsilon _{OCl}}}l{ \propto _1}\;TOTCl} \\ & = \;{{TOTFre{e_2}} \over {TOTCl}} \tag{5} \end{align*} \end{document}

TOTFree₂ can be calculated using Equation (5) and the measured (by an appropriate method such as colorimetric DPD) free or total chlorine concentration in Solution #1 as shown by Equation (6). \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} TOTFre { e_2 } = { \frac { Ab { s_2 } } { Ab { s_1 } } } \left( { TOTCl } \right) = { \frac { Ab { s_2 } } { Ab { s_1 } } } \left( { TOTFre { e_1 } } \right) \tag { 6 } \end{align*} \end{document}

Equation (6) does not require knowledge of the HOCl and OCl⁻ molar absorptivities. Rather, the only requirement for Equation (6) is that at each temperature, HOCl and OCl⁻ have a constant molar absorptivity (i.e., that concentration is proportional to absorbance) as has been determined by Morris (1966).

Simple model

Equation (7) represents a third-order polynomial equation derived from the simple model in terms of total free chlorine (TOTFree) concentration (see Derivation of Equation (7) from the Simple Model section in Supplementary Data). In Equation (7), α_o is the concentration fraction of TOTFree that is HOCl [Supplementary Eq. (S1)]; \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\left\{ {{{ \rm{H}}^ + }} \right\} $$ \end{document} is hydrogen ion activity that equals 10^−pH; K_6c is the concentration-based cyanuric acid first equilibrium constant calculated from K₆ [Supplementary Eq. (S4)]; and K_7a and K_9a are chlorinated cyanurate equilibrium constants. Equation (7) was used along with experimentally determined TOTFree concentrations [Eq. (6)] to estimate the temperature dependence of K_7a and K_9a as subsequently discussed (see Equilibrium Constant Temperature Dependence section and Parameter Estimation section). \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} \\ & \left( { { \frac { { { \rm { K } } _ { 6 { \rm { c } } } } { \propto _0 } ^2 } { { { \rm { K } } _ { 7 { \rm { a } } } } { { \rm { K } } _ { 9 { \rm { a } } } } \left\{ { { { \rm { H } } ^ + } } \right\} } } } \right) { \rm { TOTFre } } { { \rm { e } } ^3 } + \left( { 2 { \rm { TOTCy } } { \frac { { { \rm { K } } _ { 6 { \rm { c } } } } { \propto _0 } ^2 } { { { \rm { K } } _ { 7 { \rm { a } } } } { { \rm { K } } _ { 9 { \rm { a } } } } \left\{ { { { \rm { H } } ^ + } } \right\} } } } \right.\\ & \quad\left. + { \frac { { { \rm { K } } _ { 6 { \rm { c } } } } { \propto _0 } } { { { \rm { K } } _ { 9 { \rm { a } } } } \left\{ { { { \rm { H } } ^ + } } \right\} } } - { \rm { TOTCl } } { \frac { { { \rm { K } } _ { 6 { \rm { c } } } } { \propto _0 } ^2 } { { { \rm { K } } _ { 7 { \rm { a } } } } { { \rm { K } } _ { 9 { \rm { a } } } } \left\{ { { { \rm { H } } ^ + } } \right\} } } \right.) { \rm { TOTFre } } { { \rm { e } } ^2 } + \\ & \left ( {{ \rm{TOTCy}}{{{{ \rm{K}}_{6{ \rm{c}}}}{ \propto _0}} \over {{{ \rm{K}}_{9{ \rm{a}}}} \left\{ {{{ \rm{H}}^ + }} \right\} }} + 1 + {{{{ \rm{K}}_{6{ \rm{c}}}}} \over { \left\{ {{{ \rm{H}}^ + }} \right\} }} - { \rm{TOTCl}}{{{{ \rm{K}}_{6{ \rm{c}}}}{ \propto _0}} \over {{{ \rm{K}}_{9{ \rm{a}}}} \left\{ {{{ \rm{H}}^ + }} \right\} }}} \right) \\ & \quad{ \rm{TOTFree}} - { \rm{TOTCl}} \left( {1 + {{{{ \rm{K}}_{6{ \rm{c}}}}} \over { \left\{ {{{ \rm{H}}^ + }} \right\} }}} \right) = 0 \tag{7} \end{align*} \end{document}

Ionic strength and activity corrections for K₆ and K

To determine the concentration-based equilibrium constants K_c and K_6c from K and K₆ [Supplementary Eqs. (S1) and (S4)], the range of μ for the experimental conditions was calculated using Equation (8) based on known chemical additions, the measured pH and temperature, and a charge balance to account for the acid or base addition required to obtain the measured pH (Wahman, 2018b). For the μ calculations, relationships for carbonate and water ionization constants with temperature were taken from Goldberg et al. (2002) and Stumm and Morgan (1996), respectively. \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} { \rm { \mu } } = \frac { 1 } { 2 } \mathop \sum \limits_ { { \rm { all \;ions } } } { { \rm { c } } _ { \rm { i } } } { { \rm { z } } _ { \rm { i } } } ^2 \tag { 8 } \end{align*} \end{document}

Subsequently, μ was used to calculate the single charge ion (z = 1) activity coefficient (γ₁) from the Davies equation using Equation (9) (Benjamin, 2002). In Equation (9), the term A_DH is determined from Equation (10) at the specified T_K (Benjamin, 2002), and in Equation (10), the dielectric constant for water (ɛ_water) is determined at the specified temperature in °C (T_c) from Equation (11) (Malmberg and Maryott, 1956). \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} \log { { \rm { \gamma } } _ { \rm { z } } } = { \rm { \; } } - { { \rm { A } } _ { { \rm { DH } } } } { { \rm { z } } ^2 } \left( { { \frac { \sqrt { \rm { \mu } } } { 1 + \sqrt { \rm { \mu } } } } - 0.2 { \rm { \mu } } } \right) \tag { 9 } \end{align*} \end{document} \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} {{ \rm{A}}_{{ \rm{DH}}}} = 1.82{ \times {10^6}{ \left( {{{ \rm{ \varepsilon }}_{{ \rm{water}}}}{{ \rm{T}}_{ \rm{K}}}} \right) ^{ - 3 / 2}}} \tag{10} \end{align*} \end{document} \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} {{ \rm{ \varepsilon }}_{{ \rm{water}}}} & = 87.740 - 0.40008{{ \rm{T}}_{ \rm{c}}} + 9.398{ \times }{10^{ - 4}}{{ \rm{T}}_{ \rm{c}}}^2 \\ & \quad- 1.410{\times }{10^{ - 6}}{{ \rm{T}}_{ \rm{c}}}^3 \tag{11} \end{align*} \end{document}

The determined range of μ (0.0040 to 0.0045 M) and corresponding γ₁ (0.928 to 0.935) values covered a narrow range. Because γ₁ was identical to two significant figures over its entire range, γ₁ was fixed at 0.93 for all experimental conditions to simplify analysis, representing at worst a 0.5% deviance from the actual γ₁.

Equilibrium constant temperature dependence

The basic relationship between an equilibrium constant (K_eq) and T_K is shown by Equation (12) (Stumm and Morgan, 1996). In Equation (12), ΔH⁰ is the standard enthalpy change of the reaction, whereas R is the universal gas constant (8.314 × 10⁻³ kJ/[mol·K]. If ΔH⁰ is assumed to be independent of T_K, Equation (13) results. Equation (13) can be further simplified to Equation (14) or Equation (15). Using Equation (14), specific relationships for K_7a and K_9a with T_K are represented by Equations (16) and (17). Equations (1), (2), (7), (16), and (17) were used for parameter estimation, and A_7a, A_9a, B_7a, and B_9a represent constants estimated during parameter estimation with \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$ { { \rm { A } } _ { 7 { \rm { a } } } } = { \frac { \Delta { { \rm { H } } ^0 } _ { 7 { \rm { a } } } } { \rm { R } } } $$ \end{document} and \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$ { { \rm { A } } _ { 9 { \rm { a } } } } = { \frac { \Delta { { \rm { H } } ^0 } _ { 9 { \rm { a } } } } { \rm { R } } } $$ \end{document} . \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} { \frac { { \rm { d } } \ln { { \rm { K } } _ { { \rm { eq } } } } } { { \rm { d } } { { \rm { T } } _ { \rm { K } } } } } = { \frac { \Delta { { \rm { H } } ^0 } } { { \rm { R } } { { \rm { T } } _ { \rm { K } } } ^2 } } \tag { 12 } \end{align*} \end{document} \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} \ln { \frac { { { \rm { K } } _ { { \rm { eq } } , { \rm { \; } } { { \rm { T } } _ { { \rm { K } } 2 } } } } } { { { \rm { K } } _ { { \rm { eq } } , { { \rm { T } } _ { { \rm { K } } 1 } } } } } } { \rm { \; } } = \mathop \smallint \limits_ { { { \rm { T } } _ { { \rm { K } } 1 } } } ^ { { { \rm { T } } _ { { \rm { K } } 2 } } } { \frac { \Delta { { \rm { H } } ^0 } } { { \rm { R } } { { \rm { T } } _ { \rm { K } } } ^2 } } { \rm { d } } { { \rm { T } } _ { \rm { K } } } = { \frac { \Delta { { \rm { H } } ^0 } } { \rm { R } } } \left( { \frac { 1 } { { { { \rm { T } } _ { { \rm { K } } 1 } } } } - \frac { 1 } { { { { \rm { T } } _ { { \rm { K } } 2 } } } } } \right) \tag { 13 } \end{align*} \end{document} \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} \ln { { \rm { K } } _ { { \rm { eq } } , { { \rm { T } } _ { \rm { K } } } } } = - { \frac { \Delta { { \rm { H } } ^0 } } { { \rm { R } } { { \rm { T } } _ { \rm { K } } } } } + { \rm { constant } } = - { \frac { \rm { A } } { { { \rm { T } } _ { \rm { K } } } } } + { \rm { B } } \tag { 14 } \end{align*} \end{document} \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} { \rm { p } } { { \rm { K } } _ { { \rm { eq } } , { { \rm { T } } _ { \rm { K } } } } } = { \frac { \rm { A } } { \left( { \ln 10 } \right) \left( { { { \rm { T } } _ { \rm { K } } } } \right) } } - { \frac { \rm { B } } { \ln 10 } } \tag { 15 } \end{align*} \end{document} \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} \ln { { \rm { K } } _ { 7 { \rm { a } } , { { \rm { T } } _ { \rm { K } } } } } = - { \frac { { { \rm { A } } _ { 7 { \rm { a } } } } } { { { \rm { T } } _ { \rm { K } } } } } + { { \rm { B } } _ { 7 { \rm { a } } } } \tag { 16 } \end{align*} \end{document} \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} \begin{align*} \ln { { \rm { K } } _ { 9 { \rm { a } } , { { \rm { T } } _ { \rm { K } } } } } = - { \frac { { { \rm { A } } _ { 9 { \rm { a } } } } } { { { \rm { T } } _ { \rm { K } } } } } + { { \rm { B } } _ { 9 { \rm { a } } } } \tag { 17 } \end{align*} \end{document}

Parameter estimation

Estimation of A_7a, A_9a, B_7a, and B_9a and corresponding determinations of \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $${{ \rm{K}}_{7{ \rm{a}} , {{ \rm{T}}_{ \rm{K}}}}}$$ \end{document} , \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $${{ \rm{K}}_{9{ \rm{a}} , {{ \rm{T}}_{ \rm{K}}}}}$$ \end{document} , \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\Delta {{ \rm{H}}^0}_{7{ \rm{a}}}$$ \end{document} , and \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\Delta {{ \rm{H}}^0}_{9{ \rm{a}}}$$ \end{document} were accomplished in R, a freeware language and environment for statistical computing and graphics, and R is available for download at www.r-project.org (R Core Team, 2017). R's capabilities were extended through the use of the polynom package, a free user-contributed package downloaded through the R software, to numerically solve Equation (7) for the TOTFree concentration (Venables et al., 2016). In the R software, the function optim() was used to vary estimates of A_7a, A_9a, B_7a, and B_9a by minimizing the residual sum of squares between all (n = 208) the experimentally measured TOTFree concentrations and the TOTFree concentrations required to solve Equation (7). To provide an estimate of uncertainty in the parameter estimates, bootstrapping with sample replacement (n = 10,000) was used.

Results and Discussion

Absorbance verification experiments

Figure 3 shows results demonstrating the verification of a linear relationship (R² = 0.999) between concentration and absorbance on the current equipment at 292 nm over the range of temperatures evaluated in this research, using OCl⁻ as the demonstration chemical. For each temperature, similar slopes (345 to 356 M⁻¹ cm⁻¹ molar absorptivity) were determined, verifying that the current equipment can reproduce similar results of Morris (1966). Based on variability in the lowest standard (0.25 mg Cl₂/L), the detection limit corresponded to an absorbance of 0.002 with a 10-cm cell (APHA et al., 2005). To provide a measure of safety, the minimum absorbance values used for any sample used in parameter estimation was set to 0.020 with a 10-cm cell.

FIG. 3.

Linearity verification between concentration and absorbance at 292 nm and 5°C (A), 20°C (B), or 35°C (C). Note that the ordinate has been normalized to a 1-cm path length and 1 mg Cl₂/L = 1.41 × 10⁻⁵ M.

Equilibrium constant determination experiments

Source data

Experimental data obtained from the equilibrium constant determination experiments are summarized in Supplementary Table S1, consisting of 208 data points. A summary of the experimental conditions is provided in Table 2. For each temperature, 24 possible samples containing both free chlorine and cyanuric acid were created for analysis. The lower number of samples shown in Table 2 at each temperature is a result of samples whose absorbance were <0.020 and removed from subsequent analyses. As temperature decreased, a greater number of samples did not achieve the 0.020 absorbance threshold, which is a result of lower concentrations of free chlorine being present as temperature decreased.

Full model (without temperature dependence) representation of experimental data

To evaluate the need to incorporate temperature into the water chemistry model, a comparison of the experimentally determined free chlorine concentrations (Supplementary Table S1) with the simulated free chlorine concentrations using the full model (Table 1) was conducted (Fig. 4A). If there was no impact with temperature, the experimentally measured and model-simulated concentrations would fall on the 1:1 line displayed in Fig. 4A. It is clear in Fig. 4A that as temperatures decreased below 25°C, the full model (without temperature dependence) overestimated the free chlorine concentrations. To understand the discrepancy between the measured free chlorine concentrations and the full model-simulated free chlorine concentration, Supplementary Fig. S2 displays the difference on a concentration (Supplementary Fig. S2A) and percentage basis (Supplementary Fig. S2B). Differences of up to 1.2 mg Cl₂/L are evident with percentages reaching nearly 60% of the simulated value. Taken together, these results indicated that temperature must be considered to properly represent the experimental data.

FIG. 4.

Comparison of model-simulated versus experimentally determined free chlorine concentrations for the full model (A) and simple model (B). Full model simulations conducted using equilibrium constants shown in Table 1. Simple model simulations conducted using equilibrium constants determined in the current research.

Simple model estimated parameters

Using the data presented in Supplementary Table S1, A_7a, A_9a, B_7a, and B_9a were estimated using Equations (7), (16), and (17) to determine K_7a and K_9a at various temperatures. No apparent pattern is evident in the simple model residuals (experimentally measured minus model-simulated free chlorine concentrations) associated with the parameter estimation (Supplementary Fig. S3). A comparison of the simple model-simulated and experimentally measured free chlorine concentrations is summarized in Fig. 4B. If the simple model-simulated and experimentally measured free chlorine concentrations were identical, the points shown in Fig. 4B would form a line with a slope of one passing through the origin. The linear regression results in a line with a slope that is not significantly different than 1 (0.99 ± 0.010, 95% confidence interval [CI]), an intercept not significantly different from 0 (0.018 ± 0.021, 95% CI), and excellent R² (0.994), indicating that the simple model is representing the experimental data extremely well and that the assumption of constant ΔH⁰ [Eq. (13)] was justified.

A summary of the estimates of A_7a, A_9a, B_7a, and B_9a and corresponding pK_7a and pK_9a values ranging from 5 to 35°C is provided in Table 3 along with comparisons with the 25°C O'Brien et al. (1974) estimates adjusted to μ = 0 M (refer to Supplementary Data, Adjustment of O'Brien et al. (1974) Equilibrium Constants). From the estimated values of A_7a and A_9a, \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\Delta { \rm{H}}_{7{ \rm{a}}}^0$$ \end{document} = 38.8 ± 6.0 kJ/mol (95% CI) and \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\Delta { \rm{H}}_{9{ \rm{a}}}^0$$ \end{document} = 42.7 ± 3.0 kJ/mol (95% CI). Supplementary Fig. S4 shows histograms of the bootstrap results for A_7a, A_9a, B_7a, and B_9a that were used to determine their reported uncertainty in Table 3, indicating normally distributed estimates. In addition, Fig. 5 shows the paired estimates of A_7a and B_7a (Fig. 5A) and A_9a and B_9a (Fig. 5B) from the bootstrap results, detailing the correlation between each estimated pair and corresponding joint uncertainty that is expected based on the linear relationship [Eqs. (16) and (17)]. The correlation would need to be incorporated if a future uncertainty analysis on the simple model is conducted to determine safety factors so that free chlorine concentration is not overestimated by the simple model to ensure conservative estimates when estimating disinfectant efficacy in dichlor or trichlor systems.

FIG. 5.

Paired bootstrap estimates (n = 10,000) of parameters for K_7a (A) and K_9a (B).

Table 3.

Summary of Estimated Equilibrium Constants K_7a and K_9a and Associated Parameters from 5 to 35°C Determined in the Current Research Unless Otherwise Noted

Temperature (°C)	Value ±95% CI (μ = 0 M)
Temperature (°C)	pK_7a	pK_9a
5	5.14 ± 0.061	6.37 ± 0.041
10	5.01 ± 0.050	6.23 ± 0.033
15	4.89 ± 0.047	6.09 ± 0.026
20	4.77 ± 0.051	5.96 ± 0.021
25	4.65 ± 0.059	5.83 ± 0.020
25 ^a	4.51 ± 0.058^a	5.68 ± 0.035^a
30	4.54 ± 0.071	5.71 ± 0.022
35	4.43 ± 0.084	5.59 ± 0.027

Parameter (units)	Value ±95% CI (μ = 0 M)
Parameter (units)	K_7a	K_9a
A (K)	4,671 ± 725	5,133 ± 363
B	4.95 ± 2.52	3.79 ± 1.22
ΔH⁰ (kJ/mol)	38.8 ± 6.0	42.7 ± 3.0

O'Brien et al. (1974) after adjusting to 0 M ionic strength and converting the reported standard deviations to 95% confidence intervals using the experimental data from O'Brien (1972).

CI, confidence interval.

Using the reported standard deviations from O'Brien et al. (1974) and after adjusting their estimates determined at μ = 0.020 M to μ = 0 M with the experimental data provided in O'Brien (1972) (pK_7a = 4.45–4.57 and pK_9a = 5.64–5.71; 95% CI), the current 25°C simple model estimates (pK_7a = 4.59–4.71 and pK_9a = 5.81–5.85; 95% CI) are comparable yet significantly different than the full model estimates of O'Brien et al. (1974). The differences between the current estimates of pK_7a and pK_9a and those of O'Brien et al. (1974) may be a result of the different experimental conditions, methods used to estimate equilibrium constants, or models implemented (i.e., simple versus full). This research conducted experiments under conditions representative of drinking water (pH 7.5 to 9.0, <4 mg Cl₂/L TOTCl, and <4 mg Cy/L TOTCy), used a simplified model that only included chemical species important under drinking water conditions, and simultaneously estimated two equilibrium constants (K_7a and K_9a) using all the data collected from 5 to 35°C.

To estimate K_9a, O'Brien (1972) and O'Brien et al. (1974) used pHs from 9.4 to 10.6, chemical concentrations one to two orders of magnitude greater than the current research (149 mg Cl₂/L TOTCl and 378 mg Cy/L TOTCy), and relied on two estimated equilibrium constants from the full model (Table 1; K₈ and K₁₀) to estimate K_9a. In addition, O'Brien (1972) discussed an apparent linear trend with the experimental pH for individual K_9a estimates used to determine the reported K_9a and associated uncertainty, noting that further experimental work should be conducted to investigate the pH trend as it should not exist. The O'Brien (1972) linear trend with pH (R² = 0.96) is illustrated in Supplementary Fig. S5, where the O'Brien (1972) individual pK_9a estimates increased as the experimental pH under which the pK_9a estimates were determined decreased. If the O'Brien (1972) pK_9a estimates (5.75 and 5.77) at the lowest experimental pH (9.45) are more representative of the actual pK_9a value, then the current research estimated pK_9a (5.81–5.85) corresponds better with these previous results. For K_7a, O'Brien (1972) and O'Brien et al. (1974) did not conduct specific experiments to estimate K_7a. Rather, K_7a was calculated from three estimated equilibrium constants (Fig. 1 and Table 1; K₂, K_3a, and K₄) that were determined with TOTCl and TOTCy concentrations like those used to estimate K_9a. Regardless of the differences between the current estimates and those of O'Brien (1972) and O'Brien et al. (1974), the simple model presented herein represented the entire data set extremely well.

Practical implications

With the addition of K_7a and K_9a's temperature dependence to the known temperature dependence of K and K₆, a relevant drinking water chemistry model for the free chlorine and cyanuric acid system (i.e., systems using dichlor or trichlor) has been determined over a practical range of drinking water temperatures (5 to 35°C), allowing estimation of free chlorine concentrations. Figure 6 summarizes the temperature impact on the percentage free chlorine for a 2 mg Cl₂/L dichlor addition (1.8 mg Cy/L TOTCy) at 5, 15, 25, and 35°C. The percentage free chlorine decreases with decreasing temperature at all pH values simulated, highlighting the importance of accounting for temperature so as not to overestimate free chlorine concentration, which is the true disinfectant. For example and at pH 7.0, the percentage free chlorine decreases from 45% (0.90 mg Cl₂/L free chlorine) at 25°C to 31% (0.60 mg Cl₂/L free chlorine) at 5°C. If temperature was not considered, a system operating at 5°C would overestimate their free chlorine concentration by 50%, which could have significant implications for understanding disinfection efficacy (e.g., compliance with U.S. primary disinfection requirements).

FIG. 6.

Simple model simulation of percentage free chlorine concentration at various temperatures for a 2 mg Cl₂/L dichlor addition (i.e., 2 mg Cl₂/L total chlorine). Simulations conducted using equilibrium constants determined in the current research. Percentage free chlorine = 100 × actual free chlorine/total chlorine.

Summary

A simplified version of the free chlorine and cyanuric acid equilibrium water chemistry model (simple model) was used to estimate the unknown temperature dependence (5 to 35°C) of the two remaining equilibrium constants (K_7a and K_9a) required for the simple model. At 0 M ionic strength (μ), the estimated equilibrium constants were \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\ln { { \rm { K } } _ { 7 { \rm { a } } } } = - { \frac { 4 , 671 } { { { \rm { T } } _ { \rm { K } } } } } + 4.95$$ \end{document} or \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$ { \rm { p } } { { \rm { K } } _ { 7 { \rm { a } } } } = { \frac { 2 , 028 } { { { \rm { T } } _ { \rm { K } } } } } - 2.15$$ \end{document} and \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\ln { { \rm { K } } _ { 9 { \rm { a } } } } = - { \frac { 5 , 133 } { { { \rm { T } } _ { \rm { K } } } } } + 3.79$$ \end{document} or \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$ { \rm { p } } { { \rm { K } } _ { 9 { \rm { a } } } } = { \frac { 2 , 229 } { { { \rm { T } } _ { \rm { K } } } } } - 1.65$$ \end{document} . The corresponding estimated standard enthalpy change of the reactions was \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\Delta { \rm{H}}_{7{ \rm{a}}}^0$$ \end{document} = 38.8 ± 6.0 kJ/mol and \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland, xspace}\usepackage{amsmath, amsxtra}\usepackage{upgreek}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document} $$\Delta { \rm{H}}_{9{ \rm{a}}}^0$$ \end{document} = 42.7 ± 3.0 kJ/mol (95% CI). With the incorporation of temperature dependence, the completed simple model is now capable of simulating free chlorine concentrations to allow assessment of disinfection efficacy, including U.S. primary disinfection or numeric free chlorine distribution system residual (i.e., greater than detectable) requirements.

To apply the simple model in drinking water practice, one would follow the methodology outlined by Wahman (2018a) with the addition of temperature as another required measurement. Therefore, for a water sample, one needs to know the pH, TOTCl, TOTCy, and temperature as inputs to the simple model. As discussed in Wahman (2018a), pH, TOTCl, and now temperature could be measured directly on the water sample whereas the TOTCy concentration could be calculated from the total mass of dichlor or trichlor that was dosed into the water system, accounting for chemical purity. Because chlorine demand may have occurred, TOTCy cannot be calculated from the measured TOTCl in the water sample and dichlor or trichlor's chemical formula. Doing so would lead to an underestimation of TOTCy and an overestimation of the actual free chlorine concentration. With pH, TOTCl, TOTCy, and temperature known, the simple model provides an estimate of the sample's actual free chlorine concentration.

Footnotes

Acknowledgments

The authors acknowledge the technical assistance provided by Meridith Bartley of Neptune and Company for the implementation of the bootstrapping analysis in R and Richard Falk for fruitful discussions on relevant chemistry.

Disclaimer

This work has been subjected to the United States Environmental Protection Agency's (Agency's) review and has been approved for publication. The views expressed in this article are those of the authors and do not necessarily represent the views or policies of the Agency. Any mention of trade names, products, or services does not imply an endorsement by the Agency. The Agency does not endorse any commercial products, services, or enterprises. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Author Disclosure Statement

No competing financial interests exist.

References

APHA, AWWA, and WEF. (2005). Standard Methods for the Examination of Water and Wastewater. Washington, DC: American Public Health Association: American Water Works Association: Water Environment Federation.

Brady

A.P.

, Sancier

K.M.

, and Sirine

(1963). Equilibria in Solutions of Cyanuric Acid and its Chlorinated Derivatives. J. Am. Chem. Soc. 85, 3101.

Benjamin

M.M.

(2002). Water Chemistry. New York, NY: McGraw–Hill.

Canelli

(1974). Chemical, bacteriological, and toxicological properties of cyanuric acid and chlorinated isocyanurates as applied to swimming pool disinfection—A review. Am. J. Public Health, 64, 155.

Canelli

(1980). Evaluation of a portable bare–electrode amperometric analyzer for determining free chlorine in potable and swimming–pool waters. Water Res. 14, 1533.

Clasen

, and Edmondson

(2006). Sodium dichloroisocyanurate (NaDCC) tablets as an alternative to sodium hypochlorite for the routine treatment of drinking water at the household level. Int. J. Hyg. Environ. Health, 209, 173.

Engel

J.P.

, Rubin

A.J.

, and Sproul

O.J.

(1983). Inactivation of Naegleria gruberi cysts by chlorinated cyanurates. Appl. Environ. Microbiol. 46, 1157.

Goldberg

R.N.

, Kishore

, and Lennen

R.M.

(2002). Thermodynamic quantities for the ionization reactions of buffers. J. Phys. Chem. Ref. Data, 31, 231.

Hach Company. (2015). Method 10260, Free Chlorine. USEPA DPD Method, Chemkey Reagents, 0.12 to 4.6 mg/L as Cl₂. Loveland, CO: Hach Company.

10.

Hach Company. (2017). Method 10260, Total Chlorine. USEPA DPD Method, Chemkey Reagents, 0.04 to 10.0 mg/L as Cl₂. Hach Company, Loveland, Colorado.

11.

Hammond

B.G.

, Barbee

S.J.

, Inoue

, Ishida

, Levinskas

G.J.

, Stevens

M.W.

, Wheeler

A.G.

, and Cascieri

(1986). A review of toxicology studies on cyanurate and its chlorinated derivatives. Environ. Health Perspect. 69, 287.

12.

Jensen

J.N.

, and Johnson

J.D.

(1990). Interferences by monochloramine and organic chloramines in free available chlorine methods. 2. N,N–Diethyl–p–phenylenediamine. Environ. Sci. Technol., 24, 985.

13.

Kuechler

T.C.

(2009). Use of chlorinated isocyanurates for drinking water chlorination. Proc. Water Environ. Fed. 2009, 799.

14.

Malmberg

C.G.

, and Maryott

A.A.

(1956). Dielectric constant of water from 0° to 100°C. J. Res. Natl. Bur. Stand. 56, 1.

15.

Matte

, Solastiouk

, Merlin

, and Deglise

(1990). Étude cinétique de la N–chloration de l'acide cyanurique en phase aqueuse. Can. J. Chem. 68, 307.

16.

Morris

J.C.

(1966). The acid ionization constant of HOCl from 5 to 35°C. J. Phys. Chem. 70, 3798.

17.

Murphy

J.L.

, Arrowood

M.J.

, Lu

, Hlavsa

M.C.

, Beach

M.J.

, and Hill

V.R.

(2015). Effect of cyanuric acid on the inactivation of Cryptosporidium parvum under hyperchlorination conditions. Environ. Sci. Technol. 49, 7348.

18.

NSF International. Search for NSF Certified Drinking Water Treatment Chemicals. Available at: http://info.nsf.org/Certified/PwsChemicals/ (accessed January 11, 2018).

19.

O'Brien

J.E.

, 1972. Hydrolytic and Ionization Equilibria of Chlorinated Isocyanurate in Water. Dissertation, Harvard University, Cambridge, Massachusetts.

20.

O'Brien

J.E.

, Morris

J.C.

, and Butler

J.N.

(1974). Equilibria in aqueous solutions of chlorinated isocyanurate. In Rubin

A.J.

, Ed. Chemistry of Water Supply, Treatment, and Distribution. Ann Arbor, MI: Ann Arbor Science Publishers, Inc., pp. 333–358.

21.

R Core Team. (2017). R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.

22.

Stumm

, and Morgan

J.J.

(1996). Aquatic Chemistry Chemical Equilibria and Rates in Natural Waters. Hoboken, NJ: Wiley.

23.

Summers

R.S.

, Hooper

S.M.

, Shukairy

H.M.

, Solarik

, and Owen

(1996). Assessing the DBP yield: Uniform formation conditions. J. Am. Water Works Assoc. 88, 80.

24.

USEPA. Memorandom: Compliance Determination for Inactivation Requirements of the National Primary Drinking Water Regulations when a Public Water Systems Uses Dichlor and Trichlor for Primary Disinfection (April 18, 2018). Available at: https://asdwa.files.wordpress.com/2018/04/dichlor-trichlor-memo-2018-04-18.pdf (accessed August 27, 2018).

25.

Utah Division of Drinking Water. Plan Review Policy for Dichlor and Trichlor Installations (DDW–Eng–0032 April 27, 2017). Available at: https://documents.deq.utah.gov/drinking-water/engineering/new-facilities-new-public-water-systems/DDW-2017-001930.pdf (accessed August 27, 2018).

26.

Venables

, Hornik

, and Maechler

(2016). polynom: A Collection of Functions to Implement a Class for Univariate Polynomial Manipulations. R package version 1.3–9. https://CRAN.R-project.org/package=polynom

27.

Wahman

D.G.

(2017). Free Chlorine and Cyanuric Acid Simulator Application Description–Version 0.50, EPA/600/S–17/165. Washington, DC: United States Environmental Protection Agency.

28.

Wahman

D.G.

(2018a). Chlorinated cyanurates: Review of water chemistry and associated drinking water implications. J. Am. Water Works Assoc. 110, E1.

29.

Wahman

D.G.

(2018b). First acid ionization constant of the drinking water relevant chemical cyanuric acid from 5 to 35°C. Environ. Sci. Water Res. Technol. 4, 1522.

30.

Wahman

D.G.

, Alexander

, and Dugan

A.G.

(2017). Chlorinated cyanurates: Method interferences and application implications. In 2017 Water Quality Technology Conference. Portland, Oregon: American Water Works Association.

31.

Wahman

D.G.

, and Pressman

J.G.

(2015). Distribution system residuals—Is “detectable” still acceptable for chloramines? J. Am. Water Works Assoc. 107, 53.

32.

Wajon

J.E.

, and Morris

J.C.

(1980). The analysis of free chlorine in the presence of nitrogenous organic compounds. Environ. Int. 3, 41.

33.

Whittle

G.P.

(1970). Recent advances in determining free chlorine. In National Specialty Conference on Disinfection. Amherst, MA. American Society of Civil Engineers, New York, NY.

34.

Wojtowicz

J.A.

(1996). Relative bactericidal effectiveness of hypochlorous acid and chloroisocyanurates. J. Swimming Pool Spa Indus. 2, 34.

35.

World Health Organization. (2007). Sodium Dichloroisocyanurate in Drinking-Water. Background Document for Development of WHO Guidelines of Drinking-Water Quality. Geneva, Switzerland: World Health Organization.

36.

World Health Organization. (2011). Guidelines for Drinking-Water Quality. Geneva, Switzerland: World Health Organization.

37.

Yamashita

, Sakae

, Ishihara

, Isomura

, and Inoue

(1988). Virucidal effect of chlorinated water containing cyanuric acid. Epidemiol. Infect. 101, 631.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

1.55 MB

A Drinking Water Relevant Water Chemistry Model for the Free Chlorine and Cyanuric Acid System from 5°C to 35°C

Abstract

Abstract

Introduction

Materials and Methods

Reagent and glassware preparation

Analytical methods

Absorbance verification experiments

Equilibrium constant determination experiments

Equilibrium constant experiments

Simple model

Ionic strength and activity corrections for K6 and K

Equilibrium constant temperature dependence

Parameter estimation

Results and Discussion

Absorbance verification experiments

Equilibrium constant determination experiments

Source data

Full model (without temperature dependence) representation of experimental data

Simple model estimated parameters

Practical implications

Summary

Footnotes

Acknowledgments

Disclaimer

Author Disclosure Statement

References

Supplementary Material

Ionic strength and activity corrections for K₆ and K