Mapping tropical forest aboveground biomass using airborne SAR tomography

Ramachandran, N.; Department of Civil Engineering, India; email: [email protected] Researc

The BIOMASS level 2 prototype processor : design and experimental results of above-ground biomass estimation

BIOMASS is ESA’s seventh Earth Explorer mission, scheduled for launch in 2022. The satellite will be the first P-band SAR sensor in space and will be operated in fully polarimetric interferometric and tomographic modes. The mission aim is to map forest above-ground biomass (AGB), forest height (FH) and severe forest disturbance (FD) globally with a particular focus on tropical forests. This paper presents the algorithms developed to estimate these biophysical parameters from the BIOMASS level 1 SAR measurements and their implementation in the BIOMASS level 2 prototype processor with a focus on the AGB product. The AGB product retrieval uses a physically-based inversion model, using ground-canceled level 1 data as input. The FH product retrieval applies a classical PolInSAR inversion, based on the Random Volume over Ground Model (RVOG). The FD product will provide an indication of where significant changes occurred within the forest, based on the statistical properties of SAR data. We test the AGB retrieval using modified airborne P-Band data from the AfriSAR and TropiSAR campaigns together with reference data from LiDAR-based AGB maps and plot-based ground measurements. For AGB estimation based on data from a single heading, comparison with reference data yields relative Root Mean Square Difference (RMSD) values mostly between 20% and 30%. Combining different headings in the estimation process significantly improves the AGB retrieval to slightly less than 20%. The experimental results indicate that the implemented retrieval scheme provides robust results that are within mission requirements

Temporal Characteristics of Boreal Forest Radar Measurements

Radar observations of forests are sensitive to seasonal changes, meteorological variables and variations in soil and tree water content. These phenomena cause temporal variations in radar measurements, limiting the accuracy of tree height and biomass estimates using radar data. The temporal characteristics of radar measurements of forests, especially boreal forests, are not well understood. To fill this knowledge gap, a tower-based radar experiment was established for studying temporal variations in radar measurements of a boreal forest site in southern Sweden. The work in this thesis involves the design and implementation of the experiment and the analysis of data acquired. The instrument allowed radar signatures from the forest to be monitored over timescales ranging from less than a second to years. A purpose-built, 50 m high tower was equipped with 30 antennas for tomographic imaging at microwave frequencies of P-band (420-450 MHz), L-band (1240-1375 MHz) and C-band (5250-5570 MHz) for multiple polarisation combinations. Parallel measurements using a 20-port vector network analyser resulted in significantly shorter measurement times and better tomographic image quality than previous tower-based radars. A new method was developed for suppressing mutual antenna coupling without affecting the range resolution. Algorithms were developed for compensating for phase errors using an array radar and for correcting for pixel-variant impulse responses in tomographic images. Time series results showed large freeze/thaw backscatter variations due to freezing moisture in trees. P-band canopy backscatter variations of up to 10 dB occurred near instantaneously as the air temperature crossed 0⁰C, with ground backscatter responding over longer timescales. During nonfrozen conditions, the canopy backscatter was very stable with time. Evidence of backscatter variations due to tree water content were observed during hot summer periods only. A high vapour pressure deficit and strong winds increased the rate of transpiration fast enough to reduce the tree water content, which was visible as 0.5-2 dB backscatter drops during the day. Ground backscatter for cross-polarised observations increased during strong winds due to bending tree stems. Significant temporal decorrelation was only seen at P-band during freezing, thawing and strong winds. Suitable conditions for repeat-pass L-band interferometry were only seen during the summer. C-band temporal coherence was high over timescales of seconds and occasionally for several hours for night-time observations during the summer. Decorrelation coinciding with high transpiration rates was observed at L- and C-band, suggesting sensitivity to tree water dynamics.The observations from this experiment are important for understanding, modelling and mitigating temporal variations in radar observables in forest parameter estimation algorithms. The results also are also useful in the design of spaceborne synthetic aperture radar missions with interferometric and tomographic capabilities. The results motivate the implementation of single-pass interferometric synthetic aperture radars for forest applications at P-, L- and C-band

The NASA AfriSAR campaign: Airborne SAR and lidar measurements of tropical forest structure and biomass in support of current and future space missions

International audienceIn 2015 and 2016, the AfriSAR campaign was carried out as a collaborative effort among international space and National Park agencies (ESA, NASA, ONERA, DLR, ANPN and AGEOS) in support of the upcoming ESA BIOMASS, NASA-ISRO Synthetic Aperture Radar (NISAR) and NASA Global Ecosystem Dynamics Initiative (GEDI) missions. The NASA contribution to the campaign was conducted in 2016 with the NASA LVIS (Land Vegetation and Ice Sensor) Lidar, the NASA L-band UAVSAR (Uninhabited Aerial Vehicle Synthetic Aperture Radar). A central motivation for the AfriSAR deployment was the common AGBD estimation requirement for the three future spaceborne missions, the lack of sufficient airborne and ground calibration data covering the full range of ABGD in tropical forest systems, and the intercomparison and fusion of the technologies. During the campaign, over 7000 km2 of waveform Lidar data from LVIS and 30,000 km2 of UAVSAR data were collected over 10 key sites and transects. In addition, field measurements of forest structure and biomass were collected in sixteen 1-hectare sized plots. The campaign produced gridded Lidar canopy structure products, gridded aboveground biomass and associated uncertainties, Lidar based vegetation canopy cover profile products, Polarimetric Interferometric SAR and Tomographic SAR products and field measurements. Our results showcase the types of data products and scientific results expected from the spaceborne Lidar and SAR missions; we also expect that the AfriSAR campaign data will facilitate further analysis and use of waveform lidar and multiple baseline polarimetric SAR datasets for carbon cycle, biodiversity, water resources and more applications by the greater scientific community

