A Macroecological Analysis of SERA Derived Forest Heights and Implications for Forest Volume Remote Sensing

Individual trees have been shown to exhibit strong relationships between DBH, height and volume. Often such studies are cited as justification for forest volume or standing biomass estimation through remote sensing. With resolution of common satellite remote sensing systems generally too low to resolve individuals, and a need for larger coverage, these systems rely on descriptive heights, which account for tree collections in forests. For remote sensing and allometric applications, this height is not entirely understood in terms of its location. Here, a forest growth model (SERA) analyzes forest canopy height relationships with forest wood volume. Maximum height, mean, H100, and Lorey's height are examined for variability under plant number density, resource and species. Our findings, shown to be allometrically consistent with empirical measurements for forested communities world-wide, are analyzed for implications to forest remote sensing techniques such as LiDAR and RADAR. Traditional forestry measures of maximum height, and to a lesser extent H100 and Lorey's, exhibit little consistent correlation with forest volume across modeled conditions. The implication is that using forest height to infer volume or biomass from remote sensing requires species and community behavioral information to infer accurate estimates using height alone. SERA predicts mean height to provide the most consistent relationship with volume of the height classifications studied and overall across forest variations. This prediction agrees with empirical data collected from conifer and angiosperm forests with plant densities ranging between 102–106 plants/hectare and heights 6–49 m. Height classifications investigated are potentially linked to radar scattering centers with implications for allometry. These findings may be used to advance forest biomass estimation accuracy through remote sensing. Furthermore, Lorey's height with its specific relationship to remote sensing physics is recommended as a more universal indicator of volume when using remote sensing than achieved using either maximum height or H100

Comparison between TanDEM-X- and ALS-based estimation of aboveground biomass and tree height in boreal forests

Interferometric Synthetic Aperture Radar (InSAR) data from TerraSAR-X add-on for Digital Elevation Measurement (TanDEM-X) were used to estimate aboveground biomass (AGB) and tree height with linear regression models. These were compared to models based on airborne laser scanning (ALS) data at two Swedish boreal forest test sites, Krycklan (64 degrees N19 degrees E) and Remningstorp (58 degrees N13 degrees E). The predictions were validated using field data at the stand-level (0.5-26.1 ha) and at the plot-level (10 m radius). Additionally, the ALS metrics percentile 99 (p99) and vegetation ratio, commonly used to estimate AGB and tree height, were estimated in order to investigate the feasibility of replacing ALS data with TanDEM-X InSAR data. Both AGB and tree height could be estimated with about the same accuracy at the stand-level from both TanDEM-X- and ALS-based data. The AGB was estimated with 17.2% and 14.6% root mean square error (RMSE) and the tree height with 7.6% and 4.1% RMSE from TanDEM-X data at the stand-level at the two test sites Krycklan and Remningstorp. The Pearson correlation coefficients between the TanDEM-X height and the ALS height p99 were r=.98 and r=.95 at the two test sites. The TanDEM-X height contains information related to both tree height and forest density, which was validated from several estimation models

Relationship between Lidar-Derived Canopy Densities and the Scattering Phase Center of High-Resolution TanDEM-X Data

Abstract: The estimation of forestry parameters is essential to understanding the three-dimensional structure of forests. In this respect, the potential of X-band synthetic aperture radar (SAR) has been recognized for years. Many studies have been conducted on deriving tree heights with SAR data, but few have paid attention to the effects of the canopy structure. Canopy density plays an important role since it provides information about the vertical distribution of dominant scatterers in the forest. In this study, the position of the scattering phase center (SPC) of interferometric X-band SAR data is investigated with regard to the densest vegetation layer in a deciduous and coniferous forest in Germany by applying a canopy density index from high-resolution airborne laser scanning data. Two different methods defining the densest layer are introduced and compared with the position of the TanDEM-X SPC. The results indicate that the position of the SPC often coincides with the densest layer, with mean differences ranging from −1.6 m to +0.7 m in the deciduous forest and +1.9 m in the coniferous forest. Regarding relative tree heights, the SAR signal on average penetrates up to 15% (3.4 m) of the average tree height in the coniferous forest. In the deciduous forest, the difference increases to 18% (6.2 m) during summer and 24% (8.2 m) during winter. These findings highlight the importance of considering not only tree height but also canopy density when delineating SAR-based forest heights. The vertical structure of the canopy influences the position of the SPC, and incorporating canopy density can improve the accuracy of SAR-derived forest height estimations

Microwave models of snow characteristics for remote sensing

