Search

Extensive efforts have been made to observe the accumulation and melting of seasonal snow. However, making accurate observations of snow water equivalent (SWE) at global scales is challenging. Active radar systems show promise, provided the dielectric properties of the snowpack are accurately constrained. The dielectric constant (k) determines the velocity of a radar wave through snow, which is a critical component of time-of-flight radar techniques such as ground penetrating radar and interferometric synthetic aperture radar (InSAR). However, equations used to estimate k have been validated only for specific conditions with limited in situ validation for seasonal snow applications. The goal of this work was to further understand the dielectric permittivity of seasonal snow under both dry and wet conditions. We utilized extensive direct field observations of k, along with corresponding snow density and liquid water content (LWC) measurements. Data were collected in the Jemez Mountains, NM; Sandia Mountains, NM; Grand Mesa, CO; and Cameron Pass, CO from February 2020 to May 2021. We present empirical relationships based on 146 snow pits for dry snow conditions and 92 independent LWC observations in naturally melting snowpacks. Regression results had r2 values of 0.57 and 0.37 for dry and wet snow conditions, respectively. Our results in dry snow showed large differences between our in situ observations and commonly applied equations. We attribute these differences to assumptions in the shape of the snow grains that may not hold true for seasonal snow applications. Different assumptions, and thus different equations, may be necessary for varying snowpack conditions in different climates, suggesting that further testing is necessary. When considering wet snow, large differences were found between commonly applied equations and our in situ measurements. Many previous equations assume a background (dry snow) k that we found to be inaccurate, as previously stated, and is the primary driver of resulting uncertainty. Our results suggest large errors in SWE (10–15%) or LWC (0.05–0.07 volumetric LWC) estimates based on current equations. The work presented here could prove useful for making accurate observations of changes in SWE using future InSAR opportunities such as NISAR and ROSE-L.

Urbanization increases the amount of impervious surfaces, making accurate information on spatial and temporal expansion trends essential; the challenge is to develop a cost- and labor-effective technique that is compatible with the assessment of multiple geographical locations in developing countries. Several studies have identified the potential of remote sensing and multiple source information in impervious surface quantification. Therefore, this study aims to fuse datasets from the Sentinel 1 and 2 Satellites to map the impervious surfaces of nine Pakistani cities and estimate their growth rates from 2016 to 2020 utilizing the random forest algorithm. All bands in the optical and radar images were resampled to 10 m resolution, projected to same coordinate system and geometrically aligned to stack into a single product. The models were then trained, and classifications were validated with land cover samples from Google Earth’s high-resolution images. Overall accuracies of classified maps ranged from 85% to 98% with the resultant quantities showing a strong linear relationship (R-squared value of 0.998) with the Copernicus Global Land Services data. There was up to 9% increase in accuracy and up to 12 % increase in kappa coefficient from the fused data with respect to optical alone. A McNemar test confirmed the superiority of fused data. Finally, the cities had growth rates ranging from 0.5% to 2.5%, with an average of 1.8%. The information obtained can alert urban planners and environmentalists to assess impervious surface impacts in the cities.

This research focused on the quantitative assessment of the surface deformation velocities and rates and their natural and man-made controlling factors at Tengiz Oilfield in Kazakhstan using the Small Baseline Subset remote sensing technique followed by 3D and 2D decompositions and cosine corrections to derive vertical and horizontal movements from line-of-sight (LOS) measurements. In the present research we applied time-series of Sentinel-1 satellite images acquired during 2018–2020. All ground deformation derivatives showed the continuous subsidence at the Tengiz oilfield with increasing velocity. 3D and 2D decompositions of LOS measurements to vertical movement showed that the Tengiz Oil Field 2018–2020 continuously subsided with the maximum annual vertical deformation velocity around 70 mm. Based on the LOS measurements, the maximum annual subsiding velocity was observed to be 60 mm. Cosine corrections of LOS measurements to vertical movement, however, revealed a maximum annual vertical deformation velocity of 77 mm. The vertical deformation confirmed typical patterns of subsidence caused by oil extraction. Detected east-west and north-south horizontal movements at the Tengiz field clearly indicated that the study area crossed by seismic faults is affected by natural tectonic processes. The overall RMSE of 3D decomposed vertical deformation in relationship to LOS measurements and cosine corrections were in the range of 10–13 mm and 6–8 mm, correspondingly. The results of the present research will support operators of oil and gas fields and also other types of infrastructure to evaluate the actual differences of InSAR ground deformation measurements against the required standards and the precision of measurements depending on the operational needs, timeframes and availability of radar imagery.

In the Ferret Valley (NW Italy), anthropic activities coexist close to the Grandes Jorasses massif’s glaciological complex. In the past, break-off events have caused damage to people and infrastructure. These events concerned two specific sectors: the Montitaz Lobe (Planpincieux Glacier) and the Whymper Serac (Grandes Jorasses Glacier). Since 2010, permanent and discontinuous survey campaigns have been conducted to identify potential failure precursors, investigate the glacier instability processes, and explore different monitoring approaches. Most of the existing terrestrial apparatuses that measure the surface kinematics have been adopted in the Grandes Jorasses area. The monitoring sites in this specific area are characterized by severe weather, complex geometry, logistic difficulties, and rapid processes dynamics. Such exceptional conditions highlighted the limitations and potentialities of the adopted monitoring approaches, including robotic total station (RTS), GNSS receivers, digital image correlation applied to time-lapse imagery, and terrestrial radar interferometry (TRI). We examined the measurement uncertainty of each system and their monitoring performances. We discussed their principal limitations and possible use for warning purposes. In the Grandes Jorasses area, the use of a time-lapse camera appeared to be a versatile and cost-effective solution, which, however is not suitable for warning applications, as it does not guarantee data continuity. RTS and GNSS have warning potentialities, but the target installation and maintenance in remote environments remain challenging. TRI is the most effective monitoring system for early warning purposes in such harsh conditions, as it provides near-real-time measurements. However, radar equipment is very costly and requires extreme logistic effort. In this framework, we present data integration strategies to overcome the abovementioned limits and we demonstrate that these strategies are optimal solutions to obtain data continuity and robustness.

International audience; In the past decades, Bragg scattering has been considered to be an important scattering mechanism of microwave backscattering from sea surfaces. However, as reported in many recent literatures, non-Bragg scattering (which is often attributed to wave breaking) also makes a significant impact on radar scattering, especially for Horizontal–Horizontal (HH) polarized radar signals. To date, we know far less about non-Bragg scattering than Bragg scattering. Herein, this paper carries out an investigation on non-Bragg scattering and its effect on radar echoes at moderate incidence angles, both for oil-free and oil-covered sea surfaces. This paper firstly presents a systematic comparison of several sea spectra commonly used for the simulation of microwave scattering from sea surfaces. It is found that none of them perform well for the description of Bragg waves. Then, the “pure” Bragg wave spectra are inverted in the framework of the two-scale model (TSM) and geophysical model functions (GMFs). The normalized radar cross sections (NRCS) related to total scattering, non-Bragg scattering, and “pure” Bragg scattering in C, X, and Ku-bands are simulated under various conditions (i.e., incidence angles, wind speeds, and wind directions). Quantitative assessments of the relative contributions of non-Bragg scattering to total scattering are conducted. We also perform a survey on the non-Bragg scattering from the oil-covered sea surface. This article provides some new insights for a better understanding of the non-Bragg microwave scattering from rough sea surfaces at moderate incidence angles.

Strong variability of Ka-band radar backscattering from short wind waves on the surface of water covered with surfactant films in the presence of internal waves (IW) was studied in wave tank experiments. It has been demonstrated that modulation of Ka-band radar return due to IW strongly depends on the relationship between the phase velocity of IW and the velocity of drifting surfactant films. An effect of the strong increase in surfactant concentration was revealed in convergent zones, associated with IW orbital velocities in the presence of a “resonance” surface steady current, the velocity of which was close to the IW phase velocity. A phenomenological model of suppression and modulations in the spectrum of small-scale wind waves due to films and IW was elaborated. It has been shown that backscatter modulation could not be explained by the modulation of free (linear) millimeter-scale Bragg waves, but was associated with the modulation of bound (parasitic) capillary ripples generated by longer, cm–dm-scale waves—a “cascade” modulation mechanism. Theoretical analysis based on the developed model was found to be consistent with experiments. Field observations which qualitatively illustrated the effect of strong modulation of Ka-band radar backscatter due to IW in the presence of resonance drift of surfactant films are presented.

