# Green Roof Model (EcoRoof)[LINK]

The input object Material:RoofVegetation provides a model for green roofs (aka ecoroofs or vegetated roofs) that are becoming increasingly common for both new and retrofit buildings. There is widespread recognition and a growing literature of measured data that suggest green roofs can reduce building energy consumption. Currently, however, there are few design tools available to assist developers and architects in assessing the likely magnitude of energy savings associated with various implementation options (e.g., soil type/depth, irrigation options, plant type). As a result there is a significant need for a quantitative and physically-based building energy simulation tool that represents the effects of green roof constructions. Such a tool would facilitate more rapid spread of green roof technologies and make it possible to account for green roof benefits in state energy codes and related energy efficiency standards such as LEED.

In response to the need for green roof design tools a computational model of the heat transfer processes involved on a vegetated roof has been developed. This model accounts for:

long wave and short wave radiative exchange within the plant canopy,

plant canopy effects on convective heat transfer,

evapotranspiration from the soil and plants, and

heat conduction (and storage) in the soil layer

The ability to track moisture-dependent thermal properties is not implemented yet due to stability issues in the Conduction Transfer Function (CTF) scheme. Finite difference solution scheme currently is not supported in EnergyPlus Green Roof (EcoRoof) Model.

As implemented in EnergyPlus the green roof module allows the user to specify “ecoroof” as the outer layer of a rooftop construction using a “Material:RoofVegetation” object. The user can then specify various aspects of the green roof construction including growing media depth, thermal properties, plant canopy density, plant height, stomatal conductance (ability to transpire moisture), and soil moisture conditions (including irrigation).

The model formulation includes the following:

simplified moisture balance that allows precipitation, irrigation, and moisture transport between two soil layers (top and root zone).

soil and plant canopy energy balance based on the Army Corps of Engineers’ FASST vegetation models (Frankenstein and Koenig), drawing heavily from BATS (Dickenson et al.) and SiB (Sellers et al.).

soil surface (Tg) and foliage (Tf) temperature equations are solved simultaneously each time step, inverting the CTF to extract heat flux information for the energy balance calculation.

The detailed energy balance analysis and resulting equations, being rather complicated, are summarized here. The interested reader is referred to the FASST documentation cited herein for the complete development. The end result is a set of two simultaneous equations for temperature—one for the soil surface and the other for the foliage.

## Green Roof Model Description[LINK]

As with a traditional roof, the energy balance of an green roof is dominated by radiative forcing from the sun. This solar radiation is balanced by sensible (convection) and latent (evaporative) heat flux from soil and plant surfaces combined with conduction of heat into the soil substrate. This energy balance is illustrated in Figure 1. The variables introduced in this figure are defined in the equations that follow.

The energy budget analysis follows the Fast All Season Soil Strength (FASST) model developed by Frankenstein and Koenig for the US Army Corps of Engineers. FASST was developed, in part, to determine the ability of soils to support manned and unmanned vehicles and personnel movement. In order to accomplish this, however, FASST tracks the energy and moisture balance (including ice and snow) within a vegetated soil. It is a one-dimensional model that draws heavily from other plant canopy models including BATS (Dickinson et al.) and SiB (Sellers et al.). We have implemented FASST here with only a few modifications to adapt it for use with a relatively thin soil layer. The sign convention used assumes all heat fluxes are positive when energy is absorbed into the layer.

In the following discussion this energy budget is divided into a budget for the foliage layer (Ff) and a budget for the ground surface (Fg). The various parameterizations for latent and sensible heat flux are described in some detail and then the equation set is reduced to the simultaneous solution of two equations involving the temperatures of the foliage and ground surface.

### Energy budget in the foliage layer[LINK]

The foliage energy balance is given by:

Ff=σf[IS↓(1−αf)+εfIir↓−εfσTf4]+σfεgεfσε1(Tg4−Tf4)+Hf+Lf

Where σf is calculated as a function of Leaf-Area-Index (LAI):

σf=0.9−0.7exp−0.75LAI

In addition to convective and sensible heat transfer this equation accounts for both the short and longwave radiation absorbed by the vegetation, including the effects of multiple reflections. The sensible and latent heat flux terms (Hf and Lf−) are somewhat complicated and therefore discussed in some detail below.

#### Sensible heat flux in the foliage layer[LINK]

The sensible heat transfer between the leaf surface and near-canopy air (Hf) is influenced by the temperature difference between them, wind speed, and Leaf Area Index (LAI). The Leaf Area Index is the dimensionless ratio of the projected leaf area for a unit ground area (Oke). In contrast fractional vegetative cover (sf) is the ratio of shaded ground surface to total ground surface area. The sensible heat flux is given by:

Hf=(1.1∗LAIρafCp,aCfWaf)∗(Taf−Tf)

In this equation the constant 1.1 accounts for heat transfer from the stems, twigs and limbs (Deardorff). The properties of air near the foliage are modeled using the average from the foliage and instrument conditions:

ρaf=0.5(ρa+ρf)

where ra is the density of air at the instrument height and rf is the density of air at the leaf temperature. The air temperature within the foliage is estimated by:

Taf=(1−σf)(Ta)+σf(0.3Ta+0.6Tf+0.1Tg)

*where, Ta* is the air temperature at the instrument height in Kelvin *Tf*, is leaf temperature in Kelvin and *Tg*, is the ground surface temperature in Kelvin. The foliage wind speed is estimated as:

Waf=0.83σfW√Cfhn+(1−σf)W

Here *W* is the larger of 2.0 m/s or the actual wind speed above the canopy (Hughes et al.) and Cfhn is the transfer coefficient at near-neutral atmospheric stability conditions:

Cfhn=Kv2⋅(ln(Za−ZdZfo))−2

where *Kv,* is von Karmen’s constant (0.4), Za is the instrument height, Zd is the zero displacement height in meters (height above soil within which the wind speed is effectively zero), and Zfo is the foliage roughness length scale (m). The formulations for zero displacement height, roughness length are based on Balick et al.:

Zd=0.701Z0.979f

Zo=0.131Z0.997f

Finally, the bulk transfer coefficient as defined by Deardorff is given by:

Cf=0.01∗(1+0.3(m/s)Waf(m/s))

#### Latent heat flux in the foliage layer[LINK]

The process of water loss through plant respiration is known as transpiration. It is controlled by the closing and opening of stomata - the intercellular openings between to epidermal (guard) cells (Gates). The resistance to the diffusion of water vapor from these spaces into the atmosphere is called stomatal resistance. It depends on factors such as light intensity, soil moisture content and vapor pressure difference between inside leaf and the outside atmosphere. It is measured in units of s/m and is formulated as:

rs=rs,minLAI⋅f1⋅f2⋅f3

Here, rs,min is the minimum stomatal resistance. The actual stomatal resistance at any time is proportional to this minimum resistance and inversely proportional to LAI. The stomatal resistance is further modified by fractional multiplying factors that relate to incoming solar radiation and atmospheric moisture. As found in Frankenstein and Koenig the inverses of the multiplying factors f1, f2, and f3 are given by:

1f1=min[1,0.004∗Is↓+0.0050.81∗(0.004∗Is↓+1)]1f2=⎧⎪⎨⎪⎩0whenθr>¯¯¯θ¯θ−θrθmax−θrwhenθr≤¯¯¯θ≤θmax1f3=exp[−gd(ef,sat−ea)]

Here, Qr, is the residual moisture content (defined as the amount of moisture in soil when plants begin to wilt), Qmax is the maximum moisture content (defined as the maximum amount of moisture a particular type of soil can hold and above which run off occurs), and ¯¯¯¯Θ is the average soil moisture in the root zone. The residual moisture content is typically around 0.01 m3/m3 (Frankenstein and Koenig). The maximum moisture content depends upon the soil, but generally varies from 0.3 to 0.6 m3/m3 (Guymon et al.). In the expression for f3, gd is a plant specific characteristic that is only non-zero for trees, ef,sat is the saturated vapor pressure at the leaf temperature, and ea is the air vapor pressure.

Resistance to moisture exchange offered by the boundary layer formed on the leaf surface is known as aerodynamic resistance. It is measured in units of (s/m) and is influenced by wind speed, surface roughness and stability of the atmosphere (Oke). It is formulated as:

ra=1cfWaf

The combined effect of aerodynamic and stomatal resistances to vapor diffusion is integrated into a foliage surface wetness factor:

r′′=rara+rs

This surface wetness factor is simply a ratio of the aerodynamic resistance to the total resistance. When the aerodynamic resistance is small the wetness factor approaches zero (leaf surfaces remain dry as surface moisture is readily evaporated). As the aerodynamic resistance increases in importance relative to stomatal resistance the wetness factor approaches 1.0 (moisture readily travels to the leaf surfaces, but is not easily evaporated).

The latent heat flux is then given by:

Lf=lf∗LAIρafCfWafr′′(qaf−qf,sat)

Here *lf* , is the latent heat of vaporization (J/kg), qf,sat is the saturation mixing ratio at the leaf surface temperature, and qaf is the mixing ratio of the air within the canopy. As developed in Frankenstein and Koenig the mixing ratio within the canopy can be determined from:

qaf=[(1−σf)qa+σf(0.3qa+0.6qf,satr′′+0.1qf,satMg)1−σf[0.6(1−r′′)+0.1(1−Mg)]]

where the factor *Mg* (ranging from 0 to 1) is the ratio of volumetric moisture content to the porosity of the soil (Koenig). The latent heat of vaporization (*lf*) is the amount of energy required to convert a unit mass of water to vapor. It is measured in units of J/kg and is inversely proportional to the temperature. From Henderson-Sellers it is estimated as:

lf=1.91846∗106[TfTf−33.91]2

### Soil Energy budget[LINK]

The energy budget at the soil surface is mainly influenced by the soil thermal properties, the amount of foliage coverage (sf) and the amount of moisture in the soil. If the soil surface is densely covered the diurnal range of surface temperature is small. In the soil energy budget the heat released or gained due to phase changes of soil water, precipitation heat flux and heat flux due to vertical transport of water in the soil are ignored. Future refinements to this model will incorporate these phenomena. The sign convention followed here is the same as above (heat flux into the soil is positive). The overall energy balance at the soil surface (as given in Frankenstein and Koenig) is:

Fg=(1−σf)[I↓s(1−αg)+εgI↓ir−εgT4g]−σfεgεfσε1(T4g−T4f)+Hg+Lg+K∗∂Tg∂z

As with the energy equation for the foliage this equation represents sensible heat flux (Hg), latent heat flux (Lg) and the multiple reflections associated with long and short wave radiation. The final term on the right side gives the conduction of heat into the soil substrate.

#### Sensible heat flux in the soil layer[LINK]

Sensible heat flux between the soil surface and air in its vicinity is dependent on the temperature difference between them and the wind speed within the canopy. It is given as

Hg=ρagCp,aCghWaf(Taf−Tg)

where Cgh is the bulk transfer coefficient and rag is the density of air near the soil surface (kg/m3) given by:

pag=pa+pg2

Here rg is the density of air at the ground surface temperature

