SWAT+ Documentation
  • Introduction to SWAT+
  • Watershed Configuration
    • Spatial Objects
  • Calibration
  • SWAT+ Input Files
    • Input File Format
    • Master File (file.cio)
    • Simulation Settings
      • object.cnt
      • time.sim
        • day_start
        • yrc_start
        • day_end
        • yrc_end
        • step
      • print.prt
        • nyskip
        • interval
        • aa_int_cnt
        • hydcon
        • object
      • object.prt
        • obj_typ
        • obj_typ_no
        • hyd_typ
    • Climate
      • weather-sta.cli
        • name (weather_sta.cli)
        • wgn
        • pcp
        • tmp
        • slr
        • hmd
        • wnd
        • pet
        • atmo_dep
      • weather-wgn.cli
        • name (weather_wgn.cli)
        • yrs_pcp
        • tmp_max_ave
        • tmp_min_ave
        • tmp_max_sd
        • tmp_min_sd
        • pcp_ave
        • pcp_sd
        • pcp_skew
        • wet_dry
        • wet_wet
        • pcp_days
        • pcp_hhr
        • slr_ave
        • dew_ave
        • wnd_ave
      • pcp.cli and pcp data files
      • tmp.cli and tmp data files
      • hmd.cli and hmd data files
      • slr.cli and slr data files
      • wnd.cli and wnd data files
      • atmo.cli
        • timestep
        • num_aa
    • Basin
      • codes.bsn
        • pet
        • crack
        • swift_out
        • rte_cha
        • nostress
        • carbon
        • lapse
        • uhyd
        • tiledrain
        • wtable
        • soil_p
        • gampt
        • qual2e
        • gwflow
      • parameters.bsn
        • sw_init
        • surq_lag
        • adj_pkrt
        • prf
        • orgn_min
        • n_uptake
        • p_uptake
        • n_perc
        • p_perc
        • p_soil
        • p_avail
        • rsd_decomp
        • pest_perc
        • msk_co1
        • msk_co2
        • msk_x
        • nperco_lchtile
        • evap_adj
        • denit_exp
        • denit_frac
        • adj_uhyd
        • cn_froz
        • dorm_hr
        • plaps
        • tlaps
        • n_fix_max
        • rsd_decay
        • urb_init_abst
        • uhyd_alpha
        • splash
        • rill
        • surq_exp
        • cov_mgt
        • co2
    • Landscape Units
      • ls_unit.def
        • id (ls_unit.def)
        • name (ls_unit.def)
        • elem_tot
        • elements
      • ls_unit.ele
        • id (ls_unit.ele)
        • name (ls_unit.ele)
        • obj_typ
        • obj_typ_no
        • bsn_frac
    • Routing Units
      • rout_unit.rtu
        • id (rout_unit.rtu)
        • name (rout_unit.rtu)
        • define
        • topo
        • field
      • rout_unit.def
        • id (rout_unit.def)
        • name (rout_unit.def)
        • elem_tot
        • elements
      • rout_unit.ele
        • id (rout_unit.ele)
        • name (rout_uni.ele)
        • obj_typ
        • obj_id
    • Hydrologic Response Units
      • hru-data.hru
        • id (hru-data.hru)
        • name (hru-data.hru)
        • topo
        • hydro
        • soil
        • lu_mgt
        • soil_plant_init
        • surf_stor
        • snow
        • field
      • hru-lte.hru
        • id (hru-lte.hru)
        • name (hru-lte.hru)
        • area
        • cn2
        • cn3_swf
        • t_conc
        • soil_dp
        • perco_co
        • slp
        • slp_len
        • et_co
        • aqu_sp_yld
        • alpha_bf
        • revap
        • rchg_dp
        • sw_init
        • aqu_init
        • aqu_sh_flo
        • aqu_dp_flo
        • snow_h2o
        • lat
        • soil_text
        • trop_flag
        • grow_start
        • grow_end
        • plnt_typ
        • stress
        • pet_flag
        • irr_flag
        • irr_src
        • t_drain
        • usle_k
        • usle_c
        • usle_p
        • usle_ls
    • Hydrology
      • topography.hyd
        • name (topography.hyd)
        • slp
        • slp_len
        • lat_len
        • depos
      • hydrology.hyd
        • name (hydrology.hyd)
        • lat_time
        • lat_sed
        • can_max
        • esco
        • epco
        • orgn_enrich
        • orgp_enrich
        • cn3_swf
        • bio_mix
        • perco
        • lat_orgn
        • lat_orgp
        • pet_co
        • latq_co
      • snow.sno
        • name (snow.sno)
        • fall_tmp
        • melt_tmp
        • melt_max
        • melt_min
        • tmp_lag
        • snow_h2o
        • cov50
        • snow_init
      • field.fld
        • name (field.fld)
        • len
        • wd
    • Soils
      • soils.sol
        • name (soils.sol)
        • hyd_grp
        • dp_tot
        • anion_excl
        • perc_crk
        • texture
        • bd
        • awc
        • soil_k
        • carbon
        • clay
        • silt
        • sand
        • rock
        • alb
        • usle_k
        • ec
        • caco3
        • ph
      • soil_plant.ini
        • name (soil_plant.ini)
        • sw_frac
        • nutrients
        • pest
        • salt
      • nutrients.sol
        • name (nutrients.sol)
        • exp_co
        • lab_p
        • nitrate
        • fr_hum_act
        • hum_c_n
        • hum_c_p
    • Landuse and Management
      • landuse.lum
        • name (landuse.lum)
        • cal_grp
        • plnt_com
        • mgt
        • cn2
        • cons_prac
        • urban
        • urb_ro
        • ov_mann
        • tile
        • sep
        • vfs
        • grww
        • bmp
      • plant.ini
        • name (plant.ini)
        • plnt_cnt
        • rot_yr_ini
        • plnt_name
        • lc_status
        • lai_init
        • bm_init
        • phu_init
        • plnt_pop
        • yrs_init
        • rsd_init
      • management.sch
        • name (management.sch)
        • op_typ
        • op_data1
        • op_data2
        • op_data3
      • cntable.lum
        • name (cntable.lum)
        • cn_a
        • cn_b
        • cn_c
        • cn_d
      • cons_practice.lum
        • name (cons_practice.lum)
        • usle_p
        • slp_len_max
      • ovn_table.lum
        • name (ovn_table.lum)
        • ovn_mean
        • ovn_min
        • ovn_max
    • Decision Tables
      • "name".dtl
        • var
        • obj
        • obj_num
        • lim_var
        • lim_op
        • lim_const
        • alt
        • act_typ
        • obj
        • obj_num
        • name
        • option
        • const
        • const2
        • fp
        • outcome
    • Management Practices
      • harv.ops
        • name (harv.ops)
        • harv_typ
        • harv_idx
        • harv_eff
        • harv_bm_min
      • graze.ops
        • name (graze.ops)
        • fertname
        • bm_eat
        • bm_tramp
        • man_amt
        • grz_bm_min
      • irr.ops
        • name (irr.ops)
        • irr_eff
        • surq_rto
        • irr_amt
        • irr_dep
        • irr_salt
        • irr_no3n
        • irr_po4
      • chem_app.ops
        • name (chem_app.ops)
        • app_eff
        • inject_dp
        • surf_frac
      • fire.ops
        • name (fire.ops)
        • chg_cn2
        • frac_burn
      • sweep.ops
        • name (sweep.ops)
        • swp_eff
        • frac_curb
    • Structural Practices
      • tiledrain.str
        • name (tiledrain.str)
        • dp
        • t_fc
        • lag
        • rad
        • dist
        • drain
        • pump
        • lat_ksat
      • filterstrip.str
        • name (filterstrip.str)
        • fld_vfs
        • con_vfs
        • cha_q
      • grassedww.str
        • name (grassedww.str)
        • mann
        • sed_co
        • dp
        • wd
        • len
        • slp
      • bmpuser.str
        • name (bmpuser.str)
        • sed_eff
        • ptlp_eff
        • solp_eff
        • ptln_eff
        • soln_eff
        • bact_eff
      • septic.str
        • name (septic.str)
        • typ
        • yr
        • operation
        • residents
        • area
        • t_fail
        • dp_bioz
        • thk_bioz
        • cha_dist
        • sep_dens
        • bm_dens
        • bod_decay
        • bod_conv
        • fc_lin
        • fc_exp
        • fecal_decay
        • tds_conv
        • mort
        • resp
        • slough1
        • slough2
        • nit
        • denit
        • p_sorp
        • p_sorp_max
        • solp_slp
        • solp_int
    • Databases
      • plants.plt
        • name (plants.plt)
        • plnt_typ
        • gro_trig
        • nfix_co
        • days_mat
        • bm_e
        • harv_idx
        • lai_pot
        • frac_hu1
        • lai_max1
        • frac_hu2
        • lai_max2
        • hu_lai_decl
        • dlai_rate
        • can_ht_max
        • rt_dp_max
        • tmp_opt
        • tmp_base
        • frac_n_yld
        • frac_p_yld
        • frac_n_em
        • frac_n_50
        • frac_n_mat
        • frac_p_em
        • frac_p_50
        • frac_p_mat
        • harv_idx_ws
        • usle_c_min
        • stcon_max
        • vpd
        • frac_stcon
        • ru_vpd
        • co2_hi
        • bm_e_hi
        • plnt_decomp
        • lai_min
        • bm_tree_acc
        • yrs_mat
        • bm_tree_max
        • ext_co
        • leaf_tov_min
        • leaf_tov_max
        • bm_dieoff
        • rt_st_beg
        • rt_st_end
        • plnt_pop1
        • frac_lai1
        • plnt_pop2
        • frac_lai2
        • frac_sw_gro
        • aeration
        • rsd_pctcov
        • rsd_covfac
      • urban.urb
        • name (urban.urb)
        • frac_imp
        • frac_dc_imp
        • curb_den
        • urb_wash
        • dirt_max
        • t_halfmax
        • conc_totn
        • conc_totp
        • conc_no3
        • urb_cn
      • tillage.til
        • name (tillage.til)
        • mix_eff
        • mix_dp
      • fertilizer.frt
        • name (fertilizer.frt)
        • min_n
        • min_p
        • org_n
        • org_p
        • nh3_n
      • pesticide.pes
        • name (pesticide.pes)
        • soil_ads
        • frac_wash
        • hl_foliage
        • hl_soil
        • solub
        • aq_reac
        • aq_volat
        • mol_wt
        • aq_resus
        • aq_settle
        • ben_act_dep
        • ben_bury
        • ben_reac
      • septic.sep
        • name (septic.sep)
        • q_rate
        • bod
        • tss
        • nh4_n
        • no3_n
        • no2_n
        • org_n
        • min_p
        • org_p
        • fcoli
    • Aquifers
      • aquifer.aqu
        • id (aquifer.aqu)
        • name (aquifer.aqu)
        • aqu_init
        • gw_flo
        • dep_bot
        • dep_wt
        • no3_n
        • sol_p
        • carbon
        • flo_dist
        • flo_max
        • alpha_bf
        • revap
        • rchg_dp
        • spec_yld
        • hl_no3n
        • flo_min
        • revap_min
      • initial.aqu
        • name (initial.aqu)
        • org_min
        • pest
        • path
        • hmet
        • salt
    • GWFlow
    • Channels
      • channel-lte.cha
        • id (channel-lte.cha)
        • name (channel-lte.cha)
        • ini
        • hyd
        • nut
      • initial.cha
        • name (initial.cha)
        • org_min
        • pest
        • salt
      • hyd-sed-lte.cha
        • name (hyd-sed-lte.cha)
        • wd
        • dp
        • slp
        • len
        • mann
        • k
        • erod_fact
        • cov_fact
        • sinu
        • eq_slp
        • d50
        • clay
        • carbon
        • dry_bd
        • side_slp
        • bed_load
        • fps
        • fpn
        • n_conc
        • p_conc
        • p_bio
      • nutrients.cha
        • name (nutrients.cha)
        • plt_n
        • plt_p
        • alg_stl
        • ben_disp
        • ben_nh3n
        • ptln_stl
        • ptlp_stl
        • cst_stl
        • ben_cst
        • cbn_bod_co
        • air_rt
        • cbn_bod_stl
        • ben_bod
        • bact_die
        • cst_decay
        • nh3n_no2n
        • no2n_no3n
        • ptln_nh3n
        • ptlp_solp
        • q2e_lt
        • q2e_alg
        • chla_alg
        • alg_n
        • alg_p
        • alg_o2_prod
        • alg_o2_resp
        • o2_nh3n
        • o2_no2n
        • alg_grow
        • alg_resp
        • slr_act
        • lt_co
        • const_n
        • const_p
        • lt_nonalg
        • alg_shd_l
        • alg_shd_nl
        • nh3_pref
    • Reservoirs and Ponds
      • reservoir.res
        • id (reservoir.res)
        • name (reservoir.res)
        • init
        • hyd
        • rel
        • sed
        • nut
      • initial.res
        • name (initial.res)
        • org_min
        • pest
        • salt
      • hydrology.res
        • name (hydrology.res)
        • yr_op
        • mon_op
        • area_ps
        • vol_ps
        • area_es
        • vol_es
        • k
        • evap_co
        • shp_co1
        • shp_co2
      • sediment.res
        • name (sediment.res)
        • nsed
        • d50
        • carbon
        • bd
        • sed_stl
        • stl_vel
      • nutrients.res
        • name (nutrients.res)
        • mid_start
        • mid_end
        • mid_n_stl
        • n_stl
        • mid_p_stl
        • p_stl
        • chla_co
        • secchi_co
        • theta_n
        • theta_p
        • n_min_stl
        • p_min_stl
      • weir.res
        • name (weir.res)
        • linear_c
        • exp_k
        • width
        • height
    • Wetlands
      • wetland.wet
        • id (wetland.wet)
        • name (wetland.wet)
        • init
        • hyd
        • rel
        • sed
        • nut
      • hydrology.wet
        • name (hydrology.wet)
        • hru_ps
        • dp_ps
        • hru_es
        • dp_es
        • k
        • evap
        • vol_area_co
        • vol_dp_a
        • vol_dp_b
        • hru_frac
    • Nutrient initialization
      • om_water.ini
        • name (om_water.ini)
        • vol
        • sed
        • part_n
        • part_p
        • no3
        • solp
        • chl_a
        • nh3
        • no2
        • cbn_bod
        • dis_ox
        • sand
        • silt
        • clay
        • sm_ag
        • l_ag
        • gvl
        • tmp
    • Constituents
      • constituents.cs
      • pest_water.ini
      • pest_hru.ini
      • salt_water.ini
      • salt_hru.ini
    • Point Sources and Inlets
      • recall.rec
        • id (recall.rec)
        • rec_typ
        • file
      • 'filename'.rec
    • Connectivity
      • 'object'.con
        • id
        • name
        • gis_id
        • area
        • lat
        • lon
        • elev
        • 'obj'
        • wst
        • cst
        • ovfl
        • rule
        • out_tot
        • obj_typ
        • obj_id
        • hyd_typ
        • frac
      • aqu_cha.lin
      • chan_surf.lin
    • Water Allocation
      • water_allocation.wro
        • rule_typ
        • cha_ob
        • ob_typ (source)
        • limit_mon
        • ob_typ (demand)
        • withdr
        • amount
        • right
        • rcv_ob
        • rcv_num
        • srcs
        • src
        • frac
        • comp
    • Calibration
      • codes.sft
      • wb_parms.sft
      • water_balance.sft
      • cal_parms.cal
      • calibration.cal
  • SWAT+ Output Files
    • Output File Format
    • Debugging Outputs
    • Soil
    • Management
    • Flow Duration Curve
    • Water Balance
    • Nutrient Balance
    • Losses
    • Plant and Weather
    • Channel
    • Aquifer
    • Reservoir
    • Recall
    • Hydrographs
    • Routing Unit
    • Pesticides
    • Object Outputs
  • Theoretical Documentation
    • ☁️Section 1: Climate
      • 🟰Chapter 1:1 Equations: Energy
        • ☀️1:1.1 Sun-Earth Relationships
          • 1:1.1.1 Distance between Earth and Sun
          • 1:1.1.2 Solar Declination
          • 🌄1:1.1.3 Solar Noon, Sunrise, Sunset, and Daylength
        • 1:1.2 Solar Radiation
          • 👽1:1.2.1 Extraterrestrial Radiation
          • 1:1.2.2 Solar Radiation under Cloudless Skies
          • 1:1.2.3 Daily Solar Radiation
          • 1:1.2.4 Hourly Solar Radiation
          • 1:1.2.5 Daily Net Radiation
        • 1:1.3 Temperature
          • 1:1.3.1 Daily Air Temperature
          • 1:1.3.2 Hourly Air Temperature
          • 1:1.3.3 Soil Temperature
          • 1:1.3.4 Water Temperature
        • 1:1:4 Wind Speed
      • Chapter 1:2 Atmospheric Water
        • 1:2.1 Precipitation
        • 1:2.2 Maximum Half-Hour Rainfall
        • 1:2.3 Water Vapor
        • 1:2.4 Snow Cover
        • 1:2.5 Snow Melt
          • 1:2.5.1 Snow Pack Temperature
            • Snow Melt Equation
      • Chapter 1:3 Weather Generator
        • 1:3.1 Precipitation
          • 1:3.1.1 Occurrence of Wet or Dry Day
          • 1:3.1.2 Amount of Precipitation
        • 1:3.2 Maximum Half-Hour Rainfall
          • 1:3.2.1 Monthly Maximum Half-Hour Rain
          • 1:3.2.2 Daily Maximum Half-Hour Rain Value
        • 1:3.3 Distribution of Rainfall Within Day
          • 1:3.3.1 Normalized Intensity Distribution
          • 1:3.3.2 Generated Time to Peak Intensity
          • 1:3.3.3 Total Rainfall and Duration
        • 1:3.4 Solar Radiation & Temperature
          • 1:3.4.1 Daily Residuals
          • 1:3.4.2 Generated Values
          • 1:3.4.3 Adjustment for Clear/Overcast Conditions
            • Maximum Temperature
            • Solar Radiation
        • 1:3.5 Relative Humidity
          • 1:3.5.1 Mean Monthly Relative Humidity
          • 1:3.5.2 Generated Daily Value
          • 1:3.5.3 Adjustment for Clear/Overcast Conditions
        • 1:3.6 Wind Speed
      • Chapter 1:4 Climate Customization
        • 1:4.1 Elevation Bands
        • 1:4.2 Climate Change
    • Section 2: Hydrology
      • Surface Runoff
        • Runoff Volume: SCS Curve Number Procedure
        • SCS Curve Number
        • Soil Hydrologic Groups
        • Antecedent Soil Moisture Condition
        • Retention Parameter
        • Slope Adjustments
        • Runoff Volume: Green & Ampt Infiltration Method
        • Peak Runoff Rate
        • Time of Concentration
          • Overland Flow Time of Concentration
          • Channel Flow Time of Concentration
        • Runoff Coefficient
        • Rainfall Intensity
        • Modified Rational Formula
        • Surface Runoff Lag
        • Transmission Losses
      • Chapter 2:2 Evapotranspiration
        • 2:2.1 Canopy Storage
        • 2:2.2 Potential Evapotranspiration
          • 2:2.2.1 Penman-Monteith Method
            • Soil Heat Flux
            • Aerodynamic Resistance
            • Canopy Resistance
            • Combined Term
          • 2:2.2.2 Priestley-Taylor Method
          • 2:2.2.3 Hargreaves Method
          • 2:2.2.4 Reading Measured or Estimated PET
        • 2:2.3 Actual Evapotranspiration
          • 2:2.3.1 Evaporation of Intercepted Rainfall
          • 2:2.3.2 Transpiration
          • 2:2.3.3 Sublimation and Evaporation from the Soil
            • Sublimation
            • Soil Water Evaporation
          • 2:2.3.4 Evaporation from Ponded Water
      • Soil Water
        • Soil Structure
        • Percolation
        • Bypass Flow
        • Perched Water Table
        • Lateral Flow
          • Lateral Flow Lag
      • Groundwater
        • Groundwater Systems
        • Shallow Aquifer
          • Recharge
          • Partitioning of Recharge Between Shallow and Deep Aquifer
          • Groundwater/Base Flow
          • Revap
          • Pumping
          • Groundwater Height
        • Deep Aquifer
    • Section 3: Nutrients/Pesticides
      • Nitrogen
        • Nitrogen Cycle in the Soil
          • Initialization of Soil Nitrogen Levels
        • Mineralization & Decomposition / Immobilization
          • Humus Mineralization
          • Residue Decomposition & Mineralization
        • Nitrification & Ammonia Volatilization
        • Denitrification
        • Atmospheric Deposition
          • Nitrogen in Rainfall
          • Nitrogen Dry Deposition
        • Fixation
        • Upward Movement of Nitrate in Water
        • Leaching
        • Nitrate in the Shallow Aquifer
      • Phosphorus
        • Phosphorus Cycle
          • Initialization of Soil Phosphorus Levels
        • Mineralization & Decomposition / Immobilization
          • Humus Mineralization
          • Residue Decomposition & Mineralization
        • Sorption of Inorganic P
        • Leaching
        • Phosphorus in the Shallow Aquifer
      • Pesticides
        • Wash-off
        • Degradation
        • Leaching
      • Bacteria
        • Wash-off
        • Bacteria Die-off/Re-growth
        • Leaching
      • Carbon
        • Sub-model Description
        • Changes from previous version
        • Analytical Solutions
    • Section 4: Erosion
      • Sediment
        • MUSLE
          • Soil Erodibility Factor
          • Cover and Management Factor
          • Support Practice Factor
          • Topographic Factor
          • Coarse Fragment Factor
        • USLE
          • Rainfall Erodibility Index
        • Snow Cover Effects
        • Sediment Lag in Surface Runoff
        • Sediment in Lateral & Groundwater Flow
      • Nutrient Transport
        • Nitrate Movement
        • Organic N in Surface Runoff
          • Enrichment Ratio
        • Soluble Phosphorus Movement
        • Organic & Mineral P Attached to Sediment in Surface Runoff
          • Enrichment Ratio
        • Nutrient Lag in Surface Runoff and Lateral Flow
      • Pesticide Transport
        • Phase Distribution of Pesticide
        • Movement of Soluble Pesticide
        • Transport of Sorbed Pesticide
          • Enrichment Ratio
        • Pesticide Lag in Surface Runoff and Lateral Flow
      • Bacteria Transport
        • Bacteria in Surface Runoff
        • Bacteria Attached to Sediment in Surface Runoff
          • Enrichment Ratio
        • Bacteria Lag in Surface Runoff
      • Water Quality Parameters
        • Algae
        • Carbonaceous Biological Oxygen Demand
          • Enrichment Ratio
        • Dissolved Oxygen
          • Oxygen Saturation Concentration
    • Section 5: Land Cover/Plant
      • Growth Cycle
        • Heat Units
          • Heat Unit Scheduling
        • Dormancy
        • Plant Types
      • Optimal Growth
        • Potential Growth
          • Biomass Production
            • Impact of Climate on Radiation-Use Efficiency
            • Modification of Biomass Calculation for Trees
          • Canopy Cover and Height
          • Root Development
          • Maturity
        • Water Uptake by Plants
          • Impact of Low Soil Water Content
          • Actual Water Uptake
        • Nutrient Uptake by Plants
          • Nitrogen Uptake
            • Nitrogen Fixation
          • Phosphorus Uptake
        • Crop Yield
      • Actual Growth
        • Growth Constraints
          • Water Stress
          • Temperature Stress
          • Nitrogen Stress
          • Phosphorus Stress
        • Actual Growth
          • Biomass Override
        • Actual Yield
          • Harvest Index Override
          • Harvest Efficiency
    • Section 6: Management Practices
      • General Management
        • Planting/Beginning of Growing Season
        • Harvest Operation
        • Grazing Operation
        • Harvest & Kill Operation
        • Kill/End of Growing Season
        • Tillage
          • Biological Mixing
        • Fertilizer Application
        • Auto-Application of Fertilizer
        • Continuous Application of Fertilizer
        • Pesticide Application
      • Water Management
        • Irrigation
          • Manual Application of Irrigation
          • Auto-Application of Irrigation
        • Tile Drainage
        • Impounded/Depressional Areas
        • Water Transfer
        • Consumptive Water Use
        • Point Source Loadings
      • Urban Areas
        • Characteristics of Urban Areas
        • Surface Runoff from Urban Areas
        • USGS Regression Equations
        • Build Up/Wash Off
          • Street Cleaning
      • Septic Systems
        • Biozone Algorithm
          • Buildup of Live Bacterial Biomass
          • Fate and Transport of Bacterial Biomass
          • Field Capacity
          • Clogging Effect on Hydraulic Conductivity
          • Soil Moisture and Percolation
          • Nitrogen, BOD, Fecal Coliform
          • Phosphorus Removal
          • Model Assumptions
        • Integration of Biozone Algorithm
          • Simulating Active and Failing Systems
      • Filter Strips and Grassed Waterways
        • Filter Strips
          • Empirical Model Development
          • Sediment Reduction Model
          • Nutrient Reduction Models
            • Total Nitrogen
            • Nitrate Nitrogen
            • Total Phosphorus
            • Soluble Phosphorus
          • VFS SWAT+ Model Structure
        • Grassed Waterways
    • Section 7: Main Channel Processes
      • Water Routing
        • Channel Characteristics
        • Flow Rate and Velocity
        • Variable Storage Routing Method
        • Muskingum Routing Method
        • Transmission Losses
        • Evaporation Losses
        • Bank Storage
        • Channel Water Balance
      • Sediment Routing
        • Landscape Contribution to Subbasin Routing Reach
        • Sediment Routing In Stream Channels
          • Simplified Bagnold Equation (Default method)
          • Physics Based Approach for Channel Erosion
        • Channel Erodibility Factor
        • Channel Cover Factor
        • Channel Downcutting and Widening
      • In-Stream Nutrient Processes
        • Algae
          • Chlorophyll a
          • Algal Growth
            • Local Specific Growth Rate of Algae
            • Local Respiration Rate of Algae
            • Local Settling Rate of Algae
        • Nitrogen Cycle
          • Organic Nitrogen
          • Ammonium
          • Nitrite
          • Nitrate
        • Phosphorus Cycle
          • Organic Phosphorus
          • Inorganic/Soluble Phosphorus
        • Carbonaceous Biological Oxygen Demand
        • Oxygen
          • Oxygen Saturation Concentration
          • Reaeration
            • Reaeration By Fickian Diffusion
            • Reaeration By Turbulent Flow Over A Dam
      • In-Stream Pesticide Transformations
        • Pesticide In The Water
          • Solid-Liquid Partitioning
          • Degradation
          • Volatilization
          • Settling
          • Outflow
        • Pesticide In The Sediment
          • Solid-Liquid Partitioning
          • Degradation
          • Resuspension
          • Diffusion
          • Burial
        • Mass Balance
      • Bacteria Routing
        • Bacteria Decay
        • Bacteria Sediment
      • Heavy Metal Routing
    • Section 8: Water Bodies
      • Impoundment Water Routing
        • Reservoirs
          • Surface Area
          • Precipitation
          • Evaporation
          • Seepage
          • Outflow
            • Measured Daily Outflow
            • Measured Monthly Outflow
            • Average Annual Release Rate For Uncontrolled Reservoir
            • Target Release For Controlled Reservoir
        • Ponds/Wetlands
          • Surface Area
          • Precipitation
          • Inflow
          • Evaporation
          • Seepage
          • Outflow
            • Pond Outflow
            • Wetland Outflow
        • Depressions/Potholes
          • Surface Area
          • Precipitation
          • Inflow
          • Evaporation
          • Seepage
          • Outflow
            • Overflow
            • Release Operation
            • Tile Flow
      • Sediment In Water Bodies
        • Mass Balance
        • Settling
        • Sediment Outflow
      • Nutrients In Water Bodies
        • Nutrient Transformations
        • Total Balance
        • Eutrophication
          • Phosphorus/Chlorophyll a Correlations
      • Pesticides In Water Bodies
        • Pesticide In The Water
          • Solid-Liquid Partitioning
          • Degradation
          • Volatilization
          • Settling
          • Outflow
        • Pesticide In The Sediment
          • Solid-Liquid Partitioning
          • Degradation
          • Resuspension
          • Diffusion
          • Burial
        • Mass Balance
      • Bacteria In Water Bodies
        • Bacteria Decay
    • References
