medfate (version 4.7.0)

soil: Soil initialization

Description

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

Usage

soil(x, VG_PTF = "Toth")
# S3 method for soil
summary(object, model = "SX", ...)

Value

Function soil returns a data frame of class soil with the following columns:

• widths: 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).

• nitrogen: Sum of total nitrogen (ammonia, organic and reduced nitrogen) for each layer (in g/kg).

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

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

• Ksat: Saturated soil conductivity for each layer (estimated using function soil_saturatedConductivitySX.

• 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.

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

• Temp: State variable with temperature (in ºC) of each layer.

Arguments

x

A data frame of soil parameters (see an example in defaultSoilParams).

VG_PTF

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

object

An object of class soil.

model

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

...

Additional parameters to summary.

Author

Miquel De Cáceres Ainsa, CREAF

Details

Function summary 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.

If available, the user can specify columns VG_alpha, VG_n, VG_theta_res, VG_theta_sat and K_sat, to override Van Genuchten parameters an saturated conductivity estimated from pedotransfer functions when calling function soil.

References

Carsel, R.F., and Parrish, R.S. 1988. Developing joint probability distributions of soil water retention characteristics. Water Resources Research 24: 755–769.

Tóth, B., Weynants, M., Nemes, A., Makó, A., Bilas, G., and Tóth, G. 2015. New generation of hydraulic pedotransfer functions for Europe. European Journal of Soil Science 66: 226–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–459.

See Also

soil_redefineLayers, soil_psi2thetaSX, soil_psi2thetaVG, spwb, defaultSoilParams

Examples

# Default parameters
df_soil <- defaultSoilParams()

# Initializes soil
s = soil(df_soil)
s

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

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

# Add columns 'VG_theta_sat' and 'VG_theta_res' with custom values
df_soil$VG_theta_sat <- 0.400 
df_soil$VG_theta_res <- 0.040 

# Reinitialize soil (should override estimations)
s2 = soil(df_soil)
s2
summary(s2, model="VG")