Advances in Meteorology
Volume 2016 (2016), Article ID 1835487, 9 pages
http://dx.doi.org/10.1155/2016/1835487

Research Article

Relationship between Evapotranspiration and Land Surface Temperature under Energy- and Water-Limited Conditions in Dry and Cold Climates

Zhigang Sun,1,2 Qinxue Wang,³ Ochirbat Batkhishig,⁴ and Zhu Ouyang1,2

¹Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
²Yucheng Comprehensive Experimental Station, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
³Center for Regional Environmental Research, National Institute for Environmental Studies, Tsukuba 305-8506, Japan
⁴Institute of Geography, Mongolian Academy of Sciences, 14192 Ulaanbaatar, Mongolia

Received 24 May 2015; Revised 29 September 2015; Accepted 7 October 2015

Academic Editor: Ke Zhang

Copyright © 2016 Zhigang Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Remotely sensed land surface temperature- (LST-) dependent evapotranspiration (ET) models and vegetation index- (VI-) LST methods may not be suitable for ET estimation in energy-limited cold areas. In this study, the relationship of ET to LST was simulated using the process-based Simultaneous Heat and Water (SHAW) model for energy- and water-limited conditions in Mongolia, to understand the differences in ET processes under these two limiting conditions in dry and cold climates. Simulation results from the SHAW model along with ground observational data showed that ET and LST have a positive relationship when air temperature () is less than or equal to the temperature () above which plants transpire and have a negative relationship when is greater than under the energy-limited condition. However, ET and LST maintain a negative relationship with changes in under the water-limited condition. The differences in the relationship between ET and LST under the energy-limited and water-limited conditions could be attributed to plant transpiration and energy storage in moist/watered soil and plants. This study suggests that different strategies should be used to estimate ET under the energy-limited condition in dry and cold climates.

1. Introduction

Terrestrial evapotranspiration (ET), defined as the loss of water from the land surface to the atmosphere, is a key process in water cycles [1, 2]. ET closely relates to greenhouse gas efflux and production [3], plant growth [4–6], and droughts [7, 8]. Therefore, it is essential to understand the ET phenomenon over wide ranges of space and time for water resources management and climate change studies.

Regional or global ET maps are often obtained through satellite remote sensing [9, 10]. Among remote sensing-based ET estimation methods, remotely sensed land surface temperature- (LST) dependent methods, such as the vegetation index (VI) and LST scattered plot (VI-LST), and models that calculate ET as a residual of the land surface energy balance equation have been widely used. The LST-dependent ET models include SEBAL (remote sensing Surface Energy Balance Algorithm for Land) [11], Sim-ReSET (Simple Remote Sensing EvapoTranspiration model) [9], SEBI (Surface Energy Balance Index) and S-SEBI (Simplified-SEBI) [12], SEBS (Surface Energy Balance System) [13], and METRIC (Mapping EvapoTranspiration at high Resolution with Internalized Calibration) [14, 15], among others. The VI-LST methods use the VI-LST triangle/trapezoidal feature space to derive an index, and this index is used to partition sensible heat and latent heat fluxes from available energy (net radiation minus soil heat flux) [16–19]. Both LST-dependent ET models and VI-LST methods assume that ET can cool land surfaces under the condition of homogeneous atmospheric forcing. In other words, land surfaces with larger ET have lower LST while air temperature (), air humidity, wind speed, and solar radiation remain homogeneous over land surfaces of interest. However, this assumption may be invalid in high latitudes (>50°) and cold areas [2], where ET commonly exhibits a positive relationship with LST [2, 20]. A positive relationship between LST and VI has also been observed in Alaska tundra ecosystems [21] and in North America above 45°N [22]. These findings indicate that both LST-dependent ET models and VI-LST methods may not be suitable for ET estimation in the energy-limited cold areas. A study to estimate water deficit (defined as the ratio of actual to potential ET) using the VI-LST method in southern Spain also suggests that the VI-LST method should not be used to calculate ET under energy-limited conditions [19].

Reasons for the incompetence of LST-dependent ET models and VI-LST methods in energy-limited cold areas have been seldom investigated so far. In this study, the relationship between ET and LST was investigated using the process-based Simultaneous Heat and Water (SHAW) model for energy- and water-limited conditions in Mongolia, to understand the difference in ET processes under such limiting conditions in dry and cold climates.

2. Methodology

2.1. SHAW Model

