Ground Range Detected Research Articles

Integrating Multisource Data and Machine Learning for Supraglacial Lake Detection: Implications for Environmental Management and Sustainable Development Goals in High Mountainous Regions

The accurate detection and monitoring of supraglacial lakes in high mountainous regions are crucial for understanding their dynamic nature and implications for environmental management and sustainable development goals. In this study, we propose a novel approach that integrates multisource data and machine learning techniques for supra-glacial lake detection in the Passu Batura glacier of the Hunza Basin, Pakistan. We extract pertinent features or parameters by leveraging multisource datasets such as radar backscatter intensity VH and VV parameters from Sentinel-1 Ground Range Detected (GRD) data, near-infrared (NIR), NDWI_green, NDWI_blue parameters from Sentinel-2 Multi-spectral Instrument(MSI) data, and surface slope, aspect, and elevation parameters from topographic data. The entire dataset is partitioned into training and testing sets, with machine learning models including the artificial neural network (ANN), the support vector machine (SVM), logistic regression (LR), random forest (RF), and K-nearest neighbour (KNN) trained on the training data (70%). Accuracy assessment employs testing data and involves the evaluation of metrics such as ROC curves and confusion matrices. The best-performing model, ANN, is validated against manually digitized lake polygons derived from Sentinel-2 and Google Earth Pro imagery. Furthermore, the digitized lake polygons are used to analyze glacial lake dynamics from 2016 to 2022. Key findings of this research presented that the NDWI_green, Sigma0_VH, and elevation are the most significant predictors in detecting supra-glacial lakes. Among the various trained and evaluated models, the Artificial Neural Network (ANN) achieved the highest performance (accuracy: 95%, AUC: 0.99) and accurately mapped supra-glacial lakes regardless of their small size. The findings have significant implications for understanding glacial lake behavior in the context of climate change and informing future research and monitoring efforts.

Early Season Crop Acreage Estimation for Pigeon Pea using SAR Data: A Case Study of Kalaburagi District, Karnataka

Crop acreage estimates using space-based observations during mid-season and pre-harvest periods are operationally provided for major Kharif and Rabi crops in India. These are critical inputs for planning and decision-making concerning food security, storage, pricing, and procurement. This work explores the potential of Satellite C-band Synthetic Aperture Radar (SAR) data in VV and VH polarization, acquired during monsoon season, for early crop acreage assessment of Pigeon Pea. The present study evaluates two classification approaches based on multi-temporal SAR data to generate a spatial distribution map of the pigeon pea crop in the Kalaburagi region of Karnataka. The temporal profiles of sigma nought (backscatter coefficient) and Entropy vs. Alpha Angle were studied with reference to crop phenology. Random forest algorithm has been used to classify multi-temporal SAR Ground Range Detected (GRD) data from June to November 2019. Temporal SAR data from June to early November with VV+VH polarization was optimal for tur area assessment with 85% accuracy. This study also examines phase and intensity information from 3-date SAR Single Look Complex (SLC) data with an unsupervised Wishart Classification approach for improved crop acreage estimation. The results emphasize the effectiveness of dual-polarimetric SAR data for timely crop inventory of Pigeon Pea during the Kharif season.

Dynamic flood mapping by a normalized probabilistic classification method using satellite radar amplitude time series

ABSTRACT Owing to the vast development of Synthetic Aperture Radar (SAR), especially the improvement of spatio-temporal resolution, observing and quantifying the complex and dynamic flood process becomes increasingly feasible. Utilizing the Sentinel-1 Ground Range Detected (GRD) dataset, we proposed an improved probabilistic flood mapping method combining image Pareto Scaling (PS) normalization and Bayesian probability estimation. We validated our method during a flood event in Xinjiang County, Shaanxi Province of China in October 2021 using a high spatial resolution (0.1 m) Unmanned Aerial Vehicle (UAV) image. The overall reliability of the new method agrees 95% to the UAV measurements and achieves the highest accuracy (85.2%) when compared to the Sentinel-1 dual-polarized water index (SDWI) threshold method and the Z-score method. Our results distinguished four flood stages: flood emergence, peak, receding, and disappearance, which provide valuable insights into the dynamic change process of floods. Notably, we observed that pixels with different flood probabilities exhibited distinct temporal characteristics. The extremely high probability pixel experienced rapid fluctuations, while the medium probability pixel showed more gradual changes over time. We believe our proposed method can enhance our understanding of flood-prone areas and their dynamics so that decision-makers can develop targeted mitigation measures and response plans.

