Laserfiche WebLink
288 D. C. Goodrich et al. /Agricultural and Forest Meteorology 105 (2000) 181 -309 <br />measurement periods. Due to limited availability of <br />heat pulse velocity probes, sap flux measurements <br />were not continuously made for the entire grow- <br />ing season. The measurements alternated between <br />the perennial Lewis Springs site and the ephemeral <br />Escapule site (Schaeffer et al., 2000). To estimate <br />C/W transpiration between, and outside the synoptic <br />measurement periods, a model driven by more easily <br />measured continuous variables is required. <br />4.2. Penman — Monteith (P —M) model for C/W <br />transpiration <br />The well -known P —M model (Monteith and <br />Unsworth, 1990) was selected to model C/W transpi- <br />ration throughout the riparian growing season. The <br />P —M equation for evaporation is given by <br />_ AA + pacpD/ra <br />1 <br />�E A +y(1 +rclra) <br />where XE is the evaporation (W m -2), A the slope <br />of the saturation vapor pressure /temperature curve <br />(kPa °C -1), A the available energy (Wm-2), Pa the <br />density of moist air (kgm3), cp = 1013Jkg °C -1 <br />the specific heat capacity of dry air under constant <br />pressure, D the vapor pressure deficit (kPa), ra the <br />aerodynamic resistance (s m -1), y the psychrometric <br />constant (kPa °C -1), and r, the bulk canopy resistance <br />(s m -1). In the above equation, A, pa, D, and y can be <br />approximated by formulae described by Shuttleworth <br />(1993), based on measurements of air temperature, <br />Ta ( °C), relative humidity, RH ( %), and atmospheric <br />pressure, P (kPa). These quantities were measured <br />at the nearby mesquite site (Scott et al., 2000), and <br />were found to reasonably approximate conditions in- <br />side the C/W canopy. This assessment was based on <br />comparisons to a limited set of measurements made <br />from a 12.5 in tower inside the C/W canopy. Measure- <br />ments from this tower were available in 1997 from <br />DOY 190 to 290. For this period of time, the root <br />mean square error (RMSE) and R2 between the vapor <br />pressure deficit computed from measurements at the <br />C/W tower and the mesquite tower were 0.26 kPa and <br />0.96, respectively. For air temperature, the RMSE <br />was 1.21 °C, with R2 = 0.97 between the two towers. <br />The available energy to the canopy is given by <br />A = S J (1 — a) + Lnet — St (2) <br />where SJ is the incoming solar radiation (W m -2), a <br />the canopy albedo, Lnet the net long -wave radiation <br />(W m -2), and St the temporary storage of energy into <br />the tree itself (trunk and limbs) and the energy used in <br />the photosynthesis process (W m -2). Because the bulk <br />of the canopy is typically 10-20 in above the ground, <br />soil heat flux contributions to the available energy for <br />the canopy were considered negligible. The incom- <br />ing solar radiation was measured over the mesquite <br />site. Canopy albedo was estimated to be 0. 18, a value <br />which has been measured over broadleaf oak trees <br />(Bras, 1990). St was estimated to be 5% of the incom- <br />ing solar radiation based on work by Moore and Fisch <br />(1986) who found that St ranged between 0 and 10% <br />of the net radiation available to a tropical forest. The <br />net long -wave radiation contribution to the available <br />energy was calculated from a formula provided by <br />Shuttleworth (1993, p. 4.7). Because the P —M model <br />only applies to tree transpiration, the net radiation was <br />adjusted to account for the portion not intercepted <br />by the leaves using a simple Beer's law relationship <br />from Shuttleworth and Gurney (1990). This adjust- <br />ment reduced the net radiation available to the canopy <br />by 25 %. <br />The aerodynamic resistance (ra) was assumed to <br />be the sum of the turbulent resistance between the <br />canopy and the atmosphere from turbulent eddies and <br />the boundary layer resistance (Thom, 1975). Due to <br />the relatively open nature of the cottonwood canopy, <br />the turbulent canopy resistance is assumed negligible <br />in comparison to the boundary layer resistance. Thus <br />ra is assumed to equal the boundary layer resistance <br />(rb). To estimate the boundary layer resistance, the <br />model proposed by Choudhury and Monteith (1988) <br />was used: <br />1 aatt w 3 <br />/b = Lb (1 — exp( —watt)) U ( ) <br />In this equation, L is the canopy projected leaf area <br />index estimated to be 2.0 (Schaeffer et al., 2000). The <br />quantity b was set equal to 0.0067ms -1/2. It is a <br />scaling coefficient for leaf boundary layer resistance <br />(Magnani et al., 1998). aatt is an attenuation coeffi- <br />cient for wind speed inside the canopy, w = 0.05 in <br />is a typical leaf width, and U the wind speed outside <br />the canopy (measured at 10m above the ground at <br />the mesquite site). The value for the wind attenuation <br />