All pages
Powered by GitBook
1 of 4

Loading...

Loading...

Loading...

Loading...

Sediment In Water Bodies

SWAT+ incorporates a simple mass balance model to simulate the transport of sediment into and out of water bodies. SWAT+ defines four different types of water bodies: ponds, wetlands, reservoirs and potholes. Sediment processes modeled in ponds, wetlands, reservoirs, and potholes are identical.

When calculating sediment movement through a water body, SWAT+ assumes the system is completely mixed. In a completely mixed system, as sediment enters the water body it is instantaneously distributed throughout the volume.

Sediment Outflow

The amount of sediment transported out of the water body on a given day is calculated as a function of the final concentration. The initial suspended solid concentration is:

sedflowout=concsed,fβˆ—Vflowoutsed_{flowout}=conc_{sed,f}*V_{flowout}sedflowout​=concsed,fβ€‹βˆ—Vflowout​ 8:2.3.1

where sedflowoutsed_{flowout}sedflowout​ is the amount of sediment transported out of the water body with outflow (metric tons), concsed,fconc_{sed,f}concsed,f​ is the final sediment concentration in the water body (Mg/m3^33), and VflowoutV_{flowout}Vflowout​ is the volume of outflow from the impoundment (m3^33 H2_22​O).

Settling

Incoming sediment is deposited using a modified overflow rate model (EPA 1986, cited in Haan et al., 1994). For each day, the deposition routine begins with the computation of the detention times. The actual detention time is based upon the ratio of the impoundment volume to the outflow rate.

tD=(Ct(1βˆ’DS)Vol)Qot_D=\frac{(C_t(1-DS)Vol)}{Q_o}tD​=Qo​(Ct​(1βˆ’DS)Vol)​ 8:2.2.1

Where tDt_DtD​ is detention time (sss), CtC_tCt​ is an empirical parameter to account for impoundment geometry, hydraulic response, and stratification of the suspended sediment, DSDSDS is the dead storage (the portion of the pond are that does not contribute to settling) (Griffin et al., 1985), Vol is the average impoundment volume over the time step (ft3^33), and QOQ_OQO​ is the average outflow rate over the time step (ft3^33 s%βˆ’1^{-1}βˆ’1). The detention time required for 100% of the suspended sediment to settle out of suspension is computed form the average impoundment depth (volume / area) and the settling velocity.

The trapping efficiency is calculated as

8:2.2.2

Where is trapping efficiency (fraction), is the settling velocity (m/d), is overflow velocity (m/d). is defined as

8:2.2.3

Where is reservoir outflow in m and is reservoir surface area in .

During days of no sediment inflow the amount of suspended solid settling that occurs in the water body on a given day is calculated as a function of concentration. The initial suspended solid concentration is:

8:2.2.4

where is the initial concentration of suspended solids in the water (Mg/m), is the amount of sediment in the water body at the beginning of the day (metric tons), is the amount of sediment added to the water body with inflow (metric tons), is the volume of water stored in water body or channel at the beginning of the day (m HO), and is the volume of water entering water body on given day (m HO).

Settling occurs only when the sediment concentration in the water body exceeds the equilibrium sediment concentration specified by the user, . The concentration of sediment in the water body at the end of the day is calculated:

​

if 8:2.2.5

if 8:2.2.6

where is the final sediment concentration in the water body (Mg/m), is the initial concentration of suspended solids in the water body (Mg/m), is the equilibrium concentration of suspended solids in the water body (Mg/m), is the decay constant (1/day), is the length of the time step (1 day), and is the median particle size of the inflow sediment (Β΅m). Assuming 99% of the 1 Β΅m size particles settle out of solution within 25 days, is equal to 0.184.

For ponds, wetlands, and potholes, the median particle size of the inflow sediment is calculated:

8:2.2.7

where is the median particle size of the inflow sediment (m), is percent clay in the surface soil layer in the subbasin, is the percent silt in the surface soil layer in the subbasin, is the percent sand in the surface soil layer in the subbasin. Because reservoirs are located on the main channel network and receive sediment from the entire area upstream, defaulting the sand, silt, and clay fractions to those of a single subbasin or HRU in the upstream area is not appropriate. Instead the user is allowed to set the median particle size diameter to a representative value for reservoirs.

The amount of sediment settling out of solution on a given day is then calculated:

8:2.2.8

where is the amount of sediment removed from the water by settling (metric tons), is the initial concentration of suspended solids in the water body (Mg/m), is the final sediment concentration in the water body (Mg/m), and is the volume of water in the impoundment (m HO).

Table 8:2-1: SWAT+ input variables that pertain to sediment settling.

Variable Name
Definition
Input File

: Equilibrium sediment concentration in water body (mg/L)

.hru

CLAY

: Percent clay in the surface soil layer in the subbasin

.sol

SILT

: Percent silt in the surface soil layer in the subbasin

.sol

SAND

: Percent sand in the surface soil layer in the subbasin

.sol

RES_D50

