NASA-ISRO Synthetic Aperture Radar Research Articles

Evaluating D-InSAR Performance to Detect Small Water Level Fluctuations in Two Small Lakes in Sweden

Abstract Monitoring lake water level fluctuations is critical for managing water resources, predicting the impacts of climatic change, and preserving ecosystem services lakes provide. However, traditional gauging stations are insufficient to monitor all lakes worldwide due to the large number of existing lakes, the challenges of installation and maintenance, and the remote locations of some. Although satellite altimetry is an alternative for measuring water levels, it cannot monitor small lakes effectively. This study evaluates the potential of Differential Interferometric Synthetic Aperture Radar (D-InSAR) for tracking minor water level changes in small lakes, a method more typically used in wetland studies. We investigate two Swedish lakes using Sentinel-1A and Sentinel-1B data from 2019, generating six-day interferograms and filtering out those with in-situ water level changes exceeding one phase cycle. Our results show that D-InSAR can detect small water level changes with Lin's correlations up to 0.63 and 0.89 and RMSE values of approximately 9 and 4 mm, respectively. These results evidence the potential of future L-band SAR missions with larger wavelengths, such as the NASA-ISRO SAR (NISAR) of the National Aeronautics and Space Administration (NASA) and the Indian Space Research Organisation (ISRO), to track water level changes in lakes and aid water tracking missions such as the SWOT (Surface Water and Ocean Topography).

Uncertainty estimates in the NISAR high-resolution soil moisture retrievals from multi-scale algorithm

It is important to know the amount of systematic and random uncertainties in any state variable to improve its geophysical application potential. The expected high-resolution (200 [m]) soil moisture product from the NASA-ISRO Synthetic Aperture Radar (NISAR) mission is no exception. Thus, knowing the quality of the soil moisture retrievals through the estimation of various error sources is imperative. The estimation error sources in soil moisture retrievals can be obtained by various methods. In situ measurements provide a reliable estimate of the uncertainty of soil moisture retrievals. However, in situ measurements are available only for limited locations, as they are typically very tedious and expensive to obtain. Thus, an analytical approach has been developed to obtain an estimate of the uncertainty in the soil moisture retrievals that vary in space and time across grid-cells. This uncertainty estimation is specifically developed for the multi-scale algorithm of the upcoming NISAR mission, which will provide soil moisture retrievals at 200 [m] resolution. The multi-scale algorithm for the NISAR mission disaggregates the coarser resolution soil moisture (∼9 [km]) to high-resolution (∼200 [m]) using NISAR L-band SAR measurements. However, uncertainty in high-resolution soil moisture retrievals might be introduced due to errors in input datasets (e.g., coarse resolution soil moisture, instrument error of SAR, etc.) and multi-scale algorithm parameters. Therefore, this study carried out a detailed sensitivity analysis of input datasets and algorithm parameters using the proposed approach. The sensitivity analysis shows that error in the input coarse resolution soil moisture is one of the primary drivers of uncertainty in the high-resolution soil moisture retrievals. The other portion of the uncertainty comes from errors in the algorithm parameters, and noise in SAR co-pol and cross-pol backscatter observations. Furthermore, the approach was tested on the UAVSAR L-band data time-series that had been simulated to closely match the expected characteristics of NISAR (e.g., spatial resolution and noise). The uncertainty estimates in UAVSAR-based high-resolution retrievals were compared with the SMAPVEX-12 in situ measurements. The uncertainties estimated for different crops were found to be close to the ubRMSE metric, which is also lower than the NISAR mission accuracy goal (0.06 [m3/m3]).

Improved soil moisture estimation and detection of irrigation signal by incorporating SMAP soil moisture into the Indian Land Data Assimilation System (ILDAS)

