Lorey's height regression for ICESAT-GLAS waveforms in hyrcanian deciduous forests of Iran

IGARSS 2015, Milan, ITA, 26-/07/2015 - 31/07/2015International audienceSince Lidar technology provides the most direct measurements of 3D of phenomena, it plays a critical role in a variety of applications. Forest canopy height as a main factor in forest biomass estimation is costly and time consuming to be measured on the ground. This study aims to estimate Lorey's height “Hlorey” using GLAS data based on regression models. Different metrics like waveform extent “Wext”, trail-edge extent “Htrail” and lead-edge extent “Hlead” were extracted from waveforms and a terrain index “TI” was also calculated using a digital elevation model. Hlorey estimated using multiple regression models were compared to field measurements data. A 5-fold cross validation method was used to validate the results. Best model with lowest AIC (297.440) was resulted using combination of Wext and TI (R_a^2=0.72; RMSE= 5.04m). The results show capability of ICESat-GLAS to estimate Lorey's height in sloped area with a simple regression model. It is prospected to reach better result using other statistical methods and also improvement of processing techniques for LiDAR waveforms in the case of sloped terrai

Regional scale rain-forest height mapping using regression-kriging of spaceborneand airborne lidar data: application on French Guiana

IGARSS 2015, Milan, ITA, 26-/07/2015 - 31/07/2015International audienceLiDAR remote sensing has been shown to be a good technique for the estimation of forest parameters such as canopy heights and aboveground biomass. Whilst airborne LiDAR data are in general very dense but only available over small areas due to the cost of their acquisition, spaceborne LiDAR data acquired from the Geoscience Laser Altimeter System (GLAS) have a coarser acquisition density associated with a global cover. It is therefore valuable to analyze the integration relevance of canopy heights estimated from LiDAR sensors with ancillary data such as geological, meteorological, and phenological variables in order to propose a forest canopy height map with good precision and high spatial resolution.In this study, canopy heights extracted from both airborne and spaceborne LiDAR, were first extrapolated from available environmental data. The estimated canopy height maps using random forest (RF) regression from the airborne or GLAS calibration datasets showed similar precisions (RMSE better than 6.5 m). In order to improve the precision of the canopy height estimates regression-kriging (kriging of RF regression residuals) was used. Results indicated an improvement in the RMSE (decrease from 6.5 to 4.2 m) for the regression-kriging maps from the GLAS dataset, and from 5.8 to 1.8 m for the regression-kriging map from the airborne LiDAR dataset

Capability of GLAS/ICESat data to estimate forest canopy height and volume in mountainous forests of Iran

International audienceThe importance of measuring biophysical properties of forest for ecosystem health monitoring and forest management encourages researchers to find precise, yet low cost methods especially in mountainous and large area. In the present study Geoscience Laser Altimeter System (GLAS) on board ICESat was used to estimate three biophysical characteristics of forests located in north of Iran: 1) maximum canopy height (Hmax), 2) Lorey's height (HLorey), and 3) Forest volume (V). A large number of Multiple Linear Regressions (MLR) and also Random Forest (RF) regressions were developed using different set of variables: waveform metrics, Principal Components (PCs) produced from Principal Component Analysis (PCA) and Wavelet Coefficients (WCs) generated from wavelet transformation. To validate and compare different models, statistical criteria were calculated based on a five-fold cross validation. The best model concerning the maximum canopy height was an MLR with an RMSE of 5.0 m which combined two metrics extracted from waveforms (waveform extent "Wext" and height at 50% of waveform energy "H50"), and one from the Digital Elevation Model (Terrain Index: TI). The mean absolute error (MAPE) of maximum canopy height estimates is about 16.4%. For Lorey's height, a simple MLR model including two metrics (Wext and TI) represents the highest performance (RMSE=5.1 m, MAPE=24.0%). Totally, MLR models showed better performance rather than RF models, and accuracy of height estimations using waveform metrics was greater than those based on PCs or WCs. Concerning forest volume, employing regression models to estimate volume directly from GLAS data led to a better result (RMSE=128.8 m3/ha) rather than volume-HLorey relationship (RMSE=167.8 m3/ha)