Powered by GitBook
On this page
Export as PDF
  1. Theoretical Documentation
  2. Section 1: Climate
  3. Chapter 1:1 Equations: Energy
  4. 1:1.3 Temperature

1:1.3.3 Soil Temperature

Last updated 3 months ago

Soil temperature will fluctuate due to seasonal and diurnal variations in temperature at the surface. Figure 1:1-2 plots air temperature and soil temperature at 5 cm and 300 cm below bare soil at College Station, Texas.

This figure illustrates several important attributes of temperature variation in the soil. First, the annual variation in soil temperature follows a sinusoidal function. Second, the fluctuation in temperature during the year (the amplitude of the sine wave) decreases with depth until, at some depth in the soil, the temperature remains constant throughout the year. Finally, the timing of maximum and minimum temperatures varies with depth. Note in the above graph that there is a three month difference between the recording of the minimum temperature at the surface (January) and the minimum temperature at 300 cm (March).

Carslaw and Jaeger (1959) developed an equation to quantify the seasonal variation in temperature:

Tsoil(z,dn)=T‾AA+Asurfexp(−z/dd)sin(ωtmpdn−z/dd)T_{soil}(z,d_n)=\overline T_{AA} +A_{surf}exp(-z/dd)sin(\omega_{tmp}d_n-z/dd)Tsoil​(z,dn​)=TAA​+Asurf​exp(−z/dd)sin(ωtmp​dn​−z/dd) 1:1.3.2