The bulk transfer coefficient is given as the linear combination of bulk transfer coefficient near ground (Cfhn) and near foliage-atmosphere interface (Cghn) multiplied by the stability factor (Gh) and is formulated as:

Cgh=Γk[(1−σf)Cghn+σfCfhn]

The ground and foliage bulk transfer coefficients, in turn, are given by:

Cghn=r−1ch⎡⎢
⎢
⎢
⎢⎣Kvln(ZaZgo)⎤⎥
⎥
⎥
⎥⎦2

And

Cfhn=⎡⎢
⎢
⎢
⎢⎣Kvln(Za−ZdZfo)⎤⎥
⎥
⎥
⎥⎦2

whereZgo and Zfo are the ground and foliage roughness lengths, rch is turbulent Schmidt number (0.63), and Kv is the von Karman constant (0.4).

The condition of the atmosphere (Gh) is determined as stable or unstable based on the sign of the bulk Richardson number:

Rib=2gZa(Taf−Tg)(Taf+Tg)W2af

The atmospheric stability factor is then given by Businger and Lumley and Panofsky as:

Γh=⎧⎪⎨⎪⎩1.0(1.0−16.0Rib)0.5forRib<01.0(1.0−5.0Rib)forRib>0

#### Latent heat flux in the soil layer[LINK]

Removal of water vapor from the soil surface depends on the difference between the mixing ratio of the soil surface and air and the wind speed within the canopy. The resulting latent heat flux is then given by:

LR=CgelgWafρag(qaf−qg)

Here Cge is the bulk transfer coefficient, lg is the latent heat of vaporization at the ground surface temperature, qaf is the mixing ratio at the foliage-atmosphere interface, and qf is the mixing ratio at the ground surface, given by:

qg=Mgqg,sat+(1−Mg)qaf

The bulk transfer coefficient for latent heat exchange is analogous to that for sensible heat exchange and is given by:

Cge=Γe[(1−σf)Cgen+σfCfhn]

where Cgen is the near ground bulk transfer coefficient for Latent heat flux and Ge is the latent heat exchange stability correction factor (assumed to be the same as Gh).

## Linearization[LINK]

In order to solve the foliage and soil heat budget equations, the 4th order terms Tf4 and Tg4 and mixing ratio terms qg,sat and qf,sat are linearized as given by Deardorff:

[T(n+1)f]4=[Tnf]4+4[Tnf]3[Tn+1f−Tnf]

[T(n+1)g]4=[Tng]4+4[Tng]3[Tn+1g−Tng]

Here Tfn+1 and Tgn+1 are the current time step leaf and ground surface temperatures in Kelvin. Tfn and Tgn are the corresponding temperatures at the previous time step.

The saturation mixing ratio at the ground and leaf surface temperatures are given as:

qg,sat(Tn+1g)=qsat(Tng)+(∂qsat∂T)Tng∗(Tn+1g−Tng)

qf,sat(Tn+1f)=qsat(Tnf)+(∂qsat∂T)Tnf∗(Tn+1f−Tnf)

where qsat(Tgn) is the saturation mixing ratio at the previous time step and is formulated as given in Garratt:

qsat(Tng)=0.622e∗(Tng)P−e∗(Tng)

Here the saturation vapor pressure e* (Pa) is evaluated at the ground temperature from the previous time step (Tgn) as:

e∗=611.2exp[17.67(Tng−273.15Tng−29.65)]

The derivative of saturation mixing ratio at the previous time step is given by:

dq∗dTng=[0.622∗P(P−0.378∗e∗)2](de∗dTng)

Here, the derivative of the saturation vapor pressure can be calculated from the Clausius-Clapeyron equation:

de∗dTng=lg∗e∗(Tng)(Rv∗(Tng)2)

Where Rv is the gas constant for water vapor and lg is the latent heat of vaporization at the soil surface temperature.

The corresponding saturation mixing ratio relations for the leaf surfaces can be obtained by replacing Tg with Tf in the above relations.

## Final Equations[LINK]

After linearization the final equations are of the form:

Cf1+Cf2Tg+Cf3Tf=0

Cg1+Cg2Tg+Cg3Tf=0

The coefficients in these equations result from the direct combination of the equations from the above development. The interested reader is directed to the papers by Frankenstein and Koenig for the complete and somewhat complicated expressions.

This final set of equations is then solved simultaneously to obtain T*g* and T*f* . One key difference in our implementation of the FASST algorithm is that the conduction terms in the equations for Cg1 and Cg2 are solved by inverting the Conduction Transfer Functions (CTF) within the EnergyPlus solution scheme.

## Green Roof Nomenclature[LINK]

C1, C2, C3 = coefficients in linearized temperature equations

*Ceg = * latent heat flux bulk transfer coefficient at ground layer

*Cf * = bulk heat transfer coefficient

Chg = sensible heat flux bulk transfer coefficient at ground layer

Chnf = near-neutral transfer coefficient at foliage layer

Chng = near-neutral transfer coefficient at ground layer

Cp,a = specific heat of air at constant pressure (1005.6 J/kg k)

e∗ = saturation vapor pressure (Pa)

f1 = multiplying factor for radiation effect on stomatal resistance

f2 = multiplying factor for moisture effect on stomatal resistance

f3 = additional multiplying factor for stomatal resistance

Ff = net heat flux to foliage layer (W/m2)

Fg = net heat flux to ground surface (W/m2)

gd = plant specific characteristic related to stomatal resistance

Hf = foliage sensible heat flux (W/m2)

Hg = ground sensible heat flux (W/m2)

I↓s = total incoming short wave radiation (W/m2)

I↓iγ = total incoming longwave radiation (W/m2)

Kv = von Karmen constant (0.4)