ResNet Deep Learning‐Based Inversion Method for Sea Surface Wind Field

This paper presents a method to invert the sea surface wind field from Sentinel‐1 synthetic aperture radar (SAR) interferometric wide swath (IW) band mode ground range detected (GRD) vertical polarization (VV) data based on the deep residual network (ResNet). Compared with the traditional algorithm (i.e., using geophysical model functions (GMFs)), the wind field inversion based on the ResNet deep learning method solves the problem of separate wind direction and wind speed inversion. The ResNet model is developed based on the results of 212 VV polarized IW images and cross calibrated multiplatform (CCMP) data collected in 2021 in the offshore waters of China. The root mean square error, standard deviation, and bias of the inverse wind direction of the ResNet model are 12.5°, 9.8°, and 7.4°, respectively, and the root mean square error, standard deviation, and bias of the inverse wind speed are 2.33, 1.495, and −1.375. The ResNet is better than the conventional local gradient and C‐band model models, further confirming the advantages shown by the deep learning algorithm in nonlinear fitting and its potential application in sea surface wind field inversion.

Flood susceptibility mapping using Sentinel 1 and frequency ratio technique in Jinjiram River watershed, India.

Assam is one of the most flood-prone states in India, and the state frequently experiences catastrophic floods that cause significant damage in terms of loss of life and property. Flood susceptibility is considered the most essential and crucial input for managing floodplains and fostering local and regional development. This study focuses on the generation of flood susceptibility maps using the Frequency Ratio (FR) technique and microwave remote sensing inputs in the Jinjiram watershed which experienced disastrous flooding in 2020. The study has been carried out by taking into consideration of different morphological, lithological, and hydrological factors. In this study the flood inventory map was created by extracting the time series SAR (Synthetic Aperture Radar) Sentinel 1 GRD (Ground Range Detected) images of flooded areas for the past 5years, from 2016 to 2020. A total of 72 inventory samples were identified of which 70% of total flooded samples were chosen for training and 30% for model testing at random basis. Applying these FR methods, the study determines a range of flood susceptibility which was then divided into five classes, from very low to very high. The Receiver Operating Characteristics (ROC) analysis was used to evaluate the accuracy of flood susceptibility maps generated using FR models. The AUC of ROC in flood susceptibility mapping is 0.81167 achieved, corresponding to a prediction accuracy of 81.17%. The findings can be used to calculate risk, develop flood control measures and infrastructural policies, and formulate sustainable water management policies for the watershed.

Evaluating SAR Radiometric Terrain Correction Products: Analysis-Ready Data for Users

Operational applications for Synthetic Aperture Radar (SAR) are under development around the world, driven by the free-and-open access of SAR C-band observations that Sentinel-1 of Copernicus has provided since 2014. Radiometric Terrain Correction (RTC) data are key entry-level products for multiple applications ranging from ecosystem to hazard monitoring. Various open-source software packages exist to create RTC products from Single Look Complex (SLC) or Ground Range Detected (GRD) level SAR data, including the Interferometric SAR Computing Environment (ISCE), and the Sentinel-1 Toolbox from the European Space Agency (SNAP 8). Despite the growing availability of RTC software solutions, little work has been performed to identify differences between RTC products generated using different software packages. This work evaluates several Sentinel-1 RTC products and two other Sentinel-1 Analysis Ready Data (ARD) to address the following questions: (1) Which software provides the most accurate RTC product? and (2) how appropriate for analysis are other non-RTC products that are readily available? The RTCs are produced with GAMMA, ISCE-2, and SNAP 8. The other two ARD products evaluated consisted of an angular-based radiometric slope correction produced in Google Earth Engine (GEE) following Vollrath et al., and the Sentinel-1 GRD product. Products are evaluated across 10 sites in a single image approach for (1) radiometric calibration, (2) geometric corrections, and for (3) geolocation quality. In addition, time-series stacks over two sites representing varied terrain and ecosystems are evaluated. The GAMMA-derived RTC product implemented by the Alaska Satellite Facility (ASF) is used as a reference for some of the time-series metrics. The results provide direct guidance and recommendations about the quality of the RTC and ARD products obtained from open source methods. The results indicate that it is not recommended to use the GRD product with no radiometric or geometric corrections for any applications given low performance in multiple metrics. The radiometric calibration and geometric corrections have overall good performance for all open-source solutions, only the non-RTC products (Vollrath et al. and GRD) portray some significant variances in steep terrain. The geolocation assessment indicated that the GRD product has the most significant displacement errors, followed by SNAP 8 with Digital Elevation Model (DEM) matching, and ISCE-2. RTCs created without DEM-matching performed better for both GAMMA and SNAP 8. The time-series results indicate that SNAP 8 products align more closely to GAMMA products than other open-source software in terms of radiometric and geometric quality. This understanding of software performance for SAR image processing is key to designing the affordable and scalable solutions needed for the operational application of SAR Sentinel-1 data.

