Radar Backscatter Research Articles

Estimation of lunar surface roughness using Chandrayaan-2 full-polarimetric DFSAR data

Polarimetric radar backscattering coefficients depend on the impinging frequency, dielectric constant, incidence angle, polarization, and surface roughness. Therefore, surface roughness can be estimated using high-resolution polarimetric radar datasets. For this purpose, we examine the utility of model-based and eigenvalue-based decomposition approaches in this work. While using model-based decomposition, the dielectric constant is estimated at the outset. Thereafter, the rms height is derived from a scalar multiplier that models surface roughness. We also propose a modified single-bounce eigenvalue relative difference (SERD) and establish that it is a better indicator of surface roughness than the circular polarization coherence and the original SERD for non-reflection-symmetric lunar terrain. The Chandrayaan-2 Dual Frequency Synthetic Aperture Radar (DFSAR) datasets acquired over three Apollo mission landing sites are used for demonstration.

Automated detection and tracking of medium-large icebergs from Sentinel-1 imagery using Google Earth Engine

Monitoring Antarctic icebergs helps us understand the interaction between ocean, atmosphere, and sea ice in the Southern Ocean. Although gigantic icebergs have been the objects of many remote sensing studies, medium icebergs in the Southern Ocean have been rarely monitored or traced. In this study, we develop an iceberg detection and tracking tool particularly for medium and large icebergs (0.4–10 km2), by using Python programming language and Sentinel-1 (S1) imagery, based on Google Earth Engine (GEE). To detect icebergs, we employ the simple non-iterative clustering (SNIC) and region adjacency graph (RAG) merging for object-based image segmentation and train/test the support vector machine (SVM) model with 6432 labeled segments of iceberg or non-icebergs from 40 images of S1 (2019–2021). Radar backscatter features and morphological features of those segments are used as the inputs of the SVM model. After icebergs are detected in two image scenes of different dates, we track the displacements of detected icebergs by comparing their 1-D shape signals. Our model shows ∼99% of accuracy in detecting icebergs and ∼ 93–98% of accuracy in tracking icebergs depending on the day difference between image scenes. When using our tool for the monitoring of icebergs in the Amundsen Sea, we find that the iceberg fraction varies from 2% to 8% in 2021 and most of icebergs move westward with a speed of <0.2 km/day.

The influence of two-dimensional currents on the sea surface radar backscattering coefficient

ABSTRACT In this paper, the current field modulation is described using a solution of the wave action equation. Using an integrated equation model (IEM), a current field modulated scattering model is established. The sea surface backscattering coefficient decreases with increasing current speed. At high wind speeds, the change of backscattering coefficient under 1 m ⋅ s−1 current field modulation is about 1 dB. An increase of 60° in the angle between wind and current direction leads to a decrease of 1.5 dB in the modulation coefficient. Finally, the model is verified using satellite data.

Tesserae: Surface differences across Venus’s “continents”

The heavily deformed upland tesserae are some of the most ancient geologic units on Venus and, as such, record the longest history of surface evolution. Our geologic understanding of these landforms is based largely on radar images from the Magellan mission, in which gross morphology and small-scale properties can be difficult to deconvolve. Here we use Magellan radar backscatter data for ridge slope surfaces in 22 highland areas to understand whether the tesserae can be subdivided in ways that differentiate surface property variations. Significant variations occur in the mean backscatter of ridge slopes, and we divide the tesserae into two groups with echoes lower (n = 15) or higher (n = 7) than an average tessera radar scattering behavior. While both few-kilometers-scale slopes and centimeter-scale roughness can modulate the radar returns, at least seven out of 15 tesserae with lower echoes are correlated with fine-grained impact crater ejecta deposits that smooth the surface. We propose that distal ejecta deposition plays a major role in creating the observed range of tessera radar properties and obscuring aspects of their original formation and in situ weathering. Our twofold classification system provides a new way of assessing the physical characteristics of tesserae from the Magellan data. Upcoming missions must consider both their original morphology and post-emplacement processes if we are to unlock the geologic record preserved in tesserae.

C and K band microwave penetration into snow on sea ice studied with off-the-shelf tank radars

