Function plot produces temporal graphs of the results of the soil plant water balance model (see spwb), plant water balance model (see pwb), the forest growth model (see growth) or the forest dynamics model (see fordyn).
# 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,
summary.freq = 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,
summary.freq = 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,
summary.freq = NULL, ...)
# S3 method for fordyn
plot(x, type="StandBasalArea",
ylim=NULL, xlab=NULL, ylab=NULL, ...)An object of class spwb, pwb, growth or fordyn.
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 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).
Additional parameters for function plot (not used).
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:
"GrossPhotosynthesis": Gross photosynthesis rate per dry weight.
"MaintenanceRespiration": Maintenance respiration cost per dry weight.
"CarbonBalance": Carbon balance per dry weight.
"SugarLeaf": Sugar concentration in leaves.
"StarchLeaf": Starch concentration in leaves.
"SugarSapwood": Sugar concentration in sapwood.
"StarchSapwood": Starch concentration in sapwood.
"SugarTransport": Phloem sugar transport rate.
"RootExudation": Root exudation rate per dry weight.
"LeafPI0": Leaf osmotic concentration at full turgor.
"StemPI0": Sapwood osmotic concentration at full turgor.
"SapwoodArea": Sapwood area per individual.
"LeafArea": Leaf area per individual.
"FineRootArea": Fine root area per individual.
"SapwoodBiomass": Sapwood dry biomass per individual.
"LeafBiomass": Leaf dry biomass per individual.
"FineRootBiomass": Fine root dry biomass per individual.
"LabileBiomass": Labile C biomass per individual.
"TotalLivingBiomass": Total dry biomass per individual.
"SAgrowth": Sapwood area growth rate.
"LAgrowth": Leaf area growth rate.
"FRAgrowth": Fine root area growth rate.
"HuberValue": Ratio of leaf area to sapwood area.
"RootAreaLeafArea": Ratio of fine root area to leaf area.
The following are only available for fordyn:
"StandBasalArea":Stand basal area of living trees.
"StandLAI":Stand leaf area index.
"StandDensity":Stand density of living trees.
"SpeciesBasalArea":Basal area of living trees by species.
"SpeciesLAI":Leaf area index by species.
"SpeciesDensity":Density of living trees by species.
"CohortBasalArea":Basal area of living trees by plant cohort.
"CohortLAI":Leaf area index by plant cohort.
"CohortDensity":Density of living trees by plant cohort.
# 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)
# }
Run the code above in your browser using DataLab