Land surface models have facilitated the estimation of soil moisture over a range of spatiotemporal scales. However, limitations in model parameterization and under-representation of anthropogenic processes restrict their ability to estimate local-scale soil moisture variability, especially over irrigated areas. Assimilation of satellite-based soil moisture retrievals into land surface models can be a viable approach to overcome these constraints, specially over highly irrigated countries such as India, where such applications are rare. Additionally, large-scale validation of modeled soil moisture has been limited over India till now due to lack of a representative station network. By assimilating Soil Moisture Active Passive (SMAP)-based estimates into the state-of-the-art Indian Land Data Assimilation System (ILDAS) and combining with a new soil moisture station network of more than 200 stations, this study demonstrates improved soil moisture estimations and capture of irrigation signals over the region. The Noah-MP land surface model is forced by multiple local and global meteorological datasets and Ensemble Kalman Filter (EnKF) is used for assimilation of soil moisture. Comparison of open-loop and data assimilated soil moisture against station soil moisture data shows relative spatial mean improvement of 0.0178 in correlation and 0.0029 m3/m3 in RMSE. Further statistical comparison with in-situ data has also shown better results over most of the stations, as evident from improved correlations and reduced unbiased RMSE after assimilation. Finally, the climatology of soil moisture over the different irrigation fractions reveals that data assimilated outputs over irrigated grid cells tend to have higher soil moisture during dry winter season, demonstrating the ability to capture irrigation signals. These findings quantify the value of data assimilation in improving soil moisture estimates and the ability to capture unmodeled processes such as irrigation, which lays the science groundwork for upcoming space missions such as NASA ISRO Synthetic Aperture Radar (NISAR).

Snow water equivalent retrieval over Idaho – Part 1: Using Sentinel-1 repeat-pass interferometry

Abstract. Snow water equivalent (SWE) is identified as the key element of the snowpack that impacts rivers' streamflow and water cycle. Both active and passive microwave remote sensing methods have been used to retrieve SWE, but there does not currently exist a SWE product that provides useful estimates in mountainous terrain. Active sensors provide higher-resolution observations, but the suitable radar frequencies and temporal repeat intervals have not been available until recently. Interferometric synthetic aperture radar (InSAR) has been shown to have the potential to estimate SWE change. In this study, we apply this technique to a long time series of 6 d temporal repeat Sentinel-1 C-band data from the 2020–2021 winter. The retrievals show statistically significant correlations both temporally and spatially with independent in situ measurements of SWE. The SWE change measurements vary between −5.3 and 9.4 cm over the entire time series and all the in situ stations. The Pearson correlation and RMSE between retrieved SWE change observations and in situ stations measurements are 0.8 and 0.93 cm, respectively. The total retrieved SWE in the entire 2020–2021 time series shows an SWE error of less than 2 cm for the nine in situ stations in the scene. Additionally, the retrieved SWE using Sentinel-1 data is well correlated with lidar snow depth data, with correlation of more than 0.47. Low temporal coherence is identified as the main reason for degrading the performance of SWE retrieval using InSAR data. We also show that the performance of the phase unwrapping algorithm degrades in regions with low temporal coherence. A higher frequency such as L-band improves the temporal coherence and SWE ambiguity. SWE retrieval using C-band Sentinel-1 data is shown to be successful, but faster revisit is required to avoid low temporal coherence. Global SWE retrieval using radar interferometry will have a great opportunity with the upcoming L-band 12 d repeat-pass NASA-ISRO Synthetic Aperture Radar (NISAR) data and the future 6 d repeat-pass Radar Observing System for Europe in L-band (ROSE-L) data.

NISAR S‐SAR Digital Beamforming Architecture and Implementation