lf = latent heat of vaporization at foliage temperature (J/kg)

lg = latent heat of vaporization at ground temperature (J/kg)

L*f* = foliage latent heat flux (W/m2)

Lg = ground latent heat flux (W/m2)

*LAI* = leaf area index (m2/m2)

Mg = moisture saturation factor

qa = mixing ratio for air

qaf = mixing ratio for air within foliage canopy

qf,sat = saturation mixing ratio at foliage temperature

qg,sat = saturation mixing ratio at ground temperature

ra = aerodynamic resistance to transpiration (s/m)

rs = foliage leaf stomatal resistance (s/m)

rs,min = minimal leaf stomatal resistance (s/m)

r“ = surface wetness factor

Rib = bulk Richardson number

Rv = gas constant for water vapor (461.53 J/kgK)

*Ta* = the air temperature at the instrument height (Kelvin)

Taf = air temperature with in the canopy (Kelvin)

T*f* = leaf temperature (Kelvin)

Tg = ground surface temperature (Kelvin)

W = wind speed above canopy (m/s)

Waf = wind speed with in the canopy (m/s)

z = height or depth (m)

Za = instrument height (m)

Zd = displacement height (m)

Zof = foliage roughness length scale (m)

### Greek letters[LINK]

αf = albedo (short wave reflectivity) of the canopy

αg = albedo (short wave reflectivity) of ground surface

ε1 = εg + εf - εfεg

εf = emissivity of canopy

εg = emissivity of the ground surface

Γh = stability factor

ρa = density of air at instrument height (kg/m3)

ρf = density of air at foliage temperature (kg/m3)

ρaf = density of air at foliage temperature (kg/m3)

ρag = density of air at ground surface temperature (kg/m3)

σ = the Stefan-Boltzmann constant (5.699*108 W/m2 ∘K4)

σf = fractional vegetation coverage

### Subscripts and superscripts[LINK]

a = air

af = air within the foliage layer

e = latent heat flux term

f = foliage surface

g = ground surface

h = sensible heat flux term

n = current time step

n+1 = future time step

ir = infrared (or long-wave)

sat = saturation value

S = short-wave

ASHRAE. 2005. 2005 ASHRAE Handbook – Fundamentals. Chapter 16, Air flow Around Buildings, Atlanta: American Society of Heating Refrigerating and Air-Conditioning Engineers, Inc.

Balick, L. R., R. K. Scoggins, and L. E. Link. 1981. Inclusion of a simple vegetation layer in terrain temperature models for thermal IR signature prediction. IEEE Geoscience and Remote Sensing GE-19(3), pp.143-152.

Businger, J. A. 1966. In ‘Arctic Heat Budget and Atmospheric Circulation’, Symposium Proceedings, pp. 305-332. The Rand Corporation.

Deardorff, J.W. 1978. “Efficient Prediction of ground surface temperature and moisture with inclusion of a layer of vegetation”, Journal Geophysical Research, pp. 1889-1902.

Dickinson, R.E., A. Henderson-Sellers, P.J. Kennedy, and M.F. Wilson. 1986. Biosphere-Atmosphere Transfer Scheme (BATS) for the NCAR community climate model. NCAR Technical Note, TN-275+STR.

ECMWF. 2002. European Centre for Medium-Range Weather Forecasts, Integrated Forecast System. Documentation, CY25R1 (Operational implementation 9 April 2002). .

Frankenstein, S., and G. Koenig. 2004.FASST Vegetation Models. U. S. Army Engineer Research and Development Center, Cold regions Research and Engineering Laboratory, ERDC/CRREL Technical Report TR-04-25.

Frankenstein, S., and G. Koenig. 2004.Fast All-season Soil Strength (FASST). U.S. Army Engineer Research and Development Center, Cold regions Research and Engineering Laboratory, ERDC/CRREL Special Report SR-04-1.

Garratt, J.R. 1992. The Atmospheric Boundary Layer, Cambridge university press.

Gates, D.M. 1980. Biophysical Ecology. New York: Springer-Verlag

Guymon, G.L., R.L. Berg, and T.V. Hromadka. 1993. Mathematical Model of Frost Heave and Thaw Settlement in Pavements. U.S. Army Cold Regions Research and Engineering Laboratory, CRREL Report 93-2.

Henderson-Sellers, B. 1984. “A New Formula for Latent Heat of Vaporization of water as function of temperature”, Quarterly Journal Royal Meteorological Society, 10 pp. 1186-1190.

Hughes, P.A., T.J.L. McComb, A.B. Rimmer, and K.E. Turver. 1993. “A mathematical model for the prediction of temperature of man-made and natural surfaces”, International Journal of Remote Sensing 14 (7), pp. 1383-1412.

Koenig, G.G. 1994. Smart Weapons Operability Enhancement (SWOE) Joint Test and Evaluation (JT and E) Program: Final Report. Dr. James P. Welch, Joint Test Director, SWOE JT and E, SWOE Report 94-10, Annex D.

Lumley, J. L. and Panofsky, H. A. 1964. ‘The structure of Atmospheric Turbulence’. Interscience Monographs and Texts in Physics and Astronomy, Vol. XII. Wiley, New York.

Oke, T.R. 1987. Boundary Layer Climates, University Press, Cambridge

Sellers, P.J., Y. Mintz, Y.C. Sud, and A. Dalcher. 1986. A simple biosphere model (SiB) for use within general circulation models. Journal of Atmospheric Science, 43 (6), pp. 505-532.

## Green Roof Model (EcoRoof)[LINK]

## Overview[LINK]

