Synthetic Aperture Radar Backscatter Research Articles

Mapping high marsh and salt pannes/flats along the northern Gulf of Mexico coast

Coastal wetlands are predicted to undergo extensive transformation due to climate and land use change. Baseline maps of coastal wetlands can be used to help assess changes. Found in the upper portion of the estuarine zone, high marsh and salt pannes/flats provide ecosystem goods and services and are particularly important to fish and wildlife. We developed the first map of high marsh and salt pannes/flats along the northern Gulf of Mexico using regional models that included spectral indices related to greenness and wetness from optical satellite imagery, elevation data, irregularly flooded wetland probability information, and synthetic aperture radar backscatter. We found the greatest relative coverage of high marsh along the Texas coast (30% to 65%) and the Florida Panhandle (40%), whereas the greatest relative coverage of salt pannes/flats was along the lower Texas coast (74%) and the middle Texas coast (15%). As part of this effort, we also developed a map that highlighted irregularly flooded wetlands dominated by Juncus roemerianus (black needlerush) for part of the study area. Both maps had an overall accuracy of around 80%. Our results advance the understanding of estuarine marsh zonation and provide a baseline for assessing future transformations.

Synthetic aperture radar polarised backscattering behaviour in partially inundated agricultural fields

ABSTRACT A methodology for classifying rainfed paddy fields based on their hydrological conditions is lacking. This study analysed the behaviour of synthetic aperture radar (SAR) backscatter coefficients at each polarisation and cross-polarisation ratio index in a rainfed rice paddy field with a partially waterlogged surface. The SAR polarisations used were the VH and VV of the C-band SAR Sentinel-1 and the HV and HH of the L-band SAR Palsar-2. The relationship between backscatter coefficient and terrains including topographic categories and local relative elevation, which affect waterlogging conditions in rainfed paddy fields, were evaluated. The VV of C-band, and HH and HV of L-band showed different patterns in the time series variation according to the topographic categories. We observed that the combination of L-band HH and HV could be used to assess waterlogging conditions because they differ depending on variations in local elevation. The study results suggest that future studies can evaluate the microtopography and associated local hydrological environment using the combination of L-band HH and HV backscatter coefficient.

Mapping Antarctic crevasses and their evolution with deep learning applied to satellite radar imagery

Abstract. The fracturing of glaciers and ice shelves in Antarctica influences their dynamics and stability. Hence, data on the evolving distribution of crevasses are required to better understand the evolution of the ice sheet, though such data have traditionally been difficult and time-consuming to generate. Here, we present an automated method of mapping crevasses on grounded and floating ice with the application of convolutional neural networks to Sentinel-1 synthetic aperture radar backscatter data. We apply this method across Antarctica to images acquired between 2015 and 2022, producing a 7.5-year record of composite fracture maps at monthly intervals and 50 m spatial resolution and showing the distribution of crevasses around the majority of the ice sheet margin. We develop a method of quantifying changes to the density of ice shelf fractures using a time series of crevasse maps and show increases in crevassing on Thwaites and Pine Island ice shelves over the observational period, with observed changes elsewhere in the Amundsen Sea dominated by the advection of existing crevasses. Using stress fields computed using the BISICLES ice sheet model, we show that much of this structural change has occurred in buttressing regions of these ice shelves, indicating a recent and ongoing link between fracturing and the developing dynamics of the Amundsen Sea sector.

Impact of SAR-based vegetation attributes on the SMAP high-resolution soil moisture product

