## 2.7.1. Canopy Water[¶](#canopy-water "Permalink to this headline")
------------------------------------------------------------------

Liquid precipitation is either intercepted by the canopy, falls directly to the snow/soil surface (throughfall), or drips off the vegetation (canopy drip). Solid precipitation is treated similarly, with the addition of unloading of previously intercepted snow. Interception by vegetation is divided between liquid and solid phases \\(q\_{intr,\\,liq}\\) and \\(q\_{intr,\\,ice}\\) (kg m\-2 s\-1)

(2.7.2)[¶](#equation-7-2 "Permalink to this equation")\\\[q\_{intr,\\,liq} = f\_{pi,\\,liq} \\ q\_{rain}\\\]

(2.7.3)[¶](#equation-7-3 "Permalink to this equation")\\\[q\_{intr,\\,ice} = f\_{pi,\\,ice} \\ q\_{sno}\\\]

where \\(f\_{pi,\\,liq}\\) and \\(f\_{pi,\\,ice}\\) are the fractions of intercepted precipitation of rain and snow, respectively

(2.7.4)[¶](#equation-7-2b "Permalink to this equation")\\\[f\_{pi,\\,liq} = \\alpha\_{liq} \\ tanh \\left(L+S\\right)\\\]

(2.7.5)[¶](#equation-7-3b "Permalink to this equation")\\\[f\_{pi,\\,ice} =\\alpha\_{sno} \\ \\left\\{1-\\exp \\left\[-0.5\\left(L+S\\right)\\right\]\\right\\} \\ ,\\\]

and \\(L\\) and \\(S\\) are the exposed leaf and stem area index, respectively (section [2.2.1.4](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Ecosystem/CLM50_Tech_Note_Ecosystem.html#phenology-and-vegetation-burial-by-snow)), and the \\(\\alpha\\)'s scale the fractional area of a leaf that collects water ([Lawrence et al. 2007](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#lawrenceetal2007)). Default values of \\(\\alpha\_{liq}\\) and \\(\\alpha\_{sno}\\) are set to 1. Throughfall (kg m\-2 s\-1) is also divided into liquid and solid phases, reaching the ground (soil or snow surface) as

(2.7.6)[¶](#equation-7-4 "Permalink to this equation")\\\[q\_{thru,\\, liq} = q\_{rain} \\left(1 - f\_{pi,\\,liq}\\right)\\\]

(2.7.7)[¶](#equation-7-5 "Permalink to this equation")\\\[q\_{thru,\\, ice} = q\_{sno} \\left(1 - f\_{pi,\\,ice}\\right)\\\]

Similarly, the liquid and solid canopy drip fluxes are

(2.7.8)[¶](#equation-7-6 "Permalink to this equation")\\\[q\_{drip,\\, liq} =\\frac{W\_{can,\\,liq}^{intr} -W\_{can,\\,liq}^{max } }{\\Delta t} \\ge 0\\\]

(2.7.9)[¶](#equation-7-7 "Permalink to this equation")\\\[q\_{drip,\\, ice} =\\frac{W\_{can,\\,sno}^{intr} -W\_{can,\\,sno}^{max } }{\\Delta t} \\ge 0\\\]

where

(2.7.10)[¶](#equation-7-8 "Permalink to this equation")\\\[W\_{can,liq}^{intr} =W\_{can,liq}^{n} +q\_{intr,\\, liq} \\Delta t\\ge 0\\\]

and

(2.7.11)[¶](#equation-7-9 "Permalink to this equation")\\\[W\_{can,sno}^{intr} =W\_{can,sno}^{n} +q\_{intr,\\, ice} \\Delta t\\ge 0\\\]

are the the canopy liquid water and snow water equivalent after accounting for interception, \\(W\_{can,\\,liq}^{n}\\) and \\(W\_{can,\\,sno}^{n}\\) are the canopy liquid and snow water from the previous time step, and \\(W\_{can,\\,liq}^{max }\\) and \\(W\_{can,\\,snow}^{max }\\) (kg m\-2 or mm of H2O) are the maximum amounts of liquid water and snow the canopy can hold. They are defined by

(2.7.12)[¶](#equation-7-10 "Permalink to this equation")\\\[W\_{can,\\,liq}^{max } =p\_{liq}\\left(L+S\\right)\\\]

(2.7.13)[¶](#equation-7-11 "Permalink to this equation")\\\[W\_{can,\\,sno}^{max } =p\_{sno}\\left(L+S\\right).\\\]

The maximum storage of liquid water is \\(p\_{liq}=0.1\\) kg m\-2 ([Dickinson et al. 1993](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#dickinsonetal1993)), and that of snow is \\(p\_{sno}=6\\) kg m\-2, consistent with reported field measurements ([Pomeroy et al. 1998](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#pomeroyetal1998)).

Canopy snow unloading from wind speed \\(u\\) and above-freezing temperatures are modeled from linear fluxes and e-folding times similar to [Roesch et al. (2001)](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#roeschetal2001)

(2.7.14)[¶](#equation-7-12 "Permalink to this equation")\\\[q\_{unl,\\, wind} =\\frac{u W\_{can,sno}}{1.56\\times 10^5 \\text{ m}}\\\]

(2.7.15)[¶](#equation-7-13 "Permalink to this equation")\\\[q\_{unl,\\, temp} =\\frac{W\_{can,sno}(T-270 \\textrm{ K})}{1.87\\times 10^5 \\text{ K s}} > 0\\\]

(2.7.16)[¶](#equation-7-14 "Permalink to this equation")\\\[q\_{unl,\\, tot} =\\min \\left( q\_{unl,\\, wind} +q\_{unl,\\, temp} ,W\_{can,\\, sno} \\right)\\\]

The canopy liquid water and snow water equivalent are updated as

(2.7.17)[¶](#equation-7-15 "Permalink to this equation")\\\[ W\_{can,\\, liq}^{n+1} =W\_{can,liq}^{n} + q\_{intr,\\, liq} - q\_{drip,\\, liq} \\Delta t - E\_{v}^{liq} \\Delta t \\ge 0\\\]

and

(2.7.18)[¶](#equation-7-16 "Permalink to this equation")\\\[W\_{can,\\, sno}^{n+1} =W\_{can,sno}^{n} + q\_{intr,\\, ice} - \\left(q\_{drip,\\, ice}+q\_{unl,\\, tot} \\right)\\Delta t - E\_{v}^{ice} \\Delta t \\ge 0\\\]

where \\(E\_{v}^{liq}\\) and \\(E\_{v}^{ice}\\) are partitioned from the stem and leaf surface evaporation \\(E\_{v}^{w}\\) (Chapter [2.5](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Fluxes/CLM50_Tech_Note_Fluxes.html#rst-momentum-sensible-heat-and-latent-heat-fluxes)) based on the vegetation temperature \\(T\_{v}\\) (K) (Chapter [2.5](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Fluxes/CLM50_Tech_Note_Fluxes.html#rst-momentum-sensible-heat-and-latent-heat-fluxes)) and its relation to the freezing temperature of water \\(T\_{f}\\) (K) ([Table 2.2.7](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Ecosystem/CLM50_Tech_Note_Ecosystem.html#table-physical-constants))

(2.7.19)[¶](#equation-7-17 "Permalink to this equation")\\\[\\begin{split}E\_{v}^{liq} = \\left\\{\\begin{array}{lr} E\_{v}^{w} & T\_v > T\_{f} \\\\ 0 & T\_v \\le T\_f \\end{array}\\right\\}\\end{split}\\\]

(2.7.20)[¶](#equation-7-18 "Permalink to this equation")\\\[\\begin{split}E\_{v}^{ice} = \\left\\{\\begin{array}{lr} 0 & T\_v > T\_f \\\\ E\_{v}^{w} & T\_v \\le T\_f \\end{array}\\right\\}.\\end{split}\\\]

The total rate of liquid and solid precipitation reaching the ground is then

(2.7.21)[¶](#equation-7-19 "Permalink to this equation")\\\[q\_{grnd,liq} =q\_{thru,\\, liq} +q\_{drip,\\, liq}\\\]

(2.7.22)[¶](#equation-7-20 "Permalink to this equation")\\\[q\_{grnd,ice} =q\_{thru,\\, ice} +q\_{drip,\\, ice} +q\_{unl,\\, tot} .\\\]

Solid precipitation reaching the soil or snow surface, \\(q\_{grnd,\\, ice} \\Delta t\\), is added immediately to the snow pack (Chapter [2.8](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Snow_Hydrology/CLM50_Tech_Note_Snow_Hydrology.html#rst-snow-hydrology)). The liquid part, \\(q\_{grnd,\\, liq} \\Delta t\\) is added after surface fluxes (Chapter [2.5](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Fluxes/CLM50_Tech_Note_Fluxes.html#rst-momentum-sensible-heat-and-latent-heat-fluxes)) and snow/soil temperatures (Chapter [2.6](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Soil_Snow_Temperatures/CLM50_Tech_Note_Soil_Snow_Temperatures.html#rst-soil-and-snow-temperatures)) have been determined.

The wetted fraction of the canopy (stems plus leaves), which is required for surface flux (Chapter [2.5](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Fluxes/CLM50_Tech_Note_Fluxes.html#rst-momentum-sensible-heat-and-latent-heat-fluxes)) calculations, is ([Dickinson et al.1993](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#dickinsonetal1993))

(2.7.23)[¶](#equation-7-21 "Permalink to this equation")\\\[\\begin{split}f\_{wet} = \\left\\{\\begin{array}{lr} \\left\[\\frac{W\_{can} }{p\_{liq}\\left(L+S\\right)} \\right\]^{{2\\mathord{\\left/ {\\vphantom {2 3}} \\right.} 3} } \\le 1 & \\qquad L+S > 0 \\\\ 0 &\\qquad L+S = 0 \\end{array}\\right\\}\\end{split}\\\]

while the fraction of the canopy that is dry and transpiring is

(2.7.24)[¶](#equation-7-22 "Permalink to this equation")\\\[\\begin{split}f\_{dry} = \\left\\{\\begin{array}{lr} \\frac{\\left(1-f\_{wet} \\right)L}{L+S} & \\qquad L+S > 0 \\\\ 0 &\\qquad L+S = 0 \\end{array}\\right\\}.\\end{split}\\\]

Similarly, the snow-covered fraction of the canopy is used for surface alebdo when intercepted snow is present (Chapter [2.3](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Surface_Albedos/CLM50_Tech_Note_Surface_Albedos.html#rst-surface-albedos))

(2.7.25)[¶](#equation-7-23 "Permalink to this equation")\\\[\\begin{split}f\_{can,\\, sno} = \\left\\{\\begin{array}{lr} \\left\[\\frac{W\_{can,\\, sno} }{p\_{sno}\\left(L+S\\right)} \\right\]^{{3\\mathord{\\left/ {\\vphantom {3 20}} \\right.} 20} } \\le 1 & \\qquad L+S > 0 \\\\ 0 &\\qquad L+S = 0 \\end{array}\\right\\}.\\end{split}\\\]