Interferometric Wide Swath Mode Research Articles

The Characteristics of Submesoscale Eddies near the Coastal Regions of Eastern Japan: Insights from Sentinel-1 Imagery

A long-term time series of 319 Sentinel-1 SAR Imagery with Interferometric Wide Swath (IW) mode was used to study the characteristics of submesoscale eddies over Japanese coastal regions from 2015 to 2021, including spatiotemporal eddy properties and possible mechanisms of their formation. The results showed that around 98% of the 1499 eddies identified from the SAR snapshots were submesoscale eddies (horizontal scales of O1–20 km) with a ratio of around 78% cyclones to around 22% anticyclones. Around 8% of the submesoscale eddies were found in these SAR images in winter since the submesoscale current-induced signals are masked by the stronger wind speed, compared with other seasons. Typical features of submesoscale eddies are summarized, providing a preliminary qualitative analysis of potential generation mechanisms specific to the eddy characteristics in this region. This study suggests that Sentinel-1 images are capable of providing insights into the observed submesoscale eddies near the coastal regions of eastern Japan, thereby contributing to the improved understanding of the generation of submesoscale eddies.

Modelling and Mapping Oil Pollution across the Inland and Coastal Waters: A Case of Selected Segment of Gulf of Guinea

Artisanal fisheries in the Niger Delta and Nigeria’s coastal waters support the livelihood of its inhabitants, thus management of water pollution has public health significance for the regions and beyond. Oil Spills and contaminated runoff from the hinterland always find their way to the coast via the extensive network of rivers and creeks across the delta. This is a major cause of marine pollution in the region. With this understanding there is a need to map and model the dynamics of oil pollution across the coastal waters of the Niger Delta, thereby providing evidence and data support for their management. This study explored the potential of mapping oil slicks along the inland and coastal water of the Niger Delta using the SENTINEL-1 Synthetic Aperture Radar (SAR). The Interferometric Wide Swath Mode imagery Ground Range Detected Imagery for the 2nd of January 2019 was acquired for the study area, processed using Sentinel Application Platform (SNAP) and the oil spill detection algorithm was applied on the Amplitude Vertical-Vertical Band. After preprocessing, the non-polluted areas were found to have high backscatter (amplitude >18.0). For the impacted areas, lower values were recorded - amplitude decibel value of between 15.9 and 18.0. A total of 329 polygons of oil slick were identified, covering an area of 46.34 km2. The impacted surface area ranges between 0.000064 – 15 km2. The farthest impacted surface can be found around 18.25km from the Bonny coastline. Results showed that a significant proportion of the oil slicks are concentrated along the upper reaches of the Bonny River spanning Okrika and Ogu/Bolo LGA. Towards the nearshore, the oil slicks were disaggregated into several units. Contiguous units were found closer to the banks and coastline. Out of the 329 polygons, 14 were identified along the New Calabar River, emptying into the Bonny Estuary. The impacted surfaces along the New Calabar River were smaller in size compared to those found along the Bonny River. The results identified the pattern and distribution of oil slicks across a segment of the Gulf of Guinea dominated by oil and gas industry activities. This study highlighted the potential of remotely sensed data to support coastal pollution management efforts. It further highlighted the potential for spatiotemporal monitoring of oil slicks along coastal and inland waters. Challenges for deploying such a platform include the high computing power required and the lack of in-situ monitoring data for validation.

Mitigating the impact of dense vegetation on the Sentinel-1 surface soil moisture retrievals over Europe

ABSTRACT The C-band Synthetic Aperture Radar (SAR) on board of the Sentinel-1 satellites have a strong potential to retrieve Surface Soil Moisture (SSM). Using a change detection model to Sentinel-1 backscatter, an SSM product at a kilometre scale resolution over Europe could be established in the Copernicus Global Land Service (CGLS). Over areas with dense vegetation and high biomass. The geometry and water content influence the seasonality of the backscatter dynamics and hamper the SSM retrieval quality from Sentinel-1. This study demonstrates the effect of woody vegetation on SSM retrievals and proposes a masking method at the native resolution of Sentinel-1’s Interferometric Wide (IW) swath mode. At a continental 20 m grid, four dense vegetation masks are implemented over Europe in the resampling of the backscatter to a kilometre scale. The resulting backscatter is then used as input for the TUWien (TUW) change detection model and compared to both in-situ and modelled SSM. This paper highlights the potential of high-resolution vegetation datasets to mask for non-soil moisture-sensitive pixels at a sub-kilometre resolution. Results show that both correlation and seasonality of the retrieved SSM are improved by masking the dense vegetation at a 20 m resolution.

