Skip to main content
Log in

Evaluation of TSEB turbulent fluxes using different methods for the retrieval of soil and canopy component temperatures from UAV thermal and multispectral imagery

  • Original Paper
  • Published:
Irrigation Science Aims and scope Submit manuscript

Abstract

The thermal-based Two-Source Energy Balance (TSEB) model partitions the evapotranspiration (ET) and energy fluxes from vegetation and soil components providing the capability for estimating soil evaporation (E) and canopy transpiration (T). However, it is crucial for ET partitioning to retrieve reliable estimates of canopy and soil temperatures and net radiation, as the latter determines the available energy for water and heat exchange from soil and canopy sources. These two factors become especially relevant in row crops with wide spacing and strongly clumped vegetation such as vineyards and orchards. To better understand these effects, very high spatial resolution remote-sensing data from an unmanned aerial vehicle were collected over vineyards in California, as part of the Grape Remote sensing and Atmospheric Profile and Evapotranspiration eXperiment and used in four different TSEB approaches to estimate the component soil and canopy temperatures, and ET partitioning between soil and canopy. Two approaches rely on the use of composite \(T_\mathrm{rad}\), and assume initially that the canopy transpires at the Priestley–Taylor potential rate. The other two algorithms are based on the contextual relationship between optical and thermal imagery partition \(T_\mathrm{rad}\) into soil and canopy component temperatures, which are then used to drive the TSEB without requiring a priori assumptions regarding initial canopy transpiration rate. The results showed that a simple contextual algorithm based on the inverse relationship of a vegetation index and \(T_\mathrm{rad}\) to derive soil and canopy temperatures yielded the closest agreement with flux tower measurements. The utility in very high-resolution remote-sensing data for estimating ET and E and T partitioning at the canopy level is also discussed.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9

Similar content being viewed by others

Notes

  1. The use of trade, firm, or corporation names in this article is for the information and convenience of the reader. Such use does not constitute official endorsement or approval by the US Department of Agriculture or the Agricultural Research Service of any product or service to the exclusion of others that may be suitable