One of the key problems of microwave remote sensing is the development of theoretical microwave models for terrain such as soil, vegetation, snow, forest, etc., due to the complexity of modeling of microwave interaction with the terrain. In this thesis this problem is approached from the new point of view of both empirical models and rigorous theoretical models. New information concerning radar remote sensing of snow-covered terrain and permittivity of snow has been produced. A C-band semi-empirical backscattering model is presented for the forest-snow-ground system. The effective permittivity of random media such as snow, vegetation canopy, soil, etc., describes microwave propagation and attenuation in the media and is a very important parameter in modeling of microwave interaction with the terrain. Good permittivity models are needed in microwave emission and scattering models of terrain. In this thesis, the strong fluctuation theory is applied to calculate the effective permittivity of wet snow. Numerical results for the effective permittivity of wet snow are illustrated. The results are compared with the semi-empirical and the theoretical models. A comparison with experimental data at 6, 18 and 37 GHz is also presented. The results indicate that the model presented in this work gives reasonably good accuracy for calculating the effective permittivity of wet snow. Microwave emission and scattering theoretical models of wet snow are developed based on the radiative transfer and strong fluctuation theory. It is shown that the models agree with the experimental data.reviewe

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

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

Spatial and temporal statistics of SAR and InSAR observations for providing indicators of tropical forest structural changes due to forest disturbance

Tropical forests are extremely important ecosystems which play a substantial role in the global carbon budget and are increasingly dominated by anthropogenic disturbance through deforestation and forest degradation, contributing to emissions of greenhouse gases to the atmosphere. There is an urgent need for forest monitoring over extensive and inaccessible tropical forest which can be best accomplished using spaceborne satellite data. Currently, two key processes are extremely challenging to monitor: forest degradation and post-disturbance re-growth. The thesis work focuses on these key processes by considering change indicators derived from radar remote sensing signal that arise from changes in forest structure. The problem is tackled by exploiting spaceborne Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) observations, which can provide forest structural information while simultaneously being able to collect data independently of cloud cover, haze and daylight conditions which is a great advantage over the tropics. The main principle of the work is that a connection can be established between the forest structure distribution in space and signal variation (spatial statistics) within backscatter and Digital Surface Models (DSMs) provided by SAR. In turn, forest structure spatial characteristics and changes are used to map forest condition (intact or degraded) or disturbance. The innovative approach focuses on looking for textural patterns (and their changes) in radar observations, then connecting these patterns to the forest state through supporting evidence from expert knowledge and auxiliary remote sensing observations (e.g. high resolution optical, aerial photography or LiDAR). These patterns are descriptors of the forest structural characteristics in a statistical sense, but are not estimates of physical properties, such as above-ground biomass or canopy height. The thesis tests and develops methods using novel remote sensing technology (e.g. single-pass spaceborne InSAR) and modern image statistical analysis methods (wavelet-based space-scale analysis). The work is developed on an experimental basis and articulated in three test cases, each addressing a particular observational setting, analytical method and thematic context. The first paper deals with textural backscatter patterns (C-band ENVISAT ASAR and L-band ALOS PALSAR) in semi-deciduous closed forest in Cameroon. Analysis concludes that intact forest and degraded forest (arising from selective logging) are significantly different based on canopy structural properties when measured by wavelet based space-scale analysis. In this case, C-band data are more effective than longer wavelength L-band data. Such a result could be explained by the lower wave penetration into the forest volume at shorter wavelength, with the mechanism driving the differences between the two forest states arising from upper canopy heterogeneity. In the second paper, wavelet based space-scale analysis is also used to provide information on upper canopy structure. A DSM derived from TanDEM-X acquired in 2014 was used to discriminate primary lowland Dipterocarp forest, secondary forest, mixed-scrub and grassland in the Sungai Wain Protection Forest (East Kalimantan, Indonesian Borneo) which was affected by the 1997/1998 El Niño Southern Oscillation (ENSO). The Jeffries- Matusita separability of wavelet spectral measures of InSAR DSMs between primary and secondary forest was in some cases comparable to results achieved by high resolution LiDAR data. The third test case introduces a temporal component, with change detection aimed at detecting forest structure changes provided by differencing TanDEM-X DSMs acquired at two dates separated by one year (2012-2013) in the Republic of Congo. The method enables cancelling out the component due to terrain elevation which is constant between the two dates, and therefore the signal related to the forest structure change is provided. Object-based change detection successfully mapped a gradient of forest volume loss (deforestation/forest degradation) and forest volume gain (post-disturbance re-growth). Results indicate that the combination of InSAR observations and wavelet based space-scale analysis is the most promising way to measure differences in forest structure arising from forest fires. Equally, the process of forest degradation due to shifting cultivation and post-disturbance re-growth can be best detected using multiple InSAR observations. From the experiments conducted, single-pass InSAR appears to be the most promising remote sensing technology to detect forest structure changes, as it provides three-dimensional information and with no temporal decorrelation. This type of information is not available in optical remote sensing and only partially available (through a 2D mapping) in SAR backscatter. It is advised that future research or operational endeavours aimed at mapping and monitoring forest degradation/regrowth should take advantage of the only currently available high resolution spaceborne single-pass InSAR mission (TanDEM-X). Moreover, the results contribute to increase knowledge related to the role of SAR and InSAR for monitoring degraded forest and tracking the process of forest degradation which is a priority but still highly challenging to detect. In the future the techniques developed in the thesis work could be used to some extent to support REDD+ initiatives