Ground-based weather radar data plays an essential role in monitoring severe convective weather. The detection of such weather systems in time is critical for saving people’s lives and property. However, the limited spatial coverage of radars over the ocean and mountainous regions greatly limits their effective application. In this study, we propose a novel framework of a deep learning-based model to retrieve the radar composite reflectivity factor (RCRF) maps from the Fengyun-4A new-generation geostationary satellite data. The suggested framework consists of three main processes, i.e., satellite and radar data preprocessing, the deep learning-based regression model for retrieving the RCRF maps, as well as the testing and validation of the model. In addition, three typical cases are also analyzed and studied, including a cluster of rapidly developing convective cells, a Northeast China cold vortex, and the Super Typhoon Haishen. Compared with the high-quality precipitation rate products from the integrated Multi-satellite Retrievals for Global Precipitation Measurement, it is found that the retrieved RCRF maps are in good agreement with the precipitation pattern. The statistical results show that retrieved RCRF maps have an R-square of 0.88-0.96, a mean absolute error of 0.3-0.6 dBZ, and a root-mean-square error of 1.2-2.4 dBZ.

Abandoned agricultural land (AAL) is a European problem and phenomenon when agricultural land is gradually overgrown with shrubs and forest. This wood biomass has not yet been systematically inventoried. The aim of this study was to experimentally prove and validate the concept of the satellite-based estimation of woody above-ground biomass (AGB) on AAL in the Western Carpathian region. The analysis is based on Sentinel-1 and -2 satellite data, supported by field research and airborne laser scanning. An improved AGB estimate was achieved using radar and optical multi-temporal data and polarimetric coherence by creating integrated predictive models by multiple regression. Abandonment is represented by two basic AAL classes identified according to overgrowth by shrub formations (AAL1) and tree formations (AAL2). First, an allometric model for AAL1 estimation was derived based on empirical material obtained from blackthorn stands. AAL2 biomass was quantified by different procedures related to (1) mature trees, (2) stumps and (3) young trees. Then, three satellite-based predictive mathematical models for AGB were developed. The best model reached R2 = 0.84 and RMSE = 41.2 t·ha−1 (35.1%), parametrized for an AGB range of 4 to 350 t·ha−1. In addition to 3214 hectares of forest land, we identified 992 hectares of shrub–tree formations on AAL with significantly lower wood AGB than on forest land and with simple shrub composition.

The effective normalized radar cross section (NRCS) of breaking waves, σwb, is empirically derived based on joint synchronized Ka-band radar and video records of the sea surface from a research tower. The σwb is a key parameter that, along with the breaker footprint fraction, Q, defines the contribution of non-polarized backscattering, NP =σwbQ, to the total sea surface NRCS. Combined with the right representation of the regular Bragg and specular backscattering components, the NP component is fundamental to model and interpret sea surface radar measurements. As the first step, the difference between NRCS values for breaking and non-breaking conditions is scaled with the optically-observed Q and compared with the geometric optics model of breaker backscattering. Optically-derived Q might not be optimal to represent the effect of breaking waves on the radar measurements. Alternatively, we rely on the breaking crest length that is firmly detected by the video technique and the empirically estimated breaker decay (inverse wavelength) scale in the direction of breaking wave propagation. A simplified model of breaker NRCS is then proposed using the geometric optics approach. This semi-analytical model parameterizes the along-wave breaker decay from reported breaker roughness spectra, obtained in laboratory experiments with mechanically-generated breakers. These proposed empirical breaker NRCS estimates agree satisfactorily with observations.

A melting layer detection algorithm is developed for the NCAR 94 GHz airborne cloud radar (HIAPER CloudRadar, HCR). The detection method is based on maxima in the linear depolarization ratio and a discontinuity in the radial velocity field. A melting layer field is added to the radar data, which provides detected, interpolated, and estimated altitudes of the melting layer and the altitude of the 0 °C isotherm detected in model temperature data. The icing level is defined as the lowest melting layer, and the cloud data are flagged as either above (cold) or below (warm) the icing level. Analysis of the detected melting layer shows that the offset between the 0 °C isotherm and the actual melting layer varies with cloud type: in heavy convection sampled in the tropics, the melting layer is found up to 500 m below the 0 °C isotherm, while in shallow clouds, the offset is much smaller or sometimes vanishes completely. A relationship between the offset and the particle fall speed both above and below the melting layer is established. Special phenomena, such as a lowering of the melting layer towards the center of storms or split melting layers, were observed.

Monitoring rice production is essential for securing food security against climate change threats, such as drought and flood events becoming more intense and frequent. The current practice to survey an area of rice production manually and in near real-time is expensive and involves a high workload for local statisticians. Remote sensing technology with satellite-based sensors has grown in popularity in recent decades as an alternative approach, reducing the cost and time required for spatial analysis over a wide area. However, cloud-free pixels of optical imagery are required to produce accurate outputs for agriculture applications. Thus, in this study, we propose an integration of optical (PROBA-V) and radar (Sentinel-1) imagery for temporal mapping of rice growth stages, including bare land, vegetative, reproductive, and ripening stages. We have built classification models for both sensors and combined them into 12-day periodical rice growth-stage maps from January 2017 to September 2018 at the sub-district level over Java Island, the top rice production area in Indonesia. The accuracy measurement was based on the test dataset and the predicted cross-correlated with monthly local statistics. The overall accuracy of the rice growth-stage model of PROBA-V was 83.87%, and the Sentinel-1 model was 71.74% with the Support Vector Machine classifier. The temporal maps were comparable with local statistics, with an average correlation between the vegetative area (remote sensing) and harvested area (local statistics) is 0.50, and lag time 89.5 days (n = 91). This result was similar to local statistics data, which correlate planting and the harvested area at 0.61, and the lag time as 90.4 days, respectively. Moreover, the cross-correlation between the predicted rice growth stage was also consistent with rice development in the area (r > 0.52, p < 0.01). This novel method is straightforward, easy to replicate and apply to other areas, and can be scaled up to the national and regional level to be used by stakeholders to support improved agricultural policies for sustainable rice production.

The paper investigates the penetration properties of an airborne Ku-band frequency modulated continuous waveform (FMCW) profiling radar named Tomoradar and a satellite near-infrared lidar into the boreal forest of Finland. We achieve the accumulative energy distributions based on the Tomoradar waveforms and the satellite lidar waveforms generated from the high-density airborne lidar data within Tomoradar footprints. By comparing two groups of the height percentiles and energy percentiles derived from the accumulative energy distributions, we evaluate the relationship of penetrations between the Ku-band microwave and near-infrared laser according to the coefficients of the determination (COD), and the root mean square errors (RMSE) of linear regression analyses. The quantitative analysis results demonstrate that the height and energy percentiles derived from Tomoradar waveforms correlate well with those from satellite lidar waveforms with the mean correlation coefficients of more than 0.78 and 0.85. The linear regression models for the height and energy percentile produce excellent fits with the mean CODs of 0.95 and 0.90 and the mean RMSEs of 1.25 m and 0.03, respectively. Less than 15% of height percentiles and 87.54% of the energy percentiles in the sixth stratum near the ground derived from Tomoradar waveforms surpass those from satellite lidar waveforms. Hence, the Ku-band microwave can penetrate deeper into the forest than the near-infrared laser at the same spatial scale. In addition, quadratic fitting models are established to describe the differences of the height percentile (DHP) and the energy percentile (DEP) to expound the canopy height and closure contributions numerically. The facts that the CODs of the DHP and DEP individually are more than 0.96 and 0.89 and the fitting residual histograms approximate to normal distributions reveal the reliabilities of the proposed fitting models. Thus, the penetration analyses are valid for the explorations on the FMCW radar applications and the data fusion of the Ku-band radar and near-infrared lidar in the forest investigations.

Ground Penetrating Radar (GPR) has proved to be a successful technique for the detection of landmines and Improvised Explosive Devices (IEDs) buried in the ground. In the last years, novel architectures for safe and fast detection, such as those based on GPR systems onboard Unmanned Aerial Vehicles (UAVs), have been proposed. Furthermore, improvements in GPR hardware and signal processing techniques have resulted in a more efficient detection. This contribution presents an experimental validation of a hybrid Forward-Looking–Down-Looking GPR architecture. The main goal of this architecture is to combine advantages of both GPR architectures: reduction of clutter coming from the ground surface in the case of Forward-Looking GPR (FLGPR), and greater dynamic range in the case of Down-Looking GPR (DLGPR). Compact radar modules working in the lower SHF frequency band have been used for the validation of the hybrid architecture, which involved realistic targets.