NASA's Soil Moisture Active Passive (SMAP) mission employs L-band radiometer observations with a spatial resolution of ~33 km to produce global soil moisture products at various spatial resolutions and regular intervals. Among these products is the SMAP-Sentinel active-passive high-resolution soil moisture, available at spatial resolutions of 1 km and 3 km. The SMAP-Sentinel product combines data from the SMAP L-band radiometer and the Copernicus Sentinel-1A and Sentinel-1B C-band Synthetic Aperture Radar (SAR) observations. The SMAP-Sentinel high-resolution soil moisture product enables various hydrological and agricultural applications that were previously beyond the capability of coarser resolution radiometer-only products. However, the assessment of SMAP-Sentinel product over SMAP core validation sites (CVS) has revealed high errors and biases, particularly in agricultural regions. We suspect that the discrepancy may be related to the vegetation attribute parameters (e.g., vegetation optical depth, τ) used in the SMAP-Sentinel baseline Single-Channel Algorithm (SCA) microwave emission model (i.e., τ-ω model) for soil moisture retrievals, which might be out of sync due to reliance on NDVI climatology in the τ estimation. In this study, we leveraged and evaluated the advantage of using vegetation attribute information inherently available through the SMAP ratiometer and Sentinel- 1A/1B SAR backscatter as an alternative to NDVI climatology-derived τ in the SCA to enhance soil moisture retrievals. The SAR cross-polarization (vh) backscatter potentially contains vegetation attribute information that can be correlated with the coarse resolution τ derived from SMAP observations. Therefore, we hypothesized that including Sentinel- 1A/1B vh observations to derive high-resolution (1 km and 3 km) τ, and subsequently incorporating it into the SCA, can significantly improve the accuracy of the SMAP-Sentinel soil moisture product. We established a statistical relationship between Sentinel- 1A/1B SAR backscatter and the SMAP Dual-Channel Algorithm (DCA)-based τ retrieval at ~33 km, for a specific region and different landcover classes at a global extent, using seven years of observations (2015–2022). Further, we developed an algorithm to downscale the SMAP DCA τ of ~33 km to 1 km and 3 km resolutions based on the developed statistical relationships, enabling us to capture the actual temporal variability in vegetation attributes (i.e., τ) at higher resolutions. The downscaled τ values were then incorporated into the SCA to retrieve SMAP-Sentinel high-resolution (1 km and 3 km) soil moisture. The validation analysis of high-resolution soil moisture retrievals demonstrates that the inclusion of the downscaled τ approach significantly improved the performance of the SMAP-Sentinel product by reducing both bias and error standard deviation compared to the retrievals based on τ derived from the NDVI climatology.

Characterising the ice sheet surface in Northeast Greenland using Sentinel-1 SAR data

AbstractOver half of the recent mass loss from the Greenland ice sheet, and its associated contribution to global sea level rise, can be attributed to increased surface meltwater runoff, with the remainder a result of dynamical processes such as calving and ice discharge. It is therefore important to quantify the distribution of melting on the ice sheet if we are to adequately understand past ice sheet change and make predictions for the future. In this article, we present a novel semi-empirical approach for characterising ice sheet surface conditions using high-resolution synthetic aperture radar (SAR) backscatter data from the Sentinel-1 satellite. We apply a state-space model to nine sites within North-East Greenland to identify changes in SAR backscatter, and we attribute these to different surface types with reference to optical satellite imagery and meteorological data. A set of decision-making rules for labelling ice sheet melting states are determined based on this analysis and subsequently applied to previously unseen sites. We show that our method performs well in (1) recognising some of the ice sheet surface types such as snow and dark ice and (2) determining whether the surface is melting or not melting. Sentinel-1 SAR data are of high spatial resolution; thus, in developing a method to identify the state of the surface from these data, we improve our capability to understand the variation of ice sheet melting across time and space.

Total-variation regularized U-Net for wildfire burned area mapping based on Sentinel-1 C-Band SAR backscattering data