Coupling potential of ICESat/GLAS and SRTM for the discrimination of forest landscape types in French Guiana

The Shuttle Radar Topography Mission (SRTM) has produced the most accurate nearly global elevation dataset to date. Over vegetated areas, the measured SRTM elevations are the result of a complex interaction between radar waves and tree crowns. In this study, waveforms acquired by the Geoscience Laser Altimeter System (GLAS) were combined with SRTM elevations to discriminate the five forest landscape types (LTs) in French Guiana. Two differences were calculated: (1) penetration depth, defined as the GLAS highest elevations minus the SRTM elevations, and (2) the GLAS centroid elevations minus the SRTM elevations. The results show that these differences were similar for the five LTs, and they increased as a function of the GLAS canopy height and of the SRTM roughness index. Next, a Random Forest (RF) classifier was used to analyze the coupling potential of GLAS and SRTM in the discrimination of forest landscape types in French Guiana. The parameters used in the RF classification were the GLAS canopy height, the SRTM roughness index, the difference between the GLAS highest elevations and the SRTM elevations and the difference between the GLAS centroid elevations and the SRTM elevations. Discrimination of the five forest landscape types in French Guiana was possible, with an overall classification accuracy of 81.3% and a kappa coefficient of 0.75. All forest LTs were well classified with an accuracy varying from 78.4% to 97.5%. Finally, differences of near coincident GLAS waveforms, one from the wet season and one from the dry season, were analyzed. The results showed that the open forest LT (LT12), in some locations, contains trees that lose leaves during the dry season. These trees allow LT12 to be easily discriminated from the other LTs that retain their leaves using the following three criteria: (1) difference between the GLAS centroid elevations and the SRTM elevations, (2) ratio of top energy in the wet season to top energy in the dry season, or (3) ratio of ground energy in the wet season to ground energy in the dry season

Aboveground biomass mapping in French Guiana by combining remote sensing, forest inventories and environmental data

International audienceMapping forest aboveground biomass (AGB) has become an important task, particularly for the reporting of carbon stocks and changes. AGB can be mapped using synthetic aperture radar data (SAR) or passive optical data. However, these data are insensitive to high AGB levels (>150 Mg/ha, and >300 Mg/ha for P-band), which are commonly found in tropical forests. Studies have mapped the rough variations in AGB by combining optical and environmental data at regional and global scales. Nevertheless, these maps cannot represent local variations in AGB in tropical forests. In this paper, we hypothesize that the problem of misrepresenting local variations in AGB and AGB estimation with good precision occurs because of both methodological limits (signal saturation or dilution bias) and a lack of adequate calibration data in this range of AGB values. We test this hypothesis by developing a calibrated regression model to predict variations in high AGB values (mean >300 Mg/ha) in French Guiana by a methodological approach for spatial extrapolation with data from the optical geoscience laser altimeter system (GLAS), forest inventories, radar, optics, and environmental variables for spatial inter-and extrapolation. Given their higher point count, GLAS data allow a wider coverage of AGB values. We find that the metrics from GLAS footprints are correlated with field AGB estimations (R 2 =0.54, RMSE=48.3 Mg/ha) with no bias for high values. First, predictive models, including remote-sensing, environmental variables and spatial correlation functions, allow us to obtain "wall-to-wall" AGB maps over French Guiana with an RMSE for the in situ AGB estimates of ~51 Mg/ha and R²=0.48 at a 1-km grid size. We conclude that a calibrated regression model based on GLAS with dependent environmental data can produce good AGB predictions even for high AGB values if the calibration data fit the AGB range. We also demonstrate that small temporal and spatial mismatches between field data and GLAS footprints are not a problem for regional and global calibrated regression models because field data aim to predict large and deep tendencies in AGB variations from environmental gradients and do not aim to represent high but stochastic and temporally limited variations from forest dynamics. Thus, we advocate including a greater variety of data, even if less precise and shifted, to better represent high AGB values in global models and to improve the fitting of these models for high values

