A first order decay function is used to calculate changes in bacteria concentrations (Bowie et al., 1985).
bactlprch,iβ=bactlprch,iβ1ββexp(βΞΌlprch,dieβ) 7:5.1.1
bactprch,iβ=bactprch,iβ1ββexp(βΞΌprch,dieβ) 7:5.1.2
where bactlprch,iβ is the amount of less persistent bacteria present in the reach on day i (#cfu/100mL), bactlprch,iβ1β is the amount of less persistent bacteria present in the reach on day iβ1 (#cfu/100mL), ΞΌlprch,dieβ is the rate constant for die-off of less persistent bacteria in streams (1/day), bactprch,iβ is the amount of persistent bacteria present in the reach on day i (#cfu/100mL), bactprch,iβ1β is the amount of persistent bacteria present in the reach on day iβ1 (#cfu/100mL), and ΞΌprch,dieβ is the rate constant for die-off of persistent bacteria in streams (1/day).
The die-off rate constants are adjusted for temperature using the equations:β
7:5.1.3
7:5.1.4
where , is the rate constant for die-off of less persistent bacteria in streams (1/day), is the rate constant for die-off of persistent bacteria in streams (1/day), is the rate constant for die-off of less persistent bacteria in streams at 20C (1/day), is the rate constant for die-off of persistent bacteria in streams at 20C (1/day), is the temperature adjustment factor for bacteria die-off/re-growth, and is the water temperature (C).
Table 7:5-1: SWAT+ input variables that pertain to bacteria die-off in the stream.
In sediment channel routing, the maximum concentration of sediment that can be transported by the water, concsed,ch,mxβ, (ton/m3 or kg/L) is compared to the concentration of sediment in the reach at the beginning of the time step, concsed,ch,iβ (Neitsch et al., 2005).
If concsed,ch,iβ < concsed,ch,mxβ, resuspension is the dominant process in the reach segment and the net amount of sediment reentrained is calculated:
seddegβ=(concsed,ch,mxββconcsed,ch,iβ)βVchββKchββCchβ 7:5.2.1
where seddegβ is the amount of sediment reentrained in the reach segment (metric tons), concsed,ch,mxβ is the maximum concentration of sediment that can be transported by the water (ton sediment/m3 H2βO or kg sediment/L H2βO), concsed,ch,iβ is the initial sediment concentration in the reach (ton sediment/m3 HO or kg sediment/L HO), is the volume of water in the reach segment (m HO), is the channel erodibility factor (cm/hr/Pa), and is the channel cover factor. When sediment resuspends, both bacteria in sediment solution and on sediment particles are released, and the net amount of bacteria released from streambed is calculated:
7:5.2.2
where is the amount of bacteria released from streambed in the reach segment (# cfu), is the amount of sediment reentrained in the reach segment (metric tons), and is the concentration of bacteria in streambed in the reach segment (# cfu/ton sediment). Bacteria concentration in streambed is calculated by the empirical regression equation, logarithmic sine function of the days of year:
7:5.2.3
where is the concentration of bacteria in streambed (# cfu/ton sediment), day is the days of year, and through are the regression coefficients in streambed bacteria concentration equation.
If > , deposition is the dominant process in the reach segment and the net amount of sediment deposited is calculated:
7:5.2.4
where is the amount of sediment deposited in the reach segment (metric tons), is the initial sediment concentration in the reach (ton sediment/m HO or kg sediment/L HO), is the maximum concentration of sediment that can be transported by the water (ton sediment/m HO or kg sediment/L HO), and is the volume of water in the reach segment (m HO). When suspended sediment deposits, bacteria on settling sediment particles are deposited, and the net amount of bacteria settled from stream water is calculated (Bai and Lung, 2005):
7:5.2.5
where is the amount of bacteria settled from stream water in the reach segment (# cfu), is the amount of bacteria in the stream water in the reach segment at the beginning of the time period (# cfu), is the linear partitioning coefficient of bacteria between the suspended sediment and water (m HO/ton sediment or L HO/kg sediment), is the amount of sediment deposited in the reach segment (metric tons), is the volume of water in the reach segment (m HO), and is the initial sediment concentration in the reach (ton sediment/m HO or kg sediment/L HO). The linear partitioning coefficient is calculated from the empirical regression equation (Pachepsky et al., 2006):
7:5.2.6
where is the linear partitioning coefficient of bacteria onto the suspended sediment (m HO/ton sediment or L HO/kg sediment) and is the percentage of clay in suspended sediment in stream water in the reach segment (%). clay normally varies between 2 and 50%.
Once the amount of bacteria released and settled has been calculated, the final amount of sediment in the reach is determined:
7:5.2.7
where is the amount of bacteria in the stream water in the reach segment (# cfu), is the amount of bacteria in the stream water in the reach segment at the beginning of the time period (# cfu), is the amount of bacteria released from streambed in the reach segment (# cfu), and is the amount of bacteria settled from stream water in the reach segment (# cfu).
The final bacteria concentration in the reach is calculated:
7:5.2.8
where is the concentration of bacteria in the stream water in the reach segment (# cfu/100 mL), is the amount of bacteria in the stream water in the reach segment (# cfu), and is the volume of water in the reach segment (m HO).
WDPRCH
ΞΌprch,die,20β: Die-off factor for persistent bacteria in streams at 20Β°C (1/day)
.bsn
WDLPRCH
ΞΌlprch,die,20β: Die-off factor for less persistent bacteria in streams at 20Β°C (1/day)
.bsn
THBACT
ΞΈbactβ: Temperature adjustment factor for bacteria die-off/growth
.bsn
SWAT+ calculates loading of pathogens and indicator bacteria for pathogens from land areas in the watershed. In the stream, bacteria die-off is the only process modeled.