Very High Resolution Tomographic SAR Inversion for Urban Infrastructure Monitoring — A Sparse and Nonlinear Tour

The topic of this thesis is very high resolution (VHR) tomographic SAR inversion for urban infrastructure monitoring. To this end, SAR tomography and differential SAR tomography are demonstrated using TerraSAR-X spotlight data for providing 3-D and 4-D (spatial-temporal) maps of an entire high rise city area including layover separation and estimation of deformation of the buildings. A compressive sensing based estimator (SL1MMER) tailored to VHR SAR data is developed for tomographic SAR inversion by exploiting the sparsity of the signal. A systematic performance assessment of the algorithm is performed regarding elevation estimation accuracy, super-resolution and robustness. A generalized time warp method is proposed which enables differential SAR tomography to estimate multi-component nonlinear motion. All developed methods are validated with both simulated and extensive processing of large volumes of real data from TerraSAR-X

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

Radar Imaging in Challenging Scenarios from Smart and Flexible Platforms

undefine

Modeling of Subsurface Scattering from Ice Sheets for Pol-InSAR Applications

Remote sensing is a fundamental tool to measure the dynamics of ice sheets and provides valuable information for ice sheet projections under a changing climate. There is, however, the potential to further reduce the uncertainties in these projections by developing innovative remote sensing methods. One of these remote sensing techniques, the polarimetric synthetic aperture radar interferometry (Pol-InSAR), is known since decades to have the potential to assess the geophysical properties below the surface of ice sheets, because of the penetration of microwave signals into dry snow, firn, and ice. Despite this, only very few studies have addressed this topic and the development of robust Pol-InSAR applications is at an early stage. Two potential Pol-InSAR applications are identified as the motivation for this thesis. First, the estimation and compensation of the penetration bias in digital elevation models derived with SAR interferometry. This bias can lead to errors of several meters or even tens of meters in surface elevation measurements. Second, the estimation of geophysical properties of the subsurface of glaciers and ice sheets using Pol-InSAR techniques. There is indeed potential to derive information about melt-refreeze processes within the firn, which are related to density and affect the mass balance. Such Pol-InSAR applications can be a valuable information source with the potential for monthly ice sheet wide coverage and high spatial resolution provided by the next generation of SAR satellites. However, the required models to link the Pol-InSAR measurements to the subsurface properties are not yet established. The aim of this thesis is to improve the modeling of the vertical backscattering distribution in the subsurface of ice sheets and its effect on polarimetric interferometric SAR measurements at different frequencies. In order to achieve this, polarimetric interferometric multi-baseline SAR data at different frequencies and from two different test sites on the Greenland ice sheet are investigated. This thesis contributes with three concepts to a better understanding and to a more accurate modeling of the vertical backscattering distribution in the subsurface of ice sheets. First, the integration of scattering from distinct subsurface layers. These are formed by refrozen melt water in the upper percolation zone and cause an interesting coherence undulation pattern, which cannot be explained with previously existing models. This represents a first link between Pol-InSAR data and geophysical subsurface properties. The second step is the improved modeling of the general vertical backscattering distribution of the subsurface volume. The advantages of more flexible volume models are demonstrated, but interestingly, the simple modification of a previously existing model with a vertical shift parameter lead to the best agreement between model and data. The third contribution is the model based compensation of the penetration bias, which is experimentally validated. At the investigated test sites, it becomes evident that the model based estimates of the surface elevations are more accurate than the interferometric phase center locations, which are conventionally used to derive surface elevations of ice sheets. This thesis therefore improves the state of the art of subsurface scattering modeling for Pol-InSAR applications, demonstrates the model-based penetration bias compensation, and makes a further research step towards the retrieval of geophysical subsurface information with Pol-InSAR

Study of the speckle noise effects over the eigen decomposition of polarimetric SAR data: a review

This paper is focused on considering the effects of speckle noise on the eigen decomposition of the co- herency matrix. Based on a perturbation analysis of the matrix, it is possible to obtain an analytical expression for the mean value of the eigenvalues and the eigenvectors, as well as for the Entropy, the Anisotroopy and the dif- ferent a angles. The analytical expressions are compared against simulated polarimetric SAR data, demonstrating the correctness of the different expressions.Peer ReviewedPostprint (published version