Assessment and mapping of landslides in steep mountainous terrain using PS-InSAR: A case study of Karimabad Valley in Chitral

Assessing and mapping landslides in rugged, mountainous terrain presents considerable challenges due to the high cost, risk, and complexity involved. High-resolution assessments of landslides are essential for mitigating landslide-relevant damage and understanding landslide kinematics. Recent studies have utilized satellite-based Synthetic Aperture Radar Interferometry (InSAR) to achieve millimeter-scale detection of ground surface deformation. Among these InSAR techniques, Persistent Scatterers Interferometry (PS-InSAR) is particularly useful in identifying landslide-prone areas, offering a means to map and detect slope movements and characterize the stability of potentially unstable slopes. This study focuses on the slow-moving unstable hillslopes of Karimabad Valley in Chitral. We conducted an analysis using Sentinel-1 SLC data over a three-year period, obtained in interferometric wide swath mode with vertical-vertical polarization. Our results, derived from both ascending and descending satellite passes, reveal that the hillslopes of Karimabad Valley exhibit an average velocity of −25 to +28 mm/yr. Based on the deformation pattern identified, we detected seven deformation zones, with particular concern for the Beshgram village area. To better understand terrain steepness, four topographic profile lines (AA', BB', CC', and DD') were extracted and correlated with the deformation time series. Furthermore, the deformation time series were examined in relation to rainfall, earthquake, and snowfall data. These analyses highlight the complex interplay of multiple factors contributing to the observed variations in deformation, including the region's topography characterized by steep slopes, the freezing and thawing of snow, and seismic events.

Elevation Change of CookE2 Subglacial Lake in East Antarctica Observed by DInSAR and Time-Segmented PSInSAR

In this study, elevation change and surface morphology of CookE2, one of the most active subglacial lakes in East Antarctica, were analyzed by using Differential Interferometric Synthetic Aperture Radar (DInSAR) and a newly adapted Time-Segmented Persistent Scatterer Interferometric Synthetic Aperture Radar (TS-PSInSAR) techniques. Firstly, several DInSAR pairs were used to study the surface morphology of the subglacial lake during the rapid discharge event in 2007 and the subsequent recharge in 2010 by using ALOS PALSAR data and the continuous recharge from 2018 to 2020 by using Sentinel-1 SAR data. For time-series observation from 2018 to 2020, however, simple integration of DInSAR deviates largely from the satellite altimeter data because errors from the horizontal flow of the surrounding ice field or atmospheric phase accumulate. Conventional PSInSAR deviates from the altimeter data if the LOS displacement exceeds 300 mm, i.e., approximately 1/4 of the slant range resolution of the Sentinel-1 SAR in Interferometric Wide-swath (IW) mode, during the time window. Therefore, a series of Time-Segmented PSInSAR with a 4-month time window could accurately distinguish 1.10 ± 0.01 m/year of highly linear (R2 = 0.99) surface rise rate of CookE2 and 0.63 m/year of horizontal deformation rate of the surrounding ice field from 2018 to 2020.

Multiparametric sea state fields from synthetic aperture radar for maritime situational awareness