where Tsoil(z,dn)T_{soil}(z,d_n) Tsoil​(z,dn​) is the soil temperature (°C\degree C°C) at depth zzz(mm) and day of the year dnd_ndn​ , T‾AA\overline T_{AA} TAA​ is the average annual soil temperature (°C\degree C°C), AsurfA_{surf}Asurf​ is the amplitude of the surface fluctuations (°C\degree C°C), dddddd is the damping depth (mm) and ωtmp\omega_{tmp}ωtmp​ is the angular frequency. When z=0z=0z=0 (soil surface), equation 1:1.3.2 reduces to Tsoil(0,dn)=T‾AA+Asurfsin(ωtmpdn).T_{soil}(0,d_n)=\overline T_{AA} + A_{surf}sin(\omega_{tmp}d_n).Tsoil​(0,dn​)=TAA​+Asurf​sin(ωtmp​dn​). As zzz →∞\rightarrow \infty→∞, equation 1:1.3.2 becomes Tsoil(∞,dn)=T‾AAT_{soil}(\infty,d_n)=\overline T_{AA}Tsoil​(∞,dn​)=TAA​.

In order to calculate values for some of the variables in this equation, the heat capacity and thermal conductivity of the soil must be known. These are properties not commonly measured in soils and attempts at estimating values from other soil properties have not proven very effective. Consequently, an equation has been adopted in SWAT+ that calculates the temperature in the soil as a function of the previous day’s soil temperature, the average annual air temperature, the current day’s soil surface temperature, and the depth in the profile.