AbstractISRO(Indian Space Research Organization) and JPL(Jet Propulsion Laboratory), NASA are jointly developing a state‐of‐the‐art L&S band SAR, planned to be launched through GSLV during 2023‐24. The conventional SAR systems has a limitation of providing higher spatial resolution as well as wider swath coverage simultaneously[6]. The NASA‐ISRO Synthetic Aperture Radar (NISAR) is based on a novel SAR concept, which will be utilized to image wide swath as high resolution of stripmap SAR, also called SweepSAR[1]. NISAR uses Digital Beam Forming on receive along with SweepSAR concept to get the advantage of lesser transmit power, immunity to range ambiguities with narrower receive beams, improved SNR and noise equivalent sigma naught as compared to conventional SAR system and directly radiating arrays with digital beamforming[5].In this paper, digital beamforming (DBF)[2] technique being used in NISAR S‐SAR receive chain is presented. The paper focuses on the DBF architecture, methodology, mathematical equations, practical implementation challenges and strategies used for mitigation of such challenges in order to achieve desired performance. Subsequently, DBF Hardware implementation algorithm and improvement in S‐SAR performance as a result of DBF technique is discussed.

Estimation of sea ice drift and concentration during melt season using C-band dual-polarimetric Sentinel-1 data

Characterizing sea ice changes over the Arctic can give proper insight into its extent in the current context of climate change. Sea ice motion's effect on ice area is influenced by seasonal conditions and the prevailing climate. However, the motion of thinner ice during the melt season can result in polynyas and cause significant ice loss that lasts the entire year. Sea ice drift and concentration during melt season are rarely discussed. The present study attempts to estimate sea ice drift during the melt season (July–August) using Differential Synthetic Aperture Radar Interferometry (DInSAR) with the usability of C-band Sentinel-1 data. The displacements/drifts were estimated using the phase shifts between two complex images with a temporal extent of 12 days. The study also proposed estimating sea ice concentration using polarimetric descriptors (power, purity, and randomness). The results revealed drifts of about 14–15 cm in the northeastern and southern parts followed by 45–50 cm in the middle parts and 20–25 cm along the coasts and open areas. Based on three polarimetric descriptors, thin ice and other ice fragments over open areas have shown less sea ice concentration as compared to the more concentrated pancake and pack ice with more compactness over the ice crystals. The study can be useful for computing drift in the melt season and carrying out spatiotemporal time-series analysis of seasonal sea ice drift/deformation studies. The algorithm and results can be used as a base study for the upcoming L- and S-band carrying NASA-ISRO Synthetic Aperture Radar (NISAR) mission.

Comparing NISAR (Using Sentinel-1), USDA/NASS CDL, and Ground Truth Crop/Non-Crop Areas in an Urban Agricultural Region.

A general limitation in assessing the accuracy of land cover mapping is the availability of ground truth data. At sites where ground truth is not available, potentially inaccurate proxy datasets are used for sub-field-scale resolution investigations at large spatial scales, i.e., in the Contiguous United States. The USDA/NASS Cropland Data Layer (CDL) is a popular agricultural land cover dataset due to its high accuracy (>80%), resolution (30 m), and inclusions of many land cover and crop types. However, because the CDL is derived from satellite imagery and has resulting uncertainties, comparisons to available in situ data are necessary for verifying classification performance. This study compares the cropland mapping accuracies (crop/non-crop) of an optical approach (CDL) and the radar-based crop area (CA) approach used for the upcoming NASA-ISRO Synthetic Aperture Radar (NISAR) L- and S-band mission but using Sentinel-1 C-band data. CDL and CA performance are compared to ground truth data that includes 54 agricultural production and research fields located at USDA's Beltsville Agricultural Research Center (BARC) in Maryland, USA. We also evaluate non-crop mapping accuracy using twenty-six built-up and thirteen forest sites at BARC. The results show that the CDL and CA have a good pixel-wise agreement with one another (87%). However, the CA is notably more accurate compared to ground truth data than the CDL. The 2017-2021 mean accuracies for the CDL and CA, respectively, are 77% and 96% for crop, 100% and 94% for built-up, and 100% and 100% for forest, yielding an overall accuracy of 86% for the CDL and 96% for CA. This difference mainly stems from the CDL under-detecting crop cover at BARC, especially in 2017 and 2018. We also note that annual accuracy levels varied less for the CA (91-98%) than for the CDL (79-93%). This study demonstrates that a computationally inexpensive radar-based cropland mapping approach can also give accurate results over complex landscapes with accuracies similar to or better than optical approaches.

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.

