Search

Accurate evapotranspiration (ET) estimation is crucial for sustainable water management in the diverse and water-scarce Mediterranean region. This study compares three prominent models (Simulator of Terrestrial Ecohydrological Processes and Systems (STEPS), Soil-Canopy-Observation of Photosynthesis and Energy fluxes (SCOPE), and Two-Source Energy Balance (TSEB)) with established global datasets (Moderate Resolution Imaging Spectroradiometer 8-day global terrestrial product (MOD16A2), Global Land Evaporation Amsterdam Model (GLEAM), and TerraClimate) at multiple spatial and temporal scales and validates model outcomes with eddy covariance based ground measurements. Insufficient ground-based observations limit comprehensive model validation in the eastern Mediterranean part (Turkey and Balkans). The results reveal significant discrepancies among models and datasets, highlighting the challenges of capturing ET variability in this complex region. Differences are attributed to variations in ecosystem type, energy balance calculations, and water availability constraints. Ground validation shows that STEPS performs well in some French and Italian forests and crops sites but struggles with seasonal ET patterns in some locations. SCOPE mostly overestimates ET due to detailed radiation flux calculations and lacks accurate water limitation representation. TSEB faces challenges in capturing ET variations across different ecosystems at a coarser 10 km resolution. No single model and global dataset accurately represent ET across the entire region. Model performance varies by region and ecosystem. As GLEAM and TSEB excel in semi-arid Savannahs, STEPS and SCOPE are better in grasslands, croplands, and forests in few locations (5 out of 18 sites) which indicates these models need calibration for other locations and ecosystem types. Thus, a region-specific model calibration and validation, sensitive to extremely humid and arid conditions can improve ET estimation across the diverse Mediterranean region.

Crop phenology data offer crucial information for crop yield estimation, agricultural management, and assessment of agroecosystems. Such information becomes more important in the context of increasing year-to-year climatic variability. The dataset provides in-situ crop phenology data (first leaves emergence and harvest date) of major European crops (wheat, corn, sunflower, rapeseed) from seventeen field study sites in Bulgaria and two in France. Additional information such as the sowing date, area of each site, coordinates, method and equipment used for phenophase data estimation, and photos of the France sites are also provided. The georeferenced ground-truth dataset provides a solid base for a better understanding of crop growth and can be used to validate the retrieval of phenological stages from remote sensing data.

Abstract Although soil moisture is a key factor of hydrologic and climate applications, global continuous high resolution soil moisture datasets are still limited. Here we use physics-informed machine learning to generate a global, long-term, spatially continuous high resolution dataset of surface soil moisture, using International Soil Moisture Network (ISMN), remote sensing and meteorological data, guided with the knowledge of physical processes impacting soil moisture dynamics. Global Surface Soil Moisture (GSSM1 km) provides surface soil moisture (0–5 cm) at 1 km spatial and daily temporal resolution over the period 2000–2020. The performance of the GSSM1 km dataset is evaluated with testing and validation datasets, and via inter-comparisons with existing soil moisture products. The root mean square error of GSSM1 km in testing set is 0.05 cm3/cm3, and correlation coefficient is 0.9. In terms of the feature importance, Antecedent Precipitation Evaporation Index (APEI) is the most important significant predictor among 18 predictors, followed by evaporation and longitude. GSSM1 km product can support the investigation of large-scale climate extremes and long-term trend analysis.

