Landsat and Sentinel-2 data archives provide ever-increasing amounts of satellite data for studying land cover and land use change (LCLUC) over the past four decades. However, the availability of cloud-, shadow-, and snow-free observations varies spatially and temporally due to climate and satellite data acquisition schemes. Spatio-temporal heterogeneity poses a major issue for some time-series analysis approaches, but can be addressed with pixel-based compositing that generates temporally equidistant cloud-free or near-cloud free synthetic images. Although much consideration is given to methods identifying the ‘best’ pixel value for each composite, determining the aggregation period receives less attention and is often done arbitrary, or based on expert intuition. Here, we evaluated data compositing windows ranging from five days to one year for 1984-2021 Landsat and 2015-2021 Sentinel‑2 time series across Europe. We considered separate and joint use of both data archives and analyzed spatio-temporal availability of composites during each calendar year and pixel-specific growing season. We reported mean annual composites’ availability investigating differences among biogeographical regions, checked feasibility of pan‑European analyses for three LCLUC applications based on annual, monthly and 10-day composites, and analyzed the shortest feasible compositing window ensuring ≥50% temporal data availability and interpolation of the remaining composites for individual years and across a variety of medium- and long‑term time windows. Our results highlighted low data coverage in the 1980s, 1990s, and in 2012, as well as spatial variability in data availability driven by climate and orbit overlaps, which altogether impact spatio-temporal consistency of medium- and long-term time series, limiting feasibility of some LCLUC analyses. We demonstrated that prior to 2011 monthly composites ensured overall 50-62% data coverage for each calendar year, and ~75% afterwards, with further increase to ~82% when Landsat and Sentinel-2 were combined. Temporal consistency of monthly composites was overall low and temporal interpolation augmenting up to 50% missing data each year and across a time window of interest, ensured feasibility of analyses. Applications based on shorter than monthly composites were challenging without joining Landsat and Sentinel‑2 archives after 2015, and beyond the Mediterranean biogeographical region. Using pixel-specific growing season data typically boosted data availability in most geographies and diminished most of the latitudinal differences, but feasibility of complete time series with sub-monthly compositing windows was still restricted to the most recent years, and required data interpolation. Overall, our analyses provided a detailed assessment of Landsat and Sentinel-2 data availability over Europe, and based on selected application examples, highlighted often lacking spatio-temporal consistency of time series with sub-monthly compositing windows and long-time periods, which might hinder feasibility of some LCLUC applications.

Landsat and Sentinel-2 acquisitions are among the most widely used medium-resolution optical data adopted for terrestrial vegetation applications, such as land cover and land use mapping, vegetation condition and phenology monitoring, and disturbance and change mapping. When combined, both data archives provide over 40 years, and counting, of continuous and consistent observations. Although the spatio-temporal availability of both data archives is well-known at the scene level, information on the actual availability of cloud-, snow-, and shade-free observations at the pixel level is lacking and should be explored individually for each study to correctly parametrize subsequent analyses. However, data exploration is time and resource consuming, thus is rarely performed a-priori. Consequently, the spatio-temporal heterogeneity of usable data is often inadequately accounted for in the analysis design, risking ill-advised selection of algorithms and hypotheses, and thus inferior quality of final results. Here we present precomputed data on the daily 1982‑2023 availability of usable Landsat and Sentinel-2 acquisitions across the globe, sampled at a pixel-level in a regular 0.18°-point grid. The dataset comprises separate Landsat- and Sentinel‑2‑specific data records, and is accompanied by a pixel-specific growing season information to facilitate data exploration. The dataset was derived based on freely available 1982-2023 Landsat surface reflectance (Collection 2) and Sentinel-2 top-of-the-atmosphere reflectance (pre‑Collection-1 and Collection-1) scenes from 2015 through 2023, following the methodology developed in the recent study on data availability over Europe [1]. Growing season information was derived based on 2001‑2019 time series of the yearly 500 m MODIS land cover dynamics product (MCD12Q2; Collection 6) [1]. As such, the dataset presents a unique overview of the spatio-temporal availability of usable daily Landsat and Sentinel-2 data at the global scale, hence offering much needed a-priori information aiding identification of appropriate methods and challenges for terrestrial vegetation analyses at local to global scale.