AbstractSnow cover on sea ice poses a challenge for radar measurements as microwave penetration into snow is not yet fully understood. In this study, the aim is to investigate microwave penetration into snow on Arctic sea ice using commercial C (6 GHz) and K (26 GHz) band tank radars. Nadir-looking radar measurements collected at nine study locations over first-year and multiyear landfast sea ice in the Lincoln Sea in May 2018 are analysed together with detailed measurements of the physical properties of the snow cover to determine the dominant scattering horizons at both frequencies. They are evaluated for the feasibility to determine snow depth. The results show that in 39% of the measurements and only on first-year ice a major fraction of the C band radar backscatter originated closer to the snow–ice interface potentially enabling snow depth retrieval. At K band, 81% of the radar returns originated from the snow surface. Partly confirming the findings of previous studies, however, the analysis was potentially hampered by relatively warm air temperatures (up to$-0.9^\circ$C) during the study period as well as stratigraphic features and inconclusive microwave interaction with the saline basal layers found in the snow cover on first-year ice.

How Does Wind Influence Near-Nadir and Low-Incidence Ka-Band Radar Backscatter and Coherence from Small Inland Water Bodies?

While many studies have been conducted regarding wind-driven Ka-band scattering on the ocean and sea surfaces, few have identified the impacts of Ka-band scattering on small inland water bodies, and fewer have identified the influence of wind on coherence over water. These previous studies have been limited in spatial scale, covering only large water bodies &gt;25 km2. The recently launched Surface Water and Ocean Topography (SWOT) mission is the first Ka-band InSAR satellite designed for mapping water surface elevations and open water areas for rivers as narrow as 100 m and lakes as small as 0.0625 km2. Because measurements of these types are novel, there remains some uncertainty about expected backscatter amplitudes given wind-driven water surface roughness variability. A previous study using the airborne complement to SWOT, AirSWOT, found that low backscatter and low coherence values were indicative of higher errors in the water surface elevation products, recommending minimum thresholds for backscatter and coherence for filtering the data to increase the accuracy of averaged data for lakes and rivers. We determined that the global average wind speed over lakes is 4 m/s, and after comparing AirSWOT backscatter and coherence data with ERA-5 wind speeds, we found that the minimum required speed to retrieve high backscatter and coherence is 3 m/s. We examined 11,072 lakes across Canada and Alaska, with sizes ranging from 350 m2 to 156 km2, significantly smaller than what could be measured with previous Ka-band instruments in orbit. We found that small lakes (0.0625–0.25 km2) have significantly lower backscatter (3–5 dB) and 0.20–0.25 lower coherence than larger lakes (&gt;1 km2). These results suggest that approximately 75% of SWOT observable lake areas around the globe will have consistently high-accuracy water surface elevations, though seasonal wind variability should remain an important consideration. Despite very small lakes presenting lower average backscatter and coherence, this study asserts that SWOT will be able to accurately resolve the water surface elevations and water surface extents for significantly smaller water bodies than have been previously recorded from satellite altimeters. This study additionally lays the foundation for future high-resolution inland water wind speed studies using SWOT data, when the data become available, as the relationships between wind speed and Ka-band backscatter reflect those of traditional scatterometers designed for oceanic studies.

Global clustering of recent glacier surges from radar backscatter data, 2017–2022

AbstractUsing global Sentinel-1 radar backscatter data, we systematically map the locations of glaciers with surge-type activity during 2017–22. Patterns of pronounced increases or decreases in the strongest backscatter between two winter seasons often indicate large changes in glacier crevassing, which we treat here as a sign of surge-type activity. Validations against velocity time series, terminus advances and crevassing found in optical satellite images confirm the robustness of this approach. We find 115 surge-type events globally between 2017 and 2022, around 100 of which on glaciers already know as surge-type. Our data reveal a pronounced spatial clustering in three regions, (i) Karakoram, Pamirs and Western Kunlun Shan (~50 surges), (ii) Svalbard (~25) and (iii) Yukon/Alaska (~9), with only a few other scattered surges elsewhere. This spatial clustering is significantly more pronounced than the overall global clustering of known surge-type glaciers. The 2017–22 clustering may point to climatic forcing of surge initiation.

Grassland mowing event detection using combined optical, SAR, and weather time series

