Large area forest stem volume mapping using synergy of spaceborne interferometric radar and optical remote sensing: a case study of northeast chin

More than a decade of investigations on the use of the interferometric ERS-1/2 tandem coherence for forest applications have increased the understanding of the behaviour of C-band repeat-pass coherence over forested terrain. It has been shown that under optimal imaging conditions, ERS-1/2 tandem coherence can be used for stem volume retrieval with accuracies in the range of ground surveys. Large-area applications of ERS-1/2 tandem coherence are rare though. One of the main limitations concerning large-area exploitation of the existing ERS-1/2 tandem archives for forest stem volume retrieval is related to the considerable dependence of repeat-pass coherence upon the meteorological (rain, temperature, wind speed) and environmental (soil moisture variations, snow metamorphism) acquisition conditions. Conventional retrieval algorithms require accurate forest inventory data for a dense grid of forest sites to tune models that relate coherence to stem volume to the local conditions. Accurate forest inventory data is, however, a rare commodity that is often not freely available. In this thesis, a fully automated algorithm was developed, based on a synergetic use of the MODIS Vegetation Continuous Field product (Hansen et al., 2002), that allowed the training of the Interferometric Water Cloud Model IWCM (Askne et al., 1997) without further need for forest inventory data. With the new algorithm it was possible to train the IWCM on a frame-by-frame basis and thus to account for the spatial and temporal variability of the meteorological and environmental acquisition conditions. The new algorithm was applied to a multi-seasonal ERS-1/2 tandem dataset covering Northeast China that was acquired between 1995 and 1998 with baselines up to 400 m

Hemiboreaalsete metsade kaardistamine interferomeetrilise tehisava-radari andmetelt