The equation used to calculate daily average soil temperature at the center of each layer is:

Tsoil(z,dn)=ℓTsoil(z,dn1)+[1.0−ℓ][df[T‾AAair−Tssurf]+Tssurf]T_{soil}(z,d_n)=\ell T_{soil}(z,d_n1)+[1.0-\ell] [df [\overline T_{AAair}-T_{ssurf}]+T_{ssurf}]Tsoil​(z,dn​)=ℓTsoil​(z,dn​1)+[1.0−ℓ][df[TAAair​−Tssurf​]+Tssurf​] 1:1.3.3

The depth factor is calculated using the equation:

df=zdzd+exp(−0.867−2.078zd)df=\frac{zd}{zd+exp(-0.867-2.078 zd)}df=zd+exp(−0.867−2.078zd)zd​ 1:1.3.4

where zdzdzd is the ratio of the depth at the center of the soil layer to the damping depth:

zd=zddzd=\frac{z}{dd}zd=ddz​ 1:1.3.5

where zzz is the depth at the center of the soil layer (mm) and dddddd is the damping depth (mm).

From the previous three equations (1:1.3.3, 1:1.3.4 and 1:1.3.5) one can see that at depths close to the soil surface, the soil temperature is a function of the soil surface temperature. As the depth increases, soil temperature is increasingly influenced by the average annual air temperature, until at the damping depth, the soil temperature is within 5% of T‾AAair\overline T_{AAair}TAAair​.