Detection of Floating Debris in the Lake Using Statistical Properties of Synthetic Aperture Radar Pulses

This study developed the European Space Agency (ESA) Setinel-1 Ground Range Detected (GRD) time series analysis model for monitoring floating debris in lake areas through Google Earth Engine Application Programming Interface. The study aims to monitor floating debris caused by heavy rainfall efficiently. Regarding water resources and water quality management, floating debris from multipurpose dams requires continuous monitoring from the initial generation stage. In the study, a Synthetic Aperture Radar (SAR) time series analysis model that is easy to identify water bodies was developed due to low accessibility in large areas. Although SAR satellite images could be used to observe inland water environments, debris detection on water surface surfaces has yet to be studied. For the first time, this study detected floating debris patches in a wide range of lakes from GRD imagery acquired by ESA’s Sentinel-1 satellite. It demonstrated the potential to distinguish them from naturally occurring materials such as invasive floating plants. In this study, the case of Daecheong Dam, in which predicted floating debris was detected after heavy rain using Sentinel-1 GRD data, is presented. It could quickly detect various floating debris flowing into dams used as a source of drinking water and serve as a reference for establishing a collection plan.

ReCuSum: A polyvalent method to monitor tropical forest disturbances

Change detection methods based on Earth Observations are increasingly used to monitor rainforest in the intertropical band. Until recently, deforestation monitoring was mainly based on remotely sensed optical images which often face limitations in humid tropical areas due to frequent cloud coverage. This leads to late detections of disturbance events. Since the launch of Sentinel-1 acquiring images with a revisit time of 12 days and a spatial resolution of 5 × 20 m in Brazil, Synthetic Aperture Radar (SAR) images have been increasingly used to monitor deforestation. In this study, we propose a multitemporal version of the change detection method we previously applied to timeseries of Sentinel-1 SAR images, to monitor deforestation/degradation in the Congo rainforest. Our approach is based on a Cumulative Sum (CuSum) method combined with a spatial recombination of Critical Thresholds (CuSum cross-Tc). The newly developed multitemporal CuSum method (ReCuSum) was applied to a time-series of 82 dual polarization (VH, VV) Ground Range Detected (GRD) Sentinel-1 images acquired in the Para State, in the Brazilian Amazonia, from 29/09/2016 to 01/07/2019. The ReCuSum method consists of iteratively applying the CuSum cross-Tc to monitor multiple changes in a time-series by splitting the time-series at each date of detected change and by independently iterating over the time periods resulting from the splits. The number of changes in the time-series was then analysed according to the vegetation type (Forest, non-forest vegetation) determined by visual inspection of optical Sentinel-2 image and PlanetScope monthly mosaic. This showed a difference between non-forest vegetation and forested areas. A threshold based on the number of changes (Tnbc) was then developed to differentiate forest from non-forest disturbances. The ability to monitor non-forest vegetation was analysed: the CuSum cross-Tc detected up to 90% of the total non-forest vegetation area over the study region in the past period. After removing past disturbances and past non-forest vegetation, then removing the pixels covered with non-forest vegetations based on Tnbc, the application of the ReCuSum led to a precision of 81%, a recall of 68%, a kappa coefficient of 0.72 and a F1- score of 0.74 in forest disturbance monitoring. According to these results, ReCuSum applied to Sentinel-1 time-series of images can be used for efficient forest disturbance monitoring and for generating a forest / non-forest map after the application of newly developed post-processing steps. Sentinel-1 imagery can be used for both Forest / Non-forest mapping and for forest disturbance detection. ReCuSum was released as an open-source GIT project available at: https://forgemia.inra.fr/bertrand.ygorra/cusum-deforestation_monitoring.

