Function hydrology_soilWaterBalance estimates water balance of soil layers given water inputs/outputs, including the simulation of water movement within the soil.
hydrology_soilWaterBalance(
soil,
soilFunctions,
rainfallInput,
rainfallIntensity,
snowmelt,
sourceSink,
runon = 0,
lateralFlows = NULL,
waterTableDepth = NA_real_,
infiltrationMode = "GreenAmpt1911",
infiltrationCorrection = 5,
soilDomains = "buckets",
nsteps = 24L,
max_nsubsteps = 3600L,
modifySoil = TRUE
)Returns a named vector with different elements, depending on soilDomains.
If soilDomains == "buckets":
Snowmelt: Snowmelt input (mm).
Source/sinks: Sum of source/sink input across layers (mm).
Infiltration: Water infiltrated into the soil (mm).
InfiltrationExcess: Excess infiltration in the topmost layer (mm) leading to an increase in runoff.
SaturationExcess: Excess saturation in the topmost layer (mm) leading to an increase in runoff.
Runoff: Surface runoff generated by saturation excess or infiltration excess (mm).
DeepDrainage: Water draining from the bottom layer (mm). This quantity is corrected to close the water balance.
CapillarityRise: Water entering the soil via capillarity rise (mm) from the water table, if waterTableDepth is supplied.
If soilDomains == "single" the named vector contains the following additional elements:
Correction: Amount of water (mm) added to deep drainage to correct the water balance.
VolumeChange: Change in soil water volume (mm).
Substep: Time step of the moisture solving (seconds).
If soilDomains == "dual" the named vector contains the following additional elements:
Lateral flows: Sum of water circulating between micropores and macropores, positive when filling micropores (mm).
InfiltrationMatrix: Water infiltrated into the soil matrix (mm).
InfiltrationMacropores: Water infiltrated into the soil macropore domain (mm).
InfiltrationExcessMatrix/InfiltrationExcessMacropores: Excess infiltration in the topmost layer (mm) leading to an increase in runoff.
SaturationExcessMatrix/SaturationExcessMacropores: Excess saturation in the topmost layer (mm) leading to an increase in runoff.
DrainageMatrix: Water draining from the bottom layer of the matrix domain (mm). This quantity is corrected to close water balance in the micropore domain.
DrainageMacropores: Water draining from the bottom layer of the macropore domain (mm). This quantity is corrected to close the water balance in the macropore domain.
CorrectionMatrix: Amount of water (mm) added to deep drainage of soil matrix to correct the water balance.
CorrectionMacropores: Amount of water (mm) added to deep drainage of macropores to correct the water balance.
MatrixVolumeChange: Change in soil water volume in the soil matrix domain (mm).
MacroporeVolumeChange: Change in soil water volume in the macropore domain (mm).
Object of class soil.
Soil water retention curve and conductivity functions, either 'SX' (for Saxton) or 'VG' (for Van Genuchten).
Amount of water from rainfall event (after excluding interception), in mm.
Rainfall intensity, in mm/h.
Amount of water originated from snow melt, in mm.
Local source/sink term for each soil layer (from soil evaporation or plant transpiration/redistribution) as mm/day.
Surface water amount running on the target area from upslope (in mm).
Lateral source/sink terms for each soil layer (interflow/to from adjacent locations) as mm/day.
Water table depth (in mm). When not missing, capillarity rise will be allowed if lower than total soil depth.
Infiltration model, either "GreenAmpt1911" or "Boughton1989"
Correction for saturated conductivity, to account for increased infiltration due to macropore presence
Either "buckets" (multi-bucket domain), "single" (for single-domain Richards) or "dual" (for dual-permeability model).
Number of time steps per day
Maximum number of substeps per time step
Boolean flag to indicate that the input soil object should be modified during the simulation.
Miquel De Cáceres Ainsa, CREAF
María González Sanchís, UPV-CTFC
The multi-bucket model adds/substracts water to each layer and if content is above field capacity the excess percolates to the layer below. If there is still an excess for the bottom layer, the model will progressively fill upper layers (generating saturation excess if the first layer becomes saturated). Every day the some layers are over field capacity, the model simulates deep drainage.
The single-domain model simulates water flows by solving Richards's equation using the predictor-corrector method, as described in Bonan et al. (2019).
The dual-permeability model is an implementation of the model MACRO 5.0 (Jarvis et al. 1991; Larsbo et al. 2005).
Both the multi-bucket and the single-domain model can apply a correction to the infiltration rate to account for macroporosity in infiltration. In the dual-permeability model extra infiltration through macropores is simulated explicitly.
Bonan, G. (2019). Climate change and terrestrial ecosystem modeling. Cambridge University Press, Cambridge, UK.
Jarvis, N.J., Jansson, P‐E., Dik, P.E. & Messing, I. (1991). Modelling water and solute transport in macroporous soil. I. Model description and sensitivity analysis. Journal of Soil Science, 42, 59–70.
Larsbo, M., Roulier, S., Stenemo, F., Kasteel, R. & Jarvis, N. (2005). An Improved Dual‐Permeability Model of Water Flow and Solute Transport in the Vadose Zone. Vadose Zone Journal, 4, 398–406.
spwb, hydrology_waterInputs, hydrology_infiltration
# Define soil parameters
spar <- defaultSoilParams(4)
# Initializes soil hydraulic parameters
examplesoil <- soil(spar)
# Water balance in a multi-bucket model
hydrology_soilWaterBalance(examplesoil, "VG", 10, 5, 0, c(-1,-1,-1,-1),
soilDomains = "buckets", modifySoil = FALSE)
# Water balance in a single-domain model (Richards equation)
hydrology_soilWaterBalance(examplesoil, "VG", 10, 5, 0, c(-1,-1,-1,-1),
soilDomains = "single", modifySoil = FALSE)
# Water balance in a dual-permeability model (MACRO)
hydrology_soilWaterBalance(examplesoil, "VG", 10, 5, 0, c(-1,-1,-1,-1),
soilDomains = "dual", modifySoil = FALSE)
Run the code above in your browser using DataLab