Function plot produces temporal graphs of the results of the soil plant water balance model (see spwb), plant water balance model (see pwb) or the forest growth model (see growth), whereas function summary summarizes the model's output in different temporal steps (i.e. weekly, annual, ...).
# S3 method for spwb
plot(x, type="PET_Precipitation",
cohorts = NULL, bySpecies = FALSE,
dates = NULL, subdaily = FALSE,
xlim = NULL, ylim=NULL, xlab=NULL, ylab=NULL, ...)
# S3 method for pwb
plot(x, type="PlantTranspiration",
cohorts = NULL, bySpecies = FALSE,
dates = NULL, subdaily = FALSE,
xlim = NULL, ylim=NULL, xlab=NULL, ylab=NULL, ...)
# S3 method for growth
plot(x, type="PET_Precipitation",
cohorts = NULL, bySpecies = FALSE,
dates = NULL, subdaily = FALSE,
xlim = NULL, ylim=NULL, xlab=NULL, ylab=NULL, ...)
# S3 method for pwb
summary(object, freq="years", output="WaterBalance", FUN=sum, bySpecies = FALSE, ...)
# S3 method for spwb
summary(object, freq="years", output="WaterBalance", FUN=sum, bySpecies = FALSE, ...)
# S3 method for growth
summary(object, freq="years", output="WaterBalance", FUN=sum, bySpecies = FALSE, ...)An object of class spwb.
The information to be plotted (see details)
An integer, boolean or character vector to select the plant cohorts to be plotted.
Allows aggregating output by species, before calculating summaries or drawing plots (only has an effect with some values of type). Aggregation can involve a sum (as for plant lai or transpiration) or a LAI-weighted mean (as for plant stress or plant water potential).
A Date vector with a subset of dates to be plotted.
Whether subdaily results should be shown, only for simulations using transpirationMode = "Sperry" and having set subdailyResults = TRUE in the simulation control object. If subdaily = TRUE, then the valid strings for type are listed in plot.spwb_day.
Range of values for x.
Range of values for y.
x-axis label.
y-axis label.
Frequency of summary statistics (see cut.Date).
The data table to be summarized. Accepted values are "DailyBalance", "PlantStress", "PlantPsi", "PlantTranspiration", "PlantPhotosynthesis" and "SoilWaterBalance", "Temperature" and "EnergyBalance".
The function to summarize results (e.g., sum, mean, ...)
Additional parameters for function plot or summary.
The following plots are currently available for spwb (most of them also for pwb):
"PET_Precipitation": Potential evapotranspiration and Precipitation.
"PET_NetRain": Potential evapotranspiration and Net rainfall.
"Snow": Snow precipitation and snowpack dynamics.
"Export": Water exported through deep drainage and surface runoff.
"Evapotranspiration": Plant transpiration and soil evaporation.
"SoilPsi": Soil water potential.
"SoilRWC": Soil relative water content (in percent of field capacity).
"SoilTheta": Soil moisture water content (in percent volume).
"SoilVol": Soil water volumetric content (in mm).
"PlantExtraction": Water extracted by plants from each soil layer.
"HydraulicRedistribution":Water added to each soil layer coming from other soil layers, transported through the plant hydraulic network (only for transpirationMode = "Sperry").
"WTD": Water table depth.
"LAI": Expanded and dead leaf area index of the whole stand.
"PlantLAI": Plant cohort leaf area index (expanded leaves).
"SoilPlantConductance":Average instantaneous overall soil plant conductance (calculated as the derivative of the supply function).
"PlantStress": Plant cohort average daily drought stress.
"PlantPsi": Plant cohort water potential (only for transpirationMode = "Granier").
"LeafPsi": Midday leaf water potential (only for transpirationMode = "Sperry").
"StemPsi": Midday (upper) stem water potential (only for transpirationMode = "Sperry").
"RootPsi": Midday root crown water potential (only for transpirationMode = "Sperry").
"PlantTranspiration": Plant cohort transpiration.
"TranspirationPerLeaf": Plant cohort transpiration per leaf area.
"PlantGrossPhotosynthesis": Plant cohort photosynthesis.
"GrossPhotosynthesisPerLeaf": Plant cohort photosynthesis per leaf area.
"PlantNetPhotosynthesis": Plant cohort net photosynthesis (only for transpirationMode = "Sperry").
"NetPhotosynthesisPerLeaf": Plant cohort net photosynthesis per leaf area (only for transpirationMode = "Sperry").
"PlantWUE": Plant cohort daily water use efficiency (gross photosynthesis over transpiration; only for transpirationMode = "Sperry").
"PlantAbsorbedSWR": Plant cohort absorbed short wave radiation (only for transpirationMode = "Sperry").
"AbsorbedSWRPerLeaf": Plant cohort absorbed short wave radiation per leaf area (only for transpirationMode = "Sperry").
"PlantNetLWR": Plant cohort net long wave radiation (only for transpirationMode = "Sperry").
"NetLWRPerLeaf": Plant cohort net long wave radiation per leaf area (only for transpirationMode = "Sperry").
"AirTemperature": Minimum/maximum/mean daily temperatures above canopy (only for transpirationMode = "Sperry").
"CanopyTemperature": Minimum/maximum/mean daily temperatures inside canopy (only for transpirationMode = "Sperry").
"SoilTemperature": Minimum/maximum/mean daily temperatures inside the first soil layer (only for transpirationMode = "Sperry").
"CanopyEnergyBalance": Canopy energy balance components (only for transpirationMode = "Sperry").
"SoilEnergyBalance": Soil energy balance components (only for transpirationMode = "Sperry").
The following are only available for growth:
"PlantRespiration": Plant cohort respiration.
"PlantRespirationPerLeaf": Plant cohort respiration per leaf area.
"PlantRespirationPerIndividual": Respiration per individual.
"PlantCBalance": Plant cohort carbon balance.
"PlantCBalancePerLeaf": Plant cohort carbon balance per leaf area.
"PlantCBalancePerIndividual": Carbon balance per individual.
"PlantCstorageFast": Amount of fast-dynamics carbon reserves.
"PlantCstorageSlow": Amount of slow-dynamics carbon reserves.
"PlantSA": Amount of sapwood area in an individual.
"PlantSAgrowth": Amount of newly-created sapwood area.
"PlantRelativeSAgrowth": Amount of newly-created sapwood area per sapwood area.
"PlantLAIlive": Plant cohort leaf area index of live leaves.
"PlantLAIdead": Plant cohort leaf area index of dead leaves.
De C<U+00E1>ceres M, Mart<U+00ED>nez-Vilalta J, Coll L, Llorens P, Casals P, Poyatos R, Pausas JG, Brotons L. (2015) Coupling a water balance model with forest inventory data to predict drought stress: the role of forest structural changes vs. climate changes. Agricultural and Forest Meteorology 213: 77-90 (doi:10.1016/j.agrformet.2015.06.012).
# NOT RUN {
#Load example daily meteorological data
data(examplemeteo)
#Load example plot plant data
data(exampleforestMED)
#Default species parameterization
data(SpParamsMED)
#Initialize soil with default soil params (2 layers)
examplesoil = soil(defaultSoilParams(2))
#Initialize control parameters
control = defaultControl("Granier")
#Initialize input
x = forest2spwbInput(exampleforestMED,examplesoil, SpParamsMED, control)
#Call simulation function
S1<-spwb(x, examplemeteo, latitude = 41.82592, elevation = 100)
#Plot results
plot(S1)
#Monthly summary (averages) of soil status
summary(S1, freq="months",FUN=mean, output="Soil")
# }
Run the code above in your browser using DataLab