Väitekirja elektrooniline versioon ei sisalda publikatsioone.Käesolev doktoritöö uurib tehisavaradari (SAR) kasutusvõimalusi metsa kõrguse hindamiseks hemiboreaalsete metsade vööndis. Uurimistöö viidi läbi Tartu Üli¬kooli, Tartu Observatooriumi, Aalto Ülikooli, Euroopa Kosmoseagentuuri (ESA) kaugseire keskuse ESRIN ja Reach-U koostöös. Uurimistöös kasutatud satelliidi¬andmed on pärit Saksa Kosmosekeskuse (DLR) kõrglahutusega bistaatilise X-laineala tehisavaradari TanDEM-X satelliidipaarilt. Sagedasti uuenevad satelliidiandmed, nende globaalne katvus ja kõrge ruumi¬line lahutus võimaldavad tehisavaradari abil kaardistada metsi ning nendes toimu¬vaid muutusi suurtel maa-aladel. Radari abil on võimalik saada kõrge lahutusvõimega pilte, mis on tundlikud taimestikule, maapinna karedusele ja dielektrilistele omadustele. Sünkroonis lendava radaripaari samaaegselt tehtud pildid elimineerivad võimalikud ajalised muutused taimestikus ning tänu sellele on radariandmetest võimalik tuletada metsade vertikaalset struktuuri ja kõrgust. Uurimistöös käsitletakse tehisavaradari interferomeetrilise koherentsuse tund¬likkust metsa kõrguse suhtes ning analüüsitakse, millised keskkonna ja klimaati¬lised tingimused ning satelliidi orbiidiga seotud parameetrid mõjutavad radari¬piltidelt erinevate puuliikide kõrguse hindamise täpsust. Lisaks keskendub väitekiri interferomeetrilisele koherentsusele tuginevate mudelite analüüsi¬misele ning nende täpsuse hindamisele operatiivse metsa kõrguse kaardistamise raken-duseks. Vaatluse alla on võetud kolm testala, mis asuvad Soomaa rahvuspargis, Võrtsjärve idakaldal Rannus ja Peipsiveere looduskaitsealal ning katavad kokku 2291 hektarit metsa. 23 TanDEM-X satelliidipildi koherentsuspilte võrreldakse samadel testaladel aerolaserskaneerimise (LiDAR) abil mõõdetud puistute kõrgu¬sega, mis on omakorda jagatud kolme rühma (kuused, männid ja laia¬lehised segametsad). RVoG (Random Volume over Ground) taimekatte mudel ning sellest tule¬tatud lihtsamad pooleempiirilised mudelid sobituvad olemasolevate TanDEM-X koherentsuse ning LiDARi metsa puistute kõrgusandmetega hästi. Töö tule¬mused kinnitavad, et tulevikus on suurte ja erinevatest metsatüüpidest koosne¬vate metsade kõrguse kosmosest kaardistamisel otstarbekas kasutusele võtta esmalt just soovitatud lihtsamad ja universaalsemad mudelid.This thesis presents research in the field of radar remote sensing and contributes to the forest monitoring application development using space-borne synthetic aperture radar (SAR). Satellite data is particularly useful for large-scale forestry applications making high revisit monitoring of the state of forests worldwide possible. The sensitivity of SAR to the dielectric and geometrical properties of the targets, penetration capacity and coherent imaging properties make it a unique tool for mapping and monitoring forest biomes. SAR satellites are also capable of retrieving additional information about the structure of the forest, tree height and biomass estimates as an essential input for monitoring the changes in the carbon stocks. Interferometric SAR (InSAR) is an advanced SAR imaging technique that allows the retrieval of forest parameters while working in nearly all weather conditions, independently of daylight and cloud cover. This research concen¬trates on assessing the impact of different variables affecting hemiboreal forest height estimation from space-borne X-band interferometric SAR coherence data. In particular, the research analyses the changes in coherence dynamics related to seasonal conditions, tree species and imaging properties using a large collection of interferometric SAR images from different seasons over a four-year period. The study is carried out over three test sites in Estonia using the extensive multi-temporal dataset of 23 TanDEM-X images, covering 2291 hectares of forests to describe the relation between the interferometric SAR coherence mag¬nitude and forest parameters. The work demonstrates how the correlation of interferometric coherence and Airborne LiDAR Scanning (ALS)-derived forest height varies for pine and deciduous tree species, for summer (leaf-on) and winter (leaf-off) conditions and for flooded forest floor. A simple semi-empirical modelling approach is proposed as being suitable for wide area forest mapping with limited a priori information under a range of seasonal and environ¬¬mental conditions. A Random Volume over Ground (RVoG) model and three semi-empirical models are compared and validated against a large dataset of coherence magnitude and ALS-measured data over hemiboreal forests in Estonia. The results show that all proposed models perform well in describing the relationship between hemiboreal forest height and interferometric coherence, allowing in future to derive forest stand height with an accuracy suitable for a wide range of applications

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

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

Growing stock volume estimation in temperate forsted areas using a fusion approach with SAR Satellites Imagery

Forest monitoring plays a central role in the context of global warming mitigation and in the assessment of forest resources. To meet these challenges, significant efforts have been made by scientists to develop new feasible remote sensing techniques for the retrieval of forest parameters. However, much work remains to be done in this area, in particular in establishing global assessments of forest biomass. In this context, this Ph.D. Thesis presents a complete methodology for estimating Growing Stock Volume (GSV) in temperate forested areas using a fusion approach based on Synthetic-Aperture Radar (SAR) satellite imagery. The investigations which were performed focused on the Thuringian Forest, which is located in Central Germany. The satellite data used are composed of an extensive set of L-band (ALOS PALSAR) and X-band (TerraSAR-X, TanDEM-X, Cosmo-SkyMed) images, which were acquired in various sensor configurations (acquisition modes, polarisations, incidence angles). The available ground data consists of a forest inventory delivered by the local forest offices. Weather measurements and a LiDAR DEM complete the datasets. The research showed that together with the topography, the forest structure and weather conditions generally limited the sensitivity of the SAR signal to GSV. The best correlations were obtained with ALOS PALSAR (R2 = 0.61) and TanDEM-X (R2 = 0.72) interferometric coherences. These datasets were chosen for the retrieval of GSV in the Thuringian Forest and led with regressions to an root-mean-square error (RMSE) in the range of 100─200 m3ha-1. As a final achievement of this thesis, a methodology for combining the SAR information was developed. Assuming that there are sufficient and adequate remote sensing data, the proposed fusion approach may increase the biomass maps accuracy, their spatial extension and their updated frequency. These characteristics are essential for the future derivation of accurate, global and robust forest biomass maps