Regional scale rain-forest height mapping using regression-kriging of spaceborne and airborne lidar data : application on French Guiana

LiDAR data has been successfully used to estimate forest parameters such as canopy heights and biomass. Major limitation of LiDAR systems (airborne and spaceborne) arises from their limited spatial coverage. In this study, we present a technique for canopy height mapping using airborne and spaceborne LiDAR data (from the Geoscience Laser Altimeter System (GLAS)). First, canopy heights extracted from both airborne and spaceborne LiDAR were extrapolated from available environmental data. The estimated canopy height maps using Random Forest (RF) regression from airborne or GLAS calibration datasets showed similar precisions (~6 m). To improve the precision of canopy height estimates, regression-kriging was used. Results indicated an improvement in terms of root mean square error (RMSE, from 6.5 to 4.2 m) using the GLAS dataset, and from 5.8 to 1.8 m using the airborne LiDAR dataset. Finally, in order to investigate the impact of the spatial sampling of future LiDAR missions on canopy height estimates precision, six subsets were derived from the initial airborne LiDAR dataset. Results indicated that using the regression-kriging approach a precision of 1.8 m on the canopy height map was achievable with a flight line spacing of 5 km. This precision decreased to 4.8 m for flight line spacing of 50 km

Estimating aboveground woody biomass change in Kalahari woodland: combining field, radar, and optical data sets

Maps that accurately quantify aboveground vegetation biomass (AGB) are essential for ecosystem monitoring and conservation. Throughout Namibia, four vegetation change processes are widespread, namely, deforestation, woodland degradation, the encroachment of the herbaceous and grassy layers by woody strata (woody thickening), and woodland regrowth. All of these vegetation change processes affect a range of key ecosystem services, yet their spatial and temporal dynamics and contributions to AGB change remain poorly understood. This study quantifies AGB associated with the different vegetation change processes over an 8-year period, for a region of Kalahari woodland savannah in northern Namibia. Using data from 101 forest inventory plots collected during two field campaigns (2014–2015), we model AGB as a function of the Advanced Land Observing Satellite Phased Array L-band synthetic aperture radar (PALSAR and PALSAR-2) and dry season Landsat vegetation index composites, for two periods (2007 and 2015). Differences in AGB between 2007 and 2015 were assessed and validated using independent data, and changes in AGB for the main vegetation processes are quantified for the whole study area (75,501 km2). We find that woodland degradation and woody thickening contributed a change in AGB of −14.3 and 2.5 Tg over 14% and 3.5% of the study area, respectively. Deforestation and regrowth contributed a smaller portion of AGB change, i.e. −1.9 and 0.2 Tg over 1.3% and 0.2% of the study area, respectively

Testing Different Methods of Forest Height and Aboveground Biomass Estimations From ICESat/GLAS Data in Eucalyptus Plantations in Brazil

International audienceThe Geoscience Laser Altimeter System (GLAS) has provided a useful dataset for estimating forest heights in many areas of the globe. Most of the studies on GLAS waveforms have fo- cused on natural forests and only a few were conducted over forest plantations. This work set out to estimate the stand-scale dominant height and aboveground biomass of intensively managed Euca- lyptus plantations in Brazil using the most commonly used models developed for natural forests. These forest plantations are valuable case studies, with large and numerous stands that are very uni- form, in which field measurements are precise compared to nat- ural forests. The height of planted Eucalyptus forest stands esti- mated from waveforms acquired by GLAS were compared with in situ measurements in order to determine the model that produced the best forest height estimates. For our slightly sloping study site (slope < 7°), the direct method defined as the difference between the signal begin and the ground peak provided forest height es- timates with an accuracy of 2.2 m. The use of statistical models based on waveform metrics and digital elevation models provided slightly better results (1.89 m accuracy) in comparison with the di- rect method and the most relevant metrics proved to be the trailing edge extent and the waveform extent. Moreover, a power law model was used to fit in situ aboveground biomass to in situ forest height. The results using this model with GLAS-derived heights showed an accuracy for biomass of 16.1 Mg/ha

