“…where S in is the incoming short wave radiation (W m 2 ), L in is the incoming long wave radiation (W m 2 ), L out is the outgoing long wave radiation (W m 2 ), r is the surface albedo, which can be retrieved from Band 1, 2, 3, 4, 5, and 7 of Landsat5 TM (Chen and Ohring, 1984;Koepke et al, 1985;Bastiaanssen et al, 1998a;Bastiaanssen, 2000;Wang et al, 2000), ε is the surface emissivity, which can be estimated by the Normal Difference Vegetation Index (NDVI) retrieved from the red and near infrared bands using semi-empirical relationships ( Van de Griend and Owe, 1993;Bastiaanssen et al, 1998a), T sur is the surface temperature (K), which can be retrieved by the thermal infrared band, Band 4853 6 of Landsat5 TM (Bastiaanssen, 2000;Chander and Markham, 2003). The atmospheric emissivity ε a can be estimated from atmospheric temperature and humidity (Brutsaert, 1975), is the Stefan-Boltzman constant (5Ð67 ð 10 8 W m 2 K 4 ).…”