Land cover and forest mapping in boreal zone using polarimetric and interferometric SAR data

Remote sensing offers a wide range of instruments suitable to meet the growing need for consistent, timely and cost-effective monitoring of land cover and forested areas. One of the most important instruments is synthetic aperture radar (SAR) technology, where transfer of advanced SAR imaging techniques from mostly experimental small test-area studies to satellites enables improvements in remote assessment of land cover on a global scale. Globally, forests are very suitable for remote sensing applications due to their large dimensions and relatively poor accessibility in distant areas. In this thesis, several methods were developed utilizing Earth observation data collected using such advanced SAR techniques, as well as their application potential was assessed. The focus was on use of SAR polarimetry and SAR interferometry to improve performance and robustness in assessment of land cover and forest properties in the boreal zone. Particular advances were achieved in land cover classification and estimating several key forest variables, such as forest stem volume and forest tree height. Important results reported in this thesis include: improved polarimetric SAR model-based decomposition approach suitable for use in boreal forest at L-band; development and demonstration of normalization method for fully polarimetric SAR mosaics, resulting in improved classification performance and suitable for wide-area mapping purposes; establishing new inversion procedure for robust forest stem volume retrieval from SAR data; developing semi-empirical method and demonstrating potential for soil type separation (mineral soil, peatland) under forested areas with L-band polarimetric SAR; developing and demonstrating methodology for simultaneous retrieval of forest tree height and radiowave attenuation in forest layer from inter-ferometric SAR data, resulting in improved accuracy and more stable estimation of forest tree height

Biomass Representation in Synthetic Aperture Radar Interferometry Data Sets