Error Propagation Analysis for Remotely Sensed Aboveground Biomass

Edited version available. Full version will remain embargoed due to copyright. AS DCAbstract Above-Ground Biomass (AGB) assessment using remote sensing has been an active area of research since the 1970s. However, improvements in the reported accuracy of wide scale studies remain relatively small. Therefore, there is a need to improve error analysis to answer the question: Why is AGB assessment accuracy still under doubt? This project aimed to develop and implement a systematic quantitative methodology to analyse the uncertainty of remotely sensed AGB, including all perceptible error types and reducing the associated costs and computational effort required in comparison to conventional methods. An accuracy prediction tool was designed based on previous study inputs and their outcome accuracy. The methodology used included training a neural network tool to emulate human decision making for the optimal trade-off between cost and accuracy for forest biomass surveys. The training samples were based on outputs from a number of previous biomass surveys, including 64 optical data based studies, 62 Lidar data based studies, 100 Radar data based studies, and 50 combined data studies. The tool showed promising convergent results of medium production ability. However, it might take many years until enough studies will be published to provide sufficient samples for accurate predictions. To provide field data for the next steps, 38 plots within six sites were scanned with a Leica ScanStation P20 terrestrial laser scanner. The Terrestrial Laser Scanning (TLS) data analysis used existing techniques such as 3D voxels and applied allometric equations, alongside exploring new features such as non-plane voxel layers, parent-child relationships between layers and skeletonising tree branches to speed up the overall processing time. The results were two maps for each plot, a tree trunk map and branch map. An error analysis tool was designed to work on three stages. Stage 1 uses a Taylor method to propagate errors from remote sensing data for the products that were used as direct inputs to the biomass assessment process. Stage 2 applies a Monte Carlo method to propagate errors from the direct remote sensing and field inputs to the mathematical model. Stage 3 includes generating an error estimation model that is trained based on the error behaviour of the training samples. The tool was applied to four biomass assessment scenarios, and the results show that the relative error of AGB represented by the RMSE of the model fitting was high (20-35% of the AGB) in spite of the relatively high correlation coefficients. About 65% of the RMSE is due to the remote sensing and field data errors, with the remaining 35% due to the ill-defined relationship between the remote sensing data and AGB. The error component that has the largest influence was the remote sensing error (50-60% of the propagated error), with both the spatial and spectral error components having a clear influence on the total error. The influence of field data errors was close to the remote sensing data errors (40-50% of the propagated error) and its spatial and non-spatial Overall, the study successfully traced the errors and applied certainty-scenarios using the software tool designed for this purpose. The applied novel approach allowed for a relatively fast solution when mapping errors outside the fieldwork areas.HCED iraq, Middle Technical Universit

Remote sensing environmental change in southern African savannahs : a case study of Namibia

