Function plot plots 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", bySpecies = FALSE,
xlim = NULL, ylim=NULL, xlab=NULL, ylab=NULL, ...)
# S3 method for pwb
plot(x, type="PlantTranspiration", bySpecies = FALSE,
xlim = NULL, ylim=NULL, xlab=NULL, ylab=NULL, ...)
# S3 method for growth
plot(x, type="PET_Precipitation", bySpecies = 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:
"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.
"PlantTranspirationPerLeaf": Plant cohort transpiration per leaf area.
"PlantPhotosynthesis": Plant cohort photosynthesis.
"PlantPhotosynthesisPerLeaf": Plant cohort photosynthesis per leaf area.
"PlantWUE": Plant cohort daily water use efficiency (photosynthesis over transpiration).
"PlantAbsorbedSWR": Plant cohort absorbed short wave radiation (only for transpirationMode = "Sperry").
"PlantAbsorbedSWRPerLeaf": Plant cohort absorbed short wave radiation per leaf area (only for transpirationMode = "Sperry").
"PlantAbsorbedLWR": Plant cohort absorbed long wave radiation (only for transpirationMode = "Sperry").
"PlantAbsorbedLWRPerLeaf": Plant cohort absorbed 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").
"PlantRespiration": Plant cohort respiration (only for plot.growth).
"PlantRespirationPerLeaf": Plant cohort respiration per leaf area (only for plot.growth).
"PlantRespirationPerIndividual": Respiration per individual (only for plot.growth).
"PlantCBalance": Plant cohort carbon balance (only for plot.growth).
"PlantCBalancePerLeaf": Plant cohort carbon balance per leaf area (only for plot.growth).
"PlantCBalancePerIndividual": Carbon balance per individual (only for plot.growth).
"PlantCstorageFast": Amount of fast-dynamics carbon reserves (only for plot.growth).
"PlantCstorageSlow": Amount of slow-dynamics carbon reserves (only for plot.growth).
"PlantSA": Amount of sapwood area in an individual (only for plot.growth).
"PlantSAgrowth": Amount of newly-created sapwood area (only for plot.growth).
"PlantRelativeSAgrowth": Amount of newly-created sapwood area per sapwood area (only for plot.growth).
"PlantLAIlive": Plant cohort leaf area index of live leaves (only for plot.growth).
"PlantLAIdead": Plant cohort leaf area index of dead leaves (only for plot.growth).
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).
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.
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(exampleforest)
#Default species parameterization
data(SpParamsMED)
#Initialize soil with default soil params (2 layers)
examplesoil = soil(defaultSoilParams(2))
#Initialize control parameters
control = defaultControl()
#Initialize input
x = forest2spwbInput(exampleforest,examplesoil, SpParamsMED, control)
#Call simulation function
S1<-spwb(x, examplesoil, examplemeteo, 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