This work makes an attempt to explain the origin, features and potential applications of the elevation bias of the synthetic aperture radar interferometry (InSAR) datasets over areas covered by vegetation. The rapid development of radar-based remote sensing methods, such as synthetic aperture radar (SAR) and InSAR, has provided an alternative to the photogrammetry and LiDAR for determining the third dimension of topographic surfaces. The InSAR method has proved to be so effective and productive that it allowed, within eleven days of the space shuttle mission, for acquisition of data to develop a three-dimensional model of almost the entire land surface of our planet. This mission is known as the Shuttle Radar Topography Mission (SRTM). Scientists across the geosciences were able to access the great benefits of uniformity, high resolution and the most precise digital elevation model (DEM) of the Earth like never before for their a wide variety of scientific and practical inquiries. Unfortunately, InSAR elevations misrepresent the surface of the Earth in places where there is substantial vegetation cover. This is a systematic error of unknown, yet limited (by the vertical extension of vegetation) magnitude. Up to now, only a limited number of attempts to model this error source have been made. However, none offer a robust remedy, but rather partial or case-based solutions. More work in this area of research is needed as the number of airborne and space-based InSAR elevation models has been steadily increasing over the last few years, despite strong competition from LiDAR and optical methods. From another perspective, however, this elevation bias, termed here as the “biomass impenetrability”, creates a great opportunity to learn about the biomass. This may be achieved due to the fact that the impenetrability can be considered a collective response to a few factors originating in 3D space that encompass the outermost boundaries of vegetation. The biomass, presence in InSAR datasets or simply the biomass impenetrability, is the focus of this research. The report, presented in a sequence of sections, gradually introduces terminology, physical and mathematical fundamentals commonly used in describing the propagation of electromagnetic waves, including the Maxwell equations. The synthetic aperture radar (SAR) and InSAR as active remote sensing methods are summarised. In subsequent steps, the major InSAR data sources and data acquisition systems, past and present, are outlined. Various examples of the InSAR datasets, including the SRTM C- and X-band elevation products and INTERMAP Inc. IFSAR digital terrain/surface models (DTM/DSM), representing diverse test sites in the world are used to demonstrate the presence and/or magnitude of the biomass impenetrability in the context of different types of vegetation – usually forest. Also, results of investigations carried out by selected researchers on the elevation bias in InSAR datasets and their attempts at mathematical modelling are reviewed. In recent years, a few researchers have suggested that the magnitude of the biomass impenetrability is linked to gaps in the vegetation cover. Based on these hints, a mathematical model of the tree and the forest has been developed. Three types of gaps were identified; gaps in the landscape-scale forest areas (Type 1), e.g. forest fire scares and logging areas; a gap between three trees forming a triangle (Type 2), e.g. depending on the shape of tree crowns; and gaps within a tree itself (Type 3). Experiments have demonstrated that Type 1 gaps follow the power-law density distribution function. One of the most useful features of the power-law distributed phenomena is their scale-independent property. This property was also used to model Type 3 gaps (within the tree crown) by assuming that these gaps follow the same distribution as the Type 1 gaps. A hypothesis was formulated regarding the penetration depth of the radar waves within the canopy. It claims that the depth of penetration is simply related to the quantisation level of the radar backscattered signal. A higher level of bits per pixels allows for capturing weaker signals arriving from the lower levels of the tree crown. Assuming certain generic and simplified shapes of tree crowns including cone, paraboloid, sphere and spherical cap, it was possible to model analytically Type 2 gaps. The Monte Carlo simulation method was used to investigate relationships between the impenetrability and various configurations of a modelled forest. One of the most important findings is that impenetrability is largely explainable by the gaps between trees. A much less important role is played by the penetrability into the crown cover. Another important finding is that the impenetrability strongly correlates with the vegetation density. Using this feature, a method for vegetation density mapping called the mean maximum impenetrability (MMI) method is proposed. Unlike the traditional methods of forest inventories, the MMI method allows for a much more realistic inventory of vegetation cover, because it is able to capture an in situ or current situation on the ground, but not for areas that are nominally classified as a “forest-to-be”. The MMI method also allows for the mapping of landscape variation in the forest or vegetation density, which is a novel and exciting feature of the new 3D remote sensing (3DRS) technique. Besides the inventory-type applications, the MMI method can be used as a forest change detection method. For maximum effectiveness of the MMI method, an object-based change detection approach is preferred. A minimum requirement for the MMI method is a time-lapsed reference dataset in the form, for example, of an existing forest map of the area of interest, or a vegetation density map prepared using InSAR datasets. Preliminary tests aimed at finding a degree of correlation between the impenetrability and other types of passive and active remote sensing data sources, including TerraSAR-X, NDVI and PALSAR, proved that the method most sensitive to vegetation density was the Japanese PALSAR - L-band SAR system. Unfortunately, PALSAR backscattered signals become very noisy for impenetrability below 15 m. This means that PALSAR has severe limitations for low loadings of the biomass per unit area. The proposed applications of the InSAR data will remain indispensable wherever cloud cover obscures the sky in a persistent manner, which makes suitable optical data acquisition extremely time-consuming or nearly impossible. A limitation of the MMI method is due to the fact that the impenetrability is calculated using a reference DTM, which must be available beforehand. In many countries around the world, appropriate quality DTMs are still unavailable. A possible solution to this obstacle is to use a DEM that was derived using P-band InSAR elevations or LiDAR. It must be noted, however, that in many cases, two InSAR datasets separated by time of the same area are sufficient for forest change detection or similar applications