Current and upcoming Sun-Induced chlorophyll Fluorescence (SIF) satellite products (e.g., GOME, TROPOMI, OCO, FLEX) have medium-to-coarse spatial resolutions (i.e., 0.3–80 km) and integrate radiances from different sources into a single ground surface unit (i.e., pixel). However, intrapixel heterogeneity, i.e., different soil and vegetation fractional cover and/or different chlorophyll content or vegetation structure in a fluorescence pixel, increases the challenge in retrieving and quantifying SIF. High spatial resolution Sentinel-2 (S2) data (20 m) can be used to better characterize the intrapixel heterogeneity of SIF and potentially extend the application of satellite-derived SIF to heterogeneous areas. In the context of the COST Action Optical synergies for spatiotemporal SENsing of Scalable ECOphysiological traits (SENSECO), in which this study was conducted, we proposed direct (i.e., spatial heterogeneity coefficient, standard deviation, normalized entropy, ensemble decision trees) and patch mosaic (i.e., local Moran’s I) approaches to characterize the spatial heterogeneity of SIF collected at 760 and 687 nm (SIF760 and SIF687, respectively) and to correlate it with the spatial heterogeneity of selected S2 derivatives. We used HyPlant airborne imagery acquired over an agricultural area in Braccagni (Italy) to emulate S2-like top-of-the-canopy reflectance and SIF imagery at different spatial resolutions (i.e., 300, 20, and 5 m). The ensemble decision trees method characterized FLEX intrapixel heterogeneity best (R2 > 0.9 for all predictors with respect to SIF760 and SIF687). Nevertheless, the standard deviation and spatial heterogeneity coefficient using k-means clustering scene classification also provided acceptable results. In particular, the near-infrared reflectance of terrestrial vegetation (NIRv) index accounted for most of the spatial heterogeneity of SIF760 in all applied methods (R2 = 0.76 with the standard deviation method; R2 = 0.63 with the spatial heterogeneity coefficient method using a scene classification map with 15 classes). The models developed for SIF687 did not perform as well as those for SIF760, possibly due to the uncertainties in fluorescence retrieval at 687 nm and the low signal-to-noise ratio in the red spectral region. Our study shows the potential of the proposed methods to be implemented as part of the FLEX ground segment processing chain to quantify the intrapixel heterogeneity of a FLEX pixel and/or as a quality flag to determine the reliability of the retrieved fluorescence.

Most applications of remote sensing in agricultural crop monitoring use multispectral imaging techniques, but with upcoming hyperspectral missions, the opportunity arises to better estimate pigment absorption and crop structure by exploiting the full solar reflective spectrum. In this study, we demonstrate how hyperspectral time series can be used with the Soil Canopy Observation of Photochemistry and Energy fluxes (SCOPE) model to estimate crop yield variability among fields, crop varieties and nitrogen treatments generically, i.e. without a calibration with in situ, data. Field experiments were conducted in two potato fields in the Netherlands between May and September 2019. The fields were planted with five varieties of potato, under three nitrogen fertilization treatments. By fitting the model to the full VNIR-SWIR spectrum of measured hyperspectral reflectance, we retrieved the model input parameters of Leaf Area Index (LAI), leaf chlorophyll content (Cab) and leaf water content (Cw) and simulated the photosynthesis throughout the season using data of local Automatic Weather Stations (AWS). Statistical analysis of measured and retrieved traits of LAI, Cab and canopy water content showed that two fields responded differently to the treatments, exhibiting fewer classes than were expected based on the experimental design. Potato yield, which was estimated as the sum of photosynthesis flux multiplied by the harvest index of 0.64, correlated with the measured tuber dry weight with R2 0.36 and RMSE 2.5 t ha−1. This study demonstrates that even in the absence of crop or variety specific information, hyperspectral reflectance and local weather data ingested into SCOPE can explain a substantial part of the observed variability in yield among fields.