The European Union’s Common Agricultural Policy (CAP) and the Habitats Directive aim to improve biodiversity in agricultural landscapes. Both policies require enormous monitoring, which can be facilitated by remote sensing. Use intensity, measured by mowing frequency is an important indicator of biodiversity in permanent grasslands. The frequency and timing of mowing can be determined using satellite remote sensing because photosynthetically active biomass changes rapidly in response to mowing. However, the rapid regrowth of grasses requires very dense satellite time series for reliable detection. Radar time series can complement optical time series and fill in cloud-related gaps to overcome this problem. Additional weather data can support the detection of grassland mowing events, as mowing events are associated with specific meteorological conditions. However, previous studies have not fully exploited both potentials or different machine learning approaches for mowing event detection.This study presents a new transferable two-step approach to detect grassland mowing events using combined optical and SAR data and additional weather data. First, we filled cloud-related gaps in optical time series using a supervised machine learning regression with optical and SAR data. We then classified time series sequences of optical, SAR and weather data into mown and unmown using four different machine learning algorithms. We used time series of NDVI and EVI (combined Sentinel-2 and Landsat 8), SAR backscatter, six-day interferometric coherence, backscatter radar vegetation index, backscatter cross-ratio (Sentinel-1), and temperature and precipitation sums. Our test sites are distributed across Germany and cover the entire gradient of grassland use intensities.Mowing events could be detected with F1 values of up to 89%, first cut with up to 94%. Our results show no structural advantage of infilling time series with machine learning over linearly interpolated time series. The combined Sentinel-2 and Landsat-8 time series provided dense time series with mostly median gaps less than 20 days, which proved sufficient to reliably detect mowing events. SAR data were not essential for mowing event detection in our study, but weather data improved classification results for models trained on all areas and years. However, when the model was transferred to unknown years or areas that were not used for training, SAR data improved detection accuracy, whereas weather data degrade it. Models trained on all years but not all study sites detected mowing events with an accuracy of up to F1 = 76%. Models trained with all regions but not all years detected mowing events in untrained years with F1 up to 80%.

Quantifying Multifrequency Ocean Altimeter Wind Speed Error Due to Sea Surface Temperature and Resulting Impacts on Satellite Sea Level Measurements

Surface wind speed measurements from a satellite radar altimeter are used to adjust altimeter sea level measurements via sea state bias range correction. We focus here on previously neglected ocean radar backscatter and subsequent wind speed variations due to sea surface temperature (SST) change that may impact these sea level estimates. The expected error depends on the radar operating frequency and may be significant at the Ka band (36 GHz) frequency chosen for the new Surface Water and Ocean Topography (SWOT) satellite launched in December 2022. SWOT is expected to revolutionize oceanography by providing wide-swath Ka band observations and enhanced spatial resolution compared to conventional Ku band (14 GHz) altimetry. The change to the Ka band suggests a reconsideration of SST impact on wind and sea level estimates, and we investigate this in advance of SWOT using existing long-term Ku and Ka band satellite altimeter datasets. This study finds errors up to 1.5 m/s in wind speed estimation and 1.0 cm in sea level for AltiKa altimeter data. Future SWOT data analyses may require consideration of this dependence prior to using its radar backscatter data in its sea level estimation.

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.

94 GHz Radar Backscatter Characteristics of Alpine Glacier Ice

94 GHz Radar Backscatter Characteristics of Alpine Glacier Ice

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.

Wind redistribution of snow impacts the Ka- and Ku-band radar signatures of Arctic sea ice

Abstract. Wind-driven redistribution of snow on sea ice alters its topography and microstructure, yet the impact of these processes on radar signatures is poorly understood. Here, we examine the effects of snow redistribution over Arctic sea ice on radar waveforms and backscatter signatures obtained from a surface-based, fully polarimetric Ka- and Ku-band radar at incidence angles between 0∘ (nadir) and 50∘. Two wind events in November 2019 during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition are evaluated. During both events, changes in Ka- and Ku-band radar waveforms and backscatter coefficients at nadir are observed, coincident with surface topography changes measured by a terrestrial laser scanner. At both frequencies, redistribution caused snow densification at the surface and the uppermost layers, increasing the scattering at the air–snow interface at nadir and its prevalence as the dominant radar scattering surface. The waveform data also detected the presence of previous air–snow interfaces, buried beneath newly deposited snow. The additional scattering from previous air–snow interfaces could therefore affect the range retrieved from Ka- and Ku-band satellite altimeters. With increasing incidence angles, the relative scattering contribution of the air–snow interface decreases, and the snow–sea ice interface scattering increases. Relative to pre-wind event conditions, azimuthally averaged backscatter at nadir during the wind events increases by up to 8 dB (Ka-band) and 5 dB (Ku-band). Results show substantial backscatter variability within the scan area at all incidence angles and polarizations, in response to increasing wind speed and changes in wind direction. Our results show that snow redistribution and wind compaction need to be accounted for to interpret airborne and satellite radar measurements of snow-covered sea ice.