The retrieval of precipitation (snowfall and rainfall) from satellite sensors on a global basis is essential in aiding our knowledge and understanding of the Earth System and for many societal applications. Measurements from surface-based instruments are essentially limited to populated regions, necessitating the use of satellite-based observations to provide estimates of precipitation across the whole of the Earth’s surface. The temporal and spatial variability of precipitation requires adequate sampling, especially at finer resolutions. It is, therefore, necessary to exploit all available data from precipitation-capable satellites to ensure the proper representation of precipitation. To date, the estimation of precipitation using passive microwave observations has been largely concentrated upon the conically scanning imaging instruments, with relatively few techniques exploiting the observations made from the cross-track sounders. This paper describes the development of the Precipitation Retrieval and Profiling Scheme (PRPS) to retrieve precipitation from cross-track sensors, together with its performance against surface radar data and other satellite precipitation retrievals.

A previous study has used the stable and peculiar echoes backscattered by a single “bright scatterer” (BS) during five winter days to characterize the hardware of C-band, the dual-polarization radar located at Monte Lema (1625 m altitude) in Southern Switzerland. The BS is the 90 m tall metallic tower on Cimetta (1633 m altitude, 18 km range). In this note, the statistics of the echoes from the BS were derived from other ten dry days with normal propagation conditions in winter 2015 and January 2019. The study confirms that spectral signatures, such as spectrum width, wideband noise and Doppler velocity, were persistently stable. Regarding the polarimetric signatures, the large values (with small dispersion) of the copolar correlation coefficient between horizontal and vertical polarization were also confirmed: the average value was 0.9961 (0.9982) in winter 2015 (January 2019); the daily standard deviations were very small, ranging from 0.0007 to 0.0030. The dispersion of the differential phase shift was also confirmed to be quite small: the daily standard deviation ranged from a minimum of 2.5° to a maximum of 5.3°. Radar reflectivities in both polarizations were typically around 80 dBz and were confirmed to be among the largest values observed in the surveillance volume of the Monte Lema radar. Finally, another recent 5-day data set from January 2020 was analyzed after the replacement of the radar calibration unit that includes low noise amplifiers: these five days show poorer characteristics of the polarimetric signatures and a few outliers affecting the spectral signatures. It was shown that the “historical” polarimetric and spectral signatures of a bright scatterer could represent a benchmark for an in-depth comparison after hardware replacements.

In this paper, we build on past efforts with regard to the implementation of an efficient feature tracking algorithm for the mass processing of satellite images. This generic open-source feature tracking routine can be applied to any type of imagery to measure sub-pixel displacements between images. The routine consists of a feature tracking module (autoRIFT) that enhances computational efficiency and a geocoding module (Geogrid) that mitigates problems found in existing geocoding algorithms. When applied to satellite imagery, autoRIFT can run on a grid in the native image coordinates (such as radar or map) and, when used in conjunction with the Geogrid module, on a user-defined grid in geographic Cartesian coordinates such as Universal Transverse Mercator or Polar Stereographic. To validate the efficiency and accuracy of this approach, we demonstrate its use for tracking ice motion by using ESA’s Sentinel-1A/B radar data (seven pairs) and NASA’s Landsat-8 optical data (seven pairs) collected over Greenland’s Jakobshavn Isbræ glacier in 2017. Feature-tracked velocity errors are characterized over stable surfaces, where the best Sentinel-1A/B pair with a 6 day separation has errors in X/Y of 12 m/year or 39 m/year, compared to 22 m/year or 31 m/year for Landsat-8 with a 16-day separation. Different error sources for radar and optical image pairs are investigated, where the seasonal variation and the error dependence on the temporal baseline are analyzed. Estimated velocities were compared with reference velocities derived from DLR’s TanDEM-X SAR/InSAR data over the fast-moving glacier outlet, where Sentinel-1 results agree within 4% compared to 3–7% for Landsat-8. A comprehensive apples-to-apples comparison is made with regard to runtime and accuracy between multiple implementations of the proposed routine and the widely-used “dense ampcor" program from NASA/JPL’s ISCE software. autoRIFT is shown to provide two orders of magnitude of runtime improvement with a 20% improvement in accuracy.

Profile radar allows direct characterization of the vertical forest structure. Short-wavelength, such as Ku or X band, microwave data provide opportunities to detect the foliage. In order to exploit the potential of radar technology in forestry applications, a helicopter-borne Ku-band profile radar system, named Tomoradar, has been developed by the Finnish Geospatial Research Institute. However, how to use the profile radar waveforms to assess forest canopy parameters remains a challenge. In this study, we proposed a method by matching Tomoradar waveforms with simulated ones to estimate forest canopy leaf area index (LAI). Simulations were conducted by linking an individual tree-based forest gap model ZELIG and a three-dimension (3D) profile radar simulation model RAPID2. The ZELIG model simulated the parameters of potential local forest succession scene, and the RAPID2 model utilized the parameters to generate 3D virtual scenes and simulate waveforms based on Tomoradar configuration. The direct comparison of simulated and collected waveforms from Tomoradar could be carried out, which enabled the derivation of possible canopy LAI distribution corresponding to the Tomoradar waveform. A 600-m stripe of Tomoradar data (HH polarization) collected in the boreal forest at Evo in Finland was used as a test, which was divided into 60 plots with an interval of 10 m along the trajectory. The average waveform of each plot was employed to estimate the canopy LAI. Good results have been found in the waveform matching and the uncertainty of canopy LAI estimation. There were 95% of the plots with the mean relative overlapping rate (RO) above 0.7. The coefficients of variation of canopy LAI estimates were less than 0.20 in 80% of the plots. Compared to lidar-derived canopy effective LAI estimation, the coefficient of determination was 0.46, and the root mean square error (RMSE) was 1.81. This study established a bridge between the Ku band profile radar waveform and the forest canopy LAI by linking the RAPID2 and ZELIG model, presenting the uncertainty of forest canopy LAI estimation using Tomoradar. It is worth noting that since the difference of backscattering contribution is caused by both canopy structure and tree species, similar waveforms may correspond to different canopy LAI, inducing the uncertainty of canopy LAI estimation, which should be noticed in forest parameters estimation with empirical methods.

Over the past several decades, as sensor technology has improved, the spatial resolution of satellite images has been steadily improving [...]

The availability of free and high temporal resolution satellite data and advanced SAR techniques allows us to analyze ground displacement cost-effectively. Our aim was to properly define subsidence and uplift areas to delineate a geothermal field and perform time-series analysis to identify temporal trends. A Persistent Scatterer Interferometry (PSI) algorithm was used to estimate vertical displacement in the Brady geothermal field located in Nevada by analyzing 70 Sentinel-1A Synthetic-Aperture Radar (SAR) images, between January 2017 and December 2019. To classify zones affected by displacement, an unsupervised Self-Organizing Map (SOM) algorithm was applied to classify points based on their behavior in time, and those clusters were used to determine subsidence, uplift, and stable regions automatically. Finally, time-series analysis was applied to the clustered data to understand the inflection dates. The maximum subsidence is –19 mm/yr with an average value of –6 mm/yr within the geothermal field. The maximum uplift is 14 mm/yr with an average value of 4 mm/yr within the geothermal field. The uplift occurred on the NE of the field, where the injection wells are located. On the other hand, subsidence is concentrated on the SW of the field where the production wells are located. The coupling of the PSInSAR and the SOM algorithms was shown to be effective in analyzing the direction and pattern of the displacements observed in the field.

Soil salinization is a global problem, which seriously damages the ecological environment and considerably reduces agricultural productivity, especially in arid regions. Synthetic aperture radar (SAR) has been widely used in remote sensing due to its weather and sunlight independence. Polarimetric SAR has great potential for large-scale mapping and monitoring salt-affected soils. In this study, we investigate the characteristics of saline soil in extremely arid regions using dual-band quadrature-polarimetric (quad-pol) SAR images acquired by GF-3 (C-band) and ALOS-2 (L-band). Firstly, the effectiveness of the modified dielectric mixing model and integral equation model (IEM) in describing saline soil is evaluated. Secondly, the potential relationships between polarimetric parameters and salinity are discussed in both the C- and L-band, respectively, such as co-polarization ratio, scattering entropy H, and scattering angle α. Finally, a linear regression model for monitoring salt content is established. The main contributions of this article are as follows: (1) Simulation results suggest that the radar backscattering coefficient is a weak function of salinity at low water content, but our experimental data show that soil salinity significantly contributes to the radar backscattering coefficient, which indicates the modified dielectric mixing model and IEM model is not applicable in extremely arid areas. (2) A negative correlation between the co-polarization ratio and salinity is observed, and the correlation coefficients are 0.64 (C-band) and 0.71 (L-band). Besides, scattering entropy and scattering angle exhibit a positive correlation with salinity in the C-band with correlation coefficients 0.686 and 0.669, respectively, whereas a negative correlation is found in the L-band with correlation coefficients 0.682 and 0.680, respectively. This can be attributed to the different penetration depths and sensitivity to the surface roughness of the electromagnetic waves at two frequencies. (3) A regression model for salinity estimating based on radar backscattering coefficient, co-polarization ratio, and scattering entropy is established, with a determination coefficient (R2) of 0.79 and a root mean square error (RMSE) of 6.56%, allowing us to determine soil salinity from quad-pol SAR images without using backscattering models. Therefore, our results can be a reference for future soil salinity monitoring and inversion.

