arrow-left

All pages
gitbookPowered by GitBook
1 of 51

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

Loading...

swift_out

SWIFT input file

Code
Option

0

Do not write SWIFT input file

1

Write SWIFT input file (swift_hru.inp)

nostress

Plant stress

Code
Option

0

All plant stresses applied

1

Turn off all plant stress

2

Turn off nutrient plant stress only

gampt

Surface runoff method

Code
Option

0

Curve Number

1

Green & Ampt

qual2e

Instream nutrient routing method

Code
Option

0

Instream nutrient routing using QUAL2E

1

Instream nutrient routing using QUAL2E with simplified nutrient transformations

p_perc

Phosphorus percolation coefficient

The phosphorus percolation coefficient is the ratio of the solution phosphorus concentration in the top 10 mm of soil to the concentration of phosphorus in the percolate.

rsd_decomp

Residue decomposition coefficient

The fraction of residue that will decompose in a day assuming optimal moisture, temperature, C:N ratio, and C:P ratio.

pest_perc

Pesticide percolation coefficient

This parameter controls the amount of pesticide removed from the surface layer in runoff and lateral flow relative to the amount removed via percolation. This parameter is only important if pesticide transport is simulated.

circle-info

The lower pest_perc, the lower the concentration of the pesticide in runoff. If pest_perc = 1.0, the surface runoff will have the same pesticide concentration as the percolate.

msk_co1

Coefficient to control the impact of the storage time constant for normal flow on the overall storage time constant for the channel

Normal flow is defined as the streamflow when the channel is at bankfull depth. This parameter is only important if channel routing is simulated using the Muskingum routing method.

msk_x

Weighting factor control relative importance of inflow rate and outflow rate in determining storage on reach

The weighting factor is a function of the wedge storage. This parameter is only important if channel routing is simulated using the Muskingum routing method.

circle-info

For reservoir-type storage, there is no wedge and msk_x should be 0.0. For a full-wedge, msk_x should be 0.5. For rivers, msk_x will fall between 0.0 and 0.3 with a mean value near 0.2.

denit_frac

Denitrification threshold water content

This parameter defines the fraction of field capacity water content above which denitrification takes place. Denitrification is the bacterial reduction of nitrate (NO3) to N2 or N2O gases under anaerobic (reduced) conditions. Because SWAT+ does not track the redox status of the soil layers, the presence of anaerobic conditions in a soil layer is defined by this variable. If the soil water content calculated as a fraction of field capacity is ≥ denit_frac, then anaerobic conditions are assumed to be present and denitrification is modeled. If the soil water content calculated as a fraction of field capacity is < denit_frac, then aerobic conditions are assumed to be present and denitrification is not modeled.

cn_froz

Parameter for frozen soil adjustment on infiltration/runoff

plaps

Precipitation lapse rate

A positive value denotes an increase in precipitation with an increase in elevation while a negative value denotes a decrease in precipitation with an increase in elevation.

adj_uhyd

Adjustment factor for subdaily unit hydrograph basetime

tlaps

Temperature lapse rate

A positive value denotes a decrease in temperature with an increase in elevation.

uhyd_alpha

Alpha coefficient for gamma function unit hydrograph

This parameter is required if = 1.

uhyd

n_fix_max

Maximum daily nitrogen fixation

cov_mgt

Scaling parameter for cover and management factor for overland flow erosion

rte_cha

Channel water routing method

There are two channel water routing methods available in SWAT+:

Code
Option

0

Variable Storage method

1

Muskingum method

circle-info

The user must be careful to define , and in when the Muskingum method is chosen.

msk_co1
msk_co2
msk_x
parameters.bsn

pet

Potential evapotranspiration (PET) method

Numerous methods exist to calculate potential evapotranspiration. Three of the most widely-used ones are included in SWAT+, the Priestley-Taylor, Penman/Monteith, and Hargreaves equations. The codes for the three methods are listed in the table below. If a method other than Priestley-Taylor, Penman/Monteith, or Hargreaves is recommended for the area, in which the watershed is located, the user can calculate daily PET values with the recommended method and import them into SWAT+ (using pet in weather-sta.cli).

Code
Option

0

Priestley-Taylor method

1

Penman/Monteith method

2

Hargreaves method

carbon

Carbon routine

Code
Option

0

Static soil carbon

1

C-FARM one carbon pool model

2

Century model

Papers Xuesong, Armen

lapse

Precipitation and temperature lapse rate control

Code
Option

0

Do not adjust precipitation and temperature for elevation

1

Adjust precipitation and temperature for elevation

crack

Crack flow code

Code
Option

0

Do not compute crack flow in soil

