In large subbasins with a time of concentration greater than 1 day, only a portion of the surface runoff and lateral flow will reach the main channel on the day it is generated. SWAT+ incorporates a storage feature to lag a portion of the surface runoff and lateral flow release to the main channel. Pesticides in the surface runoff and lateral flow are lagged as well.
Once the pesticide load in surface runoff and lateral flow is determined, the amount of pesticide released to the main channel is calculated:
pstsurfβ=(pstsurfβ²β+pstsurstor,iβ1β)β(1βexp[tconcββsurlagβ]) 4:3.4.1
pstlatβ=(pstlatβ²β+pstlatstor,iβ1β)β(1βexp[TTlatββ1β]) 4:3.4.2
pstsedβ=(pstsedβ²β+pstsedstor,iβ1β)β(1βexp[tconcββsurlagβ]) 4:3.4.3
where is the amount of soluble pesticide discharged to the main channel in surface runoff on a given day (kg /ha), is the amount of surface runoff soluble pesticide generated in HRU on a given day (kg /ha), is the surface runoff soluble pesticide stored or lagged from the previous day (kg /ha), is the amount of soluble pesticide discharged to the main channel in lateral flow on a given day (kg /ha), is the amount of lateral flow soluble pesticide generated in HRU on a given day (kg /ha), is the lateral flow pesticide stored or lagged from the previous day (kg /ha), is the amount of sorbed pesticide discharged to the main channel in surface runoff on a given day (kg /ha), is the sorbed pesticide loading generated in HRU on a given day (kg /ha), is the sorbed pesticide stored or lagged from the previous day (kg /ha), is the surface runoff lag coefficient, is the time of concentration for the HRU (hrs) and is the lateral flow travel time (days).
Table 4:3-4: SWAT+ input variables that pertain to pesticide lag calculations.
SURLAG
surlag: surface runoff lag coefficient
.bsn
LAT_TTIME
TTlagβ: Lateral flow travel time (days)
.hru
As surface runoff flows over the soil surface, part of the waterβs energy is used to pick up and transport soil particles. The smaller particles weigh less and are more easily transported than coarser particles. When the particle size distribution of the transported sediment is compared to that of the soil surface layer, the sediment load to the main channel has a greater proportion of clay sized particles. In other words, the sediment load is enriched in clay particles. The sorbed phase of pesticide in the soil is attached primarily to colloidal (clay) particles, so the sediment load will also contain a greater proportion or concentration of pesticide than that found in the soil surface layer.
The enrichment ratio is defined as the ratio of the concentration of sorbed pesticide transported with the sediment to the concentration in the soil surface layer. SWAT+ will calculate an enrichment ratio for each storm event, or allow the user to define a particular enrichment ratio for sorbed pesticide that is used for all storms during the simulation. To calculate the enrichment ratio, SWAT+ uses a relationship described by Menzel (1980) in which the enrichment ratio is logarithmically related to sediment concentration. The equation used to calculate the pesticide enrichment ratio, Ξ΅pst:sedβ, for each storm event is:
Ξ΅pst:sedβ=0.78β(concsed,surqβ)β0.2468 4:3.3.5
where is the concentration of sediment in surface runoff (Mg sed/m HO). The concentration of sediment in surface runoff is calculated:
4:3.3.6
where is the sediment yield on a given day (metric tons), is the HRU area (ha), and is the amount of surface runoff on a given day (mm HO).
Table 4:3-3: SWAT+ input variables that pertain to sorbed pesticide loading.
Pesticide in the soil environment can be transported in solution or attached to sediment. The partitioning of a pesticide between the solution and soil phases is defined by the soil adsorption coefficient for the pesticide. The soil adsorption coefficient is the ratio of the pesticide concentration in the soil or solid phase to the pesticide concentration in the solution or liquid phase:
4:3.1.1
where is the soil adsorption coefficient ((mg/kg)/(mg/L) or m/ton), is the concentration of the pesticide sorbed to the solid phase (mg chemical/kg solid material or g/ton), and is the concentration of the pesticide in solution (mg chemical/L solution or g/ton). The definition of the soil adsorption coefficient in equation 4:3.1.1 assumes that the pesticide sorption process is linear with concentration and instantaneously reversible.
Because the partitioning of pesticide is dependent upon the amount of organic material in the soil, the soil adsorption coefficient input to the model is normalized for soil organic carbon content. The relationship between the soil adsorption coefficient and the soil adsorption coefficient normalized for soil organic carbon content is:
Pesticide in the soluble phase may be transported with surface runoff, lateral flow or percolation. The change in the amount of pesticide contained in a soil layer due to transport in solution with flow is a function of time, concentration and amount of flow:
4:3.2.1
where is the amount of pesticide in the soil layer (kg pst/ha), is the pesticide concentration in solution (mg/L or g/ton), and is the amount of mobile water on a given day (mm HO). The amount of mobile water in the layer is the amount of water lost by surface runoff, lateral flow or percolation:
for top 10 mm 4:3.2.2
SOL_BD
Οbβ: Bulk density(Mg/m3)
.sol
PSTENR
Ξ΅pst:sedβ: Pesticide enrichment ratio
.chm
Kpβ=Kocββ100orgCβ 4:3.1.2
where Kpβ is the soil adsorption coefficient ((mg/kg)/(mg/L)), Kocβ is the soil adsorption coefficient normalized for soil organic carbon content ((mg/kg)/(mg/L) or m3/ton), and orgC is the percent organic carbon present in the soil.
Table 4:3-1: SWAT+ input variables that pertain to pesticide phase partitioning.
SOL_CBN
: Amount of organic carbon in the layer (%)
.sol
SKOC
: Soil adsorption coefficient normalized for soil organic carbon content (ml/g or (mg/kg)/(mg/L) or L/kg)
pest.dat
where wmobileβ is the amount of mobile water in the layer (mm H2βO), Qsurfβ is the surface runoff generated on a given day (mm H2βO), Qlat,lyβ is the water discharged from the layer by lateral flow (mm H2βO), and wperc,lyβ is the amount of water percolating to the underlying soil layer on a given day (mm H2βO).
The total amount of pesticide in the soil layer is the sum of the adsorbed and dissolved phases:
psts,lyβ=0.01β(CsolutionββSATlyβ+CsolidphaseββΟbββdepthlyβ) 4:3.2.4
where psts,lyβ is the amount of pesticide in the soil layer (kg pst/ha), Csolutionβ is the pesticide concentration in solution (mg/L or g/ton), SATlyβ is the amount of water in the soil layer at saturation (mm H2βO), Csolidphaseβ is the concentration of the pesticide sorbed to the solid phase (mg/kg or g/ton), Οbβ is the bulk density of the soil layer (Mg/m3), and depthlyβ is the depth of the soil layer (mm). Rearranging equation 4:3.1.1 to solve for Csolidphaseβ and substituting into equation 4:3.2.4 yields:
psts,lyβ=0.01β(CsolutionββSATlyβ+CsolutionββKpββΟbββdepthlyβ) 4:3.2.5
which rearranges to
Csolutionβ=0.01β(SATlyβ+KpββΟbββdepthlyβ)psts,lyββ 4:3.2.6
Combining equation 4:3.2.6 with equation 4:3.2.1 yields
dtdpsts,lyββ=(SATlyβ+KpββΟbββdepthlyβ)psts,lyββwmobileββ 4:3.2.7
Integration of equation 4:3.2.7 gives
psts,ly,tβ=psts,ly,oββexp[(SATlyβ+KpββΟbββdepthlyβ)βwmobileββ] 4:3.2.8
where psts,ly,tβ is the amount of pesticide in the soil layer at time t (kg pstt/ha), psts,ly,oβ is the initial amount of pesticide in the soil layer (kg pst/ha), wmobileβ is the amount of mobile water in the layer (mm H2βO), SATlyβ is the amount of water in the soil layer at saturation (mm H2βO), Kpβ is the soil adsorption coefficient ((mg/kg)/(mg/L)), Οbβ is the bulk density of the soil layer (Mg/m3), and depthlyβ is the depth of the soil layer (mm).
To obtain the amount of pesticide removed in solution with the flow, the final amount of pesticide is subtracted from the initial amount of pesticide:
pstflowβ=psts,ly,oββ(1βexp[(SATlyβ+KpββΟbββdepthlyβ)βwmobileββ]) 4:3.2.9
where pstflowβ is the amount of pesticide removed in the flow (kg pst/ha) and all other terms were previously defined.
For the top 10 mm that interacts with surface runoff, the pesticide concentration in the mobile water is calculated:
concpst,flowβ=min[pstflowβ/[wperc,surfβ+Ξ²pstβ(Qsurfβ+Qlat,surfβ)]],pstsolβ/100. 4:3.2.10
while for lower layers
concpst,flowβ=min[[pstflowβ/wmobileβ],pstsolβ/100.] 4:3.2.11
where concpst,flowβ is the concentration of pesticide in the mobile water (kg pst/ha-mm H2βO), pstflowβ is the amount of pesticide removed in the flow (kg pst/ha), Ξ²pstβ is the pesticide percolation coefficient, Qsurfβ is the surface runoff generated on a given day (mm H2βO), Qlat,lyβ is the water discharged from the layer by lateral flow (mm H2βO), wperc,lyβ is the amount of water percolating to the underlying soil layer on a given day (mm H2βO), wmobileβ is the amount of mobile water in the layer (mm H2βO), and pstsolβ is the solubility of the pesticide in water (mg/L).
Pesticide moved to the underlying layer by percolation is calculated:
pstperc,lyβ=concpst,flowββwperc,lyβ 4:3.2.12
where pstperc,lyβ is the pesticide moved to the underlying layer by percolation (kg pst/ha), concpst,flowβ is the concentration of pesticide in the mobile water for the layer (kg pst/mm H2βO), and wperc,lyβ is the amount of water percolating to the underlying soil layer on a given day (mm H2βO).
Pesticide removed in lateral flow is calculated:
pstlat,surfβ=Ξ²pstββconcpst,flowββQlat,surfβ for top 10 mm 4:3.2.13
pstlat,lyβ=concpst,flowββQlat,lyβ for lower layers 4:3.2.14
where pstlat,lyβ is the pesticide removed in lateral flow from a layer (kg pst/ha), Ξ²pstβ is the pesticide percolation coefficient, concpst,flowβ is the concentration of pesticide in the mobile water for the layer (kg pst/mm H2βO), and Qlat,lyβ is the water discharged from the layer by lateral flow (mm H2βO). The pesticide percolation coefficient allows the user to set the concentration of pesticide in runoff and lateral flow from the top 10 mm to a fraction of the concentration in percolate.
Pesticide removed in surface runoff is calculated:
pstsurfβ=Ξ²pstββconcpst,flowββQsurfβ 4:3.2.15
where pstsurfβ is the pesticide removed in surface runoff (kg pst/ha), Ξ²pstβ is the pesticide percolation coefficient, concpst,flowβ is the concentration of pesticide in the mobile water for the top 10 mm of soil (kg pst/mm H2βO), and Qsurfβ is the surface runoff generated on a given day (mm H2βO).
Table 4:3-2: SWAT+ input variables that pertain to pesticide transport in solution.
SOL_BD
: Soil bulk density (Mg m)
.sol
WSOL
: Solubility of the pesticide in water (mg/L)
pest.dat
PERCOP
: Pesticide percolation coefficient
.bsn
The transport of pesticide from land areas into streams and water bodies is a result of soil weathering and erosion processes. Excessive loading of pesticides in streams and water bodies can produce toxic conditions that harm aquatic life and render the water unfit for human consumption. This chapter reviews the algorithms governing movement of soluble and sorbed forms of pesticide from land areas to the stream network. Pesticide transport algorithms in SWAT+ were taken from EPIC (Williams, 1995).
Pesticide attached to soil particles may be transported by surface runoff to the main channel. This phase of pesticide is associated with the sediment loading from the HRU and changes in sediment loading will impact the loading of sorbed pesticide. The amount of pesticide transported with sediment to the stream is calculated with a loading function developed by McElroy et al. (1976) and modified by Williams and Hann (1978).
pstsedβ=0.001βCsolidphaseββareahruβsedββΞ΅pst:sedβ 4:3.3.1
where pstsedβ is the amount of sorbed pesticide transported to the main channel in surface runoff (kg pst/ha), Csolidphaseβ is the concentration of pesticide on sediment in the top 10 mm (g pst/ metric ton soil), sed is the sediment yield on a given day (metric tons), areahruβ is the HRU area (ha), and Ξ΅pst:sedβ is the pesticide enrichment ratio.
The total amount of pesticide in the soil layer is the sum of the adsorbed and dissolved phases:
4:3.3.2
where is the amount of pesticide in the soil layer (kg /ha), is the pesticide concentration in solution (mg/L or g/ton), is the amount of water in the soil layer at saturation (mm HO), is the concentration of the pesticide sorbed to the solid phase (mg/kg or g/ton), is the bulk density of the soil layer (Mg/m), and is the depth of the soil layer (mm). Rearranging equation 4:3.1.1 to solve for and substituting into equation 4:3.3.2 yields:
4:3.3.3
which rearranges to
4:3.3.4
where is the concentration of the pesticide sorbed to the solid phase (mg/kg or g/ton), is the soil adsorption coefficient ((mg/kg)/(mg/L) or /ton) is the amount of pesticide in the soil layer (kg /ha), is the amount of water in the soil layer at saturation (mm HO), is the bulk density of the soil layer (Mg/m), and is the depth of the soil layer (mm).