Previous studies have shown that Synthetic Aperture Radar (SAR) is able to detect burned areas, serving as a key data source for monitoring active wildfires in situations where optical sensors are hindered by dense smoke or cloud cover. Radar remote sensing provides rich and useful data on historical burned areas, which are critical for large-scale wildfire burned area mapping. This study aims to unveil the potential inherent in Sentinel-1 C-Band SAR data and to investigate the impact of reference masks on wildfire burned area mapping. SAR images frequently exhibit disruptive speckle noise, which manifests as isolated pixels or small regions, thereby hindering the seamless identification of burned areas through SAR data. To mitigate noise and enhance the connectivity of SAR-based burned area delineation, we propose a novel approach: a Total-Variation (TV) regularized U-Net model, tailored to learn from noisy pseudo masks derived from Sentinel-1 C-Band SAR backscattering data. To validate its efficacy, we assembled a dataset encompassing 16 geographically diverse wildfire events, incorporating SAR-based kMap, a SAR-based binary label (designated as SARREF), and an optical-based binary label (referred to as OptREF). The kMap, serving as a relative indicator of SAR backscattering change, is defined as the divergence between post-fire SAR backscattering and the temporal average, normalized by the temporal standard deviation of pre-fire SAR time series. Utilizing SAR-based kMap as input and SARREF or OptREF as labels, we systematically explored the impact of TV regularization on SAR-based burned area mapping. Our experimental findings highlight several key insights: 1) Vanilla U-Net considerably elevates the F1 score from 0.61 to 0.66 when compared to threshold-based approaches; 2) The TV-regularized U-Net exhibits a notable capacity for enhancing both accuracy and the connectivity of SAR-based burned areas. When trained with noisy SARREF, it achieves a substantial F1 score of 0.69 and an IoU score of 0.53. 3) In scenarios involving OptREF, the vanilla U-Net attains peak performance with the highest F1 (0.75) and IoU (0.61) scores, with the expected outcome of TV regularization yielding no additional enhancement.

A novel semi-empirical model for crop leaf area index retrieval using SAR co- and cross-polarizations

The retrieval of continuous leaf area index (LAI) in space and time from remote sensing is beneficial for cropland monitoring and management. Synthetic Aperture Radar (SAR) with the advantages of all-weather operation and fine spatial resolutions has been utilized in various agricultural applications. Although the water cloud model (WCM) has been extensively used for LAI estimation over croplands, it is modified for application to crop canopies with large gaps. The requirement of prior information on soil moisture is also a hinderance for the model application. In this study, WCM is theoretically modified to consider the correlation of active microwave propagation through the canopy in downward and upward directions. Such a modification is particularly important for sparse vegetation with large gaps between crop plants in the early stage of the growing season. A parameter with explicit physical meaning, i.e. the full-vegetation backscattering coefficient, was defined to simplify our model scheme and make the model applicable for the whole growing season. In addition, a physics-based SAR data processing scheme is developed to entangle the influences of LAI and soil moisture on SAR backscatter by taking advantage of the multiple polarizations of RADARSAT-2 (R2) SAR data. In this way, LAI was estimated using modified WCM without the prior knowledge of soil moisture. To evaluate LAI retrieved from the R2 datasets (R2 LAI), ground-based LAI measurements were made at the experimental area of SMAPVEX16-MB in Canada with twenty approximately 800 m х 800 m plots in soybeans and corn. R2 LAI was well correlated to these ground-based LAI (n = 15, R2 = 0.63, RMSE = 0.99 m2·m−2 for soybeans; n = 5, R2 = 0.66, RMSE = 1.20 m2·m−2 for corn). LAI was also retrieved from optical data acquired by Sentinel-2/MSI (S2), denoted as S2 LAI. The R2 LAI and S2 LAI are well correlated and achieved coefficients of determination (R2) of 0.76 and 0.71 and root mean square errors (RMSE) of 1.1 and 1.4 m2·m−2 for soybeans and corn, respectively. The seasonal variations of R2 LAI and S2 LAI are generally similar except at the end of the growing season when S2 LAI is considerably larger than R2 LAI. S2 LAI followed the trend of the reduction in leaf chlorophyll content, while R2 LAI reduced only slightly due to the decrease in leaf water content near the end the growing season. R2 LAI represents the total standing leaf area useful for surface energy balance estimation, while S2 LAI responds to green leaf area useful for crop productivity modeling.

Sentinel-1 time series data for sunflower (Helianthus annuus) phenology monitoring