Measuring Vertical Urban Growth of Patna Urban Agglomeration Using Persistent Scatterer Interferometry SAR (PSInSAR) Remote Sensing

In the present study, the vertical and horizontal growth of Patna Urban Agglomeration was evaluated using the Persistent Scatterer Interferometry Synthetic Aperture Radar (PSInSAR) technique during 2015–2018. The vertical urban growth assessment of the city landscape was assessed using microwave time series (30 temporal) datasets of Single Look Complex (SLC) Sentinel-1A interferometric Synthetic Aperture Radar using SARPROZ software (ver. 2020). This study demonstrated that peripheral city regions experienced higher vertical growth (~4 m year−1) compared to the city core regions, owing to higher urban development opportunities leading to significant land use alterations, the development of high-rise buildings, and infrastructural development. While the city core of Patna observed an infill and densification process, as it was already saturated and highly densified. The rapidly urbanizing city in the developing region witnessed a considerable horizontal urban expansion as estimated through the normalized difference index for built-up areas (NDIB) and speckle divergence (SD) using optical Sentinel 2A and microwave Sentinel-1A ground range detected (GRD) satellite data, respectively. The speckle divergence-based method exhibited high urban growth (net growth of 11.28 km2) with moderate urban infill during 2015–2018 and reported a higher accuracy as compared to NDIB. This study highlights the application of SAR remote sensing for precise urban area delineation and temporal monitoring of urban growth considering horizontal and vertical expansion through processing a long series of InSAR datasets that provide valuable information for informed decision-making and support the development of sustainable and resilient cities.

Automatic lineaments detection using radar and optical data with an emphasis on geologic and tectonic implications: a case study of Kabul Block, eastern Afghanistan

The earth’s surface linear features, expressing geological lineaments, play a key role in identifying hydrothermal alteration and mineralization zones, as well as in understanding tectonic settings of a region. The objective of this investigation is to utilize a method for extracting lineaments automatically, which will be integrated and applied to identify geological-based lineaments by making use of remote sensing data. The study will further examine the structural pattern and temporal-spatial evolution of the lineaments and establish their connection with the primary active faults present in the Kabul Block. Multi-sensor data from radar (DEM (Digital Elevation Model)-5m, Sentinel-1B GRD (Ground Range Detected)) and optical sensors (Sentinel-2 MSI (Multispectral Instrument) and ASTER) were processed using spatial and spectral filters before automatic lineament extraction. LINE-module algorithm was applied to various illuminated hill-shades of DEM-5m, filtered HH (Horizonal-Horizonal) and VH (Vertical-Horizontal) of Sentinel-1GRD and PC1 (Principle Component) of Sentinel-2A MSI, and ASTER data to detect linear surface features. Extensive testing was conducted to verify the accuracy of extracted lineaments and to exclude any artificial lineaments in the study area. The radar and optical data results were compared while taking into consideration the geological and tectonic settings of the study area to select the most appropriate extracted lineaments. The DEM-5m and Sentinel-1B GRD showed the best result for identifying lineaments, and these were found to be highly correlated with previously available data in the Kabul Block. The final results of DEM-5m and Sentinel-1B GRD were further analyzed. The extracted lineaments were found to trend predominantly in a NW-SE and NE-SW directions, which is consistent with the results of other data. Temporal evolution and spatial distribution reveal that a high density of the lineaments is associated with Paleogene and Quaternary formations, while a low density is observed in Proterozoic, Paleozoic, and Mesozoic formations. The west and southwest edges of the Kabul Block are controlled by compressive stress trending NNW-SSE, while the north and southeast margins are influenced by NE-SW and ∼ N-S trending compressive stress.

