SAGA-GIS Tool Library Documentation (v7.3.0)

Tool WETNESS


This tool calculates a topographic wetness index (TWI) following Montgomery & Dietrich (1994) that can be used to estimate the degree of saturation of unconsolidated, permeable materials above (more or less) impermeable bedrock. In contrast to the common TOPMODEL (Beven & Kirkby, 1979) - based TWI, this index differs in such that it considers hydraulic conductivity to be constant in a soil mantle overlying relatively impermeable bedrock. Also, it uses the sine of the slope rather than its tangens, which is more correct and significantly matters for steeper slopes that give raise to landslides. For computation, a slope (in radians) and a catchment area (in m2) grid are required. Additionally, information on groundwater recharge (m/hr), material hydraulic conductivity (m/hr), and depth to potential shear plane (m) are required that can be specified either globally or through grids. The tool produces a continuous wetness index (-) where cells with WI values > 1 (overland flow) set to 1, and optionally creates a classified WI grid rendering three saturation classes:.
0): Low moisture (WI smaller 0.1)
1): Partially wet (0.1 smaller WI smaller 1)
2): Saturation zone (WI larger 1)

References:
Beven, K.J., Kirkby, M.J. (1979) A physically-based variable contributing area model of basin hydrology. Hydrology Science Bulletin, 24, 43-69..

Montgomery D. R., Dietrich, W. E. (1994) A physically based model for the topographic control on shallow landsliding. Water Resources Research, 30, 1153-1171..

Parameters

 NameTypeIdentifierDescriptionConstraints
InputDEMGrid (input)DEMA DEM-
Min hydraulic conductivity grid (m/hr) (*)Grid (optional input)CminA grid representing minimum material hydraulic conductivity (in m/hr)-
Max hydraulic conductivity grid (m/hr) (*)Grid (optional input)CmaxA grid representing maximum material hydraulic conductivity (in m/hr)-
Min groundwater recharge grid (m/hr) (*)Grid (optional input)DminA grid representing minimum groundwater recharge (in m/hr)-
Max groundwater recharge grid (m/hr) (*)Grid (optional input)DmaxA grid representing maximum groundwater recharge (in m/hr)-
Min material depth grid (m) (*)Grid (optional input)EminA grid representing minimum depth to potential shear plane (in m)-
Max material depth grid (m) (*)Grid (optional input)EmaxA grid representing maximum depth to potential shear plane (in m)-
OutputWI valuesGrid (output)FResulting wetness index (-) grid-
WI classes (*)Grid (optional output)GClassified wetness (-) grid-
OptionsGrid systemGrid systemPARAMETERS_GRID_SYSTEM--
Min global material conductivity (m/hr)Floating pointfCminConstant value if no raster setDefault: 2.700000
Max global material conductivity (m/hr)Floating pointfCmaxConstant value if no raster setDefault: 2.700000
Min global groundwater recharge (m/hr)Floating pointfDminConstant value if no raster setDefault: 0.001000
Max global groundwater recharge (m/hr)Floating pointfDmaxConstant value if no raster setDefault: 0.001000
Min global material depth (m)Floating pointfEminConstant value if no raster setDefault: 1.000000
Max global material depth (m)Floating pointfEmaxConstant value if no raster setDefault: 1.000000
Parameter sampling runsIntegerfHNumber of sampling cyclesDefault: 1
Catchment Area CalculationChoiceMETHOD-Available Choices:
[0] Deterministic 8
[1] Rho 8
[2] Braunschweiger Reliefmodell
[3] Deterministic Infinity
[4] Multiple Flow Direction
[5] Multiple Triangular Flow Directon
Default: 4
PreprocessingBooleanPREPROC-Default: 0
(*) optional

Command-line

Usage: saga_cmd ta_slope_stability 3 [-DEM <str>] [-Cmin <str>] [-Cmax <str>] [-Dmin <str>] [-Dmax <str>] [-Emin <str>] [-Emax <str>] [-fCmin <double>] [-fCmax <double>] [-fDmin <double>] [-fDmax <double>] [-fEmin <double>] [-fEmax <double>] [-fH <num>] [-F <str>] [-G <str>] [-METHOD <str>] [-PREPROC <str>]
  -DEM:<str>     	DEM
	Grid (input)
  -Cmin:<str>    	Min hydraulic conductivity grid (m/hr) 
	Grid (optional input)
  -Cmax:<str>    	Max hydraulic conductivity grid (m/hr) 
	Grid (optional input)
  -Dmin:<str>    	Min groundwater recharge grid (m/hr) 
	Grid (optional input)
  -Dmax:<str>    	Max groundwater recharge grid (m/hr) 
	Grid (optional input)
  -Emin:<str>    	Min material depth grid (m)
	Grid (optional input)
  -Emax:<str>    	Max material depth grid (m)
	Grid (optional input)
  -fCmin:<double>	Min global material conductivity (m/hr)
	Floating point
	Default: 2.700000
  -fCmax:<double>	Max global material conductivity (m/hr)
	Floating point
	Default: 2.700000
  -fDmin:<double>	Min global groundwater recharge (m/hr)
	Floating point
	Default: 0.001000
  -fDmax:<double>	Max global groundwater recharge (m/hr)
	Floating point
	Default: 0.001000
  -fEmin:<double>	Min global material depth (m)
	Floating point
	Default: 1.000000
  -fEmax:<double>	Max global material depth (m)
	Floating point
	Default: 1.000000
  -fH:<num>      	Parameter sampling runs
	Integer
	Default: 1
  -F:<str>       	WI values
	Grid (output)
  -G:<str>       	WI classes
	Grid (optional output)
  -METHOD:<str>  	Catchment Area Calculation
	Choice
	Available Choices:
	[0] Deterministic 8
	[1] Rho 8
	[2] Braunschweiger Reliefmodell
	[3] Deterministic Infinity
	[4] Multiple Flow Direction
	[5] Multiple Triangular Flow Directon
	Default: 4
  -PREPROC:<str> 	Preprocessing
	Boolean
	Default: 0