Investigating the Stability and Distribution of Surface Ice in Mercury’s Northernmost Craters

Observations made by Earth-based radar telescopes and the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft provided compelling evidence for water ice in Mercury's polar craters. In our investigation, we constructed higher-resolution (125 m pixel−1) digital elevation models (DEMs) for four of the largest northernmost craters, Kandinsky, Tolkien, Chesterton, and Tryggvadóttir. The DEMs were leveraged to model solar illumination and the thermal environment, products that were used to identify permanently shadowed regions and simulate surface temperatures. From these models, we predicted the regions of surface stability for ice and volatile organic compounds. These predictions were then compared against the Arecibo radar, Mercury Laser Altimeter (MLA), and Mercury Dual Imaging System (MDIS) data. Our radar analysis shows that areas of high radar backscatter are correlated with areas predicted to host surface ice. Additionally, we identify radar backscatter heterogeneities within the deposits that could be associated with variations in ice purity, mantling of the ice, or ice abundances. The MDIS analysis did not reveal conclusive evidence for ice or volatiles at the surface, while MLA results support the presence of water ice at the surface in these craters. However, evidence for boundaries between the surface ice and low-reflectance volatile organic compounds, as suggested could be present by our models, was inconclusive owing to the limited MESSENGER data in these regions. BepiColombo’s upcoming orbital mission at Mercury has the opportunity to obtain new measurements of these high-latitude craters and test our predictions for the distribution of surface volatiles in these environments.

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.

Convolutional Neural Network-Based Deep Learning Approach for Automatic Flood Mapping Using NovaSAR-1 and Sentinel-1 Data

The accuracy of most SAR-based flood classification and segmentation derived from semi-automated algorithms is often limited due to complicated radar backscatter. However, deep learning techniques, now widely applied in image classifications, have demonstrated excellent potential for mapping complex scenes and improving flood mapping accuracy. Therefore, this study aims to compare the image classification accuracy of three convolutional neural network (CNN)-based encoder–decoders (i.e., U-Net, PSPNet and DeepLapV3) by leveraging the end-to-end ArcGIS Pro workflow. A specific objective of this method consists of labelling and training each CNN model separately on publicly available dual-polarised pre-flood data (i.e., Sentinel-1 and NovaSAR-1) based on the ResNet convolutional backbone via a transfer learning approach. The neural network results were evaluated using multiple model training trials, validation loss, training loss and confusion matrix from test datasets. During testing on the post-flood data, the results revealed that U-Net marginally outperformed the other models. In this study, the overall accuracy and F1-score reached 99% and 98% on the test data, respectively. Interestingly, the segmentation results showed less use of manual cleaning, thus encouraging the use of open-source image data for the rapid, accurate and continuous monitoring of floods using the CNN-based approach.

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.

Effects of Directional Wave Spectra on the Modeling of Ocean Radar Backscatter at Various Azimuth Angles by a Modified Two-Scale Method

Knowledge of the ocean backscatter at various azimuth angles is critical to the radar detection of the ocean environment. In this study, the modified two-scale model (TSM), which introduces a correction term in the conventional TSM, is improved based on the empirical model, CMOD5.n. Then, the influences of different directional wave spectra on the prediction of azimuthal behavior of ocean radar backscatter are investigated by comparing the simulated results with CMOD5.n and the Advanced Scatterometer (ASCAT) measurements. The results show that the overall performance of the single spectra of D, A, E, and H18 and the composite spectra of AH18 and AEH18 in predicting ocean backscatter are different at different wind speeds and incidence angles. Generally, the AH18 spectrum has better performance at low and moderate wind speeds, while the A spectrum works better at high wind speed. Nevertheless, the wave spectra have little effect on the prediction of the azimuthal fluctuation of scattering, which is highly dependent on the directional spreading function. The relative patterns of azimuthal undulation produced by different spreading functions are rather different at different wind speeds, but similar under different incidence angles. The Gaussian spreading function generally has better performance in predicting the azimuthal fluctuation of scattering.