The study of peatland is challenging due to the water saturation and evergreen mixed vegetation that ranges from simple forms of plants such as mosses to higher forms of plants such as cranberries, grasses, etc. The changing water level through the growing season makes the peatland vegetation very dynamic. In this work, we have used ground-level remote-sensing signals to understand the dynamic nature of peatland vegetation. We have also estimated the leaf area index (LAI) and Sun-Induced fluorescence (SIF) through the Soil Canopy Observation of Photosynthesis and Energy fluxes (SCOPE) model. The estimated LAI and SIF were compared with the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Near-Infrared Reflectance of vegetation (NIRv), and measured SIF. The modeled LAI was observed to be significantly correlated with NDVI, EVI, and NIRv, whereas a good correlation was observed between measured and modeled SIF. Along with showing the dynamic behavior of peatland vegetation, the study indicates that SCOPE in its inverted form can be used to estimate reflectance-based LAI for peatland, which can be more reliable to present biomass and productivity of peatland ecosystem in comparison to transmittance-based LAI measurement for such ecosystem. The good correlation between measured and modeled SIF at 760 nm indicates that a reliable SIF value can be estimated through the SCOPE model for a complex ecosystem such as peatland, which can be very helpful in the absence of high-resolution hyperspectral data (usually used for SIF measurements).

In this study, we demonstrate that the Google Earth Engine (GEE) dataset of Sentinel-3 Ocean and Land Color Instrument (OLCI) level-1 deviates from the original Copernicus Open Access Data Hub Service (DHUS) data by 10–20 W m−2 sr−1μμm−1 per pixel per band. We compared GEE and DHUS single pixel time series for the period from April 2016 to September 2020 and identified two sources of this discrepancy: the ground pixel position and reprojection. The ground pixel position of OLCI product can be determined in two ways: from geo-coordinates (DHUS) or from tie-point coordinates (GEE). We recommend using geo-coordinates for pixel extraction from the original data. When the Sentinel Application Platform (SNAP) Pixel Extraction Tool is used, an additional distance check has to be conducted to exclude pixels that lay further than 212 m from the point of interest. Even geo-coordinates-based pixel extraction requires the homogeneity of the target area at a 700 m diameter (49 ha) footprint (double of the pixel resolution). The GEE OLCI dataset can be safely used if the homogeneity assumption holds at 2700 m diameter (9-by-9 OLCI pixels) or if the uncertainty in the radiance of 10% is not critical for the application. Further analysis showed that the scaling factors reported in the GEE dataset description must not be used. Finally, observation geometry and meteorological data are not present in the GEE OLCI dataset, but they are crucial for most applications. Therefore, we propose to calculate angles and extraterrestrial solar fluxes and to use an alternative data source—the Copernicus Atmosphere Monitoring Service (CAMS) dataset—for meteodata.

Sentinel-3 satellite has provided simultaneous observations in the optical (visible, near infrared (NIR), shortwave infrared (SWIR)) and thermal infrared (TIR) domains since 2016, with a revisit time of 1−2 days. The high temporal resolution and spectral coverage make the data of this mission attractive for vegetation monitoring. This study explores the possibilities of using the Soil Canopy Observation, Photochemistry and Energy fluxes (SCOPE) model together with Sentinel-3 to exploit the two sensors onboard of Sentinel-3 (the ocean and land color instrument (OLCI) and sea and land surface temperature radiometer (SLSTR)) in synergy. Sobol’ variance based global sensitivity analysis (GSA) of top of atmosphere (TOA) radiance produced with a coupled SCOPE-6S model was conducted for optical bands of OLCI and SLSTR, while another GSA of SCOPE was conducted for the land surface temperature (LST) product of SLSTR. The results show that in addition to ESA level-2 Sentinel-3 products, SCOPE is able to retrieve leaf area index (LAI), leaf chlorophyll content (Cab), leaf water content (Cw), leaf senescent material (Cs), leaf inclination distribution (LAD). Leaf dry matter content (Cdm) and soil brightness, despite being important, were not confidently retrieved in some cases. GSA of LST in TIR domain showed that plant biochemical parameters—maximum carboxylation rate (Vcmax) and stomata conductance-photosynthesis slope (Ball-Berry m)—can be constrained if prior information on near-surface weather conditions is available. We conclude that the combination of optical and thermal domains facilitates the constraint of the land surface energy balance using SCOPE.