Synthetic Aperture Radar Backscatter Research Articles

The potential of optical and SAR time-series data for the improvement of aboveground biomass carbon estimation in Southwestern China’s evergreen coniferous forests

ABSTRACT Accurate assessments of forest biomass carbon are invaluable for managing forest resources, evaluating effects on ecological protection, and achieving goals related to climate change and sustainable development. Currently, the integration of optical and synthetic aperture radar (SAR) data has been extensively utilized in estimating forest aboveground biomass carbon (AGC), while it is limited by using single-phase remote sensing images. Time-series data, which capture the interannual dynamic growth and seasonal variations of photosynthetic phenology in forests, can sufficiently describe forest growth characteristics. However, there remains a gap in research focusing on utilizing satellite-based time-series data for AGC estimation, especially for SAR sensors. This study investigated the potential of satellite-based optical and SAR time-series data for estimating AGC. Here, we undertook nine quantitative experiments of AGC estimation from Landsat 8 and Sentinel-1 and tested several regression algorithms (including multiple linear regression (MLR), random forests (RF), artificial neural network (ANN), and extreme gradient boosting (XGBoost)) to explore the contributions of spatiotemporal features to AGC estimation. The results suggested that the XGBoost algorithm was suitable for AGC estimation with explanatory solid power and stable performance. The temporal features representing forest growth trends and periodic change characteristics (such as coefficients of continuous wavelet transform) were more valuable for AGC estimation than spatial features for both sensor types, accounting for around 40% ~50% of the variance compared to 17% ~25%. The combination of optical and SAR time-series data produced the best performance (R2 = 0.814, RMSE = 18.789 Mg C/ha, rRMSE = 26.235%), compared with when utilizing optical or SAR time-series data alone (optical: R2 of 0.657 and rRMSE of 35.317%; SAR: R2 of 0.672 and rRMSE of 34.701%). Feature importance analysis also verified that temporal features of optical vegetation indices, SWIR 1/2 bands, and SAR backscatter from VV polarization were the most critical variables for AGC estimation. Furthermore, incorporating temporal features into the modeling is illustrated to be effective in reducing saturation effects within high-biomass forests. This study demonstrated the superiority of time-series data for forest carbon estimation. While the applicability of this methodology has only been investigated in evergreen coniferous forests, it may provide a viable approach needed to make full use of increasingly better and free satellite time-series data to estimate forest AGC with high accuracy, supporting policy making of forest management and sustainable development.

Interacting Sentinel-2A, Sentinel 1A, and GF-2 Imagery to Improve the Accuracy of Forest Aboveground Biomass Estimation in a Dry-Hot Valley

Carbon absorption and storage in forests is one of the important ways to mitigate climate change. Therefore, it is essential to use a variety of remote-sensing resources to accurately estimate forest aboveground biomass (AGB) in dry-hot valley regions. In this study, satellite images from the Sentinel-1A, Sentinel-2A, and Gaofen-2 satellites were utilized to estimate the forest AGB in Yuanmou County, Yunnan Province, China. Different combinations of image data, based on selected variables of stepwise regression and their performance in constructing linear stepwise regression (LSR) and random forest (RF) models, were explored. The results showed that: (1) after adding the polarized values of the synthetic aperture radar backscatter coefficients, the combination fitting effect was significantly improved; (2) the fitting effect of the Sentinel-1A + Sentinel-2A + Gaofen-2 data combination was superior to the other combinations, indicating that the effective extraction of forest horizon and vertical information can improve the estimation effect of the forest AGB; and (3) the RF model exhibited superior fitting performance compared to the LSR model across all permutations of remotely sensed image datasets, with R2 values of 0.71 and 0.65, and RMSE values of 30.67 and 33.79 Mg/ha, respectively. These findings lay the groundwork for enhancing the precision of AGB estimation in dry-hot valley areas by integrating Sentinel-2A, Sentinel-1A, and GF-2 imagery, providing valuable insights for future research and applications.

