Search

The Biomass Proxy is a new cloud-free vegetation monitoring product that offers timely and analysis-ready data indicative of above-ground crop biomass dynamics at 10m spatial resolution. The Biomass Proxy links the consistent and continuous temporal signal of the Sentinel-1 Cross Ratio (CR), a vegetation index derived from Synthetic Aperture Radar backscatter, with the spatial information of the Sentinel-2 Normalized Difference Vegetation Index (NDVI), a vegetation index derived from optical observations. A global scaling relationship between CR and NDVI forms the basis of a novel fusion methodology based on static and dynamic combinations of temporal and spatial responses of CR and NDVI at field level. The fusion process is used to mitigate the impact on product quality of low satellite revisit periods due to acquisition design or persistent cloud coverage, and to respond to rapid changes in a timely manner to detect environmental and management events. The resulting Biomass Proxy provides time series that are continuous, unhindered by clouds, and produced uniformly across all geographical regions and crops. The Biomass Proxy offers opportunities including improved crop growth monitoring, event detection, and phenology stage detection. [ABSTRACT FROM AUTHOR]

The monitoring of Natura 2000 habitats (Habitat Directive 92/43/EEC) is a key activity ensuring the sufficient protection of European biodiversity. Reporting on the status of Natura 2000 habitats is required every 6 years. Although field mapping is still an indispensable source of data on the status of Natura 2000 habitats, and very good field-based data exist in some countries, keeping the field-based habitat maps up to date can be an issue. Remote sensing techniques represent an excellent alternative. Here, we present a new method for detecting habitats that were likely misclassified during the field mapping or that have changed since then. The method identifies the possible habitat mapping errors as the so-called "attribute outliers", i.e., outlying observations in the feature space of all relevant (spectral and other) characteristics of an individual habitat patch. We used the Czech Natura 2000 Habitat Layer as field-based habitat data. To prepare the feature space of habitat characteristics, we used a fusion of Sentinel-1 and Sentinel-2 satellite data along with a Digital Elevation Model. We compared outlier ratings using the robust Mahalanobis distance and Local Outlier Factor using three different thresholds (Tukey rule, histogram-based Scott's rule, and 95% quantiles in χ2 distribution). The Mahalanobis distance thresholded by the 95% χ2 quantile achieved the best results, and, because of its high specificity, appeared as a promising tool for identifying erroneously mapped or changed habitats. The presented method can, therefore, be used as a guide to target field updates of Natura 2000 habitat maps or for other habitat/land cover mapping activities where the detection of misclassifications or changes is needed. [ABSTRACT FROM AUTHOR]

The ash dieback pandemic, caused by the invasive fungus Hymenoscyphus fraxineus, represents one of Europe's biggest threats to preserving natural biodiversity. To ensure the suppression of forest damage caused by fungi, timely recognition of the symptoms of ash dieback and further continuous monitoring on an adequate spatial scale are essential. Visual crown damage assessment is currently the most common method used for identifying ash dieback, but it lacks the spatial and temporal coverage required for effective disease suppression. Remote sensing technologies, with the capabilities of fast and repetitive retrieval of information over a large spatial scale, could present efficient supplementary methods for ash damage detection and disease monitoring. In this study, we provided a synthesis of the existing remote sensing methods and applications that considers ash dieback disease, and we described the lifecycle of the disease using the major symptoms that remote sensing technologies can identify. Unfortunately, although effective methods of monitoring biotic damage through remote sensing have been developed, ash dieback has only been addressed in two research studies in the United Kingdom and Germany. These studies were based on single-date hyperspectral and very-high-resolution imagery in combination with machine learning, using previously specified ground-truth information regarding crown damage status. However, no study exists using high-resolution imagery such as Sentinel-2 or radar Sentinel-1, although some preliminary project results show that these coarser sources of information could be applicable for ash dieback detection and monitoring in cases of Fraxinus angustifolia, which forms pure, more homogenous stands in Southern Europe. [ABSTRACT FROM AUTHOR]

Corn (Zea mays L.) nitrogen (N) management requires monitoring plant N concentration (Nc) with remote sensing tools to improve N use, increasing both profitability and sustainability. This work aims to predict the corn Nc during the growing cycle from Sentinel-2 and Sentinel-1 (C-SAR) sensor data fusion. Eleven experiments using five fertilizer N rates (0, 60, 120, 180, and 240 kg N ha−1) were conducted in the Pampas region of Argentina. Plant samples were collected at four stages of vegetative and reproductive periods. Vegetation indices were calculated with new combinations of spectral bands, C-SAR backscatters, and sensor data fusion derived from Sentinel-1 and Sentinel-2. Predictive models of Nc with the best fit (R2 = 0.91) were calibrated with spectral band combinations and sensor data fusion in six experiments. During validation of the models in five experiments, sensor data fusion predicted corn Nc with lower error (MAPE: 14%, RMSE: 0.31 %Nc) than spectral band combination (MAPE: 20%, RMSE: 0.44 %Nc). The red-edge (704, 740, 740 nm), short-wave infrared (1375 nm) bands, and VV backscatter were all necessary to monitor corn Nc. Thus, satellite remote sensing via sensor data fusion is a critical data source for predicting changes in plant N status. [ABSTRACT FROM AUTHOR]

Selective logging is a major cause of forest degradation in the tropics, but its precise scale, location and timing are not known as wide-area, automated remote sensing methods are not yet available at this scale. This limits the abilities of governments to police illegal logging, or monitor (and thus receive payments for) reductions in degradation. Sentinel-1, a C-band Synthetic Aperture Radar satellite mission with a 12-day repeat time across the tropics, is a promising tool for this due to the known appearance of shadows in images where canopy trees are removed. However, previous work has relied on optical satellite data for calibration and validation, which has inherent uncertainties, leaving unanswered questions about the minimum magnitude and area of canopy loss this method can detect. Here, we use a novel bi-temporal LiDAR dataset in a forest degradation experiment in Gabon to show that canopy gaps as small as 0.02 ha (two 10 m × 10 m pixels) can be detected by Sentinel-1. The accuracy of our algorithm was highest when using a timeseries of 50 images over 20 months and no multilooking. With these parameters, canopy gaps in our study site were detected with a false alarm rate of 6.2%, a missed detection rate of 12.2%, and were assigned disturbance dates that were a good qualitative match to logging records. The presence of geolocation errors and false alarms makes this method unsuitable for confirming individual disturbances. However, we found a linear relationship ( r 2 = 0.74 ) between the area of detected Sentinel-1 shadow and LiDAR-based canopy loss at a scale of 1 hectare. By applying our method to three years' worth of imagery over Gabon, we produce the first national scale map of small-magnitude canopy cover loss. We estimate a total gross canopy cover loss of 0.31 Mha, or 1.3% of Gabon's forested area, which is a far larger area of change than shown in currently available forest loss alert systems using Landsat (0.022 Mha) and Sentinel-1 (0.019 Mha). Our results, which are made accessible through Google Earth Engine, suggest that this approach could be used to quantify the magnitude and timing of degradation more widely across tropical forests. [ABSTRACT FROM AUTHOR]

Fire is one of the significant drivers of vegetation loss and threat to Amazonian landscapes. It is estimated that fires cause about 30% of deforested areas, so the severity level is an important factor in determining the rate of vegetation recovery. Therefore, the application of remote sensing to detect fires and their severity is fundamental. Radar imagery has an advantage over optical imagery because radar can penetrate clouds, smoke, and rain and can see at night. This research presents algorithms for mapping the severity level of burns based on change detection from Sentinel-1 backscatter data in the southeastern Peruvian Amazon. Absolute, relative, and Radar Forest Degradation Index (RDFI) predictors were used through singular polarization length (dB) patterns (Vertical, Vertical-VV and Horizontal, Horizontal-HH) of vegetation and burned areas. The Composite Burn Index (CBI) determined the algorithms' accuracy. The burn severity ratios used were estimated to be approximately 40% at the high level, 43% at the moderate level, and 17% at the low level. The validation dataset covers 384 locations representing the main areas affected by fires, showing the absolute and relative predictors of cross-polarization (k = 0.734) and RDFI (k = 0.799) as the most concordant in determining burn severity. Overall, the research determines that Sentinel-1 cross-polarized (VH) data has adequate accuracy for detecting and quantifying burns. [ABSTRACT FROM AUTHOR]

Data assimilation techniques allow researchers to optimally merge remote sensing observations in ecohydrological models, guiding them for improving land surface fluxes predictions. Presently, freely available remote sensing products, such as those of Sentinel 1 radar, Landsat 8 sensors, and Sentinel 2 sensors, allow the monitoring of land surface variables (e.g., radar backscatter for soil moisture and the normalized difference vegetation index (NDVI) and for leaf area index (LAI)) at unprecedentedly high spatial and time resolutions, appropriate for heterogeneous ecosystems, typical of semiarid ecosystems characterized by contrasting vegetation components (grass and trees) competing for water use. A multiscale assimilation approach that assimilates radar backscatter and grass and tree NDVI in a coupled vegetation dynamic–land surface model is proposed. It is based on the ensemble Kalman filter (EnKF), and it is not limited to assimilating remote sensing data for model predictions, but it uses assimilated data for dynamically updating key model parameters (the ENKFdc approach), including saturated hydraulic conductivity and grass and tree maintenance respiration coefficients, which are highly sensitive parameters of soil–water balance and biomass budget models, respectively. The proposed EnKFdc assimilation approach facilitated good predictions of soil moisture, grass, and tree LAI in a heterogeneous ecosystem in Sardinia for a 3-year period with contrasting hydrometeorological (dry vs. wet) conditions. Contrary to the EnKF-based approach, the proposed EnKFdc approach performed well for the full range of hydrometeorological conditions and parameters, even assuming extremely biased model conditions with very high or low parameter values compared with the calibrated ("true") values. The EnKFdc approach is crucial for soil moisture and LAI predictions in winter and spring, key seasons for water resources management in Mediterranean water-limited ecosystems. The use of ENKFdc also enabled us to predict evapotranspiration and carbon flux well, with errors of less than 4% and 15%, respectively; such results were obtained even with extremely biased initial model conditions. [ABSTRACT FROM AUTHOR]

This paper shows the work carried out to obtain a methodology capable of monitoring the Common Agricultural Policy (CAP) aid line for the protection of steppe birds, which aims to improve the feeding and breeding conditions of these species and contribute to the improvement of their overall biodiversity population. Two methodologies were initially defined, one based on remote sensing (BirdsEO) and the other on Machine Learning (BirdsML). Both use Sentinel-1 and Sentinel-2 data as a basis. BirdsEO encountered certain impediments caused by the land's slope and the crop's height. Finally, the methodology based on Machine Learning offered the best results. It evaluated the performance of up to 7 different Machine Learning classifiers, the most optimal being RandomForest. Fourteen different datasets were generated, and the results they offered were evaluated, the most optimal being the one with more than 150 features, including a time series of 8 elements with Sentinel-1, Sentinel-2 data and derived products, among others. The generated model provided values higher than 97% in metrics such as accuracy, recall and Area under the ROC Curve, and 95% in precision and recall. The methodology is transformed into a tool that continuously monitors 100% of the area requesting aid, continuously over time, which contributes positively to optimizing the use of administrative resources and a fairer distribution of CAP funds. [ABSTRACT FROM AUTHOR]

Surface snowmelt causes changes in mass and energy balance, and endangers the stabilities of the ice shelves in the Antarctic Peninsula (AP). The dynamic changes of the snow and ice conditions in the AP were observed by Sentinel-1 images with a spatial resolution of 40min this study. Snowmelt detected by the special sensor microwave/imager (SSM/I) is used to study the relationship between summer snowmelt and winter synthetic aperture radar (SAR) backscatter. Radar glacier zones (RGZs) classifications were conducted based on their differences in liquid snow content, snow grain size, and the relative elevations. We developed a practical method based on the simulations of a microwave scattering model to classify RGZs by using Sentinel-1 images in the AP. The summer snowmelt detected by SSM/I and Sentinel-1 data are compared between 2014 and 2015. The SSM/I-derived melting days is used to validate the winter dry snow line (DSL). RGZs derived from Sentinel-1 images suggest that snowmelt expanded from inland of the Larsen C Ice Shelf to the coastal area, whereas an opposite direction was found in the George VI Ice Shelf. The long melting season in the grounding zone of the Larsen C Ice Shelf may result from the adiabatically-dried föhn winds on the east side of the AP. As the uppermost limit of summer snowmelt, DSL was mapped based on the winter Sentinel-1 mosaic of the AP. Compared with the SSM/I-derived melting days, the winter DSL mainly distributed in the areas melted for one to three days in summer. DSL elevations on the Palmer Land increased from south to north. [ABSTRACT FROM AUTHOR]

The 40 km resolution SMOS (Soil Moisture and Ocean Salinity) soil moisture, previously disaggregated at a 1 km resolution using the DISPATCH (DISaggregation based on Physical And Theoretical scale CHange) method based on MODIS optical/thermal data, is further disaggregated to 100 m resolution using Sentinel-1 backscattering coefficient (σº). For this purpose, three distinct radar-based disaggregation methods are tested by linking the spatio-temporal variability of σº and soil moisture data at the 1 km and 100 m resolution. The three methods are: (1) the weight method, which estimates soil moisture at 100 m resolution at a certain time as a function of σº ratio (100 m to 1 km resolution) and the 1 km DISPATCH products of the same time; (2) the regression method which estimates soil moisture as a function of σº where the regression parameters (e.g., intercept and slope) vary in space and time; and (3) the Cumulative Distribution Function (CDF) method, which estimates 100 m resolution soil moisture from the cumulative probability of 100 m resolution backscatter and the maximum to minimum 1 km resolution (DISPATCH) soil moisture difference. In each case, disaggregation results are evaluated against in situ measurements collected between 1 January 2016 and 11 October 2016 over a bare soil site in central Morocco. The determination coefficient (R²) between 1 km resolution DISPATCH and localized in situ soil moisture is 0.31. The regression and CDF methods have marginal effect on improving the DISPATCH accuracy at the station scale with a R² between remotely sensed and in situ soil moisture of 0.29 and 0.34, respectively. By contrast, the weight method significantly improves the correlation between remotely sensed and in situ soil moisture with a R² of 0.52. Likewise, the soil moisture estimates show low root mean square difference with in situ measurements (RMSD = 0.032 m³ m-3). [ABSTRACT FROM AUTHOR]

Soil moisture (SM) data are required at high spatio-temporal resolution—typically the crop field scale every 3–6 days—for agricultural and hydrological purposes. To provide such high-resolution SM data, many remote sensing methods have been developed from passive microwave, active microwave and thermal data. Despite the pros and cons of each technique in terms of spatio-temporal resolution and their sensitivity to perturbing factors such as vegetation cover, soil roughness and meteorological conditions, there is currently no synergistic approach that takes advantage of all relevant (passive, active microwave and thermal) remote sensing data. In this context, the objective of the paper is to develop a new algorithm that combines SMAP L-band passive microwave, MODIS/Landsat optical/thermal and Sentinel-1 C-band radar data to provide SM data at the field scale at the observation frequency of Sentinel-1. In practice, it is a three-step procedure in which: (1) the 36 km resolution SMAP SM data are disaggregated at 100 m resolution using MODIS/Landsat optical/thermal data on clear sky days, (2) the 100 m resolution disaggregated SM data set is used to calibrate a radar-based SM retrieval model and (3) the so-calibrated radar model is run at field scale on each Sentinel-1 overpass. The calibration approach also uses a vegetation descriptor as ancillary data that is derived either from optical (Sentinel-2) or radar (Sentinel-1) data. Two radar models (an empirical linear regression model and a non-linear semi-empirical formulation derived from the water cloud model) are tested using three vegetation descriptors (NDVI, polarization ratio (PR) and radar coherence (CO)) separately. Both models are applied over three experimental irrigated and rainfed wheat crop sites in central Morocco. The field-scale temporal correlation between predicted and in situ SM is in the range of 0.66–0.81 depending on the retrieval configuration. Based on this data set, the linear radar model using PR as a vegetation descriptor offers a relatively good compromise between precision and robustness all throughout the agricultural season with only three parameters to set. The proposed synergistical approach combining multi-resolution/multi-sensor SM-relevant data offers the advantage of not requiring in situ measurements for calibration. [ABSTRACT FROM AUTHOR]

This research focused on the quantitative assessment of the surface deformation velocities and rates and their natural and man-made controlling factors at Tengiz Oilfield in Kazakhstan using the Small Baseline Subset remote sensing technique followed by 3D and 2D decompositions and cosine corrections to derive vertical and horizontal movements from line-of-sight (LOS) measurements. In the present research we applied time-series of Sentinel-1 satellite images acquired during 2018–2020. All ground deformation derivatives showed the continuous subsidence at the Tengiz oilfield with increasing velocity. 3D and 2D decompositions of LOS measurements to vertical movement showed that the Tengiz Oil Field 2018–2020 continuously subsided with the maximum annual vertical deformation velocity around 70 mm. Based on the LOS measurements, the maximum annual subsiding velocity was observed to be 60 mm. Cosine corrections of LOS measurements to vertical movement, however, revealed a maximum annual vertical deformation velocity of 77 mm. The vertical deformation confirmed typical patterns of subsidence caused by oil extraction. Detected east-west and north-south horizontal movements at the Tengiz field clearly indicated that the study area crossed by seismic faults is affected by natural tectonic processes. The overall RMSE of 3D decomposed vertical deformation in relationship to LOS measurements and cosine corrections were in the range of 10–13 mm and 6–8 mm, correspondingly. The results of the present research will support operators of oil and gas fields and also other types of infrastructure to evaluate the actual differences of InSAR ground deformation measurements against the required standards and the precision of measurements depending on the operational needs, timeframes and availability of radar imagery. [ABSTRACT FROM AUTHOR]

The present paper focuses on a sensitivity analysis of Sentinel-1 backscattering signatures from oil palm canopies cultivated in Gabon, Africa. We employed one Sentinel-1 image per year during the 2015–2021 period creating two separated time series for both the wet and dry seasons. The first images were almost simultaneously acquired to the initial growth stage of oil palm plants. The VH and VV backscattering signatures were analysed in terms of their corresponding statistics for each date and compared to the ones corresponding to tropical forests. The times series for the wet season showed that, in a time interval of 2–3 years after oil palm plantation, the VV/VH ratio in oil palm parcels increases above the one for forests. Backscattering and VV/VH ratio time series for the dry season exhibit similar patterns as for the wet season but with a more stable behaviour. The separability of oil palm and forest classes was also quantitatively addressed by means of the Jeffries–Matusita distance, which seems to point to the C-band VV/VH ratio as a potential candidate for discrimination between oil palms and natural forests, although further analysis must still be carried out. In addition, issues related to the effect of the number of samples in this particular scenario were also analysed. Overall, the outcomes presented here can contribute to the understanding of the radar signatures from this scenario and to potentially improve the accuracy of mapping techniques for this type of ecosystems by using remote sensing. Nevertheless, further research is still to be done as no classification method was performed due to the lack of the required geocoded reference map. In particular, a statistical assessment of the radar signatures should be carried out to statistically characterise the observed trends. [ABSTRACT FROM AUTHOR]

Land-Use/Land-Cover (LULC) products are a common source of information and a key input for spatially explicit models of ecosystem service (ES) supply and demand. Global, continental, and regional, readily available, and free land-cover products generated through Earth Observation (EO) data, can be potentially used as relevant to ES mapping and assessment processes from regional to national scales. However, several limitations exist in these products, highlighting the need for timely land-cover extraction on demand, that could replace or complement existing products. This study focuses on the development of a classification workflow for fine-scale, object-based land cover mapping, employed on terrestrial ES mapping, within the Greek terrestrial territory. The processing was implemented in the Google Earth Engine cloud computing environment using 10 m spatial resolution Sentinel-1 and Sentinel-2 data. Furthermore, the relevance of different training data extraction strategies and temporal EO information for increasing the classification accuracy was also evaluated. The different classification schemes demonstrated differences in overall accuracy ranging from 0.88% to 4.94% with the most accurate classification scheme being the manual sampling/monthly feature classification achieving a 79.55% overall accuracy. The classification results suggest that existing LULC data must be cautiously considered for automated extraction of training samples, in the case of new supervised land cover classifications aiming also to discern complex vegetation classes. The code used in this study is available on GitHub and runs on the Google Earth Engine web platform. [ABSTRACT FROM AUTHOR]