Establishing the sensitivity of Synthetic Aperture Radar to above-ground biomass in wooded savannas

Radar for biomass estimation has been widely investigated for temperate, boreal and tropical forests, yet tropical savanna woodlands, which generally form non-continuous cover canopies or sparse woodlands, have been largely neglected in biomass studies. This thesis evaluates the capability of Synthetic Aperture Radar (SAR) for estimating the above-ground biomass of the woody vegetation in a savanna in Belize, Central America. This is achieved by evaluating (i) polarimetric Synthetic Aperture Radar (SAR) backscatter and (ii) single-pass shortwave interferometric SAR (InSAR) as indicators of above-ground biomass. Specifically, the effect on SAR backscatter of woody vegetation structure such as canopy cover, basal area, vegetation height and above-ground biomass is evaluated. Since vegetation height is often correlated to above-ground biomass, the effectiveness of vegetation height retrieval from InSAR is evaluated as an indicator of above-ground biomass. The study area, situated in Belize, is representative of Central American savannas. Radar data used are AIRSAR fully polarimetric L- and P-band SAR, and AIRSAR C-band InSAR, Intermap Technologies STAR-3i X-band InSAR, and Shuttle Radar Topography Mission (SRTM) C-band InSAR. The field data comprise accurately georeferenced three-dimensional measurements for 1,133 trees and shrubs and 75 palmetto clumps and thickets in a transect of 800 m x 60 m which spans the main savanna vegetation strata of the study area. An additional 2,464 ground points were observed. Results show that savanna woodlands present a challenge for radar remote sensing methods due to the sparse and heterogeneous nature of savanna woodlands. Long-wave SAR backscatter is dominated not only by high biomass areas, but also by areas of leafy palmetto which have low vegetative biomass. Retrieved woodland canopy heights from X- and C-band InSAR are indicative of the general patterns of tree height, although retrieved heights are underestimated. The amount of underestimation is variable across the different canopy conditions. Of these two methods, the shortwave InSAR data give a better indication of the spatial distribution of the above-ground biomass of the woody vegetation in the savannas than SAR backscatter. These results have implications for new and planned future global biomass estimation missions, such as ALOS PALSAR, ESA’s planned P-band BIOMASS and TanDEM-X. Without appropriate mediation, SAR backscatter methods might overestimate above-ground biomass of the woody vegetation of savannas while InSAR height retrieval methods might underestimate biomass estimates. Some possible mediating approaches are discussed

Radar backscatter modelling of forests using a macroecological approach

This thesis provides a new explanation for the behaviour of radar backscatter of forests using vegetation structure models from the field of macroecology. The forests modelled in this work are produced using allometry-based ecological models with backscatter derived from the parameterisation of a radiative transfer model. This work is produced as a series of papers, each portraying the importance of macroecology in defining the forest radar response. Each contribution does so by incorporating structural and dynamic effects of forest growth using one of two allometric models to expose variations in backscatter as a response to vertical and horizontal forest profiles. The major findings of these studies concern the origin of backscatter saturation effects from forest SAR surveys. In each work the importance of transition from Rayleigh to Optical scattering, combined with the scaling effects of forest structure, is emphasised. These findings are administered through evidence including the transition’s emergence as the region of dominant backscatter in a vertical profile (according to a dominant canopy scattering layer), also through the existence of a two trend backscatter relationship with volume in the shape of the typical “saturation curve” (in the absence of additional attenuating factors). The importance of scattering regime change is also demonstrated through the relationships with volume, basal area and thinning. This work’s findings are reinforced by the examination of the relationships between forest height and volume, as collective values, providing evidence to suggest the non-uniqueness of volume-toheight relationships. Each of the studies refer to growing forest communities not single trees, so that unlike typical studies of radar remote sensing of forests the impact of the macroecological structural aspects are more explicit. This study emphasises the importance of the overall forest structure in producing SAR backscatter and how backscatter is not solely influenced by electrical properties of scatteres or the singular aspects of a tree but also by the collective forest parameters defining a dynamically changing forest