References

  • Aboutalebi A, Torres-Rua A, Nieto H, Kustas W (2018) Assessment of different methods for the shadows detection in high-resolution imagery and shadows impact on calculation of NDVI, LAI, and evapotranspiration. Irrig Sci (this issue)

  • Agam N, Alfieri J, Kustas W, Jones S, McKee L, Prueger J (2018) Spatial variability in soil heat flux from a detailed soil heat flux plates array. Irrig Sci (this issue)

  • Alfieri J, Kustas W, Gao F, Prueger J, Nieto H, Hipps L (2018) Influence of wind direction on the effective surface roughness of vineyards. Irrig Sci (this issue)

  • Anderson MC, Norman JM, Kustas WP, Houborg R, Starks PJ, Agam N (2008) A thermal-based remote sensing technique for routine mapping of land-surface carbon, water and energy fluxes from field to regional scales. Remote Sens Environ 112(12):4227–4241. https://doi.org/10.1016/j.rse.2008.07.009

    Article  Google Scholar 

  • Anderson MC, Kustas WP, Norman JM, Hain CR, Mecikalski JR, Schultz L, González-Dugo MP, Cammalleri C, d’Urso G, Pimstein A, Gao F (2011) Mapping daily evapotranspiration at field to continental scales using geostationary and polar orbiting satellite imagery. Hydrol Earth Syst Sci 15(1):223–239. https://doi.org/10.5194/hess-15-223-2011

    Article  Google Scholar 

  • Bellvert J, Zarco-Tejada P, Marsal J, Girona J, González-Dugo V, Fereres E (2016) Vineyard irrigation scheduling based on airborne thermal imagery and water potential thresholds. Aust J Grape Wine Res 22(2):307–315. https://doi.org/10.1111/ajgw.12173

    Article  Google Scholar 

  • Berni JAJ, Zarco-Tejada PJ, Suarez L, Fereres E (2009) Thermal and narrowband multispectral remote sensing for vegetation monitoring from an unmanned aerial vehicle. IEEE Trans Geosci Remote Sens 47(3):722–738. https://doi.org/10.1109/TGRS.2008.2010457

    Article  Google Scholar 

  • Brutsaert W (1999) Aspects of bulk atmospheric boundary layer similarity under free-convective conditions. Rev Geophys 37(4):439–451

    Article  Google Scholar 

  • Brutsaert W (2005) Hydrology: an introduction. Cambridge University Press, Cambridge

    Book  Google Scholar 

  • Cammalleri C, Anderson MC, Gao F, Hain CR, Kustas WP (2013) A data fusion approach for mapping daily evapotranspiration at field scale. Water Resour Res 49(8):4672–4686. https://doi.org/10.1002/wrcr.20349

    Article  Google Scholar 

  • Cammalleri C, Anderson M, Gao F, Hain C, Kustas W (2014) Mapping daily evapotranspiration at field scales over rainfed and irrigated agricultural areas using remote sensing data fusion. Agric For Meteorol 186:1–11. https://doi.org/10.1016/j.agrformet.2013.11.001

    Article  Google Scholar 

  • Campbell G (1990) Derivation of an angle density function for canopies with ellipsoidal leaf angle distributions. Agric For Meteorol 49(3):173–176. https://doi.org/10.1016/0168-1923(90)90030-A

    Article  Google Scholar 

  • Campbell GS (1986) Extinction coefficients for radiation in plant canopies calculated using an ellipsoidal inclination angle distribution. Agric For Meteorol 36(4):317–321

    Article  Google Scholar 

  • Campbell GS, Norman JM (1998) An introduction to environmental biophysics, 2nd edn. Springer, New York

    Book  Google Scholar 

  • Chehbouni A, Nouvellon Y, Lhomme J, Watts C, Boulet G, Kerr Y, Moran M, Goodrich D (2001) Estimation of surface sensible heat flux using dual angle observations of radiative surface temperature. Agric For Meteorol 108(1):55–65

    Article  Google Scholar 

  • Colaizzi P, Evett S, Howell T, Li F, Kustas W, Anderson M (2012a) Radiation model for row crops: I. Geometric view factors and parameter optimization. Agron J 104(2):225–240

    Article  Google Scholar 

  • Colaizzi P, Agam N, Tolk J, Evett S, Howell T, Gowda P, O’Shaughnessy S, Kustas W, Anderson M (2014) Two-source energy balance model to calculate E, T, and ET: Comparison of Priestley-Taylor and Penman-Monteith formulations and two time scaling methods. Trans ASABE 57:479–498

    Google Scholar 

  • Colaizzi PD, Kustas WP, Anderson MC, Agam N, Tolk JA, Evett SR, Howell TA, Gowda PH, O’Shaughnessy SA (2012b) Two-source energy balance model estimates of evapotranspiration using component and composite surface temperatures. Adv Water Resour 50:134–151. https://doi.org/10.1016/j.advwatres.2012.06.004

    Article  Google Scholar 

  • Colaizzi PD, Schwartz RC, Evett SR, Howell TA, Gowda PH, Tolk JA (2012c) Radiation model for row crops: II. Model evaluation. Agron J 104(2):241–255. https://doi.org/10.2134/agronj2011.0083

    Article  Google Scholar 

  • Coops NC, Hilker T, Wulder MA, St-Onge B, Newnham G, Siggins A, Trofymow JAT (2007) Estimating canopy structure of Douglas-fir forest stands from discrete-return LiDAR. Trees 21(3):295. https://doi.org/10.1007/s00468-006-0119-6

    Article  Google Scholar 

  • Detto M, Montaldo N, Albertson JD, Mancini M, Katul G (2006) Soil moisture and vegetation controls on evapotranspiration in a heterogeneous mediterranean ecosystem on Sardinia. Italy. Water Resour Res 42(8):w08419. https://doi.org/10.1029/2005WR004693

    Article  Google Scholar 

  • Gao F, Kustas WP, Anderson MC (2012) A data mining approach for sharpening thermal satellite imagery over land. Remote Sens 4(11):3287. https://doi.org/10.3390/rs4113287

    Article  Google Scholar 

  • Goudriaan J (1977) Crop micrometeorology: a simulation stud. Center for Agricultural Publications and Documentation, Wageningen, Tech. rep

  • Guzinski R, Anderson MC, Kustas WP, Nieto H, Sandholt I (2013) Using a thermal-based two source energy balance model with time-differencing to estimate surface energy fluxes with day-night MODIS observations. Hydrol Earth Syst Sci 17(7):2809–2825. https://doi.org/10.5194/hess-17-2809-2013

    Article  Google Scholar 

  • Guzinski R, Nieto H, Jensen R, Mendiguren G (2014) Remotely sensed land-surface energy fluxes at sub-field scale in heterogeneous agricultural landscape and coniferous plantation. Biogeosciences 11(18):5021–5046. https://doi.org/10.5194/bg-11-5021-2014

    Article  Google Scholar 

  • Guzinski R, Nieto H, Stisen S, Fensholt R (2015) Inter-comparison of energy balance and hydrological models for land surface energy flux estimation over a whole river catchment. Hydrol Earth Syst Sci 19(4):2017–2036. https://doi.org/10.5194/hess-19-2017-2015

    Article  Google Scholar 

  • Hillel D (1998) Environmental soil physics. Academic Press, New York

    Google Scholar 

  • Hoffmann H, Nieto H, Jensen R, Guzinski R, Zarco-Tejada P, Friborg T (2016) Estimating evaporation with thermal UAV and two-source energy balance models. Hydrol Earth Syst Sci 20(2):697–713. https://doi.org/10.5194/hess-20-697-2016

    Article  Google Scholar 

  • Knipper KR, Kustas WP, Anderson MC, Alfieri JG, Prueger JH, Hain CR, Gao F, Yang Y, McKee LG, Nieto H, Hipps LE, Alsina MM, Sanchez6 L (2018) Evapotranspiration estimates derived using thermal-based satellite remote sensing and data fusion for irrigation management in california vineyards. Irrig Sci (this issue)

  • Kondo J, Ishida S (1997) Sensible heat flux from the earth’s surface under natural convective conditions. J Atmos Sci 54(4):498–509

    Article  Google Scholar 

  • Kustas W, Anderson M (2009) Advances in thermal infrared remote sensing for land surface modeling. Agric For Meteorol 149(12):2071–2081 (Environmental Biophysics - Tribute to John Norman)

  • Kustas W, Norman J (2000) A two-source energy balance approach using directional radiometric temperature observations for sparse canopy covered surfaces. Agron J 92(5):847–854

    Article  Google Scholar 

  • Kustas W, Alfieri J, Nieto H, Gao F, Anderson M (2018) Utility of the two-source energy balance model TSEB in vine and inter-row flux partitioning over the growing season. Irrig Sci. https://doi.org/10.1007/s00271-018-0586-8 (this issue)

    Article  PubMed  Google Scholar 

  • Kustas WP, Norman JM (1997) A two-source approach for estimating turbulent fluxes using multiple angle thermal infrared observations. Water Resour Res 33(6):1495–1508. https://doi.org/10.1029/97WR00704

    Article  Google Scholar 

  • Kustas WP, Norman JM (1999) Evaluation of soil and vegetation heat flux predictions using a simple two-source model with radiometric temperatures for partial canopy cover. Agric For Meteorol 94(1):13–29. https://doi.org/10.1016/S0168-1923(99)00005-2

    Article  Google Scholar 

  • Kustas WP, Alfieri JG, Anderson MC, Colaizzi PD, Prueger JH, Evett SR, Neale CM, French AN, Hipps LE, Chávez JL, Copeland KS, Howell TA (2012) Evaluating the two-source energy balance model using local thermal and surface flux observations in a strongly advective irrigated agricultural area. Adv Water Resour 50:120–133. https://doi.org/10.1016/j.advwatres.2012.07.005

    Article  Google Scholar 

  • Kustas WP, Nieto H, Morillas L, Anderson MC, Alfieri JG, Hipps LE, Villagarcía L, Domingo F, García M (2016) Revisiting the paper “using radiometric surface temperature for surface energy flux estimation in mediterranean drylands from a two-source perspective”. Remote Sens Environ 184:645–653

    Article  Google Scholar 

  • Kustas WP, Anderson MC, Alfieri JG, Knipper K, Torres-Rua A, Parry CK, Nieto H, Agam N, White A, Gao F, McKee L, Prueger JH, Hipps LE, Los S, Alsina M, Sanchez L, Sams B, Dokoozlian N, McKee M, Jones S, McElrone A, Heitman JL, Howard AM, Post K, Melton F, Hain C (2018) The grape remote sensing atmospheric profile and evapotranspiration eXperiment (GRAPEX). Bull Am Meteorol Soc. https://doi.org/10.1175/BAMS-D-16-0244.1 (in press)

    Article  Google Scholar 

  • Li F, Kustas WP, Prueger JH, Neale CM, Jackson TJ (2005) Utility of remote sensing-based two-source energy balance model under low-and high-vegetation cover conditions. J Hydrometeorol 6(6):878–891

    Article  Google Scholar 

  • Mallet C, Bretar F (2009) Full-waveform topographic lidar: state-of-the-art. ISPRS J Photogramm Remote Sens 64(1):1–16. https://doi.org/10.1016/j.isprsjprs.2008.09.007

    Article  Google Scholar 

  • Massman W (1987) A comparative study of some mathematical models of the mean wind structure and aerodynamic drag of plant canopies. Bound-Layer Meteorol 40(1):179–197

    Article  Google Scholar 

  • Massman W, Forthofer J, Finney M (2017) An improved canopy wind model for predicting wind adjustment factors and wildland fire behavior. Can J For Res 47(5):594–603. https://doi.org/10.1139/cjfr-2016-0354

    Article  Google Scholar 

  • Morillas L, García M, Nieto H, Villagarcía L, Sandholt I, González-Dugo M, Zarco-Tejada P, Domingo F (2013) Using radiometric surface temperature for surface energy flux estimation in Mediterranean drylands from a two-source perspective. Remote Sens Environ 136:234–246. https://doi.org/10.1016/j.rse.2013.05.010

    Article  Google Scholar 

  • Nieto H, Kustas W, Gao F, Alfieri J, Torres A, Hipps L (2018) Impact of different within-canopy wind attenuation formulations on modelling sensible heat flux using TSEB. Irrig Sci. https://doi.org/10.1007/s00271-018-0611-y

    Article  PubMed  Google Scholar 

  • Norman J, Kustas W, Prueger J, Diak G (2000) Surface flux estimation using radiometric temperature: a dual-temperature-difference method to minimize measurement errors. Water Resour Res 36(8):2263–2274

    Article  Google Scholar 

  • Norman JM, Kustas WP, Humes KS (1995) Source approach for estimating soil and vegetation energy fluxes in observations of directional radiometric surface temperature. Agric For Meteorol 77(3–4):263–293. https://doi.org/10.1016/0168-1923(95)02265-Y

    Article  Google Scholar 

  • Ortega-Farías S, Ortega-Salazar S, Poblete T, Kilic A, Allen R, Poblete-Echeverría C, Ahumada-Orellana L, Zuñiga M, Sepúlveda D (2016) Estimation of energy balance components over a drip-irrigated olive orchard using thermal and multispectral cameras placed on a helicopter-based unmanned aerial vehicle (UAV). Remote Sens 8(8):638. https://doi.org/10.3390/rs8080638

    Article  Google Scholar 

  • Parry C, Nieto H, Guillevic P, Agam N, Kustas B, Alfieri J, McKee L, McElrone A (2018) An intercomparison of radiation partitioning models in vineyard row structured canopies. Irrig Sci (this issue)

  • Pope G, Treitz P (2013) Leaf area index (LAI) estimation in boreal mixedwood forest of Ontario, Canada using light detection and ranging (LiDAR) and WorldView-2 imagery. Remote Sens 5(10):5040–5063. https://doi.org/10.3390/rs5105040

    Article  Google Scholar 

  • Priestley CHB, Taylor RJ (1972) On the assessment of surface heat flux and evaporation using large-scale parameters. Mon Weather Rev 100(2):81–92. https://doi.org/10.1175/1520-0493(1972)100%3c0081:OTAOSH%3e2.3.CO;2

    Article  Google Scholar 

  • Santanello J Jr, Friedl M (2003) Diurnal covariation in soil heat flux and net radiation. J Appl Meteorol 42(6):851–862

    Article  Google Scholar 

  • Sauer TJ, Norman JM, Tanner CB, Wilson TB (1995) Measurement of heat and vapor transfer coefficients at the soil surface beneath a maize canopy using source plates. Agric For Meteorol 75(1–3):161–189. https://doi.org/10.1016/0168-1923(94)02209-3

    Article  Google Scholar 

  • Semmens KA, Anderson MC, Kustas WP, Gao F, Alfieri JG, McKee L, Prueger JH, Hain CR, Cammalleri C, Yang Y, Xia T, Sanchez L, Alsina MM, Vélez M (2016) Monitoring daily evapotranspiration over two california vineyards using Landsat 8 in a multi-sensor data fusion approach. Remote Sens Environ 185:155–170. https://doi.org/10.1016/j.rse.2015.10.025

    Article  Google Scholar 

  • Song L, Liu S, Zhang X, Zhou J, Li M (2015) Estimating and validating soil evaporation and crop transpiration during the HiWATER-MUSOEXE. IEEE Geosci Remote Sens Lett 12(2):334–338. https://doi.org/10.1109/LGRS.2014.2339360

    Article  Google Scholar 

  • Song L, Liu S, Kustas WP, Zhou J, Xu Z, Xia T, Li M (2016) Application of remote sensing-based two-source energy balance model for mapping field surface fluxes with composite and component surface temperatures. Agric For Meteorol 230231:8–19. https://doi.org/10.1016/j.agrformet.2016.01.005

    Article  Google Scholar 

  • Torres-Rua A (2017) Vicarious calibration of sUAS microbolometer temperature imagery for estimation of radiometric land surface temperature. Sensors 17:1499. https://doi.org/10.3390/s17071499

    Article  Google Scholar 

  • Torres-Rua A, Nieto H, Parry C, Alarab M, McKee M (2018) Inter-comparison of thermal sensors from ground aircraft, and satellite. Irrig Sci (this issue)

  • Willmott CJ (1981) On the validation of models. Phys Geogr 2(2):184–194. https://doi.org/10.1080/02723646.1981.10642213

    Article  Google Scholar 

  • Xia T, Kustas WP, Anderson MC, Alfieri JG, Gao F, McKee L, Prueger JH, Geli HME, Neale CMU, Sanchez L, Alsina MM, Wang Z (2016) Mapping evapotranspiration with high-resolution aircraft imagery over vineyards using one- and two-source modeling schemes. Hydrol Earth Syst Sci 20(4):1523–1545. https://doi.org/10.5194/hess-20-1523-2016

    Article  Google Scholar 

  • Zhao K, Popescu S (2009) Lidar-based mapping of leaf area index and its use for validating GLOBCARBON satellite LAI product in a temperate forest of the southern USA. Remote Sens Environ 113(8):1628–1645. https://doi.org/10.1016/j.rse.2009.03.006

    Article  Google Scholar 

  • Zipper SC, Loheide SP II (2014) Using evapotranspiration to assess drought sensitivity on a subfield scale with HRMET, a high resolution surface energy balance model. Agric For Meteorol 197:91–102. https://doi.org/10.1016/j.agrformet.2014.06.009

    Article  Google Scholar 