The Simultaneous Heat and Water (SHAW) model, a one-dimensional model that simulates heat, water, and solute transfer within the atmosphere-plant-soil system, was used to simulate ET and LST in this study. The SHAW model, developed by USDA Agricultural Research Service, has a detailed solution to simulate soil freezing and thawing and a sophisticated approach to simulate transpiration from a plant canopy and evaporation from soil [23, 24]. The SHAW model has been documented that it is capable to simulate heat and water movement through plant cover, snow, residue, and soil for investigating climate and anthropogenic effects on soil freezing, snowmelt, runoff, soil temperature, water, evaporation, and transpiration [23, 24]. Therefore, the calibrated SHAW model could be competent at numerically simulating the relationship between ET and LST.

Vegetation, climate, and soil data are required to run the SHAW model. These inputs include initial conditions for snow, soil temperature, and water content profiles; daily or hourly weather data; and parameters describing the vegetative cover, snow, residue, and soil. General site information is also needed to initialize the model, including the slope, aspect, latitude, and surface roughness parameters. Detailed information on these input parameters is listed in Table 1.

Table 1: Main inputs to the SHAW model.

Date type
Data item
Unit
Data source

Weather data
Air temperature
°C
AMS

Air humidity
%
AMS

Precipitation
Inch
AMS

Wind speed
mph
AMS

Solar radiation
W/m²
AMS

Soil data
Bulk density
kg/m³
Soil sample

Organic matter
%
Soil sample

Sand
%
Soil sample

Silt
%
Soil sample

Clay
%
Soil sample

Soil temperature
°C
AMS

Soil moisture
m³/m³
AMS

Vegetation data
Vegetation canopy height
m
Vegetation sample

Dry matter
kg/m²
Vegetation sample

Leaf area index
—
Vegetation sample

Root depth
m
Vegetation sample

General site information
Slope
%
Field survey

Aspect
Degree
Field survey

Latitude
Degree
GPS

Elevation
m
GPS

2.2. Data

Weather Data. Hourly air temperature (), wind speed, air humidity, precipitation, and solar radiation data were collected from an auto-meteorological station (AMS; 47.75°N 107.33°E) located near Ulaanbaatar in Mongolia. A dry and cold climate dominates the region around the AMS, with an annual mean air temperature of −2.8°C and annual total precipitation of 260mm. Detailed information regarding the AMS can be found in Sun et al. [25].

Soil Data. Data for organic matter content, bulk density, and soil component ratio of sand, silt, and clay were obtained from soil samples that were analyzed in the lab. Soil temperature and moisture profile data were obtained from the AMS.

Vegetation Data. Grass height, dry matter, leaf area index (LAI), and root depth data were obtained by measuring grass samples.

2.3. Strategies for ET and LST Simulations

In this study, nine scenarios of and four scenarios of topsoil water content (SWC, 0–25cm) were combined to simulate ET and LST for energy- and water-limited conditions. was set at −1–15°C to match the local cold climate, and SWC was set at 0.1–0.4 to cover the soil moisture range from dry to wet. Combinations of and SWC could thereby simulate a wide range of energy and water conditions under which ET and LST are simulated using the SHAW model (Figure 1). The scenario of SWC = 0.1 represents the water-limited condition where ET is mainly limited by soil water shortage. The scenario of SWC = 0.4 with low values represents the energy-limited condition where ET is mainly limited by low incoming solar radiation. Since low incoming solar radiation leads to low , the proxy of low along with saturated soil water content could be used to define the energy-limited condition. One combination of and SWC along with other inputs was used to run the model. For each model run, the SWC data were used to initialize the moisture condition for the soil layer of 0–25cm, and data along with other weather data representing typical weather conditions at the AMS site in early spring were selected from the AMS dataset. A typical record of and other weather data was repeated 84 times to simulate the condition of homogeneous atmospheric forcing. Therefore, ET and LST values from each run of the SHAW model represent the water and heat conditions of land surfaces under homogeneous atmospheric forcing. Outputted ET and LST values were plotted and analyzed using the linear regression method, and the slope and coefficient of determination () of the linear regression were used to understand the relationship of ET to LST under different energy and moisture conditions.

Figure 1: Strategies to simulate ET and LST using the SHAW model.

Before numerical simulations, intensive daily observational data during the period of August 2008 to August 2013 were used to calibrate the SHAW model by forcing the simulated soil temperature and water content to be consistent with observed values, respectively.

3. Results and Discussion

When the SHAW model was calibrated, simulations of nine scenarios of (−1 to 15°C) and four scenarios of topsoil water content (SWC = 0.1, 0.2, 0.3, and 0.4) were conducted to estimate ET and LST. For each scenario, the relationship of ET and LST for SWC = 0.2, 0.3, and 0.4 is similar, but for SWC = 0.1 the relationship is quite different. Below are detailed results and discussion regarding model calibration and the relationship of ET and LST under energy- and water-limited conditions.