The potential of time series of Sentinel-1 (S1) Synthetic Aperture Radar (SAR) data for monitoring crops and their phenological stages has long been recognized. However, sunflower (Helianthus annuus) is one of the crops whose interaction with microwave signal is still not fully understood. Here, we aim to analyze and interpret time series of S1 data for sunflower phenology monitoring. The SAR backscattering coefficients in VH (vertical-horizontal) and VV (vertical-vertical) polarizations as well as their ratio (VH/VV) in both ascending and descending orbits were retrieved for three sites located in three different agro-ecological regions in Ukraine. We observed that sunflower backscattering response differs for the ascending and descending orbits for the VV polarization and VH/VV polarization ratio due to the directional behavior of the flower head. Our analysis did not find major change in VH polarization for the two orbits. We also compared the S1-based sunflower phenology with that of Normalized Difference Vegetation Index (NDVI) derived from optical data acquired by Sentinel-2 (S2). Our analysis revealed a considerable correlation r2 > 0.4 between S1 and NDVI for VH and VV polarizations in both ascending and descending orbits; however, we observed lower correlation (r2 < 0.20) for the VH/VV polarization ratio in descending orbit (VH/VVdes). Based on the characteristics of sunflower phenology, we identified two sunflower phenological stages: beginning of flowering (BoF) and end of flowering (EoF). The difference between S1-based dates using VH/VVdes and reference dates was 4.1 and 4.4 days for the BoF and EoF stages, respectively; the use of VH/VVdes+asc yielded weaker results: 7.4 and 6.1 days. These results suggest that the VH/VVdes is a preferable parameter for identifying sunflower phenological stages such as BoF and EoF.

A multi-scale algorithm for the NISAR mission high-resolution soil moisture product

This study proposes a multi-scale soil moisture algorithm for the upcoming NASA-ISRO SAR (NISAR) mission to estimate high-resolution (200 [m]) soil moisture (the water content of the soil). The algorithm takes advantage of the high-resolution (∼10 [m]) synthetic aperture radar (SAR) backscatter and coarse resolution modeled/reanalysis soil moisture products (∼ 9 [km]) to create a high-resolution (200 [m]) soil moisture product at a global extent. The end goal of the algorithm is to remove dependencies on any complex modeling, tedious retrieval steps, or multiple ancillary data needs, and subsequently decrease the degrees of freedom to achieve optimal accuracy in soil moisture retrievals. The use of modeled/reanalysis soil moisture products with high temporal resolution gives an added advantage in reducing the temporal mismatch between the two different inputs used in the algorithm. In this study, the proposed algorithm is tested using L-band UAVSAR backscatter (σ°) data and Advanced Land Observing Satellite −2 (ALOS-2) SAR σ° as a substitute for the NISAR L-band SAR observations. The algorithm uses the L-band SAR σ° to disaggregate coarse resolution (∼9 [km]) reanalysis soil moisture of the European Centre for Medium-Range Weather Forecast (ECMWF) to a high-resolution of ∼200 [m] soil moisture product. The potential of the algorithm is demonstrated over three sites in different hydroclimatic regions of the world, such as India, the USA, and Canada. The high-resolution soil moisture estimates were compared with the in-situ soil moisture measurements available for three sites (North India, Southern California, and Carman, Manitoba, Canada). In North India, in-situ measurements are from paddy crops with high vegetation water content, the unbiased root-mean-square-error (ubRMSE) for the high-resolution soil moisture retrievals was found to be 0.036 [m3/m3] with a bias of −0.051 [m3/m3]. For the southern California site, the validation statistics shows low ubRMSE of 0.027 [m3/m3] and a low bias of 0.016 [m3/m3]. At the Carman, Manitoba test site of Canada, where in-situ soil measurements were available for multiple crop types, the comparison shows that the ubRMSE for all the crop types lies below 0.05 [m3/m3] with an average bias of <0.07 [m3/m3]. The result confirms that the proposed algorithm meets the NISAR mission's accuracy goals, i.e., 0.06 [m3/m3] ubRMSE over areas with vegetation water content (VWC) below 5 [kg/m2]. Rigorous validation work will still need to be carried out in the future based on the availability of L-band SAR datasets and after the launch of the NISAR satellite.