The input object Material:RoofVegetation provides a model for green roofs (aka ecoroofs or vegetated roofs) that are becoming increasingly common for both new and retrofit buildings. There is widespread recognition and a growing literature of measured data that suggest green roofs can reduce building energy consumption. Currently, however, there are few design tools available to assist developers and architects in assessing the likely magnitude of energy savings associated with various implementation options (e.g., soil type/depth, irrigation options, plant type). As a result there is a significant need for a quantitative and physically-based building energy simulation tool that represents the effects of green roof constructions. Such a tool would facilitate more rapid spread of green roof technologies and make it possible to account for green roof benefits in state energy codes and related energy efficiency standards such as LEED.

In response to the need for green roof design tools a computational model of the heat transfer processes involved on a vegetated roof has been developed. This model accounts for:

long wave and short wave radiative exchange within the plant canopy,

plant canopy effects on convective heat transfer,

evapotranspiration from the soil and plants, and

heat conduction (and storage) in the soil layer

The ability to track moisture-dependent thermal properties is not implemented yet due to stability issues in the Conduction Transfer Function (CTF) scheme. Finite difference solution scheme currently is not supported in EnergyPlus Green Roof (EcoRoof) Model.

As implemented in EnergyPlus the green roof module allows the user to specify “ecoroof” as the outer layer of a rooftop construction using a “Material:RoofVegetation” object. The user can then specify various aspects of the green roof construction including growing media depth, thermal properties, plant canopy density, plant height, stomatal conductance (ability to transpire moisture), and soil moisture conditions (including irrigation).

The model formulation includes the following:

simplified moisture balance that allows precipitation, irrigation, and moisture transport between two soil layers (top and root zone).

soil and plant canopy energy balance based on the Army Corps of Engineers’ FASST vegetation models (Frankenstein and Koenig), drawing heavily from BATS (Dickenson et al.) and SiB (Sellers et al.).

soil surface (Tg) and foliage (Tf) temperature equations are solved simultaneously each time step, inverting the CTF to extract heat flux information for the energy balance calculation.

The detailed energy balance analysis and resulting equations, being rather complicated, are summarized here. The interested reader is referred to the FASST documentation cited herein for the complete development. The end result is a set of two simultaneous equations for temperature—one for the soil surface and the other for the foliage.

## Green Roof Model Description[LINK]

As with a traditional roof, the energy balance of an green roof is dominated by radiative forcing from the sun. This solar radiation is balanced by sensible (convection) and latent (evaporative) heat flux from soil and plant surfaces combined with conduction of heat into the soil substrate. This energy balance is illustrated in Figure 1. The variables introduced in this figure are defined in the equations that follow.

The Energy Balance for a Green Roof.

The energy budget analysis follows the Fast All Season Soil Strength (FASST) model developed by Frankenstein and Koenig for the US Army Corps of Engineers. FASST was developed, in part, to determine the ability of soils to support manned and unmanned vehicles and personnel movement. In order to accomplish this, however, FASST tracks the energy and moisture balance (including ice and snow) within a vegetated soil. It is a one-dimensional model that draws heavily from other plant canopy models including BATS (Dickinson et al.) and SiB (Sellers et al.). We have implemented FASST here with only a few modifications to adapt it for use with a relatively thin soil layer. The sign convention used assumes all heat fluxes are positive when energy is absorbed into the layer.

In the following discussion this energy budget is divided into a budget for the foliage layer (Ff) and a budget for the ground surface (Fg). The various parameterizations for latent and sensible heat flux are described in some detail and then the equation set is reduced to the simultaneous solution of two equations involving the temperatures of the foliage and ground surface.

## Energy budget in the foliage layer[LINK]

The foliage energy balance is given by:

Ff=σf[IS↓(1−αf)+εfIir↓−εfσTf4]+σfεgεfσε1(Tg4−Tf4)+Hf+Lf

Where σf is calculated as a function of Leaf-Area-Index (LAI):

σf=0.9−0.7exp−0.75LAI

In addition to convective and sensible heat transfer this equation accounts for both the short and longwave radiation absorbed by the vegetation, including the effects of multiple reflections. The sensible and latent heat flux terms (Hf and Lf−) are somewhat complicated and therefore discussed in some detail below.

## Sensible heat flux in the foliage layer[LINK]

The sensible heat transfer between the leaf surface and near-canopy air (Hf) is influenced by the temperature difference between them, wind speed, and Leaf Area Index (LAI). The Leaf Area Index is the dimensionless ratio of the projected leaf area for a unit ground area (Oke). In contrast fractional vegetative cover (sf) is the ratio of shaded ground surface to total ground surface area. The sensible heat flux is given by:

Hf=(1.1∗LAIρafCp,aCfWaf)∗(Taf−Tf)

In this equation the constant 1.1 accounts for heat transfer from the stems, twigs and limbs (Deardorff). The properties of air near the foliage are modeled using the average from the foliage and instrument conditions:

ρaf=0.5(ρa+ρf)

where ra is the density of air at the instrument height and rf is the density of air at the leaf temperature. The air temperature within the foliage is estimated by:

Taf=(1−σf)(Ta)+σf(0.3Ta+0.6Tf+0.1Tg)

where, Tais the air temperature at the instrument height in KelvinTf, is leaf temperature in Kelvin andTg, is the ground surface temperature in Kelvin. The foliage wind speed is estimated as:Waf=0.83σfW√Cfhn+(1−σf)W

Here

Wis the larger of 2.0 m/s or the actual wind speed above the canopy (Hughes et al.) and Cfhn is the transfer coefficient at near-neutral atmospheric stability conditions:Cfhn=Kv2⋅(ln(Za−ZdZfo))−2

where