3.1. Calibration Results for the SHAW Model

The comparison of simulated soil temperature and soil water content against observed values at four depths was used to evaluate the performance of the calibrated SHAW model. Figure 2 shows that simulated soil temperature matches observed one well in magnitude and seasonal variation. A linear analysis was conducted for each comparison. The coefficient of determination () is greater than 0.9, and the mean absolute error (MAE) is less than 1.4°C at four depths. Figure 3 shows that simulated soil water content also matches observed one well in magnitude. The mean absolute error (MAE) is less than 0.015 at four depths. The simulated soil water content can capture the soil water recharge from precipitation as shown in Figure 3(a). Soil temperature and soil water content are two key variables of water and thermal processes in soil. Therefore, the consistence of soil temperature and soil water content between simulated and observed values could indicate that the calibrated SHAW model could perform well to investigate the relationship of ET and LST.

(a)

(b)

(c)

(d)

Figure 2: Comparison of simulated and observed soil temperature at the depths of 0.5m, 1m, 1.5m, and 2m, respectively. A linear analysis was conducted for each comparison, where is the coefficient of determination and MAE is the mean absolute error.

(a)

(b)

(c)

(d)

Figure 3: Comparison of simulated and observed soil water content at the depths of 0.5m, 1m, 1.5m, and 2m, respectively. A linear analysis was conducted for each comparison, where is the coefficient of determination and MAE is the mean absolute error.

3.2. Relationship of ET and LST under the Energy-Limited Condition in a Cold and Dry Climate

The relationship of ET and LST under the soil moisture condition of SWC = 0.4 exhibited two patterns for nine simulations when varied from −1 to 15°C. As shown in Figure 4, simulations with from −1 to 5°C show that ET increases as LST increases, but this trend becomes progressively weaker with increases in . The slope of linear regression analysis gradually decreases from 0.0059 to 0.0037, and the value of also decreases from 0.8356 to 0.1153. The relationship of ET and LST abruptly changes at °C, and the slope of linear regression analysis becomes −0.0329 (). This negative trend becomes stronger as increases from 7°C to 15°C; the slope of linear regression analysis decreases from −0.0329 to −0.0636, and the value of increases from 0.1525 to 0.2946.

(a)°C

(b)°C

(c)°C

(d)°C

(e)°C

(f)°C

(g)°C

(h)°C

(i)°C

Figure 4: ET and LST under nine scenarios of air temperature and four scenarios of soil water content.

Figures 4(a)–4(e) show that higher land surface temperature enhances ET when air is cold ( ≤ 5°C) and soil water is not limited. These results suggest that some caution should be noted when ET is estimated from LST-dependent ET models and VI-LST methods. These methods for ET estimation assume that ET cools terrestrial surfaces through vegetation transpiration and soil (or water surface) evaporation and that moist/watered surfaces have relatively low values of LST. However, these assumptions are invalid under the energy-limited condition (°C, SWC = 0.4).

3.3. Relationship of ET and LST under the Water-Limited Condition in a Cold and Dry Climate

Compared with the energy-limited condition, the relationship of ET and LST under the water-limited condition (SWC = 0.1) is relatively straightforward. As shown in Figure 4, the trend of decreasing ET with increasing LST becomes more distinct as is increased from −1 to 15°C. The slope of linear regression analysis gradually decreases from −0.0022 to −0.0097, and the value of increases from 0.1323 to 0.6552. Results from Figure 4 indicate that ET has an apparent cooling effect even when is less than or equal to 5°C under the water-limited condition (SWC = 0.1). Under such a condition, the assumption of ET cooling appears to be valid, and ET can be estimated from LST-dependent ET models and VI-LST methods.

3.4. Effect of Plant Transpiration on the Relationship of ET and LST under the Energy- and Water-Limited Conditions in a Cold and Dry Climate

Plant transpiration has a clear impact on the relationship of ET and LST under the energy-limited condition. As shown in Figure 5, the relationship of ET and LST becomes negative when is greater than 5°C, which for this study was set as the temperature above which plants start to transpire. It indicates that soil evaporation has no clear cooling function when °C, whereas plant transpiration above 5°C can cool land surfaces under the energy-limited condition.

Figure 5: Slope of linear regression for ET and LST with nine scenarios of air temperature and two scenarios of soil water content.

