GEOtop_VisSoilWaterRet_gg: Visualisation Soil Water Retention Curve with ggplot

In JBrenn/AnalyseGEOtop: Analyse GEOtop simulations

Description

Visulisation of the Soil Water Retention Curve van Genuchten Model (1980). For soil water calculation from soil water pressure head (hPa or cm water column) the package soilwater is used. Density functions for provided soil depths are visualized as side plots.

Usage

Geotop_VisSoilWaterRet_gg(alpha, n, theta_sat, theta_res, accurate=10,
                        observed=NULL, add_ref_curves=TRUE)

Arguments

alpha	alpha van Genuchten parameter [1/m]; inverse of a length - scale parameter
n	N van Genuchten parameter [-]; shape parameter
theta_sat	saturated soil water content [m³/m³]
theta_res	residual soil water content [m³/m³]
accurate	default = 10: water potential is given as seq(1,10^7,10), for accurate plotting set argument 1
observed	numeric data.frame, providing observed soil moisture data [m³/m³] (first column, named SWC) and water potential [hPa] (second column, named SWP) and soil depth [cm] (third coumn, named depth)
add_ref_curves	boolean, default = TRUE: reference soil water retention curves are added for clay, loam and sand; van Genuchten parameterisation according to ROSETTA manual
colors	colors for water retention curves and density functions, default = NULL: grey color range is used

Details

If multiple vanGenuchten models should be visualised, simply provide alpha, n, theta_sat, and theta_res values as vectors.

Value

gg table object

Author(s)

Johannes Brenner, [e-mail]

References

Van Genuchten, M. (1980). A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal, 44, 892–898.

Examples

## Not run:
library(AnalyseGeotop)

data("Kaltern_SoilWatRet")

colors <- c("red","#E69F00","blue","darkgreen")

gg <- Geotop_VisSoilWaterRet_gg(alpha = parameter$alpha, n = parameter$n, 
                                theta_sat = parameter$thetaS, 
                                theta_res = rep(.05,length(parameter$alpha)), 
                                accurate = 10, 
                                add_ref_curves = T, observed = observed, 
                                colors = colors)
## End(**Not run**)                                

JBrenn/AnalyseGEOtop documentation built on May 7, 2019, 6:48 a.m.