Download references

Acknowledgements

Partial funding provided by E.&J. Gallo Winery made possible the acquisition and processing of the high-resolution manned aircraft and UAV imagery collected during GRAPEX IOPs. In addition, thanks are given to the Utah Water Research Laboratory for the use of the AggieAir UAV platform, support personnel and partial funding. In addition, we would like to thank the staff of Viticulture, Chemistry and Enology Division of E.&J. Gallo Winery for the collection and processing of field data during GRAPEX IOPs. Finally, this project would not have been possible without the cooperation of Mr. Ernie Dosio of Pacific Agri Lands Management, along with the Borden/ McMannis vineyard staff, for logistical support of GRAPEX field and research activities. USDA is an equal opportunity provider and employer.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Héctor Nieto.

Additional information

Communicated by N. Agam.

Part of this research was conducted thanks to the MC-COFUND Talentia Program.

Appendices

Appendix

A TSEB model

The basic equation of the energy balance at the surface can be expressed following Eq. 4:

$$\begin{aligned} R_{\mathrm{n}}\approx&H+\lambda E+G \end{aligned}$$
(4a)
$$\begin{aligned} R_{\mathrm{n,S}}\approx&H_{S}+\lambda E_{S}+G \end{aligned}$$
(4b)
$$\begin{aligned} R_{\mathrm{n,C}}\approx&H_{C}+\lambda E_{C} \end{aligned}$$
(4c)