Soil moisture retrieval from Sentinel-1 using a first-order radiative transfer model—A case-study over the Po-Valley

Soil moisture is an important variable controlling many land surface processes and is used to quantify precipitation, drought, flooding, irrigation and other factors that influence decision making and risk-assessment. This paper presents the retrieval of high resolution (∼1 km) soil moisture data from Sentinel-1 C-band Synthetic Aperture Radar (SAR) backscatter measurements using a new bistatic radiative transfer modeling framework (RT1) previously only tested for scatterometer data. The model is applied over a diverse set of landcover types across the entire Po-Valley in Italy over a 4-year time-period from 2016 to 2019. The performance of the soil moisture retrievals is analyzed with respect to the ERA5-Land reanalysis dataset. The model parameterisation and retrieval method are chosen such as to constitute a trade-off between a physically plausible and a computationally feasible modeling approach. The results demonstrate the potential of RT1 for the retrieval of high-resolution soil moisture data from SAR time series.

Characterizing and Mapping Volcanic Flow Deposits on Mount St. Helens via Dual-Band SAR Imagery

Mapping volcanic flow deposits can be achieved by considering backscattering characteristics as a metric of surface roughness. In this study, we developed an approach to extract a measure of surface roughness from dual-band airborne Synthetic Aperture Radar (ASAR) backscattering data to characterize and map various volcanic flow deposits—namely, debris avalanches, lahars, lava flows, and pyroclastic density currents. We employed ASAR and Indian Space Research Organization (ISRO) airborne SAR datasets, from a joint project (ASAR-ISRO), acquired in December 2019 at 2 m spatial resolution, to assess the role and importance of incorporating dual-band data, i.e., L-band and S-band, into surface roughness models. Additionally, we derived and analyzed surface roughness from a digital surface model (DSM) generated from unoccupied aircraft systems (UAS) acquisitions using Structure from Motion (SfM) photogrammetry techniques. These UAS-derived surface roughness outputs served as meter-scale calibration products to validate the radar roughness data over targeted areas. Herein, we applied our method to a region in the United States over the Mount St. Helens volcano in the Cascade Range of Washington state. Our results showed that dual-band systems can be utilized to characterize different types of volcanic deposits and range of terrain roughness. Importantly, we found that a combination of radar wavelengths (i.e., 9 and 24 cm), in tandem with high-spatial-resolution backscatter measurements, yields improved surface roughness maps, compared to single-band, satellite-based approaches at coarser resolution. The L-band (24 cm) can effectively differentiate small, medium, and large-scale structures, namely, blocks/boulders from fine-grained lahar deposits and hummocks from debris avalanche deposits. Additionally, variation in the roughness estimates of lahar and debris avalanche deposits can be identified and quantified individually. In contrast, the S-band (9 cm) can distinguish different soil moisture conditions across variable terrain; for example, identify wet active channels. In principle, this dual-band approach can also be employed with time series of various other SAR data of higher coherence (such as satellite SAR), using different wavelengths and polarizations, encompassing a wider range of surface roughness, and ultimately enabling additional applications at other volcanoes worldwide and even beyond volcanology.

Discrimination and monitoring of rice cultural types using dense time series of Sentinel-1 SAR data