In light of recent Arctic change, there is a need to better understand sea ice dynamic processes at the floe scale to evaluate sea ice stability, deformation, and fracturing. This work investigates the use of the Gamma portable radar interferometer (GPRI) to characterize sea ice displacement and surface topography. We find that the GPRI is best suited to derive lateral surface deformation due to mm-scale horizontal accuracy. We model interferometric phase signatures from sea ice displacement and evaluate possible errors related to noise and antenna motion. We compare the analysis with observations acquired during a drifting ice camp in the Beaufort Sea. We used repeat-scan and stare-mode interferometry to identify two-dimensional shear and to track continuous uni-directional convergence. This paper demonstrates the capacity of the GPRI to derive surface strain on the order of 10−7 and identify different dynamic regions based on sub-mm changes in displacement. The GPRI is thus a promising tool for sea ice applications due to its high accuracy that can potentially resolve pre- and post-fracture deformation relevant to sea ice stability and modeling.

Spaceborne and airborne polarimetric synthetic-aperture radar interferometry (PolInSAR) data have been extensively used for forest parameter retrieval. The PolInSAR models have proven their potential in the accurate measurement of forest vegetation height. Spaceborne monostatic multifrequency data of different SAR missions and the Global Ecosystem Dynamics Investigation (GEDI)-derived forest canopy height map were used in this study for vegetation height retrieval. This study tested the performance of PolInSAR complex coherence-based inversion models for estimating the vegetation height of the forest ranges of Doon Valley, Uttarakhand, India. The inversion-based forest height obtained from the three-stage inversion (TSI) model had higher accuracy than the coherence amplitude inversion (CAI) model-based estimates. The vegetation height values of GEDI-derived canopy height map did not show good relation with field-measured forest height values. It was found that, at several locations, GEDI-derived forest height values underestimated the vegetation height. The statistical analysis of the GEDI-derived estimates with field-measured height showed a high root mean square error (RMSE; 5.82 m) and standard error (SE; 5.33 m) with a very low coefficient of determination (R2; 0.0022). An analysis of the spaceborne-mission-based forest height values suggested that the L-band SAR has great potential in forest height retrieval. TSI-model-based forest height values showed lower p-values, which indicates the significant relation between modelled and field-measured forest height values. A comparison of the results obtained from different SAR systems is discussed, and it is observed that the L-band-based PolInSAR inversion gives the most reliable result with low RMSE (2.87 m) and relatively higher R2 (0.53) for the linear regression analysis between the modelled tree height and the field data. These results indicate that higher wavelength PolInSAR datasets are more suitable for tree canopy height estimation using the PolInSAR inversion technique.

Satellite sensor systems for soil moisture measurements have been continuously evolving. The Soil Moisture Active Passive (SMAP) mission represents one of the latest advances in this regard. Thus far, much of our knowledge of the accuracy of SMAP soil moisture over the Great Lakes region of North America has originated from evaluation studies using in situ data from the U.S. Department of Agriculture (USDA) Natural Resources Conservation Service Soil Climate Analysis Network and/or the U.S. Climate Reference Network, which provide only several in situ sensor stations for this region. As such, these results typically underrepresent the accuracy of SMAP soil moisture in this region, which is characterized by a relatively large soil moisture variability and is one of the least studied regions. In this work, SMAP Level 2‒4 soil moisture products: SMAP/Sentinel-1 L2 Radiometer/Radar Soil Moisture (SPL2SMAP_S), SMAP Enhanced L3 Radiometer Soil Moisture (SPL3SMP_E), and SMAP L4 Surface and Root-Zone Soil Moisture Analysis Update (SPL4SMAU) are evaluated over the southern portion of the Great Lakes region using in situ measurements from Michigan State University’s Enviro-weather Automated Weather Station Network. The unbiased root-mean-square error (ubRMSE) values for both SPL4SMAU surface and root zone soil moisture estimates are below 0.04 m3 m−3 at the 36-km scale, with an average ubRMSE of 0.045 m3 m−3 (0.037 m3 m−3) for the surface (root-zone) soil moisture against the sparse network. The ubRMSE values for SPL3SMP_E a.m. (i.e., descending overpasses) soil moisture retrievals are close to or below 0.04 m3 m−3 at the 36-km scale, with an average ubRMSE of ~0.06 m3 m−3 against the sparse network. The average ubRMSE values are ~0.05‒0.06 m3 m−3 for high-resolution SPL2SMAP_S soil moisture retrievals against the sparse network, with the skill of the baseline algorithm-based soil moisture retrievals exceeding that of the optional algorithm-based counterparts. Clearly, the skill of SPL4SMAU surface soil moisture exceeds that of the SPL3SMP_E and SPL2SMAP_S soil moisture retrievals.

Currently, shore-based HF radars are widely used for coastal observations, and airborne radars are utilized for monitoring the ocean with a relatively large coverage offshore. In order to take the advantage of airborne radars, the theoretical mechanism of airborne HF/VHF radar for ocean surface observation has been studied in this paper. First, we describe the ocean surface wave height with the linear and nonlinear parts in a reasonable mathematical form and adopt the small perturbation method (SPM) to compute the HF/VHF radio scattered field induced by the sea surface. Second, the normalized radar cross section (NRCS) of the ocean surface is derived by tackling the field scattered from the random sea as a stochastic process. Third, the NRCS is simulated using the SPM under different sea states, at various radar operating frequencies and incident angles, and then the influences of these factors on radar sea echoes are investigated. At last, a comparison of NRCS using the SPM and the generalized function method (GFM) is done and analyzed. The mathematical model links the sea echoes and the ocean wave height spectrum, and it also offers a theoretical basis for designing a potential airborne HF/VHF radar for ocean surface remote sensing.

The skills of CODAR SeaSonde coastal high-frequency radars (HFR) located in the West Iberian Peninsula on measuring wave parameters are compared to in situ (buoy) and satellite altimeters (SA) wave observations. Significant wave heights (SWH), wave periods, and wave directions are compared over a time window of 36-months, from January 2017 to December 2019. The ability of HFR systems to capture extreme wave events is also assessed by comparing SWH measurements during the Emma storm, which hit the Iberian Peninsula in March 2018. The analysis presented in this study shows a slight overestimation of the SWH by the HFR systems. Comparisons with in situ observations revealed correlation coefficients (R) of the order of 0.69–0.87, biases below 0.60 m, root-mean-squared errors (RMSE) between 0.89 m to 1.18 m, and a slope regression between 1.01 and 1.26. Using buoy observations as reference ground truth, the comparisons with SA revealed Rs higher than 0.94, biases under 0.19 m, and RMSEs between 0.17 m and 0.42 m. Since in situ observations do not overlap all the HFR range cells (RC), and its correlation coefficients with SA have shown good agreement (R > 0.94), Sentinel-3 SA (SRAL) SWH measurements are further used for the validation of the HFR systems SWH observations. The comparison between the HFR and the SA collocated SWH observations allowed the evaluation of the ability of the radars to retrieve wave data as a function of the distance to the coast, particularly during extreme wave events. The comparison of the lower frequency (4.86 MHz) HFR coastal radars with the SA measurements showed an R of 0.94–0.99, a negative but reduced bias (−0.37), and an RMSE of 0.53 m. The higher frequency HFR systems (12–13.5 MHz) showed R between 0.53 and 0.82, and a clear overestimation of the SWH by the HFR sites.

