Introduction
The Amazon rainforest is one of the world's most extensive natural ecosystems, influencing the Earth's water, energy, and carbon cycles (Malhi, 2012), and is also a major source of global terrestrial evapotranspiration () or latent heat flux () (Costa et al., 2010; Harper et al., 2014). An intensification of the Amazon hydrological cycle was observed in the past 2 decades (Cox et al., 2000; Huntingford et al., 2008; Gloor et al., 2013). Recent Amazonian droughts have gained particular attention due to the sensitivity of the tropical forest to climate change (Hilker et al., 2014). If persistent precipitation extremes become more prevalent (Hilker et al., 2014), the Amazon rainforest may increasingly become a net source of carbon as a result of both the suppression of net biome exchange by drought and carbon emissions from fires (Gatti et al., 2014). Changes in land cover due to conversion of tropical forest to pastures significantly alter the energy partitioning by decreasing and increasing sensible heat fluxes () over pasture sites (e.g., Priante Filho et al., 2004). This will ultimately lead to severe consequences for the water balance in the region, with changes to river discharge already observed in some parts of the basin (Davidson et al., 2012). Evaluating the response to changing climate and land use in the Amazon Basin is critical to understanding the stability of the tropics within the Earth system (Lawrence and Vandecar, 2015). The control of can be viewed as complex supply–demand interactions, where net radiation and soil moisture represent the supply and the atmospheric vapor pressure deficit represents the demand. This supply–demand interaction accelerates the biophysical feedbacks in , and understanding these biophysical feedbacks is necessary to assess the terrestrial biosphere response to water availability. Therefore, quantifying the critical role of biophysical variables on will add substantial insight to assessments of the resilience of the Amazon Basin under global change.
The aerodynamic and canopy conductances ( and hereafter) (unit m s) are the two most important biophysical variables regulating the evaporation () and transpiration ( flux components of (Monteith and Unsworth, 2008; Dolman et al., 2014; Raupach, 1995; Colaizzi et al., 2012; Bonan et al., 2014). While controls the bulk aerodynamic transfer of energy and water through the near-surface boundary layer, represents the restriction on water vapor flow through the aggregated conductance from stomata of the leaves, in the case of a vegetated surface. In the case of partial vegetation cover, also includes the soil surface conductance for evaporation. At a small / ratio, the vapor pressure deficit close to the canopy source/sink height () approximates the atmospheric vapor pressure deficit () due to aerodynamic mixing and/or low transpiration. This results in a strong canopy–atmosphere coupling, and such conditions are prevalent under soil moisture deficits, which prevails under conditions of soil moisture deficit. By contrast, a large / ratio influences the gradients of vapor pressure deficit just above the canopy, such that tends towards zero and thus remains different from (Jarvis and McNaughton, 1986). This situation reflects a weak canopy–atmosphere coupling, and such situations prevail under predominantly wet conditions and/or poor aerodynamic mixing due to wetness-induced low aerodynamic roughness. The Penman–Monteith (PM) equation is a physically based scheme for quantifying such biophysical controls on canopy-scale and from terrestrial ecosystems, treating the vegetation canopy as “big-leaf” (Monteith, 1965, 1981). Despite its development based on biophysical principles controlling water vapor exchange, quantifying the and controls on through the PM equation suffers from the continued longstanding uncertainty in the aggregated stomatal and aerodynamic behavior within the soil–plant–atmosphere continuum (Matheny et al., 2014; Prihodko et al., 2008).
One of the major sources of uncertainties in modeling is associated with the empirical (and uncertain) parameterizations of near-surface boundary-layer dynamics, which is invariably confounded by space–time variability in atmospheric stability (van der Tol et al., 2009; Shuttleworth, 1989; Gibson et al., 2011). For example, the Monin–Obukhov similarity theory (MOST) used for modeling appears to be only valid over uniform, extensive, and flat surfaces (Monteith and Unsworth, 2008; van der Tol et al., 2009; Holwerda et al., 2012), and its application to complex “real” canopy systems is problematic due to chaotic interactions between turbulence, canopy roughness, and topography (Raupach and Finnigan, 1995; Shuttleworth, 2007; Holwerda et al., 2012). Similarly, varies in space and time due to variations in plant species, photosynthetic capacity, soil moisture variability, and environmental drivers (Monteith and Unsworth, 2008; van der Tol et al., 2009). Despite the existence of several semi-mechanistic and empirical parameterizations for (e.g., Ball et al., 1987; Leuning, 1995; Tuzet et al., 2003; Medlyn et al., 2011), the adaptive tendencies of plant canopies severely compromise the efficacy of such approaches (Matheny et al., 2014), limiting their applicability over most landscapes. Thus, debate on the most appropriate model for the canopy conductance has endured for decades.
Previous studies in the Amazon Basin focused on developing an observational understanding of the biogeochemical cycling of energy, water, carbon, trace gases, and aerosols in Amazonia (Andreae et al., 2002; Malhi et al., 2002; da Rocha et al., 2009), model-based understanding of surface ecophysiological behavior and seasonality of (Baker et al., 2013; Christoffersen et al., 2014), modeling the environmental controls on (Hasler and Avissar, 2007; Costa et al., 2010), understanding the seasonality of photosynthesis and of (da Rocha et al., 2004; Restrepo-Coupe et al., 2013), and the impact of land use on hydrometeorology (Roy and Avissar, 2002; von Randow et al., 2012). However, the combination of climatic and ecohydrological disturbances will significantly affect stomatal functioning, the partitioning of , and carbon–water–climate interactions of tropical vegetation (Cox et al., 2000; Mercado et al., 2009). Hence, investigation of the effects of drought and land cover changes on conductances, and , is a topic requiring urgent attention (Blyth et al., 2010) both because of the cursory way it is handled in the current generation of parametric models (Matheny et al., 2014) and because of the centrality of and in controlling modeled flux behaviors (Villagarcía et al., 2010). The persistent risk of deforestation is likely to alter the radiation interception, surface temperature, surface moisture, associated meteorological conditions, and vegetation biophysical states of different plant functional types (PFTs). Conversion from forest to pasture is expected to change the / ratio of these ecosystems and impact the evapotranspiration components. Besides inverting the PM equation using field measurements of , to date either photosynthesis-dependent modeling or leaf-scale experiments were performed to directly quantify (Ball et al., 1987; Meinzer et al., 1993, 1997; Monteith, 1995; Jones, 1998; Motzer et al., 2005). However, an analytical or physical retrieval for and is required not only to better understand the role of the canopy in regulating evaporation and transpiration, but also to enable our capability to characterize the conductances using remote observations, across large spatial domains where in situ observations are not available. This paper aims to leverage this emerging opportunity by exploring data from the Large-scale Biosphere-Atmosphere Experiment in Amazonia (LBA) eddy covariance (EC) observations (e.g., de Goncalves et al., 2013; Restrepo-Coupe et al., 2013) using a novel analytical modeling technique, the Surface Temperature Initiated Closure (STIC) (Mallick et al., 2014, 2015), in order to quantify the biophysical control on and over several representative PFTs of the Amazon Basin.
STIC (STIC1.0 and STIC1.1) provides a unique framework for simultaneously estimating and , surface energy balance fluxes, and and . It is based on finding analytical solutions for and by physically integrating radiometric surface temperature ( information (along with radiative fluxes and meteorological variables) into the PM model (Mallick et al., 2014, 2015). The direct estimates of canopy-scale conductances and obtained through STIC are independent of any land surface parameterization. This contrasts with the multi-layer canopy models that explicitly parameterize the leaf-scale conductances and perform bottom-up scaling to derive the canopy-scale conductances (Baldocchi et al., 2002; Drewry et al., 2010). A primary advantage of the approach on which STIC is based is the ability to directly utilize remotely sensed to estimate , thereby providing a capability to estimate over large spatial scales using a remotely sensed variable that is central to many ongoing and upcoming missions. This study presents a detailed examination of the performance of STIC to better understand land–atmosphere interactions in one of the most critical global ecosystems and addresses the following science questions and objectives.
-
How realistic are canopy-scale conductances when estimated analytically (or non-parametrically) without involving any empirical leaf-scale parameterization?
-
What are the controls of canopy-scale and on evaporation and transpiration in the Amazon Basin, as evaluated using STIC?
-
How do the STIC-based canopy-scale conductances compare with known environmental constraints?
-
Is the biophysical response of consistent with the leaf-scale theory (Jarvis and McNaughton, 1986; McNaughton and Jarvis, 1991; Monteith, 1995)?
Methodology
Theory
The retrievals of , , and are based on finding a “closure” of the PM equation (Eq. 1 below) using the STIC framework (Fig. A1) (Mallick et al., 2015). STIC is a physically based single-source surface energy balance scheme that includes internally consistent estimation of and (Mallick et al., 2014, 2015). Originally designed for application to thermal remote sensing data from Earth observation sensors, the STIC framework exploits observations of radiative () and environmental variables, including net radiation (), ground heat flux (), air temperature (), relative humidity (), or vapor pressure () at a reference level above the surface.
Variables and symbols and their description used in the present study.
Variables | Description |
---|---|
and | |
symbols | |
Evapotranspiration (evaporation transpiration) as latent heat flux (W m) | |
Sensible heat flux (W m) | |
Net radiation (W m) | |
Ground heat flux (W m) | |
Net available energy (W m) | |
Air temperature (C) | |
Dew-point temperature (C) | |
Radiometric surface temperature (C) | |
Relative humidity (%) | |
Atmospheric vapor pressure at the level of measurement (hPa) | |
Atmospheric vapor pressure deficit at the level of measurement (hPa) | |
Wind speed (m s) | |
Friction velocity (m s) | |
Dew-point temperature at the source/sink height (C) | |
Aerodynamic temperature or source/sink height temperature (C) | |
“Effective” vapor pressure of evaporating front near the surface (hPa) | |
Saturation vapor pressure of the surface (hPa) | |
Saturation vapor pressure at the source/sink height (hPa) | |
Atmospheric vapor pressure at the source/sink height (hPa) | |
Atmospheric vapor pressure deficit at the source/sink height (hPa) | |
Equilibrium latent heat flux (W m) | |
Imposed latent heat flux (W m) | |
Evaporation as flux (W m) | |
Transpiration as flux (W m) | |
Evapotranspiration (evaporation transpiration) as depth of water (mm) | |
Potential evaporation as flux (W m) | |
Potential transpiration as flux (W m) | |
Wet environment evaporation as flux (W m) | |
Potential evaporation as flux according to Penman (W m) | |
Potential evaporation as flux according to Penman–Monteith (W m) | |
Potential evaporation as flux according to Priestley–Taylor (W m) | |
Potential evaporation as depth of water (mm) | |
Potential evaporation as depth of water according to Penman (mm) | |
Potential evaporation as depth of water according to Penman–Monteith (mm) | |
Potential evaporation as depth of water according to Priestley–Taylor (mm) | |
Wet environment evaporation as depth of water (mm) | |
Aerodynamic conductance (m s) | |
Stomatal/surface conductance (m s) | |
Momentum conductance (m s) | |
Quasi-laminar boundary-layer conductance (m s) | |
Maximum stomatal/surface conductance (m s) ( ) | |
Surface moisture availability (0–1) | |
Slope of saturation vapor pressure vs. temperature curve (hPa K (estimated at ) | |
Slope of the saturation vapor pressure and temperature between ( | |
vs. ( ) (approximated at (hPa K) | |
Slope of the saturation vapor pressure and temperature between ( | |
vs. ( (hPa K), estimated according to Mallick et al. (2015) | |
Slope of the saturation vapor pressure and temperature between ( | |
vs. ( ) (approximated at ) (hPa K) | |
Ratio between ( ) and ( ) | |
Latent heat of vaporization of water (J kg K) | |
Reference height (m) |
Continued.
Variables | Description |
---|---|
and | |
symbols | |
Effective source–sink height of momentum (m) | |
Roughness length (m) | |
Displacement height (m) | |
Psychrometric constant (hPa K) | |
Density of air (kg m) | |
Specific heat of dry air (MJ kg K) | |
Evaporative fraction (unitless) | |
Bowen ratio (unitless) | |
Priestley–Taylor parameter (unitless) | |
Decoupling coefficient (unitless) | |
Sc | Schmidt number (unitless) |
Pr | Prandtl number (unitless) |
von Kármán constant (0.4) |
The foundation of the development of STIC is based on the goal of finding an analytical solution of the two unobserved “state variables” ( and ) in the PM equation while exploiting the radiative ( and , meteorological (, ), and radiometric surface temperature () as external inputs. The fundamental assumption in STIC is the first-order dependence of and on the aerodynamic temperature () and soil moisture (through ). This assumption allows a direct integration of into the PM equation while simultaneously constraining the conductances through . Although the signal is implicit in , which appears in the numerator of the PM equation (Eq. 1), it may be noted that has a relatively weak dependence on (compared to the sensitivity of to soil moisture and ). Given that is a direct signature of the soil moisture availability, inclusion of in the PM equation also works to add water-stress controls in . Until now the explicit use of in the PM model was hindered due to the unavailability of any direct method to integrate into this model, and, furthermore, due to the lack of physical models expressing biophysical states of vegetation as a function of . Therefore, the majority of the PM-based modeling approaches strongly rely on surface reflectance and meteorology while exploiting the empirical leaf-scale parameterizations of the biophysical conductances (Prihodko et al., 2008; Bonan et al., 2014; Ershadi et al., 2015).
The PM equation is commonly expressed as where is the air density (kg m), is the specific heat of air (J kg K), is the psychrometric constant (hPa K), is the slope of the saturation vapor pressure vs. air temperature (hPa K), is the saturation deficit of the air (hPa) or vapor pressure deficit at the reference level, and is the net available energy (W m) (the difference between and ). The units of all the surface fluxes and conductances are in W m and m s, respectively. For a dense canopy, in the PM equation represents the canopy surface conductance. Although it is not equal to the canopy stomatal conductance, it contains integrated information on the stomata. For a heterogeneous landscape, in the PM equation is an aggregated surface conductance containing information on both canopy and soil. Traditionally, the two unknown “state variables” in Eq. (1) are and , and the STIC methodology is based on formulating “state equations” for these conductances that satisfy the PM model (Mallick et al., 2014, 2015). The PM equation is “closed” upon the availability of canopy-scale measurements of the two unobserved biophysical conductances, and if we assume the empirical models of and to be reliable. However, neither nor can be measured at the canopy scale or at larger spatial scales. Furthermore, as shown by some recent studies (Matheny et al., 2014; Van Dijk et al., 2015), a more appropriate and model is currently not available. This implies that a true “closure” of the PM equation is only possible through an analytical estimation of the conductances.
State equations
By integrating with standard surface energy balance (SEB) theory and vegetation biophysical principles, STIC formulates multiple “state equations” that eliminate the need for exogenous parametric submodels for and , associated aerodynamic variables, and land–atmosphere coupling. The state equations of STIC are as follows and their detailed derivations are described in Appendix A1. Here, is the temperature (C) at the source/sink height (or at the roughness length () or in-canopy air stream), is the atmospheric vapor pressure (hPa) at the source/sink height, is the saturation vapor pressure (hPa) at the source/sink height, is the evaporative fraction (the ratio of and ), is the Priestley–Taylor parameter (unitless) (Priestley and Taylor, 1972), and is a unitless quantity that describes the relative wetness (or moisture availability) of the surface. controls the transition from potential to actual evaporation and hence is critical for providing a constraint against which the conductances can be estimated ( estimation is explained in Appendix A2). Given values of , , , and or , the four state equations (Eqs. 2–5) can be solved simultaneously to derive analytical solutions for the four state variables. This also produces a “closure” of the PM model, which is independent of empirical parameterizations for both and . However, the analytical solutions to the above state equations have four accompanying unknowns, (surface moisture availability), (vapor pressure at the source/sink height), (saturation vapor pressure at the source/sink height), and the Priestley–Taylor coefficient (), and as a result there are four equations with eight unknowns. Consequently an iterative solution is needed to determine the four unknown variables (as described in Appendix A2), which is a further modification of the STIC1.1 framework (Mallick et al., 2015). The present version of STIC is designated as STIC1.2 and its uniqueness is the physical integration of into a combined structure of the PM and Shuttleworth–Wallace (SW, hereafter) (Shuttleworth and Wallace, 1985) models to estimate the source/sink height vapor pressures (Appendix A2). In addition to physically integrating observations into a combined PM–SW framework, STIC1.2 also establishes a feedback loop describing the relationship between and , coupled with canopy–atmosphere components relating to and . To estimate , the radiometric surface temperature () is extensively used in a physical retrieval framework, thus treating as an external input. In Eq. (5), the Priestley–Taylor coefficient () appeared due to the use of the advection–aridity (AA) hypothesis (Brutsaert and Stricker, 1979) for deriving the state equation of (Supplement S1). However, instead of optimizing as a “fixed parameter”, we have developed a physical equation of (Eq. A15) and numerically estimated as a “variable”. The derivation of the equation for is described in Appendix A2. The fundamental differences between STIC1.2 and earlier versions are described in Table A1.
In STIC1.2, is a function of , and they are not assumed equal ( ). The analytical expression of is dependent on and the estimation of is based on . To further elaborate this point on the inequality of and , we show an intercomparison of retrieved vs. for forest and pasture (Fig. A2). This indicates the distinct difference of the retrieved from for the two different biomes.
Partitioning
The terrestrial latent heat flux is an aggregate of both transpiration () and evaporation () (sum of soil evaporation and interception evaporation from the canopy). During rain events the land surface becomes wet and tends to approach the potential evaporation (), while surface drying after rainfall causes to approach the potential transpiration rate () in the presence of vegetation, or zero without any vegetation. Hence, at any time is a mixture of these two end-member conditions depending on the degree of surface moisture availability or wetness () (Bosveld and Bouten, 2003; Loescher et al., 2005). Considering the general case of evaporation from an unsaturated surface at a rate less than the potential, is the ratio of the actual to potential evaporation rate and is considered as an index of evaporation efficiency during a given time interval (Boulet et al., 2015). Partitioning of into and was performed according to Mallick et al. (2014) as follows: The estimates of in the current method consist of an aggregated contribution from both interception and soil evaporation, and no further attempt is made to separate these two components. In the Amazon forest, soil evaporation has a negligible contribution, while the interception evaporation contributes substantially to the total evaporative fluxes, and therefore the partitioning of into and is crucial. After estimating , was estimated according to the Penman equation (Penman, 1948) and was estimated as the residual in Eq. (6).
In this study, we use the term “canopy conductance” instead of “stomatal conductance” given that the term “stomata” is applicable at the leaf scale only. As stated earlier, for a heterogeneous surface, should principally be a mixture of the canopy surface (integrated stomatal information) and soil conductances. However, given the high vegetation density of the Amazon Basin, the soil surface exposure is negligible, and hence we assume to be the canopy-scale aggregate of the stomatal conductance. Similarly, a different exists for soil–canopy, sun–shade, and dry–wet conditions (Leuning, 1995), which are currently integrated into a lumped (given the big-leaf nature of STIC). From the big-leaf perspective, it is generally assumed that the aerodynamic conductance of water vapor and heat are equal (Raupach, 1998). However, to obtain partitioned aerodynamic conductances, explicit partitioning of is needed, which is beyond the scope of the current paper.
Evaluating and
Due to the lack of direct canopy-scale measurements, a rigorous evaluation of cannot be performed. To evaluate the STIC retrievals of (), we adopted three different methods.
- a.
By using the measured friction velocity () and wind speed () at the EC towers and using the equation of Baldocchi and Ma (2013) () in which was expressed as the sum of turbulent conductance and canopy (quasi-laminar) boundary-layer conductance aswhere is the von Kármán constant, 0.4; Sc is the Schmidt number; Pr is the Prandtl number, and their ratio is generally considered to be unity. Here the conductances of momentum, sensible, and latent heat fluxes are assumed to be identical (Raupach, 1998).
- b.
By inverting observations for wet conditions, hence assuming and estimating () as
- c.
By inverting the aerodynamic equation of and estimating a hybrid () from observed and STIC as (),Like , direct verification of STIC () could not be performed, as canopy-scale observations are not possible with current measurement techniques. Although leaf-scale measurements are relatively straightforward, these values are not comparable to values retrieved at the canopy scale. However, assuming -based as the baseline aerodynamic conductance, we have estimated canopy-scale by inverting the PM equation () (Monteith, 1995) to evaluate by exploiting in conjunction with the available , , , and measurements from the EC towers.
Decoupling coefficient and biophysical controls
The decoupling coefficient or “Omega” () is a dimensionless coefficient ranging from 0.0 to 1.0 (Jarvis and McNaughton, 1986) and is considered as an index of the degree of stomatal control on transpiration relative to the environment. The equation of is as follows: Introducing into the Penman–Monteith (PM) equation for results in where is the equilibrium latent heat flux, which depends only on and would be obtained over an extensive surface of uniform moisture availability (Jarvis and McNaughton, 1986; Kumagai et al., 2004). is the imposed latent heat flux, which is “imposed” by the atmosphere on the vegetation surface through the effects of vapor pressure deficit (triggered under limited soil moisture availability), and becomes proportional to .
When the ratio is very small (i.e., water-stress conditions), stomata principally control the water loss, and a change in will result in a nearly proportional change in transpiration. Such conditions trigger a strong biophysical control on transpiration. In this case the value approaches zero and vegetation is believed to be fully coupled to the atmosphere. In contrast, for a high / ratio (i.e., high water availability), changes in will have little effect on the transpiration rate, and transpiration is predominantly controlled by . In this case the value approaches unity, and vegetation is considered to be poorly coupled to the atmosphere.
Given that both and are the independent estimates in STIC1.2, the concept of was used to understand the degree of biophysical control on , which indicates the extent to which the transpiration fluxes are approaching the equilibrium limit. However, the biophysical characterization of and through STIC1.2 significantly differs from previous approaches (Ma et al., 2015; Chen et al., 2011; Kumagai et al., 2004), and the fundamental differences are centered on the specifications of and (as described in Table A2). While the estimation of in previous approaches was based on and , the estimation of was based on inversion of observed based on the PM equation (e.g., Stella et al., 2013). However, none of these approaches allow independent quantification of biophysical controls of , as is constrained by itself.
Datasets
Eddy covariance and meteorological quantities
We used the LBA (Large-Scale Biosphere-Atmosphere Experiment in Amazonia) data for quantifying the biophysical controls on the evaporative flux components. LBA was an international research initiative conducted during 1995–2005 to study how Amazonia functions as a regional entity within the larger Earth system, and how changes in land use and climate will affect the hydrological and biogeochemical functioning of the Amazon ecosystem (Andreae et al., 2002).
Overview of the LBA tower sites. Here, (–) refers to (S) and (W) for latitude and longitude, respectively.
Biome | PFT | Site | LBA | Data | Latitude | Longitude | Tower | Annual |
---|---|---|---|---|---|---|---|---|
code | availability | () | () | height | rainfall | |||
period | (m) | (mm) | ||||||
Forest | Tropical | Manaus | K34 | Jun 1999 to | 2.609 | 60.209 | 50 | 2329 |
rainforest | KM34 | Sep 2006 | ||||||
(TRF) | ||||||||
Forest | Tropical | Santarem | K67 | Jan 2002 to | 2.857 | 54.959 | 63 | 1597 |
moist | KM67 | Jan 2006 | ||||||
forest | ||||||||
(TMF) | ||||||||
Forest | Tropical | Santarem | K83 | Jul 2000 to | 3.018 | 54.971 | 64 | 1656 |
moist | KM83 | Dec 2004 | ||||||
forest | ||||||||
(TMF) | ||||||||
Forest | Tropical | Reserva | RJA | Mar 1999 to | 10.083 | 61.931 | 60 | 2354 |
dry forest | Biológica | Oct 2002 | ||||||
(TDF) | Jarú | |||||||
Pasture | Pasture | Santarem | K77 | Jan 2000 to | 3.012 | 54.536 | 18 | 1597 |
(PAS) | KM77 | Dec 2001 | ||||||
Pasture | Pasture | Fazenda | FNS | Mar 1999 to | 10.762 | 62.357 | 8.5 | 1743 |
(PAS) | Nossa | Oct 2002 | ||||||
Senhora |
Examples of monthly averages of the diurnal time series of canopy-scale (a) and (b) estimated for two different biomes (forest and pasture) in the Amazon Basin (LBA sites K34 and FNS). The time series of four different estimates and their corresponding estimates are shown here.
[Figure omitted. See PDF]
A network of eddy covariance (EC) towers was operational during the LBA
experiment, such that data from nine EC towers were obtained from the ORNL
Distributed Archive Active Centre
(
Results
Evaluating , , and surface energy balance fluxes
Examples of monthly averages of the diurnal cycles of the four different estimates and their corresponding estimates over two different PFTs (K34 for forest and FNS for pasture) reveal that and tend to be generally higher for the forest than their counterparts, varying from 0 to 0.06 m s and 0 to 0.04 m s respectively (Fig. 1a and b). The magnitude of varied between 0 and 0.025 m s for the pasture (Fig. 1a), while values were less than half of those estimated over the forest (0–0.01 m s) (Fig. 1b). The conductances showed a marked diurnal variation expressing their overall dependence on net radiation, vapor pressure deficit, and surface temperature. Despite the absolute differences between the conductances from the different retrieval methods, their diurnal patterns were comparable.
The canopy-scale evaluation of is illustrated in Fig. 2a (and Table 3), combining data from the four PFTs. Estimated values range between zero and 0.1 m s and show modest correlation ( 0.44) ( range between 0.22 [0.18] and 0.55 [0.12]) between and , with regression parameters ranging from 0.81 (0.023) to 1.07 (0.047) for the slope and 0.0019 (0.0006) to 0.0006 (0.0006) m s for the offset (Table 3). The root mean squared deviation (RMSD) varied between 0.007 (TDF) and 0.013 m s (TRF). Statistical comparisons between and revealed relatively low RMSD and high correlation between them (RMSD 0.007 m s and 0.77) as compared to the error statistics between and (RMSD 0.011 m s and 0.50) (Fig. 2b and c). The residuals between and are plotted as a function of and in Fig. 2d with the aim of ascertaining whether significant biases are introduced by ignoring wind and shear information within STIC1.2. As illustrated in Fig. 2d, there appears to be a weak systematic relationship between the residual difference with either or ( 0.26 and 0.17). However, a considerable relationship was found between wind- and shear-driven (i.e., ) vs. , , and ( 0.83, 0.48, and 0.42) (Fig. 2e and f), which indicates that these three energy and water constraints can explain 69, 23, and 17 % variance of , respectively.
Canopy-scale evaluation of hourly is presented in Fig. 3a (and Table 3), combining data from the four PFTs. Estimated values range between zero and 0.06 m s for and show reasonable correlation ( 0.39) ( range between 0.14 [0.04] and 0.58 [0.12]) between and , with regression parameters ranging between 0.30 (0.022) and 0.85 (0.025) for the slope and between 0.0024 (0.0003) and 0.0097 (0.0007) m s for the offset (Table 3). The RMSD varied between 0.007 (PAS) and 0.012 m s (TRF and TDF). Given that significantly controls , we also examined whether biases in are introduced by ignoring wind and shear information within STIC. The scatter plots between the residual difference ( ) vs. both and (Fig. 3b) showed residuals to be evenly distributed across the entire range of and , and no systematic pattern was evident.
(a) Comparison between STIC-derived () with an estimated aerodynamic conductance based on friction velocity () and wind speed () according to Baldocchi and Ma (2013) (), (b) comparison between with an inverted () based on EC observations of and , (c) comparison between with a hybrid ( based on EC observations of and estimated over the LBA EC sites, (d) comparison between residual differences vs. and , and (e, f) relationship between wind- and shear-derived vs. , , and over the LBA EC sites.
[Figure omitted. See PDF]
(a) Comparison between STIC-derived () and computed by inverting the PM model () over the LBA EC sites, where was used as aerodynamic input in conjunction with tower measurements of , radiation and meteorological variables, and (b) residual differences vs. wind speed () and friction velocity () over the LBA EC sites.
[Figure omitted. See PDF]
Comparison between STIC-derived (a) and (b) over four different PFTs in the Amazon Basin (LBA tower sites). MAPD is the percent error defined as the mean absolute deviation between the predicted and observed variables divided by the mean observed variable.
[Figure omitted. See PDF]
The reliability of STIC1.2-based and retrievals was further verified by evaluating and estimates (Fig. 4). Both the predicted and are generally in good agreement with the observations, with substantial correlation () ( from 0.61 to 0.94), reasonable RMSD of 33 and 37 W m, and mean absolute percent deviation (MAPD) of 14 and 32 % between the observed and STIC fluxes (Fig. 4), respectively. Regression parameters varied between 0.96 (0.008) and 1.14 (0.010) for the slope and between 16 (2) and 2 (2) W m for the offset for (Table 4), whereas for , these were 0.60 (0.025) to 0.89 (0.035) for the slope and 9 (1) to 29 (2) W m for the offset (Table 3), respectively. The RMSD in varied from 20 to 31 W m and from 23 to 34 W m for (Table 3).
The evaluation of the conductances and surface energy fluxes indicates some efficacy for the STIC-derived fluxes and conductance estimates that represent a weighted average of these variables over the source area around the EC tower.
Canopy coupling, transpiration, and evaporation
From Fig. 5a an overall weak to moderate relationship ( 0.31 to 0.42) is apparent between the coupling (i.e., 1 ) and , where is negatively related to the coupling for all the PFTs, thus indicating the influence of weak to moderate biophysical controls on throughout the year in addition to radiative controls. The biophysical control was substantially enhanced in TRF ( increased from 0.36 to 0.53 and 0.60) (47 to 67 % increase) and TMF ( increased from 0.31 to 0.53 and 0.58) (70 to 85 % increase) during the dry seasons (July–September) (Fig. 5a). A profound increase in the biophysical control on during the dry season was also found in TDF (52 % increase) and PAS (37 % increase) (Fig. 5a). The negative relationship ( 0.29 to 0.45) between (1 ) and (Fig. 5b) in all four PFTs indicated the role of aerodynamic control on . The aerodynamic control was also enhanced during the dry seasons, as shown by the increased negative correlation ( 0.50 to 0.69) (Fig. 5b) between (1 ) and .
Correlation of coupling (1 ) with (a) transpiration ( and (b) evaporation () and over four different PFTs by combining data for all the years, only during dry seasons for all the years, and during drought year 2005. Data for 2005 were not available for TDF and PAS. (c)–(e) Examples of the diurnal pattern of (black lines), (grey dotted lines) and (grey solid lines) estimated over two ecohydrologically contrasting biomes (K34 for forest and FNS for pasture) in the Amazon Basin (LBA tower sites) during wet and dry seasons.
[Figure omitted. See PDF]
Illustrative examples of the diurnal variations of , , and for two different PFTs with different annual rainfall (2329 mm in rainforest, K34, and 1597 mm in pasture, FNS) for 3 consecutive days during both dry and wet seasons are shown in Fig. 5c–f. This shows the morning rise of and a near-constant afternoon in the wet season (Fig. 5c and 5d), thus indicating no biophysical controls on and during this season. By contrast, during the dry season, the morning rise in is followed by a decrease during noontime (15 to 25 % increase in coupling in forest and pasture) (Fig. 5e and f) due to dominant biophysical control, which is further accompanied by a transient increase from mid-afternoon till late afternoon, and steadily declines thereafter. Interestingly, coupling was relatively higher in pasture during the dry seasons, and the reasons are detailed in the following section and discussion.
and vs. transpiration and evaporation
Scatter plots between and vs. and showed a triangular pattern that became wider with increasing conductances (Fig. 6). To explain this behavior of vs. and , we further examined the entire mechanism of conductance– interactions through two-dimensional scatters between and conductances for two consecutive diurnal cycles during wet and dry seasons over rainforest and pasture sites with different annual rainfall (e.g., K34 as a wet site and FNS as a dry site, annual rainfall 2329 and 1597 mm) (Fig. 7). Our results confirm the occurrence of diurnal hysteresis between and and explain the reason for the shape of the curves obtained in Fig. 6. During the wet season, a distinct environmental control is detectable on and in the morning hours (Fig. 7a and b) in both the PFTs where and increased as a result of increasing , , and . From the late morning to afternoon, a near-constant (forest) or negligible (pasture) increase in is observed despite substantial reduction of both and (25 to 50 % decrease), after which starts decreasing. This behavior of was triggered due to the concurrent changes in (15 to 50 % change), (20 to 60 % change), and (5 to 14 % change), which indicates the absence of any dominant biophysical regulation on during the wet season (Fig. 7a and b). On the contrary, in the dry season, although the morning rise in is steadily controlled by the integrated influence of environmental variables, a modest to strong biophysical control is found for both PFTs during the afternoon, where substantially decreased with decreasing conductances (Fig. 7c and d). This decrease in is mainly caused by the reduction in as a result of increasing and (as seen later in Fig. 8a and c). In the dry season, the area under the hysteretic relationship between , , and environmental variables was substantially wider in pasture (Fig. 7d) than for the rainforest (Fig. 7c), which is attributed to a greater hysteresis area between and in pasture as a result of reduced water supply. The stronger hysteresis effects in pasture during the dry season (Fig. 7d) ultimately led to the stronger relationship between coupling and (as seen in Fig. 5a).
Scatter plots of transpiration () and evaporation () vs. and over four different PFTs in the Amazon Basin (LBA tower sites).
[Figure omitted. See PDF]
Illustrative examples of the occurrence of diurnal hysteresis of transpiration () during wet and dry seasons with canopy and environmental controls over two different sites with different annual rainfall (2329 and 1597 mm, respectively) in the Amazon Basin (LBA tower sites K34 and FNS).
[Figure omitted. See PDF]
(a) Response of retrieved to atmospheric vapor pressure deficit () for different classes of net radiation (), (b) response of retrieved to transpiration for different classes of , (c) response of retrieved to radiometric surface temperature () for different classes , and (d) relationship between retrieved and radiometric surface temperature and air temperature difference ( ) in the Amazon Basin (LBA tower sites).
[Figure omitted. See PDF]
Factors affecting variability of and
The sensitivity of stomatal conductance to vapor pressure deficit is a key governing factor of transpiration (Ocheltree et al., 2014; Monteith, 1995). We examined whether the feedback or feed-forward response hypothesis (Monteith, 1995; Farquhar, 1987) between , , and is reflected in our canopy-scale retrievals. Combining data of all the PFTs, we found an exponential decline of in response to increasing regardless of the variations of net radiation (Fig. 8a). High is consistent with high humidity and low evaporative demand. Five negatively logarithmic scatters fit the data with values of 0.38 (0 150 W m), 0.63 (150 300 W m), 0.73 (300 450 W m), 0.78 (450 600 W m), and 0.87 ( 600 W m). The sensitivity of to was at the maximum in the high range beyond 600 W m and the sensitivity progressively declined with declining magnitude of (0–150 W m).
Scatter plots between and for different levels of revealed a linear pattern between them for a wide range of (20 0 hPa) (Fig. 8b). Following Monteith (1995), isopleths of are delineated by the solid lines passing through on the axis and through on the axis. Isobars of (dotted lines) pass through the origin because approaches zero as approaches zero. Figure 8b shows substantial reduction of with increasing without any increase in , like an inverse hyperbolic pattern to (Monteith, 1995; Jones, 1998). For all the PFTs, an active biological (i.e., stomatal) regulation maintained almost constant when was changed from low to high values (Fig. 8b). At high (above 10 hPa), after an initial increase in with , approached a maximum limit and remained nearly independent of (Fig. 8b). Among all the levels, the maximum control of on variability (62 to 80 %) was found at high atmospheric water demand (i.e., 30 hPa 20 hPa). The scatter plots between and (Fig. 8c) for different levels of revealed an exponential decline in with increasing and atmospheric water demand. When retrieved was plotted against the radiometric surface temperature and air temperature difference ( ), an exponential decline in was found in response to increasing ( ) (Fig. 8d). High is persistent with low ( ) irrespective of the variations in (with the exception of very low ). Four negatively logarithmic scatters fit the vs. ( ) relationship with values of 0.28 (150 300 W m), 0.55 (3000 450 W m), 0.64 (450 600 W m), and 0.77 ( 600 W m).
(a) Scatter plots between source/sink height (or in-canopy) vapor pressure deficit () and atmospheric vapor pressure deficit () for two different classes of / ratios over four PFTs, which clearly depicts a strong coupling between and for low / ratios. (b) Histogram distribution of / ratios over the four PFTs in the Amazon Basin (LBA tower sites). (c) Scatter plots between / ratio vs. surface air temperature difference ( ) for the four PFTs during the wet season and dry season in the Amazon Basin (LBA tower sites).
[Figure omitted. See PDF]
Discussion
Evaluating , , and surface energy balance fluxes
The aerodynamic conductance retrieved with STIC1.2 showed acceptable correlation and valid estimates of when compared against an empirical model that uses and to derive (Figs. 1 and 2a) and two other inversion-/hybrid-based estimates. The differences between and were mainly attributed to the structural differences and empirical nature of the parameterization for the near-surface boundary-layer conductance ((2/)(Sc/Pr)) in , which results in some discrepancies between and , particularly in the pasture (Fig. 2a). The extent to which the structural discrepancies between and relate to actual differences in the conductances for momentum vs. heat is beyond the scope of this paper, and a detailed investigation using data on atmospheric profiles of wind speed, temperature, etc. are needed to actually quantify such differences. Momentum transfer is associated with pressure forces and is not identical to heat and mass transfer (Massman, 1999). In STIC1.2, is directly estimated and is a robust representative of the conductances to heat (and water vapor) transfer, whereas estimates based on and are more representative of the momentum transfer. Therefore, the difference between the two different estimates (Fig. 2) can be largely attributed to the actual difference in the conductances for momentum and heat (water vapor). The turbulent conductance equation (/) in is also very sensitive to the uncertainties in the sonic anemometer measurement (Contini et al., 2006; Richiardone et al., 2012). However, the evidence of a weak systematic relationship between the residuals and (Fig. 2d) and the capability of the thermal (, radiative (), and meteorological (, variables in capturing the variability of (Fig. 2e and f) indicates the diagnostic potential estimates to explain the wind-driven variability. Excluding might introduce errors in cases where wind is the only source of variations in and surface fluxes (Mallick et al., 2015). In general, the accuracies in commonly used parametric estimates based on and surface roughness parameters several meters distant from canopy foliage are limited due to the uncertainties concerning the attenuation of close to the vegetation surface (Meinzer et al., 1997; Prihodko et al., 2008). The magnitude of near the foliage can be substantially lower than that measured considerably away at some reference location above or within the canopy (Meinzer et al., 1997). Notwithstanding the inequalities of estimated with different methods, it is challenging to infer the accuracy of the different estimates. It is imperative to mention that is one of the main anchors in the PM–SW model because it not only appears in the numerator and denominator of these models, but also provides feedback to , , and (the seminal paper of Jarvis and McNaughton, 1986). Therefore, the estimates of in the PM–SW framework are very sensitive to parameterization of and stable estimates might be possible if estimation is unambiguous (Holwerda et al., 2012; Van Dijk et al., 2015). Given the lack of consensus in the community on the “true” and from the nature of surface flux validation results (Fig. 4), it appears that tends to be the appropriate aerodynamic conductance that satisfies the PM–SW equation. Discrepancies between and originated from the differences in estimates between the two methods.
Comparative statistics for the STIC- and tower-derived hourly and for a range of PFTs in the Amazon Basin (LBA tower sites). Values in parentheses are 1 standard deviation (standard error for correlation).
vs. | vs. | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
PFTs | RMSD | Slope | Offset | RMSD | Slope | Offset | ||||
(m s) | (m s) | (m s) | (m s) | |||||||
TRF | 0.013 | 0.41 | 1.07 | 0.0031 | 0.012 | 0.14 | 0.39 | 0.0097 | 1159 | |
(0.03) | (0.047) | (0.0008) | (0.04) | (0.039) | (0.0007) | |||||
TMF | 0.012 | 0.55 | 0.81 | 0.0006 | 0.009 | 0.55 | 0.85 | 0.0032 | 1927 | |
(0.12) | (0.023) | (0.0006) | (0.12) | (0.025) | (0.0005) | |||||
TDF | 0.007 | 0.49 | 0.89 | 0.0019 | 0.012 | 0.33 | 0.30 | 0.0050 | 787 | |
(0.15) | (0.041) | (0.0006) | (0.19) | (0.022) | (0.0005) | |||||
PAS | 0.012 | 0.22 | 1.03 | 0.0059 | 0.007 | 0.58 | 0.65 | 0.0024 | 288 | |
(0.18) | (0.083) | (0.0007) | (0.12) | (0.025) | (0.0003) | |||||
Mean | 0.012 | 0.44 | 0.76 | 0.0047 | 0.010 | 0.39 | 0.63 | 0.0046 | 4161 | |
(0.10) | (0.016) | (0.003) | (0.08) | (0.016) | (0.0003) |
: number of data points; RMSD: root mean square deviation between predicted () and observed () variables .
Comparative statistics for the STIC- and tower-derived hourly and for a range of PFTs in the Amazon Basin (LBA tower sites). Values in parentheses are 1 standard deviation (standard error for correlation).
PFTs | RMSD | Slope | Offset | RMSD | Slope | Offset | ||||
---|---|---|---|---|---|---|---|---|---|---|
(W m) | (W m) | (W m) | (W m) | |||||||
TRF | 28 | 0.96 | 1.10 | 16 | 34 | 0.52 | 0.60 | 29 | 1159 | |
(0.007) | (0.008) | (2) | (0.030) | (0.025) | (2) | |||||
TMF | 20 | 0.98 | 1.08 | 11 | 23 | 0.71 | 0.61 | 20 | 1927 | |
(0.004) | (0.004) | (1) | (0.019) | (0.014) | (1) | |||||
TDF | 26 | 0.96 | 0.96 | 7 | 30 | 0.66 | 0.89 | 20 | 787 | |
(0.009) | (0.008) | (2) | (0.032) | (0.035) | (3) | |||||
PAS | 31 | 0.96 | 1.14 | 2 | 33 | 0.88 | 0.67 | 9 | 288 | |
(0.009) | (0.010) | (2) | (0.016) | (0.011) | (1) | |||||
Mean | 33 | 0.94 | 1.04 | 1 | 37 | 0.61 | 0.58 | 24 | 4161 | |
(0.005) | (0.005) | (1) | (0.021) | (0.009) | (2) |
Despite the good agreement between the measured and predicted and (Fig. 4, Table 4), the larger error in was associated with the higher sensitivity of to the errors in (due to poor emissivity correction) (Mallick et al., 2015). Since the difference between and is considered to be the primary driving force of (van der Tol et al., 2009), the modeled errors in are expected to arise due to the uncertainties associated with .
Canopy coupling, , and vs. transpiration and evaporation
The correlation analysis between 1 and revealed the extent of biophysical and radiative controls on (Fig. 5). The degree of biophysical control is a function of the ratio of to . Minor biophysical control on was apparent for forest and pasture during the wet seasons (Fig. 5c and d) as a result of a high / ratio along with increasing . Such conditions stimulate local humidification of air surrounding the canopy and uncoupling of the in-canopy vapor pressure deficit () from that in the air above (i.e., ) (Meinzer et al., 1997; Motzer et al., 2005) (Fig. 9a), which implies that becomes largely independent of . By contrast, an enhanced biophysical control on was apparent during the dry season and drought year 2005 during the period of reduced water supply, particularly over PAS (Figs. 5e, f, and 7). Such a condition leads to a relatively dry canopy surface and substantially high compared to , thus resulting in low / ratios regardless of their absolute values (Meinzer et al., 1993; McNaughton and Jarvis, 1991). Here, a fractional change in results in an equivalent fractional change in . This impedes transpiration from promoting local equilibrium of and minimizing (or maximizing) the gradient between and atmospheric vapor pressure deficit () (i.e., or ) (Eq. A10) (Fig. 9a), thereby resulting in strong coupling between and (Meinzer et al., 1993; Jarvis and McNaughton, 1986). Besides, a supplemental biophysical control on might have been imposed as a consequence of a direct negative feedback of and on (McNaughton and Jarvis, 1991; Jarvis, 1986). An increase in (or ) beyond a certain limit decreases (Figs. 7 and 8), resulting in a low and narrow increase in , despite a steady increase in and . The combination of a negative feedback response between and with the overall radiative–aerodynamic coupling significantly dampens the variation of transpiration in PAS and TDF in the dry season, thus featuring increased biophysical control in these PFTs. These results are in agreement with von Randow et al. (2012), who found enhanced biophysical control on for the pasture during the dry season. For the wet season, evidence of minor biophysical control indicates the dominance of -driven equilibrium evaporation in these PFTs (Hasler and Avissar, 2007; da Rocha et al., 2009; Costa et al., 2010). In the TRF and TMF, 94 and 99 % of the retrieved / ratios fall above 0.5, and only 1 and 6 % of the retrieved / ratios fall below the 0.5 range (Fig. 9b). In contrast, 90 and 73 % of the / ratios range above 0.5, and 10 to 27 % of the / ratios were below 0.5 for TDF and PAS, respectively (Fig. 9b). This shows that, although radiation control is prevailing in all the sites, biophysical control is relatively stronger in TDF and PAS as compared to the other sites. For large / ratios, the conditions within the planetary boundary layer (PBL) become decoupled from the synoptic scale (McNaughton and Jarvis, 1991) and the net radiative energy becomes the important regulator of transpiration. For small / ratios (e.g., in the dry season), the conditions within the PBL are strongly coupled to the atmosphere above by rapid entrainment of air from the capping inversion and by some ancillary effects of sensible heat flux on the entrainment (McNaughton and Jarvis, 1991). These findings substantiate the earlier theory of McNaughton and Jarvis (1991), who postulated that large / ratios result in minor biophysical control on canopy transpiration due to the negative feedback on the canopy from the PBL. The negative relationship between 1 and (Fig. 5b) over all the PFTs is due to the feedback of on . However, over all the PFTs, a combined control of and environmental variables on again highlighted the impact of realistically estimated on (Holwerda et al., 2012).
It is important to mention that forests are generally expected to be better coupled to the atmosphere as compared to the pastures, which is related to generally higher of the forests (due to high surface roughness). This implies that forests exhibit stronger biophysical control on . However, due to the broad leaves of the rainforests (larger leaf area index) and higher surface wetness (due to higher rainfall amounts), the wet surface area is much larger in the forests than in the pastures. This results in much higher values for forests than for pastures during the wet season ( ), and / 1. Consequently, no significant difference in coupling was found between them during the wet season (Fig. 5c and d). Despite the absolute differences in and between forest and pasture, the high surface wetness is largely offsetting the expected difference between them. Although the surface wetness is substantially lower during the dry season, the high water availability in the forests due to the deeper root systems helps in maintaining a relatively high compared to the pastures. Hence, despite (forest) (pasture) during the dry season, substantially lower values for the pasture result in a lower / ratio for the pasture compared to the forest, thus causing more biophysical control on during the dry season. The relatively better relationship between coupling vs. in PAS and TDF during the dry season was also attributed to high surface air temperature differences ( ) in these PFTs that resulted in low / ratios (Fig. 9c).
Factors affecting and variability
The stomatal feedback-response hypothesis (Monteith, 1995) also became apparent at the canopy scale (Fig. 8a and b), which states that a decrease in with increasing is caused by a direct increase in (Monteith, 1995; Matzner and Comstock, 2001; Streck, 2003), and responds to the changes in the air humidity by sensing rather than . This feedback mechanism is found because of the influence of on both and , which in turn changes by influencing the air humidity (Monteith, 1995). The change in is dominated by an increase in the net available energy, which is partially offset by an increase in . After the net energy input in the canopy exceeds a certain threshold, starts decreasing even if increases. High increases the water potential gradient between guard cells and other epidermal cells or reduces the bulk leaf water potential, thus causing stomatal closure (Monteith, 1995; Jones, 1998; Streck, 2003). The control of soil water on transpiration also became evident from the scatter plots between vs. and for different levels (Figs. 8b, c, and 7). Denmead and Shaw (1962) hypothesized that reduced and stomatal closure occurs at moderate to higher levels of soil moisture (high ) when the atmospheric demand of water vapor increases (high . The water content in the immediate vicinity of the plant root depletes rapidly at high , which decreases the hydraulic conductivity of soil, and the soil is unable to efficiently supply water under these conditions. For a given evaporative demand and available energy, transpiration is determined by the / ratio, which is further modulated by the soil water availability. These combined effects tend to strengthen the biophysical control on transpiration (Leuzinger and Kirner, 2010; Migletta et al., 2011). The complex interaction between , , and (Fig. 8c) explains why different parametric models produce divergent results.
Although and estimates are interdependent on and (as shown in Figs. 6–8), the figures reflect the credibility of the conductances as well as transpiration estimates by realistically capturing the hysteretic behavior between biophysical conductances and water vapor fluxes, which is frequently observed in natural ecosystems (Zhang et al., 2014; Renner et al., 2016) (also Zuecco et al., 2016). These results are also compliant with the theories postulated earlier from observations that the magnitude of hysteresis depends on the radiation–vapor pressure deficit time lag, while the soil moisture availability is a key factor modulating the hysteretic transpiration–vapor pressure deficit relation as soil moisture declines (Zhang et al., 2014; O'Grady et al., 1999; Jarvis and McNaughton, 1986). This shows that despite being independent of any predefined hysteretic function, the interdependent conductance–transpiration hysteresis is still captured in STIC1.2.
Figure 8d is in accordance with the existing theory that under conditions of extremely high atmospheric turbulence (i.e., high ), a close coupling exists between the surface and the atmosphere, which causes and to converge (i.e., 0). When is low, the difference between and increases due to poor vertical mixing of the air.
Conclusions
By integrating the radiometric surface temperature () into a combined structure of the PM–SW model, we have estimated the canopy-scale biophysical conductances and quantified their control on the terrestrial evapotranspiration components in a simplified SEB modeling perspective that treats the vegetation canopy as “big-leaf”. The STIC1.2 biophysical modeling scheme is independent of any leaf-scale empirical parameterization for stomata and associated aerodynamic variables.
Stomata regulate the coupling between terrestrial carbon and water cycles, which implies that their behavior under global environmental change is critical to predicting vegetation functioning (Medlyn et al., 2011). The combination of variability in precipitation (Hilker et al., 2014) and land cover change (Davidson et al., 2012) in the Amazon Basin is expected to increase the canopy–atmosphere coupling of pasture or forest systems under drier conditions by altering the ratio of the biological and aerodynamic conductances. An increase in biophysical control will most likely be an indicator of shifting the transpiration from an energy-limited to a water-limited regime (due to the impact of , , and on the / ratio), with further consequences for the surface water balance and rainfall recycling. At the same time, a transition from forest to pasture or agriculture lands will substantially reduce the contribution of interception evaporation in the Amazon; hence, it will affect the regional water cycle. This might change the moisture regime of the Amazonian Basin and affect the moisture transport to other regions. In this context, STIC1.2 provides a new quantitative and internally consistent method for interpreting the biophysical conductances and is able to quantify their controls on the water cycle components in response to a range of climatic and ecohydrological conditions (excluding rising atmospheric CO) across a broad spectrum of PFTs. It could also provide the basis for improving existing land surface parameterizations for simulating vegetation water use at large spatial scales.
It should also be noted that although the case study described here provides general insights into the biophysical controls of and associated feedback between , , , and in the framework of the PM–SW equation, there is a tendency to overestimation of due to the embedded evaporation information in the current single-source composition of STIC1.2. For accurate characterization of canopy conductance, explicit partitioning of into transpiration and evaporation (both soil and interception) is one of the further scopes for improving STIC1.2, and this assumption needs to be tested further.
Data availability
The LBA eddy covariance datasets are available through
Description of STIC1.2
Derivation of “state equations” in STIC 1.2
Neglecting horizontal advection and energy storage, the surface energy balance equation is written as follows: Figure A1 shows that, while is controlled by a single aerodynamic resistance () (or 1/); is controlled by two resistances in series, the surface resistance () (or 1/) and the aerodynamic resistance to vapor transfer ( ). For simplicity, it is implicitly assumed that the aerodynamic resistance of water vapor and heat are equal (Raupach, 1998), and both the fluxes are transported from the same level from near surface to the atmosphere. The sensible and latent heat flux can be expressed in the form of aerodynamic transfer equations (Boegh et al., 2002; Boegh and Soegaard, 2004) as follows: where and are the air temperature and vapor pressure at the source/sink height (i.e., aerodynamic temperature and vapor pressure) or at the so-called roughness length (), where wind speed is zero. They represent the vapor pressure and temperature of the quasi-laminar boundary layer in the immediate vicinity of the surface level (Fig. A1), and can be obtained by extrapolating the logarithmic profile of down to . is the saturation vapor pressure at (hPa).
By combining Eqs. (A1)–(A3) and solving for , we get the following equation.
Combining the aerodynamic expressions of in Eq. (A3) and solving for , we can express in terms of , , , and .
While deriving the expressions for and , two more unknown variables are introduced ( and ), thus there are two equations and four unknowns. Therefore, two more equations are needed to close the system of equations.
An expression for is derived from the Bowen ratio () (Bowen, 1926) and evaporative fraction () (Shuttleworth et al., 1989) equation.
Schematic representation of the one-dimensional description of STIC1.2. In STIC1.2, a feedback is established between the surface-layer evaporative fluxes and source/sink height mixing and coupling, and the connection is shown in dotted arrows between , , , , and . Here, and are the aerodynamic and canopy (or surface in the case of partial vegetation cover) resistances, and are the aerodynamic and canopy conductances (reciprocal of resistances), is the saturation vapor pressure of the surface, is the saturation vapor pressure at the source/sink height, is the source/sink height temperature (i.e., aerodynamic temperature) that is responsible for transferring the sensible heat (), is the source/sink height vapor pressure, is the vapor pressure at the surface, is the roughness length, is the radiometric surface temperature, is the source/sink height dew-point temperature, is the surface moisture availability or evaporation coefficient, and are net radiation and ground heat flux, , , and are temperature, vapor pressure, and vapor pressure deficit at the reference height (), is the latent heat flux, and is the sensible heat flux, respectively.
[Figure omitted. See PDF]
Aerodynamic temperature obtained from STIC1.2 () vs. radiometric surface temperature () over two different biomes in the Amazon Basin. The regression equation of line of best fit is 0.67 (0.10) 10.59 (2.79) with 0.65.
[Figure omitted. See PDF]
This expression for introduces another new variable (); therefore, one more equation that describes the dependence of on the conductances ( and ) is needed to close the system of equations. In order to express in terms of and , we had adopted the advection–aridity (AA) hypothesis (Brutsaert and Stricker, 1979) with a modification introduced by Mallick et al. (2015). The AA hypothesis is based on a complementary connection between the potential evaporation (), sensible heat flux (), and ; and leads to an assumed link between and . However, the effects of surface moisture (or water stress) were not explicit in the AA equation, and Mallick et al. (2015) implemented a moisture constraint in the original AA hypothesis while deriving a “state equation” of (Eq. A8 below). A detailed derivation of the “state equation” for is described in Supplement S1 (also see Mallick et al., 2014, 2015). Estimation of , , , and is described in Appendix A2.
Iterative solution of , , , and in STIC 1.2
In STIC1.0 and 1.1 (Mallick et al., 2014, 2015), no distinction was made between the surface and source/sink height vapor pressures. Therefore, was approximated as the saturation vapor pressure at and was empirically estimated from based on the assumption that the vapor pressure at the source/sink height ranges between extreme wet–dry surface conditions. However, the level of and should be consistent with the level of the aerodynamic temperature () from which the sensible heat flux is transferred (Lhomme and Montes, 2014). The predictive use of the PM model could be hindered by neglecting of the feedbacks between the surface-layer evaporative fluxes and source/sink height mixing and coupling (McNaughton and Jarvis, 1984) and their impact on the canopy-scale conductances. Therefore, in STIC1.2, we have used physical expressions for estimating and followed by estimating and as described below. The fundamental differences between the STIC1.0, 1.1, and 1.2 modeling philosophies are described in Table A1.
An estimate of is obtained by inverting the aerodynamic transfer equation of .
Following Shuttleworth and Wallace (1985) (SW), the vapor pressure deficit () ( ) and vapor pressure () at the source/sink height are expressed as follows.
A physical equation of is derived by expressing the evaporative fraction () as a function of the aerodynamic equations of [ ] and [ ] as follows.
Combining Eqs. (A14) and (A8) (eliminating ), we can derive a physical equation of .
Following Venturini et al. (2008), can be expressed as the ratio of the vapor pressure difference to the vapor press deficit between the surface and atmosphere as follows. where is the dew-point temperature at source/sink height and is the air dew-point temperature; and are the psychrometric slopes of the saturation vapor pressure and temperature between the ( ) vs. ( ) and ( ) vs. ( ) relationship (Venturini et al., 2008), and is the ratio between ( ) and ( ). Despite driving the sensible heat flux, the comprehensive dry–wet signature of the underlying surface due to soil moisture variations is directly reflected in (Kustas and Anderson, 2009). Therefore, using in the denominator of Eq. (A16) tends to give a direct signature of the surface moisture availability (). In Eq. (A16), computation is challenging because both and are unknown. By decomposing the aerodynamic equation of , can be expressed as follows.
In the earlier STIC versions, was approximated at , was approximated at , was estimated from , , , and related saturation vapor pressures (Mallick et al., 2014, 2015), and was estimated from Eq. (A16) (estimation of and was stand-alone earlier). However, since depends on and , an iterative procedure is applied to estimate and as described in the following paragraph.
(a) Convergence of the iteration method for retrieving the source/sink height (or in-canopy) vapor pressures ( and ) and Priestley–Taylor coefficient (). (b) Convergence of the iteration method for retrieving the surface wetness () and source/sink height dew-point temperature (). The initial values of , , and were determined with 1.26. The process is then iterated by updating , , , , , and in subsequent iterations with the previous estimates of , , and .
[Figure omitted. See PDF]
In STIC1.2, an initial value of is assigned as 1.26 and initial estimates of and are obtained from and as 6.13753 and . Initial and were estimated as described in Mallick et al. (2014, 2015). With the initial estimates of these variables, first estimates of the conductances, , , and are obtained. The process is then iterated by updating (using Eq. A9), (using Eq. A10), (using Eq. A11), (using Eq. A17 with estimated at ), (using Eq. A16), and (using Eq. A15), with the first estimates of , , and , and recomputing , , , , and in the subsequent iterations with the previous estimates of , , , and until the convergence of is achieved. Stable values of , , , , , and are obtained within 25 iterations. Illustrative examples of the convergence of , , , , and are shown in Fig. A3.
To summarize, the computational steps of the conductances and evaporative fluxes in STIC are
-
Step 1: analytical solution of the conductances, and , by solving the “state equations” (Eqs. 2–5);
-
Step 2: initial estimates of the conductances ( and ), , , , and ;
-
Step 3: simultaneous iteration of , , , , , and , and final estimation of the conductances ( and ), , , , and ; and
-
Step 4: partitioning into and .
Differences in the modeling philosophies of source/sink height vapor pressures (, ) and dew-point temperature (), surface wetness (), and between STIC1.0, STIC1.1, and STIC1.2.
Variable | Principles | ||
---|---|---|---|
estimation | STIC1.0 | STIC1.1 | STIC1.2 |
(Mallick et al., 2014) | (Mallick et al., 2015) | (this study – Mallick et al., 2016) | |
Saturation | was approximated as the | Same as STIC1.0 | is estimated through numerical |
vapor | saturation vapor pressure at | iteration by inverting the aerodynamic | |
pressure at | . | equation of (as described in Appendix A2). | |
source/sink | |||
height () | |||
Actual vapor | was empirically estimated | Same as STIC1.0 | is estimated as , where |
pressure at | from based on the | was iteratively estimated by combining | |
source/sink | assumption that the vapor | PM with the Shuttleworth–Wallace | |
height () | pressure at the source/sink | approximation (as described in Appendix A2). | |
height ranges between | |||
extreme wet–dry surface | |||
conditions. | |||
Dew-point | Same as STIC1.0 | is estimated through numerical | |
temperature at | and are the slopes of | iteration by inverting the aerodynamic | |
source/sink | saturation vapor pressures at | equation of (as described in Appendix A2). | |
height () | temperatures, approximated at | ||
and , respectively. | |||
Surface | As a stand-alone equation, | Same as STIC1.0 | A feedback of into is introduced |
moisture | without any feedback to . | and is iteratively estimated after | |
availability | estimating (as described in Appendix A2). | ||
() | |||
Priestley– | As fixed parameter (1.26). | A physical equation | A physical equation of is derived as a |
Taylor | of is derived as a | function of the conductances and is | |
parameter () | function of the | numerically estimated as a variable | |
conductances and | (Eq. A15) (as described in Appendix A2). | ||
is numerically | |||
estimated as a | |||
variable. |
Fundamental differences in the modeling principles between STIC1.2 and previous approaches for characterizing the biophysical controls on components.
Biophysical | Modeling principles | |
---|---|---|
states | Parametric modeling | STIC1.2 |
(Ma et al., 2015; Chen et al., 2011; | ||
Kumagai et al., 2004) | ||
Either is assumed to be the momentum | Analytically retrieved by solving “” state | |
conductance () or estimated as a sum of and | equations and “” unknowns, with explicit | |
quasi-laminar boundary-layer conductance (). | convective feedback and without any wind | |
speed () information. | ||
1/ 1/ 1/ | In a hallmark paper by Choudhury and | |
Monteith (1986), it is clearly stated that | ||
Nusselt number, leaf dimension, thermal | “aerodynamic conductance determined by wind | |
conductivity of air in the boundary layer, , kinematic | speed and roughness is assumed to be | |
viscosity, Reynolds number | unaffected by buoyancy. Strictly, the | |
aerodynamic conductance should be replaced | ||
If is available from EC tower, it is directly | by a term that accounts for radiative as well | |
used, otherwise is parametrized using Monin– | as convective heat transfer”. The role of is | |
Obukhov Similarity Theory (MOST). | associated with the role of convection | |
(Choudhury and Monteith, 1986) according to | ||
Disadvantages: (1) MOST is only valid for an | the surface energy balance principle as | |
extended, uniform, and flat surface (Foken, | reflected in the derivation of Eq. (A4). Wind | |
2006). MOST tends to fail over rough surfaces | is generated as a result of the differences in | |
due to breakdown of the similarity relationships | atmospheric pressure that is a result of | |
for heat and water vapor transfer in the roughness | uneven surface radiative heating. Therefore, | |
sub-layer, which results in an underestimation of | the aerodynamic conductance (and wind as | |
the “true” by a factor of 1–3 (Thom et al., 1975; | well) is an effect of net radiative heating and | |
Chen and Schwerdtfeger, 1989; Simpson et al., | there should be a physical relationship between | |
1998; Holwerda et al., 2012). (2) In the state-of-the- | these two. | |
art modeling, the parametric sub-models | ||
are stand-alone and empirical, and do not provide | Advantages: (1) STIC1.2 consists of a | |
any feedback to , aerodynamic temperature | feedback describing the relationship between | |
(), and aerodynamic vapor pressures ( and | and , coupled with canopy–atmosphere | |
). (3) Additional challenges in grid-scale or | components relating to and . | |
spatial-scale estimation are the requirements | (2) Supports the findings of Villani et al. (2003) | |
of numerous site-specific parameters (e.g., | which stated that during unstable surface-layer | |
vegetation height, measurement height, | conditions the major source of net available | |
vegetation roughness, leaf size, soil roughness), | energy is located at the canopy top and drives | |
and coefficients needed to correct the | the convective motion in the layers above. | |
atmospheric stability conditions (Raupach, | ||
1998). | ||
(a) If measurements are available from the EC | Analytically retrieved by solving “” state | |
towers, is estimated by inverting the PM | equations and “” unknowns where physical | |
equation. None of these approaches allow | feedbacks of , soil moisture, and vapor | |
independent quantification of biophysical | pressure deficit are embedded (as explained in | |
controls of as is constrained by itself. | STIC1.2 equations in the Appendix). | |
(b) Sometimes is modeled either by coupled | ||
leaf-scale photosynthesis models (Ball et al., | ||
1987; Leuning, 1995) or is estimated from | ||
stand-alone empirical models (Jarvis, 1976) |
The Supplement related to this article is available online at
Acknowledgements
This study was funded by the Luxembourg Institute of Science and Technology (LIST). The developed modeling framework contributes to the Catchments As Organized Systems (CAOS) Phase-2 research group (FOR 1598) funded by the German Science Foundation (DFG) and to the HiWET (High-resolution modelling and monitoring of Water and Energy Transfers in wetland ecosystems) consortium funded by BELSPO and FNR. We sincerely thank Andrew Jarvis (Lancaster University, UK), Monica Garcia (Technical University of Denmark, Denmark), and Georg Wohlfahrt (University of Innsbruck, Austria) for very helpful discussions and edits of the manuscript. We are grateful to all Brazilian and international collaborators and all the funding agencies that have contributed to the Large-scale Biosphere Atmosphere Experiment in Amazônia (LBA). The authors are indebted to Pavel Kabat, Antônio Ocimar Manzi, David R. Fitzjarrald, Julio Tota, Humberto Ribeiro da Rocha, Michael Goulden, Maarten J. Waterloo, and Luiz Martinelli for planning, coordinating, conducting, and evaluating the eddy covariance, meteorological, and leaf–gas exchange measurements at the LBA sites. We are particularly grateful to all field technicians whose hard work was the key ingredient to establishing the quality of the datasets used in this paper. The authors declare no conflict of interest. Darren T. Drewry acknowledges the support of the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. Edited by: P. Saco Reviewed by: four anonymous referees
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2016. This work is published under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Canopy and aerodynamic conductances (
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details






1 Department of Environmental Research and Innovation, Luxembourg Institute of Science and Technology (LIST), L4422, Belvaux, Luxembourg
2 Department of Science and Environment, Roskilde University, Roskilde, Denmark
3 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, 91109, USA; Joint Institute for Regional Earth System Science and Engineering, University of California, Los Angeles, California, USA
4 Instituto Nacional de Pesquisas Espaciais (INPE), Centro de Ciência do Sistema Terrestre, São José dos Campos, SP, Brazil
5 Wageningen Environmental Research (ALTERRA), Wageningen, the Netherlands
6 Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA), Belém, PA, Brazil
7 Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ, USA
8 Department of Biology, University of Utah, Salt Lake City, UT, USA
9 Faculdade de Filosofia Ciências e Letras de Ribeirão Preto, Universidade de São Paulo (USP), São Paulo, SP, Brazil
10 Centro Nacional de Monitoramento e Alertas de Desastres Naturais, São Paulo, SP, Brazil
11 Department of Earth and Planetary Science, Harvard University, Cambridge, MA, USA