Strategies to Measure Soil Moisture Using Traditional Methods, Automated Sensors, Remote Sensing, and Machine Learning Techniques: Review, Bibliometric Analysis, Applications, Research Findings, and Future Directions

This review provides a detailed synthesis of various in-situ, remote sensing, and machine learning approaches to estimate soil moisture. Bibliometric analysis of the published literature on soil moisture shows that Time-Domain Reflectometry (TDR) is the most widely used in-situ instrument, while remote sensing is the most preferred application, and the random forest is the widely applied algorithm to simulate surface soil moisture. We have applied ten most widely used machine learning models on a publicly available dataset (in-situ soil moisture measurement and satellite images) to predict soil moisture and compared their results. We have briefly discussed the potential of using the upcoming NASA-ISRO Synthetic Aperture Radar (NISAR) mission images to estimate soil moisture. Finally, this review discusses the capabilities of physics-informed and automated machine learning (AutoML) models to predict surface soil moisture at higher spatial and temporal resolutions. This review will assist researchers in investigating the applications of soil moisture in the broad domain of earth sciences.

Comparison of Different Dielectric Models to Estimate Penetration Depth of L- and S-Band SAR Signals into the Ground Surface

We evaluate the penetration depth of synthetic aperture radar (SAR) signals into the ground surface at different frequencies. We applied dielectric models (Dobson empirical, Hallikainen, and Dobson semi-empirical) on the ground surface composed of different soil types (sandy, loamy, and clayey). These models result in different penetration depths for the same set of sensors and soil properties. The Dobson semi-empirical model is more sensitive to the soil properties, followed by the Hallikainen and Dobson empirical models. We used the Dobson semi-empirical model to study the penetration depth of the upcoming NASA-ISRO synthetic aperture radar (NISAR) mission operated at the L-band (1.25 GHz) and the S-band (3.22 GHz) into the ground. We observed that depending upon the soil types, the penetration depth of the SAR signals ranges between 0 to 10 cm for the S-band and 0 to 25 cm for the L-band.

Evaluating the Robustness of NISAR's Cropland Product to Time of Observation, Observing Mode, and Dithering

AbstractCropland mapping is important for monitoring agricultural practices, cropland distribution and for supporting food security programs. Radar remote sensing will likely provide a means of cropland mapping, which can be efficient, accurate, and globally applied. Algorithms used for spaceborne radar cropland mapping are still undergoing calibration and validation efforts, necessitating ground and aircraft‐based campaigns ahead of launch. The upcoming NASA ISRO SAR (NISAR) mission will map croplands at 1 ha spatial resolution using a coefficient of variation approach. NASA Jet Propulsion Laboratory recently released simulated NISAR data over cropland calibration sites including morning (AM), evening (PM), and dithered with and without gaps. These freely available simulated data allow for making accurate a priori estimates of how NISAR’s cropland science algorithm is impacted by factors such as orbit (i.e., the incidence angle, time of day), bandwidth (i.e., whether the 5, 20, or 40 MHz channels are used), and varying the system's pulse repetition frequency—a process known as dithering. This study investigates these science mission questions over agricultural regions in Mississippi and Arkansas, USA. Results show that compositing data from both orbits (AM + PM) yielded comparable classification accuracy to using only one of AM or PM data. The study also found that dithering only slightly decreased cropland mapping capability (&lt;0.3% accuracy reduction).

Additive Manufacturing of Compliant Mechanisms for Deployable Aerospace Structures

