After Landsat 8 was launched in 2013, it was observed that for Thermal Infrared sensor (TIRS) bands, radiance from outside of an instrument’s field-of-view produced a non-uniform ghost signal across the focal plane that varied depending on the out-of-scene content (i.e., the stray light effect). A new stray light correction algorithm (SLCA) is currently operational and has been implemented into the United States Geological Survey (USGS) ground system since February 2017. The SLCA has also been applied to reprocess historical Landsat 8 scenes. After approximately two years of SLCA implementation, more land surface temperature (LST) validation studies are required to check the effect of correction in the estimation of LST from different retrieval algorithms. For this purpose, three different LST estimation method algorithms (i.e., the radiative transfer equation (RTE), single-channel algorithm (SCA), and split-window algorithm (SWA)) have been assessed. The study site is located on the campus of the University of Balearic Islands on the island of Mallorca (Spain) in the western Mediterranean Sea. The site is considered a heterogeneous area that is composed of different types of surfaces, such as buildings, asphalt roads, farming areas, sloped terrains, orange fields, almond trees, lawns, and some natural vegetation regions. Data from 21 scenes, which were acquired by the Landsat 8-TIRS sensor and extracted from a 100 × 100 m² pixel, were used to retrieve the LST with different algorithms; then, they were compared with in situ LST measurements from a broadband thermal infrared radiometer located on the same Landsat 8 pixel. The results show good performances of the three methods, with the SWA showing the lowest observed RMSE (within 1.6–2 K), whereas the SCA applied to the TIRS band 10 (10 µm) was also appropriate, with a RMSE ranging within 2.0–2.3 K. The LST estimates using the RTE algorithm display the highest observed RMSE values (within 2.0–3.6 K) of all of the compared methods, but with an almost unbiased value of −0.1 K for the case of techniques applied to band 10 data. The SWAs are the preferred method to estimate the LST in our study area. However, further validation studies around the world are required.

With a wide range of satellite-derived vegetation bio-geophysical products now available to users, validation efforts are required to assess their accuracy and fitness for purpose. Substantial progress in the validation of such products has been made over the last two decades, but quantification of the uncertainties associated with in situ reference measurements is rarely performed, and the incorporation of uncertainties within upscaling procedures is cursory at best. Since current validation practices assume that reference data represent the truth, our ability to reliably demonstrate compliance with product uncertainty requirements through conformity testing is limited. The Fiducial Reference Measurements for Vegetation (FRM4VEG) project, initiated by the European Space Agency, is aiming to address this challenge by applying metrological principles to vegetation and surface reflectance product validation. Following FRM principles, and in accordance with the International Standards Organisation’s (ISO) Guide to the Expression of Uncertainty in Measurement (GUM), for the first time, we describe an end-to-end uncertainty evaluation framework for reference data of two key vegetation bio-geophysical variables: the fraction of absorbed photosynthetically active radiation (FAPAR) and canopy chlorophyll content (CCC). The process involves quantifying the uncertainties associated with individual in situ reference measurements and incorporating these uncertainties within the upscaling procedure (as well as those associated with the high-spatial-resolution imagery used for upscaling). The framework was demonstrated in two field campaigns covering agricultural crops (Las Tiesas–Barrax, Spain) and deciduous broadleaf forest (Wytham Woods, UK). Providing high-spatial-resolution reference maps with per-pixel uncertainty estimates, the framework is applicable to a range of other bio-geophysical variables including leaf area index (LAI), the fraction of vegetation cover (FCOVER), and canopy water content (CWC). The proposed procedures will facilitate conformity testing of moderate spatial resolution vegetation bio-geophysical products in future validation exercises.

Almost fifty years have passed since the idea to retrieve a value for Evapotranspiration (ET) using remote sensing techniques was first considered. Numerous ET models have been proposed, validated and improved along these five decades, as the satellites and sensors onboard were enhanced. This study reviews most of the efforts in the progress towards providing a trustworthy value of ET by means of thermal remote sensing data. It starts with an in-depth reflection of the surface energy balance concept and of each of its terms, followed by the description of the approaches taken by remote sensing models to estimate ET from it in the last thirty years. This work also includes a chronological review of the modifications suggested by several researchers, as well as representative validations studies of such ET models. Present limitations of ET estimated with remote sensors onboard orbiting satellites, as well as at surface level, are raised. Current trends to face such limitations and a future perspective of the discipline are also exposed, for the reader’s inspiration.

The Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) and emissivity separation (MODTES) algorithm is the basis of the MOD21 product, which provides 1-km LSTs and emissivities for bands 29 (8.55 μm), 31 (11 μm), and 32 (12 μm). The MODTES algorithm uses the TES method with the water vapor scaling (WVS) method for refined atmospheric correction. The performance of the MODTES algorithm was tested with a set of MODIS data concurrent with ground LST and emissivity measurements. The test site is a large area of homogeneous full-cover rice crops (graybody), with high atmospheric water vapor. The data included LSTs measured along transects with multiple calibrated radiometers and emissivity measurements in different bands within 8-13 μm. We applied the full MODTES algorithm and that without the WVS method for comparison. For the data used here, MODTES minus in situ LSTs yielded a median difference of 0.5 K and a robust standard deviation (RSD) of 0.4 K. When the WVS method was not applied, we obtained a median bias of 1.3 K and an RSD of 0.9 K. MODTES retrieved emissivities were (median ± RSD) 0.966 ± 0.011 for band 29, 0.976 ± 0.004 for band 31, and 0.981 ± 0.003 for band 32 (0.961 ± 0.029, 0.971 ± 0.013, and 0.978 ± 0.007 without WVS), underestimating the ground values by 0.019, 0.009, and 0.001, respectively. Emissivities from the full MODTES algorithm presented less dispersion and were closer to the ground data, particularly for lower water vapor cases. These results show that MODTES performs within the predicted algorithm uncertainty and the importance of the WVS method for reducing residual errors in atmospheric correction.