In the last few decades, researchers have developed a plethora of hyperspectral Earth Observation (EO) remote sensing techniques, analysis and applications. While hyperspectral exploratory sensors are demonstrating their potential, Sentinel-2 multispectral satellite remote sensing is now providing free, open, global and systematic high resolution visible and infrared imagery at a short revisit time. Its recent launch suggests potential synergies between multi- and hyper-spectral data. This study, therefore, reviews 20 years of research and applications in satellite hyperspectral remote sensing through the analysis of Earth observation hyperspectral sensors’ publications that cover the Sentinel-2 spectrum range: Hyperion, TianGong-1, PRISMA, HISUI, EnMAP, Shalom, HyspIRI and HypXIM. More specifically, this study (i) brings face to face past and future hyperspectral sensors’ applications with Sentinel-2’s and (ii) analyzes the applications’ requirements in terms of spatial and temporal resolutions. Eight main application topics were analyzed including vegetation, agriculture, soil, geology, urban, land use, water resources and disaster. Medium spatial resolution, long revisit time and low signal-to-noise ratio in the short-wave infrared of some hyperspectral sensors were highlighted as major limitations for some applications compared to the Sentinel-2 system. However, these constraints mainly concerned past hyperspectral sensors, while they will probably be overcome by forthcoming instruments. Therefore, this study is putting forward the compatibility of hyperspectral sensors and Sentinel-2 systems for resolution enhancement techniques in order to increase the panel of hyperspectral uses.

The frequent acquisitions of fine spatial resolution imagery (10 m) offered by recent multispectral satellite missions, including Sentinel-2, can resolve single agricultural fields and thus provide crop-specific phenology metrics, a crucial information for crop monitoring. However, effective phenology retrieval may still be hampered by significant cloud cover. Synthetic aperture radar (SAR) observations are not restricted by weather conditions, and Sentinel-1 thus ensures more frequent observations of the land surface. However, these data have not been systematically exploited for phenology retrieval so far. In this study, we extracted crop-specific land surface phenology (LSP) from Sentinel-1 and Sentinel-2 of major European crops (common and durum wheat, barley, maize, oats, rape and turnip rape, sugar beet, sunflower, and dry pulses) using ground-truth information from the “Copernicus module” of the Land Use/Cover Area frame statistical Survey (LUCAS) of 2018. We consistently used a single model-fit approach to retrieve LSP metrics on temporal profiles of CR (Cross Ratio, the ratio of the backscattering coefficient VH/VV from Sentinel-1) and NDVI (Normalized Difference Vegetation Index from Sentinel-2). Our analysis revealed that LSP retrievals from Sentinel-1 are comparable to those of Sentinel-2, particularly for winter crops. The start of season (SOS) timings, as derived from Sentinel-1 and -2, are significantly correlated (average r of 0.78 for winter and 0.46 for summer crops). The correlation is lower for end of season retrievals (EOS, r of 0.62 and 0.34). Agreement between LSP derived from Sentinel-1 and -2 varies among crop types, ranging from r = 0.89 and mean absolute error MAE = 10 days (SOS of dry pulses) to r = 0.15 and MAE = 53 days (EOS of sugar beet). Observed deviations revealed that Sentinel-1 and -2 LSP retrievals can be complementary; for example for winter crops we found that SAR detected the start of the spring growth while multispectral data is sensitive to the vegetative growth before and during winter. To test if our results correspond reasonably to in-situ data, we compared average crop-specific LSP for Germany to average phenology from ground phenological observations of 2018 gathered from the German Meteorological Service (DWD). Our study demonstrated that both Sentinel-1 and -2 can provide relevant and at times complementary LSP information at field- and crop-level.

Land cover and land use maps derived from satellite remote sensing imagery are critical to support biodiversity and conservation, especially over large areas. With its 10 m to 20 m spatial resolution, Sentinel-2 is a promising sensor for the detection of a variety of landscape features of ecological relevance. However, many components of the ecological network are still smaller than the 10 m pixel, i.e., they are sub-pixel targets that stretch the sensor’s resolution to its limit. This paper proposes a framework to empirically estimate the minimum object size for an accurate detection of a set of structuring landscape foreground/background pairs. The developed method combines a spectral separability analysis and an empirical point spread function estimation for Sentinel-2. The same approach was also applied to Landsat-8 and SPOT-5 (Take 5), which can be considered as similar in terms of spectral definition and spatial resolution, respectively. Results show that Sentinel-2 performs consistently on both aspects. A large number of indices have been tested along with the individual spectral bands and target discrimination was possible in all but one case. Overall, results for Sentinel-2 highlight the critical importance of a good compromise between the spatial and spectral resolution. For instance, the Sentinel-2 roads detection limit was of 3 m and small water bodies are separable with a diameter larger than 11 m. In addition, the analysis of spectral mixtures draws attention to the uneven sensitivity of a variety of spectral indices. The proposed framework could be implemented to assess the fitness for purpose of future sensors within a large range of applications.