In the past 10 years, complex deployable structures have become common on CubeSats and large-scale spacecraft. As new missions are pursued, there is an increased need for more mass and volume efficient deployments. Over the same period, metal additive manufacturing (AM) has enabled new forms of spaceflight hardware. However, AM of compliant mechanisms has not been fully leveraged for deployable aerospace structures. The Surface Water Ocean Topography and NASA-ISRO Synthetic Aperture Radar mission missions launching in 2022 both utilize large deployable masts. Each mast deployment is driven by numerous spring and damper mechanisms. Because of volume constraints, the spring mechanisms designed utilize high aspect ratio rectangular cross section torsion springs that represent the state of the art of manufacturing. This extreme spring design resulted in manufacturing difficulties and hardware failures during ground mechanism testing. Upon re-examining the mechanism design, AM enables torque performance, mass, and complexity improvements. AM allows for torsion spring cross sections not otherwise possible with traditional spring manufacturing methods. Prototype springs of various cross sections were printed in maraging steel and tested. Results confirmed design analysis, and doubling of the spring constant was achieved when compared to the traditional springs. The use of AM also allows springs to be built monolithically with surrounding structure. Design, manufacturing, and test findings will be discussed along with future implications for deployable aerospace structures.

Establishment of SAR calibration site at Antarctica: Pre-NISAR calibration activity

Upcoming Synthetic Aperture Radar (SAR) missions with larger swath like NASA-ISRO Synthetic Aperture Radar (NISAR) (~250 km) requires a large, flat and homogeneous low-background site, free from human-made structures for point targets deployment for calibration. Finding out such a large, flat, homogeneous area devoid of perceived sources of radar clutter is a challenging task. In this regard, Antarctica is a potential site for calibration as it fulfills many of the criteria required for the ideal calibration site and suitable for setting up a corner reflector (CR) network for SAR calibration. This network will help in the calibration study and aid the studies related to Interferometric SAR (InSAR) viz. ice velocity estimation. As a part of this activity, in-house designed and developed CRs were installed near Indian research base stations Maitri and Bharati at Antarctica during the 38th Indian Scientific Expedition to Antarctica (ISEA) during 2018–2019. Corner reflectors were designed keeping in mind the harsh environment of Antarctica, such as high katabatic winds, blizzards, snowfall and low temperature. Temporal and seasonal analysis of radar backscatter data using available Indian Radar Imaging Satellite (RISAT-1) and European SENTINEL-1 SAR satellite over Maitri and Bharati was carried out to find out the suitable CR deployment site. Super-hydrophobic microwave transparent cover was designed and installed over CRs to protect snow accumulation. Hydrophobic radar absorbing materials were installed over central-mount to decrease the background noise due to it. Various tests such as structural analysis of total CR system, thermal analysis of materials and Radar Cross Section (RCS) characterization of CR were done before sending it to Antarctica. Integrated Wide Swath (IW) mode data of Sentinel-1 satellite was used for analyzing time-series response of CR in SAR image.

Land cover classification of spaceborne multifrequency SAR and optical multispectral data using machine learning