Global navigation satellite systems reflectometry (GNSS-R) is a relatively novel remote sensing technique, but it can be understood as a multi-static radar using satellite navigation signals as signals of opportunity. The scattered signals over sea ice, flooded areas, and even under dense vegetation show a detectable coherent component that can be separated from the incoherent component and processed accordingly. This work derives an analytical formulation of the response of a GNSS-R instrument to a step function in the reflectivity using well-known principles of electromagnetic theory. The evaluation of the spatial resolution then requires a numerical evaluation of the proposed equations, as the width of the transition depends on the reflectivity values of two regions. However, it is found that results are fairly constant over a wide range of reflectivities, and they only vary faster for very high or very low reflectivity gradients. The predicted step response is then satisfactorily compared to airborne experimental results at L1 (1575.42 MHz) and L5 (1176.45 MHz) bands, acquired over a water reservoir south of Melbourne, in terms of width and ringing, and several examples are provided when the transition occurs from land to a rough ocean surface, where the coherent scattering component is no longer dominant. This work was funded by the Spanish MCIU and EU ERDF project (RTI2018-099008-B-C21/AEI/10.13039/501100011033) “Sensing with pioneering opportunistic techniques” and grant to ”CommSensLab-UPC” Excellence Research Unit Maria de Maeztu (MINECO grant MDM-2016-600).

The role of wave breaking in microwave backscattering from the sea surface is a problem of great importance for the development of theories and methods on ocean remote sensing, in particular for oil spill remote sensing. Recently it has been shown that microwave radar return is determined by both Bragg and non-Bragg (non-polarized) scattering mechanisms and some evidence has been given that the latter is associated with wave breaking, in particular, with strong breaking such as spilling or plunging. However, our understanding of mechanisms of the action of strong wave breaking on small-scale wind waves (ripples) and thus on the radar return is still insufficient. In this paper an effect of suppression of radar backscattering after strong wave breaking has been revealed experimentally and has been attributed to the wind ripple suppression due to turbulence generated by strong wave breaking. The experiments were carried out in a wind wave tank where a frequency modulated wave train of intense meter-decimeter-scale surface waves was generated by a mechanical wave maker. The wave train was compressed according to the gravity wave dispersion relation (“dispersive focusing”) into a short-wave packet at a given distance from the wave maker. Strong wave breaking with wave crest overturning (spilling) occurred for one or two highest waves in the packet. Short decimeter-centimeter-scale wind waves were generated at gentle winds, simultaneously with the long breaking waves. A Ka-band scatterometer was used to study microwave backscattering from the surface waves in the tank. The scatterometer looking at the area of wave breaking was mounted over the tank at a height of about 1 m above the mean water level, the incidence angle of the microwave radiation was about 50 degrees. It has been obtained that the radar return in the presence of short wind waves is characterized by the radar Doppler spectrum with a peak roughly centered in the vicinity of Bragg wave frequencies. The radar return was strongly enhanced in a wide frequency range of the radar Doppler spectrum when a packet of long breaking waves arrived at the area irradiated by the radar. After the passage of breaking waves, the radar return strongly dropped and then slowly recovered to the initial level. Measurements of velocities in the upper water layer have confirmed that the attenuation of radar backscattering after wave breaking is due to suppression of short wind waves by turbulence generated in the breaking zone. A physical analysis of the effect has been presented.

Canada’s successful space-based earth-observation (EO) radar program has earned widespread and expanding user acceptance following the launch of RADARSAT-1 in 1995. RADARSAT-2, launched in 2007, while providing data continuity for its predecessor’s imaging capabilities, added new polarimetric modes. Canada’s follow-up program, the RADARSAT Constellation Mission (RCM), launched in 2019, while providing continuity for its two predecessors, includes an innovative suite of polarimetric modes. In an effort to make polarimetry accessible to a wide range of operational users, RCM uses a new method called hybrid compact polarization (HCP). There are two essential elements to this approach: (1) transmit only one polarization, circular; and (2) receive two orthogonal polarizations, for which RCM uses H and V. This configuration overcomes the conventional dual and full polarimetric system limitations, which are lacking enough polarimetric information and having a small swath width, respectively. Thus, HCP data can be considered as dual-pol data, while the resulting polarimetric classifications of features in an observed scene are of comparable accuracy as those derived from the traditional fully polarimetric (FP) approach. At the same time, RCM’s HCP methodology is applicable to all imaging modes, including wide swath and ScanSAR, thus overcoming critical limitations of traditional imaging radar polarimetry for operational use. The primary image data products from an HCP radar are different from those of a traditional polarimetric radar. Because the HCP modes transmit circularly polarized signals, the data processing to extract polarimetric information requires different approaches than those used for conventional linearly polarized polarimetric data. Operational users, as well as researchers and students, are most likely to achieve disappointing results if they work with traditional polarimetric processing tools. New tools are required. Existing tutorials, older seminar notes, and reference papers are not sufficient, and if left unrevised, could succeed in discouraging further use of RCM polarimetric data. This paper is designed to provide an initial response to that need. A systematic review of studies that used HCP SAR data for environmental monitoring is also provided. Based on this review, HCP SAR data have been employed in oil spill monitoring, target detection, sea ice monitoring, agriculture, wetland classification, and other land cover applications.

Following the mission science plan of EPS/Metop-SG C-band scatterometer for 2023&ndash, 2044, we consider the potential application of the sea ice/water discrimination method based on the minimum statistical distance of the measured normalized radar cross sections (NRCS) to the geophysical model functions (GMF) of the sea ice and water, respectively. The application of the method is considered for the classical spacecraft scatterometer geometry with three fixed fan-beam antennas and the hypothetical prospective scatterometer geometry with the five fixed fan-beam antennas. Joint vertical (VV) and horizontal (HH) transmit and receive polarization are considered for the spaceborne scatterometer geometries. We show explicitly that the hypothetical five fixed fan-beam antenna geometry combined with the dual VV and HH polarization for all antennas provides better estimates of the sea wind speed and direction as well as sea ice/water discrimination during single spacecraft pass. The sea ice/water discrimination algorithms developed for each scatterometer geometry and dual VV/HH polarization are presented. The obtained results can be used to optimize the design of new spaceborne scatterometers and will be beneficial to the forthcoming satellite missions.

The total variation (TV) method has been applied to realizing airborne scanning radar super-resolution imaging while maintaining the outline of the target. The iterative reweighted norm (IRN) approach is an algorithm for addressing the minimum Lp norm problem by solving a sequence of minimum weighted L2 norm problems, and has been applied to solving the TV norm. However, during the solving process, the IRN method is required to update the weight term and result term in each iteration, involving multiplications and the inversion of large matrices. Consequently, it suffers from a huge calculation load, which seriously restricts the application of the TV imaging method. In this work, by analyzing the structural characteristics of the matrix involved in iteration, an efficient method based on suitable matrix blocking is proposed. It transforms multiplications and the inversion of large matrices into the computation of multiple small matrices, thereby accelerating the algorithm. The proposed method, called IRN-FTV method, is more time economical than the IRN-TV method, especially for high dimensional observation scenarios. Numerical results illustrate that the proposed IRN-FTV method enjoys preferable computational efficiency without performance degradation.

Stable and efficient ground moving target tracking and refocusing is a hard task in synthetic aperture radar (SAR) data processing. Since shadows in video-SAR indicate the actual positions of moving targets at different moments without any displacement, shadow-based methods provide a new approach for ground moving target processing. This paper constructs a novel framework to refocus ground moving targets by using shadows in video-SAR. To this end, an automatic-registered SAR video is first obtained using the video-SAR back-projection (v-BP) algorithm. The shadows of multiple moving targets are then tracked using a learning-based tracker, and the moving targets are ultimately refocused via a proposed moving target back-projection (m-BP) algorithm. With this framework, we can perform detecting, tracking, imaging for multiple moving targets integratedly, which significantly improves the ability of moving-target surveillance for SAR systems. Furthermore, a detailed explanation of the shadow of a moving target is presented herein. We find that the shadow of ground moving targets is affected by a target’s size, radar pitch angle, carrier frequency, synthetic aperture time, etc. With an elaborate system design, we can obtain a clear shadow of moving targets even in X or C band. By numerical experiments, we find that a deep network, such as SiamFc, can easily track shadows and precisely estimate the trajectories that meet the accuracy requirement of the trajectories for m-BP.

Summarization: This work presents the latest calibration results for the Copernicus Sentinel-3A and -3B and the Jason-3 radar altimeters as determined by the Permanent Facility for Altimetry Calibration (PFAC) in west Crete, Greece. Radar altimeters are used to provide operational measurements for sea surface height, significant wave height and wind speed over oceans. To maintain Fiducial Reference Measurement (FRM) status, the stability and quality of altimetry products need to be continuously monitored throughout the operational phase of each altimeter. External and independent calibration and validation facilities provide an objective assessment of the altimeter’s performance by comparing satellite observations with ground-truth and in-situ measurements and infrastructures. Three independent methods are employed in the PFAC: Range calibration using a transponder, sea-surface calibration relying upon sea-surface Cal/Val sites, and crossover analysis. Procedures to determine FRM uncertainties for Cal/Val results have been demonstrated for each calibration. Biases for Sentinel-3A Passes No. 14, 278 and 335, Sentinel-3B Passes No. 14, 71 and 335, as well as for Jason-3 Passes No. 18 and No. 109 are given. Diverse calibration results by various techniques, infrastructure and settings are presented. Finally, upgrades to the PFAC in support of the Copernicus Sentinel-6 ‘Michael Freilich’, due to launch in November 2020, are summarized. Presented on: Remote Sensing

