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

Estimation of Forest Height Using Spaceborne Repeat-Pass L-Band InSAR Correlation Magnitude over the US State of Maine

This paper describes a novel, simple and efficient approach to estimate forest height over a wide region utilizing spaceborne repeat-pass InSAR correlation magnitude data at L-band. We start from a semi-empirical modification of the RVoG model that characterizes repeat-pass InSAR correlation with large temporal baselines (e.g., 46 days for ALOS) by taking account of the temporal change effect of dielectric fluctuation and random motion of scatterers. By assuming (1) the temporal change parameters and forest backscatter profile/extinction coefficient follow some mean behavior across each inteferogram; (2) there is minimal ground scattering contribution for HV-polarization; and (3) the vertical wavenumber is small, a simplified inversion approach is developed to link the observed HV-polarized InSAR correlation magnitude to forest height and validated using ALOS/PALSAR repeat-pass observations against LVIS lidar heights over the Howland Research Forest in central Maine, US (with RMSE \u3c 4 m at a resolution of 32 hectares). The model parameters derived from this supervised regression are used as the basis for propagating the estimates of forest height to available interferometric pairs for the entire state of Maine, thus creating a state-mosaic map of forest height. The present approach described here serves as an alternative and complementary tool for other PolInSAR inversion techniques when full-polarization data may not be available. This work is also meant to be an observational prototype for NASA’s DESDynI-R (now called NISAR) and JAXA’s ALOS-2 satellite missions

The European Space Agency BIOMASS mission: Measuring forest above-ground biomass from space

The primary objective of the European Space Agency's 7th Earth Explorer mission, BIOMASS, is to determine the worldwide distribution of forest above-ground biomass (AGB) in order to reduce the major uncertainties in calculations of carbon stocks and fluxes associated with the terrestrial biosphere, including carbon fluxes associated with Land Use Change, forest degradation and forest regrowth. To meet this objective it will carry, for the first time in space, a fully polarimetric P-band synthetic aperture radar (SAR). Three main products will be provided: global maps of both AGB and forest height, with a spatial resolution of 200 m, and maps of severe forest disturbance at 50 m resolution (where “global” is to be understood as subject to Space Object tracking radar restrictions). After launch in 2022, there will be a 3-month commissioning phase, followed by a 14-month phase during which there will be global coverage by SAR tomography. In the succeeding interferometric phase, global polarimetric interferometry Pol-InSAR coverage will be achieved every 7 months up to the end of the 5-year mission. Both Pol-InSAR and TomoSAR will be used to eliminate scattering from the ground (both direct and double bounce backscatter) in forests. In dense tropical forests AGB can then be estimated from the remaining volume scattering using non-linear inversion of a backscattering model. Airborne campaigns in the tropics also indicate that AGB is highly correlated with the backscatter from around 30 m above the ground, as measured by tomography. In contrast, double bounce scattering appears to carry important information about the AGB of boreal forests, so ground cancellation may not be appropriate and the best approach for such forests remains to be finalized. Several methods to exploit these new data in carbon cycle calculations have already been demonstrated. In addition, major mutual gains will be made by combining BIOMASS data with data from other missions that will measure forest biomass, structure, height and change, including the NASA Global Ecosystem Dynamics Investigation lidar deployed on the International Space Station after its launch in December 2018, and the NASA-ISRO NISAR L- and S-band SAR, due for launch in 2022. More generally, space-based measurements of biomass are a core component of a carbon cycle observation and modelling strategy developed by the Group on Earth Observations. Secondary objectives of the mission include imaging of sub-surface geological structures in arid environments, generation of a true Digital Terrain Model without biases caused by forest cover, and measurement of glacier and icesheet velocities. In addition, the operations needed for ionospheric correction of the data will allow very sensitive estimates of ionospheric Total Electron Content and its changes along the dawn-dusk orbit of the mission

Temporal Characteristics of P-band Tomographic Radar Backscatter of a Boreal Forest