Accurately characterizing land surface changes with Earth Observation requires geo-located ground truth. In the European Union (EU), a tri-annual surveyed sample of land cover and land use has been collected since 2006 under the Land Use/Cover Area frame Survey (LUCAS). A total of 1351293 observations at 651780 unique locations for 106 variables along with 5.4 million photos were collected during five LUCAS surveys. Until now, these data have never been harmonised into one database, limiting full exploitation of the information. This paper describes the LUCAS point sampling/surveying methodology, including collection of standard variables such as land cover, environmental parameters, and full resolution landscape and point photos, and then describes the harmonisation process. The resulting harmonised database is the most comprehensive in-situ dataset on land cover and use in the EU. The database is valuable for geo-spatial and statistical analysis of land use and land cover change. Furthermore, its potential to provide multi-temporal in-situ data will be enhanced by recent computational advances such as deep learning.

A novel methodology is proposed to robustly map oil seed rape (OSR) flowering phenology from time series generated from the Copernicus Sentinel-1 (S1) and Sentinel-2 (S2) sensors. The time series are averaged at parcel level, initially for a set of 229 reference parcels for which multiple phenological observations on OSR flowering have been collected from April 21 to May 19, 2018. The set of OSR parcels is extended to a regional sample of 32,355 OSR parcels derived from a regional S2 classification. The study area comprises the northern Brandenburg and Mecklenburg-Vorpommern (N) and the southern Bavaria (S) regions in Germany. A method was developed to automatically compute peak flowering at parcel level from the S2 time signature of the Normalized Difference Yellow Index (NDYI) and from the local minimum in S1 VV polarized backscattering coefficients. Peak flowering was determined at a temporal accuracy of 1 to 4 days. A systematic flowering delay of 1 day was observed in the S1 detection compared to S2. Peak flowering differed by 12 days between the N and S. Considerable local variation was observed in the N-S parcel-level flowering gradient. Additional in-situ phenology observations at 70 Deutscher Wetterdienst (DWD) stations confirm the spatial and temporal consistency between S1 and S2 signatures and flowering phenology across both regions. Conditions during flowering strongly determine OSR yield, therefore, the capacity to continuously characterize spatially the timing of key flowering dates across large areas is key. To illustrate this, expected flowering dates were simulated assuming a single OSR variety with a 425 growing degree days (GDD) requirement to reach flowering. This GDD requirement was calculated based on parcel-level peak flowering dates and temperatures accumulated from 25-km gridded meteorological data. The correlation between simulated and S2 observed peak flowering dates still equaled 0.84 and 0.54 for the N and S respectively. These Sentinel-based parcel-level flowering parameters can be combined with weather data to support in-season predictions of OSR yield, area, and production. Our approach identified the unique temporal signatures of S1 and S2 associated with OSR flowering and can now be applied to monitor OSR phenology for parcels across the globe.

Accurate maps of surface water extent are of paramount importance for water management, satellite data processing and climate modeling. Several maps of water bodies based on remote sensing data have been released during the last decade. Nonetheless, none has a truly (90 ${}^{\circ }$ N/90 ${}^{\circ }$ S) global coverage while being thoroughly validated. This paper describes a global, spatially-complete (void-free) and accurate mask of inland/ocean water for the 2000–2012 period, built in the framework of the European Space Agency (ESA) Climate Change Initiative (CCI). This map results from the synergistic combination of multiple individual SAR and optical water body and auxiliary datasets. A key aspect of this work is the original and rigorous stratified random sampling designed for the quality assessment of binary classifications where one class is marginally distributed. Input and consolidated products were assessed qualitatively and quantitatively against a reference validation database of 2110 samples spread throughout the globe. Using all samples, overall accuracy was always very high among all products, between $98\%$ and $100\%$ . The CCI global map of open water bodies provided the best water class representation (F-score of $89\%$ ) compared to its constitutive inputs. When focusing on the challenging areas for water bodies’ mapping, such as shorelines, lakes and river banks, all products yielded substantially lower accuracy figures with overall accuracies ranging between $74\%$ and $89\%$ . The inland water area of the CCI global map of open water bodies was estimated to be 3.17 million km ${}^2$ ± 0.24 million km ${}^2$ . The dataset is freely available through the ESA CCI Land Cover viewer.