with \(R_\mathrm{n}\) being the net radiation, H is the sensible heat flux, \(\lambda E\) is the latent heat flux or evapotranspiration, and G is the soil heat flux. “C” and “S” subscripts refer to canopy and soil layers, respectively. The symbol “\(\approx\)” appears, since there are additional components of the energy balance that are usually neglected, such as heat advection, storage of energy in the canopy layer, or energy for the fixation of CO\(_2\) (Hillel 1998)

The key in TSEB models is the partition of sensible heat flux into the canopy and soil layers, which depends on the soil and canopy temperatures (\(T_\mathrm{S}\) and \(T_\mathrm{C}\), respectively). If we assume that there is an interaction between the fluxes of canopy and soil, due to an expected heating of the in-canopy air by heat transport coming from the soil, the resistances network in TSEB can be considered to be in series. In that case, H can be estimated as in Eq. 5 (Norman et al. 1995, Eqs. A1–A3):

$$\begin{aligned} \begin{aligned} H=H_{C}+H_{S}& =\rho _{air} C_{p}\frac{T_\mathrm{AC}-T_{A}}{R_{a}}\\ & =\rho _{air} C_{p}\left[ \frac{T_\mathrm{C}-T_\mathrm{AC}}{R_{x}}+\frac{T_\mathrm{S}-T_\mathrm{AC}}{R_{s}}\right] , \end{aligned} \end{aligned}$$
(5)

