A Surface Temperature Initiated Closure (STIC) for surface energy balance fluxes

The use of Penman-Monteith (PM) equation in thermal remote sensing based surface energy balance modeling is not prevalent due to the unavailability of any direct method to integrate thermal data into the PM equation and due to the lack of physical models expressing the surface (or stomatal) and boundary layer conductances (g(S) and g(B)) as a function of surface temperature. Here we demonstrate a new method that physically integrates the radiometric surface temperature (T-S) into the PM equation for estimating the terrestrial surface energy balance fluxes (sensible heat, H and latent heat, lambda E). The method combines satellite T-S data with standard energy balance closure models in order to derive a hybrid closure that does not require the specification of surface to atmosphere conductance terms. We call this the Surface Temperature Initiated Closure (STIC), which is formed by the simultaneous solution of four state equations. Taking advantage of the psychrometric relationship between temperature and vapor pressure, the present method also estimates the near surface moisture availability (M) from T-S, air temperature (T-A) and relative humidity (R-H), thereby being capable of decomposing lambda E into evaporation (lambda E-E) and transpiration (lambda E-T). STIC is driven with T-S, T-A, R-H, net radiation (R-N), and ground heat flux (G). T-S measurements from both MODIS Terra (MOD11A2) and Aqua (MYD11A2) were used in conjunction with FLUXNET R-N, G, T-A, R-H, lambda E and H measurements corresponding to the MODIS equatorial crossing time. The performance of STIC has been evaluated in comparison to the eddy covariance measurements of lambda E and H at 30 sites that cover a broad range of biomes and climates. We found a RMSE of 37.79 (11%) (with MODIS Terra T-S) and 44.27 W m(-2) (15%) (with MODIS Aqua T-S) in lambda E estimates, while the RMSE was 37.74(9%) (with Terra) and 44.72 W m(-2) (8%) (with Aqua) in H. STIC could efficiently capture the lambda E dynamics during the dry down period in the semi-arid landscapes where lambda E is strongly governed by the subsurface soil moisture and where the majority of other lambda E models generally show poor results. Sensitivity analysis revealed a high sensitivity of both the fluxes to the uncertainties in T-S. A realistic response and modest relationship was also found when partitioned lambda E components (lambda E-E and lambda E-T) were compared to the observed soil moisture and rainfall. This is the first study to report the physical integration of T-S into the PM equation and finding analytical solution of the physical (g(B)) and physiological conductances (g(S)). The performance of STIC over diverse biomes and climates points to its potential to benefit future NASA and NOAA missions having thermal sensors, such as HyspIRI, GeoSTAR and GOES-R for mapping multi-scale lambda E and drought. (C) 2013 Elsevier Inc. All rights reserved.