Savannah biomes cover a fifth of Earth’s surface, harbour many of the world’s most iconic species and most of its livestock and rangeland, while sustaining the livelihoods of an important proportion of its human population. They provide essential ecosystem services and functions, ranging from forest, grazing and water resources, to global climate regulation and carbon sequestration. However, savannahs are highly sensitive to human activities and climate change. Across sub-Saharan Africa, climatic shifts, destructive wars and increasing anthropogenic disturbances in the form of agricultural intensification and urbanization, have resulted in widespread land degradation and loss of ecosystem services. Yet, these threatened ecosystems are some of the least studied or protected, and hence should be given high conservation priority. Importantly, the scale of land degradation has not been fully explored, thereby comprising an important knowledge gap in our understanding of ecosystem services and processes, and effectively impeding conservation and management of these biodiversity hotspots. The primary drivers of land degradation include deforestation, triggered by the increasing need for urban and arable land, and concurrently, shrub encroachment, a process in which the herbaceous layer, a defining characteristic of savannahs, is replaced with hardy shrubs. These processes have significant repercussions on ecosystem service provision, both locally and globally, although the extents, drivers and impacts of either remain poorly quantified and understood. Additionally, regional aridification anticipated under climate change, will lead to important shifts in vegetation composition, amplified warming and reduced carbon sequestration. Together with a growing human population, these processes are expected to compound the risk of land degradation, thus further impacting key ecosystem services. Namibia is undergoing significant environmental and socio-economic changes. The most pervasive change processes affecting its savannahs are deforestation, degradation and shrub encroachment. Yet, the extent and drivers of such change processes are not comprehensively quantified, nor are the implications for rural livelihoods, sustainable land management, the carbon cycle, climate and conservation fully explored. This is partly due to the complexities of mapping vegetation changes with satellite data in savannahs. They are naturally spatially and temporally variable owing to erratic rainfall, divergent plant functional type phenologies and extensive anthropogenic impacts such as fire and grazing. Accordingly, this thesis aims to (i) quantify distinct vegetation change processes across Namibia, and (ii) develop methodologies to overcome limitations inherent in savannah mapping. Multi-sensor satellite data spanning a range of spatial, temporal and spectral resolutions are integrated with field datasets to achieve these aims, which are addressed in four journal articles. Chapters 1 and 2 are introductory. Chapter 3 exploits the Landsat archive to track changes in land cover classes over five decades throughout the Namibian Kalahari woodlands. The approach addresses issues implicit in change detection of savannahs by capturing the distinct phenological phases of woody vegetation and integrating multi-sensor, multi-source data. Vegetation extent was found to have decreased due to urbanization and small-scale arable farming. An assessment of the limitations leads to Chapter 4, which elaborates on the previous chapter by quantifying aboveground biomass changes associated with deforestation and shrub encroachment. The approach centres on fusing multiple satellite datasets, each acting as a proxy for distinct vegetation properties, with calibration/validation data consisting of concurrent field and LiDAR measurements. Biomass losses predominate, demonstrating the contribution of land management to ecosystem carbon changes. To identify whether biomass is declining across the country, Chapter 5 focuses on regional, moderate spatial resolution time-series analyses. Phenological metrics extracted from MODIS data are used to model observed fractional woody vegetation cover, a proxy for biomass. Trends in modelled fractional woody cover are then evaluated in relation to the predominant land-uses and precipitation. Negative trends slightly outweighed positive trends, with decreases arising largely in protected, urban and communal areas. Since precipitation is a fundamental control on vegetation, Chapter 6 investigates its relation to NDVI, by assessing to what extent observed trends in vegetation cover are driven by rainfall. NDVI is modelled as a function of precipitation, with residuals assumed to describe the fraction of NDVI not explained by rainfall. Mean annual rainfall and rainfall amplitude show a positive trend, although extensive “greening” is unrelated to rainfall. NDVI amplitude, used as a proxy for vegetation density, indicates a widespread shift to a denser condition. In Chapter 7, trend analysis is applied to a Landsat time-series to overcome spatial and temporal limitations characteristic of the previous approaches. Results, together with those of the previous chapters, are synthesized and a synopsis of the main findings is presented. Vegetation loss is predominantly caused by demand for urban and arable land. Greening trends are attributed to shrub encroachment and to a lesser extent conservation laws, agroforestry and rangeland management, with precipitation presenting little influence. Despite prevalent greening, degradation processes associated with shrub encroachment, including soil erosion, are likely to be widespread. Deforestation occurs locally while shrub encroachment occurs regionally. This thesis successfully integrates multi-source data to map, measure and monitor distinct change processes across scales