With the continuous expansion of the high-speed railway network in China, long-span railway bridges carrying multiple tracks demand reliable and fast testing procedures and techniques. Bridge dynamic behavior analysis is a critical process in ensuring safe operation of structures. In this study, we present some experimental results of the vibration monitoring of a four-track high-speed railway bridge with a metro–track on each side: the Nanjing–Dashengguan high-speed railway bridge (NDHRB). The results were obtained using a terrestrial microwave radar interferometer named IBIS-S. The radar measurements were interpreted with the support of lidar point clouds. The results of the bridge dynamic response under different loading conditions, including high-speed trains, metro and wind were compared with the existing bridge structure health monitoring (SHM) system, underlining the high spatial (0.5 m) and temporal resolutions (50 Hz–200 Hz) of this technique for railway bridge dynamic monitoring. The detailed results can help engineers capturing the maximum train-induced bridge displacement. The bridge was also monitored by the radar from a lateral position with respect to the bridge longitudinal direction. This allowed us to have a more exhaustive description of the bridge dynamic behavior. The different effects induced by the passage of trains through different tracks and directions were distinguished. In addition, the space deformation map of the wide bridge deck under the eccentric load of trains, especially along the lateral direction (30 m), can help evaluating the running stability of high-speed trains.

Remote sensing data and methods are increasingly being embedded into assessments of volcanic processes and risk. This happens thanks to their capability to provide a spectrum of observation and measurement opportunities to accurately sense the dynamics, magnitude, frequency, and impacts of volcanic activity in the ultraviolet (UV), visible (VIS), infrared (IR), and microwave domains. Launched in mid-2018, the Special Issue “Remote Sensing of Volcanic Processes and Risk” of Remote Sensing gathers 19 research papers on the use of satellite, aerial, and ground-based remote sensing to detect thermal features and anomalies, investigate lava and pyroclastic flows, predict the flow path of lahars, measure gas emissions and plumes, and estimate ground deformation. The strong multi-disciplinary character of the approaches employed for volcano monitoring and the combination of a variety of sensor types, platforms, and methods that come out from the papers testify the current scientific and technology trends toward multi-data and multi-sensor monitoring solutions. The research advances presented in the published papers are achieved thanks to a wealth of data including but not limited to the following: thermal IR from satellite missions (e.g., MODIS, VIIRS, AVHRR, Landsat-8, Sentinel-2, ASTER, TET-1) and ground-based stations (e.g., FLIR cameras); digital elevation/surface models from airborne sensors (e.g., Light Detection And Ranging (LiDAR), or 3D laser scans) and satellite imagery (e.g., tri-stereo Pléiades, SPOT-6/7, PlanetScope); airborne hyperspectral surveys; geophysics (e.g., ground-penetrating radar, electromagnetic induction, magnetic survey); ground-based acoustic infrasound; ground-based scanning UV spectrometers; and ground-based and satellite Synthetic Aperture Radar (SAR) imaging (e.g., TerraSAR-X, Sentinel-1, Radarsat-2). Data processing approaches and methods include change detection, offset tracking, Interferometric SAR (InSAR), photogrammetry, hotspots and anomalies detection, neural networks, numerical modeling, inversion modeling, wavelet transforms, and image segmentation. Some authors also share codes for automated data analysis and demonstrate methods for post-processing standard products that are made available for end users, and which are expected to stimulate the research community to exploit them in other volcanological application contexts. The geographic breath is global, with case studies in Chile, Peru, Ecuador, Guatemala, Mexico, Hawai’i, Alaska, Kamchatka, Japan, Indonesia, Vanuatu, Réunion Island, Ethiopia, Canary Islands, Greece, Italy, and Iceland. The added value of the published research lies on the demonstration of the benefits that these remote sensing technologies have brought to knowledge of volcanoes that pose risk to local communities; back-analysis and critical revision of recent volcanic eruptions and unrest periods; and improvement of modeling and prediction methods. Therefore, this Special Issue provides not only a collection of forefront research in remote sensing applied to volcanology, but also a selection of case studies proving the societal impact that this scientific discipline can potentially generate on volcanic hazard and risk management.

Space-based radar sensors have transformed Earth observation since their first use by Seasat in 1978. Radar instruments are less affected by daylight or weather conditions than optical counterparts, suitable for continually monitoring the global biosphere. The current trends in synthetic aperture radar (SAR) platform design are distinct from traditional approaches in that miniaturized satellites carrying SAR are launched in multiples to form a SAR constellation. A systems engineering perspective is presented in this paper to track the transitioning of space-based SAR platforms from large satellites to small satellites. Technological advances therein are analyzed in terms of subsystem components, standalone satellites, and satellite constellations. The availability of commercial satellite constellations, ground stations, and launch services together enable real-time SAR observations with unprecedented details, which will help reveal the global biomass and their changes owing to anthropogenic drivers. The possible roles of small satellites in global biospheric monitoring and the subsequent research areas are also discussed.

An extremely heavy rainfall event hit Guangdong province, China, from 27 August to 1 September 2018. There were two different extreme rain regions, respectively, at the Pearl River estuary and eastern Guangdong, and a record-breaking daily precipitation of 1056.7 mm was observed at Gaotan station on 30 August. This paper utilizes a suite of observations from soundings, a gauge network, disdrometers, and polarimetric radars to gain insights to the two rainfall centers. The large-scale meteorological forcing, rainfall patterns, and microphysical processes, as well as radar-based precipitation signatures are investigated. It is concluded that a west-moving monsoon depression played a critical role in sustaining the moisture supply to the two extreme rain regions, and the combined orographic enhancement further contributed to the torrential rainfall over Gaotan station. The raindrop size distributions (DSD) observed at Zhuhai and Huidong stations, as well as the observed polarimetric radar signatures indicate that the rainfall at Doumen region was characterized by larger raindrops but a lower number concentration compared with that at Gaotan region. In addition, the dual-polarization radars are used to quantify precipitation intensity during this extreme event, providing timely information for flood warning and emergency management decision-making.

Ground-based radar interferometry (GBSAR) is a useful method to control the stability of engineering objects and elements of geographical spaces at risk of deformation or displacement. To secure accurate and credible measurement results, it is crucial to consider atmospheric conditions as they influence the corrections to distance measurements. These conditions are especially important considering the radar bandwidth used. Measurements for the stability of engineering objects are not always performed in locations where meteorological monitoring is prevalent; however, information about the range of variability in atmospheric corrections is always welcome. The authors present a hybrid method to estimate the probable need of atmospheric corrections, which allows partly eliminating false positive alarms of deformations as caused by atmospheric fluctuations. Unlike the numerous publications on atmospheric reductions focused on the current state of the atmosphere, the proposed solution is based on applying a classic machine learning algorithm designed for the SARIMAX (Seasonal Autoregressive Integrated Moving Average with covariate at time) time series data model for satellite data shared by NASA (National Aeronautics and Space Administration) during the Landsat MODIS (Moderate Resolution Imaging Spectroradiometer) mission before performing residual estimation during the monitoring phase. Example calculations (proof of concept) were made for ten-year satellite data covering a region for experimental flood bank stability observations as performed using the IBIS-L (Image by Interferometric Survey—Landslide) radar and for target monitoring data (ground measurements).

Analysis of synthetic aperture radar (SAR) images in L-band of short-period internal waves (IWs), and classification of their radar signatures is presented by means of a polarimetric data set from ALOS-PALSAR mission. We choose the polarimetric feature named standard deviation(std) of the co-polarized phase difference (CPD) to identify fundamental differences in SAR signatures of internal waves, and divided them into three different classes, according to their backscattered modulation depths and morphology as well as the std CPD, namely: double-signed, single-negative, and single-positive signatures, for IW normalized image transects that display, respectively, signatures in the form of bright/dark, dark, and bright bands that correspond to positive/negative, negative, or positive variations of radar backscatter. These radar power types of signatures have a counterpart in the std CPD normalized transects, and in this paper we discuss those correlations and decorrelations. We focus in the single-negative type of signature, that is dark bands on gray background, and show that the std CPD is greatly enhanced over the troughs and rear slopes of those IWs. It is suggested that such behavior is consistent with the presence of surface slicks owing to enhanced surfactant concentration. Furthermore, those single-negative SAR signatures appear at locations where and when biological productivity is enhanced. It is found that the modulation depths associated to the std CPD is higher than the one associated to the HH-polarized radar backscatter for single-negative signatures propagating in the range direction, while the reverse occurs for the other types of signatures.

