arrow-left

All pages
gitbookPowered by GitBook
1 of 5

Loading...

Loading...

Loading...

Loading...

Loading...

Nitrogen Cycle

In aerobic water, there is a stepwise transformation from organic nitrogen to ammonia, to nitrite, and finally to nitrate. Organic nitrogen may also be removed from the stream by settling. This section summarizes the equations used to simulate the nitrogen cycle in the stream.

Organic Nitrogen

The amount of organic nitrogen in the stream may be increased by the conversion of algal biomass nitrogen to organic nitrogen. Organic nitrogen concentration in the stream may be decreased by the conversion of organic nitrogen to NH4+^+_44+​ or the settling of organic nitrogen with sediment. The change in organic nitrogen for a given day is:

Ξ”orgNstr=(Ξ±1βˆ—Οaβˆ—algaeβˆ’Ξ²N,3βˆ—orgNstrβˆ’Οƒ4βˆ—orgNstr)βˆ—TT\Delta orgN_{str}=(\alpha_1 * \rho_a*algae-\beta_{N,3}*orgN_{str}-\sigma_4*orgN_{str})*TTΞ”orgNstr​=(Ξ±1β€‹βˆ—Οaβ€‹βˆ—algaeβˆ’Ξ²N,3β€‹βˆ—orgNstrβ€‹βˆ’Οƒ4β€‹βˆ—orgNstr​)βˆ—TT 7:3.2.1

where Ξ”orgNstr\Delta_{orgN_{str}}Ξ”orgNstr​​ is the change in organic nitrogen concentration (mg N/L), Ξ±1\alpha_1Ξ±1​ is the fraction of algal biomass that is nitrogen (mg N/mg alg biomass), ρa\rho_aρa​ is the local respiration or death rate of algae (dayβˆ’1^{-1}βˆ’1 or hrβˆ’1^{-1}βˆ’1), algaealgaealgae is the algal biomass concentration at the beginning of the day (mg alg/L), Ξ²N,3\beta_{N,3}Ξ²N,3​ is the rate constant for hydrolysis of organic nitrogen to ammonia nitrogen (dayβˆ’1^{-1}βˆ’1 or hrβˆ’1^{-1}βˆ’1), orgNstrorgN_{str}orgNstr​ is the organic nitrogen concentration at the beginning of the day (mg N/L), Οƒ4\sigma_4Οƒ4​ is the rate coefficient for organic nitrogen settling (dayβˆ’1^{-1}βˆ’1 or hrβˆ’1^{-1}βˆ’1), and is the flow travel time in the reach segment (day or hr). The fraction of algal biomass that is nitrogen is user-defined. Equation 7:3.1.17 describes the calculation of the local respiration rate of algae. The calculation of travel time is reviewed in Chapter 7:1.

The user defines the local rate constant for hydrolysis of organic nitrogen to NH at 20C. The organic nitrogen hydrolysis rate is adjusted to the local water temperature using the relationship:

7:3.2.2

where is the local rate constant for hydrolysis of organic nitrogen to NH (day or hr), is the local rate constant for hydrolysis of organic nitrogen to NH at 20C (day or hr), and is the average water temperature for the day or hour (C).

The user defines the rate coefficient for organic nitrogen settling at 20C. The organic nitrogen settling rate is adjusted to the local water temperature using the relationship:

7:3.2.3

where is the local settling rate for organic nitrogen (day or hr), is the local settling rate for organic nitrogen at 20C (day or hr), and is the average water temperature for the day or hour (C).