This paper introduces a method for estimating a series of sea state parameters from satellite-borne synthetic aperture radar (SAR). The method was realized in a near real time (NRT) application which allows for the processing of data from different satellites and modes. The algorithm estimates the total significant wave height Hs, dominant and secondary swell and windsea wave heights, first, and second moment wave periods, the mean wave period and the period of wind sea. The algorithm was applied to the Sentinel-1 (S1) C-band Interferometric Wide Swath Mode (IW), Extra Wide (EW) and Wave Mode (WV) Level-1 (L1) products and also extended to X-band TerraSAR-X (TS-X) StripMap (SM) products. The scenes are processed in raster and result in continuous sea state fields, with the exception of S1 WV, where averaged values for sea state parameters for along-orbit imagettes of 20 km × 20 km are presented.The developed empirical algorithm consists of two parts: a first CWAVE_EX (extended CWAVE) part, based on a linear regression approach, and a subsequent machine learning part using the support vector machine (SVM) technique. A series of new data preparation steps (i.e. filtering, denoising) and new features estimated from SAR images are also introduced. The algorithm was tuned and validated using two independent global wave model hindcasts, WaveWatch-3 and MFWAM as well as National Data Buoy Center (NDBC) measurements. The achieved root mean squared errors (RMSE) for CWAVE_EX for the total Hs are 0.60 m for low-resolution modes S1 IW (10 m pixel) and EW (40 m pixel) and 0.35 m for S1 WV and TS-X SM (pixel spacing ca.1–4 m) in comparison to model predictions. The RMSEs of the retrieved wave periods are in the range of 0.45–0.90 s for all of the satellites and models considered. Similarly, the dominant and secondary swell, and wind sea wave height RMSEs are in the range of 0.35–0.80 m. The SVM postprocessing improves the accuracy of the initial results of CWAVE_EX for Hs and reaches an RMSE of 0.25 m for S1 WV. Comparisons to 64 NDBC buoys, collocated at distances shorter than 50 km to S1 WV imagettes worldwide, result in an RMSE of 0.41 m. All results and the methods presented are novel in terms of the accuracy achieved, combining the classical approach with machine learning techniques, and performing an automatic NRT processing of multiparametric sea state fields from L1 data with automatic switching for different satellites and modes.The complete archive of S1 WV L1 Single Look Complex products from December 2014 until February 2021 was processed to create a sea state parameter database and validated using model hindcast and buoy measurements. The derived parameters are available to the public within the scope of the European Space Agency's Climate Change Initiative.

High-Resolution Mapping of Winter Cereals in Europe by Time Series Landsat and Sentinel Images for 2016–2020

Winter cereals, including wheat, rye, barley, and triticale, are important food crops, and it is crucial to identify the distribution of winter cereals for monitoring crop growth and predicting yield. The production and plating area of winter cereals in Europe both contribute 12.57% to the total global cereal production and plating area in 2020. However, the distribution maps of winter cereals with high spatial resolution are scarce in Europe. Here, we first used synthetic aperture radar (SAR) data from Sentinel-1 A/B, in the Interferometric Wide (IW) swath mode, to distinguish rapeseed and winter cereals; we then used a time-weighted dynamic time warping (TWDTW) method to discriminate winter cereals from other crops by comparing the similarity of seasonal changes in the Normalized Difference Vegetation Index (NDVI) from Landsat and Sentinel-2 images. We generated winter cereal maps for 2016–2020 that cover 32 European countries with 30 m spatial resolution. Validation using field samples obtained from the Google Earth Engine (GEE) platform show that the producer’s and user’s accuracies are 91% ± 7.8% and 89% ± 10.3%, respectively, averaged over 32 countries in Europe. The winter cereal map agrees well with agricultural census data for planted winter cereal areas at municipal and country levels, with the averaged coefficient of determination R2 as 0.77 ± 0.15 for 2016–2019. In addition, our method can identify the distribution of winter cereals two months before harvest, with an overall accuracy of 88.4%, indicating that TWDTW is an effective method for timely crop growth monitoring and identification at the continent level. The winter cereal maps in Europe are available via an open-data repository.

Global seasonal Sentinel-1 interferometric coherence and backscatter data set

This data set is the first-of-its-kind spatial representation of multi-seasonal, global C-band Synthetic Aperture Radar (SAR) interferometric repeat-pass coherence and backscatter signatures. Coverage comprises land masses and ice sheets from 82° Northern to 79° Southern latitudes. The data set is derived from multi-temporal repeat-pass interferometric processing of about 205,000 Sentinel-1 C-band SAR images acquired in Interferometric Wide-Swath Mode from 1-Dec-2019 to 30-Nov-2020. The data set encompasses three sets of seasonal (December-February, March-May, June-August, September-November) metrics produced with a pixel spacing of three arcseconds: 1) Median 6-, 12-, 18-, 24-, 36-, and 48-days repeat-pass coherence at VV or HH polarizations, 2) Mean radiometrically terrain corrected backscatter (γ0) at VV and VH, or HH and HV polarizations, and 3) Estimated parameters of an exponential coherence decay model. The data set has been produced to obtain global, spatially detailed information on how decorrelation affects interferometric measurements of surface displacement and is rich in spatial and temporal information for a variety of mapping applications.

Sentinel-1 time series for mapping snow cover depletion and timing of snowmelt in Arctic periglacial environments: case study from Zackenberg and Kobbefjord, Greenland

