Skip to contents

Compute all hydraulic data using recalibrate new Morel equations

Source: R/TnetHydraulic.R
TnetHydraulic_MorelNewRec.Rd

Calculate rivers width (\(Bm\)), depth (\(H\)) and water travel time (\(TPS\)) with Discharge as input data, using Morel formulas exactly like TnetHydraulic_MorelNew(). The difference is that some river segments cant have the mean term recalibrate using observed width at a specific discharge quantile.

Usage

TnetHydraulic_MorelNewRec(
  path_data,
  shapefile,
  export_files = c("Bm", "H", "TPS"),
  shapefile_largeur,
  quantile_calc
)

Arguments

path_data

Path to the folder containing discharge data created by Tnet*_computeQ

shapefile

Path to the shapefile with all info on segments. Columns needed are detailed in the Shapefile columns section.

export_files

Vector of files that will be exported. It can contain "Bm", "H", "CV" and "TPS", but only "Bm", "H" and "TPS" are needed for T-NET

shapefile_largeur

Path to a shapefile with all segments to recalibrate and the width. Columns needed are detailed in the Shapefile columns section.

quantile_calc

Quantile of discharge to use for observed river width in shapefile_largeur.

Value

All hydraulic data as NetCDF files in the folder path_data

Equations used

River segment width (\(Bm\))

On normal segments

$$Bm = ad \cdot (Qmean)^{bd} \cdot \left(\frac{Qaval}{Qmean}\right)^{b}$$

with \(ad = 7.482\), \(bd = 0.477\), \(b = 0.148\)

On recalibrate segments

$$Bm = largeur \cdot \left(\frac{Qaval}{Qquantile}\right)^{b}$$

with \(b = 0.148\)

where:

  • \(largeur\) is the observed width of the river segment in meters. Provided in shapefile_largeur.

  • \(Qquantile\) is the flow quantile at which the \(largeur\) is observed. Provided in quantile_calc.

River segment depth (\(H\))

$$H = cd \cdot (Qmean)^{fd} \cdot \left(\frac{Qaval}{Qmean}\right)^{f}$$

with \(cd = 0.340\), \(fd = 0.259\), \(f = 0.292\)

Water speed in the river segment (\(CV\))

$$CV = \frac{Qaval}{H \cdot Bm}$$

Water travel time in river segment (\(TPS\))

$$TPS = \frac{Longueur\_m}{CV} \div 3600$$

where:

  • \(Longueur\_m\) is the length of the river segment in meters.

  • \(Qaval\) is the flow rate at the segment's outlet.

  • \(Qmean\) is the mean flow rate on the entire period.



Shapefiles columns

In shapefile

\(gid\_new\)ID of the Topage segment
\(Longueur\_m\)Length of the river segment (in meter)

In shapefile_largeur

\(gid\_new\)ID of the Topage segment
\(largeur\)Width to recalibrate the river segment (in meter)

References

Maxime Morel, Doug J. Booker, Frédéric Gob, Nicolas Lamouroux, Intercontinental predictions of river hydraulic geometry from catchment physical characteristics, Journal of Hydrology, https://doi.org/10.1016/j.jhydrol.2019.124292.

See also

TNET_initializeSim(), TNET_computeQ(), TnetHydraulic_MorelNew()

Examples


TnetHydraulic_MorelNewRec(path_data = path/to/data,
                      shapefile = path/to/shapefile.shp,
                      calcul_Bm = 'EKrecal',
                      shapefile_largeur = path/to/shapefile_width.shp,
                      quantile_calc = 0.99)
#> Error in TnetHydraulic_MorelNewRec(path_data = path/to/data, shapefile = path/to/shapefile.shp,     calcul_Bm = "EKrecal", shapefile_largeur = path/to/shapefile_width.shp,     quantile_calc = 0.99): unused argument (calcul_Bm = "EKrecal")


############################################################
## Using this function with TNET_initializeSim() function ##
############################################################
infoSimu <- TNET_initializeSim(...)
#> Error: '...' used in an incorrect context

#Do discharge computation needed using a Tnet*_computeQ() function

TnetHydraulic_MorelNewRec(path_data = infoSimu$hydroPath, 
                      shapefile = infoSimu$TOPAGE_shape,
                      shapefile_largeur = infoSimu$Largeur_shape,
                      quantile_calc =  infoSimu$Bm_quantile)
#> Error: object 'infoSimu' not found