This study compares the utility of multifrequency SAR and Optical multispectral data for land-cover classification of Mumbai city and its nearby regions with a special focus on water body mapping. The L-band ALOS-2 PALSAR-2, X-band TerraSAR-X, C-band RISAT-1, and Sentinel-2 datasets have been used in this work. This work is done as a retrospective study for the dual-frequency L and S-band NASA-ISRO Synthetic Aperture Radar (NISAR) mission. The ALOS-2 PALSAR-2 data has been pre-processed before the implementation of machine learning algorithms for image segmentation. Multi-looking is performed on ALOS-2 PALSAR-2 data to generate square pixels of size 5.78 m and then target decomposition is applied to generate a false-color composite RGB image. While in the case of TerraSAR-X and RISAT-1 datasets, no multi-looking was performed and direct target decomposition was applied to generate false-color composite RGB images. Similarly, for the optical dataset that has a resolution of 10 m, a true color composite, and a false color composite RGB image are generated. For the comparative study between ALOS-2 PALSAR-2 and Sentinel-2 dataset, the RGB images are divided into smaller chunks of size 500*500 pixels each to create a training and testing dataset. Ten image patches were taken from the large dataset, out of which eight patches were used to train the machine learning models Random Forest (RF), K Nearest Neighbor (KNN), and Support Vector Machine (SVM), and two patches were kept for testing and validation purpose. For training the machine learning models, feature vectors are generated using the Gabor filter, Scharr filter, Gaussian filter, and Median filter. For patch 1, the mIOU for true-color composite based Optical image varies from 0.2323 to 0.2866 with the RF classifier performing the best and the mIOU for false-color composite based Optical image varies from 0.4130 to 0.4941 with the RF classifier performing the best while for ALOS-2 PALSAR-2 data, the mIOU varies from 0.4033 to 0.4663 with the RF classifier outperforming the KNN and the SVM classifiers. For patch 2, the mIOU for true-color composite based Optical data varies from 0.3451 to 0.4517 with KNN performing the best and the mIOU for false-color composite based Optical image varies from 0.5156 to 0.5832 with the RF classifier performing the best while for ALOS-2 PALSAR-2 data, the mIOU varies from 0.4600 to 0.5178 with the RF classifier outperforming the KNN and the SVM classifiers. The gap between the performance of ALOS-2 PALSAR-2 data and Sentinel-2 optical data is observed when the IOU of water class is compared, with IOUw for the true-color composite based optical image at a maximum of 0.2525 and for false-color composite based optical image at a maximum of 0.7366 while for ALOS-2 PALSAR-2 data a maximum IOUw of 0.7948 is achieved. The better performance of SAR data as compared to true-color composite based optical image data is due to the misclassification of ground and water classes into urban and forest in the case of the true-color composite based optical dataset which can be attributed to the high similarity between water and forest classes in the case of true-color composite based optical data whereas both these classes are easily separable in case of SAR data. This issue is however resolved by using the false color composite based optical image dataset for the classification task which performs slightly better than ALOS-2 PALSAR-2 data in the overall classification task. However, the SAR data works best in water body detection as notable from the high IOU for water class in the case of SAR data. In addition to the comparative analysis between Sentinel-2 optical and ALOS-2 PALSAR-2 data, land-cover classification has been performed on X-band TerraSAR-X and C-band RISAT-1 data on a single patch and it has been found that the RF classifier performs the best, recording the mIOU 0.5815 for X-band TerraSAR-X data, mIOU of 0.4031 for the C-band RISAT-1 data, and mIOU of 0.6153 for the L-band ALOS-2 data.

Processing Techniques for Nadir Echo Suppression in Staggered Synthetic Aperture Radar

Synthetic aperture radar (SAR) is a class of high-resolution imaging radar particularly suitable for satellite remote sensing with diverse applications, such as biomass and ice monitoring, generation of digital elevation models, and measuring of subsidence. Staggered SAR is a novel mode of operation under consideration for the next-generation SAR missions such as Tandem-L and NASA-ISRO SAR. It uses digital beamforming and a continuous variation of the pulse repetition interval (PRI) to achieve high azimuth resolution over a much wider, continuous swath than is traditionally possible. This PRI variation renders infeasible use of existing methods to avoid nadir echoes, which might impair the quality of staggered SAR images. This work proposes processing techniques that mitigate the impact of nadir echoes in staggered SAR through localization and thresholding-and-blanking of these echoes in range-compressed data and recovery of part of the underlying useful signal through interpolation. The performance of these processing techniques is evaluated through simulations using real TerraSAR-X data. The proposed technique can be implemented as an optional stage in the processing chain of future staggered SAR missions and leads to improved image quality at a reasonable additional computational cost.

A New Algorithm for Land-Cover Classification Using PolSAR and InSAR Data and Its Application to Surface Roughness Mapping Along the Gulf Coast