where \(\rho _{air}\) is the density of the air (\(\hbox {kg} hbox{m}^{-3}\)), \(C_{p}\) is the heat capacity of the air at constant pressure (\({\hbox {J} \hbox {kg}^{-1} \hbox {K}^{-1}}\)), and \(T_\mathrm{AC}\) is the air temperature at the canopy interface, equivalent to the aerodynamic temperature \(T_0\), computed with Eq. 6 (Norman et al. 1995, Eq. 4):

$$\begin{aligned} T_\mathrm{AC}=\frac{\frac{T_A}{R_a}+\frac{T_\mathrm{C}}{R_x}+\frac{T_\mathrm{S}}{R_\mathrm{s}}}{\frac{1}{R_a} +\frac{1}{R_x}+\frac{1}{R_\mathrm{s}}}. \end{aligned}$$
(6)

Here, \(R_{a}\) is the aerodynamic resistance to heat transport (\({\hbox {s} \hbox {m}^{-1}}\)), \(R_{s}\) is the resistance to heat flow in the boundary layer immediately above the soil surface (\({\hbox {s} \hbox {m}^{-1}}\)), and \(R_{x}\) is the boundary layer resistance of the canopy of leaves (\({\hbox {s} \hbox {m}^{-1}}\)). The mathematical expressions of these resistances are detailed in Eq. 7 and in Norman et al. (1995) and Kustas and Norman (2000) and discussed in Kustas et al. (2016):

