On the Sensitivity of TanDEM-X-Observations to Boreal Forest Structure

The structure of forests is important to observe for understanding coupling to global dynamics of ecosystems, biodiversity, and management aspects. In this paper, the sensitivity of X-band to boreal forest stem volume and to vertical and horizontal structure in the form of forest height and horizontal vegetation density is studied using TanDEM-X satellite observations from two study sites in Sweden: Remningstorp and Krycklan. The forest was analyzed with the Interferometric Water Cloud Model (IWCM), without the use of local data for model training, and compared with measurements by Airborne Lidar Scanning (ALS). On one hand, a large number of stands were studied, and in addition, plots with different types of changes between 2010 and 2014 were also studied. It is shown that the TanDEM-X phase height is, under certain conditions, equal to the product of the ALS quantities for height and density. Therefore, the sensitivity of phase height to relative changes in height and density is the same. For stands with a phase height >5 m we obtained an root-mean-square error, RMSE, of 8% and 10% for tree height in Remningstorp and Krycklan, respectively, and for vegetation density an RMSE of 13% for both. Furthermore, we obtained an RMSE of 17% for estimation of above ground biomass at stand level in Remningstorp and in Krycklan. The forest changes estimated with TanDEM-X/IWCM and ALS are small for all plots except clear cuts but show similar trends. Plots without forest management changes show a mean estimated height growth of 2.7% with TanDEM-X/IWCM versus 2.1% with ALS and a biomass growth of 4.3% versus 4.2% per year. The agreement between the estimates from TanDEM-X/IWCM and ALS is in general good, except for stands with low phase height

Boreal Forest Properties from TanDEM-X Data Using Interferometric Water Cloud Model and Implications for a Bistatic C-Band Mission

Data from TanDEM-X in single-pass and bistatic interferometric mode together with the interferometric water cloud model (IWCM) can provide estimates of forest height and stem volume (or the related above-ground biomass) of boreal forests with high accuracy. We summarize results from two boreal test sites using two approaches, i.e., 1) based on model calibration using reference insitu stands, and 2) based on minimization of a cost function. Both approaches are based on inversion of IWCM, which models the complex coherence and backscattering coefficient of a homogeneous forest layer, which includes gaps where free-space wave propagation is assumed. A digital terrain model of the ground is also needed. IWCM is used to estimate forest height or stem volume, since the two variables are assumed to be related through an allometric equation. A relationship between the fractional area of gaps, the area-fill, and stem volume is also required to enable model inversion. The accuracy of the stem volume estimate in the two sites varies between 16% and 21% for height of ambiguity <100 m. The results clearly show the importance of using summer-time acquisitions. Based on the TanDEM-X results at X-band, C-band data from the ERS-1/ERS-2 tandem mission are revisited to investigate the potential of a future bistatic C-band interferometric mission. Out of nine ERS-1/ERS-2 pairs, only one pair was found to be acquired at summer temperatures, without precipitation and with high coherence. A simulated bistatic phase height is shown to give approximately the same sensitivity to stem volume as TanDEM-X

Simultaneous Estimation of Sub-canopy Topography and Forest Height with Single-baseline Single-polarization TanDEM-X Interferometric Data Combined with ICESat-2 Data

To address the challenge of retrieving sub-canopy topography using single-baseline single-polarization TanDEM-X InSAR data, we propose a novel InSAR processing framework. Our methodology begins by employing the SINC model to estimate the penetration depth (PD). Subsequently, we establish a linear relationship between PD and phase center height (PCH) to generate a wall-to-wall PCH product. To achieve this, space-borne LiDAR data are employed to capture the elevation bias between actual ground elevation and InSAR-derived elevation. Finally, the sub-canopy topography is derived by subtracting the PCH from the conventional InSAR-based DEM. Moreover, this approach enables the simultaneous estimation of forest height from single-baseline TanDEM-X data by combining the estimated PD and PCH components. The approach has been validated against Airborne Lidar Scanning data over four diverse sites encompassing different forest types, terrain conditions, and climates. The derived sub-canopy topography in the boreal and hemi-boreal forest sites (Krycklan and Remningstorp) demonstrated notable improvement in accuracy. Additionally, the winter acquisitions outperformed the summer ones in terms of inversion accuracy. The achieved RMSEs for the winter scenarios were 2.45 m and 3.83 m, respectively, representing a 50% improvement over the InSAR-based DEMs. And the forest heights are also close to the ALS measurements, with RMSEs of 2.70 m and 3.33 m, respectively. For the Yanguas site in Spain, characterized by rugged terrain, sub-canopy topography in forest areas was estimated with an accuracy of 4.27m, a 35% improvement over the original DEM. For the denser tropical forest site, only an average elevation bias could be corrected.This work is funded by the National Key R&D Program of China (No. 2022YFB3902605), the National Natural Science Foundation of China (Nos. 42227801, 42030112, 42204024, 42104016, 42330717), the Spanish Ministry of Science and Innovation (State Agency of Research, AEI) and the European Funds for Regional Development under Project PID2020-117303GB-C22/AEI/10.13039/501100011033, the Natural Science Foundation for Excellent Young Scholars of Hunan Province (No. 2023JJ20061), and in part by the China Scholarship Council Foundation to the Joint Ph.D. Studies at University of Alicante (No. 202106370125)

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

Deep Learning Model Transfer in Forest Mapping Using Multi-Source Satellite SAR and Optical Images

