.TH r.uslek 1grass "" "GRASS 8.3.1" "GRASS GIS User's Manual"
.SH NAME
\fI\fBr.uslek\fR\fR  \- Computes USLE Soil Erodibility Factor (K).
.SH KEYWORDS
raster, hydrology, soil, erosion
.SH SYNOPSIS
\fBr.uslek\fR
.br
\fBr.uslek \-\-help\fR
.br
\fBr.uslek\fR \fBpsand\fR=\fIname\fR \fBpclay\fR=\fIname\fR \fBpsilt\fR=\fIname\fR \fBpomat\fR=\fIname\fR \fBoutput\fR=\fIname\fR  [\-\-\fBoverwrite\fR]  [\-\-\fBhelp\fR]  [\-\-\fBverbose\fR]  [\-\-\fBquiet\fR]  [\-\-\fBui\fR]
.SS Flags:
.IP "\fB\-\-overwrite\fR" 4m
.br
Allow output files to overwrite existing files
.IP "\fB\-\-help\fR" 4m
.br
Print usage summary
.IP "\fB\-\-verbose\fR" 4m
.br
Verbose module output
.IP "\fB\-\-quiet\fR" 4m
.br
Quiet module output
.IP "\fB\-\-ui\fR" 4m
.br
Force launching GUI dialog
.SS Parameters:
.IP "\fBpsand\fR=\fIname\fR \fB[required]\fR" 4m
.br
Name of soil sand fraction raster map [0.0\-1.0]
.IP "\fBpclay\fR=\fIname\fR \fB[required]\fR" 4m
.br
Name of soil clay fraction raster map [0.0\-1.0]
.IP "\fBpsilt\fR=\fIname\fR \fB[required]\fR" 4m
.br
Name of soil silt fraction raster map [0.0\-1.0]
.IP "\fBpomat\fR=\fIname\fR \fB[required]\fR" 4m
.br
Name of soil organic matter raster map [0.0\-1.0]
.IP "\fBoutput\fR=\fIname\fR \fB[required]\fR" 4m
.br
Name for output USLE K factor raster map [t.ha.hr/ha.MJ.mm]
.SH DESCRIPTION
\fIr.uslek\fR calculates the USLE K factor, that is the Soil Erodibility
Factor. It takes input of soil texture classes (sand, clay, silt) and
organic matter, all in range of [0.0\-1.0]. The FAO World Soil CD
documentation was used to produce the conversion system between soil
textures and soil classes. The soil classes are in number of 12 and
apparently come from a USDA publication of 1951 (p.209). Once the soil
classes have been identified (by vector cross\-products tests), a general
conversion table was applied for transforming soil classes into K factor.
.SH NOTES
\fIr.watershed\fR provides for USLE L, S, LS factors.
.SH EXAMPLE
Calculate the USLE K factor (synthetic values, North Carolina sample dataset):
.br
.nf
\fC
# NC data: simulation
g.region rural_1m \-p
r.mapcalc \(dqsand  = 0.11\(dq
r.mapcalc \(dqclay  = 0.47\(dq
r.mapcalc \(dqsilt  = 0.42\(dq
r.mapcalc \(dqorgan = 0.0136\(dq
r.uslek psand=sand pclay=clay psilt=silt pomat=organ output=usle_k
r.univar usle_k
# Result: uslek = 0.25
\fR
.fi
.PP
Example with small variations (condition: d_sand + d_clay + d_silt => 1.0)
.br
.nf
\fC
# NC data: simulation
r.mapcalc \-s \(dqsand_r = rand(0.09,0.13)\(dq
r.mapcalc \-s \(dqclay_r = rand(0.45,0.49)\(dq
r.mapcalc \-s \(dqsilt_r = rand(0.40,0.44)\(dq
r.mapcalc \-s \(dqorgan  = rand(0.01,0.70)\(dq
r.mapcalc \(dqtemp_sum = sand_r + clay_r + silt_r\(dq
r.mapcalc \(dqsand = sand_r / temp_sum\(dq
r.mapcalc \(dqclay = clay_r / temp_sum\(dq
r.mapcalc \(dqsilt = silt_r / temp_sum\(dq
r.mapcalc \(dqtest = sand + clay + silt\(dq
r.info \-r test
g.remove \-f rast name=sand_r,clay_r,silt_r,temp_sum,test
r.uslek psand=sand pclay=clay psilt=silt pomat=organ output=usle_k
r.univar usle_k
# Result: uslek =  variable values
\fR
.fi
.SH SEE ALSO
\fI
r.usler,
r.watershed
\fR
.SH AUTHOR
Yann Chemin, SIC\-ISDC, Turkmenistan
.SH SOURCE CODE
.PP
Available at:
r.uslek source code
(history)
.PP
Accessed: Sunday Dec 17 17:22:07 2023
.PP
Main index |
Raster index |
Topics index |
Keywords index |
Graphical index |
Full index
.PP
© 2003\-2023
GRASS Development Team,
GRASS GIS 8.3.1 Reference Manual