*2.2. Satellite Images*

̊ Terra/MODIS images with near nadir observations of the study site (field of view <25◦ ) were selected to minimize the bowtie effect [37]. Six different dates were used for this work (Table 1). MODIS VNIR and TIR data were extracted from the MOD09GQ, and MOD11\_L2 products, respectively, downloaded from the NASA Earthdata Search tool. MOD09GQ offers surface Red and NIR reflectivity values at a 250 m spatial resolution. MOD11\_L2 product provides LST, atmospherically corrected with a split-window algorithm, at a 1000 m spatial resolution [38].

Sentinel-2A and Sentinel-2B images concurrent or within ±1-day timing difference with MODIS were used (see dates in Table 1). S2 Level-2A products were downloaded from the Copernicus Open Access Hub, and they contain 10-m surface reflectance values. Bands 4 and 8 were used to compose the Normalized Difference Vegetation Index (NDVI). Figure 2 shows an example of the spatial distribution of the NDVI over the study site. Note the wide range in NDVI values available in the area during the experiment. Plot in Figure 3 shows the NDVI values for the different crop fields in the selected dates.



**Figure 2.** Sentinel-2 Normalized Difference Vegetation Index (NDVI) image of "Las Tiesas" farm, corresponding to date 25 July 2018.

**Figure 3.** Distribution of the Sentinel-2 NDVI values for the different sites labeled in Figure 1, and the images/dates listed in Table 1.

Landsat-7/ETM+ overpasses were also available for 3 of the selected dates (9, 16 and 25 July). Although ETM+ TIR band has a resolution of 60 m, a cubic convolution resampling to 30 m is applied for user distribution. Thus, these images, with a 30-m spatial resolution, were used as a reference for an extended validation of the disaggregated LST. For the VNIR bands the Landsat Surface Reflectance (CDR) product was used, whereas the original TIR data in band 6 were corrected from atmospheric and emissivity effects following the method proposed by Galve et al. [39].