In this paper, precipitation estimates derived from the Italian ground radar network (IT GR) are used in conjunction with Spinning Enhanced Visible and InfraRed Imager (SEVIRI) measurements to develop an operational oriented algorithm (RAdar INfrared Blending algorithm for Operational Weather monitoring (RAINBOW)) able to provide precipitation pattern and intensity. The algorithm evaluates surface precipitation over five geographical boxes (in which the study area is divided). It is composed of two main modules that exploit a second-degree polynomial relationship between the SEVIRI brightness temperature at 10.8 µm TB10.8 and the precipitation rate estimates from IT GR. These relationships are applied to each acquisition of SEVIRI in order to provide a surface precipitation map. The results, based on a number of case studies, show good performance of RAINBOW when it is compared with ground reference (precipitation rate map from interpolated rain gauge measurements), with high Probability of Detection (POD) and low False Alarm Ratio (FAR) values, especially for light to moderate precipitation range. At the same time, the mean error (ME) values are about 0 mmh−1, while root mean square error (RMSE) is about 2 mmh−1, highlighting a limited variability of the RAINBOW estimations. The precipitation retrievals from RAINBOW have been also compared with the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Support to Operational Hydrology and Water Management (H SAF) official microwave (MW)/infrared (IR) combined product (P-IN-SEVIRI). RAINBOW shows better performances than P-IN-SEVIRI, in terms of both detection and estimates of precipitation fields when they are compared to the ground reference. RAINBOW has been designed as an operational product, to provide complementary information to that of the national radar network where the IT GR coverage is absent, or the quality (expressed in terms of Quality Index (QI)) of the RAINBOW estimates is low. The aim of RAINBOW is to complement the radar and rain gauge network supporting the operational precipitation monitoring.

The characterization of carbon stocks and dynamics at the national level is critical for countries engaging in climate change mitigation and adaptation strategies. However, several tropical countries, including Kenya, lack the essential information typically provided by a complete national forest inventory. Here we present the most detailed and rigorous national-scale assessment of aboveground woody biomass carbon stocks and dynamics for Kenya to date. A non-parametric random forest algorithm was trained to retrieve aboveground woody biomass carbon (AGBC) for the year 2014 &plusmn, 1 and forest disturbances for the 2014&ndash, 2017 period using in situ forest inventory plot data and satellite Earth Observation (EO) data. The ecosystem carbon cycling of Kenya&rsquo, s forests and wooded grassland were assessed using a model-data fusion framework, CARDAMOM, constrained by the woody biomass datasets from this study as well as time series information on leaf area, fire events and soil organic carbon. Our EO-derived AGBC stocks were estimated as 140 Mt C for forests and 199 Mt C for wooded grasslands. The total AGBC loss during the study period was estimated as 1.89 Mt C with a dispersion below 1%. The CARDAMOM analysis estimated woody productivity to be three times larger in forests (mean = 1.9 t C ha&minus, 1 yr&minus, 1) than wooded grasslands (0.6 t C ha&minus, 1), and the mean residence time of woody C in forests (16 years) to be greater than in wooded grasslands (10 years). This study stresses the importance of carbon sequestration by forests in the international climate mitigation efforts under the Paris Agreement, but emphasizes the need to include non-forest ecosystems such as wooded grasslands in international greenhouse gas accounting frameworks.