Details about the spatial extent of rice cultivation in different seasons are essential for ensuring global food security as well as addressing environmental issues and climate change. This study aims to introduce a novel approach to discriminate between different rice cultural types and subtypes and identify their optimum stages to represent crop growth profiles using Synthetic Aperture Radar (SAR) data. Dense time series Sentinel-1 SAR backscatter data were analyzed at 12 days intervals throughout the growth period of summer and kharif rice. Statistical analysis using ANOVA at 5% level of significance was conducted to identify the critical stages of discrimination for all cultural types and subtypes. Knowledge-based decision tree algorithm was employed to classify all major types of rice grown during summer and kharif season. A comprehensive analysis revealed significant variations in backscatter profiles among different rice types. The most critical period for effectively distinguishing between various rice cultural types and subtypes was found to be from transplanting to tillering stage. Decision tree algorithm was able to discriminate between rice cultural types and subtypes with high overall accuracy and Kappa coefficient, viz., 94.74% and 0.94 for summer rice and 91.80% and 0.90 for kharif rice, respectively. The findings of this study revealed that C-band co and cross-polarized Sentinel-1 data have the potential to detect differences between different cultural types and subtypes of rice and, therefore, will be helpful in crop growth monitoring and agro-ecological change detection in the cultivation pattern.

Synergistic Integration of Time Series Optical and SAR Satellite Data for Mariculture Extraction

Mariculture is an important part of aquaculture, and it is important to address global food security and nutrition issues. However, seawater environmental conditions are complex and variable, which causes large uncertainties in the remote sensing spectral features. At the same time, mariculture types are distinct because of the different types of aquaculture (cage aquaculture and raft aquaculture). These factors bring great challenges for mariculture extraction and mapping using remote sensing. In order to solve these problems, an optical remote sensing aquaculture index named the marine aquaculture index (MAI) is proposed. Based on this spectral index, using time series Sentinel-1 and Sentinel-2 satellite data, a random forest classification scheme is proposed for mapping mariculture by combining spectral, textural, geometric, and synthetic aperture radar (SAR) backscattering. The results revealed that (1) MAI can emphasize the difference between mariculture and seawater; (2) the overall accuracy of mariculture in the Bohai Rim is 94.10%, and the kappa coefficient is 0.91; and (3) the area of cage aquaculture and raft aquaculture in the Bohai Rim is 16.89 km2 and 1206.71 km2, respectively. This study details an effective method for carrying out mariculture monitoring and ensuring the sustainable development of aquaculture.

SAR Based Sea Surface Complex Wind Fields Estimation: An Analysis over the Northern Adriatic Sea

Nowadays, sea surface analysis and monitoring increasingly use remote sensing, with particular interest in Synthetic Aperture Radar (SAR). Several SAR techniques exist in literature to understand the marine phenomena affecting the sea surface. In this work, we focus on the Doppler Centroid Anomaly (DCA), which accounts for the Doppler shift induced by sea surface movements. Starting from SAR raw data, we develop a processing chain to elaborate them and output the surface velocity map using DCA. The DCA technique has often been presented in the marine literature for estimating sea surface velocity, but more recently it has also been used to detect near-surface wind fields. This paper deals with estimating the sea surface wind field using Doppler information and SAR backscatter, combined with wind information provided by ECMWF and geophysical wind and Doppler model functions. We investigate the application of the approach in the coastal area of the northern Adriatic Sea (Northeast Italy). The test site is interesting, both for its particular orography, as it is a semi-enclosed basin largely surrounded by mountains, and for its complex meteorological phenomena, such as the Bora wind. Results obtained combining SAR backscatter and DCA information show an improvement in wind field estimation.

A Globally Applicable Method for NDVI Estimation from Sentinel-1 SAR Backscatter Using a Deep Neural Network and the SEN12TP Dataset

