Utilization of bistatic TanDEM-X data to derive land cover information

Forests have significance as carbon sink in climate change. Therefore, it is of high importance to track land use changes as well as to estimate the state as carbon sink. This is useful for sustainable forest management, land use planning, carbon modelling, and support to implement international initiatives like REDD+ (Reducing Emissions from Deforestation and Degradation). A combination of field measurements and remote sensing seems most suitable to monitor forests. Radar sensors are considered as high potential due to the weather and daytime independence. TanDEM-X is a interferometric SAR (synthetic aperture radar) mission in space and can be used for land use monitoring as well as estimation of biophysical parameters. TanDEM-X is a X-band system resulting in low penetration depth into the forest canopy. Interferometric information can be useful, whereas the low penetration can be considered as an advantage. The interferometric height is assumable as canopy height, which is correlated with forest biomass. Furthermore, the interferometric coherence is mainly governed by volume decorrelation, whereas temporal decorrelation is minimized. This information can be valuable for quantitative estimations and land use monitoring. The interferometric coherence improved results in comparison to land use classifications without coherence of about 10% (75% vs. 85%). Especially the differentiation between forest classes profited from coherence. The coherence correlated with aboveground biomass in a R² of about 0.5 and resulted in a root mean square error (RSME) of 14%. The interferometric height achieved an even higher correlation with the biomass (R²=0.68) resulting in cross-validated RMSE of 7.5%. These results indicated that TanDEM-X can be considered as valuable and consistent data source for forest monitoring. Especially interferometric information seemed suitable for biomass estimation

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

Forest Height Inversion by Combining Single-Baseline TanDEM-X InSAR Data with External DTM Data

Forest canopy height estimation is essential for forest management and biomass estimation. In this study, we aimed to evaluate the capacity of TanDEM-X interferometric synthetic aperture radar (InSAR) data to estimate canopy height with the assistance of an external digital terrain model (DTM). A ground-to-volume ratio estimation model was proposed so that the canopy height could be precisely estimated from the random-volume-over-ground (RVoG) model. We also refined the RVoG inversion process with the relationship between the estimated penetration depth (PD) and the phase center height (PCH). The proposed method was tested by TanDEM-X InSAR data acquired over relatively homogenous coniferous forests (Teruel test site) and coniferous as well as broadleaved forests (La Rioja test site) in Spain. Comparing the TanDEM-X-derived height with the LiDAR-derived height at plots of size 50 m × 50 m, the root-mean-square error (RMSE) was 1.71 m (R2 = 0.88) in coniferous forests of Teruel and 1.97 m (R2 = 0.90) in La Rioja. To demonstrate the advantage of the proposed method, existing methods based on ignoring ground scattering contribution, fixing extinction, and assisting with simulated spaceborne LiDAR data were compared. The impacts of penetration and terrain slope on the RVoG inversion were also evaluated. The results show that when a DTM is available, the proposed method has the optimal performance on forest height estimation.This work was supported in part by the National Natural Science Foundation of China under Grant 41820104005, Grant 42030112, and Grant 41904004, Hunan Natural Science Foundation under Grant 2021JJ30808, and in part by the Spanish Ministry of Science and Innovation, Agencia Estatal de Investigacion, under Projects PID2020-117303GB-C22/AEI/10.13039/501100011033 and PROWARM (PID2020-118444GA-I00/AEI/10.13039/501100011033)

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

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

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

Applications of SAR Interferometry in Earth and Environmental Science Research

This paper provides a review of the progress in regard to the InSAR remote sensing technique and its applications in earth and environmental sciences, especially in the past decade. Basic principles, factors, limits, InSAR sensors, available software packages for the generation of InSAR interferograms were summarized to support future applications. Emphasis was placed on the applications of InSAR in seismology, volcanology, land subsidence/uplift, landslide, glaciology, hydrology, and forestry sciences. It ends with a discussion of future research directions

Assessment of L-Band SAOCOM InSAR coherence and its comparison with C-Band: A case study over managed forests in Argentina

