medfate (version 0.8.2)

soil: Soil initialization

Description

Initializes soil parameters and state variables for its use in simulations.

Usage

soil(SoilParams, VG_PTF = "Toth", W = as.numeric(c(1)), SWE = 0)
# S3 method for soil
print(x, model="SX", ...)

Arguments

SoilParams

A list of soil parameters (see defaultSoilParams).

VG_PTF

Pedotransfer functions to obtain parameters for the van Genuchten-Mualem equations. Either "Carsel" (Carsel & Parrish 1988) or "Toth" (Toth et al. 2015).

W

A numerical vector with the initial relative water content of each soil layer.

SWE

Initial snow water equivalent of the snow pack on the soil surface (mm).

x

An object of class soil.

model

Either 'SX' or 'VG' for Saxton or Van Genuchten pedotransfer models.

...

Additional parameters to print.

Value

An list of class soil with the following elements:

• SoilDepth: Soil depth (in mm).

• W: State variable with relative water content of each layer (in as proportion relative to FC).

• Temp: State variable with temperature (in <U+00BA>C) of each layer.

• Ksoil: Kappa parameter for infiltration.

• Gsoil: Gamma parameter for infiltration.

• dVec: Width of soil layers (in mm).

• sand: Sand percentage for each layer (in percent volume).

• clay: Clay percentage for each layer (in percent volume).

• om: Organic matter percentage for each layer (in percent volume).

• usda_Type: USDA texture type.

• VG_alpha, VG_n, VG_theta_res, VG_theta_sat: Parameters for van Genuchten's pedotransfer functions, for each layer, corresponding to the USDA texture type.

• macro: Macroporosity for each layer (estimated using Stolf et al. 2011).

• rfc: Percentage of rock fragment content for each layer.

Details

Function print prompts a description of soil characteristics and state variables (water content and temperature) according to a water retention curve (either Saxton's or Van Genuchten's). Volume at field capacity is calculated assuming a soil water potential equal to -0.033 MPa. Parameter Temp is initialized as missing for all soil layers.

References

Carsel, R.F., & Parrish, R.S. 1988. Developing joint probability distributions of soil water retention characteristics. Water Resources Research 24: 755<U+2013>769.

T<U+00F3>th, B., Weynants, M., Nemes, A., Mak<U+00F3>, A., Bilas, G., & T<U+00F3>th, G. 2015. New generation of hydraulic pedotransfer functions for Europe. European Journal of Soil Science 66: 226<U+2013>238.

Stolf, R., Thurler, A., Oliveira, O., Bacchi, S., Reichardt, K., 2011. Method to estimate soil macroporosity and microporosity based on sand content and bulk density. Rev. Bras. Ciencias do Solo 35, 447<U+2013>459.

See Also

soil_psi2thetaSX, soil_psi2thetaVG, spwb, defaultSoilParams

Examples

# NOT RUN {
# Initializes soil
s = soil(defaultSoilParams())

# Prints soil characteristics according to Saxton's water retention curve
print(s, model="SX")

# Prints soil characteristics according to Van Genuchten's water retention curve
print(s, model="VG")

# }