1

Compute crack flow in soil

If crack = 1, the crack volume potential is controlled by perc_crk in soils.sol.

codes.bsn

This file contains control codes for the simulation of basin-level processes.

Field
Description
Type

Basin

General watershed attributes are defined in two basin input files, and . These attributes control a diversity of physical processes at the watershed level. Users can use the default values set by the interfaces or change them to better reflect what is happening in a given watershed.

circle-info

Even if nutrients are not being studied in a watershed, some attention must be paid to basin nutrient variables, because nutrient cycling impacts plant growth, which in turn affects the hydrologic cycle. Variables governing bacteria or pesticide transport need to be initialized only if these processes are being modeled in the watershed.

wtable

Water table depth algorithm code

Code
Option

0

Simulate shallow water table using original water table depth routine (fill to upper limit)

1

Simulate shallow water table using DRAINMOD water table depth routine

codes.bsn
parameters.bsn

tiledrain

Tile drainage equation code

Code
Option

0

Simulate tile flow using drawdown days equation

1

Simulate tile flow using DRAINMOD equations

Papers Daniel, Tássia

sw_init

Initial soil water storage expressed as a fraction of field capacity water content

All soils in the watershed will be initialized to the same fraction. If sw_init = 0.0, the model will calculate it as a function of average annual precipitation.

circle-info

We recommend using a warm-up period of at least 1 year, i.e. start the simulation at least 1 year prior to the period of interest. This allows the model to get the water cycling properly before any comparisons between measured and simulated data are made. If a warm-up period is incorporated, the value for sw_init will not impact model results.

integer

event

Currently not used

integer

Crack flow

integer

Writing of input file for SWIFT

integer

sed_det

Currently not used

integer

Channel routing

integer

deg_cha

Currently not used

integer

wq_cha

Currently not used

integer

Turning off of plant stress

integer

cn

Currently not used

integer

c_fact

Currently not used

integer

Carbon routine

integer

Precipitation and temperature lapse rate control

integer

Unit Hydrograph method

integer

sed_cha

Currently not used

integer

Tile drainage equation code

integer

Water table depth algorithms

integer

Soil phosphorus model

integer

Surface runoff method

integer

atmo_dep

Currently not used

string

stor_max

Currently not used

integer

qual2e

Instream nutrient routing method

integer

Flood routing

integer

pet_file

Currently not used

string

wq_file

Currently not used

string

pet

Potential Evapotranspiration (PET) method

uhyd

Unit Hydrograph method

Code
Option

The Unit Hydrograph method is only relevant when simulating at a sub-daily timestep and = 0.

crack
swift_out
rte_cha
nostress
carbon
lapse
uhyd
tiledrain
wtable
soil_p
gampt
gwflow

0

Triangular Unit Hydrograph

1

Gamma Function Unit Hydrograph

gampt

gwflow

Flood routing

Code
Option

0

GWFlow module not active

1

GWFlow module active

soil_p

Soil phosphorus model

Code
Option

0

Original soil phosphorus model

1

New soil phosphorus model

Vadas & White (2010)

adj_pkrt

Peak rate adjustment factor for sediment routing in the subbasin (tributary channels)

Sediment routing is a function of peak flow rate and mean daily flow. Because SWAT originally could not directly calculate the sub-daily hydrograph due to the use of precipitation summarized on a daily basis, this variable was incorporated to allow adjustment for the effect of the peak flow rate on sediment routing. This factor is used in the MUSLE equation and impacts the amount of erosion generated in the HRUs.

surq_lag

Surface runoff lag coefficient

In large routing units with a time of concentration greater than 1 day, only a portion of the surface runoff will reach the main channel on the day it is generated. SWAT+ incorporates a surface runoff storage feature to lag a portion of the surface runoff release to the main channel.

This parameter controls the fraction of the total available water that will be allowed to enter the reach on any one day. For a given time of concentration, as surq_lag decreases in value more water is held in storage. The delay in release of surface runoff will smooth the streamflow hydrograph simulated in the reach.

orgn_min

Rate factor for humus mineralization of active organic nutrients (N and P)

Please refer to the Theoretical Documentation for a description of the use of this parameter in the mineralization calculations.

n_uptake

Nitrogen uptake distribution parameter

Root density is greatest near the surface, and plant nitrogen uptake in the upper portion of the soil will be greater than in the lower portion. The depth distribution of nitrogen uptake is controlled by n_uptake, the nitrogen uptake distribution parameter.