Abstract. Snow cover (SC) and timing of snowmelt are key regulators of a wide range of Arctic ecosystem functions. Both are strongly influenced by the amplified Arctic warming and essential variables to understand environmental changes and their dynamics. This study evaluates the potential of Sentinel-1 (S-1) synthetic aperture radar (SAR) time series for monitoring SC depletion and snowmelt with high spatiotemporal resolution to capture their understudied small-scale heterogeneity. We use 97 dual-polarized S-1 SAR images acquired over northeastern Greenland and 94 over southwestern Greenland in the interferometric wide swath mode from the years 2017 and 2018. Comparison of S-1 intensity against SC fraction maps derived from orthorectified terrestrial time-lapse imagery indicates that SAR backscatter can increase before a decrease in SC fraction is observed. Hence, the increase in backscatter is related to changing snowpack properties during the runoff phase as well as decreasing SC fraction. We here present a novel empirical approach based on the temporal evolution of the SAR signal to identify start of runoff (SOR), end of snow cover (EOS) and SC extent for each S-1 observation date during melt using backscatter thresholds as well as the derivative. Comparison of SC with orthorectified time-lapse imagery indicates that HV polarization outperforms HH when using a global threshold. The derivative avoids manual selection of thresholds and adapts to different environmental settings and seasonal conditions. With a global configuration (threshold: 4 dB; polarization: HV) as well as with the derivative, the overall accuracy of SC maps was in all cases above 75 % and in more than half of cases above 90 %. Based on the physical principle of SAR backscatter during snowmelt, our approach is expected to work well in other low-vegetation areas and, hence, could support large-scale SC monitoring at high spatiotemporal resolution (20 m, 6 d) with high accuracy.

Characteristics of suspended sediment in Sentinel-1 synthetic aperture radar observations

ABSTRACT The main purpose of this study is to quantitatively examine the characteristics of suspended sediment on synthetic aperture radar (SAR) observations of sea roughness. In this work, seven Sentinel-1 (S-1) SAR images in 2019 acquired in the Interferometric Wide-swath (IW) mode at vertical-vertical (VV) polarization channel were collected, where Yangtze River freshwater dominants and suspended sediment concentration (SSC) are higher in continental shelf of East China Sea (ECS). These images were collocated with waves simulated from WAVEWATCH-III (WW3) and SSC inverted from the Haiyang-1C (HY-1C) Coastal Zone Imager (CZI) image. The comparison of WW3-simulated SWH yielded a 0.45-m root-mean-square error (RMSE) with a 0.78 correlation (COR) when comparing the simulations with the measurements from Haiyang-2B (HY-2B). The SAR-measured normalized radar cross-section (NRCS) is linearly related with the HY-1 C SSC at SSC > 1200 mg l–1, which is caused by the change of dielectric constant of sea water and viscosity. Moreover, the NRCS is simulated by the three-scale radar backscattering model. It is generally found that HY-1 C SSC enhances the NRCS in difference between simulations and observations up to 8 dB at SSC of 2000 mg l–1, which could distort the accuracy of SAR wind retrieval.

Automated Rain Detection by Dual-Polarization Sentinel-1 Data

Rain Signatures on C-band Synthetic Aperture Radar (SAR) images acquired over ocean are common and can dominate the backscattered signal from the ocean surface. In many cases, the inability to decipher between ocean and rain signatures can disturb the analysis of SAR scenes for maritime applications. This study relies on Sentinel-1 SAR acquisitions in the Interferometric Wide swath mode and high-resolution measurements from ground-based weather radar to document the rain impact on the radar backscattered signal in both co- and cross-polarization channels. The dark and bright rain signatures are found in connection with the timeliness of the rain cells. In particular, the bright patches are demonstrated by the hydrometeors (graupels, hails) in the melting layer. In general, the radar backscatter under rain increases with rain rate for a given sea state and decreases when the sea state strengthens. The rain also has a stronger impact on the radar signal in both polarizations when the incidence angle increases. The complementary sensitivity of the SAR signal of rain in both channels is then used to derive a filter to locate the areas in SAR scenes where the signal is not dominated by rain. The filter optimized to match the rain observed by the ground-based weather radar is more efficient when both polarization channels are considered. Case studies are presented to discuss the advantages and limitations of such a filtering approach.