Vegetation monitoring is important for many applications, e.g., agriculture, food security, or forestry. Optical data from space-borne sensors and spectral indices derived from their data like the normalised difference vegetation index (NDVI) are frequently used in this context because of their simple derivation and interpretation. However, optical sensors have one major drawback: cloud coverage hinders data acquisition, which is especially troublesome for moderate and tropical regions. One solution to this problem is the use of cloud-penetrating synthetic aperture radar (SAR) sensors. Yet, with very different image characteristics of optical and SAR data, an optical sensor cannot be easily replaced by SAR sensors. This paper presents a globally applicable model for the estimation of NDVI values from Sentinel-1 C-band SAR backscatter data. First, the newly created dataset SEN12TP consisting of Sentinel-1 and -2 images is introduced. Its main features are the sophisticated global sampling strategy and that the images of the two sensors are time-paired. Using this dataset, a deep learning model is trained to regress SAR backscatter data to NDVI values. The benefit of auxiliary input information, e.g., digital elevation models, or land-cover maps is evaluated experimentally. After selection of the best model configuration, another experimental evaluation on a carefully selected hold-out test set confirms that high performance, low error, and good level of spatial detail are achieved. Finally, the potential of our approach to create dense NDVI time series of frequently clouded areas is shown. One limit of our approach is the neglect of the temporal characteristics of the SAR and NDVI data, since only data from a single date are used for prediction.

Deep Learning-Based Framework for Soil Moisture Content Retrieval of Bare Soil from Satellite Data

Machine learning (ML) is a branch of artificial intelligence (AI) that has been successfully applied in a variety of remote sensing applications, including geophysical information retrieval such as soil moisture content (SMC). Deep learning (DL) is a subfield of ML that uses models with complex structures to solve prediction problems with higher performance than traditional ML. In this study, a framework based on DL was developed for SMC retrieval. For this purpose, a sample dataset was built, which included synthetic aperture radar (SAR) backscattering, radar incidence angle, and ground truth data. Herein, the performance of five optimized ML prediction models was evaluated in terms of soil moisture prediction. However, to boost the prediction performance of these models, a DL-based data augmentation technique was implemented to create a reconstructed version of the available dataset. This includes building a sparse autoencoder DL network for data reconstruction. The Bayesian optimization strategy was employed for fine-tuning the hyperparameters of the ML models in order to improve their prediction performance. The results of our study highlighted the improved performance of the five ML prediction models with augmented data. The Gaussian process regression (GPR) showed the best prediction performance with 4.05% RMSE and 0.81 R2 on a 10% independent test subset.

Gaussian process regression-based forest above ground biomass retrieval from simulated L-band NISAR data

Accurate mapping of forest above ground biomass (AGB) is essential for understanding changes in the rate of ecosystem processes such as biomass productivity, litter productivity, actual litter decomposition, and potential litter decomposition during secondary succession. They also play a vital role in evaluating forest carbon pools. This study presents a Matérn kernel-based Gaussian process regression (GPR) approach for biomass estimation integrating Synthetic Aperture Radar (SAR) backscatter intensities with LiDAR measurements. We use the backscatter intensities from five scenes of L-band Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) simulated NISAR data collected during the 2019 AM–PM NASA airborne campaign. The biomass map derived from LVIS (Land, Vegetation, and Ice Sensor) Lidar point clouds over the Lenoir landing site, Alabama, was used with these simulated NISAR data. We have utilized the GPR model to estimate the AGB of the entire forest alongside the major forest classes in the study area, namely deciduous, evergreen, and woody wetlands. To examine the dependency of the model on acquisition conditions, we calibrated and validated our proposed model utilizing scene combinations. The experiment shows that multiple-scene retrieval delivered improved AGB estimates compared to single-scene retrieval. We evaluated the efficacy of the GPR model for three AGB ranges, i.e., (i) 8to100Mgha−1, (ii) 8to230Mgha−1, and (iii) 8to470Mgha−1 at 20m, 30m, 50m and 100m spatial resolution, respectively. The results indicate that the RMSE incurred by the GPR model for all these AGB ranges seemed to reduce as we increased spatial extent or reduced spatial resolution from 20m to 100m. Further, we demonstrate the performance of the optimized GPR model to retrieve AGB within the biomass range of 8Mgha−1 to 100Mgha−1 for all the above-mentioned forest types. Finally, this research highlights the advantages of the Matérn kernel-based GPR model over other statistical regression models towards improved AGB mapping.

Toward improved urban flood detection using Sentinel-1: dependence of the ratio of post- to preflood double scattering cross sections on building orientation