Assessing riverbank erosion in Bangladesh using time series of Sentinel-1 radar imagery in the Google Earth Engine

Abstract. Riverbank erosion occurs along many of the Earth's river systems, affecting riverine populations by destroying agricultural land and housing. In this study, we detected past events of riverbank erosion along the Jamuna River in Bangladesh using time series of Sentinel-1 satellite radar imagery, ground-range-detected (GRD) data with a 12 d revisit cycle, available in the Google Earth Engine (GEE). Eroded land is detected by performing a land cover classification and by detecting land cover changes from vegetated areas before the monsoon to sand or water after the monsoon. Further, settlements are detected as persistent scatterers and classified as eroded if they are located on eroded land. We found that with Sentinel-1 data, erosion locations can be determined already 1 month after the end of the monsoon and hence potentially earlier than using optical satellite images which depend on cloud-free daylight conditions. Further, we developed an interactive GEE-based online tool allowing the user to explore where riverbank erosion has destroyed land and settlements along the Jamuna in five monsoon seasons (2015–2019). The source code of our implementation is publicly available, providing the opportunity to reproduce the results, to adapt the algorithm and to transfer our results to assess riverbank erosion in other geographical settings.

Fusion of VNIR Optical and C-Band Polarimetric SAR Satellite Data for Accurate Detection of Temporal Changes in Vegetated Areas

In this paper, we propose a processing chain jointly employing Sentinel-1 and Sentinel-2 data, aiming to monitor changes in the status of the vegetation cover by integrating the four 10 m visible and near-infrared (VNIR) bands with the three red-edge (RE) bands of Sentinel-2. The latter approximately span the gap between red and NIR bands (700 nm–800 nm), with bandwidths of 15/20 nm and 20 m pixel spacing. The RE bands are sharpened to 10 m, following the hyper-sharpening protocol, which holds, unlike pansharpening, when the sharpening band is not unique. The resulting 10 m fusion product may be integrated with polarimetric features calculated from the Interferometric Wide (IW) Ground Range Detected (GRD) product of Sentinel-1, available at 10 m pixel spacing, before the fused data are analyzed for change detection. A key point of the proposed scheme is that the fusion of optical and synthetic aperture radar (SAR) data is accomplished at level of change, through modulation of the optical change feature, namely the difference in normalized area over (reflectance) curve (NAOC), calculated from the sharpened RE bands, by the polarimetric SAR change feature, achieved as the temporal ratio of polarimetric features, where the latter is the pixel ratio between the co-polar and the cross-polar channels. Hyper-sharpening of Sentinel-2 RE bands, calculation of NAOC and modulation-based integration of Sentinel-1 polarimetric change features are applied to multitemporal datasets acquired before and after a fire event, over Mount Serra, in Italy. The optical change feature captures variations in the content of chlorophyll. The polarimetric SAR temporal change feature describes depolarization effects and changes in volumetric scattering of canopies. Their fusion shows an increased ability to highlight changes in vegetation status. In a performance comparison achieved by means of receiver operating characteristic (ROC) curves, the proposed change feature-based fusion approach surpasses a traditional area-based approach and the normalized burned ratio (NBR) index, which is widespread in the detection of burnt vegetation.

A Generalized Geodesic Distance-Based Approach for Analysis of SAR Observations Across Polarimetric Modes

Present and future sensors are diversifying from traditional quad polarimetric mode of synthetic aperture radar acquisition. Thus, an approach that is interpretative in nature and applicable across polarimetric modes is required. In this context, the geodesic distance (GD)-based approach within the polarimetric synthetic aperture radar (PolSAR) literature is seen as an eigenvalue-decomposition free approach to interpret and analyze quad PolSAR data. This approach is highly adaptive toward applications due to its ability to compare the SAR observation with a known scatterer/model, or with another SAR observation in general providing a means for direct interpretation. In this work, we show that the GD (originally defined for the quad polarization mode) is generalizable across any arbitrary SAR polarimetric mode while retaining its simple form for ready computation. We show that the GD-based approach provides level ground for comparison of different polarimetric modes given a fixed application. We demonstrate it using change detection as the chosen application. In addition, we show how the behavior of the three roll-invariant GD-based parameters change under different polarimetric modes (e.g., quad, dual, and compact polarization modes). We also discuss how the GD-based approach can also be adapted to ground range detected (GRD) product data, which is presently available from Sentinel-1 and widely used in many applications. However, in this case, we show that only one of the three GD-derived parameters can be defined. We believe this work will make the GD-based approach important for present and future PolSAR applications cutting across sensors and its available polarimetric modes.