HIGH-RESOLUTION URBAN MAPPING BY FUSION OF SAR AND OPTICAL DATA

Abstract. Mapping the exact extent of urban areas is a critical prerequisite in many remote sensing applications, such as hazard evaluation and change detection. The usage of Synthetical Aperture Radar (SAR) data has gained popularity due to the unique characteristics of the backscattered radio signal from human-made targets. The Sentinel-1 (S1) constellation, with a global revisit time of 6–12 days in Interferometric Wide Swath (IW) mode and free and open access to the data, allows the development of new applications to monitor urban sites. However, S1 is rarely considered when fine resolution is required due to the large pixel size and the need for spatial averaging to obtain robust estimators. We propose a method to improve Sentinel-1 urban classification performance by exploiting one Multi-Spectral (MS) image acquired by Sentinel-2 (S2). MS data is used for tracing the precise natural boundaries in a scene through superpixels segmentation. A machine learning approach is then applied to interpret the thematic context of each segment from short temporal stacks of coregistered SAR data. We use a short sensing period (around two months), so rapid changes can be traces. The proposed fusion of S1 and S2 data was tested in the area of Milan (Italy), with a total accuracy of about 90%. The ability to follow high-resolution details in a mixed environment is demonstrated, opening the possibility of efficiently tracing the human footprint.

Extreme Wind Speeds Retrieval Using Sentinel-1 IW Mode SAR Data

With the improvement in microwave radar technology, spaceborne synthetic aperture radar (SAR) is widely used to observe the tropical cyclone (TC) wind field. Based on European Space Agency Sentinel-1 Interferometric Wide swath (IW) mode imagery, this paper evaluates the correlation between vertical transmitting–horizontal receiving (VH) polarization signals and extreme ocean surface wind speeds (>40 m/s) under strong TC conditions. A geophysical model function (GMF) Sentinel-1 IW mode wind retrieval model after noise removal (S1IW.NR) was proposed, according to the SAR images of nine TCs and collocated stepped frequency microwave radiometer (SFMR) and soil moisture active passive (SMAP) radiometer wind speed measurements. Through curve fitting and regression correction, the new GMF exploits the relationships between VH-polarization normalized radar cross section, incident angle, and wind speed in each sub-swath and covers wind speeds up to 74 m/s. Based on collocated SAR and SFMR measurements of four TCs, the new GMF was validated in the wind speed range from 2 to 53 m/s. Results show that the correlation coefficient, bias, and root mean squared error were 0.89, −0.89 m/s, and 4.13 m/s, respectively, indicating that extreme winds can be retrieved accurately by the new model. In addition, we investigated the relationship between the S1IW.NR wind retrieval bias and the SFMR-measured rain rate. The S1IW.NR model tended to overestimate wind speeds under high rain rates.

Mapping Arctic Lake Ice Backscatter Anomalies Using Sentinel-1 Time Series on Google Earth Engine

Seepage of geological methane through sediments of Arctic lakes might contribute conceivably to the atmospheric methane budget. However, the abundance and precise locations of such seeps are poorly quantified. For Lake Neyto, one of the largest lakes on the Yamal Peninsula in Northwestern Siberia, temporally expanding regions of anomalously low backscatter in C-band SAR imagery acquired in late winter and spring have been suggested to be related to seepage of methane from hydrocarbon reservoirs. However, this hypothesis has not been verified using in-situ observations so far. Similar anomalies have also been identified for other lakes on Yamal, but it is still uncertain whether or how many of them are related to methane seepage. This study aimed to document similar lake ice backscatter anomalies on a regional scale over four study regions (the Yamal Peninsula and Tazovskiy Peninsulas; the Lena Delta in Russia; the National Petroleum Reserve Alaska) during different years using a time series based approach on Google Earth Engine (GEE) that quantifies changes of σ0 from the Sentinel-1 C-band SAR sensor over time. An algorithm for assessing the coverage that takes the number of acquisitions and maximum time between acquisitions into account is presented, and differences between the main operating modes of Sentinel-1 are evaluated. Results show that better coverage can be achieved in extra wide swath (EW) mode, but interferometric wide swath (IW) mode data could be useful for smaller study areas and to substantiate EW results. A classification of anomalies on Lake Neyto from EW Δσ0 images derived from GEE showed good agreement with the classification presented in a previous study. Automatic threshold-based per-lake counting of years where anomalies occurred was tested, but a number of issues related to this approach were identified. For example, effects of late grounding of the ice and anomalies potentially related to methane emissions could not be separated efficiently. Visualizations of Δσ0 images likely reflect the temporal expansions of anomalies and are expected to be particularly useful for identifying target areas for future field-based research. Characteristic anomalies that clearly resemble the ones observed for Lake Neyto could be identified solely visually in the Yamal and Tazovskiy study regions. All data and algorithms produced in the framework of this study are openly provided to the scientific community for future studies and might potentially aid our understanding of geological lake seepage upon the progression of related field-based studies and corresponding evaluations of formation hypotheses.