The importance of the nitrogen uptake distribution parameter lies in its control over the maximum amount of nitrate removed from the upper layers. Because the top 10 mm of the soil profile interacts with surface runoff, the nitrogen uptake distribution parameter will influence the amount of nitrate available for transport in surface runoff. The model allows lower layers in the root zone to fully compensate for lack of nitrate in the upper layers, so there should not be significant changes in nitrogen stress with variation in the value used for n_uptake.

p_soil

Phosphorus soil partitioning coefficient

The phosphorus soil partitioning coefficient p_soil is the ratio of the soluble phosphorus concentration in the top 10 mm of soil to the concentration of soluble phosphorus in surface runoff.

The primary mechanism of phosphorus movement in the soil is by diffusion. Diffusion is the migration of ions over small distances (1-2 mm) in the soil solution in response to a concentration gradient. Due to the low mobility of solution phosphorus, surface runoff will only partially interact with the solution P stored in the top 10 mm of soil.

parameters.bsn

This file contains basin-level parameters.

Field
Description
Type
Unit
Default
Range

lai_noevap

Currently not used

real

sw_init

Initial soil water storage expressed as a fraction of field capacity water content

real

0

0-1

surq_lag

Surface runoff lag coefficient

real

4

1-24

adj_pkrt

Peak rate adjustment factor for sediment routing in the subbasin (tributary channels)

real

1

0.5-2

prf

Peak rate adjustment factor for sediment routing in the main channel

real

1

0-2

lin_sed

Currently not used

real

exp_sed

Currently not used

real

orgn_min

Rate factor for humus mineralization of active organic nutrients (N and P)

real

0.0003

0.001-0.003

n_uptake

Nitrogen uptake distribution parameter

real

20

0-100

p_uptake

Phosphorus uptake distribution parameter

real

20

0-100

n_perc

Nitrate percolation coefficient

real

0.2

0-1

p_perc

Phosphorus percolation coefficient

real

10m^3/M

10

10-17.5

p_soil

Phosphorus soil partitioning coefficient

real

m^3/Mg

175

100-200

p_avail

Phosphorus availability index

real

0.4

0.01-0.7

rsd_decomp

Residue decomposition coefficient

real

0.05

0.02-0.1

pest_perc

Pesticide percolation coefficient

real

0.5

0-1

msk_co1

Coefficient to control the impact of the storage time constant for normal flow on the overall storage time constant for the channel

real

0.75

0-10

msk_co2

Coefficient to control the impact of the storage time constant for low flow on the overall storage time constant for the channel

real

0.25

0-10

msk_x

Weighting factor control relative importance of inflow rate and outflow rate in determining storage on reach

real

0.2

0-0.3

nperco_lchtile

Nitrogen concentration coefficient for tile flow and leaching from bottom layer

real

0

0-1

evap_adj

Reach evaporation adjustment factor

real

0.6

0.5-1

scoef

Currently not used

real

denit_exp

Denitrification exponential rate coefficient

real

1.4

0-3

denit_frac

Denitrification threshold water content

real

1.3

0-1

man_bact

Currently not used

real

adj_uhyd

Adjustment factor for subdaily unit hydrograph basetime

real

0

0-1

cn_froz

Parameter for frozen soil adjustment on infiltration/runoff

real

0.000862

0-0

dorm_hr

Time threshold used to define dormancy

real

hrs

0

0-24

plaps

Precipitation lapse rate

real

mm/km

tlaps

Temperature lapse rate

real

deg C/km

n_fix_max

Maximum daily nitrogen fixation

real

kg/ha

20

1-20

rsd_decay

Minimum daily residue decay

real

fraction

0.01

0-0.05

rsd_cover

Currently not used

real

urb_init_abst

Maximum initial abstraction for urban areas

real

5

0-10

petco_pmpt

Currently not used

real

uhyd_alpha

Alpha coefficient for gamma function unit hydrograph

real

5

0.5-10

splash

Splash erosion coefficient

real

1

0.9-3.1

rill

Rill erosion coefficient

real

0.7

0.5-2

surq_exp

Exponential coefficient for overland flow

real

1.2

1-3

cov_mgt

Scaling parameter for cover and management factor for overland flow erosion

real

0.03

0.001-0.45

cha_d50

Currently not used

real

co2

CO2 concentration at start of simulation

real

ppm

1.57

1-5

day_lag_max

Currently not used

real

igen

Currently not used

integer

p_uptake

Phosphorus uptake distribution parameter

This parameter controls plant uptake of phosphorus from the different soil horizons in the same way that n_uptake controls nitrogen uptake.