Surface Water Mapping from SAR Images Using Optimal Threshold Selection Method and Reference Water Mask

Water resources are an important component of ecosystem services. During long periods of cloudiness and precipitation, when a ground-based sample is not available, the water bodies are detected from satellite SAR (synthetic-aperture radar) data using threshold methods (e.g., Otsu and Kittler–Illingworth). However, such methods do not enable to obtain the correct threshold value for the backscattering coefficient (σ0) of relatively small water areas in the image. The paper proposes and substantiates a method for the mapping of the surface of water bodies, which makes it possible to correctly identify water bodies, even in “water”/“land” class imbalance situations. The method operates on a principle of maximum compliance of the resulting SAR water mask with a given reference water mask. Therefore, the method enables the exploration of the possibilities of searching and choosing the optimal parameters (polarization and speckle filtering), which provide the maximum quality of SAR water mask. The method was applied for mapping natural and industrial water bodies in the Pohjois-Pohjanmaa region (North Ostrobothnia), Finland, using Sentinel-1A and -1B ground range detected (GRD) data (ascending and descending orbits) in 2018–2021. Reference water masks were generated based on optical spectral indices derived from Sentinel-2A and -2B data. The polarization and speckle filtering parameters were chosen since they provide the most accurate σ0 threshold (on average for all observations above 0.9 according to the Intersection over Union criterion) and are resistant to random fluctuations. If a reference water mask is available, the proposed method is more accurate than the Otsu method. Without a reference mask, the σ0 threshold is calculated as an average of thresholds obtained from previous observations. In this case, the proposed method is as good in accuracy as the Otsu method. It is shown that the proposed method enables the identification of surface water bodies under significant class imbalance conditions, such as when the water surface covers only a fraction of a percent of the area under study.

Probabilistic SAR-based water segmentation with adapted Bayesian convolutional neural network

Geospatial resources, including satellite-based synthetic aperture radar (SAR) and optical data, have been instrumental in providing time-sensitive information about the extent and impact of natural hazard events, such as floods, to support emergency response and hazard management efforts. In effect, finite resources can be better optimized to support the needs of often extensively affected areas. However, the derivation of SAR-based flood information is not without its challenges and inaccurate flood detection can result in poor decision-making and non-trivial, adverse consequences. Reliable uncertainty quantification of flood extent estimates addresses this risk. In this context, our study presents the results of two probabilistic convolutional neural networks (CNNs) adapted for SAR-based water segmentation with freely available Sentinel-1 interferometric wide (IW) swath ground range detected (GRD) data. In particular, the performance of a variational inference-based Bayesian convolutional neural network (BCNN) is evaluated against that of a Monte Carlo dropout network (MCDN). In previous studies where uncertainty information has been generated along with segmentation results, MCDN has been more commonly applied as an approximation of Bayesian deep learning. Differences between the two approaches are highlighted with the application of a set of extended performance metrics. Both segmentation and uncertainty outputs are evaluated at data-, model-, tile- and scene-levels. The methods are demonstrated on the binary segmentation of a reference water dataset. We show that while different probabilistic techniques return comparable segmentation accuracies, they are differentiated based on their performance in assigning reliable probabilities. In particular, MCDNs characterized by more restrictive architectures generally lead to overconfident prediction intervals, whereas BCNNs have greater flexibility to learn the mean and the spread of the parameter posterior. Furthermore, we demonstrate that examining the (inaccurate, certain) metric is a better indicator of reliable uncertainty quantification and the BCNN is recommended to quantify uncertainties associated with SAR-based segmentation outputs. This information is especially valuable where the cost of inaccurate detections (false-positive and false-negative) is high. These preliminary findings highlight how the manner in which probabilities are properly assigned and their inclusion are instrumental and complementary to the production of flood masks and should be considered as a standard in the natural hazards domain.