We introduce a multiscale superpixel approach that leverages repeat-pass interferometric coherence and sparse AGB estimates from a simulated spaceborne lidar in order to extend the NISAR mission’s applicable range of aboveground biomass (AGB) in tropical forests. Airborne and spaceborne L-band radar and full-waveform airborne lidar data are used to simulate the NISAR and GEDI mission, respectively. In addition to UAVSAR data, we use spaceborne ALOS-2/PALSAR-2 imagery with 14-day temporal baseline, which is comparable to NISAR’s 12-day baseline. Our reference AGB maps are derived from the airborne LVIS data during the AfriSAR campaign for three sites (Mondah, Ogooue, and Lope). Each tropical site has mean AGB of at least 125 Mg/ha in addition to areas with AGB exceeding 700 Mg/ha. Spatially sampling from these LVIS-derived AGB reference maps, we approximate GEDI AGB estimates. To evaluate our methodology, we perform several different analyses. First, we partition each study site into low (≤100 Mg/ha) and high (>100 Mg/ha) AGB areas, in conformity with the NISAR mission requirement to provide AGB estimates for forests between 0 and 100 Mg/ha with a RMSE below 20 Mg/ha. In the low AGB areas, this RMSE requirement is satisfied in Lope and Mondah and it fell short of the requirement in Ogooue by less 3 Mg/ha with UAVSAR and 6 Mg/ha with PALSAR-2. We note that our maps have finer spatial resolution (50 m) than NISAR requires (1 hectare). In the high AGB areas, the normalized RMSE increases to 51% (i.e.

The snow depth on Antarctic sea ice is critical to estimating the sea ice thickness distribution from laser altimetry data, such as from Operation IceBridge or ICESat-2. Snow redistributed by wind collects around areas of deformed ice and forms a wide variety of feature on sea ice, the morphology of these features may provide some indication of the mean snow depth. Here, we apply a textural segmentation algorithm to classify and group similar textures to infer the distribution of snow using snow surface freeboard measurements from Operation IceBridge campaigns over the Weddell Sea. We find that texturally-similar regions have similar snow/ice ratios, even when they have different absolute snow depth measurements. This allows for the extrapolation of nadir-looking snow radar data using two-dimensional surface altimetry scans, providing a two-dimensional estimate of the snow depth with &sim, 22% error. We show that at the floe scale (&sim, 180 m), snow depth can be directly estimated from the snow surface with &sim, 20% error using deep learning techniques, and that the learned filters are comparable to standard textural analysis techniques. This error drops to &sim, 14% when averaged over 1.5 km scales. These results suggest that surface morphological information can improve remotely-sensed estimates of snow depth, and hence sea ice thickness, as compared to current methods. Such methods may be useful for reducing uncertainty in Antarctic sea ice thickness estimates from ICESat-2.

Across-track acceleration is a major source of estimation error of along-track velocity in synthetic-aperture radar (SAR) ground moving-target indication (GMTI). This paper presents the theory and a method of compensating across-track acceleration to improve the accuracy of along-track velocity estimated from single-channel SAR single-look complex data. A unique feature of the proposed method is the utilisation of phase derivatives in the Doppler frequency domain, which is effective for azimuth-compressed signals. The performance of the method was evaluated through experimental data acquired by TerraSAR-X and speed-controlled and measured vehicles. The application results demonstrate a notable improvement in along-track velocity estimates. The amount of along-track velocity correction is particularly significant when a target has irregular motion with a low signal-to-clutter ratio. A discontinuous velocity jump rather than a constant acceleration was also observed and verified through comparison between actual data and simulations. By applying this method, the capability of single-channel SAR GMTI could be substantially improved in terms of accuracy of velocity, and moving direction. However, the method is effective only if the correlation between the actual Doppler phase derivatives and a model derived from the residual Doppler rate is sufficiently high. The proposed method will be applied to X-band SAR systems of KOMPSAT-5 and -6.

The traveling public judges the quality of a road mostly by its roughness and/or ride quality. Hence, mapping, monitoring, and maintaining adequate pavement smoothness is of high importance to State Departments of Transportation in the US. Current methods rely mostly on in situ measurements and are, therefore, time consuming and costly when applied at the network scale. This paper studies the applicability of satellite radar remote sensing data, specifically, high-resolution Synthetic Aperture Radar (SAR) data acquired at X-band, to the network-wide mapping of pavement roughness of roads in the US. Based on a comparison of high-resolution X-band Cosmo-SkyMed images with road roughness data in the form of International Roughness Index (IRI) measurements, we found that X-band radar brightness generally increases when pavement roughness worsens. Based on these findings, we developed and inverted a model to distinguish well maintained road segments from segments in need of repair. Over test sites in Augusta County, VA, we found that our classification scheme reaches an overall accuracy of 92.6%. This study illustrates the capacity of X-band SAR for pavement roughness mapping and suggests that incorporating SAR into DOT operations could be beneficial.

This article belongs to the Special Issue Applications of GNSS Reflectometry for Earth Observation., As opposed to monostatic radars where incoherent backscattering dominates, in bistatic radars, such as Global Navigation Satellite Systems Reflectometry (GNSS-R), the forward scattered signals exhibit both an incoherent and a coherent component. Current models assume that either one or the other are dominant, and the calibration and geophysical parameter retrieval (e.g., wind speed, soil moisture, etc.) are developed accordingly. Even the presence of the coherent component of a GNSS reflected signal itself has been a matter of discussion in the last years. In this work, a method developed to separate the leakage of the direct signal in the reflected one is applied to a data set of GNSS-R signals collected over the ocean by the Microwave Interferometer Reflectometer (MIR) instrument, an airborne dual-band (L1/E1 and L5/E5a), multi-constellation (GPS and Galileo) GNSS-R instrument with two 19-elements antenna arrays with 4 beam-steered each. The presented results demonstrate the feasibility of the proposed technique to untangle the coherent and incoherent components from the total power waveform in GNSS reflected signals. This technique allows the processing of these components separately, which increases the calibration accuracy (as today both are mixed and processed together), allowing higher resolution applications since the spatial resolution of the coherent component is determined by the size of the first Fresnel zone (300–500 meters from a LEO satellite), and not by the size of the glistening zone (25 km from a LEO satellite). The identification of the coherent component enhances also the location of the specular reflection point by determining the peak maximum from this coherent component rather than the point of maximum derivative of the incoherent one, which is normally noisy and it is blurred by all the glistening zone contributions., This work was founded by the Spanish Ministry of Science, Innovation and Universities, “Sensing with Pioneering Opportunistic Techniques”, grant RTI2018-099008-B-C21, and the grant for recruitment of early-stage research staff FI-DGR 2015 of the AGAUR—Generalitat de Catalunya (FEDER), Spain, and Unidad de Excelencia María de Maeztu MDM-2016-060.

As the first Geostationary Satellite with the LMI (Lightning Mapping Imager) instrument aboard running over the eastern hemisphere, FY-4A (Feng-Yun-4A) can better indicate severe convection and compensate for the limitations of radar observation in temporal and spatial resolution. In order to realize the application of FY-4A lightning data in numerical weather prediction (NWP) models, a logarithmic relationship between FY-4A lightning density and maximum radar reflectivity is presented to convert FY-4A lightning data into maximum FY-4A proxy reflectivity. Then, according to the profiles of radar reflectivity, the maximum FY-4A proxy reflectivity is extended to 3D FY-4A proxy reflectivity. Finally, the 3D FY-4A proxy reflectivity is assimilated in RMAPS-ST (Rapid-refresh Multi-scale Analysis and Prediction System—Short Term) to compare with radar assimilation. Four groups of continuous cycling data assimilation and forecasting experiments are carried out for a severe rainfall case. The results demonstrate that cycling assimilation of 3D FY-4A proxy reflectivity can adjust the moisture condition effectively, and indirectly affects the temperature and wind fields, then makes the thermal and dynamic analysis more reasonable. The Fractions Skill Scores (FSSs) show that the rainfall forecasts are improved significantly within 6 h by assimilating 3D FY-4A proxy reflectivity, which is similar to the parallel experiments in assimilating radar reflectivity. In addition, other cycling data assimilation experiments are carried out in mountainous areas without radar data. The improvement of precipitation forecasts in mountainous areas further proves that the application of assimilating 3D FY-4A proxy reflectivity can be considered a useful substitute where observed radar data are missing. Through the two severe rainfall cases, this method could be framed as an example of how to use lightning for data assimilation.

This article presents the possibility of helping navigators direct the movement of an object, while safely passing through other objects, using an artificial neural network and optimization methods. It has been shown that the best trajectory of an object in terms of optimality and security, from among many possible options, can be determined by the method of dynamic programming with the simultaneous use of an artificial neural network, by depicting the encountered objects as moving in forbidden domains. Analytical considerations are illustrated with examples of simulation studies of the developed calculation program on real navigational situations at sea. This research took into account both the number of objects encountered and the different shapes of domains assigned to the objects encountered. Finally, the optimal value of the safe object trajectory time was compared on the setpoint value of the safe passing distance of objects in given visibility conditions at sea, and the degree of discretization of calculations was determined by the density of the location of nodes along the route of objects.

High-frequency (HF) radars are efficient tools for measuring vast areas and gathering ocean parameters in real-time. However, the accuracy of their wave estimates is under analysis. This paper presents a new methodology for analyzing and validating the wave data estimated by two CODAR SeaSonde radars located on the Galician coast (NW Spain). Approximately one and a half years of wave data (January, 2014–April, 2015) were obtained for ten range cells employing two different sampling times used by the radar software. The resulting data were screened by an updated method, and their abundance and quality were described for each radar range cell and different wave regime; the latter were defined using the spectral significant wave height (Hm0) and mean wave direction (Dm) estimated by two buoys and three SIMAR points (SImulación MARina in Spanish, from the wave reanalysis model by Puertos del Estado (PdE)). The correlation between the results and the particularities of the different sea states (broadband or bimodal), the wind and the operation of the devices are discussed. Most HF radar wave parameters’ errors occur for waves from the NNE and higher than 6 m. The best agreement between the Vilán radar and the Vilano-Sisargas buoy wave data was obtained for the dominant wave regime (from the northwest) and the southwest wave regime. However, relevant contradictions regarding wave direction were detected. The possibilities of reducing the wave parameters’ processing time by one hour and increasing the numbers of range cells of the radars have been validated This research was funded by Interreg Atlantic Area project MyCOAST (EAPA 285/2016) and INTERREG V-A Spain-Portugal (POCTEP) project RADAR_ON_RAIA co-funded by the European Regional Development Fund (ERDF) (EU). V.P.-M. and A.B. acknowledge financial support by CRETUS strategic partnership (ED431E2018/01), co-funded by the ERDF (EU) http://www.usc.es/cretus/ SI

Airborne Synthetic Aperture Radar (SAR) systems are gaining increasing interest within the remote sensing community due to their operational flexibility and observation capabilities. Among these systems, those exploiting the Frequency-Modulated Continuous-Wave (FMCW) technology are compact, lightweight, and comparatively low cost. For these reasons, they are becoming very attractive, since they can be easily mounted onboard ever-smaller and highly flexible aerial platforms, like helicopters or unmanned aerial vehicles (UAVs). In this work, we present the imaging and topographic capabilities of a novel Italian airborne SAR system developed in the frame of cooperation between a public research institute (IREA-CNR) and a private company (Elettra Microwave S.r.l.). The system, which is named AXIS (standing for Airborne X-band Interferometric SAR), is based on FMCW technology and is equipped with a single-pass interferometric layout. In the work we first provide a description of the AXIS system. Then, we describe the acquisition campaign carried out in April 2018, just after the system completion. Finally, we perform an analysis of the radar data acquired during the campaign, by presenting a quantitative assessment of the quality of the SLC (Single Look Complex) SAR images and the interferometric products achievable through the system. The overall analysis aims at providing first reference values for future research and operational activities that will be conducted with this sensor.

The effect of cloud parallax shift occurs in satellite imaging, particularly for high angles of satellite observations. This study demonstrates new methods of parallax effect correction for clouds observed by geostationary satellites. The analytical method that could be found in literature, namely the Vicente et al./Koenig method, is presented at the beginning. It approximates a cloud position using an ellipsoid with semi-axes increased by the cloud height. The error values of this method reach up to 50 meters. The second method, which is proposed by the author, is an augmented version of the Vicente et al./Koenig approach. With this augmentation, the error can be reduced to centimeters. The third method, also proposed by the author, incorporates geodetic coordinates. It is described as a set of equations that are solved with the numerical method, and its error can be driven to near zero by adjusting the count of iterations. A sample numerical solution procedure with application of the Newton method is presented. Also, a simulation experiment that evaluates the proposed methods is described in the paper. The results of an experiment are described and contrasted with current technology. Currently, operating geostationary Earth Observation (EO) satellite resolutions vary from 0.5 km up to 8 km. The pixel sizes of these satellites are much greater than for maximal error of the least precise method presented in this paper. Therefore, the chosen method will be important when the resolution of geostationary EO satellites reaches 50 m. To validate the parallax correction, procedure data from on-ground radars and the Meteosat Second Generation (MSG) satellite, which describes stormy events, was compared before and after correction. Comparison was performed by correlating the logarithm of the cloud optical thickness (COT) with radar reflectance in dBZ (radar reflectance – Z in logarithmic form).