Synthetic aperture radar backscatter is influenced by bubbles at the ice/water interface of an Antarctic lake

Synthetic aperture radar is an invaluable tool for monitoring lake ice. This study utilizes synthetic aperture radar to analyze the 2019-2023 backscatter time series of perennially ice-covered Lake Untersee in East Antarctica. We observed stark seasonal backscatter variations, averaging –9.6 dB from December to March and –3.7 dB from May to November. These fluctuations correspond to the abundance of sub-centimeter bubbles at the ice/water interface. Notably, the backscatter increase in April-May aligns closely with variations in ice thickness across the lake. Our findings suggest that ice cover thickness influences the timing and duration of ice accretion at the bottom, the accumulation of dissolved gases and bubbles, and the resultant changes in surface roughness at the ice/water interface. These factors collectively impact the backscatter response. This study enhances our understanding of the interactions between subsurface processes and synthetic aperture radar backscatter, shedding light on the seasonal dynamics of perennially ice-covered lakes.

CROP HEIGHT ESTIMATION OF WHEAT USING SENTINEL-1 SATELLITE IMAGERY: PRELIMINARY RESULTS

Abstract. Wheat is one of the primary crop productions for half of the global production. It is the most exported cereal, which reached up to 40% according to the records of the Food and Agriculture Organization of the United Nations (FAO) in 2020. The relationship between the crop parameters and remote sensing tools has become essential for the temporal monitoring and estimations of crop growth. In this study, we investigate the relationship between crop height and Synthetic Aperture Radar (SAR) backscatter. Plant height was collected twice, in the early (13 April 2023) and late stages (21 June 2023) of wheat, for a total of 70 samples. Sentinel-1 SAR data was also attempted to be synchronized with ground measurements. For this purpose, 8 images were obtained, 4 in both ascending (8 and 15 April, 16 and 26 June) and descending (9 and 14 April, 20 and 25 June) directions. The basic steps of image pre-processing, calibration, filtering, and topography correction, were applied to each image. By evaluating the correlation coefficients between plant height and images, data with low correlations were excluded and the first three data, dated 8 April and 19 June, were used for plant height estimation. For prediction purposes, Linear Regression (LR) and Random Forest (RF) methods were evaluated with three different training data sets. The highest correlation and minimum error value were achieved by VH polarization with random forest (r = 0.868, RMSE = 5.377 cm) in the early stage and LR (r = 0.616, RMSE = 14.451 cm) in the late stage. In general, higher training (80%) and lower test data (20%) produced better results.

The Biomass Proxy: Unlocking Global Agricultural Monitoring through Fusion of Sentinel-1 and Sentinel-2

The Biomass Proxy is a new cloud-free vegetation monitoring product that offers timely and analysis-ready data indicative of above-ground crop biomass dynamics at 10m spatial resolution. The Biomass Proxy links the consistent and continuous temporal signal of the Sentinel-1 Cross Ratio (CR), a vegetation index derived from Synthetic Aperture Radar backscatter, with the spatial information of the Sentinel-2 Normalized Difference Vegetation Index (NDVI), a vegetation index derived from optical observations. A global scaling relationship between CR and NDVI forms the basis of a novel fusion methodology based on static and dynamic combinations of temporal and spatial responses of CR and NDVI at field level. The fusion process is used to mitigate the impact on product quality of low satellite revisit periods due to acquisition design or persistent cloud coverage, and to respond to rapid changes in a timely manner to detect environmental and management events. The resulting Biomass Proxy provides time series that are continuous, unhindered by clouds, and produced uniformly across all geographical regions and crops. The Biomass Proxy offers opportunities including improved crop growth monitoring, event detection, and phenology stage detection.

Estimating vegetation water content from Sentinel-1 C-band SAR data over savanna and grassland ecosystems