Plant transpiration has less impact on the relationship of ET and LST under the water-limited condition. As shown in Figure 5, the relationship of ET and LST remains negative under the water-limited condition, and this negative relationship becomes more pronounced when °C (Figure 6). This indicates that soil evaporation has a cooling function even when °C and that plant transpiration above 5°C can enhance the cooling effect under the energy-limited condition.

Figure 6: of linear regression for ET and LST with nine scenarios of air temperature and two scenarios of soil water content.

Summarizing the inferences drawn from Figures 5 and 6, plant transpiration shows a cooling effect under both energy- and water-limited conditions, while soil evaporation has a cooling effect under the water-limited condition only. These results suggest that care should be taken with ET estimation from LST-dependent ET models and VI-LST methods at temperatures below which plants do not transpire (°C in this study).

3.5. Comparison of ET and LST Relationships between Energy- and Water-Limited Conditions

Under the water-limited condition, the relationship of ET and LST is always negative across the range from −1 to 15°C. With the assumption of an ET cooling function under the dry and cold condition, therefore, LST-dependent ET models and VI-LST methods can be used for ET estimation. Under the energy-limited condition, however, the ET estimating methods mentioned above do not work and the assumption of an ET cooling function is invalid when °C. The reason for this might be that the energy storage of moist/watered soil and plants is larger than that of dry soil and plants, and it is comparable to available energy for moist/watered soil and plants. This suggests that ET estimation methods based on the land surface energy balance equation might be useful under the energy-limited condition if the energy storage by moist/watered soil and plants is accounted for in the land surface energy balance equation.

4. Conclusions

The energy- and water-limited ET processes should be considered and included in the methods used for ET estimation in cold and dry climates. In this study, the process-based SHAW model along with ground observational data was used to simulate ET and LST in order to understand the relationship of ET and LST under energy- and water-limited conditions. Simulation results indicated that ET and LST have a positive relationship when °C and a negative relationship when °C under the energy-limited condition (SWC = 0.4). However, ET and LST maintain a negative relationship under the water-limited condition (SWC = 0.1). Plant transpiration and energy storage in moist/watered soil and plants can potentially explain the differences in the relationship of ET and LST simulated under the energy-limited and water-limited conditions. Plant transpiration likely affects the relationship of ET and LST due to its relatively strong cooling effect under the energy-limited condition. The energy storage of moist/watered soil and plants, which is comparable to the available energy, has a relatively large contribution to ET when °C. This study suggests that different strategies should be used to estimate ET under the energy-limited condition in a dry and cold climate.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was supported by the 100 Talents Program of the Chinese Academy of Sciences. This research was also supported by the Grant-in-Aid for Encouragement, Center for Regional Environmental Research, National Institute for Environmental Studies, Japan, and by the Environment Research and Technology Development Fund, the Ministry of Environment, Japan.

References

S. Jasechko, Z. D. Sharp, J. J. Gibson, S. J. Birks, Y. Yi, and P. J. Fawcett, “Terrestrial water fluxes dominated by transpiration,” Nature, vol. 496, no. 7445, pp. 347–350, 2013.

K. Wang and R. E. Dickinson, “A review of global terrestrial evapotranspiration: observation, modeling, climatology, and climatic variability,” Reviews of Geophysics, vol. 50, no. 2, pp. 1–54, 2012.

J. Balogh, S. Fóti, K. Pintér et al., “Soil CO₂ efflux and production rates as influenced by evapotranspiration in a dry grassland,” Plant and Soil, vol. 388, no. 1-2, pp. 157–173, 2015.

A. Harper, I. T. Baker, A. S. Denning, D. A. Randall, D. Dazlich, and M. Branson, “Impact of evapotranspiration on dry season climate in the Amazon forest,” Journal of Climate, vol. 27, no. 2, pp. 574–591, 2014.

J. B. Wu, Y. L. Jing, D. X. Guan et al., “Controls of evapotranspiration during the short dry season in a temperate mixed forest in Northeast China,” Ecohydrology, vol. 6, no. 5, pp. 775–782, 2013.

M. C. R. Alberto, J. R. Quilty, R. J. Buresh et al., “Actual evapotranspiration and dual crop coefficients for dry-seeded rice and hybrid maize grown with overhead sprinkler irrigation,” Agricultural Water Management, vol. 136, pp. 1–12, 2014.

M. C. Anderson, C. Hain, B. Wardlow, A. Pimstein, J. R. Mecikalski, and W. P. Kustas, “Evaluation of drought indices based on thermal remote sensing of evapotranspiration over the continental United States,” Journal of Climate, vol. 24, no. 8, pp. 2025–2044, 2011.