Toward a Comprehensive Dam Monitoring: On-Site and Remote-Retrieved Forcing Factors and Resulting Displacements (GNSS and PS–InSAR)

Many factors can influence the displacements of a dam, including water level variability and environmental temperatures, in addition to the dam composition. In this work, optical-based classification, thermal diachronic analysis, and a quasi-PS (Persistent Scatter) Interferometric SAR technique have been applied to determine both forcing factors and resulting displacements of the crest of the Castello dam (South Italy) over a one-year time period. The dataset includes Sentinel-1A images acquired in Interferometric Wide swath mode using the Terrain Observation with Progressive Scans SAR (TOPSAR); Landsat 8 Thermal Infrared Sensor (TIRS) thermal images, and Global Navigation Satellite System (GNSS) for interpreting the motion of the top of the dam retrieved via interferometry. Results suggest that it is possible to monitor both dam water level and temperature periodic forcing factors and resulting displacements via a synergistic use of different satellite images.

Soil Moisture Estimation using Sentinel-1 SAR Data and Land Surface Temperature in Panchmahal District, Gujarat State

This paper presents the potential for soil moisture (SM) retrieval using Sentinel-1 C-band Synthetic Aperture Radar (SAR) data acquired in Interferometric Wide Swath (IW) mode along with Land Surface Temperature (LST) estimated from analysis of LANDSAT-8 digital thermal data. In this study Sentinel-1 data acquired on 27 February 2020 was downloaded from Copernicus website and LANDSAT-8 OLI data acquired on 24 February 2020 from the website https://earthexplorer.usgs.gov/.The soil samples were collected from 70 test fields in different villages of three talukas for estimating soil moisture content using the gravimetric method. The Sentinel-1 SAR microwave data was analysed using open source tools of Sentinel Application Platform (SNAP) software for estimation of backscattering coefficient. Land surface temperature estimated using Landsat-8 thermal data. The Landsat-8, Thermal infrared sensor Band-10 data and operational land imager Band-4 and Band-5 data were used in estimating LST. The Soil Moisture Index (SMI) for all field test sites was computed using the LST values. The regression analysis using σ0VV and σ0VH polarization with soil moisture indicated that σ0VV polarization was more sensitive to soil moisture content as compared to σ0VH polarization. The multiple regression analysis using field measured soil moisture (MS %) as dependent variable, and σ0VV and SMI as independent variable was carried which resulted in the coefficient of determination (R2) of 0.788, 0.777 and 0.778 for Godhra, Goghamba and Kalol talukas, respectively. These linear regression equations were used to compute the predicted soil moisture in three talukas. The maps of spatial distribution of soil moisture in three talukas were generated using the respective regression equations of three talukas.

The potential of sentinel-1 InSAR coherence for grasslands monitoring in Eastern Cape, South Africa

This study evaluates the potential of the repeat-pass synthetic aperture radar interferometry (InSAR) coherence for grasslands monitoring in Eastern Cape Province, in South Africa. We used 12 pairs of InSAR data in the Interferometric Wide swath mode (IW), in VV-polarization (with different spatial baseline B⊥ and different temporal repeat intervals between the subsequent acquisitions), covering grasslands in the study area during two growing seasons; drought season (2015/2016) and non-drought season (2016/2017). Results show that non-drought season has higher variability in the coherence between the different regions of interest (ROIs) than the drought season and higher variability in the backscattering coefficient. Based on our results, we conclude that InSAR-based data have shown potential for grasslands monitoring, and the regularity of Sentinel-1 data is well suited for InSAR applications. However, the potential to separate grasslands from other vegetation categories appears to be limited with longer temporal baselines of 12 days intervals using only Sentinel-1A data. This potential will be improved with the 6 days intervals that will become available with Sentinel-1B in the study area.