Kv,is von Karmen’s constant (0.4), Za is the instrument height, Zd is the zero displacement height in meters (height above soil within which the wind speed is effectively zero), and Zfo is the foliage roughness length scale (m). The formulations for zero displacement height, roughness length are based on Balick et al.:Zd=0.701Z0.979f

Zo=0.131Z0.997f

Finally, the bulk transfer coefficient as defined by Deardorff is given by:

Cf=0.01∗(1+0.3(m/s)Waf(m/s))

## Latent heat flux in the foliage layer[LINK]

The process of water loss through plant respiration is known as transpiration. It is controlled by the closing and opening of stomata - the intercellular openings between to epidermal (guard) cells (Gates). The resistance to the diffusion of water vapor from these spaces into the atmosphere is called stomatal resistance. It depends on factors such as light intensity, soil moisture content and vapor pressure difference between inside leaf and the outside atmosphere. It is measured in units of s/m and is formulated as:

rs=rs,minLAI⋅f1⋅f2⋅f3

Here, rs,min is the minimum stomatal resistance. The actual stomatal resistance at any time is proportional to this minimum resistance and inversely proportional to LAI. The stomatal resistance is further modified by fractional multiplying factors that relate to incoming solar radiation and atmospheric moisture. As found in Frankenstein and Koenig the inverses of the multiplying factors f1, f2, and f3 are given by:

1f1=min[1,0.004∗Is↓+0.0050.81∗(0.004∗Is↓+1)]1f2=⎧⎪⎨⎪⎩0whenθr>¯¯¯θ¯θ−θrθmax−θrwhenθr≤¯¯¯θ≤θmax1f3=exp[−gd(ef,sat−ea)]

Here, Qr, is the residual moisture content (defined as the amount of moisture in soil when plants begin to wilt), Qmax is the maximum moisture content (defined as the maximum amount of moisture a particular type of soil can hold and above which run off occurs), and ¯¯¯¯Θ is the average soil moisture in the root zone. The residual moisture content is typically around 0.01 m3/m3 (Frankenstein and Koenig). The maximum moisture content depends upon the soil, but generally varies from 0.3 to 0.6 m3/m3 (Guymon et al.). In the expression for f3, gd is a plant specific characteristic that is only non-zero for trees, ef,sat is the saturated vapor pressure at the leaf temperature, and ea is the air vapor pressure.

Resistance to moisture exchange offered by the boundary layer formed on the leaf surface is known as aerodynamic resistance. It is measured in units of (s/m) and is influenced by wind speed, surface roughness and stability of the atmosphere (Oke). It is formulated as:

ra=1cfWaf

The combined effect of aerodynamic and stomatal resistances to vapor diffusion is integrated into a foliage surface wetness factor:

r′′=rara+rs

This surface wetness factor is simply a ratio of the aerodynamic resistance to the total resistance. When the aerodynamic resistance is small the wetness factor approaches zero (leaf surfaces remain dry as surface moisture is readily evaporated). As the aerodynamic resistance increases in importance relative to stomatal resistance the wetness factor approaches 1.0 (moisture readily travels to the leaf surfaces, but is not easily evaporated).

The latent heat flux is then given by:

Lf=lf∗LAIρafCfWafr′′(qaf−qf,sat)

Here

lf, is the latent heat of vaporization (J/kg), qf,sat is the saturation mixing ratio at the leaf surface temperature, and qaf is the mixing ratio of the air within the canopy. As developed in Frankenstein and Koenig the mixing ratio within the canopy can be determined from:qaf=[(1−σf)qa+σf(0.3qa+0.6qf,satr′′+0.1qf,satMg)1−σf[0.6(1−r′′)+0.1(1−Mg)]]

where the factor

Mg(ranging from 0 to 1) is the ratio of volumetric moisture content to the porosity of the soil (Koenig). The latent heat of vaporization (lf) is the amount of energy required to convert a unit mass of water to vapor. It is measured in units of J/kg and is inversely proportional to the temperature. From Henderson-Sellers it is estimated as:lf=1.91846∗106[TfTf−33.91]2

## Soil Energy budget[LINK]

The energy budget at the soil surface is mainly influenced by the soil thermal properties, the amount of foliage coverage (sf) and the amount of moisture in the soil. If the soil surface is densely covered the diurnal range of surface temperature is small. In the soil energy budget the heat released or gained due to phase changes of soil water, precipitation heat flux and heat flux due to vertical transport of water in the soil are ignored. Future refinements to this model will incorporate these phenomena. The sign convention followed here is the same as above (heat flux into the soil is positive). The overall energy balance at the soil surface (as given in Frankenstein and Koenig) is:

Fg=(1−σf)[I↓s(1−αg)+εgI↓ir−εgT4g]−σfεgεfσε1(T4g−T4f)+Hg+Lg+K∗∂Tg∂z

As with the energy equation for the foliage this equation represents sensible heat flux (Hg), latent heat flux (Lg) and the multiple reflections associated with long and short wave radiation. The final term on the right side gives the conduction of heat into the soil substrate.

## Sensible heat flux in the soil layer[LINK]

Sensible heat flux between the soil surface and air in its vicinity is dependent on the temperature difference between them and the wind speed within the canopy. It is given as

Hg=ρagCp,aCghWaf(Taf−Tg)

where Cgh is the bulk transfer coefficient and rag is the density of air near the soil surface (kg/m3) given by:

pag=pa+pg2

Here rg is the density of air at the ground surface temperature

The bulk transfer coefficient is given as the linear combination of bulk transfer coefficient near ground (Cfhn) and near foliage-atmosphere interface (Cghn) multiplied by the stability factor (Gh) and is formulated as:

Cgh=Γk[(1−σf)Cghn+σfCfhn]

The ground and foliage bulk transfer coefficients, in turn, are given by:

Cghn=r−1ch⎡⎢ ⎢ ⎢ ⎢⎣Kvln(ZaZgo)⎤⎥ ⎥ ⎥ ⎥⎦2

And

Cfhn=⎡⎢ ⎢ ⎢ ⎢⎣Kvln(Za−ZdZfo)⎤⎥ ⎥ ⎥ ⎥⎦2

whereZgo and Zfo are the ground and foliage roughness lengths, rch is turbulent Schmidt number (0.63), and Kv is the von Karman constant (0.4).

The condition of the atmosphere (Gh) is determined as stable or unstable based on the sign of the bulk Richardson number:

Rib=2gZa(Taf−Tg)(Taf+Tg)W2af

The atmospheric stability factor is then given by Businger and Lumley and Panofsky as:

Γh=⎧⎪⎨⎪⎩1.0(1.0−16.0Rib)0.5forRib<01.0(1.0−5.0Rib)forRib>0

## Latent heat flux in the soil layer[LINK]

Removal of water vapor from the soil surface depends on the difference between the mixing ratio of the soil surface and air and the wind speed within the canopy. The resulting latent heat flux is then given by:

LR=CgelgWafρag(qaf−qg)

Here Cge is the bulk transfer coefficient, lg is the latent heat of vaporization at the ground surface temperature, qaf is the mixing ratio at the foliage-atmosphere interface, and qf is the mixing ratio at the ground surface, given by:

qg=Mgqg,sat+(1−Mg)qaf

The bulk transfer coefficient for latent heat exchange is analogous to that for sensible heat exchange and is given by:

Cge=Γe[(1−σf)Cgen+σfCfhn]

where Cgen is the near ground bulk transfer coefficient for Latent heat flux and Ge is the latent heat exchange stability correction factor (assumed to be the same as Gh).

## Linearization[LINK]

In order to solve the foliage and soil heat budget equations, the 4th order terms Tf4 and Tg4 and mixing ratio terms qg,sat and qf,sat are linearized as given by Deardorff:

[T(n+1)f]4=[Tnf]4+4[Tnf]3[Tn+1f−Tnf]

[T(n+1)g]4=[Tng]4+4[Tng]3[Tn+1g−Tng]

Here Tfn+1 and Tgn+1 are the current time step leaf and ground surface temperatures in Kelvin. Tfn and Tgn are the corresponding temperatures at the previous time step.

The saturation mixing ratio at the ground and leaf surface temperatures are given as:

qg,sat(Tn+1g)=qsat(Tng)+(∂qsat∂T)Tng∗(Tn+1g−Tng)

qf,sat(Tn+1f)=qsat(Tnf)+(∂qsat∂T)Tnf∗(Tn+1f−Tnf)

where qsat(Tgn) is the saturation mixing ratio at the previous time step and is formulated as given in Garratt:

qsat(Tng)=0.622e∗(Tng)P−e∗(Tng)

Here the saturation vapor pressure e* (Pa) is evaluated at the ground temperature from the previous time step (Tgn) as:

e∗=611.2exp[17.67(Tng−273.15Tng−29.65)]

The derivative of saturation mixing ratio at the previous time step is given by:

dq∗dTng=[0.622∗P(P−0.378∗e∗)2](de∗dTng)

Here, the derivative of the saturation vapor pressure can be calculated from the Clausius-Clapeyron equation:

de∗dTng=lg∗e∗(Tng)(Rv∗(Tng)2)

Where Rv is the gas constant for water vapor and lg is the latent heat of vaporization at the soil surface temperature.

The corresponding saturation mixing ratio relations for the leaf surfaces can be obtained by replacing Tg with Tf in the above relations.

## Final Equations[LINK]

After linearization the final equations are of the form:

Cf1+Cf2Tg+Cf3Tf=0

Cg1+Cg2Tg+Cg3Tf=0

The coefficients in these equations result from the direct combination of the equations from the above development. The interested reader is directed to the papers by Frankenstein and Koenig for the complete and somewhat complicated expressions.

This final set of equations is then solved simultaneously to obtain T

gand Tf. One key difference in our implementation of the FASST algorithm is that the conduction terms in the equations for Cg1 and Cg2 are solved by inverting the Conduction Transfer Functions (CTF) within the EnergyPlus solution scheme.## Green Roof Nomenclature[LINK]

C1, C2, C3 = coefficients in linearized temperature equations

Ceg =latent heat flux bulk transfer coefficient at ground layerCf= bulk heat transfer coefficientChg = sensible heat flux bulk transfer coefficient at ground layer

Chnf = near-neutral transfer coefficient at foliage layer

Chng = near-neutral transfer coefficient at ground layer

Cp,a = specific heat of air at constant pressure (1005.6 J/kg k)

e∗ = saturation vapor pressure (Pa)

f1 = multiplying factor for radiation effect on stomatal resistance

f2 = multiplying factor for moisture effect on stomatal resistance

f3 = additional multiplying factor for stomatal resistance

Ff = net heat flux to foliage layer (W/m2)

Fg = net heat flux to ground surface (W/m2)

gd = plant specific characteristic related to stomatal resistance

Hf = foliage sensible heat flux (W/m2)

Hg = ground sensible heat flux (W/m2)

I↓s = total incoming short wave radiation (W/m2)

I↓iγ = total incoming longwave radiation (W/m2)

Kv = von Karmen constant (0.4)

lf = latent heat of vaporization at foliage temperature (J/kg)

lg = latent heat of vaporization at ground temperature (J/kg)

L

