Simplified Normalization of C-Band Synthetic Aperture Radar Data for Terrestrial Applications in High Latitude Environments

Quantifying Irrigated Winter Wheat LAI in Argentina Using Multiple Sentinel-1 Incidence Angles

Synthetic aperture radar (SAR) data provides an appealing opportunity for all-weather day or night Earth surface monitoring. The European constellation Sentinel-1 (S1) consisting of S1-A and S1-B satellites offers a suitable revisit time and spatial resolution for the observation of croplands from space. The C-band radar backscatter is sensitive to vegetation structure changes and phenology as well as soil moisture and roughness. It also varies depending on the local incidence angle (LIA) of the SAR acquisition's geometry. The LIA backscatter dependency could therefore be exploited to improve the retrieval of the crop biophysical variables. The availability of S1 radar time-series data at distinct observation angles holds the feasibility to retrieve leaf area index (LAI) evolution considering spatiotemporal coverage of intensively cultivated areas. Accordingly, this research presents a workflow merging multi-date S1 smoothed data acquired at distinct LIA with a Gaussian processes regression (GPR) and a cross-validation (CV) strategy to estimate cropland LAI of irrigated winter wheat. The GPR-S1-LAI model was tested against in situ data of the 2020 winter wheat campaign in the irrigated valley of Colorador river, South of Buenos Aires Province, Argentina. We achieved adequate validation results for LAI with R C V 2 = 0.67 and RMSE _CV = 0.88 m² m^-2. The trained model was further applied to a series of S1 stacked images, generating temporal LAI maps that well reflect the crop growth cycle. The robustness of the retrieval workflow is supported by the associated uncertainties along with the obtained maps. We found that processing S1 smoothed imagery with distinct acquisition geometries permits accurate radar-based LAI modeling throughout large irrigated areas and in consequence can support agricultural management practices in cloudprone agri-environments.

Spatial Dispersion and Non-Negative Matrix Factorization of SAR Backscattering as Tools for Monitoring Snow Depth Evolution in Mountain Areas: A Case Study at Central Pyrenees (Spain)

Accurate knowledge of snow cover extent, depth (SD), and water equivalent is essential for studying the global water cycle, climate, and energy–mass exchange in the Earth–atmosphere system, as well as for water resources management. Ratio between SAR cross- and co-polarization backscattering (σVH/σVV) was used to monitor SD during snowy months in mountain areas; however, published results refer to short periods and show lack of correlation during non-snowy months. We analyze Sentinel-1A images from a study area in Central Pyrenees to generate and investigate (i) time series of σVH/σVV spatial dispersion, (ii) spatial distribution of pixelwise σVH/σVV temporal standard deviation, and (iii) fundamental modes of σVH/σVV evolution by non-negative matrix factorization. The spatial dispersion evolution and the first mode are highly correlated (correlation coefficients larger than 0.9) to SD evolution during the whole seven-year-long period, including snowy and non-snowy months. The local incidence angle strongly affects how accurately σVH/σVV locally follows the first mode; thus, areas where it predominates are orbit-dependent. When combining ascending- and descending-orbit images in a single data matrix, the first mode becomes predominant almost everywhere snow pack persists during winter. Capability of our approach to reproduce SD evolution makes it a very effective tool.

Land Cover-Specific Local Incidence Angle Correction: A Method for Time-Series Analysis of Forest Ecosystems

This study deals with a local incidence angle correction method, i.e., the land cover-specific local incidence angle correction (LC-SLIAC), based on the linear relationship between the backscatter values and the local incidence angle (LIA) for a given land cover type in the monitored area. Using the combination of CORINE Land Cover and Hansen et al.’s Global Forest Change databases, a wide range of different LIAs for a specific forest type can be generated for each scene. The algorithm was developed and tested in the cloud-based platform Google Earth Engine (GEE) using Sentinel-1 open access data, Shuttle Radar Topography Mission (SRTM) digital elevation model, and CORINE Land Cover and Hansen et al.’s Global Forest Change databases. The developed method was created primarily for time-series analyses of forests in mountainous areas. LC-SLIAC was tested in 16 study areas over several protected areas in Central Europe. The results after correction by LC-SLIAC showed a reduction of variance and range of backscatter values. Statistically significant reduction in variance (of more than 40%) was achieved in areas with LIA range >50° and LIA interquartile range (IQR) >12°, while in areas with low LIA range and LIA IQR, the decrease in variance was very low and statistically not significant. Six case studies with different LIA ranges were further analyzed in pre- and post-correction time series. Time-series after the correction showed a reduced fluctuation of backscatter values caused by different LIAs in each acquisition path. This reduction was statistically significant (with up to 95% reduction of variance) in areas with a difference in LIA greater than or equal to 27°. LC-SLIAC is freely available on GitHub and GEE, making the method accessible to the wide remote sensing community.