The damping depth, dddddd, is calculated daily and is a function of the maximum damping depth, bulk density and soil water. The maximum damping depth, ddmaxdd_{max}ddmax​, is calculated:

ddmax=1000+2500ρbρb+686exp(−5.63ρb)dd_{max} = 1000+\frac{2500\rho_b}{\rho_b+686exp(-5.63\rho_b)}ddmax​=1000+ρb​+686exp(−5.63ρb​)2500ρb​​ 1:1.3.6

where ddmaxdd_{max}ddmax​ is the maximum damping depth (mm), and bbb is the soil bulk density (Mg/m3Mg/m^3Mg/m3). The impact of soil water content on the damping depth is incorporated via a scaling factor,, that is calculated with the equation:

φ=SW(0.356−0.144ρb)ztot\varphi=\frac{SW}{(0.356-0.144\rho_b) z_{tot}}φ=(0.356−0.144ρb​)ztot​SW​ 1:1.3.7

where SWSWSW is the amount of water in the soil profile expressed as depth of water in the profile (mm H2OH_{2}OH2​O), ρb\rho_bρb​ is the soil bulk density (Mg/m3Mg/m^3Mg/m3), and ztotz_{tot}ztot​ is the depth from the soil surface to the bottom of the soil profile (mm).

The daily value for the damping depth, dddddd, is calculated:

dd=ddmaxexp[ln(500ddmax)(1−φ1+φ)2]dd=dd_{max} exp[ln(\frac{500}{dd_{max}}) (\frac{1-\varphi}{1+\varphi})^2]dd=ddmax​exp[ln(ddmax​500​)(1+φ1−φ​)2] 1:1.3.8

where ddmaxdd_{max}ddmax​ is the maximum damping depth (mm), and is the scaling factor for soil water. The soil surface temperature is a function of the previous day’s temperature, the amount of ground cover and the temperature of the surface when no cover is present. The temperature of a bare soil surface is calculated with the equation:

Tbare=T‾av+εsr(Tmx−Tmn)2T_{bare}=\overline T_{av}+\varepsilon_{sr} \frac{(T_{mx}-T_{mn})}{2}Tbare​=Tav​+εsr​2(Tmx​−Tmn​)​ 1:1.3.1.9

where TbareT_{bare}Tbare​ is the temperature of the soil surface with no cover (°C\degree C°C), T‾av\overline T_{av}Tav​ is the average temperature on the day (°C\degree C°C), TmxT_{mx}Tmx​ is the daily maximum temperature (°C\degree C°C), TmnT_{mn}Tmn​ is the daily minimum temperature (°C\degree C°C), and εsr\varepsilon_{sr}εsr​ is a radiation term. The radiation term is calculated with the equation:

εsr=Hday(1−α)−1420\varepsilon_{sr}=\frac{H_{day} (1-\alpha)-14}{20}εsr​=20Hday​(1−α)−14​ 1:1.3.10

where HdayH_{day}Hday​ is the solar radiation reaching the ground on the current day (MJm−2d−1MJ m^{-2}d^{-1}MJm−2d−1), and α\alphaα is the albedo for the day. Any cover present will significantly impact the soil surface temperature. The influence of plant canopy or snow cover on soil temperature is incorporated with a weighting factor, bcvbcvbcv, calculated as:

bcv=max{CVCV+exp(7.563−1.297X10−4∗CV),SNOSNO+exp(6.055−0.3002∗SNO)}bcv=max \{{{\frac{CV}{CV+exp(7.563-1.297X10^-4*CV)}}}, \frac{SNO}{SNO+exp(6.055-0.3002*SNO)}\}bcv=max{CV+exp(7.563−1.297X10−4∗CV)CV​,SNO+exp(6.055−0.3002∗SNO)SNO​} 1:1.3.11

where CVCVCV is the total aboveground biomass and residue present on the current day (kg ha−1^{-1}−1) and SNO is the water content of the snow cover on the current day (mm H2OH_2OH2​O ). The weighting factor, bcvbcvbcv, is 0.0 for a bare soil and approaches 1.0 as cover increases.

