HydroCode (version 1.0.3)

rhineHYBASlev08: HydroBASINS Level 08 Data of the Rhine Catchment

Description

This spatial dataset is a subset of the data layer hybas_eu_lev08_v1c from the HydroBASINS v1.c database (Lehner and Grill, 2013). The original data format was SHP (ESRI Shape File). The subset consists of subbasins within the Rhine basin, i.e. MAIN_BAS = 2080023010.

Format

An Object of class SpatialPolygonsDataFrame with 289 features. The attribute table comprises 13 fields:

HYBAS_ID

factor, unique basin identifier

NEXT_DOWN

factor, Hybas_id of the next downstream polygon

NEXT_SINK

factor, Hybas_id of the next downstream sink

MAIN_BAS

factor, Hybas_id of the most downstream sink, i.e. the outlet of the main river basin

DIST_SINK

numeric, Distance from polygon outlet to the next downstream sink along the river network, in kilometers

DIST_MAIN

numeric, Distance from polygon outlet to the most downstream sink, i.e. the outlet of the main river basin along the river network, in kilometers

SUB_AREA

numeric, Area of the individual polygon (i.e. sub-basin), in square kilometers

UP_AREA

numeric, Total upstream area, in square kilometers, calculated from the headwaters to the polygon location (including the polygon)

PFAF_IF

int, the Pfafstetter code

ENDO

int, Indicator for endorheic (inland) basins without surface flow connection to the ocean: 0 = not part of an endorheic basin; 1 = part of an endorheic basin; 2 = sink (i.e. most downstream polygon) of an endorheic basin.

COAST

int, indicator for lumped coastal basins: 0 = no; 1 = yes

ORDER

int, Indicator of river order (classical ordering system): order 1 represents the main stem river from sink to source; order 2 represents all tributaries that flow into a 1 st order river; order 3 represents all tributaries that flow into a 2 nd order river; etc.; order 0 is used for conglomerates of small coastal watersheds

SORT

int, Indicator showing the record number (sequence) in which the original polygons are stored in the shapefile (i.e. counting upwards from 1 in the original shapefile).

The PROJ.4 string is +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84+towgs84=0,0,0

References

B. Lehner and G. Grill, 2013, Global river hydrography and network routing: baseline data and new approaches to study the world’s large river systems, Hydrological Processes, 27, 2171--2186, tools:::Rd_expr_doi("10.1002/hyp.9740").

Examples

if (FALSE) {
## Catchment of gauge Mainz / Rhein
## HYBAS uses a modified version of Pfafstetter
op.old <- options()
options("hydrocode.alloddorzero" = TRUE)
ok <- is.downstream(23263070, rhineHYBASlev08$PFAF_ID)
MzBas <- subset(rhineHYBASlev08, subset = ok)
sp::plot(rhineHYBASlev08)
sp::plot(MzBas, add = TRUE, col = "red")

## get catchment size
## AEO : 98206 km2 according to website
sum(MzBas$SUB_AREA)

## reset previous settings
options(op.old)
}