During a flooding event, the ability of the terrain to dissipate water flow energy depends on its land-cover type and the associated surface roughness. In this study, we developed a new land-cover classification algorithm using repeat-pass polarimetric synthetic aperture radar (PolSAR) and interferometric synthetic aperture radar (InSAR) data. Through a two-level hierarchical approach, we classified nine land-cover types with distinct surface roughness coefficients (Manning’s <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$n$ </tex-math></inline-formula> ). We demonstrated the performance of this algorithm using available L-band ALOS PALSAR scenes acquired between April 2007 and April 2011 over the Houston area. The radar-based surface roughness estimates show a good agreement with those independently derived from NOAA’s 22-class Coastal Change Analysis Program (C-CAP) 2010 land-cover classification data. Our algorithm is robust, and the randomly selected training sets only account for 0.3% of the total multilooked radar pixels (30-m spacing). Furthermore, we were able to accurately map surface roughness over the New Orleans area using available ALOS PALSAR scenes without selecting any new training sets. We note that NOAA’s C-CAP data are currently used for estimating surface roughness in the operational storm-surge models, and a new version is typically released every five to six years. With the launch of the L-band NASA-ISRO Synthetic Aperture Radar (NISAR) mission in the near future, our algorithm can be used to fill the temporal gaps of the existing C-CAP-based surface roughness database and improve the accuracy of near real-time hydrodynamic modeling.

An Area-Based Projection Algorithm for SAR Radiometric Terrain Correction and Geocoding

This article describes a projection algorithm between radar and map coordinates based on the representation of radar samples as area elements (AEs) rather than point elements. Each AE on the map grid (geographic grid) is associated with a number of radar grid samples that intersect completely or partially the AE. The association enables the geocoding (i.e., the map projection of radar imagery) with adaptive multilooking, accurately accounting for all radar samples contributing to the geocoded elements according to topography and radar geometry. By using averaging rather than interpolation, the proposed projection does not suffer from interpolation overfitting. The area-based geocoding also enables the generation of the geocoded polarimetric covariance matrix (GCOV) and geocoded synthetic aperture radar (SAR) interferograms with adaptive multilooking. Analogously, the slant-range projection of geocoded data is improved by projecting geographic grid pixels onto the radar grid according to their corresponding location based on the radar geometry without leaving gaps. This approach is used to reduce the computation time of previously published radiometric terrain correction (RTC) algorithms, performing 3.6–6.5 times faster over multilooked data and up to 26.3 times faster over single-look data. We demonstrate the strengths of the proposed area projection (AP) algorithm for RTC and geocoding using Uninhabited Aerial Vehicle SAR (UAVSAR), Sentinel-1B, and ALOS-2/PALSAR-2 data, and evaluate the results in the context of the upcoming NASA-ISRO SAR (NISAR) mission.

Airborne and Spaceborne Lidar Reveal Trends and Patterns of Functional Diversity in a Semi-Arid Ecosystem

