SAGA-GIS Tool Library Documentation (v8.1.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, inputDEMA DEM-
Min hydraulic conductivity grid (m/hr) (*)Grid, input, optionalCminA grid representing minimum material hydraulic conductivity (in m/hr)-
Max hydraulic conductivity grid (m/hr) (*)Grid, input, optionalCmaxA grid representing maximum material hydraulic conductivity (in m/hr)-
Min groundwater recharge grid (m/hr) (*)Grid, input, optionalDminA grid representing minimum groundwater recharge (in m/hr)-
Max groundwater recharge grid (m/hr) (*)Grid, input, optionalDmaxA grid representing maximum groundwater recharge (in m/hr)-
Min material depth grid (m) (*)Grid, input, optionalEminA grid representing minimum depth to potential shear plane (in m)-
Max material depth grid (m) (*)Grid, input, optionalEmaxA grid representing maximum depth to potential shear plane (in m)-
OutputWI valuesGrid, outputFResulting wetness index (-) grid-
WI classes (*)Grid, output, optionalGClassified 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, input, optional
  -Cmax:<str>    	Max hydraulic conductivity grid (m/hr) 
	Grid, input, optional
  -Dmin:<str>    	Min groundwater recharge grid (m/hr) 
	Grid, input, optional
  -Dmax:<str>    	Max groundwater recharge grid (m/hr) 
	Grid, input, optional
  -Emin:<str>    	Min material depth grid (m)
	Grid, input, optional
  -Emax:<str>    	Max material depth grid (m)
	Grid, input, optional
  -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, output, optional
  -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