The objective of this work is to analyze the behavior of short temporal baseline interferomet ric coherence in forested areas for L-band spaceborne SAR data. Hence, an exploratory assessment of the impacts of temporal and spatial baselines on coherence, with emphasis on how these effects vary between SAOCOM-1 L-band and Sentinel-1 C-band data is presented. The interferometric coherence is analyzed according to different imaging parameters. In the case of SAOCOM-1, the impacts of the variation of the incidence angle and the ascending and descending orbits over forested areas are also assessed. Finally, short-term 8-day interferometric coherence maps derived from SAOCOM-1 are especially addressed, since this is the first L-band spaceborne mission that allows us to acquire SAR images with such a short temporal span. The analysis is reported over two forest-production areas in Argentina, one of which is part of the most important region in terms of forest plantations at the national level. In the case of SAOCOM, interferometric configurations are characterized by a lack of control on the spatial baseline, so a zero-baseline orbital tube cannot be guaranteed. Nevertheless, this spatial baseline variability is crucial to exploit volume decorrelation for forest monitoring. The results from this exploratory analysis demonstrates that SAOCOM-1 short temporal baseline interferograms, 8 to 16 days, must be considered in order to mitigate temporal decorrelation effects and to be able to experiment with different spatial baseline configurations, in order to allow appropriate forest monitoring.This research was funded by the project INTERACT PID2020-114623RB-C32 funded by the Spanish MCIN /AEI /10.13039 /501100011033.Peer ReviewedPostprint (published version

Puistute takseertunnuste hindamine aerolidari mõõtmisandmete põhjal hemiboreaalsetes metsades

A Thesis for applying for the degree of Doctor of Philosophy in Forestry.Forest management and planning requires up-to-date data, which commonly is acquired using field experts and ground measurements. Nowadays, more and more of data about forest stands is measured using remotely sensing methods. Most common methods include aerial photography and laser scanning from airplanes, also spectral measurements from satellites or even drone images and applications. This doctoral thesis focuses on developing applications and methods for utilising the airborne laser scanning (ALS) data that is freely available for the whole Estonia. The ALS measurements are carried out by the Estonian Land Board on a routine basis twice a year – in spring and summer. The first variable that was studied in this thesis was forest height. Based on the thesis, the most reliable method for forest height assessment was using the ALS point-cloud 80th height percentile (HP80). The small circular plot (radius of 15…30 m) and stand based studies showed high correlations with the field-measured forest heights and with great confidence it can be said, that ALS-based forest height estimations are close or even with higher accuracy, than field inspected. The second studied variable was standing wood volume. The ALS-based methods and models that were developed throughout this thesis used the idea, that standing wood volume is based on forest height and density. For this the HP80 and a threshold-based point count ratio was used (canopy cover - CC). ALS-based CC (CCALS) estimates were studied and compared with digital hemispherical photo based measurements. The results showed similar errors as were shown in other similar studies, with around 10-15% root mean square error (RMSE). The strongest correlation was shown using all echoes above a 1.3 metre threshold. Combining the CCALS and HP80 showed standing wood volume estimates with a similar error as we would receive from field measurements (<20%). The freely available multitemporal ALS data showed promising results for forest height growth monitoring and detecting small-scale disturbances. CCALS was shown to have strong predictive value, when compared with a four year difference in thinned and unthinned stands. The nation-wide ALS data can also be combined with forest height predictions from satellites, providing a faster update compared to the ALS data. Promising results were shown using the interferometric synthetic aperture radar (InSAR). Stand species maps generated using self-learning algorithms and satellite based spectral data can be used for developing species specific models of standing wood volume prediction. By combining these different datasets we can construct a nation-wide forest resource to help make better decisions for forest management and targeting fieldwork.Metsades majandamisotsuste langetamiseks ja metsamajanduslike tööde planeerimiseks on metsaomanikel vaja andmeid. Harjumuspäraselt on andmete kogumiseks tehtud metsas maapealseid mõõtmisi. Viimastel aastakümnetel on metsade inventeerimiseks üha enam aga kasutatud mittekontaktseid mõõtmisi - lennukitelt tehtavad aerofotosid, laserskaneerimist, satelliitidelt tehtavaid kiirgusmõõtmisi või viimastel aastatel ka droonidelt tehtud pilte. Antud doktoritöö on võtnud fookusesse aerolaserskaneerimise (ALS) andmete põhjal Eesti metsadesse sobilike rakenduste väljatöötamise. ALS mõõtmisi teeb Eesti Maa-amet rutiinsete lendude käigus kaks korda aastas, nii kevadel kui ka suvel. Aastast 2008 alustatud mõõtmiste tulemusel on Eesti üks väheseid riike maailmas, kus on vabalt kasutada mitmekordselt kogu riiki kattev ALS andmestik. Doktoritöö tulemusel töötati välja metsa kõrguse hindamiseks sobilikud meetodid, kasutades selleks punktipilvede kõrgusprotsentiile. Tugevamaid seoseid metsas proovitükkidel mõõdetud kõrgustega näitas punktipilve 80-protsentiil (HP80) ja uuringute põhjal võib väita, et metsa kõrguse määramine suvistelt aerolidari andmetelt on ligilähedane täpsustele, mida saadakse metsas kohapeal mõõtes. Teine oluline tunnus, mida metsade majandamise planeerimisel silmas peetakse, on kasvava metsa tagavara. Teadustöö põhjal töötati välja mudelite kujud ja metoodika, mille abil prognoositud tagavara oli sarnase veapiiriga, mis on lubatud metsas hinnanguid tegevatele taksaatoritele (<20%). Väljatöötatud ALS-põhine mudeli kuju järgib loogikat, et metsa tagavara on otseselt seotud mõõdetud kõrguse ja metsa tihedusega. Tihenduse hindamiseks aerolidari andmetelt kasutatakse nivoopõhist punktide suhtearvu ehk nn katvushinnangut (CCALS). Katvushinnangu täpsuse valideerimiseks ja tihedas metsas sobiva prognoosimeetodi väljatöötamiseks tehti välimõõtmisi kasutades poolsfäärikaameraid. Poolsfääripiltide põhjal tehtud valideerimise tulemused andsid sarnaseid veahinnanguid, mida on ka varasemates teadusuuringutes esitletud (RMSE = 10…15%). Kahe sarnasest fenoloogilisest perioodist ALS andmestiku lahutamisel uuriti ka muutuste tuvastamise võimalikkust. Uuringud andsid paljulubavaid tulemusi metsade kõrguskasvu hindamiseks ja CCALS osutus ka oluliseks tunnuseks väiksemate häiringute, nagu näiteks harvendusraie, tuvastamiseks. Kogu riiki katva ALS andmestiku kombineerimisel erinevate satelliitandmetega või näiteks spektraalsete mõõtmiste põhjal tehtud puistu liigiliste koosseisu kaartidega on võimalik antud töös välja pakutud meetodite abil anda igal aastal kogu Eesti metsaressursside ülevaade. Samuti on võimalik koostada vaid kaugseirevahendeid ja proovitükkidel lähendatud mudeleid kasutades eraldiste põhised takseerkirjeldused, mida siis taksaatorid saavad näiteks kasutada oma välitööde kavandamisel.  Publication of this thesis is supported by the Estonian University of Life Sciences