The Priestley–Taylor (PT) approximation for computing evapotranspiration was initially developed for conditions of a horizontally uniform saturated surface sufficiently extended to obviate any significant advection of energy. Nevertheless, the PT approach has been effectively implemented within the framework of a thermal-based two-source model (TSM) of the surface energy balance, yielding reasonable latent heat flux estimates over a range in vegetative cover and climate conditions. In the TSM, however, the PT approach is applied only to the canopy component of the latent heat flux, which may behave more conservatively than the bulk (soil + canopy) system. The objective of this research is to investigate the response of the canopy and bulk PT parameters to varying leaf area index (LAI) and vapor pressure deficit (VPD) in both natural and agricultural vegetated systems, to better understand the utility and limitations of this approximation within the context of the TSM. Micrometeorological flux measurements collected at multiple sites under a wide range of atmospheric conditions were used to implement an optimization scheme, assessing the value of the PT parameter for best performance of the TSM. Overall, the findings suggest that within the context of the TSM, the optimal canopy PT coefficient for agricultural crops appears to have a fairly conservative value of ∼1.2 except when under very high vapor pressure deficit (VPD) conditions, when its value increases. For natural vegetation (primarily grasslands), the optimal canopy PT coefficient assumed lower values on average (∼0.9) and dropped even further at high values of VPD. This analysis provides some insight as to why the PT approach, initially developed for regional estimates of potential evapotranspiration, can be used successfully in the TSM scheme to yield reliable heat flux estimates over a variety of land cover types.
The two-source model (TSM) of Norman et al. (1995) is a land surface parameterization of the radiative and turbulent energy exchanges between the soil, vegetation, and lower atmosphere in which radiometric land surface temperature, derived from remote sensing images collected in the thermal wave band, serves as the key boundary condition. The TSM formulation is based on the two-source energy balance equation of
where the subscripts c and s denote the canopy and soil components, respectively; Rn is net radiation; G is soil heat flux; and H and λE are sensible and latent heat fluxes, respectively. This modeling framework follows earlier development of two-source schemes based on the Penman–Monteith big-leaf model (Monteith 1981) to deal with sparse canopy cover conditions (Shuttleworth and Wallace 1985; Shuttleworth and Gurney 1990). Given an estimate of fractional vegetation cover, the TSM partitions the observed surface temperature into soil and canopy contributions (illustrated in Fig. 1):
where TRAD(θ) is the surface radiometric temperature at look angle θ, f (θ) is fractional vegetation cover seen by the sensor, and TC and TS are the derived canopy and soil temperatures, respectively. The component (soil and vegetation) sensible heat fluxes (HS and HC) are then computed along the gradients in temperature, regulated by transport resistances (see Fig. 1). Extinction of net radiation within the canopy (Rnc) is approximated with an analytical formalism described by Campbell and Norman (1998), based primarily on leaf absorptivity and leaf area index (LAI), whereas G is parameterized as a fraction (0.31) of the net radiation above the soil surface (RnS), following Choudhury et al. (1994). The canopy transpiration component of the latent heat flux λEC is approximated using the Priestley–Taylor approach:
where αc is the PT parameter applied to the canopy (details provided later), Δ is the slope of saturation vapor pressure versus temperature, and γ is the psychrometric constant. Lastly, the soil evaporation term is computed as the residual of the overall energy balance equation [Eq. (1)].
The use of the PT approach within the TSM scheme to estimate canopy transpiration was motivated by its simplicity and the apparent robustness of its predictions. It provides information required to efficiently solve the system of equations defining the component soil and canopy heat fluxes and temperatures, which are critical to constraining fluxes from the two-source system. If these component temperatures are known a priori—for example, through direct measurements (impractical over wide spatial and temporal scales) or through dual-angle thermal infrared decomposition (e.g., Otterman et al. 1992; François 2002)—then the soil and canopy latent heat fluxes can be computed directly as a residual to the component energy budgets, eliminating the need for the PT approximation (Sanchez et al. 2008). However, if only the composite temperature is available (which is generally the case), then some initial assumption about either soil or canopy evaporation must be made to inform the partitioning of temperatures and heat fluxes between these two sources.
The performance of the TSM has been tested in both modes, that of using the PT approach and using observations of component temperatures derived from detailed infrared thermometer (IRT) observations (Kustas and Norman 1999; Sanchez et al. 2008) and from dual-angle measurements (Kustas and Norman 1997). In these tests, model fluxes from the PT mode provided similar or better agreement with observations than did model fluxes using observed values of TS and TC. The generally better performance using the PT approach was most pronounced in cases in which dual-angle measurements were employed to derive TS and TC, which often has much greater uncertainty in estimating reliable component temperatures. Within the context of the TSM, the PT mode appears to be a simple and economical approach, and it provides robust assessments of component and system evaporative fluxes over a wide range of vegetation cover and climatic conditions (e.g., Kustas and Norman 1999, 2000a,b; French et al. 2003, 2005; Li et al. 2006).
Still, legitimate questions can be raised regarding the validity of implementing the PT approximation in the TSM, as a parameterization for the canopy transpiration flux component (e.g., Sanchez et al. 2008). Studies have indicated that the PT coefficient associated with the combined soil–plant–atmosphere system is not constant, but it does show dependency on LAI, vapor pressure deficit (VPD), and soil moisture (e.g., Jury and Tanner 1976; Flint and Childs 1991; Pereira 2004; Diaz-Espejo et al. 2005; Baldocchi and Xu 2007). The TSM, however, applies the PT approach only to the canopy component of the latent heat flux, which may behave more conservatively than that associated with the bulk system. The objective of this research was to investigate the behavior of the canopy PT parameter diagnosed by the TSM, and its utility and limitations over a range in vegetation cover and moisture conditions.
2. Scientific background
a. Priestley–Taylor approach applied to the bulk system
The Priestley and Taylor (1972) approach is a simplification of Penman’s formulation of evapotranspiration ET (Penman 1948), assuming that the equilibrium term λEeq (Slatyer and McIlroy 1961) is significantly larger than the aerodynamic term λEa. On the basis of this assumption, the Priestley–Taylor formulation for ET is expressed as
where α is the so-called Priestley–Taylor coefficient.
Numerous studies have found that the value of α can be quite variable (e.g., Jury and Tanner 1976; Flint and Childs 1991; Castellvi et al. 2001; Pereira 2004; Diaz-Espejo et al. 2005; Baldocchi and Xu 2007). When the soil has ample water, the latent heat flux is determined by the atmospheric evaporative demand, and α increases with increasing λEa/λE—that is, under conditions of strong horizontal advection of sensible heat flux (Priestley and Taylor 1972; Jury and Tanner 1975; Tanner and Jury 1976). This has been experimentally documented by Singh and Taillefer (1986), Diaz-Espejo et al. (2005), and Li and Yu (2007). Once soil moisture decreases, and the soil and/or vegetation can no longer transpire at the potential rate, the surface resistance to λE increases and α decreases (e.g., De Bruin 1983; Flint and Childs 1991; Raupach 2000).
Each of the previously mentioned studies treated the soil–biosphere–atmosphere as a bulk system, based on the big-leaf assumption, which may be reasonable for dense vegetation, but is less appropriate for partial canopies. Soil evaporation becomes rapidly limited by the soil hydraulic conductivity when the water content of the surface layer (∼0–5 cm) drops below a certain threshold, whereas transpiration continues at an energy-limited rate until the moisture content over a large part of the root zone drops below a critical level (Tanner and Jury 1976). For vegetated surfaces, 50%–80% of the available soil water in the plant root zone can be extracted at the potential rate. For bare soil, however, potential evaporation is limited to only ∼40% of the available soil water in the near-surface (0–5 cm) layer (Flint and Childs 1991). With α for the bulk system effectively representing a combination of soil evaporation and canopy transpiration, only under the wettest soil moisture conditions can α be expected to remain at ∼1.26 for partial canopy cover conditions. At low vegetation cover fraction (or LAI), it is more likely that α will be less than 1.26, because the soil evaporation component of ET rapidly decreases as the soil surface dries out (Stannard 1993). Conversely, α has been reported to exceed 1.26 and reach values of ∼2, particularly under well-watered closed-canopy conditions (e.g., Baldocchi 1994; Baldocchi et al. 1997).
There are both theoretical and experimental studies indicating that the value of α varies significantly with LAI, VPD, and soil moisture (e.g., Jury and Tanner 1976; Flint and Childs 1991; Baldocchi 1994; Baldocchi et al. 1997; Pereira 2004; Diaz-Espejo et al. 2005; Baldocchi and Xu 2007). However, a recent study by Pereira et al. (2007) used the Priestley–Taylor approach only for the canopy elements (with the traditional value of αc = 1.26) for estimating transpiration from well-watered isolated fruit trees and obtained excellent agreement with sap-flow measurements over a wide range of LAI. This implies that applying the Priestley–Taylor method to the canopy component may result in a coefficient that is more conservative than α for the whole (bulk) soil–biosphere–atmosphere system αbulk when the canopy is unstressed. This study investigates the behavior of the value of α when applied to solely the vegetation component αc, compared to αbulk.
b. Two-source formulations for α
1) Tanner and Jury
One of the first attempts to account for the effects of vegetation amount on evaporation from the bulk system by separating evaporation and transpiration components was formulated by Ritchie (1972). Tanner and Jury (1976) adapted Ritchie’s methodology by applying the Priestley–Taylor approach for the soil and vegetated canopy components. They assumed that under energy-limited conditions (i.e., ample soil water), the soil evaporation λEs and transpiration λEc can be expressed as
where τ, a canopy transmission factor representing the ratio of net all-wave radiation at the soil surface (Rns) to that above the canopy (Rn), is given by
In Eqs. (5) and (6), β is a semiempirical coefficient dependent on leaf angle distribution and soil and canopy temperatures (Kustas and Norman 1999), and αs and αc are modified PT coefficients applied to λEs and λEc, respectively, and are given by
where α = 1.26 for energy limited (wet) conditions with a full canopy cover, τ0 is a critical value of τ below which the canopy is sufficiently closed that soil evaporation approaches equilibrium (αs ≈ 1 for wet soil conditions). Equation (7) is applied under energy-limited conditions. The choice of τ0 may vary between 0.2 and 0.5, and its exact value was found not to be crucial (Tanner and Jury 1976). This approach was recently applied to cornfields in Wisconsin (Morgan et al. 2003) and soybeans in Iowa (Sauer et al. 2007).
Stannard (1993) suggested that for estimating λE at the soil–biosphere–atmosphere level, the “system” or “bulk” α should be modified to account for separate contributions from the canopy and soil:
where α = 1.26, C8 and C9 are parameters determined using nonlinear regression, and is the amount of water evaporated since the last rainfall. For details about the individual terms in Eq. (8), see Stannard (1993). In this formulation, either a large LAI, or a recent rainfall, are sufficient to drive αbulk toward its maximum value of 1.26, whereas either a small LAI or an extended dry spell results in a minimum value of αbulk = (α − C8) ≈ 0.35. Note that in determining values of C8 and C9, Stannard (1993) used data from a sparsely vegetated site with a sandy soil overlying a shallow water table (0.8–1.6-m depth). Therefore, although the soil surface may have dried out in the absence of precipitation, the vegetation had ample moisture at all times.
3) The two-source model
Detailed descriptions of the principles and formulation of the two-source energy balance model are given by Norman et al. (1995) and Kustas and Norman (1999, 2000a). In the model, the transpiration component of the latent heat flux is approximated using the Priestley–Taylor equation [Eq. (3)], in which Rnc is the amount of net radiation intercepted by the canopy, given in the original form of the TSM (Norman et al. 1995) as
Given an initial estimate of λEc using αc = 1.26 (assuming potential or unstressed transpiration), the TSM computes λEs as the residual of the overall energy balance. If λEs is found to be less than 0, suggesting condensation of water vapor onto the soil, which is unlikely to occur during daytime dry convective conditions, then an additional iterative stage is executed, during which αc is gradually reduced until λEs ≥ 0, reflecting some degree of canopy stress.
Although current versions of the TSM use a more physically based scheme for partitioning net radiation, accounting for differences in longwave and shortwave radiation attenuation in the canopy layer (Kustas and Norman 2000a), the simplified form in Eq. (10) has been used in this study to develop a conceptual/analytical relationship between αbulk and αs and αc. For this reason, λEs is also expressed here in terms of the Priestley–Taylor equation
where G is estimated as a fraction of Rns, that is, G = cRns = cRn exp(−βLAI). Note that this differs from the TSM algorithm, which solves for λEs as a residual from the soil surface energy balance equation, thus not requiring an estimate of αs for the model operation.
The original Priestley–Taylor expression for the total λE is
However, λE can also be expressed as the sum of the λEc and λEs as follows:
As in Stannard’s formulation, the dependency of αbulk on LAI in Eq. (14) is explicit.
Note that although the TSM has some ability to adjust αc in response to perceived canopy stress, it is unlikely that this adjustment will result in the “optimal” or the most appropriate αc value in every case for the following reasons: 1) the model has an iterative mechanism to reduce αc but not to increase it, which may be required under strong advective conditions; 2) some degree of canopy stress (αc < 1.26) could be realized before soil evaporation goes to zero (λEs > 0), but there is no mechanism in the model to capture this condition.
a. Data description
Micrometeorological measurements collected at multiple sites under a wide range of atmospheric conditions were used in this study. Data from agricultural fields were collected in Maryland [Optimizing Production Inputs for Economic and Environmental Enhancement (OPE3)], Iowa [Soil Moisture Experiment 2002 (SMEX02)], Texas (Bushland) and Arizona (Maricopa), representing gradually changing climatic zones from the temperate climate in Maryland to the Sonoran Desert in Arizona. The crops grown at these sites were corn, soybean, alfalfa, sorghum, wheat, and cotton. Data were also collected at rangeland sites (both grassland and shrubland) in Oklahoma [Southern Great Plains 1997 (SGP97)], Arizona [Audubon, Monsoon ’90, and Soil Moisture Experiment 2004 (SMEX04)], Mississippi (Goodwin Creek), and Montana (Fort Peck). The exact location of the sites, the dates of data collection, and mean LAI and VPD at each site are listed in Table 1; a further description of the site conditions and instrumentation configurations can be found in the references listed in Table 1.
The sensible and latent heat fluxes were measured using several micrometeorological techniques. In Maricopa, both Bowen ratio and eddy covariance (EC) methods were used (Kustas et al. 1989; Kustas 1990). In Bushland, λE was measured by large weighing lysimeters (Howell et al. 1995a), and H was determined as the residual in the energy balance equation. In Monsoon ’90, both EC and flux variance approaches were used (Kustas and Goodrich 1994; Stannard et al. 1994). For all other sites, H and λE were measured only using the EC method. For comparison with fluxes computed by the TSM (which assumes closure of the energy budget), a closure correction was applied to all EC flux datasets by assigning the residual flux to λE (the residual method; Li et al. 2005). At all sites, measurements of net radiation, soil heat flux, and radiometric surface temperature were also collected, along with meteorological data (e.g., air temperature, VPD, and wind speed).
At these sites, LAI was measured using several different methods. In OPE3 (Sanchez et al. 2008), SMEX04 (Li et al. 2008), SMEX02 (Anderson et al. 2005), and SGP97 (Norman et al. 2000), LAI data were acquired with LAI-2000 plant canopy analyzer (LICOR Inc.; Company and trade names are given for the benefit of the reader and imply no endorsement of the product by the USDA-ARS). During Monsoon ’90, LAI was estimated using the line-intercept transect method (Canfield 1941) as detailed by (Weltz et al. 1994). LAI data from Maricopa were estimated by agronomic measurements described by Kustas et al. (1989) and Kustas (1990). The LAI estimates from the Bushland sites were obtained by descriptive plant sampling (1.5 m2 samples) in the 5-ha fields surrounding each lysimeter. In Audubon, Goodwin Creek, and Fort Peck, LAI was estimated from onsite normalized difference vegetation index (NDVI) observations (Wilson and Meyers 2007). For all sites, linear interpolation was used to estimate LAI between measurement dates.
For many sites, with the exception of OPE3, Maricopa, Audubon, Goodwin Creek, and Fort Peck, data from multiple flux towers (multiple lysimeters in Bushland) were used, and each flux tower was considered as a separate dataset. For sites with multiple years of flux acquisition, each year was considered as a distinct “flux site dataset” in the analysis. In total, 21 different flux sites were included in the study (see Table 1).
Half-hourly data at midday (1000–1400 local time) were analyzed, corresponding with periods of peak flux and typical times of daytime satellite data acquisition when the TSM is generally applied. Data were further screened to exclude cloudy conditions, as suggested by Black (1979) and Flint and Childs (1991), because under such conditions, even relatively dry soils can supply enough water to meet potential evapotranspiration.
b. Data analysis
1) Scheme for optimizing αc
Direct and reliable measurements of vegetation and soil components of the surface energy balance are very difficult to obtain at canopy and field scale. Therefore, an approach was employed to estimate optimal values of αc using the representation of flux partitioning inherent in the TSM. The optimization scheme determined the optimal αc for each data point by iteratively running the TSM over a range in initial αc values, varying between 0.5 and 2.5 in increments of 0.1. For each time step at each flux site, the TSM was executed 21 times with the various αc input values; after each execution, the modeled λE was compared to the measured flux. The optimal αc was defined as the αc value obtained from the TSM computation that resulted in the minimal error with the λE measurements. In addition, a bulk observed value of α was computed from the flux measurements using Eq. (12).
An alternative approach would have been to optimize αc on average over the course of each day. Although this would reduce noise in the derived values of optimal αc (in comparison with the hourly optimization employed here), it would significantly reduce the range in VPD represented in the dataset, concealing the high-VPD responses discussed later.
2) Independent variables and data binning
To emphasize trends in αc and αbulk with varying environmental factors, the resulting dataset for each flux site was averaged over bins in LAI and VPD. LAI and VPD were selected as driving variables for investigation because both potentially affect αbulk but are independent of the Priestley–Taylor formulation [cf. Eq. (12)]. Note that soil moisture content was not measured at all sites, and therefore it could not be used as an independent variable in this study.
Data were binned into LAI ranges in increments of 1 (0–1, 1–2 … 4–5). Within each bin, averages, standard deviations (STD), and standard errors (SE = STD/n) in LAI, αc, and αbulk were computed. In the same way, data were binned into VPD ranges in 1-kPa increments. The standard errors were computed to account for differing number of observations available for the various sites in the different bins. Linear regressions describing trends in αc and αbulk versus LAI and VPD were weighted by the standard errors for each binned data point (Bevington 1969).
In this study, data binning was required for trend analysis to suppress noise in optimized α while maintaining a wide range in the driving variables, VPD and LAI. Noise in the hourly optimized α values reflect uncertainties in the observations, enhanced by the variety of observational methods employed in these various studies, and the fact that these experiments were not originally designed to specifically address the topic of this study. Potentially, most problematic are the LAI and surface radiometric temperature data—key inputs to the TSM in terms of model sensitivity. An additional and significant source of uncertainty is the measurement of latent heat flux from the various micrometeorological methods and issues related to forcing energy balance closure (e.g., Twine et al. 2000). Without binning, the primary trends discussed here are obscured.
In Fig. 2 the range in LAI and VPD sampled by the observational dataset used in the current study is illustrated. Each point represents the mean LAI for the corresponding 1-unit increments described earlier and the corresponding mean VPD value. The bars represent the STD values in LAI and VPD. Several sites, where the data used were collected during multiple years and/or the data cover seasonal growth cycle, will have more than one LAI/VPD data point. Note that this study did not include data from sites with both high LAI and high VPD. In nature, such conditions are not likely to exist, because high VPD usually represents dry conditions where natural vegetation tends to be sparse. Sites with both high LAI and persistently high VPD are typically found only in small riparian areas, where micrometeorological measurements are difficult (or impossible) to conduct, or where mature crops are being artificially irrigated. In the current study, high VPD values (greater than 3 kPa) were only measured in sites with LAI < 0.5 (most with natural vegetation).
4. Results and discussion
a. Effect of LAI on α
The bulk Priestley–Taylor coefficient does not consider the effect of the fractional vegetation cover on the soil–canopy partitioning of latent heat flux, yet this is the form in which the PT method is typically applied (e.g., Kustas et al. 1996). Figure 3a shows a set of points representing the distribution of observed αbulk (derived directly from the flux measurements) versus the observed LAI for the combined set of flux sites. Each point represents an average for each flux site in the respective LAI bin, with X-axis “error bars” representing the STD of LAI (to indicate the LAI variability within each LAI bin) and Y-axis error bars representing the SE in αbulk (to normalize the standard deviation by the number of observations), which tends to be quite small. Also illustrated in Fig. 3a is the theoretical dependency of αbulk on LAI as described by a set of analytical curves, computed using Eq. (14). These curves assume four different surface soil moisture conditions, as represented by the different values of αs. A value of αs = 1.3 represents a wet bare soil surface, whereas αs = 0 represents a completely dry soil surface (i.e., no evaporation from the soil), and αs = 0.5 and 1.0, represent intermediate conditions. The effect of vegetation cover on αs at nonzero LAI for the wet soil case (energy-limited case) was accounted for using Eq. (7) with τ0 = 0.3 and β = 0.5. For the other cases where αs ≤ 1, αs was assumed to be constant, which is a reasonable assumption when soil evaporation reaches stage two drying (Hillel 1971). Four of the curves use αc = 1.3, assuming the vegetation is transpiring at the potential rate (unstressed). A fifth curve represents a dry soil surface (αs = 0) and αc of 0.9, which is the average optimized αc value found for the natural vegetation dataset (see later). Clearly, the drier the soil, the stronger the effect of LAI on the value of αbulk, especially for LAI < 2.
In comparing the analytical curves of αbulk with the data points in Fig. 3a, most of the agricultural crop data points at low LAI are clustered around the unstressed canopy curves with αs = 0.5 and 1.0, indicating intermediate surface soil moisture conditions at these sites (at high LAI, the curves are virtually indistinguishable). The αbulk values for the majority of the sites with natural vegetation at higher cover (i.e., LAI > 2) are clustered around the curve representing dry soil and vegetation transpiring at less than the potential rate. These findings agree well with the fact that agricultural crops tend to be grown in areas with adequate moisture achieved either by natural rainfall or by irrigation, whereas the water supply for natural vegetation is largely dependent on precipitation, which can be highly variable in space and time and is often inadequate to sustain potential ET.
Figure 3b shows the mean final αc for each LAI bin, as derived by the TSM optimization scheme. In contrast to αbulk, which showed dependency on LAI (p < 0.01), αc in this dataset is independent of LAI (p > 0.05). The large scatter in the optimal αc at low LAI is attributed to the fact that αc is not well constrained in the TSM optimization procedure using flux data over sparse vegetation. Under low vegetation cover conditions, Rnc is relatively small, thus large changes in αc are required to effect a significant change in λE. However, in these cases, the contribution of canopy transpiration to the overall latent heat flux is small, and errors in computations of the transpiration flux will not induce large errors to the overall flux estimates.
For the agricultural crop dataset used in this study, the average value of optimal αc was 1.2, close to the standard value of αbulk = 1.3 that is commonly used in most bulk applications of the Priestley–Taylor approximation. This suggests that in most cases, the agricultural crops studied here were not under stress and were transpiring at close to the potential rate. The average optimal αc value obtained for the natural vegetation dataset was 0.9. This lower value may be because most of the natural vegetation sites studied here were located in semiarid and subhumid climates, where vegetation growth is often water limited. Under such conditions, natural vegetation can adapt to frequent moisture deficits by controlling the rate of transpiration. Consequently, a lower value of αc would be expected for these natural vegetation sites, unless sufficient root-zone moisture is always available (e.g., Stannard 1993). The average value of optimal αc for the combined dataset (including both crops and natural vegetation) is 1.1.
The general response of αbulk and αc to LAI depicted in Figs. 3a and 3b was also reproduced with a multilayer soil–plant–atmosphere model, Cupid. Detailed descriptions of Cupid can be found in Norman and Campbell (1983), Norman and Arkebauer (1991), and Kustas et al. (2004). In brief, Cupid computes both radiative and turbulent exchanges across the soil–canopy–air interface while considering a wide variety of physiological and environmental processes simultaneously.
Kustas et al. (2004, 2007) describe a numerical experiment in which the Cupid model was first validated with the Monsoon ’90 flux dataset, and then a set of hypothetical simulations was performed assuming a range of soil moisture and LAI conditions. Figure 4 shows mean values of αc, αs and αbulk derived from Cupid simulations of λEc, λEs, and λE using Eqs. (3), (11), and (12), respectively, assuming unstressed canopy conditions. Both αs and αbulk increase with increasing LAI, whereas αc remains relatively constant (see also Kustas and Anderson 2009). The Cupid-derived average value of αc of ∼1.4 is slightly higher than 1.26 and reflects the fact that the unstressed simulations were run using local meteorological forcing under arid (advective) conditions (see Kustas et al. 2004). The fact that the detailed Cupid model reproduces the PT coefficient behavior observed with the simple TSM optimization scheme provides further evidence that αc will be relatively conservative with LAI when the canopy is unstressed.
b. Effect of VPD on α
1) Agricultural crops
The dependency of αc and αbulk on VPD for flux sites with agricultural crops is illustrated in Figs. 5a and 5b, respectively. As in Fig. 3, each point represents an average for each site in the respective VPD bin, with X-axis error bars representing the STD of VPD and Y-axis error bars representing the SE in α. The trend line (solid line) in Fig. 5a is derived from a statistically weighted linear regression (r 2 = 0.68) and indicates a dependency of αc on VPD. From a modeling standpoint, however, significant improvement in TSM latent heat estimates resulting from replacing the nominal value of αc = 1.3 with the optimal αc is obtained only under the high VPD range, as detailed in section 4c below. This suggests that a modification to αc as utilized in the TSM for agricultural crops is required only at high VPD, as indicated by the piecewise linear relationship (i.e., dashed line) in Fig. 6a. Additional data, particularly collected under strong advective conditions and high VPD, will be required to determine a robust quantitative relationship between αc and VPD for use in the TSM.
Interestingly, αbulk in Fig. 5b shows little dependency on VPD, with the statistically weighted linear regression (not shown) yielding a slope not significantly different than 0 (p > 0.05). This should not be interpreted, however, as an indication that αbulk is conservative over this wide range in VPD. As mentioned in section 3b(2), the flux dataset used here did not include data from sites with both high LAI and high VPD (Fig. 2), as might be encountered in smaller irrigated agricultural systems in the western United States. For well-watered agricultural crops, αbulk is expected to behave similarly to αc in Fig. 5a (i.e., increasing at high VPD). This is represented by the piecewise linear curves in Fig. 5b derived from the analytical model for αbulk in Eq. (14), for a range of LAI cases (5, 2, 1, and 0.5), with αs = 0.5 and αc from Fig. 5a. If the dataset contained data from mature irrigated crops in arid climates, then the upper curves would be populated. The trends at the higher LAI values represent a theoretical upper bound in αbulk because sheltering effects of dense vegetation are likely to cause dampening of the αc response to increasing VPD, analogous to the effects on αs described by Tanner and Jury (1976). These curves reveal how the increase in canopy transpiration but decrease in soil evaporation can cancel out under low canopy cover conditions, such that a relatively small change in the overall αbulk is computed.
2) Natural vegetation
For the natural vegetation sites, both αc (Fig. 6a) and αbulk (Fig. 6b) appear to decrease with increasing VPD, although the scatter is large. The response of αc to VPD may relate to physiological characteristics of natural vegetation growing in arid and semiarid environments. Although an increase in VPD enhances transpiration by producing stronger humidity gradient between the leaf and the atmosphere, it also initiates a negative feedback on stomatal conductance, which leads to a reduction in transpiration (Baldocchi and Xu 2007). Different adaptation mechanisms in natural vegetation lead to a stronger response of stomatal resistance to high VPD to preserve water (e.g., Baldocchi and Xu 2007; Galmes et al. 2007), causing a decrease in transpiration and thus a decrease in αc.
The decrease in αbulk with VPD is in part due to the behavior of αc, but it also reflects the dependence of αs on VPD, which is expected to be significant in rain-fed systems. High VPD will induce high soil evaporation rates and a rapid decrease in soil moisture, moving quickly from the first stage to the second stage of evaporation (Hillel 1971). The decrease in soil evaporation affects the overall decrease in latent heat flux, reflected in the decrease in αbulk with VPD. The lowest group of points, in the VPD range between 3 and 6 kPa, falling below the αbulk curve with αs = 0 indicates that both the surface and profile soil moisture conditions are likely affecting the overall decrease in αbulk. This is in contrast to the agricultural crops, which are likely to continue transpiring at the potential rate up to high values of VPD as long as ample water exists, as illustrated by the αbulk curves (dashed lines) in Fig. 6b, derived from analysis of the α–VPD relationship for the agricultural crops.
c. Errors in TSM using αc = 1.3
For operational applications of the TSM (e.g., Anderson et al. 2007a,b), it is important to understand the errors in λE incurred by using the conventional initialization value of αc = 1.3 for various land use and climatic conditions. In Fig. 7, the mean-absolute-error (MAE) statistic between modeled and measured λE recommended by Willmott and Matsurura (2005) is used to quantify the performance of the TSM as applied to the agricultural and natural vegetation datasets, segregated by VPD (above and below 4 kPa) and by LAI bins. For each case, the TSM was initialized with αc = 1.3, and the model was allowed to downregulate the PT coefficient in response to perceived canopy stress, as per the standard mode of model operation.
For the agricultural vegetation dataset under low-to-moderate atmospheric demand (VPD < 4 kPa), MAE values of 25 W m−2 are obtained on average over all LAI bins using αc = 1.3, whereas the natural vegetation sites yielded a higher average error of 38 W m−2. If a lower value of αc = 1.0 is used to initialize runs for the natural vegetation sites, as suggested by optimization results from section 4a, then this error reduces to 24 W m−2. In either case, the errors accrued by using a fixed nominal value of αc are comparable or superior to model–measurement errors encountered in many other studies applying remote sensing–based energy balance methods (Kalma et al. 2008; Kustas and Norman 2000b), which are typically on the order of 50 W m−2. Note that errors using αc = 1.3 tend to decrease with increasing LAI.
Under high VPD conditions, the choice of nominal input αc is more critical for obtaining good flux estimates from the TSM over both agricultural crops and natural vegetation. For these cases, using an input value of αc = 1.3 produces MAE values of 54 and 67 W m−2 for agricultural and natural vegetation datasets, respectively—somewhat higher than the 50 W m−2 errors typically reported for surface energy balance. This suggests that a VPD-dependent function of αc may be beneficial in reducing model errors under conditions of very high atmospheric demand. For natural vegetation, a lower value of input αc may be warranted, reflecting the relatively conservative water-use tendencies of undomesticated plants. Additional data under high VPD conditions (>4 kPa) is required to thoroughly evaluate this hypothesis.
This study examined the utility of the PT approach in providing an initial value for canopy transpiration, required for solving the surface energy balance using the TSM scheme. In the TSM, the PT approximation is applied only to the canopy component of the system latent heat flux, reflecting the ability of the vegetation to transpire, and hence αc is expected to be relatively conservative with respect to changes in LAI. The PT coefficient for the bulk (soil + canopy) system αbulk shows a different behavior, typically decreasing at low vegetation cover, density, or LAI with the increasing influence of soil evaporation component of ET.
For agricultural crops under low-to-moderate VPD conditions (0 < VPD < 4), initializing the TSM with a constant αc = 1.3 did not significantly degrade model performance in estimating λE compared to using the optimal αc (Fig. 7). For the natural vegetation, data using αc = 1.3 yielded somewhat larger errors, which might be reduced if the TSM were initialized with a lower value of αc. Nevertheless, in either case, the model errors for runs using a nominal value of αc = 1.3 are comparable to model-measurement errors encountered in many other studies applying remote sensing–based energy balance methods (Kalma et al. 2008). This is a benefit, because in practice, the optimal value of αc will not be known a priori.
On the other hand, under strong advective conditions with high VPD (>4 kPa) and at sites with sufficient available soil moisture (in this study represented only by agricultural crops), the estimated value of αc increased significantly with VPD (cf. Fig. 5a). Under such conditions, agricultural crops will increase transpiration to meet atmospheric demand as long as an ample root-zone water supply exists. For these cases, significantly larger errors in modeled λE resulted for both vegetation types (particularly agricultural crops) when using a nominal αc = 1.3 as compared to the optimal αc (Fig. 7). For sites with natural vegetation, αc did not increase with increasing VPD (Fig. 6a), as was observed for crop sites (Fig. 5a). This may reflect limitations in available soil moisture at these sites, which tended to be in semiarid and subhumid climates, and a greater tendency for undomesticated vegetation to reduce water loss/transpiration under high VPD conditions.
Overall, the findings suggest that, except under very high VPD conditions, initialization of the TSM with the conventional value of αc = 1.3 does not significantly degrade the model performance, especially for agricultural crops. This analysis provides some insight as to why the PT approach, initially developed for regional estimates of potential evapotranspiration, can be used successfully in the TSM scheme, yielding reliable heat flux estimates over typical agricultural crops, as well as grassland- and shrubland-type natural vegetation. Future work will focus on evaluating the PT approach under very high VPD conditions and for deciduous and coniferous forests.
This research was supported by Vaadia-BARD Postdoctoral Fellowship Award FI-371-2005 from BARD, the United States–Israel Binational Agricultural Research and Development Fund.
* Current affiliation: Gilat Research Center, Agricultural Research Organization of Israel, D. N. Negev, Israel.
Corresponding author address: Nurit Agam, Gilat Research Center, Agricultural Research Organization of Israel, D. N. Negev 85280, Israel. Email: email@example.com