Temporal variations in synthetic aperture radar (SAR) backscatter over forests are of concern for any SAR mission with the goal of estimating forest parameters from SAR data. In this article, a densely sampled, two-year long time series of P-band (420 to 450 MHz) boreal forest backscatter, acquired by a tower-based radar, is analyzed. The experimental setup provides time series data at multiple polarizations. Tomographic capabilities allow the separation of backscatter at different heights within the forest. Temporal variations of these multi-polarimetric, tomographic radar observations are characterized and quantified. The mechanisms studied are seasonal variations, effects of freezing conditions, diurnal variations, effects of wind and the effects of rainfall on backscatter. An emphasis is placed on upper-canopy backscatter, which has been shown to be a robust proxy for forest biomass. The canopy backscatter was most sensitive to freezing conditions but was more stable than ground-level backscatter and full-forest backscatter during non-frozen conditions. The analysis connects tree water transport mechanisms and P-band radar backscatter for the first time. The presented results are useful for designing boreal forest parameter estimation algorithms, using data from P-band SARs, that are robust to temporal variations in backscatter. The results also present new forest remote sensing opportunities using P-band radars

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

Temporal survey of P- A nd L-band polarimetric backscatter in boreal forests

Environmental conditions and seasonal variations affect the backscattered radar signal from a forest. This potentially causes errors in a biomass retrieval scheme using data from the synthetic aperture radar (SAR) data. A better understanding of these effects and the electromagnetic scattering mechanisms in forests is required to improve biomass estimation algorithms for current and upcoming P- A nd L-band SAR missions. In this paper, temporal changes in HH-, VV-, and HV-polarized P- A nd L-band radar backscatter and temporal coherence from a boreal forest site are analyzed in relation to environmental parameters. The radar data were collected from a stand of mature Norway spruce ( Picea abies (L.) Karst.) with an above-ground biomass of approximately 250 tons/ha at intervals of 5 min from January to August 2017 using the BorealScat tower-based scatterometer. It was observed that subzero temperatures during the winters cause large variations (4 to 10 dB) in P- A nd L-band backscatter, for which the HH/VV backscatter ratio offered some mitigation. High wind speeds were also seen to cause deviations in the average backscatter at P-band due to decreased double-bounce scattering. Severe temporal decorrelation was observed at L-band over timescales of days or more, whereas the P-band temporal coherence remained high (> 0.9) for at least a month neglecting windy periods. Temporal coherence at P-band was highest during night times when wind speeds are low

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

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

Research progress on geosynchronous synthetic aperture radar

Based on its ability to obtain two-dimensional (2D) high-resolution images in all-time and all-weather conditions, spaceborne synthetic aperture radar (SAR) has become an important remote sensing technique and the study of such systems has entered a period of vigorous development. Advanced imaging modes such as radar interferometry, tomography, and multi-static imaging, have been demonstrated. However, current in-orbit spaceborne SARs, which all operate in low Earth orbits, have relatively long revisit times ranging from several days to dozens of days, restricting their temporal sampling rate. Geosynchronous SAR (GEO SAR) is an active research area because it provides significant new capability, especially its much-improved temporal sampling. This paper reviews the research progress of GEO SAR technologies in detail. Two typical orbit schemes are presented, followed by the corresponding key issues, including system design, echo focusing, main disturbance factors, repeat-track interferometry, etc, inherent to these schemes. Both analysis and solution research of the above key issues are described. GEO SAR concepts involving multiple platforms are described, including the GEO SAR constellation, GEO-LEO/airborne/unmanned aerial vehicle bistatic SAR, and formation flying GEO SAR (FF-GEO SAR). Due to the high potential of FF-GEO SAR for three-dimensional (3D) deformation retrieval and coherence-based SAR tomography (TomoSAR), we have recently carried out some research related to FF-GEO SAR. This research, which is also discussed in this paper, includes developing a formation design method and an improved TomoSAR processing algorithm. It is found that GEO SAR will continue to be an active topic in the aspect of data processing and multi-platform concept in the near future