Deep learning (DL) models are gaining popularity in forest variable prediction using Earth observation (EO) images. However, in practical forest inventories, reference datasets are often represented by plot- or stand-level measurements, while high-quality representative wall-to-wall reference data for end-to-end training of DL models are rarely available. Transfer learning facilitates expansion of the use of deep learning models into areas with sub-optimal training data by allowing pretraining of the model in areas where high-quality teaching data are available. In this study, we perform a “model transfer” (or domain adaptation) of a pretrained DL model into a target area using plot-level measurements and compare performance versus other machine learning models. We use an earlier developed UNet based model (SeUNet) to demonstrate the approach on two distinct taiga sites with varying forest structure and composition. The examined SeUNet model uses multi-source EO data to predict forest height. Here, EO data are represented by a combination of Copernicus Sentinel-1 C-band SAR and Sentinel-2 multispectral images, ALOS-2 PALSAR-2 SAR mosaics and TanDEM-X bistatic interferometric radar data. The training study site is located in Finnish Lapland, while the target site is located in Southern Finland. By leveraging transfer learning, the SeUNet prediction achieved root mean squared error (RMSE) of (Formula presented.) m and R2 of 0.882, considerably more accurate than traditional benchmark methods. We expect such forest-specific DL model transfer can be suitable also for other forest variables and other EO data sources that are sensitive to forest structure.</p

National Forest Biomass Mapping Using the Two-Level Model

This article uses the two-level model (TLM) to predict above-ground biomass (AGB) from TanDEM-X synthetic aperture radar (SAR) data for Sweden. The SAR data were acquired between October 2015 and January 2016 and consisted of 420 scenes. The AGB was estimated from forest height and canopy density estimates obtained from TLM inversion with a power law model. The model parameters were estimated separately for each satellite scene. The prediction accuracy at stand-level was evaluated using field inventoried references from entire Sweden 2017, provided by a forestry company. AGB estimation performance varied throughout the country, with smaller errors in the north and larger in the south, but when the errors were expressed in relative terms, this pattern vanished. The error in terms of root mean square error (RMSE) was 45.6 and 27.2 t/ha at the plot- and stand-level, respectively, and the corresponding biases were -8.80 and 11.2 t/ha. When the random errors related to using sampled field references were removed, the RMSE decreased about 24% to 20.7 t/ha at the stand-level. Overall, the RMSE was of similar order to that obtained in a previous study (27-30 t/ha), where one linear regression model was used for all scenes in Sweden. It is concluded that, using the power law model with parameters estimated for each scene, the scene-wise variations decreased

Estimation of change in forest variables using synthetic aperture radar

Large scale mapping of changes in forest variables is needed for both environmental monitoring, planning of climate actions and sustainable forest management. Remote sensing can be used in conjunction with field data to produce wall-to-wall estimates that are practically impossible to produce using traditional field surveys. Synthetic aperture radar (SAR) can observe the forest independent of sunlight, clouds, snow, or rain, providing reliable high frequency coverage. Its wavelength determines the interaction with the forest, where longer wavelengths interact with larger structures of the trees, and shorter wavelengths interact mainly with the top part of the canopy, meaning that it can be chosen to fit specific applications. This thesis contains five studies conducted on the Remningstorp test site in southern Sweden. Studies I – III predicted above ground biomass (AGB) change using long wavelength polarimetric P- (in I) and L-band (in I – III) SAR data. The differences between the bands were small in terms of prediction quality, and the HV polarization, just as for AGB state prediction, was the polarization channel most correlated with AGB change. A moisture correction for L-band data was proposed and evaluated, and it was found that certain polarimetric measures were better for predicting AGB change than all of the polarization channels together. Study IV assessed the detectability of silvicultural treatments in short wavelength TanDEM-X interferometric phase heights. In line with earlier studies, only clear cuts were unambiguously distinguishable. Study V predicted site index and stand age by fitting height development curves to time series of TanDEM-X data. Site index and age were unbiasedly predicted for untreated plots, and the RMSE would likely decrease with longer time series. When stand age was known, SI was predicted with an RMSE comparable to that of the field based measurements. In conclusion, this thesis underscores SAR data's potential for generalizable methods for estimation of forest variable changes

Estimation of Forest Biomass From Two-Level Model Inversion of Single-Pass InSAR Data

A model for aboveground biomass estimation from single-pass interferometric synthetic aperture radar (InSAR) data is presented. Forest height and canopy density estimates Delta h and eta(0), respectively, obtained from two-level model (TLM) inversion, are used as biomass predictors. Eighteen bistatic VV-polarized TanDEM-X (TDM) acquisitions are used, made over two Swedish test sites in the summers of 2011, 2012, and 2013 (nominal incidence angle: 41 degrees; height-of-ambiguity: 32-63 m). Remningstorp features a hemiboreal forest in southern Sweden, with flat topography and where 32 circular plots have been sampled between 2010 and 2011 (area: 0.5 ha; biomass: 42-242 t/ha; height: 14-32 m). Krycklan features a boreal forest in northern Sweden, 720-km north-northeast from Remningstorp, with significant topography and where 31 stands have been sampled in 2008 (area: 2.4-26.3 ha; biomass: 23-183 t/ha; height: 7-21 m). A high-resolution digital terrain model has been used as ground reference during InSAR processing. For the aforementioned plots and stands and if the same acquisition is used for model training and validation, the new model explains 65%-89% of the observed variance, with root-mean-square error (RMSE) of 12%-19% (median: 15%). By fixing two of the three model parameters, accurate biomass estimation can also be done when different acquisitions or different test sites are used for model training and validation, with RMSE of 12%-56% (median: 17%). Compared with a simple scaling model computing biomass from the phase center elevation above ground, the proposed model shows significantly better performance in Remningstorp, as it accounts for the large canopy density variations caused by active management. In Krycklan, the two models show similar performance