$$\begin{aligned} R_{a}&=\frac{\ln \left( \frac{z_{T}-d_0}{z_{0H}}\right) -\varPsi _{h}\left( \frac{z_{T}-d_0}{L}\right) +\varPsi _{h}\left( \frac{z_{0H}}{L}\right) }{\kappa '\,u_{*}}\end{aligned}$$
(7a)
$$\begin{aligned} R_{s}&=\frac{1}{c\left( T_\mathrm{S}-T_{A}\right) ^{1/3}+b\,u_{s}} \end{aligned}$$
(7b)
$$\begin{aligned} R_{x}&=\frac{C'}{\mathrm {LAI}}\left( \frac{l_w}{U_{d_0+z_{0M}}}\right) ^{1/2}, \end{aligned}$$
(7c)

where \(u_{*}\) is the friction velocity (\({\hbox {m} \hbox {s}^{-1}}\)) computed as

$$\begin{aligned} u_{*}=\frac{\kappa '\,u}{\left[ \ln \left( \frac{z_{u}-d_0}{z_{0M}}\right) -\varPsi _{m} \left( \frac{z_{u}-d_0}{L}\right) +\varPsi _m\left( \frac{z_{0M}}{L}\right) \right] }. \end{aligned}$$
(8)

In Eq. , \(z_{u}\) and \(z_{T}\) are the measurement heights for wind speed u (\({\hbox {m} \hbox {s}^{-1}}\)) and air temperature \(T_{A}\) (K), respectively. \(d_0\) is the zero-plane displacement height, and \(z_{0M}\) and \(z_{0H}\) are the roughness length for momentum and heat transport, respectively (all those magnitudes expressed in m), with \(z_{0H}=z_{0M}\exp \left( {-kB^{-1}}\right)\). In the series version of TSEB, \(z_{0H}\) is assumed equal to \(z_{0M}\), since the term \(R_x\) already accounts for the different efficiencies between heat and momentum transport (Norman et al. 1995), and therefore, \(kB^{-1}=0\). The value of \(\kappa '=0.4\) is the von Karman’s constant. The \(\varPsi _{m}\left( \zeta \right)\) terms in Eqs. 7a and are the adiabatic correction factors for momentum. The formulations of these two factors are described in Brutsaert (1999) and Brutsaert (2005). These corrections depend on the atmospheric stability, which is expressed using the Monin–Obukhov length L (m):

$$\begin{aligned} L=\frac{-u_{*}^{3}\rho _{air}}{k\,g\left[ ^H/_{\left( T_{A}C_{p}\right) } +0.61E\right] }, \end{aligned}$$
(9)

where H is the bulk sensible heat flux (\({\hbox {W} \hbox {m}^{-2}}\)), E is the rate of surface evaporation (\({\hbox {kg} \hbox {s}^{-1}}\)), and g is the acceleration of gravity (\({\hbox {m} \hbox {s}^{-2}}\))

The coefficients b , c in Eq. 7b depend on turbulent length scale in the canopy, soil-surface roughness, and turbulence intensity in the canopy, which are discussed in Sauer et al. (1995), Kondo and Ishida (1997) and Kustas et al. (2016). \(C'\) is assumed to be 90\(\hbox { s}^{^1/_2} \hbox { m}^{-1}\) and \(l_w\) is the average leaf width (m)

B Modifications to TSEB model for row crops

B.1 Radiation transmission in row crops

The clumping index for row crops is defined as the factor that modifies the leaf area index of a real canopy (F) in a fictitious homogeneous canopy with \(\mathrm{LAI}_{\mathrm{eff}}=\varOmega F\) such as its gap fraction is the same as the gap fraction of the real-world canopy (\(G\left( \theta ,\phi \right)\)):

$$\begin{aligned} \varOmega \left( \theta ,\phi \right) =\frac{-\log \left[ G\left( \theta ,\phi \right) \right] }{\kappa _{be}\left( \theta \right) F}, \end{aligned}$$
(10)

where \(\kappa _{be}\left( \theta \right)\) is the beam extinction coefficient through a plant with an ellipsoidal inclination distribution (Campbell 1986, 1990), \(\theta\) is the zenith incidence angle, and \(\psi\) is the relative azimuth angle between the incidence beam and the row direction

Our modelled real canopy consists of a horizontally infinite long prism with a total height \(h_c\) (i.e., the canopy height) and a width \(w_c\) (i.e., canopy width) that is placed above the ground at \(h_b\) (i.e., the height of the first living branch). This canopy contains finite-sized leaves randomly placed (no clumping within the canopy) oriented according to a ellipsoidal leaf angle distribution function (Campbell 1990) with a total leaf area index F (Fig. 3).

Then, the real canopy gap fraction consists of the sunlit part of the bare soil that is not shaded by the canopy plus the gaps caused by the solar beam passing through the crop canopy (Eq. 11):

$$\begin{aligned} G\left( \theta ,\phi \right) =f_{sc}\left( \theta ,\phi \right) \exp \left[ -\kappa _{be}\left( \theta \right) F\right] +\left[ 1-f_{sc}\left( \theta ,\phi \right) \right] . \end{aligned}$$
(11)

The solar canopy view factor \(f_{sc}\left( \theta ,\phi \right)\) is the fraction of soil that is cast by shadows (Colaizzi et al. 2012a) and in our case is estimated as

$$\begin{aligned} f_{sc}\left( \theta ,\phi \right) =\frac{w_{c}+\left( h_c-h_b\right) tan\theta |\sin \phi |}{L}, \end{aligned}$$
(12)

where L is the row separation (m). For a vertical projection (\(\theta =0\)), Eq. 12 reduces to \(w_{c}/L\), the fractional cover.

B.2 Massman et al. (2017) wind attenuation profile

Compared to previously used canopy wind profiles such as Goudriaan (1977) or Massman (1987), the additional key input required in Massman et al. (2017) wind attenuation model is the relative canopy foliage distribution, computed as in Eq. 13:

$$\begin{aligned} ha\left( \xi \right) =PAI\frac{f_a\left( \xi \right) }{\int _0^1f_a\left( \xi '\right) d\xi '}, \end{aligned}$$
(13)

where \(\tfrac{f_a\left( \xi \right) }{\int _0^1f_a\left( \xi '\right) d\xi '}\) is the relative canopy shape (i.e., \(\sum {\tfrac{f_a\left( \xi \right) }{\int _0^1f_a\left( \xi '\right) d\xi '}}=1\), and \(\xi =z/h_c\)) and PAI is the plant (leaves+stems) area index. Massman et al. (2017) modelled \(f_a\left( \xi \right)\) as a combination of asymmetric Gaussian curves, but \(f_a\left( \xi \right)\) can also be estimated as a continuous curve obtained from canopy structure measurements or three-dimensional cloud points, such as in Nieto et al. (2018).

The canopy wind speed profile is then the product of two terms: one logarithmic profile (\(U_b\)) that is dominant near the ground and a second a hyperbolic cosine profile (\(U_t\)) that dominates near the top of the canopy, where the canopy foliage distribution plays a major role. Ancillary input in \(U_t\) is the the drag coefficient of the individual foliage elements (\(C_d\)), which is usually considered equal to 0.2 (Massman et al. 2017; Goudriaan 1977). Massman et al. (2017) model has as well the ability to consider variations of the drag coefficient due to either wind sheltering between foliage elements, or vertical variations independently of wind blocking. This effect can usually be disregarded in most canopies (Massman et al. 2017), so was it in this study.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nieto, H., Kustas, W.P., Torres-Rúa, A. et al. Evaluation of TSEB turbulent fluxes using different methods for the retrieval of soil and canopy component temperatures from UAV thermal and multispectral imagery. Irrig Sci 37, 389–406 (2019). https://doi.org/10.1007/s00271-018-0585-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00271-018-0585-9

Navigation