Reliably mapping low-intensity forest disturbance using satellite radar data

In the last decades tropical forests have experienced increased fragmentation due to a global growing demand for agricultural and forest commodities. Satellite remote sensing offers a valuable tool for monitoring forest loss, thanks to the global coverage and the temporal consistency of the acquisitions. In tropical regions, C-band Synthetic Aperture Radar (SAR) data from the Sentinel-1 mission provides cloud-free and open imagery on a 6- or 12-day repeat cycle, offering the unique opportunity to monitor forest disturbances in a timely and continuous manner. Despite recent advances, mapping subtle forest losses, such as those due to small-scale and irregular selective logging, remains problematic. A Cumulative Sum (CuSum) approach has been recently proposed for forest monitoring applications, with preliminary studies showing promising results. Unfortunately, the lack of accurate in-situ measurements of tropical forest loss has prevented a full validation of this approach, especially in the case of low-intensity logging. In this study, we used high-quality field measurements from the tropical Forest Degradation Experiment (FODEX), combining unoccupied aerial vehicle (UAV) LiDAR, Terrestrial Laser Scanning (TLS), and field-inventoried data of forest structural change collected in two logging concessions in Gabon and Peru. The CuSum algorithm was applied to VV-polarized Sentinel-1 ground range detected (GRD) time series to monitor a range of canopy loss events, from individual tree extraction to forest clear cuts. We developed a single change metric using the maximum of the CuSum distribution, retrieving location, time, and magnitude of the disturbance events. A comparison of the CuSum algorithm with the LiDAR reference map resulted in a 78% success rate for the test site in Gabon and 65% success rate for the test site in Peru, for disturbances as small as 0.01 ha in size and for canopy height losses as fine as 10 m. A correlation between the change metric and above ground biomass (AGB) change was found with R2 = 0.95, and R2 = 0.83 for canopy height loss. From the regression model we directly estimated local AGB loss maps for the year 2020, at 1 ha scale and in percentages of AGB loss. Comparison with the Global Forest Watch (GFW) Tree Cover Loss (TCL) product showed a 61% overlap between the two maps when considering only deforested pixels, with 504 ha of deforestation detected by CuSum vs. 348 ha detected by GFW. Low intensity disturbances captured by the CuSum method were largely undetected by GFW and by the SAR-based Radar for Detecting Deforestation (RADD) Alert System. The results of this study confirm this approach as a simple and reproducible change detection method for monitoring and quantifying fine-scale to high intensity forest disturbances, even in the case of multi-storied and high biomass forests.

Monitoring of Plastic Islands in River Environment Using Sentinel-1 SAR Data

Plastics in the river environment are of major concern due to their potential pathways into the ocean, their persistence in the environment, and their impacts on human and marine health. It has been documented that plastic concentrations in riparian environments are higher following major rain events, where plastic can be moved through surface runoff. Considering the hazard that plastic waste poses to the environment, monitoring techniques are needed to aid in locating, monitoring, and remediating plastic waste within these systems. Dams are known to trap sediments and pollutants, such as metals and Polychlorinated Biphenyls (PCBs). While there is an established background on the monitoring of dams using the synoptic coverage provided by satellite imaging to observe water quality and volume, the detection of marine debris in riparian systems remains challenging, especially in cloudy conditions. Herein, we exploit the use of Synthetic Aperture Radar (SAR) to understand its capabilities for monitoring marine debris. This research focuses on detecting plastic islands within the Drina River system in Bosnia and Herzegovina and Serbia. Here, the results show that the monitoring of these plastic accumulations is feasible using Sentinel-1 SAR data. A quantitative analysis of detection performance is presented using traditional and state-of-the-art change detectors. The analysis of these detectors indicates that detectors that can utilise the coherent data from Single Look Complex (SLC) acquisitions are perform better when compared with those that only utilise incoherent data from Ground Range-Detected (GRD) acquisitions, with true positive detection ratings of ~95% with 0.1% false alarm rates seen in the best-performing detector. We also found that that the cross-pol VH channel provides better detection than those based on single-pol VV polarisation.

Field-scale soil moisture estimation using sentinel-1 GRD SAR data