Phosphorus removed from the soil by plants is taken from the solution phosphorus pool. The importance of the phosphorus uptake distribution parameter lies in its control over the maximum amount of solution P removed from the upper layers. Because the top 10 mm of the soil profile interacts with surface runoff, the phosphorus uptake distribution parameter will influence the amount of labile phosphorus available for transport in surface runoff. The model allows lower layers in the root zone to fully compensate for lack of solution P in the upper layers, so there should not be significant changes in phosphorus stress with variation in the value used for p_uptake.

evap_adj

Reach evaporation adjustment factor

The evaporation coefficient is a calibration parameter for the user that was created to allow reach evaporation to be dampened in arid regions. The original equation tends to overestimate evaporation in these areas.

prf

Peak rate factor

Sediment routing is a function of peak flow rate and mean daily flow. Because SWAT originally could not directly calculate the sub-daily hydrograph, this variable was incorporated to allow adjustment for the effect of the peak flow rate on sediment routing. This variable impacts channel degradation.

n_perc

Nitrate percolation coefficient

This parameter controls the amount of nitrate removed from the surface layer in runoff relative to the amount removed via percolation.

circle-info

The smaller n_perc, the lower the concentration of nitrate in runoff. If n_perc = 1.0, the surface runoff will have the same nitrate concentration as the percolate.

msk_co2

Coefficient to control the impact of the storage time constant for low flow on the overall storage time constant for the channel

Normal flow is defined as the streamflow when the channel is at 0.1*bankfull depth. This parameter is only important if channel routing is simulated using the Muskingum routing method.

nperco_lchtile

Nitrogen concentration coefficient for tile flow and leaching from bottom layer

splash

Splash erosion coefficient

rill

Rill erosion coefficient

Multiplier to USLE_K for soils susceptible to rill erosion.

rsd_decay

Minimum daily residue decay

urb_init_abst

Maximum initial abstraction for urban areas

This parameter is only relevant when using Green & Ampt.

surq_exp

Exponential coefficient for overland flow

dorm_hr

Time threshold used to define dormancy

The maximum day length minus dorm_hr is equal to when dormancy occurs.

co2

CO2 concentration in the atmosphere at start of simulation

p_avail

Phosphorus availability index

Many studies have shown that after an application of soluble P fertilizer, solution P concentration decreases rapidly with time due to reaction with the soil. This initial "fast" reaction is followed by a much slower decrease in solution P that may continue for several years (Barrow and Shaw, 1975; Munns and Fox, 1976; Rajan and Fox, 1972; Sharpley, 1982). In order to account for the initial rapid decrease in solution P, SWAT+ assumes a rapid equilibrium exists between solution P and an "active" mineral pool. The subsequent slow reaction is simulated by the slow equilibrium assumed to exist between the "active" and "stable" mineral pools. The algorithms governing movement of inorganic phosphorus between these three pools are taken from Jones et al. (1984).

Equilibration between the solution and active mineral pool is governed by the phosphorus availability index. This index specifies the fraction of fertilizer P which is in solution after an incubation period, i.e. after the rapid reaction period.

A number of methods have been developed to measure the phosphorus availability index. Jones et al. (1984) recommends a method outlined by Sharpley et al. (1984) in which various amounts of phosphorus are added in solution to the soil as K2HPO4. The soil is wetted to field capacity and then dried slowly at 25°C. When dry, the soil is rewetted with deionized water. The soil is exposed to several wetting and drying cycles over a 6-month incubation period. At the end of the incubation period, solution phosphorus is determined by extraction with anion exchange resin.

The P availability index is then calculated as:

where pai is the phosphorus availability index, is the amount of phosphorus in solution after fertilization and incubation, is the amount of phosphorus in solution before fertilization, and is the amount of soluble P fertilizer added to the sample.

hashtag
References

Barrow, N.J. and T.C. Shaw. 1975. The slow reactions between soil and anions. 2. Effect of time and temperature on the decrease in phosphate concentration in soil solution. Soil Science 119(2): 167-177.

Jones et al. 1984

Munns and Fox. 1976

Rajan and Fox. 1972

Sharpley. 1982

Sharpley et al. 1984

denit_exp

Denitrification exponential rate coefficient

This coefficient allows the user to control the rate of denitrification.

Impact of denit_exp value on amount of nitrogen lost to denitrification assuming initial nitrate content in layer is 200 kg/ha, temperature of layer is 10ºC, and organic carbon content of layer is 2%.
pai=Psolution,f−Psolution,ifertminPpai=\frac{P_{solution,f}-P_{solution,i}}{fert_{minP}}pai=fertminP​Psolution,f​−Psolution,i​​
Psolution,fP_{solution,f}Psolution,f​
Psolution,iP_{solution,i}Psolution,i​
fertminPfert_{minP}fertminP​
Linkarrow-up-right