High-resolution synthetic aperture radar (SAR) sensors are now commonly used for flood detection. Automated detection tends to be limited to rural areas owing to the complicated backscattering mechanisms occurring in urban areas. Flooding can be identified in urban areas by searching for increased SAR backscatter in a postflood image due to double scattering between water and adjacent buildings, compared with a preflood image where double scattering is between unflooded ground and buildings. For co-polarized data, if φ is the angle between the building wall and the satellite direction of travel, double scattering is strongest for φ = 0 deg and falls off as φ increases. Theoretical studies estimating the ratio of flooded-to-unflooded double scatterer (DS) radar cross section (RCS) using X-band SAR data, found that the ratio was high at φ = 0 deg but only small at φ > 10 deg. Ostensibly, this implies that few DSs might be detected in an urban area. However, experiments on real images have called into question the veracity of the modeling. We describe an empirical study to examine the relationship between the flooded-to-unflooded DS RCS ratio and φ in Sentinel-1 (S-1) C-band data. We use high-resolution light detection and ranging and aerial photographs so that φ can be measured accurately and is based on S-1 images from flood events that occurred in the United Kingdom during the storms of winter 2019 to 2020. Results indicate that vertical–vertical polarization is better than vertical–horizontal at distinguishing flooded from unflooded DS; that the theoretical model used underestimates the number of DS with high RCS ratios in the φ range 10 deg to 30 deg; and that sufficient DS ground heights can be determined to estimate an accurate local average flood level, although in high density housing there are less of these due to the presence of adjacent buildings.

Rapid pre-explosion increase in dome extrusion rate at La Soufrière, St. Vincent quantified from synthetic aperture radar backscatter

The extrusion rate of a lava dome is a critical parameter for monitoring silicic eruptions and forecasting their development. Satellite radar backscatter can provide unique information about dome growth during a volcanic eruption when other datasets (e.g., optical, thermal, ground-based measurements, etc.) may be limited. Here, we present an approach for estimating volcanic topography from individual backscatter images. Using data from multiple SAR sensors we apply the method to the dome growth during the 2021 eruption at La Soufrière, St. Vincent. We measure an average extrusion rate of 1.8 m3s−1 between December 2020 and March 2021 before an acceleration in extrusion rate to 17.5 m3s−1 in the 2 days prior to the explosive eruption on 9 April 2021. We estimate a final dome volume of 19.4 million m3, extrapolated from the SAR sensors, with approximately 15% of the total extruded volume emplaced in the last 2 days. A possible explanation for the acceleration in extrusion rate could be the combined emptying of a conduit and reservoir of older material before the ascent of gas-rich magma in April 2021.

A Step-Wise Workflow for SAR Remote Sensing of Perennial Heaving Mound/Crater on the Yamal Peninsula, Western Siberia

Climate change in the Arctic region is more significant than in other parts of our planet. One of the manifestations of these changes is crater creation with blowouts of a gas, ice and frozen soil mixture. In this context, dynamics studies of long-term heaving mounds that turn into craters as a result are relevant. A workflow for detecting and assessing anomalous dynamics of heaving mounds in the Arctic regions is proposed. Areas with anomalous increase of ALOS-2 PALSAR-2 synthetic aperture radar (SAR) backscattering intensity are detected in the first stage. These increases take place due to sudden changes in local terrain slopes when the scattering surface (mound slope) turns toward the radar. Radar backscattering intensity also rises due to depolarization at newly formed frost cracks. Validation of the detected anomaly is carried out at the second stage through a comparison of multi-temporal digital elevation models obtained from bistatic radar interferometry TerraSAR-X/TanDEM-X data. At the final stage, the deformations are assessed within the detected areas using differential SAR interferometry (DInSAR) technique by ALOS-2 PALSAR-2 data. The magnitude of the heaving along the line of sight (LOS) was 22–24 cm in the period from January 2019 to January 2020. In general, effectiveness for detecting the perennial heaving mounds and the rate assessment of their increase were demonstrated in the suggested workflow.