Quality checked and gap-filled daily MODIS observations of surface reflectance and land surface temperature at global eddy co-variance sites for the time period 2000-2020. Two product versions: one features all MODIS pixels within 2km radius around a given site, and a second version consists of an average time series that represents the area within 1km2 around a site. All data layers have a complementary layer with gap-fill information. MODIS data comprise all sites in the Fluxnet La Thuile, Berkeley and ICOS Drought 2018 data releases. FluxnetEO MODIS reflectance products: enhanced vegetation index (EVI), normalized difference vegetation index (NDVI), generalized NDVI (kNDVI), near infra-red reflectance of vegetation (NIRv), normalized difference water index (NDWI) with band 5, 6, or 7 as reference, the scaled wide dynamic range vegetation index (sWDRVI), surface reflectance in MODIS bands 1-7. Based on the NASA MCD43A4 and MCD43A2 collection 6 products with a pixel size of 500m. FluxnetEO MODIS land surface temperature: Terra and Aqua, day and night, at native viewing zenith angle as well as corrected to viewing zenith angles of 0 and 40degrees (Ermida et al., 2018, RSE). Based on NASA MOD11A1 and MYD11A1 collection 6 at a pixel size of 1km. Supplementary data to Walther, S., Besnard, S., Nelson, J.A., El-Madany, T. S., Migliavacca, M., Weber, U., Ermida, S. L., Brümmer, C., Schrader, F., Prokushkin, A., Panov, A., Jung, M. , 2021. A view from space on global flux towers by MODIS and Landsat: The FluxnetEO dataset, in preparation for Biogeosciences Discussions. The data are separated in zip-files by continents and groups of countries.