The equation used to calculate the soil surface temperature is:

Tssurf=bcvTsoil(1,dn−1)+(1−bcv)TbareT_{ssurf}=bcv T_{soil}(1,d_n-1)+(1-bcv) T_{bare}Tssurf​=bcvTsoil​(1,dn​−1)+(1−bcv)Tbare​ 1:1.3.12

where TssurfT_{ssurf}Tssurf​ is the soil surface temperature for the current day (°C\degree C°C), bcvbcvbcv is the weighting factor for soil cover impacts, Tsoil(1,dn−1)T_{soil}(1,d_n-1)Tsoil​(1,dn​−1) is the soil temperature of the first soil layer on the previous day (°C\degree C°C), and TbareT_{bare}Tbare​ is the temperature of the bare soil surface (°C\degree C°C). The influence of ground cover is to place more emphasis on the previous day’s temperature near the surface.

Table 1:1-7: SWAT+ input variables that pertain to soil temperature.

Definition
Source Name
Input Name
Input File

tmpmx

tmp_max_ave

tmpmn

tmp_min_ave

z

dp

bd

bd

Moist soil albedo of the top layer

alb

alb

max temp

tmpmax

min temp

tmpmin

where Tsoil(z,dn)T_{soil}(z,d_n)Tsoil​(z,dn​) is the soil temperature (°C\degree C°C) at depth zzz (mm) and day of the year dnd_ndn​, ℓ\ell ℓ is the lag coefficient (ranging from 0.0 to 1.0) that controls the influence of the previous day's temperature on the current day's temperature , Tsoil(z,dn−1)T_{soil}(z,d_n-1)Tsoil​(z,dn​−1) is the soil temperature (°C\degree C°C) in the layer from the previous day, dfdfdf is the depth factor that quantifies the influence of depth below surface on soil temperature , T‾AAair\overline T_{AAair}TAAair​ is the average annual temperature (°C\degree C°C), and TssurfT_{ssurf}Tssurf​ is the soil surface temperature on the day. SWAT+ sets the lag coefficient, ℓ,\ell,ℓ, to 0.80. The soil temperature from the previous day is known and the average annual air temperature is calculated from the long-term monthly maximum and minimum temperatures reported in the weather generator input () file. This leaves the depth factor, dfdfdf, and the soil surface temperature, TssurfT_{ssurf}Tssurf​, to be defined.

SWAT+ input variables that directly impact soil temperature calculations are listed in Table 1:1-7. There are several other variables that initialize residue and snow cover in the subbasins or HRUs (snow_init in snow.sno and rsd_init in plant.ini). The influence of these variables will be limited to the first few months of simulation. Finally, the timing of management operations in the file will affect ground cover and consequently soil temperature.

Average maximum air temperature for month ()

Average minimum air temperature for month ()

: Depth from soil surface to bottom of layer (mm)

: Moist bulk density (Mg m or g cm)

: Daily maximum temperature ()

: Daily minimum temperature ()

☁️
🟰
°C\degree C°C
°C\degree C°C
zzz
bbb
−3^{-3}−3
−3^{-3}−3
TmxT_{mx}Tmx​
°C\degree C°C
TmnT_{mn}Tmn​
°C\degree C°C
weather-wgn.cli
management.sch
weather-wgn.cli
weather-wgn.cli
soils.sol
soils.sol
soils.sol
.tmp
.tmp
Figure 1:1-2: Four-year average air and soil temperature at College Station, Texas.