Soil moisture is a critical land variable that controls the energy and mass balance in land–atmosphere interactions. Spaceborne Synthetic Aperture Radar (SAR) sensors offer an efficient way to map and monitor soil moisture because of their sensitivity towards the dielectric and geometric properties of the target. In addition, SAR acquisitions are weather-independent, providing a significant advantage over optical imaging during periods of cloud cover. However, vegetation cover makes these processes more complex and influences the interaction of SAR backscatter resulting from combined soil matrix and vegetation cover. Therefore, using SAR data, it is necessary to compensate for vegetation contribution in total backscatter while estimating soil moisture over the vegetated soil surface. This study presents a technique that utilizes a vegetation index derived from SAR data to generate high-resolution soil moisture maps. It is noteworthy that this proposed soil moisture retrieval method uses only the dual-polarimetric Ground Range Detected (GRD) SAR product, i.e., only backscatter intensities. Hence, the proposed method has a high potential for operational soil moisture monitoring globally. We validated over 34 soil moisture stations of the Texas Soil Observation Network (TxSON) using time-series Sentinel-1 SAR data. The Root Mean Square Error (RMSE) values for estimated volumetric soil moisture are within the range of 0.048 m3 m−3 to 0.055 m3 m−3 with the Pearson correlation coefficient r>0.79. The code to generate DpRVIc in Google Earth Engine is available at: https://github.com/Narayana-Rao/dual_pol_descriptors.

Soil moisture retrieval over croplands using dual-pol L-band GRD SAR data

Vegetation cover significantly influences the hydrometeorological processes of land surfaces. The heterogeneity of vegetation cover makes these processes more complex and impacts the interaction between water held in the soil matrix and vegetation cover. The backscatter measured by Synthetic Aperture Radar (SAR) is sensitive to target dielectric and morphological properties. In addition, SAR acquisitions are weather-independent, providing a major advantage over optical imaging during periods of cloud cover. Most often, vegetation properties are measured using vegetation indices, including simple ratios of backscatter intensities. However, this approach can miss the scattered wave purity from vegetation targets. The scattered wave characteristic is an essential parameter as it is susceptible to the complex random structure of vegetation. In this study, we propose a novel dual-polarimetric SAR vegetation descriptor based on the co-pol purity component of the wave. We use this descriptor within a semi-empirical vegetation model to estimate soil moisture. It is noteworthy that this proposed method to retrieve soil moisture uses only the dual-polarimetric Ground Range Detected (GRD) SAR product, i.e., only backscatter intensities. Therefore, the proposed method has the potential for operational scale monitoring applications at a global scale. This study validated different crop types (viz. canola, oats, forage, maize, soybeans, sugar beet, wheat, winter wheat, winter barley) over two test sites using airborne SAR data. The Root Mean Square Error (RMSE) values for soil moisture are within the range of 4.3% to 7.7% with the Pearson correlation coefficient r greater than 0.6. This novel descriptor will facilitate the operational use of dual-polarimetric GRD SAR data for vegetation monitoring and soil moisture estimation.The code to generate DpRVIc in Google Earth Engine is available at: https://github.com/Narayana-Rao/dual_pol_descriptors.

An Operational Analysis Ready Radar Backscatter Dataset for the African Continent

Digital Earth Africa is now providing an operational Sentinel-1 normalized radar backscatter dataset for Africa. This is the first free and open continental scale analysis ready data of this kind that has been developed to be compliant with the CEOS Analysis Ready Data for Land (CARD4L) specification for normalized radar backscatter (NRB) products. Partnership with Sinergise, a European geospatial company and Earth observation data provider, has ensured this dataset is produced efficiently in the cloud infrastructure and can be sustained in the long term. The workflow applies radiometric terrain correction (RTC) to the Sentinel-1 ground range detected (GRD) product, using the Copernicus 30 m digital elevation model (DEM). The method is used to generate data for a range of sites around the world and has been validated as producing good results. This dataset over Africa is made available publicly as a AWS public dataset and can be accessed through the Digital Earth Africa platform and its Open Data Cube API. We expect this dataset to support a wide range of applications, including natural resource monitoring, agriculture, and land cover mapping across Africa.