A. J. Teuling, A. F. Van Loon, S. I. Seneviratne et al., “Evapotranspiration amplifies European summer drought,” Geophysical Research Letters, vol. 40, no. 10, pp. 2071–2075, 2013.

Z. G. Sun, Q. X. Wang, B. Matsushita, T. Fukushima, Z. Ouyang, and M. Watanabe, “Development of a simple remote sensing EvapoTranspiration model (Sim-ReSET): algorithm and model test,” Journal of Hydrology, vol. 376, no. 3-4, pp. 476–485, 2009.

Z.-L. Li, R. L. Tang, Z. M. Wan et al., “A review of current methodologies for regional evapotranspiration estimation from remotely sensed data,” Sensors, vol. 9, no. 5, pp. 3801–3853, 2009.

W. G. M. Bastiaanssen, M. Menenti, R. A. Feddes, and A. A. M. Holtslag, “A remote sensing surface energy balance algorithm for land (SEBAL): 1. Formulation,” Journal of Hydrology, vol. 212, no. 1–4, pp. 198–212, 1998.

G. J. Roerink, Z. Su, and M. Menenti, “S-SEBI: a simple remote sensing algorithm to estimate the surface energy balance,” Physics and Chemistry of the Earth, Part B: Hydrology, Oceans and Atmosphere, vol. 25, no. 2, pp. 147–157, 2000.

Z. Su, “The Surface Energy Balance System (SEBS) for estimation of turbulent heat fluxes,” Hydrology and Earth System Sciences, vol. 6, no. 1, pp. 85–99, 2002.

R. G. Allen, M. Tasumi, A. Morse et al., “Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC)—applications,” Journal of Irrigation and Drainage Engineering, vol. 133, no. 4, pp. 395–406, 2007.

R. G. Allen, M. Tasumi, and R. Trezza, “Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC)—model,” Journal of Irrigation and Drainage Engineering, vol. 133, no. 4, pp. 380–394, 2007.

L. Jiang and S. Islam, “Estimation of surface evaporation map over southern Great Plains using remote sensing data,” Water Resources Research, vol. 37, no. 2, pp. 329–340, 2001.

T. Carlson, “An overview of the ‘triangle method’ for estimating surface evapotranspiration and soil moisture from satellite imagery,” Sensors, vol. 7, no. 8, pp. 1612–1629, 2007.

R. L. Tang, Z.-L. Li, and B. H. Tang, “An application of the Ts-VI triangle method with enhanced edges determination for evapotranspiration estimation from MODIS data in arid and semi-arid regions: implementation and validation,” Remote Sensing of Environment, vol. 114, no. 3, pp. 540–551, 2010.

M. Garcia, N. Fernández, L. Villagarcía, F. Domingo, J. Puigdefábregas, and I. Sandholt, “Accuracy of the Temperature-Vegetation Dryness Index using MODIS under water-limited vs. Energy-limited evapotranspiration conditions,” Remote Sensing of Environment, vol. 149, pp. 100–117, 2014.

R. R. Nemani, C. D. Keeling, H. Hashimoto et al., “Climate-driven increases in global terrestrial net primary production from 1982 to 1999,” Science, vol. 300, no. 5625, pp. 1560–1563, 2003.

A. Hope, R. Engstrom, and D. Stow, “Relationship between AVHRR surface temperature and NDVI in Arctic tundra ecosystems,” International Journal of Remote Sensing, vol. 26, no. 8, pp. 1771–1776, 2005.

A. Karnieli, N. Agam, R. T. Pinker et al., “Use of NDVI and land surface temperature for drought assessment: merits and limitations,” Journal of Climate, vol. 23, no. 3, pp. 618–633, 2010.

G. N. Flerchinger, “Sensitivity of soil freezing simulated by the SHAW model,” Transactions of the ASAE, vol. 34, no. 6, pp. 2381–2389, 1991.

G. N. Flerchinger, T. G. Caldwell, J. Cho, and S. P. Hardegree, “Simultaneous heat and water (SHAW) model: model use, calibration, and validation,” Transactions of the ASABE, vol. 55, no. 4, pp. 1395–1411, 2012.

Z. G. Sun, Q. X. Wang, Q. G. Xiao, O. Batkhishig, and M. Watanabe, “Diverse responses of remotely sensed grassland phenology to interannual climate variability over frozen ground regions in Mongolia,” Remote Sensing, vol. 7, no. 1, pp. 360–377, 2015.