: Median particle size of sediment in a reservoir

.res

trappeff=Vsetl/Vovfltrappeff=V_{setl}/V_{ovfl}trappeff=Vsetl​/Vovfl​
trappefftrappefftrappeff
VsetlV_{setl}Vsetl​
VovflV_{ovfl}Vovfl​
VovflV_{ovfl}Vovfl​
Vovfl=(Qo/SAres)(10,000)V_{ovfl}=\frac{(Q_o/SA_{res})}{(10,000)}Vovfl​=(10,000)(Qo​/SAres​)​
QoQ_oQo​
3^33
SAresSA_{res}SAres​
hahaha
concsed,i=(sedwb,i+sedflowin)(Vstored+Vflowin)conc_{sed,i}=\frac{(sed_{wb,i}+sed_{flowin})}{(V_{stored}+V_{flowin})}concsed,i​=(Vstored​+Vflowin​)(sedwb,i​+sedflowin​)​
concsed,iconc_{sed,i}concsed,i​
3^33
sedwb,ised_{wb,i}sedwb,i​
sedflowinsed_{flowin}sedflowin​
VstoredV_{stored}Vstored​
3^33
2_22​
VflowinV_{flowin}Vflowin​
3^33
2_22​
concsed,eqconc_{sed,eq}concsed,eq​
concsed,f=(concsed,iβˆ’concsed,eq)βˆ—exp[βˆ’ksβˆ—tβˆ—d50]+concsed,eqconc_{sed,f}=(conc_{sed,i}-conc_{sed,eq})*exp[-k_s*t*d_{50}]+conc_{sed,eq}concsed,f​=(concsed,iβ€‹βˆ’concsed,eq​)βˆ—exp[βˆ’ksβ€‹βˆ—tβˆ—d50​]+concsed,eq​
concsed,i>concsed,eqconc_{sed,i}>conc_{sed,eq}concsed,i​>concsed,eq​
concsed,f=concsed,iconc_{sed,f}=conc_{sed,i}concsed,f​=concsed,i​
concsed,i≀concsed,eqconc_{sed,i} \le conc_{sed,eq}concsed,i​≀concsed,eq​
concsed,fconc_{sed,f}concsed,f​
3^33
concsed,iconc_{sed,i}concsed,i​
3^33
concsed,eqconc_{sed,eq}concsed,eq​
3^33
ksk_sks​
ttt
d50d_{50}d50​
ksk_sks​
d50=exp(0.41βˆ—mc100+2.71βˆ—msilt100+5.7βˆ—ms100)d_{50}=exp(0.41*\frac{m_c}{100}+2.71*\frac{m_{silt}}{100}+5.7*\frac{m_s}{100})d50​=exp(0.41βˆ—100mc​​+2.71βˆ—100msilt​​+5.7βˆ—100ms​​)
d50d_{50}d50​
ΞΌ\muΞΌ
mcm_cmc​
msiltm_{silt}msilt​
msm_sms​
sedstl=(concsed,iβˆ’concsed,f)βˆ—Vsed_{stl}=(conc_{sed,i}-conc_{sed,f})*Vsedstl​=(concsed,iβ€‹βˆ’concsed,f​)βˆ—V
sedstlsed_{stl}sedstl​
concsed,iconc_{sed,i}concsed,i​
3^33
concsed,fconc_{sed,f}concsed,f​
3^33
VVV
3^33
2_22​

RES_NSED

concsed,eqconc_{sed,eq}concsed,eq​: Equilibrium sediment concentration in water body (mg/L)

.res

PND_NSED

concsed,eqconc_{sed,eq}concsed,eq​: Equilibrium sediment concentration in water body (mg/L)

.pnd

WET_NSED

concsed,eqconc_{sed,eq}concsed,eq​: Equilibrium sediment concentration in water body (mg/L)

.pnd

POT_NSED

Mass Balance

The mass balance equation for sediment in a water body is:

8:2.1.1

where is the amount of sediment in the water body at the end of the day (metric tons), is the amount of sediment in the water body at the beginning of the day (metric tons), is the amount of sediment added to the water body with inflow (metric tons), is the amount of sediment removed from the water by settling (metric tons), is the amount of sediment transported out of the water body with outflow (metric tons).

concsed,eqconc_{sed,eq}concsed,eq​
mcm_cmc​
msiltm_{silt}msilt​
msm_sms​
d50d_{50}d50​
sedwb=sedwb,i+sedflowinβˆ’sedstlβˆ’sedflowoutsed_{wb}=sed_{wb,i}+sed_{flowin}-sed_{stl}-sed_{flowout}sedwb​=sedwb,i​+sedflowinβ€‹βˆ’sedstlβ€‹βˆ’sedflowout​
sedwbsed_{wb}sedwb​
sedwb,ised_{wb,i}sedwb,i​
sedflowinsed_{flowin}sedflowin​
sedstlsed_{stl}sedstl​
sedflowoutsed_{flowout}sedflowout​