Studying vegetation water content (VWC) dynamics is essential for understanding plant growth, water and carbon cycles, and ecosystem stability. However, acquiring field-based VWC estimates, consistently through space and time, is challenging due to time and resource constraints. This study investigates the potential of Sentinel-1 C-band Synthetic Aperture Radar (SAR) data for estimating VWC in natural ecosystems in central Brazil. We assessed (i) how well Sentinel-1 SAR data can capture variations in VWC over three different vegetation types (i.e. dry and waterlogged grasslands, and savannas) and (ii) how the studied vegetation types respond to seasonal dry periods in terms of water content. Field data from 82 plots, distributed across the three vegetation types and revisited in four different seasons, were used to calibrate and validate a model for VWC estimation. The calibrated model, with an R 2 of 0.52 and RMSE of 0.465 kg m−2, was then applied to Sentinel-1 SAR backscatter data to generate monthly VWC maps for grassland and savanna ecosystems at 30 m spatial resolution between April 2015 and September 2023. These maps, combined with rainfall and potential evapotranspiration data, provided insights into how the studied vegetation types respond to water shortage during the dry season at the community scale. More specifically, savannas showed to be better able to retain higher levels of water content during the dry season, probably due to a higher water holding capacity of the woody component together with its deep-root system ability to access deeper groundwater. This research demonstrates the potential of Sentinel-1 SAR data for monitoring VWC in natural ecosystems, allowing for future studies to assess ecosystems’ response to drought events and changes in their functioning, ultimately supporting land management decisions.

X- and Ku-Band SAR Backscattering Signatures of Snow-Covered Lake Ice and Sea Ice

In this work, backscattering signatures of snow-covered lake ice and sea ice from X- and Ku-band synthetic aperture radar (SAR) data are investigated. The SAR data were acquired with the ESA airborne SnowSAR sensor in winter 2012 over Lake Orajärvi in northern Finland and over landfast ice in the Bay of Bothnia of the Baltic Sea. Co-incident with the SnowSAR acquisitions, in situ snow and ice data were measured. In addition, time series of TerraSAR-X images and ice mass balance buoy data were acquired for Lake Orajärvi in 2011–2012. The main objective of our study was to investigate relationships between SAR backscattering signatures and snow depth over lake and sea ice, with the ultimate objective of assessing the feasibility of retrieval of snow characteristics using X- and Ku-band dual-polarization (VV and VH) SAR over freshwater or sea ice. This study constitutes the first comprehensive survey of snow backscattering signatures at these two combined frequencies over both lake and sea ice. For lake ice, we show that X-band VH-polarized backscattering coefficient (σo) and the Ku-band VV/VH-ratio exhibited the highest sensitivity to the snow depth. For sea ice, the highest sensitivity to the snow depth was found from the Ku-band VV-polarized σo and the Ku-band VV/VH-ratio. However, the observed relations were relatively weak, indicating that at least for the prevailing snow conditions, obtaining reliable estimates of snow depth over lake and sea ice would be challenging using only X- and Ku-band backscattering information.

Evaluation of limiting factors for SAR backscatter based cut detection of alpine grasslands

Several studies utilized C-band Synthetic Aperture Radar (SAR) backscatter time series to detect cut events of grasslands. They identified several potential factors hindering the detection: Vegetation characteristics, precipitation, and the timing of salvage of the harvested grass. This study uses a comprehensive in situ database to assess the impact of those factors on the detection rate of cut events by performing a cut detection based on Sentinel-1 backscatter time series and relating the accuracy to the potentially limiting factors. The results can be summarized in the following key findings: (i) The detection rate decreases significantly with grass heights below 35 cm and a biomass of less than 2100 kg/ha. As the grass of the first growth is typically characterized by greater height and higher biomass, first cuts achieved a higher accuracy with 85% compared to re-growth cuts with 65%. (ii) False positive cut events were related to higher precipitation amounts, but adding precipitation data to the model led only to a slight increase of the accuracy of re-growth cuts, but a decrease of the overall accuracy. (iii) No relation was found between the timing of salvage and the backscatter behaviour. These insights contribute to a better utilization of C-band backscatter for vegetation analysis and agricultural applications, including cut detection. Further research with dense in situ measurements, including Vegetation Water Content (VWC) is required to fully understand the behaviour of C-band backscatter over managed grasslands.