Learning Convolutional Sparse Coding on Complex Domain for Interferometric Phase Restoration.

Interferometric phase restoration has been investigated for decades and most of the state-of-the-art methods have achieved promising performances for InSAR phase restoration. These methods generally follow the nonlocal filtering processing chain, aiming at circumventing the staircase effect and preserving the details of phase variations. In this article, we propose an alternative approach for InSAR phase restoration, that is, Complex Convolutional Sparse Coding (ComCSC) and its gradient regularized version. To the best of the authors' knowledge, this is the first time that we solve the InSAR phase restoration problem in a deconvolutional fashion. The proposed methods can not only suppress interferometric phase noise, but also avoid the staircase effect and preserve the details. Furthermore, they provide an insight into the elementary phase components for the interferometric phases. The experimental results on synthetic and realistic high- and medium-resolution data sets from TerraSAR-X StripMap and Sentinel-1 interferometric wide swath mode, respectively, show that our method outperforms those previous state-of-the-art methods based on nonlocal InSAR filters, particularly the state-of-the-art method: InSAR-BM3D. The source code of this article will be made publicly available for reproducible research inside the community.

Thermal Noise Removal From Polarimetric Sentinel-1 Data

This study proposes, for the first time, an approach to remove thermal noise from the wave coherency matrix, <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\mathrm {C_{2}}$ </tex-math></inline-formula> , estimated from single-look complex dual-polarization Interferometric Wide Swath mode Sentinel-1 synthetic aperture radar data. The approach is straightforward; it exploits the ThermalNoiseRemoval module, provided by the European Space Agency (ESA) in its Sentinel Application Platform (SNAP) software, to remove thermal noise from the channel intensities. Then, noise correction on the complex data is applied, in order to estimate the noise-free <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\mathrm {C_{2}}$ </tex-math></inline-formula> matrix. As a further novelty, the proposed approach can be implemented in SNAP, through the use of a processing graph that is here provided. The method is applied on a dense time series of Sentinel-1 data, collected on an agricultural area located near Seville, Spain. The impact of thermal noise on the estimation of the eigendecomposition parameters of <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\mathrm {C_{2}}$ </tex-math></inline-formula> , i.e., entropy ( <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$H_{2}$ </tex-math></inline-formula> ), average alpha angle ( <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\overline {\alpha _{2}}$ </tex-math></inline-formula> ), and anisotropy ( <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$A_{2}$ </tex-math></inline-formula> ), is assessed for different land-cover types, namely river, rice, forest, and urban areas. Monte Carlo simulations are implemented to assess the performance of the proposed approach in estimating <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$H_{2}$ </tex-math></inline-formula> , <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\overline {\alpha _{2}}$ </tex-math></inline-formula> , and <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$A_{2}$ </tex-math></inline-formula> . Results show that the proposed noise removal method improves the estimation of these parameters for the considered land-cover classes.

Deformations Prior to the Brumadinho Dam Collapse Revealed by Sentinel-1 InSAR Data Using SBAS and PSI Techniques

Differential Interferometric SAR (DInSAR) has been used to monitor surface deformations in open pit mines and tailings dams. In this paper, ground deformations have been detected on the area of tailings Dam-I at the Córrego do Feijão Mine (Brumadinho, Brazil) before its catastrophic failure occurred on 25 January 2019. Two techniques optimized for different scattering models, SBAS (Small BAseline Subset) and PSI (Persistent Scatterer Interferometry), were used to perform the analysis based on 26 Sentinel-1B images in Interferometric Wide Swath (IW) mode, which were acquired on descending orbits from 03 March 2018 to 22 January 2019. A WorldDEM Digital Surface Model (DSM) product was used to remove the topographic phase component. The results provided by both techniques showed a synoptic and informative view of the deformation process affecting the study area, with the detection of persistent trends of deformation on the crest, middle, and bottom sectors of the dam face until its collapse, as well as the settlements on the tailings. It is worth noting the detection of an acceleration in the displacement time-series for a short period near the failure. The maximum accumulated displacements detected along the downstream slope face were −39 mm (SBAS) and −48 mm (PSI). It is reasonable to consider that Sentinel-1 would provide decision makers with complementary motion information to the in situ monitoring system for risk assessment and for a better understanding of the ongoing instability phenomena affecting the tailings dam.