Reaeration occurs by diffusion of oxygen from the atmosphere into the stream and by the mixing of water and air that occurs during turbulent flow.
Reareation will occur when water falls over a dam, weir, or other structure in the stream. To simulate this form of reaeration, a βstructureβ command line is added in the watershed configuration file (.fig) at every point along the stream where flow over a structure occurs.
The amount of reaeration that occurs is a function of the oxygen deficit above the structure and a reaeration coefficient:
ΞOxstrβ=DaββDbβ=Daβ(1βrea1β) 7:3.5.10
where ΞOxstrβ is the change in dissolved oxygen concentration (mg O2β/L), Daβ is the oxygen deficit above the structure (mg O2β/L), Dbβ is the oxygen deficit below the structure (mg O2β/L), and rea is the reaeration coefficient.
The oxygen deficit above the structure, , is calculated:
7:3.5.11
where is the equilibrium saturation oxygen concentration (mg O/L), and is the dissolved oxygen concentration in the stream (mg O/L).
Butts and Evans (1983) documents the following relationship that can be used to estimate the reaeration coefficient:
7:3.5.12
where is the reaeration coefficient, is an empirical water quality factor, is an empirical dam aeration coefficient, is the height through which water falls (m), and is the average water temperature (C).
The empirical water quality factor is assigned a value based on the condition of the stream:
= 1.80 in clean water
= 1.60 in slightly polluted water
= 1.00 in moderately polluted water
= 0.65 in grossly polluted water
The empirical dam aeration coefficient is assigned a value based on the type of structure:
= 0.70 to 0.90 for flat broad crested weir
= 1.05 for sharp crested weir with straight slope face
= 0.80 for sharp crested weir with vertical face
= 0.05 for sluice gates with submerged discharge
Table 7:3-5: SWAT+ input variables used in in-stream oxygen calculations.
The user defines the reaeration rate at 20Β°C. The reaeration rate is adjusted to the local water temperature using the relationship:
ΞΊ2β=ΞΊ2,20ββ1.024(Twaterββ20) 7:3.5.4
where ΞΊ2β is the reaeration rate (dayβ1 or hrβ1), ΞΊ2,20β is the reaeration rate at 20Β°C (dayβ1 or hrβ1), and Twaterβ is the average water temperature for the day or hour (Β°C).
Numerous methods have been developed to calculate the reaeration rate at 20Β°C, ΞΊ2,20β. A few of the methods are listed below. Brown and Barnwell (1987) provide additional methods.
Using field measurements, Churchill, Elmore and Buckingham (1962) derived the relationship:
7:3.5.5
where is the reaeration rate at 20C (day), is the average stream velocity (m/s), and is the average stream depth (m).
OβConnor and Dobbins (1958) incorporated stream turbulence characteristics into the equations they developed. For streams with low velocities and isotropic conditions,
7:3.5.6
where is the reaeration rate at 20C (day), is the molecular diffusion coefficient (m/day), is the average stream velocity (m/s), and is the average stream depth (m). For streams with high velocities and nonisotropic conditions,
7:3.5.7
where is the reaeration rate at 20C (day), is the molecular diffusion coefficient (m/day), is the slope of the streambed (m/m), and is the average stream depth (m). The molecular diffusion coefficient is calculated
7:3.5.8
where is the molecular diffusion coefficient (m/day), and is the average water temperature (C).
Owens et al. (1964) developed an equation to determine the reaeration rate for shallow, fast moving streams where the stream depth is 0.1 to 3.4 m and the velocity is 0.03 to 1.5 m/s.
7:3.5.9
where is the reaeration rate at 20C (day), is the average stream velocity (m/s), and is the average stream depth (m).
: Local algal respiration rate at 20C (day)
.wwq
RK1
: CBOD deoxygenation rate at 20C (day)
.swq
RK4
:Sediment oxygen demand rate at 20C(mg O/(m.day))
.swq
AI5
: Rate of oxygen uptake per unit NH oxidation (mg O/mg N)
.wwq
AI6
: Rate of oxygen uptake per unit NO oxidation (mg O/mg N)
.wwq
AERATION_COEF
: Reaeration coefficient
.fig
RK2
ΞΊ2,20β: Reaeration rate at 20Β°C (dayβ1)
.swq
AI3
Ξ±3β: Rate of oxygen production per unit algal photosynthesis (mg O2β/mg alg)
.wwq
AI4
Ξ±4β: Rate of oxygen uptake per unit algal respiration (mg O2β/mg alg)
.wwq
RHOQ