Better localized predictions with Out-of-Scope information and Explainable AI: One-Shot SAR backscatter nowcast framework with data from neighboring region

Synthetic Aperture Radar (SAR) provides 10-m weather-independent global Earth surface observations for various tasks such as land cover land use mapping, water body delineation, and vegetation change monitoring. However, the application of SAR imagery has been limited to retrospective tasks by a “first event then observation” rule. Recent studies have proven the feasibility of a one-shot nowcast of backscatters of SAR imagery using meteorological driving forces, soil moisture, geomorphic factors, and previous SAR images collected for the target area. Although the approach is promising, spatial connectivity, more specifically, the influence of the status of surrounding areas of the target location has yet to be considered. To fill that gap, this study proposed two nowcasting frameworks that can integrate precipitation and soil moisture data collected from surrounding areas through spatial aggregation (SA) and by processing spatial series (SS), respectively. The catastrophic 2019 Central US Flooding was used as a case study with the goal of predicting backscatters of SAR imagery captured during the event. The results from SA, SS, and a framework that only considers localized input (S0) are compared against each other as well as with the benchmark performance created with persistence assumption. Results show that S0, SA, and SS outperform the benchmark. More importantly, considering data from upstream neighboring areas further improves prediction accuracy compared to S0 where only inputs from the target region are considered. The gradients of results considering/not considering neighboring data indicates that neighboring data can alter the model’s attention on each feature of the localized input matrix. The difference in gradients between SA and SS indicates the way the neighboring information is integrated also matters. The methodology proposed by this study can serve as a building block for more active usage of SAR imagery in forward-looking tasks such as early flood warning and response.

Forest Aboveground Biomass Estimation in Subtropical Mountain Areas Based on Improved Water Cloud Model and PolSAR Decomposition Using L-Band PolSAR Data

Forest aboveground biomass (AGB) retrieval using synthetic aperture radar (SAR) backscatter has received extensive attention. The water cloud model (WCM), because of its simplicity and physical significance, has been one of the most commonly used models for estimating forest AGB using SAR backscatter. Nevertheless, forest AGB estimation using the WCM is usually based on simplified assumptions and empirical fitting, leading to results that tend to overestimate or underestimate. Moreover, the physical connection between the model and the polarimetric synthetic aperture radar (PolSAR) is not established, which leads to the limitation of the inversion scale. In this paper, based on the fully polarimetric SAR data from the Advanced Land Observing Satellite-2 (ALOS-2) Phased Array-type L-band Synthetic Aperture Radar (PALSAR-2), the relative contributions of the three major scattering mechanisms were first analyzed in a hilly area of southern China. On this basis, the traditional WCM was extended by considering the secondary scattering mechanism. Then, to establish the direct relationship between the vegetation scattering mechanism and forest AGB, a new relationship equation between the PolSAR decomposition model and the improved water cloud model (I-WCM) was constructed without the help of external data. Finally, a nonlinear iterative method was used to estimate the forest AGB. The results show that volume scattering is the dominant mechanism, accounting for more than 60%. Double-bounce scattering accounts for the smallest fraction, but still about 10%, which means that the contribution of the double-bounce scattering component is not negligible in forested areas because of the strong penetration capability of the long-wave SAR. The modified method provides a correlation coefficient R2 of 0.665 and a root mean square error (RMSE) of 21.902, which is an improvement of 36.42% compared to the traditional fitting method. Moreover, it enables the extraction of forest parameters at the pix scale using PolSAR data without the need for low-resolution external data and is thus helpful for high-resolution mapping of forest AGB.

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.