TTTTTT
4+^+_44+​
Β°\degreeΒ°
Ξ²N,3=Ξ²N,3,20βˆ—1.047(Twaterβˆ’20)\beta_{N,3}=\beta_{N,3,20}*1.047^{(T_{water}-20)}Ξ²N,3​=Ξ²N,3,20β€‹βˆ—1.047(Twaterβ€‹βˆ’20)
Ξ²N,3\beta_{N,3}Ξ²N,3​
4+^+_44+​
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
Ξ²N,3,20\beta_{N,3,20}Ξ²N,3,20​
4+^+_44+​
Β°\degreeΒ°
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
TwaterT_{water}Twater​
Β°\degreeΒ°
Β°\degreeΒ°
Οƒ4=Οƒ4,20βˆ—1.024(Twaterβˆ’20)\sigma_4=\sigma_{4,20}*1.024^{(T_{water}-20)}Οƒ4​=Οƒ4,20β€‹βˆ—1.024(Twaterβ€‹βˆ’20)
Οƒ4\sigma_4Οƒ4​
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
Οƒ4,20\sigma_{4,20}Οƒ4,20​
Β°\degreeΒ°
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
TwaterT_{water}Twater​
Β°\degreeΒ°

Nitrate

The amount of nitrate (NO3βˆ’NO_3^-NO3βˆ’β€‹) in the stream may be increased by the oxidation of NO2βˆ’NO_2^-NO2βˆ’β€‹. The nitrate concentration in the stream may be decreased by the uptake of NO3βˆ’NO_3^-NO3βˆ’β€‹ by algae. The change in nitrate for a given day is:

Ξ”NO3str=(Ξ²N,2βˆ—NO2strβˆ’(1βˆ’frNH4)βˆ—Ξ±1βˆ—ΞΌaβˆ—algae)βˆ—TT\Delta NO3_{str}=(\beta_{N,2}*NO2_{str}-(1-fr_{NH4})*\alpha_1*\mu_a*algae)*TTΞ”NO3str​=(Ξ²N,2β€‹βˆ—NO2strβ€‹βˆ’(1βˆ’frNH4​)βˆ—Ξ±1β€‹βˆ—ΞΌaβ€‹βˆ—algae)βˆ—TT 7:3.2.10

where Ξ”NO3str\Delta NO3_{str}Ξ”NO3str​ is the change in nitrate concentration (mg N/L), Ξ²N,2\beta_{N,2}Ξ²N,2​ is the rate constant for biological oxidation of nitrite to nitrate (dayβˆ’1^{-1}βˆ’1 or hrβˆ’1^{-1}βˆ’1), NO2strNO2_{str}NO2str​ is the nitrite concentration at the beginning of the day (mg N/L), frNH4fr_{NH4}frNH4​ is the fraction of algal nitrogen uptake from ammonium pool, Ξ±1\alpha_1Ξ±1​ is the fraction of algal biomass that is nitrogen (mg N/mg alg biomass), ΞΌa\mu _aΞΌa​ is the local growth rate of algae (dayβˆ’1^{-1}βˆ’1 or hrβˆ’1^{-1}βˆ’1), algaealgaealgae is the algal biomass concentration at the beginning of the day (mg alg/L), and TTTTTT is the flow travel time in the reach segment (day or hr). The local rate constant for biological oxidation of nitrite to nitrate is calculated with equation 7:3.2.9 while the fraction of algal nitrogen uptake from ammonium pool is calculated with equation 7:3.2.7. Section 7:3.1.2.1 describes the calculation of the local growth rate of algae. The calculation of travel time is reviewed in Chapter 7:1.

Table 7:3-2: SWAT+ input variables used in in-stream nitrogen calculations.

Variable Name
Definition
File Name

Ammonium

The amount of ammonium (NH4+^+_44+​) in the stream may be increased by the mineralization of organic nitrogen and diffusion of ammonium from the streambed sediments. The ammonium concentration in the stream may be decreased by the conversion of NH4+^+_44+​ to NO2βˆ’^-_22βˆ’β€‹ or the uptake of NH4+^+_44+​ by algae. The change in ammonium for a given day is:

Ξ”NH4str=(Ξ²N,3βˆ—orgNstrβˆ’Ξ²N,1βˆ—NH4str+Οƒ3(1000βˆ—depth)βˆ’frNH4βˆ—Ξ±1βˆ—ΞΌaβˆ—algae)βˆ—TT\Delta NH4_{str}=(\beta_{N,3}*orgN_{str}-\beta_{N,1}*NH4_{str}+\frac{\sigma_3}{(1000*depth)}-fr_{NH4}*\alpha_1*\mu_a*algae)*TTΞ”NH4str​=(Ξ²N,3β€‹βˆ—orgNstrβ€‹βˆ’Ξ²N,1β€‹βˆ—NH4str​+(1000βˆ—depth)Οƒ3β€‹β€‹βˆ’frNH4β€‹βˆ—Ξ±1β€‹βˆ—ΞΌaβ€‹βˆ—algae)βˆ—TT

7:3.2.4

where Ξ”NH4str\Delta NH4_{str}Ξ”NH4str​ is the change in ammonium concentration (mg N/L), Ξ²N,3\beta_{N,3}Ξ²N,3​ is the rate constant for hydrolysis of organic nitrogen to ammonia nitrogen (dayβˆ’1^{-1}βˆ’1 or hrβˆ’1^{-1}βˆ’1), orgNstrorgN_{str}orgNstr​ is the organic nitrogen concentration at the beginning of the day (mg N/L), Ξ²N,1\beta_{N,1}Ξ²N,1​ is the rate constant for biological oxidation of ammonia nitrogen (day or hr), is the ammonium concentration at the beginning of the day (mg N/L), is the benthos (sediment) source rate for ammonium (mg N/m-day or mg N/m-hr), is the depth of water in the channel (m), is the fraction of algal nitrogen uptake from ammonium pool, is the fraction of algal biomass that is nitrogen (mg N/mg alg biomass), is the local growth rate of algae (day or hr), is the algal biomass concentration at the beginning of the day (mg alg/L), and is the flow travel time in the reach segment (day or hr). The local rate constant for hydrolysis of organic nitrogen to NH is calculated with equation 7:3.2.2. Section 7:3.1.2.1 describes the calculation of the local growth rate of algae. The calculation of depth and travel time is reviewed in Chapter 7:1.

The rate constant for biological oxidation of ammonia nitrogen will vary as a function of in-stream oxygen concentration and temperature. The rate constant is calculated:

7:3.2.5

where is the rate constant for biological oxidation of ammonia nitrogen (day or hr), is the rate constant for biological oxidation of ammonia nitrogen at 20C (day or hr), is the dissolved oxygen concentration in the stream (mg O/L), and is the average water temperature for the day or hour (C). The second term on the right side of equation 7:3.2.5,, is a nitrification inhibition correction factor. This factor inhibits nitrification at low dissolved oxygen concentrations.

The user defines the benthos source rate for ammonium at 20C. The benthos source rate for ammonium nitrogen is adjusted to the local water temperature using the relationship:

7:3.2.6

where is the benthos (sediment) source rate for ammonium (mg N/m-day or mg N/m2-hr), is the benthos (sediment) source rate for ammonium nitrogen at 20C (mg N/m-day or mg N/m-hr), and is the average water temperature for the day or hour (C).

The fraction of algal nitrogen uptake from ammonium pool is calculated:

7:3.2.7

where is the fraction of algal nitrogen uptake from ammonium pool, is the preference factor for ammonia nitrogen, is the ammonium concentration in the stream (mg N/L), and is the nitrate concentration in the stream (mg N/L).

βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
NH4strNH4_{str}NH4str​
Οƒ3\sigma_3Οƒ3​
2^22
2^22
depthdepthdepth
frNH4fr_{NH4}frNH4​
Ξ±1\alpha_1Ξ±1​
ΞΌa\mu _aΞΌa​
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
algaealgaealgae
TTTTTT
4+^+_44+​
Ξ²N,1=Ξ²N,1,20βˆ—(1βˆ’exp[βˆ’0.6βˆ—Oxstr])βˆ—1.083(Twaterβˆ’20)\beta_{N,1}=\beta_{N,1,20}*(1-exp[-0.6*Ox_{str}])*1.083^{(T_{water}-20)}Ξ²N,1​=Ξ²N,1,20β€‹βˆ—(1βˆ’exp[βˆ’0.6βˆ—Oxstr​])βˆ—1.083(Twaterβ€‹βˆ’20)
Ξ²N,1\beta_{N,1}Ξ²N,1​
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
Ξ²N,1,20\beta_{N,1,20}Ξ²N,1,20​
Β°\degreeΒ°
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
OxstrOx_{str}Oxstr​
2_22​
TwaterT_{water}Twater​
Β°\degreeΒ°
(1βˆ’exp[βˆ’0.6βˆ—Oxstr])(1-exp[-0.6*Ox_{str}])(1βˆ’exp[βˆ’0.6βˆ—Oxstr​])
Β°\degreeΒ°
Οƒ3=Οƒ3,20βˆ—1.074(Twaterβˆ’20)\sigma_3=\sigma_{3,20}*1.074^{(T_{water}-20)}Οƒ3​=Οƒ3,20β€‹βˆ—1.074(Twaterβ€‹βˆ’20)
Οƒ3\sigma_3Οƒ3​
2^22
2^22
Οƒ3,20\sigma_{3,20}Οƒ3,20​
Β°\degreeΒ°
2^22
2^22
TwaterT_{water}Twater​
Β°\degreeΒ°
frNH4=fNH4βˆ—NH4str(fNH4βˆ—NH4str+(1βˆ’fNH4)βˆ—NO3str)fr_{NH4}=\frac{f_{NH4}*NH4_{str}}{(f_{NH4}*NH4_{str}+(1-f_{NH4})*NO3_{str})}frNH4​=(fNH4β€‹βˆ—NH4str​+(1βˆ’fNH4​)βˆ—NO3str​)fNH4β€‹βˆ—NH4str​​
frNH4fr_{NH4}frNH4​
fNH4f_{NH4}fNH4​
NH4strNH4_{str}NH4str​
NO3strNO3_{str}NO3str​

: Local settling rate for organic nitrogen at 20C (day)

.swq

BC1

: Rate constant for biological oxidation of ammonia nitrogen at 20C (day)

.swq

RS3

: Benthos (sediment) source rate for ammonium nitrogen at 20C (mg N/m-day or mg N/m-hr)

.swq

P_N

: Preference factor for ammonia nitrogen

.wwq

BC2

: Rate constant for biological oxidation of nitrite to nitrate at 20C (day or hr)

.swq

AI1

Ξ±1\alpha_1Ξ±1​: Fraction of algal biomass that is nitrogen (mg N/mg alg biomass)

.wwq

RHOQ

ρa,20\rho_{a,20}ρa,20​: Local algal respiration rate at 20Β°\degreeΒ°C (dayβˆ’1^{-1}βˆ’1)

.wwq

BC3

Ξ²N,3,20\beta_{N,3,20}Ξ²N,3,20​: Local rate constant for hydrolysis of organic nitrogen to NH4+^+_44+​ at 20Β°\degreeΒ°C (dayβˆ’1^{-1}βˆ’1 or hrβˆ’1^{-1}βˆ’1)

.swq

RS4

Nitrite

The amount of nitrite (NO2βˆ’NO_2^-NO2βˆ’β€‹) in the stream will be increased by the conversion of NH4+NH_4^+NH4+​ to NO2βˆ’NO_2^-NO2βˆ’β€‹ and decreased by the conversion of NO2βˆ’NO_2^-NO2βˆ’β€‹ to NO3βˆ’NO_3^-NO3βˆ’β€‹. The conversion of NO2βˆ’NO_2^-NO2βˆ’β€‹ to NO3βˆ’NO_3^-NO3βˆ’β€‹ occurs more rapidly than the conversion of NH4+NH_4^+NH4+​ to NO2βˆ’NO_2^-NO2βˆ’β€‹, so the amount of nitrite present in the stream is usually very small. The change in nitrite for a given day is:

Ξ”NO2str=(Ξ²N,1βˆ—NH4strβˆ’Ξ²N,2βˆ—NO2str)βˆ—TT\Delta NO2_{str}=(\beta_{N,1}*NH4_{str}-\beta_{N,2}*NO2_{str})*TTΞ”NO2str​=(Ξ²N,1β€‹βˆ—NH4strβ€‹βˆ’Ξ²N,2β€‹βˆ—NO2str​)βˆ—TT 7:3.2.8

where Ξ”NO2str\Delta NO2_{str}Ξ”NO2str​ is the change in nitrite concentration (mg N/L), Ξ²N,1\beta_{N,1}Ξ²N,1​ is the rate constant for biological oxidation of ammonia nitrogen (dayβˆ’1^{-1}βˆ’1 or hrβˆ’1^{-1}βˆ’1), NH4strNH4_{str}NH4str​ is the ammonium concentration at the beginning of the day (mg N/L), Ξ²N,2\beta_{N,2}Ξ²N,2​ is the rate constant for biological oxidation of nitrite to nitrate (dayβˆ’1^{-1}βˆ’1 or hrβˆ’1^{-1}βˆ’1), NO2strNO2_{str}NO2str​ is the nitrite concentration at the beginning of the day (mg N/L), and is the flow travel time in the reach segment (day or hr). The local rate constant for biological oxidation of ammonia nitrogen is calculated with equation 7:3.2.5. The calculation of travel time is reviewed in Chapter 7:1.

The rate constant for biological oxidation of nitrite to nitrate will vary as a function of in-stream oxygen concentration and temperature. The rate constant is calculated:

7:3.2.9

where is the rate constant for biological oxidation of nitrite to nitrate (day or hr), is the rate constant for biological oxidation of nitrite to nitrate at 20C (day or hr), is the dissolved oxygen concentration in the stream (mg O/L), and is the average water temperature for the day or hour (C). The second term on the right side of equation 7:3.2.9, , is a nitrification inhibition correction factor. This factor inhibits nitrification at low dissolved oxygen concentrations.

TTTTTT
Ξ²N,2=Ξ²N,2,20βˆ—(1βˆ’exp[βˆ’0.6βˆ—Oxstr])βˆ—1.047(Twaterβˆ’20)\beta_{N,2}=\beta_{N,2,20}*(1-exp[-0.6*Ox_{str}])*1.047^{(T_{water}-20)}Ξ²N,2​=Ξ²N,2,20β€‹βˆ—(1βˆ’exp[βˆ’0.6βˆ—Oxstr​])βˆ—1.047(Twaterβ€‹βˆ’20)
Ξ²N,2\beta_{N,2}Ξ²N,2​
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
Ξ²N,2,20\beta_{N,2,20}Ξ²N,2,20​
Β°\degreeΒ°
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1
OxstrOx_{str}Oxstr​
2_22​
TwaterT_{water}Twater​
Β°\degreeΒ°
(1βˆ’exp[βˆ’0.6βˆ—Oxstr])(1-exp[-0.6*Ox_{str}])(1βˆ’exp[βˆ’0.6βˆ—Oxstr​])
Οƒ4,20\sigma_{4,20}Οƒ4,20​
Β°\degreeΒ°
βˆ’1^{-1}βˆ’1
Ξ²N,1,20\beta_{N,1,20}Ξ²N,1,20​
Β°\degreeΒ°
βˆ’1^{-1}βˆ’1
Οƒ3,20\sigma_{3,20}Οƒ3,20​
Β°\degreeΒ°
2^22
2^22
fNH4f_{NH4}fNH4​
Ξ²N,2,20\beta_{N,2,20}Ξ²N,2,20​
Β°\degreeΒ°
βˆ’1^{-1}βˆ’1
βˆ’1^{-1}βˆ’1