Assessing functional diversity and its abiotic controls at continuous spatial scales are crucial to understanding changes in ecosystem processes and services. Semi-arid ecosystems cover large portions of the global terrestrial surface and provide carbon cycling, habitat, and biodiversity, among other important ecosystem processes and services. Yet, the spatial trends and patterns of functional diversity in semi-arid ecosystems and their abiotic controls are unclear. The objectives of this study are two-fold. We evaluated the spatial pattern of functional diversity as estimated from small footprint airborne lidar (ALS) with respect to abiotic controls and fire in a semi-arid ecosystem. Secondly, we used our results to understand the capabilities of large footprint spaceborne lidar (GEDI) for future applications to semi-arid ecosystems. Overall, our findings revealed that functional diversity in this ecosystem is mainly governed by elevation, soil, and water availability. In burned areas, the ALS data show a trend of functional recovery with time since fire. With 16 months of data (April 2019-August 2020), GEDI predicted functional traits showed a moderate correlation (r = 41–61%) with the ALS predicted traits except for the plant area index (PAI) (r = 11%) of low height vegetation (&amp;lt;5 m). We found that the number of GEDI footprints relative to the size of the fire-disturbed areas (=&amp;lt; 2 km2) limited the ability to estimate the full effects of fire disturbance. However, the consistency of diversity trends between ALS and GEDI across our study area demonstrates GEDI’s potential of capturing functional diversity in similar semi-arid ecosystems. The capability of spaceborne lidar to map trends and patterns of functional diversity in this semi-arid ecosystem demonstrates its exciting potential to identify critical biophysical and ecological shifts. Furthermore, opportunities to fuse GEDI with complementary spaceborne data such as ICESat-2 or the upcoming NASA-ISRO Synthetic Aperture Radar (NISAR), and fine scale airborne data will allow us to fill gaps across space and time. For the first time, we have the potential to monitor carbon cycle dynamics, habitats and biodiversity across the globe in semi-arid ecosystems at fine vertical scales.

Estimating Forest Stand Height in Savannakhet, Lao PDR Using InSAR and Backscatter Methods with L-Band SAR Data

Forest stand height (FSH), or average canopy height, serves as an important indicator for forest monitoring. The information provided about above-ground biomass for greenhouse gas emissions reporting and estimating carbon storage is relevant for reporting for Reducing Emissions from Deforestation and Forest Degradation (REDD+). A novel forest height estimation method utilizing a fusion of backscatter and Interferometric Synthetic Aperture Radar (InSAR) data from JAXA’s Advanced Land Observing Satellite Phased Array type L-band Synthetic Aperture Radar (ALOS PALSAR) is applied to a use case in Savannakhet, Lao. Compared with LiDAR, the estimated height from the fusion method had an RMSE of 4.90 m and an R2 of 0.26. These results are comparable to previous studies using SAR estimation techniques. Despite limitations of data quality and quantity, the Savannakhet, Lao use case demonstrates the applicability of these techniques utilizing L-band SAR data for estimating FSH in tropical forests and can be used as a springboard for use of L-band data from the future NASA-ISRO SAR (NISAR) mission.

Moving Toward L‐Band NASA‐ISRO SAR Mission (NISAR) Dense Time Series: Multipolarization Object‐Based Classification of Wetlands Using Two Machine Learning Algorithms

AbstractGiven the key role wetlands play in climate regulation and shoreline stabilization, identifying their spatial distribution is essential for the management, restoration, and protection of these invaluable ecosystems. The increasing availability of high spatial and temporal resolution optical and synthetic aperture radar (SAR) remote sensing data coupled with advanced machine learning techniques have provided an unprecedented opportunity for mapping complex wetlands’ ecosystems. A recent partnership between the National Aeronautics and Space Administration (NASA) and the Indian Space Research Organization (ISRO) resulted in the design of the NASA‐ISRO SAR (NISAR) mission. In this study, the capability of L‐band simulated NISAR data for wetland mapping in Yucatan Lake, Louisiana, is investigated using two object‐based machine learning approaches: Support vector machine (SVM) and random forest (RF). L‐band Unmanned Aerial Vehicle SAR (UAVSAR) data are exploited as a proxy for NISAR data. Specifically, we evaluated the synergistic use of different polarimetric features for efficient delineation of wetland types, extracting 84 polarimetric features from more than 10 polarimetric decompositions. High spatial resolution National Agriculture Imagery Program imagery is applied for image segmentation using the mean‐shift algorithm. Overall accuracies of 74.33% and 81.93% obtained by SVM and RF, respectively, demonstrate the great possibility of L‐band prototype NISAR data for wetland mapping and monitoring. In addition, variable importance analysis using the Gini index for RF classifier suggests that H/A/ALPHA, Freeman‐Durden, and Aghababaee features have the highest contribution to the overall accuracy.