Towards a Consistent Wind Data Record for the CFOSAT Scatterometer

Since launch, the Ku-band rotating fan-beam scatterometer onboard the China–France Oceanography Satellite (CFOSAT) has provided valuable sea surface wind measurements for more than four years. The performance of CFOSAT scatterometer (CSCAT)-derived wind vectors is generally good in terms of root-mean-square error, while the absolute calibration error remains an issue in the current CSCAT product. In this paper, the temporal variation in CSCAT winds is overviewed by analyzing the collocated CSCAT and numerical weather prediction (NWP) model winds. Then, the reasons for the inconsistency of CSCAT-retrieved winds are discussed. The results show that the imperfect calibration of radar backscatter coefficients is likely the main problem of CSCAT wind processing. Consequently, a running-window-based (i.e., weekly) ocean calibration is proposed to evaluate the consistency of CSCAT radar backscatters, and in turn, to recalibrate CSCAT backscattering measurements before the reprocessing of CSCAT wind data. Although the proposed method is not feasible for the near-real-time processing of CSCAT data, it significantly mitigates the temporal variations in CSCAT wind speed bias, resulting in a more consistent CSCAT wind data record that may be beneficial to meteorological quantitative applications.

A global long-term, high-resolution satellite radar backscatter data record (1992–2022+): merging C-band ERS/ASCAT and Ku-band QSCAT

Abstract. Satellite radar backscatter contains unique information on land surface moisture, vegetation features, and surface roughness and has thus been used in a range of Earth science disciplines. However, there is no single global radar data set that has a relatively long wavelength and a decades-long time span. We here provide the first long-term (since 1992), high-resolution (∼8.9 km instead of the commonly used ∼25 km resolution) monthly satellite radar backscatter data set over global land areas, called the long-term, high-resolution scatterometer (LHScat) data set, by fusing signals from the European Remote Sensing satellite (ERS; 1992–2001; C-band; 5.3 GHz), Quick Scatterometer (QSCAT, 1999–2009; Ku-band; 13.4 GHz), and the Advanced SCATterometer (ASCAT; since 2007; C-band; 5.255 GHz). The 6-year data gap between C-band ERS and ASCAT was filled by modelling a substitute C-band signal during 1999–2009 from Ku-band QSCAT signals and climatic information. To this end, we first rescaled the signals from different sensors, pixel by pixel. We then corrected the monthly signal differences between the C-band and the scaled Ku-band signals by modelling the signal differences from climatic variables (i.e. monthly precipitation, skin temperature, and snow depth) using decision tree regression. The quality of the merged radar signal was assessed by computing the Pearson r, root mean square error (RMSE), and relative RMSE (rRMSE) between the C-band and the corrected Ku-band signals in the overlapping years (1999–2001 and 2007–2009). We obtained high Pearson r values and low RMSE values at both the regional (r≥0.92, RMSE ≤ 0.11 dB, and rRMSE ≤ 0.38) and pixel levels (median r across pixels ≥ 0.64, median RMSE ≤ 0.34 dB, and median rRMSE ≤ 0.88), suggesting high accuracy for the data-merging procedure. The merged radar signals were then validated against the European Space Agency (ESA) ERS-2 data, which provide observations for a subset of global pixels until 2011, even after the failure of on-board gyroscopes in 2001. We found highly concordant monthly dynamics between the merged radar signals and the ESA ERS-2 signals, with regional Pearson r values ranging from 0.79 to 0.98. These results showed that our merged radar data have a consistent C-band signal dynamic. The LHScat data set (https://doi.org/10.6084/m9.figshare.20407857; Tao et al., 2023) is expected to advance our understanding of the long-term changes in, e.g., global vegetation and soil moisture with a high spatial resolution. The data set will be updated on a regular basis to include the latest images acquired by ASCAT and to include even higher spatial and temporal resolutions.