1 Introduction 2 pH Calculation - US EPA
1 Introduction
The pH, alkalinity and total inorganic carbon algorithms that are incorporated into WASP come directly from QUAL2K & QUAL2Kw (Chapra, Pelletier and Tao, 2008). Every effort was made to insure that the implementation of the pH calculation within WASP is consistent with QUAL2K/.
The following sections were mostly taken from the QUAL2K documentation and where appropriate the equations were modified to match the assumptions and kinetic implementations in WASP.
2 pH Calculation
The following equilibrium, mass balance and electroneutrality equations define a freshwater dominated by inorganic carbon (Stumm and Morgan 1996),
K 1
[HCO3 ][H ]
[H
2
CO
* 3
]
(1)
K 2
[CO
2 3
][H
]
[HCO
3
]
(2)
Kw [H ][OH ]
(3)
c T
[H
2CO*3
]
[HCO3
]
[CO
2 3
]
(4)
Alk [HCO3 ] 2[CO32 ] [OH ] [H ]
(5)
where K1, K2 and Kw are acidity constants, Alk = alkalinity [eq L1], H2CO3* = the sum of dissolved carbon dioxide and carbonic acid, HCO3 = bicarbonate ion, CO32 = carbonate ion, H+ = hydronium ion, OH = hydroxyl ion, and cT = total inorganic carbon concentration [mole L1]. The brackets [ ] designate molar concentrations.
Note that the alkalinity is expressed in units of eq/L for the internal calculations. For input and output, it is expressed as mgCaCO3/L. The two units are related by
Alk(mgCaCO3/L) 50,000 Alk(eq/L)
(6)
The equilibrium constants are corrected for temperature by Harned and Hamer (1933):
pK w
=
4787.3 Ta
7.1321log10 (Ta ) 0.010365Ta
22.80
(7)
Plummer and Busenberg (1982):
logK1 = 356.3094 0.06091964Ta 21834.37 / Ta 126.8339logTa 1,684,915/ Ta2
(8)
Plummer and Busenberg (1982):
logK2 = 107.8871 0.03252849Ta 5151.79/ Ta 38.92561logTa 563,713.9 / Ta2
(9)
The nonlinear system of five simultaneous equations (1 through 5) can be solved numerically for the five unknowns: [H2CO3*], [HCO3], [CO32], [OH], and {H+}. An efficient solution method can be derived by combining Eqs. (1), (2) and (4) to define the quantities (Stumm and Morgan 1996)
0
[H ]2
[H ]2 K1[H ] + K1K2
(10)
1
[H
]2
K1[H ] K1[H ]
+
K1 K 2
(11)
2
[H ]2
K1 K 2 K1[H ] + K1K2
(12)
where 0, 1, and 2 = the fraction of total inorganic carbon in carbon dioxide, bicarbonate, and carbonate, respectively. Equations (3), (11), and (12) can be substituted into Eq. (4) to yield,
Alk = (1
2 2 )cT
Kw [H ] [H ]
(13)
Thus, solving for pH reduces to determining the root, {H+}, of
f ([H ]) = (1
2 2 )cT
Kw [H ]
[H ] Alk
(14)
where pH is then calculated with
pH log10[H ]
(15)
The root of Eq. (14) is determined with a numerical method. The user can choose bisection, Newton-Raphson or Brent's method (Chapra and Canale 2006, Chapra 2007) as specified on the QUAL2K sheet. The Newton-Raphson is the fastest but can sometimes diverge. In contrast, the bisection method is slower, but more reliable. Because it balances speed with reliability, Brent's method is the default.
2.1 Total Inorganic Carbon (cT)
Total inorganic carbon concentration increases due to fast carbon oxidation and plant respiration. It is lost via plant photosynthesis. Depending on whether the water is undersaturated or oversaturated with CO2, it is gained or lost via reaeration,
S cT
rcco FastCOxid rcca PhytoResp rcca
Bot AlgResp H
(16)
Bot AlgP hoto
rcca PhytoPhoto rcca
H
CO2Reaer
where
CO2Reaer kac(T )[CO2 ]s 0cT
(17)
where kac(T) = the temperature-dependent carbon dioxide reaeration coefficient [/d], and [CO2]s = the saturation concentration of carbon dioxide [mole/L].
The stoichiometric coefficients are computed as1
rcca
rca
gC mgA
moleC 12 gC
m3 1000L
(18)
rcco
1 roc
gC gO 2
moleC m3 12 gC 1000L
(19)
2.2 Carbon Dioxide Saturation
The CO2 saturation is computed with Henry's law,
[CO2 ]s K H pCO2
(20)
where KH = Henry's constant [mole (L atm)1] and pCO2 = the partial pressure of carbon dioxide in the atmosphere [atm]. Note that the partial pressure is input as a constant or as an
1 The conversion, m3 = 1000 L is included because all mass balances express volume in m3, whereas total inorganic carbon is expressed as mole/L.
environmental time function in units of ppm. The program internally converts ppm to atm using the conversion: 106 atm/ppm.
The value of KH can be computed as a function of temperature by (Edmond and Gieskes 1970)
pK H
=
2385.73 Ta
0.0152642Ta
14.0184
(21)
The partial pressure of CO2 in the atmosphere has been increasing, largely due to the combustion of fossil fuels (Error! Reference source not found.). Values in 2007 are approximately 103.416 atm (= 383.7 ppm).
Figure 1 Concentration of carbon dioxide in the atmosphere as recorded at Mauna Loa Observatory, Hawaii.
The CO2 reaeration coefficient can be computed from the oxygen reaeration rate by
kac (20)
32 0.25 44
0.923 k a
(20)
(22)
2.3 Effect of Control Structures: CO2
As was the case for dissolved oxygen, carbon dioxide gas transfer in streams can be influenced by the presence of control structures. WASP assumes that carbon dioxide behaves similarly to dissolved oxygen. Thus, the inorganic carbon mass balance for the element immediately downstream of the structure is written as
dcT,i dt
Qi1 Vi
c'T
,i1
Qi Vi
cT ,i
Qab,i Vi
cT ,i
Ei' Vi
cT ,i1 cT ,i
WcT,i Vi
ScT,i
(23)
where c'T,i1 = the concentration of inorganic carbon entering the element [mgO2/L], where
c'T ,i1
( 1
2
)cT ,i1
CO2,s,i1
CO2,s,i1 rd
2cT ,i1
(24)
where rd is calculated from the dam reaeration algorithm for dissolved oxygen.
2.4 Alkalinity (Alk)
As summarized in the present model accounts for changes in alkalinity due to several mechanisms:
Table 1 Processes that effect alkalinity.
Process
Utilize Create
Alkalinity change
Nitrif
NH4
NO3
Decrease
Denitr
NO3
Increase
OPHydr
SRP
Decrease
ONHydr
NH4
Increase
PhytoPhoto NH4
Decrease
NO3
Increase
SRP
Increase
PhytoResp
NH4
Increase
SRP
Decrease
PhytoUpN
NH4
Decrease
NO3
Increase
PhytoUpP
SRP
Increase
PhytoExcrN
NH4
Increase
PhytoExcrP
SRP
Decrease
BotAlgUpN NH4
Decrease
NO3
Increase
BotAlgUpP SRP
Increase
................
................
In order to avoid copyright disputes, this page is only a partial summary.
To fulfill the demand for quickly locating and searching documents.
It is intelligent file search solution for home and business.
Related download
- soil ph and liming iowa state university
- a guide to ph measurement mettler toledo
- cooling tower mathematics acid feed requirements
- injecting soda ash sodium carbonate to increase ph
- allowable adjustments to european pharmacopeia ep methods
- ph adjustment
- dosages for adjusting alkalinity
- study on ph adjustment systems and recommendations for
- acid base and abg interpretation made simple
- 1 introduction 2 ph calculation us epa
Related searches
- us epa methods
- 1 or 2 374 374 1 0 0 0 1 168 1 1 default username and password
- 1 or 2 711 711 1 0 0 0 1 168 1 1 default username and password
- 1 or 2 693 693 1 0 0 0 1 168 1 1 default username and password
- 1 or 2 593 593 1 0 0 0 1 or 2dvchrbu 168 1 1 default username and password
- 1 or 2 910 910 1 0 0 0 1 168 1 1 default username and password
- 192 1 or 2 33 33 1 0 0 0 1 1 1 default username and password
- 1 or 2 364 364 1 0 0 0 1 168 1 1 admin username and password
- 1 or 2 633 633 1 0 0 0 1 168 1 1 admin username and password
- 192 1 or 2 735 735 1 0 0 0 1 1 1 default username and password
- 1 or 2 297 297 1 0 0 0 1 168 1 1 username and password verizon
- 1 or 2 948 948 1 0 0 0 1 168 1 1 admin username and password