Tcanopy: Compute below-canopy air temperature

In ilyamaclean/microctools: Various worker functions for microclimc package

Compute below-canopy air temperature

Description

Function to compute below canopy air temperatures

Usage

Tcanopy(
  climdata,
  z,
  lat,
  long,
  hgt,
  pai,
  paia = NA,
  x = 1,
  gsmax = 0.33,
  lw = 0.05,
  slope = 0,
  aspect = 0,
  alb = 0.23,
  galb = 0.15,
  zu = 2,
  Bowen = 1.5,
  upper = 25,
  umin = 0.5,
  merid = 0,
  dst = 0,
  clump = 0,
  method = "C",
  lims = TRUE
)

Arguments

climdata	a dataframe of hourly weather data formated and with units as per the internal dataset climdata
z	the height (below canopy) for which temperature estimates are required (m)
lat	the latitude of the location (decimal degrees)
long	the longitude of the location (decimal degrees)
hgt	canopy height (m)
pai	the total one sided area of canopy elements per unit ground area (see details)
paia	optionally, the total one sided area of canopy elements per unit ground area above z (see details)
x	leaf distribution angle coefficient
gsmax	maximum stomatal conductance of leaves (mol / m^2 / s)
lw	average leaf width (m)
slope	the slope of the underlying ground surface (decimal degrees). Ignored if method = S.
aspect	the aspect of the underlying ground surface (decimal degrees). Ignored if method = S.
alb	the albedo of the canopy surface (either the combined ground and canopy albedo if method = S or just the canopy albedo)
galb	ground surface albedo
zu	the height above ground of wind speeds in climdata (m)
Bowen	Optional parameter specifying the Bowen Ratio of the surface. Used to improve estimates of temperature in application of the Penman-Monteith equation
upper	optional upper limit to temperature offset (difference between reference and canopy surface temperature cannot exceed this value). Ignored if lims = FALSE
umin	optional minimum wind speed for computing conductances (avoids conductances being too low)
merid	an optional numeric value representing the longitude (decimal degrees) of the local time zone meridian (0 for GMT)
dst	an optional numeric value representing the time difference from the timezone meridian (hours, e.g. +1 for BST if merid = 0).
clump	clumpiness factor (0-1, see details)
method	if S treats the vegetation surface as a flat surface (see details)
lims	optional logical indicating whether to limit temperatures by upper and dewpoint temperature

Details

if pai is unknown it can be estimated as -ln(1-fractional canopy cover). If paia is unspecified vertically uniform foliage distribution is assumed. if clump = 0 the canopy is assumed entirely uniform and radiation transmission is as for a turbid medium. As clump approaches 1, the canopy is assumed to be increasingly patchy, such that a greater proportion of reaches the ground without being obscured by leaves. If method = S, the radiation intercepted by the canopy is assumed to be that for a flat surface. If method is not S, the radiation absorption by the ground and canopy surface are computed separately accounting for the inclination of the ground surface and the distribution of leaf angles.

Value

a list of the following elements:

tleaf a vector of leaf surface temperatures (deg C)

tn a vector of air temperatures (deg C)

rh a vector of relative humidities (percentage)

Examples

# Compute temperatures below canopy
Tln<-Tcanopy(climdata,0.25,50.2178,-5.32656,hgt=0.5,pai=2,x=1,gsmax=0.33,alb=0.23)
plot(Tln$tn,type="l") # air temperature
plot(Tln$tleaf,type="l") # leaf temperature

ilyamaclean/microctools documentation built on Jan. 25, 2023, 5:29 a.m.