Towards Circumpolar Mapping of Arctic Settlements and Infrastructure Based on Sentinel-1 and Sentinel-2

Infrastructure expands rapidly in the Arctic due to industrial development. At the same time, climate change impacts are pronounced in the Arctic. Ground temperatures are, for example, increasing as well as coastal erosion. A consistent account of the current human footprint is needed in order to evaluate the impact on the environments as well as risk for infrastructure. Identification of roads and settlements with satellite data is challenging due to the size of single features and low density of clusters. Spatial resolution and spectral characteristics of satellite data are the main issues regarding their separation. The Copernicus Sentinel-1 and -2 missions recently provided good spatial coverage and at the same time comparably high pixel spacing starting with 10 m for modes available across the entire Arctic. The purpose of this study was to assess the capabilities of both, Sentinel-1 C-band Synthetic Aperture Radar (SAR) and the Sentinel-2 multispectral information for Arctic focused mapping. Settings differ across the Arctic (historic settlements versus industrial, locations on bedrock versus tundra landscapes) and reference data are scarce and inconsistent. The type of features and data scarcity demand specific classification approaches. The machine learning approaches Gradient Boosting Machines (GBM) and deep learning (DL)-based semantic segmentation have been tested. Records for the Alaskan North Slope, Western Greenland, and Svalbard in addition to high-resolution satellite data have been used for validation and calibration. Deep learning is superior to GBM with respect to users accuracy. GBM therefore requires comprehensive postprocessing. SAR provides added value in case of GBM. VV is of benefit for road identification and HH for detection of buildings. Unfortunately, the Sentinel-1 acquisition strategy is varying across the Arctic. The majority is covered in VV+VH only. DL is of benefit for road and building detection but misses large proportions of other human-impacted areas, such as gravel pads which are typical for gas and oil fields. A combination of results from both GBM (Sentinel-1 and -2 combined) and DL (Sentinel-2; Sentinel-1 optional) is therefore suggested for circumpolar mapping.

Slope Hazard Monitoring Using High-Resolution Satellite Remote Sensing: Lessons Learned from a Case Study

In this study, a highway slope monitoring project for a section of US highway I-77 in Virginia was carried out with the InSAR technique. This paper attempts to provide insights into the complete and practical solution for the monitoring project, including two parts: what to consider for selecting the optimal satellites and configurations for the given area of interest (AoI) and budget; and how to best process the selected data for the monitoring purposes. To answer the first question, the simulated geometric distortion map, cumulative change detection map, intensity map, interferograms and coherence maps from all available historical datasets were generated. The satellite configuration that gives the best coherence and least geometric distortion with the given budget was selected for the monitoring project. For this project, it was the X-band COSMO stripmap with 3 m resolution and eight-days revisit time. To answer the second question, a multi-temporal InSAR (MTInSAR) was applied to retrieve the settlement time series of the slopes along the highway. Several special techniques were applied to increase the level of confidence, i.e., dividing AoI into smaller and independent areas, using a non-linear approach, etc. Finally, fieldwork was carried out for the interpretation and validation of the results. The AoI was overall stable, though some local changes were detected by the SAR signal which were validated by the fieldwork.

Seasonal Progression of Ground Displacement Identified with Satellite Radar Interferometry and the Impact of Unusually Warm Conditions on Permafrost at the Yamal Peninsula in 2016

Ground subsidence monitoring by Synthetic Aperture Radar interferometry (InSAR) over Arctic permafrost areas is largely limited by long revisit intervals, which can lead to signal decorrelation. Recent satellite missions such as COSMO-Skymed (X-band) and Sentinel-1 (C-band) have comparably short time intervals of a few days. We analyze dense records of COSMO-Skymed from 2013 and 2016 and of Sentinel-1 from 2016, 2017, and 2018 for the unfrozen period over central Yamal (Russia). These years were distinct in environmental conditions and 2016 in particular was unusually warm. We evaluate the InSAR-derived displacement with in situ subsidence records, active-layer thickness measurements, borehole temperature records, meteorological data, C-band scatterometer records, and a land-cover classification based on Sentinel-1 and -2 data. Our results indicate that a comparison of seasonal thaw evolution between years is feasible after accounting for the early thaw data gap in InSAR time series (as a result of snow cover) through an assessment with respect to degree-days of thawing. Average rates of subsidence agree between in situ and Sentinel-1 (corrected for viewing geometry), with 3.9 mm and 4.3 mm per 100 degree-days of thaw at the test site. X-band and C-band records agree well with each other, including seasonal evolution of subsidence. The average displacement is more than twice in magnitude at the active-layer monitoring test site in 2016 compared to the other years. We further demonstrate that InSAR displacement can not only provide information on the magnitude of ground thaw but also on soil properties through analyses of seasonal evolution in extreme years.