f= foliage latent heat flux (W/m2)Lg = ground latent heat flux (W/m2)

LAI= leaf area index (m2/m2)Mg = moisture saturation factor

qa = mixing ratio for air

qaf = mixing ratio for air within foliage canopy

qf,sat = saturation mixing ratio at foliage temperature

qg,sat = saturation mixing ratio at ground temperature

ra = aerodynamic resistance to transpiration (s/m)

rs = foliage leaf stomatal resistance (s/m)

rs,min = minimal leaf stomatal resistance (s/m)

r“ = surface wetness factor

Rib = bulk Richardson number

Rv = gas constant for water vapor (461.53 J/kgK)

Ta= the air temperature at the instrument height (Kelvin)Taf = air temperature with in the canopy (Kelvin)

T

f= leaf temperature (Kelvin)Tg = ground surface temperature (Kelvin)

W = wind speed above canopy (m/s)

Waf = wind speed with in the canopy (m/s)

z = height or depth (m)

Za = instrument height (m)

Zd = displacement height (m)

Zof = foliage roughness length scale (m)

## Greek letters[LINK]

αf = albedo (short wave reflectivity) of the canopy

αg = albedo (short wave reflectivity) of ground surface

ε1 = εg + εf - εfεg

εf = emissivity of canopy

εg = emissivity of the ground surface

Γh = stability factor

ρa = density of air at instrument height (kg/m3)

ρf = density of air at foliage temperature (kg/m3)

ρaf = density of air at foliage temperature (kg/m3)

ρag = density of air at ground surface temperature (kg/m3)

σ = the Stefan-Boltzmann constant (5.699*108 W/m2 ∘K4)

σf = fractional vegetation coverage

## Subscripts and superscripts[LINK]

a = air

af = air within the foliage layer

e = latent heat flux term

f = foliage surface

g = ground surface

h = sensible heat flux term

n = current time step

n+1 = future time step

ir = infrared (or long-wave)

sat = saturation value

S = short-wave

## References[LINK]

ASHRAE. 2005. 2005 ASHRAE Handbook – Fundamentals. Chapter 16, Air flow Around Buildings, Atlanta: American Society of Heating Refrigerating and Air-Conditioning Engineers, Inc.

Balick, L. R., R. K. Scoggins, and L. E. Link. 1981. Inclusion of a simple vegetation layer in terrain temperature models for thermal IR signature prediction. IEEE Geoscience and Remote Sensing GE-19(3), pp.143-152.

Businger, J. A. 1966. In ‘Arctic Heat Budget and Atmospheric Circulation’, Symposium Proceedings, pp. 305-332. The Rand Corporation.

Deardorff, J.W. 1978. “Efficient Prediction of ground surface temperature and moisture with inclusion of a layer of vegetation”, Journal Geophysical Research, pp. 1889-1902.

Dickinson, R.E., A. Henderson-Sellers, P.J. Kennedy, and M.F. Wilson. 1986. Biosphere-Atmosphere Transfer Scheme (BATS) for the NCAR community climate model. NCAR Technical Note, TN-275+STR.

ECMWF. 2002. European Centre for Medium-Range Weather Forecasts, Integrated Forecast System. Documentation, CY25R1 (Operational implementation 9 April 2002). .

Frankenstein, S., and G. Koenig. 2004.FASST Vegetation Models. U. S. Army Engineer Research and Development Center, Cold regions Research and Engineering Laboratory, ERDC/CRREL Technical Report TR-04-25.

Frankenstein, S., and G. Koenig. 2004.Fast All-season Soil Strength (FASST). U.S. Army Engineer Research and Development Center, Cold regions Research and Engineering Laboratory, ERDC/CRREL Special Report SR-04-1.

Garratt, J.R. 1992. The Atmospheric Boundary Layer, Cambridge university press.

Gates, D.M. 1980. Biophysical Ecology. New York: Springer-Verlag

Guymon, G.L., R.L. Berg, and T.V. Hromadka. 1993. Mathematical Model of Frost Heave and Thaw Settlement in Pavements. U.S. Army Cold Regions Research and Engineering Laboratory, CRREL Report 93-2.

Henderson-Sellers, B. 1984. “A New Formula for Latent Heat of Vaporization of water as function of temperature”, Quarterly Journal Royal Meteorological Society, 10 pp. 1186-1190.

Hughes, P.A., T.J.L. McComb, A.B. Rimmer, and K.E. Turver. 1993. “A mathematical model for the prediction of temperature of man-made and natural surfaces”, International Journal of Remote Sensing 14 (7), pp. 1383-1412.

Koenig, G.G. 1994. Smart Weapons Operability Enhancement (SWOE) Joint Test and Evaluation (JT and E) Program: Final Report. Dr. James P. Welch, Joint Test Director, SWOE JT and E, SWOE Report 94-10, Annex D.

Lumley, J. L. and Panofsky, H. A. 1964. ‘The structure of Atmospheric Turbulence’. Interscience Monographs and Texts in Physics and Astronomy, Vol. XII. Wiley, New York.

Oke, T.R. 1987. Boundary Layer Climates, University Press, Cambridge

Sellers, P.J., Y. Mintz, Y.C. Sud, and A. Dalcher. 1986. A simple biosphere model (SiB) for use within general circulation models. Journal of Atmospheric Science, 43 (6), pp. 505-532.

Documentation content copyright © 1996-2022 The Board of Trustees of the University of Illinois and the Regents of the University of California through the Ernest Orlando Lawrence Berkeley National Laboratory. All rights reserved. EnergyPlus is a trademark of the US Department of Energy.

This documentation is made available under the EnergyPlus Open Source License v1.0.