Skip to content

r.subdayprecip.design

Computes subday design precipitation totals.

r.subdayprecip.design map=name return_period=name [,name,...] rainlength=integer [area_size=float] [--verbose] [--quiet] [--qq] [--ui]

Example:

r.subdayprecip.design map=name return_period=name rainlength=integer

grass.script.run_command("r.subdayprecip.design", map, return_period, rainlength, area_size=20, verbose=False, quiet=False, superquiet=False)

Example:

gs.run_command("r.subdayprecip.design", map="name", return_period="name", rainlength=integer)

Parameters

map=name [required]
    Vector map of location under analysis
    Or data source for direct OGR access
return_period=name [,name,...] [required]
    Rainfall raster maps of required return period
rainlength=integer [required]
    Design rainfall length in minutes
    Allowed values: 0-1439
area_size=float
    Maximum area size to be processed (in km2, -1 for no limit)
    Default: 20
--help
    Print usage summary
--verbose
    Verbose module output
--quiet
    Quiet module output
--qq
    Very quiet module output
--ui
    Force launching GUI dialog

map : str, required
    Vector map of location under analysis
    Or data source for direct OGR access
    Used as: input, vector, name
return_period : str | list[str], required
    Rainfall raster maps of required return period
    Used as: input, raster, name
rainlength : int, required
    Design rainfall length in minutes
    Allowed values: 0-1439
area_size : float, optional
    Maximum area size to be processed (in km2, -1 for no limit)
    Default: 20
verbose: bool, optional
    Verbose module output
    Default: False
quiet: bool, optional
    Quiet module output
    Default: False
superquiet: bool, optional
    Very quiet module output
    Default: False

DESCRIPTION

r.subdayprecip.design computes subday design precipitation totals based on Hradek's method of reduction of daily maximums to chosen duration.

The tool uses methods of zonal statistics to compute average values of 24 hours precipitation amounts of specified return period for a provided area or a spot. Rasters of daily maxima were derived from statistics published by Samaj et al. in 1985, which were based on precipitation series observed in 1901-1980. Calculated average value (in millimeters) is then reduced to the chosen length of design rain event.

NOTES

Subday design precipitation series are important for hydrological modelling and soil erosion problems in a small catchment scale when designing common measures for promoting water retention, landscape drainage systems, flood mitigation measures etc.

First automatization has been implemented by well-known method which is based on reduction of 24 hours design precipitation to shorter time. GIS can used for spatial supervised classification of point values of specified return periods (2, 10, 20, 50 and 100 years).

Basins (in orange) with orthophoto on background
Figure: Basins (in orange) with orthophoto on background.

Return periods (2, 10, 20, 50 years) in the area of the Czech Republic
Figure: Return periods (2, 10, 20, 50 years) in the area of the Czech Republic.

IV.order basins colored by mean H_N2T60 value
Figure: IV.order basins colored by mean H_N2T60 value (in millimeters)

EXAMPLE

r.subdayprecip.design map=basin return_period=N2,N5,N10,N20 rainlength=60

REFERENCES

Acknowledgement

This work has been supported by the research project QJ1520265 - "Variability of Short-term Precipitation and Runoff in Small Czech Drainage Basins and its Influence on Water Resources Management".

SEE ALSO

v.rast.stats, v.what.rast

AUTHOR

Martin Landa, GeoForAll (OSGeoREL) Lab, Czech Technical University in Prague, Czech Republic
The module is inspired by Python script developed for Esri ArcGIS platform by M. Tomasu in 2013.

SOURCE CODE

Available at: r.subdayprecip.design source code (history)
Latest change: Friday Feb 21 23:33:10 2025 in commit 7c35290